林业科学  2011, Vol. 47 Issue (3): 59-67   PDF    
0

文章信息

周炎花, 乔小燕, 马春雷, 乔婷婷, 金基强, 姚明哲, 陈亮
Zhou Yanhua, Qiao Xiaoyan, Ma Chunlei, Qiao Tingting, Jin Jiqiang, Yao Mingzhe, Chen Liang
广西茶树地方品种遗传多样性和遗传结构的EST-SSR分析
Genetic Diversity and Structure of Tea Landraces from Guangxi ased on EST-SSR Analysis
林业科学, 2011, 47(3): 59-67.
Scientia Silvae Sinicae, 2011, 47(3): 59-67.

文章历史

收稿日期:2010-01-07
修回日期:2010-04-29

作者相关文章

周炎花
乔小燕
马春雷
乔婷婷
金基强
姚明哲
陈亮

广西茶树地方品种遗传多样性和遗传结构的EST-SSR分析
周炎花1,2, 乔小燕3, 马春雷1, 乔婷婷1, 金基强1, 姚明哲1, 陈亮1    
1. 中国农业科学院茶叶研究所 国家茶树改良中心 杭州 310008;
2. 漳州天福茶职业技术学院 漳州 363202;
3. 广东省农业科学院茶叶研究所 广州 510640
摘要: 利用设计合成的60对扩增稳定的EST-SSR核心引物分析51份广西茶树地方品种的遗传多样性和遗传结构,共检测到232个等位基因,每对引物检测到2~7个,平均3.88个。群体的观测杂合度(Ho)变化范围为0.02 (TM149)~0.92(TM186),平均为0.37;期望杂合度(He)的变化范围为0.13(TM144,TM147)~0.73(TM189),平均为0.47;位点的多态信息含量(PIC)介于0.12~0.68之间,平均为0.41。地区间遗传相似系数在0.83~0.94之间,说明不同地区地方品种间的亲缘关系较近。对于茶与白毛茶群体而言,平均遗传分化系数为0.02 (P < 0.001),即有2%的遗传变异存在于群体间,而有98%的遗传变异存在于群体内。另外,种群内的近交系数Fis平均为0.21,说明群体内近交现象较严重。基于Nei’s遗传距离,利用PowerMarker 3.25软件进行Neighbor-Joining聚类分析,将51份地方品种聚为3类:桂林、来宾和贺州的大部分品种,钦州和梧州的全部品种聚在第Ⅰ类; 第Ⅱ类包含多数地区的品种; 第Ⅲ类主要由崇左、防城港和南宁的品种组成。该结果与基于模型的遗传结构分析结果基本一致。
关键词:茶树    遗传多样性    遗传结构    EST-SSR    
Genetic Diversity and Structure of Tea Landraces from Guangxi ased on EST-SSR Analysis
Zhou Yanhua1,2, Qiao Xiaoyan3, Ma Chunlei1, Qiao Tingting1, Jin Jiqiang1, Yao Mingzhe1, Chen Liang1    
1. National Center for Tea Improvement Tea Research Institute, the Chinese Academy of Agricultural Sciences Hangzhou 310008;
2. Zhangzhou Tenfu Tea College Zhangzhou 363202;
3. Tea Research Institute of the Guangdong Academy of Agricultural Sciences Guangzhou 510640
Abstract: Sixty core EST-SSR primers, among which 59 were newly developed from the ESTs of tea plant, were selected and used to analyze genetic diversity and genetic structure of 51 landraces of tea (Camellia sinensis and C. sinensis var. pubilimba) distributed in Guangxi, China. A dendrogram based on genetic distance was constructed from software PowerMarker 3.25 using Neighbor-Joining clustering method. Two hundred and thirty-two alleles were amplified, and the number of alleles per primer ranged from 2 to 7, with an average of 3.88. The observed heterozygosity (Ho) ranged from 0.02 to 0.92 but expected heterozygosity (He) varied from 0.13 to 0.73, with average of 0.37 and 0.47, respectively. The average polymorphism information content (PIC) was 0.41. The genetic similarity coefficient among eight different regions was from 0.83 to 0.94, indicating there were close relationships among tea landraces in the different regions. Most of the diversity (98%) was found within populations, only 2% was between populations of C. sinensis and C. sinensis var. pubilimba. It was suggested that serious inbreeding occurred within the groups as shown by the coefficient of inbreeding (Fis=0.21). The 51 landraces from the different regions were grouped into three. Most tea landraces from Guilin, Laibin and Hezhou, and all from Qinzhou and Wuzhou were clustered into groupⅠof the dendrogram, while those from Chongzuo, Fangchenggang and Nanning were assigned into group Ⅲ. Other landraces were assigned into group Ⅱ. The results from Neighbor-Joining clustering method agreed with that from the model-based analysis of genetic structure.
Key words: Camellia sinensis    Camellia sinensis var. pubilimba    genetic diversity    genetic structure    EST-SSR    

茶树起源于中国西南,滇、桂、黔毗邻地区为起源的中心(闵天禄,1992),广西正处在这个起源中心的边缘,其气候夏长冬暖、炎热多雨,极适宜于各种类型的山茶科(Theaceae)植物的生长,茶树种质资源非常丰富。在全国茶叶生产区划中,广西被列为最适生长区(陈宗懋,2007; 陈亮等,1996)。在国家种质杭州茶树圃保存的2 665份种质资源中,来自广西的种质资源为149份(陈亮等,2004),仅次于云南、福建、浙江,位列第4。广义的茶树包含山茶属茶组[Camellia L.Sect.Thea (L.) Dyer]所有种和变种,目前茶组植物分类系统主要有席勒系统(Sealy,1958)、张宏达系统(张宏达,1981a; 1981b; 1984; 1990)、闵天禄系统(闵天禄,1992)等。广西茶树种质资源除桂西有一些大厂茶(C.tachangensis Chang)和厚轴茶(C.crassicolumna Chang)外,基本上属于茶(C.sinensis (L.) O.Kuntze)和白毛茶(C.sinensis var.pubilimba Chang) (陈亮等,2006)。广西拥有很多历史悠久的地方品种及名茶,如‘苍梧六堡茶’、‘贺州开山白毛茶’、‘武宣金龙茶’、‘昭平象棋茶’、‘资源云雾茶’、‘钟山雷电茶’、‘三江牙己茶’等(陈宗懋,2007)。据统计,到2010年底全国通过国家审(认、鉴)定的茶树品种有123个,而广西仅有‘桂红3号’、‘桂红4号’、‘桂绿1号’、‘尧山秀绿’和‘桂香18’5个国家级无性系良种; 通过省级审(鉴)定的品种有150多个(Chen et al., 2007),广西仅有‘桂热1号’、‘桂热2号’2个。茶树地方品种资源在新品种选育中有很大的潜力,而对其进行准确和有效的遗传多样性和遗传结构分析,是茶树遗传改良利用的前提。

SSR (Simple Sequence Repeat)标记以其多态性高、重复性好、可靠性高、操作简单等优点,已成为大豆(Glycine max)、番茄(Solanum lycopersicum)、小麦(Triticum aestivum)、玉米(Zea mays)、水稻(Oryza sativa)等作物遗传育种研究的最重要方法之一(吴晓雷等,2001; He et al., 2003; Gupta et al., 2003; 朴红梅等,2007; Gao et al., 2008),在茶树遗传研究中得到越来越广泛的应用(Freeman et al., 2004; 金基强等,2007; Kaundun et al., 2002; Ohsako et al., 2008; Zhao et al., 2008; Sharma et al., 2009)。刘振等(2008; 2009)姚明哲等(2009)通过EST-SSR (Expressed Sequence Tag based SSR)技术对我国云南、四川、贵州、福建等省份的茶树种质资源进行了比较详尽的遗传多样性研究,而对毗邻的广西茶树地方品种资源未见系统的研究报道。本研究利用EST-SSR标记分析了广西茶树地方品种的遗传多样性和遗传结构,比较并分析了茶与白毛茶变种的遗传差异,为更好地利用和改良广西茶树地方品种提供参考。

1 材料与方法 1.1 材料

共选用了原产于广西有代表性的51份茶树地方品种(表 1),分类上均为茶(C.sinensis (L.) O. Kuntze),其中白毛茶(C.sinensis var. pubilimba Chang)31份。用于DNA提取的一芽二叶鲜叶采自中国农业科学院茶叶研究所的“国家种质杭州茶树圃”,样品采集后经液氮迅速冷冻处理,保存于-70 ℃冰箱中备用。

表 1 广西51份茶树地方品种 Tab.1 Fifty-one accessions of tea landraces from Guangxi
1.2 试验方法 1.2.1 DNA提取

采用改进的SDS法(陈亮等,1997)分别提取51个地方品种的基因组DNA,0.8%琼脂糖凝胶电泳检测DNA质量,ND-1000 UVVis分光光度计(Nanodrop Technologies,Inc.,Wilmington,DE,USA)检测DNA浓度后,稀释部分DNA原液到10 ng·μL-1的工作浓度,用于PCR扩增,剩余原液-20 ℃保存备用。

1.2.2 EST-SSR引物的选择

所用的60对ESTSSR核心引物重复性好、扩增条带清晰稳定、具多态性,其中59对由本实验室设计、筛选(乔婷婷等,2010),1对为Freeman等(2004)发表,作为参照引物,标记为Camsin M5 (表 2)。所有引物均由上海生工生物工程技术服务有限公司合成。

表 2 60个微卫星标记的特征及在广西地方品种中的扩增结果和F统计量 Tab.2 The characteristics, amplification information of tea landraces in Guangxi and F-statistics of 60 EST-SSRs
1.2.3 PCR扩增及其产物检测

SSR-PCR反应体系与反应程序参照刘振等(2008)的方法,在BioRad PTC-220型PCR仪(MJ Research,Inc.,Waltham,MA,USA)上进行PCR扩增。扩增产物中加入2 μL的上样缓冲液(loading buffer),用10%的聚丙烯酰胺凝胶,在CBS垂直电泳系统(C.B. S.Scientific Company,Inc.,San Diego,CA,USA) 150 V电压下电泳90 min。电泳后参照Charters和Wilkinson(2000)的方法银染显色。Bio-Rad Gel Doc 2000 (Bio-Rad Laboratories,Inc.,Philadelphia,PA,USA)凝胶成像仪拍照记录。

1.3 数据统计

EST-SSR为共显性标记,同一引物扩增产物中电泳迁移率一致的条带被认为具有同源性,其显带的记录方法为:电泳图上目的片段范围内,清晰、重复的条带,根据分子质量大小,从大到小依次记作A,B,C,D……,1条为纯合,2条为杂合,分别记录为aa,ab,bb,ac,ad,cd等,按照等位基因的位置不同而不同,从而获得原始的数据。

1.3.1 遗传多样性及遗传分化分析

利用Popgene version 1.32软件(http://cc.oulu.fi/~jaspi/popgen/popdown.htm)计算在物种水平和群体水平的等位基因数(N)、Shannon信息指数(I)、观测杂合度(Ho)、期望杂合度(He),以及遗传分化的固定指数(FitFstFis),即总近交系数Fit、种群分化系数Fst和种群内的近交系数Fis。利用Powermarker 3.25软件(http://www.powermarker.net)计算Nei基因多样性(H)和多态信息含量(PIC) (Liu et al., 2005)。

1.3.2 聚类分析

利用PowerMarker 3.25软件,基于Nei’s遗传距离(Nei,1973)进行Neighbor-Joining聚类分析,构建无根树状聚类图。

1.3.3 遗传结构分析

采用Structure 2.0软件(http://pritch.bsd.uchicago.edu/software)分析群体的遗传结构(Pritchard et al., 2000),假定试材的群体数目(K)为2~9,将迭代周期(length of burnin period)和MCMC (Markov Chain Monte Carlo)各设为10 000次,重复3次计算,重复得到一致的结果,然后根据似然值最大的原则选取1个合适的K值。

2 结果与分析 2.1 遗传多样性分析

应用60个EST-SSR标记在广西51个茶树地方品种资源中,共检测到233个等位基因,平均每个位点约有3.88个等位基因。如表 2所示,TM154位点的等位基因最多,有7个; 位点TM142,TM144,TM149,TM180的等位基因数最少,为2个。对于整个群体,群体的观测杂合度(Ho)变化范围为0.02 (TM149)~0.92 (TM186),平均为0.37;期望杂合度(He)的变化范围为0.13(TM144,TM147)~0.73 (TM189),平均为0.47。Shannon信息指数(I)为0.18(TM144)~1.42 (TM169),平均为0.82。所使用的60个微卫星标记的多态性信息含量(PIC)为0.12~0.68,平均为0.41。根据Botstein等(1980)的理论,有8个标记为低度多态位点(0 < PIC < 0.25),33个标记为中度多态位点(0.25 < PIC < 0.5),其他19个为高度多态位点(PIC > 0.5),TM169位点PIC最高为0.68,TM144,TM147位点PIC最低为0.12,其他位点的PIC值见表 2

2.2 聚类分析

广西51份地方品种的亲缘关系聚类如图 1所示,全部品种聚为3类:第Ⅰ类共有24份,包含了桂林(7份)、来宾(3份)和贺州(2份)的大部分品种,钦州(4份)和梧州(3份)的全部品种,还有少部分来自百色(1份)、崇左(1份)、防城港(1份)、贵港(1份)和柳州(1份); 第Ⅱ类共有13份,来自崇左(3份)、桂林(3份)、贺州(2份)、来宾(2份)、南宁(2份)和柳州(1份); 第Ⅲ类包含14份品种,主要来自崇左(5份)、防城港(7份),还有部分来自南宁(2份)。

图 1 51份茶树资源EST-SSR数据的NJ聚类分析 Figure 1 Dendrogram of 51 accessions tea germplasms based on EST-SSR data using Neighbor-Joining clustering method 编号见表 1 See Tab. 1 for details.

另外,茶与白毛茶在3大类群中的分布见图 2,茶与白毛茶在各类群中穿插分布,茶主要分布在第Ⅲ类; 白毛茶主要分布在第Ⅰ,Ⅱ类,且第Ⅱ类以白毛茶为主。

图 2 茶与白毛茶在NJ聚类图中的分布 Figure 2 Distribution of tea cultivars of C. sinensis and C. sinensis var. pubilimba in population Ⅰ, Ⅱ and Ⅲ
2.3 遗传结构分析

利用Structure 2.0软件对广西51份茶树地方品种进行遗传结构分析,当K=3时,能够取得最大的lnP(D)值和比较稳定的α值。即广西茶树地方品种可分为3个类群(图 3) :类群A主要由桂林、钦州、贺州和梧州的组成; 类群B包含了多数地区的品种; 类群C包含了防城港、崇左和南宁的大部分品种。各地区在遗传结构图中的具体分布如表 3所示。茶与白毛茶在各类群中的分布如图 4所示,茶主要分布于类群C,白毛茶主要布于类群A和B。基于模型和遗传距离的分析结果基本一致,广西南部(崇左、防城港和南宁)的茶树地方品种和北部(桂林、钦州、梧州和来宾)的品种分别聚在一起,茶与白毛茶在各类群中穿插分布。

图 3 基于模型的遗传结构 Figure 3 The model-based analysis of genetic structure 编号见表 1 See Tab. 1 for details.
表 3 各地区地方品种在遗传结构图中的分布 Tab.3 Distribution of accessions from different regions in population structure
图 4 茶与白毛茶在数学模型聚类图中的分布 Figure 4 Distribution of tea cultivars of C. sinensis and C. sinensis var. pubilimba in genetic structure
2.4 地区间遗传多样性和遗传距离分析

为了进一步明确地区间的遗传差异,对广西8个主要地区共47份茶树地方品种的平均观测杂合度和平均期望杂合度进行分析(表 4)。就各地区而言,其平均期望杂合度在0.42(梧州)~0.48(南宁)之间,平均观测杂合度0.34(钦州)~0.40(来宾)之间,Shannon信息指数在0.55(梧州)~0.72(崇左,防城港)之间,说明地区间的遗传差异较小。8个地区间的亲缘关系见图 5,来源于梧州和钦州的地方品种聚为一类; 南宁、崇左及防城港的先聚为一小类,再与桂林和来宾聚为一类; 贺州单独为一类。

表 4 不同地区间遗传多样性分析 Tab.4 Genetic diversity analysis in different regions
图 5 基于Nei遗传距离的各地区地方品种间亲缘关系的NJ聚类 Figure 5 Dendrogram of phylogenetic relationships among cities of Guangxi based on Nei's genetic distance, using Neighbor-Joining method
2.5 茶与白毛茶遗传多样性和遗传分化分析

茶与白毛茶的遗传多样性比较分析结果表明:茶的I,PIC和Nei基因多样性(H)分别为0.80,0.41,0.46,均略高于白毛茶(图 6)。另外,将茶品种与白毛茶品种分为2个群体,根据Popgene的计算,60对EST-SSR标记在2个群体的F-统计量结果(表 2)可知: 60个标记中TM127,129,137,139等12个标记显示出杂合子剩余,种群内的近交系数Fis平均为0.21,说明群内近交现象较严重。另外,种群分化系数为-0.76 (TM186)~0.89 (TM149),Fst平均为0.02,即有2%的变异来源于群体间,有98%的变异来源于群体内。

图 6 茶与白毛茶资源的遗传多样性比较 Figure 6 Genetic diversity comparison between C. sinensis and C. sinensis var. pubilimba
3 结论与讨论 3.1 广西茶树地方品种的遗传多样性

等位基因数是反映群体遗传多样性的指标之一,在保护遗传学研究中,尤其更强调等位基因数对研究种群遗传变异的影响。在茶树遗传研究中,其有效等位基因数有待进一步确定。Ohsako等(2008)采用6个微卫星标记对日本京都8个群体共146个地方品种的遗传多样性进行研究,检测到的等位基因数为3.83~4.67。Sharma等(2009)采用61对微卫星标记检测来自茶树及近缘种的34份种质资源,检测到的等位基因数为2~16,平均为5.31。Zhao等(2008)采用茶树及其变种共检测到等位基因数平均为4.46,而刘振等(2008; 2009)采用微卫星标记检测我国西南地区及福建茶树资源的平均等位基因数分别为4.74和3.30。本研究所采用的60个微卫星标记,平均等位基因为3.88,与上述的研究结果相近,说明本研究采用的EST-SSR标记对揭示广西茶树遗传多样性水平是可靠的。

另外,遗传杂合度是衡量群体内物种遗传变异的重要指标,其值高则遗传变异大,反之则变异小。而遗传变异的大小是衡量选择潜力的依据,变异大则选择潜力大,反之则选择潜力小。采用60个标记对广西51个茶树地方品种进行检测,发现其平均观测杂合度(Ho)为0.37,平均期望杂合度(He)为0.47,Shannon信息指数(I)平均为0.82,多态信息含量(PIC)平均为0.41。根据左志明等(1994)对桂西山区茶树种质资源考察发现,桂西山区各种类型的茶树种质资源非常丰富。本试验由于客观条件等原因,采用样品中包含桂西地方品种较少(尤其是野生资源),对试验结果可能有一定影响。另外,本文主要针对广西茶树地方品种进行研究,研究发现各地区地方品种间的遗传多样性比较接近,这在一定程度上降低了整体的多态性水平。

3.2 广西不同地区茶树地方品种间的亲缘关系

基于Nei's (Nei,1973)遗传距离,利用PowerMarker 3.25软件进行NJ聚类分析,构建51份地方品种的无根树状聚类图。如图 1所示,全部品种聚为3类:桂林和来宾的大部分品种、钦州和梧州的全部品种聚在第Ⅰ类; 第Ⅱ类包含了多数地区的品种; 第Ⅲ类主要由崇左、防城港的品种组成。另外,基于模型的遗传结构分析表明,当K=3时,能够取得最大的lnP(D)值和比较稳定的α值,即广西茶树地方品种可分为3个类群(图 3) :类群A主要由桂林、钦州、贺州和梧州的品种组成; 类群B包含了多数地区的品种; 类群C包含了防城港、崇左和南宁的大部分品种。基于模型和遗传距离的分析结果基本一致,广西南部(崇左、防城港和南宁)和北部(桂林、钦州、梧州和来宾)的茶树地方品种分别聚在一类,即来自防城港、崇左和南宁的地方品种亲缘关系较近,来自桂林、钦州、贺州、来宾和梧州的品种亲缘关系较近。

3.3 茶与白毛茶的分化情况

一个自然群体由于各种原因往往会发生变异,形成各种亚群,亚群间的遗传差异就是亚群间的遗传分化。本研究采用F-统计量(FitFstFis),即总近交系数Fit、种群分化系数Fst、种群内的近交系数Fis来衡量遗传分化。Fis值变化在-1~1之间,值越小,杂合子越多; 值越大,纯合子越高; 值为0则说明群体是随机交配的。本研究采用的60个标记中只有TM127,129,137,139等12个标记表现出杂合子剩余,其余的位点表现杂合子缺失,说明茶与白毛茶群体间的近交现象比较严重,导致杂合子缺失。另外,种群分化系数Fst为0.02,即有2%的变异来源于群体间,有98%的变异来源于群体内。

根据张宏达(1981a; 1981b; 1984; 1990)、闵天禄(1992)陈亮等(2000)对茶树系统分类的研究,白毛茶都是作为茶的一个变种,在表型方面二者表现出较高的相似性,其区别主要在于白毛茶花萼和嫩枝外部密被茸毛,叶片长多内折,叶面光滑,叶质较硬多革质。本研究结果表明,二者之间的遗传分化不显著,可能因为茶与白毛茶本身的相似度较高,本研究采用的样品量或标记量不足于揭示二者之间的差异,在今后的研究中应适当增加标记数量。另外,基于遗传距离和基于模型的聚类分析结果一致,均表明茶与白毛茶在不同类群中穿插分布现象严重,可能是由于地区间基因的频繁交流,导致遗传背景复杂。因此,白毛茶的分类归属问题值得进一步深入研究。

EST-SSR方法揭示了广西茶树地方品种遗传多样性较丰富,但各地区间(不包含桂西)的遗传相似性较高,区域差异对品种变异影响不显著,茶与白毛茶群体间的遗传分化较低。因此在今后的茶树选育过程中,在充分利用地方品种的同时,还应着重引入外源的遗传变异。

参考文献(References)
陈亮, 虞富莲, 左志明. 1996. 广西茶树种质资源考察研究初报[J]. 广西农业科学, (2): 77-80.
陈亮, 陈大明, 高其康, 等. 1997. 茶树基因组DNA提纯与鉴定[J]. 茶叶科学, 17(2): 172-176.
陈亮, 虞富莲, 童启庆. 2000. 关于茶组植物分类与演化的讨论[J]. 茶叶科学, 20(2): 89-94.
陈亮, 杨亚军, 虞富莲. 2004. 中国茶树种质资源研究的主要进展和展望[J]. 植物遗传资源学报, 5(4): 389-392.
陈亮, 虞富莲, 杨亚军, 等. 2006. 茶树种质资源与遗传改良[M]. 北京: 中国农业科学技术出版社.
陈宗懋. 2007. 中国茶经[M]. 上海: 上海文化出版社.
金基强, 崔海瑞, 龚晓春, 等. 2007. 用EST-SSR标记对茶树种质资源的研究[J]. 遗传, 29(1): 103-108.
刘振, 王新超, 赵丽萍, 等. 2008. 基于EST-SSR的西南茶区茶树资源遗传多样性和亲缘关系分析[J]. 分子植物育种, 6(1): 100-110.
刘振, 姚明哲, 王新超, 等. 2009. 基于EST-SSR的福建地区茶树资源遗传多样性和亲缘关系分析[J]. 中国农业科学, 42(5): 1720-1727.
闵天禄. 1992. 山茶属茶组植物的订正[J]. 云南植物研究, 14(2): 115-132.
朴红梅, 王玉民, 李万良, 等. 2007. SSR标记在玉米育种中的研究与应用[J]. 玉米科学, 15(S1): 157-159.
乔婷婷, 马春雷, 周炎花, 等. 2010. 浙江省茶树地方品种与选育品种遗传多样性和群体结构的EST-SSR分析[J]. 作物学报, 36(5): 744-753.
吴晓雷, 贺超英, 陈受宜, 等. 2001. 用分子标记研究大豆属种间亲缘进化关系[J]. 遗传学报, 28(4): 359-366.
姚明哲, 刘振, 陈亮, 等. 2009. 利用EST-SSR分析江北茶区茶树资源的遗传多样性和遗传结构[J]. 茶叶科学, 29(3): 243-250.
张宏达. 1981a. 茶树的系统分类[J]. 中山大学学报:自然科学版, (1): 87-99.
张宏达. 1981b. 山茶属植物的系统研究[J]. 中山大学学报:自然科学版(论丛), (1): 1-124.
张宏达. 1984. 茶叶植物资源的修订[J]. 中山大学学报:自然科学版, (1): 1-12.
张宏达. 1990. 中国山茶科植物新种[J]. 中山大学学报:自然科学版, 29(2): 85-93.
左志明, 陈亮, 虞富莲. 1994. 桂西山区茶树种质资源考察初报[J]. 广西农业科学, (2): 63-65.
Botstein D, White R L, Skolnick M, et al. 1980. Construction of a genetic linkage maps in man using restriction fragment length polymorphisms[J]. The American Journal of Human Genetics, 32(3): 314-331.
Charters Y M, Wilkinson M J. 2000. The use of self-pollinated progenies as 'in-groups' for the genetic characterization of cocoa germplasm[J]. Theoretical and Applied Genetics, 100(1): 160-166. DOI:10.1007/PL00002903
Chen Liang, Zhou Zhixiu, Yang Yajun. 2007. Genetic improvement and breeding of tea plant (Camellia sinensis) in China: from individual selection to hybridization and molecular breeding[J]. Euphytica, 154(1/2): 239-248.
Freeman S, West J, James C, et al. 2004. Isolation and characterization of highly polymorphic microsatellites in tea (Camellia sinensis)[J]. Molecular Ecology Notes, 4(3): 324-326. DOI:10.1111/men.2004.4.issue-3
Gao L, Innan H. 2008. Nonindependent domestication of the two rice subspecies, Oryza sativa spp. indica and ssp. japonica, demonstrated by multilocus microsatellites[J]. Genetics, 179(2): 965-976. DOI:10.1534/genetics.106.068072
Gupta P K, Rustgi S, Sharma S, et al. 2003. Transferable EST-SSR markers for the study of polymorphism and genetic diversity in bread wheat[J]. Molecular Genetics and Genomics, 270(4): 315-323. DOI:10.1007/s00438-003-0921-4
He C, Poysa V, Yu K. 2003. Development and characterization of simple sequence repeat (SSR) markers and their use in determining relationships among Lycopersicon esculentum cultivars[J]. Theoretical and Applied Genetics, 106(2): 363-373. DOI:10.1007/s00122-002-1076-0
Kaundun S S, Matsumoto S. 2002. Heterologous nuclear and chloroplast microsatellite amplification and variation in tea, Camellia sinensis[J]. Genome, 45(6): 1041-1048. DOI:10.1139/g02-070
Liu K, Muse S V. 2005. Integrated analysis environment for genetic marker data[J]. Bioinformatics, 21(9): 2128-2129. DOI:10.1093/bioinformatics/bti282
Nei M. 1973. Analysis of gene diversity in subdivided populations[J]. Proceedings of the National Academy of Sciences of the United States of America, 70(12): 3321-3323. DOI:10.1073/pnas.70.12.3321
Ohsako T, Ohgushi T, Motosugi H, et al. 2008. Microsatellite variability within and among local landrace populations of tea, Camellia sinensis (L.) O. Kuntze, in Kyoto, Japan[J]. Genetic Resources and Crop Evolution, 55(7): 1047-1053. DOI:10.1007/s10722-008-9311-4
Pritchard J K, Stephens M, Donnelly P. 2000. Inference of population structure using multilocus genotype data[J]. Genetics, 155(2): 945-959.
Sealy J R. 1958. A revision of the genus Camellia[M]. London: Royal Horticultural Society.
Sharma R K, Bhardwaj P, Negi R, et al. 2009. Identification, characterization and utilization of unigene derived microsatellite markers in tea (Camellia sinensis L.[J]. BMC Plant Biology, 9: 53. DOI:10.1186/1471-2229-9-53
Zhao Liping, Liu Zhen, Chen Liang, et al. 2008. Generation and characterization of 24 novel EST derived microsatellites from tea plant (Camellia sinensis) and cross-species amplification in its closely related species and varieties[J]. Conservation Genetics, 9(5): 1327-1331. DOI:10.1007/s10592-007-9476-y