基于EST-SSR标记的贵州野生刺梨居群遗传多样性分析

删除或更新信息，请邮件至freekaoyan#163.com(#换成@)

本站小编 Free考研考试/2021-12-26

张怀山, 鄢秀芹, 鲁敏, 王道平, 安华明. 基于EST-SSR标记的贵州野生刺梨居群遗传多样性分析[J]. , 2017, 50(6): 1098-1108 https://doi.org/10.3864/j.issn.0578-1752.2017.06.011
ZHANG HuaiShan, YAN XiuQin, LU Min, WANG DaoPing, AN HuaMing. Analysis of the Genetic Diversity of Wild Rosa roxburghii Populations in Guizhou Province Based on EST-SSR Marker[J]. Scientia Acricultura Sinica, 2017, 50(6): 1098-1108 https://doi.org/10.3864/j.issn.0578-1752.2017.06.011

0 引言

【研究意义】刺梨（Rosa roxburghii Tratt,2n=2x=14）为蔷薇科（Rosaceae）蔷薇属（Rosa）植物,是中国特有的果树资源。因其果实含有丰富的营养物质和药用保健成分,尤其是极高的维生素C而倍受关注。刺梨作为贵州省重点发展的特色产业,种植面积已超过300 km²。同时,贵州省还含有丰富的野生刺梨资源,多分布于海拔1 000—1 600 m的山区、丘陵地带^[1-2]。这些丰富的野生资源无疑是刺梨品种选育和产业可持续发展的重要基础和必要前提。然而,近年来对野生刺梨资源进行调查收集时发现,贵州野生刺梨分布面积和群体数量正在急剧减少,这使得开展刺梨资源的收集、遗传多样性评价、保护性开发等方面的工作势在必行。【前人研究进展】遗传多样性,包括群体遗传结构的研究,有利于对种质资源提出科学合理的保护策略,尤其是对于现今正遭受人类直接干扰和破坏的物种。SSR标记是研究生物遗传多样性的重要方法,已在蔷薇科植物群体遗传学、种质资源、分类学与种系发生学等方面得到广泛应用。ZHANG等^[3]利用8对SSR标记对中国新疆伊犁地区4个种下居群109个新疆野苹果实生株系,进行群体遗传结构研究,结果表明巩留县遗传多样性最为丰富,故在制定原位种质保护计划时应优先考虑巩留县居群。ZONG等^[4]利用14对SSR标记分析了川梨的遗传多样性和遗传分化,这14个位点在幼苗种群中显示较高的多态性,并确定出3个群体将成为重点保护和调查使用的对象。LIU等^[5]利用14对基因组SSR引物研究了中国浙江省8个居群共77份豆梨资源的遗传多样性和种群结构,结果显示地理距离是目前豆梨遗传结构形成的一个关键因素,并确定FY群体是最适合进行原地保护的种群。陈娇等^[6]采用10对SSR引物对四川野生中国樱桃5个居群共133株的遗传多样性水平及居群遗传结构进行了研究,并提出野生中国樱桃的保护利用策略。MENG等^[7]利用SSR标记及单拷贝核基因对云南特有物种大花香水月季的遗传结构进行了研究,为保护这一濒危物种提供了重要的规划和决策。然而,刺梨种质资源分子水平的研究起步较晚,迄今为止,相关研究仅采用了少数几个刺梨材料作为参比试材研究蔷薇属不同种间亲缘进化的关系^[8-13],以及开展了刺梨种内遗传多样性研究,但仅局限于30个刺梨样品的RAPD分析^[14-15]。此外,笔者课题组在开发EST-SSR标记的过程中,曾涉及16份刺梨资源作为研究材料^[16-17]。由此可知,刺梨种质资源在分子水平上的遗传多样性研究基础还相当薄弱,表现为样品数量不足、代表性差,且缺乏群体遗传结构的研究。【本研究切入点】SSR标记根据其来源可分为基因组SSR和EST-SSR,目前对刺梨遗传背景的研究刚刚起步^[18],基因组信息较缺乏,可用的遗传标记较少,EST-SSR作为一种新型的分子标记,不仅具有传统SSR标记多态性高、共显性与重复性好等特点,更重要的是开发成本低;而且由于EST-SSR来源于表达的基因组区域,可直接反映相关基因的多样性,在不同物种间也具有良好的通用性^[19]。本研究基于课题组前期对刺梨果实转录组测序获得的数据,筛选有效的刺梨EST-SSR引物,以研究刺梨种质资源的遗传多样性和群体遗传结构。【拟解决的关键问题】从批量开发并合成的102对引物中,筛选出多态性最好的10对引物用于研究;基于EST-SSR标记技术,分析12个居群共255份贵州野生刺梨的遗传多样性及其居群遗传结构,以期为刺梨种质资源的科学保护和可持续利用提供理论依据。

1 材料与方法

试验于2014—2016年在贵州省果树工程技术研究中心进行。

1.1 试验材料

选择刺梨资源分布最为丰富的12个县作为居群代表采样点（表1、图1）,根据现有分布规模决定居群样本数量,每个居群取样12—38个不等,共255株,居群内单个样本间距离至少50 m以上,每个采样株均进行GPS定位,并详细记录各居群经度、纬度、海拔等地理条件（表1）。每株采取幼嫩叶片10—20 g,将采样叶片就地迅速保存于自封袋,置于冰盒,带回实验室用液氮处理后存于-70℃超低温冰箱,用于DNA提取。
Table 1
表1
表1贵州12个野生刺梨居群采样点概况
Table 1Description of the 12 sampling sites of wild R. roxburghii in Guizhou province

居群代号 Population code	居群来源地 Sampling localities	样本数量 Sample size	海拔 Altitude (m)	纬度 Latitude	经度 Longitude
SQ	石阡县 Shiqian county	18	657—892	27°28.798′—27°35.507′	108°01.646′—108°12.007′
FQ	福泉市 Fuquan city	12	957—1000	26°48.027′—26°48.793′	107°34.976′—107°36.185′
HS	惠水县 Huishui county	14	946—1032	25°44.113′—26°04.586′	106°34.440′—106°59.407′
LD	罗甸县 Luodian county	13	836—934	25°37.137′—25°43.969′	106°32.206′—106°40.331′
PX	盘县 Pan county	38	1445—1870	25°10.433′—25°54.893′	104°32.206′—104°50.843′
XY	兴义市 Xingyi city	28	1528—1920	24°54.607′—25°57.106′	104°43.413′—104°50.288′
AL	安龙县 Anlong county	20	1258—1334	25°11.480′—25°14.121′	105°19.107′—105°29.241′
ZY	遵义市 Zunyi city	19	778—981	27°26.875′—27°68.885′	106°54.649′—107°18.528′
SY	绥阳县 Suiyang county	15	734—1405	27°93.375′—29°94.346′	106°58.852′—107°19.976′
HX	花溪区 Huaxi district	35	1000—1500	26°17.367′—26°41.885′	106°42.898′—106°67.528′
PB	平坝县 Pingba county	30	1243—1400	26°21.673′—26°41.524′	106°12.623′—106°33.598′
QL	晴隆县 Qinglong county	13	1499—1583	25°81.843′—25°83.478′	105°19.752′—105°21.905′

新窗口打开

显示原图|下载原图ZIP|生成PPT
图1贵州省12个野生刺梨居群分布情况
-->Fig. 1Description of the 12 sampling sites distribution of wild R. roxburghii in Guizhou province
-->

1.2 DNA提取和SSR-PCR

1.2.1 DNA的提取试验材料基因组DNA提取参照POREBSKI等^[20]的CTAB法。
1.2.2 引物筛选从本课题组前期已完成的批量开发的EST-SSR引物中,随机选取102对交由上海英骏生物技术公司合成。选取栽培种‘贵农5号’（Rosa roxburghii Tratt cv. GuiNong 5）、2份野生刺梨资源进行初步引物筛选。选取‘贵农5号’、15份不同来源的野生刺梨资源对筛选出的有扩增条带的引物进行多态性检测,最终确定出10对多态性较好的引物用于本研究（表2）。
Table 2
表2
表2用于本研究的EST-SSR引物信息
Table 2Polymorphic SSR primers used for the analyses of all 255 R. roxburghii samples

引物编号 Primer No.	引物序列(5′→3′) Primer sequences(5′→3′)	重复单元 Repeat motif	退火温度 Anneal temperature (℃)	预期产物 Expected size (bp)	扩增带数 Number of bands	多态性带数 Number of polymorphic bands	多态率 Percentage of Polymorphism (%)
04	F TTGTGGTTATAGTTCAGCCCCTA R CAATCACGGAAATCATACATTCA	CCT (3×5)	60	95	3	2	66.67
07	F AATATACACGAACAACAACCATCG R GGAACATGACCCCTTTTCTTATT	CTACA (5×4)	60	147	8	8	100.00
14	F GTGGATGTGTCAAATCTAATGGC R GGAGGGAAGAAGTAGTGAAGAACA	CT (2×6)	60	152	2	1	50.00
17	F GATGCTTTTCATTCTGCTTCAAC R ATTTTTACCGTACTCTGGGTGCT	TC (2×6)	60	143	3	1	33.33
19	F ATTTGTTTTTCGTTTTTCTTCCC R TTTGGTCATTCATTCTTCTCCTC	CAT (3×5)	60	128	2	1	50.00
24	F GAGTCATGTTGAATGATATTGGC R TTTCCTCTTTTCTTCTTTTTCCC	TGGA (4×13)	59	131	2	2	100.00
56	F AGCAATAGAGAACTAGAATATCAGAGA R TTTTCCCCCTTCAGTTGTAAGAT	GT (2×10)	58	116	3	2	66.67
58	F TGAGAACCAAAAGAAACGACTACA R GTCACAAATTCACGAGCCATATT	AT (2×8)	60	146	3	3	100.00
65	F GAGAGAGTAGCCTTTGATGAGGA R GAATGTCGTCAGGGAGCAGTAG	CTC (3×7)	60	146	4	4	100.00
70	F TGATTGGTAGTGGAAGCGATATT R ACACATCAGTCCTTCAACAGTCC	TCA (3×5)	60	108	3	3	100.00
	平均 Average				3.3	2.7	76.67

新窗口打开
1.2.3 PCR扩增 PCR反应体系为20 µL,其中包括
2×Mix 10 µL;10 µmol·L^-1的Primer-F、Primer-R各0.5 µL;ddH₂O 8 µL;50 ng·µL^-1的DNA模板1 µL。扩增程序为94℃下预变性3 min;然后进行35个循环,每个循环包括94℃变性40 s,55℃退火40 s（退火温度因不同引物而异）,72℃延伸1 min;最后72℃延伸10 min。
1.2.4 电泳检测扩增产物利用8%非变性聚丙烯酰胺凝胶,于DYCZ-30C型垂直电泳槽中电泳分离,150 V电压下电泳90 min。电泳后参照BASSAM等^[21]的方法进行银染显色,并在BIO-RAD凝胶成像系统中拍照记录、分析。

1.3 数据分析处理

根据分子量大小对扩增结果读带,以二倍体形式记录,从大到小依次记为A、B、C…^[22]。应用POPGENE 1.31软件^[23]对各居群在多基因座位上的Hardy-Weinberg平衡进行检测,计算观测等位基因数（N_a）、有效等位基因数（N_e）、多态性位点比率（P）、Nei’s基因多样性指数（Nei’s）、Shannon信息指数（I）,观测杂合度（H_o）和期望杂合度（H_e）、基因流（N_em）、居群内近交系数（F_it）、总近交系数（F_is）、居群间分化系数（F_st）,Nei’s标准遗传距离（GD）和遗传一致度（GI）。根据Nei’s标准遗传距离,利用NTSYS-pc2.10e 软件中的UPGMA方法对各居群进行聚类分析,并进行Nei’s遗传距离与地理距离的Mantel 相关性检验。

2 结果

2.1 EST-SSR标记的多态性

以‘贵农5号’和2份野生刺梨材料的基因组DNA为模板,对随机选取并合成的102对引物进行PCR扩增、筛选。结果表明,其中的71对引物产生理想的PCR产物,有效扩增率为69.61%,通过筛选确定为有效的EST-SSR引物。进一步选取‘贵农5号’和15份野生刺梨材料对71对有效引物进行扩增、多态性评价。结果31对引物呈现出多态性,占有效引物的43.66%。从31对具有多态性的引物中选取10对效果良好、结果稳定的引物（表2）,对255份刺梨材料进行扩增,进而分析刺梨种质资源的遗传多样性。10对引物在255份样品上均能扩增出清晰的条带,共扩增出33个条带,不同引物的扩增条带在2—8条之间,平均3.3条,其中有27条为多态性条带,多态率为76.67%（表2）。图2为引物24在44份材料中的扩增情况。

显示原图|下载原图ZIP|生成PPT
图2引物24对刺梨基因组DNA的扩增图谱 M：DNA标准分子量
-->Fig. 2SSR amplification pattern of R. roxburghii genomic DNA by primer 24 DNA molecular weight marker
-->

2.2 居群遗传多样性

在居群水平上,10对EST-SSR引物在12个自然居群中共扩增出33个等位基因,平均每个微卫星位点扩增3.3个（N_a=3.3）。除晴隆（QL）居群在位点17、19上,福泉（FQ）居群在位点19上只扩增出单一条带,未表现出多态性,10个居群的多态位点百分率P为100%（表 3）,表明贵州野生刺梨具有较丰富的遗传多样性。各居群观测等位基因数N_a为2.20—3.70;有效等位基因数N_e为1.84—2.29;Shannon信息指数I为0.62—0.88,观测杂合度H_o为0.64—0.76;预期杂合度 H_e为0.44—0.54,且H_o极显著高于H_e（P=0.0001）;惠水（HS）和花溪（HX）居群拥有最高的基因多样性指数（Nei’s为0.52）,晴隆（QL）居群的基因多样性指数最低,为0.40（表3）。
Table 3
表3
表3贵州12个刺梨居群在10个微卫星位点的多样性指数
Table 3The genetic diversity index of 12 populations of R. roxburghii at 10 microsatellite loci

位点 Locus	遗传参数 Parameter	居群代号 Population code												平均 Mean
位点 Locus	遗传参数 Parameter	SQ	FQ	HS	LD	PX	XY	AL	ZY	SY	HX	PB	QL	平均 Mean
04	N_a	4.00	3.00	3.00	2.00	3.00	4.00	3.00	4.00	2.00	4.00	4.00	2.00	3.17
	N_e	2.48	2.17	2.14	2.00	2.10	2.51	2.55	2.23	2.00	2.40	2.21	2.00	2.23
	I	1.05	0.84	0.82	0.69	0.80	1.06	1.00	0.91	0.69	1.00	0.89	0.69	0.87
	H_o	1.00	1.00	1.00	1.00	1.00	1.00	1.00	1.00	1.00	1.00	1.00	1.00	1.00
	H_e	0.62	0.56	0.55	0.54	0.53	0.61	0.62	0.57	0.52	0.59	0.56	0.56	0.57
	Nei’s	0.60	0.54	0.53	0.50	0.52	0.60	0.61	0.55	0.50	0.58	0.55	0.50	0.55
07	N_a	3.00	3.00	6.00	4.00	7.00	6.00	6.00	5.00	3.00	8.00	6.00	3.00	5.00
	N_e	2.51	2.46	3.50	3.27	1.93	3.44	3.54	2.25	2.97	3.23	3.10	1.85	2.84
	I	1.00	0.99	1.44	1.28	1.02	1.40	1.47	1.05	1.09	1.43	1.35	0.80	1.19
	H_o	0.59	0.58	0.64	0.57	0.58	0.64	0.74	0.67	0.57	0.74	0.83	0.60	0.65
	H_e	0.62	0.62	0.74	0.75	0.49	0.72	0.74	0.57	0.69	0.70	0.69	0.51	0.65
	Nei’s	0.60	0.59	0.71	0.69	0.48	0.71	0.72	0.56	0.66	0.69	0.68	0.46	0.63
14	N_a	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00
	N_e	1.78	1.88	1.96	1.96	2.00	1.88	1.95	1.98	1.91	1.99	2.00	2.00	1.94
	I	0.63	0.66	0.68	0.68	0.69	0.66	0.68	0.69	0.67	0.69	0.69	0.69	0.68
	H_o	0.65	0.75	0.86	0.86	0.97	0.75	0.84	0.89	0.79	0.91	1.00	1.00	0.86
	H_e	0.45	0.49	0.51	0.53	0.51	0.48	0.50	0.51	0.49	0.50	0.51	0.56	0.50
	Nei’s	0.44	0.47	0.49	0.49	0.50	0.47	0.49	0.49	0.48	0.50	0.50	0.50	0.49
17	N_a	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	1.00	1.92
	N_e	1.71	1.38	1.69	1.32	1.63	1.11	1.23	1.60	1.51	1.29	1.36	1.00	1.40
	I	0.61	0.45	0.60	0.41	0.58	0.21	0.34	0.56	0.52	0.38	0.43	0.00	0.42
	H_o	0.59	0.33	0.57	0.29	0.53	0.11	0.21	0.50	0.43	0.26	0.31	0.00	0.34
	H_e	0.43	0.29	0.42	0.26	0.39	0.10	0.19	0.39	0.35	0.23	0.27	0.00	0.28
	Nei’s	0.42	0.28	0.41	0.24	0.39	0.10	0.19	0.38	0.34	0.22	0.26	0.00	0.27
19	N_a	2.00	1.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	1.00	1.83
	N_e	1.26	1.00	1.41	1.51	1.11	1.11	1.17	1.53	1.07	1.47	1.49	1.00	1.26
	I	0.36	0.00	0.47	0.52	0.21	0.21	0.28	0.53	0.15	0.50	0.51	0.00	0.31
	H_o	0.24	0.00	0.36	0.43	0.11	0.11	0.16	0.44	0.07	0.40	0.41	0.00	0.23
	H_e	0.21	0.00	0.30	0.36	0.10	0.10	0.15	0.36	0.07	0.32	0.33	0.00	0.19
	Nei’s	0.21	0.00	0.29	0.34	0.10	0.10	0.15	0.35	0.07	0.32	0.33	0.00	0.19
24	N_a	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00	2.00
	N_e	1.90	1.99	1.69	2.00	1.76	1.91	1.82	2.00	1.42	1.88	1.75	1.72	1.82
	I	0.67	0.69	0.60	0.69	0.62	0.67	0.64	0.69	0.47	0.66	0.62	0.61	0.64
	H_o	0.41	0.42	0.29	0.71	0.58	0.50	0.37	0.78	0.36	0.69	0.48	0.60	0.52
	H_e	0.49	0.52	0.42	0.54	0.44	0.49	0.46	0.51	0.30	0.47	0.44	0.47	0.46
	Nei’s	0.47	0.50	0.41	0.50	0.43	0.48	0.45	0.50	0.29	0.47	0.43	0.42	0.45
56	N_a	4.00	4.00	3.00	3.00	3.00	4.00	4.00	4.00	3.00	4.00	3.00	3.00	3.50
	N_e	2.71	2.80	2.67	2.51	2.50	2.62	2.81	2.40	2.35	2.70	2.64	2.38	2.59
	I	1.12	1.15	1.04	0.99	1.00	1.10	1.16	1.01	0.94	1.09	1.03	0.94	1.05
	H_o	1.00	1.00	1.00	1.00	0.92	0.96	1.00	0.94	0.93	1.00	1.00	1.00	0.98
	H_e	0.65	0.67	0.65	0.65	0.61	0.63	0.66	0.60	0.60	0.64	0.63	0.64	0.64
	Nei’s	0.63	0.64	0.62	0.60	0.60	0.62	0.64	0.58	0.57	0.63	0.62	0.58	0.61
位点 Locus	遗传参数 Parameter	居群代号 Population code												平均 Mean
位点 Locus	遗传参数 Parameter	SQ	FQ	HS	LD	PX	XY	AL	ZY	SY	HX	PB	QL	平均 Mean
58	N_a	4.00	3.00	4.00	2.00	4.00	3.00	4.00	3.00	3.00	4.00	4.00	3.00	3.42
	N_e	2.18	2.38	2.56	1.96	2.50	2.26	2.53	2.31	2.43	2.19	2.39	2.63	2.36
	I	0.98	0.96	1.09	0.68	1.09	0.89	1.07	0.94	0.98	0.91	1.03	1.03	0.97
	H_o	0.76	0.92	0.93	0.86	0.82	0.96	1.00	0.89	0.93	0.86	0.90	1.00	0.90
	H_e	0.56	0.61	0.63	0.53	0.61	0.57	0.62	0.58	0.61	0.55	0.59	0.69	0.60
	Nei’s	0.54	0.58	0.61	0.49	0.60	0.56	0.61	0.57	0.59	0.54	0.58	0.62	0.57
65	N_a	3.00	2.00	5.00	2.00	3.00	3.00	4.00	4.00	4.00	5.00	3.00	2.00	3.33
	N_e	2.03	1.70	2.14	1.69	1.77	2.47	2.13	1.91	1.56	2.41	1.82	1.47	1.93
	I	0.78	0.60	1.04	0.60	0.70	0.99	0.92	0.83	0.71	1.04	0.70	0.50	0.78
	H_o	0.71	0.58	0.50	0.57	0.55	0.79	0.58	0.56	0.36	0.74	0.48	0.40	0.57
	H_e	0.52	0.43	0.55	0.44	0.44	0.61	0.54	0.49	0.37	0.59	0.46	0.36	0.48
	Nei’s	0.51	0.41	0.53	0.41	0.44	0.60	0.53	0.48	0.36	0.59	0.45	0.32	0.47
70	N_a	3.00	4.00	3.00	4.00	4.00	3.00	4.00	4.00	4.00	4.00	4.00	3.00	3.67
	N_e	2.64	2.34	2.56	2.65	2.91	2.78	3.13	2.85	2.74	2.67	2.44	2.38	2.67
	I	1.03	0.98	1.00	1.10	1.16	1.06	1.19	1.14	1.13	1.11	1.00	0.94	1.07
	H_o	1.00	1.00	0.93	0.86	0.89	0.93	0.74	0.94	1.00	0.94	0.86	1.00	0.92
	H_e	0.64	0.60	0.63	0.67	0.66	0.65	0.70	0.67	0.66	0.63	0.60	0.64	0.65
	Nei’s	0.62	0.57	0.61	0.62	0.66	0.64	0.68	0.65	0.64	0.62	0.59	0.58	0.62
Mean	N_a	2.90	2.60	3.30	2.50	3.20	3.10	3.30	3.30	2.70	3.70	3.30	2.20	3.30
	N_e	2.12	2.01	2.23	2.09	2.02	2.21	2.29	2.11	2.00	2.22	2.12	1.84	2.11
	I	0.82	0.73	0.88	0.77	0.79	0.82	0.88	0.84	0.74	0.88	0.83	0.62	0.80
	H_o	0.69	0.66	0.71	0.71	0.69	0.67	0.66	0.76	0.64	0.75	0.73	0.66	0.69
	H_e	0.52	0.48	0.54	0.53	0.48	0.50	0.52	0.52	0.47	0.52	0.51	0.44	0.50
	Nei’s	0.50	0.46	0.52	0.49	0.47	0.49	0.51	0.51	0.45	0.52	0.50	0.40	0.49

N_a: Observed number of alleles; N_e: Effective number of alleles; I: Shannon’s information index; Nei’s: Nei’s gene diversity; H_o: Observed heterozygosity; H_e: Expected heterozygosity; P: Percentage of polymorphic lociN_a：观测等位基因数;N_e：有效等位基因数;I：Shannon信息指数;Nei’s：基因多样性指数;H_o：观测杂合度;H_e：预期杂合度;P：多态位点百分率

新窗口打开

2.3 居群遗传结构

采用卡方检验对各居群的多基因座位进行检测发现,除罗甸（LD）居群外,11个刺梨居群均在1/3以上的位点显著偏离Hardy-Weinberg平衡（表4）,在04位点甚至全部显著偏离。12个刺梨居群的基因多样性指数Nei’s均大于0.4;Shannon信息指数I除晴隆（QL）外均大于0.70,最高为0.88（表3）。以上结果表明,12个刺梨居群内具有较高的遗传多样性。
Table 4
表4
表4Hardy-Weinberg平衡检测的P 值
Table 4The P values for Hardy-Weinberg Equilibrium test

位点 Locus	居群代号 Population code
位点 Locus	SQ	FQ	HS	LD	PX	XY	AL	ZY	SY	HX	PB	QL
04	0.01**	0.01**	0.00**	0.01**	0.00**	0.00**	0.00**	0.01**	0.00**	0.00**	0.00**	0.05*
07	0.00**	0.00**	0.13	0.08	0.00**	0.00**	0.01**	0.12	0.00**	0.01**	0.00**	0.90
14	0.06	0.05*	0.01**	0.07	0.00**	0.00**	0.00**	0.00**	0.02*	0.00**	0.00**	0.05*
17	0.11	0.56	0.17	0.76	0.03*	0.81	0.66	0.19	0.36	0.41	0.35	0.00**
19	0.64	0.00**	0.47	0.57	0.76	0.81	0.76	0.26	1.00	0.15	0.18	0.00**
24	0.51	0.48	0.20	0.35	0.04*	0.87	0.36	0.03*	0.47	0.01**	0.55	0.45
56	0.01**	0.09	0.00**	0.11	0.00**	0.00**	0.01**	0.04*	0.02*	0.00**	0.00**	0.26
58	0.43	0.05*	0.14	0.07	0.00**	0.00**	0.01**	0.01**	0.02*	0.00**	0.01**	0.26
65	0.36	0.19	0.00**	0.37	0.33	0.02*	0.53	0.97	0.40	0.37	0.50	0.71
70	0.00**	0.09	0.02*	0.68	0.00**	0.00**	0.01**	0.02*	0.04*	0.00**	0.00**	0.26

*presented that the locus was away from Hardy-Weinberg Equilibrium significantly (at 0.05 level);**presented that the locus was away from Hardy-Weinberg Equilibrium highly significantly (at 0.01 level)*为显著偏离Hardy-Weinberg平衡（0.05水平）,**为极显著偏离Hardy-Weinberg平衡（0.01水平）

新窗口打开
贵州12个刺梨居群在10个微卫星位点上的F-statistics分析结果见表5。12个刺梨居群在各位点上F_it、F_is为负值,F_st值在位点65最大,为0.0773;在位点4最小,为0.0092。F_st平均值为 0.0403,居群各位点基因流N_em值均大于2,且平均值为5.9484。
Table 5
表5
表5各位点固定指数及基因流
Table 5The fixed index and gene flow of 10 microsatellite loci

位点 Locus	居群内近交系数 F_it	总近交系数 F_is	居群间分化系数 F_st	基因流 N_em
4	-0.8055	-0.8223	0.0092	26.9562
7	-0.0206	0.0255	0.0450	5.3047
14	-0.7475	-0.7672	0.0112	22.1621
17	-0.2072	-0.2785	0.0558	4.2330
19	-0.1279	-0.2109	0.0686	3.3947
24	-0.0737	-0.1561	0.0713	3.2581
56	-0.5711	-0.5992	0.0176	13.9487
58	-0.5176	-0.5713	0.0342	7.0641
65	-0.1198	-0.2136	0.0773	2.9840
70	-0.4195	-0.4829	0.0428	5.5964
平均Mean	-0.3784	-0.4363	0.0403	5.9484

新窗口打开
贵州刺梨居群间GI为0.9068—0.9926,平均值0.9586,遗传一致度非常高（表6）;平均Nei’s遗传距离（GD）0.0431,花溪（HX）和平坝（PB）居群Nei’s遗传距离最近,仅为0.0074,兴义（XY）居群和晴隆（QL）居群之间的Nei’s遗传距离最大,也仅为0.0978。
Table 6
表6
表6居群间的Nei’s遗传距离（GD）和遗传一致度（GI）
Table 6Nei’s genetic distance and genetic identity between populations

	SQ	FQ	HS	LD	PX	XY	AL	ZY	SY	HX	PB	QL
SQ	****	0.9573	0.9378	0.9464	0.9443	0.9768	0.9502	0.9658	0.9134	0.9639	0.9453	0.9084
FQ	0.0436	****	0.9544	0.9519	0.9689	0.9572	0.9786	0.9733	0.9698	0.9637	0.9542	0.9405
HS	0.0643	0.0467	****	0.9624	0.9802	0.9213	0.9697	0.9656	0.9734	0.9758	0.9798	0.9465
LD	0.0551	0.0493	0.0384	****	0.9543	0.9445	0.9699	0.9644	0.9422	0.9822	0.9784	0.9386
PX	0.0573	0.0316	0.0200	0.0468	****	0.9246	0.9727	0.9807	0.9611	0.9736	0.9781	0.9712
XY	0.0234	0.0438	0.0819	0.0571	0.0784	****	0.9622	0.9478	0.9164	0.9634	0.9378	0.9068
AN	0.0511	0.0217	0.0308	0.0305	0.0276	0.0386	****	0.9672	0.9672	0.9836	0.9784	0.9532
ZY	0.0348	0.0270	0.0350	0.0363	0.0195	0.0536	0.0333	****	0.9511	0.9795	0.9780	0.9548
SY	0.0906	0.0307	0.0270	0.0595	0.0397	0.0873	0.0334	0.0501	****	0.9597	0.9631	0.9420
HX	0.0368	0.0370	0.0245	0.0595	0.0268	0.0373	0.0166	0.0208	0.0411	****	0.9926	0.9590
PB	0.0562	0.0469	0.0204	0.0219	0.0221	0.0642	0.0219	0.0222	0.0376	0.0074	****	0.9704
QL	0.0960	0.0613	0.0550	0.0633	0.0292	0.0978	0.0479	0.0463	0.0598	0.0418	0.0301	****

Nei’s genetic distance (below diagonal) and genetic identity (above diagonal)Nei’s遗传距离（下三角）、遗传一致度（上三角）

新窗口打开

2.4 居群聚类结果

基于贵州刺梨居群间Nei’s遗传距离^[24],采用UPGMA 法进行聚类分析,进一步直观分析居群间遗传关系,如聚类图（图3）所示,在Nei’s遗传距离0.06处可将12个贵州刺梨居群分为2个组。石阡（SQ）和兴义（XY）居群聚为第一组,其他居群聚为第二组,其中地理位置相近的花溪（HX）和平坝（PB）居群首先聚在一起,显示了最近的亲缘关系。对居群间的Nei’s遗传距离和地理距离之间的Mantel（r=0.2498,P=0.9512）检测也表明,在采集地范围,刺梨的遗传变异分布和地理位置具有显著的相关性,尤其当居群间地理距离小于500 km时,其相关关系如图4所示。

显示原图|下载原图ZIP|生成PPT
图3基于 Nei's 遗传距离的居群UPGMA 聚类图
-->Fig. 3UPGMA dendrogram based on Nei’s genetic distance
-->

显示原图|下载原图ZIP|生成PPT
图412个刺梨居群间Nei’s遗传距离与地理距离的相关关系
-->Fig. 4The correlation between Nei’s genetic distance and geographic distance for 12 populations of R. roxburghii
-->

3 讨论

3.1 刺梨遗传多样性

本研究利用10对EST-SSR引物,对采自贵州省12个自然居群的255份野生刺梨种质资源进行遗传多样性分析,10个微卫星位点在各居群中均表现出不同程度的多态性,除福泉（FQ）和晴隆（QL）居群外,各居群的多态位点百分率P 均为100%,表明这10个微卫星位点具有高度的多态性,可作为有效的遗传标记用于刺梨居群遗传多样性和遗传结构分析。
观测等位基因数（N_a）是衡量SSR位点多态性和居群变异程度高低的重要指标。10对EST-SSR引物在12个居群（255份野生刺梨种质资源）中共扩增出33个等位基因,平均每个微卫星位点扩增3.3个（N_a=3.3）,与MENG等^[7]采用7对SSR标记对濒危物种大花香水月季27个居群的研究结果接近（N_a=3.9）,低于孙萍等^[25]利用 10对SSR引物对清凉峰地区36份三叶海棠的扩增结果（N_a=7.1）,低于LIU等^[5]利用14对SSR引物对中国浙江省77份豆梨资源的扩增结果（N_a=9.5）,也低于陈娇等^[6]采用10对SSR引物对5个中国樱桃野生居群共133株的研究结果（N_a=7.8）,这可能是由于EST-SSR来源于表达的基因组区域,所以多态性比基因组SSR标记稍低^[26]。
平均有效等位基因数（N_e）、平均Shannon信息指数（I）、平均期望杂合度（H_e）平均基因多样性指数（Nei’s）等是评价遗传多样性的重要指标。本试验研究结果（N_e=2.19,I=0.87,H_o=0.70,H_e=0.51,Nei’s=0.51）,低于同属于蔷薇科的中国樱桃野生居群^[6]（Nei’s=0.70,I=1.53）,三叶海棠^[25]（H_e=0.70,I=1.46,N_e=3.95）;却高于野杏^[27]（Nei’s=0.29,I=0.46）和新疆野苹果^[3]（H_e=0.26,I=0.41,N_e=1.43）,表明相较于蔷薇科其他植物而言,贵州野生刺梨群体遗传多样性处于适中到较高水平。这可能与刺梨的生活习性（如繁殖方式等）有关,刺梨可自花授粉结实,也可异花授粉结实,因此其遗传多样性倾向于界于杂交物种和自交物种之间。

3.2 刺梨遗传分化

贵州野生刺梨居群间分化系数F_st平均值为0.0403,低于浙江野生豆梨^[5]（F_st=0.136）,云南野生大花香水月季^[7]（F_st=0.240）及云南野生川梨^[28]（F_st=0.404）,表明刺梨居群间属低度分化^[29],居群间的遗传变异只占总变异的较少部分（4.03%）;居群间基因流N_em=5.9484,表明居群间基因交流频繁,从而有效抑制了由遗传漂变而引起的遗传分化^[30]。同时,频繁的基因交流必然导致居群间遗传的均质化,使各居群间表现出较高的遗传一致度（GI>0.9）及较小的遗传距离（GD<0.0978）。本试验所收集材料多集中在道路两旁,人类采集活动频繁,可能是促进刺梨居群间基因交流的重要原因。此外,也有相当部分试验材料采自人迹罕至的山坡、河沟边,因此流水和采食刺梨果实的鸟兽,也可能是刺梨远距离传播基因交流的主要方式。
贵州野生刺梨12个居群在多个位点上偏离Hardy-Weinberg平衡,且观测杂合度高于期望杂合度,表明贵州野生刺梨居群内不仅遗传多样性高,还存在杂合度过剩的现象。推测可能主要存在两方面的原因：一方面,由于刺梨遗传背景较为狭窄,起源亲本数量有限,奠基者效应（Founder effect）会导致连锁不平衡现象,继而造成杂合子过剩;另一方面,可能是虫媒和风媒造成的高异花授粉率,导致刺梨野生自然居群内有性繁殖的自交率低（F_is、F_it<0）,而频繁的基因交流（N_em=5.9484）使种内杂交产生的杂合子可通过无性繁殖方式进行有效固定,从而导致杂合度过剩。对居群间的Nei’s遗传距离和地理距离之间的Mantel检测表明,刺梨的遗传变异分布和地理位置具有显著的相关性,暗示地理距离可能是影响刺梨遗传分化的主要因素。

3.3 刺梨保护策略

由于土地开垦以及城市化进程等导致的土地侵占、环境破坏以及生境破碎化等,加上农业耕作时的乱砍滥伐,使得野生刺梨分布面积和群体数量急剧减少。以前在田间地头、房前屋后随处可见的刺梨资源,现在必须要深入到荒山坡地才能觅其踪影。可以预见,随着城市化发展日新月异和生态环境的日趋恶化,刺梨种质资源大量丢失的现象将日趋严重。因此,刺梨资源的收集、评价、保护等方面的工作迫在眉睫。本试验结果表明,贵州野生刺梨12个居群中以花溪（HX）、惠水（HS）、平坝（PB）、安龙（AL）、遵义（ZY）等5个居群拥有高于平均值的观测等位基因数（N_a）、有效等位基因数（N_e）、Nei’s基因多样性指数（Nei’s）、Shannon信息指数（I）,观测杂合度（H_o）和期望杂合度（H_e）（表3）,遗传多样性较高,应当采取就地保存的策略进行种质保存,其中平坝（PB）和花溪（HX）居群遗传距离最小,且地理距离也较近,可优先选择遗传多样性更高的花溪（HX）居群进行保护;而对于人为破坏最为严重的石阡（SQ）、罗甸（LD）、盘县（PX）等3个居群则应采取迁地保存的策略。下一步将扩大野生资源的采集范围和SSR引物的筛选范围,并选择更多的分子标记技术,更加系统地开展不同区域间遗传多样性研究,为刺梨和其他蔷薇科植物的资源分类、分子标记辅助育种、遗传图谱构建等研究提供更多专一有效的分子标记。

4 结论

基于EST-SSR标记对贵州野生刺梨12个居群255份种质资源的遗传多样性分析表明,10对EST-SSR引物具有高度的多态性,可作为有效的遗传标记用于刺梨遗传多样性和遗传结构评价。贵州省野生刺梨自然居群内的遗传多样性较高,居群内存在杂合度过剩的现象,且绝大多数的遗传变异发生在居群内（居群间分化系数F_st平均值0.0403）,居群间具有基因交流频繁（基因流N_em平均值为5.9484）、遗传一致度高（GI为0.9068—0.9926）、Nei’s遗传距离小（最大遗传距离为0.0978）等特点。基于以上结果,应优先考虑对遗传多样性较高的花溪（HX）、惠水（HS）、安龙（AL）、遵义（ZY）等4个居群实施就地保存策略。
The authors have declared that no competing interests exist.