• 中国中文核心期刊
  • 中国科学引文数据库（CSCD）核心库来源期刊
  • 中国科技论文统计源期刊（CJCR）
  • 第二届国家期刊奖提名奖

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

蕨类植物芒萁幼孢子体转录组高通量测序及特征分析

刘丽婷 ,  温强 ,  黄小春 ,  刘琪璟

downloadPDF
文章导航 > 林业科学研究  > 2016 >  29(4): 500-507
引用本文:
shu
Citation:
shu

蕨类植物芒萁幼孢子体转录组高通量测序及特征分析

  • 1.

    北京林业大学 北京 100083

  • 2.

    江西省林业科学院省植物生物技术重点实验室 江西 南昌 330013

  • 基金项目:

    国家林业局948项目（213-4-62）;江西省重大财政专项青年人才培养计划（芒萁根系微生物环境与成片发育研究）。

  • 中图分类号: S718.46

De novo Sequencing and Characterization of Juvenile Sporophyte Transcriptome of a Fern, Dicranopteris dichotoma

  • 1.

    Beijing Forestry University, Beijing 100083, China

  • 2.

    Jiangxi Provincial Biotech Key Lab for Plant, Jiangxi Academy of Forestry, Nanchang 330032, Jiangxi, China

  • CLC number: S718.46

  • 摘要: ［目的］ 采用高通量测序技术Illumina MiSeq250获得蕨类植物芒萁的孢子体转录组数据，以期为芒萁的生长、发育、代谢调控、微进化机制分析等提供重要的分子信息。［方法］ 应用生物信息学方法对测序获得的大量单基因簇（Unigene）进行基因功能注释、代谢途径及微卫星分析等。［结果］ 本研究共获得18 463 296条序列读取片段（reads），总碱基数为4.62Gbp序列信息，经序列组装最终得到63 169个Unigene，平均单条Unigene长度为863 bp，N50为1 587 bp，其中分布在200500 bp长度区间Unigene占总数的55.4%。数据库中的序列同源性比较表明，26 826个Unigene与其他物种的已知基因具有不同程度的同源性。芒萁转录组中的Unigene根据GO功能大致可分为细胞组成、分子功能和生物学过程3大类47个分支，其中有大量的Unigene与细胞进程、绑定活性、代谢过程和催化活性相关。将Unigene与COG数据库进行比对，根据其功能大致可分为26类。以KEGG数据库作为参考，依据代谢途径可将Unigene定位到257个代谢途径分支。此外，利用MISA软件检索26碱基微卫星，共找到13 286个SSR。在不同长度微卫星中，三核苷重复数量最多，占总数的40.41%。在各重复基序类型中出现频率最高的为AG/CT（14.45%）与AAG/CTT（12.39%）。利用重复基序开发的多态性SSR标记，可应用于芒萁不同个体的基因型分型鉴定。［结论］ 本研究获得了较高质量的芒萁转录组数据库，揭示了芒萁孢子体生长发育过程中表达基因的功能总体特征，可为芒萁进一步的功能基因挖掘和分子标记规模化开发奠定基础。
    Abstract: [Objective] The sporophyte transcriptome of Dicranopteris dichotoma was sequenced by Illumina MiSeq 250 to provide molecular information of its growth, development, metabolism, and the micro evolutionary mechanism. [Method] The functional annotations, metabolic pathways and microsatellite analysis of some Unigenes were conducted using bioinformatics methods. [Result] A total of 18 463 296 reads containing 4.62 Gbp of sequence information were generated. A total of 63 169 unigenes were formed by initial sequence splicing, with an average read length of 863 bp and N50 value of 1 587 bp. 26 826 unigenes were annotated using BLASTX searches against the Nr, Nt and SwissProt databases. The unigenes of the transcriptome of D. dichotoma were roughly divided into cellular components, molecule function and biological processes categories of 47 branches by gene ontology, of which related with cellular process cell, binding, metabolism processes and catalytic activities. Further annotated based on COG category, Unigenes could be grouped into 26 functional categories. KEGG pathway analysis showed that Unigenes could be divided into 276 classes based on their metabolic function. Meanwhile, 13 286 SSRs (simple sequence repeats) were mined with repeat motif of 2 to 6 bp by MISA. The trinucleotide repeats were most dominant, accounting for a total of 40.41%. AG/CT (14.45%) and AAG/CTT (12.39%) were the most common repeat motifs. Polymorphic SSR markers were developed from repeat motifs, which could be used for genotyping of different individuals of D. dichotoma. [Conclusion] A higher quality of transcriptome database was obtained in this study, which could reveal the general characteristics of gene expression in the process of growth and development, and lay the foundation for further gene function mining and the large-scale development of molecular markers of D. dichotoma.
  • [1]

    Xu X L, Li Q K, Wang J Y, et al. Inorganic and Organic Nitrogen Acquisition by a Fern Dicranopteris dichotoma in a Subtropical Forest in South China[J]. PLOS ONE, 2014, 9(5): e9005.
    [2] 董 鸣.克隆植物生态学[M]. 北京:科学出版社. 2011,5-6.
    [3] 岩槻邦男. 日本野生植物-蕨类(日文) [M]. 东京: 平凡社. 1999,311.
    [4] 张明如, 何 明, 温国胜, 等. 芒萁种群特征及其对森林更新影响评述[J]. 内蒙古农业大学学报, 2010, 31(4): 303-308.
    [5] 李小飞, 陈志彪, 陈志强, 等.南方红壤侵蚀区芒萁生长特征及其对环境因子的响应[J]. 水土保持通报, 2013, 33(3): 33-37.
    [6] 刘迎春, 刘琪璟, 汪宏清, 等.芒萁生物量分布特征[J].生态学杂志, 2008, 27(5): 705-711.
    [7] 侯晓龙, 刘明新, 蔡丽平, 等.安溪崩岗侵蚀区不同植被配置模式与恢复效果研究[J]. 亚热带水土保持, 2010, 22(1): 5-10.
    [8]

    Wei C Y, Wang C, Sun X, et al. Arsenic accumulation by ferns:a field survey in southern China[J]. Environmental Geochemistry and Health, 2007, 29(3):169-177.
    [9] 刘足根, 杨国华, 杨 帆, 等. 赣南钨矿区土壤重金属含量与植物富集特征[J].生态学杂志, 2008, 27(8): 1345-1350.
    [10] 李小飞, 陈志彪, 陈志强.南方稀土采矿恢复地土壤稀土元素含量及植物吸收特征[J]. 生态学杂志, 2013, 32( 8): 2126-2132.
    [11] 王立丰, 季红兵, 田维敏.重稀土矿区芒萁稀土元素精细地位及光抑制对其光合活性的影响[J]. 中国稀土学报, 2010, 28(3): 379-386.
    [12] 骆永明. 金属污染土壤的植物修复[J]. 土壤, 1999, 33(5): 261-265.
    [13] 李交昆, 龚育龙, 唐璐璐, 等. 金属型植物的研究进展[J]. 生命科学研究, 2011, 15(6):560-564.
    [14] 邓敏捷, 董焱鹏, 赵振利, 等. 基于Illumina 高通量测序的泡桐转录组研究[J]. 林业科学, 2013, 49(6):30-36.
    [15]

    Wang Z W, Jiang C, Wen Q, et al. Deep sequencing of the Camellia chekiangoleosa transcriptome revealed candidate genes for anthocyanin biosynthesis[J]. Gene, 2014, 538(1):1-7.
    [16]

    Barker M S, Wolf P G. Unfurling fern biology in the genomics age[J]. Bioscience, 2010, 60: 177-185.
    [17]

    Barker M S. Evolutionary genomic analyses of ferns reveal that high chromosome numbers are a product of high retention and fewer rounds of polyploidy relative to angiosperms[J]. Amer Fern J, 2009, 99: 136-137.
    [18]

    Nakazato T, Barker M S, Rieseberg L H, et al. Evolution of the nuclear genome of ferns and lycophytes [M]// Ranker T A, Haufler C H.Biology and Evolution of Ferns and Lycophytes. Cambridge University Press. 2008, 175-198.
    [19]

    Aya K, Kobayashi M, Tanaka J, et al. De novo transcriptome assembly of a fern, Lygodium japonicum, and a web resource database, Ljtrans DB[J]. Plant & Cell Physiology, 2015, 56(1):e5.
    [20]

    Weng J K, Tanurdzic M, Chapple C. Functional analysis and comparative genomics of expressed sequence tags from the lycophyte Selaginella moellendorffii[J]. BMC Genomics, 2005, 6: 85-97.
    [21]

    Yamauchi D, Sutoh K, Kanegae H, et al. Analysis of expressed sequence tags in prothallia of Adiantum capillus-veneris[J]. Journal of Plant Research, 2005, 118:223-227.
    [22]

    Der J D, Barker M S, Wickett N J, et al. De novo characterization of the gametophyte transcriptome in bracken fern, Pteridium aquilinum[J]. BMC Genomics, 2011, 12(1):99-113.
    [23]

    Bushart T J, Cannon A E, Haque U L A, et al. RNA-seq analysis identifies potential modulators of gravity response in spores of Ceratopteris (Parkeriaceae): evidence for modulation by calcium pumps and apyrase activity[J]. Amer J Bot, 2013, 100: 161-174.
    [24] 贾新平, 孙晓波, 邓衍明, 等. 鸟巢蕨转录组高通量测序及分析[J]. 园艺学报, 2014, 41(11): 2329-2341.
    [25] 温 强, 叶金山, 雷小林, 等. 油茶ISSR反应体系建立及优化[J]. 中南林学院学报, 2006, 26(6): 22-26.
    [26]

    Grabherr M G, Haas B J, Yassour M, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome[J]. Nature Biotechnology, 2011, 29, 644-652.
    [27]

    Weber J L. Infomativeness of human (dC-dA)n-(dG-dT)n polymorphisms[J]. Genomics, 1990, 7: 524-530.
    [28]

    Trapnell C, Williams B A, Pertea G, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation[J]. Nature Biotechnology, 2010, 28, 511-515.
    [29] 全先庆, 张洪涛, 单 雷, 等. 植物金属硫蛋白及其重金属解毒机制研究进展[J]. 遗传, 2006, 28(3): 375-382.
    [30] 刘洪亮, 郑丽明, 刘青青, 等. 非模式生物转录组研究[J]. 遗传, 2013, 35(8): 955-970.
    [31] 梁 烨, 陈双燕, 刘公社. 新一代测序技术在植物转录组研究中的应用[J]. 遗传, 2011, 33(12): 1317-1326.
    [32] 雷淑云, 张发起, Khan G, 等. 利用高通量测序分析青藏高原地区青杨的SSR和SNP特征[J]. 林业科学研究, 2015, 28(1): 37-43.
    [33] 温 强, 徐林初, 江香梅, 等. 基于454测序的油茶DNA序列微卫星观测与分析[J]. 林业科学, 2013, 49(8): 43-50.
  • [1] 刘雄芳 , 李太强 , 李正红 , 万友名 , 刘秀贤 , 张序 , 安静 , 马宏 . 云南干热河谷地区余甘子转录组分析. 林业科学研究, 2018, 31(5): 1-8. doi: 10.13275/j.cnki.lykxyj.2018.05.001
    [2] 李太强 , 刘雄芳 , 万友名 , 李正红 , 李钰莹 , 刘秀贤 , 马宏 . 滇东南濒危植物长梗杜鹃转录组微卫星特征分析. 林业科学研究, 2017, 30(4): 533-541. doi: 10.13275/j.cnki.lykxyj.2017.04.001
    [3] 杜明凤 , 丁贵杰 . 基于马尾松干旱转录组的抗旱功能SSR位点分析. 林业科学研究, 2018, 31(5): 9-19. doi: 10.13275/j.cnki.lykxyj.2018.05.002
    [4] 陈奶寿 , 张秋芳 , 陈坦 , 杨玉盛 , 谢锦升 . 退化红壤恢复过程中芒萁的N、P化学计量特征. 林业科学研究, 2016, 29(5): 735-742.
    [5] 董国鑫 , 黄雪蔓 , 王一 , 任立宁 . 两种芒萁覆盖度下毛竹林土壤酶活性及生态化学计量特征比较研究. 林业科学研究, 2024, 37(1): 183-193. doi: 10.12403/j.1001-1498.20230137
    [6] 王书珍 , 张羽佳 , 黄诗颖 , 罗炎炎 , 金正强 , 李志良 , 金卫斌 . 基于锦绣杜鹃花蕾转录组的SSR标记开发及应用. 林业科学研究, 2019, 32(3): 97-104. doi: 10.13275/j.cnki.lykxyj.2019.03.013
    [7] 贺君星 , 马庆国 , 裴东 , 张俊佩 . 基于核桃参考基因组的SSR位点鉴定分析和单态性标记开发. 林业科学研究, 2022, 35(6): 89-100. doi: 10.13275/j.cnki.lykxyj.2022.006.010
    [8] 王帅 , 邵芬娟 , 李论 , 芦强 , 邱德有 . 穗花杉的转录组测序及其转录组特性分析. 林业科学研究, 2017, 30(5): 759-764. doi: 10.13275/j.cnki.lykxyj.2017.05.008
    [9] 马婧 , 邓楠 , 褚建民 , 纪敬 , 史胜青 , 江泽平 , 成铁龙 . 泡泡刺高通量转录组鉴定及其黄酮类代谢途径初步分析. 林业科学研究, 2016, 29(1): 61-66.
    [10] 饶龙兵 , 杨汉波 , 郭洪英 , 段红平 , 陈益泰 . 基于桤木属转录组测序的SSR分子标记的开发. 林业科学研究, 2016, 29(6): 875-882.
    [11] 郝心怡 , 王哲舒 , 范志斌 , 王丽娟 . 杨树响应腐皮镰刀菌侵染的转录组学分析. 林业科学研究, 2024, 37(): 1-13. doi: 10.12403/j.1001-1498.20230458
    [12] 李珊珊 , 曾艳飞 , 何彩云 , 张建国 . 基于沙棘转录组序列开发EST-SSR分子标记. 林业科学研究, 2017, 30(1): 69-74. doi: 10.13275/j.cnki.lykxyj.2017.01.0010
    [13] 张振 , 张含国 , 周宇 , 张磊 , 于宏影 , 张莉 . 红松4个组织的转录组数据分析与次生代谢产物的表达差异初探. 林业科学研究, 2015, 28(4): 597-603.
    [14] 封润霞 , 赵婕 , 张苏芳 , 王建军 , 魏建荣 , 刘建凤 . 绒毛白蜡（Fraxinus velutina Torr）韧皮部响应白蜡窄吉丁（Agrilus planipennis Fairmaire）危害的转录组变化. 林业科学研究, 2021, 34(1): 47-55. doi: 10.13275/j.cnki.lykxyj.2021.01.006
    [15] 邓楠 , 史胜青 , 常二梅 , 刘建锋 , 兰倩 , 江泽平 . 基于中麻黄萌发种子转录组的黄酮类化合物合成途径基因的挖掘. 林业科学研究, 2014, 27(6): 758-763.
    [16] 于淑惠 , 亓倩 , 孙涛 , 王雪庆 , 杨璞 , 冯颖 . 白蜡虫雄虫真蛹转录组分析. 林业科学研究, 2016, 29(3): 413-417.
    [17] 田林 , 尹丹丹 , 成铁龙 , 夏新莉 , 尹伟伦 . 盐胁迫下比拉底白刺差异表达基因的转录组分析. 林业科学研究, 2020, 33(1): 1-10. doi: 10.13275/j.cnki.lykxyj.2020.01.001
    [18] 陈沫 , 何沙娥 , 陈少雄 , 欧阳林男 , 张程 , 张维耀 . 基于转录水平解析尾巨桉径向生长对种植密度的响应. 林业科学研究, 2021, 34(5): 1-12. doi: 10.13275/j.cnki.lykxyj.2021.005.001
    [19] 刘文汝 , 冀胜鑫 , 甄志先 . 感染国槐带化病植株内源激素及转录组分析. 林业科学研究, 2022, 35(1): 141-149. doi: 10.13275/j.cnki.lykxyj.2022.01.016
    [20] 郭晋敏 , 杨升 , 刘星 , 王金旺 , 王文卿 , 陈秋夏 . 秋茄低温胁迫转录组分析及脱落酸信号途径基因挖掘. 林业科学研究, 2023, 36(2): 39-49. doi: 10.12403/j.1001-1498.20220417
  • 加载中
WeChat 点击查看大图
计量
  • 文章访问数:  3019
  • HTML全文浏览量:  407
  • PDF下载量:  921
  • 被引次数: 0
出版历程
  • 收稿日期:  2015-11-25

蕨类植物芒萁幼孢子体转录组高通量测序及特征分析

  • 1. 北京林业大学 北京 100083
  • 2. 江西省林业科学院省植物生物技术重点实验室 江西 南昌 330013
  • 收稿日期: 2015-11-25
基金项目:  国家林业局948项目（213-4-62）;江西省重大财政专项青年人才培养计划（芒萁根系微生物环境与成片发育研究）。

摘要: ［目的］ 采用高通量测序技术Illumina MiSeq250获得蕨类植物芒萁的孢子体转录组数据，以期为芒萁的生长、发育、代谢调控、微进化机制分析等提供重要的分子信息。［方法］ 应用生物信息学方法对测序获得的大量单基因簇（Unigene）进行基因功能注释、代谢途径及微卫星分析等。［结果］ 本研究共获得18 463 296条序列读取片段（reads），总碱基数为4.62Gbp序列信息，经序列组装最终得到63 169个Unigene，平均单条Unigene长度为863 bp，N50为1 587 bp，其中分布在200500 bp长度区间Unigene占总数的55.4%。数据库中的序列同源性比较表明，26 826个Unigene与其他物种的已知基因具有不同程度的同源性。芒萁转录组中的Unigene根据GO功能大致可分为细胞组成、分子功能和生物学过程3大类47个分支，其中有大量的Unigene与细胞进程、绑定活性、代谢过程和催化活性相关。将Unigene与COG数据库进行比对，根据其功能大致可分为26类。以KEGG数据库作为参考，依据代谢途径可将Unigene定位到257个代谢途径分支。此外，利用MISA软件检索26碱基微卫星，共找到13 286个SSR。在不同长度微卫星中，三核苷重复数量最多，占总数的40.41%。在各重复基序类型中出现频率最高的为AG/CT（14.45%）与AAG/CTT（12.39%）。利用重复基序开发的多态性SSR标记，可应用于芒萁不同个体的基因型分型鉴定。［结论］ 本研究获得了较高质量的芒萁转录组数据库，揭示了芒萁孢子体生长发育过程中表达基因的功能总体特征，可为芒萁进一步的功能基因挖掘和分子标记规模化开发奠定基础。

English Abstract

    全文HTML

参考文献 (33)

目录

    /

    返回文章
    返回

    京ICP备09064894号-3

    版权所有：《林业科学研究》编辑部

    地址：北京万寿山后中国林科院电话：010-62889680；62889702Email：lykxyj@caf.ac.cn

    本系统由 北京仁和汇智信息技术有限公司 设计开发 技术支持： info@rhhz.net   百度统计