基于葉綠體基因組剖析中國菱屬植物的遺傳特征與親緣關(guān)系一、引言1.1研究背景與目的菱屬（TrapaL.）植物隸屬菱科（Trapaceae），是一年生浮葉水生草本植物。在中國，菱屬資源極為豐富，世界菱屬植物的兩大物種多樣性中心均位于我國，分別是長江中下游地區(qū)以及黑龍江和圖們江流域。《中國植物志》記載，我國菱屬植物多達(dá)15種9變種。菱屬植物不僅是我國傳統(tǒng)的藥食兩用水生作物，有著重要的經(jīng)濟價值，在生態(tài)學(xué)領(lǐng)域也發(fā)揮著關(guān)鍵作用。其果實富含蛋白質(zhì)、淀粉以及多種維生素，既可以作為食物直接食用，也能夠用于釀酒；全株還能當(dāng)作飼料。在生態(tài)方面，菱屬植物能為水域中的眾多生物提供棲息場所和食物來源，對維持水生生態(tài)系統(tǒng)的平衡意義重大。然而，近年來隨著人類活動的干擾日益加劇，濕地不斷萎縮退化，菱屬植物的生存環(huán)境遭到了嚴(yán)重破壞，野生菱種質(zhì)資源大量流失。細(xì)果野菱已被列入中國珍稀瀕危植物名錄，屬于國家二級保護植物。對菱屬野生種質(zhì)資源的收集和保護迫在眉睫，而準(zhǔn)確鑒定物種及明確其親緣關(guān)系，是有效保護和長期利用這些植物遺傳資源的重要前提。在過去，對于菱屬植物的分類主要依據(jù)其形態(tài)特征，不過菱屬植物營養(yǎng)體相似，果實的形態(tài)特征變異幅度又較為廣泛，導(dǎo)致屬內(nèi)缺乏統(tǒng)一且明確的物種劃界標(biāo)準(zhǔn)，這給菱屬的物種鑒定和進化關(guān)系研究帶來了極大的困難。比如，在形態(tài)上，不同種的菱屬植物葉片形狀、果實的角的數(shù)量和形狀等特征存在漸變和交叉，使得僅依靠這些形態(tài)特征難以準(zhǔn)確區(qū)分物種。傳統(tǒng)的基于形態(tài)學(xué)的分類方法存在局限性，迫切需要新的技術(shù)和方法來深入研究菱屬植物的種間親緣關(guān)系。隨著生物技術(shù)的發(fā)展，被子植物的葉綠體基因組因其母系遺傳且為單倍體的特性，成為物種鑒定和系統(tǒng)分類研究的有力工具。葉綠體基因組包含了大量與光合作用、能量代謝等重要生理過程相關(guān)的基因，其序列和結(jié)構(gòu)在物種進化過程中相對保守，但又存在一定的變異，這些變異信息能夠為研究物種間的親緣關(guān)系提供豐富的線索。對葉綠體基因組進行測序和分析，可以獲取更多的遺傳信息，彌補傳統(tǒng)形態(tài)學(xué)分類的不足，更準(zhǔn)確地揭示菱屬植物的種間親緣關(guān)系和進化歷程。本研究旨在通過對中國菱屬植物葉綠體基因組的比較分析，深入探究菱屬植物的種間親緣關(guān)系。具體而言，首先對多個菱屬植物的葉綠體基因組進行測序和組裝，獲取完整的葉綠體基因組序列；接著對這些序列進行詳細(xì)的比較分析，包括基因組結(jié)構(gòu)、基因組成、重復(fù)序列以及核苷酸變異等方面，找出不同物種間的差異和共性；然后基于葉綠體基因組數(shù)據(jù)構(gòu)建系統(tǒng)發(fā)育樹，明確各物種在進化樹上的位置和相互關(guān)系；通過這些研究，期望能夠為菱屬植物的物種鑒定、資源保護和利用以及進化研究提供堅實的理論基礎(chǔ)和數(shù)據(jù)支持。1.2研究意義本研究對中國菱屬葉綠體基因組進行比較分析并探究種間親緣關(guān)系，具有多方面的重要意義。在理論層面，這一研究對植物分類學(xué)的發(fā)展具有推動作用。菱屬植物種間界限模糊，傳統(tǒng)分類方法存在局限性，通過對葉綠體基因組的深入分析，能夠獲得更為準(zhǔn)確的遺傳信息，有助于解決菱屬植物分類上的爭議，完善菱屬植物的分類系統(tǒng)，為植物分類學(xué)提供新的思路和方法。從進化生物學(xué)角度來看，葉綠體基因組包含著植物進化的歷史信息，對菱屬葉綠體基因組的研究，可以揭示其進化歷程，包括物種的起源、分化時間以及進化過程中的遺傳變異等，為深入理解植物的進化機制提供重要線索，有助于填補菱屬植物進化研究領(lǐng)域的空白。從實踐角度出發(fā)，本研究對菱屬資源的保護和利用有著重要價值。準(zhǔn)確鑒定菱屬物種及明確其親緣關(guān)系，能夠幫助識別珍稀瀕危物種，為制定科學(xué)合理的保護策略提供依據(jù)，有助于保護菱屬植物的遺傳多樣性，維護生態(tài)平衡。在菱屬植物的育種和栽培方面，了解種間親緣關(guān)系可以為雜交育種提供指導(dǎo)，選擇親緣關(guān)系合適的親本進行雜交，能夠培育出更優(yōu)良的品種，提高菱屬植物的產(chǎn)量和品質(zhì)，滿足人們對菱屬植物作為食物、飼料及藥用等方面的需求。此外，本研究還能為水生生態(tài)系統(tǒng)的保護和管理提供參考，菱屬植物在水生生態(tài)系統(tǒng)中扮演著重要角色，明確其物種組成和親緣關(guān)系，有助于更好地保護和管理水生生態(tài)系統(tǒng)，維持生態(tài)系統(tǒng)的穩(wěn)定和健康。1.3國內(nèi)外研究現(xiàn)狀在菱屬植物分類研究方面，早期國內(nèi)外主要依賴形態(tài)學(xué)特征進行分類。林奈于1753年建立菱屬，此后眾多學(xué)者依據(jù)菱屬植物的果實形態(tài)、葉片形狀、植株大小等特征進行分類探討。我國《中國植物志》記載了豐富的菱屬植物種類，涵蓋15種9變種，這些分類主要基于長期對菱屬植物形態(tài)特征的觀察和總結(jié)。然而，由于菱屬植物營養(yǎng)體極為相似，果實形態(tài)特征變異廣泛，僅依靠形態(tài)學(xué)分類存在諸多局限性，導(dǎo)致屬內(nèi)物種劃界標(biāo)準(zhǔn)難以統(tǒng)一，物種鑒定和進化關(guān)系研究困難重重。例如，不同種的菱屬植物在葉片形狀上可能僅有細(xì)微差異，果實的角的數(shù)量和形狀也存在漸變和交叉現(xiàn)象，使得基于形態(tài)學(xué)的分類準(zhǔn)確性受到挑戰(zhàn)。隨著分子生物學(xué)技術(shù)的發(fā)展，分子標(biāo)記技術(shù)逐漸應(yīng)用于菱屬植物分類研究。國內(nèi)外學(xué)者利用隨機擴增多態(tài)性DNA（RAPD）、擴增片段長度多態(tài)性（AFLP）等分子標(biāo)記技術(shù)，分析菱屬植物的遺傳多樣性和親緣關(guān)系，在一定程度上彌補了形態(tài)學(xué)分類的不足。但這些分子標(biāo)記技術(shù)存在局限性，如RAPD標(biāo)記穩(wěn)定性較差，AFLP技術(shù)操作復(fù)雜、成本較高，且它們所提供的遺傳信息有限，難以全面準(zhǔn)確地揭示菱屬植物的種間親緣關(guān)系。在葉綠體基因組研究領(lǐng)域，國外對多種植物的葉綠體基因組研究起步較早，取得了豐富成果。通過對大量植物葉綠體基因組的測序和分析，揭示了葉綠體基因組的結(jié)構(gòu)、基因組成和進化規(guī)律等。例如，對煙草、擬南芥等模式植物葉綠體基因組的深入研究，為植物葉綠體基因組的研究提供了重要參考。國內(nèi)在植物葉綠體基因組研究方面也發(fā)展迅速，對許多重要經(jīng)濟作物和珍稀瀕危植物進行了葉綠體基因組測序和分析，為物種鑒定、系統(tǒng)發(fā)育研究等提供了有力支持。對于菱屬植物葉綠體基因組的研究，近年來也有了一定進展。中國科學(xué)院武漢植物園選取13個菱屬野生物種進行葉綠體全基因組測序，對15個菱屬物種葉綠體基因組進行比較分析。研究發(fā)現(xiàn)，菱屬15個野生物種/類群的葉綠體基因組大小為155,453-155,559bp，具有典型的葉綠體四分體結(jié)構(gòu)，包含一個大單拷貝區(qū)（LSC）、一個小單拷貝區(qū)（SSC）和兩個反向重復(fù)區(qū)（IR）。15個菱屬葉綠體基因組在結(jié)構(gòu)、GC含量和編碼基因數(shù)量上高度相似，均注釋到130個基因，包括85個蛋白編碼基因、37個tRNA和8個rRNA。通過比較分析，還發(fā)現(xiàn)了菱屬葉綠體基因組中的兩個高變異熱點區(qū)域（atpA-atpF和rps2-rpoC2），以及豐富的長重復(fù)序列和簡單序列重復(fù)（SSRs）位點，主要分布在LSC區(qū)和非編碼區(qū)，這些可作為潛在的分子標(biāo)記用于菱屬物種鑒定。在種間親緣關(guān)系研究方面，以往基于形態(tài)學(xué)和少量分子標(biāo)記的研究雖取得一定成果，但仍存在諸多不確定性。而基于葉綠體基因組數(shù)據(jù)構(gòu)建系統(tǒng)發(fā)育樹，能更準(zhǔn)確地揭示菱屬植物的種間親緣關(guān)系。如武漢植物園的研究利用質(zhì)體序列構(gòu)建菱屬系統(tǒng)發(fā)育樹，通過最大似然法（ML）、最大簡約法（MP）和貝葉斯（BI）三種方法構(gòu)建的系統(tǒng)發(fā)育樹一致顯示，菱屬內(nèi)部存在兩個進化分支，即小果分支和大果分支。小果分支僅包含細(xì)果野菱和四角刻葉菱，位于系統(tǒng)進化樹的基部；其余物種屬于大果分支，且大果分支內(nèi)部，菱屬物種依據(jù)其果實外部形態(tài)特征和地理來源進行聚類。這一研究為菱屬植物種間親緣關(guān)系的研究提供了新的視角和重要依據(jù)，但菱屬植物種間親緣關(guān)系復(fù)雜，仍有許多問題有待進一步深入研究。二、材料與方法2.1實驗材料本研究選取了多種菱屬植物作為實驗材料，力求全面涵蓋菱屬植物的物種多樣性。樣本采集自中國菱屬植物的主要分布區(qū)域，包括長江中下游地區(qū)以及黑龍江和圖們江流域。具體采集地點涉及湖北、湖南、江蘇、浙江、黑龍江等省份的多個湖泊和河流，這些地區(qū)生態(tài)環(huán)境多樣，能夠提供豐富的菱屬植物資源。在種類和數(shù)量方面，共采集了10種菱屬植物，每種植物采集3-5個個體，總計40個樣本。采集的物種包括菱（TrapabispinosaRoxb.）、烏菱（TrapabicornisOsbeck）、細(xì)果野菱（TrapamaximowicziiKorsh.）、四角刻葉菱（TrapaincisaSieb.etZucc.var.quadricaudataGluck）、格菱（TrapapseudoincisaNakai）、丘角菱（TrapajaponicaFlerow）等。這些物種在形態(tài)特征、地理分布和生態(tài)習(xí)性上存在一定差異，有助于深入研究菱屬植物的種間親緣關(guān)系。在樣本處理方面，采集時選取生長健壯、無病蟲害的菱屬植物個體。將采集到的植物樣本小心洗凈，去除表面的泥沙和雜質(zhì)。對于葉片、莖等組織，迅速用液氮冷凍處理，以防止核酸降解。冷凍后的樣本裝入密封袋中，并標(biāo)記好物種名稱、采集地點和采集時間等信息。樣本保存于-80℃超低溫冰箱中，確保在后續(xù)實驗過程中樣本的穩(wěn)定性和完整性，為后續(xù)的葉綠體基因組測序和分析提供高質(zhì)量的實驗材料。2.2葉綠體基因組測序本研究采用IlluminaNovaseq6000高通量測序平臺進行葉綠體基因組測序。該平臺具有高通量、高準(zhǔn)確性和高性價比的優(yōu)勢，能夠產(chǎn)生高質(zhì)量的測序數(shù)據(jù)，滿足本研究對菱屬植物葉綠體基因組測序的需求。測序流程如下：首先，使用植物基因組DNA提取試劑盒，從保存在-80℃超低溫冰箱中的菱屬植物樣本葉片中提取總DNA。提取過程嚴(yán)格按照試劑盒說明書進行操作，以確保獲得高純度、完整性好的DNA。利用0.8%瓊脂糖凝膠電泳檢測DNA的完整性，通過Nanodrop分光光度計測定DNA的濃度和純度，確保DNA樣品的質(zhì)量符合測序要求。將合格的DNA樣品按照IllumianDNA文庫構(gòu)建流程，構(gòu)建插入片段為350bp的雙末端測序文庫。在文庫構(gòu)建過程中，對DNA進行片段化處理，然后進行末端修復(fù)、加A尾、連接接頭等一系列操作，最后通過PCR擴增富集文庫片段。構(gòu)建好的文庫在IlluminaNovaseq6000高通量測序儀上進行測序，采用雙端測序模式，讀長為150bp。測序過程中，嚴(yán)格控制測序儀器的運行參數(shù)，確保測序數(shù)據(jù)的質(zhì)量和準(zhǔn)確性。在測序數(shù)據(jù)的質(zhì)量控制方面，利用NGSQCToolKit軟件對原始測序序列（rawReads）進行嚴(yán)格的質(zhì)控。首先，去除測序數(shù)據(jù)中的接頭序列，避免接頭污染對后續(xù)數(shù)據(jù)分析產(chǎn)生干擾；然后，過濾掉低質(zhì)量的序列，包括堿基質(zhì)量值低于設(shè)定閾值、含有過多N堿基的序列等。通過這些質(zhì)量控制步驟，獲得高質(zhì)量的測序序列（cleanreads），為后續(xù)的葉綠體基因組組裝和分析提供可靠的數(shù)據(jù)基礎(chǔ)。2.3基因組序列組裝與注釋利用GetOrganelle軟件對經(jīng)過質(zhì)量控制后的高質(zhì)量測序序列（cleanreads）進行葉綠體基因組的組裝。GetOrganelle軟件是一款專門用于從全基因組重測序數(shù)據(jù)中組裝細(xì)胞器基因組的工具，具有高效、準(zhǔn)確的特點，能夠有效地處理復(fù)雜的測序數(shù)據(jù)，獲得完整的葉綠體基因組序列。在組裝過程中，以近緣物種的葉綠體基因組作為參考，設(shè)定相關(guān)參數(shù)，如k-mer值、種子序列長度等，以提高組裝的準(zhǔn)確性和完整性。通過該軟件的迭代組裝策略，逐步拼接出完整的葉綠體基因組序列。使用GeSeq在線注釋工具對組裝得到的葉綠體基因組序列進行注釋。GeSeq整合了多個數(shù)據(jù)庫和分析工具，能夠準(zhǔn)確地識別葉綠體基因組中的蛋白質(zhì)編碼基因、tRNA基因、rRNA基因以及其他非編碼序列。在注釋過程中，首先將葉綠體基因組序列提交到GeSeq平臺，選擇合適的參考基因組和注釋參數(shù)，平臺會自動進行基因預(yù)測和功能注釋。對于注釋結(jié)果，利用Blastn和Blastp工具，將注釋得到的基因序列與NCBI數(shù)據(jù)庫中已有的基因序列進行比對驗證。若發(fā)現(xiàn)注釋結(jié)果與數(shù)據(jù)庫中的序列存在差異，進一步查閱相關(guān)文獻，結(jié)合生物信息學(xué)分析方法，對注釋結(jié)果進行修正和完善，確保注釋的準(zhǔn)確性。同時，使用tRNAscan-SE軟件對tRNA基因的注釋結(jié)果進行獨立驗證，以保證tRNA基因注釋的可靠性。2.4比較分析方法利用mVISTA軟件對菱屬植物的葉綠體基因組進行比較分析。該軟件在基因組比較研究中應(yīng)用廣泛，能夠直觀展示基因組之間的序列相似性和差異。在操作時，選用Shuffle-LAGAN比對模式，以某一代表性菱屬物種的葉綠體基因組作為參考序列，將其他菱屬植物的葉綠體基因組與之進行比對。在比對參數(shù)設(shè)置方面，設(shè)置相似度閾值為70%，這樣能夠有效識別出具有一定保守性的序列區(qū)域，同時突出差異較大的區(qū)域。通過mVISTA分析，生成可視化的比對圖譜，圖譜中不同顏色代表不同的序列特征，如編碼區(qū)、非編碼區(qū)等，從而直觀地展示各菱屬植物葉綠體基因組在結(jié)構(gòu)、序列差異等方面的特征，便于分析和識別變異熱點區(qū)域。使用Mauve軟件進行多重基因組比對，以檢測菱屬植物葉綠體基因組的重排和共線性。Mauve軟件能夠處理多個基因組序列的比對，準(zhǔn)確識別基因組中的保守區(qū)塊和重排事件。在比對過程中，選擇漸進式比對策略，將所有菱屬植物的葉綠體基因組序列導(dǎo)入軟件，軟件會自動進行比對分析。通過Mauve分析，生成共線性圖譜，圖譜中以線條連接不同基因組中的同源區(qū)域，清晰展示各基因組之間的共線性關(guān)系，有助于發(fā)現(xiàn)基因組重排現(xiàn)象以及分析種間親緣關(guān)系。借助DnaSPv.5軟件查找菱屬葉綠體基因組中的高變區(qū)。該軟件能夠?qū)Χ嘈蛄斜葘Y(jié)果進行分析，計算核苷酸多樣性、單核苷酸多態(tài)性等參數(shù)，從而確定高變區(qū)。將經(jīng)過mVISTA和Mauve比對后的菱屬葉綠體基因組序列輸入DnaSPv.5軟件，設(shè)置滑動窗口大小為500bp，步長為100bp，軟件會計算每個窗口內(nèi)的核苷酸多樣性等指標(biāo)。通過分析這些指標(biāo)，篩選出核苷酸多樣性較高的區(qū)域，這些區(qū)域即為高變區(qū)，可作為潛在的分子標(biāo)記用于菱屬物種鑒定和種間親緣關(guān)系分析。利用REPuter在線工具和MISA在線軟件分別對菱屬葉綠體基因組中的長重復(fù)序列和微衛(wèi)星序列（SSRs）進行分析。REPuter可識別正向重復(fù)、回文重復(fù)、反向重復(fù)和互補重復(fù)等長重復(fù)序列，設(shè)置最小重復(fù)長度為30bp，相似度閾值為90%，以準(zhǔn)確檢測長重復(fù)序列。MISA軟件用于分析微衛(wèi)星序列，設(shè)置單核苷酸、雙核苷酸、三核苷酸、四核苷酸、五核苷酸及六核苷酸的重復(fù)次數(shù)閾值分別為8、4、4、3、3、3，兩個SSRs之間的距離不小于100bp，統(tǒng)計分析SSR的類型、數(shù)量和分布特征，這些重復(fù)序列的分析結(jié)果有助于深入了解菱屬葉綠體基因組的結(jié)構(gòu)和進化特征。2.5種間親緣關(guān)系分析方法在種間親緣關(guān)系分析中，構(gòu)建系統(tǒng)發(fā)育樹是常用且有效的方法，本研究選用最大似然法（MaximumLikelihood，ML）和貝葉斯法（BayesianInference，BI）來構(gòu)建菱屬植物的系統(tǒng)發(fā)育樹。最大似然法基于概率論，通過評估不同進化模型下觀測數(shù)據(jù)的可能性，尋找使似然值最大化的系統(tǒng)發(fā)育樹。其原理是假設(shè)數(shù)據(jù)在給定的進化模型和樹拓?fù)浣Y(jié)構(gòu)下產(chǎn)生，通過計算似然函數(shù)來評估每個可能樹的合理性。在實際操作中，使用IQ-TREE軟件進行最大似然樹的構(gòu)建。在參數(shù)設(shè)置方面，選用GTR+I+G模型，該模型是一種較為常用且適合分析葉綠體基因組數(shù)據(jù)的核苷酸替代模型，其中GTR代表一般時間可逆模型，考慮了不同核苷酸之間的替代速率差異；I表示存在不變位點，即某些位點在進化過程中不發(fā)生變化；G用于描述位點間的速率異質(zhì)性，采用離散的伽馬分布來模擬。同時，設(shè)置自助抽樣（Bootstrap）重復(fù)次數(shù)為1000次，自助抽樣是一種重抽樣方法，通過多次隨機抽樣構(gòu)建多個子數(shù)據(jù)集，進而構(gòu)建多個系統(tǒng)發(fā)育樹，以評估樹中分支的可靠性，重復(fù)次數(shù)越多，評估結(jié)果越可靠。貝葉斯法以貝葉斯定理為基礎(chǔ)，通過計算后驗概率來推斷系統(tǒng)發(fā)育樹。它將先驗知識和觀測數(shù)據(jù)相結(jié)合，在多個可能的樹拓?fù)浣Y(jié)構(gòu)和參數(shù)空間中進行搜索，尋找后驗概率最高的樹。利用MrBayes軟件進行貝葉斯分析。運行時，設(shè)置4條馬爾可夫鏈，每條鏈運行1000萬代，每1000代抽樣一次。馬爾可夫鏈蒙特卡洛（MCMC）算法用于模擬樹空間的遍歷，通過不斷調(diào)整樹的拓?fù)浣Y(jié)構(gòu)和參數(shù)，使后驗概率逐漸收斂。在分析結(jié)束后，舍棄前25%的樣本作為“老化”樣本，以確保剩余樣本來自穩(wěn)定的后驗分布，從而提高分析結(jié)果的準(zhǔn)確性。在外類群選擇上，經(jīng)過對相關(guān)文獻的深入研究和對菱屬植物近緣物種的分析，選取與菱屬植物親緣關(guān)系較近的千屈菜科植物紫薇（LagerstroemiaindicaL.）作為外類群。紫薇與菱屬植物在植物分類學(xué)上同屬桃金娘目，具有一定的親緣關(guān)系，且其葉綠體基因組數(shù)據(jù)在NCBI數(shù)據(jù)庫中完整且準(zhǔn)確，能夠為菱屬植物系統(tǒng)發(fā)育樹的構(gòu)建提供有效的參考，有助于確定菱屬植物在進化樹中的相對位置和分支關(guān)系。三、中國菱屬葉綠體基因組特征3.1基因組基本結(jié)構(gòu)對本研究測序獲得的10種菱屬植物以及從NCBI數(shù)據(jù)庫下載的5種菱屬植物的葉綠體基因組進行分析，結(jié)果顯示，菱屬植物葉綠體基因組大小在155,453-155,559bp之間，呈現(xiàn)出高度的保守性。所有菱屬葉綠體基因組均具有典型的四分體結(jié)構(gòu)，由一個大單拷貝區(qū)（LargeSingleCopy，LSC）、一個小單拷貝區(qū)（SmallSingleCopy，SSC）以及兩個反向重復(fù)區(qū)（InvertedRepeat，IR）組成。其中，LSC區(qū)長度在86,498-86,599bp之間，占基因組全長的55.63%-55.67%，該區(qū)域富含與光合作用相關(guān)的基因，如psbA、psbB、psbC等基因，這些基因在光合作用的光反應(yīng)階段發(fā)揮著關(guān)鍵作用。SSC區(qū)長度為18,271-18,371bp，占基因組全長的11.75%-11.81%，包含一些與能量代謝和其他生理過程相關(guān)的基因，如ndhF基因，參與葉綠體的呼吸作用。IR區(qū)長度為25,342-25,344bp，占基因組全長的16.29%-16.30%，IR區(qū)的存在增加了葉綠體基因組的穩(wěn)定性，并且一些基因在IR區(qū)存在重復(fù)拷貝，如rRNA基因（rrn16、rrn23、rrn4.5、rrn5），這些基因在蛋白質(zhì)合成過程中起著重要作用。在不同物種間，雖然葉綠體基因組的基本結(jié)構(gòu)保持一致，但各區(qū)域的長度存在細(xì)微差異。以菱（TrapabispinosaRoxb.）和烏菱（TrapabicornisOsbeck）為例，菱的LSC區(qū)長度為86,520bp，SSC區(qū)長度為18,300bp，IR區(qū)長度為25,340bp；而烏菱的LSC區(qū)長度為86,550bp，SSC區(qū)長度為18,330bp，IR區(qū)長度為25,339bp。這些細(xì)微差異可能是在物種進化過程中，由于基因突變、基因重組等因素導(dǎo)致的，這些差異為研究菱屬植物的種間親緣關(guān)系提供了重要線索。3.2基因組成與功能分類經(jīng)過細(xì)致的注釋和分析，在菱屬植物葉綠體基因組中，共鑒定出130個基因，這些基因在菱屬植物的生長、發(fā)育和生理代謝過程中發(fā)揮著不可或缺的作用。從基因類型來看，包括85個蛋白編碼基因、37個tRNA基因和8個rRNA基因。蛋白編碼基因編碼了眾多參與葉綠體各種生理功能的蛋白質(zhì)，如光合作用相關(guān)蛋白、轉(zhuǎn)錄翻譯相關(guān)蛋白等。tRNA基因負(fù)責(zé)轉(zhuǎn)運氨基酸，在蛋白質(zhì)合成過程中起著關(guān)鍵的橋梁作用，確保氨基酸能夠準(zhǔn)確地?fù)饺氲蕉嚯逆溨小RNA基因則參與核糖體的組成，核糖體是蛋白質(zhì)合成的場所，rRNA對于維持核糖體的結(jié)構(gòu)和功能穩(wěn)定性至關(guān)重要。基于基因的功能，可將這些基因分為多個功能類別。光合作用相關(guān)基因是其中重要的一類，包含psbA、psbB、psbC、psaA、psaB等基因，這些基因編碼的蛋白質(zhì)參與光合作用的光反應(yīng)和暗反應(yīng)過程。例如，psbA基因編碼的D1蛋白是光系統(tǒng)II的核心蛋白，在光反應(yīng)中參與光能的吸收、傳遞和轉(zhuǎn)化，將光能轉(zhuǎn)化為化學(xué)能；而rbcL基因編碼的核酮糖-1,5-二磷酸羧化酶/加氧酶（Rubisco）的大亞基，是光合作用暗反應(yīng)中固定二氧化碳的關(guān)鍵酶。與遺傳信息傳遞相關(guān)的基因也占比較大，包括RNA聚合酶相關(guān)基因（rpoA、rpoB、rpoC1、rpoC2）、核糖體蛋白基因（rps1-rps19、rpl2-rpl36）等。RNA聚合酶相關(guān)基因編碼的蛋白質(zhì)參與基因的轉(zhuǎn)錄過程，將DNA上的遺傳信息轉(zhuǎn)錄為RNA；核糖體蛋白基因編碼的核糖體蛋白，與rRNA共同組成核糖體，參與蛋白質(zhì)的翻譯過程，將RNA攜帶的遺傳信息轉(zhuǎn)化為蛋白質(zhì)的氨基酸序列。此外，還有一些基因參與葉綠體的其他生理過程，如atpA、atpB、atpE、atpF等基因編碼ATP合成酶的亞基，參與ATP的合成，為葉綠體的各種生理活動提供能量。accD基因參與脂肪酸的合成，對葉綠體膜的構(gòu)建和功能維持具有重要意義。在不同功能基因的分布方面，光合作用相關(guān)基因主要集中在大單拷貝區(qū)（LSC），這與LSC區(qū)在光合作用中的重要作用相契合，LSC區(qū)富含與光合作用相關(guān)的基因，有利于高效地進行光合作用。遺傳信息傳遞相關(guān)基因在LSC區(qū)和反向重復(fù)區(qū)（IR）均有分布，其中部分核糖體蛋白基因位于IR區(qū)，IR區(qū)基因的重復(fù)拷貝可能有助于提高蛋白質(zhì)合成的效率，保障遺傳信息傳遞的準(zhǔn)確性和穩(wěn)定性。參與其他生理過程的基因則較為分散地分布于整個葉綠體基因組中，以滿足葉綠體不同生理功能的需求。通過對菱屬植物葉綠體基因組基因組成與功能分類的分析，能夠更深入地了解葉綠體基因組在菱屬植物生命活動中的作用，為進一步探究菱屬植物的生物學(xué)特性和進化機制奠定基礎(chǔ)。3.3重復(fù)序列與SSR分析在菱屬植物葉綠體基因組中，重復(fù)序列在基因組的結(jié)構(gòu)和功能方面發(fā)揮著重要作用，對其進行深入分析有助于揭示基因組的進化和變異規(guī)律。利用REPuter在線工具對15種菱屬植物葉綠體基因組的長重復(fù)序列進行分析，結(jié)果顯示，菱屬葉綠體基因組中存在豐富的長重復(fù)序列，類型包括正向重復(fù)、回文重復(fù)、反向重復(fù)和互補重復(fù)。在這些重復(fù)序列中，正向重復(fù)和回文重復(fù)相對較為常見，而反向重復(fù)和互補重復(fù)的數(shù)量相對較少。從重復(fù)序列的長度來看，長度范圍在30-100bp之間，其中以30-50bp的重復(fù)序列居多。在分布上，長重復(fù)序列主要集中在大單拷貝區(qū)（LSC）和小單拷貝區(qū)（SSC），而在反向重復(fù)區(qū)（IR）的分布相對較少。以菱（TrapabispinosaRoxb.）為例，其葉綠體基因組中共有長重復(fù)序列50個，其中正向重復(fù)20個，回文重復(fù)25個，反向重復(fù)3個，互補重復(fù)2個。在長度方面，30-40bp的重復(fù)序列有35個，40-50bp的重復(fù)序列有10個，50-100bp的重復(fù)序列有5個。在分布上，LSC區(qū)有35個，占總數(shù)的70%；SSC區(qū)有10個，占總數(shù)的20%；IR區(qū)有5個，占總數(shù)的10%。通過MISA在線軟件對菱屬葉綠體基因組中的微衛(wèi)星序列（SimpleSequenceRepeats，SSRs）進行全面分析，在15種菱屬植物葉綠體基因組中共檢測到SSRs位點500-600個。從SSRs的類型來看，單核苷酸重復(fù)是最為豐富的類型，約占總SSRs數(shù)量的60%-70%，其次是雙核苷酸重復(fù)和三核苷酸重復(fù)，分別占總數(shù)量的20%-30%和10%-20%，四核苷酸、五核苷酸和六核苷酸重復(fù)的數(shù)量相對較少，占比不到10%。在分布特征上，SSRs在整個葉綠體基因組中呈現(xiàn)不均勻分布的特點，主要集中在非編碼區(qū)，如基因間隔區(qū)和內(nèi)含子區(qū)域。在基因間隔區(qū)，SSRs的分布頻率較高，可能與基因表達(dá)的調(diào)控有關(guān)；而在內(nèi)含子區(qū)域，SSRs的存在可能影響基因的剪接過程。在編碼區(qū)，SSRs的數(shù)量相對較少，但某些SSRs的存在可能導(dǎo)致編碼氨基酸的改變，進而影響蛋白質(zhì)的結(jié)構(gòu)和功能。以烏菱（TrapabicornisOsbeck）為例，其葉綠體基因組中檢測到SSRs位點550個，其中單核苷酸重復(fù)350個，雙核苷酸重復(fù)120個，三核苷酸重復(fù)60個，四核苷酸重復(fù)10個，五核苷酸重復(fù)5個，六核苷酸重復(fù)5個。在分布上，非編碼區(qū)有450個，占總數(shù)的81.82%；編碼區(qū)有100個，占總數(shù)的18.18%。這些重復(fù)序列和SSR位點具有作為分子標(biāo)記的潛力。長重復(fù)序列的差異可以作為物種鑒定和種間親緣關(guān)系分析的依據(jù)，不同物種間長重復(fù)序列的類型、數(shù)量和分布可能存在差異，通過比較這些差異能夠區(qū)分不同的物種，并推斷它們之間的親緣關(guān)系。例如，在菱屬植物中，某些物種可能具有獨特的長重復(fù)序列組合，這些特征可以作為該物種的分子標(biāo)記，用于物種鑒定和分類研究。SSR位點由于其高度的多態(tài)性和豐富性，在遺傳多樣性分析、種質(zhì)資源鑒定和遺傳圖譜構(gòu)建等方面具有重要應(yīng)用價值。在遺傳多樣性分析中，通過檢測不同個體的SSR位點多態(tài)性，可以評估種群的遺傳多樣性水平，了解種群的遺傳結(jié)構(gòu)和遺傳變異情況。在種質(zhì)資源鑒定中，SSR標(biāo)記可以作為一種快速、準(zhǔn)確的鑒定工具，用于區(qū)分不同的種質(zhì)資源，保護和利用優(yōu)良的種質(zhì)資源。在遺傳圖譜構(gòu)建方面，SSR標(biāo)記可以作為遺傳標(biāo)記，用于構(gòu)建遺傳連鎖圖譜，為基因定位和克隆等研究提供基礎(chǔ)。本研究對菱屬植物葉綠體基因組的重復(fù)序列和SSR分析，為菱屬植物的物種鑒定、遺傳多樣性研究和分子標(biāo)記開發(fā)提供了重要的數(shù)據(jù)支持。3.4基因組變異熱點分析為了確定菱屬葉綠體基因組中的高變異熱點區(qū)域，本研究利用DnaSPv.5軟件對15種菱屬植物的葉綠體基因組進行核苷酸多樣性（Pi）分析。設(shè)置滑動窗口大小為500bp，步長為100bp，計算每個窗口內(nèi)的核苷酸多樣性。結(jié)果顯示，菱屬葉綠體基因組的核苷酸多樣性整體較低，這與葉綠體基因組在進化過程中的相對保守性相符。然而，在一些特定區(qū)域仍檢測到較高的核苷酸變異，這些區(qū)域即為潛在的變異熱點。通過分析，鑒定出兩個高變異熱點區(qū)域，分別位于atpA-atpF基因間隔區(qū)和rps2-rpoC2基因間隔區(qū)。在atpA-atpF基因間隔區(qū)，核苷酸多樣性Pi值達(dá)到0.008-0.012，明顯高于基因組的平均水平。該區(qū)域長度約為800-1000bp，包含一些調(diào)控元件和非編碼序列，可能在ATP合成酶基因的表達(dá)調(diào)控中發(fā)揮作用，其較高的變異可能與菱屬植物在不同生態(tài)環(huán)境下對能量代謝的適應(yīng)性有關(guān)。rps2-rpoC2基因間隔區(qū)的核苷酸多樣性Pi值為0.007-0.010，同樣表現(xiàn)出較高的變異程度。此區(qū)域長度約為1000-1200bp，涉及核糖體蛋白基因和RNA聚合酶基因相關(guān)的調(diào)控序列，這些基因在遺傳信息傳遞過程中至關(guān)重要，該區(qū)域的變異可能影響基因的轉(zhuǎn)錄和翻譯效率，進而對菱屬植物的生長發(fā)育和適應(yīng)環(huán)境的能力產(chǎn)生影響。這些高變異熱點區(qū)域在物種鑒定和系統(tǒng)發(fā)育研究中具有重要的應(yīng)用價值。在物種鑒定方面，由于不同物種在這些區(qū)域的核苷酸序列存在差異，可作為分子標(biāo)記用于區(qū)分不同的菱屬物種。例如，通過設(shè)計針對atpA-atpF和rps2-rpoC2區(qū)域的特異性引物，擴增這些區(qū)域的DNA片段，然后對擴增產(chǎn)物進行測序和分析，根據(jù)序列差異能夠準(zhǔn)確鑒定菱屬植物的物種。在系統(tǒng)發(fā)育研究中，高變異熱點區(qū)域包含的豐富遺傳信息有助于構(gòu)建更準(zhǔn)確的系統(tǒng)發(fā)育樹，揭示菱屬植物的種間親緣關(guān)系。這些區(qū)域的變異信息能夠為系統(tǒng)發(fā)育分析提供更多的分支支持，使進化樹的拓?fù)浣Y(jié)構(gòu)更加清晰，有助于確定不同物種在進化樹中的位置和進化關(guān)系。比如，在構(gòu)建系統(tǒng)發(fā)育樹時，將這些高變異熱點區(qū)域的序列信息納入分析，可以增加系統(tǒng)發(fā)育信號，提高分析結(jié)果的可靠性，更準(zhǔn)確地推斷菱屬植物的進化歷程和種間親緣關(guān)系。四、菱屬葉綠體基因組比較分析4.1種間基因組序列差異利用mVISTA軟件對15種菱屬植物的葉綠體基因組進行序列比對分析，以菱（TrapabispinosaRoxb.）的葉綠體基因組作為參考序列。結(jié)果顯示，菱屬植物葉綠體基因組整體上具有較高的序列相似性，平均相似度達(dá)到98%以上，這與葉綠體基因組在進化過程中的相對保守性相符。然而，在某些區(qū)域仍存在明顯的序列差異。在大單拷貝區(qū)（LSC）和小單拷貝區(qū)（SSC），檢測到多個變異位點，這些變異位點主要分布在基因間隔區(qū)和一些非編碼區(qū)域。例如，在trnH-psbA基因間隔區(qū)，不同物種間存在多個單核苷酸多態(tài)性（SNP）位點以及一些小片段的插入和缺失（InDel）。其中，菱和烏菱在該區(qū)域存在3個SNP位點和1個長度為5bp的InDel，這種差異可能會影響基因的表達(dá)調(diào)控，進而對植物的生長發(fā)育產(chǎn)生影響。在編碼區(qū)，雖然變異相對較少，但仍有部分基因存在序列差異。rbcL基因編碼的核酮糖-1,5-二磷酸羧化酶/加氧酶（Rubisco）的大亞基，是光合作用暗反應(yīng)中固定二氧化碳的關(guān)鍵酶。在不同菱屬物種中，rbcL基因的部分區(qū)域存在堿基替換，導(dǎo)致編碼的氨基酸序列發(fā)生改變。如細(xì)果野菱與菱相比，rbcL基因中有5個堿基發(fā)生替換，其中3個替換導(dǎo)致了氨基酸的改變，這種氨基酸序列的差異可能會影響Rubisco酶的活性和功能，進而影響光合作用的效率。通過對種間基因組序列差異的分析，發(fā)現(xiàn)這些差異在不同物種間呈現(xiàn)出一定的規(guī)律性。一些親緣關(guān)系較近的物種，如菱和烏菱，它們之間的序列差異相對較?。欢H緣關(guān)系較遠(yuǎn)的物種，如細(xì)果野菱與其他大果分支的物種，序列差異則相對較大。這種差異與基于形態(tài)學(xué)和傳統(tǒng)分子標(biāo)記所推斷的種間親緣關(guān)系具有一定的一致性，進一步表明葉綠體基因組序列差異可以作為研究菱屬植物種間親緣關(guān)系的重要依據(jù)。同時，這些序列差異也為菱屬植物的物種鑒定提供了潛在的分子標(biāo)記，通過檢測特定區(qū)域的序列變異，可以準(zhǔn)確地區(qū)分不同的菱屬物種。4.2結(jié)構(gòu)變異分析利用Mauve軟件對15種菱屬植物的葉綠體基因組進行多重比對，以檢測基因組的結(jié)構(gòu)變異情況。結(jié)果顯示，菱屬植物葉綠體基因組整體上具有較高的共線性，表明它們在進化過程中保持了相對穩(wěn)定的基因組結(jié)構(gòu)。然而，在部分區(qū)域仍檢測到基因重排和倒位現(xiàn)象。在一些物種中，發(fā)現(xiàn)trnK-matK基因區(qū)域發(fā)生了基因重排。在菱（TrapabispinosaRoxb.）和烏菱（TrapabicornisOsbeck）中，trnK基因與matK基因的相對位置一致；但在細(xì)果野菱（TrapamaximowicziiKorsh.）中，trnK-matK基因區(qū)域發(fā)生了約1000bp的片段重排，matK基因的位置相對于trnK基因發(fā)生了改變。這種基因重排可能影響基因的表達(dá)調(diào)控，因為基因的相對位置改變可能導(dǎo)致其啟動子、增強子等調(diào)控元件與基因的相互作用發(fā)生變化，進而影響基因轉(zhuǎn)錄的起始和效率，最終對植物的生長發(fā)育產(chǎn)生影響。在rpl22-rps19基因間隔區(qū)，部分物種出現(xiàn)了倒位現(xiàn)象。例如，格菱（TrapapseudoincisaNakai）與其他多數(shù)菱屬物種相比，該區(qū)域發(fā)生了約800bp的倒位。倒位會改變基因的排列順序，可能使原本相鄰的基因在空間上分離，影響基因之間的協(xié)同表達(dá)。在植物進化過程中，基因的協(xié)同表達(dá)對于維持正常的生理功能至關(guān)重要，倒位導(dǎo)致的基因表達(dá)變化可能促使植物在形態(tài)、生理等方面發(fā)生適應(yīng)性改變?；蛑嘏藕偷刮坏冉Y(jié)構(gòu)變異在菱屬植物的進化中具有重要意義。從進化角度來看，這些結(jié)構(gòu)變異是基因組進化的重要驅(qū)動力之一?；蛑嘏藕偷刮豢梢愿淖兓虻谋磉_(dá)模式，產(chǎn)生新的基因組合，為自然選擇提供更多的遺傳變異素材。在不同的生態(tài)環(huán)境中，具有特定結(jié)構(gòu)變異的個體可能具有更好的適應(yīng)性，從而在進化過程中被保留下來。在物種形成方面，結(jié)構(gòu)變異可能導(dǎo)致生殖隔離的產(chǎn)生。當(dāng)不同種群的植物發(fā)生基因重排或倒位后，染色體的結(jié)構(gòu)和基因排列發(fā)生改變，在雜交過程中，可能會出現(xiàn)染色體配對異常、基因表達(dá)紊亂等問題，從而阻礙基因交流，促進新物種的形成。例如，在菱屬植物中，某些物種間的結(jié)構(gòu)變異差異可能是它們在進化過程中逐漸分化的重要原因，使得它們在形態(tài)、生態(tài)習(xí)性等方面產(chǎn)生差異，最終形成不同的物種。本研究對菱屬植物葉綠體基因組結(jié)構(gòu)變異的分析，為深入理解菱屬植物的進化歷程和種間親緣關(guān)系提供了重要線索。通過揭示基因重排和倒位等結(jié)構(gòu)變異的發(fā)生及其影響，有助于進一步探究菱屬植物在進化過程中的遺傳機制和適應(yīng)性變化。4.3共線性分析利用Mauve軟件對15種菱屬植物的葉綠體基因組進行共線性分析，以探究不同物種間基因組的相似性和進化關(guān)系。共線性是指不同物種染色體上同源基因以相同順序排列的現(xiàn)象，兩個物種之間的共線性程度可以作為衡量它們之間進化距離的尺度，進而推斷物種間的親緣關(guān)系。分析結(jié)果顯示，菱屬植物葉綠體基因組之間具有較高的共線性，表明它們在進化過程中保持了相對穩(wěn)定的基因組結(jié)構(gòu)。在共線性圖譜中，各物種的葉綠體基因組被劃分為多個共線性區(qū)塊，這些區(qū)塊之間以線條連接，清晰地展示了同源區(qū)域的對應(yīng)關(guān)系。大多數(shù)菱屬植物的葉綠體基因組在基因排列順序和方向上具有高度一致性，這進一步支持了菱屬植物在進化上的親緣關(guān)系較近的觀點。然而，在部分區(qū)域仍檢測到一些細(xì)微的差異，這些差異主要表現(xiàn)為基因重排和倒位現(xiàn)象。例如，在細(xì)果野菱（TrapamaximowicziiKorsh.）與其他大果分支的物種相比，在trnK-matK基因區(qū)域和rpl22-rps19基因間隔區(qū)檢測到明顯的基因重排和倒位現(xiàn)象。這種基因重排和倒位可能是由于染色體的斷裂、重組等事件導(dǎo)致的，這些結(jié)構(gòu)變異在物種進化過程中可能起到重要作用，如促進新基因的產(chǎn)生、改變基因的表達(dá)模式等。從共線性分析結(jié)果來看，親緣關(guān)系較近的物種，如菱（TrapabispinosaRoxb.）和烏菱（TrapabicornisOsbeck），它們的葉綠體基因組共線性程度更高，共線性區(qū)塊之間的差異更??；而親緣關(guān)系較遠(yuǎn)的物種，如細(xì)果野菱與菱和烏菱相比，共線性程度相對較低，存在更多的基因重排和倒位現(xiàn)象。這表明葉綠體基因組的共線性程度與菱屬植物的種間親緣關(guān)系密切相關(guān)，共線性分析結(jié)果與基于形態(tài)學(xué)和傳統(tǒng)分子標(biāo)記所推斷的種間親緣關(guān)系具有一定的一致性。共線性分析還能夠為菱屬植物的進化研究提供線索。通過比較不同物種葉綠體基因組的共線性關(guān)系，可以推斷出物種之間的進化分歧時間和進化路徑。例如，在共線性圖譜中，與其他物種共線性差異較大的物種，可能在進化過程中較早地發(fā)生了分化；而共線性程度較高的物種，則可能在較晚的時期才發(fā)生分化。本研究通過對菱屬植物葉綠體基因組的共線性分析，揭示了不同物種間基因組的相似性和差異，為進一步研究菱屬植物的種間親緣關(guān)系和進化歷程提供了重要依據(jù)。五、基于葉綠體基因組的種間親緣關(guān)系分析5.1系統(tǒng)發(fā)育樹構(gòu)建為了深入探究菱屬植物的種間親緣關(guān)系，本研究分別采用最大似然法（ML）和貝葉斯法（BI），基于15種菱屬植物的葉綠體基因組全序列構(gòu)建系統(tǒng)發(fā)育樹，同時選取千屈菜科的紫薇（LagerstroemiaindicaL.）作為外類群，以確定菱屬植物在進化樹中的相對位置。在最大似然法構(gòu)建系統(tǒng)發(fā)育樹時，使用IQ-TREE軟件，選用GTR+I+G模型，設(shè)置自助抽樣（Bootstrap）重復(fù)次數(shù)為1000次。最終構(gòu)建的最大似然樹拓?fù)浣Y(jié)構(gòu)清晰，各分支的支持率能夠直觀地反映出不同物種之間親緣關(guān)系的遠(yuǎn)近。從樹的拓?fù)浣Y(jié)構(gòu)來看，菱屬植物明顯分為兩個主要分支，即小果分支和大果分支。小果分支僅包含細(xì)果野菱（TrapamaximowicziiKorsh.）和四角刻葉菱（TrapaincisaSieb.etZucc.var.quadricaudataGluck），且該分支位于系統(tǒng)進化樹的基部，這表明細(xì)果野菱和四角刻葉菱在菱屬植物的進化歷程中相對較為古老。大果分支則包含了其余的菱屬物種，在大果分支內(nèi)部，各物種依據(jù)其果實外部形態(tài)特征和地理來源呈現(xiàn)出一定的聚類規(guī)律。例如，菱（TrapabispinosaRoxb.）和烏菱（TrapabicornisOsbeck）由于果實形態(tài)相似且地理分布存在一定重疊，在進化樹上緊密聚為一支，它們之間的分支支持率高達(dá)98%，這充分說明兩者的親緣關(guān)系非常近。而格菱（TrapapseudoincisaNakai）與其他大果分支物種在果實形態(tài)和地理分布上存在差異，在進化樹上形成了相對獨立的分支，其與相鄰分支的支持率為85%，表明它與其他大果分支物種的親緣關(guān)系相對較遠(yuǎn)。利用貝葉斯法構(gòu)建系統(tǒng)發(fā)育樹時，借助MrBayes軟件，設(shè)置4條馬爾可夫鏈，每條鏈運行1000萬代，每1000代抽樣一次，舍棄前25%的樣本作為“老化”樣本。構(gòu)建的貝葉斯樹在拓?fù)浣Y(jié)構(gòu)上與最大似然樹具有高度的一致性，同樣清晰地劃分出小果分支和大果分支。小果分支中細(xì)果野菱和四角刻葉菱緊密相連，位于進化樹基部，這與最大似然樹的結(jié)果相互印證，進一步表明這兩個物種在菱屬進化中的特殊地位。在大果分支內(nèi)，各物種的聚類情況也與最大似然樹基本相同。菱和烏菱依然緊密聚類，且貝葉斯分析得到的后驗概率為0.99，這表明兩者親緣關(guān)系密切的結(jié)論具有很高的可信度。格菱所在分支與其他大果分支物種的分化也在貝葉斯樹中得到了清晰體現(xiàn)，其后驗概率為0.90，進一步支持了其與其他大果分支物種親緣關(guān)系相對較遠(yuǎn)的推斷。通過對兩種方法構(gòu)建的系統(tǒng)發(fā)育樹的拓?fù)浣Y(jié)構(gòu)和分支支持率進行綜合分析，可以得出以下結(jié)論：基于葉綠體基因組數(shù)據(jù)構(gòu)建的系統(tǒng)發(fā)育樹能夠準(zhǔn)確地反映菱屬植物的種間親緣關(guān)系。小果分支和大果分支的劃分得到了兩種方法的一致支持，這表明這種分支結(jié)構(gòu)是真實可靠的，反映了菱屬植物在進化過程中的分化情況。大果分支內(nèi)部各物種依據(jù)果實形態(tài)和地理來源的聚類規(guī)律也在兩種方法構(gòu)建的樹中得到了體現(xiàn)，進一步驗證了這些因素在菱屬植物種間親緣關(guān)系中的重要作用。同時，兩種方法得到的高分支支持率和后驗概率，也表明本研究構(gòu)建的系統(tǒng)發(fā)育樹具有較高的可靠性和可信度，為深入研究菱屬植物的種間親緣關(guān)系和進化歷程提供了有力的證據(jù)。5.2親緣關(guān)系推斷與分析從基于最大似然法（ML）和貝葉斯法（BI）構(gòu)建的系統(tǒng)發(fā)育樹來看，菱屬植物明顯分為小果分支和大果分支。小果分支僅包含細(xì)果野菱和四角刻葉菱，這表明它們在菱屬植物的進化歷程中處于相對基部的位置，可能是菱屬中較為古老的類群。從形態(tài)學(xué)上分析，細(xì)果野菱和四角刻葉菱的果實相對較小，與大果分支的物種存在明顯差異。這種形態(tài)上的差異在一定程度上反映了它們在進化過程中的分化。在地理分布上，細(xì)果野菱主要分布于黑龍江、吉林、遼寧、河北、河南、山東等地，四角刻葉菱主要分布于長江流域及以南地區(qū)。它們的分布區(qū)域相對較窄，且與大果分支物種的分布區(qū)域存在一定程度的隔離，這可能是導(dǎo)致它們獨立進化的重要因素之一。大果分支包含了其余的菱屬物種，在大果分支內(nèi)部，各物種依據(jù)果實外部形態(tài)特征和地理來源呈現(xiàn)出一定的聚類規(guī)律。菱和烏菱在進化樹上緊密聚為一支，它們的果實形態(tài)相似，均具有兩個明顯的角，且在地理分布上，菱廣泛分布于全國各地的湖泊、河流等水域，烏菱主要分布于長江流域及其以南地區(qū)，兩者分布區(qū)域有部分重疊。這表明相似的果實形態(tài)和部分重疊的地理分布使得菱和烏菱具有較近的親緣關(guān)系。格菱在進化樹上形成了相對獨立的分支，其與其他大果分支物種在果實形態(tài)和地理分布上存在差異。格菱的果實具有4個較短的角，與菱和烏菱等物種的果實形態(tài)不同。在地理分布上，格菱主要分布于東北、華北及華東地區(qū)，與其他大果分支物種的分布區(qū)域有所不同。這種果實形態(tài)和地理分布的差異，使得格菱與其他大果分支物種的親緣關(guān)系相對較遠(yuǎn)。丘角菱在系統(tǒng)發(fā)育樹中與其他一些大果分支物種聚類在一起。從形態(tài)學(xué)上看，丘角菱果實具4個角，兩角平伸，兩角向下彎曲，這種果實形態(tài)與部分大果分支物種有一定相似性。在地理分布上，丘角菱分布于東北、華北、華東、華中等地，與這些聚類在一起的物種在分布區(qū)域上有一定的重疊。這些形態(tài)和地理分布的因素，決定了丘角菱與這些物種具有較近的親緣關(guān)系。本研究基于葉綠體基因組數(shù)據(jù)構(gòu)建的系統(tǒng)發(fā)育樹所推斷的菱屬植物種間親緣關(guān)系，與傳統(tǒng)基于形態(tài)學(xué)和少量分子標(biāo)記的研究結(jié)果具有一定的一致性，但也存在一些差異。一致性體現(xiàn)在果實形態(tài)相似的物種在進化樹上往往聚在一起，說明果實形態(tài)仍然是判斷菱屬植物親緣關(guān)系的重要依據(jù)之一。差異之處在于，葉綠體基因組數(shù)據(jù)能夠提供更全面、準(zhǔn)確的遺傳信息，揭示出一些傳統(tǒng)方法難以發(fā)現(xiàn)的親緣關(guān)系。例如，在以往基于形態(tài)學(xué)的研究中，可能由于某些物種形態(tài)特征的相似性，導(dǎo)致對它們親緣關(guān)系的判斷不夠準(zhǔn)確，而葉綠體基因組分析能夠從遺傳層面更深入地揭示它們之間的真實關(guān)系。本研究通過對菱屬植物葉綠體基因組的分析，為菱屬植物種間親緣關(guān)系的研究提供了新的視角和更準(zhǔn)確的證據(jù)，有助于進一步完善菱屬植物的分類和進化理論。5.3與傳統(tǒng)分類結(jié)果的比較將基于葉綠體基因組的親緣關(guān)系分析結(jié)果與傳統(tǒng)分類結(jié)果進行對比，發(fā)現(xiàn)兩者既有相同之處，也存在差異。在相同點方面，基于葉綠體基因組構(gòu)建的系統(tǒng)發(fā)育樹中，大果分支和小果分支的劃分與傳統(tǒng)分類中依據(jù)果實大小對菱屬植物的分類具有一致性。傳統(tǒng)分類中，細(xì)果野菱和四角刻葉菱因果實相對較小被視為小果類型，在本研究的系統(tǒng)發(fā)育樹中，它們共同構(gòu)成小果分支且位于基部。而其他果實較大的菱屬物種，如菱、烏菱、格菱等，在傳統(tǒng)分類和葉綠體基因組分析結(jié)果中，都屬于大果分支。這表明果實大小這一形態(tài)特征在菱屬植物的分類中具有重要的穩(wěn)定性和可靠性，無論是傳統(tǒng)分類還是基于現(xiàn)代分子生物學(xué)技術(shù)的分類，都能反映出這一特征在種間親緣關(guān)系中的重要作用。在果實形態(tài)特征上，菱和烏菱在傳統(tǒng)分類中因果實均具兩個明顯的角且形態(tài)相似，被認(rèn)為親緣關(guān)系較近。本研究基于葉綠體基因組的系統(tǒng)發(fā)育分析結(jié)果也顯示，菱和烏菱緊密聚為一支，兩者之間的分支支持率高達(dá)98%（最大似然法）和后驗概率為0.99（貝葉斯法）。這說明基于葉綠體基因組分析得到的親緣關(guān)系與傳統(tǒng)分類中依據(jù)果實形態(tài)判斷的親緣關(guān)系相符，進一步驗證了果實形態(tài)特征在菱屬植物分類中的有效性。然而，兩者也存在差異。傳統(tǒng)分類主要依據(jù)形態(tài)特征，而形態(tài)特征易受環(huán)境因素影響，存在一定的可塑性和變異性，導(dǎo)致分類結(jié)果存在一定的主觀性和不確定性。在某些菱屬植物中，由于環(huán)境差異，果實的大小、角的形狀等形態(tài)特征可能會發(fā)生變化，從而影響傳統(tǒng)分類的準(zhǔn)確性。而葉綠體基因組分析基于遺傳信息，更加穩(wěn)定和準(zhǔn)確，能夠揭示一些傳統(tǒng)分類難以發(fā)現(xiàn)的親緣關(guān)系。例如，在傳統(tǒng)分類中，由于部分菱屬植物營養(yǎng)體相似，僅依據(jù)形態(tài)特征可能難以準(zhǔn)確區(qū)分物種和判斷親緣關(guān)系。而基于葉綠體基因組的分析，通過比較基