Genomic variation reveals demographic history and biological adaptation of the ancient relictual, lotus (Nelumbo Adans.)

Abstract

Lotus (Nelumbo Adans.), a relict plant, is testimony to long-term sustained ecological success, but the underlying genetic changes related to its survival strategy remain unclear. Here, we assembled the high-quality lotus genome, investigated genome variation of lotus mutation accumulation (MA) lines and reconstructed the demographic history of wild Asian lotus. We identified and validated 43 base substitutions fixed in MA lines, implying a spontaneous mutation rate of 1.4 × 10⁻⁹ bases/generation in lotus shoot stem cells. The past history of the lotus revealed that the ancestors of the lotus in eastern and southern Asia could be traced back to ∼20 million years ago and twice experienced significant bottlenecks and population splits. We further identified selected genes among three lotus groups in different habitats, suggesting that 453 differed genes between the tropical and temperate group and 410 differed genes between two subgroups from Northeastern China and the Yangtze River–Yellow River Basin might play important roles in natural selection in the lotus’s adaptation and resilience. Our findings not only improve understanding of the evolutionary history of the lotus and the genetic basis of its survival advantages, but also provide valuable data for addressing various questions in evolution and protection for relict plants.

References

[1] Werth AJ, Shear WA . The evolutionary truth about living fossils. Am Sci. 2014; 102: 434-40.

[2] Scott L, Love AC . Rethinking living fossils. Bioscience. 2018; 68: 760-70.

[3] Simpson MG . Diversity and classification of flowering plants: eudicots. In: Plant Systematics . San Diego: Academic Press, 2010, 275-448.

[4] Sciences, W. B. G. C. A. o. In: Chinese Lotus . Beijing: Science Press, 1987.

[5] Chen D, Ji Q, Ke W . Study on origin, geographical distribution and classification of lotus. J Changjiang Vegetables. 2017; 18: 36-40.

[6] Liu L, Li Y, Min J et al. Analysis of the cultivar names and characteristics of global lotus (Nelumbo) . Hans J Agric Sci. 2019; 9: 163-181.

[7] Connecticut Invasive Plants Council . Connecticut Invasive Plant List . 2018. Retrieved from https://cipwg.uconn.edu/wp-content/uploads/sites/244/2018/12/CT-Invasive-Plant-List-2018Scientific-Name.pdf.

[8] Ming R, VanBuren R, Liu Y et al. Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn.) . Genome Biol. 2013; 14: R41.

[9] Wang Y, Fan G, Liu Y et al. The sacred lotus genome provides insights into the evolution of flowering plants. Plant J. 2013; 76: 557-67.

[10] Gui S, Peng J, Wang X et al. Improving Nelumbo nucifera genome assemblies using high-resolution genetic maps and BioNano genome mapping reveals ancient chromosome rearrangements. Plant J. 2018; 94: 721-34.

[11] Zhao M, Yang J-X, Mao T-Y et al. Detection of highly differentiated genomic regions between lotus (Nelumbo nucifera Gaertn.) with contrasting plant architecture and their functional relevance to plant architecture . Front Plant Sci. 2018; 9: 1219-33.

[12] Fu Y, Liu F, Li S et al. Genetic diversity of the wild Asian lotus (Nelumbo nucifera) from northern China . Hortic Plant J. 2021; 7: 488-500.

[13] Li H, Yang X, Wang Q et al. Distinct methylome patterns contribute to ecotypic differentiation in the growth of the storage organ of a flowering plant (sacred lotus). Mol Ecol. 2021; 30: 2831-45.

[14] Shi T, Rahmani RS, Gugger PF et al. Distinct expression and methylation patterns for genes with different fates following a single whole-genome duplication in flowering plants. Mol Biol Evol. 2020; 37: 2394-413.

[15] Ossowski S, Schneeberger K, Lucas-Lledó JI et al. The rate and molecular spectrum of spontaneous mutations in Arabidopsis thaliana . Science. 2010; 327: 92-6.

[16] PD, RW, Halligan DL et al. Estimation of the spontaneous mutation rate per nucleotide site in a Drosophila melanogaster full-sib family. Genetics. 2014; 196: 313-33.

[17] Sablowski R . Plant and animal stem cells: conceptually similar, molecularly distinct? Trends Cell Biol. 2004; 14: 605-16.

[18] Burian A, de Reuille PB, Kuhlemeier C . Patterns of stem cell divisions contribute to plant longevity. Curr Biol. 2016; 26: 1385-94.

[19] Yang N, Xu X-W, Wang R-R et al. Contributions of Zea mays subspecies mexicana haplotypes to modern maize . Nat Commun. 2017; 8: 1874-84.

[20] Xu S, Stapley J, Gablenz S et al. Low genetic variation is associated with low mutation rate in the giant duckweed. Nat Commun. 2019; 10: 1243-9.

[21] Lynch M . Evolution of the mutation rate. Trends Genet. 2010; 26: 345-52.

[22] Plomion C, Aury J-M, Amselem J et al. Oak genome reveals facets of long lifespan. Nat Plants. 2018; 4: 440-52.

[23] Yang S, Wang L, Huang J et al. Parent-progeny sequencing indicates higher mutation rates in heterozygotes. Nature. 2015; 523: 463-7.

[24] Dobzhansky T . Genetics and the origin of species. Nature. 1959; 184: 587-8.

[25] Liu Z, Zhu H, Zhou J et al. Resequencing of 296 cultivated and wild lotus accessions unravels its evolution and breeding history. Plant J. 2020; 104: 1673-84.

[26] Zhao Y-P, Fan G, Yin P-P et al. Resequencing 545 ginkgo genomes across the world reveals the evolutionary history of the living fossil. Nat Commun. 2019; 10: 4201-11.

[27] Bai WN, Yan P-C, Zhang B-W et al. Demographically idiosyncratic responses to climate change and rapid Pleistocene diversification of the walnut genus Juglans (Juglandaceae) revealed by whole-genome sequences . New Phytol. 2018; 217: 1726-36.

[28] Zhou X, Meng X, Liu Z et al. Population genomics reveals low genetic diversity and adaptation to hypoxia in snub-nosed monkeys. Mol Biol Evol. 2016; 33: 2670-81.

[29] Yang Y, Ma T, Wang Z et al. Genomic effects of population collapse in a critically endangered ironwood tree Ostrya rehderiana . Nat Commun. 2018; 9: 5449-58.

[30] Li Y, Jones L, McQueen-Mason S . Expansins and cell growth. Curr Opin Plant Biol. 2003; 6: 603-13.

[31] Zhang Y, Zhang B, Yan D et al. Two Arabidopsis cytochrome P450 monooxygenases, CYP714A1 and CYP714A2, function redundantly in plant development through gibberellin deactivation. Plant J. 2011; 67: 342-53.

[32] Cheng CY, Krishnakumar V, Chan AP et al. Araport11: a complete reannotation of the Arabidopsis thaliana reference genome. Plant J. 2017; 89: 789-804.

[33] Goring DR, Rothstein SJ . The S-locus receptor kinase gene in a self-incompatible Brassica napus line encodes a functional serine/threonine kinase. Plant Cell. 1992; 4: 1273-81.

[34] Barrett SCH . Influences of clonality on plant sexual reproduction. Proc Natl Acad Sci USA. 2015; 112: 8859-66.

[35] Holsinger KE . Reproductive systems and evolution in vascular plants. Proc Natl Acad Sci USA. 2000; 97: 7037-42.

[36] Bricker TM, Roose JL, Zhang P et al. The PsbP family of proteins. Photosynth Res. 2013; 116: 235-50.

[37] Ping L, Liu H, Yang H et al. Translocation of drought-responsive proteins from the chloroplasts. Cell. 2020; 9: 259-76.

[38] Sun F, Yuan J, Guan Y . Asymmetric change of maximum and minimum temperature in the Northeast China from 1959-2002. Sci Geogr Sin. 2008; 28: 532-6.

[39] Francisco L, María L-CR, Rafael C et al. A novel cold-inducible gene from Arabidopsis, RCI3, encodes a peroxidase that constitutes a component for stress tolerance. Plant J. 2002; 32: 13-24.

[40] Otto H, Reche PA, Bazan F et al. In silico characterization of the family of PARP-like poly(ADP-ribosyl)transferases (pARTs). BMC Genomics. 2005; 6: 139-62.

[41] Tchou J, Kasai H, Shibutani S et al. 8-Oxoguanine (8-hydroxyguanine) DNA glycosylase and its substrate specificity. Proc Natl Acad Sci USA. 1991; 88: 4690-4.

[42] PD, REJ, Barton MK . Interactions between the cell cycle and embryonic patterning in Arabidopsis uncovered by a mutation in DNA polymerase epsilon. Plant Cell. 2005; 17: 3362-77.

[43] Feng X, Koh DW . Roles of poly(ADP-ribose) glycohydrolase in DNA damage and apoptosis. Int Rev Cell Mol Biol. 2013; 304: 227-81.

[44] Ding L, Cao J, Duan Y et al. Proteomic and physiological responses of Arabidopsis thaliana exposed to salinity stress and N-acyl-homoserine lactone. Physiol Plant. 2016; 158: 414-34.

[45] Alfenito MR, Souer E, Goodman CD et al. Functional complementation of anthocyanin sequestration in the vacuole by widely divergent glutathione S-transferases. Plant Cell. 1998; 10: 1135-49.

[46] Zhang Y, Xu YH, Yi HY et al. Vacuolar membrane transporters OsVIT1 and OsVIT2 modulate iron translocation between flag leaves and seeds in rice. Plant J. 2012; 72: 400-10.

[47] Takehisa H, Sato Y, Antonio BA et al. Global transcriptome profile of rice root in response to essential macronutrient deficiency. Plant Signal Behav. 2013; 8: e24409.

[48] Dong C, Zheng X, Li G et al. Molecular cloning and expression of two cytosolic copper-zinc superoxide dismutases genes from Nelumbo nucifera . Appl Biochem Biotechnol. 2011; 163: 679-91.

[49] Xue J, Xu G . Pu Lan Dian antique lotus. Life World. 2015; 6: 20-5.

[50] Mather JR, Wright HE, Kutzbach JEJ et al. Global climates since the last glacial maximum. Rev Palaeobot Palynol. 1993; 91: 425-6.

[51] Yang D. On basic features of climate and environment in quaternary glacial period in eastern China. Mar Geol Quat Geol 1990; 1: 71-79.

[52] Nadachowska-Brzyska K, Li C, Smeds L et al. Temporal dynamics of avian populations during Pleistocene revealed by whole-genome sequences. Curr Biol. 2015; 25: 1375-80.

[53] Allen GC, Flores-Vergara MA, Krasynanski S et al. A modified protocol for rapid DNA isolation from plant tissues using cetyltrimethylammonium bromide. Nat Protoc. 2006; 1: 2320-5.

[54] Chen S, Zhou Y, Chen Y et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018; 34: 884-90.

[55] Michael TP, Jupe F, Bemm F et al. High contiguity Arabidopsis thaliana genome assembly with a single nanopore flow cell. Nat Commun. 2018; 9: 541-9.

[56] Shi J, Ma X, Zhang J et al. Chromosome conformation capture resolved near complete genome assembly of broomcorn millet. Nat Commun. 2019; 10: 464-73.

[57] Ruan J, Li H . Fast and accurate long-read assembly with wtdbg2. Nat Methods. 2020; 17: 155-8.

[58] Walker BJ, Abeel T, Shea T et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014; 9: e112963.

[59] Simão FA, Waterhouse RM, Ioannidis P et al. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015; 31: 19-23.

[60] Ou S, Chen J, Jiang N . Assessing genome assembly quality using the LTR assembly index (LAI). Nucleic Acids Res. 2018; 46: e126.

[61] Langmead B, Salzberg SL . Fast gapped-read alignment with bowtie 2. Nat Methods. 2012; 9: 357-9.

[62] Servant N, Varoquaux N, Lajoie BR et al. HiC-pro: an optimized and flexible pipeline for Hi-C data processing. Genome Biol. 2015; 16: 259.

[63] Robinson JT, Turner D, Durand NC et al. Juicebox.Js provides a cloud-based visualization system for Hi-C data. Cell Syst. 2018; 6: 256-8.

[64] Burton JN, Adey A, Patwardhan RP et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat Biotechnol. 2013; 31: 1119-25.

[65] Xuewen W, Le W . GMATA: an integrated software package for genome-scale SSR mining, marker development and viewing. Front Plant Sci. 2016; 7: 1350-61.

[66] Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999; 27: 573-80.

[67] Han Y, Wessler SR . MITE-hunter: a program for discovering miniature inverted-repeat transposable elements from genomic sequences. Nucleic Acids Res. 2010; 38: e199.

[68] Zhao X, Hao W . LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons. Nucleic Acids Res. 2007; 35: 265-8.

[69] Ellinghaus D, Kurtz S, Willhoeft U . LTRharvest, an efficient and flexible software for de novo detection of LTR retrotransposons . BMC Bioinformatics. 2008; 9: 18.

[70] Shujun O, Ning J . LTR_retriever: a highly accurate and sensitive program for identification of long terminal repeat retrotransposons. Plant Physiol. 2018; 176: 1410-22.

[71] Abrusán G, Grundmann N, DeMester L et al. TEclass-a tool for automated classification of unknown eukaryotic transposable elements. Bioinformatics. 2009; 25: 1329-30.

[72] Stanke M, Waack S . Gene prediction with a hidden Markov model and a new intron submodel. Bioinformatics. 2003; 19: 215-25.

[73] Keilwagen J, Wenk M, Erickson JL et al. Using intron position conservation for homology-based gene prediction. Nucleic Acids Res. 2016; 44: e89.

[74] Haas BJ, Delcher AL, Mount SM et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res 2003; 31: 5654-66.

[75] Haas BJ, Salzberg SL, Zhu W et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the Program to Assemble Spliced Alignments. Genome Biol. 2008; 9: 7-15.

[76] Chen K, Li K, Zhou Y et al. SNP genotyping based on next generation sequencing data. J Donghua Univ (Nat Sci). 2017; 43: 370-6.

[77] Li H, Handsaker B, Wysoker A et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009; 25: 2078-9.

[78] Yang J, Lee SH, Goddard ME et al. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011; 88: 76-82.

[79] Li H, Durbin R . Inference of human population history from individual whole-genome sequences. Nature. 2011; 475: 493-6.

[80] Pickrell JK, Pritchard JK . Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 2012; 8: e1002967.

[81] Excofffier L, Marchi N, Marques DA et al. Fastsimcoal2: demographic inference under complex evolutionary scenarios. Bioinformatics. 2021; 468: 4882-5.

[82] Yu G, Wang LG, Han Y et al. ClusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16: 284-7.