用分子标记揭示植物随机大居群中亚居群的遗传结构——茅栗自然居群空间遗传结构的SSR分析

doi:10.17521/cjpe.2006.0021

摘要/Abstract

摘要：

采用微卫星标记对茅栗(Castanea sequinii)随机大居群以及其中各亚居群的遗传结构进行了空间自相关分析,以探讨植物自然居群内遗传变异的分布特征及其形成机制。通过9对微卫星引物所产生的119个多态位点,测定了大别山区域内茅栗居群以及各亚居群的空间自相关系数Moran's I值。结果表明:大别山分布的野生茅栗为一个缺乏空间结构的随机大居群,茅栗亚居群之间频繁的花粉流削弱了地理隔离导致的遗传漂变或分化作用在维系居群随机遗传结构中具有的重要作用。但是,在接近亚居群大小的地域范围内(0.228 km)具有一定的空间结构,即小地域尺度中的亚居群存在着空间遗传结构。取样的3个亚居群在小格局范围内都存在一定的空间结构,遗传变异基本上呈非随机分布,在短距离内(61 m)3个亚居群一致表现出不同程度的显著正相关,而随着距离的增加,Moran's I值虽然在不同亚居群间存在一定差异变化,但是总体而言趋向预期值,即不存在空间结构,说明其遗传变异在亚居群内只是在短距离内形成一定的空间结构。研究认为有限的种子散播以及微生境选择等因素可能是产生这种小格局的遗传结构的主要原因。上述研究结果有助于进一步了解植物随机大居群的进化历史和生态过程,同时也为栗属植物中国特有种的保育策略提供了科学依据。

关键词: 茅栗, 微卫星, 随机大居群, 亚居群, 居群遗传结构, 空间自相关, 遗传多样性, 基因流

Abstract:

Castanea squinii, an endemic tree widely distributed in China, has great value both for chestnut breeding and in natural forest ecosystems. The spatial genetic structure within and among populations is an important part of the evolutionary and ecological genetic dynamics of natural populations and can provide insights into effective conservation of genetic resources. In the present study, the spatial genetic structure of a natural population of C. sequinii in the Dabie Mountain region, including three subpopulations, was investigated using SSR markers. Nine prescreened SSR loci generated 29-33 alleles and were used for spatial autocorrelation analysis. Based on Moran's I coefficient, a panmictic population of C. sequinii was revealed that lacked spatial genetic structure in the Dabie Mountain region. These results suggested that high pollen flow among subpopulations counteracted genetic drift or genetic differentiation and played an important role in maintaining a random and panmictic population structure in the region. However, spatial genetic structure was detected in the geographic range of each subpopulation's scale (0.228 km) indicating that spatial genetic structure occurred at a fine scale within each subpopulation. All three subpopulations showed significant fine-scale structure. The genetic variation was not randomly distributed within 61 m in each subpopulation (Moran's I positive values). However, although Moran's I values varied among the different subpopulations, Moran's I in all the three subpopulations reached the expected values with an increase in distance suggesting a general patchy distribution in the subpopulations. The fine-scale structure seemed to be influenced by restricted seed dispersal and microenvironment selection. These results have important implications for understanding the evolutionary history and ecological process of natural populations of C. squinii and provide baseline data for formulating conservation strategies of Castanea species.

Key words: Castanea sequinii, SSR, Panmictic large population, Subpopulation, Genetic structure, Spatial autocorrelation, Genetic diversity, Gene flow

王英, 康明, 黄宏文. 用分子标记揭示植物随机大居群中亚居群的遗传结构——茅栗自然居群空间遗传结构的SSR分析. 植物生态学报, 2006, 30(1): 147-156. DOI: 10.17521/cjpe.2006.0021

WANG Ying, KANG Ming, HUANG Hong-Wen. SUBPOPULATION GENETIC STRUCTURE IN A PANMICTIC POPULATION AS REVEALED BY MOLECULAR MARKERS: A CASE STUDY OF CASTANEA SEQUINII USING SSR MARKERS. Chinese Journal of Plant Ecology, 2006, 30(1): 147-156. DOI: 10.17521/cjpe.2006.0021

导出引用管理器 EndNote|Ris|BibTeX

链接本文: https://www.plant-ecology.com/CN/10.17521/cjpe.2006.0021

https://www.plant-ecology.com/CN/Y2006/V30/I1/147

图/表 7

参考文献 42

[1]	Asuka Y, Tomaru N, Nisimura N, Tsumura Y, Yamamoto S (2004). Heterogeneous genetic structure in a Fagus crenata population in an old-growth beech forest revealed by microsatellite markers. Molecular Ecology, 13,1241-1250. URL PMID
[2]	Ausubel FM, Brent R, Kingston RE, Moore DD, Seidman JG, Smith JA, Stuhl K (1987). Current Protocols in Molecular Biology. John Wiley & Sons, New York.
[3]	Bjornstad ON, Falck W, Stenseth NC (1995). A geographic gradient in small rodent density fluctuations: a statistical modeling approach. Proceedings of the Royal Society of London B, 262,127-133.
[4]	Chung MG (1996). Spatial genetic structure among Korean populations of Hosta minor and H. capitata (Liliaceae). Botanical Bulletin of Academia Sinica, 37,25-30.
[5]	Chung MG, Epperson BK (2000). Spatial genetic structure of allozyme polymorphisms in a population of Eurya japonica (Theaceae). Silvae Genetica, 49,1-4.
[6]	Chung MY, Nason J, Chung MG, Kim KJ, Park CW, Sun BY, Pak JH (2002). Landscape-level spatial genetic structure in Quercus acutissima (Fagaceae). American Journal of Botany, 89,1229-1236.
[7]	Clegg MT, Brown AHD (1983). The founding of plant populations. In: Schonewald-Cox CM, Chambers SM, MacBryde B, Thomas WL eds. Genetics and Conservation. Benjamin-Cummings, Menlo Park, CA,216-228.
[8]	Cottrell JC, Munro RC, Tabbener HE, Milner AD, Forrest GI, Lowe AJ (2003). Comparison of fine-scale genetic structure using nuclear microsatellites within two British oakwoods differing in population history. Forest Ecology and Management, 176,287-303.
[9]	Ellstrand NC, Elam DR (1993). Population genetic consequences of small population size: implications for plant conservation. Annual Review Ecology and Systematics, 34,217-242.
[10]	Epperson BK, Allard RW (1989). Spatial structure of two-locus genotypes under isolation by distance. Genetics, 140,365-375. URL PMID
[11]	Epperson BK (1990). Spatial patterns of genetic variation within plant populations. In: Brown AHD, Clegg MT, Kahler AL, Weir BS eds. Plant Population Genetics,Breeding and Genetic Resources. Sinauer Associates, MA,229-253.
[12]	Epperson BK (1992). Spatial structure of genetic variation within populations of forest trees. New Forests, 6,257-278.
[13]	Epperson BK (1993). Recent advances in correlation studies of spatial patterns of genetic variation. Evolutionary Biology, 27,95-155.
[14]	Escudero A, Iriondo JM, Torres ME (2003). Spatial analysis of genetic diversity as a tool for plant conservation. Biological Conservation, 113,351-365.
[15]	Geburek T (1993). Are gene randomly distributed over space in mature populations of sugar maple ( Acer saccharum Marsh.)? Annals of Botany, 71,217-222.
[16]	Huang HW (黄宏文) (1998). Review of current research of the world Castanea species and importance of germplasm conservation of China native Castanea species. Journal of Wuhan Botanical Research (武汉植物学研究), 16,171-176. (in Chinese)
[17]	Lang P (郎萍), Huang HW (黄宏文) (1999). Genetic diversity and geographic variation in natural populations of the endemic Castanea species in China. Acta Botanica Sinica (植物学报), 41,651-657. (in Chinese with English abstract)
[18]	Levin DA, Kerster HW (1974). Gene flow in seed plants. Evoutionary Biology, 7,139-220.
[19]	Li CC, Horvitz DG (1953). Some methods of estimating the inbreeding coefficient. American Journal of Human Genetics, 5,107-117.
[20]	Li ZZ (李作洲), Lang P (郎萍), Huang HW (黄宏文) (2002). Spatial structure of allozyme frequencies in Castanea mollissima BL. Journal of Wuhan Botanical Research (武汉植物学研究), 20,165-170.. (in Chinese with English abstract)
[21]	Liu EH, Godt MJW (1983). The differentiation of populations over short distance. In: Schonewald-Cox CM, Chambers MS, MacBryde B, Thomas W eds. Genetics and Conservation. The Benjamin/Cummings Publishing Company Inc, California,78-95.
[22]	Marinoni D, Akkak A, Bounous G, Edwards KJ, Botta R (2003). Development and characterization of microsatellite markers in Castanea sativa(Mill.). Molecular Breeding, 11,127-136.
[23]	Marquardt PE, Epperson BK (2004). Spatial and population genetic structure of microsatellites in white pine. Molecular Ecology, 13,3305-3315.
[24]	McCauley DE, Stevens JE, Peroni PA, Raveill JA (1996). The spatial distribution of chloroplast DNA and allozyme polymorphisms within a population of Silene alba (Caryophyllaceae). American Journal of Botany, 83,727-731.
[25]	Peakall R, Smouse PE (2001). GenAlEx V5: Genetic Analysis in Excel. Population Genetic Software for Teaching and Research. Australian National University, Canberra, Australia.
[26]	Pigliucci M, Villani F, Benedettelli S (1990). Geographic and climatic factors associated with the spatial structure of gene frequencies in Castanea sativa Mill.from Turkey. Journal of Genetics, 69,141-149.
[27]	Robert R, Sokal FMLS (1978). Spatial autocorrelation in biology. Biological Journal of the Linnean Society, 10,199-228.
[28]	Rogers SO, Bendich AJ (1985). Extraction of DNA from milligram amount of fresh, herbarium and mummified plant tissues. Plant Molecular Biology, 5,69-76. URL PMID
[29]	Rutter PA, Miller G, Payne JA (1990). Chestnuts. In: Moore JN, Ballington JR Jr eds. Genetic Resources of Temperate Fruit and Nut Crops. The International Society for Horticultural Science, Wageningen, The Netherlands,761-788.
[30]	Sanguinetti CJ, Dias NE, Simpson AJG (1994). Rapid silver staining and recovery of PCR products separated on polyacrylamide gels. Biotechniques, 17,915-919.
[31]	Schaal BA, Levin DA (1976). The demographic genetics of Liatris cylindracea Michx. The American Naturalist, 110,191-206.
[32]	Slatkin M, Barton H (1989). A comparison of three indirect methods for estimating average levels of gene flow. Evolution, 43,1349-1368. DOI URL PMID
[33]	Sokal RR, Oden NL (1978a). Spatial autocorrelation in bilology. 1. Methodology. Biological Journal of the Linnean Society, 10,199-228.
[34]	Sokal RR, Oden NL (1978b). Spatial autocorrelation in bilology. 2. Some biological implications and four apolications of evolutionary and ecological interest. Biological Journal of the Linnean Society, 10,229-249. DOI URL
[35]	Sokal RR, Jacquez GM, Wooten M (1989). Spatial autocorrelation analysis of migration and selection. Genetics, 121,845-855. URL PMID
[36]	Streiff R, Labbe T, Bacilieri R, Steinkellner H, Glossl H, Kremer A (1998). Within-population genetic structure in Quercus robur L. & Quercus petraea (Matt.) Leibl. assessed with isozymes and microsatellites. Molecular Ecology, 7,317-328. DOI URL
[37]	Vekemans X, Hardy OJ (2004). New insights from fine-scale spatial genetic structure analyses in plant populations. Molecular Ecology, 13,921-935. URL PMID
[38]	Vornam B, Decarli N, Gailing O (2004). Spatial distribution of genetic variation in a natural beech stand ( Fagus sylvatica L.) based on microsatellite markers. Conservation Genetics, 5,561-570. DOI URL
[39]	Waternberg D (1989). SAAP 4.3.—A Spatial Autocorrelation Analysis Program. Exeter Software, Setauket, NY, USA.
[40]	Weir BS, Cockerham CC (1984). Estimating F-statistics for the analysis of population structure. Evolution, 38,1358-1370. DOI URL PMID
[41]	Williams CG, Hamrick JL (1996). Elite populations for conifer breeding and gene conservation. Canadian Journal of Forest Research, 26,453-461. DOI URL
[42]	Workman PL, Niswander JD (1970). Population studies on southwestern Indian tribes.II.Local genetic differentiation in the Papago. American Journal of Human Genetics, 22,24-49. URL PMID

引物编号 Locus code	片段大小范围 Size range (bp)	等位基因数 No.of alleles	预期杂合度 H_e	观察杂合度 H_o	近交系数 F_IS	居群分化系数 F_ST	哈迪-温格平衡检验 Hardy-weinberg equilibrium test
CsCAT1	160~201	19	0.879	0.538	0.359	0.094	NS
CsCAT2	180~226	13	0.864	0.774	0.034	0.078	NS
CsCAT3	197~215	6	0.637	0.594	-0.023	0.089	NS
CsCAT4	217~261	13	0.818	0.817	-0.047	0.054	NS
CsCAT5	217~242	7	0.753	0.607	0.181	0.031	NS
CsCAT7	190~217	6	0.761	0.538	0.256	0.038	NS
CsCAT8	190~210	6	0.729	0.658	0.117	0.018	NS
CsCAT14	150~160	6	0.754	0.770	-0.151	0.049	NS
CsCAT15	127~142	12	0.790	0.355	0.436	0.119	NS
平均值Mean		9.78	0.778	0.625	0.129	0.063	NS

引物编号 Locus code	片段大小范围 Size range (bp)	等位基因数 No.of alleles	预期杂合度 H_e	观察杂合度 H_o	近交系数 F_IS	居群分化系数 F_ST	哈迪-温格平衡检验 Hardy-weinberg equilibrium test
CsCAT1	160~201	19	0.879	0.538	0.359	0.094	NS
CsCAT2	180~226	13	0.864	0.774	0.034	0.078	NS
CsCAT3	197~215	6	0.637	0.594	-0.023	0.089	NS
CsCAT4	217~261	13	0.818	0.817	-0.047	0.054	NS
CsCAT5	217~242	7	0.753	0.607	0.181	0.031	NS
CsCAT7	190~217	6	0.761	0.538	0.256	0.038	NS
CsCAT8	190~210	6	0.729	0.658	0.117	0.018	NS
CsCAT14	150~160	6	0.754	0.770	-0.151	0.049	NS
CsCAT15	127~142	12	0.790	0.355	0.436	0.119	NS
平均值Mean		9.78	0.778	0.625	0.129	0.063	NS

居群 Population	平均等位基因数A	平均有效等位基因数A_e	预期杂合度 H_e	观察杂合度 H_o	特有等位基因数A_pe	近交系数 F_IS	取样个体数 Sampling size
随机大居群 Panmictric population	9.778(1.561 5)	4.954(0.595 8)	0.776(0.024 4)	0.625		0.190	197
广水 Guangshui	6.111(0.696 1)	3.578(0.332 9)	0.705(0.021 4)	0.634	3	0.109	43
麻城 Macheng	8.444(1.168 0)	4.537(0.556 7)	0.750(0.033 9)	0.620	20	0.169	118
岳西 Yuexi	6.000(0.500 0)	3.638(0.330 7)	0.704(0.030 2)	0.622	6	0.104	36

居群 Population	平均等位基因数A	平均有效等位基因数A_e	预期杂合度 H_e	观察杂合度 H_o	特有等位基因数A_pe	近交系数 F_IS	取样个体数 Sampling size
随机大居群 Panmictric population	9.778(1.561 5)	4.954(0.595 8)	0.776(0.024 4)	0.625		0.190	197
广水 Guangshui	6.111(0.696 1)	3.578(0.332 9)	0.705(0.021 4)	0.634	3	0.109	43
麻城 Macheng	8.444(1.168 0)	4.537(0.556 7)	0.750(0.033 9)	0.620	20	0.169	118
岳西 Yuexi	6.000(0.500 0)	3.638(0.330 7)	0.704(0.030 2)	0.622	6	0.104	36

位点 Locus	10个距离等级(km)的Moran's I值 Moran's I for 10 distance classes (km)										累积机率 Cum. Prob
位点 Locus	1	2	3	4	5	6	7	8	9	10
CsCAT1-180	0.32^**	0.26^**	0.31^**	0.16^**	-0.03	-0.13^**	-0.15^**	-0.55^**	-0.57^**	0.15^**	0.000^*
CsCAT1-188	0.01	0.02	-0.02	-0.07	0.02	-0.02	-0.00	-0.00	-0.01	0.01	0.045^*
CsCAT1-190	0.07^**	-0.09^**	0.04^*	-0.06^**	-0.04^*	0.01	0.00	-0.02	0.03	-0.02	0.001^*
CsCAT2-209	0.15^**	0.05^**	0.06^**	-0.01	-0.23^**	-0.13^**	0.02	-0.01	0.17^**	-0.11^**	0.000^*
CsCAT3-201	0.09^**	-0.02	0.03	-0.04	-0.17^**	0.04^*	-0.00	0.01	0.05^*	-0.04	0.000^*
CsCAT3-207	0.04^**	0.01	-0.01	0.01	0.07^**	-0.00	0.05^**	-0.04^*	-0.08^**	-0.13^**	0.000^*
CsCAT4-230	0.06^**	0.04^*	-0.02	0.02	-0.00	-0.03	-0.12^**	-0.02	-0.03	0.04^**	0.000^*
CsCAT4-238	0.10^**	0.05^**	-0.07^**	-0.03	0.02	0.04^**	0.04^*	-0.11^**	-0.03	-0.08^**	0.000^*
CsCAT4-242	0.07^**	0.08^**	0.01	0.15^**	0.00	-0.11^**	-0.20^**	-0.07^**	-0.11^**	0.11^**	0.000^*
CsCAT4-253	0.05^**	-0.05^*	-0.01	-0.05^*	-0.01	0.01	0.01	-0.02	-0.00	-0.00	0.056
CsCAT5-227	0.04^**	-0.03	0.02	-0.01	0.01	-0.03	-0.06^**	-0.01	-0.04	0.04^*	0.037^*
CsCAT5-232	0.12^**	0.07^**	0.07^**	0.04^*	0.07^**	-0.01	-0.05^*	-0.20^**	-0.13^**	-0.11^**	0.000^*
CsCAT7-205	0.16^**	0.08^**	-0.08^**	0.00	-0.06^**	0.01	-0.08^**	-0.18^**	0.00	0.09^**	0.000^*
CsCAT7-209	0.05^**	0.05^*	-0.04^*	-0.01	-0.08^**	0.02	0.00	-0.08^**	-0.01	0.05^**	0.000^*
CsCAT7-215	0.02	-0.03	-0.03	0.02	-0.02	0.02	-0.03	0.01	-0.01	0.00	1.000
CsCAT7-217	-0.00	-0.02	0.02	0.00	0.02	-0.00	0.03	-0.01	-0.03	-0.07^*	0.176
CsCAT8-190	0.02	-0.02	0.01	-0.01	-0.03	-0.05	0.02	-0.01	-0.02	0.02	0.780
CsCAT8-201	0.02	0.04^*	0.04^*	0.02	-0.05	-0.03	0.05^*	-0.05	-0.11^**	-0.00	0.001^*
CsCAT8-204	-0.04	0.03	-0.02	-0.01	0.01	-0.03	0.01	-0.01	0.00	-0.01	0.798
CsCAT8-207	0.04^*	0.01	-0.00	0.05^**	-0.04^*	-0.02	0.04^**	-0.07^**	-0.09^**	0.01	0.000^*
CsCAT14-150	-0.02	0.00	-0.01	-0.03	0.03^*	-0.03	0.02	-0.01	0.00	-0.00	0.400
CsCAT14-153	0.03^*	0.01	-0.01	0.05^**	0.03^*	-0.04	-0.05^*	-0.01	0.02	-0.10^**	0.000^*
CsCAT14-155	0.14^**	-0.01	0.05^**	0.19^**	0.07^**	-0.12^**	-0.11^**	-0.02	0.09^**	-0.41^**	0.000^*
CsCAT14-157	0.10^**	0.08^**	0.14^**	0.11^**	-0.13^**	-0.26^**	-0.08^**	-0.02	-0.02	0.02	0.000^*
CsCAT15-133	0.19^**	0.16^**	-0.01	-0.04	-0.16^**	-0.15^**	-0.04	0.09^**	0.01	-0.09^**	0.000^*
CsCAT15-134	0.16^**	0.11^**	0.09^**	0.12^**	0.04^*	-0.00	-0.12^**	-0.28^**	-0.06^**	-0.08^**	0.000^*
CsCAT15-135	0.02	0.02	0.05^*	-0.00	-0.06^**	-0.06^**	-0.02	-0.02	-0.05^*	0.06^**	0.022^*
平均Average	0.08	0.03	0.02	0.02	-0.03	-0.04	-0.03	-0.06	-0.04	-0.03	-
最大值Maximum	0.32	0.26	0.31	0.19	0.07	0.04	0.05	0.09	0.17	0.15	-
最小值Minimum	-0.04	-0.09	-0.08	-0.07	-0.23	-0.26	-0.20	-0.55	-0.57	-0.41	-