碧桃花瓣转录组微卫星特征分析

为开发碧桃转录组微卫星信息,利用454高通量测序技术,对其花瓣转录组序列进行SSR位点发掘,结果发现含SSR的序列4 705条,共得到5 668个SSR,平均每3.49 kb出现1个SSR。微卫星序列主要以三碱基重复为主,约占总数的42.66%。笔者共发现516种碱基重复基元,所占比例最高的为(AG/CT)n(18.34%),其次是(AAG/CTT)n(12.42%)。微卫星多为重复长度小于20 bp的短序列,长度大于20 bp的微卫星仅占总数的12.13%。研究还发现碧桃花瓣微卫星的频率和长度呈显著负相关(P0.05),相关系数为-0.246。     

湿地松转录组SSR分析及EST-SSR标记开发

【目的】湿地松是南方地区优先推广的优质产脂树种。但湿地松分子遗传基础薄弱,基因组序列信息匮乏,影响了湿地松基因组学的深入研究。目前,湿地松分子研究所用的SSR标记主要来自其他近缘种或利用公共数据库中有限的基因序列资源开发的SSR标记,其多态性和通用性较差。为解决这一问题, 笔者根据湿地松转录组测序数据开发EST-SSR位点,并揭示其在转录组序列中的分布类型及特征,为湿地松分子标记辅助育种奠定基础。【方法】利用MISA软件对转录组序列进行SSRs查找和分布特征分析。查找标准参数设置为: 单核苷酸重复>10次,二核苷酸重复>6次,三、四、五、六核苷酸重复>5次。根据SSR位点两端的保守区域,利用Primer3.0设计并随机挑选120对SSR引物,通过琼脂糖电泳和毛细管电泳对来自美国和吉安的113份家系个体进行遗传多样性分析,确定引物多态性。【结果】79 574条unigenes序列中搜索到3 818个SSR位点,出现频率为4.80%,平均18.27 kb出现1个SSR位点,3 373个unigenes含有SSR位点,SSR发生频率(含SSR位点的序列数/搜索序列总数)为4.24%,其中2 980条序列含1个SSR位点,含1个以上SSR位点的序列有393条。在检测到的3 818个SSR标记中,单核苷酸分布最多,其次是二核苷酸和三核苷酸,SSR数量分别占总数的63.54%、19.15%和16.27%,而四、五、六核苷酸重复类型所占比例较小,分别为0.52%、0.13%和0.31%。SSR重复单元的重复次数分布在5~22次之间,除单核苷酸重复外的1 391个SSR中,重复5次的数量最多,为498个(35.80%);重复6次和7次的次之,分别为417个(29.98%)和198个(14.23%);重复10次以上的仅有38个(2.73%)。在检测到的731个二核苷酸重复SSR中,最常出现的重复单元为AT/AT,数量为491个(12.86%),AG/CT和AC/GT两种类型的重复单元出现的次数次之,分别为156(4.09%)和81个(2.12%)。在检测到的621个三核苷酸重复中,AAT/ATT是出现频率最高的单元,共139个(3.64%),其次是AAG/CTT,共122个(3.20%)。3 818个SSR中有24.59%的位置未知,其余的SSR则分布在非编码区域(untranslated region,UTR)或者编码区(coding sequence,CDS)上,分布数量表现为3'UTR>5'UTR>CDS。参试的120对SSR引物,有92对扩增成功(76.78%),其中24对呈现多态(20%)。24对引物(13个二核苷酸重复、7个三核苷酸重复和4个四核苷酸重复)共检测出81个等位基因,每个位点的等位基因数为2~9,平均为3.38个。多态性信息含量(polymorphism information content,PIC)为0.103~0.726,平均为0.349。【结论】通过对湿地松转录组数据的挖掘,共获得 3 818个SSR位点,主要重复单元为AT/AT和AAT/ATT,可扩增出多态性位点的引物重复单元以二、三核苷酸重复为主。基于湿地松转录组序列的SSR标记开发是可行的。     

基于转录组测序的波纹巴非蛤微卫星标记研究

应用MISA软件对波纹巴非蛤转录组测序数据进行分析,共获得7378个微卫星标记,从中选取95个位点设计引物,并在一个野生波纹巴非蛤群体中进行分析验证,结果显示:95个位点中能够扩增出清晰条带的位点有30个,其中15个位点表现为单态,15个位点呈现多态.使用15个多态位点进行波纹巴非蛤群体的遗传多样性分析,结果表明:这15个位点的等位基因数为2~9,平均等位基因数为4.733,平均有效等位基因数为2.438,平均观测杂合度和平均期望杂合度分别为0.270、0.516,平均多态信息含量为0.457.     

雌核发育棕点石斑鱼EST-SSR的开发验证

通过人工雌核发育棕点石斑鱼(Epinephelus fuscogutatus)的转录组序列,共获得了24 359个微卫星位点,棕点石斑鱼转录组的微卫星重复基序具有不同的分布特征,其中以单碱基(34.05%)、二碱基(37.58%)和三碱基(22.75%)为主要的EST-SSR基序重复类型.根据微卫星引物设计原则,随机筛选了93个EST-SSR位点来进行引物合成和多态性检测,最终开发了48个具有多态性的微卫星标记,并在不同棕点石斑鱼群体中进行了初步的验证.结果显示:每个位点等位基因数(Na)为2～22,平均等位基因数为9;观测杂合度(H_0)为0.111～1.000,期望杂合度(He)为0.636～0.940;多态信息含量(PIC)为0.538～0.922个,表明48个EST-SSR位点均表现出高多态性(PIC 0.5).结果表明:基于雌核发育棕点石斑鱼转录组数据,开发微卫星标记是可行的.本研究所开发的具有多态性的微卫星标记可应用于棕点石斑鱼及其近缘种的遗传多样性分析和分子标记辅助育种研究.     

高通量转录组测序的数据分析与基因发掘

高通量转录组测序(RNA-seq)是在转录组水平上进行深度测序的一项技术,为真核生物转录组学的研究开创了新平台,但同时测序所得到的海量数据的生物信息学分析成为科研工作者的一大挑战。对转录组测序技术进行了阐述,重点介绍了转录组测序后的数据分析,以及在真核生物尤其是非模式物种中的基因发掘方法。     

云南松转录组SSR的分布及其序列特征

分析云南松转录组中的SSR位点信息的分布及其序列特征,为SSR引物的开发做前期工作.从80000条转录组Unigene序列中,搜索到2455个SSR,出现频率为3.07%,分布的平均距离29 kb.SSR以三、二核苷酸占优势,分别占总数42.73%和27.25%,而以四核苷酸的较少,仅占总数的0.98%.重复片段长度12~25bp,其中近90%的集中在12~20bp,平均16.09bp.AT/AT与AGC/CTG和AAG/CTT为二、三核苷酸优势重复基元,分别占基元总数的16.66%、9.65%和8.39%.从2899对引物中随机挑选32对引物用于检测,24对成功扩增.云南松转录组SSR位点出现的频率较高,分布距离较近,基元类型丰富,重复次数较高,具有较高的多态性潜力.云南松转录组产生的Unigene是开发SSR的有效来源,获得的SSR引物可为云南松种质资源的遗传变异研究提供更多的标记.     

基于转录组序列的小麦ETS-SSR标记筛选与染色体定位

通过小麦转录组序列分析,获得121 210个非重复序列(Unigenes),在10 672(8.8%)个Unigenes中搜索出1-6个重复基元的11 650条EST-SSR信息位点,筛选并设计出308条SSR引物,选取前30对引物进行合成.有效性扩增检测结果表明,20(66.7%)对引物有清晰条带.利用缺体-四体系共定位了14对引物,分别位于13条染色体的21个位点上.研究表明,利用小麦转录组中EST-SSR信息开发新的SSR标记是可行的,开发的新ESTSSR标记可有效用于小麦基因的定位和遗传多样性的分析等.     

拟南芥花粉单细胞转录组分析

为了深入探究花粉基因表达的规律,应用单细胞技术对拟南芥花粉细胞进行了单细胞转录组测序.花粉单细胞RNA-seq结果表明,与叶片转录组相比,成熟花粉中大部分基因表达较低,可能与花粉转录沉默有关.花粉中富集果胶代谢、细胞骨架和钙信号转导等通路,与细胞壁形成和花粉管生长相关.花粉还富集组蛋白甲基化因子,说明花粉中存在表观遗传学调控.应用RT-PCR鉴定出6个花粉特异性表达且功能未知的基因,其中大多数为十字花科特有.总之,本研究精确定量了花粉细胞的基因表达谱,揭示出拟南芥花粉转录组的复杂性和独特性.     

鸟王茶转录组中SSR位点信息分析

为开发鸟王茶分子标记技术,对使用高通量测序所获得的鸟王茶转录组原始数据,经过软件Trinity拼接,共获得161 561条Unigenes,使用软件MISA进行SSR位点搜索,共得到51 453个SSR,分布于36 691条Unigenes中,总长度为554 050bp,发生频率为31.85%,平均分布频率是1/1.79kb,平均长度为10.77bp。在鸟王茶转录组SSR中,二核苷酸重复和三核苷酸重复是主要的重复类型,分别占总SSRs的45.05%和13.86%。鸟王茶转录组SSR种类共有223种重复基元类型,其中优势重复基元是A/T、AG/CT、AC/GT它们所占比例分别是:36%、32%、13%,这3种重复基元,占总SSRs的比例是81%。     

基于色木枫转录组数据的核低拷贝基因引物开发

基于色木枫的转录组数据初步设计166对引物,并根据琼脂糖凝胶电泳与Sanger测序结果筛选出其中的26对核低拷贝基因引物．在东亚63个种群的182个个体中大批量扩增所开发引物,单倍型多样性范围为0.45~0.97,核苷酸多样性为(1.90~18.03)×10?3;9~15对引物可在金钱槭、茶条槭、青榨槭和鸡爪槭中扩增;STRUCTURE聚类结果与之前核微卫星研究结果类似,且展示了更多的谱系地理结构,结果表明引物多态性高、拓展性强、具有可靠性．本研究的转录组数据为槭属的研究提供重要遗传信息,所用方法被证明可以有效开发大量核低拷贝基因引物,已开发出的26对引物将为色木枫的进化历史和槭属系统发育研究提供丰富有效的分子标记．      

灰树花菌（Griflola frondosa）子实体原基转录组测序及分析

采用Illumina测序平台对灰树花菌（Griflola frondosa）原基进行转录组测序,获得了大量转录组信息.结显示果,测序得到38 189个非冗余Unigene,在Nr数据库中有31 454个Unigene被注释;与KEGG数据库比对,共有1 832个Unigene被注释,分属20条生物通路;在COG数据库中注释可划分为25类;与GO数据库进行注释可分为分子功能、细胞组分和生物过程3个主类和60个二级亚类.     

三叶木通雌雄花芽转录组测序及分析

为寻找三叶木通雌雄花性别分化的相关基因,从基因的表达水平上揭示三叶木通雌雄分化的分子机制.以三叶木通雌雄分化期的花芽为材料,采用高通量测序技术,对自然状态下三叶木通的雌性花芽与雄性花芽进行转录组测序,通过生物信息学对雌雄花芽的转录本进行分析,筛选差异表达基因.将雌花与雄花Unigene进行表达量分析,筛选表达量高且可信...     

转录组测序技术在绵山羊育种中的应用研究进展

转录组测序技术已广泛应用于生命科学的各个领域.羊作为被最早驯化的家养动物之一,在人类日常生活中有着不可取代的作用.文章主要从绵山羊肉品质、毛绒生长、泌乳等方面应用转录组测序进行了阐述,为今后绵山羊的分子育种工作提供参考.     

基于PacBio平台的七叶一枝花全长转录组测序

以根状茎、茎、叶、花、果荚及种子6个部位为样品,利用PacBio技术对七叶一枝花进行了全长转录组测序,并对测序数据进行了生物信息学分析.研究共获得52537个平均长度为2607 bp的高质量isoforms,其中40725条isoforms得到注释.GO注释中,4596个isoforms分为生物学过程、细胞组成和分子功...     

双棘黄姑鱼染色体组型分析

以头肾细胞为材料，采用PHA体内注射。空气干燥法制片，分析了双棘黄姑(Nibea diacanthus)的染色体组型。结果表明：其染色体组型为：2n=48=24t，NF=48。在已进行过染色体组型分析的石首鱼类中，除大黄鱼外，染色体组型均为2n=48t。     

冷驯化下的牛皮杜鹃转录组变化分析

冷胁迫是限制植物生长和分布,降低作物产量的重要环境因子.冷驯化是温带植物对抗冷胁迫的有效手段.本研究以具有耐寒能力的高山植物牛皮杜鹃(Rhododendron chrysanthum Pall.)作为研究对象,探究冷驯化前后牛皮杜鹃的基因表达和代谢通路的变化.基于转录组数据的测序结果:常温组和驯化组的总reads分别为...     

绵羊骨骼肌转录组高通量测序从头组装和特征分析

本实验旨在对绵羊肌肉转录组进行组装分析,以丰富绵羊的基因组并为进一步对绵羊的相关基因的分子遗传学研究奠定基础。选取Illumina Hiseq 2000测序平台对杜泊羊和小尾寒羊的骨骼肌转录组文库进行了高通量测序,并对测序数据进行从头组装分析。结果表明:两个测序文库共得到103058824个长为2×90bp的高质量的测序序列,经从头组装共产生了145524个unigenes。两个文库间有5718个unigenes差异表达,经GO注释分析后发现它们主要分布在细胞、细胞过程和结合条目中。将两个文库的unigenes进一步组装得到70348个平均长为863bp的all-unigenes,其中35201个被注释到Nr数据库,12219个被注释到COG数据库,并与KEGG数据库中的258个功能通路中的蛋白质和氨基酸新陈代谢通路密切相关。     

基于转录组数据对糙苏三萜皂苷合成途径的研究

糙苏为我国民间常用的传统草药,具有消肿、续筋接骨、祛风化痰、通络止痉、安胎等作用,三萜皂苷类化合物为其主要的生物活性成分.对糙苏组织转录组测序,挖掘其三萜皂苷生物合成途径相关基因,探究糙苏三萜皂苷化合物生物合成的分子基础.利用Trinity软件对测序结果获得的Unigenes数据进行过滤组装,并结合KEGG、GO数据库对Unigenes进行注释,预测糙苏三萜皂苷代谢的生物合成途径关键基因.结果获得了参与萜类化合物骨架生物合成的Unigenes数量为73条,39条Unigenes编码19个与萜类生物合成的关键酶,且有23条基因参与萜类物质骨架代谢途径(途径编号为ko00900),19条基因参与了倍半萜类物质代谢途径(途径编号为ko00909),涉及萜类物质骨架合成代谢途径相关酶共有15种,以及三萜代谢途径相关酶有10种.本研究获取了糙苏的三萜皂苷合成相关候选基因核苷酸和氨基酸序列,为糙苏三萜皂苷物质体外表达和生物合成奠定理论基础.     

一款基于转录组差异基因表达分析的软件包——findDEG

【目的】随着二代测序技术的不断发展,转录组测序技术在许多物种里已被广泛地应用于基因差异表达分析和基因注释研究。现有的多种基因差异表达分析软件,分析步骤多而且复杂,不同分析方法其结果差别也较大,这给研究者分析实际数据带来了不少困难。为了简化基因差异表达分析的过程,利用现有的软件开发一个集成的软件包。【方法】针对Trinity、TopHat+Cufflinks和HISAT2+StringTie 3种比较成熟的基因差异表达分析流程,考虑研究对象有无参考基因组序列、样本数据是否有重复、单端还是双端测序、不同基因表达量的计算方法以及不同的基因差异表达显著性检验方法等因素,将多种转录组测序数据分析软件整合起来形成一个集成的软件包。【结果】 使用Perl语言开发了一个名为findDEG软件包用于转录组测序数据的基因差异表达分析。软件包共分为3个模块,即Trinity、TopHat+Cufflinks和HISAT2+StringTie模块。Trinity模块提供3种计算转录本表达量方法和4种差异表达基因显著性检验方法,TopHat+Cufflinks模块可供用户选择新版或旧版的Cufflinks分析方案,HISAT2+StringTie模块则只有一种分析方案。该软件包可以自由下载使用,其网址为http://www.bioseqdata.com/findDEG/findDEG.htm。采用新版和旧版的Cufflinks分析方案以及一种Trinity组合方法,分别对小叶杨在正常和干旱胁迫条件下的转录组数据进行了分析。结果两种Cufflinks方法分别识别出了53和33个差异表达基因,其中25个是相同的; Trinity方法识别了高达1 641个差异表达基因,其中与Cufflinks两种方法相同的分别有14和3个。【结论】 新开发的软件包findDEG有十多种基因差异表达分析方案可供选择,采用一键的方式进行数据计算分析,避免了中间环节参数输入和结果利用等操作步骤,使用方便。     

祁连山黄参叶片转录组测序及生物信息学分析

在《本草纲目》中黄参被誉为“小人参”,以黄参叶片为试材,采用高通量测序平台BGISEQ-500进行转录组测序,利用转录组分析软件进行组装、注释。结果表明：1)利用组装软件,获得99 981个Unigene,总长度是113 850 816 bp,平均长度是1 138 bp, N50的长度是1 874 bp, GC含量是39.93%。2)将Unigene比对到7大功能数据库进行注释,分别有49 390(NT:49.40%)、48 281(SwissProt:48.29%)、61 116(KOG:61.13%)以及55 859(Pfam:55.87%)个Unigene获得功能注释。3)比对到NR数据库共有66 451条,黄参与胡萝卜Daucus carota subsp.sativus有较高同源性,与其他物种的同源性较低。4)基因本体(gene ontology, GO)数据库注释显示,有78 040条Unigene得到注释,按功能分为生物过程、细胞组分、分子功能三大类,分别有15、11、14个亚类,其中执行生物过程的类区较多。5)51 479条Unigene富集在KEGG数据库的20条代谢...