Predicting photosynthetic pathway from anatomy using machine learning.

Plants with Crassulacean acid metabolism (CAM) have long been associated with a specialized anatomy, including succulence and thick photosynthetic tissues. Firm, quantitative boundaries between non-CAM and CAM plants have yet to be established - if they indeed exist. Using novel computer vision software to measure anatomy, we combined new measurements with published data across flowering plants. We then used machine learning and phylogenetic comparative methods to investigate relationships between CAM and anatomy. We found significant differences in photosynthetic tissue anatomy between plants with differing CAM phenotypes. Machine learning-based classification was over 95% accurate in differentiating CAM from non-CAM anatomy, and had over 70% recall of distinct CAM phenotypes. Phylogenetic least squares regression and threshold analyses revealed that CAM evolution was significantly correlated with increased mesophyll cell size, thicker leaves, and decreased intercellular airspace. Our findings suggest that machine learning may be used to aid the discovery of new CAM species and that the evolutionary trajectory from non-CAM to strong, obligate CAM requires continual anatomical specialization.

Reconciling continuous and discrete models of C4 and CAM evolution.

BACKGROUND: A current argument in the CAM biology literature has focused on the nature of the CAM evolutionary trajectory: whether there is a smooth continuum of phenotypes between plants with C3 and CAM photosynthesis or whether there are discrete steps of phenotypic evolutionary change such as has been modelled for the evolution of C4 photosynthesis. A further implication is that a smooth continuum would increase the evolvability of CAM, whereas discrete changes would make the evolutionary transition from C3 to CAM more difficult. SCOPE: In this essay, I attempt to reconcile these two viewpoints, because I think in many ways this is a false dichotomy that is constraining progress in understanding how both CAM and C4 evolved. In reality, the phenotypic space connecting C3 species and strong CAM/C4 species is both a continuum of variably expressed quantitative traits and yet also contains certain combinations of traits that we are able to identify as discrete, recognizable phenotypes. In this sense, the evolutionary mechanics of CAM origination are no different from those of C4 photosynthesis, nor from the evolution of any other complex trait assemblage. CONCLUSIONS: To make progress, we must embrace the concept of discrete phenotypic phases of CAM evolution, because their delineation will force us to articulate what aspects of phenotypic variation we think are significant. There are some current phenotypic gaps that are limiting our ability to build a complete CAM evolutionary model: the first is how a rudimentary CAM biochemical cycle becomes established, and the second is how the 'accessory' CAM cycle in C3+CAM plants is recruited into a primary metabolism. The connections to the C3 phenotype we are looking for are potentially found in the behaviour of C3 plants when undergoing physiological stress - behaviour that, strangely enough, remains essentially unexplored in this context.

The CAM lineages of planet Earth.

BACKGROUND AND SCOPE: The growth of experimental studies of crassulacean acid metabolism (CAM) in diverse plant clades, coupled with recent advances in molecular systematics, presents an opportunity to re-assess the phylogenetic distribution and diversity of species capable of CAM. It has been more than two decades since the last comprehensive lists of CAM taxa were published, and an updated survey of the occurrence and distribution of CAM taxa is needed to facilitate and guide future CAM research. We aimed to survey the phylogenetic distribution of these taxa, their diverse morphology, physiology and ecology, and the likely number of evolutionary origins of CAM based on currently known lineages. RESULTS AND CONCLUSIONS: We found direct evidence (in the form of experimental or field observations of gas exchange, day-night fluctuations in organic acids, carbon isotope ratios and enzymatic activity) for CAM in 370 genera of vascular plants, representing 38 families. Further assumptions about the frequency of CAM species in CAM clades and the distribution of CAM in the Cactaceae and Crassulaceae bring the currently estimated number of CAM-capable species to nearly 7 % of all vascular plants. The phylogenetic distribution of these taxa suggests a minimum of 66 independent origins of CAM in vascular plants, possibly with dozens more. To achieve further insight into CAM origins, there is a need for more extensive and systematic surveys of previously unstudied lineages, particularly in living material to identify low-level CAM activity, and for denser sampling to increase phylogenetic resolution in CAM-evolving clades. This should allow further progress in understanding the functional significance of this pathway by integration with studies on the evolution and genomics of CAM in its many forms.

Atmospheric CO2 decline and the timing of CAM plant evolution.

BACKGROUND AND AIMS: CAM photosynthesis is hypothesized to have evolved in atmospheres of low CO2 concentration in recent geological time because of its ability to concentrate CO2 around Rubisco and boost water use efficiency relative to C3 photosynthesis. We assess this hypothesis by compiling estimates of when CAM clades arose using phylogenetic chronograms for 73 CAM clades. We further consider evidence of how atmospheric CO2 affects CAM relative to C3 photosynthesis. RESULTS: Where CAM origins can be inferred, strong CAM is estimated to have appeared in the past 30 million years in 46 of 48 examined clades, after atmospheric CO2 had declined from high (near 800 ppm) to lower (<450 ppm) values. In turn, 21 of 25 clades containing CAM species (but where CAM origins are less certain) also arose in the past 30 million years. In these clades, CAM is probably younger than the clade origin. We found evidence for repeated weak CAM evolution during the higher CO2 conditions before 30 million years ago, and possible strong CAM origins in the Crassulaceae during the Cretaceous period prior to atmospheric CO2 decline. Most CAM-specific clades arose in the past 15 million years, in a similar pattern observed for origins of C4 clades. CONCLUSIONS: The evidence indicates strong CAM repeatedly evolved in reduced CO2 conditions of the past 30 million years. Weaker CAM can pre-date low CO2 and, in the Crassulaceae, strong CAM may also have arisen in water-limited microsites under relatively high CO2. Experimental evidence from extant CAM species demonstrates that elevated CO2 reduces the importance of nocturnal CO2 fixation by increasing the contribution of C3 photosynthesis to daily carbon gain. Thus, the advantage of strong CAM would be reduced in high CO2, such that its evolution appears less likely and restricted to more extreme environments than possible in low CO2.

Spatial resolution of an integrated C4+CAM photosynthetic metabolism.

C4 and CAM photosynthesis have repeatedly evolved in plants over the past 30 million years. Because both repurpose the same set of enzymes but differ in their spatial and temporal deployment, they have long been considered as distinct and incompatible adaptations. Portulaca contains multiple C4 species that perform CAM when droughted. Spatially explicit analyses of gene expression reveal that C4 and CAM systems are completely integrated in Portulaca oleracea, with CAM and C4 carbon fixation occurring in the same cells and CAM-generated metabolites likely incorporated directly into the C4 cycle. Flux balance analysis corroborates the gene expression findings and predicts an integrated C4+CAM system under drought. This first spatially explicit description of a C4+CAM photosynthetic metabolism presents a potential new blueprint for crop improvement.

Replicated radiation of a plant clade along a cloud forest archipelago.

Replicated radiations, in which sets of similar forms evolve repeatedly within different regions, can provide powerful insights into parallel evolution and the assembly of functional diversity within communities. Several cases have been described in animals, but in plants we lack well-documented cases of replicated radiation that combine comprehensive phylogenetic and biogeographic analyses, the delimitation of geographic areas within which a set of 'ecomorphs' evolved independently and the identification of potential underlying mechanisms. Here we document the repeated evolution of a set of leaf ecomorphs in a group of neotropical plants. The Oreinotinus lineage within the angiosperm clade Viburnum spread from Mexico to Argentina through disjunct cloud forest environments. In 9 of 11 areas of endemism, species with similar sets of leaf forms evolved in parallel. We reject gene-flow-mediated evolution of similar leaves and show, instead, that species with disparate leaf forms differ in their climatic niches, supporting ecological adaptation as the driver of parallelism. Our identification of a case of replicated radiation in plants sets the stage for comparative analyses of such phenomena across the tree of life.

Gene co-expression reveals the modularity and integration of C4 and CAM in Portulaca.

C4 photosynthesis and Crassulacean acid metabolism (CAM) have been considered as largely independent adaptations despite sharing key biochemical modules. Portulaca is a geographically widespread clade of over 100 annual and perennial angiosperm species that primarily use C4 but facultatively exhibit CAM when drought stressed, a photosynthetic system known as C4 + CAM. It has been hypothesized that C4 + CAM is rare because of pleiotropic constraints, but these have not been deeply explored. We generated a chromosome-level genome assembly of Portulaca amilis and sampled mRNA from P. amilis and Portulaca oleracea during CAM induction. Gene co-expression network analyses identified C4 and CAM gene modules shared and unique to both Portulaca species. A conserved CAM module linked phosphoenolpyruvate carboxylase to starch turnover during the day-night transition and was enriched in circadian clock regulatory motifs in the P. amilis genome. Preservation of this co-expression module regardless of water status suggests that Portulaca constitutively operate a weak CAM cycle that is transcriptionally and posttranscriptionally upregulated during drought. C4 and CAM mostly used mutually exclusive genes for primary carbon fixation, and it is likely that nocturnal CAM malate stores are shuttled into diurnal C4 decarboxylation pathways, but we found evidence that metabolite cycling may occur at low levels. C4 likely evolved in Portulaca through co-option of redundant genes and integration of the diurnal portion of CAM. Thus, the ancestral CAM system did not strongly constrain C4 evolution because photosynthetic gene networks are not co-regulated for both daytime and nighttime functions.

CAM photosynthesis in desert blooming Cistanthe of the Atacama, Chile.

When plants of the Atacama desert undergo episodic blooms, among the most prominent are succulent-leaved Cistanthe (Montiaceae). We demonstrate that two Cistanthe species, the perennial Cistanthe sp. aff. crassifolia and the annual/biannual Cistanthe sp. aff. longiscapa, can exhibit net CO2 uptake and leaf acidification patterns typical of crassulacean acid metabolism (CAM). In C. sp. aff. crassifolia leaves, CAM expression was facultative. CAM-type nocturnal net CO2 uptake and acid accumulation occurred in drought-stressed but not in well-watered plants. By contrast, CAM expression in C. sp. aff. longiscapa was largely constitutive. Nocturnal acid accumulation was present in leaves of well-watered and in droughted plants. Following water-deficit stress, net nocturnal CO2 uptake was induced and the level of acid accumulated increased. Neither nocturnal CO2 uptake nor acid accumulation was reduced when the plants were re-watered. δ13C values of a further nine field-collected Cistanthe species are consistent with a contribution of CAM to their carbon pools. In the Portulacinae, a suborder with eight CAM-containing families, Cistanthe becomes the sixth genus with CAM within the family Montiaceae, and it is likely that the ancestor of all Portulacineae also possessed CAM photosynthesis. In the stochastic rainfall landscape of the Atacama, carbon uptake in the dark is a water-use efficient mechanism that increases the carbon pool available for seed production or dormancy. The next rain event may be years away.

Developmental and biophysical determinants of grass leaf size worldwide.

One of the most notable ecological trends-described more than 2,300  years ago by Theophrastus-is the association of small leaves with dry and cold climates, which has recently been recognized for eudicotyledonous plants at a global scale1-3. For eudicotyledons, this pattern has been attributed to the fact that small leaves have a thinner boundary layer that helps to avoid extreme leaf temperatures4 and their leaf development results in vein traits that improve water transport under cold or dry climates5,6. However, the global distribution of leaf size and its adaptive basis have not been tested in the grasses, which represent a diverse lineage that is distinct in leaf morphology and that contributes 33% of terrestrial primary productivity (including the bulk of crop production)7. Here we demonstrate that grasses have shorter and narrower leaves under colder and drier climates worldwide. We show that small grass leaves have thermal advantages and vein development that contrast with those of eudicotyledons, but that also explain the abundance of small leaves in cold and dry climates. The worldwide distribution of leaf size in grasses exemplifies how biophysical and developmental processes result in convergence across major lineages in adaptation to climate globally, and highlights the importance of leaf size and venation architecture for grass performance in past, present and future ecosystems.

Modeling Phylogenetic Biome Shifts on a Planet with a Past.

The spatial distribution of biomes has changed considerably over deep time, so the geographical opportunity for an evolutionary lineage to shift into a new biome may depend on how the availability and connectivity of biomes has varied temporally. To better understand how lineages shift between biomes in space and time, we developed a phylogenetic biome shift model in which each lineage shifts between biomes and disperses between regions at rates that depend on the lineage's biome affinity and location relative to the spatial distribution of biomes at any given time. To study the behavior of the biome shift model in an empirical setting, we developed a literature-based representation of paleobiome structure for three mesic forest biomes, six regions, and eight time strata, ranging from the Late Cretaceous (100 Ma) through the present. We then fitted the model to a time-calibrated phylogeny of 119 Viburnum species to compare how the results responded to various realistic or unrealistic assumptions about paleobiome structure. Ancestral biome estimates that account for paleobiome dynamics reconstructed a warm temperate (or tropical) origin of Viburnum, which is consistent with previous fossil-based estimates of ancestral biomes. Imposing unrealistic paleobiome distributions led to ancestral biome estimates that eliminated support for tropical origins, and instead inflated support for cold temperate ancestry throughout the warmer Paleocene and Eocene. The biome shift model we describe is applicable to the study of evolutionary systems beyond Viburnum, and the core mechanisms of our model are extensible to the design of richer phylogenetic models of historical biogeography and/or lineage diversification. We conclude that biome shift models that account for dynamic geographical opportunities are important for inferring ancestral biomes that are compatible with our understanding of Earth history.[Ancestral states; biome shifts; historical biogeography; niche conservatism; phylogenetics].

Joint Phylogenetic Estimation of Geographic Movements and Biome Shifts during the Global Diversification of Viburnum.

Phylogeny, molecular sequences, fossils, biogeography, and biome occupancy are all lines of evidence that reflect the singular evolutionary history of a clade, but they are most often studied separately, by first inferring a fossil-dated molecular phylogeny, then mapping on ancestral ranges and biomes inferred from extant species. Here we jointly model the evolution of biogeographic ranges, biome affinities, and molecular sequences, while incorporating fossils to estimate a dated phylogeny for all of the 163 extant species of the woody plant clade Viburnum (Adoxaceae) that we currently recognize in our ongoing worldwide monographic treatment of the group. Our analyses indicate that while the major Viburnum lineages evolved in the Eocene, the majority of extant species originated since the Miocene. Viburnum radiated first in Asia, in warm, broad-leaved evergreen (lucidophyllous) forests. Within Asia, we infer several early shifts into more tropical forests, and multiple shifts into forests that experience prolonged freezing. From Asia, we infer two early movements into the New World. These two lineages probably first occupied warm temperate forests and adapted later to spreading cold climates. One of these lineages (Porphyrotinus) occupied cloud forests and moved south through the mountains of the Neotropics. Several other movements into North America took place more recently, facilitated by prior adaptations to freezing in the Old World. We also infer four disjunctions between Asia and Europe: the Tinus lineage is the oldest and probably occupied warm forests when it spread, whereas the other three were more recent and in cold-adapted lineages. These results variously contradict published accounts, especially the view that Viburnum radiated initially in cold forests and, accordingly, maintained vessel elements with scalariform perforations. We explored how the location and biome assignments of fossils affected our inference of ancestral areas and biome states. Our results are sensitive to, but not entirely dependent upon, the inclusion of fossil biome data. It will be critical to take advantage of all available lines of evidence to decipher events in the distant past. The joint estimation approach developed here provides cautious hope even when fossil evidence is limited. [Biogeography; biome; combined evidence; fossil pollen; phylogeny; Viburnum.].

Evolutionary dynamics of genome size in a radiation of woody plants.

PREMISE: Plant genome size ranges widely, providing many opportunities to examine how genome size variation affects plant form and function. We analyzed trends in chromosome number, genome size, and leaf traits for the woody angiosperm clade Viburnum to examine the evolutionary associations, functional implications, and possible drivers of genome size. METHODS: Chromosome counts and genome size estimates were mapped onto a Viburnum phylogeny to infer the location and frequency of polyploidization events and trends in genome size evolution. Genome size was analyzed with leaf anatomical and physiological data to evaluate the influence of genome size on plant function. RESULTS: We discovered nine independent polyploidization events, two reductions in base chromosome number, and substantial variation in genome size with a slight trend toward genome size reduction in polyploids. We did not find strong relationships between genome size and the functional and morphological traits that have been highlighted at broader phylogenetic scales. CONCLUSIONS: Polyploidization events were sometimes associated with rapid radiations, demonstrating that polyploid lineages can be highly successful. Relationships between genome size and plant physiological function observed at broad phylogenetic scales may be largely irrelevant to the evolutionary dynamics of genome size at smaller scales. The view that plants readily tolerate changes in ploidy and genome size, and often do so, appears to apply to Viburnum.

Diversification in evolutionary arenas-Assessment and synthesis.

Understanding how and why rates of evolutionary diversification vary is a key issue in evolutionary biology, ecology, and biogeography. Evolutionary rates are the net result of interacting processes summarized under concepts such as adaptive radiation and evolutionary stasis. Here, we review the central concepts in the evolutionary diversification literature and synthesize these into a simple, general framework for studying rates of diversification and quantifying their underlying dynamics, which can be applied across clades and regions, and across spatial and temporal scales. Our framework describes the diversification rate (d) as a function of the abiotic environment (a), the biotic environment (b), and clade-specific phenotypes or traits (c); thus, d ~ a,b,c. We refer to the four components (a-d) and their interactions collectively as the "Evolutionary Arena." We outline analytical approaches to this framework and present a case study on conifers, for which we parameterize the general model. We also discuss three conceptual examples: the Lupinus radiation in the Andes in the context of emerging ecological opportunity and fluctuating connectivity due to climatic oscillations; oceanic island radiations in the context of island formation and erosion; and biotically driven radiations of the Mediterranean orchid genus Ophrys. The results of the conifer case study are consistent with the long-standing scenario that low competition and high rates of niche evolution promote diversification. The conceptual examples illustrate how using the synthetic Evolutionary Arena framework helps to identify and structure future directions for research on evolutionary radiations. In this way, the Evolutionary Arena framework promotes a more general understanding of variation in evolutionary rates by making quantitative results comparable between case studies, thereby allowing new syntheses of evolutionary and ecological processes to emerge.

Crassulacean acid metabolism.

Gilman and Edwards introduce crassulacean acid metabolism and highlight how recent advances in molecular biology are deepening our knowledge of CAM evolution.

Lineage-based functional types: characterising functional diversity to enhance the representation of ecological behaviour in Land Surface Models.

Process-based vegetation models attempt to represent the wide range of trait variation in biomes by grouping ecologically similar species into plant functional types (PFTs). This approach has been successful in representing many aspects of plant physiology and biophysics but struggles to capture biogeographic history and ecological dynamics that determine biome boundaries and plant distributions. Grass-dominated ecosystems are broadly distributed across all vegetated continents and harbour large functional diversity, yet most Land Surface Models (LSMs) summarise grasses into two generic PFTs based primarily on differences between temperate C3 grasses and (sub)tropical C4 grasses. Incorporation of species-level trait variation is an active area of research to enhance the ecological realism of PFTs, which form the basis for vegetation processes and dynamics in LSMs. Using reported measurements, we developed grass functional trait values (physiological, structural, biochemical, anatomical, phenological, and disturbance-related) of dominant lineages to improve LSM representations. Our method is fundamentally different from previous efforts, as it uses phylogenetic relatedness to create lineage-based functional types (LFTs), situated between species-level trait data and PFT-level abstractions, thus providing a realistic representation of functional diversity and opening the door to the development of new vegetation models.

C4 and crassulacean acid metabolism within a single leaf: deciphering key components behind a rare photosynthetic adaptation.

Although biochemically related, C4 and crassulacean acid metabolism (CAM) systems are expected to be incompatible. However, Portulaca species, including P. oleracea, operate C4 and CAM within a single leaf, and the mechanisms behind this unique photosynthetic arrangement remain largely unknown. Here, we employed RNA-seq to identify candidate genes involved exclusively or shared by C4 or CAM, and provided an in-depth characterization of their transcript abundance patterns during the drought-induced photosynthetic transitions in P. oleracea. Data revealed fewer candidate CAM-specific genes than those recruited to function in C4 . The putative CAM-specific genes were predominantly involved in night-time primary carboxylation reactions and malate movement across the tonoplast. Analysis of gene transcript-abundance regulation and photosynthetic physiology indicated that C4 and CAM coexist within a single P. oleracea leaf under mild drought conditions. Developmental and environmental cues were shown to regulate CAM expression in stems, whereas the shift from C4 to C4 -CAM hybrid photosynthesis in leaves was strictly under environmental control. Moreover, efficient starch turnover was identified as part of the metabolic adjustments required for CAM operation in both organs. These findings provide insights into C4 /CAM connectivity and compatibility, contributing to a deeper understanding of alternative ways to engineer CAM into C4 crop species.

Celebrating a New Division of Botany at SICB: An Introduction to the Integrative Plant Biology Symposium.

The Society for Integrative and Comparative Biology (SICB) should, in theory, be a home for scientists working across the entire Tree of Life. In practice, SICB has remained principally a society that supports integrative zoological research. Here we highlight a broad collection of what we consider to the best in integrative and comparative plant biology, gathered together for a special symposium at the 2019 SICB meeting. This symposium and special issue mark the initiation of a new Division of Botany within SICB, which we hope will usher in a new era of SICB where botanists and zoologists engage, collaborate, and celebrate together in this especially creative period of integrative and comparative biology.

The Evolution of CAM Photosynthesis in Australian Calandrinia Reveals Lability in C3+CAM Phenotypes and a Possible Constraint to the Evolution of Strong CAM.

Australian Calandrinia has radiated across the Australian continent during the last 30 Ma, and today inhabits most Australian ecosystems. Given its biogeographic range and reports of facultative Crassulacean acid metabolism (CAM) photosynthesis in multiple species, we hypothesized (1) that CAM would be widespread across Australian Calandrinia and that species, especially those that live in arid regions, would engage in strong CAM, and (2) that Australian Calandrinia would be an important lineage for informing on the CAM evolutionary trajectory. We cultivated 22 Australian Calandrinia species for a drought experiment. Using physiological measurements and δ13C values we characterized photosynthetic mode across these species, mapped the resulting character states onto a phylogeny, and characterized the climatic envelopes of species in their native ranges. Most species primarily utilize C3 photosynthesis, with CAM operating secondarily, often upregulated following drought. Several phylogenetically nested species are C3, indicating evolutionary losses of CAM. No strong CAM was detected in any of the species. Results highlight the limitations of δ13C surveys in detecting C3+CAM phenotypes, and the evolutionary lability of C3+CAM phenotypes. We propose a model of CAM evolution that allows for lability and reversibility among C3+CAM phenotypes and C3 and suggest that an annual life-cycle may preclude the evolution of strong CAM.