蕨類植物葉綠體基因ycf94的分子進化研究

李子菲, 蘇應娟, 王艇*

蕨類植物葉綠體基因的分子進化研究

李子菲1, 蘇應娟2,3*, 王艇1*

(1. 華南農業大學生命科學學院, 廣州 510642; 2. 中山大學生命科學學院, 廣州 510275; 3. 中山大學深圳研究院, 深圳 518057)

基因是近年來在葉綠體基因組中新發現的一個基因，在蕨類植物中表現高度保守。該研究共選取94種蕨類植物，在系統發育背景下，對基因的結構特征、密碼子偏好性、進化速率和適應性進化進行分析。結果表明,基因的密碼子偏好性較弱，偏好使用以 A/U 結尾的密碼子，且不同物種間的偏好性存在一定差異。密碼子偏好性的形成主要受到突變壓的影響，同時也存在其他因素的作用；基于鳳尾蕨科和其他蕨類中基因的結構特征存在區別，對兩者的分子替換速率進行了比較，表明顛換率、非同義替換率和值間存在顯著差異；僅檢測出1個正選擇位點74A，強烈的負選擇作用表明基因的結構和功能基本趨于穩定。這為蕨類系統發育分析提供了新依據，并提供了解析基因功能的線索。

蕨類植物；；密碼子偏好性；進化速率；適應性進化

葉綠體基因組具有長度短、拷貝數高、保守性強的特點，長期被廣泛用于植物系統發育和分子進化研究[1]。葉綠體基因組作為雙鏈環狀閉合DNA,大小一般為110～116 kb，具有典型的四區結構，包含了大單拷貝區(large single copy, LSC)、小單拷貝區(small single copy, SSC)和2個反向重復區(inverted repeats, IR)。目前普遍認為陸地植物葉綠體基因組的基因組成是大致相似的，大約有80個蛋白編碼基因(CDS, protein-coding gene)、30個tRNA基因和4個rRNA基因[2]。其中也包含一些保守的開放閱讀框(open reading frame, ORF)，其所編碼的蛋白相關功能仍在摸索中，通常稱為假定葉綠體開放閱讀框(hypothetical chloroplast open reading frame)基因()。

蕨類作為最早脫離水體適應陸生環境、進化出維管組織的植物類群，在整個陸地植物的起源和演化歷史研究中有著重要的地位[3]。蕨類植物現存約有12 000種，它在形態上具有豐富的多樣性，廣泛分布在熱帶和亞熱帶地區。同時，蕨類的質體基因組有獨特的結構，相比于與種子植物的質體基因組結構基本一致的蕨類祖先類型“蓮座蕨型”，“鐵線蕨型”分別在跨IR/LSC區和LSC區中發生了約20和3 kb的兩次重疊倒位，而“鐵線蕨型”占據了蕨類90%的種類[4]。2016年發表的PPG I (the pteridophyte phylogeny group Ⅰ)系統是一個基本建立在分子系統發育研究之上的現代蕨類植物分類系統[5]，其中描述的分類較為全面且被接受。目前，有關蕨類葉綠體基因組和分類系統等各類信息仍在不斷完善中。

Song等[6]于2018年在薄囊蕨類中發現1個新的假定ORF，并將其命名為。研究表明，位于K和16之間，在薄囊蕨類中表現高度保守。同時在蕨類、石松類和苔蘚類的質體中都發現其同源性，而種子植物中不存在該序列。在陸地植物中質體的進化歷史是動態的，其基因組存在基因丟失的現象，基因會從質體基因組轉移到核基因組[7]。或許存在于陸生植物的祖先質體中，隨后在種子植物中發生基因丟失。

根據序列的保守程度，以及其編碼跨膜蛋白和具有多個RNA編輯位點的特征[6]，可以表明是一個具有功能意義的蛋白編碼基因[8]，但其起源和功能目前尚不明確。本研究以94種蕨類植物的基因為研究對象，通過分析的堿基組成和密碼子偏好性，并在系統發育背景下，分析其進化速率和選擇壓力，來探究的基因結構特征和在蕨類植物中的進化樣式，挖掘在葉綠體基因組系統發育關系分析中的可用性，為蕨類植物分類提供新的依據，為進一步研究的起源和功能等提供基礎信息。

1 材料和方法

1.1 提取序列數據

根據的基因序列作為參考序列，用NCBI-BLAST在線比對(E value<1e-6)在GenBank數據庫上選取了94種蕨類植物葉綠體全基因組，涵蓋了3目19科60屬(附錄1)。采用Genious prime軟件[9]提取基因，并使用MEGA X軟件通過MUSCLE (密碼子)模塊進行序列比對并對結果進行手工校正。

1.2 構建系統發育樹

進行94種蕨類植物基于K+L串聯數據集的多基因系統發育分析，構建3種系統發育樹。PhyloSuite v1.2.2軟件[10]導入94條序列并提取K和L基因，在MATTF插件上序列比對后進行基因串聯。用ModelFinder插件來選擇最佳的核苷酸替代模型GTR+F+I+G4，以6種合囊蕨目作為外類群, 利用MrBayes方法重建貝葉斯(bayesian inference, BI)樹。并且用MEGA分別構建鄰接(neighbor-joining, NJ)樹和最大似然(maximum likehood, ML)樹。綜合分析3種樹，結合PPGI系統發育關系進行手工調整，得出系統發育樹。

1.3 密碼子偏好性分析

以94種蕨類植物的基因為對象，使用在線軟件CUSP計算密碼子第一、第二、第三位的GC含量(分別用GC1、GC2、GC3表示)和總GC含量；CodonW v1.4.2軟件統計同義密碼子相對使用度(RSCU, relative synonymous codon usage)和有效密碼子(ENC, effective number of codon)。使用軟件SPSS 24對數據進行Pearson相關性分析。

GC3-GC12分析 密碼子的使用受到堿基的突變壓和自然選擇的雙重作用[11]，通過分析GC12和GC3的相關性，研究密碼子3個位置上的堿基組成是否相似和影響密碼子偏好性的主要因素。當GC12和GC3之間顯著相關，說明密碼子3個位置上的堿基組成無差異，密碼子的使用主要受突變的影響；當GC12和GC3之間相關性不顯著, 說明密碼子第一、第二和第三位堿基組成不相似, 密碼子的使用主要受自然選擇的影響。

ENC-plot分析 分析GC3和ENC的相關性,研究堿基組成對密碼子偏好性的影響。以GC3為橫坐標、ENC值為縱坐標繪制散點圖，并繪制ENC期望值標準曲線(ENC=2+GC3+29/[GC32+(1–GC3)2])進行比較。當基因散點分布在標準曲線附近，說明密碼子偏好性受突變影響較大；而散點位于曲線較遠下方位置時，密碼子偏好性受自然選擇影響較大。

最優密碼子 ENC值反映了同義密碼子非均衡使用的偏好程度，取值范圍為20～61。通常基因表達量越高，密碼子偏好性越強，ENC值越小[12]。分析94種蕨類植物基因主要偏好密碼子, 構建以ENC值為標準的高低表達庫(選擇ENC值最低和最高的10個基因)，計算RSCU值(ΔRSCU= RSCUHigh–RSCULow)，選擇RSCU>0.08且RSCU>1的作為最優密碼子。

1.4 分子進化速率分析

采用HyPhy v2.2.4軟件[13]基于核苷酸模型HKY85和密碼子模型MG94×HKY85_3×4計算各物種分支基因的轉換率(transition rate, trst)、顛換率(transversion rate, trsv)、轉換率和顛換率的比值(trst/trsv)、非同義替換率(nonsynonymous substi- tution rate, dN)、同義替換率(synonymous substitu- tion rate, dS)，以及非同義替換速率和同義替換速率的比值(=dN/dS)。通過軟件SPSS 24使用Mann- Whitney秩和檢驗對數據進行相關性分析。

1.5 適應性進化分析

使用PAML v4.9軟件[14]中的codeml程序檢測基因中是否存在氨基酸正選擇位點。采用檢測譜系之間選擇壓力的3種模型：分支模型(Branch model)、位點模型(Site model)和分支-位點模型(Branch-site model)。

分支模型中，單比率模型(one-ratio, M0)假定所有分支的值相同；自由比率模型(free-ratio, F)假定所有分支的值均不相同；二比率模型(two-ratio, M2)假定前景支與背景支的值不相同。通過似然比檢驗進行模型比較，M0和檢驗不同分支間值的差異性，M0和M2檢驗前景支和背景支間值的差異性。

位點模型假定不同位點的值不同，其中存在4對比較模型：M0 (one-ratio)和M3 (discrete)檢測位點間值是否一致；M1a (neutral)和M2a (selection), M7 (beta)和M8 (beta&)，M8和M8a (beta&=1)用于檢測是否存在正選擇位點。

分支-位點模型假定了分支間和位點間值都存在差異。對模型Model A和零假設進行比較，檢測前景支中的部分位點是否受到正選擇作用。

2 結果和分析

2.1 ycf94的基因特征分析

取樣的94種蕨類植物基因長度為222～156 bp。其中可以大致分為2種類型：以起始密碼子為T(ACG)開頭發生C to U RNA編輯的序列共28條，其中包括全部的鳳尾蕨科27條，大多為219 bp (共20條)；起始密碼子為M(AUG)的序列共66條, 包含除鳳尾蕨科外的18科，大多為216 bp (共50條)，其中6條合囊蕨目的大小均為177 bp。

2.2 94種蕨類植物的系統發育樹

利用K+L串聯基因構建NJ、ML和BI3種系統發育樹并對其進行比較，用蕨類分類系統PPGI進行修正。

在ML和BI中水龍骨亞目和鐵角蕨亞目為姐妹群且都是單系，與PPGI系統一致；而NJ中水龍骨亞目屬于多系。BI支持鐵角蕨科和腸蕨科為姐妹群, ML支持鐵角蕨科位于鐵角蕨亞目的基部位置。BI與PPGI系統最為接近，所以最后以BI為標準。構建出的系統發育樹(附錄2)主要為后續研究分析提供系統發育背景。

2.3 密碼子偏好性分析

2.3.1 密碼子堿基組成分析

附錄3為94種蕨類基因的密碼子堿基組成，其中GC含量為23.73%～47.75%，GC1、GC2和GC3分別為19.72%～47.22%、23.94%～47.22%和20.34%～56.94%，說明密碼子總體的堿基組成偏向A/U。ENC值為61～34.89，在近緣物種中基因的密碼子偏好性存在差異。其中有19條序列的ENC=61即每個密碼子都被均衡使用。ENC值以35作為判斷密碼子偏好性強弱的標準[15]，大部分的密碼子偏好性較弱。

對基因的密碼子相關參數進行熱圖統計，大多數的GC3要高于GC1、GC2和GC值；合囊蕨目的GC各項數值要低于其他蕨類，鳳尾蕨科的GC1和GC2總體上也低于其他蕨類(圖1)。

由表1可見，GC1、GC2、GC3和GC相互間均為極顯著相關，表明基因密碼子三位的堿基組成相似；而ENC值與各項GC值相關性不顯著，堿基組成對密碼子偏好性的影響較小。

表1 ycf94密碼子各參數的相關性分析

*:<0.05; **:<0.01

圖1 ycf94基因的密碼子相關參數熱圖統計

2.3.2 GC3-GC12分析

由圖2可見，大部分散點位于對角線之下，除天星蕨()外，大部分物種的基因的GC3值大于GC12。相關性分析結果表明，雙尾檢驗=0.00 (<0.01)，相關系數為0.739, GC3和GC12相關性極顯著，說明密碼子偏好性主要受到突變的影響。

2.3.3 ENC-plot分析

由圖3可見，散點的分布較為分散，且一部分落在距離曲線較遠的位置，則該部分基因密碼子偏好性受選擇影響較大。有40條序列的ENC比值為–0.1～0.1，表明ENC值與期望值相差較小，該部分基因的密碼子偏好性受突變影響較大。

圖2 中性繪圖

圖3 ENC-plot繪圖

2.3.4最優密碼子

從表2可見，共篩選出高頻密碼子(RSCU>1) 26個，高表達密碼子(△RSCU>0.08) 24個，而同時滿足這兩個條件的最優密碼子有14個：UUU、GUU、GUA、GUG、UCU、CCC、ACC、GCC、UAU、AAA、GAU、AGG、GGC、GGA。其中8個以A/U結尾，6個以G/C結尾。

2.4 進化速率分析

對鳳尾蕨科和其余蕨類的替換速率進行差異顯著性分析，數據表明(附錄4)，鳳尾蕨科除了同義替換率外的進化速率均值都大于其他蕨類(圖4)，秩和檢驗中顛換率，非同義替換率和的值分別是0.011、0.017和0.036 (<0.05)，說明鳳尾蕨科和其他蕨類基因的trsv、dN和存在顯著差異, 而trst、trsv/trst和dS無顯著差異。

表2 ycf94基因的密碼子RSCU統計以及最優密碼子

*: 最優密碼子; **: 高表達密碼子; ***: 高頻密碼子。

*: Optimal codon; **: High expression codon; ***: High frequency codon.

圖4 鳳尾蕨科與其他蕨類ycf94進化速率均值比較及Mann-Whitney檢驗

2.5 適應性進化分析

表3和4是PAML得出的基因適應性進化分析的結果。對各模型下計算出的對數似然值進行模型間的似然比檢驗，分析顯著性選擇更為適合的模型。結果表明，分支模型下M0和檢驗結果顯著(=0.024)，表明不同分支間值存在顯著差異。二比率模型中設定了以鳳尾蕨科為前景支，其他蕨類為背景支，M0和M2檢驗結果極顯著(=6.188× 10–6)，表明鳳尾蕨科和其他蕨類的值存在極顯著差異。

位點模型下，結果顯示位點間值存在顯著差異(=0.000)，且3對模型比較都顯示基因中存在正選擇位點。以BEB后驗概率大于95%為標準，模型M2a和M8均鑒定出1個正選擇位點74A。

在分支-位點模型中，同樣以鳳尾蕨科為前景支，似然比檢測結果顯示零假設模型和Model A間不存在差異(=1.000)，故不接受模型A的結果。

表3 不同模型下的參數估計值和對數似然值

**: 后驗概率<99%

**: Posterior probability <99%

表4 似然比檢驗

*:<0.05; **:<0.01

3 結論和討論

基因在蕨類中表現出高度保守性，但在種子植物中沒有發現同源序列這個重要特征，可能成為系統發育分析中的一個新的分子標記。同時，研究發現在蕨類中存在2種結構樣式：起始密碼子ACG發生RNA編輯的、基因長度多為219 bp的序列和以起始密碼子為AUG、基因長度多為216 bp的序列。其中全部的鳳尾蕨科基因序列都表現為前一種模式，以及三叉蕨科的也屬于同一結構，其余的序列均為后者。因此，的起始密碼子發生C to U RNA編輯雖然不屬于鳳尾蕨科的獨有結構，但RNA編輯作為表觀遺傳現象之一擁有可遺傳性，推測鳳尾蕨科的共同祖先中基因起始密碼子存在RNA編輯，該基因結構是鳳尾蕨科的共衍征。

編譯氨基酸過程中同義密碼子的使用是存在偏好性的，而物種在進化中會形成特有的密碼子使用模式，不同物種間的密碼子偏好性不同，而近緣物種間的密碼子使用模式有可能相似[16]。這種偏好性的差異對揭示物種間的進化關系，探究基因特征和蛋白質功能有重要意義。本研究中，基因的密碼子偏好性較弱，同義密碼子的使用傾向平均，其中有19條序列密碼子沒有使用偏好，說明的表達水平可能相對較低，非高表達基因。統計表明堿基組成整體偏向A/T，密碼子第三位也以A/U為主，這與葉綠體基因組中密碼子大多偏好A/U結尾這一特征相符合[17]，且GC3大多高于GC1和GC2。研究結果中94種蕨類的基因ENC值從34到61不等，其密碼子偏好性存在一定差異。通常，在沒有其他壓力作用下，同義密碼子的使用概率可能因為基因長度的限制而造成統計上的誤差，這種物種間的差異性可能與之相關。同時，GC3與GC12表現出顯著的相關性，密碼子不同位置下堿基組成的差異較小。而ENC與GC各項數值都不相關, 說明堿基組成對密碼子偏好性的影響較小。通過GC3-GC12和ENC-plot分析得出結論，基因的密碼子偏好性主要受突變影響，同時也受到選擇壓力和基因長度等其他因素的作用。

基于鳳尾蕨科和其他蕨類中基因的結構特征存在區別，對兩者的分子替換速率進行比較, 研究表明，鳳尾蕨科除了同義替換率外的進化速率均值都大于其他蕨類，且顛換率、非同義替換率和值間兩者存在顯著差異。鳳尾蕨科具有陸生、水生、旱生和附生等多個生態位[18]，有可能是生境的差異使在某一生態位上受到較高的選擇壓力，導致鳳尾蕨科的進化速率升高。

在基因的適應性進化分析中，分支模型下似然比檢驗支持的自由比率模型中檢測出15條>1的分支，說明該部分分支可能有正選擇存在，其中有7條分支屬于鳳尾蕨科，其余分支包括了蹄蓋蕨科、鱗毛蕨科、水龍骨科和合囊蕨科。同時似然比檢驗表示二比率模型優于單比率模型，鳳尾蕨科和其他蕨類的值存在極顯著差異，這與HyPhy的分子進化速率計算結果一致。但是前景支值(= 0.578 8)并不高于1，尚不能認為它是受正選擇作用的，選擇壓力約束放松可能是較為合理的解釋。基因共編碼80個氨基酸位點，位點模型M2a和M8均檢測出一個正選擇位點74A，占總序列1.25%; 基于M8模型檢測負選擇位點63個(<0.05), 占總序列78.75%。當一個基因受到正選擇壓力時，蛋白功能可能發生改變來適應外界環境變化；而當一個基因受到負選擇壓力，則說明該基因傾向保持原有的功能[19]。總的來說基因大部分位點受到負選擇作用，推斷其進化高度保守，在蕨類中基因的結構和功能基本趨于穩定。下一步分析氨基酸位點在蛋白質結構域中的分布，對正選擇位點的定位可以探究基因功能，進一步了解適應性進化發生的原因和作用。

本研究以94種蕨類植物的基因為對象,揭示了基因的基本結構特征，為葉綠體基因組的基因組成補充信息。在探討的堿基組成與密碼子偏好性的研究中，基因的密碼子偏好較弱,偏好使用A和U結尾的密碼子，其密碼子偏好性主要受到突變的影響，并且不同物種之間密碼子偏好性具有一定差異。同時，通過對基因的進化分析，發現該基因表現高度保守，為基因的功能研究提供參考。

[1] RAUBESON L A, JANSEN R K. Chloroplast genomes of plants [M]// HENRY R. Plant Diversity and Evolution: Genotypic and Phenotypic Variation in Higher Plants. Cambridge: CABI, 2005: 45–68. doi: 10. 079/9780851999043.0045.

[2] JANSEN R K, RUHLMAN T A. Plastid genomes of seed plants [M]// Genomics of Chloroplasts and Mitochondria: Advances in Photosyn- thesis and Respiration. New York: Springer, 2012: 103–126. doi: 10. 1007/978-94-007-2920-9_5.

[3] ZHANG X C, WEI R, LIU H M, et al. Phylogeny and classification of the extant lycophytes and ferns from China [J]. Chin Bull Bot, 2013, 48(2): 119–137. [張憲春, 衛然, 劉紅梅, 等. 中國現代石松類和蕨類的系統發育與分類系統[J]. 植物學報, 2013, 48(2): 119–137. doi: 10.3724/SP.J.1259.2013.00119.]

[4] PPG Ⅰ. A community-derived classification for extant lycophytes and ferns [J]. J Syst Evol, 2016, 54(6): 563–603. doi: 10.1111/jse.12229.

[5] WOLF P G, ROPER J M, DUFFY A M. The evolution of chloroplast genome structure in ferns [J]. Genome, 2010, 53(9): 731–738. doi: 10. 1139/G10-061.

[6] SONG M, KUO L Y, HUIET L, et al. A novel chloroplast gene reported for flagellate plants [J]. Am J Bot, 2018, 105(1): 117–121. doi: 10. 1002/ajb2.1010.

[7] CULLIS C A, VORSTER B J, VAN DER VYVER C, et al. Transfer of genetic material between the chloroplast and nucleus: How is it related to stress in plants? [J]. Ann Bot, 2009, 103(4): 625–633. doi: 10.1093/ aob/mcn173.

[8] HSU P Y, BENFEY P N. Small but mighty: Functional peptides encoded by small ORFs in plants [J]. Proteomics, 2018, 18(10): 1700038. doi: 10.1002/pmic.201700038.

[9] KEARSE M, MOIR R, WILSON A, et al. Geneious Basic: An integrated and extendable desktop software platform for the organi- zation and analysis of sequence data [J]. Bioinformatics, 2012, 28(12): 1647–1649. doi: 10.1093/bioinformatics/bts199.

[10] ZHANG D, GAO F L, JAKOVLI? I, et al. PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies [J]. Mol Ecol Resour, 2020, 20(1): 348–355. doi: 10.1111/1755-0998.13096.

[11] SUEOKA O. Total cross section measurements for positron and electron scattering on benzene molecules [J]. J Phys B: At Mol Opt Phys, 1988, 21(20): L631. doi: 10.1088/0953-4075/21/20/003.

[12] WRIGHT F. The ‘effective number of codons’ used in a gene [J]. Gene, 1990, 87(1): 23–29. doi: 10.1016/0378-1119(90)90491-9.

[13] POND S L K, FROST S D W, MUSE S V. HyPhy: Hypothesis testing using phylogenies [J]. Bioinformatics, 2005, 21(5): 676–679. doi: 10. 1093/bioinformatics/bti079.

[14] YANG Z H. PAML 4: Phylogenetic analysis by maximum likelihood [J]. Mol Biol Evol, 2007, 24(8): 1586–1591. doi: 10.1093/molbev/msm088.

[15] JIANG Y, DENG F, WANG H L, et al. An extensive analysis on the global codon usage pattern of baculoviruses [J]. Arch Virol, 2008, 153 (12): 2273–2282. doi: 10.1007/s00705-008-0260-1.

[16] ROMERO H, ZAVALA A, MUSTO H. Codon usage inis the result of strand-specific mutational biases and a complex pattern of selective forces [J]. Nucl Acids Res, 2000, 28(10): 2084–2090. doi: 10.1093/nar/28.10.2084.

[17] ZHOU M, LONG W, LI X. Patterns of synonymous codon usage bias in chloroplast genomes of seed plants [J]. For Stud China, 2008, 10(4): 235–242. doi: 10.1007/s11632-008-0047-1.

[18] SCHUETTPELZ E, SCHNEIDER H, HUIET L, et al. A molecular phylogeny of the fern family Pteridaceae: Assessing overall relation- ships and the affinities of previously unsampled genera [J]. Mol Phylo- genet Evol, 2007, 44(3): 1172–1185. doi: 10.1016/j.ympev.2007.04.011.

[19] PRINCE V E, PICKETT F B. Splitting pairs: The diverging fates of duplicated genes [J]. Nat Rev Genet, 2002, 3(11): 827–837. doi: 10. 1038/nrg928.

附錄1 蕨類物種信息

附錄2 94種蕨類植物的系統發育樹

附錄3基因的密碼子堿基組成與ENC值

附錄4基因分子進化速率

http://jtsb.ijournals.cn/ajax/common/download_attache_file.aspx?seq_id=20230218234211003&file_no

Molecular Evolution of Chloroplast Genein Ferns

LI Zifei1, SU Yingjuan2,3*, WANG Ting1*

(1. College of Life Sciences, South China Agricultural University, Guangzhou 510642, China; 2. School of Life Sciences, Sun Yat-Sen University,Guangzhou 510275, China; 3. Research Institute of Sun Yat-Sen University in Shenzhen, Shenzhen 518057, Guangdong, China)

As a newly discovered chloroplast gene in recent years,gene, with length of about 200 bp, is conserved among ferns. However, the origin and function ofrequire further investigation. The study was performed on thegene sequence of 94 species of ferns, in the phylogenetic background, analyzing codon usage bias, evolution rate and selection pressure. Results showed that the codon bias was weak ingene which the 3rd position of codon prefers A/U, variously in closely species. The codon bias was mainly produced by gene mutation. In addition, based on the differences in the structural characteristics ofbetween Pteridaceae and other ferns, their evolution rate was compared, suggesting that there were significant difference in transversion rate, nonsynonymous substitution rate and omega. Besides, only one position selection site 74A was detected. The strong negative selection pressure indicated that the function and structure ofwere mostly stabilized. The results provide a new clue of phylogenetic analysis in ferns and functional studies ofgene.

Fern;gene; Codon usage bias; Evolutionary rate; Adaptive evolution

10.11926/jtsb.4780

2023-02-18

2023-03-13

國家自然科學基金項目(31770587)資助

This work was supported by the National Natural Science Foundation of China (Grant No. 31770587).

李子菲(1998年生)，女，碩士研究生，研究方向為植物系統發育與分子進化。E-mail: fuyjun420@qq.com

通訊作者 Corresponding author.E-mail: suyj@mail.sysu.edu.cn; tingwang@scau.edu.cn