瓠瓜幼叶转录组功能基因测序及分析

为探讨瓠瓜转录组功能基因信息,对福州芋瓠瓠瓜叶片进行了Illumina Hiseqxten高通量无参转录组测序分析,共获得664 252 268个reads片段,经序列组装共计获得87 518个Unigene,总长度高达91 405 320 bp,76.03%的Unigene长度主要集中在1~1 000 bp。功能注释结果显示,有55 725个Unigene在Nr数据库获得注释,注释到葫芦科的Unigene数量最多,有26 354个,占总Unigene的47.29%。GO数据库注释到18 278个Unigene,可以分为分子功能、细胞组分和生物学过程共三大类的55个亚类。COG数据库中有41 635个Unigene获得注释,分布在信息存储与处理、细胞过程和信号传递、新陈代谢、无特征基因四大类的25个功能区域,其中参与氨基酸运输和代谢的Unigene数量最多。KEGG数据库中有24 770个Unigene获得注释,涉及到220个KEGG代谢途径。SSR位点检索结果显示,包含SSR序列的8 617条Unigene中存在11 029个SSR位点,发生频率为9.846%。本研究结果为进一步深入开展瓠瓜功能基因研究和SSR分子标记开发奠定了重要基础。     

大豆抗胞囊线虫4号小种转录组测序及分析

大豆胞囊线虫(soybean cyst nematode,SCN)(Heterodera glycines)是一种世界性的大豆(Glycine max)病害.本研究以抗SCN4号生理小种的黑豆为研究对象,对大豆胞囊线虫侵染后的两个发育时期进行转录组测序和生物信息学分析,宏观了解大豆的抗病机制.通过Illumina HiSeqTM2500测序共获得1.96× 1010 bp的数据.经过分析后,接种后第9天(第一时期)和接种后第17天(第二时期)分别得到2 180个和4 210个差异表达基因,同时,对差异表达基因进行基因本体(Gene Ontology,GO)注释,蛋白质直系同源簇(Clusters of Orthologous Groups of proteins,COG)注释,结果显示,GO注释和COG注释分别将差异表达基因分为了58和25个功能类别.另外,京都基因与基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,KEGG)富集分析结果表明,可将两个时期的差异表达基因分别定位到83和105个具体的代谢途径分支.其中,参与到激素信号转导途径,苯丙氨酸代谢,能量代谢等过程中的差异基因最多,这些数据为进一步筛选抗病关键基因及功能验证等研究提供了依据.     

百香果低温胁迫转录组及茉莉酸代谢基因分析

为研究百香果低温胁迫响应机制,以紫果百香果(Passiflora edulia Sims)为试验材料在0℃下低温胁迫处理,以常温处理为对照组(CK),采用Illumina HiSeq测序平台进行转录组测序,并对茉莉酸代谢相关基因进行挖掘。结果显示,共获得百香果转录组数据45.30 Gb,组装得到39 521条Unigene和5 311 个差异基因;GO分类显示注释的Unigene分为细胞组件、分子功能及生物过程三大类,其中差异基因数量最多为生物过程大类的代谢过程,包括甾醇生物合成、类黄酮糖脂化、酪氨酸代谢、L-苯丙氨酸生物合成、软木脂生物合成、芥子油苷代谢及长链脂肪-酰基辅酶A代谢等。KEGG途径富集分析结果显示,核糖体途径、淀粉与蔗糖代谢途径、植物激素信号转导途径及植物与病原体互作途径为百香果响应低温胁迫的重要代谢途径。实时定量PCR(qRT-PCR)分析表明,百香果低温胁迫后其茉莉酸代谢途径相关基因AOC、AOS、JAR1、MYC2、PYL和JAZ均上调表达,该结果与测序获得FPKM值变化趋势较为相似,说明测序结果较为准确。但低温胁迫后,COI1的表达水平呈下调趋势。研究发现,在百香果中茉莉酸对低温胁迫的响应机制大体上与模式植物一致,关于COI1和MYC2等基因的调控方式还有待进一步功能验证。本研究结果为进一步明确百香果抗寒机制提供了科学参考。     

青梗花椰菜和白梗花椰菜转录组分析

为探讨花椰菜花梗颜色的遗传基础,采用高通量测序技术(Illumina Hi Seq 4000)对青梗和白梗花椰菜进行转录组测序,共获得52 253 822个序列读取片段(reads),对测序的结果从头组装获得66 450个单基因(Unigene),N50为1 285 bp,平均长度为717.40 bp。转录物注释结果显示,45 390个基因有同源比对信息;21 060个基因无匹配序列信息,可能为花椰菜特异的基因序列。对所得差异基因(DEGS)进行不同数据库注释,GO注释到2 540个DEGS,分为细胞组分、分子功能及生物学过程3大类54个功能组;COG注释到的753个Unigene功能系统分为24类;以KEGG数据库为参考,将946个DEGS定位到119个代谢途径分支,注释到6个差异基因与类胡萝卜合成途径有关,9个DEGS与类黄酮生物合成途径有关,1个DEGS与叶绿素代谢相关,且与花椰菜的花梗颜色形成密切相关。16个差异基因通过qRT-PCR验证,结果与RNA-Seq测序结果一致。本研究结果进一步揭示花椰菜花梗颜色的形成机理,为花椰菜品种遗传改良奠定了一定的理论基础。     

甲基营养型芽胞杆菌WM7低温响应转录谱分析

为从分子水平揭示甲基营养型芽胞杆菌WM7低温响应机制,采用高通量测序技术对低温和常温培养的菌体进行转录组测序分析,对部分差异表达显著基因进行RT-qPCR验证。结果表明,共检测到2 029个差异表达基因,其中1 519个低温上调,510个低温下调。GO和KEGG分析显示,1 199个差异表达基因归属细胞组分、分子功能和生物过程3个部分,683个基因参与到225个代谢途径中,其中涉及基因最多的是磷酸戊糖途径、嘌呤和氨基酸合成途径。菌株低温响应主要与脂肪酸代谢、ABC转运体、转录调控、信号转导、碳氮代谢等途径相关。此外,检测到大量差异表达的未知功能基因,多个基因差异表达倍数均大于50倍。本研究从基因表达方面阐述了甲基营养型芽胞杆菌WM7低温响应机制,为全面揭示其低温生长机制奠定了理论基础,提供了基因信息。     

四倍体青花菜小孢子培养及胚胎发育途径研究

以同源四倍体青花菜为试材,进行了游离小孢子培养的胚胎诱导条件及胚胎发育途径的研究。结果表明:基因型是影响四倍体青花菜小孢子胚胎诱导率的重要因子,基因型A36的出胚率最高,达到28.5个·蕾-1;热激胁迫(32.5℃)诱导小孢子从配子体发育途径向胚胎发生途径转换,其小孢子胚胎发育途径为典型的单核期小孢子经对称分裂产生胚胎的途径(B途径);4℃低温预处理1d或2d结合32.5℃热激1 d后,对称分裂和多细胞团增加,供试基因型小孢子出胚率显著提高。本研究将有助于游离小孢子胚胎发生机理的研究,以及其他芸薹属小孢子胚胎体系的建立及优化。     

紫色黄秋葵转录组功能基因测序及分析

黄秋葵为药食兼用型植物,富含花青素、多糖、黄酮及萜类等活性成分。为探讨黄秋葵次级代谢过程中这些活性物质合成的遗传基础,以紫色黄秋葵嫩茎叶为材料,采用Illumina Hi Seq 2500高通量转录组测序技术获得23 026 500个短读序,经组装获得42 484个Unigene,平均长度877 bp。31 931个和21 926个Unigene分别在Nr和Swiss Prot数据库有同源比对信息;Unigene与COG、KOG数据库进行比对,根据其功能各分为24类和25类;通过Pfam数据库,所注释到的20 031个Unigene分属于7 277类蛋白结构域;GO注释到18 602个Unigene分为细胞组分、分子功能及生物学过程等三大类的51个功能组,其中大量Unigene与细胞部分、细胞、催化活性、结合活性、代谢进程及细胞进程相关;根据KEGG数据库,可将7 619个Unigene定位到119个代谢途径分支,其中有1、3、35、61、13和70个Unigene分属于花色素苷、黄酮、类黄酮、N-多糖、二萜类和萜类骨架生物合成路径,它们可能参与这些活性物质的生物合成。本研究结果为后续深入开展黄秋葵花青素、黄酮类等物质代谢合成途径及相关功能基因研究奠定了基础。     

北苍术对持续干旱的生理和转录组响应分析

为了探究北苍术响应干旱胁迫的生理和分子机制，本试验采用穴盘育苗法，分析持续自然干旱3～18 d北苍术响应干旱胁迫的生理特征；利用高通量转录组测序技术，分析干旱胁迫下北苍术的分子响应特征。结果表明，持续干旱6 d起，北苍术植株相对含水量极显著低于对照，可溶性糖、脯氨酸和丙二醛含量、抗氧化酶活性均显著或极显著高于对照；复水15 d时，北苍术过氧化氢酶（CAT）活性和植株相对含水量恢复到与对照无显著差异，但其他指标均仍极显著高于对照。主成分分析结果表明，可溶性糖和总蛋白含量、超氧化物歧化酶（SOD）和CAT活性可较好地反应北苍术的抗旱能力，干旱胁迫下北苍术隶属函数综合评价D值为0.765～0.865。转录组分析共获得14 538个差异表达基因（DEGs）,GO功能注释最多的Unigene与细胞组分相关，KEGG分类中注释最多的Unigene是糖代谢；注释到萜类与聚酮类代谢途径的Unigene共71条，其中12条编码倍半萜和三萜类生物合成途径中的7个关键酶。综合分析可知，干旱胁迫影响北苍术相关功能基因表达，进而影响其生理特征。本试验为北苍术的抗旱机理和主要活性成分生物合成途径研究提供了参考。     

春兰与大花蕙兰杂交种红花呈色机理初探

为探究兰属杂交种红花的呈色机理,选择绿花春兰和黄花大花蕙兰的红花杂交种黄金梅为试验材料,通过小花蕾期和始花期花瓣的转录组测序分析,筛选影响花色的关键结构基因及调节基因,实时荧光定量PCR(qPCR)验证基因的差异表达,分光光度计法检测花瓣色素(类黄酮、叶绿素A、叶绿素B、花青苷)的含量。结果表明,通过测序共获得113 780条单基因(Unigene),其中200~300 bp的Unigene占比42.68%;44 088 个Unigene获得注释信息;与小花蕾期花瓣比较,始花期花瓣上调基因有1 855个,下调基因有2 494 个;差异表达基因(DEGs)被GO数据库注释划分为47个功能组;1 219个DEGs被KEGG数据库注释,涉及166条代谢途径;根据KEGG代谢通路及基因注释结果,筛选出与花色相关的54个关键结构基因和21个转录因子;花青苷合成基因DFR、ANS、3,5GT及转录因子Zm1、Hv1、MYB305表达量在始花期上调, 在NCBI库BLAST同源基因发现黄金梅DFR、ANS、3GT基因与大花蕙兰的基因同源性最高,5,3GT、3,5GT、Zm1基因与蝴蝶兰和石斛兰的基因同源性最高,转录因子Hv1、MYB305与建兰的基因同源性最高。qPCR参试基因表达量变化与转录组测序结果趋势一致。各阶段花瓣总黄酮的含量显著高于类胡萝卜素和叶绿素含量;始花期花瓣中花青苷含量升高。本研究结果为兰属杂交种花色基因的功能分析提供了参考。     

利用建兰转录数据开发genic-SSR标记

建兰(Cymbidium ensifolium)因基因组信息缺乏导致其分子标记的开发和利用受到一定限制。为了开发可靠、有效的分子标记,本研究对建兰转录组进行了大规模测序,并对此数据进行简单重复序列(simple sequence repeat,SSR)搜索。根据基元、长度和所在位置选取80对genic-SSR进行引物设计;在12个建兰品种中进行扩增验证,结果发现,61对引物成功扩增,其中39对引物有多态性,其多态信息含量(polymorphism information content,PIC)范围为0.141~0.746,平均为0.348;对成功扩增的61个genic-SSR所在的独立基因(unigene)序列进行包括非冗余(non-redundant,NR)数据库、基因本体论(gene ontology,GO)、真核细胞的直系同源组(eukaryotic orthologous groups,KOGs)和代谢途径(kyoto encyclopedia of genes and genomes,KEGG)的注释,发现52个genic-SSR所在的unigenes比对上了NR的序列、38个unigenes归入了25个GO词汇、18个归入了9个KOG分类,15对genic-SSR涉及到了14条KEGG途径。以上结果表明,来自于建兰转录组的39对引物不但具有一定的多态性,而且通过注释赋予了其许多基因功能的信息,是图谱构建和基因定位的理想分子标记。由于这些引物具有一定的通用性,还可用于整个兰属遗传多样性和群体结构的分析。     

60Co-γ射线辐射水稻小RNA测序变异分析

为探讨γ射线辐射引起水稻损伤效应的小RNA变异,采用0、300和400 Gy的⁶⁰Co-γ射线辐射处理水稻高粳糯种子,利用Illumina HiSeq^TM 2500高通量测序平台对辐射后的高粳糯三叶期叶片进行小RNA测序。结果显示,0 Gy(对照)组共获得7 395 578个小RNA,254个已知miRNA成熟序列,26个新miRNA成熟序列;300 Gy剂量组(Gy3)共获得10 315 701个小RNA,265个已知miRNA成熟序列,29个新miRNA成熟序列;400 Gy剂量组(Gy4)共获得6 469 869个小RNA,261个已知miRNA成熟序列,29个新miRNA成熟序列。GO富集分析显示,Gy4 vs CK组差异miRNAs对应的候选基因富集到生物过程、细胞组分和分子功能三大类共58个小类,其中涉及分子功能的基因最多。Gy3 vs Gy4组差异miRNAs对应的候选基因富集到生物过程和分子功能两大类共7个小类,以离子结合、过渡金属离子结合和氧化还原酶活性占主要比例。KEGG富集分析显示前20条显著富集途径,Gy3 vs CK组差异miRNAs对应的候选基因富集到植物病原互作最高,富集最显著的是植物昼夜节律。Gy4 vs CK组差异miRNAs对应的候选基因组富集到代谢途径最高,富集最显著的是甘油磷脂代谢和醚脂类代谢,Gy3 vs Gy4组差异miRNAs对应的候选基因富集到甘油磷脂代谢、真核生物的核糖体发生和内吞作用三个途径最高,富集最显著的是醚脂类代谢。本研究从RNA水平为辐射诱变水稻引起当代的苗期损伤效应提供了新见解。     

利用EST序列鉴定葡萄应答外源赤霉素的基因

为了利用NCBI上大量的葡萄(Vitis viniferaL.)EST序列进行葡萄应答外源赤霉素基因的信息挖掘,通过本地化Blast、Blast2go等工具以及其他生物信息学工具和数据库,对NCBI上经过赤霉素(gibberellins,GA)诱导后的葡萄EST序列进行处理,获得了45条仅在赤霉素处理后表达的无冗余EST序列。Blastx注释结果显示,45条序列中有25条序列注释为假定蛋白或未知蛋白(约占55.6%),14条序列无注释信息(约占31.1%),6条序列(13.1%)注释有推定功能,其中包括整合酶蛋白(EE077049)、SPX结构域基因1(EE085000)、丝氨酸羧肽酶类44蛋白(EE091188)、CHY1(EE092187)、假尿苷合酶/转运蛋白(EE106096)、钾离子通道蛋白(EE108944)。Blast2go分析显示,共有29条序列通过基因本体论(GO)注释180个GO号,分属于分子功能(molecula rfunction)、细胞成分(cellular component)、生物过程(biological process)三大类的不同水平。初步推断,在外源赤霉素对葡萄进行处理后,应答基因主要具有结合活性(46.34%)和催化活性(39.02%)两方面的分子功能。其作用空间分布为整个细胞(44.68%)或者细胞器(40.43%)中,并且主要通过参与细胞生理过程(28.57%)和代谢过程(25.71%)等一系列复杂的生理生化反应完成整个应答过程。为进一步研究赤霉素处理后导致葡萄的基因表达情况提供了基本资料。     

干旱胁迫下野生大豆ABC转运蛋白转录组测序分析

为获取干旱胁迫状态下ABC转运蛋白家族相关基因,以抗旱性强的野生大豆为试验材料,浇灌20%PEG 6000模拟干旱胁迫,分别在干旱处理0、 6、12、24和48 h后取叶片提取RNA进行转录组测序;以GO、KEGG和NR三大数据库注释的结果为基础,再以ABC转运蛋白为关键词进行筛选。结果表明,共获得337条ABC转运蛋白相关序列,包括8大亚家族(ABCA、ABCB、ABCC、ABCD、ABCE、ABCF、ABCG和ABCI)。差异表达分析结果表明,共有182条表现差异表达。GO注释到的差异表达基因多数为嘌呤核苷酸分解过程,KEGG注释到的差异表达基因分布在5个亚家族(ABCA、ABCB、ABCC、ABCD和ABCG)。对NR数据库有注释的168条差异表达基因进行进化分析,结果表明ABCE、ABCF、ABCG家族与ABCI11为一类,ABCB、ABCC家族与ABCI10为一类,ABCA、ABCD家族与ABCI19为一类。本研究为进一步研究野生大豆干旱胁迫下ABC转运蛋白抗旱基因提供了一定的理论依据。     

番鸭攻击行为的转录组学研究

为从分子遗传学角度探讨番鸭(Cairna moschata)出现攻击行为的原因,本实验应用The Medial Recorder 2.0软件对选取的120只生长状况良好的番鸭进行为期1个月(每天24 h)的行为观察记录,并采用瞬时观察法和连续观察法筛选出表现攻击行为的番鸭(实验组)和未出现攻击行为且表现正常的番鸭(对照组)。对筛选出的对照和实验两组番鸭,本实验采用Illumina HiSeq 2000高通量转录组测序技术对实验和对照两组番鸭(每组3次重复)的下丘脑组织的差异表达基因进行筛选,利用基因本体(Gene Ontology,GO)和京都基因与基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,KEGG)数据库对筛选的差异表基因进行功能注释,并进行差异表达基因聚类分析。结果显示,在出现攻击行为番鸭的下丘脑组织中,被注释到GO数据库中的差异表达基因有626个,其中上调基因137个,下调基因489个,参与调控番鸭攻击行为的差异表达基因有26个,其中上调基因4个,下调基因22个,KEGG数据库注释结果显示,这些差异基因显著富集的信号通路为:吞噬体、溶酶体、沙门氏菌感染、MAPK信号通路、孕激素介导的卵母细胞成熟、mTOR信号通路等6个通路。结果还显示在监管行动的细胞骨架和神经活性配体-受体相互作用通路上差异表达基因富集不显著,但富集程度很大,并在该通路上筛选到了番鸭攻击行为差异表达相关基因。通过GO和KEGG数据库的注释,本实验初步筛选出了33个与番鸭攻击行为相关联的候选基因及其调控机制,并发现ERBB信号通路可能是参与番鸭攻击行为调节的关键通路。该实验结果为从分子遗传学角度探讨番鸭攻击行为的发生机制奠定基础。     

大豆孢囊线虫大豆和烟草群体的转录组比较分析

大豆孢囊线虫(Heterodera glycines,soybean cyst nematode,SCN)为植物专性内寄生线虫,主要危害大豆(Glycine max)等豆科植物,烟草(Nicotiana tabacum)是其非寄主。但近年来,本课题组在山东地区发现一个特殊的SCN群体(简称SCN_T),其可以寄生烟草,但对大豆的致病性很差。为揭示SCN_T与普通大豆孢囊线虫群体(SCN)对同一寄主致病性存在显著差异的分子机制,本研究采用Illumina平台Hi Seq~(TM)2500高通量测序技术对两个群体的2龄幼虫(second stage juvenile,J2)进行转录组测序和生物信息学分析,并采用q RT-PCR对转录组测序结果进行验证。结果表明,SCN和SCN_T共检测出1 628个差异表达基因(differentially expressed genes,DEGs),其中1 347个基因在SCN_T中上调表达,281个下调表达。对DEGs进行基因本体(Gene Ontology,GO)分析发现,被注释到分子功能的核糖体结构组成的DEGs数目最多,高达284个;其次是生物学进程的翻译过程,为211个;京都基因与基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,KEGG)富集分析表明,639个DEGs定位到277个代谢途径分支,其中参与到核糖体和氧化磷酸化途径的差异基因最多。另外,与线虫生长发育和寄生相关的代谢途径以及一些编码食道腺细胞分泌蛋白的DEGs同样被显著富集。q RT-PCR分析数据与转录组测序结果相关性较好(R~2=0.98),说明转录组测序数据可靠。本研究首次研究了不同致病力群体SCN和SCN_T在转录水平上的差异,为进一步解析SCN致病性差异形成的分子机制提供科学依据。     

应用Solexa测序技术分析棉花不同抗病品种的数字表达谱

棉花枯萎病是一种危害严重的世界性病害,构建棉花枯萎病不同抗病品种的基因表达谱将为进一步研究棉花抗枯萎病的分子调控机制及筛选抗枯萎病相关基因奠定基础.本研究选用高抗枯萎病的陆地棉品种中棉所12和高感枯萎病的陆地棉品种新陆早7号为材料,利用Solexa高通量测序技术分析棉花幼苗根部组织中的基因表达情况.Solexa高通量测序后构建了多达350万个reads的标签文库,根据已知序列信息,显著差异表达基因1 080个,其中上调基因302个,下调基因778个.基因功能分析表明,已注释的显著差异表达基因可划分为分子功能、生物过程和细胞组件3个功能注释本体,并进一步细分为37个功能类别.鉴定出112条代谢通路,注释了1 238个显著差异表达基因,共涉及氨基酸代谢、多糖的生物合成和代谢、脂类代谢、能量代谢及信号转导等方面.在植物激素信号转导通路分析中,涉及到26个已知功能基因,4个基因编码AP2/ERF转录因子,4个基因与LRR类受体丝氨酸/苏氨酸蛋白激酶或蛋白激酶有关,5个基因与蛋白磷酸酶2C有关,而这些基因都参与植物的抗病反应.随机抽取5个基因进行实时荧光定量PCR验证,其结果与测序结果一致.本研究初步揭示了棉花枯萎病抗感品种在转录组水平上的表达差异,获得了大量差异表达基因.     

大麦白化颖壳突变体的转录组学分析

为了探究大麦白化颖壳突变体形成机理,以白化颖壳突变体(AL)和野生型植株(WT)为试验材料,对2份材料抽穗期的颖壳进行转录组测序分析,通过Illumina Hi SeqTM2500高通量转录组测序技术分别从WT和AL中获得2.7 Gb和4.1 Gb有效数据,拼接组装后得到64 634条通用基因(unigenes),其中长度大于1 kb的有23 696条。通过对2份材料间差异基因的筛选,共得到672个差异表达基因(DEGs),其中139个上调表达,533个下调表达。GO功能富集分析发现,在分子功能类型中注释的差异基因主要与叶绿体相关;KEGG显著富集分析发现差异基因广泛涉及碳代谢、光合作用等调控途径。对其中6个相关基因进行了RT-qPCR验证,表达趋势与RNA-Seq测序结果一致。表明AL的形成可能是由于叶绿体形成和发育相关基因的表达受到抑制,进一步导致了光合作用的下降,该研究结果对深入探究大麦的白化形成机理具有重要的参考价值。     

鸭早期胚胎发育相关基因消减文库的构建与分析

动物胚胎发育受到多个基因的表达调控,为探索鸭早期胚胎发育的分子机制,本实验构建了鸭早期胚胎发育相关基因的消减文库,筛选与鸭早期胚胎发育相关的候选基因.采集两个不同时段的鸭(Anas platyrhynchos)胚(6和20 h)分别作为驱动组(driver)和实验组(tester),通过由双链特异性核酸酶(designatedduplexspecific nuclease,DSN)介导的消减杂交方法,构建鸭早期胚胎发育相关基因的消减文库,并进行生物信息学分析.结果表明,在构建好的文库中随机挑选216个阳性克隆进行测序,最终共得到54条有效序列,平均长度为1 084 bp.将这些序列在GenBank中进行BLAST比对,其中44条有同源基因,基因注释(geneontology,GO)显示,这些基因与代谢、建立细胞定位、生物调节等各种生物途径相关,它们在鸭早期胚胎发育过程中的不同方面(如神经系统发育、免疫调节、转录调节等)都发挥着非常重要的作用.这些差异表达基因的分离与鉴定,有助于家禽早期胚胎发育的相关研究.     

碱胁迫下榆叶梅差异基因表达分析

为了解榆叶梅(Prumus triloba Lindl.)在短期碱胁迫下的分子响应机制,本研究采用Illumina HiSeq^TM2500高通量测序技术对其在0、1、3、6、12 h碱胁迫下的叶片进行了de novo转录组测序,以此鉴定与碱胁迫相关的差异表达基因(DEGs)。结果表明,从5 960万条raw reads中获得了约5 300万条高质量clean reads,经从头组装后获得了124 786条榆叶梅unigenes。基于8 948个 unigenes差异表达(DEGs), GO富集分析发现,其中有28个基因可能在早期碱胁迫响应机制中起着重要作用; KEGG富集分析发现,不同碱处理时间下的KEGG途径(pathways)具有显著差异。对7个碱性胁迫反应相关基因进行RT-qPCR验证,所得基因的表达模式与RNA-Seq数据结果一致,证实了RNA-Seq结果的可靠性。本研究首次对碱胁迫下的榆叶梅叶片进行了转录组分析,获得的序列数据极大地丰富了榆叶梅可利用的基因资源,深入了解了榆叶梅碱胁迫的分子响应机制,同时为研究以榆叶梅为砧木的中国李(Prunus salicina)在碱胁迫下的适应机制奠定了分子研究基础。