基于轉錄組測序初步揭示天麻生長代謝的分子機制

劉云霞，狄永國，仇全雷，肖舒卉，譚彧文，陳麗梅，徐慧妮，李昆志*

1.昆明理工大學生命科學與技術學院，云南 昆明 650500

2.昭通市昭陽區農業局 種植業科，云南 昭通 657000

3.西藏波密高原藏天麻產業開發有限公司，西藏 林芝 860300

天麻是名貴傳統中藥材，具有息風止痙、平抑肝陽、祛風通絡之功效，主治小兒驚風、癲癇抽搐、破傷風、頭痛眩暈、風濕痹痛等癥[1-3]。天麻為藥食兩用的蘭科植物，具有極高的藥用和食用價值，擁有巨大的潛在市場，天麻需求量大。野生天麻活性成分含量高、藥效好，但產量低，不能滿足實際需求，且由于人工過度采集瀕臨滅絕，天麻供給主要靠人工種植。因此，天麻栽培成為天麻領域的研究熱點。

天麻生長經歷種子、原球莖、米麻、白麻和箭麻5個生長發育階段[4]。Kusano[5]首次報道了天麻與木材腐爛病原菌蜜環菌存在真菌菌根關系。隨后的研究表明，天麻種子萌發需要萌發菌如小菇屬Mycena dendrobiiL.FanetS.X.Guo 侵入種胚細胞，為天麻種子萌發提供營養物質，促進天麻種子萌發長出幼小的原球莖[6]。隨著原球莖進一步生長，其前端分生組織處可直接分化出營養繁殖莖米麻。天麻可消化萌發菌和蜜環菌入侵的菌絲作為其營養來源[7-9]，米麻從侵入的萌發菌和蜜環菌獲得營養發育成營養繁殖莖白麻[10]。在接下來的生長中白麻的直接營養物質是蜜環菌，蜜環菌侵入白麻后，白麻消化侵入的蜜環菌菌絲獲得營養后進一步長出商品麻箭麻[11]。最近，研究者通過天麻轉錄組學分析，闡明了天麻種子萌發的機制[12-14]。譚彧文等[15]通過轉錄組學的方法發現侵入到天麻體內的蜜環菌其胞外酶和抗氧化酶的差異表達基因（differentially expressed genes，DEGs）98.31%表達下調，說明侵入天麻體內的蜜環菌活性減弱，不能進一步入侵天麻，且不受到天麻脅迫作用，使蜜環菌和天麻處于共生狀態，初步揭示了天麻與蜜環菌的共生機制。袁媛等[16]測出天麻的全基因組序列發現，大部分與光合作用和抗病菌相關的基因如核酸結合部位（nucleotide-binding site，NBS）類基因發生了大量丟失，同時發現天麻的線粒體基因組擴大且單子葉甘露糖結合凝集素抗真菌蛋白（gastrodia antifungal protein，GAFP）基因的數量增加，且發現天麻與蜜環菌共生關系建立的重要信號是獨角金內酯，這些結果展現了完全異養植物天麻是如何通過實現廣泛的基因收縮甚至丟失、擴張以及基因的新功能化來完成其獨特的生長特征，更深一步地闡述了天麻與蜜環菌的共生關系。文歡等[17]也通過轉錄組學的方法發現了天麻內存在較完整的天麻苯丙烷類產物合成代謝通路。到目前為此，天麻各生長發育階段的生理形態特征比較清楚，蜜環菌能入侵營養繁殖莖白麻與其共生，但共生后進一步生長為箭麻的生長代謝特征仍不清楚。

本實驗通過轉錄組分析，在轉錄組水平的基礎上研究箭麻和共生天麻生長代謝過程中相關差異基因表達水平變化，揭示兩者間的基因表達水平的差異，以期闡釋箭麻和共生天麻間的生長代謝特征，為進一步研究天麻不同生長發育階段的生理生化代謝特征奠定基礎,同時為天麻栽培提供理論指導。

1 材料與儀器

1.1 材料

選取昭通產烏天麻作為研究材料，經筆者鑒定為烏天麻Gastrodia.e lataBl.(F.gl auca) S.Chow。將蜜環菌Armillaria.mellea(Vahl) P.Kumm.接種于蘋果樹枝一段時間待長出蜜環菌菌絲后接種于營養繁殖莖白麻，成為共生天麻，接種3 個月長出箭麻后分別取2 種樣品，立即放入液氮中冷凍后儲存入-80 ℃冰箱中用于總RNA 提取。

1.2 儀器

Bio-rad CFX96 型熒光定量PCR 儀（美國伯樂公司）；DHP-9502 型電熱恒溫培養箱（上海一恒科學儀器有限公司）；heraeus mutifuge X1R 型實驗室用離心機（賽默飛世爾科技有限公司）。

2 方法

2.1 RNA 提取及質量檢測

采用Trizol Reagent（Invitrogen）提取箭麻和共生天麻的RNA。將約0.1 g 樣品放入用液氮預冷的研缽中，磨碎后加入1 mL 的Trizol 提取液于研缽中繼續研磨至研磨液呈紅色透明狀，于室溫靜置5 min 后移入1.5 mL 離心管并加入0.2 mL 氯仿振蕩混勻；隨后在12 000 r/min、4 ℃離心15 min 后吸取上清液至新的離心管中；在加入0.25 mL 異丙醇和0.25 mL 氯化鈉與檸檬酸鈉混合高鹽液，混勻-20 ℃放置30 min后在12 000 r/min、4 ℃條件下離心30 min，棄上清，沉淀用1 mL（-20 ℃ 預冷）的75%乙醇清洗，隨后在7500 r/min、4 ℃條件下離心5 min，棄乙醇，重復清洗2 次；然后將沉淀在室溫條件下自然晾干，最后用20 μL 焦碳酸二乙酯（diethyl pyrocarbonate，DEPC）處理水溶解RNA。用分光光度計測量樣品RNA 的純度和濃度，檢測RNA 樣品的完整性。

2.2 轉錄組測序樣品制備及文庫建立

提取樣品總RNA 后，將破碎緩沖液加入到帶有Oligo（dT）磁珠富集的真核生物mRNA 中將其打斷成短片段。以mRNA 為模板，用六堿基隨機引物合成第1 條cDNA 鏈，隨后加入緩沖液、RNase H、DNA 聚合酶I 和dNTPs 以合成第2 條cDNA 鏈。在經過QiaQuick PCR 試劑盒純化并加EB 緩沖液洗脫之后做末端修復、加poly（A）并連接測序接頭。最后用瓊脂糖凝膠電泳進行片段大小選擇，并進行PCR 擴增，建好的測序文庫用Illumina HiSeqTM2000 進行高通量測序分析。

2.3 測序數據處理及分析

經過統計箭麻和共生天麻分別獲得平均28.97 Mb 和67.73 Mb 條原始測序序列（raw reads）。對raw reads過濾低質量數據和測序引物、接頭等人工序列，得到clean reads。對過濾處理得到的clean reads 進行Trinity 組裝生成轉錄本（transcript）和單基因簇（unigene）。再將Unigene 與Swiss-Prot（瑞士-波特蛋白質序列數據庫）、非冗余蛋白庫（non-redundant protein sequence database，NR）、京都基因和基因組百科全書（kyoto encyclopedia of genes and genomes，KEGG）和蛋白相鄰類的聚簇（cluster of orthologous groups of proteins，KOG）數據庫比對，獲得Unigene 功能注釋信息。

2.4 表達差異分析

使用RPKM 值表示Unigene 的表達豐度。滿足false discovery rate（FDR）＜0.05 且log2|FC|≥2 條件的基因為差異表達基因（differentially expressed genes，DEGs）。隨后對差異表達基因做基因本體論（Gene Ontology，GO）功能富集分析、KOG 注釋和KEGG 注釋與分類。

2.5 實時熒光定量PCR (qRT-PCR) 分析

為了驗證轉錄組數據中基因差異表達的正確性，選取轉錄組基因文庫內與生長代謝過程中氮代謝相關的谷氨酰胺合成酶（glutamine synthetase，GS）、天冬酰胺合成酶（asparagine synthetase，ASNS）、谷氨酸脫氫酶（glutamate dehydrogenase，GDH）、精氨酸酶（arginase，Argase）；碳代謝相關的蔗糖合成酶（sucrose synthase，SuSy）、可溶性淀粉合成酶（soluble starch synthase，SSS）、ADP-葡萄糖焦磷酸化酶（ADP-glucose pyrophosphorylase，AGPase）、淀粉酶（amylase，AMS）和與能量代謝相關的丙酮酸激酶（pyruvate kinase，PK）、己糖激酶（hexokinase，HK）、丙酮酸脫氫酶（pyruvate dehydrogenase，PDH）、蘋果酸脫氫酶（malate dehydrogenase，MDH）、異檸檬酸脫氫酶（isocitrate dehydrogenase，IDH）、ATP 酶（adenosine triphosphatase，ATPase）14個具有代表性的基因，用Primer 6 設計引物，以β-actin作為內參基因進行qRT-PCR 相對表達分析[16]，根據其cDNA 片段設計特異性引物（表1）。提取RNA后，用HiScript II 1stStrand cDNA Synthesis Kit 試劑盒進行反轉錄。將反轉錄產物稀釋10 倍，用ChamQTM Universal SYBR ? qPCR Master Mix 做熒光定量。熒光定量在20 μL 體系含有10 μL 2×ChamQTM Universal SYBR qPCR Master Mix，2 μL cDNA 模板和0.4 μL 每種基因特異性引物的反應混合物中進行。最后使用CFX96TMReal-Time System 及其相對定量軟件進行2次生物重復3次技術重復。反應參數為95 ℃、30 s，并以95 ℃、10 s，58 ℃、30 s 循環40 次。將cDNA 文庫分別標準化為參考基因β-actin。通過2-ΔΔCt方法計算基因相對表達量。

3 結果與分析

3.1 測序數據組裝結果

3.1.1 數據過濾統計表 測序后在G 樣本中平均獲得28.36 Mb 的高質量clean reads，占raw reads 的97.92%，在GA-G 樣本平均獲得66.27 Mb 的高質量的clean reads，占raw reads的97.84%。4個樣本clean reads Q30均大于92%，平均GC 含量為47.28%（表2）。對片段進行拼接得72 244 個Unigene 轉錄物。

3.1.2 組裝質量統計 對片段進行拼接得到72 244個Unigene 轉錄物。組裝結果質量評估可從N50數值來評估。將所有Unigene 從長到短排序，并依次累加長度。當累加片段長度達到總片段長度（所有Unigene 的長度）的50%時，對應那個片段的長度和數量，即為Unigene N50長度和數量。Unigene N50越長，數量越少，說明組裝質量越好。在本實驗中，由組裝結果得N50 長度為1733 bp，其中Unigene 序列最長為15 844 bp，最短為201 bp，平均長度為916 bp，長度在1000 bp 以上的Unigene有18 698 個，占總數的25.88%（表3）。

表1 qRT-PCR 擴增特異引物Table 1 Specific primers of qRT-PCR amplification

表2 數據過濾統計表Table 2 Data filter statistics table

表3 Unigene 長度分布統計Table 3 Length statistics of Unigene distribution

3.1.3 轉錄物功能注釋及分類 為了解轉錄物Unigene 序列信息，將組裝好的72 244 個Unigene序列比對到Swissprot、Nr、KEGG 和KOG 4 個數據庫中，有26 312 個Unigenes 獲得注釋信息，所占比例為36.42%。其中Unigenes 注釋成功最多的為Nr 和Swissprot 數據庫，分別為26 217 和17 886個。注釋在KOG 和KEGG 數據庫中的Unigenes 個數分別為16 625 和9 712（表4）。

表4 Unigene 功能注釋結果統計Table 4 Unigene functional annotation statistical results

3.2 基因差異表達分析

3.2.1 主成分分析（ principal component analysis，PCA） 在RNA 組學研究中利用PCA，將樣本所包含的上萬個維度的信息（上萬個基因的表達量），降維為數個維度的綜合指標（主成分），以便于進行樣本間的比較，同時保證原始數據中包含信息盡可能多地被保留[18]。使用R 語言包（http://www.r-project.org/）對共生天麻和箭麻所有樣品進行PCA，如圖1 所示。根據共生天麻和箭麻在第1 主成分（PC1）和第2 主成分（PC2）2 個綜合指標中的數值大小，做二維坐標圖。從圖1 可以看出共生天麻-1 和共生天麻-2 分布在右側，箭麻-1 和箭麻-2 分布在左側，說明2 種樣品具有一定的差異性。PCA 分析中PC1（95.1%）和PC2（4.2%）對樣品中所有基因表達量總體方差貢獻率為99.3%，表明通過PCA 分析可以很好地區分箭麻和共生天麻。

圖1 共生天麻和箭麻樣品的PCAFig.1 Principal component analysis of symbiotic G.elata and mature tuber

3.2.2 分組間差異基因火山圖 比較箭麻與共生天麻的測序結果，以FDR＜0.05 且|log2FC|>1 為判斷標準，檢測到共有12 498 條基因發生顯著差異表達，其中上調表達基因數為9 000，占差異表達基因個數的72.01%；下調表達基因數為3 498占差異表達基因數的27.99%。差異表達直觀展示火山圖見圖2。

圖2 差異表達基因火山圖Fig.2 Differential expressed gene volcanic map

3.2.3 差異表達基因的GO 分類 如表5 所示，有10 473 條差異表達基因成功進行了GO 分類。在GO 分類的3 個大類中有4723 條差異表達基因與生物學過程有關，占比為45.10%，其中富集最多的是代謝過程亞類，占比對到該功能類別總數的23.65%；有2508 條差異表達基因與分子功能相關，占比23.95%，其中富集最多的是催化活性亞類，共比對到該功能類別總數46.93%；與細胞成分有關的差異表達基因有3242 條，占比30.96%，富集最多的是細胞和細胞成分亞類，共比對到該功能類別總數的46.33%。

3.2.4 差異表達基因的KEGG 功能分析 箭麻與共生天麻樣本間的差異基因進行通路富集分析發現，有803 個差異基因比對到20 個KEGG 通路中。有較多基因歸類為苯丙素生物合成（phenylpropanroid biosynthesis）、淀粉與蔗糖代謝（starch and sucrose metabolism）、植物激素信號傳導（plant hormone signal transduction）、碳代謝（carbon metabolism）、糖酵解/葡萄糖酵素（glycolysis/gluconeogenesis）和氮代謝（nitrogen metabolism）這幾類相關的通路之中，包含了生長過程中氮代謝、碳代謝和能量代謝（如糖酵解）等途徑，說明蜜環菌侵入天麻與天麻共生長出箭麻過程中，箭麻與共生天麻生長代謝過程中相關酶基因差異表達，兩者的代謝強弱有很大差別（表6）。

表5 箭麻和共生天麻差異表達基因GO 注釋Table 5 GO annotation of differentially expressed genes in mature tuber and symbiotic G.elata.

表6 箭麻和共生天麻差異基因的顯著富集的代謝途徑Table 6 Significant enrichment pathway of differentially expressed genes in mature tuber and symbiotic G.elata.

3.2.5 差異表達基因的KOG 功能注釋 在實驗中共有26 980 個DEGs 注釋到KOG 的25 個分類中，其中7796 個Unigenes 注釋為一般功能預測（general functional predictio only，R），占比28.90%；2782個Unigenes 注釋為翻譯后修飾、蛋白質分解及伴侶分子（ posttranslational modification,protein turnover，chaperenes，O），占比10.31%；2295 個Unigenes 注釋為信號轉導（signal transduction mechanism，T），占比8.51%；1519 個Unigenes注釋為RNA 的加工和修飾（RNA processing and modification，A），占比5.63%；774 個Unigenes注釋為能量的產生和轉化（energy production and conversion，C），占比2.87%；816 個Unigenes 注釋為碳水化合物轉運與代謝（carbohydrate transport and metabolism，G），占比3.02%，628 個Unigenes注釋為次生代謝產物的合成、轉運和代謝（secondary metabolites biosynthesis，transport and catabolism，Q），占比2.32%（表7）。

3.2.6 箭麻和共生天麻中生長代謝關鍵差異基因表達分析 根據差異基因GO 分類、KEGG 功能分析和KOG 功能注釋，對箭麻和共生天麻生長代謝過程的相關基因進行篩選。在箭麻和共生天麻轉錄組測序數據注釋的結果中篩選到與氮代謝、碳代謝和能量代謝相關酶的基因序列。在氮代謝途徑中，選取了差異表達的基因GS、ASNS、GDH、Argase；在碳代謝中選取了差異表達的基因SuSy、SSS、AGPase、AMS；在能量代謝中選取了差異表達的基因PK、HK、PDH、MDH、IDH、ATPase（表8）。從這些基因的表達趨勢來看，與合成氮，碳和能量的酶基因大部分呈上調，除分解能量的酶基因上調外，分解氮，碳的酶基因呈下調。說明蜜環菌侵染天麻與天麻共生時共生天麻各種代謝旺盛，合成各種有機物和能量物質。

表7 差異表達基因的KOG 功能注釋Table 7 DEGs KOG annotation

3.2.7 箭麻和共生天麻中與各種代謝相關基因相對表達水平qRT-PCR 分析 選取的14 個與生長代謝過程相關關鍵差異酶基因的相對定量qRT-PCR分析結果見圖3。由2-ΔΔCt值計算箭麻和共生天麻中的相對表達量，統計分析結果顯示箭麻和共生天麻的這些基因表達水平與轉錄組分析結果一致。即與箭麻相比，共生天麻中合生碳、氮，能量物質的顯著差異酶基因呈上調表達，除分解ATP的ATPase基因呈上調表達外，分解氮、碳的GDH、Argase和AMS基因呈下調表達。因此，本課題組認為在轉錄組分析中用RPKM（Reads Per Kilobase per Million mapped reads，每百萬reads 中來自于某基因每千堿基長度的reads 數）法推斷的基因表達水平是可靠的。

表8 箭麻和共生天麻中與生長代謝相關關鍵差異基因的表達情況Table 8 Expression of key differentially expressed genes in growth and metabolism of mature tuber and symbiotic G.elata.

圖3 箭麻和共生天麻中生長代謝相關基因相對表達水平分析Fig.3 Analysis of relative expression level of related genes in growth and metabolism in symbiotic G.elata and mature tuber

4 討論

4.1 共生天麻和箭麻的轉錄組測序分析

本研究利用Illumina HiseqTM2000 測序技術箭麻和共生天麻進行測序，箭麻和共生天麻分別獲得4.23 Gb 和9.90 Gb 有效數。通過數據組裝獲得72 244 條Unigenes，N50長度為1773 bp，平均長度為916 bp 的非冗余Unigene，且Q30值均大于92%，說明組裝數據良好可信。同時將組裝好的72 244 個Unigene 序列比對到Swissprot、Nr、KEGG和KOG 4 個數據庫中有26 312 個Unigene 獲得注釋信息，所占比例為36.42%。將箭麻和共生天麻進行PCA，發現箭麻和共生天麻在一定程度上有所差異，通過PCA 可以很好地區分箭麻和共生天麻。在FDR＜0.05 和|log2FC|>1 的條件下，比較兩組測序數據的結果，發現12 498 條DEGs，上調表達9000條，下調表達3498條。GO數據庫將10 473條差異表達基因分為生物學過程、分子功能和細胞組成3 大類，其中代謝過程、細胞過程催化活性等功能小類富集較多。箭麻與共生天麻樣本中有803 個DEGs 比對到20 個KEGG 通路中，與箭麻和共生天麻的生長代謝過程相關的淀粉與蔗糖代謝、碳代謝、糖酵解/葡萄糖酵素等表現出顯著差異。將DEGs 注釋到KOG 的25 個分類中，含有能量的產生和轉化，碳水化合物轉運與代謝和次生代謝產物的合成、轉運和代謝功能分類。根據DEGs 的GO 分類，KEGG 功能分析和KOG 功能注釋篩選出箭麻和共生天麻生長代謝（氮代謝、碳代謝和能量代謝）過程中的關鍵差異酶基因。qRT-PCR分析結果與轉錄組數據結果基本一致。在差異表達的基因中，篩選到與天麻氮、碳和能量代謝相關基因，合成氮和碳營養物質的差異基因表達上調，分解氮碳營養物質的差異基因表達下調，能量相關的差異基因大多數表達上調。與箭麻相比較，共生天麻有利于氮、碳營養物質和能量物質合成。

4.2 共生天麻和箭麻中生長代謝過程相關關鍵差異基因分析

氮、碳和能量的積累在植物生長過程中至關重要，因此選取了氮代謝，碳代謝和能量代謝過程中的關鍵酶基因。氮素同化在植物生長和發育過程中是一個十分重要的生理過程。植物從外界吸收的無機銨態氮和硝態氮必須先同化為谷氨酰胺和谷氨酸等有機物才能被植物利用。氨同化過程中谷氨酰胺合成酶是一個關鍵酶，它在ATP 供能條件下催化NH4＋同化為谷氨酰胺，然后在谷氨酸合酶作用下將其轉化生成谷氨酸。許多研究表明正常條件下高等植物氨同化的主要途徑是谷氨酰胺合成酶和谷氨酸合酶構成的循環反應途徑[19]。谷氨酸脫氫酶參與催化谷氨酸分解為α-酮戊二酸和銨離子及其逆反應合成谷氨酸的可逆反應。天冬酰胺和精氨酸是植物體內N/C 比值較高的氨基酸，是植物氮代謝過程中有機氮運輸和儲存的主要形式。近年來進一步研究表明，蘋果樹在越冬期間其體內可溶性和蛋白態氨基酸中，精氨酸含量均很高，這說明蘋果樹中的氮素在越冬期間主要以精氨酸形式貯藏[20-22]。qRT-PCR 分析表明共生天麻GDH和ASNS基因都明顯上調，分解谷氨酸的GDH基因和分解精氨酸的Argase基因都呈下調，這些基因不同表達水平有利于共生天麻積累便于在植物體內儲存運輸的谷氨酰胺，谷氨酸、天冬酰胺和精氨酸，為共生天麻和箭麻的基礎代謝和生長代謝提供氮素營養。

淀粉是高等植物在生長、發育和繁殖過程中重要碳水化合物營養物質的貯藏形式，并且以往研究表明蔗糖主要參與淀粉等的合成。植物蔗糖代謝過程的關鍵酶是蔗糖合酶，它催化蔗糖的合成和分解，但研究表明它主要催化蔗糖分解為淀粉合成所需的葡萄糖和果糖[23]。并有研究表明蔗糖合酶作為關鍵代謝酶在馬鈴薯塊莖淀粉合成中為其提供前體物質二磷酸尿核苷葡萄糖（uridine diphosphate glucose，UDPG）[24]，馬鈴薯中過量表達該酶可引起二磷酸腺苷葡萄糖（adenosine diphosphate glucose，ADPG）和UDPG 的含量增加，淀粉產量增加了近1 倍，同時馬鈴薯干質量較野生型也有所增加[25]。植物合成淀粉的另一個限速酶是AGPase，研究表明AGPase 是淀粉生物合成過程中的關鍵酶，也是淀粉合成的限速酶，抑制ADP-葡萄糖焦磷酸化酶活性將會部分或全部終止淀粉的合成[26-27]。可溶性淀粉合成酶也是淀粉合成過程中的關鍵酶。它以腺苷二磷酸葡萄糖為葡萄糖的供體，催化形成淀粉顆粒。淀粉酶是能將淀粉分子中的糖苷鍵水解獲得葡萄糖和寡糖等產物的一類酶。轉錄組分析和qRT-PCR 結果都顯示相較于箭麻，共生天麻中SuSy、AGPase和SSS基因均上調，AMS基因下調，這些基因的差異表達有利于共生天麻淀粉的合成和累積，供共生天麻和箭麻的生長發育需要。

糖酵解和檸檬酸循環是產生能量的主要途徑。丙酮酸激酶和己糖激酶是糖酵解途徑的限速酶[28]。一分子葡萄糖通過糖酵解途徑分解為2 分子丙酮酸可以產生2 分子ATP。蘋果酸脫氫酶和異檸檬酸脫氫酶是檸檬酸循環中主要限速酶[28]。一分子丙酮酸通過檸檬酸循環可以產生12.5 分子ATP。通過糖酵解和檸檬酸循環分解一分子葡萄糖可產生32 分子ATP。ATP 酶能催化ATP 水解為ADP 和磷酸根離子，并釋放能量的反應。基于轉錄組數據和 qRT-PCR 分析，共生區天麻HK、PDH、IDH、MDH基因都表達上調，說明共生區天麻糖酵解和檸檬酸循環都處于快速代謝產生較多的ATP 儲存能量。同時ATP基因也呈上調表達，它催化ATP基因分解釋放較多的能量以滿足共生天麻與蜜環菌共生所需要的能量和箭麻生長代謝所需要的能量。

5 結論

近年來從分子水平上對天麻及天麻與蜜環菌共生的研究頗多，最新研究已從生物學角度揭示天麻與蜜環菌共生關系及其形成機制[16]。但是天麻（白麻）與蜜環菌共生后進一步生長出新天麻（箭麻）的生長代謝特征還不清楚。本研究通過轉錄組比較分析初步闡述共生天麻和箭麻生長代謝特征。轉錄組分析表明，共生天麻中與氮、碳和能量代謝相關差異基因與箭麻相比較大多數合成氮、碳和能量物質表達上調，分解谷氨酸和精氨酸的谷氨酸脫氫酶和精氨酸酶以及分解淀粉的淀粉酶表達下調，有利于共生區天麻產生較多的氮和碳營養儲存物質和較多的能量物質供其進一步生長為箭麻。可以看出白麻與蜜環菌共生狀況是否良好對白麻能否進一步生長出箭麻有很大影響。在生產中科學地掌握白麻與蜜環菌共生時適宜的生長環境是白麻可以進一步生長為箭麻的關鍵，這是天麻栽培過程中需要注意的事項。然而，天麻蜜環菌共生時的最適宜環境仍需要進行下一步的探究。與此同時，箭麻和共生天麻與碳代謝、氮代謝及能量代謝相關的基因中，同一基因家族成員基因時空表達不一致，今后應研究清楚基因家族成員在箭麻和共生天麻中的亞細胞定位和時空表達特性，闡明其功能與作用，為人工栽培天麻提供進一步理論指導。

利益沖突所有作者均聲明不存在利益沖突