外源CaCl2缓解NaCl胁迫下酸枣幼苗叶和根转录组测序分析

【目的】研究CaCl₂缓解酸枣盐胁迫的分子机制。通过转录组测序分析,为研究CaCl₂缓解酸枣盐害分子机制提供参考。【方法】研究以水培酸枣幼苗为材料,150 mmol/L NaCl和150 mmol/L NaCl + 10 mmol/L CaCl₂分别处理6 h,利用高通量测序技术对酸枣幼苗叶片和根系进行了转录组测序与组装,对获得的所有Unigene在GO、KEGG数据库中进行功能注释。【结果】 与对照组相比,在外源CaCl₂处理下酸枣幼苗叶片上调和下调的差异表达基因分别7 365个和1 247个,根系分别有762个和4 861个。【结论】GO数据库比对,酸枣幼苗叶片富集了43个GO条件,根系富集了42个GO条件,叶片和根系均是催化活性被富集到差异基因最多;KEGG数据库比对,差异基因主要涉及植物病原互作、MAPK信号通路-植物、植物激素信号转导、苯丙醇生物合成、戊糖和葡萄糖醛酸相互转化等。     

盐胁迫下转 DcCIPK24 拟南芥的转录组比较分析

【 目 的】 探 究 铁 皮 石 斛（Dendrobium catenatum Lindl.）DcCIPK24 的 下 游 响 应 基 因， 为 研 究DcCIPK24 提高耐盐性的分子调控途径提供指导。【方法】利用 Illumina NovaSeq 6000 测序平台对 200 mmol/L NaCl 处理和对照组的转 DcCIPK24 拟南芥、野生型拟南芥进行测序，分析差异表达基因的富集通路，并对通路中的差异基因进行荧光定量 PCR 验证。【结果】盐胁迫下，从 2 个转基因拟南芥株系与野生型拟南芥的比较组合中共鉴定出 96 个差异表达基因；GO 注释分析发现 96 个差异基因主要被富集在分子功能相关通路；KEGG 富集分析发现差异基因主要被注释在氨基糖和核苷酸糖代谢、植物 MAPK 信号通路和甘油酯代谢通路中，选取上调表达的差异基因 AtGPAT5、AtSIRK、AtMEE25、AtUGE5 和下调表达基因 AtAMT1-2、AtATTI3、AtCYP71A25、AtLHCB4.3 进行实时荧光定量 PCR 验证，结果表明盐胁迫下 8 个差异基因表达趋势与转录组数据基本一致。【结论】盐胁迫下，DcCIPK24 主要通过氨基糖和核苷酸糖代谢途径、植物 MAPK 信号途径和甘油酯代谢通路响应耐盐。     

新疆野苹果幼苗叶片响应NaCl胁迫的miRNAs及靶基因筛选

【目的】 筛选新疆野苹果响应NaCl胁迫相关差异miRNA及靶基因。【方法】 以新疆野苹果组培苗为材料,对150 mmol/L NaCl处理6和48 h时的新疆野苹果幼苗叶片进行small RNA测序。【结果】 NaCl处理6 h时3个miRNA差异表达,其中2个上调表达,1个下调表达,miRNA对应的靶基因数目为64个;NaCl处理48 h时13个miRNA差异表达,其中4个上调表达,9个下调表达,miRNA对应的靶基因数目为108个。对差异靶基因进行GO和KEGG分类注释,NaCl处理6 h时,GO主要富集包括:肌醇磷酸2激酶活性、寡糖基转移酶活性、蛋白丝氨酸/苏氨酸磷酸酶抑制剂活性、木聚糖生物合成过程等;NaCl处理48 h时,GO主要富集包括:氧氧化还原酶活性、DNA结合、质外体、次生代谢过程等。KEGG共同富集在N-glycan生物合成通路。【结论】 mdm-miR168b、mdm-miR159c、mdm-miR827、mdm-miR390f、mdm-miR171i、mdm-miR399j基因在NaCl处理6和48 h时表达量存在差异,其中mdm-miR390f及其靶基因HF10976、mdm-miR171i及其靶基因HF33844、mdm-miR399j及其靶基因HF11095表达量呈负相关,3个miRNA参与了新疆野苹果对NaCl胁迫的响应过程。     

2种不同生长规格墨瑞鳕肌肉组织转录组测序分析

【目的】掌握不同生长规格墨瑞鳕个体间差异表达基因的表达特点,为其功能相关基因深度挖掘及分子遗传育种提供科学依据。【方法】挑选同一养殖条件下极大个体和极小个体的墨瑞鳕,构建肌肉组织cDNA文库后,采用Illumina HiSeqTM 4000测序平台对存在生长差异的墨瑞鳕肌肉组织进行转录组测序分析,获得的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等数据库中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行GO功能注释及KEGG信号通路富集分析,并采用MISA进行SSR鉴定分析。【结果】从墨瑞鳕肌肉组织中共测序获得39749条Unigenes,其长度范围在301~55230 bp,平均长度为1705 bp。注释到Nt、Nr、Swiss-Prot、Pfam数据库的Unigenes分别有27046、20824、18268和17772条,在7个数据库中均得到注释的Unigenes共计6742条,占Unigenes总数的16.96%。根据差异表达基因筛选条件P<0.05且|log₂ Fold Change|>1,共筛选出722个差异表达基因,其中上调基因308个、下调基因414个。差异表达基因GO功能注释分析结果表明,注释基因数目较多的GO功能条目包括细胞过程、代谢过程、膜、细胞器及结合等;KEGG信号通路富集分析发现,差异表达基因被成功富集到234条信号通路上,主要涉及磷脂酰肌醇3激酶/蛋白激酶信号通路、MAPK信号通路、胰岛素信号通路及FoxO信号通路等。在39749条Unigenes中鉴定筛选出22120个SSRs,占Unigene总数的55.65%,SSR的平均间距为3063 bp。【结论】基于转录组测序分析获得的墨瑞鳕肌肉组织差异表达基因以发挥结合、细胞过程及代谢过程等功能为主,且主要富集在PI3K-Akt信号通路、核糖体信号通路、FoxO信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。     

海马齿根系响应盐胁迫的转录组分析

【目的】对盐胁迫下海马齿根系进行转录组测序分析,挖掘海马齿根系耐盐相关基因,为揭示海马齿耐盐的分子机制提供参考。【方法】利用Illumina测序技术对0 mmol/L NaCl （对照组）和400 mmol/L NaCl胁迫处理（盐胁迫处理组）下的海马齿根系进行转录组测序分析,从中筛选出差异表达基因,选取13个基因进行实时荧光定量PCR （qRTPCR）检测,以验证转录组数据的可靠性。【结果】在海马齿根系转录组中共鉴定出305145个转录本,平均长度为622 bp,其中,对照组有146177个长度>300 bp的转录本,盐胁迫处理组有72173个长度>300 bp的转录本;共有65535条Unigenes在Nr、GO、Swiss-Prot、COG和KEGG五大数据库注释成功,占Unigenes总数的52.36%。对照组和盐胁迫处理组共有65535个差异Unigenes,其中,有182个热休克蛋白基因。对照组和盐胁迫处理组间共有24042个差异表达基因,从中选取13个基因进行qRT-PCR检测,结果显示,9个基因表达上调,其余4个基因表达下调,与转录组测序结果一致。24042个差异表达基因中,共有10106个显著差异基因富集到129条代谢通路,其中富集程度排名前10的代谢途径为核糖体、次级代谢生物合成、RNA转运、内吞作用、剪接体、甘油磷脂代谢、内质网加工、吞噬、醚脂类代谢和植物-病原体相互作用,参与盐胁迫相关的硫代谢、脯氨酸积累、活性氧（ROS）代谢、与盐胁迫相关的钙信号通路和过氧化氢代谢等途径的差异基因上调。【结论】在盐胁迫下海马齿差异表达基因如小分子量热激蛋白基因、抗氧化酶相关基因及与离子交换相关基因发挥了重要调控作用。     

育成期高能饲喂下开产与未开产蛋鸡肝脏miRNA表达谱差异分析

【目的】对育成期高能饲喂下开产与未开产蛋鸡的肝脏进行miRNA高通量测序分析,探明高能饮食状态下蛋鸡肝脏中影响其开产的miRNA,为提高产蛋性能打下基础。【方法】以代谢能水平为12.14 MJ/kg的高能饲粮饲喂育成期蛋鸡,通过Illumina NextSeq500高通量测序平台对开产与未开产蛋鸡肝脏进行small RNA测序,使用DESeq分析miRNA表达量,并以实时荧光定量PCR进行验证;采用miRanda对差异表达miRNA进行靶基因和靶位点预测,同时以超几何分布对差异表达的靶基因进行GO功能注释分析和KEGG信号通路富集分析。【结果】 6个样本（3个开产蛋鸡样本,3个未开产蛋鸡样本）能注释到的miRNA均超过300个,约占miRBase中已鉴定鸡miRNA的30.00%,且前体鉴定结果显示各样本注释到的miRNA有部分定位在同一前体上。与未开产组样本相比,开产组样本有9个miRNA表达上调、3个miRNA表达下调。其中,gga-miR-132c-5p、gga-miR-132b-5p、gga-miR-2184a-5p、gga-miR-132c-3p、gga-miR-132b-3p及gga-miR-2184a-3p属于同一miRNA基因簇,且gga-miR-1682、gga-miR-132b-3p和gga-miR-2184a-3p等3个上调miRNA在未开产蛋鸡样本中不表达。通过miRanda共预测得到129个差异表达潜在靶基因,以miR-203a预测得到的靶基因最多,为32个;其次是gga-miR-2184a-5p、gga-miR-148a-5p和gga-miR-12211-5p,分别预测到30、25和23个靶基因;而miR-132c-5p预测得到的靶基因最少,仅为9个。129个靶基因显著注释到8条GO功能条目上,生物学过程（Biological process）主要注释到脂质代谢过程、细胞脂质代谢过程、脂质生物合成过程、磷脂生物合成过程、磷脂代谢过程、甘油磷脂生物合成过程和解剖学结构发育,分子功能（Molecular funcion）仅注释到1条GO功能条目,即利钠肽受体活性,未发现涉及细胞组分（Cellular component）的GO功能条目;在注释到的8条GO功能条目中有6条与脂质代谢相关,涉及的靶基因有22个,占总潜在靶基因的17.10%。KEGG信号通路富集分析结果显示共显著富集到7条KEGG信号通路,其中脂肪酸降解、脂肪酸生物合成、酮体合成与降解及PPAR信号通路等4条信号通路与脂质代谢相关。【结论】育成期高能饲喂下开产与未开产蛋鸡中共存在12个差异表达miRNA,涉及129个差异表达潜在靶基因,且主要富集在肝脏脂质代谢相关过程和信号通路,说明miRNA是通过调控脂质代谢及其相关基因表达而影响蛋鸡开产。     

基于转录-代谢联合分析哈茨木霉ACCC32527 对NaCl胁迫的分子调节

【目的】 通过NaCl胁迫下哈茨木霉（Trichoderma harzianum）ACCC32527差异转录组和代谢组数据分析,挖掘关键功能型差异表达基因（DEG）及次级代谢产物。【方法】 采用RNA-seq和GC-TOF-MS技术分别完成0（T1）、0.4（T2）、0.6（T3）mol?L ^-1 NaCl胁迫下哈茨木霉ACCC32527的差异表达基因和次级代谢产物检测,利用相关软件及数据库（GO、COG、KEGG pathway等）完成对差异表达基因（|log₂fold change)|>1 & FDR<0.01）和次级代谢产物（P-value≤0.05 & VIP>1）的筛选、注释和分类,并进行RT-qPCR验证。 【结果】 NaCl胁迫处理后,ACCC32527的生长受到抑制,繁殖速率明显降低,并在该条件下分别获得T1 vs T2和T1 vs T3比较组637、1 570个DEG,获得DEG的16种表达模式,包括9种上调表达模式（950 gene）、3种下调表达模式（662 gene）和4种不规则表达模式（309 gene）,共有的GO功能分类为催化活性、结合、细胞器、细胞膜部分、细胞膜、细胞部分、细胞、单组织过程和代谢过程,且占主要比例,约为61%—94%,其中处理前后基因比例差异较大的分类为催化活性。KEGG代谢通路富集分析发现,不同NaCl浓度处理下的DEG富集到的代谢通路种类及数量存在较大差异,分别有225和535个差异基因富集于20条KEGG通路,包括代谢通路、抗生素的合成和次级代谢产物合成, 其中T2和T3处理下核糖体的生物合成途径分别有12个和18个相关基因受到抑制。从NaCl胁迫下差异转录组中筛选出活性氧清除、离子转运和细胞壁及胞外结构等共73个相关基因,并且多数为上调表达基因。另外,差异代谢组学分析表明,T3处理下,胞内小分子代谢物出现明显变化,累积量增加的代谢产物有30种,包括氨基酸及其衍生物、糖类及其衍生物和醇类,而减少的代谢物有53种,涉及脂肪酸和有机酸等。其中,甘油是已知的真菌重要渗透保护物,T3处理下ACCC32527胞内甘油水平提高,可参与细胞的渗透调节,维持渗透压平衡。RT-qPCR验证了7个差异表达显著基因,虽然在表达量上存在一定的差异,但表达趋势与RNA-seq分析结果基本一致,表明转录组测序结果的可靠性。【结论】 NaCl胁迫引起哈茨木霉ACCC32527大量基因及次级代谢产物变化,获得了与NaCl胁迫调节相关基因和代谢产物,这些机制共同作用减少盐离子对细胞的毒害作用,研究结果可为研究木霉菌的耐盐机理提供重要信息。     

卵形鲳鲹卵巢不同发育时期的转录组测序分析

【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的cDNA文库,采用Illumina HiSeqTM 2500进行转录组测序,经过滤、质量控制及拼接组装后获得的Unigenes在七大数据库（Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO）中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行功能注释及信号通路富集分析,并采用MISA和GATK3进行SSR鉴定及SNP分析。【结果】卵形鲳鲹卵巢组织转录组测序获得的325156432条Raw reads,经过滤筛选得到317206752条Clean reads,拼接组装后得到59554条Unigenes;69.65%的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等七大数据库中注释成功,其中有24599条Unigenes被注释到GO数据库,15997条Unigenes被注释到KEGG数据库。在卵形鲳鲹卵巢组织的2个发育时期共鉴定获得56115个基因,经差异表达分析后获得17737个差异基因,其中8169个基因在卵巢Ш期上调表达、9568个基因在卵巢Ш期下调表达。GO功能注释分析发现,卵形鲳鲹卵巢差异表达基因主要注释在细胞过程、氮化合物代谢过程、初级代谢过程、核、核部分、离子结合及水解酶活性等条目上;而KEGG信号通路富集分析结果显示,17737个差异表达基因显著富集在318条代谢途径上,其中前20条KEGG信号通路包括2-氧代羧酸代谢、PI3K-Akt信号通路、甲状腺激素信号通路、磷脂酶D信号通路、Fc εRI信号通路和细胞周期等。卵形鲳鲹卵巢转录组（59554条Unigenes）中共存在30133个SSRs和82490个SNPs。【结论】GnRHR、FSHR、FSHβ、CYP11A、SIRT3和PEG3等差异表达基因及PI3K-Akt信号通路和VEGF信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。     

高温胁迫下半夏叶片的转录组分析

【目的】通过分析高温（38℃）胁迫下半夏（Pinellia ternate）叶片转录组数据,挖掘其响应高温胁迫相关的基因,探究高温胁迫影响半夏的分子机制。【方法】通过Illumina HiSeq测序对高温胁迫处理后的半夏叶片进行转录组分析,并用qRT-PCR技术对DEGs进行验证。【结果】转录组测序共获得38.48 Gb Clean data,与对照组（25℃）相比,高温处理后1 138个DEGS上调表达,578个DEGs下调表达。GO富集分析表明,DEGs主要集中在转录因子复合体、内肽酶活性和血红素结合;KEGG富集分析表明,内质网蛋白质加工代谢通路的DEGs最多。基因表达水平分析表明,植物激素信号转导相关基因（PtARP1、PtGRP1、PtABF）、7个谷胱甘肽S转移酶基因PtGST和23个热激蛋白基因PtHSP在高温胁迫下表达量上调。qRT-PCR分析结果表明,9个DEGs在高温胁迫下的差异表达与转录组分析结果一致。【结论】通过对高温胁迫下半夏叶片进行转录组分析,挖掘高温胁迫相关基因,为进一步研究半夏高温胁迫响应机制提供理论参考。     

低氧—复氧胁迫下脊尾白虾鳃组织差异表达基因趋势分析

【目的】探究脊尾白虾（Exopalaemon carinicauda）鳃组织在渐变式低氧—复氧胁迫下的分子调控机制,为今后开展脊尾白虾耐低氧品系（种）的选育提供理论指导。【方法】通过模拟自然低氧环境的形成过程,分别于低氧处理0（对照）、3和6 h及复氧后1和8 h采集脊尾白虾鳃组织,利用Illumina HiSeqTM 4000测序平台进行转录组测序分析,经过滤和Trinity组装获得Unigenes,选取Nr、Swiss-Prot、KEGG和COG/KOG等数据库进行注释分析,在Omicsmart平台上完成差异表达基因筛选及其表达趋势分析,然后进行GO功能注释分析和KEGG信号通路富集分析,并随机选取5个差异表达基因进行实时荧光定量PCR验证。【结果】脊尾白虾鳃组织样本转录组测序数据经过滤后的Cleanreads进行Trinity组装共获得93227条Unigenes,其长度范围在201~35402 bp,平均长度为834 bp,N50长度为1352 bp。通过组间两两比较分析,共鉴定出4750个差异表达基因,其中上调差异表达基因3557个、下调差异表达基因2829个;超过50%的差异表达基因被显著富集到6种基因表达趋势模式中（P<0.01）,具体表现为:Profile 0模式富集到415个基因,Profile 5模式富集到201个基因,Profile 11模式富集到371个基因,Profile 13模式富集到841个基因,Profile 17模式富集到387个基因,Profile 18模式富集到411个基因。6种基因表达趋势模式中的差异表达基因被注释到代谢进程、细胞进程、单一有机体进程、细胞、细胞零件、大分子复合物及催化活性等GO功能条目上;而KEGG信号通路富集分析结果显示,以Profile 13模式中的差异表达基因富集到最多信号通路（86条）,其中呈显著富集的有8条,分别为核糖体、碳代谢、氧化磷酸化、氨基酸生物合成、内质网蛋白质加工、糖酵解/糖异生、谷胱甘肽代谢和蛋白输出。【结论】脊尾白虾鳃组织在受低氧胁迫早期通过合成蛋白质及提高代谢能力来抵御低氧环境,随着低氧胁迫时间的延长,物质合成和能量代谢活动均显著下降;但在复氧后随着复氧时间的延长,其蛋白质合成和能量代谢水平又逐渐升高恢复至常氧水平。     

青花菜花蕾转录组测序分析及蜡粉合成相关基因挖掘

【目的】对青花菜花蕾进行转录组测序分析，并挖掘与蜡粉合成相关基因，为探明青花菜花球表面蜡粉形成的分子机制提供理论参考。【方法】分别提取野生型和蜡粉缺失型青花菜花球总RNA，采用Illumina HiSeqTM 2500平台进行转录组测序，获得高质量Clean reads，采用Trinity进行序列组装后获得青花菜Unigene库，将获得的Unigene序列与Nr、Nt、KEGG、Pfam、KOG/COG、Swiss-Prot和GO数据库比对，获得基因功能注释信息；使用DESeq2进行差异表达分析。【结果】共获得44.68 Gb Clean data，De novo组装得到41244条Unigenes，N50长度为1847 bp。从所获得的Unigenes中筛选出8685个差异表达基因（DEGs）（上调基因5747个，下调基因2938个），共有8038个基因被注释到不同数据库，其中，5220个基因注释到Pfam数据库； 2066个基因注释到COG数据库，3866个基因注释到KOG数据库； 2580个差异表达基因被注释到75个转录因子家族中，注释最多的是MYB家族（235个）；GO数据库中6095个差异表达基因注释到细胞组分、分子功能和生物学过程三大类的52个功能分类； KEGG数据库中，1671个差异表达基因富集到138条代谢通路，其中13个差异表达基因与脂肪酸合成有关，7个差异表达基因与蜡粉生物合成途径有关。【结论】转录因子MYB家族在调控青花菜蜡粉合成中发挥重要作用。蜡粉合成过程中相关酶基因的差异表达是调控青花菜蜡粉合成的关键，尤其是野生型和蜡粉缺失突变体中特异性表达的差异表达基因，可作为后续研究青花菜花球表面蜡粉形成分子机制的对象。     

抗氟苯虫酰胺小菜蛾差异表达基因及其通路

【目的】随着氟苯虫酰胺等双酰胺类杀虫剂的广泛应用,小菜蛾（Plutella xylostella）对该类药剂的抗性愈发突出。研究旨在从转录组水平研究小菜蛾对氟苯虫酰胺的抗性机制,明确小菜蛾对氟苯虫酰胺产生抗性的关键基因及通路,为揭示小菜蛾对氟苯虫酰胺的抗性机制打下基础。【方法】以室内筛选的小菜蛾抗氟苯虫酰胺品系（Rh28）、田间抗性种群（Rz36）和敏感品系（S）为材料,通过转录组测序技术（RNA-Seq）,获得抗性小菜蛾品系/种群在转录水平上的差异表达基因及抗性显著上调表达基因,采用基因本体（GO）显著性富集分析及京都基因与基因组百科全书数据库（KEGG）富集分析,确定差异表达基因具有的主要生物学功能及其参与的主要生化代谢和信号转导途径。【结果】通过将不同抗性品系/种群与敏感品系的RNA-Seq数据对比,获得数量不等的差异表达基因。将差异基因富集的生物学过程（biological process）和信号通路（pathway）进行GO分析,发现其主要集中于对刺激物的反应（response to stimulus）、催化活性（catalytic activity）和代谢过程（metabolic process）等条目中;而KEGG分析发现,差异基因主要集中于代谢通路（metabolic pathway）中,基因数量占比最高,为17.84%。通过统计分析进一步获得了抗性显著上调表达基因,并对其中218个基因开展了GO分析,其仍然主要集中在代谢过程、对刺激物的反应和生物调节（biological regulation）等条目中;通过对抗性显著上调基因进行KEGG分析,发现其仍主要集中在代谢通路中,其中谷胱甘肽-S-转移酶、细胞色素P450解毒酶及参与昆虫多种生理反应的热激蛋白（Hsp）超基因家族等基因的表达量均显著上调。通过聚类热图分析发现,显著表达上调的基因较多集中在内质网上蛋白质的加工（protein processing in endoplasmic reticulum）、酪氨酸代谢（tyrosine metabolism）、咖啡因代谢（caffeine metabolism）和Wnt信号通路（Wnt signaling pathway）中。【结论】获得了小菜蛾抗氟苯虫酰胺差异表达基因及抗性显著上调表达基因,其主要富集于代谢过程、应激反应及对刺激物的反应等条目中,这些基因的相互协同调控是小菜蛾对氟苯虫酰胺产生抗性的重要机制。     

不同桑树品种响应干旱胁迫的比较转录组学分析

【目的】发掘响应干旱胁迫的关键抗旱基因,从转录水平上揭示桑树的抗旱分子机制,为后续开展桑树分子抗旱性育种工作提供科学依据。【方法】以干旱敏感性品种德果1号和耐旱性品种湖桑32号为研究对象,通过盆栽种植方式开展干旱胁迫及复水处理试验,采集24个桑叶样本提取总RNA后构建cDNA文库,在Illumina HiSeqTM 4000测序平台上进行高通量测序,并结合生物信息学对相关基因进行注释分析。【结果】经转录组测序,各样本的Clean reads范围为45723096~67280168条,有效碱基数（Clean bases）集中在6.86~10.09 Gb,GC含量在44.84%~46.48%（平均为45.61%）,Q30在90.36%~93.07%。差异表达基因（DEGs）筛选结果显示,6个差异分组（A2 vs A1,B2 vs B₁,A3 vsA1,B3 vs B₁,A4 vs A1,B4 vs B₁）分别筛选出3510、3399、5677、5507、5124和2734个差异表达基因;经干旱胁迫处理后,桑叶功能组基因中呈下调表达的差异表达基因明显多于呈上调表达的差异表达基因,说明桑树在生长过程中存在不同的功能基因以控制桑叶生长发育。不同抗旱性桑叶转录组测序数据中以涉及生物学过程的差异表达基因最多,且主要集中在小分子代谢过程（Small molecule metabolic process）、跨膜运输（Transmembrane transport）及碳水化合物代谢过程（Carbohydrate metabolic process）等方面;与分子功能相关的差异表达基因次之,主要涉及转移酶活性（Transferase activity）和水解酶（Hydrolase activity）等。干旱胁迫下不同抗旱性桑叶差异表达基因主要富集到59条KEGG信号通路上,可划分为代谢、遗传信息处理、环境信息处理、细胞过程和生物系统五大信号通路;其中,抗旱性桑树品种通过提高能量代谢、碳水化合物代谢、增强光合作用及脂质代谢来更好地适应干旱胁迫。综合GO功能注释分析和KEGG信号通路富集分析,得到以下可能与干旱相关的基因: LOC21410404、LOC21404884、LOC21409623、LOC21401352、LOC21398977、LOC21388561、LOC21401447、LOC21398764、LOC21385254、LOC21384661、LOC-21410404及LOC21408971。【结论】不同桑树品种的耐旱性存在显著差异,其中抗旱性桑树品种是通过提高能量代谢、碳水化合物代谢、脂质代谢及光合作用共同应对干旱胁迫。     

苹果叶片不定芽再生过程的差异表达基因鉴定与分析

【目的】筛选分析‘GL-3’苹果叶片不定芽再生过程中的差异表达基因（differentially expressed gene,DEG）,进一步解析苹果叶片不定芽再生的潜在分子机制,为提高苹果的遗传转化效率提供理论参考。【方法】‘GL-3’苹果继代组培苗叶片外植体接种在再生培养基上,分别于0、3、7、14和21 d后取样并提取RNA,构建mRNA文库后采用Illumina Nova seq平台进行测序。筛选出各时间点的DEGs,根据GO（Gene ontology）和KEGG（Kyoto encyclopedia of genes and genomes）注释结果以及官方分类,使用R软件中的phyper函数对筛选到的DEGs进行GO和KEGG富集分析;利用BLAST软件进行基因比对注释;重点分析植物再生相关的激素、酶、转录因子、多胺等DEGs;采用qRT-PCR对DEGs进行定量验证。【结果】再生培养基上培养3、7、14和21 d的苹果叶片外植体与对照组相比,分别筛选到5 250、4 937、6 852、6 493个DEGs,4个时间点共有的DEGs有3 027个。DEGs的GO功能富集显示,4个时间点筛选到的共有DEGs中上调表达的DEGs主要与氧化还原过程、细胞外围、蛋白激酶活性和有机环化合物结合等功能有关;下调表达的DEGs主要与单细胞代谢过程、钙离子结合、光合膜和类囊体部分等功能有关。DEGs的KEGG通路富集分析显示,4个时间点筛选到的共有DEGs中上调表达的DEGs主要富集在磷酸戊糖途径、植物激素信号转导、植物-病原菌相互作用和内质网蛋白质加工等途径中;下调表达的DEGs主要富集在α-亚麻酸代谢、苯丙烷生物合成、碳代谢和光合作用等途径中。对与植物离体叶片再生相关的激素、酶、转录因子和多胺等相关DEGs的表达模式进行分析发现,这些DEGs大部分呈上调表达趋势。经qRT-PCR验证后,所检测基因的表达趋势与转录组测序结果一致。【结论】通过对苹果叶片不定芽再生过程中不同时间点的基因表达谱进行检测和对比分析,获得了大量与苹果叶片不定芽再生相关的基因,研究结果为深入探讨苹果离体叶片再生机理提供了理论依据。