枯萎病菌诱导感、抗陆地棉品种的基因表达谱分析

选用陆地棉抗枯萎病品种中棉所12和感枯萎病品种新陆早7号,利用Solexa高通量测序技术分析棉花幼苗受到枯萎病菌侵染后3 h和6 h根部组织的基因表达谱。中棉所12受病菌侵染后3 h与0 h相比,6 h与0 h相比以及6 h与3 h相比,分别鉴定出4447、5481和2559个差异表达基因;新陆早7号分别鉴定出8615、6727和2078个差异表达基因;2个品种比较,在3 h和6 h分别鉴定出1879个和500个差异表达基因。基因功能分析表明,差异表达基因划分为生物学过程、细胞组分和分子功能注释本体,并进一步细分为48个功能类别。代谢通路分析显示,中棉所12和新陆早7号在枯萎病菌侵染后的6 h与0 h相比鉴定出的代谢通路(Pathway)最多,均有126条;在枯萎病菌侵染后6 h,2个品种相比鉴定出的代谢通路最少,仅有89条。各个比对组的代谢通路涉及的基因可被划分为次生代谢产物的生物合成、多糖的生物合成和代谢、环境适应和免疫系统等13类。在环境适应和免疫系统分类中存在着植物与病原菌相互作用代谢通路,涉及到996个已知功能的差异表达基因,有444个基因上调表达,552个基因下调表达。其中,差异表达基因最多的是WRKY转录因子家族基因,其次为丝氨酸/苏氨酸蛋白激酶基因,而DNA损伤修复/耐受性蛋白,JAZ1、RAR1和RPM1互作蛋白,S位点的特异性糖蛋白S6前体以及钙牵引蛋白等6类基因参与该代谢通路的数目最少。最后随机抽取6个基因进行实时荧光定量PCR验证,其结果与测序结果一致。     

海岛棉枯萎病抗性与类黄酮代谢相关基因表达量的相关

枯萎病是危害海岛棉生产的重要因素之一, 研究枯萎病抗性分子机制将为培育抗病海岛棉品种、解决枯萎病对海岛棉的危害问题提供坚实的基础。本研究在前期转录组测序的基础上, 对海岛棉枯萎病抗性差异表达基因进行分析(DEG, Differentially Expressed Gene); 以7个抗病性表现不同的海岛棉品种为材料, 利用qRT-PCR的方法研究抗病差异表达基因在不同接菌时间点的表达量差异, 并分析基因表达量与病情指数的相关性。结果表明, DEG分析得出类黄酮代谢通路相关基因与海岛棉枯萎病抗性有关。qRT-PCR分析显示抗病材料中类黄酮代谢通路关键基因的表达量显著高于感病材料。在接菌后多个时间点, 类黄酮代谢通路中的关键基因TT7、CHI和DFR在抗病材料中的表达量显著或极显著高于感病材料, 其中CHI和DFR基因的表达量与病情指数呈显著负相关。综上所述, 类黄酮代谢通路相关基因对海岛棉枯萎病抗性均有影响, 且CHI、TT7和DFR基因是关键基因。     

海岛棉枯萎病抗性与类黄酮代谢途径基因表达量的相关性

枯萎病是危害海岛棉生产的重要因素之一,研究枯萎病抗性分子机制将为培育抗病海岛棉品种、解决枯萎病对海岛棉的危害问题提供坚实的基础。本研究在前期转录组测序的基础上,对海岛棉枯萎病抗性差异表达基因进行分析(Differentially Expressed Gene,DEG);以7个抗病性表现不同的海岛棉品种为材料,利用qRT-PCR方法研究抗病差异表达基因在接种0~40 h的表达量差异,分析基因表达量与病情指数的相关性。结果表明,DEG分析得出类黄酮代谢通路相关基因与海岛棉枯萎病抗性有关。qRT-PCR分析显示抗病材料中类黄酮代谢通路关键基因的表达量显著高于感病材料。在接菌后多个时间点,类黄酮代谢通路中的关键基因TT7、CHI和DFR在抗病材料中的表达量显著或极显著高于感病材料,其中CHI和DFR基因的表达量与病情指数呈显著负相关。综上所述,类黄酮代谢通路相关基因对海岛棉枯萎病抗性均有影响,且CHI、TT7和DFR基因是关键基因。     

山药品种苏蓣8号抗炭疽病基因的转录组测序与表达分析

为深入研究炭疽病抗病机制提供候选基因，从分子水平探索炭疽菌侵染对山药炭疽病基因转录水平的影响。苏蓣8号是高感炭疽病品系024经组培诱变所育成的高抗炭疽病山药新品种，为挖掘苏蓣8号炭疽病抗性基因，采用胶孢炭疽菌菌株4-2分别接种苏蓣8号和品系024,以未接种苏蓣8号的叶片为对照，利用转录组测序技术分析抗感材料在接种炭疽菌后不同时间点的差异表达基因。结果表明，抗感材料接种24,48,72,96 h共同差异表达基因个数分别为197,132,187,313个，去除各接种时间点共同差异表达的基因数，共获得711个差异表达基因。GO功能富集显示，差异基因主要与响应生物刺激、防卫反应、细胞壁代谢和氧化还原过程等有关。KEGG富集分析表明，生长素、茉莉酸和乙烯等防御相关的植物激素信号转导途径多个基因表达出现差异，其中5个生长素早期响应基因、茉莉酸和脱落酸合成关键酶基因均上调表达，乙烯响应因子ERF036下调表达，可能负向调控山药炭疽菌的侵染；参与次生代谢产物修饰的多个细胞色素P450基因、泛素连接酶及植物甾醇合成酶、防御素和凝集素基因上调表达；WRKY类、MYB类和TIFY类转录因子正向或负向调控抗病...     

基于高通量测序技术鉴定松辽黑猪与长白猪背最长肌差异基因

旨在利用高通量测序技术对松辽黑猪和长白猪背最长肌组织样本进行mRNA测序和差异分析,筛选影响猪肌肉生长、肉质和脂肪沉积的关键基因,从而为猪肉品质研究提供新的参考信息。采集6头松辽黑猪和6头长白猪的背最长肌组织样本,提取其RNA,采用Illumina HiSeq 2500高通量测序技术对mRNA进行测序,并对获得的reads进行比对、注释和差异表达分析,用NOISeq筛选出差异表达基因并进行相关生物学功能富集分析。结果表明,从2个猪种间共筛选出了664个差异表达基因,其中,364个基因在松辽黑猪中高表达,300个基因在长白猪中高表达。通过对差异表达基因进行生物学功能分析,筛选出LPIN1、FADS1、FADS2、PLIN2、PPARGC1A、PRKAG2和ACSL1等基因参与脂类代谢和肌肉发育相关的调控,相关通路为脂肪酸代谢、PPAR信号通路、AMPK信号通路、胰岛素信号通路和脂肪细胞因子信号通路等。     

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

为探明藏黄牛对低氧适应的相关差异表达基因并分析差异表达基因功能和代谢通路,以藏黄牛和三江黄牛心脏组织为研究对象,提取总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细胞受体信号通路参与低氧适应性调控。采用高通量测序技术对藏黄牛和三江黄牛心脏组织进行转录组测序分析,获得大量差异表达基因,为研究动物低氧适应性机制奠定基础。     

水稻芽期响应镉胁迫的转录组分析

为研究不同类型水稻对镉(Cd)胁迫差异的原因,以镉敏感籼稻恢复系昌恢891(CH891)和镉耐受粳稻品种02428为材料,选取种子萌发后,连续镉处理3 d的和未进行镉处理的幼芽,使用高通量Illumina HiSeq 4000测序技术进行RNA-Seq分析,共得到539 524 490个有效读长,与参考基因组的比对率在94.81%～96.82%,GC含量均在49%以上。通过比较分析,共检测到7 204个差异表达基因(DEGs),其中在CH891中特异表达的基因(SCR3)有849个,02428中特异表达的基因(RCR3)有676个,2个品种均有表达但表达量存在差异的基因(CCR3)有770个。KEGG和GO富集通路分析表明,异黄酮生物合成、磷酸肌醇代谢、黄酮类生物合成和花青素生物合成途径在SCR3中显著富集,细胞核、自噬调节和胰岛素抗性途径在RCR3中显著富集,α-亚麻酸代谢、吞噬体、柠檬烯和蒎烯降解和脂肪酸降解途径在CCR3中显著富集。结果表明,Cd处理后,次级代谢过程主要在CH891中富集,蛋白质代谢主要在02428中富集,综上所述,CH891和02428响应的代谢路径存在差异,控...     

甜瓜与白粉病菌非亲和互作的数字基因表达谱分析

为探讨甜瓜与白粉病菌非亲和互作在转录水平的分子机制,本研究选用高抗甜瓜品种‘Yuntian930’,利用Solexa高通量测序技术分析甜瓜幼苗接种白粉病菌前0 h和接菌后24 h、48 h和72 h的基因表达谱。与未接种相比,接种后3个时间点分别鉴定出219个、1 784个和2 371个差异表达基因,且上调表达基因数多于下调表达基因数。GO功能分析表明,已注释的差异表达基因显著富集到代谢过程、生物学过程、磷代谢过程、磷酸化及蛋白磷酸化修饰、光合膜、类囊体、氧化还原酶活性和序列特异结合位点等生物学过程。Pathway分析显示,差异表达基因富集到多条物质代谢、次级代谢物合成、甘氨酸/丝氨酸/苏氨酸代谢、光和系统、植物激素信号转导等通路中。比较各接种时间的差异表达基因显示,接种后48 h鉴定出的差异基因和Pathway最多,并发现有7个涉及多胺代谢的差异表达基因,其中6个呈上调表达。这一结果进一步验证了多胺代谢途径基因能有效响应白粉病菌胁迫反应,从而提高对白粉病菌的抵抗能力。利用q RT-PCR验证4个多胺代谢相关基因在接种白粉病菌后的差异表达,其结果与DGE一致,证实了DGE测序结果的可靠性。     

假禾谷镰孢侵染小麦后3种植物激素相关基因的差异表达分析

假禾谷镰孢(Fusarium pseudograminearum)是我国近年新发现的小麦病原真菌,本研究目的是揭示植物激素信号传导途径中相关基因表达对假禾谷镰孢侵染的响应。利用假禾谷镰孢野生菌株WZ2-8A侵染小麦品种周麦24,对接种后5 d和15 d样品进行转录组测序,分析差异表达基因,并对候选基因进行q RT-PCR验证。假禾谷镰孢侵染后,小麦幼苗生长受到明显抑制,根长、株高、根重和地上部分鲜重均显著降低。转录组分析结果表明,植物激素信号传导途径中有29个基因差异表达,涉及到生长素、细胞分裂素和脱落酸3种植物激素。在接种后5 d有11个差异表达基因,其中2个上调,9个下调;在接种后15 d共有25个差异表达基因,其中8个上调,17个下调。在生长素信号传导途径中,生长素输入载体AUX1差异表达,影响生长素的极性运输,从而影响小麦根部的细胞伸长。在细胞分裂素信号传导途径中,起正调控作用的B-ARR上调表达,推测其促进细胞分裂素的信号传导,从而抑制细胞分裂,与生长素协同作用,造成小麦的长势减弱。在脱落酸途径中,脱落酸受体PYR/PYL下调表达;起到负调控作用的PP2C相关基因均上调表达。脱落酸使植物对真菌和细菌的抗性起到负调控作用,其信号传导途径与茉莉酸/乙烯途径相互拮抗,脱落酸信号传导途径的阻遏可能会使茉莉酸/乙烯途径信号通路打开。q RT-PCR结果基本能够和转录组测序结果相拟合,说明在假禾谷镰孢侵染胁迫下,脱落酸的信号传导被抑制可能是中抗品种周麦24对假禾谷镰孢产生一定的抗性的生理基础。     

花生油脂合成相关基因的表达谱分析

为研究不同发育时期花生籽仁油脂合成过程中基因表达调控模式,本研究以高油酸、中油花生品系F18和低油酸、低油花生品种‘鲁花6号’为材料,对下针后10、30、40、60 DAP(days after pegging)的花生种子进行表达谱芯片测序。结果表明,130、3556、2783个基因分别在30、40、60 DAP时期差异表达。GO注释和KEGG富集结果显示,差异表达的基因主要富集在脂肪酸合成和光合等代谢进程中,其中FAB2、FAD2、WRI1等主要参与油酸的积累,参与光合作用的基因均为捕光叶绿素a/b结合蛋白,全部上调表达。代谢通路图结果表明,籽仁发育的40 DAP和60 DAP时期,脂肪酸合成途径的基因均上调表达。研究结果为花生油脂代谢的分子机制提供理论基础,同时也为花生品质改良贡献了基因资源。     

利用WGCNA鉴定棉花抗黄萎病相关基因共表达网络

加权基因共表达网络分析(Weighted Gene Co-expression Network Analysis,WGCNA)作为一种典型的系统生物学方法,可用来鉴定协同表达的基因模块,探索模块与目标性状之间的生物学相关性,并挖掘网络中的核心基因。黄萎病(Verticillium wilt)可严重降低棉花(Gossypium spp.)的产量及品质,因此研究抗病基因及抗病机制,对于棉花抗病育种工作具有重要意义。本研究以黄萎病菌侵染不同时间点的海岛棉(Gossypium barbadense)幼苗根尖转录组数据为基础,分析其差异表达基因,并选取在不同样本之间表达水平变异最大的前50%基因(共35,647个)进行WGCNA。结果表明,在黄萎病菌侵染条件下,共有22,850个基因差异表达,其中4685个为所有侵染时间点共有的差异基因。共鉴定到18个基因共表达模块,其中5个为抗黄萎病相关特异性模块(black、mediumpurple3、darkolivegreen、plum3模块分别与侵染2 h、6 h、48 h、72 h时间点正相关;mediumpurple2与侵染2 h时间点负相关)。GO和KEGG富集分析表明,特异性模块可以富集到有生物学意义的富集结果,如刺激响应、钙离子结合、黄酮类化合物合成代谢等抗逆相关通路。通过计算模块内基因的连通性,挖掘出了网络中的核心基因,功能预测表明,这些基因可能在逆境胁迫响应中发挥重要作用。本研究为进一步理解棉花抗黄萎病的分子机制、选育棉花抗病新品种提供了理论指导。     

海岛棉枯萎病抗性与糖基转移酶类基因表达量的相关性

枯萎病是危害海岛棉生产的重要因素之一，研究枯萎病抗性分子机制是培育抗病海岛棉品种的分子基础。基于对前期转录组测序数据分析，本研究对参与海岛棉抗枯萎病调控的糖基转移酶基因进行了初步研究，以8个不同抗病性的海岛棉品种为材料，在枯萎病菌处理后，利用实时荧光定量PCR检测9个棉花抗枯萎病相关基因进行表达量分析，结果表明，GB＿A03G0575在海岛棉枯萎病侵染过程中随着时间的推移，表达量会先减少再增加，推断其可能在受到病菌胁迫时植物生长会受到影响。GB＿A13G1825和GB＿A11G1787基因在抗病材料中的表达量会比感病材料的多，最大表达量出现的时间，感、抗材料基本一致，推测这些基因对抗病有一定影响，表明在抗病过程中三个基因可能起着比较重要的作用。本研究结果为海岛棉抗枯萎病育种及抗病机理提供候选基因资源。     

水稻植物激素响应低温胁迫反应的转录组分析

植物激素在调节和控制植物生长发育、代谢过程、应对逆境等方面具有重要的作用。为明确水稻苗期植物激素对低温的应答机制,本研究以野生型粳稻品种‘中花11’为研究材料,经低温处理后,利用RNA-Seq技术分析植物激素调控水稻幼苗对低温胁迫的应答模式。通过转录组测序,共获得9.9×10~7条干净的序列,筛选出2 044个差异表达基因(DEGs)。首先,GO富集分析结果表明差异表达基因主要富集在生物学过程、细胞组分和分子功能三个方面。进一步的KEGG通路分析表明,差异表达基因主要富集在代谢途径、次生代谢生物合成、植物激素信号转导等途径。其中植物激素信号转导中低温应答的差异表达基因为31个,分别为25个上调表达基因和6个下调表达基因;主要参与了茉莉酸和脱落酸信号途径。这些研究结果对水稻苗期植物激素的低温应答机制完善和栽培调控具有重要的参考价值。     

缺磷胁迫对黑籽南瓜幼苗根系生长和根系分泌物的影响

采用营养液培养方法,以1/2单位日本园试配方为基本营养液,在温室内培养黑籽南瓜(Cucurbita ficifolia B.)幼苗,测定正常供磷处理和缺磷胁迫下黑籽南瓜幼苗在根系生长、根系形态和根系分泌物等方面的差异,以期为南瓜磷高效种质创新和耐低磷机制研究提供理论依据.结果表明,缺磷胁迫后21 d,黑籽南瓜幼苗地上部生物量显著降低,根系生物量与正常供磷处理相比没有明显差异,根冠比增大;总根长和根表面积显著增加,根数、根横径和根体积与正常供磷处理相比没有明显的差异;根系吸磷量和磷转运率下降,根系吸磷量占总吸磷量的百分数增加,根系磷利用效率显著增高.总根长和根表面积增加可能与其磷吸收效率密切相关,是黑籽南瓜幼苗对缺磷胁迫的适应性变化.与正常供磷处理相比,缺磷胁迫后21 d,黑籽南瓜幼苗根系分泌的化合物种类组成类似,主要由酚类、芳香烃类、酯类、胺类、烯烃类和烷烃类化合物组成.但缺磷胁迫后根系分泌物中的化学成分多于正常供磷处理,邻苯二甲酸二(1丁基2异丁基)酯等8种化合物是缺磷胁迫下黑籽南瓜根系分泌物的特有成分.缺磷胁迫后黑籽南瓜的自毒作用增强.     

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

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

肉牛大理石花纹等级相关lncRNA筛选与特征分析

牛肉大理石花纹特征是由多基因调控的重要肉质性状,是衡量高档牛肉品质的条件之一。LncRNA是一种重要的非编码RNA,在调节脂肪沉积的相关生物学过程中发挥重要作用。为了解lncRNA在肉牛不同大理石花纹表型中的潜在生物学作用,采集高(A5)、低(A1)2种极端等级大理石花纹的杂交和牛背最长肌组织,利用高通量测序方法检测与肉质大理石花纹相关的差异lncRNA,预测差异lncRNA反式作用靶基因,并对其进行功能富集分析。结果表明,在2个等级大理石花纹牛肉组织中共筛选出45个显著差异lncRNA,其中在高等级大理石花纹组织中分别有21,24个显著上调及下调的lncRNA,包括MSTRG.9509.3、ENSBTAT00000077612、ENSBTAT00000072841等与脂肪沉积相关的lncRNA。对反式作用靶基因的功能富集分析表明,预测靶基因显著富集到脂多糖结合、脂肪分解等脂肪沉积相关的生物学过程。差异lncRNA可能通过反式作用靶向脂肪沉积相关基因(FABP4、PLIN1、LIPE等)参与调节肉牛大理石花纹表型特征。检测到的牛肉大理石花纹相关差异lncRNA、预测靶基因及其富集通路将...     

外源DNA导入冬瓜条件的优化及其运用

针对冬瓜(Benincasa hispida(Thunb.)Cogn.)缺乏抗性种质资源的现状,本文开展了冬瓜抗枯萎病种质资源创新研究.采用花粉管通道法将高抗枯萎病的黑籽南瓜DNA导入优质冬瓜材料,以外源DNA剂量、pH值、注入时间、注入方法以及冬瓜品种为因素,通过正交设计,设计处理组合,筛选和优化外源DNA导入冬瓜条件.结果表明:以子房注射法、授粉后24 h、导入浓度为1.2 mg/mL DNA溶液(pn 7.5)的处理,对冬瓜坐果率的影响最小.利用此方法将高抗冬瓜枯萎病的黑籽南瓜总DNA导入冬瓜,经病圃田间筛选和6代自交纯化已获得2份稳定的冬瓜新材料.经苗期人工接种和病圃田间自然抗病鉴定,其抗病性明显提高.所获抗病材料的品质与受体无明显差异.其中1份抗枯萎病材料已运用于育种实践,培育出了抗病、优质粉皮冬瓜"广利1号",目前该品种每年推广面积700 hm2以上.     

开发黑籽南瓜发展高效农业

黑籽南瓜是葫芦科南瓜属中种子为黑色的栽培种,具有较高的经济价值,我国仅在云贵高原乌蒙山区有零星种植。近年来,由于黄瓜和西瓜用黑籽南瓜作砧木嫁接苗栽培,可有效地防治枯萎病等瓜类土传病害,同时可大幅度提高产量和改善品质,其黑籽南瓜的经济价值培增。我省具有黑籽南瓜生长的独特生态条件和资源优势,抓住时机开发我省黑籽南瓜资源,对我省“三高”农业发展具有重要的意义。同时,由于黑籽南瓜分布在我省高寒贫困的乌蒙山区,还可以为我省高寒贫困山区人民脱贫致富开辟新途径。     

转录组学挖掘黑果枸杞中花青素生物合成代谢关键bHLH调控基因

黑果枸杞果实中富含花青素,为探究黑果枸杞和白果枸杞的基因表达差异,挖掘控制黑果枸杞黑果性状的主效基因,本研究利用RNA-Seq技术对黑果枸杞和白果枸杞进行转录组测序与生物信息学分析。结果表明,黑果枸杞和白果枸杞中共检测到25 279个基因差异表达,其中12 381个基因在黑果枸杞表达上调,12 898个基因表达下调。KEGG富集结果显示,共有571个差异表达基因参与植物次生代谢产物合成,尤其是在黄酮类和花青素生物合成代谢通路中差异最为明显。黑果枸杞中共有8个参与花青素合成代谢的结构基因的表达量均高于白果枸杞,说明花青素的生物合成可能与黑果枸杞黑果性状相关。调控花青素生物合成的bHLH转录因子CL8159.Contig5＿All的表达水平比白果枸杞增强了2 139.57倍,这个基因很有可能是控制黑果枸杞黑果性状的候选基因。本研究筛选到控制黑果枸杞花青素合成代谢的候选基因,为黑果枸杞的分子遗传机理及进一步的基因功能研究提供了依据。     

大丽轮枝菌侵染陆地棉早期的转录组分析

【目的】旨在深入挖掘大丽轮枝菌致病相关基因。【方法】运用转录组测序技术比较分析了强致病力大丽轮枝菌Vd991在侵染陆地棉品种Coker 312早期阶段的基因表达情况。【结果】对Vd991侵染前后的转录组数据进行差异表达分析,共筛选到979个差异表达基因(Differentially expressed genes,DEGs),其中376个上调,603个下调。通过基因注释和功能分类发现,这些DEGs涉及细胞增殖与生长、产孢、碳水化合物结合、蛋白激酶活性调节、裂解酶活性、水解酶活性等功能。对上调DEGs进行基因本体(Gene ontology,GO)功能富集分析,发现共有277个GO词条达到显著富集水平(P0.05),涉及生长速率正向调控、核苷酸合成、解旋酶活性等功能;京都基因与基因组百科全书(Kyoto encyclopedia of genes and genomes,KEGG)通路分析表明,上调DEGs参与53个代谢通路,包括嘧啶代谢、嘌呤代谢等。【结论】在早期侵染过程中,大丽轮枝菌细胞增殖相关基因表现活跃。上述分析为进一步揭示大丽轮枝菌的致病机理提供了有价值的参考。