为缓解不断增加的能源需求与有限化石燃料储备以及环境恶化之间的矛盾,寻找更加生态、可持续再生的替代能源已成为未来新能源的发展方向。利用能源植物生产生物柴油与生物酒精类燃料是目前研究的热点。其中,小桐子作为一种新兴的能源植物,因其种子含油量高可制备生物柴油、不与粮食作物争地及改良生态环境等方面的特性而备受关注。由于起源热带,冷害成为限制小桐子产量及种植面积进一步扩大的主要因素。因此,深入研究小桐子的抗冷性形成机制并利用基因工程手段对其进行改良,培育耐冷新品种,是目前亟待解决的种质问题。本研究中,我们通过Illumina Hiseq 2000高通量测序技术对小桐子低温锻炼条件下的转录组(Transcriptome)进行了测序,共获得了55,112,142条clean reads片段,包含4,960,092,780nt。经过从头拼接,共得到106,749条序列重叠群(Contig),平均长度为338bp;进一步组装得到45,251条Unigene序列,平均长度为789bp。通过Blast序列比对与ESTscan软件预测,发现有34,274条Unigenes(75.74%)包含完整或部分的蛋白质编码区(33,362条通过Blast比对得到,912条通过ESTscan软件预测得到)。有35,791条Unigenes注释到了如Nr、Swiss-Prot、GO、COG、KEGG等不同蛋白质数据库。其中,27,293条Unigenes注释到GO数据库中的61个功能组中;11,887条Unigenes注释到了COG数据库中25个功能蛋白家族中;18,787条Unigenes注释到了KEGG数据库中的128条主要的代谢途径中。另外,45,251条Unigenes中共鉴定到6665个SSR位点,分布于5922条Unigenes中,平均长度12bp,平均发生频率为13.09%,SSR种类主要以单、二及三核苷酸重复为主,共占94.49%。基于转录组测序数据,我们对小桐子12℃低温锻炼12h、24h、48h的叶片材料进行了数字基因表达谱(Digital Gene Expression Profiling,DGE)的测序分析。与对照相比,小桐子经过12℃低温锻炼共发现4185个差异表达基因,在锻炼12h、24h、48h分别有1915、2085、3213个基因表达表现出差异性,在三个锻炼时间点都表现出差异表达的基因有553个。进一步分析表明,经过12℃低温锻炼12h、24h、48h,分别有1211、1302、2516个基因表达上调,同时分别有704、783、697个基因表达下调。上调表达的基因主要涉及电子传递与氧化还原平衡、水解酶与小分子运输蛋白家族、DNA或RNA结合蛋白类、各类转录因子如AP2/ERF与ABC等。对于主要代谢途径分析表明,小桐子在低温锻炼条件下,淀粉的分解代谢加剧,积累的渗透调节可溶性糖可能主要以肌醇、肌醇半乳糖苷以及蔗糖为主。而不依赖aba的cbf冷信号转导途径中,低温锻炼条件下各层次信号转导膜蛋白、激酶、cbf及其它多种转录因子都上调表达,说明cbf信号系统在小桐子抗冷性形成中起着重要作用。同时,小桐子在低温锻炼12h、24h、48h下,新诱导表达基因分别为238、160、330个,其中在不同锻炼时间点都诱导表达的新基因有61个,主要编码转移酶家族、dna结合蛋白家族、受体蛋白家族及氧化还原酶家族等。结合小桐子抗冷相关转录组与数字基因表达谱的测序数据分析结果,本工作还揭示了22个小桐子抗冷相关基因在小桐子抗冷性形成过程中的表达信息,同时也揭示了它们的组织表达特异性。22个基因在根中基本上都有基础表达,但随着低温锻炼时间的延长,其表达量总体表现为先下调后上调的表达模式,一般在12℃低温锻炼24h时达到最大表达量,其中以nsltpa、ga20ox、admt、dxs、pip2、dhn2、lea5表现最为显著。在茎中,所筛选的22个基因有18个对低温锻炼应答不明显,与对照一样都表现为本底表达;gsts、apx1、lea5表达则表现出低温诱导性;而pdi基因表达量甚微。在叶中,22个基因表达差异性最大,大部分基因表现为随低温锻炼时间的延长逐渐上调表达,而且在12℃低温锻炼24h时达到最大表达量。另外,部分低温锻炼新表达基因在对照条件下没有本底表达,这也与数字基因表达谱分析结果一致,如dnaj、chs、lrr-rlk、admt、dxs、apx1等基因。为进一步理解小桐子抗冷性形成机制,我们从小桐子中克隆了低温锻炼下差异表达变化较大的gs3、nsltpa基因,低温新诱导表达的chs、ga20ox基因,抗氧化和维持大分子功能相关的pdi与msra基因,以及cbf基因家族的两个成员cbf1与cbf2基因。氨基酸序列分析表明,gs3包含由7个平行β-折叠及8个氨基酸残基组成的保守活性中心,而nsltpa包含由4对二硫键(由8个cys残基形成)维系的4段α-螺旋构成的基序结构,都属于多基因家族,且分化上存在跳跃现象。chs属于pks家族,存在典型的查耳酮合酶及对苯乙烯合酶活性位点、产物结合位点与丙二酰-coa结合位点;而ga20ox属于2-odd家族,存在双层β-桶结构及h-d-h基序;两者在进化上都较为保守,与其它植物的直系同源蛋白都存在至少70%以上的序列相似性。pdi属于硫氧还蛋白家族,而我们克隆的小桐子pdi蛋白属于典型的l型成员,具有-a-b-b′-a′-结构域排列方式与两个-cghc-催化核心,并具备-kdel-内质网驻留信号,锚定在内质网膜上,发挥二硫键氧化、还原及异构化的功能。msra属于msr家族a亚家族,与该家族中的msrb、frmsr在氨基酸一级序列、高级结构等方面没有同源性,具有a亚家族蛋白的典型-gcgwp-保守序列以及3个α-螺旋包裹的6个β-折叠片层组成的α/β卷状结构;克隆的CBF1基因氨基酸序列中存在一个典型AP2基序,负责结合COR基因的CRT/DRE顺式作用元件,但没有鉴定到保守的Val-Glu-Glu残基,而是仅发现了中间部位的Glu。对其启动子序列分析发现,CBF1基因启动子区域包含典型的TATA框、CAAT框,以及ICE结合元件、MYB结合元件、MYC结合元件,表明CBF1受多种激素等信号的诱导,并参与多种逆境胁迫的应答。同时,重组GS3、nsLTPA、CHS及PDI酵母在低温下的长势都明显好于对照酵母菌,暗示这些基因与小桐子抗冷性直接相关,对于小桐子抗冷性的基因工程改良具有较大应用潜力。总之,我们在本研究中通过高通量测序技术得到了小桐子低温锻炼条件下的转录组数据,通过数字基因表达谱分析揭示了小桐子抗冷相关基因的表达模式及参与的主要代谢途径,并且实验鉴定了22个抗冷相关基因在根、茎及叶中的表达模式。另外,基于以上结果综合筛选,克隆了8个关键基因的全长cDNA,并对其进行了生物信息学分析及功能鉴定。从意义上而言,本研究初步阐明了小桐子低温锻炼下抗冷性形成的分子机制,得到的Unigenes序列极大地丰富了目前NCBI中小桐子的核酸数据库,为后续抗冷相关基因的克隆并通过基因工程手段对小桐子进行改良提供了理论与实践基础。