房山紫堇葉綠體基因組特征及其系統(tǒng)進(jìn)化關(guān)系

孫 哲， 李澳旋， 杜曉蓉， 宋 蕓， 喬永剛*

(1. 山西農(nóng)業(yè)大學(xué)生命科學(xué)學(xué)院， 山西 太谷 030801； 2. 中獸醫(yī)藥現(xiàn)代化山西省重點(diǎn)實(shí)驗(yàn)室， 山西 太谷 030801)

紫堇屬是罌粟科最大的一個(gè)屬，在全球約有440種，為草本植物，無乳汁[1]。房山紫堇(CorydalisfangshanensisW. T. Wang ex S. Y. He)又稱石黃連、土黃連，是罌粟科紫堇屬的一種多生長于海拔500～1 600 m左右石灰?guī)r山坡的植物[2]，開白色小花，生命力頑強(qiáng)，可防風(fēng)固沙、保持水土[3]；其根可入藥，具有清熱解毒的功效[4]。因獨(dú)特的形態(tài)，特殊的生長環(huán)境和稀少的數(shù)量，房山紫堇被列為北京市重點(diǎn)保護(hù)野生植物。開展房山紫堇的保護(hù)和科學(xué)研究，可以使其在園林綠化、水土保持等方面發(fā)揮重要作用，也有利于豐富草地生態(tài)系統(tǒng)。

葉綠體是植物光合作用的場所[5]，隨著一些生物技術(shù)的開發(fā)和應(yīng)用，葉綠體的基因組大小得到確定，被子植物的葉綠體基因組大小為120～170 kb不等，通常編碼110～130個(gè)基因[6]。高等植物的葉綠體基因組結(jié)構(gòu)高度保守，可以確保在進(jìn)化群體之間基因的直接同源性，包括反向重復(fù)區(qū)(Inverted repeats,IRs)和大單拷貝區(qū)(Large single copy,LSC)、小單拷貝區(qū)(Small single copy,SSC)，植物葉綠體基因組的大小差異由IRa和IRb決定[7]。葉綠體基因組所包含的遺傳信息沒有核基因組多，但其結(jié)構(gòu)穩(wěn)定并且保守[8]，可以為高等植物群體遺傳學(xué)、系統(tǒng)發(fā)育和基因組進(jìn)化提供依據(jù)[9-10]。

紫堇屬植物在形態(tài)上很難區(qū)分，被認(rèn)為是罌粟科中最難分類的一個(gè)屬[11]，葉綠體基因組可以為紫堇屬植物的系統(tǒng)發(fā)育提供參考，國內(nèi)外已經(jīng)有紫堇(CorydalisedulisMaxim.)、陜西紫堇(CorydalisshensianaLiden)[12]、毛黃堇(CorydalistomentellaFranch.)、石生黃堇(CorydalissaxicolaBunting)[13]、斑花黃堇(CorydalisconspersaMaxim.)[14]等紫堇屬植物葉綠體基因組的相關(guān)報(bào)道。本研究對(duì)房山紫堇葉綠體基因組進(jìn)行測序，對(duì)其GC含量、SSR位點(diǎn)等進(jìn)行分析，并與紫堇屬其他物種的基因組進(jìn)行多重比較，以期為房山紫堇的系統(tǒng)發(fā)育和紫堇屬植物的親緣關(guān)系提供數(shù)據(jù)支撐，同時(shí)豐富紫堇屬植物的葉綠體基因組數(shù)據(jù)庫。

1 材料與方法

1.1 試驗(yàn)材料

房山紫堇種子采集自山西省陵川縣丈河村(北緯35°33′～35°38′，東經(jīng)113°08′～113°17′)，在實(shí)驗(yàn)室育苗后，移栽到藥用植物資源圃。于2020年10月采集健康無病害的房山紫堇植株的新鮮葉片，用酒精棉擦拭后，快速放入液氮，保存于—80℃冰箱備用。

1.2 試驗(yàn)方法

1.2.1DNA提取和測序 用CTAB法提取保存的房山紫堇葉片DNA[15]，DNA質(zhì)量和濃度檢測合格后，將DNA片段化(超聲波打斷)，經(jīng)過片段純化、末端修復(fù)等途徑獲得測序文庫，質(zhì)檢合格后在陜西博瑞德生物科技有限公司用Illumina HiSeq進(jìn)行測序。

1.2.2基因的組裝和注釋 用NOVOPlasty軟件[16]對(duì)測序所得的房山紫堇葉綠體基因組進(jìn)行組裝，參考序列為Coreanomeconhylomeconoides(NC_031446.1)[17]，獲得的葉綠體基因組用CPGAVAS軟件[18]進(jìn)行基因注釋，結(jié)果提交至GenBank，序列號(hào)為MZ772033，用OGDRAW在線軟件[19]將結(jié)果可視化，GC含量用MEGA7[20]軟件進(jìn)行分析。

1.2.3簡單重復(fù)序列(Simple sequence repeats，SSR)位點(diǎn)分析 利用MISA軟件[21]對(duì)房山紫堇的SSR位點(diǎn)進(jìn)行分析，設(shè)置重復(fù)次數(shù)：單核苷酸>10，雙核苷酸>6，三核苷酸>5，四核苷酸>5，五核苷酸>5，六核苷酸>5；2個(gè)SSR之間的最小距離為100 bp。

1.2.4密碼子組成及偏好性 運(yùn)用軟件CodonW1.4.2[22]分析其葉綠體基因組的有效密碼子數(shù)(Effective number of codon,Nc)、相對(duì)同義密碼子使用情況(Relative synonymous codon usage,RSCU)、密碼子適應(yīng)指數(shù)(Codon adaption index,CAI)、密碼子偏好性指數(shù)(Codon bias index,CBI)、最優(yōu)密碼子使用頻率(Frequency of optimal codons,FOP)等密碼子組成參數(shù)，RSCU>1的密碼子具有偏好性[22]。

1.2.5多重比較和IR收縮與擴(kuò)張分析 利用mVISTA軟件[23]將房山紫堇與石生黃堇、紫堇、毛黃堇葉綠體基因組進(jìn)行多重比較，使用IRscope軟件[24]分析房山紫堇與石生黃堇、紫堇、毛黃堇的IR收縮與擴(kuò)張情況。

1.2.6系統(tǒng)發(fā)育分析 以房山紫堇基因組序列和已經(jīng)公布的9個(gè)紫堇屬植物基因組序列為研究對(duì)象，選取苦參(SophoraflavescensAlt.)和苦豆子(SophoraalopecuroidesL.)為外群[25]，利用MAFFT軟件[26]對(duì)12個(gè)物種的葉綠體全基因組進(jìn)行比對(duì)，并使用MEGA7軟件構(gòu)建ML系統(tǒng)發(fā)育樹，自展值Bootstrap value設(shè)為1 000，使用Kimura 2-parameter模型。

2 結(jié)果與分析

2.1 葉綠體基因組特征

基于邊合成邊測序技術(shù)，使用Illumina HiSeq高通量測序平臺(tái)對(duì)房山紫堇葉綠體基因組進(jìn)行測序，共得到10.45 GB原始數(shù)據(jù)，各樣品Q30堿基百分比均不小于90%。房山紫堇的葉綠體基因組為191 128 bp，GC含量為40.2%;呈現(xiàn)典型的四分體結(jié)構(gòu)(圖1)，大單拷貝區(qū)(LSC)為97 302 bp、小單拷貝區(qū)(SSC)為9 836 bp，反向重復(fù)區(qū)(IRa，IRb)為83 990 bp，GC含量分別為39.1%，35.2%，42.1%。

圖1 房山紫堇葉綠體基因組示意圖Fig.1 Chloroplast genome map of Corydalis fangshanensis

房山紫堇葉綠體基因組共注釋到127個(gè)基因(表1)，包括84個(gè)蛋白質(zhì)編碼基因、35個(gè)tRNA基因以及8個(gè)rRNA基因，占比分別為66.1%，27.6%，6.3%。根據(jù)不同的功能可將這些基因分為4類，包括atpA，atpB，atpE，atpF，atpH等44個(gè)與光合作用有關(guān)的基因；rpl14，rpl16等73個(gè)與自我復(fù)制相關(guān)的基因；還有infA等其他基因以及一些功能未知的基因。其中17個(gè)基因具有內(nèi)含子，rps12，trnA-UGC，trnE-UUC，ndhB基因的內(nèi)含子多于1。

表1 房山紫堇葉綠體基因功能分類Table 1 Chloroplast gene functional classification of Corydalis fangshanensis

2.2 SSR位點(diǎn)分析

對(duì)房山紫堇單基因序列進(jìn)行分析，結(jié)果顯示(表2)，共檢測到40個(gè)SSR位點(diǎn)，其中單核苷酸有31個(gè)，雙核苷酸有8個(gè)，三核苷酸有1個(gè)。單核苷酸最多，包括7個(gè)A(占比17.5%)，1個(gè)C(占比2.5%)，23個(gè)T(占比57.5%)；雙核苷酸總共8個(gè)，包括4個(gè)AT(占比10%)和4個(gè)TA(占比10%)；三核苷酸只有1個(gè)CAA(占比2.5%)。

表2 房山紫堇SSR位點(diǎn)分析Table 2 SSR loci analysis of Corydalis fangshanensis

2.3 密碼子偏好性分析

對(duì)房山紫堇葉綠體基因組的密碼子組成分析結(jié)果顯示(表3)：房山紫堇基因組的有效密碼子數(shù)為56.57，說明其密碼子偏好性較弱；GC含量為41.2%，GC3s含量為38.5%，均小于50%，說明密碼子偏向使用堿基A和堿基T。

表3 房山紫堇葉綠體基因組密碼子使用的總體特征Table 3 Overall characterization of codon usage in the chloroplast genome of Corydalis fangshanensis

對(duì)房山紫堇葉綠體基因組的相對(duì)同義密碼子進(jìn)行分析(圖2)，RSCU>1的密碼子有35個(gè)，其中以A/T堿基結(jié)尾的有28個(gè)，占80%，表明這些密碼子偏向以A/T結(jié)尾。

圖2 房山紫堇葉綠體基因組的相對(duì)同義密碼子使用情況Fig.2 Relative synonymous codon usage in chloroplast genomes of Corydalis fangshanensis

2.4 比較基因組學(xué)和IR收縮與擴(kuò)張分析

以測序所得的房山紫堇基因組序列為參照，將其與石生黃堇、紫堇、毛黃堇進(jìn)行相似性比較分析，結(jié)果如圖3所示，4種紫堇屬植物的編碼區(qū)較為保守，在ycf2，trnN-GUU，ndhF等基因所在的位置存在明顯的變異。邊界分析結(jié)果如圖4所示，房山紫堇rpl2基因位于LSC區(qū)；ndhA基因位于IRb與SSC的交界處，在IRb區(qū)267 bp，SSC區(qū)域1 902 bp；ycf1位于SSC區(qū)域，距離SSC與IRa的交界407 bp；rpl23位于IR區(qū)域，距離IRa與LSC的交界35 bp，石生黃堇與房山紫堇的邊界收縮情況最為相似。

圖3 房山紫堇與石生黃堇、紫堇、毛黃堇的多重比較結(jié)果Fig.3 The results of multiple comparisons between Corydalis fangshanensis and Corydalis saxicola,Corydalis edulis and Corydalis tomentella

圖4 房山紫堇與石生黃堇、紫堇、毛黃堇的邊界比較Fig.4 Boundary comparison of Corydalis fangshanensis and Corydalis saxicola,Corydalis edulis and Corydalis tomentella

2.5 系統(tǒng)發(fā)育分析

為了鑒定房山紫堇在紫堇屬植物中的進(jìn)化位置，使用MEGA7構(gòu)建ML系統(tǒng)進(jìn)化樹，結(jié)果顯示(圖5)，紫堇屬植物可以分為兩大類，房山紫堇與同屬植物紫堇、毛黃堇、石生黃堇聚為一大類，與紫堇[10]的親緣關(guān)系最近；灰綠黃堇(CorydalisaduncaMaxim.)[27]與其余5種紫堇屬植物聚為另一大類；外源植物苦參與苦豆子聚為一類，與前人的研究一致[16]。

圖5 基于葉綠體基因組的紫堇屬植物系統(tǒng)發(fā)育分析Fig.5 Phylogenetic analysis of Corydalis based on chloroplast genome

3 討論

紫堇屬植物物種豐富，種類繁多，僅根據(jù)形態(tài)很難對(duì)其分類。2015年山西首次報(bào)道發(fā)現(xiàn)了房山紫堇，房山紫堇全株可入藥，抗逆性強(qiáng)，對(duì)環(huán)境要求較低，可在干旱缺水的石灰?guī)r環(huán)境中生存，具有防風(fēng)固沙、水土保持的生態(tài)價(jià)值[28]。房山紫堇被列入《世界自然保護(hù)聯(lián)盟紅色名錄》(IUCN)中，中國珍稀瀕危植物信息系統(tǒng)將其列為易危(VU)植物，禁止直接采挖[29]。

目前有部分紫堇屬植物的葉綠體基因組數(shù)據(jù)得到公布[30]，大部分紫堇屬植物的葉綠體基因組有待揭示。葉綠體基因、線粒體基因和細(xì)胞核基因是系統(tǒng)發(fā)育和物種起源研究中常用的三種數(shù)據(jù)來源[31]。高度保守的葉綠體基因組[8]含有豐富的遺傳信息，研究房山紫堇的葉綠體基因組對(duì)了解房山紫堇的屬內(nèi)進(jìn)化位置、遺傳變異、系統(tǒng)發(fā)育以及開發(fā)利用等具有重要的價(jià)值。

SSR位點(diǎn)分析對(duì)研究物種鑒定、遺傳變異、系統(tǒng)進(jìn)化有重要意義，對(duì)房山紫堇進(jìn)行SSR位點(diǎn)分析，發(fā)現(xiàn)其有40個(gè)SSR位點(diǎn)，包括單核苷酸、雙核苷酸和三核苷酸，房山紫堇葉綠體基因組SSR多含有AT堿基，與前人研究結(jié)果一致[32]。Nc，RSCU，GC，GC3s均可以反映密碼子的偏好性，Nc值通常在20～61之間，靠近61偏好性弱，靠近20偏好性強(qiáng)；RSCU>1時(shí)，該密碼子為使用頻率較高的密碼子，RSCU=1時(shí)，密碼子沒有偏好性，RSCU<1時(shí)，其為低使用頻率的密碼子；GC和GC3s含量大于50%時(shí)，說明該物種葉綠體基因組偏向使用GC堿基，反之，偏向使用AT堿基。對(duì)房山紫堇的密碼子分析結(jié)果顯示：Nc值靠近61，說明房山紫堇葉綠體基因組密碼子使用性偏弱；RSCU>1的密碼子超過半數(shù)，其中以A/T堿基結(jié)尾的占80%，并且GC含量和GC3s含量均小于50%，房山紫堇偏向使用AT堿基，與段春燕[33]、宋蕓[34]等人的研究結(jié)果一致。植物葉綠體基因組的收縮和擴(kuò)張是其變異的主要機(jī)制，本研究對(duì)房山紫堇基因組與石生黃堇、紫堇、毛黃堇進(jìn)行比較分析，結(jié)果顯示石生黃堇與房山紫堇的邊界收縮情況高度相似。聚類結(jié)果顯示：紫堇屬植物序列被聚為一類，兩種豆科植物作為外源種被聚為一類[25]，與預(yù)期結(jié)果一致，房山紫堇與紫堇的親緣關(guān)系最近，在對(duì)房山紫堇進(jìn)行開發(fā)利用時(shí)，可將紫堇作為參考物種。本研究豐富了紫堇屬植物的葉綠體基因組數(shù)據(jù)，為探索紫堇屬植物的起源及親緣關(guān)系提供了依據(jù)，為房山紫堇的進(jìn)一步科學(xué)研究奠定了基礎(chǔ)。

4 結(jié)論

本研究用Illumina HiSeq對(duì)房山紫堇進(jìn)行測序，將測序結(jié)果進(jìn)行整理得出以下結(jié)論：房山紫堇的葉綠體基因組大小為191 128 bp，GC含量為40.2%，呈現(xiàn)典型的四段式結(jié)構(gòu)，注釋到127個(gè)基因；有40個(gè)SSR位點(diǎn)；葉綠體基因組的密碼子偏好性較弱，偏向使用A和T這兩種堿基；編碼區(qū)較為保守；房山紫堇與紫堇的親緣關(guān)系最近。目前對(duì)房山紫堇的報(bào)道很少，本研究為房山紫堇的物種保護(hù)、系統(tǒng)進(jìn)化、資源的開發(fā)利用提供了相關(guān)試驗(yàn)數(shù)據(jù)和理論依據(jù)。