茶树是山茶科山茶属茶组的多年生常绿木本植物[1]，且为常异花授粉植物，具有高度杂合性，其遗传基础较为复杂，种质资源的农艺性状与加工品质鉴定难度大、时间长，对茶树种质资源保护与新品种选育工作造成较大的困难。近些年，EST-SSR来自于功能基因，具有多态性高、共显性、技术简单、重复性好等特点[2-5]，在茶树遗传多样性、亲缘关系分析、群体结构分析、种质资源鉴定筛选与辅助育种中广泛应用[6-9]。对茶树种质资源遗传多样性与亲缘关系分析是野生茶树种质初步鉴定的重要手段。黄晓霞等[10]的研究结果说明千家寨不同海拔梯度上野生茶树居群遗传多样性随着海拔梯度的递增，居群的遗传多样性呈现出低-高-低的分布；毛娟等[11]发现样本量在40个时，能较好地反映云南白莺山地区茶树资源的遗传多样性；张明泽等[12]认为黔南野生茶树资源间的遗传差异较大，遗传基础较宽，具有丰富的遗传多样性。郭燕等[13]研究表明40个贵州古茶树资源有较高的遗传多样性，主要是由遗传因素和地理位置不同引起。周萌等[14]认为树龄3200年的香竹簧茶树王与野生茶树、栽培型茶树品种之问的遗传关系较远。姜晓辉等[15]研究表明广东历史名茶茶树群体罗坑与马图群体间在过去的某个时间可能有过基因交流。茶树种质资源的群体结构分析也有较多研究。王丽鸳等[16]认为龙井群体的遗传分化程度较低。姚明哲等[17]认为江北茶区主要省份间茶树种质的遗传分化程度低。乔婷婷等[18]认为茶树地方品种、选育品种(系)具有相对独立的群体结构，选育品种(系)根据亲缘关系的不同形成不同的类群。吴晓梅等[19]认为乌龙茶品种与绿茶品种间的遗传结构存在差异。本试验对福建省主要审(鉴)定品种、茶树自然杂交后代群体、地方搜集资源群体进行个体材料的亲缘关系、遗传多样性与群体结构分析，以期为茶树种质初步鉴定与分子标记辅助育种提供参考价值。

将校园打造成为具有典型北方自然特色及人文精神的绿色生态型校园,亦为我校师生营造优美舒适、刚柔并济的生态功能型校园。同时，立足气候、地址等原生特征，力求打造吉林工师的校园文化旅游名片和校园文化景观品牌。

1 材料与方法

1.1 试验材料

参试材料来源于福建省农业科学院茶叶研究所，福建省茶树种质资源及其自然杂交后代共94个(详见表1)，采集各个材料一芽二叶鲜叶后用液氮迅速冷冻处理，保存在-70℃冰箱中备用。

表1 94份供试材料的名称与来源

Table 1 Name and origin of 94 germplasms in collection

群体名称代号材料来源群体名称代号材料来源群体1A1-A6白鸡冠自然杂交后代群体9L1九龙大白茶A7白鸡冠L2大叶乌龙群体2B1-B10福云6号自然杂交后代L3绿芽佛手B11福云6号L4铁观音群体3C1-C 10金牡丹自然杂交后代L5梅占C 11金牡丹L6福鼎大白茶群体4D1-D 10茗科1号(金观音)自然杂交后代L8矮脚乌龙D 11茗科1号(金观音)L9肉桂群体5E1-E10紫玫瑰自然杂交后代L10白芽奇兰E11紫玫瑰L11黄棪群体6F1- 10悦茗香自然杂交后代F11悦茗香群体7G1- G11武夷山茶树种质资源群体8H1-H 11武夷山茶树种质资源

1.2 试验方法

1.2.1 基因组DNA的提取 采用改进的CTAB法[20]提取茶树基因组DNA。用1.0%的琼脂糖凝胶电泳进行茶树基因组分子量大小检测，用756-MC型紫外分光光度计测定茶树基因组DNA的纯度(OD260/280比值)。

1.2.2 引物合成 参照文献[21]引物序列，序列由上海Sangon公司合成，引物信息见表2。

从图1上看出，白鸡冠与其自然杂交后代4个聚类在一起，说明白鸡冠自然杂交后代遗传组成主要来自于母本白鸡冠，后代遗传相对稳定(i)；福云6号与其自然杂交后代分为2个类群，其中，第一类群B1、B2、B3与福云6号(B11)聚类在一起(ii)，B4、B5、B8、B10聚类在一起，推测B1、B2、B3与福云6号(B11)亲缘关系较近；金牡丹(C11)及其所有自然杂交后代都聚类在一起(iii)，说明其后代遗传组成主要来自于母本金牡丹；金观音(D11)与D2、D4、D5、D6、D8聚类在一起(iv)，推测其与母本金观音亲缘关系较近；9个材料紫玫瑰自然杂交后代聚类在一起(v)，却没有和亲本(E11)聚类在一起，推测父本为同一类型材料，且遗传组成主要来自于父本；7个悦茗香自然杂交后代与悦茗香(F11)聚类在一起(vi)，说明其后代遗传成分主要来自于母本悦茗香；武夷山种质资源G1、G3、G4、G6、G10、H1、H2、H3、H10、与肉桂聚类在一起(vii)，推测与肉桂亲缘关系较近；武夷山地方资源G5、H6、H8、H5、H11与大叶乌龙(L2)、梅占(L5)、矮脚乌龙(L8)、白芽奇兰(L10)聚类在一起(viii)，武夷山种质资源与这些品种遗传组成较为相似；武夷山种质资源H7与所有参试材料都没有聚类在一起，说明H7遗传组成相对特殊。

孟导权衡了一会，便答应了老贾的主意。不过吃一堑长一智，孟导还是想留点什么给自己提个醒。于是就和老贾商量能不能留五枚伪币给自己带回去。老贾起初还有些不乐意，不过想到孟导这么爽快答应了自己出的价格，也就没吱声了。

1.2.5 数据处理 采用ABI公司的GeneMapper4.0软件，选择GeneScanTM 500 ROXTM作为分子量标准，对每一个扩增出的条带记录大小。运用PopGen3.2软件计算Shannon信息指数(I)和Nei’s基因多样性指数(H)。运用NTSYSpc2.1中的Tree plot程序进行聚类分析。

1.2.4 群体结构分析 使用Struture 2.3.4[22]进行群体结构分析，先计算各材料归属第k亚群的概率Q值。使用混合模型K值从3设置到10模拟群体结构，每一个K值重复5次。将MCMC迭代设为100000次，初始burn-in次数为100000。

引物统一在反向引物(R)的5′段标记荧光(FAM/TAMRA)，由上海百力格生物科技有限公司合成。扩增产物0.5 μL，GeneScanTM 500 ROXTM 0.5 μL，HiDi 9 μL，混匀后使用ABI公司3730XL进行毛细管电泳。

表2 17对 SSR 引物及其序列

Table 2 Nucleotide sequences of 17 primer pairs

引物名称重复序列正向引物序列反向引物序列TM241(GAGAA)3ATCGGCGACGGTGGAAGTGCCAGCGGAGAGGAGAAGTM262(CT)21CGACCAGACGGTGAAATAGGCTTGTGAGCAAAATCTM341(TA)10CATGCTCCCATCCCACCTATGCTGCTCATTCAAACCAACTTM345(CT)12TTCTTGTCCTTCATGGAACTCAGGAAATCTTATGGCACGGAATM369(GAA)8CGGAGCTGGAATCTGAAGAGGGAAGGGTTGCAAATTCTGATM422(TTC)7GGACTTCGTTGCTTCCTTTGCCATTCTCGACGAATCCAGTTM426(AGA)11TGAGAGTGCTTGTCTGGGTGCAACTACCCCTTTTCCCCATTM428(CAC)7TCTCCTCCTCGATCCTCAGACCCTCTTCTTCGGATCCTTCTM440(TTTGC)3TTGACCCGAATAAAATGGGACCTCAAAACATGCTTTTCTTAATCTM442(ATACAC)3CAAGCCAAACCTTGCTGAATCTGTCCTGTGTCTGGTGGTGTM447(AAAAG)5TGTTGTTAACGGTGTTCGGAGCATTTGTTTTCTCTCTCTGCCTM461(ATTTTT)6GGCTAGGGTTTCTCCCACTTGAAGGTCGAAGCGATGTTGTTM499(AGA)5AACTGTGACACCGATTGCAGAAGTTTCACTTGCCAGCACCTM513(AG)10CAAGCGATCAACAACAATGGTTGAGAAATCAACCCCTTGGTM514TM514(TCA)5ATGTCTGGCCGTGGATTAAGTM569(TTTTTG)3TTGCGATAACATCCCAGTGATTTTCATTTGCCATCACCAATM601(GGA)5TTGCACTGGAGTGCGATAAGCATCGCCACCAAACTCTTCT

2 结果与分析

2.1 94个材料间亲缘关系分析

1.2.3 PCR扩增和产物鉴定 PCR反应体系为：ddH2O 18.8 μL，10×Buffer 2.5 μL(Mg2+)，dNTP(10 mmol·L-1)0.5 μL，上、下游primer(10 μmol·L-1)各0.5 μL，Taq polymerase 0.2 μL(0.5U)，模板DNA 1μL。PCR反应于美国ABI-9600型扩增仪上进行，热循环程序为：94℃预变性4 min，使模板DNA 充分变性，然后进入下列温度循环：94℃变性45 s，不同温度条件退火60 s，72℃延伸75 s，重复35个热循环；72℃延伸10 min，最后4℃保存。

2.2 群体间遗传多样性比较与亲缘关系分析

采用Struture 2.3.4软件对94个茶树品种(系)进行系统解析，当K=3时△K值为最大值，从图3中看出，94个茶树品种(系)分为3个类群。参考刘志斋[24]等以类群属性比率50%为标准的划分依据，94个茶树品种(系)分为S1、S2、S3类群，S1类群从品种来源看，代表性材料为黄棪、铁观音，主要材料为金观音及其自然杂交后代、金牡丹及其自然杂交后代，体现的是黄棪、铁观音等闽南高香资源类群属性；S2类群从品种来源看，代表性材料为白鸡冠、福云6号、福鼎大白茶、悦茗香，主要材料为白鸡冠、福云6号、紫玫瑰与悦茗香自然杂交后代，体现的是乌龙茶与绿茶资源的混合类群属性；S3类群从品种来源看，代表性材料为肉桂、矮脚乌龙，主要材料为武夷山茶树种质资源，体现闽北优异乌龙茶种质资源属性(详见表4)。

图1 遗传相似系数聚类图 Fig.1 Dendrogram on genetic similarity

Hardy-Weinberg遗传偏离指数D主要反映Ho(观测杂合度)和He(期望杂合度)之间的平衡关系，D值越接近0，基因型的分布就越接近于平衡状态，D值大于0说明存在杂合子过剩，D值小于0，说明存在杂合子缺失现象[23]。从表3可以看出，6个自然杂交后代群体Ho高于对应的He，D值大于0，属于亲本单一，杂合子过剩现象；武夷山2个地方种质资源群体Ho低于对应He，D值小于0，属于试验材料有限杂合子缺失现象。

从群体聚类图2看出，群体1、群体6、群体7与群体8聚类在一起，说明这4个群体亲缘关系相对较近，悦茗香自然杂交后代群体的父本可能来自于闽北武夷山茶树种质资源；群体2、群体9、群体3与群体4聚类在一起，推测福云6号自然杂交后代群体的父本可能来自于闽南铁观音、黄棪等高香茶树种质资源。

2.3 94个材料群体结构分析

从表3中看出，Shannon遗传多样性指数最高的为群体7、群体8(武夷山茶树种质资源)，且高于群体9(福建省茶树品种)；6个自然杂交后代群体遗传多样性指数最高的是福云6号自然杂交后代，推测是由于母本福云6号有云南茶树种质资源遗传基础，遗传多样性相对广泛；白鸡冠、悦茗香、紫玫瑰的自然杂交后代群体遗传多样性相当，略高于金观音、金牡丹的遗传多样性指数。

结合我国农村的实际情况，建立信息数据库。加大资金投入力度，充分发挥政府的职能，及时将最新的数据信息提供给广大农村地区[2]。要强化数据库管理，保证农业信息记录的正确性，确保入库信息真实、准确，实现农业信息共享。同时，在实际工作中要合理运用农村信息管理技术，保证农民能够在第一时间接收到相关信息，最大程度提高农村经济管理的工作效率。

表3 9个群体遗传多样性

Table 3 Comparison on genetic diversity of 9 groups

群体名称观测等位基因数有效等位基因数Shannon多样性指数观测杂合度(Ho)期望杂合度(He)基因多样度群体13.6472.6580.9950.6350.5730.532群体24.4123.0631.1860.7510.6550.625群体32.2941.9460.6770.5350.4570.436群体43.0002.1530.8260.5940.5220.499群体53.2942.5280.9530.6080.5760.549群体64.4122.5231.0380.5880.5570.532群体74.5883.2161.2550.6170.6910.657群体84.5883.0471.2190.5700.6680.638群体94.0592.9351.1360.5650.6440.612平均值3.8102.6741.0320.6070.5940.564

图2 群体遗传相似系数聚类图 Fig.2 Dendrogram on genetic similarities of groups

图3 假设K=3时群体按照所估计的成分(Q)值的结构分析结果 Fig.3 Membership coefficients of 3 clusters based on assumption of cluster number (k=3) and accessions sorted by individuals

表4 94个材料的群体结构

Table 4 Population structure of 94 germplasms

第一类群铁观音、黄棪第二类群白鸡冠、福云6号、福鼎大白茶、悦茗香第三类群肉桂、矮脚乌龙代号组成1组成2组成3代号组成1组成2组成3代号组成1组成2组成3A50.55460.26280.1827A10.01540.97210.013F30.0080.04540.9466B60.61210.01840.3695A20.0270.94840.025G10.03110.01650.9524B90.69490.29510.01A30.01060.97530.014G100.0090.01140.9796C10.9010.03040.0686A40.04970.93940.011G110.01070.01190.9774C100.98780.00640.0059A60.0710.91530.014G20.03130.03460.9341C110.96740.01950.0131A70.04390.9250.031G30.0060.01410.9799C20.98560.00840.006B10.28890.70370.007G40.01220.0140.9738C30.98690.0070.0061B100.38390.53440.082G50.03330.01330.9535C40.98730.00690.0059B110.01150.97890.01G60.00810.0090.9829C50.98690.00710.006B20.19070.77140.038G70.04150.0420.9165C60.98810.00690.005B30.11440.86780.018G80.01110.01660.9722C70.98780.0070.0052B40.1970.77790.025G90.0070.0150.978C80.98770.00710.0051B50.32380.64970.027H10.02750.01950.953

(续表4)

第一类群铁观音、黄棪第二类群白鸡冠、福云6号、福鼎大白茶、悦茗香第三类群肉桂、矮脚乌龙代号组成1组成2组成3代号组成1组成2组成3代号组成1组成2组成3C90.98540.00760.007B70.38050.57460.045H100.01010.010.9799D10.57590.41090.0131B80.40320.59020.007H110.0070.01150.9815D110.95390.03470.0114D100.36220.55420.084H20.02610.01550.9584D20.98650.00750.006E10.15740.80540.037H30.0080.01040.9816D30.98530.00790.0069E100.04040.94590.014H40.01250.08860.8989D40.98550.00850.006E20.05540.93090.014H50.01480.01790.9674D50.98550.00840.0061E30.17420.56510.261H60.01140.01590.9727D60.97340.01820.0084E40.15630.82420.02H70.0120.04350.9445D70.98460.00930.0061E50.15710.83070.012H80.0070.01420.9788D80.98590.0080.0061E60.17270.80780.02H90.00930.15720.8335D90.9840.00910.0069E70.05210.92940.019L100.46150.02040.5181E110.68540.18570.1289E80.33910.64090.02L20.11750.18890.6937F70.58810.24860.1633E90.03440.95540.01L90.14610.34850.5054L110.76330.19550.0412F10.02180.83820.14L80.46680.02560.5076L30.58060.30430.1151F100.0070.92750.066平均值0.057290.04560.8971074L40.94260.03360.0237F110.01250.70370.284L50.53790.1320.3301F20.0110.96440.025平均值0.87040.0770930.0525F40.00990.9790.011F50.01190.96450.024F60.0080.91490.077F80.0190.95090.03F90.03080.9430.026L10.01650.89940.084L60.01050.93980.05平均值0.11820.8360.05

3 讨论

3.1 武夷山种质资源初步鉴定

武夷山茶树种质资源丰富，据历史记载有1187种，其中铁罗汉、水金龟、白鸡冠、半天妖等四大名枞就是从武夷山的菜茶中选出[25]，武夷菜茶在世界茶树植物学分类学上占有重要地位，被称为武夷种、武夷亚种或武夷变种[26]。刘振[27]、朱晨[28]、姜燕华[29]分别对福建茶树种质资源进行遗传多样性分析，Shannon多样性指数平均值分别为1.16、0.368、0.877，本试验中2个武夷山茶树种质资源群体的Shannon多样性指数平均值分别为1.219、1.136，参试武夷山茶树种质资源群体遗传多态性相对最丰富；武夷山种质资源2个群体间没有明显的差异性，群体亲缘关系较近；22个资源主要分为2个类型，第1个类型9个材料与肉桂亲缘关系较近，第2个类型5个材料与大叶乌龙(L2)、梅占(L5)、矮脚乌龙(L8)、白芽奇兰(L10)构成的小群体聚类在一起，遗传组成较相似；其中，武夷山种质资源H7与所有参试材料都未聚类在一起，推测参试材料H7遗传基础较为特殊，需进一步观察与试验。

项目建设范围内主要为荒地、房屋、鱼塘和山地，地面标高为29.95～12.56m（黄海高程、下同），地势起伏较大，呈北高南低的趋势。沿线无任何市政管线，雨污水处于散排状态，通过地面径流就近排入附近鱼塘、洼地，造成污染。

3.2 6个自然杂交后代初步鉴定

白鸡冠自然杂交后代遗传基础与白鸡冠多相一致；紫玫瑰自然杂交后代未与紫玫瑰聚类在一起，推测其遗传基础受父本影响较大；悦茗香自然杂交后代遗传组成与悦茗香一致，群体与武夷山地方资源亲缘关系较为接近，推测悦茗香自然杂交后代的父本是武夷山种质资源类型；福云6号自然杂交后代分为2种类型，第1个类型3个材料与福云6号聚类在一起，遗传组成偏向绿茶，另外6个材料与铁观音、黄棪等高香乌龙茶品种聚类在一起，推测其父本为铁观音、黄棪等高香乌龙茶品种资源；金观音、金牡丹及自然杂交后代与铁观音、黄棪属于1个类型。

参考文献

[1]虞富莲.论茶树原产地和起源中心[J].茶叶科学,1986,6(1):1-8.

[2]Senior M L, Murphy J P, Goodman M M, et al. Utility of SSRs for determining genetic similarities and relationships in maize using an agarose gel system[J]. Crop Sciences,1998(38):1088-1098.

[3]Wang F G, Tian H L, Zhao J R, et al. Development and characterization of a core set of SSR markers for fingerprinting analysis of Chinese maize varieties[J]. Maydica, 2011(56):1693-1699.

[4]Würschum T, Langer S M, Longin C F H, et al. Population structure, genetic diversity and linkage disequilibrium in elite winter wheat assessed with SNP and SSR markers[J]. Theoretical and Applied Genetics, 2013,126(6):1477-1486.

[5]Adugna A. Analysis of in situ diversity and population structure in Ethiopian cultivated Sorghum bicolor (L.) landraces using phenotypic traits and SSR markers[J]. Springerplus, 2014(3):212.

[6]姜燕华,张成才,成浩,等．茶树良种场不同品种的SSR鉴定研究[J]．茶叶学报,2016,57(3):105-112.

[7]王松琳,马春雷,陈亮,等．基于SSR标记的白化和黄化茶树品种遗传多样性分析及指纹图谱构建[J]．茶叶科学,2018,38(1):58-68．

[8]雷雨,段继华,李赛君,等.茶树新品系“53-34”父本的EST-SSR标记鉴定?[J].茶叶通讯,2016,43(1):28-31．

[9]杨阳,刘振,赵洋,等.利用 EST-SSR 标记研究黄金茶群体遗传多样性及遗传分化[J]．茶叶科学,2009,29(3):236-242．

[10]黄晓霞,唐探,姜永雷,等.千家寨不同海拔野生茶树的EST-SSR遗传多样性研究[J].茶叶科学,2015,35(4):347-353．

[11]毛娟,江鸿键,陈亮,等.云南白莺山地区茶树遗传多样性研究[J].茶叶科学,2018,38(1):69-77．

[12]张明泽,姚玉仙,陈世军,等．黔南60份茶树种质资源遗传多样性的SSR分析[J].西北植物学报,2016,36(6):1117-1124.

[13]郭燕,刘声传,曹雨,等.基于SSR标记贵州古茶树资源的遗传多样性分析及指纹图谱构建[J].西南农业学报,2016,29(3):491-497.

[14]周萌,李友勇,刘本英,等.基于EST-SSR分子标记对香竹箐茶树王的遗传多样性分析[J].西南农业学报,2016,29(2):231-239.

[15]姜晓辉,陈栋,吴华玲,等．基于EST-SSR毛细管电泳荧光标记技术分析广东2个历史名茶群体遗传多样性[J].热带作物学报,2018,39(1):46-54.

[16]王丽鸳，姜燕华，段云裳,等．基于SSR 分子标记的龙井群体种的遗传多样性及遗传分化研究[J].茶叶科学,2011,31(1):40-44．

[17]姚明哲,刘振,梁月荣,等.利用EST-SSR分析江北茶区茶树资源的遗传多样性和遗传结构[J]．茶叶科学,2009,29(3):243-250．

[18]乔婷婷,马春雷,陈亮,等．浙江省茶树地方品种与选育品种遗传多样性和群体结构的EST-SSR分析[J].作物学报,2010,36(5):744?753．

[19]吴晓梅,姚明哲,陈亮,等．利用EST-SSR标记研究适制绿茶与乌龙茶品种的遗传多样性与遗传结构[J].茶叶科学,2010,30(3):195-202．

[20]金基强，崔海瑞，龚晓春,等．用EST-SSR标记对茶树种质资源的研究[J]．遗传，2007，29(1)：103-108．

[21]Ma J Q, Yao M Z, Ma C L, et al. Construction of a SSR-based genetic map and identification of QTLs for catechins content in tea plant(Camellia sinensis)[J].PLOS ONE,2014,9(3):e93131.

[22]Prichard J K, Stepphen M, Donnelly P. Inference of population structure using multilocus genotype data[J]. Genetics, 2000(155):945-959.

[23]刘志斋,吴迅,刘海利,等．基于40个核心SSR标记揭示的820份中国玉米重要自交系的遗传多样性与群体结构[J].中国农业科学,2012,45(11):2107-2138．

[24]宋炜，孟永永，马凌波,等．棘头梅童鱼七个野生群体遗传多样性的微卫星分析[J]．水产学报, 2017，41(1):31-39.

[25]罗盛财.武夷岩茶名枞录[M].北京:科学出版社,2007:15-63.

[26]白堃元.中国茶树品种志[M].上海:上海科学出版社,2001:4-6.

[27]刘振,姚明哲,陈亮,等．基于EST-SSR的福建地区茶树资源遗传多样性和亲缘关系分析[J].中国农业科学,2009,42(5):1720-1727.

[28]朱晨，赖钟雄，郭玉琼,等．福建省 53 份茶树种质资源遗传多样性ISSR分析[J]．热带作物学报,2017，38(7):1303-1310.

[29]姜燕华,段云裳,成浩,等．福建茶树品种的SSR分析及其人为选择的影响[J].浙江林业科技,2010，30(3):12-16.