摘要 以6個鵝掌楸（Liriodendron chinense）種源和3個北美鵝掌楸（Liriodendron tulipifera）種源共計23個個體作為研究群體，選取11對真實SNP位點，對這253個SNP標(biāo)記的PCR產(chǎn)物進(jìn)行正反向Sanger測序。結(jié)果表明：在11個真實SNP位點中，8個SNP位點在中國鵝掌楸和北美鵝掌楸間呈現(xiàn)出多態(tài)性，其中3個位點是檢測鵝掌楸和北美鵝掌楸種間分化的潛在標(biāo)記。11個SNP標(biāo)記在北美鵝掌楸種內(nèi)都沒有多態(tài)性。lm_ll_10205、lm_ll_10836、lm_ll_3853和lm_ll_8246 4個SNP位點在我國東部的湖南瀏陽、浙江松陽和江西廬山3個種源間存在多態(tài)性，僅有l(wèi)m_ll_10205位點在我國西部種源間存在多態(tài)性。前期開發(fā)的SNP標(biāo)記可以用于鵝掌楸種源遺傳多樣性的研究。研究種源間SNP位點的變異特征為揭示鵝掌楸的起源、進(jìn)化以及遺傳變異多樣性提供了科學(xué)數(shù)據(jù)。

關(guān)鍵詞 北美鵝掌楸；鵝掌楸；SNP標(biāo)記；SNP位點；種源

中圖分類號 S792.21 文獻(xiàn)標(biāo)識碼 A 文章編號 0517-6611（2024）21-0093-04

doi：10.3969/j.issn.0517-6611.2024.21.019

開放科學(xué)（資源服務(wù)）標(biāo)識碼（OSID）：

Application of SNP Markers of Liriodendron Genome in Provenance Evaluation

LU Ye， SHI Ji-sen

（Key Laboratory of Forest Genetics & Biotechnology of Ministry of Education of China， Co-Innovation Center for Sustainable Forestry in Southern China， Nanjing Forestry University， Nanjing， Jiangsu 210037）

Abstract In this study， a total of 23 individuals were selected from 9 provenances with representative geographical distribution （including six Liriodendron chinense and three Liriodendron tulipifera provenances） as the study population. After PCR with 11 pairs of primers from real SNPs， Sanger sequencing was performed on the 253 SNP-labeled PCR products. The results showed that among the 11 real SNP loci， 8 SNP loci showed polymorphism between Chinese and North American provenances， and 3 of them were potential markers for detecting the interspecific differentiation of Chinese and North American provenances. None of the 11 SNP markers were polymorphic within the North American Liriodendron provenances. Four SNP loci （lm_ll_10205， lm_ll_10836， lm_ll_3853 and lm_ll_8246） showed polymorphism in Liuyang （Hunan）， Songyang （Zhejiang） and Lushan （Jiangxi） species in eastern China， while only lm_ll_10205 showed polymorphism in western species of China. The results showed that the SNP markers developed in earlier studies could be used to study genetic diversity of Liriodendron provenance. The variation characteristics of SNP sites between provenances provide useful scientific data for revealing the origin， evolution and diversity of genetic variation of Liriodendron.

Key words Liriodendron tulipifera;Liriodendron chinense;SNP makers;SNPs loci;Provenance

基金項目 江蘇省農(nóng)業(yè)科技自主創(chuàng)新資金項目“耐鹽型雜交鵝掌楸新品種培育及推廣示范”（SCX（23）3628）。

作者簡介 陸葉（1982— ），女，江蘇鎮(zhèn)江人，實驗師，從事林木遺傳育種研究。*通信作者，教授，博士生導(dǎo)師，從事林木遺傳育種研究。

收稿日期 2023-12-29

木蘭科（Magnoliaceae）鵝掌楸屬（Liriodendron）現(xiàn)僅存2個種，即鵝掌楸［Liriodendron chinense（Hemsl.）Sarg.］和北美鵝掌楸（Liriodendron tulipifera Linn.），是典型的孑遺植物［1］。新生代第三紀(jì)，該屬物種廣泛分布于歐亞大陸和北美洲。第四紀(jì)冰川之后，鵝掌楸分布于我國長江流域和越南北部，分化為東部和西部兩個支系，北美鵝掌楸自然分布于北美洲東南部，遺傳上位于我國東西部鵝掌楸類群之間，更接近東部類群［2］。由于結(jié)實率低，加上人為對其適生環(huán)境的干擾，鵝掌楸在我國已被列為二類瀕危保護(hù)物種［3］。鵝掌楸屬物種不僅在林業(yè)產(chǎn)業(yè)方面具有重要的應(yīng)用前景，而且也是林木遺傳多樣性、群體結(jié)構(gòu)、系統(tǒng)進(jìn)化及雜交育種等領(lǐng)域理論研究的理想材料［4-5］。

RAD-seq（Restriction-site associated DNA sequencing）技術(shù)是利用限制性內(nèi)切酶對基因組DNA進(jìn)行酶切，然后進(jìn)行高通量測序［6］，可作為一種有效的方法，進(jìn)行高通量的SNP標(biāo)記的開發(fā)［7］。目前已有很多研究基于RAD-seq技術(shù)對物種的遺傳多樣性和遺傳結(jié)構(gòu)進(jìn)行分析。劉學(xué)鋒等［8］對168份藍(lán)花楹種質(zhì)材料進(jìn)行RAD-seq測序，獲得45 552個高質(zhì)量的SNP，分析表明，來自川渝地區(qū)的藍(lán)楹花具有相對較近的親緣關(guān)系。潘文婷等［9］用RAD-seq鑒定了北美和中國種源的鵝掌楸共4 454個SNP標(biāo)記，對鵝掌楸的遺傳多樣性和遺傳結(jié)構(gòu)進(jìn)行分析，認(rèn)為鵝掌楸屬遺傳結(jié)構(gòu)的形成與其地理隔離和片段化分布有關(guān)。

前期研究以北美鵝掌楸和鵝掌楸的雜交F1代和2個親本為材料，利用RAD-seq技術(shù)已開發(fā)的大量可靠的SNP標(biāo)記［10］，筆者從中選取了11對真實SNP位點，以6個鵝掌楸種源和3個北美鵝掌楸種源共計23個個體為研究材料，分析SNP位點在不同種源間的變異特征，旨在為研究鵝掌楸物種分子水平上的變異提供科學(xué)依據(jù)。

1 材料與方法

1.1 試驗材料

種源材料均來自國家林業(yè)和草原局世行貸款項目管理中心于1990年推廣支持體系中的“闊葉樹中心組”的科研項目建立的實驗林。在鵝掌楸全分布區(qū)內(nèi)抽樣收集了15個地理種源的鵝掌楸種子，同時分別從美國密蘇里（Missouri）、路易斯安那（Louisiana）、北卡羅來納（North Carolina，北卡）、南卡羅來納（South Carolina，南卡）、佐治亞（Georgia）引進(jìn)該地區(qū)北美鵝掌楸種源的種子，翌年分別在四川邛崍、湖北京山、湖南桃源、江西分宜和福建邵武進(jìn)行了地理種源研究試驗。2014年4月從定植于湖北省京山縣虎爪山林場的24年生的種源材料中選取了9個種源的材料，每個種源23株，共計23株作為后續(xù)試驗材料。取樣具體如下：湖南瀏陽的3株鵝掌楸（1-10、1-13、1-23），浙江松陽的2株鵝掌楸（4-1、4-24），江西廬山的3株鵝掌楸（8-3、8-17、8-24），四川敘永的3株鵝掌楸（9-5、9-19、9-22），貴州松桃的3株鵝掌楸（11-12、11-15、11-22），云南勐臘的2株鵝掌楸（12-5、12-22），北卡羅來納的3株北美鵝掌楸（15-19、15-21、15-23），南卡羅來納的2株北美鵝掌楸（18-16、18-21），密蘇里的2株北美鵝掌楸（19-16、19-21）。其中，來自湖南瀏陽、浙江松陽和江西廬山的種源為我國東南部種源，來自四川敘永、貴州松桃和云南勐臘的種源為我國西南部種源，來自北卡羅來納、南卡羅來納和密蘇里的種源為北美種源。

1.2 試驗方法

1.2.1 基因組DNA的提取。

新鮮葉片用液氮冷凍，-80 ℃保存。采用天根試劑盒（DNAsecure Plant Kit，離心柱型）提取DNA，保存在-20 ℃冰箱中。

1.2.2 引物合成。

根據(jù)前期研究開發(fā)的鵝掌楸27個真實SNP位點［10］，選取其中的11個SNP位點，利用鵝掌楸和北美鵝掌楸不同種源個體，對SNP位點在種群間的多態(tài)性進(jìn)行研究，引物序列見表1。

1.2.3 PCR擴增體系。

用11對引物和KOD FX（Toyobo，KFX 101）DNA聚合酶對23個DNA樣本進(jìn)行PCR，將具有特異性條帶的PCR產(chǎn)物進(jìn)行Sanger測序。測序由上海英駿生物技術(shù)有限公司完成。用BioXM2.6軟件比對測序結(jié)果。

PCR反應(yīng)體系（50 μL）為：2×PCR Buffer for KOD FX 25 μL，2 mmol/L dNTPs 10 μL，10 μmol/L Primer 1.5 μL，KOD FX DNA polymerase 1 μL，DNA 5 μL，ddH2O 6 μL。PCR 反應(yīng)程序如下：94 ℃預(yù)變性 2 min；98 ℃變性10 s，Tm 5 ℃退火 30 s，68 ℃延伸 1 min，36 個循環(huán)；最后4 ℃保溫。

2 結(jié)果與分析

2.1 DNA質(zhì)量檢測

對提取的DNA進(jìn)行瓊脂糖電泳檢測，如圖1所示，DNA條帶清晰完整，無明顯彌散條帶，表明DNA質(zhì)量好且純度較高。用Nanodrop檢測DNA樣品的OD260/OD280值，均在1.8～2.0。高質(zhì)量的DNA達(dá)到后續(xù)試驗要求。

2.2 PCR產(chǎn)物檢測

PCR產(chǎn)物經(jīng)過瓊脂糖凝膠電泳檢測，將符合預(yù)期的目的條帶切膠回收、純化后進(jìn)行Sanger測序驗證。圖2為部分引物在不同個體間的PCR擴增結(jié)果。

2.3 SNP位點檢測

對總計253個包含SNP標(biāo)記的PCR產(chǎn)物進(jìn)行正、反向測序，因49個樣品未能測序成功，最后獲得204個測序結(jié)果，結(jié)果見表2。

在11個真實的SNP位點中，nn_np_7542、lm_ll_11011、nn_np_5667、nn_np_4694、lm_ll_10205、lm_ll_10836、lm_ll_3853和lm_ll_8246 8個SNP位點在鵝掌楸和北美鵝掌楸間呈現(xiàn)出多態(tài)性，其中l(wèi)m_ll_7542、lm_ll_5667、lm_ll_4694 3個位點雖在鵝掌楸種內(nèi)和北美鵝掌楸種內(nèi)未檢測到多態(tài)性，但是在鵝掌楸和北美鵝掌楸種間呈現(xiàn)出多態(tài)性，是檢測鵝掌楸和北

美鵝掌楸種間分化的潛在標(biāo)記。nn_np_9274、nn_np_9109和nn_np_9757 3個SNP位點無論是在鵝掌楸種內(nèi)，還是北美鵝掌楸種內(nèi)，甚至鵝掌楸和北美鵝掌楸種間都未檢測出多態(tài)性。 11個SNP位點，在北美鵝掌楸種內(nèi)皆為單一且相同堿基，因此11個SNP位點在北美鵝掌楸種內(nèi)無多態(tài)性。

從我國種源內(nèi)部來看，除了lm_ll_11011 SNP位點缺少測序結(jié)果，lm_ll_10205、lm_ll_10836、lm_ll_3853和lm_ll_8246 4個SNP位點在我國東部的湖南瀏陽、浙江松陽和江西廬山3個種源間存在多態(tài)性，其中l(wèi)m_ll_10205位點在浙江松陽、江西廬山、四川敘永、貴州松桃、云南勐臘間呈現(xiàn)多態(tài)性，并且其在四川敘永種內(nèi)也存在多樣性。lm_ll_10836、lm_ll_3853和lm_ll_8246位點在我國西部的四川敘永、貴州松桃和云南勐臘3個種源間沒有多態(tài)性，其基因型分別為T、G、A，但都與東部的江西廬山種源呈現(xiàn)多態(tài)性。lm_ll_10836位點在浙江松陽的基因型也是T，與西部種源沒有多態(tài)性，但在湖南瀏陽種源中出現(xiàn)了基因型C。lm_ll_3853和lm_ll_8246位點，在東部種源湖南瀏陽的基因型和西部種源相同，沒有多態(tài)性。

3 討論

SNP即單核苷酸多態(tài)性（single nucleotide polymorphism）是繼RFLP、SSR標(biāo)記后的第3代分子標(biāo)記，普遍存在于生物基因組中［11］。由于SNP遺傳穩(wěn)定性高，位點豐富且分布廣泛，富有代表性、二態(tài)性和等位基因性，檢測快速，易于實現(xiàn)自動化分析等優(yōu)點，因此，SNP標(biāo)記已經(jīng)在遺傳圖譜、基因遺傳、分子進(jìn)化、功能基因組研究等領(lǐng)域得到了廣泛運用［12］。

為了研究SNP位點在不同種源之間變異情況，該研究從前期開發(fā)的SNP位點中選取11個真實的SNP位點，采集6個我國種源和3個北美種源，共23個個體作為研究群體，結(jié)果表明，SNP位點在鵝掌楸不同種源間呈現(xiàn)不同的變異。11個SNP標(biāo)記在北美鵝掌楸種內(nèi)都沒有多態(tài)性。lm_ll_7542、lm_ll_5667、lm_ll_4694在鵝掌楸和北美鵝掌楸中間呈現(xiàn)出多態(tài)性，可以作為檢測鵝掌楸和北美鵝掌楸種間分化的潛在標(biāo)記。11個SNP位點在6個鵝掌楸種源和3個北美鵝掌楸種源間的多態(tài)性高于6個鵝掌楸種源內(nèi)的多態(tài)性，這可能與鵝掌楸物種的進(jìn)化關(guān)系有關(guān)。鵝掌楸在經(jīng)歷了第四紀(jì)冰川時期后，僅存在于我國漢水流域以南地區(qū)，以及北美洲的美國東部沿海至加拿大的安大略省五大湖區(qū)的多倫多地區(qū)。北美鵝掌楸種源在進(jìn)化歷程中SNP位點未見較大的突變發(fā)生；而在鵝掌楸種源中，我國東南部種源相比于西南部種源發(fā)生了一定比例的堿基突變，在SNP位點處堿基種類較多。由此可見，前期開發(fā)出來的SNP標(biāo)記可以用于鵝掌楸的起源和遺傳多樣性研究。在后續(xù)的研究中可增加種源的個體數(shù)和SNP標(biāo)記數(shù)量，確保測序準(zhǔn)確度的同時，為鵝掌楸種源變異的研究提供科學(xué)依據(jù)。

參考文獻(xiàn)

［1］ 王章榮.鵝掌楸屬雜交育種成就與育種策略［J］.林業(yè)科技開發(fā)，2008（5）：1-4.

［2］ CHEN J H，HAO Z D，GUANG X M，et al.Liriodendron genome sheds light on angiosperm phylogeny and species-pair differentiation［J］.Nature plants，2019，5（1）：18-25.

［3］ 季孔庶，王章榮.鵝掌楸屬植物研究進(jìn)展及其繁育策略［J］.世界林業(yè)研究，2001，14（1）：8-14.

［4］ 羅西，郭秋菊，姚蘭，等.瀕危植物鵝掌楸的天然種群結(jié)構(gòu)特征［J］.中南林業(yè)科技大學(xué)學(xué)報，2021，41（7）：115-123.

［5］ 祝園青，白天道，王曉波，等.馬褂木母樹林與不同來源無性系的遺傳多樣性比較分析［J］.江西農(nóng)業(yè)大學(xué)學(xué)報，2021，43（3）：630-639.

［6］ MILLER M R，DUNHAM J P，AMORES A，et al.Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA（RAD）markers［J］.Genome research，2007，17（2）：240-248.

［7］ RODRGUEZ-EZPELETA N，BRADBURY I R，MENDIBIL I，et al.Population structure of Atlantic mackerel inferred from RAD-seq-derived SNP markers：Effects of sequence clustering parameters and hierarchical SNP selection［J］.Molecular ecology resources，2016，16（4）：991-1001.

［8］ 劉學(xué)鋒，李小梅，張國武，等.基于RAD高通量測序的藍(lán)花楹群體遺傳分析［J］.熱帶亞熱帶植物學(xué)報，2022，30（5）：613-622.

［9］ 潘文婷，孫建軍，原勤勤，等.RAD-seq技術(shù)研究鵝掌楸屬種源遺傳多樣性和遺傳結(jié)構(gòu)［J］.林業(yè)科學(xué)，2022，58（4）：74-81.

［10］ 陸葉，龍曉飛，王鵬凱，等.基于RAD-seq技術(shù)的鵝掌楸基因組SNP標(biāo)記開發(fā)［J］.南京林業(yè)大學(xué)學(xué)報（自然科學(xué)版），2019，43（4）：1-7.

［11］ DURSTEWITZ G，POLLEY A，PLIESKE J，et al.SNP discovery by amplicon sequencing and multiplex SNP genotyping in the allopolyploid species Brassica napus［J］.Genome，2010，53（11）：948-956.

［12］ 王富強，樊秀彩，張穎，等.SNP分子標(biāo)記在作物品種鑒定中的應(yīng)用和展望［J］.植物遺傳資源學(xué)報，2020，21（5）：1308-1320.