木本高山花卉马缨杜鹃的全基因组扫描分析

本研究采用高通量测序技术(Illumina Hiseq ~(TM)2000)首次对世界名贵园林花卉马缨杜鹃(Rhododendron delavayi)进行了基因组大小的测定,并评估了马缨杜鹃的GC含量、杂合率和序列重复度等信息参数,为马缨杜鹃全基因组测序策略的选择提供了有据参考。全基因组扫描结果如下:马缨杜鹃的基因组大小约为698 Mb,杂合率较高约为1.25%,GC含量约为39%,且具有一定的基因组序列重复度。研究表明:马缨杜鹃的全基因组测序分析不宜采用全基因组鸟枪法(WGS),可尝试WGS与fosmid或BAC-to-BAC相结合的测序策略,亦或采用杂合组装软件进行序列拼装,以提高全基因组序列的拼装质量。     

崖壁植物太行菊与长裂太行菊全基因组大小及特征分析北大核心CSCD

太行菊(Opisthopappus taihangensis)、长裂太行菊(O.longilobus),为太行山特有多年生崖壁草本植物,菊科(Compositae)重要野生资源,具有较高的经济与生态价值。为确定适合两物种的全基因组测序策略,该研究利用流式细胞法和高通量测序技术,分析两物种基因组大小、杂合率、重复序列及GC含量等信息。结果表明:(1)流式细胞法估算太行菊基因组大小约为2.1 Gb,长裂太行菊基因组大小约为2.4 Gb。(2)高通量测序修正后太行菊基因组大小为3.13 Gb,重复序列比例为84.35%,杂合度为0.99%,GC含量为36.56%;长裂太行菊基因组为3.18 Gb,重复序列比例为83.83%,杂合度为1.17%,GC含量为36.62%。(3)初步组装后GC含量分布及平均深度存在异常,出现分层现象,可能是两物种基因组杂合率较高所致。综上结果表明,太行菊、长裂太行菊均属于高重复、高杂合、大基因组的复杂基因组,建议使用Illumina+PacBio测序组装策略,进行全基因组测序分析。     

樟树全基因组调查

樟树是我国特有的珍贵用材和经济树种,富含多糖多酚、萜类等次生代谢物质,是香精香料、油脂化工和医药等的重要原料树种。本研究采用高通量测序技术(Illumina HiseqTM2000)首次测定了樟树基因组大小,并利用生物信息方法估计樟树杂合率、重复序列情况和GC含量等基因组信息,为全基因组测序策略的选择提供依据。主要结论如下:(1)樟树基因组大小粗略估计为760 Mb左右;(2)樟树基因组有较高的杂合率和一定的重复,杂合率约为0.65%;(3)由于樟树杂合率较高,全基因组鸟枪法策略不适合该基因组测序分析,可尝试使用BAC-to-BAC策略或fosmid策略,有利于樟树基因组的序列拼接和组装。     

罗汉果全基因组Survey分析

罗汉果是广西特有药用及甜料植物,其主要成分之一甜苷V作为天然、非糖甜味剂,具有广阔的开发前景,但罗汉果目前完全来自于栽培,适生区狭窄,连作障碍严重,加之含量低导致甜苷V生产成本居高不下,严重限制了其应用。为了减少盲目性,在大规模全基因组深度测序之前,先做低覆盖度的基因组Survey测序,评价基因组的大小及复杂程度,以确定适合该植物全基因组的测序研究策略。该研究采用第二代高通量测序技术(Illumina Hiseq TM 2000)首次测定了罗汉果基因组大小,并利用生物信息学方法估计罗汉果杂合率、重复序列和GC含量等基因组信息。结果表明:(1)获得了18.1 Gb罗汉果基因组测序数据,基因组大小估计为344.95 Mb左右,测序深度为52×;(2)从K-mer分布曲线发现罗汉果基因组有明显的杂合峰,杂合率达1.5%,基因组高杂合导致组装的结果中Contig N50和Scaffold N50的长度比预期的要短很多,还造成GC平均深度及含量分布明显异常,存在一个低深度分布区域。基因组主峰后面有微弱的重复峰,说明罗汉果存在较多的重复序列;(3)由于罗汉果存在高杂合率和重复序列较多的特点,该基因组测序分析仅采用全基因组鸟枪法(WGS)策略不合适,为了更好地对全基因组进行序列拼接和组装,可尝试结合采用Fosmid-to-Fosmid或BAC-to-BAC策略。该研究结果对于揭示罗汉果产量、有效成分含量、发育及抗病虫的分子机制,以及通过分子育种来提高甜苷V含量和降低生产成本具有重要意义,为全基因组测序策略的选择提供了依据。     

宽边黄粉蝶(Eurema hecabe)基因组Survey分析

宽边黄粉蝶(Eurema hecabe)是重要的传粉昆虫,广布于非洲热带区、东洋区、澳洲区及古北区东部,具有较高的学术价值和经济意义。为了深入研究其广泛适应性,确定适合该蝴蝶的全基因组的测序研究策略,我们首先做低覆盖度的基因组Survey测序,然后做大规模的全基因组深度测序。采用第二代高通量的测序技术作为该研究的研究方法,测定了宽边黄粉蝶基因组大小,并利用生物信息学方法估计该种的杂合率、重复序列和GC含量等基因组信息。结果表明:(1)宽边黄粉蝶的基因组大小估计为285.34 Mb,测序深度51×;(2)从K-mer分布曲线发现黄粉蝶基因组有明显的杂合峰,杂合率达1.97%,重复序列比例为35.37%。该研究结果对于揭示宽边黄粉蝶物种的起源和进化及适应性具有重要意义,为宽边黄粉蝶选择全基因组测序策略提供依据。     

甘薯属耐盐植物马鞍藤基因组大小及特征分析

马鞍藤(Ipomoea pes-caprae(L.)R. Br.)是甘薯的近缘野生种之一,为热带、亚热带海滩生长的多年生藤蔓植物,作为沿海地区的园林绿化植物之一,具有较强的耐盐性。马鞍藤基因组学的研究可为耐盐甘薯种质创制提供信息参考。本试验调查了马鞍藤基因组大小和特征概况,为后续全基因组精细图谱绘制打下基础。研究以已知基因组大小的三裂叶薯(Ipomoea triloba L.)为对照,运用流式细胞术对马鞍藤基因组大小进行初步预估;使用二代高通量测序技术(Illumina Hiseq2500)对马鞍藤基因组进行survey评估,测序深度30×,利用生物信息学方法估算马鞍藤GC含量(即鸟嘌呤和胞嘧啶所占比例)、杂合率、重复序列等基因组概况。结果表明:流式细胞术估算马鞍藤基因组大小为1012.704±17.37 Mb;经全基因组Survey测定获得马鞍藤有效数据为21.71 Gb,基因组大小经修正后估算为1041.65 Mbp;通过K-mer分布曲线估算马鞍藤基因组中重复序列所占比率为74.52%,杂合率为0.99%;经初步组装后,GC平均深度及含量分布存在异常,出现分层的现象,这可能与马鞍藤基因组杂合率较高有关。本试验首次报导甘薯属耐盐植物马鞍藤基因组大小及特征信息,为马鞍藤进一步的全基因组深度测序和甘薯耐盐基因挖掘打下基础。     

基于流式细胞术和K-mer分析的银缕梅属(Parrotia C.A.Mey.)植物基因组大小测定

金缕梅科(Hamamelidaceae)银缕梅属(Parrotia C.A.Mey.)仅包含银缕梅(Parrotia subaequalis(H.T.Chang)R.M.Hao&H.T.Wei)和波斯铁木(Parrotia persica(DC.)C.A.Mey.)两种落叶阔叶乔木,其中银缕梅是我国华东地区特有的Ⅰ级濒危珍稀保护植物,属东亚第三纪孑遗成分;其姊妹种波斯铁木则间断分布于伊朗北部,属北极第三纪孑遗植物类群。本研究首次利用流式细胞术和K-mer分析方法对银缕梅属两姊妹种的基因组大小进行了测定,建立和优化了以萝卜(Raphanus sativus L.‘Saxa’)为内标、WPB(Woody plant buffer)为细胞核解离液的两种植物单倍体基因组的DNA含量(DNA C值)流式测定的适宜体系,旨在为金缕梅科银缕梅属植物的全基因组测序、基因组学研究、种质资源开发和利用以及物种保育等提供前期基础数据参考;同时也可为金缕梅科其他属、种的基因组大小测定提供借鉴。主要研究结果如下:(1)通过流式测定银缕梅基因组大小约为971.45±13.91 Mb,波斯铁木基因组大小约为890.52±24.69 Mb;(2)K-mer分析估测银缕梅基因组大小为951.70 Mb,杂合率为1.740%,重复序列比例为77.50%;波斯铁木基因组大小为858.50 Mb,杂合率为0.695%,重复序列占74.30%;(3)银缕梅属于高杂合和高重复基因组,波斯铁木则属于微杂合和高重复基因组。本研究的结果为银缕梅属植物后续基于DNA三代高通量测序技术的全基因组测序、组装及去冗余处理等工作提供了重要的数据参考。     

蒙古沙冬青基因组调查及SSR分析

蒙古沙冬青是第三纪冰川孑遗物种,具有极强的耐逆性,是研究植物耐逆分子机制的良好材料。利用高通量测序技术和生物信息学工具预测蒙古沙冬青的基因组大小及杂合度,进行SSR分子遗传标记的初步鉴定。蒙古沙冬青基因组大小约为812 Mb,杂合率为0.506%,基因组杂合度较高。对基因组进行拼接,进行SSR分子遗传标记分析,共鉴定183 102个SSR,不同类型核苷酸重复差异较大,其中,二核苷酸重复中的AT/TA含量最高,四核苷酸重复最少,共占总数的1.48%。     

金针菇基因组测序与组装策略分析

采用二代和三代测序技术分别对金针菇单核体菌株“6-3”进行测序,应用4种组装策略进行基因组的de novo组装,对比组装效果。基因组组装的参数方面,仅使用二代测序组装的效果最差,长度大于10kb的Contig全长只有24.6Mb,Contig N50只有23kb,组装率只有59.27%。采用三代组装二代校正的组装策略效果最好,长度大于10kb的Contig全长为38.3Mb,Contig N50为2.8Mb,组装率高达92.16%。保守单拷贝基因拼接效果方面,4种组装策略获得基因组序列与BUSCO数据库里的担子菌的保守单拷贝基因比对,基因完整性均大于94%。在组装准确性方面,经过PCR扩增、Sanger测序验证,三代组装二代校正的基因组序列完整并且连续,同时序列上碱基的SNP、InDel数量最少。综上所述,三代组装二代校正得到的基因组序列具有Contig N50值大、组装率高、碱基准确性高的特点,是食用菌基因组测序较为理想的方案。     

雄烯二酮高产菌新金分枝杆菌MN4的全基因组测序及序列分析北大核心CSCD

【目的】新金分枝杆菌(Mycobacterium neoaurum)MN4是1株经诱变育种获得的高产雄烯二酮,并且能够耐受高浓度底物植物甾醇的突变菌株。为深入研究MN4菌株耐底物的机制及雄烯二酮的生物合成途径,有必要解析MN4菌株的全基因组序列信息。【方法】本研究采用高通量测序技术对MN4进行全基因组测序,然后使用相关软件对测序数据进行基因组装、基因预测与功能注释、COG聚类分析以及次级代谢产物合成基因簇预测等。【结果】新金分枝杆菌MN4基因组装获得33个Contigs,整个基因组大小为5.39 Mb,GC含量为66.9%,编码4920个蛋白基因,序列提交至Gen Bank数据库,登录号为JXYZ00000000。【结论】本研究首次报道了1株高产雄烯二酮菌株MN4的全基因组序列,分析了基因组的基本特征,初步解析了该菌株降解植物甾醇生产雄烯二酮的关键基因,将为MN4的功能基因组学研究及相关次级代谢产物的生物合成途径与异源表达研究提供基础。     

西北干旱区孑遗濒危植物四合木(Tetraena mongolica)叶绿体基因组特征研究及比较分析

四合木(Tetraena mongolica)是我国特有的蒺藜科(Zygophyllaceae)强旱生小灌木,因其起源古老、抗逆性强,所以可作为生物多样性起源和环境演变研究的理想对象,具有重要的学术研究价值。本研究采用Illumina双末端测序技术对四合木叶绿体基因组进行建库测序和分析。选取蒺藜目及牻牛儿苗目共计30个物种叶绿体基因组,与四合木进行系统发育关系分析探讨。结果表明:四合木叶绿体基因组长度为106259bp,其中反向重复区(IR区中)有7种基因,包括4种PCG基因,3种tRNA基因。叶绿体基因组共编码98种基因,包括65种蛋白编码基因、29种tRNA基因与4种rRNA基因。生物信息学表明,在四合木中共搜到92个SSR位点,其中包括74个单核苷酸重复基序,7个二核苷酸重复基序,1个三核苷酸重复基序,9个四核苷酸重复基序和1个五核苷酸基序。没有发现六核苷酸,其中单核苷酸重复在四合木的叶绿体基因组SSR中占比为80.1%。通过MEGA软件采用近邻结合法(neighbor-joining,NJ)对四合木等31个物种的叶绿体基因组进行聚类分析,发现四合木与蒺藜科三齿拉雷亚灌木为最近的姐妹种,其次为牻牛儿苗科智利白桦植物亲缘关系较近,与牻牛儿苗科天竺葵属和牻牛儿苗科高桂花属亲缘关系最远,说明四合木属于蒺藜科物种,这对于四合木的研究等具有一定的参考价值。     

卡奇霉素产生菌棘孢小单胞菌ATCC 15837的全基因组测序及序列分析

棘孢小单胞菌(Micromonospora echinospora) ATCC 15837是一种高GC含量的革兰氏阳性稀有放线菌,能够合成烯二炔类抗肿瘤抗生素卡奇霉素(calicheamicin, CLM)。目前,还没有相关研究报道棘孢小单胞菌ATCC 15837的全基因组序列,这限制了其代谢产物合成途径和比较基因组学等研究。本研究首次通过高通量测序技术对棘孢小单胞菌ATCC 15837进行全基因组测序,使用相关生物信息学软件对数据进行组装和注释等分析。使用Velvet软件进行组装拼接得到77个Contigs,GC含量为72.36%,基因组大小约为7.69 Mb。序列已提交至美国国立生物技术信息中心(NCBI)的GenBank数据库(登录号为NGNT00000000)。本研究首次报道了一株烯二炔类抗肿瘤抗生素卡奇霉素产生菌棘孢小单胞菌ATCC 15837的全基因组序列,分析了基因组基本特征,预测了该菌株的次级代谢产物生物合成基因簇,为后续的进一步代谢调控与合成生物学提供了理论基础。     

膨大剂处理对樟子松种实性状的影响

樟子松(Pinus sylvestris var.mongolica)是我国北方重要的用材、园林绿化及生态防护树种。研究膨大剂处理对樟子松种实性状的影响,可以为樟子松优质球果的高产、稳产的科学管理方式提供依据。本研究以吉林省延边自治区汪清县亲和种子园樟子松优良单株为研究对象,分析氯吡脲(PPU)和α-萘乙酸钠(α-NA)2种膨大剂处理对樟子松球果产量、球果质量、种子产量和种子质量等指标的影响。膨大剂处理后,樟子松种实性状的平均值均大于对照,表明膨大剂处理对于樟子松的种实发育有促进作用;樟子松的球果重、单果出籽量、种长及种宽在各处理间差异均达极显著水平(P<0.001),表明膨大剂处理对樟子松球果产量以及种子质量有极显著影响。膨大剂处理后单果出籽量和单果出籽重的表型变异系数分别达到28.54%、33.88%。利用球果性状对各处理进行综合评价,初步选出氯吡脲浓度10 mg·L^-1为最佳处理,该处理条件下球果长、球果宽、球果重的现实增益分别为11.81%、4.65%和24.35%;根据种子性状初步选出α-萘乙酸钠质量浓度为10 mg·L^-1的...     

Whole-genome survey and phylogenetic analysis of Gadus macrocephalus

Gadus macrocephalus (Pacific cod) is an economically important species on the northern coast of the Pacific. Although numerous studies on G. macrocephalus exist, there are few reports on its genomic data. Here, we used whole-genome sequencing data to elucidate the genomic characteristics and phylogenetic relationship of G. macrocephalus. From the 19-mer frequency distribution, the genome size was estimated to be 658.22 Mb. The heterozygosity, repetitive sequence content and GC content were approximately 0.62%, 27.50% and 44.73%, respectively. The draft genome sequences were initially assembled, yielding a total of 500,760 scaffolds (N50 = 3565 bp). A total of 789,860 microsatellite motifs were identified from the genomic data, and dinucleotide repeat was the most dominant simple sequence repeat motif. As a byproduct of whole-genome sequencing, the mitochondrial genome was assembled to investigate the evolutionary relationships between G. macrocephalus and its relatives. On the basis of 13 protein-coding gene sequences of the mitochondrial genome of Gadidae species, the maximum likelihood phylogenetic tree showed that complicated relationships and divergence times among Gadidae species. Demographic history analysis revealed changes in the G. macrocephalus population during the Pleistocene by using the pairwise sequentially Markovian coalescent model. These findings supplement the genomic data of G. macrocephalus, and make a valuable contribution to the whole-genome studies on G. macrocephalus.     

喜马红景天叶绿体基因组特征及其系统发育分析

以青藏高原特有药用植物——喜马红景天（Rhodiola himalensis）为试验材料,利用高通量测序技术对喜马红景天进行叶绿体基因组测序、组装和注释,获得完整的叶绿体基因组。结果显示：喜马红景天叶绿体基因组全长为151 074 bp,GC含量为37.8%,具有1个长单拷贝区、1个短单拷贝区和1对反向重复区的典型四分体结构,其序列长度分别为82 309、17 017、25 874 bp;叶绿体基因组共编码130个基因,其中编码蛋白的基因86个、编码tRNA的基因37个、编码rRNA的基因7个;叶绿体基因组共检测出25 513个密码子,其中编码亮氨酸（Leu）的密码子占比最大;喜马红景天IRa和IRb区的rps19和ycf1基因缺失,长单拷贝区的trnH基因收缩;喜马红景天与圣地红景天（R. sacra）亲缘关系最近;短单拷贝区域的单核苷酸多态性（SNP）变异频率最高。本研究报道了喜马红景天的叶绿体基因组,并对其进行了组装、注释和序列分析,为今后开展喜马红景天的遗传多样性研究和合理开发利用提供理论依据。     

鲁桑叶绿体基因组序列及特征分析

鲁桑(Morus multicaulis)是亚洲地区栽培的重要经济作物。以鲁桑品种日本胡橙为实验材料, 利用高通量测序技术对鲁桑叶绿体基因组进行测序, 获得NCBI登录号(KU355297), 并研究鲁桑的叶绿体基因组结构。结合前人对蒙桑(M. mongolica)、印度桑(M. indica)和川桑(M. notabilis)的研究结果, 对鲁桑的系统进化关系进行了探讨。研究结果表明: 鲁桑叶绿体基因组是一个典型的四部分结构, 全长159 154 bp, 共注释130个基因, 包含85个蛋白质编码基因(18个基因在反向重复区重复)、37个转运RNA (tRNA)基因和8个核糖体RNA (rRNA)基因。生物信息学分析表明, 在鲁桑中共搜索到82个SSR位点, 单核苷酸、二核苷酸、三核苷酸、四核苷酸和五核苷酸重复基序个数分别为63、7、2、9和1个, 并没有发现六核苷酸; 其中单核苷酸重复在鲁桑的叶绿体基因组SSR中占76.8%。采用MEGA 6.0软件, 通过最大似然法和近邻结合法对包括4个桑属物种在内的15个物种的叶绿体基因组序列进行聚类分析, 2种方法得到的聚类结果均为鲁桑和蒙桑聚在一起。研究结果对叶绿体基因组工程研究及桑属种间的分子标记开发和优良品种培育具有一定的参考价值。     

六妹羊肚菌PacBio基因组测序和DNA 6mA甲基化分析

羊肚菌属于子囊门真菌,是世界范围内最受欢迎的食用菌之一。本研究通过PacBio单分子实时测序技术对我国四川省成功栽培的六妹羊肚菌Morchella sextelata SCLS菌株进行全基因组测序,获得其高质量核基因组组装,大小为53.57Mb,重复序列含量为17.03%,包含42条重叠群（contigs）,重叠群N50高达1.82Mb,其中13条重叠群两端均含有端粒重复序列,为完整的染色体。通过链特异性RNA-seq测序和转录本拼接,并结合多种基因预测策略,预测到13 182个蛋白编码基因,包括267个碳水化合物活性酶,11个次生代谢产物合成基因簇。通过与内蒙古地区栽培的六妹羊肚菌菌株NZTD180501373基因组比较发现,六妹羊肚菌进化过程中可能发生过染色体重组事件,二者具有33 055个SNP和48 726个InDel位点差异,并且各自拥有超过6Mb的特有序列。此外,相比于梯棱羊肚菌M. importuna,六妹羊肚菌SCLS菌株中与逆转录转座酶有关的orthogroup发生了扩张,并且拥有5个成员数超过100的特有逆转录转座酶基因家族。对SCLS菌株中的DNA N⁶腺嘌呤甲基化（6mA）修饰进行了鉴定,发现SCLS菌株基因组中0.42%的腺嘌呤被6mA甲基化,其含量显著高于已报道的6种双核亚界（Dikarya）真菌。6mA甲基化位点在逆转录转座子上显著富集,表明六妹羊肚菌中6mA甲基化可能调控逆转录转座子的活性,这也是首次在真菌中报道6mA甲基化与转座子相关。