摘要：【目的】基于轉(zhuǎn)錄組學解析基因超顯性效應在煙葉生物量雜種優(yōu)勢形成中的作用機制，為高產(chǎn)雜交種選育提供理論依據(jù)?！痉椒ā恳跃哂猩锪侩s種優(yōu)勢的雜交種（K326×GDH88）及其雙親為材料，通過測定煙草親本與雜交種的生物量，確定雜種優(yōu)勢形成的關(guān)鍵時期，并采用轉(zhuǎn)錄組測序技術(shù)分析煙草煙葉生物量雜種優(yōu)勢相關(guān)差異表達基因（DEGs）及其表達模式?！窘Y(jié)果】移栽后45～52 d雜交種的生物量出現(xiàn)正向雜種優(yōu)勢，移栽后45～52 d的中親優(yōu)勢變化最大的階段，且雜交種基因的表現(xiàn)模型從加性中親模型轉(zhuǎn)變?yōu)檎蝻@性模型，故移栽后45 d是解析雜交種煙葉生物量形成關(guān)鍵時期。對移栽后45 d的樣品進行轉(zhuǎn)錄組測序分析，結(jié)果發(fā)現(xiàn)雜交種與雙親之間存在2660個DEGs，超顯性表達的DEGs有1683個，占比63.27%，可能是煙葉生物量雜種優(yōu)勢形成的主要原因。在超顯性表達的DEGs中，有28個超顯性上調(diào)表達的DEGs富集在光合作用相關(guān)通路上；10個超顯性下調(diào)表達的DEGs富集在呼吸作用相關(guān)的蔗糖和淀粉代謝通路上，光合作用和呼吸作用中超顯性表達的DEGs可能是調(diào)控雜交種生物量雜種優(yōu)勢形成的關(guān)鍵候選基因。【結(jié)論】煙草雜交種煙葉生物量雜種優(yōu)勢形成的關(guān)鍵時期是移栽后45 d。煙葉生物量雜種優(yōu)勢的形成主要受DEGs的超顯性效應影響；光合作用和呼吸作用相關(guān)的DEGs超顯性上調(diào)表達和超顯性下調(diào)表達可能是煙葉生物量雜種優(yōu)勢的形成原因。

關(guān)鍵詞：煙草；煙葉生物量；雜種優(yōu)勢；轉(zhuǎn)錄組學；超顯性效應

中圖分類號：S572.035.1文獻標志碼：A文章編號：2095-1191（2024）10-2966-14

Mechanisms of overdominance gene in the formation of biomass heterosis in tobacco leaves based on transcriptomics analysis

TAN Yong-yan1，ZENG Shuai-bo2，LU An-bin1，TANG Hua-jiang1，YANG Hao1，PENG Li-sha1，SHEN Si1，LIU Ren-xiang1*

（1Tobacco College，Guizhou University/Key Laboratory of Tobacco quality Research of Guizhou，Guiyang，Guizhou550025，China；2Tongren Company，Guizhou Tobacco Company，Tongren，Guizhou 554300，China）

Abstract：【Objective】This study aimed to elucidate the mechanisms by which overdominance gene contributed to biomass heterosis in tobacco leaves through transcriptomic analysis，thereby providing theoretical reference for the breeding of high-yield hybrids.【Method】The study utilized a hybrid（K326×GDH88）exhibiting biomass heterosis and its parental lines as experimental materials.By measuring the biomass of the tobacco parents and the hybrid，the critical period for heterosis formation was identified.Transcriptome sequencing technology was employed to analyze differen-tially expressed genes（DEGs）associated with biomass heterosis in tobacco leaves and their expression patterns.【Result】Positive heterosis in biomass was observed in the hybrid between 45 and 52 d after transplantation.The mid-parent advan-tage was the highest between 45 and 52 d post-transplantation.During this period，the gene expression model of the hybrid shifted from an additive mid-parent model to a positive dominance model，indicating that 45 d post-transplantation was the critical period for analyzing biomass formation in hybrid tobacco leaves.Transcriptome sequencing of samples col-lected 45 d after transplantation revealed 2660 DEGs between the hybrid and its parents.Of these，1683 DEGs（63.27%）exhibited overdominant expression，which was likely the primary cause of biomass heterosis in tobacco leaves.Amongtheoverdominantly expressed DEGs，28 were up-regulated and enriched in photosynthesis-related pathways，while 10 were down-regulated and enriched in sucrose and starch metabolism pathway related to respiration.These DEGs with over-dominant expression in photosynthesis and respiration were potential key candidate genes regulating biomass heterosis in hybrid.【Conclusion】The critical period for the formation of biomass heterosis in tobacco hybrid leaves is 45 d post-transplantation.The formation of biomass heterosis is primarily influenced by the overdominance effect of DEGs.Theoverdominant up-regulation and down-regulation of related DEGs in photosynthesis and respiration may be the underlying reasons for biomass heterosis in tobacco leaves.

Key words：tobacco；tobacco leaf biomass；heterosis；transcriptomics；overdominance effect

Foundation items：National Natural Science Foundation of China（32060510）；Guizhou Tobacco Company Project（2022XM02）；Bijie Tobacco Company Project（2024520500240064）

0引言

【研究意義】煙草（Nicotiana tabacum L.）是我國重要的葉用經(jīng)濟作物，其葉片生物量的大小直接決定了煙農(nóng)的種植收益。雜種優(yōu)勢主要表現(xiàn)在雜交后代的長勢、適應性、抗逆性和產(chǎn)量等性狀優(yōu)于雙親，近年來烤煙雜種優(yōu)勢利用研究成效顯著（劉仁祥等，2016），已在煙草煙堿（王國琴，2015）、K+含量（羅雯，2022）、葉片數(shù)（梁婷等，2023）、根系生物量（皮凱，2023）等性狀中開展了雜種優(yōu)勢研究，但性狀不同，其遺傳機制也不同。遺傳學家對于雜種優(yōu)勢的解釋提出了許多假說，其中顯性假說、超顯性假說和上位性假說被廣泛認可和接受（Jones，1917；Stuber et al.，1992；曹立勇等，2003），但雜種優(yōu)勢的產(chǎn)生取決于物種、性狀和遺傳背景，單一的假說無法獨立地解釋雜種優(yōu)勢形成的遺傳機制（Li etal.，2008）。因此，對具有生物量雜種優(yōu)勢的雜交組合進行轉(zhuǎn)錄組測序，以期解析雜種優(yōu)勢的形成機理，對明確煙葉生物量雜種優(yōu)勢的形成機制和提高豐產(chǎn)雜交種選配效率具有重要的研究意義?！厩叭搜芯窟M展】利用雜種優(yōu)勢可有效提高作物產(chǎn)量和品質(zhì)，已在小麥（張愛民等，2002）、水稻（唐志明等，2021）、玉米（魏鋒等，2022）等作物上得到了廣泛的應用。煙草葉片生物量受許多因素影響，包括栽培措施（陸星星等，2012）、光照條件（魏明月等，2017）、煙株個體發(fā)育（劉國俠等，2018）和品種（悅飛雪等，2019）等因素，其中受品種因素影響最大（尹劍藤等，2011）。尹劍藤等（2011）對10個煙草品種的田間性狀及煙葉組織結(jié)構(gòu)進行比較，結(jié)果發(fā)現(xiàn)品種基因型差異是導致其生物量出現(xiàn)差異的主要原因。陸星星等（2012）對4個烤煙品種在耕地條件、土壤肥力和前作等不同耕作措施下的產(chǎn)量進行差異分析，結(jié)果發(fā)現(xiàn)同一品種烤煙在不同耕作措施下產(chǎn)量存在明顯差異。劉國俠等（2018）對41個煙草品種產(chǎn)量與農(nóng)藝性狀進行相關(guān)分析，結(jié)果發(fā)現(xiàn)產(chǎn)量與節(jié)距和葉數(shù)呈顯著正相關(guān)，而與莖圍則呈顯著負相關(guān)。悅飛雪等（2019）研究發(fā)現(xiàn)，在同一鎘、鉛脅迫水平下，15個不同基因型煙草的生物量之間存在明顯差異，且在這2種脅迫下煙草生物量變化趨勢基本一致。魏明月等（2017）對遮蔭、自然光、強光3個光照強度下的盆栽煙草進行生物量及光合特性分析，結(jié)果發(fā)現(xiàn)不同光照環(huán)境對煙草生物量積累有顯著影響，高光強條件對煙草生物量的積累有利，而弱光下則受到抑制。曾帥波（2023）研究發(fā)現(xiàn)，煙葉生物量遺傳力高，雜種優(yōu)勢明顯，可利用雜種優(yōu)勢進行煙葉生物量改良?！颈狙芯壳腥朦c】隨著基因測序技術(shù)的發(fā)展，轉(zhuǎn)錄組測序技術(shù)為解析雜種優(yōu)勢形成機制提供了可能（Nejat etal.，2018），且已在玉米（Wang et al.，2022）、油菜（Xiong et al.，2022）、水稻（Huang et al.，2023）等作物中開展相關(guān)研究，為雜種優(yōu)勢的形成機理解析提供了基礎(chǔ)。本課題組前期發(fā)現(xiàn)，煙草雜交種煙葉生物量具有雜種優(yōu)勢，但對煙葉生物量雜種優(yōu)勢形成的相關(guān)機制及基因表達模式尚不明確，也未見相關(guān)研究報道?！緮M解決的關(guān)鍵問題】以具有生物量雜種優(yōu)勢的雜交種（K326×GDH88）及其雙親為材料，通過測定煙草親本與雜交種的生物量，確定雜種優(yōu)勢形成的關(guān)鍵時期，并采用轉(zhuǎn)錄組測序技術(shù)分析煙草煙葉生物量雜種優(yōu)勢相關(guān)基因及其表達模式，以期解析煙葉生物量雜種優(yōu)勢形成機制，為豐產(chǎn)雜交種的選育提供理論依據(jù)。

1材料與方法

1.1試驗材料

以本課題組前期篩選獲得的具有生物量雜種優(yōu)勢的雜交種（K326×GDH88）及其雙親為供試材料，由貴州省煙草品質(zhì)研究重點實驗室提供。RNAprep Pure Plant Kit試劑盒、FastKing gDNA Dispelling RT SuperMix試劑盒和Talent qPCR PreMix（SYBR Green）試劑盒購自天根生化科技（北京）有限公司。主要儀器設備：CFX96 Real-Time System（Bio-Rad公司）、凝膠成像系統(tǒng)（Bio-Rad公司）、冷凍離心機（Thermo-Fisher Scientific公司）、DYCP-44N電泳儀（上海新諾儀器設備有限公司）和LI-6400便攜式光合測定儀（LICOR公司）。

1.2田間試驗設計

試驗于2022年在安順市西秀區(qū)楊武鄉(xiāng)貴州大學煙草科研基地進行，田間采用隨機區(qū)組設計，設3次重復，每小區(qū)種植3行（行距1.10 m、株距0.55 m），每行15株，共計45株。土壤肥力中等且一致，田間栽培管理按優(yōu)質(zhì)烤煙種植規(guī)程進行，地塊四周設置保護行（2～3行），每行第一株和最后一株不參與數(shù)據(jù)測定。

1.3測定項目及方法

1.3.1生物量測定生物量采用殺青烘干方式測定：將全株葉片于105℃下殺青30min，75℃烘干至恒重，測定其干重。

1.3.2光合性狀測定于上午9：00－11：00在各小區(qū)挑選3株生長整齊的健壯煙株，采集其從下往上數(shù)第7～9葉，利用LI-6400便攜式光合測定儀進行凈光合速率測定。

葉綠素含量測定采用無水乙醇浸提法（努爾凱麥爾·木拉提等，2021），用紫外分光光度計分別在470、649和665 nm下測定吸光度，計算葉綠素含量（葉綠素含量=葉綠素a含量+葉綠素b含量）。

1.4轉(zhuǎn)錄組測序

1.4.1樣品采集煙葉樣品于團棵期（移栽后38 d）開始采集，每隔7 d取樣1次，共取4次。各小區(qū)選取3株生長整齊的健壯煙株，選取從下往上數(shù)第7～9葉，利用打孔器分別取出一部分樣品混合，置于離心管中，用液氮處理后存入-80℃冰箱用于后續(xù)試驗分析。

1.4.2 RNA提取、文庫構(gòu)建及測序采用TRIzol法提取葉片組織中的RNA，并利用NanoDrop 2000超微量分光光度計進行純度和濃度檢測。質(zhì)檢合格的RNA樣品用于cDNA文庫構(gòu)建，并用Agilent 2100生物分析儀進行cDNA文庫質(zhì)檢，合格后使用Illu-mina NovaSeq 6000平臺進行測序。

1.4.3差異表達基因篩選通過HISAT2（Version 2.1.0）將Clean reads與煙草基因組進行比對，并對比對結(jié)果進行質(zhì)檢。再通過RSEM結(jié)合基因組比對結(jié)果和注釋文件，獲取每個樣本基因/轉(zhuǎn)錄本的Read Counts，并進行FPKM轉(zhuǎn)換，得到標準化的基因表達量。然后使用DESeq2（Version 1.10.1）進行差異表達分析，以Plt;0.05且Fold Change≥1.5或Fold Change≤0.67為條件篩選出差異表達基因（Differen-tially expressed genes，DEGs）。

1.4.4基因表達模式分類將DEGs的表達量矩陣數(shù)據(jù)導入STEM（http：//www.cs.cmu.edu/～jernst/stem/）中，歸一化處理后進行聚類分析（Ernst and Bar-Joseph，2006）。根據(jù)聚類分析結(jié)果，對基因的表達模式進行劃分。

1.4.5 GO功能注釋和KEGG信號通路富集分析將DEGs導入GO數(shù)據(jù)庫，將DEGs分成生物學過程（Biological process）、細胞組分（Cellular component）和分子功能（Molecular function）三大類，分析其在功能上的差異。利用KEGG數(shù)據(jù)庫解析DEGs參與的信號通路，并進行分類，篩選出顯著（Plt;0.05）富集的信號通路，深入解析DEGs的生物學功能。

1.5實時熒光定量PCR檢測

根據(jù)轉(zhuǎn)錄組數(shù)據(jù)選擇6個DEGs，獲得其編碼區(qū)（CDS）全長，利用Primer 5.0設計其擴增引物（表1），委托北京擎科生物科技有限公司合成。以煙草LR25（GenBank登錄號：LI890）為內(nèi)參基因，利用Talent qPCR PreMix（SYBR Green）試劑盒對6個DEGs進行實時熒光定量PCR檢測。反應體系20.0μL：2×qPCR PreMix 10.0μL，正、反向引物（10μmol/L）各0.6μL，65.9μmol/L cDNA模板2.0μL，ddH2O補足至20.0μL。擴增程序：95℃預變性3 min；95℃5 s，退火溫度（表1）10 s，75℃15 s，進行40個循環(huán)。采用2-ΔΔCt法計算DEGs的相對表達量（Livak and Schmittgen，2001）。

1.6雜交種基因表達模式劃分及雜種優(yōu)勢計算

雜交種基因的表達模式按照Stupar（2008）的方法進行劃分。用F1表示雜種一代的數(shù)值。當F1高于高親值10%以上，說明雜交種基因超顯性表達；當F1介于親本值的-10%～10%，說明雜交基因顯性表達；當F1介于雙親平均值的-10%～10%，說明雜交種基因加性表達；不符合以上情況的定義為其他。

根據(jù)陳澤輝（2009）的方法計算煙葉生物量雜種優(yōu)勢，包括中親優(yōu)勢（MPH）和超親優(yōu)勢（HPH），計算公式如下：

MPH（%）=（F1-MP）/MP×100

HPH（%）=（F1-HP）/HP×100

式中，MP表示雙親平均值；HP表示高親值。

1.7統(tǒng)計分析

采用Excel 2019計算數(shù)據(jù)平均值、標準差及煙葉生物量雜種優(yōu)勢。利用SPSS 26.0進行單因素方差分析和Duncan’s多重比較，并利用GraphPad Prism 9.4繪制柱形圖。

2結(jié)果與分析

2.1雜交種煙葉生物量形成關(guān)鍵時期分析

由表2和表3可知，移栽后38 d雜交種（K326×GDH88）的生物量小于雙親，中親優(yōu)勢為-5%，移栽后45 d雜交種的生物量介于雙親之間，中親優(yōu)勢為-1%，這2個時期中雜交種與雙親的生物量差異不顯著（Pgt;0.05，下同），表現(xiàn)為負向生物量雜種優(yōu)勢；移栽后52 d雜交種表現(xiàn)出正向生物量雜種優(yōu)勢，中親優(yōu)勢為18%，隨后中親優(yōu)勢下降，在移栽后59 d降至4%。由圖1可知，移栽后38、45和59 d雜交種基因的表達模式為加性（中親本）模型，移栽后52 d為正向顯性模型，表明移栽后45～52 d雜交種的生物量出現(xiàn)正向雜種優(yōu)勢，移栽后45 d的中親優(yōu)勢變化最大的階段，且雜交種基因的表達模式從加性（中親本）模型轉(zhuǎn)變?yōu)檎蝻@性模型，移栽后45 d是解析雜交種煙葉生物量形成關(guān)鍵時期。因此，選取移栽后45d所取樣本進行轉(zhuǎn)錄組測序，深入分析雜交種煙葉生物量雜種優(yōu)勢的形成機制。

2.2雜交種與親本間的轉(zhuǎn)錄組差異分析結(jié)果

2.2.1測序數(shù)據(jù)質(zhì)量評估采用Illumina NovaSeq 6000測序平臺對雜交種和親本材料進行轉(zhuǎn)錄組測序，結(jié)果如表3所示。各樣本Clean data均達10.24 Gb以上，Q20和Q30分別在97.45%和92.63%以上，GC含量為42.93%～43.55%。將各樣品質(zhì)控后的原始數(shù)據(jù)與煙草K326參考基因組（https：//solgenomics.net/organism/Nicotiana_tabacum/genome）進行序列比對，對比率為95.92%～96.82%。可見，測序質(zhì)量均已滿足測序要求，可用于后續(xù)數(shù)據(jù)分析。

2.2.2雜交種與親本間的DEGs分析為了解基因差異表達對生物量雜種優(yōu)勢的影響，對雜交種及雙親的基因表達情況進行比較，結(jié)果如圖2所示。雜交種與K326比對組（KG vs K）共有918個DEGs，其中589個上調(diào)表達，329個下調(diào)表達；雜交種與GDH88比對組（KG vs G）共有2123個DEGs，其中1412個上調(diào)表達，711個下調(diào)表達；雜交種與雙親（K326和GDH88）比對組（KG vs K and G）有868個DEGs，其中503個上調(diào)表達，365個下調(diào)表達；排除雜交種與雙親之間重復的DEGs后，共篩選到2660個DEGs，推測這些DEGs與雜交種生物量雜種優(yōu)勢的形成有關(guān)。

2.2.3 DEGs的表達模式分析結(jié)果為更深入分析雜交種與親本間的DEGs，按Rapp（2009）的方法將DEGs劃分成12種表達模式（P1～P12）。由圖3-A可知，P1和P2模式為加性表達，其余模式為非加性表達；P3～P6模式為顯性表達，其中P3～P4模式偏母本顯性表達，P5～P6模式偏父本顯性表達；P7～P12模式為超顯性表達，其中P7～P9模式為下調(diào)超顯性，P10～P12模式為上調(diào)超顯性。DEGs表達模式分析結(jié)果（圖3-B和表5）顯示，在雜交種中，加性表達的DEGs有346個，占DEGs總數(shù)的13.01%，非加性表達的DEGs有2314個，占DEGs總數(shù)的82.99%。在非加性表達的DEGs中，顯性表達的DEGs有631個，占DEGs總數(shù)23.73%；超顯性表達的DEGs有1683個，占DEGs總數(shù)的63.27%，其中超顯性下調(diào)表達的DEGs有612個，占23.01%，超顯性上調(diào)表達的DEGs有1071個，占40.26%，表明雜交種和雙親間DEGs的非加性效應可能是導致生物量雜種優(yōu)勢形成的主要原因，其中基因超顯性表達的影響最大。

2.2.4雜交種超顯性表達的DEGs GO功能注釋和KEGG信號通路富集分析結(jié)果

2.2.4.1 GO功能注釋超顯性表達的DEGs GO功能注釋結(jié)果如圖4所示。1683個超顯性表達的DEGs被分成三大類，分別為分子功能、細胞組分和生物學過程。進一步深入分析發(fā)現(xiàn)，這些超顯性表達的DEGs在分子功能中被注釋到催化活性（Cata-lytic activity）、結(jié)合（Binding）、轉(zhuǎn)錄調(diào)節(jié)活性（Tran-scription regulator activity）和轉(zhuǎn)錄調(diào)節(jié)因子活性（Tran-scription factor activity）等5個條目，其中，催化活性和結(jié)合的DEGs較多，分別為664和853個，其次是轉(zhuǎn)錄調(diào)節(jié)活性和轉(zhuǎn)錄調(diào)節(jié)因子活性，分別為111和80個；在細胞組分中主要注釋到膜部分（Membrane part）、細胞部分（Cell part）、細胞器（Organelle）和膜（Mem-brane）等7個條目，其中膜部分和細胞部分的DEGs較多，分別為584和646個，其次是細胞器和膜，分別為345和201個；在生物學過程中主要注釋到刺激反應（Response to stimulus）、生物調(diào)節(jié)（Biological regu-lation）、細胞過程（Cellular process）和代謝過程（Meta-bolic process）等8個條目，其中，較多的DEGs被注釋到細胞過程和代謝過程上，分別為545和599個，其次是生物調(diào)節(jié)過程和刺激反應過程，分別為303和198個。

2.2.4.2雜交種超顯性表達的DEGs KEGG信號通路富集分析為了解這些超顯性表達的DEGs功能和代謝通路，對這些超顯性表達的DEGs進行KEGG信號通路富集分析，結(jié)果顯示，1683個超顯性表達的DEGs被富集到99個信號通路上。對富集程度排名前20條信號通路進行分析，結(jié)果如圖5所示。前20條信號通路包括光合作用—天線蛋白（Photosynthesis-antenna proteins）、ABC轉(zhuǎn)運體（ABC transporters），MAPK信號通路（MAPK signaling pathway），植物病原體相互作用（Plant-pathogen interaction），角質(zhì)、軟堿和蠟的生物合成（Cutin，suberine and wax biosyn-thesis），植物激素信號轉(zhuǎn)導（Plant hormone signal trans-duction），倍半萜和三萜的生物合成（Sesquiterpenoid and triterpenoid biosynthesis），苯丙素的生物合成（Phenylpropanoid biosynthesis）和光合作用（Photo-synthesis）等。在排名前20條信號通路中，與光合直接相關(guān)的通路如光合作用—天線蛋白通路和光合作用通路分別排名第1和第10，說明雜交種光合相關(guān)基因的超顯性表達有助于煙葉生物量雜種優(yōu)勢的形成。

2.2.5超顯性表達的DEGs在煙葉生物量雜種優(yōu)勢形成中的作用為了進一步探究超顯性表達的DEGs對生物量雜種優(yōu)勢形成的影響，將超顯性表達的DEGs劃分為正向超顯性表達和負向超顯性表達，分別對其進行KEGG信號通路富集分析。由于光合作用和呼吸作用是同化物產(chǎn)生和消耗的主要途徑，因此，后續(xù)對這2條途徑及相關(guān)信號通路進行深入分析。

超顯性上調(diào)表達的DEGs KEGG信號通路富集結(jié)果如圖6所示。有28個超顯性上調(diào)表達的DEGs被富集在光合作用相關(guān)的通路上（表6）。其中，光合作用—天線蛋白通路上有16個編碼LHHCB1、LHCB2和LHCB3蛋白的DEGs超顯性上調(diào)表達；光合作用通路上有12個超顯性上調(diào)表達的DEGs，其中有1個超顯性上調(diào)表達的DEGs編碼Psb蛋白（位于光系統(tǒng)Ⅱ）;有6個超顯性上調(diào)表達的DEGs編碼Psa蛋白（位于光系統(tǒng)Ⅰ）中；有3個超顯性上調(diào)表達的DEGs編碼Fd（位于光合電子傳遞）；有2個超顯性上調(diào)表達的DEGs編碼F型ATP酶（F-type ATPase）亞基b。這些DEGs超顯性上調(diào)表達使雜交種擁有更強的光合能力，在光能吸收和傳遞上起到重要作用，可產(chǎn)生更多的同化物，提高了雜交種的煙葉生物量。

超顯性下調(diào)表達的DEGs KEGG信號通路富集結(jié)果如圖7-A所示。有10個超顯性下調(diào)表達的DEGs富集在蔗糖和淀粉代謝相關(guān)通路，其中，蔗糖分解過程中調(diào)控酶β-呋喃果糖苷酶和淀粉分解的α-淀粉酶DEGs超顯性下調(diào)表達，減少了雜交種蔗糖和淀粉的分解，降低了果糖和葡萄糖的含量（圖7-B）。蔗糖和淀粉是植物光合作用產(chǎn)物的重要存在形式，其分解產(chǎn)物是植物呼吸作用的重要底物，這些下調(diào)超顯性表達的基因使蔗糖和淀粉的分解減弱，通過減少呼吸作用的底物含量來減弱雜交種的呼吸作用，減少了雜交種同化物的消耗，使雜交種擁有更高的生物量（圖7-B）。篩選獲得調(diào)控煙草生物量雜種優(yōu)勢的關(guān)鍵候選基因Nitab4.5_0004782g0020和Nitab4.5_0003280g0040（表7）。

2.3雜交種光合性狀分析結(jié)果

為進一步解析光合性狀在生物量雜種優(yōu)勢形成中的重要作用，對雜交種及其親本光合特性進行對比分析，結(jié)果（表8和表9）發(fā)現(xiàn)雜交種凈光合速率介于雙親之間，為19.94μmol/（m2·s），與K326無顯著差異，顯著高于GDH88（Plt;0.05，下同），雜種優(yōu)勢值為8.27%，表現(xiàn)為正向優(yōu)勢；葉綠素含量以雜交種（K326×GDH88）最高，達2.73μg/cm2，與K326無顯著差異，顯著高于GDH88，雜種優(yōu)勢值為6.24%，表現(xiàn)為正向優(yōu)勢，表明超顯性表達的光合相關(guān)基因可能提升了雜交種的葉綠素含量和凈光合速率，從而形成了煙葉生物量雜種優(yōu)勢。

2.4實時熒光定量PCR驗證結(jié)果

為了驗證轉(zhuǎn)錄組測序結(jié)果的準確性，隨機選擇6個DEGs，以LR25基因為內(nèi)參基因進行實時熒光定量PCR檢測，結(jié)果（圖8）顯示，實時熒光定量PCR檢測數(shù)據(jù)和轉(zhuǎn)錄組測序數(shù)據(jù)雖有部分差異，但基本表達趨勢一致，說明轉(zhuǎn)錄組測序數(shù)據(jù)可靠。

3討論

雜種優(yōu)勢利用作為提高作物產(chǎn)量和品質(zhì)的重要途徑之一，通過對雜種優(yōu)勢的利用可實現(xiàn)許多作物產(chǎn)量的提升（張愛民等，2002；唐志明等，2021；魏鋒等，2022）。煙草作為我國重要的葉用經(jīng)濟作物，其葉片生物量的大小直接決定煙葉產(chǎn)量和煙農(nóng)種植收益。本研究以具有煙葉生物量雜種優(yōu)勢的雜交種K326×GDH88及其親本為材料，在移栽后不同時期對各材料生物量進行動態(tài)測量，最終確定移栽后45 d為生物量雜種優(yōu)勢形成關(guān)鍵時期。雜種優(yōu)勢的形成原因較復雜，Huang等（2011）在水稻中驗證了基因的顯性效應是雜種優(yōu)勢產(chǎn)生的原因；Shang等（2016）在陸地棉產(chǎn)量性狀的雜種優(yōu)勢定位中發(fā)現(xiàn)，上位性是雜種優(yōu)勢產(chǎn)生的主要原因；Mo等（2022）從煙葉鉀離子雜種優(yōu)勢形成的研究中發(fā)現(xiàn)基因的超顯性效應發(fā)揮主要作用。本研究對雜種優(yōu)勢形成關(guān)鍵時期的樣本進行轉(zhuǎn)錄組測序分析，結(jié)果發(fā)現(xiàn)63.27%的DEGs表現(xiàn)為超顯性表達模式，說明煙葉生物量雜種優(yōu)勢主要受超顯性表達的基因影響，與Stuber等（1992）、孫其信等（1997）的結(jié)果一致。

對超顯性表達的DEGs進行KEGG信號通路富集分析，結(jié)果發(fā)現(xiàn)超顯性上調(diào)表達的DEGs主要富集在光合作用—天線蛋白通路上，編碼Lhcb1、Lhcb2和Lhcb3蛋白，這些蛋白收集光能并將光能傳遞到葉綠體的光合反應中心（Hasegawa et al.，2002），編碼這類蛋白的DEGs超顯性上調(diào)表達使雜交種具有更強的光合能力，與Joshi等（2023）沉默Lhcb1基因?qū)е鹿夂夏芰ο陆档慕Y(jié)論基本一致。同時超顯性上調(diào)表達的DEGs還富集在光合作用通路的光系統(tǒng)Ⅱ、光系統(tǒng)Ⅰ、光合電子傳遞鏈和F型ATP酶上，編碼Psb、Psa、Fd、F型ATP酶亞基b等物質(zhì)的作用，這些基因的超顯性上調(diào)表達使雜交種在電荷分離、電荷穩(wěn)定性、光保護、脅迫抗性（Sitaramametal.，2008；Paques et al.，2013；Li et al.，2021）和ATP合成（Liu et al.，2019）等方面具有更強的能力，保證了光合作用高效運行。光合相關(guān)通路上超顯性上調(diào)表達的DEGs使得雜交種具有更強的光合能力，可產(chǎn)生更多的同化物，為生物量雜種優(yōu)勢提供物質(zhì)基礎(chǔ)。

大多數(shù)植物光合產(chǎn)物的30%～60%會因為呼吸作用而被損耗（Cannell and Thornley，2000），故低呼吸速率與作物高產(chǎn)密切相關(guān)（Hauben etal.，2009）。在與呼吸作用相關(guān)的蔗糖與淀粉代謝通路中，調(diào)控蔗糖分解的β-呋喃果糖苷酶和淀粉分解的α-淀粉酶DEGs超顯性下調(diào)表達，使得雜交種蔗糖和淀粉的分解能力減弱。蔗糖和淀粉是植物光合作用產(chǎn)物的主要存儲方式，其分解產(chǎn)物果糖-6-磷酸（F6P）和葡萄糖是呼吸作用的重要底物（董文科等，2019）。雜交種蔗糖和淀粉的分解能力減弱導致雜交種同化物積累增加；同時，呼吸作用也因為底物的減少而減弱，與β-呋喃果糖苷酶（Fugate etal.，2024）和α-淀粉酶（譚彪等，2024）的活性和呼吸作用呈正相關(guān)性的結(jié)論一致。這些超顯性下調(diào)表達的DEGs使得雜交種同化物的分解減弱，呼吸消耗降低，減少了同化物的消耗，使雜交種積累更高的生物量。

4結(jié)論

煙草雜交種煙葉生物量雜種優(yōu)勢形成的關(guān)鍵時期是移栽后45 d。煙葉生物量雜種優(yōu)勢的形成主要受DEGs的超顯性效應影響；光合作用和呼吸作用相關(guān)的DEGs超顯性上調(diào)表達和超顯性下調(diào)表達可能是煙葉生物量雜種優(yōu)勢的形成原因。

參考文獻（References）：

曹立勇，占小登，莊杰云，鄭康樂，程式華.2003.水稻產(chǎn)量性狀的QTL定位與上位性分析［J］.中國農(nóng)業(yè)科學，36（11）：1241-1247.［Cao LY，Zhan X D，Zhuang J Y，Zheng K L，Cheng S H.2003.QTL mapping and epistasis analysis for yield components in a RIL population of rice（Oryza sativa L.subsp.Indica）［J］.Scientia Agricultura Sinica，36（11）：1241-1247.］doi：10.3321/j.issn：0578-1752.2003.11.002.

曾帥波.2023.雜交種煙葉生物量優(yōu)勢形成的生理及分子基礎(chǔ)研究［D］.貴陽：貴州大學.［Zeng S B.2023.Physi-ological and molecular basis research on the formation of biomass dominance in hybrid tobacco leaves［D］.Guiyang：Guizhou University.］

陳澤輝.2009.群體與數(shù)量遺傳學［M］.貴陽：貴州科技出版社.［Chen Z H.2009.Population and quantitative genetics［M］.Guiyang：Guizhou Science and Technology Press.］

董文科，馬祥，張玉娟，郭睿，張兆恒，張巖，馬暉玲.2019.低溫脅迫對不同早熟禾品種糖酵解代謝及其相關(guān)基因表達的影響［J］.草地學報，27（6）：1503-1510.［Dong W K，Ma X，Zhang Y J，Guo R，Zhang Z H，Zhang Y，Ma H L.2019.Effects of low temperature stress on glycolysis metabolism and related gene expression of different Poa pratensis varieties［J］.Acta Agrestia Sinica，27（6）：1503-1510.］doi：10.11733/j.ISSN.1007-0435.2019.06.006.

梁婷，闕遠慧，姜超英，賈傛吏，段麗麗，黃鶯，劉洋，劉仁祥.2023.烤煙品種韭菜坪2號葉片數(shù)的遺傳分析［J］.山地農(nóng)業(yè)生物學報，42（1）：70-74.［Liang T，Que Y H，Jiang C Y，Jia R L，Duan L L，Huang Y，Liu Y，Liu R X.2023.Genetic analysis of leaf number of flue-cured tobacco varie-ty Jiucaiping No.2［J］.Journal of Mountain Agriculture and Biology，42（1）：70-74.］doi：10.15958/j.cnki.sdnyswxb.2023.01.010.

劉國俠，冉法芬，李云，薛琳，李世金，朱啟法，季學軍，劉旦.2018.皖南地區(qū)不同烤煙種質(zhì)產(chǎn)量產(chǎn)值與農(nóng)藝性狀的相關(guān)分析［J］.科技通報，34（2）：46-50.［Liu G X，Ran F F，Li Y，Xue L，Li S J，Zhu Q F，Ji X J，Liu D.2018.Analysis of correlation between yield and production value and agronomic traits in different varieties of flue-cured tobacco［J］.Bulletin of Science and Technology，34（2）：46-50.］doi：10.13774/j.cnki.kjtb.2018.02.011.

劉仁祥，夏志林，崔慶偉，喻奇?zhèn)?，張慶明，聶瓊，徐秀紅，劉洋.2016.貴州烤煙育種工作進展與對策探討［J］.山地農(nóng)業(yè)生物學報，35（4）：1-6.［Liu R X，Xia Z L，Cui Q W，Yu Q W，Zhang Q M，Nie Q，Xu X H，Liu Y.2016.Prog-ress in tobacco breeding in Guizhou Province and discus-sion for solutions［J］.Journal of Mountain Agricultural and Biology，35（4）：1-6.］doi：10.15958/j.cnki.sdnyswxb.2016.04.001.

陸星星，龔蘭芳，李明福，董軍忠.2012.不同栽培條件與烤煙產(chǎn)質(zhì)量關(guān)系分析［J］.山地農(nóng)業(yè)生物學報，31（3）：231-234.［Lu X X，Gong L F，Li M F，Dong J Z.2012.Analy-sis of relationship between different cultivation conditions and yield and quality of flue-cured tobacco［J］.Journal of Moun-tain Agricultural and Biology，31（3）：231-234.］doi：10.15958/j.cnki.sdnyswxb.2012.03.017.

羅雯.2022.基于轉(zhuǎn)錄組學的煙葉鉀含量雜種優(yōu)勢形成機理研究［D］.貴陽：貴州大學.［Luo W.2022.Study on the formation mechanism of heterosis of potassium content in tobacco leaves based on transcriptome［D］.Guiyang：Gui-zhou University.］

努爾凱麥爾·木拉提，楊亞杰，帕爾哈提·阿布都克日木，瑪依努爾·吾斯曼.2021.小麥葉綠素含量測定方法比較［J］.江蘇農(nóng)業(yè)科學，49（9）：156-159.［Murat N，Yang Y J，Abudukerimu P，Osman M.2021.Comparative study on determination methods of chlorophyll content in wheat［J］.Jiangsu Agricultural Sciences，49（9）：156-159.］doi：10.15889/j.issn.1002-1302.2021.09.028.

皮凱.2023.煙草雜交種根系數(shù)量優(yōu)勢相關(guān)基因分析［D］.貴陽：貴州大學.［Pi K.2023.Analysis of genes related to heterosis of tobacco root number［D］.Guiyang：Guizhou University.］

孫其信，倪中福，陳希勇，劉志勇，黃鐵城.1997.冬小麥部分基因雜合性與雜種優(yōu)勢表達［J］.中國農(nóng)業(yè)大學學報，2（1）：64.［Sun Q X，Ni Z F，Chen X Y，Liu Z Y，Huang T C.1997.Heterozygosity and heterosis expression of some genes in winter wheat［J］.Journal of China Agricultural University，2（1）：64.］

譚彪，任慕瑤，楊正鵬，徐佳依，鄭華斌，唐啟源，王慰親.2024.GA3和KNO2引發(fā)處理對水稻芽期生長、淀粉及呼吸代謝的影響［J］.中國稻米，30（3）：40-47.［Tan B，Ren M Y，Yang Z P，Xu J Y，Zheng H B，Tang Q Y，Wang W Q.2024.Effects of GA 3 priming and KNO2 priming on the seedling establishment，starch metabolism and respiratory metabolism of rice［J］.China Rice，30（3）：40-47.］doi：10.3969/j.issn.1006-8082.2024.03.006.

唐志明，馬榮榮，王曉燕，陸永法，周華成，蔡克峰.2021.甬優(yōu)系列秈粳雜交組合產(chǎn)量優(yōu)勢與親本秈粳成分及農(nóng)藝性狀間的關(guān)系研究［J］.雜交水稻，36（5）：65-72.［Tang Z M，Ma R R，Wang X Y，Lu Y F，Zhou H C，Cai K F.2021.Relationship of yield heterosis of Yongyou series of indica-japonica hybrid combinations to the genetic distance and agronomical traits of parents［J］.Hybrid Rice，36（5）：65-72.］doi：10.16267/j.cnki.1005-3956.20200612.199.

王國琴.2015.煙草煙堿雜種優(yōu)勢表現(xiàn)及其基因差異表達分析［D］.貴陽：貴州大學.［Wang G Q.2015.Studies on tobacco nicotine trait of heterosis and related gene differen-tial expression analysis［D］.Guiyang：Guizhou University.］

魏鋒，史大坤，衛(wèi)曉軼，馬毅，洪德峰，馬俊峰，王稼苜，鄭秋道.2022.不同玉米雜交種產(chǎn)量和品質(zhì)性狀的雜種優(yōu)勢分析［J］.山東農(nóng)業(yè)科學，54（12）：25-30.［Wei F，Shi D K，Wei X Y，Ma Y，Hong D F，Ma J F，Wang J M，Zheng Q D.2022.Heterosis analysis of yield and quality charac-ters of different maize hybrids［J］.Shandong Agricultural Science，54（12）：25-30.］doi：10.14083/j.issn.1001-4942.2022.12.004.

魏明月，云菲，劉國順，宋亮.2017.不同光環(huán)境下煙草光合特性及同化產(chǎn)物的積累與分配機制［J］.應用生態(tài)學報，28（1）：159-168.［Wei M Y，Yun F，Liu G S，Song L.2017.Response of photosynthetic characteristics and accumula-tion and distribution of assimilation products in tobacco to different light environments［J］.Chinese Journal of Applied Ecology，28（1）：159-168.］doi：10.13287/j.1001-9332.201701.010.

尹劍藤，李佛琳，楊煥文，段樹蒼，陳茂建，李成杰，蔡寒玉.2011.10個烤煙新品種田間性狀及葉片組織結(jié)構(gòu)的比較［J］.中國農(nóng)學通報，27（10）：92-97.［Yin J T，Li F L，Yang H W，Duan S C，Chen M J，Li C J，Cai H Y.2011.Field property and leaf dissecting structure comparison in ten new flue-cured tobacco cultivars［J］.Chinese Agricul-tural Science Bulletin，27（10）：92-97.］

悅飛雪，李繼偉，王艷芳，劉杉，王小東，陳明燦，劉領(lǐng).2019.不同基因型烤煙對鎘、鉛富集特征評價［J］.江蘇農(nóng)業(yè)科學，47（8）：105-111.［Yue F X，Li J W，Wang Y F，Liu S，Wang X D，Chen M C，Liu L.2019.Evaluation of accumu-lation characteristics of cadmium and lead in tobacco with different genotypes［J］.Jiangsu Agricultural Science，47（8）：105-111.］doi：10.15889/j.issn.1002-1302.2019.08.024.

張愛民，劉冬成，聶秀玲，郭小麗，黃鐵城.2002.雜種小麥育種的戰(zhàn)略［J］.中國農(nóng)業(yè)科技導報，（5）：42-48.［Zhang A M，Liu D C，Nie X L，Guo X L，Huang T C.Strategy of hybrid wheat breeding［J］.Journal of Agricultural Science and Technology，（5）：42-48.］doi：10.3969/j.issn.1008-0864.2002.05.011.

Cannell M G R，Thornley J H M.2000.Modelling the compo-nents of plant respiration：Some guiding principles［J］.An-nals of Botany，85（1）：45-54.

Ernst J，Bar-Joseph Z.2006.STEM：A tool for the analysis of short time series gene expression data［J］.BMC Bioinfor-matics，7：191.doi：10.1186/1471-2105-7-191.

Fugate K K，Eide J D，Lafta A M，Tehseen M M，Chu C G，Khan M F R，F(xiàn)inger F L.2024.Transcriptomic and metabolomic changes in postharvest sugarbeetroots reveal widespread metabolic changes in storage and identify genes potentially responsible for respiratory sucrose loss［J］.Frontiers in Plant Science，15：1320705.doi：10.3389/fpls.2024.1320705.

Hasegawa K，Yukawa Y，Sugita M，Sugiura M.2002.Organi-zation and transcription of the gene family encoding chlo-rophyll a/b-binding proteins in Nicotiana sylvestris［J］.Gene，289（1-2）：161-168.doi：10.1016/S0378-1119（02）00539-5.

Hauben M，Haesendonckx B，Standaert E，Van Der Kelen K，Azmi A，Akpo H，Van Breusegem F，Guisez Y，Bots M，Lambert B，Laga B，De Block M.2009.Energy use effi-ciency is characterized by an epigenetic component that can be directed through artificial selection to increase yield［J］.Proceedings of the National Academy of Sciences ofthe United States of America，106（47）：20109-20114.doi：10.1073/pnas.0908755106.

Huang X H，Zhao Y，Wei X H，Li C Y，Wang A H，Zhao Q，Li W J，Guo Y L，Deng L W，Zhu C R，F(xiàn)an D L，Lu Y Q，Weng Q J，Liu KY，Zhou TY，JingY F，Si L Z，Dong G J，Huang T，Lu T T，F(xiàn)eng Q，Qian Q，Li J Y，Han B.2011.Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germ-plasm［J］.Nature Genetics，44（1）：32-39.doi：10.1038/ng.1018.

Huang Z B，Ye J，Zhai R R，Wu M M，Yu F M，Zhu G F，Wang Z F，Zhang X M，Ye S H.2023.Comparative transcrip-tome analysis of the heterosis of salt tolerance in inter-subspecific hybrid rice［J］.International Journal of Mole-cular Sciences，24（3）：2212.doi：10.3390/IJMS24032212.

Jones D F.1917.Dominance of linked factors as a means of accounting for heterosis［J］.Genetics，2（5）：466-479.doi：10.1093/genetics/2.5.466.

Joshi J，Amthor J S，McCarty D R，Messina C D，Wilson M A，Millar A H，Hanson A D.2023.Why cutting respiratory CO2 loss from crops is possible，practicable，and pruden-tial［J］.Modern Agriculture，1（1）：16-26.doi：10.1002/moda.1.

Li G Y，Chen T T，F(xiàn)eng B H，Peng S B，Tao L X，F(xiàn)u G F.2021.Respiration，rather than photosynthesis，determines rice yield loss under moderate high-temperature conditions［J］.Frontiers Plant Science，12：678653.doi：10.3389/fpls.2021.678653.

Li L Z，Lu K Y，Chen Z M，Mu T M，Hu Z L，Li X Q.2008.Dominance，overdominance and epistasis condition the he-terosis in two heterotic rice hybrids［J］.Genetics，180（3）：1725-1742.doi：10.1534/geneti cs.108.091942.

Liu Y M，Qu J T，Zhang L，Xu X Y，Wei G，Zhao Z F，Ren M Z，Cao M J.2019.Identification and characterization ofthe TCA cycle genes in maize［J］.BMC Plant Biology，19（1）：592.doi：10.1186/s 12870-019-2213-0.

Livak K J，Schmittgen T D.2001.Analysis of relative gene expression data using real-time quantitative PCR and the 2（-Delta DeltaC（T））［J］.Methods，25（4）：402-408.doi：10.1006/meth.2001.1262.

Mo Z J，Luo W，Pi K，Duan L L，Wang P S，Ke Y Z，Zeng S B，Jia R L，Liang T，Huang Y，Liu R X.2022.Comparative transcriptome analysis between inbred lines and hybrids provides molecular insights into K+content heterosis of tobacco（Nicotiana tabacumL.）［J］.Frontiers in Plant Scien-ce，13：940787.doi：10.3389/fpls.2022.940787.

Nejat N，Ramalingam A，Mantri N.2018.Advances in tran-scriptomics of plants［J］.Adv BiochemEngBiotechnol，164：161-185.doi：10.1007/10_2017_52.

Paques E L，García-Casas M D C，Charpentier J P.2013.Distri-bution of heartwood extractives in hybrid larches and in their related European and Japanese larch parents：Rela-tionship with wood colour parameters［J］.European Jour-nal of Forest Research，132（1）：61-69.doi：10.1007/s 10342-012-0654-1.

Rapp R A，Udall J A，Wendel J F.2009.Genomic expression dominance in allopolyploids［J］.BMC Biology，7（1）：18.doi：10.1186/1741-7007-7-18.

Shang L G，Liang Q Z，Wang Y M，Zhao Y P，Wang K B，Hua J P.2016.Epistasis together with partial dominance，over-dominance and QTL by environment interactions contri-bute to yield heterosis in upland cotton［J］.Theoretical and Applied Genetics，129（7）：1429-1446.doi：10.1007/s00122-016-2714-2.

Sitaramam V，Bhate R，Kamalraj P，Pachapurkar S.2008.Res-piration hastens maturation and lowers yield in rice［J］.Physiology and Molecular Biology of Plants，14（3）：253-271.doi：10.1007/s 1229 8-008-0024-0.

Stuber C W，Lincoln S E，Wolff D W，Helentjaris T，Lander E S.1992.Identification of genetic factors contributing to heterosis in a hybrid from two elite maize inbred lines using molecular markers［J］.Genetics，132（3）：823-839.doi：10.1093/genetics/132.3.823.

Stupar R M，Gardiner J M，Oldre A G，Haun W J，Chandler V L，Springer N M.2008.Gene expression analyses in maize inbreds and hybrids with varying levels of heterosis［J］.BMC Plant Biology，8（1）：33.doi：10.1186/1471-2229-8-33.

Wang L F，Li J，Lin Y，Dang K T，Wan J，Meng S J，Qiu X Q，Wang Q Y，Mu L Q，Ding D，Luo H B，Tang J H.2022.Comparative transcriptomics analysis at the key stage of maize ear development dissect heterosis［J］.The Plant Ge-nome，16（1）：e20293.doi：10.1002/tpg2.20293.

Xiong J，Hu K N，Shalby N，Zhuo C J，Wen J，Yi B，Shen J X，Ma C Z，F(xiàn)u T D，Tu J X.2022.Comparative transcrip-tomic analysis reveals the molecular mechanism under-lying seedling biomass heterosis in Brassica napus［J］.BMC Plant Biology，22（1）：283.doi：10.1186/S12870-022-03671-0.

（責任編輯陳燕）