烟草是我国重要的经济作物,栽培面积和产量均居世界首位,经济效益高,是我国财政收入的重要组成部分。干旱严重影响烟叶产量和品质,是烟草生产主要的非生物胁迫因素之一。烟草响应干旱胁迫是一个多水平的复杂调控反应,包括了生理生化水平的调节和基因表达网络的分子适应。干旱胁迫下的基因表达调控网络十分重要,尤其是转录因子可用来提高植物对干旱胁迫的耐受性。通过对模式植物的分子生物学和基因组研究,现已发现一些干旱胁迫响应基因和各种参与调节胁迫诱导基因的转录因子。干旱胁迫增加内源脱落酸(ABA)水平,诱导依赖ABA和不依赖ABA的转录调控网络的表达。对拟南芥的相关研究已初步建立了干旱响应转录调控网络,一些重要的启动子元件,如ABA反应元件(ABREs)和耦合元件(CE),也已通过实验验证,但迄今对烟草抗旱的分子机制尚不清楚。此外,microRNAs (miRNAs)也是重要的基因表达调控因子,其主要原理是通过抑制mRNA的表达实现转录后水平的调控。拟南芥中干旱胁迫响应的转录和转录后调控之间存在复杂的相互作用,有些miRNA(如miR159和miR169等)在干旱胁迫应答反应中发挥了重要作用,但这种作用在烟草中仍未被确定。本研究以烟草品种红花大金元为研究材料,经过生理生化指标的检测和筛选,分别对六叶期幼苗进行6h和48h的短期干旱胁迫处理,分析对照和各处理在干旱胁迫下根系mRNA和miRNA的表达动态、鉴定差异表达的miRNA所调控的靶基因并预测干旱胁迫下的基因表达调控网络,旨在揭示烟草苗期在干旱逆境条件下诱导根系mRNA和miRNA的快速响应机理,为阐述烟草苗期耐旱性的分子机制及利用分子标记辅助育种提供理论依据。本研究取得以下结果：1.为选择干旱胁迫下烟草根系基因表达分析的最佳时间点,用20%的PEG6000对烟草六叶期的幼苗进行模拟干旱处理,分别于Oh,3h,6h,12h,24h,48h和96h这7个时间点取其根系(5株烟苗根系均匀混合),其中0h的样品为CK,然后测定其超氧化物歧化酶(SOD)活性、丙二醛(MDA)和脯氨酸(Pro)的含量。结果表明,SOD活性、MDA和Pro含量显着增加,6h和48h两个时间点相对于3h和24h两个时间点,其变化更为显著(P<0.05)。结合胁迫后的表型和生理生化指标的结果,本研究干旱胁迫检测的最佳时间点是6h和48h。2.选择0h,6h和48h的根组织样品作为测序的试验材料,分别命名为NCK(对照组),N6H和N48H。分别对三组样品提取总RNA,采用高通量测序技术构建了数字表达谱(DGE)的文库。经测序,3个文库平均获得约3.37Mb个reads,3个文库中共有21128个基因表达,注释的烟草基因占43.7%,其中3个文库中基因均表达的有13101个基因,仅有1887个基因差异表达,占烟草根系总表达基因的8.9%。qRT-PCR验证结果和测序结果基本一致。利用k-means聚类算法,依照基因功能分类和基因分组对差异基因进行Gene Ontology (GO)功能显著性富集分析,共分成6个类群(P<0.05),主要包含一些编码脂肪酸代谢、酰基转移酶活性、氧化还原酶活性、乙醇代谢及初始乙醇代谢和转移酶活性等的基因；Pathway分析发现有17条通路可能被影响,主要涉及谷胱甘肽代谢、脂肪酸延长、二苯乙烯类化合物和姜辣素的生物合成、次生代谢物的生物合成和苯丙生物合成等。3.根据数据库预测烟草基因组中同源的转录因子,共确定了609个转录因子,DGE的结果表明：所有的转录因子都可以在干旱胁迫下的烟草幼苗根系中检测到。这些根系中的转录因子,有82个差异表达,分属于24个转录因子家族,主要是与抗旱性相关的MYB、NAC和ERF家族,而其他差异表达的转录因子家族(HD-ZIP,NF-YA,NAC, GRAS,TCP)主要参与了发育和分生组织保持,防御／应激信号通路(HSP,WRKY和bZIP),生长素(Aux/IAA)诱导的的激素介导或胁迫介导的信号传导等。53%转录因子在NCK的表达水平最高(G1),而只有10%的转录因子是在N6H的表达水平最高(G2),另外还有37%的转录因子表达高峰在N48H(G3)。此外,我们还确定了一些家族特异性转录因子的表达趋势。C2H2、MYB、WRKY、ERF和Dof家族的转录调节因子分别有3、3、4、3、13和3个基因在NCK中高效表达,GATA、MYB和ERF家族的转录因子在干旱胁迫过程的NCK和N6H时间点分别有1、3和3个基因的表达水平最高；MYB、NAC、MYB相关、NF-YA,HD-ZIP、和ERF家族分别有6、6、2、2、2和2个基因在N6H优先表达。4.通过与已知基因和Nicotiana benthamiana烟草的注释信息进行序列同源性比较,我们发现了276个干旱应答候选基因(DRGs),其中,约40%(110/276)的基因是WRKY、 NAC、ERF和bZIP家族的转录因子。同时,我们也分析了这些干旱应答候选基因(DRGs)的功能,发现了其中有46个是差异表达的干旱应答候选基因(DRGs)。在这些基因中,21个(46%)是分属于NAC (6),MYB (4),ERF (10)和bZIP (1)家族的转录因子；其他干旱应答候选基因(DRGs),如GRF6、ABF1、APX2, SIPK和ZPT2,在干旱胁迫时也呈现不同的表达模式。5.从干旱处理和对照的根组织样品以及正常发育的叶、茎组织中,分离所有小分子RNA,构建small RNA文库并进行高通量测序,测序经生物信息学分析,共检出122个烟草miRNAs。在全部的miRNAs中,保守的miRNAs比非保守的miRNAs的表达量要高得多,如miR166和miR168家族的的表达丰度占总测序片段的57%和16%。然而,仍有43个前人报道的miRNAs未在这些样品中检测到。相对于对照,干旱胁迫的样品中有5个已知的烟草miRNA家族表达量的相对变化(log2root-ck/root-treat)大于2,表明表达差异显著或极显著。虽然miR159, miR169, miR402和miR408的表达在其它植物受干旱胁迫的差异显著,但在本实验中根系样品处理前后表达并没发生明显的改变。6.为了解烟草响应干旱胁迫的miRNA的功能,以前人通过降解组测序得到的已鉴定的烟草miRNA的靶基因为参考,我们的研究仅有27个靶基因GSS序列可以比对到烟草参考基因组中,并与87个转录子相对应。只有两个干旱响应的miRNA家族(miR160和miR395)找到了对应的靶基因。7.在干旱和低温应急反应过程中,许多研究人员提出了响应过程中的复杂的调控网络,涉及到miRNA和其调控的靶基因。本研究在前人的基础上,结合DGE、small RNA的测序数据和其调控的靶基因也绘制了一个简易的调控网络。这个网络主要包含两条通路(依赖ABA和不依赖ABA的通路)。在依赖ABA的通路中,NCED参与了对干旱胁迫的快速和应急反应。其级联转录过程,包括了AREB/ABF,MYB,bZIP, NAC和CBF/DREB1,他们参与了在应激反应中的渐进和适应的过程。SnRK2.6蛋白激酶也参与了在ABA中的信号传导。在不依赖ABA通路中,未知蛋白被认为具有渗透传感器和ERF系统的上游元件的功能。此外,响应的miRNAs (miR160、miR162、 miR394miR395,和miR827)也都表明是其表达存在时间的特异性。这个调控网络的绘制和分析将为以后研究烟草应对干旱,低温和重金属等非生物胁迫反应的基因表达调控提供参考。