火龙果转录组测序、基因表达与功能分析

利用高通量测序技术对火龙果（Hylocereus undulatus Britt）红肉品种‘大红二号’的花芽、果实和枝条不同发育阶段的基因表达进行研究。结果显示,转录组测序共获得468.68 Gb原始数据（Raw data）,从头组装获得239 152条转录本和162 519条unigene,约53.74%的unigene得到注释。分别在43 506条和16 251条unigene中检测到600 283个SNP位点和56 147个SSR位点。基因表达分析结果表明,在火龙果不同组织Fl510、Fl513、Fl514、Fl518、F711、F715、S513、S419中分别有31、7、5、152、17、63、17、8个特异表达的unigene。通过GO和KEGG富集分析,发现了一些组织特异的GO条目和代谢通路,如在Fl510中富集的类萜骨架生物合成代谢通路等。本研究还对参与花发育的候选基因进行了鉴定和表达分析,他们包括COL基因、FT-like基因、分生组织决定基因和器官决定基因等。     

草地贪夜蛾雄性成虫和5龄幼虫的转录组比较分析

【目的】草地贪夜蛾Spodoptera frugiperda是一种新近入侵我国的重要害虫。本研究旨在对草地贪夜蛾雄性成虫和5龄幼虫两个不同发育阶段的转录组进行比较分析。【方法】利用高通量测序技术对草地贪夜蛾雄性成虫和5龄幼虫进行转录组测序和数据组装,并对转录组数据进行功能注释和比较分析。【结果】经de novo组装共获得209 002条转录本,平均长度为687.55 bp,N50为982 bp。共有46 198条(57.43%) unigene在至少一个数据库中获得功能注释,其中1 713条(2.13%) unigene在所有数据库中均能获得注释。在GO数据库中获得205 269条unigene的注释,主要包括68个功能分类;在KEGG数据库中共有3 408条unigene得到注释,涉及277个代谢通路。共鉴定到424个嗅觉相关的基因,并且在雄性成虫和5龄幼虫之间的表达存在差异。通过比较转录组分析,在雄性成虫中鉴定到9 162个上调和6 399个下调差异表达基因(DEGs);功能富集分析发现在上调DEGs中涉及信息素以及信号转导的代谢通路显著富集,而下调DEGs中涉及解毒相关的通路显著富集。【结论】这些转录组数据为探究草地贪夜蛾的生长发育、嗅觉相关功能基因以及候选分子靶标提供资源信息。     

基于高通量测序的文冠果转录组分析

应用Illumina Hi-seqTM2000高通量测序技术对文冠果花芽进行转录组分析。共获得N50为1 180bp、平均长度为686bp的unigene 58 311条。与公共数据库Nr和Swiss-Prot同源性比较后发现37 047条unigene获得基因注释,另有21 264条unigene未被注释。利用COG数据库将unigene分成25类。通过GO分类和KEGG Pathway富集性分析,将unigene分别归类于55个GO term和128个代谢途径。此外,在9 794条unigene中共搜索到12 213个SSR位点,单核苷酸重复基元出现频率最高(34.95%),其次分别为二核苷酸(32.74%)和三核苷酸(28.64%)。在获得的unigene中发掘出涉及4个开花调控途径(光周期途径、春化途径、GA途径和自主途径)多个基因的同源序列。研究结果可在一定程度上解析文冠果花芽形态分化的分子调控模式与机制。     

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

【目的】建立沙葱萤叶甲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信号及糖胺聚糖生物合成等通路可能在沙葱萤叶甲卵滞育调控中起着重要的作用。     

沙葱萤叶甲对蜕皮激素响应的转录组分析

为建立沙葱萤叶甲Gauleruca daurica成虫对蜕皮激素（20-hydroxyecdysone,20E）响应的转录组数据库,挖掘对20E响应的基因以及代谢和信号通路,并在转录组水平探讨20E调控生殖滞育的分子机制,本研究采用Illumina HiSeqTM4000高通量测序平台对20E及二甲基亚砜（dimethyl sulfoxide,DMSO）处理后的沙葱萤叶甲成虫进行了转录组测序,共获得80 313个unigene;与阴性对照DMSO相比,共获得201个差异表达基因,其中106个上调、95个下调。GO和KEGG富集分析表明,多数差异表达功能基因富集于各种代谢通路,其中核黄素代谢（riboflavin metabolism）、溶酶体（lysosome）和泛酸与乙酰辅酶A合成（pantothenate and CoA biosynthesis）通路显著富集（q<0.05）。结果表明,20E可能通过影响多种代谢通路调控沙葱萤叶甲的生殖滞育。     

苗期水稻响应褐飞虱取食的基因差异表达分析

【目的】褐飞虱Nilaparvata lugens (St?l)是水稻的重要害虫之一,主要刺吸水稻韧皮部汁液为害,对水稻生产造成严重的产量和经济损失。为研究水稻响应褐飞虱取食的分子机制,对褐飞虱取食6 h后的苗期水稻进行转录组测序及基因差异表达分析。【方法】采用illumina二代测序技术获得褐飞虱取食前后水稻组织的转录组数据,利用RSEM软件进行基因表达定量和DEseq2进行差异表达分析;从差异表达基因中随机选取20个基因采用荧光定量PCR技术进行验证;采用GeneMerge软件对差异表达基因进行KEGG和GO富集分析。【结果】褐飞虱取食后,水稻转录组中的1 104个基因出现了差异表达,其中435个基因表达上调,669个基因表达下调。荧光定量PCR结果显示,20个差异表达基因中18个基因的表达变化趋势和测序结果一致,证明了转录组分析结果可靠。GO和KEGG富集分析表明,表达上调基因主要与水稻氧化应激、海藻糖合成及次生化合物代谢有关,显著富集在14个KEGG通路和30个GO功能分类中;而表达下调基因主要参与水稻纤维素、蛋白质及脂肪酸合成过程,显著富集在29个KEGG通路和26个GO功能分类。在差异表达基因中,分别有61个转录因子和13个水杨酸和茉莉酸信号通路相关基因。【结论】褐飞虱取食激发了水稻的应激反应和保护机制,同时还降低了营养合成的过程,是飞虱为害造成水稻减产的原因之一。本研究初步揭示了苗期水稻响应褐飞虱取食的差异表达基因,为研究水稻-褐飞虱互作机制以及褐飞虱抗性水稻品种培育提供了参考和依据。     

夏季高温干旱逆境下梭梭枝条中转录组表达谱分析

本研究以荒漠地区极为抗旱、耐热的超旱生固沙植物梭梭为材料,利用Solexa测序技术对高温和适温下野生梭梭的枝条进转录组测序和数据de novo组装。将得到的unigene进行功能注释、分类及代谢通路分析。结果表明,获得长度大于200 bp差异表达基因162 504条,将得到的unigene与Nr和Swiss-Prot数据库比对后发现,分别有38 529条和23 198条unigene与其它物种的基因具有同源性。利用COG数据库可将有关的13 012条梭梭unigene分成25类,KEGG数据库分析发现共有22 249个unigene参与298种代谢通路。研究表明梭梭逆境胁迫反应是一个多基因参与、多个生物过程协同调控的过程,基因表达量的变化可能是调控的主要方式,为揭示梭梭耐旱机理以及梭梭耐旱性相关关键基因奠定了基础。     

携带与未携带松材线虫的松墨天牛转录组差异表达基因分析

【目的】本研究旨在建立携带与未携带松材线虫Bursaphelenchus xylophilus的松墨天牛Monochamus alternatus成虫各组织的转录组数据库,揭示松墨天牛对松材线虫响应的转录组整体表达特征。【方法】以携带和未携带松材线虫的松墨天牛成虫的表皮、气管和脂肪体为材料,采用Illumina HiSeq TM 2000测序平台开展转录组测序,利用Trinity软件对RNA-Seq数据进行从头组装,利用NCBI数据库进行基因注释。通过DEG-Seq分析携带和未携带松材线虫的松墨天牛成虫中的差异表达基因,对上调表达基因进行GO和KEGG代谢途径富集分析,并通过实时荧光定量PCR技术对部分差异表达的基因包括ENV,CDK1,HSP70-C,HSP70,HSP75,DUO,PABPN1和IGFP基因的表达水平进行验证。【结果】测序过滤后共获得55059条单基因簇(unigene),平均长度为1536 bp。将unigene与数据库中的序列进行BLASTX比对,成功注释24354条unigenes,其中4022条unigene在GO数据库中获得注释,6098条unigene在KEGG数据库中获得注释。经过GO与KEGG富集分析显示,松墨天牛成虫携带松材线虫后,差异表达基因主要为与其气管中压力调节、组织与DNA修复、激素响应方面相关的基因。实时荧光定量PCR分析结果显示CDK1,HSP70,HSP75,DUO,PABPN1和IGFP基因的表达量在携带松材线虫的松墨天牛气管中相比未携带松材线虫的松墨天牛气管中的显著上调。【结论】本研究初步阐明松墨天牛携带松材线虫后转录组的整体表达模式,为进一步研究松墨天牛的基因功能及对松材线虫的胁迫响应过程奠定了基础。     

紫鸭跖草根组织应答铜胁迫的转录组分析

紫鸭跖草是一种高耐铜的超积累植物,本研究首次应用RNA-Seq技术对其转录组进行分析。通过全长转录组分析组装了紫鸭跖草耐铜相关候选基因,通过转录组分析,共获得82 471条N50长度为2 299 bp的高质量unigene,为紫鸭跖草的进一步研究提供了丰富的数据。对照组(CK)、300 μmol/L胁迫组(CT1)和1 000 μmol/L胁迫组(CT2)的测序数据已保存在NCBI的SRA数据库中,登录号为SAMN11265427。CT1相比于CK,共有5 028条unigene在根组织中有显著性差异表达,约占全部unigene的6.10%,其中富集上调和下调的基因分别为3 138和1 890条;CT2相比于CT1,根中共有6 813个unigene富集差异表达,占所有unigene的8.26%。其中富集上调和下调的基因分别为2 555和4 258。随机选取10个基因进行qRT-PCR荧光定量分析,结果与Illumina测序数据一致,验证了基因数据差异表达的有效性。上述实验从分子水平上为研究铜胁迫下铜耐受的分子机制提供了理论依据。     

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

【目的】筛选柑橘木虱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验证及表达量分析,研究其在柑橘木虱应对球孢白僵菌侵染过程中的作用。     

不同大小规格墨西哥湾扇贝(Argopecten irradians concentricus)转录组分析及生长相关基因筛选

墨西哥湾扇贝(Argopecten irradians concentricus)是中国双壳贝类养殖的重要经济品种,具有较高的研究价值,然而,转录组学及基因组学信息的匮乏,严重制约了该贝的遗传改良进程。本研究基于Illumina Hi SeqTM2000平台,对5月龄大规格(平均体质量16.05 g)和小规格(平均体质量5.57 g)墨西哥湾扇贝进行转录组测序,经de nove组装,大、小规格墨西哥湾扇贝分别获得120 563和118 794条高质量unigenes,unigene的平均长度分别为944 bp和969 bp;按照FDR≤0.001且|log2Ratio|≥1的原则,共筛选出9 901个差异表达基因,其中,大规格组相对于小规格组,上调基因4 976个(50.26%),下调基因4 925个(49.74%);六大公共数据库Blast比对结果显示,差异表达基因在Nr数据库的注释率最高,其次是GO数据库;从已获功能注释的差异表达基因中筛选得到47个生长相关候选基因,涉及多个生长因子及其受体;GO功能富集结果显示,膜、L-氨基酸运输和糖基键水解酶活性分别是细胞组分、生物过程和分子功能富集unigene最多的条目;Pathway富集分析发现,3 699条差异表达基因被成功注释到252条代谢通路中,其中33条为显著性富集通路(p<0.05),富集度最高的是焦点粘附;大规格组中共鉴定出10 870个SSRs和242 551个SNPs,小规格组中共鉴定出10 857个SSRs和228 121个SNPs。研究结果为墨西哥湾扇贝分子标记批量开发、功能相关基因挖掘以及分子遗传育种提供了理论数据支撑。     

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

为了探讨牦牛适应高海拔低氧环境的基因表达特征与规律,对在高海拔（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）基因的表达量变化与转录组测序结果相符。本研究为全局和深入理解牦牛肺组织转录本表达对高海拔低氧的响应提供了有价值的切入点。     

侧孢短芽孢杆菌S2-31拮抗下细辛叶枯病菌转录组差异表达分析

通过分析侧孢短芽孢杆菌拮抗作用下叶枯病菌的转录组学特征,研究差异表达基因(DEGs)和代谢通路的富集情况,初步探索侧孢短芽孢杆菌拮抗叶枯病菌的分子机制。首先利用S2-31与叶枯病菌的对峙培养观察其拮抗作用,然后利用转录组测序探究侧孢短芽孢杆菌拮抗下和正常生长下的叶枯病菌的基因表达水平差异,并进行RT-qPCR验证,最后,利用相关数据库对DEGs进行注释和富集分析。对峙培养后,叶枯病菌菌丝出现皱缩、扭曲等现象。测序后共得到3681个DEGs,其中2224个相对上调,1457个相对下调。GO功能注释分析表明,上调和下调表达的DEGs均注释到生物过程8个亚群、细胞组分8个亚群和分子功能4个亚群。KEGG富集分析中,共有732个unigenes定位到115条生物学通路中,其中富集程度相对较高的通路有氨基糖和核苷糖代谢、糖基磷脂酰肌醇(GPI)锚定生物合成和过氧化物酶体;涉及差异基因较多的通路有淀粉和蔗糖代谢、剪接体和胞吞作用;上调表达的DEGs主要富集在代谢途径,下调表达的DEGs主要富集在遗传信息处理和细胞过程途径。从显著富集的通路中筛选出与菌丝生长发育相关的差异表达基因,结果表明侧孢短芽孢杆菌主要抑制病原菌菌丝过氧化物酶体的形成、胞吞过程和GPI锚定的合成等过程。差异表达基因的RT-qPCR验证结果与转录组测序的结果一致。在侧孢短芽孢杆菌S2-31的拮抗下抑制了叶枯病菌的生长,其转录组特征也发生显著变化,差异基因主要涉及代谢、细胞过程以及遗传信息处理等途径的相关通路。     

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

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾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途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.     

金钱鱼下丘脑转录组分析及其对雌激素在体注射的响应

为探究XX雌性和XY雄性金钱鱼(Scatophagus argus)下丘脑中的差异表达基因及雌激素(17β-Estradiol,E2)在体注射对XY雄鱼下丘脑中基因表达的影响,开展了转录组测序分析,包括测序数据质控、基因功能注释,差异基因筛选、鉴定和差异基因功能富集分析等,并利用实时定量PCR(Real-time quantitative PCR, qPCR)检测了金钱鱼下丘脑中相关基因的表达。金钱鱼下丘脑转录组共测得Clean reads 275833710,测序质量30(Q30)大于95%, GC含量值大于48%。其中, Ctrl-XX-H vs. Ctrl-XY-H组,共筛选和鉴定了91个差异表达基因,包括36个上调基因和55个下调基因。在Ctrl-XY-H vs. E2-XY-H组,筛选和鉴定了28个差异表达基因,包括11个上调基因和17个下调基因。GO和KEGG富集分析发现,差异表达基因主要显著富集在细胞、单生物过程,膜、膜组分,结合等生物学功能及泛醌和其他萜类-醌生物合成及类固醇激素生物合成通路等。转录组和qPCR结果表明, XX和XY金钱鱼下丘脑中prl、ccna1和hs...     

葡萄糖诱导绿色杜氏藻转录组及相关通路差异分析

绿色杜氏藻是研究耐盐机理的模式绿藻．葡萄糖不仅是营养物质,而且还是信号物质．目前,对绿色杜氏藻转录组、糖处理后差异表达基因和β-胡萝卜素生物合成途径关键基因表达还不清楚．本研究通过Illumina HiSeqTM 2000高通量测序,获得葡萄糖处理和未处理绿色杜氏藻转录组信息．利用P value值和差异倍数对样本进行差异表达分析,共111条转录本存在差异表达,3条为上调转录本,108条为下调转录本．利用RT-qPCR检验差异表达分析的准确性. 结果表明,转录本表达结果与转录组分析结果一致．GO功能富集结果表明,71条下调转录本与代谢相关,占所有下调转录本的65.74%．KEGG富集分析结果表明,21条KEGG通路含89条下调转录本,14条通路与代谢相关．代谢中通路最多的为能量代谢（6条）,含63条下调转录本．能量代谢中与光合作用相关的下调转录本最多,为29条．通过分析找到2条与β-胡萝卜素生物合成相关通路（MVA/MEP途径及β-胡萝卜素合成途径）,并发现通路的关键基因hmgs、dxs、dxr、psy、pds、chyb,对其进行差异表达分析,均不存在差异表达．研究表明,葡萄糖抑制了绿色杜氏藻光合作用,代谢受阻,但未影响β-胡萝卜素生物合成相关通路及关键基因．     

五指山猪和长白猪背最长肌差异基因的筛选与注释

为了筛选五指山猪和长白猪背最长肌组织差异表达基因,本研究采用RNA-seq技术和生物信息学方法对6月龄、8月龄五指山猪和长白猪的背最长肌进行转录组测序分析,并对差异表达基因进行GO和KEGG Pathway显著性富集分析。通过与猪基因组比对分析后发现,在6月龄和8月龄五指山猪和长白猪均有约44 000 000条Clean reads能够比对到猪基因组上相关基因,每个文库均获得约18 200个表达基因,而转录本的数据为25 000个。在6月龄五指山猪与长白猪中共有30个基因差异表达,在8月龄五指山猪与长白猪中获得29个差异表达的基因。这些差异表达的基因主要富集在与糖酵解代谢、MAPK信号代谢通路和胰岛素信号通路等肌肉发育信号通路,通过筛选获得了与肌肉生长发育有关的候选基因,为探索五指山猪肌肉发育的分子机制提供了理论依据。     

桑树桑黄菌丝体和子实体基因差异表达分析

通过对桑树桑黄Sanghuangporus sanghuang菌丝体和子实体2个不同生长阶段的转录组进行分析,为研究桑黄子实体生长发育相关机制奠定基础。采用Illumina测序技术,对桑树桑黄菌株S23菌丝体和子实体2个不同生长发育阶段进行了全转录组测序。将转录组测序reads比对到参考序列上,菌丝体测序样本的reads比对率为82.89%;子实体测序样本的reads比对率为83%。基因差异表达分析显示,与菌丝体相比,子实体中显著上调表达基因为2 898个,显著下调表达基因为1 965个。经过Blast nr比对发现,桑黄菌在子实体阶段表达量上升的基因主要与各种氧化酶活性、疏水蛋白等相关;表达量下降的基因主要与糖类、氨基酸结合、运输等相关。基因本体（gene ontology,GO）富集分析表明,菌丝体及子实体两个阶段与跨膜转运相关的差异表达基因富集明显。代谢通路（pathway）富集分析表明,类固醇生物合成、精氨酸生物合成、丝裂原活化蛋白激酶（mitogen-activated protein kinase,MAPK）信号通路等差异基因富集明显。     

多裂骆驼蓬叶片转录组分析

多裂骆驼蓬为西北荒漠地区常见植物,具有抗风固沙、防止水土流失、抑菌杀虫和抗肿瘤等功效。为了增加骆驼蓬属植物开发利用的强度,弥补其基因功能、代谢通路等分子生物学研究层面的空缺,该文利用Illumina高通量测序平台对多裂骆驼蓬叶片进行转录组测序,根据测序结果进行转录组数据拼接、功能注释、序列水平和表达水平等分析。结果表明:共获得7 723 653 900 bp核苷酸序列信息,拼接组装得到78 641条Unigene,预测出55 535条CDS序列。与多个数据库对比后,获得33 184个NR数据库注释信息;31 835个GO数据库注释信息; 17 206个KOG数据库Unigene功能分类信息; 7 617个KEGG数据库代谢通路注释信息。对单核苷酸多态性位点和微卫星信息进行检测分析,共发现86 113个SNP位点,6 987个SSR信息。该研究首次获得并分析了多裂骆驼蓬的转录组数据,通过基因比对、CDS预测、通路注释、SNP检测和SSR检测等方法,对该植物的基因、通路以及分子标记等方面有了初步的认识,为本种植物后续研究及资源的开发利用奠定了基础。