球孢白僵菌诱导的柑橘木虱免疫应答转录组分析

【目的】筛选柑橘木虱Diaphorina citri Kuwayama应对球孢白僵菌侵染的免疫应答及其网络调控基因,以进一步探讨柑橘木虱对球孢白僵菌的免疫防御机制。【方法】采用Illumina高通量测序平台对感染球孢白僵菌24、48、72 h与健康的柑橘木虱转录组进行了测序,利用生物信息学软件对筛选到的差异表达基因进行了功能注释、分类以及参与的信号通路分析。【结果】组装得到138 313条不可延长的非冗余Unigene,其N50和N90分别为2 532 bp和413 bp,平均长度为1 191.26 bp。在CK vs.S24h、CK vs.S48h和CK vs.S72h 3个转录组里分别获得了971、1 671和752个显著差异表达基因(DEGs),GO富集分析表明分别有405、614和542个DEGs显著富集到56、83和60个GO terms中,KEGG pathway分析结果显示分别有98、333和247个DEGs显著富集到10、27和25条代谢通路里。【结论】差异表达基因主要富集在能量代谢、离子运输、转录和翻译调控、生殖和发育调控以及免疫防御反应等相关通路,大部分编码潜在的与免疫识别及调控相关的基因,并从中筛选到5个显著上调表达的免疫相关基因,为开展柑橘木虱免疫应答虫生真菌的研究奠定理论基础。下一步将进行免疫相关基因的q-PCR验证及表达量分析,研究其在柑橘木虱应对球孢白僵菌侵染过程中的作用。     

BbRho5对球孢白僵菌生长速率的作用研究

为阐明BbRho5对球孢白僵菌生防潜能的作用,构建了Bbrho5单基因敲除菌株ΔBbrho5,以野生型菌株WT作为对照,在不同培养基上测定菌落生长速率,并测定了菌株对多菌灵胁迫耐受性及对大蜡螟幼虫体壁侵染能力。进一步获取和分析了ΔBbrho5和WT细胞内基因转录组数据。结果表明,BbRho5蛋白功能缺陷显著抑制球孢白僵菌菌丝生长速率,同时微弱影响其多菌灵胁迫抗逆性及生防能力。相较于WT,ΔBbrho5中具有770个差异表达基因(DEGs),其中上调基因395个,下调基因375个。GO分析显示,ΔBbrho5 VS WT中DEGs主要富集于氧化还原酶活力(oxidoreductase activity)和单加氧酶活力(monooxygenase activity)功能。KEGG通路富集结果显示,DEGs主要富集于氮代谢及多种氨基酸代谢通路。在氮代谢通路中富集到7个功能基因,其中有5个上调,2个下调,说明敲除菌株可能采用增强氮源利用及谷氨酸合成以应对Bbrho5缺陷引起的生长迟缓。以上研究结果揭示了球孢白僵菌中小GTP酶BbRho5对球孢白僵菌生长速率具有重要影响,且氮代谢和氨基酸代谢可能为其重要的响应代谢通路。     

BbRho5对球孢白僵菌生长速率的作用研究

为阐明BbRho5对球孢白僵菌生防潜能的作用,构建了Bbrho5单基因敲除菌株ΔBbrho5,以野生型菌株WT作为对照,在不同培养基上测定菌落生长速率,并测定了菌株对多菌灵胁迫耐受性及对大蜡螟幼虫体壁侵染能力。进一步获取和分析了ΔBbrho5和WT细胞内基因转录组数据。结果表明,BbRho5蛋白功能缺陷显著抑制球孢白僵菌菌丝生长速率,同时微弱影响其多菌灵胁迫抗逆性及生防能力。相较于WT,ΔBbrho5中具有770个差异表达基因(DEGs),其中上调基因395个,下调基因375个。GO分析显示,ΔBbrho5 VS WT中DEGs主要富集于氧化还原酶活力(oxidoreductase activity)和单加氧酶活力(monooxygenase activity)功能。KEGG通路富集结果显示,DEGs主要富集于氮代谢及多种氨基酸代谢通路。在氮代谢通路中富集到7个功能基因,其中有5个上调,2个下调,说明敲除菌株可能采用增强氮源利用及谷氨酸合成以应对Bbrho5缺陷引起的生长迟缓。以上研究结果揭示了球孢白僵菌中小GTP酶BbRho5对球孢白僵菌生长速率具有重要影响,且氮代谢和氨基酸代谢可能为其重要的响应代谢通路。     

感染球孢白僵菌的家蚕免疫应答差异表达基因分析

【目的】筛选家蚕Bombyx mori应对白僵菌Beauveria bassiana侵染的应答基因, 以进一步研究家蚕抵御真菌侵染的分子机制。【方法】采用新一代Solexa高通量测序技术对感染白僵菌及未感染白僵菌的对照组家蚕进行了测序分析, 筛选差异表达基因; 结合生物信息学工具分析差异表达基因的功能注释、 分类及涉及的信号通路等; 应用荧光定量PCR技术验证10个基因的差异表达。【结果】通过测序和生物信息学分析共获得377个差异表达基因, 其中表达上调基因236个, 下调基因141个; KEGG通路分析表明, 各通路中既有表达上调的基因, 也有下调基因; 12个上调基因、 26个下调基因参与3个显著性富集的KEGG通路, 即核糖体、 氨酰tRNA生物合成和剪接体通路。定量PCR与测序结果显示, 溶菌酶、 热激蛋白、 谷胱甘肽S-转移酶、 肽聚糖识别蛋白等与免疫应激相关的蛋白基因均呈现表达上调。【结论】本研究筛选获得的差异表达基因, 特别是上调表达的基因可能与家蚕应对白僵菌侵染的应答机制有关, 其中与免疫应激相关的蛋白基因如溶菌酶、 热激蛋白、 谷胱甘肽S 转移酶、 肽聚糖识别蛋白基因等可能直接参与了家蚕对白僵菌的免疫识别和防御, 研究结果为从分子水平阐明家蚕抵御真菌侵染的防御机制和白僵菌对家蚕的致病机理提供新的依据。     

玫烟色棒束孢(Isaria fumosorosea)侵染小菜蛾的表达谱分析

【背景】昆虫病原真菌对寄主的侵染是一个十分复杂的过程,是多基因共同作用的结果。玫烟色棒束孢(Isaria fumosorosea) IFCF01菌株对小菜蛾具有很高的致病力,然而有关玫烟色棒束孢对小菜蛾致病的相关基因少见报道。【目的】筛选玫烟色棒束孢侵染小菜蛾相关基因,为更好地利用玫烟色棒束孢防治小菜蛾提供基因靶点。【方法】采用第二代高通量测序技术RNA-Seq,对玫烟色棒束孢侵染小菜蛾2–3龄幼虫4、8、12、16、24、30、36 h的虫菌混合样品(处理组)及纯培养玫烟色棒束孢(对照组)进行测序分析并筛选差异表达基因,结合生物信息学方法分析差异基因涉及的功能模块和信号通路。【结果】玫烟色棒束孢侵染小菜蛾混合样品与纯培养玫烟色棒束孢对照组对比分析共获得28 384个差异基因,其中显著差异表达基因274个,上调表达118个,下调表达156个。筛选获得的显著差异表达基因,特别是上调表达基因可能与玫烟色棒束孢对小菜蛾的侵染有关。GO二级分类显示,差异表达基因能够注释到36个GO条目中,包含18个生物学过程、9个细胞组分和9个分子功能。KEGG通路分析显示共有171个差异表达基因(differentially expressed gene,DEG)注释到132个通路中,其中有66个DEG显著富集在14个通路中。这些显著差异表达基因中大部分为玫烟色棒束孢侵染过程中潜在致病毒力相关基因。【结论】本研究为筛选玫烟色棒束孢侵染小菜蛾致病相关基因提供重要数据库,也为阐明玫烟色棒束孢对小菜蛾的侵染机制提供基础。     

玫烟色棒束孢诱导的小菜蛾免疫响应表达谱分析

【目的】本研究旨在筛选小菜蛾Plutella xylostella应对玫烟色棒束孢Isaria fumosorosea侵染的免疫应答及其网络调控基因,以进一步探讨小菜蛾对玫烟色棒束孢的防御机制。【方法】采用第二代高通量测序技术RNA-seq,对处理后12 h的感染玫烟色棒束孢和健康(平行对照)的小菜蛾4龄幼虫转录组进行了测序,通过生物信息学分析对得到的差异基因进行了功能注释和分类,对差异基因参与的信号通路进行了分析。【结果】在感菌和健康小菜蛾幼虫转录组两个表达谱里,分别获得了12 346 987和12 315 210个clean reads,有60.93%和61.26%的reads分别能比对到参考基因库里,其中完美匹配(perfect match)的比例分别为32.15%和32.73%。共得到351个显著差异表达基因(differentially expressed unigenes,DEUs),上调表达基因有275个,下调表达基因有76个,与免疫防御反应潜在相关的基因有156个。GO(Gene Ontology)富集分析有102个DEUs分布到46个GO term里,KEGG(Kyoto Encyclopedia of Genes and Genomes)pathway富集分析结果显示,有132个DEUs显著富集在13个代谢通路(pathway)里。【结论】这些差异表达基因中,大部分编码潜在的与免疫识别及调控相关的基因,集中在能量代谢、疾病反应和防御反应等相关通路。研究结果为挖掘与玫烟色棒束孢侵染相关的小菜蛾免疫应答候选基因提供了重要数据库,也为阐明小菜蛾对玫烟色棒束孢的免疫机制奠定理论基础。     

感染球孢白僵菌的烟粉虱若虫免疫应答转录组分析

【目的】筛选烟粉虱Bemisia tabaci(Gennadius)若虫应对白僵菌Beauveria bassiana侵染的应答基因,以进一步研究烟粉虱免疫反应的分子机制。【方法】采用新一代高通量测序技术对感染和非感染白僵菌的烟粉虱4龄若虫进行了测序分析,并筛选了差异表达基因;利用生物信息学工具对转录组测序得到的基因进行了功能注释、分类以及参与的信号通路展示。【结果】组装得到非冗余Unigene 232 554个,其N50和N90分别为1 153 bp和260 bp,平均长度为67 424 bp。对所有的基因进行差异表达分析发现:以P<0.05为标准筛选得到了1 166个差异表达基因,其中上调表达基因有474个,下调表达基因有692个,其中,与免疫相关的基因有405个;GO富集分析发现:有416个GO term有富集现象,包括156个生物学过程(66 402个Unigenes),89个细胞组分(27 645个Unigenes)和154个分子功能(73 417个Unigenes);KEGG代谢通路分析将烟粉虱转录组1 145个DEG匹配到309个通路上,其中76个通路得到了富集。【结论】对405个可能参与到烟粉虱若虫对白僵菌侵染的免疫识别和防御的基因进行了测序。研究结果为从分子水平上开展应用虫生真菌防治烟粉虱的研究奠定信息学基础。     

利用GFP表达系统检测球孢白僵菌侵染昆虫过程

绿色荧光蛋白(GFP)作为一种重要的报告基因,在病原菌和寄主之间的互作研究中具有良好的应用前景。利用草丁膦抗性基因bar为筛选标记,构建了组成性表达绿色荧光蛋白基因egfp的载体pBG,并转入球孢白僵菌获得成功表达。将转化子接种菜青虫进行生物测定表明,组成性表达egfp不影响球孢白僵菌的毒力。进一步利用冰冻切片和荧光显微观察技术监测球孢白僵菌侵染寄主行为,结果显示,通过荧光观察可清楚的检测到病原菌在昆虫体表的附着、菌丝穿透体壁以及菌丝从寄主体内长出等侵染过程。由此表明,GFP可有效用于球孢白僵菌的遗传标记,进行病原菌的鉴别、侵染致病机理及相关功能基因的表达模式研究。     

冰菜盐胁迫下的转录组分析

为了解冰菜(Mesembryanthemum crystallinum)叶片抗盐相关基因组学,利用Illumina Hi-seq TM2500高通量测序技术研究冰菜叶片在400 mmol L~(–1) NaCl胁迫下转录组基因的差异表达。结果表明,从400 mmol L~(–1) NaCl胁迫和对照的冰菜叶片中共获得13.01 Gb Clean data,Q30碱基均大于90.08%。共获得123个差异表达基因(DEGs),包括73个上调基因,50个下调基因,其中功能注释的基因有96个。根据Unigene库序列进行GO、COG和KEGG注释,筛选出8个与抗盐性相关差异表达基因,植物激素代谢相关基因,脱落酸8'-羟基化酶、吲哚-3-乙酰酸酰胺合成酶和茉莉酮酸酯ZIM结构域蛋白基因均下调表达,生长素响应蛋白、细胞分裂素合酶基因则上调表达,糖代谢相关基因棉子糖合成酶基因上调表达,质膜H+-ATPase基因上调表达,脱水蛋白基因下调表达。这为冰菜耐盐基因组学和分子生物学的研究奠定基础。     

球孢白僵菌Pr1蛋白酶基因cDNA克隆与序列分析

从球孢白僵菌菌株Bb202中克隆了昆虫蛋白酶Pr1的基因。使用特异性引物,分别以基因组DNA和mRNA为模板进行PCR反应,得到了球孢白僵菌的Pr1蛋白酶基因cDNA和结构基因序列。测序结果表明球孢白僵菌Pr1蛋白酶结构基因全长大小为1606bp,含有3个内含子,分别为69、61和68bp。其ORF全长1140bp,预测编码379个氨基酸残基。成熟蛋白理论分子量为38.8kD,理论等电点为8.06。为进一步研究球孢白僵菌Pr1基因,并构建高毒力工程菌株提供了理论基础。     

Do aphids infected with entomopathogenic fungi continue to produce and respond to alarm pheromone?

The response of pea aphids, Acyrthosiphon pisum, to aphid alarm pheromone was not modified by infection with Beauveria bassiana. Approximately 50% of uninfected and infected aphids responded to synthetic alarm pheromone. The simulated attack of aphids infected with B. bassiana did not elicit a response in uninfected aphids. Preliminary air entrainment experiments of both uninfected aphids and aphids at different stages of B. bassiana (generalist pathogen) or P. neoaphidis (obligate pathogen of aphids) demonstrated that B. bassiana infected aphids produced less alarm pheromone than uninfected aphids and, conversely, P. neoaphidis infected aphids produced more alarm pheromone than uninfected aphids. These results are discussed with particular emphasis on the different life history strategies of these two pathogens. We hypothesise that the obligate, specialist pathogen, P. neoaphidis, is under greater selection pressure to increase pathogen transmission and survival resulting in modified host behaviour, than the generalist pathogen, B. bassiana.     

小菜蛾NanosO基因启动子的克隆及功能验证

【目的】克隆小菜蛾Plutella xylostella NanosO基因PxnosO启动子,并验证其具有生殖腺特异性活性,以期应用于基因功能研究或转基因昆虫的构建,为小菜蛾等农业害虫的综合治理提供新的研究思路。【方法】根据小菜蛾基因组序列信息,利用PCR技术克隆NanosO的启动子并进行序列分析。构建PxnosO-EGFP表达质粒,利用脂质体细胞转染技术将PxnosO-EGFP和IE1-EGFP表达质粒转入到小菜蛾胚胎细胞系(Px-6)和草地贪夜蛾Spodoptera frugiperda卵巢细胞系(Sf9)中,通过激光共聚焦荧光显微镜观察和qRT-PCR技术分别定性和定量分析EGFP基因的表达,验证小菜蛾NanosO启动子的活性。【结果】克隆获得小菜蛾PxnosO (Px004767)启动子区序列,长1 743 bp。对启动子序列进行分析,发现该序列不仅包含启动子共有核心元件TATA box以及上游启动子成分CAAT box和GC box等,还包含有数十个转录因子结合位点。利用细胞转染技术,在PxnosO启动子驱动下成功地在Px-6和Sf9细胞系中表达外源基因EGFP。【结论】克隆了小菜蛾NanosO基因PxnosO启动子,在细胞水平上验证其能驱动外源EGFP基因的表达,为分析PxnosO在小菜蛾不同发育时期的表达模式和PxnosO启动子在体内的功能验证奠定基础。     

Laboratory investigations on the potential of entomopathogenic fungi for biocontrol of Helicoverpa armigera (Lepidoptera: Noctuidae) larvae and pupae

The susceptibility of third instar Helicoverpa armigera to seven strains of three entomopathogenic fungal species, i.e. Metarhizium anisopliae, Beauveria bassiana and Paecilomyces fumosoroseus, was tested under laboratory conditions using the larval immersion method. High efficacies ranging from 68 to 100% corrected mortality were recorded with more profound effects in treatments with B. bassiana and P. fumosoroseus strains. The median lethal concentration (LC₅₀) for L3 was 6.0×10⁵ in M. anisopliae 79, 1.5×10⁵ in B. bassiana 124 and 4.2×10⁴ in P. fumosoroseus 14. These three strains were further used to characterize the age-dependent mortality of different larval stages (L2-L5) and the effect against pupae of H. armigera. Larval stages did not differ in their mortality but differed i in median lethal time, with shorter values recorded in the second instar. Tested fungi also caused a high reduction between 74.4 and 100% in the emergence of pupae using the soil inoculation method and the pupal immersion technique. All three fungal species, especially P. fumosoroseus, have a high potential for biocontrol of H. armigera larvae and also as a soil treatment targeting the pupae.     

果生刺盘孢CfHAC1调控应答二硫苏糖醇胁迫的转录组分析

果生刺盘孢Colletotrichum fructicola是油茶炭疽病优势病原菌。前期研究发现bZIP转录因子CfHac1参与调控该菌的生长发育和致病性。为了揭示转录因子CfHac1调控果生刺盘孢响应内质网压力和致病机理,本研究测定了ΔCfhac1突变体对内质网压力胁迫剂的敏感性,发现突变体对二硫苏糖醇（dithiothreitol,DTT）的耐受性下降,说明CfHAC1基因可能参与调控果生刺盘孢响应内质网压力胁迫过程。进一步利用高通量RNA-seq技术对该病菌野生型菌株和CfHAC1敲除突变体菌株在DTT胁迫下的转录组进行了比较分析,结果表明差异表达基因共有2 680个,其中上调表达基因有1 181个,下调表达基因有1 499个。Gene Ontology 功能分析结果显示,差异表达基因主要参与催化活性、结合、代谢过程、细胞过程、细胞成分合成、生物过程调控和应激反应等生物学过程。KEGG功能富集分析表明,上调表达基因主要被富集到核糖体、真核细胞的核糖体生物合成、RNA转运和氰基氨基酸代谢通路中;下调表达基因显著富集在内质网蛋白质加工、N-聚糖生物合成、类固醇合成和蛋白质分泌等通路中。分析发现转录因子CfHac1调控内质网胁迫应答和致病相关基因的表达。本研究提供了在全基因组水平上对CfHAC1基因与内质网压力胁迫应答之间关联的新认识,为阐明果生刺盘孢响应内质网压力胁迫和致病机制奠定了基础。     

萱草叶片响应低温胁迫的转录组分析

低温胁迫是萱草（Hemerocallis fulva）生长过程中经常会遭遇的一种非生物胁迫。比较了萱草叶片在低温处理（10、5、0 ℃）下转录组与对照（15 ℃）数据的差异,共筛选出差异表达基因2 457个,其中上调基因1 253个,下调基因1 204个。差异表达基因主要富集在细胞过程、代谢过程和催化活性等49个GO过程,代谢途径、次生代谢产物的生物合成、植物激素信号转导等42条KEGG代谢途径中。其中参与植物激素信号转导通路的差异表达基因发生了不同程度的变化,GH3.10基因上调至对照组的13.624倍,IAA1基因下调0.120倍;参与可溶性糖合成通路的差异基因发生了0.076~28.114倍不同程度的变化。随后对3个低温处理组共有的29个差异表达基因进行热图和网络调控分析,基于基因在网络调控中的位置,对ABCF5、OFPs和SWEETs等基因在冷应答的作用进行了分析。本研究结果为进一步挖掘萱草低温响应的关键基因及耐寒萱草种质开发、分子育种提供了一定的理论支撑。     

干扰NSUN2通过调控细胞周期蛋白表达抑制黑素瘤细胞增殖

为探讨太阳结构域家族蛋白酶2(NOP2/Sun RNA methyltransferase 2, Nsun2)RNA甲基化酶在黑素瘤细胞中的生物学功能,本研究以小鼠黑素瘤B16细胞为对象,构建靶向干扰Nsun2基因的shRNA慢病毒干扰载体,包装病毒后感染B16细胞。与对照组相比,干扰组细胞中Nsun2的敲低效率达到80%。EdU染色结果表明,干扰Nsun2显著抑制B16细胞DNA合成能力。转录物组测序技术系统分析干扰组与对照组细胞基因表达水平,共筛选获得1 062个差异表达基因(DEGs),其中678个表达上调,384个表达下调。DEGs主要富集在染色体、着丝粒区、蛋白质结合等GO条目。KEGG分析表明,DEGs显著富集在细胞周期、DNA复制、细胞衰老等通路。荧光定量PCR和转录物组测序结果均发现,Cdk2、Ccna2、Cdc25b等促进细胞分裂的相关基因显著下调,而Gadd45g和Gadd45a等阻滞细胞增殖的基因显著上调。本研究表明,NSUN2通过调控细胞周期和DNA复制等生物学过程,影响黑素瘤细胞增殖,为黑素瘤发生和发展的分子机制研究提供参考。     

家蚕幼虫BmToll9基因对肽聚糖和金黄色葡萄球菌的响应表达

【目的】Toll信号通路是昆虫中重要的免疫信号通路,其中Toll受体在保持Toll通路的正常免疫应答、抵抗外源病原体中起到关键的作用。本研究旨在探究肽聚糖和金黄色葡萄球菌Staphylococcus aureus对家蚕Bombyx mori Toll受体基因BmToll9-1和BmToll9-2表达的影响。【方法】将革兰氏阳性菌细胞壁主要成分肽聚糖和革兰氏阳性菌金黄色葡萄球菌分别注射感染家蚕5龄第1天幼虫,诱导其发生免疫反应;采用实时荧光定量PCR分析注射后不同感染时间点BmToll9-2和BmToll9-1基因在家蚕幼虫中肠、表皮、脂肪体和丝腺中的相对表达水平。【结果】往家蚕5龄幼虫中注射肽聚糖或金黄色葡萄球菌后,BmToll9-2基因出现了时间和组织的差异性表达。注射肽聚糖和金黄色葡萄球菌均能诱导5龄幼虫中肠BmToll9-2基因的表达上调,注射肽聚糖和金黄色葡萄球菌分别在3和6 h时对基因表达的诱导效果最好,且注射金黄色葡萄球菌比注射肽聚糖对基因表达的诱导效果更好。注射金黄色葡萄球菌能引起5龄幼虫表皮、脂肪体和丝腺中BmToll9-2基因的表达上调,分别于注射后24, 6和24 h时诱导效果最好。注射金黄色葡萄球菌亦能诱导同源的BmToll9-1基因的上调表达。【结论】家蚕幼虫BmToll9基因在肽聚糖或金黄色葡萄球菌注射处理后均能在不同组织中发生上调表达,推测BmToll9基因参与了家蚕对肽聚糖和金黄色葡萄球菌的免疫反应。     

Interaction Between the Entomopathogens Beauveria bassiana, Metarhizium anisopliae and Paecilomyces fumosoroseus and the Mycoparasites Clonostachys spp., Trichoderma harzianum and Lecanicillium lecanii

Biocontrol agents of numerous insect pests and fungal pathogens exist but virtually nothing is known about their interaction if used simultaneously. Our objective was to investigate the compatibility of the entomopathogens Beauveria bassiana, Metarhizium anisopliae and Paecilomyces fumosoroseus, and the broad host-range mycoparasites Clonostachys spp., Trichoderma harzianum and Lecanicillium lecanii. In vitro host-range tests revealed that M. anisopliae was highly susceptible to all mycoparasites tested. B. bassiana was attacked by Clonostachys rosea, and P. fumosoroseus. was resistant to mycoparasites. M. anisopliae but not P. fumosoroseus killed nymphs of Bemisia tabaci in bioassays. B. bassiana and M. anisopliae proved lethal to Cosmopolites sordidus, Diatraea saccharalis and Sitophilus oryzae. Coapplication of mycoparasites with entomopathogens did not affect their biocontrol efficacy in vivo, although the reisolation success of entomopathogens could be significantly reduced, especially from smaller insect species. Trichoderma spp. were reisolated from mycoparasite-treated insects more frequently than C. rosea. The coapplication of the highly susceptible M. anisopliae generally enhanced mycoparasite recovery. Mycofungicide preparations caused some insect mortality but less than a copper hydroxide fungicide which is still permissible in organic agriculture. We concluded that the tested entomopathogens and mycoparasites are compatible elements of integrated pest management.     

基于转录组分析重金属镉对长茎葡萄蕨藻的影响

随着人类社会工业化进程的加快,重金属镉(Cd)污染对海洋生物的毒害作用突显。为了探索Cd胁迫下长茎葡萄蕨藻生理调控的分子响应机制,本文采用RNA-Seq技术对镉胁迫(Hcd²⁺)藻体组织转录组进行研究。结果表明: 与对照相比,Hcd²⁺组共筛选出差异表达基因(DEGs)702个,其中,257个基因上调表达,445个基因下调表达。对获得的DEGs进行功能注释及富集分析,发现长茎葡萄蕨藻的多种生物学功能受到Cd²⁺胁迫的影响,最终表现为抑制生长。GO富集分析发现,长茎葡萄蕨藻体内氧自由基的产生与清除出现失衡,并产生了DNA损伤等氧化伤害。KEGG富集分析发现,长茎葡萄蕨藻中多条与光合作用相关途径受到抑制,表明Cd²⁺胁迫导致长茎葡萄蕨藻的光合反应紊乱。