李力，羅盛財，王飛權(quán)，黎巷汝，馮花，石玉濤，葉江華，劉菲，趙佳林，李舒瑩，張渤*

1. 武夷學院茶與食品學院，福建 南平 354300；2. 武夷山市農(nóng)業(yè)農(nóng)村局，福建 南平 354300；3. 武夷學院茶葉科學研究所，福建 南平 354300；4. 福建農(nóng)林大學園藝學院，福建 福州 350002

武夷茶是山茶科（Theaceae）山茶屬（Camellia）的一類成員，以“武夷”的英文音譯而被命名為Bohea或VarBohea（Thea BoheaLinnaeus. 1762）。武夷茶生長于中國福建省北部的武夷山，武夷山被認為是烏龍茶和紅茶的發(fā)源地[1-2]。武夷山地區(qū)的先民從有性系的武夷茶樹群體中篩選分離出品質(zhì)優(yōu)良的稀有單叢單株，通過無性繁殖傳承至今，成為現(xiàn)在著名的“武夷名叢”[3]。由于武夷茶樹生長在巖石土壤中，其鮮葉制作的茶葉具有“巖韻”（風味濃郁，香味持久）的品質(zhì)特點，被稱為武夷巖茶，是中國十大名茶之一，具有巨大的經(jīng)濟和文化價值[4]。

歷史上記載的武夷茶樹品種有上千種。然而，許多種質(zhì)資源不斷消失，通過對茶樹性狀的長期觀察，目前僅收集到百余個種質(zhì)資源[5-6]。此外，武夷山當?shù)氐挠N過程中存在對優(yōu)良品種過度選擇的情況，將導致群體遺傳多樣性越漸狹窄，最終影響品種性狀的綜合改良，再加上存在茶樹同物異名的現(xiàn)象，進一步阻礙了武夷茶樹種質(zhì)資源的開發(fā)。因此，對武夷茶的遺傳背景分析和品種鑒定是一項重要的工作。傳統(tǒng)的茶樹種質(zhì)資源分類基于形態(tài)或農(nóng)藝性狀，容易受到環(huán)境條件的限制[7]，這給品種的確定帶來了一定的困難[8]。國際植物新品種保護聯(lián)盟（UPOV）對作物品種建立了“差異性、統(tǒng)一性和穩(wěn)定性”測試制度[9]，在難以通過表型區(qū)分確認新品種的情況下，DNA 指紋圖譜可以幫助提高品種的識別。新修改的《中華人民共和國種子法》于2022 年3 月1 日正式施行，首次建立實質(zhì)性派生品種（Essentially derived variety，EDV）制度[10]，要求對品種的基因或基因型組合進行特異性檢測鑒定。建立真實可靠的DNA 鑒定檢測技術(shù)對新《種子法》修訂條例的實施具有重大戰(zhàn)略意義。

目前，茶樹的特異性DNA 鑒定手段仍相對落后于主要農(nóng)作物，早期的茶樹特異性DNA 鑒定方法多采用第一代或第二代分子標記，包括擴增片段長度多態(tài)性（AFLPs）、限制性片段長度多態(tài)性（RFLP）、隨機擴增多態(tài)性DNA（RAPD）、簡單序列重復序列（SSRs）等[11-14]，武夷茶樹也是如此[15-17]。雖然這些傳統(tǒng)方法各有優(yōu)勢，但所使用的標記數(shù)據(jù)量有限，不適合大樣本量的群體分析，不能充分反映物種的遺傳多樣性。相比之下，第三代分子標記技術(shù)單核苷酸多態(tài)性（Single nucleotide polymorphism，SNP）[18-20]具有一定的優(yōu)勢，已逐步成為主流的分子標記[21]。Lin 等[22]采用高通量SNP 技術(shù)構(gòu)建了4 個主要烏龍茶產(chǎn)區(qū)的DNA 指紋圖譜，揭示了4 個產(chǎn)區(qū)的遺傳關(guān)系。Liu 等[23]從茶樹的表達序列標簽（EST）數(shù)據(jù)庫中篩選出49 個SNP 位點，分析了武夷山及其鄰近地區(qū)137 份茶樹種質(zhì)的遺傳多樣性豐富程度及遺傳關(guān)系。這表明高通量SNP技術(shù)在尋找和驗證茶樹親本關(guān)系方面具有很大的發(fā)展空間。然而，茶樹的SNP 分子標記開發(fā)發(fā)展較晚，數(shù)量有限，仍需不斷完善補充。

傳統(tǒng)的SNP 篩選方法需要大量的人力資源和較長的時間，基因分型（Genotyping by sequencing，GBS）測序技術(shù)可以快速識別高密度SNP 位點，廣泛應用于作物種質(zhì)鑒定、群體結(jié)構(gòu)分析和全基因組關(guān)聯(lián)分析[24-25]。本研究收集了126 份武夷茶品種/品系，結(jié)合來自12 個省份的223 份優(yōu)異茶樹種質(zhì)資源，共349份茶樹樣本進行GBS 測序，利用高質(zhì)量SNPs對武夷茶進行遺傳多樣性及背景分析，并開發(fā)一種簡便、快速的茶葉SNP 分子標記鑒定，旨在為茶樹研究領(lǐng)域EDV 制度的有效實施提供參考。

1 材料與方法

1.1 樣本收集與DNA 提取

樣品來自武夷山市龜巖茶樹種質(zhì)資源保護基地（福建省優(yōu)異茶樹種質(zhì)資源保護區(qū)閩茶圃004，N 27°43'42.46"，E 118°0'14.40"）和武夷學院茶樹種質(zhì)資源圃（N 27°61'01.34"，E 117°96'63.51"），包括126 個武夷茶品種/品系及223 個來自12 個省份的茶樹品種/品系，共349 個（包含同一品種/品系的6 對）。樣本信息見附表1（掃文后二維碼查看附表1）。

新鮮葉片基因組 DNA 采用植物基因組DNA 提取試劑盒（康威CW0553S）提取，DNA濃度和純度采用 NanoDrop 紫外分光亮度計（Thermo Scientific，USA）測定。DNA 完整性使用0.8%瓊脂糖電泳檢測，基因組DNA 保存于-80 ℃。

1.2 基因分型測序及建庫

GBS 庫的構(gòu)建參照文獻[26]，用限制性內(nèi)切酶（EcoRⅠ和NiaⅢ）酶解質(zhì)檢合格的基因組DNA，在酶切后的片段兩端加上相應的適配接頭，利用PCR 擴增兩個接頭之間的片段，共構(gòu)建了349 個對應樣本的多重庫，每個庫DNA 樣本都有唯一的適配接頭。將樣品混合后進行電泳回收純化，純化產(chǎn)物質(zhì)檢合格后使用Illumina NovaSeq 6000 測序平臺（上海凌恩生物科技有限公司）進行雙末端測序。

測序獲得的原始序列數(shù)據(jù)（Raw read）去除低質(zhì)量序列，有效的高質(zhì)量序列（Clean read）數(shù)據(jù)經(jīng)BWA v0.7.10 軟件[27]比對茶樹基因組[28]，利用SAMTOOLS v1.9 軟件將sam 轉(zhuǎn)換為bam 文件。GATK v4.1.2.0 軟件用于檢測SNP，初始SNP 質(zhì)量過濾標準為“QD＜2.0 ||MQ ＜40.0 || FS ＞60.0 || SOR ＞3.0 ||MQRankSum ＜-12.5 || ReadPosRankSum ＜-8.0”。使用VCF v0.1.11 軟件對SNP 進行嚴格過濾[29]，經(jīng)過測序深度4X，缺失率（Miss rate）＜10%、次等位基因頻率（MAF）＞0.05，多態(tài)信息含量（PIC）＞0.15，SNP 位點兩側(cè)150 bp 無變異的過濾，最終共獲得12 937 679個SNP。從12 937 679 個SNP 中篩選出遺傳多樣性較高且均勻分布在15 個染色體上的核心SNP（去除重組率為0 的冗余位點、MAF＞0.15、Miss rate＜5%、PIC＞0.15）用于后續(xù)分析。使用PowerMaker v3.25 軟件[30]計算觀察到的雜合度（Observed heterozygosity，Ho）值和基因多樣性（Gene diversity，GD）值。

1.3 茶樹群體遺傳多樣性及結(jié)構(gòu)分析

使用Admixture 軟件分析種群結(jié)構(gòu)[31]。根據(jù)隸屬度系數(shù)（Q）的取值對所有品種進行劃分，根據(jù)交叉驗證誤差（CV 誤差）確定最優(yōu)K 值。主成分分析（PCA）使用Plink 軟件進行，使用 R v3.4 進行可視化繪圖[32]。利用MEGA 軟件的鄰接法（ Neighbor-joining method，NJ）構(gòu)建系統(tǒng)發(fā)育樹，軟件分析均采用默認參數(shù)設置，使用iTOL 在線軟件編輯系統(tǒng)發(fā)育樹[33]。

1.4 基因流分析

根據(jù)茶樹文獻資料[34-35]按茶樹原產(chǎn)地（以省/直轄市為單位）對349 個茶樹進行分組，為探索武夷山地區(qū)茶樹的歷史遺傳基因流，進一步將福建省的地區(qū)劃分為福建武夷山區(qū)（FJW）、閩東區(qū)（FJE）、閩南區(qū)（FJS）和福建未知區(qū)（FJU）。為減少誤差，少于5 個品種的地區(qū)組被排除分析。最后，共采用11地區(qū)組進行分析（附表1）。使用TreeMix v1.01軟件進行基因流分析[36]，分別設定假設有1～11 次基因流事件（-m=1～11），以箭頭表示不同組之間的基因滲入。

1.5 遺傳相似度分析

使用軟件NTSYSpc v2.11[37]進行遺傳相似度分析，計算SNP 的分布頻率和遺傳相似系 數(shù)（ Genetic similarity ， GS） ，GS=NS/(NS+ND)。式中，NS為相同SNP 基因型數(shù)量，ND為不同SNP 基因型數(shù)量[38]。缺失的基因型被視為無效基因型。

根據(jù)國際種子聯(lián)合會[39]的建議，GS 值大于0.9 可作為實質(zhì)派生關(guān)系的有力證據(jù)。因此，GS＞0.9 的品種對被認為是近似品種，具有派生關(guān)系。349 個茶樹中有6 對來自不同生長地的同一品種/品系，作為明確為實質(zhì)派生關(guān)系的陽性對照，將其中最小的GS 值作為具有無可爭議的派生關(guān)系的閾值。當茶樹樣本的GS 值大于該閾值，將這對樣本明確為具有派生關(guān)系。

1.6 簡易鑒定SNP 開發(fā)

在973 個核心SNP 的基礎上，使用Perl方法以349 個樣本成對比較的可辨性為過濾條件，選擇一組數(shù)量最少且可辨性高的SNP標記用于品種簡單快速鑒定[40-41]。在SNP 側(cè)翼兩個保守區(qū)域設計引物，隨機進行一代測序驗證。

2 結(jié)果與分析

2.1 全基因組SNP 篩選

349 份茶樹種質(zhì)資源共獲得829.05 G 高質(zhì)量測序數(shù)據(jù)，平均每個樣品獲得2.38 G 數(shù)據(jù)，Q20（堿基錯誤率在1%以下）平均占比98%，Q30（堿基錯誤率在0.1%以下）平均占比93.05%，平均獲得高質(zhì)量序列條數(shù)5 737 108 966，平均有5 668 097 099 條序列可以匹配到茶樹參考基因組，平均覆蓋率為98.83%，表明測序質(zhì)量較高。經(jīng)過嚴格過濾篩選（Miss rate＜20%、MAF＞10%），共獲得12 937 679 個SNP。測序獲得的原始數(shù)據(jù)上傳 NCBI 數(shù)據(jù)庫（項目號：PRJNA924950）。以去除重組率為0 的冗余位點、MAF＞0.15、Miss rate＜5%、PIC＞0.15 且較均勻分布在茶樹基因組15 條染色體上為篩選條件，共獲得973 個核心SNPs。

2.2 茶樹群體遺傳多樣性分析

973 個SNP 較均勻地分布在茶樹的15 條染色體上，在第1 至第15 條染色體上的個數(shù)分別為94、75、62、73、67、63、70、59、69、60、47、44、64、68 和58（圖1A）。在349 個茶樹中的Ho 值為0.355，53.04%以上Ho 值大于0.3；GD 平均值為0.319，51.28%以上GD 值大于0.3；茶樹群體PIC 平均值為0.262。973 個SNP 在武夷茶樹品種/品系中的Ho 值為0.334，GD 平均值為0.302，PIC 平均值為0.248（圖1B）。

圖1 973 個核心SNP 在茶樹基因組15 條染色體上分布情況及多樣性分析Fig. 1 Distribution and diversity analysis of 973 core SNPs on 15 chromosomes of tea plant genome

基于973 個SNP 位點對349 個茶樹進行群體結(jié)構(gòu)分析?；谀Ｐ偷娜后w結(jié)構(gòu)分析表明，5 個群體（K=5）為最佳模型（圖2A）。通過PCA 對群體結(jié)構(gòu)進行評價，得出了5 個聚類，與K=5 時的結(jié)構(gòu)分析推斷結(jié)果一致（圖2B）。NJ 系統(tǒng)發(fā)育樹與種群結(jié)構(gòu)分析及PCA結(jié)果基本吻合，但少部分個體在分支聚類上有所差異（圖2C）。亞群1～5 分別有茶樹個體159、72、51、50 個和17 個。亞群1 主要以福建省武夷山地區(qū)茶樹、浙江地區(qū)茶樹及安徽地區(qū)茶樹為代表，亞群2 主要是以福建省閩南地區(qū)鐵觀音和黃旦等茶樹為代表，亞群3 主要是以福鼎大白茶及其相關(guān)品種類型為代表，亞群4 主要是以湖南地區(qū)茶樹為代表，亞群5主要是廣東地區(qū)茶樹為代表。武夷山茶樹資源大部分集中于亞群1，共108 個，15 個處于亞群2，1 個處于亞群3，2 個處于亞群4。

2.3 基因流分析

在預估的1～11 次基因流事件中（11 個場景），共有5 個場景涉及武夷山地區(qū)茶樹的基因流動。其中2 個場景發(fā)現(xiàn)武夷山地區(qū)茶樹群體（FJW）的基因主要向浙江地區(qū)茶樹群體（ZJ）流動（圖3B 和3C），沒有發(fā)現(xiàn)外地省份的茶樹群體直接向武夷山地區(qū)茶樹群體的基因滲入情況，但3 個場景發(fā)現(xiàn)存在本省閩南地區(qū)（FJS）茶樹群體向武夷山地區(qū)（FJW）的基因滲入（圖3D～F）。

圖3 基因流分析Fig. 3 Gene flow analysis

2.4 遺傳相似性分析

基于973 個SNP 位點進行349 個品種之間遺傳相似度分析（圖4），經(jīng)過兩兩比較，GS 在0.524 0～0.996 9，平均值為0.661 8。其中，GS 在0.6～0.7 的最多（90.60%），為55 018 對；GS 在0.8～0.9 的最少，為80 對；GS＞0.9 的共有136 對。126 個武夷茶品種/品系（含6 對已知實質(zhì)性派生關(guān)系的茶樹）之間的GS 在0.617 7～0.996 9，平均值為0.686 0。GS＞0.9 的共24 對，在0.8～0.9 的共7 對，在0.7～0.8 的共831 對，在0.6～0.7 的最多（89.05%），共7 013 對。

圖4 349 個茶樹之間的遺傳相似度分析Fig. 4 Genetic similarity analysis among 349 tea resources

349 個茶樹中涉及武夷茶品種/品系GS＞0.9 的有26 對（含已知具有派生關(guān)系的6 對），被認為是近似品種，具有派生關(guān)系，以6 對來自不同生長地的同一品種/品系的兩株茶樹作為具有明確派生關(guān)系的陽性對照，其中GS 最小的為0.972 3，涉及武夷茶品種/品系遺傳相似系數(shù)大于0.972 3 的共有23 對（含已知具有派生關(guān)系的6 對），被明確為具有無爭議的派生關(guān)系（表1）。其中FJ84（JM051 肉桂）、FJ204（JM046 紅海棠）、FJ73（JM037 王母桃）、FJ200（JM035 金雞母）、FJ191（黃肉桂）5 個品種/品系之間遺傳相似系數(shù)均大于0.972 3，明確具有派生關(guān)系。大紅袍2 號株FJ91（JM062-2 大紅袍2）與FJ192（奇丹），大紅袍4 號株FJ93（JM062-4 大紅袍4）與FJ82（JM049 不知春），F(xiàn)J108（JM081 醉墨）與FJ48（JM005 雀舌），F(xiàn)J67（JM029 向天梅）與FJ71（JM034 醉貴姬），F(xiàn)J76（JM042胭脂柳）與FJ77（JM043 醉八仙），F(xiàn)J13（1113）與FJ12（1114），以及FJ86（JM054 金毛猴）與FJ149（白毛猴）相似系數(shù)均大于0.972 3，明確具有派生關(guān)系。而FJ52（JM011 石中玉）與FJ161（0205D），F(xiàn)J63（JM025 石觀音）與FJ187（SR），F(xiàn)J17（0204）與FJ18（0205）相似系數(shù)在0.90～0.97，被認為是近似品種/品系，具有一定派生關(guān)系。以上品種/品系對之間均可被視為分別屬于同一類型。

表1 126 個武夷茶樹中遺傳相似系數(shù)大于0.9 的品種/品系對及派生關(guān)系分析Table 1 The information of cultivar/strain pairs with genetic similarity greater than 0.9 and essentially derived relationship analysis in 126 Wuyi tea cultivars

2.5 簡易鑒定茶樹品種SNP 的篩選

在973 個核心SNPs 基礎上，通過349 個樣本的成對比較，選擇出具有高辨識度的21個SNPs 可用于區(qū)分349 個茶樹（圖5A），其中18 個SNPs 即可用于126 個武夷茶品種/品系建立DNA 指紋圖譜（圖5B）。進一步在21個SNP 位點兩側(cè)保守區(qū)域設計引物（表2），使用21 對引物對隨機抽取的10 個茶樹樣本進行PCR 與一代測序驗證，分型結(jié)果均顯示與位點一致，可用于茶樹種質(zhì)資源的快速鑒定（圖6）。

圖5 349 個茶樹（A）及126 個武夷茶（B）的SNP 指紋圖譜Fig. 5 SNP fingerprints of 349 tea cultivars (A) and 126 Wuyi tea cultivars (B)

表2 21 對可區(qū)分349 個茶樹樣本的引物信息Table 2 Information of 21 pairs of primers that can distinguish 349 tea resources

3 討論

3.1 GBS-SNP 及茶樹遺傳多樣性

本研究收集了126 個武夷茶品種/品系及223 個來自其他不同地區(qū)的茶樹種質(zhì)資源，共349 個茶樹樣本，利用GBS-SNPs 方法分析了武夷茶的基因遺傳多樣性與背景，共鑒定出12 937 679 個高質(zhì)量SNPs，所有樣本均通過質(zhì)量評估，平均準確率達98.8%。高通量和高質(zhì)量的SNP 表明，GBS 是進行種群遺傳多樣性分析的有效方法。本研究進一步篩選出多樣性較高（缺失率＜5%、MAF＞0.15、PIC＞0.15）并均勻分布在15 個染色體上的973 個SNP 位點進行遺傳多樣性分析。973 個SNPs 在349個茶樹樣本中雜合度Ho 平均值（0.355）和GD 平均值（0.319）均較高，表明這973 個SNP信息豐富，具有較高的基因多樣性，可用于遺傳多樣性分析。973 個SNPs 在349 個茶樹與126 個武夷茶樹中的PIC 平均值分別為0.262與0.248，茶樹群體多態(tài)性信息含量中等，這暗示著茶樹群體中存在一定數(shù)量親緣關(guān)系較近的茶樹。以12 937 679 個高質(zhì)量SNPs 和973個SNPs 構(gòu)建的群體遺傳結(jié)構(gòu)、主成分和系統(tǒng)發(fā)育關(guān)系均表明，349 個茶樹樣本可劃分為5個種群。這5 個茶樹種群主要基于茶樹之間的親緣關(guān)系進行聚類，而不是樹型或葉形等形態(tài)特征，這與以往的研究一致[42]。

圖6 隨機抽取10 個茶樹樣本進行一代測序的分型結(jié)果示例Fig. 6 The sequencing and genotyping result of the random sampling of 10 tea plant samples

3.2 武夷茶歷史遺傳背景與基因流

中國西南地區(qū)是茶樹的發(fā)源地[28]，同源茶種以不同的傳播演化路徑發(fā)展。早前有學者根據(jù)茶樹的生態(tài)型和地質(zhì)變化提出茶樹五大演化區(qū)[43]，即黔區(qū)系群、三江區(qū)系群、南嶺區(qū)系群、武夷區(qū)系群及江漢區(qū)系群，指出茶樹在經(jīng)歷自然傳播后，因地質(zhì)和氣候變化在第四紀冰期后形成隔離分布和演化。武夷山地處福建西北，在中生代白堊紀時候由海洋逐漸變?yōu)殛懙?，到第四紀末冰期地殼上升，產(chǎn)生褶皺和斷裂成山間谷地。武夷山的地質(zhì)年代，形成土壤的巖性以及存在野生茶等一系列特點，暗示武夷山茶樹是由“隔離分布”的同源演化成“武夷類群”[44]。本研究中基因流分析未發(fā)現(xiàn)其他省份茶樹直接向武夷山茶樹群體的基因流動，其結(jié)果符合該說法。此外，基因流分析顯示武夷茶主要是向浙江地區(qū)茶樹基因流動，這與種群結(jié)構(gòu)分析中大部分浙江茶樹品種與武夷山茶樹聚在一起的結(jié)果一致，也與Zhang 等[45]研究結(jié)果一致。因此，福建北部武夷山與浙江茶區(qū)在地質(zhì)和演化過程可能為同一演化區(qū)，這也暗示著浙江地區(qū)的部分茶樹種與武夷山茶樹種存在一定的基因交流，在近代的茶樹育種中，人為的跨省之間的茶樹雜交育種非常常見，這些品種是否經(jīng)由武夷山茶樹傳播發(fā)展而來還有待進一步研究。

在基因流分析中，福建省內(nèi)存在閩南地區(qū)茶樹向武夷山茶樹的基因流動。一些武夷茶品種/品系，例如紫羅蘭、正太陽、金毛猴、白毛猴、金丁香和留蘭香等在種群結(jié)構(gòu)分析中也與閩南茶樹聚在一起。根據(jù)相關(guān)文獻[46-47]記載，明朝年間，部分閩南茶農(nóng)大舉內(nèi)遷武夷山地區(qū)，部分武夷茶品種/品系可能是由閩南茶流入后經(jīng)過雜交傳播發(fā)展而來。此外，近期培育的小紅袍茶樹因在遺傳相似度分析中顯示與閩南地區(qū)的毛蟹茶樹遺傳相似系數(shù)為0.983 6，具有明顯的派生關(guān)系，未被列入本研究的126個武夷茶樹之中。盡管如此，茶樹的原產(chǎn)地溯源分析反映的是茶樹品種的遺傳背景，以便于后期的茶樹種質(zhì)資源利用與開發(fā)。

3.3 茶樹EDV 及武夷茶親緣關(guān)系

隨著育種技術(shù)的發(fā)展，針對原始品種進行修飾性育種可能出現(xiàn)大量的派生品種。盡管派生品種可提高原始品種的遺傳貢獻，但也使育種的遺傳基礎變窄，不利于作物遺傳改良，甚至可能造成對糧食安全的威脅。2022 年3 月1日，新《中華人民共和國種子法》第28 條和第90 條首次規(guī)定了EDV 制度[10]，對作物品種的遺傳基因相似度提出了要求。茶樹為異花授粉植物，其遺傳組成高度雜合，表現(xiàn)型上千姿百態(tài)，對茶樹親緣關(guān)系的分析有助于茶樹種質(zhì)資源的合理開發(fā)利用，然而茶樹品種的EDV 研究未見報道。本研究首次依據(jù)遺傳相似度對349 個茶樹進行了EDV 初步調(diào)查。在349 個茶樹的遺傳相似度分析中，共發(fā)現(xiàn)136對茶樹遺傳相似系數(shù)大于0.9，其中26 對涉及武夷茶，這說明在茶樹品種資源的大量引種、頻繁的種質(zhì)雜交以及品種選育過程中過度使用優(yōu)良品種，已出現(xiàn)一定量的茶樹EDV，造成“同物異名”的情況[48-49]。

國際種子聯(lián)合會（ International seed federation，ISF）[39]將GS 值大于0.9 作為實質(zhì)派生關(guān)系的基本閾值，基于SNP 標記遺傳相似度的派生關(guān)系閾值的設定建議高于 ISF頒布的標準，且無法對所有作物品種適用統(tǒng)一的EDV 判定標準，必須根據(jù)不同作物品種進行確定[50]。茶樹品種具體的實質(zhì)派生關(guān)系閾值需要進一步深入研究。實質(zhì)派生品種除表型與原始品種存在差異外，基因型/基因型組合與原始品種大致相同?？紤]到即使同一品種在不同的地方栽培時間較長，也會存在極少量基因突變的積累，若兩個表型存在差異的個體的遺傳相似系數(shù)大于兩個不同栽培地的相同品種的遺傳相似系數(shù)，則可以被確認為具有派生關(guān)系。本研究基于遺傳相似度分析結(jié)果發(fā)現(xiàn)，大紅袍2 號株（FJ91）與奇丹（FJ192）相似系數(shù)為0.988 7，大紅袍4 號株（FJ93）與正本不知春（FJ82）相似系數(shù)為0.986 7，均具有派生關(guān)系，屬于同一類型。4 株大紅袍茶樹與北斗（FJ33）、雀舌（FJ48）的遺傳相似系數(shù)均較低（0.7 左右），這與先前的研究結(jié)果一致[51]。醉墨（FJ108）與雀舌（FJ48）的形態(tài)特征與芽葉生長特性十分相似，在葉江華等[17]的ISSR 分子標記鑒定中，相似系數(shù)僅為0.82，被認為是不同類型品種。本研究基于SNP 的分子標記顯示兩者相似系數(shù)為0.989 7，應為同一類型的品種，這可能是因為ISSR 分子標記是多等位基因系統(tǒng)，內(nèi)部存在復雜變異，易產(chǎn)生誤差的緣故[52]。留香澗不知春（FJ42）與正本不知春（FJ82）相似系數(shù)為0.686 1，為不同類型茶樹，這與先前研究一致[17]。此外，本研究還發(fā)現(xiàn)，肉桂（FJ84）、紅海棠（FJ204）、王母桃（FJ73）、金雞母（FJ200）和黃肉桂（FJ191）5 個品種/品系之間均具有明確派生關(guān)系，屬于同一類型。金毛猴（FJ86）與白毛猴（FJ149），向天梅（FJ67）與醉貴姬（FJ71），胭脂柳（FJ76）與醉八仙（FJ77），1113（FJ13）與1114（FJ12）的遺傳相似系數(shù)均大于0.97，也具有明確派生關(guān)系，屬于同一類型。石中玉（FJ52）與0205D（FJ161），石觀音（FJ63）與SR（FJ187），0204（FJ17）與0205（FJ18）的遺傳相似系數(shù)在0.90～0.97，均為近似品種，也可認為具有一定派生關(guān)系。武夷茶種質(zhì)資源歷史悠久、種類繁多，在引種栽培過程中出現(xiàn)錯亂或混雜在所難免。本研究通過SNP 分子標記的遺傳相似度分析詮釋了武夷茶品種/品系之間的關(guān)系，不僅為茶樹育種提供了依據(jù)，而且對遺傳多樣性的保護具有重要作用。

3.4 茶樹的SNP 分子標記鑒定

SNP 分子標記技術(shù)具有基因型信息豐富、檢出效率高且準確性高等優(yōu)勢特點，目前已被廣泛地應用于作物品種的鑒定[53]。盡管近年來茶樹的SNP 分子標記技術(shù)進展快速，并有相關(guān)的報道[54-56]，但茶樹品種資源豐富，基因組龐大且復雜度高，還需要不斷地補充完善。本研究篩選出的973 個SNP 位點均勻地覆蓋在茶樹基因組上，具有較高的基因多樣性，可對349 個茶樹進行有效區(qū)分和指紋圖譜構(gòu)建。數(shù)量少且辨別度高的SNP 標記有利于降低構(gòu)建指紋圖譜成本，本研究從973 個位點中篩選21 個SNP 位點作為一組可快速100%區(qū)分349個茶樹樣本的組合，并設計相關(guān)引物進行驗證，且其中18 個位點即可區(qū)分126 個武夷茶品種/品系。這些SNP 標記在未來茶樹品種分子鑒定中可起到一定的應用，可作為篩選茶樹品種真實性、特異性、純度鑒定以及類群劃分的候選位點組合，為分子標記輔助育種中背景分析提供研究基礎。