全基因组重测序分析4个毛竹变型的秆形和秆色变异

doi:10.7668/hbnxb.20191110

摘要/Abstract

摘要： 为揭示毛竹重要变型的全基因组突变类型，选取4个有代表性的毛竹变型，采用高通量重测序技术，进行全基因组重测序，测序深度12×。分别提取毛竹4个变型与参考基因组间发生非同义突变的SNP、CDS区发生InDel与SV的基因，比较可以发现，2个秆形变异的竹种圣音竹和厚壁毛竹基因组的变异基因数目和变异类型相近，而2个秆色变异的竹种黄槽毛竹和黄皮花毛竹基因组的变异数目和变异类型相近。4个毛竹变型均存在多个基因同时存在多种类型突变。将变异基因与GO、COG、KEGG等功能数据库进行比对、注释，共获得1 739个参与细胞壁、细胞膜合成的基因，501个细胞骨架构建相关基因，1 785个转录因子相关基因，1 324个信号转导相关基因，104个赤霉素相关基因，76个激素代谢相关基因，1 938个叶绿素合成相关的基因，73个类胡萝卜素合成代谢相关基因和68个花青素合成调控基因，这些差异基因分别在毛竹秆形和秆色调控方面起着重要作用。

关键词: 毛竹, 全基因组重测序, 基因变异, 单核苷酸多态性, 结构变异, 小片段插入缺失

Abstract: For an overall understanding of the whole genome of 4 representative variations,re-sequencing was used for high-throughput sequencing with 12 depth to detect its variations. By comparing the genes with non-synonymous SNP,InDel and SV in CDS,the variation numbers and types of two variations on culm shape, Phyllostachys edulis f. tubaeformis and Phyllostachys edulis f. pachyloen were similar,and the other two variations on culm color, Phyllostachys edulis f. luteosulcata and Phyllostachys edulis f. huamozhu were similar. And multiple genes occurred multiple types of mutations simultaneously. By comparing and annotating the mutant genes with GO,COG,KEGG and other functional databases,a total of 1 739 genes involved in cell wall and cell membrane synthesis,501 cytoskeletal construction-related genes,1 785 transcription factors,1 324 signal transduction genes,104 gibberellins and 76 hormone genes,1 938 chlorophyll synthesis related genes,73 carotenoids synthetic metabolism related genes and 68 anthocyanin synthesis regulation genes were obtained. These genes play an important role in the regulation of bamboo culm shape and color.

Key words: Moso bamboo, Whole genome re-sequencing, DNA variation, Single nucleotide polymorphism, Structure variation, Small indel

中图分类号:

S795.7
Q78

牟少华, 李雪平, 李娟, 高健. 全基因组重测序分析4个毛竹变型的秆形和秆色变异[J]. 华北农学报, 2020, 35(4): 96-105. doi: 10.7668/hbnxb.20191110.

MU Shaohua, LI Xueping, LI Juan, GAO Jian. Analysis on the Culm Shape and Color Variation of 4 Moso Bamboo Variants by Whole Genome Re-sequencing[J]. ACTA AGRICULTURAE BOREALI-SINICA, 2020, 35(4): 96-105. doi: 10.7668/hbnxb.20191110.

http://www.hbnxb.net/CN/Y2020/V35/I4/96

参考文献

[1] 王嘉,曾召琼,梁建秋,于晓波,吴海英,张明荣. 基于全基因组重测序的大豆分子标记开发及籽粒蛋白质含量QTL定位[J]. 中国农业科学,2019,52(16):2743-2757. doi:10.3864/j.issn.0578-1752.2019.16. 001. Wang J,Zeng Z Q,Liang J Q,Yu X B,Wu H Y,Zhang M R. Development new molecular markers for quantitative trait locus (QTL) analysis of the seed protein content based on whole genome re-sequencing in Soybean[J]. Scientia Agricultura Sinica,2019,52(16):2743-2757.
[2] 贾小平, 张博, 全建章, 李剑峰, 王永芳, 袁玺垒. 不同光周期条件下谷子株高的全基因组关联分析[J]. 华北农学报,2019,34(4):16-23. doi:10.7668/hbnxb.201750856. Jia X P,Zhang B,Quan J Z,Li J F,Wang Y F,Yuan X L. Genome-wide association analysis of plant height in foxtail millet under different photoperiod conditions[J]. Acta Agriculturae Boreali-Sinica,2019,34(4):16-23.
[3] 张羽,François Belzile. 大豆抗菌核病的全基因组关联研究[J]. 华北农学报,2020,35(1):205-213.doi:10.7668/hbnxb.20190364. Zhang Y,François B. Genome-wide association study for Sclerotinia sclerotiorum resistance of soybean[J]. Acta Agriculturae Boreali-Sinica,2020,35(1):205-213.
[4] Wu Y B,Li G,Zhu Y J,Cheng Y C,Yang J Y,Chen H Z,Song X J,Ying J Z. Genome-wide identification of QTLs for grain protein content based on genotyping-by-resequencing and verification of qGPC1-1 in rice[J]. International Journal of Molecular Sciences,2020,21(2):408. doi:10.3390/ijms21020408.
[5] Zhao Y H,Ma J J,Li M,Deng L,Li G H,Xia H,Zhao S Z,Hou L,Li P C,Ma C,Yuan M,Ren L,Gu J Z,Guo B Z,Zhao C Z,Wang X J. Whole-genome resequencing-based QTL-seq identified AhTc1 gene encoding a R2R3-MYB transcription factor controlling peanut purple testa colour[J]. Plant Biotechnology Journal, 2020,18(1):96-105. doi:10.1111/pbi.13175.
[6] Liu T P,He J H,Dong K J,Wang X W,Wang W W,Yang P,Ren R Y,Zhang L,Zhang Z S,Yang T Y. QTL mapping of yield component traits on bin map generated from resequencing a RIL population of foxtail millet (Setaria italica)[J]. BMC Genomics,2020,21:141. doi:10.1186/s12864-020-6553-9.
[7] Jiang J F,Fan X C,Zhang Y,Tang X P,Li X M,Liu C H,Zhang Z W. Construction of a high-density genetic map and mapping of firmness in grapes (Vitis vinifera L.) based on whole-Genome resequencing[J]. International Journal of Molecular Sciences,2020,21(3):797. doi:10.3390/ijms21030797.
[8] Wang J Q,Li P,Mao D D,Zhang X M,Wang H,Li Y Y. The complete chloroplast genome sequence of Manglietia yuyuanensis (Magnoliaceae)[J]. Mitochondrial DNA Part B,2020,5(1):417-419. doi:10.1080/23802359.2019.1703568.
[9] Yang H J,Geng X Q,Zhao S M,Shi H Z. Genomic diversity analysis and identification of novel SSR markers in four tobacco varieties by high-throughput resequencing[J]. Plant Physiology and Biochemistry,2020,150:80-89. doi:10.1016/j.plaphy.2020.02.023.
[10] Ren X Y,Liu Y X,Tan Y M,Huan Y H,Liu Z Y,Jiang X L. Sequencing and functional annotation of the whole genome of Shiraia bambusicola[J]. Genome Report of Shiraia bambusicola,2020,10(1):23-35. doi:10.1534/g3.119.400694.
[11] Yu F,Song J,Liang J F,Wang S K,Lu J K. Whole genome sequencing and genome annotation of the wild edible mushroom, Russula griseocarnosa[J]. Genomics, 2020,112(1):603-614. doi:10.1016/j.ygeno.2019.04.012.
[12] 江泽慧. 世界竹藤[M]. 沈阳:辽宁科学技术出版社, 2002:3-10. Jiang Z H. Bamboo and rattan in the world[M]. Shenyang:Liaoning Science and Technology Press, 2002:3-10.
[13] 马乃训,赖广辉,张培新,张宏亮. 中国刚竹属[M]. 杭州:浙江科学技术出版社,2014:91-107. Ma N X,Lai G H,Zhang P X,Zhang H L. The genus Phyllostachys in China[M]. Hangzhou:Zhejiang Science and Technology Press,2014:91-107.
[14] Peng Z H,Lu Y,Li L B,Zhao Q,Feng Q,Gao Z M,Lu H Y,Hu T,Yao N,Liu K Y,Li Y,Fan D L,Guo Y L,Li W J,Lu Y Q,Weng Q J,Zhou C C,Zhang L,Huang T,Zhao Y,Zhu C R,Liu X G,Yang X W,Wang T,Miao K,Zhuang C Y,Cao X L,Tang W L,Liu G S,Liu Y L,Chen J,Liu Z J,Yuan L C,Liu Z H,Huang X H,Lu T T,Fei B H,Ning Z M,Han B,Jiang Z H. The draft genome of the fast-growing non-timber forest species moso bamboo (Phyllostachys heterocycla)[J]. Nature Genetics,2013,45(4):456-461. doi:10.1038/ng.2569.
[15] Zhao H S, Gao Z M, Wang L, Wang J L, Wang S H, Fei B H, Chen C H, Shi C C, Liu X C, Zhang H L, Lou Y F, Chen L F, Sun H Y, Zhou X Q, Wang S N, Zhang C, Xu H, Li L C, Yang Y H, Wei Y L, Yang W, Gao Q, Yang H M, Zhao S C, Jiang Z H. Chromosome-level reference genome and alternative splicing atlas of moso bamboo (Phylostachys edulis)[J]. Giga Science,2018,7(10):1-12. doi:10.1093/gigascience/giy115.
[16] 林翩翩,白有煌,周明兵. 毛竹细胞色素P450的基因组学分析[J]. 植物生理学报,2014,50(9):1387-1400. doi:10.13592/j.cnki.ppj.2014.0205. Lin P P,Bai Y H,Zhou M B. Genomics analysis of cytochrome P450 monooxygenase genes in Phyllostachys heterocycle[J]. Plant Physiology Journal, 2014,50(9):1387-1400.
[17] Wu H L,Lü H,Li L,Liu J,Mu S H,Li X P,Gao J. Genome-wide analysis of the AP2/ERF transcription factors family and the expression patterns of DREB genes in moso bamboo (Phyllostachys edulis)[J]. PLoS One,2015,10(5):e0126657. doi:10.1371/journal.pone.0126657.
[18] 黎帮勇, 胡尚连, 曹颖, 徐刚. 毛竹NAC转录因子家族生物信息学分析[J]. 基因组学与应用生物学,2015,34(8):1769-1777. doi:10.13417/j.gab.034.001769. Li B Y,Hu S L,Cao Y,Xu G. Bioinformatics analysis of NAC gene family in moso bamboo[J]. Genomics and Applied Biology,2015,34(8):1769-1777.
[19] Bai Q S,Hou D,Li L,Cheng Z C,Ge W,Liu J,Li X P,Mu S H,Gao J. Genome-wide analysis and expression characteristics of small auxin-up RNA (SAUR) genes in moso bamboo (Phyllostachys edulis)[J]. Genome,2017,60(4):325-336. doi:10.1139/gen-2016-0097.
[20] Sun H Y,Li L C,Lou Y F,Zhao H S,Gao Z M. Genome-wide identification and characterization of aquaporin gene family in moso bamboo (Phyllostachys edulis)[J]. Molecular Biology Reports,2016,43(5):437-450. doi:10.1007/s11033-016-3973-3.
[21] Wu M,Li Y,Chen D M,Liu H L,Zhu D Y,Xiang Y. Genome-wide identification and expression analysis of the IQD gene family in moso bamboo (Phyllostachys edulis)[J]. Scientific Reports,2016,6:24520. doi:10.1038/srep24520.
[22] Pan F,Wang Y,Liu H L,Wu M,Chu W Y, Chen D M,Xiang Y. Genome-wide identification and expression analysis of SBP-like transcription factor genes in moso bamboo (Phyllostachys edulis)[J]. BMC Genomics,2017,18(1):486. doi:10.1186/s12864-017-3882-4.
[23] Li H,Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform[J]. Bioinformatics,2009,25(14):1754-1760. doi:10.1093/bioinformatics/btp324.
[24] Chen D M, Chen Z,Wu M,Wang Y,Wang Y J,Yan H W,Xiang Y. Genome-wide identification and expression analysis of the HD-Zip gene family in moso bamboo (Phyllostachys edulis)[J]. Journal of Plant Growth Regulation,2017,2(36):323-337. doi:10.1007/s00344-016-9642-x.
[25] 刘俊, 黄容, 程占超,李雪平. 毛竹TCP基因家族全基因组鉴定与分析[J]. 基因组学与应用生物学,2018,37(12):5388-5397. doi:10.13417/j.gab.037.005388. Liu J,Huang R,Cheng Z C,Li X P. Genome-wide identification and the whole analysis of TCP gene family in moso bamboo (Phyllostachys edulis)[J]. Genomics and Applied Biology,2018,37(12):5388-5397.
[26] Xie L H,Li X Y,Hou D,Cheng Z C,Liu J,Li J,Mu S H,Gao J. Genome-wide analysis and expression profiling of the heat shock factor gene family in Phyllostachys edulis during development and in response to abiotic stresses[J]. Forests, 2019,10(2):100. doi:10.3390/f10020100.
[27] Liu J,Cheng Z C,Li X Y,Xie L H,Bai Y C,Peng L X,Li J,Gao J.Expression analysis and regulation network identification of the CONSTANS-Like gene family in moso bamboo (Phyllostachys edulis) under photoperiod treatments[J]. DNA and Cell Biology,2019,38(7):607-626.doi:10.1089/dna.2018.4611.
[28] Zidani S,Ferchichi A,Chaieb M. Genomic DNA extraction method from pearl millet (Pennisetum glaucum) leaves[J]. African Journal of Biotechnology,2005,4(8):862-866.doi:10.5897/AJB2005.000-3218
[29] McKenna A,Hanna M,Banks E,Sivachenko A,Cibulskis K,Kernytsky A,Garimella K,Altshuler D,Gabriel S,Daly M,DePristo M A. The genome analysis toolkit:A mapreduce framework for analyzing next-generation DNA sequencing data[J]. Genome Research,2010,20(9):1297-1303.doi:10.1101/gr.107524.110.
[30] Cingolani P,Platts A,Wang L L,Coon M,Nguyen T,Wang L,Land S J,Lu X Y,Ruden D M. A program for annotating and predicting the effects of single nucleotide polymorphisms,SnpEff:SNPs in the genome of Drosophila melanogaster strain w ¹¹¹⁸; iso -2; iso -3[J]. Fly,2012,6(2):80-92. doi:10.4161/fly.19695.
[31] Altschul S F,Madden T L,Schäffer A A, Zhang J H,Zhang Z,Miller W,Lipman D J. Gapped BLAST and PSI-BLAST:A new generation of protein database search programs[J]. Nucleic Acids Research,1997,25(17):3389-3402.doi:10.1093/nar/25.17.3389.
[32] Deng Y Y, Li J Q, Wu S F, Zhu Y P, Chen Y W, He F C, Deng Y, Deng Y, Li J, Wu S, Zhu Y, Chen Y, He F, Li J, Wu S, Zhu Y, He F, Deng Y Y, Li J Q, Wu S F, Zhu Y P, Chen Y W, He F C, Chen Y W, Wu S, Deng Li, Li J, Wu S, Zhu Y, Chen Y, He F, Wu S F, Deng M, Li J T, Zhu Y, He F Y. Integrated nr database in protein annotation system and its localization[J]. Comput Eng,2006,32(5):71-74. doi:10.1109/INFOCOM.2006.241.
[33] Tatusov R L,Galperin M Y,Natale D A,Koonin E V. The COG database:A tool for genome-scale analysis of protein functions and evolution[J]. Nucleic Acids Research,2000,28(1):33-36. doi:10.1093/nar/28.1.33.
[34] Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome[J]. Nucleic Acids Research,2004,32(S1):277-280. doi:10.1093/nar/gkh063.
[35] 魏强, 丁雨龙. 矢竹地下茎转录组测序及节间生长相关基因表达分析[J]. 南京林业大学学报(自然科学版),2017,41(5):42-48. doi:10.3969/j.issn.1000-2006.201601006. Wei Q,Ding Y L. Transcriptome sequencing,de novo assembly and expression analysis of several genes related to internode development in Pseudosasa japonica[J]. Journal of Nanjing Forestry University (Natural Science Edition),2017,41(5):42-48.
[36] 陈建华,晏存育,王艳梅,张新明,曾翔,曹学优,贺峭. 6种不同变异类型毛竹叶片的生理特性研究[J].中南林业科技大学学报, 2011,31(4):70-73. doi:10.14067/j.cnki.1673-923x.2011.04.037. Chen J H,Yan C Y,Wang Y M,Zhang X M,Zeng X,Cao X Y,He Q. Physiological characteristics of 6 different variation type bamboo's leaves[J]. Journal of Central South University of Forestry & Technology,2011,31(4):70-73.
[37] 阮晓赛. 毛竹种源及栽培变种遗传变异的AFLP和ISSR分析[D].杭州:浙江林学院,2008:30-31. Ruan X S. Genetic variations of Phyllostachys edulis provenance and cultivars by ISSR and AFLP[D]. Hangzhou:Zhejiang Forestry College,2008:30-31.

选择文件类型/文献管理软件名称

选择包含的内容