The genome of the warm-season turfgrass African bermudagrass (Cynodon transvaalensis)

Abstract

Cynodon species can be used for multiple purposes and have high economic and ecological significance. However, the genetic basis of the favorable agronomic traits of Cynodon species is poorly understood, partially due to the limited availability of genomic resources. In this study, we report a chromosome-scale genome assembly of a diploid Cynodon species, C. transvaalensis, obtained by combining Illumina and Nanopore sequencing, BioNano, and Hi-C. The assembly contains 282 scaffolds (~423.42 Mb, N50 = 5.37 Mb), which cover ~93.2% of the estimated genome of C. transvaalensis (~454.4 Mb). Furthermore, 90.48% of the scaffolds (~383.08 Mb) were anchored to nine pseudomolecules, of which the largest was 60.78 Mb in length. Evolutionary analysis along with transcriptome comparison provided a preliminary genomic basis for the adaptation of this species to tropical and/or subtropical climates, typically with dry summers. The genomic resources generated in this study will not only facilitate evolutionary studies of the Chloridoideae subfamily, in particular, the Cynodonteae tribe, but also facilitate functional genomic research and genetic breeding in Cynodon species for new leading turfgrass cultivars in the future.

References

[1] Forbes, I.Jr, & Burton, G. W. Chromosome numbers and meiosis in some Cynodon species and hybrids. Crop Sci. 3, 75-79 (1963).

[2] Harlan, J. R., & de Wet, J. M. J. Sources of variation in Cynodon dactylon (L). Pers. Crop Sci. 9, 774-778 (1969).

[3] Peterson, P., Columbus, T. & Pennington, S. Classification and biogeography of new world grasses: chloridoideae. Aliso 23, 580-594 (2007).

[4] Harlan, J. R., de Wet, J. M. J., Rawal, K. M., Felder, M. R. & Richardson, W. L. Cytogenetic studies in Cynodon L. C. Rich. (Gramineae). Crop Sci. 10, 288-291 (1970).

[5] Watson, L. & Dallwitz, M. J. The Grass Genera of the World (CABI, 1992).

[6] Beard, J. B. & Beard, H. J. Beard’s Turfgrass Encyclopedia for Golf Courses, Grounds, Lawns, Sports Fields (Michigan State University Press, 2005).

[7] Peterson, P. M., Romaschenko, K. & Johnson, G. A classification of the Chloridoideae (Poaceae) based on multi-gene phylogenetic trees. Mol. Phylogenet. Evol. 55, 580-598 (2010).

[8] Bethel, C. M. et al. A framework linkage map of bermudagrass (Cynodon dactylon × transvaalensis) based on single-dose restriction fragments. Theor. Appl. Genet. 112, 727-737 (2006).

[9] Khanal, S. et al. SSR-enriched genetic linkage maps of bermudagrass (Cynodon dactylon × transvaalensis), and their comparison with allied plant genomes. Theor. Appl. Genet. 130, 819-839 (2017).

[10] Turgeon, A. J. (ed.) Turfgrass Management (Prentice-Hall/Pearson Higher Education, 2011).

[11] Ball, D. M., Hoveland, C. S. & Lacefield, G. D. Southern Forages: Modern Concepts for Forage Crop Management (International Plant Nutrition Institute, 2015).

[12] Xu, J., Wang, Z. & Cheng, J. J. Bermuda grass as feedstock for biofuel production: a review. Bioresour. Technol. 102, 7613-7620 (2011).

[13] Garcia, A. & Mangaroo, A. S. Residual concentration of selected heavy metals in a sewage sludge-amended soil and uptake by coastal bermudagrass. in Sewage Sludge: Land Utilization and the Environment 187-192 (1994).

[14] Razmjoo, K. & Adavi, Z. Assessment of bermudagrass cultivars for phytoremediation of petroleum contaminated soils. Int. J. Phytoremediat. 14, 14-23 (2012).

[15] Rai, P. K., Dolly, J., Rai, D. K., Bechan, S. & Geeta, W. Antioxidant potential of oral feeding of Cynodon dactylon extract on diabetes-induced oxidative stress. J. Food Biochem. 34, 78-92 (2010).

[16] Al-Snafi, A. Chemical constituents and pharmacological effects of Cynodon dactylon-A Review. IOSR J. Pharm. 6, 17-31 (2016).

[17] Gaut, B. S. Evolutionary dynamics of grass genomes. N. Phytol. 154, 15-28 (2002).

[18] Cotton, J. L. et al. Resolving deep relationships of PACMAD grasses: a phylogenomic approach. BMC Plant Biol. 15, 178 (2015).

[19] II, G.P.W.G. New grass phylogeny resolves deep evolutionary relationships and discovers C4 origins. N. Phytol. 193, 304-312 (2012).

[20] Clayton, W. D., Vorontsova, M. S., Harman, K. T. & Williamson, H. GrassBase- the online world grass flora. http://www.kew.org/data/grasses-db.html (2006).

[21] Group, G. et al. Phylogeny and subfamilial classification of the grasses (Poaceae). Ann. Mo. Bot. Gard. 88, 373 (2002).

[22] Soreng, R. J. et al. A worldwide phylogenetic classification of the Poaceae (Gramineae) II: an update and a comparison of two 2015 classifications. J. Syst. Evol. 55, 259-290 (2017).

[23] Hittalmani, S. et al. Genome and transcriptome sequence of Finger millet (Eleusine coracana (L.) Gaertn.) provides insights into drought tolerance and nutraceutical properties. BMC Genom. 18, 465 (2017).

[24] Zhang, H. et al. Development of a goosegrass (Eleusine indica) draft genome and application to weed science research. Pest Manag. Sci. 75, 2776-2784 (2019).

[25] Tanaka, H. et al. Sequencing and comparative analyses of the genomes of zoysiagrasses. DNA Res. 23, 171-180 (2016).

[26] Carballo, J. et al. A high-quality genome of Eragrostis curvula grass provides insights into Poaceae evolution and supports new strategies to enhance forage quality. Sci. Rep. 9, 10250 (2019).

[27] Cannarozzi, G. et al. Genome and transcriptome sequencing identifies breeding targets in the orphan crop tef (Eragrostis tef). BMC Genom. 15, 581 (2014).

[28] Wu, Y. Q. et al. Genetic analyses of Chinese Cynodon accessions by flow cytometry and AFLP markers. Crop Sci. 46, 917-926 (2006).

[29] Gulsen, O. & Ceylan, A. Elucidating polyploidization of bermudagrasses as assessed by organelle and nuclear DNA markers. OMICS 15, 903-912 (2011).

[30] Beard, J. B. Origins of North American Turfgrasses in Turfgrass: Biology, Use, and Management (Agronomy Monographs) 1st edn, (eds Stier, J. C., Horgan, B. P. & Bonos, S. A.) (American Society of Agronomy, Soil Science Society of America, Crop Science Society of America, Inc., 2013).

[31] Clayton, W. D., Renvoize, S. A. & Phillips, S. M. Flora of Tropical East Africa: Gramineae (Crown Agents for Overseas Governments & Administrations, 1974).

[32] VanBuren, R. et al. Single-molecule sequencing of the desiccation-tolerant grass Oropetium thomaeum. Nature 527, 508-511 (2015).

[33] VanBuren, R., Wai, C. M., Keilwagen, J. & Pardo, J. A chromosome-scale assembly of the model desiccation tolerant grass Oropetium thomaeum. Plant Direct 2, e00096 (2018).

[34] Reasor, E. H. et al. Genotypic and phenotypic evaluation of off-type grasses in hybrid Bermudagrass [Cynodon dactylon (L.) Pers. x C. transvaalensis Burtt-Davy] putting greens using genotyping-by-sequencing and morphological characterization. Hereditas 155, 8 (2017).

[35] Du, H. et al. Sequencing and de novo assembly of a near complete indica rice genome. Nat. Commun. 8, 15324- 15324 (2017).

[36] 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).

[37] Ye, C. Y. et al. The genomes of the allohexaploid Echinochloa crus-galli and its progenitors provide insights into polyploidization-driven adaptation. Mol. Plant 13, 1298-1310 (2020).

[38] Werner, T. & Schmulling, T. Cytokinin action in plant development. Curr. Opin. Plant Biol. 12, 527-538 (2009).

[39] Kieber, J. J. & Schaller, G. E. Cytokinin signaling in plant development. Development 145, dev149344 (2018).

[40] Turgeon, A. J. Turfgrass Management, 9th edn (Pearson Prentice Hall, 2010).

[41] Werner, T., Motyka, V., Strnad, M. & Schmülling, T. Regulation of plant growth by cytokinin. Proc. Natl Acad. Sci. U. S. A. 98, 10487-10492 (2001).

[42] Prerostova, S. et al. Cytokinins: their impact on molecular and growth responses to drought stress and recovery in Arabidopsis. Front. Plant Sci. 9, 655 (2018).

[43] Gaut, B. S. & Doebley, J. F. DNA sequence evidence for the segmental allotetraploid origin of maize. Proc. Natl Acad. Sci. U. S. A. 94, 6809-6814 (1997).

[44] Yu, C. S., Chen, Y. C., Lu, C. H. & Hwang, J. K. Prediction of protein subcellular localization. Proteins 64, 643-651 (2006).

[45] Sarkar, N. K., Kundnani, P. & Grover, A. Functional analysis of Hsp70 superfamily proteins of rice (Oryza sativa). Cell Stress Chaperones 18, 427-437 (2013).

[46] Zhang, L. et al. Genome-wide analysis and expression profiling under heat and drought treatments of HSP70 gene family in soybean (Glycine max L.). Front. Plant Sci. 6, 773 (2015).

[47] Cannon, S. B., Mitra, A., Baumgarten, A., Young, N. D. & May, G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol. 4, 10 (2004).

[48] Koonin, E. V. Orthologs, paralogs, and evolutionary genomics. Annu. Rev. Genet. 39, 309-338 (2005).

[49] Moura, J. C., Bonine, C. A., de Oliveira Fernandes Viana, J., Dornelas, M. C. & Mazzafera, P. Abiotic and biotic stresses and changes in the lignin content and composition in plants. J. Integr. Plant Biol. 52, 360-376 (2010).

[50] Hamido, S. A., Guertal, E. A. & Wood, W. Seasonal variation of carbon and nitrogen emissions from turfgrass. Am. J. Clim. Change 5, 448-463 (2016).

[51] Karaca, M., Lang, D. J., Yerk-Davis, G. L., Sukumar, S. & Ainsworth, A. J. Determination of DNA content and genome in Cynodon species by flow cytometry. Crop Res. 20, 1-12 (2000).

[52] Soreng, R. et al. A world-wide phylogenetic classification of Poaceae (Gramineae). http://www.tropicos.org/projectwebportal.aspx?pagename=ClassificationNWG&projectid=10 (2014).

[53] Fang, T., Dong, H., Yu, S. & Moss, J. Q. Sequence-based genetic mapping of Cynodon dactylon Pers. reveals new insights into genome evolution in Poaceae. Commun. Biol. 3, 358 (2020).

[54] Guo, C. et al. The Coix genome provides insights into panicoideae evolution and papery hull domestication. Mol. Plant 13, 309-320 (2020).

[55] Hahn, M. W., De Bie, T., Stajich, J. E., Nguyen, C. & Cristianini, N. Estimating the tempo and mode of gene family evolution from comparative genomic data. Genome Res. 15, 1153-1160 (2005).

[56] Sharpton, T. J. et al. Comparative genomic analyses of the human fungal pathogens Coccidioides and their relatives. Genome Res. 19, 1722-1731 (2009).

[57] Brawley, S. H. et al. Insights into the red algae and eukaryotic evolution from the genome of Porphyra umbilicalis (Bangiophyceae, Rhodophyta). Proc. Natl Acad. Sci. U. S. A. 114, E6361-E6370 (2017).

[58] Feder, M. E. & Hofmann, G. E. Heat-shock proteins, molecular chaperones, and the stress response: evolutionary and ecological physiology. Annu. Rev. Physiol. 61, 243-282 (1999).

[59] Schulze, E.-D., Beck, E. & Muller-Hohenstein, K. Plant Ecology (Springer, 2005).

[60] Wang, W., Vinocur, B., Shoseyov, O. & Altman, A. Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response. Trends Plant Sci. 9, 244-252 (2004).

[61] Wahid, A., Gelani, S., Ashraf, M. & Foolad, M. R. Heat tolerance in plants: an overview. Environ. Exp. Bot. 61, 199-223 (2007).

[62] Kim, S.-R. & An, G. Rice chloroplast-localized heat shock protein 70, OsHsp70CP1, is essential for chloroplast development under high-temperature conditions. J. Plant Physiol. 170, 854-863 (2013).

[63] Sanmiya, K., Suzuki, K., Egawa, Y. & Shono, M. Mitochondrial small heat-shock protein enhances thermotolerance in tobacco plants. FEBS Lett. 557, 265-268 (2004).

[64] Hong, S. W. & Vierling, E. Mutants of Arabidopsis thaliana defective in the acquisition of tolerance to high temperature stress. Proc. Natl Acad. Sci. U. S. A. 97, 4392-4397 (2000).

[65] Hýsková, V. & Ryšlavá, H. Roles of HSP70 in plant abiotic stress. Molecular Approaches in Plant Abiotic Stress 1st edn, (eds Gaur, R. K. & Sharma, P.) 44-66 (CRC Press, Boca Raton, 2014).

[66] Ma, T. et al. Genomic insights into salt adaptation in a desert poplar. Nat. Commun. 4, 2797 (2013).

[67] Zhang, Z. et al. Adaptation to extreme Antarctic environments revealed by the genome of a sea ice green alga. Curr. Biol. 30, 3330-3341.e7 (2020).

[68] Leng, L. et al. A subclass of HSP70s regulate development and abiotic stress responses in Arabidopsis thaliana. J. Plant Res. 130, 349-363 (2017).

[69] Samuels, L., Kunst, L. & Jetter, R. Sealing plant surfaces: cuticular wax formation by epidermal cells. Annu. Rev. Plant Biol. 59, 683-707 (2008).

[70] Xue, D., Zhang, X., Lu, X., Chen, G. & Chen, Z.-H. Molecular and evolutionary mechanisms of cuticular wax for plant drought tolerance. Front. Plant Sci. 8, 621 (2017).

[71] Zhang, G. et al. The KCS gene is involved in the formation of chloroplast stromules and other physiological processes in jute (Corchorus capsularis L.). Ind. Crop Prod. 141, 111781 (2019).

[72] Lokesh, U. et al. Role of plant fatty acid elongase (3 keto acyl-CoA Synthase) gene in cuticular wax biosynthesis. J. Agric. Allied Sci. 2, 35-42 (2013).

[73] Zhou, Y., Lambrides, C. & Fukai, S. Drought resistance and soil water extraction of a perennial C4 grass: Contributions of root and rhizome traits. Funct. Plant Biol. 41, 505 (2014).

[74] Long, L. M., Patel, H. P., Cory, W. C. & Stapleton, A. E. The maize epicuticular wax layer provides UV protection. Funct. Plant Biol. 30, 75-81 (2003).

[75] Reicosky, D. A. & Hanover, J. W. Physiological effects of surface waxes: I. light reflectance for glaucous and nonglaucous Picea pungens. Plant Physiol. 62, 101-104 (1978).

[76] Christie, J. M. et al. Plant UVR8 photoreceptor senses UV-B by tryptophan-mediated disruption of cross-dimer salt bridges. Science 335, 1492-1496 (2012).

[77] Ning, D.-L. et al. Chromosomal-level assembly of Juglans sigillata genome using Nanopore, BioNano, and Hi-C analysis. Gigascience 9, giaa006 (2020).

[78] Zhang, J. et al. Genome of plant maca (Lepidium meyenii) Illuminates genomic basis for high-altitude adaptation in the central Andes. Mol. Plant 9, 1066-1077 (2016).

[79] Dolezel, J. & Bartos, J. Plant DNA flow cytometry and estimation of nuclear genome size. Ann. Bot. 95, 99-110 (2005).

[80] Marçais, G. & Kingsford, C. A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics 27, 764-770 (2011).

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

[82] de Lannoy, C., de Ridder, D. & Risse, J. A sequencer coming of age: de novo genome assembly using MinION reads. bioRxiv 142711 (2017).

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

[84] Hu, J., Fan, J., Sun, Z. & Liu, S. NextPolish: a fast and efficient genome polishing tool for long read assembly. Bioinformatics 36, 2253-2255 (2019).

[85] Zhang, H.-B., Zhao, X., Ding, X., Paterson, A. H. & Wing, R. A. Preparation of megabase-size DNA from plant nuclei. Plant J. 7, 175-184 (1995).

[86] Zhang, M. et al. Preparation of megabase-sized DNA from a variety of organisms using the nuclei method for advanced genomics research. Nat. Protoc. 7, 467-478 (2012).

[87] Weissensteiner, M. H. et al. Discovery and population genomics of structural variation in a songbird genus. Nat. Commun. 11, 3403 (2020).

[88] Belton, J. M. et al. Hi-C: a comprehensive technique to capture the conformation of genomes. Methods 58, 268-276 (2012).

[89] Shi, J. et al. Chromosome conformation capture resolved near complete genome assembly of broomcorn millet. Nat. Commun. 10, 464 (2019).

[90] Burton, J. N. et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat. Biotechnol. 31, 1119-1125 (2013).

[91] Chen, S., Zhou, Y., Chen, Y. & Gu, J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884-i890 (2018).

[92] Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357-359 (2012).

[93] Korbel, J. O. & Lee, C. Genome assembly and haplotyping with Hi-C. Nat. Biotechnol. 31, 1099-1101 (2013).

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

[95] 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).

[96] Bao, W., Kojima, K. K. & Kohany, O. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mob. DNA 6, 11 (2015).

[97] Haas, B. J. et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments. Genome Biol. 9, R7 (2008).

[98] Stanke, M., Schöffmann, O., Morgenstern, B. & Waack, S. Gene prediction in eukaryotes with a generalized hidden Markov model that uses hints from external sources. BMC Bioinform. 7, 62 (2006).

[99] Keilwagen, J. et al. Using intron position conservation for homology-based gene prediction. Nucleic Acids Res. 44, e89 (2016).

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

[101] Krzywinski, M. et al. Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639-1645 (2009).

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

[103] Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772-780 (2013).

[104] Castresana, J . Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol. Biol. Evol. 17, 540-552 (2000).

[105] Stamatakis, A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22, 2688-2690 (2006).

[106] Yang, Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586-1591 (2007).

[107] 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).

[108] Zhang, Z. et al. KaKs_Calculator: calculating Ka and Ks through model selection and model averaging. Genom, Proteom. Bioinform. 4, 259-263 (2006).

[109] Kim, D., Langmead, B. & Salzberg, S. L. HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357-360 (2015).

[110] Pertea, M. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290-295 (2015).

[111] Yu, G., Wang, L. G., Han, Y. & He, Q. Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284-287 (2012).