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

为了鉴定参与花烟草低温胁迫的转录因子,对花烟草幼苗进行了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个。本研究的结果为花烟草低温响应分子机制的研究提供了借鉴。     

黑须污蝇不同发育阶段转录组及嗅觉相关基因分析

【目的】建立黑须污蝇Wohlfahrtia magnifica转录组数据库,挖掘黑须污蝇基因数据,为更深入地研究双峰驼阴道蝇蛆病提供理论依据。【方法】采用高通量测序平台Illumina HiseqTM4000对黑须污蝇幼虫、蛹和成虫进行转录组测序,并进行生物信息学分析。利用黑须污蝇幼虫、蛹和成虫转录组数据,通过基因差异表达分析以及GO功能显著性富集分析的方法筛选出嗅觉相关基因,并通过荧光定量PCR技术对黑须污蝇幼虫、蛹和成虫中OBP99b,OBP56a和OBP99a 3个嗅觉相关基因表达水平进行验证。【结果】结果显示,每个黑须污蝇样本的转录组数据量在4.96 Gb以上,G+C含量在35.35%以上,Q20含量在97%以上。与NCBInr, GO, KEGG, Pfam, Swiss-Prot和eggNOG数据库进行对比,共注释到73 303条unigenes。其中eggNOG数据库注释到的unigenes最多,为35 066条,分布在23个蛋白功能中,其中的翻译修饰、蛋白质周转的unigenes所占比例较大;GO数据库注释到的unigenes数为29 193条,功能包括分子功能、细胞组分、生物学过程三大类50个分支,其中参与生物学过程和氧化还原过程的unigenes比例较大;KEGG数据库注释到的unigenes数为15 068条,其中参与信号转导过程的unigenes比例较大。上述数据表明该转录组测序结果较好,为后续研究奠定了基础。依据黑须污蝇幼虫、蛹、成虫转录组数据筛选到30个嗅觉相关基因,其中包含9个气味结合蛋白(OBP)基因,进一步基因注释发现OBP99b,OBP56a和OBP99a基因在黑须污蝇不同发育阶段表达差异显著(|log2FC|>1,P<0.05),这在荧光定量PCR分析结果中得到了验证。【结论】本研究获得了黑须污蝇幼虫、蛹和成虫转录组数据,筛选出黑须污蝇各发育阶段嗅觉相关基因,并验证了OBP99b,OBP56a和OBP99a基因在黑须污蝇幼虫、蛹和成虫3个发育阶段差异表达,为防治双峰驼阴道蝇蛆病提供了新思路。     

脂多糖对鲤肠上皮细胞转录组模式的调控分析

旨为探究脂多糖对鲤肠上皮细胞相关基因及功能的影响.本研究以鲤肠上皮细胞为试验对象,正常处理为对照组,脂多糖处理为实验组,对处理24 h的两组鲤肠上皮细胞进行转录组测序.转录组测序结果共获得44.77G高质量数据,其中,近81.83％的数据能够比对到鲤基因组.利用DESeq软件分析,与对照组相比,脂多糖处理可诱导产生差异...     

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

利用高通量测序技术对火龙果(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基因、分生组织决定基因和器官决定基因等。     

多杀性巴氏杆菌攻毒小鼠肺脏的转录组分析

多杀性巴氏杆菌(Pasteurella multocida)可感染多种动物宿主并引发败血症和呼吸道疾病,危害动物健康,给养殖业带来严重经济损失。本研究为探究多杀性巴氏杆菌感染期间对宿主肺脏基因表达的影响,选取羊源A型多杀性巴氏杆菌对小鼠(Mus musculus)进行攻毒实验,攻毒后72 h采集小鼠肺脏组织并提取总RNA。基于转录组测序,与对照组相比,攻毒组获得2 443个差异表达基因,其中有1 437个基因上调,1 006个基因下调。对2 443个差异表达基因进行GO富集分析,结果显示,差异基因显著富集的GO terms主要有B细胞稳态、T细胞迁移的正调控、整合素复合物和胶原蛋白结合;KEGG富集分析结果显示,差异基因显著富集于11条免疫炎症信号通路,包括细胞因子与细胞因子受体互作、趋化因子信号通路等;免疫炎症相关信号通路中,差异基因经蛋白质-蛋白质互作分析,结果显示,C3、C3ar1、C5ar1、Cxcl10、Cxcr2和Gng11基因处于网络中心,说明这些基因在多杀性巴氏杆菌感染过程中发挥了重要作用;从11条免疫炎症通路的基因中挑选10个进行实时荧光定量PCR验证,发现有8个基因...     

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

利用高通量测序技术对火龙果（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基因、分生组织决定基因和器官决定基因等。     

水分胁迫下不同抗旱性小麦品种叶片转录因子表达差异研究

通过研究不同抗旱性小麦品种中转录因子表达水平的差异,为阐明小麦抗旱机制奠定基础。依据候选基因序列设计PCR引物,以干旱胁迫后0、3、6、9、12和24 h的小麦叶片为实验材料,以26S rRNA为内参,运用荧光定量PCR技术,检测Wdreb2、Wlip19基因在干旱敏感性和干旱耐受性小麦叶片中的相对表达量。定量PCR结果显示：干旱胁迫后,Wdreb2、Wlip19基因在干旱敏感性小麦叶片中的表达明显低于干旱耐受性小麦,在不同品种叶片中的响应时间和表达趋势存在差异。研究认为,Wdreb2、Wlip19基因在不同品种小麦受到干旱胁迫后的表达差异,与该品种小麦的抗旱能力具有一定的相关性。     

里氏木霉不同糖苷水解酶基因转录水平比较分析

以前期里氏木霉RNA-seq中发现的7个糖苷水解酶基因为对象,分析其不同条件下的表达特性,以期为寻找新的纤维素降解功能酶提供证据。运用生物信息学方法,分析了7个基因可能的编码产物和结构特征。以不同的产纤维素酶菌株（QM 9414、RUT C30）为材料,采用实时荧光定量PCR,对7个糖苷水解酶基因（编号4–10）在各种碳源条件下转录情况与主要的3个纤维素酶基因cbh1,cbh2,egl1（编号1–3）进行了比较分析。信息学分析表明,7个基因编码蛋白分属于GH47（4号、5号）,GH92（6–8号）,GH16（9号）,GH31（10号）糖苷水解酶家族,具有典型的信号肽序列。cbh1,cbh2,egl1基因在纤维素酶诱导条件下,转录水平均表现显著的增加,上调倍数以QM 9414菌株表现的最高。QM 9414菌株中,cbh1,cbh2,egl1基因在纤维素条件下的上调倍数显著高于乳糖,3个基因在RUT C30菌株中的转录水平则显示乳糖条件下上调幅度更大。7个糖苷水解酶基因也存在类似的情况,而且编码α-甘露糖苷酶和内切β-葡聚糖酶的8号、9号基因上调倍数在纤维素酶诱导条件下仅次于纤维素酶基因,而以甘油为碳源条件下,8号、9号基因上调倍数高于纤维素酶基因。4号基因在上述碳源条件下,转录水平变化不大。结果表明：4号基因可能是组成型表达。基因5、6、7、8、9、10的表达呈现明显的菌株和碳源依赖性,且在纤维素酶诱导条件下基本上是和3个纤维素酶基因共转录的。     

HaSNPV的9个基因的转录谱分析

使用反转录和实时荧光定量PCR技术,我们对HaSNPv的几个预期的早期基因、早晚期基因、晚期基因、极晚期基因的转录时相进行分析,结果表明这些基因的起始转录时间与其自身的启动子类型基本是一致的.但是预期的早期基因pkip晚期才开始转录;预期的早晚期基因ha107早期不转录,仅在晚期转录.这些基因的转录水平一般都在病毒感染细胞72 h后达到最高,并且极晚期基因polyhedrin的转录水平明显高于其它基因.Iap2的转录水平仅次于polyhedrin,表明它可能是一个功能基因.与AcMNPV的p10不同,在HaSNPV/HzAM1系统中p10的转录水平并不高.     

基于转录组的不同火龙果品种抗性差异分析

不同的品种抗性不同,为进一步探究不同火龙果品种之间的抗性差异,为后续火龙果抗性育种提供参考,该研究利用Illumina HiSeq 2000测序平台对'普通白肉'(BR)和'厄瓜多尔黄龙'(EY)两个品种进行转录组测序分析,并参考GO Ontology、KEGG等公共数据库对差异表达基因进行功能分类与富集分析.结果表明...     

大花铁线莲单重被花的转录组差异分析

以即可开单被花又可开重被花的大花铁线莲重被品种"薇安"为研究对象,对"薇安"同一时期同一植株上的3种不同花瓣类型（单被、半重被、重被）采用高通量测序技术进行拼接及功能注释,筛选不同花被类型下表达量高度差异的关键基因后进行实时荧光定量PCR验证。结果显示：转录组测序共产生13.8GB原始数据,3个转录本文库两两比较共获得3 075条差异表达基因（DEG）,其中单被花与半重被样本对比（A vs B）包括649条上调DEG,605条下调DEG;半重被花与重被花对比（B vs C）包括上调DEG 1 046条和下调DEG 721条;单被花与重被花对比（A vs C）有上调DEG 1 129条和下调DEG 859条。3个不同花被下共存的差异表达基因有134条。根据基因功能注释从总DEG中筛选出26条可能与重瓣化性状相关的基因进行聚类分析,并随机挑选10个目的基因进行荧光定量PCR验证。PCR结果显示这些基因的表达量在铁线莲"薇安"同一时期同一植株的3种花被类型中均存在显著差异。最终筛选出与铁线莲重瓣化相关的关键基因有MADS-BOX类基因PMADS1、AP3、FRUITFULL、FLC;生长素反应蛋白IAA9、生长素输入载体、脱落酸8'羟化酶、吲哚乙酸诱导蛋白ARG7等。本研究为探究铁线莲重瓣花分子机制提供了基础数据和理论依据。     

银杏细胞转录组高通量测序及分析

利用Illumina的Genome Analyzer Ⅱx对银杏(Ginkgo Biloba)细胞转录组进行高通量测序,挖掘银杏内酯和紫杉醇生物合成基因,特别是新的羟基化酶基因,为今后最终完善红豆杉细胞紫杉醇生物合成途径中未知的羟基化步骤作准备.通过测序,获得了银杏细胞69 286个contig,56 387个scaffold,32 032个unigene.Unigene平均长度636bp.另外从gap分布、GC含量、基因组coverage等方面对unigene进行评估,数据显示测序质量好,可信度高.通过分析unigene的表达和功能注释等信息,发现66个属于CYP450基因家族,726个参与次生代谢物合成,其中59个与萜类合成有关,17个与二萜类合成相关.利用生物信息学方法从Michigan State University银杏成熟叶、侧根、成熟果实、无菌苗以及次生茎的转录组数据中找到了与银杏细胞CYP450高度同源的紫杉烷羟基化酶候选基因15个,为后续研究奠定了基础.     

单细胞转录组高通量测序分析新进展

细胞异质性是生物组织的普遍特征。常规转录组测序(RNA-Seq)技术需要上万个细胞,所测结果实际上是一群细胞基因表达的平均值,所以难以鉴别细胞之间基因表达的异质性。单细胞RNA-Seq技术的分辨率精确至单个细胞,为辨别异质性群体中各种细胞类型的转录组特征提供了有力的工具。近年来单细胞RNA-Seq技术发展迅速,在方法学上包括cDNA扩增方法的多样化、对灵敏度和技术噪声的定量分析、浅覆盖高通量单细胞RNA-Seq方法和原位RNA-Seq技术等;在技术应用方面应用范围从早期胚胎发育扩大到组织器官发育、免疫和肿瘤等多个领域。文章对单细胞RNA-Seq在方法学和技术应用两方面的研究进展进行了详细阐述。     

基于比较转录组分析海洋弧菌X511的褐藻胶代谢途径

【目的】从北海涠洲岛海域腐烂的马尾藻中分离得到的海洋弧菌(Vibrio X511)具有较强的利用褐藻胶能力,本文利用转录组测序的方法以研究弧菌X511的褐藻胶代谢途径。【方法】采用Illumina Hi Seq2500测序平台对菌株在褐藻胶及葡萄糖培养下的转录组进行测序;比较和分析差异转录本,利用荧光定量PCR验证测序结果;采用GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)对差异转录本进行功能和Pathway注释。【结果】经比较发现,菌株在褐藻胶培养下相对于葡萄糖的培养共有2024个差异表达基因,其中1066个基因上调,958个基因下调;某些普遍存在于代谢途径中的基因在不同培养条件下也存在差异表达;海洋弧菌X511中涉及褐藻胶利用的所有基因以及合成乙醇的关键基因其转录量均有一定程度的上调;此外,通过分析发现该菌株具有独特的褐藻胶利用方式,其中的一个代谢过程尚未在弧菌中被报道。【结论】成功解析了海洋弧菌X511的褐藻胶代谢途径,丰富了生物方法降解褐藻胶的研究,为大型海藻生物质能源的研究提供有价值的数据支持。     

基于高通量测序的麦红吸浆虫转录组分析

麦红吸浆虫Sitodiplosis mosellana(Gehin)是一种世界性的小麦害虫。为获得其转录组信息,本研究采用新一代高通量测序技术Illumina HiSeq TM 2000对麦红吸浆虫成虫转录组进行测序。共获得转录组样本数据量为27.88 G,经分析共获得59257个Unigenes,总长度49861164 bp,最短20 bp,最长29282 bp,平均长度841 bp。将Unigenes序列与NR、NT、Swiss-Prot、KEGG、GO和KOG数据库进行比对(e≤10-10),共获得95029个结果。通过GO功能分类,共有19584个Unigenes在GO数据库中细胞组分、分子功能和生物学过程等3大类50个功能组中找到对应。与KOG数据库进行比对,共有11279个麦红吸浆虫Unigenes被注释,按功能大致可分为26类。通过KEGG pathways分析,共有9110个麦红吸浆虫Unigenes被注释,分别归属于细胞进程、环境信息进程、遗传信息进程、新陈代谢和有机体系统5大类代谢途径,主要包括细胞生长与死亡、细胞运动、信号转导、能量代谢等32类代谢途径。CDS预测发现30088条序列可被编码,占全部基因的50.78%。SSR位点查找发现,在59257个Unigenes中共找到36323个SSR位点,发生率为61.30%。本研究获得的巨大的麦红吸浆虫转录组信息,为麦红吸浆虫的功能基因挖掘提供了重要的信息资源。     

不同水分条件下杜鹃花转录因子的转录组分析

为了解杜鹃花(Rhododendron pulchurum)转录因子在干旱胁迫下的表达模式,利用高通量测序技术,对不同水分条件下"白凤4号"叶的转录因子表达谱进行分析。结果表明,在不同水分处理间有34~161个差异表达的转录因子,杜鹃花响应不同水分的转录调控主要是通过ERF、b HLH和MYB基因的表达协同完成的;干旱胁迫时,特异调动了NAC的差异表达和WRKY、b ZIP、PLATZ的上调表达,干旱后复水时特异调动了GATA表达来调控。RT-q PCR验证结果表明,基因表达趋势与测序结果一致,证明了测序结果的有效性。这为明晰杜鹃花抗旱的分子机理及分子育种奠定了理论基础。     

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

应用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途径和自主途径)多个基因的同源序列。研究结果可在一定程度上解析文冠果花芽形态分化的分子调控模式与机制。     

基于高通量测序的玉米象转录组分析

玉米象Sitophilus zeamais是一种世界性的储粮害虫,但其基因信息尚不完善.本研究利用高通量测序平台Illumina HiSeq TM2000对玉米象成虫进行转录组测序,总计获得64358条unigenes,总长度39481752 bp,最短201 bp,最长29046 bp,平均长度613 bp.将unigenes序列在7大数据库Nr、GO、Swissprot、KOG、Pfam、KEGG、TrEMBL中进行注释,分别有25271、20747、16477、12060、9471、3370、25500条unigenes获得注释.通过GO功能分类,共有20747条unigenes在GO数据库中生物学过程、细胞组分、分子功能3大类68个亚类功能组中找到对应.KOG结果显示,12606条unigenes归到25个基因家族.通过KEGG代谢通路分析,共有3370条unigenes被注释,分别归属于细胞进程、环境信息进程、遗传信息进程、新陈代谢和有机体系统5大类代谢途径,共33个亚类274个功能通路.进一步的数据分析,共鉴定出146081个SNP位点和5002个简单序列重复(SSR)位点.本研究获得的玉米象转录组信息,为玉米象的功能基因挖掘提供了重要的信息资源.