烟草曲茎病毒（Tobacco curly shoot virus,TbCSV）是属于Geminiviridae科Begomovirus属的单组分病毒,伴随β卫星（Tobacco curly shoot betasatellite,TbCSB）,在我国云南和四川省广泛发生且危害严重,是引起番茄和烟草曲叶病的重要病原。我国单组份双生病毒与β卫星的进化关系及互作机制复杂,TbCSV被认为在进化中介于需要卫星的双生病毒和真正的单组分双生病毒之间,因此,开展对该病毒及其伴随卫星的研究对于解读双生病毒的起源和遗传进化机制具有重要意义。目前对TbCSV及TbCSB已开展的研究工作主要集中在病毒检测、致病性、基因功能、蛋白定位、种群变异分析、基因沉默载体构建等方面,该病毒在寄主植物体内与寄主互作的方式和作用细节,对寄主内源基因表达的影响,对寄主生物过程的影响,以及寄主植株对抵抗病毒入侵的响应这些问题尚不明确,弄清这些问题将有助于全面解析该病毒的致病机理。本论文拟利用转录组测序及生物信息学分析方法,分析TbCSV侵染对本氏烟内源基因表达的影响,以及这些差异表达基因的功能、作用方式和参与的生物途径。基于生物信息学分析结果,借助基因沉默及过表达技术,分析TbCSV及TbCSB侵染对本氏烟光合作用相关途径的影响,并明确本氏烟PR和LRR-RLK基因对TbCSV及TbCSB侵染的响应及本氏烟BR和JA途径对TbCSV及TbCSB侵染的响应。本论文的主要研究结果如下:1.病毒侵染对本氏烟内源基因表达的影响利用转录组测序技术,从TbCSV单独侵染（TbCSV或Y35A）以及TbCSV与TbCSB共侵染（TbCSV+TbCSB或Y35AB）的本氏烟中共鉴定出59814条基因,其中,3196条基因在Y35A vs CK中差异表达,占基因总数的5.3%,其中1391条上调表达,1805条下调表达;有4081条基因在Y35AB vs CK中差异表达,占基因总数的6.8%,其中1775条上调表达,2306条下调表达。在Y35AB vs CK中的差异表达基因（Differential expression genes,DEGs）共富集3074个GO term,上调基因和下调基因分别富集2215和2658个GO term,其中显著富集GO term 124个,上调表达和下调表达的基因各85和166个。Y35A vs CK中的DEGs共富集2966个GO term,上调基因和下调基因分别富集2041和2033个GO term,其中显著性富集GO term 76个,上调表达和下调表达的基因各78和66个。这些DEGs主要与生物过程、代谢过程、细胞过程、有机物代谢过程、初级代谢过程等功能相关,参与核糖体、植物激素信号转导、代谢途径、DNA复制、暗反应中的固碳作用等过程。随机筛选了12条基因用RT-qPCR进行验证,12条基因表达量的变化趋势与测序结果一致。进一步分析发现,TbCSV或TbCSV+TbCSB侵染后本氏烟光合作用相关及防御相关基因的表达均受到影响,为进一步明确因此进行了以下研究。2.病毒侵染对本氏烟光合作用的影响为明确病毒侵染对寄主光合作用的影响,在TbCSV或TbCSV+TbCSB侵染20d,提取并测定本氏烟叶绿素含量,并对其光合生理特性进行了测定。结果表明,对照植株叶绿素a的含量为579.7μg/g,TbCSV单独侵染和TbCSV+TbCSB共侵染植株叶绿素a的含量分别为458.7μg/g和475.0μg/g,显著低于对照植株。对照植株叶绿素b的含量为141.6μg/g,TbCSV单独侵染和TbCSV+TbCSB共侵染植株叶绿素b的含量分别为122.3μg/g和126.7μg/g,与对照植株相比有所下降,但差异不显著。本氏烟光合生理特性测定结果显示,与对照植株相比TbCSV或TbCSV+TbCSB侵染后本氏烟净光合速率（Pn）、气孔导度（Gs）、蒸腾速率（Tr）、气孔限制值（LS）、水分利用效率（WUE）、光能利用效率（LUE）和CO₂利用效率（CUE）显著降低,而胞间CO₂浓度（Ci）显著上升。叶绿素荧光特性测定结果显示,与对照植株相比TbCSV或TbCSV+TbCSB侵染后植株非光化学淬灭系数（qN）和激发能捕获效率（Fv’/Fm’）上升,电子传递速率（ETR）、光化学量子产量（ΦPSII）和光化学淬灭系数（qP）下降。叶绿体超微结构观察结果显示,TbCSV或TbCSV+TbCSB侵染后导致本氏烟叶绿体数量减少,TbCSV单独侵染植株叶绿体中有中等尺寸的淀粉粒积累,TbCSV+TbCSB共侵染接种植株叶绿体形态肿胀,弯曲变形,内部的片层结构因膨压被挤到一边。从测序数据中筛选了23条同时在TbCSV和TbCSV+TbCSB侵染后表达丰度较高且差异倍数较大的DEGs,其中15条基因与叶绿体相关,占筛选基因的65%。利用TRV或PVX介导23条目标基因在本氏烟中表达,结果显示NbGDCSa、NbGDCSb、NbrbcS、NbChlHa、NbChlHb、NbDXS、NbHemA、NbBchP、NbGRP和NbPU被介导表达后,本氏烟表型发生变化,主要表现为叶片的黄化、白化、褪绿、斑驳、脉明、萎蔫坏死、卷曲、突起及植株矮化等。这些结果表明,TbCSV或TbCSV+TbCSB的侵染降低了本氏烟光合作用的能力。3.本氏烟PR和LRR-RLK基因对病毒侵染的响应为了明确本氏烟PR和LRR-RLK基因对病毒侵染的响应,在前期分析显示差异表达的基因中随机筛选了11条PR基因和12条LRR-RLK基因,利用RT-qPCR检测其在TbCSV或TbCSV+TbCSB侵染4 dpi、8 dpi和12 dpi时的表达量。结果显示,NbPR2、NbPR3、NbPR11（PR基因）和Nb1g56140、Nb1g67720、Nb2g16250、NbEI-LRR（LRR-RLK基因）在TbCSV或TbCSV+TbCSB侵染后的3个时间点的表达量均上调;NbPR1、NbPR4、NbPR5、NbPR9、NbPR10a、NbPR17（PR基因）和Nb1g06840、Nb3g47570a、Nb4g26540、Nb4g30520、NbGSO1（LRR-RLK基因）在TbCSV或TbCSV+TbCSB侵染后的部分时间点表达量上调,而其余时间点与对照相比差异不显著;NbPR6、NbPR10b（PR基因）和Nb3g47570b、Nb4g36180、NbFbLRR-MAX2（LRR-RLK基因）在TbCSV或TbCSV+TbCSB侵染后有的时间点上调表达,有的时间点下调表达。随机选择NbPR2、NbPR6、NbPR9、NbPR10a和NbPR11,利用TRV分别介导在本氏烟中表达,检测TbCSV或TbCSV+TbCSB侵染4 dpi、8 dpi和12 dpi时TbCSV及TbCSB的积累量。结果显示,NbPR2、NbPR6、NbPR9、NbPR10a沉默后,在TbCSV或TbCSV+TbCSB侵染4 dpi时TbCSV的积累量均显著增加,此时TbCSB的积累量在NbPR9和NbPR10a的沉默植株中也显著增加。而NbPR11沉默后,在TbCSV或TbCSV+TbCSB侵染4 dpi和8 dpi时TbCSV的积累量均显著降低,4 dpi时TbCSB的积累量也显著降低。这些结果表明,本氏烟PR和LRR-RLK基因对TbCSV或TbCSV+TbCSB的侵染作出响应,部分PR基因对TbCSV及TbCSB的积累量具有一定的调控作用。4.本氏烟BR和JA途径对病毒侵染的响应为了明确BR和JA途径对病毒侵染的响应,筛选了6条BR途径相关基因和10条JA途径相关基因,用RT-qPCR检测其在TbCSV或TbCSV+TbCSB侵染4 dpi、8 dpi和12 dpi时的表达量。结果显示,NbBZR1、NbCYCD3-X2（BR途径）和NbDOX1、NbHPL1、NbKAT2、NbJAR1、NbJAZ1、NbJAZ2（JA途径）在TbCSV或TbCSV+TbCSB侵染后的部分时间点表达量上调,而其余时间点与对照相比差异不显著;NbCYP90C1/D1、NbBRI1、NbBZR2、NbCYCD3-X1（BR途径）和NbSAMT1、NbSAMT2、NbCOI1、NbMYC2（JA途径）在TbCSV或TbCSV+TbCSB侵染后有的时间点上调表达,有的时间点下调表达。本氏烟JA含量检测结果显示,TbCSV或TbCSV+TbCSB侵染12 dpi对照植株JA的含量为8.89 pg/g,TbCSV单独侵染植株为14.04 pg/g,TbCSV+TbCSB共侵染植株为12.33 pg/g,TbCSV或TbCSV+TbCSB侵染后JA的含量上升,但未达到显著水平。因BZR1和CYCD3基因参与调控植物的生长发育,为明确本氏烟NbBZR1和NbCYCD3-X1在生长发育和抗病防御中的作用,以PVX介导本氏烟NbBZR1和NbCYCD3-X1的表达,结果显示,NbBZR1过量表达20 d时本氏烟与对照植株相比株高明显变矮,节间距变短,顶部叶片出现向下皱缩、变形,上部部分叶片出现均匀分布的白色坏死斑。NbCYCD3-X1过量表达20 d时植株与对照相比可观察到新叶呈轻微皱缩。目标基因沉默或过表达后,TbCSV及TbCSB的积累量检测结果显示,无论NbBZR1沉默或过表达后,在TbCSV或TbCSV+TbCSB侵染4 dpi时TbCSV及TbCSB的积累量均显著降低。NbCYCD3-X1过表达后,在TbCSV或TbCSV+TbCSB侵染4 dpi时TbCSV的积累量显著降低,但TbCSB的积累量显著上升。JA途径受体基因NbCOI1沉默后,TbCSV的积累量在TbCSV单侵染8 dpi和12 dpi时显著降低,但在TbCSV+TbCSB共侵染3个时间点均显著上升,而TbCSB的积累量在4 dpi时显著降低。外源性MeJA和BL处理后,TbCSV及TbCSB的含量变化趋势在两种处理方法中一致,在TbCSV单独侵染4 dpi时TbCSV的含量显著上升,8 dpi和12 dpi时显著下降;而在TbCSV+TbCSB共侵染4 dpi和8 dpi时TbCSV及TbCSB的含量显著下降,12 dpi时TbCSV及TbCSB的含量显著上升。这些结果表明,本氏烟BR和JA途径参与了寄主抗TbCSV及其伴随卫星的侵染。5.BR和JA途径在响应病毒侵染过程中发生相互作用研究发现,BL处理本氏烟12 h后JA途径基因NbCOI1的表达量显著上升,在TbCSV单独侵染3个时间点的表达量均上升,在TbCSV+TbCSB共侵染4 dpi和12 dpi时也上升。MeJA处理6 h和12 h后NbBRI1的表达量均显著上升,在TbCSV单独侵染4 dpi和8 dpi时表达量上升,在TbCSV+TbCSB共侵染4 dpi和12 dpi时的表达量也上升。MeJA+BL共同处理后12 h NbCOI1和NbBRI1的表达量均显著上升,在TbCSV单独侵染的3个时间点,NbCOI1和NbBRI1的表达量均上升,在TbCSV+TbCSB共侵染4 dpi时,NbCOI1和NbBRI1的表达量也均上升。以上结果表明,MeJA处理后在TbCSV或TbCSV+TbCSB侵染前后,BR途径受体基因NbBRI1被诱导上调表达。同样,BL处理后在TbCSV或TbCSV+TbCSB侵染前后,JA途径受体基因NbCOI1也被诱导上调表达。同时,MeJA+BL共同处理后NbBRI1和NbCOI1均有被诱导上调表达。由此推测,MeJA和BL处理相互促进了彼此信号通路中受体基因的表达量上调,表明BR和JA途径在参与本氏烟对TbCSV或TbCSV+TbCSB侵染的调控过程中存在相互促进作用。本研究中,转录组分析表明TbCSV或TbCSV+TbCSB侵染后本氏烟内源基因的表达受到不同程度的影响,涉及本氏烟各项生理功能和生物过程。研究结果表明,TbCSV或TbCSV+TbCSB侵染后本氏烟光合作用能力降低。同时,本氏烟PR和LRR-RLK基因,以及BR和JA途径参与了本氏烟抗TbCSV及其伴随卫星的侵染,并参与调控病毒及卫星的积累量。这些研究结果可为解析TbCSV及其伴随卫星与寄主相互作用的分子机制,以及寄主抗病机制提供科学依据。