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.

Caught in the Act: Incipient Speciation at the Southern limit of Viburnum in the Central Andes.

A fundamental objective of evolutionary biology is to understand the origin of independently evolving species. Phylogenetic studies of species radiations rarely are able to document ongoing speciation; instead, modes of speciation, entailing geographic separation and/or ecological differentiation, are posited retrospectively. The Oreinotinus clade of Viburnum has radiated recently from north to south through the cloud forests of Mexico and Central America to the Central Andes. Our analyses support a hypothesis of incipient speciation in Oreinotinus at the southern edge of its geographic range, from central Peru to northern Argentina. Although several species and infraspecific taxa of have been recognized in this area, multiple lines of evidence and analytical approaches (including analyses of phylogenetic relationships, genetic structure, leaf morphology, and climatic envelopes) favor the recognition of just a single species, V. seemenii. We show that what has previously been recognized as V. seemenii f. minor has recently occupied the drier Tucuman-Bolivian forest region from Samaipata in Bolivia to Salta in northern Argentina. Plants in these populations form a well-supported clade with a distinctive genetic signature and they have evolved smaller, narrower leaves. We interpret this as the beginning of a within-species divergence process that has elsewhere in the neotropics resulted repeatedly in Viburnum species with a particular set of leaf ecomorphs. Specifically, the southern populations are in the process of evolving the small, glabrous, and entire leaf ecomorph that has evolved in four other montane areas of endemism. As predicted based on our studies of leaf ecomorphs in Chiapas, Mexico, these southern populations experience generally drier conditions, with large diurnal temperature fluctuations. In a central portion of the range of V. seemenii, characterized by wetter climatic conditions, we also document what may be the initial differentiation of the leaf ecomorph with larger, pubescent, and toothy leaves. The emergence of these ecomorphs thus appears to be driven by adaptation to subtly different climatic conditions in separate geographic regions, as opposed to parapatric differentiation along elevational gradients as suggested by Viburnum species distributions in other parts of the neotropics.

The genetics of convergent evolution: insights from plant photosynthesis.

The tree of life is resplendent with examples of convergent evolution, whereby distinct species evolve the same trait independently. Many highly convergent adaptations are also complex, which makes their repeated emergence surprising. In plants, the evolutionary history of two carbon concentrating mechanisms (CCMs) - C4 and crassulacean acid metabolism (CAM) photosynthesis - presents such a paradox. Both of these modifications of ancestral C3 photosynthesis require the integration of multiple anatomical and biochemical components, yet together they have evolved more than one hundred times. The presence of CCM enzymes in all plants suggests that a rudimentary CCM might emerge via relatively few genetic changes in potentiated lineages. Here, we propose that many of the complexities often associated with C4 and CAM photosynthesis may have evolved during a post-emergence optimization phase. The ongoing development of new model clades for young, emerging CCMs is enabling the comparative studies needed to test these ideas.

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.

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.

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.

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.

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.

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

Plant height and hydraulic vulnerability to drought and cold.

Understanding how plants survive drought and cold is increasingly important as plants worldwide experience dieback with drought in moist places and grow taller with warming in cold ones. Crucial in plant climate adaptation are the diameters of water-transporting conduits. Sampling 537 species across climate zones dominated by angiosperms, we find that plant size is unambiguously the main driver of conduit diameter variation. And because taller plants have wider conduits, and wider conduits within species are more vulnerable to conduction-blocking embolisms, taller conspecifics should be more vulnerable than shorter ones, a prediction we confirm with a plantation experiment. As a result, maximum plant size should be short under drought and cold, which cause embolism, or increase if these pressures relax. That conduit diameter and embolism vulnerability are inseparably related to plant size helps explain why factors that interact with conduit diameter, such as drought or warming, are altering plant heights worldwide.

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.

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.

Restriction-Site-Associated DNA Sequencing Reveals a Cryptic Viburnum Species on the North American Coastal Plain.

Species are the starting point for most studies of ecology and evolution, but the proper circumscription of species can be extremely difficult in morphologically variable lineages, and there are still few convincing examples of molecularly informed species delimitation in plants. Here, we focus on the Viburnum nudum complex, a highly variable clade that is widely distributed in eastern North America. Taxonomic treatments have mostly divided this complex into northern (V. nudum var. cassinoides) and southern (V. nudum var. nudum) entities, but additional names have been proposed. We used multiple lines of evidence, including RADseq, morphological, and geographic data, to test how many independently evolving lineages exist within the V. nudum complex. Genetic clustering and phylogenetic methods revealed three distinct groups-one lineage that is highly divergent, and two others that are recently diverged and morphologically similar. A combination of evidence that includes reciprocal monophyly, lack of introgression, and discrete rather than continuous patterns of variation supports the recognition of all three lineages as separate species. These results identify a surprising case of cryptic diversity in which two broadly sympatric species have consistently been lumped in taxonomic treatments. The clarity of our findings is directly related to the dense sampling and high-quality genetic data in this study. We argue that there is a critical need for carefully sampled and integrative species delimitation studies to clarify species boundaries even in well-known plant lineages. Studies following the model that we have developed here are likely to identify many more cryptic lineages and will fundamentally improve our understanding of plant speciation and patterns of species richness.

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.

Different clades and traits yield similar grassland functional responses.

Plant functional traits are viewed as key to predicting important ecosystem and community properties across resource gradients within and among biogeographic regions. Vegetation dynamics and ecosystem processes, such as aboveground net primary productivity (ANPP), are increasingly being modeled as a function of the quantitative traits of species, which are used as proxies for photosynthetic rates and nutrient and water-use efficiency. These approaches rely on an assumption that a certain trait value consistently confers a specific function or response under given environmental conditions. Here, we provide a critical test of this idea and evaluate whether the functional traits that drive the well-known relationship between precipitation and ANPP differ between systems with distinct biogeographic histories and species assemblages. Specifically, we compared grasslands spanning a broad precipitation gradient (∼200-1,000 mm/y) in North America and South Africa that differ in the relative representation and abundance of grass phylogenetic lineages. We found no significant difference between the regions in the positive relationship between annual precipitation and ANPP, yet the trait values underlying this relationship differed dramatically. Our results challenge the trait-based approach to predicting ecosystem function by demonstrating that different combinations of functional traits can act to maximize ANPP in a given environmental setting. Further, we show the importance of incorporating biogeographic and phylogenetic history in predicting community and ecosystem properties using traits.

Evolutionary trajectories, accessibility and other metaphors: the case of C4 and CAM photosynthesis.

Are evolutionary outcomes predictable? Adaptations that show repeated evolutionary convergence across the Tree of Life provide a special opportunity to dissect the context surrounding their origins, and identify any commonalities that may predict why certain traits evolved many times in particular clades and yet never evolved in others. The remarkable convergence of C4 and Crassulacean Acid Metabolism (CAM) photosynthesis in vascular plants makes them exceptional model systems for understanding the repeated evolution of complex phenotypes. This review highlights what we have learned about the recurring assembly of C4 and CAM, focusing on the increasingly predictable stepwise evolutionary integration of anatomy and biochemistry. With the caveat that we currently understand C4 evolution better than we do CAM, I propose a general model that explains and unites C4 and CAM evolutionary trajectories. Available data suggest that anatomical modifications are the 'rate-limiting step' in each trajectory, which in large part determines the evolutionary accessibility of both syndromes. The idea that organismal structure exerts a primary influence on innovation is discussed in the context of other systems. Whether the rate-limiting step occurs early or late in the evolutionary assembly of a new phenotype may have profound implications for its distribution across the Tree of Life.

Facultative crassulacean acid metabolism in a C3-C4 intermediate.

The Portulacaceae enable the study of the evolutionary relationship between C4 and crassulacean acid metabolism (CAM) photosynthesis. Shoots of well-watered plants of the C3-C4 intermediate species Portulaca cryptopetala Speg. exhibit net uptake of CO2 solely during the light. CO2 fixation is primarily via the C3 pathway as indicated by a strong stimulation of CO2 uptake when shoots were provided with air containing 2% O2. When plants were subjected to water stress, daytime CO2 uptake was reduced and CAM-type net CO2 uptake in the dark occurred. This was accompanied by nocturnal accumulation of acid in both leaves and stems, also a defining characteristic of CAM. Following rewatering, net CO2 uptake in the dark ceased in shoots, as did nocturnal acidification of the leaves and stems. With this unequivocal demonstration of stress-related reversible, i.e. facultative, induction of CAM, P. cryptopetala becomes the first C3-C4 intermediate species reported to exhibit CAM. Portulaca molokiniensis Hobdy, a C4 species, also exhibited CAM only when subjected to water stress. Facultative CAM has now been demonstrated in all investigated species of Portulaca, which are well sampled from across the phylogeny. This strongly suggests that in Portulaca, a lineage in which species engage predominately in C4 photosynthesis, facultative CAM is ancestral to C4. In a broader context, it has now been demonstrated that CAM can co-exist in leaves that exhibit any of the other types of photosynthesis known in terrestrial plants: C3, C4 and C3-C4 intermediate.

Targeted Enrichment of Large Gene Families for Phylogenetic Inference: Phylogeny and Molecular Evolution of Photosynthesis Genes in the Portullugo Clade (Caryophyllales).

Hybrid enrichment is an increasingly popular approach for obtaining hundreds of loci for phylogenetic analysis across many taxa quickly and cheaply. The genes targeted for sequencing are typically single-copy loci, which facilitate a more straightforward sequence assembly and homology assignment process. However, this approach limits the inclusion of most genes of functional interest, which often belong to multi-gene families. Here, we demonstrate the feasibility of including large gene families in hybrid enrichment protocols for phylogeny reconstruction and subsequent analyses of molecular evolution, using a new set of bait sequences designed for the "portullugo" (Caryophyllales), a moderately sized lineage of flowering plants (~ 2200 species) that includes the cacti and harbors many evolutionary transitions to C$_{\mathrm{4}}$ and CAM photosynthesis. Including multi-gene families allowed us to simultaneously infer a robust phylogeny and construct a dense sampling of sequences for a major enzyme of C$_{\mathrm{4}}$ and CAM photosynthesis, which revealed the accumulation of adaptive amino acid substitutions associated with C$_{\mathrm{4}}$ and CAM origins in particular paralogs. Our final set of matrices for phylogenetic analyses included 75-218 loci across 74 taxa, with ~ 50% matrix completeness across data sets. Phylogenetic resolution was greatly improved across the tree, at both shallow and deep levels. Concatenation and coalescent-based approaches both resolve the sister lineage of the cacti with strong support: Anacampserotaceae $+$ Portulacaceae, two lineages of mostly diminutive succulent herbs of warm, arid regions. In spite of this congruence, BUCKy concordance analyses demonstrated strong and conflicting signals across gene trees. Our results add to the growing number of examples illustrating the complexity of phylogenetic signals in genomic-scale data.

Differences in flowering time maintain species boundaries in a continental radiation of Viburnum.

PREMISE: We take an integrative approach in assessing how introgression and Pleistocene climate fluctuations have shaped the diversification of the core Lentago clade of Viburnum, a group of five interfertile species with broad areas of sympatry. We specifically tested whether flowering time plays a role in maintaining species isolation. METHODS: RAD-seq data for 103 individuals were used to infer the species relationships and the genetic structure within each species. Flowering times were compared among species on the basis of historical flowering dates documented by herbarium specimens. RESULTS: Within each species, we found a strong relationship between flowering date and latitude, such that southern populations flower earlier than northern ones. In areas of sympatry, the species flower in sequence rather than simultaneously, with flowering dates offset by ≥9 d for all species pairs. In two cases it appears that the offset in flowering times is an incidental consequence of adaptation to differing climates, but in the recently diverged sister species V. prunifolium and V. rufidulum, we find evidence that reinforcement led to reproductive character displacement. Long-term trends suggest that the two northern-most species are flowering earlier in response to recent climate change. CONCLUSIONS: We argue that speciation in the Lentago clade has primarily occurred through ecological divergence of allopatric populations, but differences in flowering time were essential to maintain separation of incipient species when they came into secondary contact. This combination of factors may underlie diversification in many other plant clades.