红脚艾蒿的转录组解析

[目的]利用高通量测序技术解析红脚艾(Artemisia rubripes Nakai)的转录组信息特征.[方法]通过高通量测序平台Illumina HiSeq 2500对红脚艾进行转录组测序,通过Trinity软件de novo组装获得Unigene,并基于序列同源性对Unigene进行功能注释,得到红脚艾的转录组信息.[结果]测序数据经过质控后共获得24126043条高质量的reads,通过de novo组装获得173093个转录本,对组装的转录本去冗余后共获得85991个Unigene,平均长度为616.87 bp,N50为925 bp.共有47216个Unigene在NR、KEGG、COG、KOG、GO数据库获得功能注释,40802个Unigene在NR数据库注释,显示红脚艾与向日葵(Helianthus annuus)的单基因匹配率最高,16846个Unigene被KEGG数据库注释到130条代谢途径中,26171个Unigene被注释到25个KOG功能分类中,23203个Unigene被GO注释到生物过程、细胞组成和分子功能三大类51个功能分类,12810个Unigene被注释到25个COG功能分类中.[结论]利用高通量测序技术获得了红脚艾转录组信息特征,这些数据将为后期开展功能基因鉴定、解析化合物次生代谢途径及其调控机制奠定研究基础.     

冬瓜高通量转录组测序及分析

【目的】获得冬瓜转录组序列、遗传变异等信息,从中挖掘冬瓜基因数据及SSR分子标记,为冬瓜后续研究提供数据支撑。【方法】以冬瓜嫩叶为材料,利用Illumina HiSeq~(TM)2000技术对冬瓜进行转录组测序,构建数据库从中获得干净序列。经De novo拼接组装后,将获得的单基因簇(Unigene)数据在非冗余蛋白数据库(nonredundant protein database,Nr)、蛋白质序列数据库(Swiss Prot protein database,Swiss Prot)、基因本体论数据库(gene ontology,GO)、蛋白质真核同源数据库(eukaryotic orthologous groups,KOG)、东京基因与基金组百科全书(Kyoto encyclopedia of genes and genomes,KEGG)、蛋白质家族域数据库(protein families database,Pfam)6个公共数据库中进行比对,最终得到冬瓜单基因簇注释信息。利用MISA软件对转录组单基因簇进行搜索,获得单基因簇中的SSR位点。【结果】从冬瓜嫩叶中得到62 021 032条高品质序列,组装后获得40 611条单基因簇,平均长度955 bp。将所有单基因簇在Nr和Swiss Prot数据库中进行比对,结果分别比对到27 474及19 573条单基因簇;在GO数据库中,所注释到的10 659条单基因簇分别匹配到生物功能、分子功能和细胞组分3个本体的47个功能组中;与KOG数据库进行注释比对,根据其功能将注释到的单基因簇划分为25类;KEGG数据库比对注释到10 799条冬瓜的单基因簇,可分为5个大类、19个亚类、125条代谢途径;在Pfam数据库中比对到17 990条单基因簇,分属于369个类群。SSR位点搜索发现,有5 086条单基因簇包含SSR序列,获得5 474个SSR位点。【结论】利用高通量测序获得大量冬瓜转录组信息,有助于从分子水平对冬瓜进行深入研究。     

基于转录组测序技术的黄唇鱼SSR分子标记筛选

[目的]筛选出可用于黄唇鱼(Bahaba flavolabiata)遗传资源评价及亲缘鉴定的SSR分子标记,为开展其保护遗传学研究打下基础.[方法]利用Illumina HiSeq 4000测序平台对野生黄唇鱼样本进行转录组测序,采用TGICL对转录组数据进行聚类去冗余以获得Unigene,并通过生物信息学方法进行功能注释分类、代谢通路和SSR特征分析等.[结果]黄唇鱼混合组织的转录组测序共获得13.43 Gb数据,经组装并去冗余后获得65047条Unigenes.采用BLAST对组装获得的Unigenes进行七大功能数据库(NR、NT、GO、COG、KEGG、Swissprot和Interpro)注释,结果共有55583条Unigenes被注释(占85.45%).通过MISA对黄唇鱼的Unigene进行SSR搜索,共检测到29797个SSR位点分布于19664条Unigenes中.SSR位点分布类型及其特征分析结果显示,最主要的重复类型为二核苷酸重复基元(11855个,占39.79%),其次是单核苷酸重复基元(9695个,占32.54%)和三核苷酸重复基元(6663个,22.36%),而四核苷酸、五核苷酸和六核苷酸重复基元的数量相对较少.二核苷酸重复基元类型中重复最多的基序是AC/GT(8355条),在三核苷酸重复基元类型中重复最多的基序是AGG/CCT(1866条).从随机选取的186对SSR引物中,最终筛选出47对扩增稳定性好且具有多态性的SSR引物,以二核苷酸和四核苷酸重复基元的扩增结果居多.[结论]通过转录组测序能有效筛选出具有多样性的黄唇鱼SSR分子标记,可为黄唇鱼的遗传多样性分析、亲缘鉴定和遗传图谱构建等后续研究提供基础资料.     

杜鹃花叶片转录组测序数据组装及功能注释

采用2代 Illumina Hi-Seq 测序技术对2年生杜鹃花白凤4号(Rhododendron pulchurum cv. BaiFeng 4)扦插苗叶片的转录组进行测序,获得了213 723 424条 Clean reads数据,经过质控和De novo。组装共获得平均长度为930 nt的53 568个Unigenes。与Nr、Swiss-Prot、KOG、KEGG四大数据库进行BLAST信息比对(E-value≤10^-5),共获得28 877个注释基因。与Nr数据库序列同源性比较发现,杜鹃花与葡萄(Vitis vinifera)具有较高的同源性,而与其他物种的同源性较低。杜鹃花转录组的Unigenes在KOG数据库比对上30 508个注释,分为25类。在GO数据库比对上17 218个Unigenes,其中与抗逆有关的Unigenes有11 928个。在KEGG数据库的132条代谢通路中富集了6 475个杜鹃花Unigenes,其中注释到植物激素生物合成代谢途径和信号转导途径的有384个。同时,有1 062个Unigenes被注释为转录因子,有8 738条SSR标记被挖掘。     

大竹蛏高通量转录组测序数据组装和分析

为进一步开发大竹蛏Solen grandis的基因资源,采用2代Illumina Hi-seq测序技术对大竹蛏的鳃组织进行了转录组测序,构建了转录组数据库,获得338 483 476条Clean Reads数据;拼接组装后获得190 856条Unigene数据,平均长度为1147 bp;与NR、NT、KO、Swiss Prot、PFAM、GO、KOG等数据库进行Blast信息比对(E-value为10-5),共获得63 337个注释基因;与NR数据库比对发现,大竹蛏转录组基因序列与长牡蛎Crassostrea gigas具有较高的同源性,为53.3%;将大竹蛏转录组的Unigene的功能通过与KOG数据库进行注释比对划分为25类;GO数据库注释可分为三类,即细胞组分、生物过程和分子功能,共包括65个分支;KEGG分析发现,大竹蛏转录组数据中按照代谢通路可分为92类,利用Blast蛋白库比对和Estscan软件进行ORF预测,获得长度大于300 nt的ORF共50 681个;通过SSR分析,共获得73 089个SSR标记。本研究中获得的转录组信息可为今后进行大竹蛏分子标记的开发和关键基因的克隆及功能分析等研究提供基础数据。     

基于RNA-Seq的杜仲转录组微卫星特征分析

对杜仲(Eucommia ulmoides)国审良种‘华仲6号’和‘华仲10号’花后70和160d的种仁共4个样本进行转录组测序,对测序数据进行组装和功能注释分类,并对转录组获得的单基因簇(unigene)进行微卫星特征分析。利用新一代高通量测序技术Illumina HiSeq~(TM)2000对杜仲样品进行转录组测序,采用软件Trinity进行组装;利用BLAST软件将unigene序列分别与Nr、GO、COG和KEGG等数据库比对分析;利用MISA软件对转录组的96 469条unigenes进行SSR搜索。结果表明:转录组测序分析,共得到72 791 399个高质量的序列读取片段(Clean reads),包含了14 702 548 161个的碱基序列(bp)信息。对reads进行序列组装,共获得96 469个平均长度为690bp的unigene,序列信息量达到了66.56 Mb。同源性分析结果显示,有49 856个与其它物种同源的unigenes得到注释,占All-unigene的51.68%。将杜仲转录组中的unigene与GO数据库进行比对分析,根据其功能可将注释到的38 983条unigene分成3大类(细胞组分、分子功能和生物学过程)56个分支;根据COG功能可将注释的14 796条unigene基因划分成25个类别;KEGG数据库作为参照,可将注释到的11 260条unigene定位到117个代谢途径分支;SSR位点搜索结果显示,96 469条unigenes中共包含9 621个完整型SSR位点,占总SSR位点的84.14%。完整型SSR位点共包含55种重复基元,其中出现频率最高的重复基序类型为单核苷酸重复中的A/T(4 597个),其次是AG/CT(2 597个)、AT/AT(439个)。     

柴胡转录组SSR的分布及序列特征分析

【目的】分析柴胡转录组中SSR的分布及序列特征,为开发多态性良好的、功能基因相关的SSR标记提供理论依据。【方法】以不同温度处理的柴胡种子为材料,经高通量转录组测序后,使用Trinity对测序结果进行de novo组装,并利用MISA对组装得到的Unigenes进行SSR位点搜索,最后统计分析SSR的分布及序列特征。【结果】从转录组数据组装获得244194条Unigenes,其N50值为1036 bp,平均长度791 bp,总长度193138105 bp。从Unigenes序列中共检测到50303个SSR位点,去除20405个复合型SSR位点,以29898个单一型SSR位点为后续分析对象。转录组SSR的出现频率为12.24%,平均分布距离6.46 kb,主要重复基元类型为二核苷酸,共17105个,占SSR总数的57.21%,其次为三核苷酸和单核苷酸,分别占SSR总数的20.75%和20.46%,四核苷酸~六核苷酸重复基元数量均较少。转录组SSR中,共有97种重复基元,其中二核苷酸和三核苷酸重复基元分别以AT/AT和ATC/GAT为主,分别占SSR总数的33.33%和3.98%;重复次数5~10次的SSR位点数量最多,共27942个,占SSR总数的93.46%。转录组SSR序列长度存在明显差异（12~76 bp）,平均长度15.28 bp。【结论】柴胡转录组的SSR位点出现频率较高,类型较丰富,具有开发出高多态性SSR分子标记的潜力,将其用于柴胡的遗传多样性分析、种质资源评价及分子标记辅助育种等研究。     

阳春砂花丝、花柱转录组测序及生物信息学分析

为进一步探究阳春砂假合蕊柱形成机制,基于RNA-Seq技术对不同生长时期阳春砂花丝、花柱进行转录组测序及生物信息学分析,获得阳春砂转录组数据并筛选出激素合成及信号转导通路基因,为后续研究阳春砂基因功能、代谢途径、花器官发育及运动调控机理等方面提供参考。结果表明,测序数据拼接组装后获得94 584条Unigene,总长度为92 501 015 bp。将获得的Unigene分别在NR、GO、KOG和KEGG等七大数据库中进行比对,共有62 174条Unigene得到注释,占全部Unigene的65.73%;NR数据库中注释到58 669条Unigene,被注释的同源序列主要来自于小果野芭蕉;有45 892条Unigene注释到KOG数据库中,涉及25个功能分类;GO数据库注释Unigene 12 050条,按照功能分为3个大类及55个亚类;KEGG数据库中,有44 000条Unigene注释到137个代谢通路中;SSR特征分析中共检测到18 895个位点,三碱基重复的数目最多,达到6 313个,占比为33.41%。2012条Unigene被注释到9种不同激素的合成与信号转导途径中,赤霉素的...     

青花菜花蕾转录组测序分析及蜡粉合成相关基因挖掘

【目的】对青花菜花蕾进行转录组测序分析，并挖掘与蜡粉合成相关基因，为探明青花菜花球表面蜡粉形成的分子机制提供理论参考。【方法】分别提取野生型和蜡粉缺失型青花菜花球总RNA，采用Illumina HiSeqTM 2500平台进行转录组测序，获得高质量Clean reads，采用Trinity进行序列组装后获得青花菜Unigene库，将获得的Unigene序列与Nr、Nt、KEGG、Pfam、KOG/COG、Swiss-Prot和GO数据库比对，获得基因功能注释信息；使用DESeq2进行差异表达分析。【结果】共获得44.68 Gb Clean data，De novo组装得到41244条Unigenes，N50长度为1847 bp。从所获得的Unigenes中筛选出8685个差异表达基因（DEGs）（上调基因5747个，下调基因2938个），共有8038个基因被注释到不同数据库，其中，5220个基因注释到Pfam数据库； 2066个基因注释到COG数据库，3866个基因注释到KOG数据库； 2580个差异表达基因被注释到75个转录因子家族中，注释最多的是MYB家族（235个）；GO数据库中6095个差异表达基因注释到细胞组分、分子功能和生物学过程三大类的52个功能分类； KEGG数据库中，1671个差异表达基因富集到138条代谢通路，其中13个差异表达基因与脂肪酸合成有关，7个差异表达基因与蜡粉生物合成途径有关。【结论】转录因子MYB家族在调控青花菜蜡粉合成中发挥重要作用。蜡粉合成过程中相关酶基因的差异表达是调控青花菜蜡粉合成的关键，尤其是野生型和蜡粉缺失突变体中特异性表达的差异表达基因，可作为后续研究青花菜花球表面蜡粉形成分子机制的对象。     

藤茶高通量转录组分析及黄酮类化合物合成相关基因挖掘

[目的]分析藤茶高通量转录组序列,从中挖掘出黄酮类化合物合成相关基因,为进一步揭示藤茶黄酮类化合物生物合成调控机制提供理论参考.[方法]分别采集藤茶的幼叶和成熟叶,提取其总RNA构建cDNA文库,采用Il-lumina HiSeqTM 4000高通量测序平台对藤茶叶片进行转录组测序,经过滤处理后运用Trinity组装,将获得的Unigene与Nr、Nt、Pfam、Swiss-Prot、GO、KO和KOG 7个数据库进行比对注释,并预测Unigenes的编码区序列(CDS);基于KEGG信号通路富集分析,发掘藤茶黄酮类化合物合成相关基因.[结果]藤茶叶片转录组测序获得82126236条原始测序序列(Raw reads),过滤处理后得到80156972条高质量序列(Clean reads),进一步组装拼接得到92472条Unige-nes,平均长度为1208 bp,N50长度为1780 bp,其中,至少在1个数据库注释的Unigenes有84217条,占Unigenes总数的91.07%,有8944条Unigenes在7个数据库均被注释,占Unigenes总数的9.67%.在GO数据库成功注释的41116条Unige-nes可分为生物学过程、细胞组分和分子功能三大类,共56个小类;在KOG数据库注释的14553条Unigenes可分成25类,其中,一般功能预测注释成功的Unigenes最多(1946条);其次是翻译后修饰、蛋白质翻转、分子伴侣(1776条),参与次生代谢物质的生物合成、转运和降解的Unigenes较少,仅有319条;KEGG信号通路富集分析发现,共有15262条Unigenes注释到128条KEGG信号通路,以注释为代谢的Unigenes最多,为8694条,其中筛选获得有98个黄酮类化合物合成相关基因,分别编码苯丙烷代谢通路的3种关键酶和类黄酮代谢通路的14种关键酶.藤茶叶片转录组Unigenes与Swiss-Prot和Nr数据库比对,获得52582条CDS序列,ESTScan 3.0.3预测获得35535条CDS序列.[结论]藤茶在细胞过程、代谢过程、单有机体过程、细胞和细胞部分、结合和催化活性能力分布的基因较丰富,在一般功能、翻译、翻译后修饰、蛋白质翻转及分子伴侣的基因表达量较高,具有较强的碳水化合物代谢能力.多种关键酶基因参与藤茶黄酮类化合物的生物合成,推测其生物合成途径存在多条分支,调控机制也较复杂.     

鲫鱼血液、肝脏、卵巢组织转录组比较分析

为进一步丰富鲫鱼基因组数据，对肝脏、血液和卵巢组织进行高通量转录组测序。采用Illumina Hiseq高通量测序平台进行转录组测序，对所获得序列进行比对、基因注释、差异表达和SNP、SSR分析。结果提示，组装得到127 801条unigene,平均长度为735 bp。与相应数据库比对，最终获得117 414(91.87%)个有注释信息的unigene,有105条unigene获得GO注释，并分类到生物学过程、细胞组分和分子功能三大类42个功能中；有20 725条unigene被归纳到KOG的26个直系同源功能类型；获得KEGG注释的unigene共参与五大类32小类代谢通路。在鲫鱼肝脏、血液和卵巢组织中分别检测到307 740、343 516、348 795个SNP位点，共检测到48 808个SSR位点。提示获得有注释信息的unigene共117 414条，此外，获得一些SNP和SSR位点。结果可为进一步克隆和挖掘鲫鱼功能基因、多态性检测及群体遗传多样性分析等方面研究奠定基础。     

基于转录组测序的铁皮石斛黄酮代谢途径及相关基因解析

【目的】分析铁皮石斛Dendrobium officinale黄酮类化合物的生物合成途径及相关基因,为铁皮石斛黄酮类化合物代谢调控、药用价值的开发研究提供参考。【方法】利用Illumina HiSeq 4000测序平台对铁皮石斛2个生长阶段茎、叶进行高通量转录组测序,对组装获得的unigenes进行功能注释和黄酮类化合物的生物合成相关基因解析。【结果】铁皮石斛黄酮类化合物代谢相关Unigenes 48个,涉及14个酶。5个CHS相关Unigenes具有CHSlike保守结构域,活性位点氨基酸残基(Cys-His-Asn)高度保守,丙二酰辅酶A结合位点和产物结合位点的部分氨基酸残基发生变异。CHI(Unigene0013781)属于类型I查尔酮异构酶,异构化6′-羟基查尔酮为5-羟基黄烷酮。表达分析表明,2个生长时期的茎和叶,CHS(Unigene0008250)、 CHI(Unigene0013781)和F3H的表达量都高于CHS(Unigene0012884)和C3′H。【结论】通过对转录组数据分析,共获得铁皮石斛黄酮类化合物代谢相关Unigenes48个,涉及14个酶;5个CHS相关Unigenes中,2个编码查尔酮合酶,3个编码联苄合酶;CHI(Unigene0013781)属于类型I查尔酮异构酶。     

基于RNA-Seq技术的黄连根茎与须根转录组分析

采用RNA-Seq技术,对黄连根茎与须根进行转录组分析,经过trinity软件组装得到130 592条Unigenes,平均长度为700 bp,其中68 976条Unigenes在NR、GO、KEGG、eggNOG、Swiss-Prot和Pfam数据库获得基因注释信息,仅占全部Unigenes的52.82%。39 572个Unigene被注释到GO数据库的生物学进程、细胞组分和分子功能三大类的52个小类中;参与KEGG的5大类代谢通路,34个代谢路径,共搜索到24 999个SSR位点,单核苷酸、两核苷酸和三核苷酸为黄连根部SSR的主要重复基序。黄连根茎与须根的差异表达基因有13 496个,在黄连须根中上调的差异基因有8 375个,下调的有5 121个。两者差异基因GO富集分析表明,跟次生代谢相关进程的前5项为次生代谢进程、次生代谢物生物合成过程、生物碱代谢过程、异喹啉生物碱生物合成过程和异喹啉生物碱代谢过程。在细胞组分分类中,膜结构、细胞壁和外封装结构为差异基因富集的前三;在分子功能分类中,跟氧化还原酶活性相关条目占差异基因富集显著的前四。在KEGG数据库中,共有3 749个差异基因定位到125个代谢途径中,其中,24个差异基因注释到异喹啉生物碱生物合成,在关联的8条KEGG通路中,筛选得到11个关键酶基因;同时,有39个差异基因注释到络氨酸代谢途径中,在关联的21条KEGG通路中,筛选得到10种关键酶基因。黄连根部转录组数据的获得为研究黄连异喹啉类生物碱的合成及黄连根茎与须根代谢差异奠定了基础。     

罗氏沼虾转录组SSR标记信息分析

对罗氏沼虾肝胰腺组织转录组测序获得的33 450条Unigenes,进行简单重复序列(simple sequence repeat,简称SSR)检测与分析。结果共检测出15 356个SSR位点,分布在12 930条Unigenes上,发生频率为38.65%,平均距离为1.91 kb。其中,优势重复基序为单核苷酸、二核苷酸、三核苷酸。SSR序列重复长度主要为12～16 bp和18 bp,占罗氏沼虾总SSR的86.33%。单核苷酸重复基序中多数为A/T,二核苷酸重复基序以GA/TC、AG/CT及AT/TA为主,三核苷酸重复基序主要为AAT/TTA、AAG/TTC、ATT/TAA、CCT/GGA、CTC/GAG和CTT/GAA。最后,对SSR所在的Unigene进行了KOG分类注释。通过对罗氏沼虾转录组SSR位点信息分析,获得多态性较为丰富的SSR位点,具有极大的开发价值,可为罗氏沼虾遗传多样性的分析、遗传图谱的构建及分子标记辅助育种提供有效工具。     

基于RNA-Seq技术的茶树花转录组分析

本研究利用RNA-seq技术对父本、母本、子代不育茶树花三个样本花的转录组进行测定,经组装分析获得403 469条高质量的Unigenes。将获得的Unigenes与SWISS-PROT、TREMBL、CDD、PFAM、NR和KOG库进行blast,共注释到307 291条Unigenes。KOG功能分类显示,有23 739个Unigenes被分为25类。KEGG通路分析表明,共识别出26 967个Unigenes涉及的pathway有328个。SSR查找发现,从403 469个Unigenes中找到46 440个含有SSR序列。这些信息为茶树不育基因筛选、不育机理研究以及分子标记开发奠定了基础。     

黄芩高通量转录组测序数据组装和分析

采用高通量测序技术对黄芩的转录组进行测序,获得了29 099 899条reads数据,拼接后得到53 353条Unigene;将所获得的Unigene与COG,GO,KEGG,Swiss-Prot,NR这5个公共数据库进行比对,结果发现,分别有10 756,21 950,8 101,20 339,29 288条Unigene可比对到以上5个数据库中;已注释的Unigene与COG数据库比对后按功能可分为25类;根据GO功能可分为三大类57个亚类;经过与KEGG数据库比对后按照代谢通路可分为116类;利用Get ORF软件进行ORF预测,获得长度大于300 nt的ORF共20 552个;通过SSR分析,共获得5 658个SSR标记。获得的转录组信息可为今后进行黄芩分子标记的开发和关键基因的克隆及功能分析等研究提供基础数据。     

基于黄地老虎转录组测序的SSR和SNP特征分析

【目的】基于黄地老虎转录组测序结果,对其SSR及SNP位点信息进行分析。【方法】对黄地老虎总RNA进行提取,使用Illumina平台对转录组进行测序,利用MISA和GATK软件分别对黄地老虎转录组的SSR和SNP位点特征信息进行分析。【结果】黄地老虎转录组数据库包含66 469条Unigene序列。利用MISA对Unigene序列进行搜索,得到SSR位点4 438个,分布于4 048条Unigene序列上,占总序列的6.09%。其中,以单碱基和三碱基重复为主,分别占总重复类型的54.73%和29.29%。A/T基元是黄地老虎SSR的优势基元。SSR基元包含4～24次重复,以5次重复（21.67%）为主。SSR位点序列长度为12～127 bp,包含3 493个中等多态性位点、820个高度多态性位点。利用GATK对Unigene序列SNP位点进行搜索,得到转换类型（237 619个）和颠换类型（133 529个）共371 148个。C/T占SNP位点总数的18.61%,比率最高;G/A位居其次（18.28%）,均属于转换类型。转换类型（64.02%）显著高于颠换类型（35.98%）。【结论...     

中国濒危植物金钱松转录组测序 及生物信息学分析

金钱松是中国特有的孑遗单种属裸子植物,现存的自然种群数量很少,多被引种栽培,也是著名的庭院观赏树种。迄今为止,其遗传背景和基因组信息并不清楚,对于金钱松的保护及其遗传结构研究迫切需要基因组资源。采用Illumina HiSeqTM2500高通量测序平台对金钱松叶片进行转录组测序,经de novo组装共获得70 761条Unigene,平均长度为699 bp,N50的长度为1 300 bp,Q20和Q30序列分别占96.59%和91.29%。通过对7个不同的蛋白质和功能域数据库进行比对和功能注释,有43 674条Unigene(61.72%)注释成功。在GO数据库中,有28 355条Unigene按功能被划分成3大类56个小类,以执行生物过程的类区所占比例最多。通过KEGG pathway分析,有14 623条Unigene注释成功,发现了显著性富集的32条代谢通路,以代谢相关的基因最多。在KOG数据库中,有15 931条Unigene被分配到26个基因功能大类中,其中以参与一般功能、转录、翻译、修饰及蛋白运输的基因最为丰富。此外,利用MISA软件对转录组序列进行EST-SSR位点搜索与分析,共检测到2 260条Unigene含有2 462个EST-SSR位点,分布频率为3.48%,其中有180条序列含有一个以上EST-SSR位点,83条序列含有复合EST-SSR位点,以三核苷酸重复基元类型最为丰富,占42.53%(1 047个EST-SSR),重复次数主要以5～8次为主。这些重要的转录组序列为进一步了解金钱松生物学过程的分子机制提供了有价值的信息,并为未来的功能基因组分析、分子标记开发和群体遗传学分析提供了丰富的资源。     

荸荠优良品种‘桂蹄3号’球茎发育过程转录组测序分析

【目的】开展荸荠球茎发育过程转录组测序研究,为研究荸荠球茎发育过程相关基因表达信息提供参考。【方法】运用高通量测序技术,对荸荠球茎不同发育时期进行转录组测序研究。【结果】组装共得到223 182条转录本和90 542条Unigene,平均长度为809 bp,N50为1119。组装完整性较高,效果较好。对所得Unigene进行不同数据库注释,共有50 583条Unigene成功注释到7个数据库(NR、GO、NT、Pfam、KEGG、KOG及Swiss-prot)。对差异表达基因分析,发现球茎膨大初期的差异表达基因数目最多,为最活跃阶段。GO功能富集分析结果表明,33 205个基因获得功能注释,分为分子功能、细胞组分和生物学过程等3大类和54个亚类。COG功能分类结果表明,17 743个基因分布于25个功能区域,其中碳水化合物代谢占重要地位。KEGG代谢通路注释结果表明有20 667个基因获得功能注释,共有116条代谢途径,其中淀粉-蔗糖代谢占主要作用。【结论】利用高通量转录组测序技术首次建立了荸荠优良品种‘桂蹄3号’球茎的转录组数据库,为进一步研究荸荠球茎淀粉生物合成相关基因的功能及形成的分子机制提供了数据基础。     

青篱柴转录组的高通量测序及分析

青篱柴具有很好的观赏价值和药用价值,但其基因组信息相对缺乏,导致其分子生物学和基因功能的研究受到限制.为了获得和挖掘青篱柴的基因数据和功能,首次利用新一代高通量测序技术平台Illumina HiSeqTM2000对青篱柴转录组进行测序,经de novo组装后共得到66 755条单基因簇(Unigene).进一步利用7大公共数据库进行Blast同源比对,注释了50 339条Unigene.研究发现,558个基因参与了青篱柴次生代谢产物的生物合成,其中在甜菜红碱、芥子油苷、类黄酮、单菌霉素和苯丙素生物合成途径中的Unigene分别有2个、18个、44个、70个和231个,这些基因很可能参与了青篱柴药用活性物质的生物合成.SSR位点搜索发现,在18 972条Unigene中含有22 542个SSR位点.研究结果将为青篱柴功能基因的克隆、基因的表达、分子标记开发和遗传多样性研究奠定基础.同时,转录因子分析结果也为青篱柴的喀斯特环境适应性机制提供了初步线索.