高粱响应丝黑穗病侵染的转录组分析

高粱是中国重要的粮食作物,用途非常广泛。丝黑穗病在中国的高粱产区均有发生,已成为影响中国高粱生产最重要的病害。本研究采用高通量测序技术对被丝黑穗病病菌侵染的抗感病高粱进行转录组测序,筛选得到差异基因620个,其中229个上调表达,391个下调表达。差异基因中有330个获得GO数据库功能注释,主要涉及生物合成过程、催化活性、细胞过程等诸多生理生化过程;KEGG数据库富集分析发现,共有105个基因被注释到52个代谢通路中,其中类黄酮生物合成、苯丙素生物合成、二苯乙烯、二庚烷和姜辣素的生物合成、植物激素信号传导、苯丙氨酸代谢、植物病原互作、类胡萝卜素生物合成、油菜素内酯生物合成等通路显著富集。有26个差异基因被鉴定为转录因子,分布在10个转录因子家族中。RT-qPCR验证了高粱中差异基因的表达量变化趋势与转录组测序分析结果一致。本研究结果表明,在抗病高粱响应病原菌侵染的过程中,大量基因被诱导或抑制表达,与病原胁迫响应有关的部分转录因子显著上调表达,抗病相关的代谢和信号传导途径被激活。本研究为进一步挖掘高粱抗病基因提供参考。     

基于iTRAQ技术的甘蔗受黑穗病菌侵染蛋白组分析

黑穗病已成为影响甘蔗产量和含糖量的重要病害。为从蛋白质水平探讨甘蔗应答黑穗病菌的分子机制,本实验选用抗黑穗病品种桂糖29号和感黑穗病品种崖城71-374,处理组用浸渍法接种黑穗病菌,对照组用无菌水模拟接菌,在接种180 d后采取甘蔗叶片,使用iTRAQ技术对蛋白质组进行研究。结果显示,桂糖29号中有定量信息蛋白1429个,差异表达蛋白290个,其中上调表达蛋白153个,下调表达蛋白137个;崖城71-374中有定量信息蛋白1576个,差异表达蛋白125个,其中上调表达蛋白55个,下调表达蛋白70个。抗病品种桂糖29号中差异表达蛋白数多于感病品种崖城71-374,且桂糖29号在KEGG富集到的代谢通路也更多,可能被侵染后抗病品种的免疫调节机制更为复杂,涉及的调控通路网更广。经对光合作用、抗氧化系统、钙信号、苯丙烷类代谢、激素相关差异表达蛋白及共有差异表达蛋白分析,发现光合作用通路、ROS、ABA、钙信号通路相关蛋白在2个品种中多为上调表达,且桂糖29号的上调表达蛋白数多于崖城71-374,可能参与甘蔗后期对黑穗病的应答。植物抗病是一个复杂的过程,需要多种功能与途径参与调控。在本实验中没有发现苯丙烷类代谢途径及一些酶(谷胱甘肽过氧化物酶、抗坏血酸过氧化物酶、超氧化物歧化酶、谷胱甘肽硫转移酶、过氧化氢酶)、激素(生长素、乙烯、赤霉素)参与甘蔗的抗病过程,可能与采样时间有关。     

2种黑籽南瓜响应枯萎病菌侵染的转录组学研究

为了筛选白皮黑籽南瓜和绿皮花纹黑籽南瓜抗枯萎病的差异表达基因,挖掘参与抗性相关的代谢通路。以尖孢镰刀菌侵染2 d后和未侵染的2种幼苗叶片作为材料,利用Illumina Hiseq 2500测序技术进行转录组测序分析。以未侵染的幼苗作为对照,白皮和绿皮花纹黑籽南瓜差异表达基因分别有2082和1116个。其中141个基因在2种黑籽南瓜中均呈现出差异表达,62个基因具有相同表达趋势。对差异表达基因进行GO功能注释,结果显示2种黑籽南瓜差异表达基因主要富集在生物学过程中。KEGG代谢通路分析发现白皮有22个差异表达基因显著富集在植物激素信号通路,绿皮花纹有30个差异表达基因富集在苯丙烷生物合成通路。本研究挖掘出黑籽南瓜抗枯萎病差异表达基因,预测了抗性通路,为瓜类抗病机制的深入研究奠定了理论基础。     

基于转录组测序筛选黄牛低氧适应性相关差异基因

为探明藏黄牛对低氧适应的相关差异表达基因并分析差异表达基因功能和代谢通路,以藏黄牛和三江黄牛心脏组织为研究对象,提取总RNA,建立测序文库后利用Illumina HiSeqTM 4000进行测序,经生物信息学分析筛选出差异表达转录本,并对差异表达基因进行GO富集和KEGG通路分析,最后采用荧光定量PCR(RT-qPCR)验证测序数据的准确性。结果表明:藏黄牛相比于三江黄牛得到显著差异表达基因608个,其中上调基因467个,下调基因141个,推测表达量较高的B2M、COX7C、NFATC1、ECE1、CAST基因是与低氧适应性相关的基因;差异转录本GO功能富集分析可知,在分子功能、细胞组分、生物过程分别富集到327,253,2 191个条目,其中富集程度最高的条目分别是细胞发育过程、结合与转录因子活性、细胞前缘。KEGG分析结果表明,差异转录本显著富集在8条通路,筛选低氧适应性候选基因NFATC1是MP-PKG信号通路与T细胞受体信号通路成员,推测其通过cGMP-PKG通路与T细胞受体信号通路参与低氧适应性调控。采用高通量测序技术对藏黄牛和三江黄牛心脏组织进行转录组测序分析,获得大量差异表达基因,为研究动物低氧适应性机制奠定基础。     

甜玉米乳熟期果皮厚度相关基因的加权基因共表达网络分析

为深入探究甜玉米果皮厚度发育的分子机制,利用华南农业大学甜玉米课题组选育的果皮厚度差异较大的甜玉米自交系M03和M08为试验材料,对其授粉后15,19,23 d的果皮进行转录组测序。联合差异表达基因分析与加权基因共表达网络(WGCNA)分析,鉴定甜玉米乳熟期果皮厚度相关的共表达基因模块,根据模块内基因的连接度鉴定核心基因,利用qRT-PCR验证核心基因表达量的可靠性。比较M03和M08不同时期果皮转录组数据共筛选出4 748个差异表达基因(DEGs),DEGs主要富集到碳水化合物代谢过程、植物-病原体相互作用、植物MAPK信号通路。WGCNA分析共构建了18个共表达模块,通过筛选得到4个具有生物学意义且与果皮厚度显著相关(相关系数的绝对值>0.60)的特异性共表达模块,分别是Turquoise、Yellow、Magenta和Pink模块。GO和KEGG功能富集分析结果均表明,特异性模块富集到的天冬氨酸、丙氨酸和谷氨酸代谢、半胱氨酸和蛋氨酸代谢以及植物MAPK信号通路在甜玉米果皮厚度发育中具有重要的作用。选取连接度最高的前20个基因进行核心基因筛选,通过基因功能注释最终筛选到包括M...     

基于转录组测序筛选鸡皮肤毛囊性状相关基因和关键通路

旨在通过对崇仁麻鸡背部皮肤毛囊组织的转录组学分析，筛选关键的差异表达基因和信号通路，探讨鸡皮肤毛囊性状的分子调控机制。采用TRIzol法提取皮肤毛囊组织总RNA,通过Illumina平台进行转录组高通量测序(RNA-Seq)。最后采用实时荧光定量的方法对转录组测序进行验证。研究结果如下：共筛选出了3 856个差异表达基因，在公鸡(BGB vs YGB)比较组合中筛选出了971个差异表达基因，其中274个表达上调，697个表达下调。在母鸡(BMB vs YMB)比较组合中筛选出了3 529个差异表达基因，其中1 477个上调，2 052个下调。在公鸡、母鸡共同的差异表达基因中筛选出了3个与皮肤毛囊性状有关的高表达基因KRT75、KRT6A和KRT14。在GO功能显著富集中，BGB vs YGB比较组合富集到了516个GO term, BMB vs YMB比较组合富集到了1 020个GO term,基因主要集中参与细胞黏附等活动。显著富集了多条KEGG通路，在BGB vs YGB比较组合中富集了8条显著的KEGG通路，在BMB vs YMB比较组合中富集了20条显著的KEGG通路，筛选出了...     

桃叶片低温响应转录组分析

抗寒是桃育种工作重要的目标性状,受多个基因的调控,转录组测序是对特定生理条件下的转录本进行测序分析,用以深入探究桃叶片低温响应的分子机制。以熊岳巨桃的叶片为试验材料对4℃低温(LT)处理下和正常叶片(RT)进行转录组测序分析。利用DESeq2进行差异表达分析,获得72个上调差异表达基因和3个下调表达基因。KEGG代谢通路分析差异表达基因的代谢通路,并进行蛋白的互作预测。结果显示,低温处理的叶片呈现卷曲、变黄和失水的生理状态,差异代谢基因主要富集在植物病原体互作、MAPK信号通路和次生物质代谢积累等途径。MYC2、SPCH和Pti4～6转录因子可能在响应冷胁迫中调节低温带来的伤害中起到重要作用,MYC2和CBF蛋白预测相互作用,且分别与多个蛋白互作。本研究表明,有多个代谢途径共同响应桃叶片的低温胁迫,还预测到了可能在冷反应调节途径中起到关键作用的相关转录因子,为探究桃响应冷胁迫的分子机制打下基础。     

烟草叶片响应镉胁迫的差异表达基因鉴定及分析

烟草具有超富集镉的能力,严重降低烟叶品质,影响其经济价值。为了阐释烟草响应镉胁迫的分子机制,本研究采集了镉浓度为0和500μmol L^-1培养条件下的烟草叶片进行转录组测序。共获得76.94 Gb有效数据(Clean data),Q30碱基百分比均达到95.43%以上;在镉胁迫的烟草叶片中,共筛选出7735个差异表达基因,其中4833个基因表达上调,2902个基因表达下调,并通过qRT-PCR分析验证了转录组数据的可靠性。对差异转录本进行GO和KEGG富集分析,GO注释表明差异基因涉及代谢过程、应激反应、细胞结构体、催化活性和转录调节活性等;KEGG富集分析表明上调差异基因主要富集在氨基酸的生物合成、碳代谢、氧化磷酸化和柠檬酸循环等通路,下调差异基因则主要富集在光合作用、次生代谢产物的生物合成、代谢途径和植物激素信号转导途径。进一步分析植物激素信号转导通路发现,共有8条植物激素途径以不同的表达方式参与烟草对镉胁迫的响应。激素喷施烟草的实验结果表明,叶片通过调控赤霉素、油菜素内酯和茉莉酸途径以应对镉胁迫;拟南芥激素信号缺失突变体验证实验表明赤霉素、油菜素内酯、茉莉酸...     

金花葵开花前后差异基因及代谢物分析

为深入研究金花葵花朵开花前后基因和代谢物的变化,利用转录组学和代谢组学技术相结合的方法对金花葵花蕾和花朵进行检测。结果显示,通过转录组分析鉴定了206 636个Unigenes,筛选出42 618个差异表达Unigenes,其中包括63个差异表达转录因子家族,24个转录调节因子家族。GO分析结果显示,差异表达基因主要富集于囊泡介导的逆行运输等生物学过程。KOG功能注释显示,差异表达基因功能以通用功能预测居多,其次是信号转导机制、翻译后修饰蛋白周转以及碳水化合物的转运和代谢。差异表达基因KEGG富集分析表明,差异基因主要富集于代谢、植物激素和信号转导、淀粉和蔗糖代谢等代谢通路。通过代谢组学检测,筛选到差异显著代谢物135个,主要包括脂类、氨基酸及衍生物和黄酮类等,KEGG富集分析表明,差异代谢物主要富集于甘油磷脂代谢、糖基磷脂酰肌醇(GPI)-锚定生物合成和次生代谢产物的生物合成-未分类等过程,其中,缬氨酸,亮氨酸和异亮氨酸的降解代谢通路同时出现在了基因和代谢物KEGG代谢通路富集Top 20中。     

猪毛菜不同部位转录组特征分析

本研究以猪毛菜的地上部分和地下部分为实验材料,通过转录组测序分析探讨猪毛菜不同组织部位的基因表达特征和代谢物质形成有关的基因调控。测序数据分析结果显示获得192 777条转录本序列。差异表达筛选得到14 424个差异基因,其中上调表达基因6 137个,下调表达基因8 287个;GO功能富集分析得到14 424个差异基因。差异基因与KEGG数据库进行比对,共有5 169个差异基因注释在138个通路上,其中参与次生代谢途径共11条通路呈显著性富集。进一步分析黄酮类化合物合成和异喹啉生物碱生物合成通路中的差异基因和其对应的相关酶。随机选取了黄酮类化合物合成通路上的8个基因,通过RT-qPCR分析验证了转录组数据的准确性。本研究利用高通量测序技术和生物信息分析获得猪毛菜不同部位的转录组信息特征,为猪毛菜的药用成分研究和基因功能鉴定提供了科学依据。     

干旱胁迫下枇杷叶片的转录组分析

分析干旱胁迫下枇杷叶片的转录组,挖掘功能基因并对其差异表达基因进行筛选和分析,为枇杷抗旱提供理论依据。利用新一代高通量测序技术测序,对测序结果进行de novo拼接、功能注释和ORF预测,将差异表达基因在COG、GO和KEGG数据库中进行比对注释。测序结果表明,获得转录本共88 530个,平均长度为740.64 bp,ORF41 748条。COG、GO和KEGG数据库将转录本分别划分为24,54个功能类别及291条代谢通路中。25 197个差异表达的基因在30条代谢通路中显著富集,与酶活性、激素合成代谢和信号转导等相关的差异基因积极响应枇杷干旱,其中双萜类、油菜素类固醇和类胡萝卜素3条生物合成途径中的差异基因呈现出较为一致的表达。采用实时荧光定量(qRT-RCR)对选取的差异基因进行验证,其中过氧化物酶、蛋白激酶byr2、丝氨酸苏氨酸蛋白激酶、脱落酸8'-羟化酶、吲哚-3-乙酰乙酸合酶相关基因上调表达,细胞色素P450 734A1基因下调表达。为枇杷提供了较为全面的基因信息和代谢途径数据,为干旱胁迫下枇杷分子调控机制的深入研究奠定了基础。     

灵武长枣嫁接和根蘖繁殖植株果实转录组差异分析

为解析嫁接繁殖和根蘖繁殖灵武长枣果实糖和酸转化差异,本研究以两种繁殖方式的成熟期果实为材料,利用高通量测序平台Illumina/Nova-Seq 2500进行了转录组测序,对转录组数据质量评估获得两样品的序列长度分别为24.48 Gb (嫁接)和28.64 Gb (根蘖),错误率为0.03%,GC含量44.11%～44.56%,Q30大于90.54%,转录组数据与冬枣基因组序列比对,所有的转录组Reads均分布于冬枣12条染色体上,基因定位的Reads数量嫁接繁殖比根蘖繁殖多;差异表达基因分析筛选到27个显著差异表达基因,其中9个显著上调,18个显著下调;GO富集分析表明15个差异基因显著富集到GO的三大功能类别分子功能、细胞组分、生物过程中;KEGG富集分析表明10个差异表达基因显著富集到12条代谢通路中,其中包括3个上调基因和7个下调基因;进一步筛选到调控糖及有机酸代谢关键差异基因BAM1、PCO和HAOX1在淀粉和蔗糖代谢、牛胆素和亚牛磺酸代谢、乙醛酸和二羧酸代谢中显著表达,并在糖和有机酸的相互转化中起着关键调控作用。本研究表明嫁接繁殖和根蘖繁殖灵武长枣果实在糖和有机酸代谢存在差异,为灵武长枣优质果品繁育提供理论基础。     

银杏三个年龄阶段雄花芽的转录组分析及差异基因筛选

本研究筛选银杏三个年龄阶段开花调控的关键基因,揭示银杏开花调控年龄途径的分子机理,为促进银杏早花、分子育种和推广种植提供科学依据。本研究利用高通量测序技术及生物信息学工具,对三个年龄阶段的银杏雄花芽进行转录组测序,对测序数据进行分析,筛选出年龄途径关键差异表达基因。转录组测序转录组测序共产生57.45 Gb原始数据,注释到8大功能数据库(GO, COG, KEGG, KOG, NR, Pfam, SwissProt, egg NOG)上的Unigene总数为35 058个。通过KEGG通路注释和GO富集分析,将Unigene分别归于55个GO类别和126个代谢途径。差异表达基因分析显示,幼年雄花芽较开花一年雄花芽有37个基因上调,75个基因下调;幼年雄花芽较开花多年雄花芽有592个基因上调,871个基因下调;开花一年雄花芽较开花多年雄花芽有961个基因上调,1 203个基因下调。本研究发掘出大量的开花相关的基因,最终筛选出SPL (gene.Gb_23724, gene.Gb_03922)、AP2 (gene.Gb_00766)、MADS-box (gene.Gb_01886, gene.Gb_15398, gene.Gb_28337, Gingko_newGene_2213)、GA调节蛋白(gene.Gb_34467, gene.Gb_28606, gene.Gb_33214)和DELLAs蛋白(gene. Gb_34644)等11个年龄途径开花调控的关键基因。     

与甘蔗间作的花生叶片转录组分析

为研究间作甘蔗对花生相关基因的表达调控及响应机制,以花生叶片作为材料,使用Illumina NovaSeq 6000高通量测序技术,对间作花生离甘蔗较近的边行(IPA_leaf)、间作花生离甘蔗较远的中间行(IPB_leaf)及单作花生(MP_leaf)叶片进行转录组的测序分析。将各样本获得的转录组数据进行两两比较,MP_leaf_vs_IPA_leaf、MP_leaf_vs_IPB_leaf和IPB_leaf_vs_IPA_leaf的差异表达基因数分别是1 806个、968个、1 302个,其中上调表达基因分别占61.13%、48.24%、65.05%。GO功能富集分析显示差异表达基因主要涉及细胞壁组织或生物合成、碳水化合物生物合成与代谢过程、细胞多糖生物合成与代谢过程、细胞葡聚糖生物合成与代谢过程、催化活性、氧化还原酶活性及过程等。KEGG富集分析表明差异表达基因主要涉及的代谢通路有苯丙素类生物合成、α-亚麻酸代谢、淀粉和蔗糖的代谢、类黄酮生物合成等。本研究通过花生叶片转录组分析,为研究间作花生叶片基因的表达调控机制提供基础数据资料。     

茶饼病诱导感、抗茶树品种的基因表达谱分析

为揭示茶树被茶饼病危害诱导的防御反应机制,本研究选择抗病和感病茶树品种为材料,通过转录组测序和数字表达谱分析茶树叶片被茶饼病危害前后的基因表达差异。结果显示,共筛选出差异表达基因974条,其中共有的为122条,抗病品种中特异364条,感病品种中特异的为488条。对差异表达基因进行分析发现,茶饼病危害主要影响了茶树体内代谢途径、内质网蛋白质加工、次生代谢产物的生物合成、植物病原相互作用、植物激素信号转导途径、淀粉与蔗糖的代谢、苯丙醇生物合成等通路中关键基因的表达水平,这些差异基因包括抗病蛋白基因(R protein)、水解酶基因、细胞壁加固基因、转录因子基因、植物激素及其信号转导基因、次生代谢和氧化酶类、转运蛋白等。利用RT-qPCR对筛选的6个基因进行验证,其表达模式与测序结果一致。本研究初步明确了茶饼病侵染对茶树基因转录水平的影响,为揭示茶树抗病的分子机制奠定了基础。     

转录组学鉴定高粱营养器官间功能差异的决定基因

本研究采用Illumina HiSeq 4000高通量测序技术,对高粱根、茎和叶的转录组进行测序以及生物信息学分析,以期鉴定出决定营养器官功能异质性的核心基因。结果表明,转录组测序获得165 413 338条原始测序序列(Raw reads),过滤后获得160195618条序列(Clean reads);共鉴定出33204个基因及230685个可变剪接事件;三种营养器官基因表达模式的相似度较低,共筛选出10 083个差异表达基因,其中在不同营养器官均有表达差异的决定基因1320个。对决定基因进行GO富集分析,得到2946个GO功能注释(p0.05),分布于光合作用、结构分子活性和细胞组成相关的GO项中;KEGG富集分析表明差异基因在光合作用、黄酮和黄酮醇生物合成以及叶酸参与的"一碳单位"代谢途径富集程度显著。本研究为深入探讨高粱营养器官间功能差异性机制及器官特异性启动子克隆等工作奠定了数据基础。     

基于转录组测序枣胚发育过程差异表达基因分析

为了探索枣胚发育的分子机理,本研究使用Illumina高通量测序平台对中度可育枣品种‘冷白玉’四个发育时期(原胚,小球形胚,心形胚,鱼雷形胚)的幼胚进行转录组测序,共获得135 743 080条高质量Reads,Q3092.88%。样品间共筛选出2 333个差异表达基因,其中有19个差异表达基因是所有样品组合共有的。原胚与小球形胚的幼胚(L2 vs L3)间的差异表达基因(differentially expressed gene, DEG)基因数(867)明显多于其他组合,小球形胚与鱼雷形胚(L3 vs L5)间的DEG基因数(503)明显多于原胚与心形胚(L2 vs L4)间的DEG基因数(498)。GO富集分析显示,四个发育时期的DEG主要共同富集到分子功能、细胞组分、生物过程等GO条目。KEGG富集分析发现,四个时期的DEGs主要共同富集到次生代谢合成、代谢过程、植物信号转导通路,主要参与了苯丙氨酸代谢、α-亚麻酸代谢、氰氨基酸代谢、光合作用-天线蛋白质等代谢途径。通过转录组测序分析了不同发育阶段的‘冷白玉’枣幼胚差异基因表达状况,为寻找枣胚发育的相关基因提供理论参考。     

菠萝蜜叶绿素缺失突变体茎次生生长转录组分析

对于木本植物而言茎的次生生长是植株发育和木材形成的基础,本研究以菠萝蜜(Artocarpus heterophyllus)正常植株为对照(control check, CK),将苗龄40 d的叶绿素缺失突变体(chlorophyll deficient mutant,CDM)和CK次生生长茎作为材料进行转录组测序,测序数据利用Trinity软件进行de novo组装共得到Unigenes 265 709个,去除低表达量(≤0.5)的转录本后筛选得到差异表达基因(differentially expressed genes,DEGs) 43 881个。利用加权基因共表达网络分析(WGCNA)找到与表型相关度最高的模块。将此模块中上下调的基因分别注释到KEGG和GO数据库,KEGG分类结果显示DEGs主要注释到代谢条目,GO分类结果中DEGs主要注释到与膜、细胞以及催化活性相关的过程。同时上调的基因在对活性氧的反应、抗坏血酸和醛酸的代谢、谷胱甘肽代谢等与胁迫相关的通路富集。这些结果可能暗示CDM茎的次生生长受到代谢和细胞过程的调控,生物或非生物因素引起的胁迫可能也参与到影响茎次生生长的机制当中。木本植物茎的次生生长涉及各种代谢通路和信号传导等复杂过程。通过转录组分析本研究为进一步了解白化病对木本植物茎次生生长的影响提供基础。     

外源6-BA调控孕穗期低温后小麦幼穗发育的转录组分析

小麦孕穗期对低温非常敏感,低温胁迫后外源喷施6-苄氨基腺嘌呤(6-BA),能够缓解低温胁迫对小麦造成的伤害,通过转录组测序技术分析6-BA提高小麦抗寒性的分子机制。选用低温敏感型品种皖麦52和低温迟钝型品种烟农19为试验材料,在孕穗期低温胁迫后喷施20 mg L–1的6-BA溶液,以喷施等量蒸馏水处理为对照。观察幼穗形态,并测定幼穗可溶性糖含量和淀粉含量。再通过转录组测序筛选并分析差异表达基因,探究差异表达基因的功能及可能参与的调控通路,采用qRT-PCR方法对测序结果进行验证。外源6-BA处理10d后,与对照相比小麦幼穗形态发育良好,可溶性糖、淀粉含量均升高。转录组结果表明,在皖麦52和烟农19中分别鉴定出22,770个和9866个差异基因,其中有661个基因在两个品种中均上调,从中筛选出ARF5、AGPL1、1-SST、SWEET15等基因,这些基因与调节植物激素水平、淀粉合成、糖代谢等有关。对筛选出的差异基因进行GO和KEGG富集分析。GO注释表明皖麦52和烟农19差异基因的功能都主要富集于细胞结构稳定性、代谢、催化活性等。KEGG富集分析表示信号转导、内源激素水平的调节、碳代谢...     

结合RNA-seq分析和QTL定位筛选甘蓝型油菜萌发期与铝毒胁迫相关的候选基因

随着土壤酸化程度的加剧,铝毒害已成为影响种子萌发质量和作物产量的重要胁迫因子之一。为研究铝毒对油菜种子萌发过程影响的分子机理,本文采用RNA-seq技术对耐铝品系18D300和敏铝品系27011进行转录组分析,共检测到9344个显著差异表达基因[|log₂(fold change)|≥1和FDR≤0.05],其中4406个DEGs基因上调, 4938个DEGs下调。GO富集分析发现,差异表达基因主要与氧化反应、细胞碳水化合物代谢过程、转运蛋白活性等相关。KEGG富集分析表明,差异表达基因主要富集于苯丙烷生物合成、淀粉和蔗糖代谢、MAPK信号通路-植物、植物-病原体相互作用、植物激素信号转导等途径。此外,通过整合RNA转录组测序和铝毒胁迫下油菜萌发期根系相关性状的QTL定位结果,共筛出44个差异表达基因(10个下调和34个上调),这些基因主要与氧化应激、渗透调节、细胞壁修饰、转运蛋白、激素信号传导等功能有关。