A genome for Cissus illustrates features underlying its evolutionary success in dry savannas

Abstract

Cissus is the largest genus in Vitaceae and is mainly distributed in the tropics and subtropics. Crassulacean acid metabolism (CAM), a photosynthetic adaptation to the occurrence of succulent leaves or stems, indicates that convergent evolution occurred in response to drought stress during species radiation. Here we provide the chromosomal level assembly of Cissus rotundifolia (an endemic species in Eastern Africa) and a genome-wide comparison with grape to understand genome divergence within an ancient eudicot family. Extensive transcriptome data were produced to illustrate the genetics underpinning C. rotundifolia’s ecological adaption to seasonal aridity. The modern karyotype and smaller genome of C. rotundifolia (n = 12, 350.69 Mb/1C), which lack further whole-genome duplication, were mainly derived from gross chromosomal rearrangements such as fusions and segmental duplications, and were sculpted by a very recent burst of retrotransposon activity. Bias in local gene amplification contributed to its remarkable functional divergence from grape, and the specific proliferated genes associated with abiotic and biotic responses (e.g. HSP-20, NBS-LRR) enabled C. rotundifolia to survive in a hostile environment. Reorganization of existing enzymes of CAM characterized as diurnal expression patterns of relevant genes further confer the ability to thrive in dry savannas.

References

[1] Wen J, Lu LM, Nie ZL et al. A new phylogenetic tribal classification of the grape family (Vitaceae). J Syst Evol. 2018; 56: 262-72.

[2] Kubitzki K, Bayer C, Stevens PF . Flowering Plants. Eudicots: Berberidopsidales, Buxales, Crossosomatales, Fabales p.p., Geraniales, Gunnerales, Myrtales P.P., Proteales, Saxifragales, Vitales, Zygophyllales, Clusiaceae Alliance, Passifloraceae Alliance, Dilleniaceae, Huaceae, Picramniaceae, Sabiaceae. Berlin: Springer; 2007.

[3] Zhang N, Wen J, Zimmer EA . Another look at the phylogenetic position of the grape order Vitales: chloroplast phylogenomics with an expanded sampling of key lineages. Mol Phylogenet Evol. 2016; 101: 216-23.

[4] Zeng LP, Zhang N, Zhang QA et al. Resolution of deep eudicot phylogeny and their temporal diversification using nuclear genes from transcriptomic and genomic datasets. New Phytol. 2017; 214: 1338-54.

[5] Jaillon O, Aury JM, Noel B et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007; 449: 463-7.

[6] DeSanto AV, Bartoli G . Crassulacean acid metabolism in leaves and stems of Cissus quadrangularis. In: Winter K, JAC S, eds. Crassulacean Acid Metabolism.Springer: Berlin, 1996, 216-29.

[7] Chu ZF, Wen J, Yang YP et al. Genome size variation and evolution in the grape family Vitaceae. J Syst Evol. 2018; 56: 273-82.

[8] Griffiths H, Males J . Succulent plants. Curr Biol. 2017; 27: R890-6.

[9] Niechayev NA, Paula NP, John CC . Understanding trait diversity associated with crassulacean acid metabolism (CAM). Curr Opin Plant Biol. 2019; 49: 74-85.

[10] Cai J, Liu X, Vanneste K et al. The genome sequence of the orchid Phalaenopsis equestris. Nat Genet. 2015; 47: 65-72.

[11] Ming R, VanBuren R, Wai CM et al. The pineapple genome and the evolution of CAM photosynthesis. Nat Genet. 2015; 47: 1435-42.

[12] Wai CM, Weise SE, Ozersky P et al. Time of day and network reprogramming during drought induced CAM photosynthesis in Sedum album. PLoS Genet. 2019; 15: e1008209.

[13] Yan L, Wang X, Liu H et al. The genome of Dendrobium officinale illuminates the biology of the important traditional Chinese orchid herb. Mol Plant. 2015; 8: 922-34.

[14] Yang XH, Hu RB, Yin HF et al. The Kalanchoë genome provides insights into convergent evolution and building blocks of crassulacean acid metabolism. Nat Commun. 2017; 8: 1899.

[15] Olivares E, Urich R, Montes G et al. Occurrence of crassulacean acid metabolism in Cissus trifoliata L. (Vitaceae). Oecologia. 1984; 61: 358-62.

[16] Sayed OH . Crassulacean acid metabolism 1975-2000, a check list. Photosynthetica. 2001; 39: 339-52.

[17] Ting IP, Sternberg LO, Deniro MJ . Variable photosynthetic metabolism in leaves and stems of Cissus quadrangularis L. Plant Physiol. 1983; 71: 677-9.

[18] Al-Bukhaiti WQ, Noman A, Mahdi AA et al. Proximate composition, nutritional evaluation and functional properties of a promising food: Arabian wax Cissus (Cissus rotundifolia Forssk) leaves. J Food Sci Technol. 2019; 56: 4844-54.

[19] Marcais G, Kingsford C . A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics. 2011; 27: 764-70.

[20] Ranallo-Benavidez TR, Jaron KS, Schatz MC . GenomeScope 2.0 and Smudgeplot for reference-free profiling of polyploid genomes. Nat Commun. 2020; 11: 1432.

[21] Galbraith DW, Harkins KR, Maddox JM et al. Rapid flow cytometric analysis of the cell cycle in intact plant tissues. Science. 1983; 220: 1049-51.

[22] Loureiro J, Rodriguez E, Dolezel J et al. Two new nuclear isolation buffers for plant DNA flow cytometry: a test with 37 species. Ann Bot. 2007; 100: 875-88.

[23] Kirov I, Divashuk M, Van Laere K et al. An easy "SteamDrop" method for high quality plant chromosome preparation. Mol Cytogenet. 2014; 7: 21.

[24] Gichuki DK, Ma L, Zhu Z et al. Genome size, chromosome number determination, and analysis of the repetitive elements in Cissus quadrangularis. PeerJ. 2019; 7: e8201.

[25] Doyle JJ, Doyle JL . A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987; 19: 11-5.

[26] Xie T, Zheng JF, Liu S et al. De novo plant genome assembly based on chromatin interactions: a case study of Arabidopsis thaliana. Mol Plant. 2015; 8: 489-92.

[27] Durand NC, Shamim MS, Machol I et al. Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst. 2016; 3: 95-8.

[28] Dudchenko O, Batra SS, Omer AD et al. De novo assembly of the Aedes aegypti genome using Hi-C yields chromosome-length scaffolds. Science. 2017; 356: 92-5.

[29] Durand NC, Robinson JT, Shamim MS et al. Juicebox provides a visualization system for Hi-C contact maps with unlimited zoom. Cell Syst. 2016; 3: 99-101.

[30] Koren S, Walenz BP, Berlin K et al. Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 2017; 27: 722-36.

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

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

[33] Pryszcz LP, Gabaldón T . Redundans: an assembly pipeline for highly heterozygous genomes. Nucleic Acids Res. 2016; 44: e113-3.

[34] Roach MJ, Schmidt SA, Borneman AR . Purge Haplotigs: allelic contig reassignment for third-gen diploid genome assemblies. BMC Bioinformatics. 2018; 19: 460.

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

[36] Jurka J, Kapitonov VV, Pavlicek A et al. Repbase update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005; 110: 462-7.

[37] Hubley R, Finn RD, Clements J et al. The Dfam database of repetitive DNA families. Nucleic Acids Res. 2016; 44: D81-9.

[38] Thiel T, Michalek W, Varshney RK et al. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003; 106: 411-22.

[39] Ou S, Jiang N . LTR_retriever: a highly accurate and sensitive program for identification of long terminal repeat retrotransposons. Plant Physiol. 2018; 176: 1410-22.

[40] Haas BJ, Papanicolaou A, Yassour M et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013; 8: 1494-512.

[41] 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.

[42] Stanke M, Keller O, Gunduz I et al. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 2006; 34: W435-9.

[43] Korf I. Gene finding in novel genomes. BMC Bioinformatics. 2004; 5: 59.

[44] Lukashin AV, Borodovsky M . GeneMark.Hmm: new solutions for gene finding. Nucleic Acids Res. 1998; 26: 1107-15.

[45] 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: R7.

[46] Camacho C, Coulouris G, Avagyan V et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009; 10: 21.

[47] Li L, Stoeckert CJ Jr, Roos DS . OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003; 13: 2178-89.

[48] Edgar RC . MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004; 32: 1792-7.

[49] Stamatakis A . RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006; 22: 2688-90.

[50] Yang Z . PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007; 24: 1586-91.

[51] Han MV, Thomas GW, Lugo-Martinez J et al. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol Biol Evol. 2013; 30: 1987-97.

[52] Albert VA, Barbazuk WB, dePamphilis CW et al. The Amborella genome and the evolution of flowering plants. Science. 2013; 342: 1241089.

[53] Wickham H . ggplot2. Wiley interdisciplinary reviews: computational statistics. 2011; 3: 180-5.

[54] Wang YP, Tang HB, DeBarry JD et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012; 40: e49.

[55] Tang H, Wang X, Bowers JE et al. Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps. Genome Res. 2008; 18: 1944-54.

[56] Zhang Z, Xiao J, Wu J et al. ParaAT: a parallel tool for constructing multiple protein-coding DNA alignments. Biochem Biophys Res Commun. 2012; 419: 779-81.

[57] Wang D, Zhang Y, Zhang Z et al. KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinformatics. 2010; 8: 77-80.

[58] Murat F, Armero A, Pont C et al. Reconstructing the genome of the most recent common ancestor of flowering plants. Nat Genet. 2017; 49: 490-6.

[59] Chen SS, Black CC . Diurnal changes in volume and specific tissue weight of crassulacean acid metabolism plants. Plant Physiol. 1983; 71: 373-8.

[60] Bolger AM, Lohse M, Usadel B . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30: 2114-20.

[61] Trapnell C, Roberts A, Goff L et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012; 7: 562-78.

[62] Xu H, Yu Q, Shi Y et al. PGD: pineapple genomics database. Hortic Res. 2018; 5: 66.

[63] Goodstein DM, Shu S, Howson R et al. Phytozome: a comparative platform for green plant genomics. Nucleic Acids Res. 2012; 40: D1178-86.

[64] Kawahara Y, de la Bastide M, Hamilton JP et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice. 2013; 6: 4.

[65] Portwood JL, Woodhouse MR, Cannon EK et al. MaizeGDB 2018: the maize multi-genome genetics and genomics database. Nucleic Acids Res. 2019; 47: D1146-54.

[66] Marchler-Bauer A, Bryant SH . CD-search: protein domain annotations on the fly. Nucleic Acids Res. 2004; 32: W327-31.

[67] El-Gebali S, Mistry J, Bateman A et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2019; 47: D427-32.

[68] Mockler TC, Michael TP, Priest HD et al. The DIURNAL project: DIURNAL and circadian expression profiling, model-based pattern matching, and promoter analysis. Cold Spring Harb Symp Quant Biol. 2007; 72: 353-63.

[69] Chen LY, Xin Y, Wai CM et al. The role of cis-elements in the evolution of crassulacean acid metabolism photosynthesis. Hortic Res. 2020; 7: 5.

[70] Grant CE, Bailey TL, Noble WS . FIMO: scanning for occurrences of a given motif. Bioinformatics. 2011; 27: 1017-8.

[71] Michael TP, McClung CR . Phase-specific circadian clock regulatory elements in Arabidopsis. Plant Physiol. 2002; 130: 627-38.

[72] Freeling M. Bias in plant gene content following different sorts of duplication: tandem, whole-genome, segmental, or by transposition. Annu Rev Plant Biol. 2009; 60: 433-53.

[73] Wan T, Liu ZM, Leitch IJ et al. The Welwitschia genome reveals a unique biology underpinning extreme longevity in deserts. Nat Commun. 2021; 12: 4247.

[74] Xiong ZY, Gaeta RT, Pires JC . Homoeologous shuffling and chromosome compensation maintain genome balance in resynthesized allopolyploid Brassica napus. Proc Natl Acad Sci USA. 2011; 108: 7908-13.

[75] Qin LY, Hu YH, Wang JP et al. Insights into angiosperm evolution, floral development and chemical biosynthesis from the Aristolochia fimbriata genome. Nat Plants. 2021; 7: 1239-53.

[76] Shepherd T, Wynne GD . The effects of stress on plant cuticular waxes. New Phytol. 2006; 171: 469-99.

[77] Dubos C, Stracke R, Grotewold E et al. MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010; 15: 573-81.

[78] Gao J, Lan T . Functional characterization of the late embryogenesis abundant (LEA) protein gene family from Pinus tabuliformis (Pinaceae) in Escherichia coli. Sci Rep. 2016; 6: 19467.

[79] Jiang CK, Rao GY . Insights into the diversification and evolution of R2R3-MYB transcription factors in plants. Plant Physiol. 2020; 183: 637-55.

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

[81] Tobias PA, Guest DI . Tree immunity: growing old without antibodies. Trends Plant Sci. 2014; 19: 367-70.

[82] Mestre P, Baulcombe DC . Elicitor-mediated oligomerization of the tobacco N disease resistance protein. Plant Cell. 2006; 18: 491-501.

[83] Warren RL, Keeling CI, Yuen MM et al. Improved white spruce (Picea glauca) genome assemblies and annotation of large gene families of conifer terpenoid and phenolic defense metabolism. Plant J. 2015; 83: 189-212.

[84] Bondino HG, Valle EM, Ten Have A . Evolution and functional diversification of the small heat shock protein/α-crystallin family in higher plants. Planta. 2012; 235: 1299-313.

[85] Guo LM, Li J, He J et al. A class I cytosolic HSP20 of rice enhances heat and salt tolerance in different organisms. Sci Rep. 2020; 10: 1383.

[86] Ahl LI, Mravec J, Jorgensen B et al. Dynamics of intracellular mannan and cell wall folding in the drought responses of succulent aloe species. Plant Cell Environ. 2019; 42: 2458-71.

[87] Hanada K, Zou C, Lehti-Shiu MD et al. Importance of lineage-specific expansion of plant tandem duplicates in the adaptive response to environmental stimuli. Plant Physiol. 2008; 148: 993-1003.

[88] Rizzon C, Ponger L, Gaut BS . Striking similarities in the genomic distribution of tandemly arrayed genes in Arabidopsis and rice. PLoS Comput Biol. 2006; 2: 989-1000.

[89] Silvera K, Neubig KM, Whitten WM et al. Evolution along the crassulacean acid metabolism continuum. Funct Plant Biol. 2010; 37: 995-1010.

[90] Borland AM, Hartwell J, Weston DJ et al. Engineering crassulacean acid metabolism to improve water-use efficiency. Trends Plant Sci. 2014; 19: 327-38.

[91] Nelson EA, Sage RF . Functional constraints of CAM leaf anatomy: tight cell packing is associated with increased CAM function across a gradient of CAM expression. J Exp Bot. 2008; 59: 1841-50.

[92] Abraham PE, Castano NH, Cowan-Turner D et al. Peeling back the layers of crassulacean acid metabolism: functional differentiation between Kalanchoë fedtschenkoi epidermis and mesophyll proteomes. Plant J. 2020; 103: 869-88.

[93] Wickell D, Kuo LY, Yang HP et al. Underwater CAM photosynthesis elucidated by Isoetes genome. Nat Commun. 2021; 12: 6348.

[94] Michael TP, Mockler TC, Breton G et al. Network discovery pipeline elucidates conserved time-of-day-specific cis-regulatory modules. PLoS Genet. 2008; 4: e14.

[95] Mustilli AC, Merlot S, Vavasseur A et al. Arabidopsis OST1 protein kinase mediates the regulation of stomatal aperture by abscisic acid and acts upstream of reactive oxygen species production. Plant Cell. 2002; 14: 3089-99.

[96] Boxall SF, Kadu N, Dever LV et al. Kalanchoë PPC1 is essential for crassulacean acid metabolism and the regulation of core circadian clock and guard cell signaling genes. Plant Cell. 2020; 32: 1136-60.

[97] Abraham PE, Yin H, Borland AM et al. Transcript, protein and metabolite temporal dynamics in the CAM plant agave. Nat Plants. 2016; 2: 16178.

[98] Soltis DE, Soltis PS, Chase MW et al. Angiosperm phylogeny inferred from 18S rDNA, rbcL, and atpB sequences. Bot J Linn Soc. 2000; 133: 381-461.

[99] Chanderbali AS, Jin LL, Xu QJ et al. Buxus and Tetracentron genomes help resolve eudicot genome history. Nat Commun. 2022; 13: 1-10.

[100] Van de Peer Y, Maere S, Meyer A . The evolutionary significance of ancient genome duplications. Nat Rev Genet. 2009; 10: 725-32.

[101] Alzohairy AM, Yousef MA, Edris S et al. Detection of LTR retrotransposons reactivation induced by in vitro environmental stresses in barley (Hordeum vulgare) via RT-qPCR. Life Sci J. 2012; 9: 5019-26.

[102] Drake PL, Froend RH, Franks PJ . Smaller, faster stomata: scaling of stomatal size, rate of response, and stomatal conductance. J Exp Bot. 2013; 64: 495-505.

[103] Roddy AB, Théroux-Rancourt G, Abbo T et al. The scaling of genome size and cell size limits maximum rates of photosynthesis with implications for ecological strategies. Int J Plant Sci. 2020; 181: 75-87.

[104] Eggli U, Nyffeler R . Living under temporarily arid conditions - succulence as an adaptive strategy. Bradleya. 2009; 27: 13-36.

[105] Newton DJ, Chan J . South Africa’s trade in Southern African succulent plants. Johannesburg: TRAFFIC East/Southern Africa, 1998.

[106] Morse SR . Water-balance in Hemizonia luzulifolia - the role of extracellular polysaccharides. Plant Cell Environ. 1990; 13: 39-48.

[107] Hundertmark M, Hincha DK . LEA (late embryogenesis abundant) proteins and their encoding genes in Arabidopsis thaliana. BMC Genomics. 2008; 9: 118.

[108] Rai A, Singh R, Shirke PA et al. Expression of rice CYP450-like gene (Os08g01480) in Arabidopsis modulates regulatory network leading to heavy metal and other abiotic stress tolerance. PLoS One. 2015; 10: e0138574.

[109] Bloom AJ, Troughton JH . High productivity and photosynthetic flexibility in a CAM plant. Oecologia. 1979; 38: 35-43.

[110] Huang H, Nusinow DA . Into the evening: complex interactions in the Arabidopsis circadian clock. Trends Genet. 2016; 32: 674-86.