Genome sequence and evolution of Betula platyphylla

Abstract

Betula L. (birch) is a pioneer hardwood tree species with ecological, economic, and evolutionary importance in the Northern Hemisphere. We sequenced the Betula platyphylla genome and assembled the sequences into 14 chromosomes. The Betula genome lacks evidence of recent whole-genome duplication and has the same paleoploidy level as Vitis vinifera and Prunus mume. Phylogenetic analysis of lignin pathway genes coupled with tissue-specific expression patterns provided clues for understanding the formation of higher ratios of syringyl to guaiacyl lignin observed in Betula species. Our transcriptome analysis of leaf tissues under a time-series cold stress experiment revealed the presence of the MEKK1–MKK2–MPK4 cascade and six additional mitogen-activated protein kinases that can be linked to a gene regulatory network involving many transcription factors and cold tolerance genes. Our genomic and transcriptome analyses provide insight into the structures, features, and evolution of the B. platyphylla genome. The chromosome-level genome and gene resources of B. platyphylla obtained in this study will facilitate the identification of important and essential genes governing important traits of trees and genetic improvement of B. platyphylla.

References

[1] Singewar, K., Moschner, C. R., Hartung, E. & Fladung, M. Species determination and phylogenetic relationships of the genus Betula inferred from multiple chloroplast and nuclear regions reveal the high methyl salicylate-producing ability of the ancestor. Trees-Struct. Funct. https://doi.org/10.1007/s00468-020-01984-x (2020).

[2] Cheng, Z.-M., Boe, A. A., Schnurr, J. P. & Herman, D. E. Betula platyphylla ‘Fargo’: an Asian White Birch for the Northern Plains. HortScience 32, 1304-1305 (1997).

[3] Osumi, K. Reciprocal distribution of two congeneric trees, Betula platyphylla var. japonica and Betula maximowicziana, in a landscape dominated by anthropogenic disturbances in northeastern Japan. J. Biogeogr. 32, 2057-2068 (2005).

[4] Rastogi, S., Pandey, M. M. & Rawat, A. K. S. Medicinal plants of the genus Betula- Traditional uses and a phytochemical-pharmacological review. J. Ethnopharmacol. 159, 62-83 (2015).

[5] Matsuki, S. & Koike, T. Comparison of leaf life span, photosynthesis and defensive traits across seven species of deciduous broad-leaf tree seedlings. Ann. Bot. 97, 813-817 (2006).

[6] Koski, V. & Rousi, M. A review of the promises and constraints of breeding silver birch (Betula pendula Roth) in Finland. Forestry 78, 187-198 (2005).

[7] Chin, C.-S. et al. Phased diploid genome assembly with single-molecule real-time sequencing. Nat. Methods 13, 1050 (2016).

[8] Ruan, J. & Li, H. Fast and accurate long-read assembly with wtdbg2. Nat. Meth. 17, 155-158 (2020).

[9] Chakraborty, M., Baldwin-Brown, J. G., Long, A. D. & Emerson, J. Contiguous and accurate de novo assembly of metazoan genomes with modest long read coverage. Nucleic Acids Res. 44, e147- e147 (2016).

[10] Lam, K.-K., LaButti, K., Khalak, A. & Tse, D. FinisherSC: a repeat-aware tool for upgrading de novo assembly using long reads. Bioinformatics 31, 3207-3209 (2015).

[11] Vaser, R., Sović, I., Nagarajan, N. & Šikić, M. Fast and accurate de novo genome assembly from long uncorrected reads. Genome Res. 27, 737-746 (2017).

[12] Walker, B. J. et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PloS ONE 9, e112963 (2014).

[13] Camacho, C. et al. BLAST+: architecture and applications. BMC Bioinformatics 10, 421 (2009).

[14] Li, H. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094-3100 (2018).

[15] Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210-3212 (2015).

[16] Sun, X. et al. SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS ONE 8, e58700 (2013).

[17] Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Preprint at https://arxiv.org/abs/1303.3997 (2013).

[18] Tang, H. et al. ALLMAPS: robust scaffold ordering based on multiple maps. Genome Biol. 16, 3 (2015).

[19] Tuskan, G. A. et al. The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science 313, 1596-1604 (2006).

[20] Jaillon, O. et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature 449, 463-467 (2007).

[21] Zhang, Q. et al. The genome of Prunus mume. Nat. Commun. 3, 1318 (2012).

[22] Initiative, A. G. Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 408, 796 (2000).

[23] Schmidt, S. M. et al. MITEs in the promoters of effector genes allow prediction of novel virulence genes in Fusarium oxysporum. BMC Genomics 14, 119 (2013).

[24] Wei, H., Albert, H. & Moore, P. Differential expression of sugarcane polyubiquitin genes and isolation of promoters from two highly-expressed members of the gene family. J. Plant Physiol. 155, 513-519 (1999).

[25] Chen, J., Hu, Q., Zhang, Y., Lu, C. & Kuang, H. P-MITE: a database for plant miniature inverted-repeat transposable elements. Nucleic Acids Res. 42, D1176-D1181 (2014).

[26] Liu, Y. & Yang, G. Tc1-like transposable elements in plant genomes. Mob. DNA 5, 17 (2014).

[27] Albert, V. A. et al. The Amborella genome and the evolution of flowering plants. Science 342, 1241089 (2013).

[28] Tang, H. et al. Synteny and collinearity in plant genomes. Science 320, 486-488 (2008).

[29] Tang, H. et al. Screening synteny blocks in pairwise genome comparisons through integer programming. BMC Bioinformatics 12, 102 (2011).

[30] Bowers, J. E., Chapman, B. A., Rong, J. K. & Paterson, A. H. Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events. Nature 422, 433-438 (2003).

[31] Salojärvi, J. et al. Genome sequencing and population genomic analyses provide insights into the adaptive landscape of silver birch. Nat. Genet. 49, 904-912 (2017).

[32] Santos, R. B., Capanema, E. A., Balakshin, M. Y., Chang, H.-M. & Jameel, H. Effect of hardwoods characteristics on kraft pulping process: emphasis on lignin structure. BioResources 6, 3623-3637 (2011).

[33] Lu, S. et al. Ptr-miR397a is a negative regulator of laccase genes affecting lignin content in Populus trichocarpa. Proc. Natl Acad. Sci. USA 110, 10848-10853 (2013).

[34] Skyba, O., Douglas, C. J. & Mansfield, S. D. Syringyl-rich lignin renders poplars more resistant to degradation by wood decay fungi. Appl. Environ. Microbiol. 79, 2560-2571 (2013).

[35] Möller, P. et al. Living at the margin of the retreating Fennoscandian Ice Sheet: the early Mesolithic sites at Aareavaara, northernmost Sweden. Holocene 23, 104-116 (2013).

[36] Qing, K. et al. The MEKK1-MKK1/MKK2-MPK4 kinase cascade negatively regulates immunity mediated by a mitogen-activated protein kinase kinase kinase in Arabidopsis. Plant Cell 24, 2225-2236 (2012).

[37] Tomoyuki, F., Daisuke, M. & Takashi, N. Membrane rigidification functions upstream of the MEKK1-MKK2-MPK4 cascade during cold acclimation in Arabidopsis thaliana. FEBS Lett. 588, 2025-2030 (2014).

[38] Zhao, M., Liu, W., Xia, X., Wang, T. & Zhang, W. H. Cold acclimation-induced freezing tolerance of Medicago truncatula seedlings is negatively regulated by ethylene. Physiol. Plant. 152, 115-129 (2014).

[39] Stam, P. Construction of integrated genetic linkage maps by means of a new computer package: Join Map. Plant J. 3, 739-744 (1993).

[40] Edgar, R. C. PILER-CR: fast and accurate identification of CRISPR repeats. BMC Bioinformatics 8, 18 (2007).

[41] Price, A. L., Jones, N. C. & Pevzner, P. A. De novo identification of repeat families in large genomes. Bioinformatics 21, i351-i358 (2005).

[42] Xu, Z. & Wang, H. LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 35, W265-W268 (2007).

[43] Han, Y. & Wessler, S. R. MITE-Hunter: a program for discovering miniature inverted-repeat transposable elements from genomic sequences. Nucleic Acids Res. 38, e199 (2010).

[44] Yang, G. MITE Digger, an efficient and accurate algorithm for genome wide discovery of miniature inverted repeat transposable elements. BMC Bioinformatics 14, 186 (2013).

[45] Stanke, M., Diekhans, M., Baertsch, R. & Haussler, D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics 24, 637-644 (2008).

[46] Burge, C. B. & Karlin, S. Finding the genes in genomic DNA. Curr. Opin. Struct. Biol. 8, 346-354 (1998).

[47] Schattner, P., Brooks, A. N. & Lowe, T. M. The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 33, W686-W689 (2005).

[48] Nawrocki, E. P. & Eddy, S. R. Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 29, 2933-2935 (2013).

[49] Jones, P. et al. InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236-1240 (2014).

[50] Li, L., Stoeckert, C. J. & Roos, D. S. OrthoMCL: Identification of ortholog groups for eukaryotic genomes. Genome Res. 13, 2178-2189 (2003).

[51] De Bie, T., Cristianini, N., Demuth, J. P. & Hahn, M. W. CAFE: a computational tool for the study of gene family evolution. Bioinformatics 22, 1269-1271 (2006).

[52] Shi, R. et al. Towards a systems approach for lignin biosynthesis in Populus trichocarpa: transcript abundance and specificity of the monolignol biosynthetic genes. Plant Cell Physiol. 51, 144-163 (2010).

[53] Raes, J., Rohde, A., Christensen, J. H., Van de Peer, Y. & Boerjan, W. Genome-wide characterization of the lignification toolbox in Arabidopsis. Plant Physiol. 133, 1051-1071 (2003).

[54] Zhang, W., Wei, R. & Chen, S. Functional characterization of cinnamoyl-coa reductase (CCR) in birch (Betula platyphylla× Betula pendula) through overexpression and suppression analysis. Physiol. Plant 154, 283-296 (2015).

[55] Yang, C. & Wei, H. Designing microarray and RNA-seq experiments for greater systems biology discovery in modern plant genomics. Mol. Plant 8, 196-206 (2015).

[56] Zhang, W. et al. Functional characterization of CCR in birch (Betula platyphylla x Betula pendula) through overexpression and suppression analysis. Physiol. Plant 154, 283-296 (2015).

[57] Nikolayeva, O. & Robinson, M. D. In Stem Cell Transcriptional Networks 45-79 (Springer, 2014).

[58] Margolin, A. A. et al. ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinformatics 7, https://doi.org/10.1186/1471-2105-7-S1-S7 (2006).

[59] Wei, H., Yordanov, Y., Kumari, S., Georgieva, T. & Busov, V. Genetic networks involved in poplar root response to low nitrogen. Plant Signal Behav. 8, e27211 (2013).