芋全长转录组测序分析及淀粉生物合成相关基因挖掘

【目的】获取芋(Colocasia esculenta)全长转录本序列和结构信息,并挖掘淀粉生物合成相关基因和转录本的结构信息,为后续阐明芋淀粉生物合成的分子机制及深入利用芋基因资源提供科学依据。【方法】以荔浦芋1号为试材,采集其3月龄植株的不同组织(叶片、叶柄、球茎、匍匐茎和根)开展混合样三代全长转录组测序,对获得的转录本进行功能注释,通过生物信息学方法分析转录本可变剪接(AS)、可变多聚腺苷酸化(APA)、长链非编码RNA(lncRNA)、转录因子(TF)和简单重复序列(SSR)等的结构信息,通过京都基因和基因组百科全书(KEGG)注释获得芋淀粉生物合成相关的基因和转录本,并对其进行AS和APA分析。【结果】通过对芋叶片、叶柄、球茎、匍匐茎和根混合样开展三代全长转录组测序,共获得全长转录本序列38 043条;通过与参考基因组进行比对,鉴定出新转录本31 058条,其中28 785条获得功能注释;通过对转录本结构分析,共鉴定出9360个AS事件、6436个基因存在poly(A)位点、25 081个SSR位点、304个lncRNA、1911个融合转录本和1608个TF。通过KEGG注释...     

福清山羊与努比亚黑山羊发情期卵巢组织RNA-Seq分析

【目的】研究福清山羊与努比亚黑山羊发情期卵巢组织转录组差异表达水平,一方面为山羊繁殖性状形成的分子机制提供理论依据,另一方面为利用努比亚黑山羊杂交改良福清山羊提供可加快遗传进展的分子标记。【方法】利用转录组测序方法对福清山羊和努比亚黑山羊发情期卵巢组织进行研究,筛选品种间的差异表达基因(differentially expressed genes,DEGs),并对DEGs功能进行注释和若干基因荧光定量PCR(quantitative real-time PCR,qRT-PCR)验证;同时通过与参考基因组比对,分析和筛选测序样品中存在的SNP/InDel。【结果】6个样品共得到46.68Gb Clean Data,DESeq分析发现了福清山羊和努比亚黑山羊发情期卵巢组织的DEGs149个(包含25个新转录本),其中表达上调53个,表达下调96个;初步认为输卵管素基因(oviductin,OVN)、类固醇合成快速调节蛋白基因(steroidogenic acute regulatory protein,STAR)、早期生长应答1基因(early growth response 1,EGR1)可作为福清山羊繁殖性能的候选基因。149个DEGs中的108个基因能被GO(gene ontology)数据库注释,30个DEGs能够被COG(Cluster of Orthologous Groups of proteins)数据库注释,91个DEGs能够被KEGG(kyoto encyclopediaof genes and genomes)数据库注释。KEGG通路分析表明DEGs共富集到21条信号通路中,6条通路显著富集。经过进一步的与参考基因组序列比对分析,共发掘1 506个新转录本。经qRT-PCR验证,所选基因(转录本)表达变化模式与转录组测序结果一致,表明测序结果可靠。【结论】在转录组水平上筛选出了福清山羊和努比亚黑山羊发情期卵巢组织的149个DEGs,发掘了1 506个新转录本,初步揭示了OVN、STAR和EGR1在山羊繁殖过程中可能发挥重要作用,为进一步探索山羊繁殖性状相关机理提供参考依据。     

福清山羊与努比亚黑山羊发情期卵巢组织RNA-Seq分析

【目的】 研究福清山羊与努比亚黑山羊发情期卵巢组织转录组差异表达水平,一方面为山羊繁殖性状形成的分子机制提供理论依据,另一方面为利用努比亚黑山羊杂交改良福清山羊提供可加快遗传进展的分子标记。【方法】 利用转录组测序方法对福清山羊和努比亚黑山羊发情期卵巢组织进行研究,筛选品种间的差异表达基因（differentially expressed genes,DEGs）,并对DEGs功能进行注释和若干基因荧光定量PCR（quantitative real-time PCR,qRT-PCR）验证;同时通过与参考基因组比对,分析和筛选测序样品中存在的SNP/InDel。【结果】 6个样品共得到46.68Gb Clean Data,DESeq分析发现了福清山羊和努比亚黑山羊发情期卵巢组织的DEGs149个（包含25个新转录本）,其中表达上调53个,表达下调96个;初步认为输卵管素基因（oviductin,OVN）、类固醇合成快速调节蛋白基因(steroidogenic acute regulatory protein,STAR)、早期生长应答1基因(early growth response 1,EGR1)可作为福清山羊繁殖性能的候选基因。149个DEGs中的108个基因能被GO（gene ontology）数据库注释,30个DEGs能够被COG（Cluster of Orthologous Groups of proteins）数据库注释,91个DEGs能够被KEGG（kyoto encyclopediaof genes and genomes）数据库注释。KEGG通路分析表明DEGs共富集到21条信号通路中,6条通路显著富集。经过进一步的与参考基因组序列比对分析,共发掘1 506个新转录本。经qRT-PCR验证,所选基因（转录本）表达变化模式与转录组测序结果一致,表明测序结果可靠。【结论】 在转录组水平上筛选出了福清山羊和努比亚黑山羊发情期卵巢组织的149个DEGs,发掘了1 506个新转录本,初步揭示了OVN、STAR和EGR1在山羊繁殖过程中可能发挥重要作用,为进一步探索山羊繁殖性状相关机理提供参考依据。     

粗山羊草（Ae.tauschii）幼苗和根全长cDNA文库构建及其EST注释与比较分析

 【目的】构建小麦D基因组供体粗山羊草根和叶的全长cDNA文库，探讨利用生物信息学方法分析发掘组织特异性表达基因的可行性。【方法】本研究利用Cap-trapper方法构建全长cDNA文库，利用高通量测序技术获得大批高质量EST序列，在基于Linux系统的本地化生物信息学分析平台上对上述EST进行了电子注释分析。【结果】成功构建了粗山羊草根和叶的全长cDNA文库，测序获得根EST序列6 973条和叶EST 6 616条。利用本地化数据分析平台对其进行功能注释。利用GO-Diff软件，发现了两个文库间表达基因的功能差异，找到组织间表达差异显著的基因和组织内特异性表达的基因。【结论】利用构建的全长cDNA测序获得大量EST序列，通过生物信息学方法挖掘差异表达或特异性表达基因的方法是可行的。     

中华蜜蜂细胞吞噬与包囊作用相关基因全长转录本鉴定及分析

【目的】系统鉴定和分析中华蜜蜂(Apis cerana cerana)吞噬与包囊作用相关基因和全长转录本,为深入开展相关基因和剪接体的功能研究奠定基础。【方法】基于前期已获得的高质量中华蜜蜂纳米孔长读段测序数据,通过Blast工具将全长转录本比对Nr数据库筛选出吞噬与包囊作用相关基因和全长转录本。利用gffcompare软件将全长转录本与东方蜜蜂(Apis cerana)参考基因组上注释的转录本进行比较,鉴定未注释的新基因和新转录本。利用TAPIS pipeline预测和分析吞噬与包囊作用相关基因的可变多聚腺苷酸化(alternative polyadenylation,APA)位点,并通过TBtools软件鉴定APA位点上游的基序(motif)。使用Astalavista软件鉴定可变剪接(alternative splicing,AS)事件,并通过IGV浏览器进行结构可视化。通过RT-PCR验证AS事件的真实性。【结果】共鉴定到中华蜜蜂吞噬与包囊作用相关的基因66个和全长转录本395条,发掘出东方蜜蜂参考基因组未注释的2个新基因和303条新转录本。对参考基因组已注释的34个基因进行了结构优化,分别延伸了18个基因的5′端和12个基因的3′端,同时延长了4个基因的5′端和3′端。共鉴定到含有1个及以上APA位点的吞噬与包囊作用相关基因47个,其中多于5个APA位点的基因最多,为32个。在APA位点上游鉴定到多个基序,一致性序列为：GRBGCNKSDAACAAYTRBGCBMRNGGBYAYTAYWCNVWNGG。共鉴定到吞噬与包囊作用相关基因的AS事件296次,其中包括131次可变3′端剪接(alternative 3′splice site,A3SS)、85次内含子保留(intron retention,IR)、70次可变5′端剪接(alternative 5′splice site,A5SS)和10次外显子跳跃(exon skipping, ES)。RT-PCR结果显示,扩增的目的片段大小符合预期,证实了随机选择的2次AS事件的真实性。【结论】系统鉴定了中华蜜蜂吞噬与包囊作用相关基因和全长转录本以及AS事件和APA位点,优化了东方蜜蜂参考基因组注释的吞噬与包囊作用相关基因的结构。     

中华蜜蜂泛素基因及其全长转录本的发掘及分析

旨在利用三代Nanopore测序技术发掘中华蜜蜂（Apis cerana cerana）泛素（ubiquitin）基因及其全长转录本。基于前期获得的中华蜜蜂工蜂4～6日龄幼虫肠道的全长转录组数据,使用BLAST工具将全长转录本的序列比对到KEGG和Nr数据库以鉴定泛素基因及其全长转录本。采用Astalavista软件分析泛素基因的可变剪接（alternative splicing,AS）事件,并通过RT-PCR加以验证。通过TAPIS pipeline软件分析泛素基因的可变多聚腺苷酸化（alternative polyadenylation,APA）位点。共鉴定到48个中华蜜蜂泛素基因和435条泛素基因相关全长转录本。共发掘到48个泛素基因的74次AS事件,包括31次内含子保留事件,19次可变3′端剪接事件,16次可变5′端剪接事件及8次外显子互斥事件。RT-PCR结果证实了3种AS事件类型的真实性。共预测到38个泛素基因含有1个及以上的APA位点,并且在APA位点的上游鉴定到多个motif,一致性序列为：KCWYTDYTMWSYGMW SCARAWCCAGAATGAYCCWYWGGHWVMWGWDRTRGC。研究结果丰富了中华蜜蜂泛素基因及其全长转录本信息,为持续深入开展相关功能研究提供参考和依据。     

基于全长转录组测序的牛脂肪组织可变剪接事件分析

为探究可变剪接影响下更为复杂的脂肪沉积调控网络,明晰可变剪接事件的发生对脂肪沉积调控机制的影响,本研究基于PacBio测序平台的第三代测序技术,对夷陵黄牛腹腔脂肪、皮下脂肪、肌间脂肪进行全长转录组测序分析。共发现15 445个基因检测到50 520个可变剪接,占牛全部基因的33.4%。对这些基因进行富集分析发现,83个GO条目显著富集,这些显著富集的通路可能参与脂肪沉积网络的调控,其中16个与脂质合成及代谢相关,27条KEGG通路显著富集,其中15条与脂质合成及代谢相关。使用KOG、KEGG、NR、Swiss Prot、GO数据库对测序序列进行注释,80 756条序列共注释到69 259个基因,94 458条CDS序列及对应的pep序列。其中GO分析中15 039条序列富集到507个生物学过程条目,7 907条序列富集到214个细胞组分条目,30 672条序列富集到559个分子功能条目。KEGG数据分析富集到34条通路共49 710条序列,其中7 002条序列参与信号转导途径。以上结果表明,牛脂肪组织存在大量可变剪接事件,并且这些发生可变剪接事件的基因对脂肪沉积的调控发挥重要作用,这可...     

基于纳米孔全长转录组数据完善东方蜜蜂微孢子虫的基因组注释

【目的】利用已获得的纳米孔全长转录组数据对现有的东方蜜蜂微孢子虫(Nosema ceranae)参考基因组的基因序列和功能注释进行完善。【方法】采用TransDecoder软件预测东方蜜蜂微孢子虫基因的开放阅读框(open reading frame,ORF)及相应的氨基酸。利用gffcompare软件将全长转录本与参考基因组注释的转录本进行比较,对基因组注释基因的非编码区向上游或下游延伸,修正基因的边界。利用MISA软件鉴定长度在500 bp以上的全长转录本的简单重复序列(simple sequence repeat,SSR)位点,包括单核苷酸重复、双核苷酸重复、三核苷酸重复、四核苷酸重复、五核苷酸重复、六核苷酸重复、混合SSR等类型。通过Blast工具将鉴定到的新基因和新转录本比对Nr、KOG、eggNOG、GO和KEGG数据库,从而获得功能注释。【结果】共预测出2 353个完整ORF,其中长度分布在0—100个氨基酸的ORF最多,占总ORF数的72.12%。共对东方蜜蜂微孢子虫的2 340个基因进行了结构优化,其中5′端延长的基因有1 182个,3′端延长的基因有1 158个。共...     

基于RNA-seq数据分析玉米抗虫响应基因的可变剪接事件

采用Illumina Hiseq~(TM)分别对黏虫取食的玉米叶片及对照组材料进行转录组测序,鉴定和分析响应黏虫取食基因的可变剪接事件。结果表明,对照组中鉴定出15 701个基因对应79 363个可变剪接事件,黏虫取食组中鉴定到11 791个基因的39 385个可变剪接事件。2组玉米基因组在不同的可变剪接类型中,均是以第1个外显子可变剪切、最后1个外显子可变剪切、单内含子滞留和可变5′或3′端剪切4种类型为主。对2组发生可变剪接事件的基因进行比对发现,黏虫处理组新增加了1 121个可变剪接基因,差异表达分析鉴定出177个表达显著差异的可变剪接基因。GO功能富集分析表明,发生可变剪接的差异基因主要富集于氧化还原酶活性、转移酶活性、DNA结合、ATP结合、代谢过程、以DNA为模板的转录及调控相关的功能,表明这些可变剪接差异表达基因参与了玉米的抗虫响应过程。     

简州大耳羊早期胎儿到新出生阶段背最长肌中可变剪接的鉴定与分析

可变剪接广泛存在于动植物的各个组织中,是形成动植物体内蛋白质多样性和复杂性的重要机制之一。以简州大耳羊为研究对象,采集母羊妊娠期第45天、60天、105天胎儿以及出生后第3天羔羊的背最长肌组织,进行转录组测序。测序数据经比对、组装后,共鉴定出来自于4 562个(16.77%)已注释基因的137 308个可变剪接事件。其中,所占比例最高的类型为外显子跳跃(30.64%),而3′端可变剪接位点、内含子保留、5′端可变剪接位点的比例分别为19.47%、18.07%和10.92%。对RBMs、HNRNPs和SRSFs 3类可变剪接调控因子中58个基因表达量分析表明,19个基因的表达量呈现显著差异(P0.05);其中15个基因在妊娠期第45天和第60天的表达量显著高于其他2个阶段。GO功能注释和富集分析结果表明,4个阶段均发生可变剪接的基因显著富集在细胞表达调控和物质合成等过程(P0.05)。值得注意的是,在妊娠期第60天,可变剪接事件、可变剪接基因的数量和可变剪接因子表达量均显著高于其他阶段,并且特异性基因显著富集在有丝分裂和大分子合成相关过程。上述结果表明,妊娠期第60天胎儿的成肌过程非常活跃,暗示可变剪接在山羊背最长肌发育过程中具有重要作用,为深入探索可变剪接在山羊肌肉发育中的作用提供理论依据。     

意大利蜜蜂Hippo信号通路相关基因及其全长转录本的鉴定与分析

采用Blast工具将已鉴定到的意大利蜜蜂(Apis mellifera ligustica，简称意蜂)全长转录本比对Nr数据库，共鉴定到意蜂Hippo信号通路中相关的49个基因及其550条全长转录本。运用Gffcompare软件将鉴定到的Hippo信号通路相关全长转录本与西方蜜蜂参考基因组(Amel_HAv3.1)上注释的转录本进行比较，分别延长西方蜜蜂参考基因组注释到Hippo信号通路的9个基因的5'' UTR和7个基因的3'' UTR。通过Astalavista软件鉴定到Hippo信号通路相关的4个基因的7次可变剪切(AS)事件，包括2次外显子跳跃、2次内含子保留和3次可变5''端剪接。通过PCR验证了1个基因的AS事件真实性。利用TAPIS pipeline鉴定到意蜂Hippo信号通路相关的24个基因含有1个及以上可变多聚腺苷酸化(APA)位点，其中含有1个APA位点的基因数量最多。通过TBtool软件鉴定APA位点上游motif的一致性序列HVNCCNBNDYVDVVHVNDBVDYCNBBDNNNHNNVNNVMNN。     

马氏珠母贝(Pinctada fucata)血细胞全长转录组测序与分析

[目的]开发马氏珠母贝(Pinctada fucata)基因组资源,挖掘功能基因。[方法]以马氏珠母贝血细胞为试材,利用单分子实时技术(single-molecule real time, SMRT)进行全长转录组测序,并对所获得unigenes进行功能注释和基因结构分析。[结果]共获得277 064条全长非嵌合序列和82 381个基因。经过比对nr、SwissProt、KEGG和KOG数据库进行注释和功能分类后,共得到59 621个注释基因。同时,在8 493条基因中共发现11 219个SSR位点,其中以二核苷酸重复基元类型最高(50.9%)。生物信息学分析鉴定得到20 013个lncRNA、2 004个转录因子、10 522个可变剪切分析位点。[结论]使用SMRT技术能够深入挖掘马氏珠母贝全长转录组数据,为进一步探讨马氏珠母贝功能基因的挖掘、免疫响应及遗传机制的研究提供可靠的基因组资源。     

不同繁殖力绵羊睾丸组织新转录本注释及功能解析

【目的】完善和优化绵羊（Ovis aries）基因组注释信息,解析绵羊睾丸组织在不同繁殖力水平的遗传机理。【方法】以高繁殖力绵羊品种湖羊和低繁殖力绵羊品种藏绵羊（各3只）为研究对象,采集其睾丸组织,提取RNA,构建cDNA 文库,进行RNA-Seq测序,对测序数据进行组装。利用Cuffcompare软件对重新组装的测序数据与绵羊参考基因组进行比对分析,优化已注释基因的结构,并挖掘新转录本,对新转录本进行GO和KEGG通路分析。使用DESeq软件挖掘差异表达新转录本,对其进行GO和KEGG通路分析。随机选取表达显著上调和下调的新转录本各5个,采用实时荧光定量PCR（RT-qPCR）对其表达水平进行验证。【结果】对6只供试羊构建了6个独立的cDNA文库,共获得266 009 566个原始读段,质量控制后得到133 004 783个过滤后读段;6个cDNA文库过滤后读段的GC含量在49.61%～51.14%,符合碱基组成规律,Q30比例≥92.96%。对4 273个已注释基因的结构进行了优化,其中5′非翻译区域（UTR）延伸基因2 027个,3′ UTR延伸基因1 907个,5′和3′ UTR同时延伸基因339个。挖掘到12 178个新转录本,GO功能分析显示,2 416个新转录本注释到生物学过程、细胞组成和分子功能;KEGG富集发现新转录本显著富集在单纯疱疹病毒Ⅰ型感染、过氧物酶体、河马信号通路。与湖羊相比,藏绵羊睾丸组织中有155个新转录本表达差异显著,其中103个表达上调,52个表达下调;GO功能分析显示,差异新转录本富集在细胞过程、膜部分、绑定等条目;KEGG富集发现差异新转录本显著富集在钙信号通路和核因子-κB炎症信号通路等。10个差异表达新转录本的RT-qPCR验证试验结果与RNA-Seq结果趋势一致。【结论】从湖羊和藏绵羊睾丸组织中共挖掘出12 178个新转录本,其中155个差异表达。这些新转录本主要通过参与细胞功能、生物调节和代谢途径,以及通过核因子-κB炎症信号和钙信号通路广泛参与调控生殖细胞的增殖与分化、睾酮的分泌、免疫应答、Ca²⁺跨膜转运等途径,进而调控绵羊睾丸的发育和精子的发生。     

中华蜜蜂酚氧化酶和丝氨酸蛋白酶相关基因 及其全长转录本发掘与分析

为深入开展中华蜜蜂相关研究提供参考信息和基础,利用已获得的纳米孔(Nanopore)测序数据对中华蜜蜂Apis cerana cerana的酚氧化酶(phenoloxidase, PO)和丝氨酸蛋白酶(serine protease, SP)相关基因和全长转录本进行发掘和分析。通过Blast工具将中华蜜蜂所有全长转录本比对KEGG和Nr数据库以筛选出PO和SP相关基因和全长转录本。采用gffcompare软件将PO和SP相关全长转录本与东方蜜蜂参考基因组(ACSNU-2. 0)已注释基因进行比较以优化结构。使用Astalavista软件鉴定PO和SP相关基因的可变剪接(alternative splicing, AS)事件,再利用IGV浏览器进行结构可视化。通过PCR验证AS事件的真实性。利用TAPIS pipeline预测和分析可变多聚腺苷酸化(alternative polyadenylation, APA)位点,并通过TBtool软件鉴定APA位点上游的基序(motif)。鉴定到中华蜜蜂PO和SP相关的42个基因与146条转录本。优化了16个参考基因组已注释PO和SP相关基因的结构,其中5′端延长和3′端延长的基因均有7个,5′端和3′端同时延长的基因有2个。鉴定到PO和SP相关的10个基因的389次AS事件,其中最丰富的AS类型是5′端可变剪接。PCR结果证实了2次AS事件的真实性。共鉴定到34个PO和SP相关基因含有1个及以上的APA位点,其中含有3个APA位点的基因数量最多;在APA位点上游鉴定到多个motif,一致性序列为：GGHKSYWSHHTRATWTCNBHDMRRYWYRTNYTVACNGCKGCDCAYTGYR。鉴定到中华蜜蜂PO和SP相关的42个基因和146条全长转录本,优化了东方蜜蜂参考基因组已注释的PO和SP相关基因结构,并发掘出PO和SP相关基因的389次AS事件和237个APA位点。     

高脂饲喂对肉鸡腹脂和肝脏组织RNA 可变剪接的影响

【目的】腹脂过度沉积会降低肉鸡抗病力和屠宰率。RNA可变剪接是影响基因功能的重要调控方式,对肉鸡腹脂沉积过程的RNA可变剪接事件进行系统分析,有助于进一步了解肉鸡腹脂沉积的分子调控机理。【方法】利用高脂饲喂肉鸡和普通饲喂肉鸡肝脏组织和腹脂组织的转录组测序结果,鉴定高脂饲喂过程中肉鸡腹脂和肝脏组织的差异表达和差异剪接基因。采用GO功能注释、KEGG通路富集分析和Metascape富集分析对显著差异剪接基因进行分析,并对显著差异剪接事件进行可视化处理,以及分析调控这些差异剪接事件的剪接因子。【结果】对高脂饲喂肉鸡与普通饲喂肉鸡的腹脂和肝脏组织进行转录组测序,在腹脂组织中发现233个显著差异表达基因和349个显著差异剪接基因,对显著差异剪接基因进行富集分析,发现显著差异剪接基因主要富集于细胞分裂、细胞生长等细胞增殖相关生物学进程,还富集于甘油磷脂代谢和胰岛素信号通路等脂肪生成相关通路。在肝脏组织中发现276个显著差异表达基因和224个显著差异剪接基因,对显著差异剪接基因进行富集分析,发现显著差异剪接基因主要富集于蛋白质乙酰化、蛋白质内部氨基酸乙酰化和m RNA加工等细胞代谢相关生物进程,还富...     

两种穗伸出差异显著的大麦最上节间全长转录组分析

【目的】对两份穗伸出差异显著的大麦(Hordeum vulgare L.)材料的最上节间(UI)进行全长转录组分析,为揭示大麦UI伸长调控机制奠定基础。【方法】以包穗大麦材料CDB0012和穗伸出正常大麦品种XL19为材料,利用Nanopore平台开展了UI全长转录组测序和生物信息分析。【结果】共获得46 246条非冗余全长转录本,鉴定了2 395个潜在新基因位点和44 805条新转录本,开展了功能注释、可变剪切、差异表达基因鉴定等分析。【结论】构建了一套高质量大麦UI全长转录组数据集。     

小麦Ta NAC基因基于可变剪切和micro RNA的转录后调控分析

【目的】以条锈菌和白粉菌胁迫的普通小麦(Triticum aestivum L.)为研究对象,分析由可变剪切(alternative splicing,AS)形成的TaNAC结构变异转录本,同时分析TaNAC基因的microRNA调控位点,为进一步解析TaNAC基因通过转录后调控参与小麦响应真菌胁迫奠定基础。【方法】普通小麦兼抗种质N9134在被白粉菌和条锈菌分别侵染后,各8个时间点取样并混合,然后从混合样本中克隆得到大量TaNAC转录本。参考中国春小麦基因组注释信息(IWGSC RefSeq v1.1)进行比对,选择由可变剪切形成的TaNAC序列结构变异转录本,分析它们的序列结构特征。利用生物信息学软件和在线工具,对这些TaNAC结构变异转录本编码产物的功能结构域、高级结构、理化性质、亚细胞定位等特征和变异情况进行比对分析。同时,利用洋葱表皮细胞瞬时表达系统验证其中1对TaNAC结构变异转录本的亚细胞定位预测结果,并选取5组TaNAC基因的可变剪切序列结构变异转录本进行酵母转录自激活试验,研究序列结构变异对TaNAC基因转录调控活性的影响。此外,利用miRBase数据库收录的小麦中已...     

含热带种质玉米自交系ZOL-1胚乳发育的转录组研究

【目的】研究含热带种质玉米自交系胚乳发育分子机制及相关调控网络。【方法】通过高通量RNA-seq测序技术,对热带血缘玉米自交系ZOL-1在授粉后5、10、15、20、25、30 d 6个时期的胚乳转录组进行比较,得到的差异表达基因进行KEGG与GO功能注释,并对FPKM平均值大于1的基因进行时间序列分析。【结果】玉米自交系ZOL-1胚乳发育在授粉后25～30 d时籽粒长、宽、厚增长速率最快,增长率分别达到15.9%、16.1%和10.2%。经过转录组测序和同参考基因组比对,过滤后reads占比至少为83.99%,测序质量和组装效果均较好。按照连续比较系统方法,在10 d＿VS＿5 d、15 d＿VS＿10 d、20 d＿VS＿15 d、25 d＿VS＿20 d和30 d＿VS＿25 d分别检测到7129、5252、4254、775和312个差异表达基因。GO功能注释表明,差异表达基因主要分布于细胞分化过程、代谢过程、催化活性和细胞器等条目中。KEGG功能富集分析显示,差异表达基因主要富集在能量代谢、碳水化合物代谢、碳代谢、氨基酸代谢、信号转导等与胚乳发育相关的通路。通过对FPKM平均值...     

控制桃粘/离核PG基因的BAC克隆筛选与序列分析

【目的】控制桃果实粘/离核性状的多聚半乳糖醛酸酶(polygalacturonase,PG)基因存在串联重复和大片段的缺失。本研究对粘核桃所处的F-M基因座序列特征进行分析,为开发相关分子标记提供依据。【方法】本研究利用构建的粘核桃单株‘87-7-1’基因组的细菌人工染色体(bacterial artificial chromosome,BAC)文库,通过PCR筛选出含F-M基因座的阳性克隆,利用单分子纳米孔技术进行全长测序,并对BAC克隆中的插入片段进行基因注释、序列比对和生物信息学分析。【结果】利用已有桃品种的基因组重测序数据设计PCR引物,以BAC文库为模板,扩增F-M基因座上/下游稳定共有的序列,获得了扩增产物上/下游条带均为阳性的目标单克隆46-B-10。全长测序结果表明,BAC克隆的插入片段全长为111 612 bp,GC含量为37.03%。利用基因同源共线性方法对Prunus_persica_v2.0桃参考基因组和BAC克隆之间的同源信息进行比对分析,确定了同源区域。而与已知的桃参考基因组(测序品种为‘Lovell’,离核)序列进行...