Ancient and recent introgression shape the evolutionary history of pollinator adaptation and speciation in a model monkeyflower radiation (Mimulus section Erythranthe)

Patterns of genomic variation reveal a single evolutionary origin of the wild allotetraploid Mimulus sookensis

Polyploidy occurs across the tree of life and is especially common in plants. Because newly formed cytotypes are often incompatible with their progenitors, polyploidy is also said to trigger "instantaneous" speciation. If a polyploid can self-fertilize or reproduce asexually, it is even possible for one individual to produce an entirely new lineage, but how often this scenario occurs is unclear. Here, we investigate the evolutionary history of the wild allotetraploid Mimulus sookensis, which was formed through hybridization between self-compatible, diploid species in the Mimulus guttatus complex. We generate a chromosome-scale reference assembly for M. sookensis and define its distinct subgenomes. Despite previous reports suggesting multiple origins of this highly selfing polyploid, we discover patterns of population genomic variation that provide unambiguous support for a single origin. One M. sookensis subgenome is clearly derived from the selfer Mimulus nasutus, which organellar variation suggests is the maternal progenitor. The ancestor of the other subgenome is less certain, but it shares variation with both Mimulus decorus and M. guttatus, two outcrossing diploids with geographic ranges that overlap broadly with M. sookensis. This study establishes M. sookensis as an example of instantaneous speciation, likely facilitated by the polyploid's predisposition to self-fertilize.

A distinct foliar pigmentation pattern formed by activator-repressor gradients upstream of an anthocyanin-activating R2R3-MYB

The emergence of novel traits is often preceded by a potentiation phase, when all the genetic components necessary for producing the trait are assembled. However, elucidating these potentiating factors is challenging. We have previously shown that an anthocyanin-activating R2R3-MYB, STRIPY, triggers the emergence of a distinct foliar pigmentation pattern in the monkeyflower Mimulus verbenaceus. Here, using forward and reverse genetics approaches, we identify three potentiating factors that pattern STRIPY expression: MvHY5, a master regulator of light signaling that activates STRIPY and is expressed throughout the leaf, and two leaf developmental regulators, MvALOG1 and MvTCP5, that are expressed in opposing gradients along the leaf proximodistal axis and negatively regulate STRIPY. These results provide strong empirical evidence that phenotypic novelties can be potentiated through incorporation into preexisting genetic regulatory networks and highlight the importance of positional information in patterning the novel foliar stripe.

Phylogenomic perspectives on speciation and reproductive isolation in a North American biodiversity hotspot: an example using California sages (Salvia subgenus Audibertia: Lamiaceae)

BACKGROUND AND AIMS

The California Floristic Province (CA-FP) is the most species-rich region of North America north of Mexico. One of several proposed hypotheses explaining the exceptional diversity of the region is that the CA-FP harbours myriad recently diverged lineages with nascent reproductive barriers. Salvia subgenus Audibertia is a conspicuous element of the CA-FP, with multiple sympatric and compatible species.

METHODS

Using 305 nuclear loci and both organellar genomes, we reconstruct species trees, examine genomic discordance, conduct divergence-time estimation, and analyse contemporaneous patterns of gene flow and mechanical reproductive isolation.

KEY RESULTS

Despite strong genomic discordance, an underlying bifurcating tree is supported. Organellar genomes capture additional introgression events not detected in the nuclear genome. Most interfertility is found within clades, indicating that reproductive barriers arise with increasing genetic divergence. Species are generally not mechanically isolated, suggesting that it is unlikely to be the primary factor leading to reproductive isolation.

CONCLUSIONS

Rapid, recent speciation with some interspecific gene flow in conjunction with the onset of a Mediterranean-like climate is the underlying cause of extant diversity in Salvia subgenus Audibertia. Speciation has largely not been facilitated by gene flow. Its signal in the nuclear genome seems to mostly be erased by backcrossing, but organellar genomes each capture different instances of historical gene flow, probably characteristic of many CA-FP lineages. Mechanical reproductive isolation appears to be only part of a mosaic of factors limiting gene flow.

Toward telomere-to-telomere cat genomes for precision medicine and conservation biology

Genomic data from species of the cat family Felidae promise to stimulate veterinary and human medical advances, and clarify the coherence of genome organization. We describe how interspecies hybrids have been instrumental in the genetic analysis of cats, from the first genetic maps to propelling cat genomes toward the T2T standard set by the human genome project. Genotype-to-phenotype mapping in cat models has revealed dozens of health-related genetic variants, the molecular basis for mammalian pigmentation and patterning, and species-specific adaptations. Improved genomic surveillance of natural and captive populations across the cat family tree will increase our understanding of the genetic architecture of traits, population dynamics, and guide a future of genome-enabled biodiversity conservation.

Macroevolution of the plant-hummingbird pollination system

Plant-hummingbird interactions are considered a classic example of coevolution, a process in which mutually dependent species influence each other's evolution. Plants depend on hummingbirds for pollination, whereas hummingbirds rely on nectar for food. As a step towards understanding coevolution, this review focuses on the macroevolutionary consequences of plant-hummingbird interactions, a relatively underexplored area in the current literature. We synthesize prior studies, illustrating the origins and dynamics of hummingbird pollination across different angiosperm clades previously pollinated by insects (mostly bees), bats, and passerine birds. In some cases, the crown age of hummingbirds pre-dates the plants they pollinate. In other cases, plant groups transitioned to hummingbird pollination early in the establishment of this bird group in the Americas, with the build-up of both diversities coinciding temporally, and hence suggesting co-diversification. Determining what triggers shifts to and away from hummingbird pollination remains a major open challenge. The impact of hummingbirds on plant diversification is complex, with many tropical plant lineages experiencing increased diversification after acquiring flowers that attract hummingbirds, and others experiencing no change or even a decrease in diversification rates. This mixed evidence suggests that other extrinsic or intrinsic factors, such as local climate and isolation, are important covariables driving the diversification of plants adapted to hummingbird pollination. To guide future studies, we discuss the mechanisms and contexts under which hummingbirds, as a clade and as individual species (e.g. traits, foraging behaviour, degree of specialization), could influence plant evolution. We conclude by commenting on how macroevolutionary signals of the mutualism could relate to coevolution, highlighting the unbalanced focus on the plant side of the interaction, and advocating for the use of species-level interaction data in macroevolutionary studies.

Phylogenomic discordance is driven by wide-spread introgression and incomplete lineage sorting during rapid species diversification within rattlesnakes (Viperidae: Crotalus and Sistrurus)

Phylogenomics allows us to uncover the historical signal of evolutionary processes through time and estimate phylogenetic networks accounting for these signals. Insight from genome-wide data further allows us to pinpoint the contributions to phylogenetic signal from hybridization, introgression, and ancestral polymorphism across the genome. Here we focus on how these processes have contributed to phylogenetic discordance among rattlesnakes (genera Crotalus and Sistrurus), a group for which there are numerous conflicting phylogenetic hypotheses based on a diverse array of molecular datasets and analytical methods. We address the instability of the rattlesnake phylogeny using genomic data generated from transcriptomes sampled from nearly all known species. These genomic data, analyzed with coalescent and network-based approaches, reveal numerous instances of rapid speciation where individual gene trees conflict with the species tree. Moreover, the evolutionary history of rattlesnakes is dominated by incomplete speciation and frequent hybridization, both of which have likely influenced past interpretations of phylogeny. We present a new framework in which the evolutionary relationships of this group can only be understood in light of genome-wide data and network-based analytical methods. Our data suggest that network radiations, like seen within the rattlesnakes, can only be understood in a phylogenomic context, necessitating similar approaches in our attempts to understand evolutionary history in other rapidly radiating species.

Hybrid Mimulus flowers attract a new pollinator

Hybridization is common in flowering plants and is believed to be an important force driving adaptation and speciation. The flowers of hybrids often exhibit new trait combinations, which, theoretically, could attract new species of pollinators. In this study, we found that the hybrids between a hummingbird-pollinated species Mimulus cardinalis and a self-pollinated species Mimulus parishii attract bumblebees (Bombus impatiens), a pollinator not attracted to either of the progenitor species. This novel attraction is explained by new combinations of floral traits in hybrids, including, most importantly, petal color, in addition to nectar concentration and corolla size. To understand how petal color variation is perceived by bumblebees, we performed reflectance spectroscopy and multispectral imaging to model the flower appearance in bee vision. This analysis showed that color variation would impact the ease of detection. We also found that YUP, the genetic locus responsible for a large portion of floral color variation and previously shown to be important in bee interactions with other Mimulus species, also played an important role in this novel attraction. These results together suggest that the attraction of new pollinators to hybrid plants could be an underexplored avenue for pollinator shift and speciation.

Karyotypic stasis and swarming influenced the evolution of viral tolerance in a species-rich bat radiation

The emergence of COVID-19 and severe acute respiratory syndrome (SARS) has prioritized understanding bats' viral tolerance. Myotis bats are exceptionally species rich and have evolved viral tolerance. They also exhibit swarming, a cryptic behavior where large, multi-species assemblages gather for mating, which has been hypothesized to promote interspecific hybridization. To resolve the coevolution of genome architecture and their unusual antiviral tolerance, we undertook a phylogenomic analysis of 60 Old World Myotis genomes. We demonstrate an extensive history of introgressive hybridization that has replaced the species phylogeny across 17%-93% of the genome except for pericentromeric regions of macrochromosomes. Introgression tracts were enriched on microchromosome regions containing key antiviral pathway genes overexpressed during viral challenge experiments. Together, these results suggest that the unusual Myotis karyotype may have evolved to selectively position immune-related genes in high recombining genomic regions prone to introgression of divergent alleles, including a diversity of interleukin loci responsible for the release of pro-inflammatory cytokines.

A lethal mitonuclear incompatibility in complex I of natural hybrids

The evolution of reproductive barriers is the first step in the formation of new species and can help us understand the diversification of life on Earth. These reproductive barriers often take the form of hybrid incompatibilities, in which alleles derived from two different species no longer interact properly in hybrids1-3. Theory predicts that hybrid incompatibilities may be more likely to arise at rapidly evolving genes4-6 and that incompatibilities involving multiple genes should be common7,8, but there has been sparse empirical data to evaluate these predictions. Here we describe a mitonuclear incompatibility involving three genes whose protein products are in physical contact within respiratory complex I of naturally hybridizing swordtail fish species. Individuals homozygous for mismatched protein combinations do not complete embryonic development or die as juveniles, whereas those heterozygous for the incompatibility have reduced complex I function and unbalanced representation of parental alleles in the mitochondrial proteome. We find that the effects of different genetic interactions on survival are non-additive, highlighting subtle complexity in the genetic architecture of hybrid incompatibilities. Finally, we document the evolutionary history of the genes involved, showing signals of accelerated evolution and evidence that an incompatibility has been transferred between species via hybridization.

Genetic basis of nectar guide trichome variation between bumblebee- and self-pollinated monkeyflowers (Mimulus): role of the MIXTA-like gene GUIDELESS

Nectar guide trichomes play crucial ecological roles in bee-pollinated flowers, as they serve as footholds and guides for foraging bees to access the floral rewards. However, the genetic basis of natural variation in nectar guide trichomes among species remains poorly understood. In this study, we performed genetic analysis of nectar guide trichome variation between two closely related monkeyflower (Mimulus) species, the bumblebee-pollinated Mimulus lewisii and self-pollinated M. parishii. We demonstrate that a MIXTA-like R2R3-MYB gene, GUIDELESS, is a major contributor to the nectar guide trichome length variation between the two species. The short-haired M. parishii carries a recessive allele due to non-synonymous substitutions in a highly conserved motif among MIXTA-like MYB proteins. Furthermore, our results suggest that besides GUIDELESS, additional loci encoding repressors of trichome elongation also contribute to the transition from bumblebee-pollination to selfing. Taken together, these results suggest that during a pollination syndrome switch, changes in seemingly complex traits such as nectar guide trichomes could have a relatively simple genetic basis, involving just a few genes of large effects.

How the switch to hummingbird pollination has greatly contributed to our understanding of evolutionary processes

The evolutionary switch to hummingbird pollination exemplifies complex adaptation, requiring evolutionary change in multiple component traits. Despite this complexity, diverse lineages have converged on hummingbird-adapted flowers on a relatively short evolutionary timescale. Here, I review how features of the genetic basis of adaptation contribute to this remarkable evolutionary lability. Large-effect substitutions, large mutational targets for adaptation, adaptive introgression, and concentrated architecture all contribute to the origin and maintenance of hummingbird-adapted flowers. The genetic features of adaptation are likely shaped by the ecological and geographic context of the switch to hummingbird pollination, with implications for future evolutionary trajectories.

Transposition, duplication, and divergence of the telomerase RNA underlies the evolution of Mimulus telomeres

Telomeres are nucleoprotein complexes with a crucial role of protecting chromosome ends. It consists of simple repeat sequences and dedicated telomere-binding proteins. Because of its vital functions, components of the telomere, for example its sequence, should be under strong evolutionary constraint. But across all plants, telomere sequences display a range of variation and the evolutionary mechanism driving this diversification is largely unknown. Here, we discovered in Monkeyflower (Mimulus) the telomere sequence is even variable between species. We investigated the basis of Mimulus telomere sequence evolution by studying the long noncoding telomerase RNA (TR), which is a core component of the telomere maintenance complex and determines the telomere sequence. We conducted total RNA-based de novo transcriptomics from 16 Mimulus species and analyzed reference genomes from 6 species, and discovered Mimulus species have evolved at least three different telomere sequences: (AAACCCT)n, (AAACCCG)n, and (AAACCG)n. Unexpectedly, we discovered several species with TR duplications and the paralogs had functional consequences that could influence telomere evolution. For instance, M. lewisii had two sequence-divergent TR paralogs and synthesized a telomere with sequence heterogeneity, consisting of AAACCG and AAACCCG repeats. Evolutionary analysis of the M. lewisii TR paralogs indicated it had arisen from a transposition-mediate duplication process. Further analysis of the TR from multiple Mimulus species showed the gene had frequently transposed and inserted into new chromosomal positions during Mimulus evolution. From our results, we propose the TR transposition, duplication, and divergence model to explain the evolutionary sequence turnovers in Mimulus and potentially all plant telomeres.

Reproductive strategies and their consequences for divergence, gene flow, and genetic diversity in three taxa of Clarkia

Differences in reproductive strategies can have important implications for macro- and micro-evolutionary processes. We used a comparative approach through a population genetics lens to evaluate how three distinct reproductive strategies shape patterns of divergence among as well as gene flow and genetic diversity within three closely related taxa in the genus Clarkia. One taxon is a predominantly autonomous self-fertilizer and the other two taxa are predominantly outcrossing but vary in the primary pollinator they attract. In genotyping populations using genotyping-by-sequencing and comparing loci shared across taxa, our results suggest that differences in reproductive strategies in part promote evolutionary divergence among these closely related taxa. Contrary to expectations, we found that the selfing taxon had the highest levels of heterozygosity but a low rate of polymorphism. The high levels of fixed heterozygosity for a subset of loci suggests this pattern is driven by the presence of structural rearrangements in chromosomes common in other Clarkia taxa. In evaluating patterns within taxa, we found a complex interplay between reproductive strategy and geographic distribution. Differences in the mobility of primary pollinators did not translate to a difference in rates of genetic diversity and gene flow within taxa - a pattern likely due to one taxon having a patchier distribution and a less temporally and spatially reliable pollinator. Taken together, this work advances our understanding of the factors that shape gene flow and the distribution of genetic diversity within and among closely related taxa.

A phylogeny of the evening primrose family (Onagraceae) using a target enrichment approach with 303 nuclear loci

BACKGROUND

The evening primrose family (Onagraceae) includes 664 species (803 taxa) with a center of diversity in the Americas, especially western North America. Ongoing research in Onagraceae includes exploring striking variation in floral morphology, scent composition, and breeding system, as well as the role of these traits in driving diversity among plants and their interacting pollinators and herbivores. However, these efforts are limited by the lack of a comprehensive, well-resolved phylogeny. Previous phylogenetic studies based on a few loci strongly support the monophyly of the family and the sister relationship of the two largest tribes but fail to resolve several key relationships.

RESULTS

We used a target enrichment approach to reconstruct the phylogeny of Onagraceae using 303 highly conserved, low-copy nuclear loci. We present a phylogeny for Onagraceae with 169 individuals representing 152 taxa sampled across the family, including extensive sampling within the largest tribe, Onagreae. Deep splits within the family are strongly supported, whereas relationships among closely related genera and species are characterized by extensive conflict among individual gene trees.

CONCLUSIONS

This phylogenetic resource will augment current research projects focused throughout the family in genomics, ecology, coevolutionary dynamics, biogeography, and the evolution of characters driving diversification in the family.

Genomic mechanisms and consequences of diverse postzygotic barriers between monkeyflower species

The evolution of genomic incompatibilities causing postzygotic barriers to hybridization is a key step in species divergence. Incompatibilities take 2 general forms-structural divergence between chromosomes leading to severe hybrid sterility in F1 hybrids and epistatic interactions between genes causing reduced fitness of hybrid gametes or zygotes (Dobzhansky-Muller incompatibilities). Despite substantial recent progress in understanding the molecular mechanisms and evolutionary origins of both types of incompatibility, how each behaves across multiple generations of hybridization remains relatively unexplored. Here, we use genetic mapping in F2 and recombinant inbred line (RIL) hybrid populations between the phenotypically divergent but naturally hybridizing monkeyflowers Mimulus cardinalis and M. parishii to characterize the genetic basis of hybrid incompatibility and examine its changing effects over multiple generations of experimental hybridization. In F2s, we found severe hybrid pollen inviability (<50% reduction vs parental genotypes) and pseudolinkage caused by a reciprocal translocation between Chromosomes 6 and 7 in the parental species. RILs retained excess heterozygosity around the translocation breakpoints, which caused substantial pollen inviability when interstitial crossovers had not created compatible heterokaryotypic configurations. Strong transmission ratio distortion and interchromosomal linkage disequilibrium in both F2s and RILs identified a novel 2-locus genic incompatibility causing sex-independent gametophytic (haploid) lethality. The latter interaction eliminated 3 of the expected 9 F2 genotypic classes via F1 gamete loss without detectable effects on the pollen number or viability of F2 double heterozygotes. Along with the mapping of numerous milder incompatibilities, these key findings illuminate the complex genetics of plant hybrid breakdown and are an important step toward understanding the genomic consequences of natural hybridization in this model system.

Strong postmating reproductive isolation in Mimulus section Eunanus

Postmating reproductive isolation can help maintain species boundaries when premating barriers to reproduction are incomplete. The strength and identity of postmating reproductive barriers are highly variable among diverging species, leading to questions about their genetic basis and evolutionary drivers. These questions have been tackled in model systems but are less often addressed with broader phylogenetic resolution. In this study we analyse patterns of genetic divergence alongside direct measures of postmating reproductive barriers in an overlooked group of sympatric species within the model monkeyflower genus, Mimulus. Within this Mimulus brevipes species group, we find substantial divergence among species, including a cryptic genetic lineage. However, rampant gene discordance and ancient signals of introgression suggest a complex history of divergence. In addition, we find multiple strong postmating barriers, including postmating prezygotic isolation, hybrid seed inviability and hybrid male sterility. M. brevipes and M. fremontii have substantial but incomplete postmating isolation. For all other tested species pairs, we find essentially complete postmating isolation. Hybrid seed inviability appears linked to differences in seed size, providing a window into possible developmental mechanisms underlying this reproductive barrier. While geographic proximity and incomplete mating isolation may have allowed gene flow within this group in the distant past, strong postmating reproductive barriers today have likely played a key role in preventing ongoing introgression. By producing foundational information about reproductive isolation and genomic divergence in this understudied group, we add new diversity and phylogenetic resolution to our understanding of the mechanisms of plant speciation.

Ancient hybridization leads to the repeated evolution of red flowers across a monkeyflower radiation

The reuse of old genetic variation can promote rapid diversification in evolutionary radiations, but in most cases, the historical events underlying this divergence are not known. For example, ancient hybridization can generate new combinations of alleles that sort into descendant lineages, potentially providing the raw material to initiate divergence. In the Mimulus aurantiacus species complex, there is evidence for widespread gene flow among members of this radiation. In addition, allelic variation in the MaMyb2 gene is responsible for differences in flower color between the closely related ecotypes of subspecies puniceus, contributing to reproductive isolation by pollinators. Previous work suggested that MaMyb2 was introgressed into the red-flowered ecotype of puniceus. However, additional taxa within the radiation have independently evolved red flowers from their yellow-flowered ancestors, raising the possibility that this introgression had a more ancient origin. In this study, we used repeated tests of admixture from whole-genome sequence data across this diverse radiation to demonstrate that there has been both ancient and recurrent hybridization in this group. However, most of the signal of this ancient introgression has been removed due to selection, suggesting that widespread barriers to gene flow are in place between taxa. Yet, a roughly 30 kb region that contains the MaMyb2 gene is currently shared only among the red-flowered taxa. Patterns of admixture, sequence divergence, and extended haplotype homozygosity across this region confirm a history of ancient hybridization, where functional variants have been preserved due to positive selection in red-flowered taxa but lost in their yellow-flowered counterparts. The results of this study reveal that selection against gene flow can reduce genomic signatures of ancient hybridization, but that historical introgression can provide essential genetic variation that facilitates the repeated evolution of phenotypic traits between lineages.

Introgression between highly divergent fungal sister species

To understand how species evolve and adapt to changing environments, it is important to study gene flow and introgression due to their influence on speciation and radiation events. Here, we apply a novel experimental system for investigating these mechanisms using natural populations. The system is based on two fungal sister species with morphological and ecological similarities occurring in overlapping habitats. We examined introgression between these species by conducting whole genome sequencing of individuals from populations in North America and Europe. We assessed genome-wide nucleotide divergence and performed crossing experiments to study reproductive barriers. We further used ABBA-BABA statistics together with a network analysis to investigate introgression, and conducted demographic modelling to gain insight into divergence times and introgression events. The results revealed that the species are highly divergent and incompatible in vitro. Despite this, small regions of introgression were scattered throughout the genomes and one introgression event likely involves a ghost population (extant or extinct). This study demonstrates that introgression can be found among divergent species and that population histories can be studied without collections of all the populations involved. Moreover, the experimental system is shown to be a useful tool for research on reproductive isolation in natural populations.

The role of recombination landscape in species hybridisation and speciation

It is now well recognised that closely related species can hybridize and exchange genetic material, which may promote or oppose adaptation and speciation. In some cases, interspecific hybridisation is very common, making it surprising that species identity is preserved despite active gene exchange. The genomes of most eukaryotic species are highly heterogeneous with regard to gene density, abundance of repetitive DNA, chromatin compactisation etc, which can make certain genomic regions more prone or more resistant to introgression of genetic material from other species. Heterogeneity in local recombination rate underpins many of the observed patterns across the genome (e.g. actively recombining regions are typically gene rich and depleted for repetitive DNA) and it can strongly affect the permeability of genomic regions to interspecific introgression. The larger the region lacking recombination, the higher the chance for the presence of species incompatibility gene(s) in that region, making the entire non- or rarely recombining block impermeable to interspecific introgression. Large plant genomes tend to have highly heterogeneous recombination landscape, with recombination frequently occurring at the ends of the chromosomes and central regions lacking recombination. In this paper we review the relationship between recombination and introgression in plants and argue that large rarely recombining regions likely play a major role in preserving species identity in actively hybridising plant species.

A genomic timescale for placental mammal evolution

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.

A clue to the evolutionary history of modern East Asian flora: insights from phylogeography and diterpenoid alkaloid distribution pattern of the Spiraea japonica complex

Each subkingdom of East Asian flora (EAF) has a unique evolutionary history, but which has rarely been described based on phylogeographic studies of EAF species. The Spiraea japonica L. complex, which is widespread in East Asia (EA), has received considerable attention because of the presence of diterpenoid alkaloids (DAs). It provides a proxy for understanding the genetic diversity and DA distribution patterns of species under various environmental conditions associated with the geological background in EA. In the present study, the plastome and chloroplast/nuclear DNA of 71 populations belonging to the S. japonica complex and its congeners were sequenced, combined with DA identification, environmental analyses, and ecological niche modelling, to investigate their phylogenetic relationships, genetic and DAs distribution patterns, biogeography, and demographic dynamics. An "ampliative" S. japonica complex was put forward, comprising all species of Sect. Calospira Ser. Japonicae, of which three evolutionary units carrying their respective unique types of DAs were identified and associated with the regionalization of EAF (referring to the Hengduan Mountains, central China, and east China). Moreover, a transition belt in central China with its biogeographic significance was revealed by genetic and DA distribution patterns from the perspective of ecological adaptation. The origin and onset differentiation of the "ampliative" S. japonica complex was estimated in the early Miocene (22.01/19.44 Ma). The formation of Japanese populations (6.75 Ma) was facilitated by the land bridge, which subsequently had a fairly stable demographic history. The populations in east China have undergone a founder effect after the Last Glacial Maximum, which may have been promoted by the expansion potential of polyploidization. Overall, the in-situ origin and diversification of the "ampliative" S. japonica complex since the early Miocene is a vertical section of the formation and development of modern EAF and was shaped by the geological history of each subkingdom.

Taxon-specific, phased siRNAs underlie a speciation locus in monkeyflowers

Taxon-specific small RNA loci are widespread in eukaryotic genomes, yet their role in lineage-specific adaptation, phenotypic diversification, and speciation is poorly understood. Here, we report that a speciation locus in monkeyflowers (Mimulus), YELLOW UPPER (YUP), contains an inverted repeat region that produces small interfering RNAs (siRNAs) in a phased pattern. Although the inverted repeat is derived from a partial duplication of a protein-coding gene that is not involved in flower pigmentation, one of the siRNAs targets and represses a master regulator of floral carotenoid pigmentation. YUP emerged with two protein-coding genes that control other aspects of flower coloration as a "superlocus" in a subclade of Mimulus and has contributed to subsequent phenotypic diversification and pollinator-mediated speciation in the descendant species.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

To stripe or not to stripe: the origin of a novel foliar pigmentation pattern in monkeyflowers (Mimulus)

The origin of phenotypic novelty is one of the most challenging problems in evolutionary biology. Although genetic regulatory network rewiring or co-option has been widely recognised as a major contributor, in most cases how such genetic rewiring/co-option happens is completely unknown. We have studied a novel foliar pigmentation pattern that evolved recently in the monkeyflower species Mimulus verbenaceus. Through genome-wide association tests using wild-collected samples, experimental crosses of laboratory inbred lines, gene expression analyses, and functional assays, we identified an anthocyanin-activating R2R3-MYB gene, STRIPY, as the causal gene triggering the emergence of the discrete, mediolateral anthocyanin stripe in the M. verbenaceus leaf. Chemical mutagenesis revealed the existence of upstream activators and repressors that form a 'hidden' prepattern along the leaf proximodistal axis, potentiating the unique expression pattern of STRIPY. Population genomics analyses did not reveal signatures of positive selection, indicating that nonadaptive processes may be responsible for the establishment of this novel trait in the wild. This study demonstrates that the origin of phenotypic novelty requires at least two separate phases, potentiation and actualisation. The foliar stripe pattern of M. verbenaceus provides an excellent platform to probe the molecular details of both processes in future studies.

Evolutionary convergence on hummingbird pollination in Neotropical Costus provides insight into the causes of pollinator shifts

The evolution of hummingbird pollination is common across angiosperms throughout the Americas, presenting an opportunity to examine convergence in both traits and environments to better understand how complex phenotypes arise. Here we examine independent shifts from bee to hummingbird pollination in the Neotropical spiral gingers (Costus) and address common explanations for the prevalence of transitions from bee to hummingbird pollination. We use floral traits of species with observed pollinators to predict pollinators of unobserved species and reconstruct ancestral pollination states on a well-resolved phylogeny. We examine whether independent transitions evolve towards the same phenotypic optimum and whether shifts to hummingbird pollination correlate with elevation or climate. Traits predicting hummingbird pollination include small flower size, brightly colored floral bracts and the absence of nectar guides. We find many shifts to hummingbird pollination and no reversals, a single shared phenotypic optimum across hummingbird flowers, and no association between pollination and elevation or climate. Evolutionary shifts to hummingbird pollination in Costus are highly convergent and directional, involve a surprising set of traits when compared with other plants with analogous transitions and refute the generality of several common explanations for the prevalence of transitions from bee to hummingbird pollination.

The strength of reproductive isolating barriers in seed plants: Insights from studies quantifying premating and postmating reproductive barriers over the past 15 years

Speciation is driven by the evolution of reproductive isolating barriers that reduce, and ultimately prevent, substantial gene flow between lineages. Despite its central role in evolutionary biology, the process can be difficult to study because it proceeds differently among groups and may occur over long timescales. Due to this complexity, we typically rely on generalizations of empirical data to describe and understand the process. Previous reviews of reproductive isolation (RI) in flowering plants have suggested that prezygotic or extrinsic barriers generally have a stronger effect on reducing gene flow compared to postzygotic or intrinsic barriers. Past conclusions have rested on relatively few empirical estimates of RI; however, RI data have become increasingly abundant over the past 15 years. We analyzed data from recent studies quantifying multiple pre- and postmating barriers in plants and compared the strengths of isolating barriers across 89 taxa pairs using standardized RI metrics. Individual prezygotic barriers were on average stronger than individual postzygotic barriers, and the total strength of prezygotic RI was approximately twice that of postzygotic RI. These findings corroborate that ecological divergence and extrinsic factors, as opposed to solely the accumulation of genetic incompatibilities, are important to speciation and the maintenance of species boundaries in plants. Despite an emphasis in the literature on asymmetric postmating and postzygotic RI, we found that prezygotic barriers acted equally asymmetrically. Overall, substantial variability in the strengths of 12 isolating barriers highlights the great diversity of mechanisms that contribute to plant diversification.

Lost in translation: Molecular basis of reduced flower coloration in a self-pollinated monkeyflower (Mimulus) species

Phenotypic evolution is usually attributed to changes in protein function or gene transcription. In principle, mutations that affect protein abundance through enhancing or attenuating protein translation also could be an important source for phenotypic evolution. However, these types of mutations remain largely unexplored in the studies of phenotypic variation in nature. Through fine-scale genetic mapping and functional interrogation, we identify a single nucleotide substitution in an anthocyanin-activating R2R3-MYB gene causing flower color variation between a pair of closely related monkeyflower (Mimulus) species, the hummingbird-pollinated Mimulus cardinalis, and self-pollinated Mimulus parishii. This causal mutation is located in the 5' untranslated region and generates an upstream ATG start codon, leading to attenuated protein translation and reduced flower coloration in the self-pollinated species. Together, our results provide empirical support for the role of mutations affecting protein translation, as opposed to protein function or transcript level, in natural phenotypic variation.

Tree House Explorer: A Novel Genome Browser for Phylogenomics

Tree House Explorer (THEx) is a genome browser that integrates phylogenomic data and genomic annotations into a single interactive platform for combined analysis. THEx allows users to visualize genome-wide variation in evolutionary histories and genetic divergence on a chromosome-by-chromosome basis, with continuous sliding window comparisons to gene annotations (GFF/GTF), recombination rates, and other user-specified, highly customizable feature annotations. THEx provides a new platform for interactive phylogenomic data visualization to analyze and interpret the diverse evolutionary histories woven throughout genomes. Hosted on Conda, THEx integrates seamlessly into new or pre-existing workflows.

Phylogenomic analyses in Phrymaceae reveal extensive gene tree discordance in relationships among major clades

PREMISE

Phylogenomic datasets using genomes and transcriptomes provide rich opportunities beyond resolving bifurcating phylogenetic relationships. Monkeyflower (Phrymaceae) is a model system for evolutionary ecology. However, it lacks a well-supported phylogeny as a basis for a stable taxonomy and for macroevolutionary comparisons.

METHODS

We sampled 24 genomes and transcriptomes in Phrymaceae and closely related families, including eight newly sequenced transcriptomes. We reconstructed the phylogeny using IQ-TREE and ASTRAL, evaluated gene tree discordance using PhyParts, Quartet Sampling, and a cloudogram, and carried out reticulation analyses using PhyloNet and HyDe. We searched for whole genome duplication (WGD) events using chromosome numbers, synonymous distances, and gene duplication events as evidence.

RESULTS

Most gene trees support the monophyly of Phrymaceae and each of its tribes. Most gene trees also support tribe Mimuleae being sister to Phrymeae + Diplaceae + Leucocarpeae, with extensive gene tree discordance among the latter three. Despite the discordance, the monophyly of Mimulus s.l. is rejected, and no individual reticulation event among the Phrymaceae tribes is well-supported. Reticulation likely occurred among Erythranthe bicolor and closely related species. No ancient WGD was detected in Phrymaceae. Instead, small-scale duplications are among potential drivers of macroevolutionary diversification of Phrymaceae.

CONCLUSIONS

We show that analysis of reticulate evolution is sensitive to taxon sampling and methods used. We also demonstrate that phylogenomic datasets using genomes and transcriptomes present rich opportunities to investigate gene family evolution and genome duplication events involved in lineage diversification and adaptation.

Phylogenomic analysis does not support a classic but controversial hypothesis of progenitor-derivative origins for the serpentine endemic Clarkia franciscana

Budding speciation involves isolation of marginal populations at the periphery of a species range and is thought to be a prominent mode of speciation in organisms with low dispersal and/or strong local adaptation among populations. Budding speciation is typically evidenced by abutting, asymmetric ranges of ecologically divergent sister species and low genetic diversity in putative budded species. Yet these indirect patterns may be unreliable, instead caused by postspeciation processes such as range or demographic shifts. Nested phylogenetic relationships provide the most conclusive evidence of budding speciation. A putative case of budding speciation in the serpentine endemic Clarkia franciscana and two closely related widespread congeners was studied by Harlan Lewis, Peter Raven, Leslie Gottlieb, and others over a 20-year period, yet the origin of C. franciscana remains controversial. Here, we reinvestigate this system with phylogenomic analyses to determine whether C. franciscana is a recently derived budded species, phylogenetically nested within one of the other two putative progenitor species. In contrast to the hypothesized pattern of relatedness among the three Clarkia species, we find no evidence for recent budding speciation. Instead, the data suggest the three species diverged simultaneously. We urge caution in using contemporary range patterns to infer geographic modes of speciation.

Hybridization alters the shape of the genotypic fitness landscape, increasing access to novel fitness peaks during adaptive radiation

Estimating the complex relationship between fitness and genotype or phenotype (i.e. the adaptive landscape) is one of the central goals of evolutionary biology. However, adaptive walks connecting genotypes to organismal fitness, speciation, and novel ecological niches are still poorly understood and processes for surmounting fitness valleys remain controversial. One outstanding system for addressing these connections is a recent adaptive radiation of ecologically and morphologically novel pupfishes (a generalist, molluscivore, and scale-eater) endemic to San Salvador Island, Bahamas. We leveraged whole-genome sequencing of 139 hybrids from two independent field fitness experiments to identify the genomic basis of fitness, estimate genotypic fitness networks, and measure the accessibility of adaptive walks on the fitness landscape. We identified 132 single nucleotide polymorphisms (SNPs) that were significantly associated with fitness in field enclosures. Six out of the 13 regions most strongly associated with fitness contained differentially expressed genes and fixed SNPs between trophic specialists; one gene (mettl21e) was also misexpressed in lab-reared hybrids, suggesting a potential intrinsic genetic incompatibility. We then constructed genotypic fitness networks from adaptive alleles and show that scale-eating specialists are the most isolated of the three species on these networks. Intriguingly, introgressed and de novo variants reduced fitness landscape ruggedness as compared to standing variation, increasing the accessibility of genotypic fitness paths from generalist to specialists. Our results suggest that adaptive introgression and de novo mutations alter the shape of the fitness landscape, providing key connections in adaptive walks circumventing fitness valleys and triggering the evolution of novelty during adaptive radiation.

Spatial variation in high temperature-regulated gene expression predicts evolution of plasticity with climate change in the scarlet monkeyflower

A major way that organisms can adapt to changing environmental conditions is by evolving increased or decreased phenotypic plasticity. In the face of current global warming, more attention is being paid to the role of plasticity in maintaining fitness as abiotic conditions change over time. However, given that temporal data can be challenging to acquire, a major question is whether evolution in plasticity across space can predict adaptive plasticity across time. In growth chambers simulating two thermal regimes, we generated transcriptome data for western North American scarlet monkeyflowers (Mimulus cardinalis) collected from different latitudes and years (2010 and 2017) to test hypotheses about how plasticity in gene expression is responding to increases in temperature, and if this pattern is consistent across time and space. Supporting the genetic compensation hypothesis, individuals whose progenitors were collected from the warmer-origin northern 2017 descendant cohort showed lower thermal plasticity in gene expression than their cooler-origin northern 2010 ancestors. This was largely due to a change in response at the warmer (40°C) rather than cooler (20°C) treatment. A similar pattern of reduced plasticity, largely due to a change in response at 40°C, was also found for the cooler-origin northern versus the warmer-origin southern population from 2017. Our results demonstrate that reduced phenotypic plasticity can evolve with warming and that spatial and temporal changes in plasticity predict one another.

Predictors of genomic differentiation within a hybrid taxon

Hybridization is increasingly recognized as an important evolutionary force. Novel genetic methods now enable us to address how the genomes of parental species are combined in hybrid lineages. However, we still do not know the relative importance of admixed proportions, genome architecture and local selection in shaping hybrid genomes. Here, we take advantage of the genetically divergent island populations of Italian sparrow on Crete, Corsica and Sicily to investigate the predictors of genomic variation within a hybrid taxon. We test if differentiation is affected by recombination rate, selection, or variation in ancestry proportions. We find that the relationship between recombination rate and differentiation is less pronounced within hybrid lineages than between the parent species, as expected if purging of minor parent ancestry in low recombination regions reduces the variation available for differentiation. In addition, we find that differentiation between islands is correlated with differences in signatures of selection in two out of three comparisons. Signatures of selection within islands are correlated across all islands, suggesting that shared selection may mould genomic differentiation. The best predictor of strong differentiation within islands is the degree of differentiation from house sparrow, and hence loci with Spanish sparrow ancestry may vary more freely. Jointly, this suggests that constraints and selection interact in shaping the genomic landscape of differentiation in this hybrid species.

Bidirectional introgression between Mus musculus domesticus and Mus spretus

Introgressed variants from other species can be an important source of genetic variation because they may arise rapidly, can include multiple mutations on a single haplotype, and have often been pretested by selection in the species of origin. Although introgressed alleles are generally deleterious, several studies have reported introgression as the source of adaptive alleles—including the rodenticide-resistant variant of Vkorc1 that introgressed from Mus spretus into European populations of Mus musculus domesticus. Here, we conducted bidirectional genome scans to characterize introgressed regions into one wild population of M. spretus from Spain and three wild populations of M. m. domesticus from France, Germany, and Iran. Despite the fact that these species show considerable intrinsic postzygotic reproductive isolation, introgression was observed in all individuals, including in the M. musculus reference genome (GRCm38). Mus spretus individuals had a greater proportion of introgression compared with M. m. domesticus, and within M. m. domesticus, the proportion of introgression decreased with geographic distance from the area of sympatry. Introgression was observed on all autosomes for both species, but not on the X-chromosome in M. m. domesticus, consistent with known X-linked hybrid sterility and inviability genes that have been mapped to the M. spretus X-chromosome. Tract lengths were generally short with a few outliers of up to 2.7 Mb. Interestingly, the longest introgressed tracts were in olfactory receptor regions, and introgressed tracts were significantly enriched for olfactory receptor genes in both species, suggesting that introgression may be a source of functional novelty even between species with high barriers to gene flow.

Demographic history and natural selection shape patterns of deleterious mutation load and barriers to introgression across Populus genome

Hybridization and resulting introgression are important processes shaping the tree of life and appear to be far more common than previously thought. However, how the genome evolution was shaped by various genetic and evolutionary forces after hybridization remains unresolved. Here we used whole-genome resequencing data of 227 individuals from multiple widespread Populus species to characterize their contemporary patterns of hybridization and to quantify genomic signatures of past introgression. We observe a high frequency of contemporary hybridization and confirm that multiple previously ambiguous species are in fact F₁ hybrids. Seven species were identified, which experienced different demographic histories that resulted in strikingly varied efficacy of selection and burdens of deleterious mutations. Frequent past introgression has been found to be a pervasive feature throughout the speciation of these Populus species. The retained introgressed regions, more generally, tend to contain reduced genetic load and to be located in regions of high recombination. We also find that in pairs of species with substantial differences in effective population size, introgressed regions are inferred to have undergone selective sweeps at greater than expected frequencies in the species with lower effective population size, suggesting that introgression likely have higher potential to provide beneficial variation for species with small populations. Our results, therefore, illustrate that demography and recombination have interplayed with both positive and negative selection in determining the genomic evolution after hybridization.

Highly Replicated Evolution of Parapatric Ecotypes

Parallel evolution of ecotypes occurs when selection independently drives the evolution of similar traits across similar environments. The multiple origins of ecotypes are often inferred based on a phylogeny that clusters populations according to geographic location and not by the environment they occupy. However, the use of phylogenies to infer parallel evolution in closely related populations is problematic because gene flow and incomplete lineage sorting can uncouple the genetic structure at neutral markers from the colonization history of populations. Here, we demonstrate multiple origins within ecotypes of an Australian wildflower, Senecio lautus. We observed strong genetic structure as well as phylogenetic clustering by geography and show that this is unlikely due to gene flow between parapatric ecotypes, which was surprisingly low. We further confirm this analytically by demonstrating that phylogenetic distortion due to gene flow often requires higher levels of migration than those observed in S. lautus. Our results imply that selection can repeatedly create similar phenotypes despite the perceived homogenizing effects of gene flow.

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize's global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low FST between mexicana ancestry tracts sampled from geographically distant maize populations. Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.

Divergence in the Aquilegia ecalcarata complex is correlated with geography and climate oscillations: Evidence from plastid genome data

Quaternary climate oscillations and geographical heterogeneity play important roles in determining species and genetic diversity distribution patterns, but how these factors affect the migration and differentiation of East Asian plants species at the population level remains poorly understood. The Aquilegia ecalcarata complex, a group that originated in the Late Tertiary and is widely distributed throughout East Asia, displays high genetic variation that is suitable for studying elaborate phylogeographic patterns and demographic history related to the impact of Quaternary climate and geography. We used plastid genome data from 322 individuals in 60 populations of the A. ecalcarata complex to thoroughly explore the impact of Quaternary climate oscillations and geography on the phylogeographic patterns and demographic history of the A. ecalcarata complex through a series of phylogenetic, divergence time estimation, and demographic history analyses. The dry, cold climate and frequent climate oscillations that occurred during the early Pleistocene and the Mid-Pleistocene transition led to the differentiation of the A. ecalcarata complex, which was isolated in various areas. Geographically, the A. ecalcarata complex can be divided into Eastern and Western Clades and five subclades, which conform to the divergence of the East Asian flora. Our results clearly show the impact of Quaternary climate and geography on evolutionary history at the population level. These findings promote the understanding of the relationship between plant genetic differentiation and climate and geographical factors of East Asia at the population level.

Extensive introgression and mosaic genomes of Mediterranean endemic lizards

The Mediterranean basin is a hotspot of biodiversity, fuelled by climatic oscillation and geological change over the past 20 million years. Wall lizards of the genus Podarcis are among the most abundant, diverse, and conspicuous Mediterranean fauna. Here, we unravel the remarkably entangled evolutionary history of wall lizards by sequencing genomes of 34 major lineages covering 26 species. We demonstrate an early (>11 MYA) separation into two clades centred on the Iberian and Balkan Peninsulas, and two clades of Mediterranean island endemics. Diversification within these clades was pronounced between 6.5-4.0 MYA, a period spanning the Messinian Salinity Crisis, during which the Mediterranean Sea nearly dried up before rapidly refilling. However, genetic exchange between lineages has been a pervasive feature throughout the entire history of wall lizards. This has resulted in a highly reticulated pattern of evolution across the group, characterised by mosaic genomes with major contributions from two or more parental taxa. These hybrid lineages gave rise to several of the extant species that are endemic to Mediterranean islands. The mosaic genomes of island endemics may have promoted their extraordinary adaptability and striking diversity in body size, shape and colouration, which have puzzled biologists for centuries.

Quantitative trait locus mapping reveals an independent genetic basis for joint divergence in leaf function, life-history, and floral traits between scarlet monkeyflower (Mimulus cardinalis) populations

PREMISE

Across taxa, vegetative and floral traits that vary along a fast-slow life-history axis are often correlated with leaf functional traits arrayed along the leaf economics spectrum, suggesting a constrained set of adaptive trait combinations. Such broad-scale convergence may arise from genetic constraints imposed by pleiotropy (or tight linkage) within species, or from natural selection alone. Understanding the genetic basis of trait syndromes and their components is key to distinguishing these alternatives and predicting evolution in novel environments.

METHODS

We used a line-cross approach and quantitative trait locus (QTL) mapping to characterize the genetic basis of twenty leaf functional/physiological, life history, and floral traits in hybrids between annualized and perennial populations of scarlet monkeyflower (Mimulus cardinalis).

RESULTS

We mapped both single and multi-trait QTLs for life history, leaf function and reproductive traits, but found no evidence of genetic co-ordination across categories. A major QTL for three leaf functional traits (thickness, photosynthetic rate, and stomatal resistance) suggests that a simple shift in leaf anatomy may be key to adaptation to seasonally dry habitats.

CONCLUSIONS

Our results suggest that the co-ordination of resource-acquisitive leaf physiological traits with a fast life-history and more selfing mating system results from environmental selection rather than functional or genetic constraint. Independent assortment of distinct trait modules, as well as a simple genetic basis to leaf physiological traits associated with drought escape, may facilitate adaptation to changing climates.