USP1通过调控E2F1对膀胱癌细胞增殖和周期的作用研究

该文探讨了泛素特异蛋白酶1(USP1)对膀胱癌细胞增殖和周期等生物学行为的作用,并进一步探索其作用机制。通过分子克隆技术构建膀胱癌T24细胞USP1过表达细胞株; CRISPRCas9技术构建膀胱癌UMUC3细胞USP1敲除细胞株; CCK8和克隆形成实验检测细胞增殖和克隆形成能力;划痕实验检测细胞迁移; PI染色流式细胞术检测细胞周期;转录组测序检测USP1敲除后基因表达差异及其相关的功能与信号通路;双荧光素酶报告基因检测USP1对信号通路的影响,并通过免疫印迹技术进行验证。结果显示, USP1过表达可以促进膀胱癌细胞增殖,敲除后显著抑制了膀胱癌细胞的增殖和克隆形成及迁移能力,促进S期细胞阻滞。转录组测序结果显示, USP1敲除后差异表达基因共4 522个,其中上调基因2 078个,下调基因2 444个。KEGG分析结果显示,这些基因涉及多个方面,包括细胞周期调控、细胞信号转导、转录翻译、蛋白折叠降解、自噬凋亡等。Hallmark数据库分析结果显示,差异表达基因与E2F信号通路密切相关。双荧光素酶报告基因显示, USP1过表达后, E2F1信号通路明显上调并且呈剂量依赖式,免疫印迹结...     

ABCA1对滋养细胞基因表达谱的影响

该文从基因学角度探讨了ABCA1对人滋养细胞功能的影响。在人滋养细胞中调控ABCA1的表达,通过表达谱芯片检测ABCA1表达改变后滋养细胞中基因及相关信号通路的改变,Western blot及qRT-PCR验证芯片结果。结果显示,上调ABCA1表达时,有197个基因表达升高,有190个基因表达下降;而下调ABCA1表达时,有335个基因表达升高, 459个基因表达下降。GO和KEGG分析表明, ABCA1表达上升或下降后可导致滋养细胞内多个信号通路发生改变。qRT-PCR及Western blot检测发现,与阴性对照组相比, ABCA1上调后细胞中S1PR1的m RNA及蛋白水平明显升高,而CCL8、CXCL10与CXCL11的mRNA及蛋白水平明显降低;下调ABCA1的表达后,细胞中S1PR1的mRNA及蛋白水平明显降低,而CCL8、CXCL10与CXCL11的mRNA及蛋白水平明显升高,与芯片的结果完全一致。这些结果表明, ABCA1可通过调控多个信号通路而影响滋养细胞功能。     

沙葱萤叶甲卵滞育的转录组学分析

【目的】建立沙葱萤叶甲Galeruca daurica滞育卵转录组数据库,挖掘卵滞育相关的基因以及代谢和信号通路,在转录组水平探讨卵滞育的分子机制。【方法】采用Illumina NovaSeq6000高通量测序平台对沙葱萤叶甲滞育卵与解除滞育卵进行转录组测序,并进行生物信息学分析;利用DESeq软件分析沙葱萤叶甲滞育卵与解除滞育卵中的差异表达基因,对差异表达基因进行KEGG通路富集分析;利用qRT PCR技术对10个差异表达基因的表达模式进行验证。【结果】基于沙葱萤叶甲滞育卵与解除滞育卵转录组测序结果,共获得53 389个unigene,其中差异表达基因2 145个,24个差异表达基因与保幼激素信号及脂肪酸生物合成和降解相关。与解除滞育卵相比,滞育卵转录组中1 297个基因上调表达,富集于124条KEGG通路,其中核糖体通路显著富集;848个基因下调表达,富集于73条KEGG通路,其中MAPK信号通路和糖胺聚糖生物合成通路显著富集。qRT-PCR结果表明,随机选取的10个差异表达基因的表达趋势与RNA-Seq转录组测序结果完全一致。【结论】保幼激素,脂肪酸生物合成和降解,核糖体,MAPK信号及糖胺聚糖生物合成等通路可能在沙葱萤叶甲卵滞育调控中起着重要的作用。     

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

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

镧对铜胁迫下水稻转录组模式的调控分析

探究稀土镧对铜胁迫下水稻转录组的影响,鉴定镧在调控水稻铜胁迫应答中的关键基因和功能。利用抗氧化酶活性筛选最适铜胁迫及镧处理浓度,进行转录组测序和差异表达分析,并用qRT-PCR验证差异表达水平。对测序数据进行差异表达分析发现3 222个基因显著上调,3 798个显著下调,qRT-PCR验证了4个细胞壁防御相关基因CHIT13、Laccase、Expansin和GH3.4。功能和通路分析获得95组GO条目和112条KEGG通路,富集于激酶和氧分子结合等分子功能、细胞壁及质膜等细胞组分、次级代谢和脂质代谢等生物学过程以及苯丙烷生物合成和植物激素信号转导通路。镧通过调节细胞壁形成和组分提高水稻铜胁迫耐受性。     

基于转录组测序的卤虫卵生和卵胎生相关基因表达研究

为揭示卤虫(Artemia)不同繁殖模式的发生机制, 文章通过构建孤雌生殖卤虫卵生和卵胎生差异转录组文库并结合生物信息学分析, 对两种繁殖模式间的差异表达基因进行查找筛选, 然后利用qRT-PCR对候选繁殖模式相关基因的表达进行分析。转录组测序显示有1452个差异表达基因, 包括601个上调基因和851个下调基因。根据差异表达基因GO功能分类结果可知, 注释到生物过程、细胞组成和分子功能的unigene分别有1243、306和530个。KEGG富集分析结果显示差异基因显著富集在抗原加工和核糖体通路中。结合转录组分析, 进一步筛选得到6个生殖相关基因, 并针对不同繁殖模式下的卤虫进行qRT-PCR, 结果表明, 6个生殖相关基因在卵生卤虫卵巢中的表达量均显著高于卵胎生卤虫。此外, 对6个候选生殖相关基因编码的蛋白质保守结构域进行预测, 发现均与之前报道的相应基因保守结构域一致。综上所述, 研究所选择的6个基因可能影响参与了卤虫的生殖过程。研究结果为孤雌卤虫繁殖模式分子机制调控的研究提供了基础信息, 有助于完善卤虫的生殖生物学理论。     

基于转录组测序的羽衣甘蓝叶色相关基因分析

叶色是羽衣甘蓝重要的观赏性状之一。本研究采用Illumina Hi-Seq2500高通量测序技术,对基因型纯合的紫叶和白叶羽衣甘蓝叶片进行转录组测序,筛选差异基因并与GO和KEGG数据库比对进行注释分析,分析羽衣甘蓝叶色形成相关基因。结果显示,获得高质量短读序共104 608 770条,筛选出紫叶相对白叶的差异表达基因1 993个,其中上调表达基因1 094个,下调表达基因899个。根据GO功能分类可分为生物过程、细胞组分和分子功能3大类64功能组。根据KEGG代谢通路分析可以分为171类,在叶色相关的类黄酮生物合成途径中黄酮醇合成酶(FLS)、二氢黄酮醇4-还原酶(DFR)上调以及类胡萝卜素生物合成途径中的类胡萝卜素β-环化酶上调与紫叶形成关系密切。本研究丰富了羽衣甘蓝的转录组信息,获得了一些差异表达基因,为进一步研究羽衣甘蓝叶色形成的遗传机制提供了有价值的信息。     

lncRNA Tmevpg1对小鼠自噬和JAK-STAT信号通路关键信号分子表达水平的影响

该文主要研究lncRNA Tmevpg1对小鼠自噬和JAK-STAT信号通路关键信号分子等相关分子表达水平的影响。该研究利用生物信息学分析构建Tmevpg1、JAK-STAT信号通路及自噬互作调控网络;雷帕霉素(rapamycin, Rapa)和氯喹(chloroquine, CQ)构建小鼠自噬模型;通过细胞共培养技术构建Tmevpg1过表达与干扰细胞模型; qRT-PCR(quantitative real-time PCR)检测小鼠脾脏和细胞模型中Tmevpg1等相关分子的RNA水平, Western blot检测自噬相关蛋白、T-bet和JAKSTAT通路关键信号分子磷酸化蛋白表达水平。结果显示, Tmevpg1、IFN-γ、T-bet、STAT1和JAK1在自噬发生的一定时间点内显著上调(P0.001);过表达Tmevpg1后ULK1表达上调(P0.05), p62下调(P0.05), IFN-γ、T-bet、JAK1以及STAT1表达未发生明显改变,而干扰Tmevpg1后IFN-γ、T-bet、JAK1以及STAT1表达减少(P0.05); Western blot结果显示小鼠自噬模型构建成功, T-bet、p-STAT1和p-JAK1表达趋势与m RNA水平一致;过表达Tmevpg1可上调ULK1和LC3-Ⅱ,下调p62,同时干扰Tmevpg1后T-bet、p-JAK1和p-STAT1表达水平下降。该项研究结果表明, lncRNA Tmevpg1和JAKSTAT信号通路对细胞自噬发挥重要调控作用。     

不同海拔饲育的麦洼牦牛肺脏组织转录组学分析

为了探讨牦牛适应高海拔低氧环境的基因表达特征与规律,对在高海拔（3 560 m）和低海拔地区（478 m）饲育4个月的2.5~3岁健康雄性麦洼牦牛肺组织进行转录组测序。转录组测序采用Illumina高通量测序平台(HiSeqTM2500/4000)进行,并以qRT-PCR验证差异表达基因的表达量。结果显示,高海拔组牦牛肺脏转录组平均每个测序样本得到约5.76亿条Clean Reads,低海拔组牦牛中得到约6.10亿条Clean Reads,比对到参考基因组上的Reads数分别占91.74%和91.28%以上,共发现了2 047个新转录本。低海拔组与高海拔组牦牛肺脏组织之间共有199个差异表达基因,其中含89个差异上调表达基因和110个差异下调表达基因。所得差异表达基因富集在297个GO条目和146个KEGG通路中,包含62个低氧适应相关的GO条目和35个低氧适应相关代谢通路。其中低氧适应相关GO条目在生物过程、细胞组成和分子功能三种类别中占比最多的分别为细胞粘附、蛋白复合物和钙离子结合。低氧适应相关KEGG通路中占比最多的为肿瘤坏死因子(TNF)信号通路,其次为低氧诱导因子1(HIF-1)信号通路。qRT-PCR验证结果显示,Ⅱ类人类白细胞抗原α链(HLA-DOA、HLA-DRA)、补体因子 (C2)和甘露糖结合凝集素相关丝氨酸蛋白酶1（MASP1）基因的表达量变化与转录组测序结果相符。本研究为全局和深入理解牦牛肺组织转录本表达对高海拔低氧的响应提供了有价值的切入点。     

RNA-seq用于获得共轭亚油酸对贵妃鸡肌内脂肪代谢调控的关键基因

通过转录物组测序获得在贵妃鸡基础日粮中添加共轭亚油酸(CLA)对肌内脂肪代谢的差异表达基因,经生物信息学分析获得相关的信号通路及可能发挥重要作用的候选基因,为CLA对肌内脂肪沉积的分子机制奠定基础。本研究选用55日龄健康的贵妃鸡为试验动物,在基础日粮中添加CLA 0%、1%和2%,预饲期1周,正饲期6周。屠宰采集胸肌组织进行转录物组测序,对测序数据进行差异表达分析,差异表达基因GO功能和差异表达基因KEGG通路富集分析,筛选出与胸肌脂类代谢相关的差异表达基因,利用qRT-PCR对差异表达基因进行验证。结果显示,共获得1 065个差异表达基因,其中上调基因703个,下调基因362个。GO富集结果显示,差异表达基因主要富集在生物过程的细胞过程、单一生物过程、生物调节和代谢过程。KEGG信号通路富集显示,差异表达基因显著富集在黏着斑、不饱和脂肪酸生物合成、脂肪酸生物合成和类固醇生物合成等信号通路中,发现11个主要与肌内脂肪代谢相关的候选基因,分别是FADS1、FADS2、ELOVL5、ACOX2、SLC27A1、FABP5、LPL、LOC107050163、ENSGALG00000030996、ENSGALG00000005043和ENSGALG00000048882。并随机选取6个基因进行qRT-PCR验证,其相对表达量变化趋势与测序结果一致。本研究筛选到CLA影响贵妃鸡胸肌脂类代谢相关的差异表达基因,并对11个主要参与脂肪代谢相关的基因进行分析,为揭示CLA调控肌内脂肪沉积的分子机制奠定基础。     

茶树种子发育过程的转录组分析

该研究以‘铁观音’茶树品种的种子为试验材料,采用转录组测序技术分析种子发育的3个时期(幼果期、膨大期、成熟期)的表达差异,探究茶树种子油脂代谢的分子机制。结果表明:(1)经转录组测序、组装后共获得30 940 581个clean reads,经数据合并拼接最终得到36 951条非冗余Unigene序列,其中28 476个Unigene可得到功能注释;在转录本中能够被注释到GO分类的Unigene有11 201条(30.3%),KEGG分析发现共有17 172个基因参与了127个代谢通路。(2)经KEGG通路筛选出14条与脂肪酸代谢相关的通路,且随着茶籽的发育,大部分脂肪酸调控途径相关基因呈下调趋势,其中上调基因数最多的有α-亚麻酸代谢途径和脂肪酸降解途径(有17个基因表达量上调),下调基因数最多的是甘油磷脂代谢途径(有58个基因表达量下调);在茶籽发育幼果期α-亚麻酸代谢途径中表达量上调的基因数超过表达量下调的基因数。(3)研究发现茶籽脂肪酸合成相关的基因涉及14个脂类调控途径,共409条差异基因;随着茶树种子发育到成熟期,上调的差异表达基因数量在减少,下调的差异表达基因数量增加,其中α-亚麻酸途径中的基因PLA2G16、DAD1、pldA、FabF、FabI表达量上调显著,随后表达量下调。(4)qRT-PCR检测结果表明,7个茶树FAD和1个ACP差异表达基因的水平与转录组测序结果基本一致;随着茶籽的发育,基因CsFAD7和Δ6-CsFAD从幼果期、果实膨大期至果实成熟期都为差异下调表达,CsFAD2、CsFAD6和Δ7-CsFAD为差异上调表达,CsFAD8、Δ8-CsFAD和CsACP在幼果期至果实膨大期差异上调表达,在果实膨大期至果实成熟期差异下调表达。     

基于转录组测序的花烟草低温胁迫响应转录因子挖掘

为了鉴定参与花烟草低温胁迫的转录因子,对花烟草幼苗进行了4℃低温处理,并在处理后12 h采集其幼苗样本,提取总RNA后,采用高通量测序技术,进行了转录组测序。在对差异表达基因进行GO及KEGG分析的基础上,对参与其中的转录因子进行了挖掘,并采用qRT-PCR的方法,对转录组测序的结果进行了验证。结果表明,低温处理后,花烟草基因表达量变化在2倍以上的基因有8388个(P<0.01),其中,上调表达4229个,下调表达4159个。这些差异表达基因的功能归类于生物过程、细胞组分及分子功能3大类69个GO条目,并显著富集在40条KEGG代谢通路中。同时,在低温胁迫下,花烟草有118个转录因子的表达发生了显著改变,其中,上调表达82个,下调表达36个。这些转录因子属于28个家族,其中,数量最多的为NAC家族19个,其次为ERF家族16个、MYB家族15个、WRKY家族15个。本研究的结果为花烟草低温响应分子机制的研究提供了借鉴。     

果生刺盘孢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基因与内质网压力胁迫应答之间关联的新认识,为阐明果生刺盘孢响应内质网压力胁迫和致病机制奠定了基础。     

马立克病病毒pUL56蛋白对DF-1细胞免疫通路相关基因的转录调控研究

马立克病病毒（Marek’s disease virus,MDV）UL56基因编码α-疱疹病毒同源蛋白pUL56。为全面了解该病毒蛋白与宿主细胞的互作关系,本研究首先将表达红色荧光蛋白标记的MDV pUL56质粒转染鸡DF-1细胞系,瞬时表达蛋白的细胞经流式细胞术分选,富集阳性细胞群进行高通量测序及mRNA转录组分析。测序结果显示,表达pUL56的DF-1细胞存在914个差异表达基因,其中462个上调,452个下调。通过GO、KEGG和Reactome通路富集分析,发现富集的显著差异表达基因分布于先天性免疫反应、细胞因子产生等相关信号通路。在此基础上,随机选取了I型干扰素信号通路和促炎症细胞因子基因进行RT-qPCR验证,IFI6、IFIH1、CD83、HSP90AA1、IL-1β和IL-8L1基因的表达量均为上调,而IL-12β表达则下调,该结果与高通量测序相符。本研究提示pUL56作为非结构蛋白,可能对宿主细胞过程有广泛的调节作用。     

利用生物信息学研究Tgfβ信号通路在人类结直肠癌发生发展中的作用

目的:系统性研究经典Tgfβ信号通路各成员,在人类早期结直肠癌组织中,转录水平上的表达差异,以探索Tgfβ信号通路在人类结直肠癌发生中作用机制。方法:利用生物信息学手段,在GEO数据库中下载符合纳入标准的早期人类结直肠癌芯片结果,通过Ultral Edit软件实现探针号和基因名的转换以建立全转录本表达数据表;通过KEGG信号通路图谱选择Tgfβ信号经典通路待测分子,统计学分析各待测分子在正常组癌症组中的表达差异。结果:在GEO数据库中下载和获得符合纳入标准的早期人类结直肠癌结直肠癌芯片12组基因芯片结果,对Tgfβ信号通路及其旁路共27个进行表达差异分析,其中表达水平上调的分子4个,下调的基因3个,具有或接近显著性意义。结论:在直肠癌发生过程中,Smads分子的表达下降阻断经典的Tgfβ信号的抗肿瘤作用,同时配体分子和受体分子的广泛上调可能刺激旁路分子的激活从而导致肿瘤的发生。     

烟草脉斑驳病毒侵染本氏烟的转录组研究

为了研究烟草脉斑驳病毒（TVMV）侵染本氏烟的抗病分子机制,本文建立了TVMV转录组数据库,挖掘抗病相关基因以及代谢和信号通路,采用Illumina Novaseq 6000高通量测序平台对侵染TVMV的本氏烟进行转录组测序,并进行生物信息学分析;利用R package:edge R等软件分析了有关抗病的差异表达基因,对差异表达基因进行基因本体论（GO）及京都基因和基因组百科全书（KEGG）基于途径的通路富集分析。基于TVMV侵染本氏烟转录组的测序结果中共获得4 593个差异表达基因,其中3 564个基因上调表达,1 029个基因下调表达。共筛选出10个抗病相关的差异表达显著基因,6个上调,4个下调。GO数据库中注释到生物学过程、细胞组分及分子功能等三大类共50个功能组。其中,在第三大类分子功能中,蛋白质结合功能以及ATP结合功能所涉及的差异表达基因最显著,分别为501个和453个差异基因。KEGG共富集20条通路,注释到核糖体途径的差异基因富集程度最高,基因最显著,差异表达基因数为307个。蛋白质结合、ATP结合、核糖体、真核生物核糖体的生物反应、DNA复制（DNA replicat...     

piwil2重编程人成纤维细胞形成肿瘤样干细胞的基因差异表达

利用基因芯片技术检测piwil2基因重编程人成纤维细胞(fibroblast, FB)形成肿瘤样干细胞的差异表达基因,分析主要差异表达信号通路。采用前期已构建的piwil2-gfp基因,以gfp作为空载对照的慢病毒转染人包皮来源的成纤维细胞,获得稳定转染的肿瘤样干细胞,人成纤维细胞作为正常对照,采用基因芯片技术检测其转录差异表达,进行pathway分析,挑选差异表达明显的支链氨基酸代谢信号通路,采用q-PCR进行验证。基因芯片结果显示,FB-piwil2组相较于FB,上调基因1 490个,下调基因1 484个,参与支链氨基酸代谢等多条信号通路;FB-piwil2组相较于FB-GFP组,上调基因3 765个,下调基因1 988个,参与支链氨基酸代谢等多条信号通路;FB-GFP组相较于FB组,上调基因2 095个,下调基因4 014个,参与破骨细胞分化等信号通路;q-PCR验证重编程形成的肿瘤样干细胞克隆的支链氨基酸代谢信号通路上调的7个基因(acad8, acadm, aldh3a2, bcat1, hibch, hmgcs1, mccc2)与基因芯片上调表达一致。由此可见,piwil2基因重编程人成纤维细胞形成的肿瘤样干细胞基因出现差异表达,与支链氨基酸代谢相关。     

SNRPA1对肝癌的成瘤性的调控机制

小核核糖核蛋白(small nuclear ribonucleoprotein polypeptide A,SNRPA1)参与mRNA加工剪接体(splicesome)的组装,与多种肿瘤的发生有关,但在肝癌发生发展过程中的分子机制尚不明确。该研究利用基因芯片技术探究SNRPA1敲减后肝癌细胞信号通路关键基因的表达动态,及其在裸鼠肿瘤发生发展中的分子机制。研究结果显示,与对照组相比,SNRPA1基因敲减后的低表达组裸鼠的荷瘤细胞的生长受到显著抑制;基因芯片分析表明,SNRPA1的敲低导致462个基因的表达下调,262个基因的表达上调。qRT-PCR分析表明,FSTL1、FGF2、JAK2、WNT5A和PPM1A基因表达均有所降低,而Western blot分析进一步证实FSTL1、JAK2、WNT5A蛋白质表达的下调。以上结果可知,SNPRA1作为促癌基因,可能通过调控多种基因的表达及信号通路对肝细胞肝癌的发生和发展进行调控。     

阿尔泰蝠蛾(鳞翅目:蝙蝠蛾科)幼虫低温胁迫转录组分析

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾Hepialus altaicola Wang幼虫基因转录差异及上调表达基因的主要功能类群和参与的主要代谢通路.[方法]采用Illumina HiSeq TM 2500测序平台对阿尔泰蝠蛾幼虫进行转录组测序、组装,利用Blast软件进行数据库比对和基因功能注释,用DEGSeq R软件包分析4℃低温处理与室内适温饲养试虫的差异表达基因,并对上调表达基因进行GO和KEGG代谢途径富集分析.[结果]经序列拼接后共获得100300个unigenes,总长度81600309 bp,平均长度813 bp,N50长度1719 bp.与7大数据库同源比对,共获得34691(34.59％)条unigenes.低温胁迫转录组分析得到11569个差异表达基因(DEGs),7158条基因上调,4411条基因下调.富集到47个GO类群,217个KEGG途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.     

阿尔泰蝠蛾(鳞翅目:蝙蝠蛾科)幼虫低温胁迫转录组分析

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾Hepialus altaicola Wang幼虫基因转录差异及上调表达基因的主要功能类群和参与的主要代谢通路.[方法]采用Illumina HiSeq TM 2500测序平台对阿尔泰蝠蛾幼虫进行转录组测序、组装,利用Blast软件进行数据库比对和基因功能注释,用DEGSeq R软件包分析4℃低温处理与室内适温饲养试虫的差异表达基因,并对上调表达基因进行GO和KEGG代谢途径富集分析.[结果]经序列拼接后共获得100300个unigenes,总长度81600309 bp,平均长度813 bp,N50长度1719 bp.与7大数据库同源比对,共获得34691(34.59％)条unigenes.低温胁迫转录组分析得到11569个差异表达基因(DEGs),7158条基因上调,4411条基因下调.富集到47个GO类群,217个KEGG途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.