Links between floral morphology and floral visitors along an elevational gradient in aPenstemonhybrid zone

Whole-genome sequence and annotation of Penstemon davidsonii

Penstemon is the most speciose flowering plant genus endemic to North America. Penstemon species' diverse morphology and adaptation to various environments have made them a valuable model system for studying evolution. Here, we report the first full reference genome assembly and annotation for Penstemon davidsonii. Using PacBio long-read sequencing and Hi-C scaffolding technology, we constructed a de novo reference genome of 437,568,744 bases, with a contig N50 of 40 Mb and L50 of 5. The annotation includes 18,199 gene models, and both the genome and transcriptome assembly contain over 95% complete eudicot BUSCOs. This genome assembly will serve as a valuable reference for studying the evolutionary history and genetic diversity of the Penstemon genus.

Ecological Diversification in an Adaptive Radiation of Plants: The Role of De Novo Mutation and Introgression

Adaptive radiations are characterized by rapid ecological diversification and speciation events, leading to fuzzy species boundaries between ecologically differentiated species. Adaptive radiations are therefore key systems for understanding how species are formed and maintained, including the role of de novo mutations versus preexisting variation in ecological adaptation and the genome-wide consequences of hybridization events. For example, adaptive introgression, where beneficial alleles are transferred between lineages through hybridization, may fuel diversification in adaptive radiations and facilitate adaptation to new environments. In this study, we employed whole-genome resequencing data to investigate the evolutionary origin of hummingbird-pollinated flowers and to characterize genome-wide patterns of phylogenetic discordance and introgression in Penstemon subgenus Dasanthera, a small and diverse adaptive radiation of plants. We found that magenta hummingbird-adapted flowers have apparently evolved twice from ancestral blue-violet bee-pollinated flowers within this radiation. These shifts in flower color are accompanied by a variety of inactivating mutations to a key anthocyanin pathway enzyme, suggesting that independent de novo loss-of-function mutations underlie the parallel evolution of this trait. Although patterns of introgression and phylogenetic discordance were heterogenous across the genome, a strong effect of gene density suggests that, in general, natural selection opposes introgression and maintains genetic differentiation in gene-rich genomic regions. Our results highlight the importance of both de novo mutation and introgression as sources of evolutionary change and indicate a role for de novo mutation in driving parallel evolution in adaptive radiations.

Ecological diversification in an adaptive radiation of plants: the role of de novo mutation and introgression

Adaptive radiations are characterized by rapid ecological diversification and speciation events, leading to fuzzy species boundaries between ecologically differentiated species. Adaptive radiations are therefore key systems for understanding how species are formed and maintained, including the role of de novo mutations vs. pre-existing variation in ecological adaptation and the genome-wide consequences of hybridization events. For example, adaptive introgression, where beneficial alleles are transferred between lineages through hybridization, may fuel diversification in adaptive radiations and facilitate adaptation to new environments. In this study, we employed whole-genome resequencing data to investigate the evolutionary origin of hummingbird-pollinated flowers and to characterize genome-wide patterns of phylogenetic discordance and introgression in Penstemon subgenus Dasanthera, a small and diverse adaptive radiation of plants. We found that magenta hummingbird-adapted flowers have apparently evolved twice from ancestral blue-violet bee-pollinated flowers within this radiation. These shifts in flower color are accompanied by a variety of inactivating mutations to a key anthocyanin pathway enzyme, suggesting that independent de novo loss-of-function mutations underlie parallel evolution of this trait. Although patterns of introgression and phylogenetic discordance were heterogenous across the genome, a strong effect of gene density suggests that, in general, natural selection opposes introgression and maintains genetic differentiation in gene-rich genomic regions. Our results highlight the importance of both de novo mutation and introgression as sources of evolutionary change and indicate a role for de novo mutation in driving parallel evolution in adaptive radiations.

Testing hypotheses of hybrid taxon formation in the shrubby beardtongues (Penstemon subgenus Dasanthera)

PREMISE

Hybridization is increasingly being identified in the genomes of species across the tree of life, leading to a general recognition that hybridization plays an important role in the generation of species diversity. While hybridization may increase species diversity directly via the formation of new taxa through hybrid speciation, it may also act indirectly via the exchange of phenotypic and genetic variance between species, which may in turn stimulate future speciation events.

METHODS

Using high-throughput sequence data, we resolved phylogenetic relationships and investigated the role of hybridization as a diversification mechanism in the shrubby beardtongues (Penstemon subgenus Dasanthera), a group of North American wildflowers that has undergone a recent and rapid adaptive radiation. Specifically, we tested four hypotheses of hybrid taxon formation resulting from hybridization between P. davidsonii and P. fruticosus.

RESULTS

Species tree inference supports the monophyly of subgenus Dasanthera and elucidates relationships between taxa distributed in the Cascades and Sierra Nevada Mountains. Results also provide evidence of gene flow between P. davidsonii and P. fruticosus and support at least one hybrid origin hypothesis (P. davidsonii var. menziesii) in a region of contemporary distributional overlap. Hybridization may have also been facilitated by historical overlap in geographic distribution caused by species' responses to climatic changes during the Pleistocene.

CONCLUSIONS

Our results support a history of hybridization between focal taxa in a rapidly radiating clade of plants and more broadly contribute to our growing understanding of the role of hybridization as a diversification mechanism in plants.

Floral scent divergence across an elevational hybrid zone with varying pollinators

Divergence in floral traits attractive to different pollinators can promote reproductive isolation in related species. When isolation is incomplete, hybridization may occur, which offers the opportunity to explore mechanisms underlying reproductive isolation. Recent work suggests that divergence in floral scent may frequently contribute to reproductive barriers, although such divergence has seldom been examined in species with generalized pollination. Here, we used two closely related Penstemon species, P. newberryi and P. davidsonii, and their natural hybrids from an elevational gradient with pollinator communities that are predicted to vary in their reliance on floral scent (i.e., primarily hummingbirds at low elevation vs. bees at high elevation). The species vary in a suite of floral traits, but scent is uncharacterized. To address whether scent varies along elevation and potentially contributes to reproductive isolation, we genetically characterized individuals collected at field and identified whether they were parental species or hybrids. We then characterized scent amount and composition. Although the parental species had similar total emissions, some scent characteristics (i.e., scent composition, aromatic emission) diverged between them and may contribute to their isolation. However, the species emitted similar compound sets which could explain hybridization in the contact area. Hybrids were similar to the parents for most scent traits, suggesting that their floral scent would not provide a strong barrier to backcrossing. Our study suggests floral scent may be a trait contributing to species boundaries even in plants with generalized pollination, and reinforces the idea that evolutionary pollinator transitions may involve changes in multiple floral traits.

Adaptation to lower latitudes and lower elevations precedes the evolution of hummingbird pollination in western North American Penstemon

PREMISE

A switch in pollinator can occur when a plant lineage enters a new habitat where the ancestral pollinator is less common, and a novel pollinator is more common. Because pollinator communities vary according to environmental tolerances and availability of resources, there may be consistent associations between pollination mode and specific regions and habitats. Such associations can be studied in lineages that have experienced multiple pollinator transitions, representing evolutionary replicates.

METHODS

Our study focused on a large clade of Penstemon wildflower species in western North America, which has repeatedly evolved hummingbird-adapted flowers from ancestral bee-adapted flowers. For each species, we estimated geographic ranges from occurrence data and inferred environmental niches from climate, topographical, and soil data. Using a phylogenetic comparative approach, we investigated whether hummingbird-adapted species occupy distinct geographic regions or habitats relative to bee-adapted species.

RESULTS

Hummingbird-adapted species occur at lower latitudes and lower elevations than bee-adapted species, resulting in a difference in their environmental niche. Bee-adapted species sister to hummingbird-adapted species are also found in relatively low elevations and latitudes, similar to their hummingbird-adapted sister species, suggesting ecogeographic shifts precede pollinator divergence. Sister species pairs-regardless of whether they differ in pollinator-show relatively little geographic range overlap.

CONCLUSIONS

Adaptation to a novel pollinator may often occur in geographic and ecological isolation from ancestral populations. The ability of a given lineage to adapt to novel pollinators may critically depend on its ability to colonize regions and habitats associated with novel pollinator communities.

Phylogeographic analysis of shrubby beardtongues reveals range expansions during the Last Glacial Maximum and implicates the Klamath Mountains as a hotspot for hybridization

Quaternary glacial cycles often altered species' geographic distributions, which in turn altered the geographic structure of species' genetic diversity. In many cases, glacial expansion forced species in temperate climates to contract their ranges and reside in small pockets of suitable habitat (refugia), where they were likely to interact closely with other species, setting the stage for potential gene exchange. These introgression events, in turn, would have degraded species boundaries, making the inference of phylogenetic relationships challenging. Using high-throughput sequence data, we employed a combination of species distribution models and hybridization tests to assess the effect of glaciation on the geographic distributions, phylogenetic relationships, and patterns of gene flow of five species of Penstemon subgenus Dasanthera, long-lived shrubby angiosperms distributed throughout the Pacific Northwest of North America. Surprisingly, we found that rather than reducing their ranges to small refugia, most Penstemon subgenus Dasanthera species experienced increased suitable habitat during the Last Glacial Maximum relative to the present day. We also found substantial evidence for gene exchange between species, with the bulk of introgression events occurring in or near the Klamath Mountains of southwestern Oregon and northwestern California. Subsequently, our phylogenetic inference reveals blurred taxonomic boundaries in the Klamath Mountains, where introgression is most prevalent. Our results question the classical paradigm of temperate species' responses to glaciation and highlight the importance of contextualizing phylogenetic inference with species' histories of introgression.

Pollinator divergence and pollination isolation between hybrids with different floral color and morphology in two sympatric Penstemon species

Differential visitation of pollinators due to divergent floral traits can lead to reproductive isolation via assortative pollen flow, which may ultimately be a driving force in plant speciation, particularly in areas of overlap. We evaluate the effects of pollinator behavioral responses to variation of intraspecific floral color and nectar rewards, on reproductive isolation between two hybrid flower color morphs (fuchsia and blue) and their parental species Penstemon roseus and P. gentianoides with a mixed-pollination system. We show that pollinators (bumblebees and hummingbirds) exhibit different behavioral responses to fuchsia and blue morphs, which could result from differential attraction or deterrence. In addition to differences in color (spectral reflectance), we found that plants with fuchsia flowers produced more and larger flowers, produced more nectar and were more visited by pollinators than those with blue flowers. These differences influenced the foraging behavior and effectiveness as pollinators of both bumblebees and hummingbirds, which contributed to reproductive isolation between the two hybrid flower color morphs and parental species. This study demonstrates how differentiation of pollination traits promotes the formation of hybrid zones leading to pollinator shifts and reproductive isolation. While phenotypic traits of fuchsia and red flowers might encourage more efficient hummingbird pollination in a mixed-pollination system, the costs of bumblebee pollination on plant reproduction could be the drivers for the repeated shifts from bumblebee- to hummingbird-mediated pollination.

Explaining Intricate Morphometric Variability with Environmental Predictors: The Case of Globularia cordifolia Species Complex

Globularia is a genus of small evergreen and perennial shrubs that are widespread in Europe. Globularia section Empetron includes a group of three species, G. cordifolia, G. meridionalis, and G. neapolitana, that have been taxonomically disputed for more than 150 years. Many morphological features have been proposed to discriminate these species. Nevertheless, evidence from both past and recent literature suggest that these differences among species are not consistent. In order to shed new light in this long-disputed group, we investigated 10 populations of the G. cordifolia species complex with both classical and geometric morphometrics and used environmental predictors in multivariate regression to explain patterns of variation. Our results showed that bract area and calyx teeth length are correlated with solar radiation and annual precipitation, whereas leaf dry mass per unit area can be explained by temperature seasonality. Leaf shape can be explained by temperature seasonality as well, although with a lower amount of explanatory and predictive power. Despite a comparatively low sample size in terms of populations, our results were based on a large number of individuals and were supported by a robust statistical approach. We can conclude that differences among the three species of Globularia could be related to the combined effects of several ecological variables and might not have taxonomical value. Our novel approach provided an ecological interpretation on a species complex that makes up a continuum of forms within the environmental framework of the Mediterranean basin.

Selection by pollinators on floral traits in generalized Trollius ranunculoides (Ranunculaceae) along altitudinal gradients

Abundance and visitation of pollinator assemblages tend to decrease with altitude, leading to an increase in pollen limitation. Thus increased competition for pollinators may generate stronger selection on attractive traits of flowers at high elevations and cause floral adaptive evolution. Few studies have related geographically variable selection from pollinators and intraspecific floral differentiation. We investigated the variation of Trollius ranunculoides flowers and its pollinators along an altitudinal gradient on the eastern Qinghai-Tibet Plateau, and measured phenotypic selection by pollinators on floral traits across populations. The results showed significant decline of visitation rate of bees along altitudinal gradients, while flies was unchanged. When fitness is estimated by the visitation rate rather than the seed number per plant, phenotypic selection on the sepal length and width shows a significant correlation between the selection strength and the altitude, with stronger selection at higher altitudes. However, significant decreases in the sepal length and width of T. ranunculoides along the altitudinal gradient did not correspond to stronger selection of pollinators. In contrast to the pollinator visitation, mean annual precipitation negatively affected the sepal length and width, and contributed more to geographical variation in measured floral traits than the visitation rate of pollinators. Therefore, the sepal size may have been influenced by conflicting selection pressures from biotic and abiotic selective agents. This study supports the hypothesis that lower pollinator availability at high altitude can intensify selection on flower attractive traits, but abiotic selection is preventing a response to selection from pollinators.

Altitudinal gradients, plant hybrid zones and evolutionary novelty

Altitudinal gradients are characterized by steep changes of the physical and biotic environment that present challenges to plant adaptation throughout large parts of the world. Hybrid zones may form where related species inhabit different neighbouring altitudes and can facilitate interspecific gene flow and potentially the breakdown of species barriers. Studies of such hybrid zones can reveal much about the genetic basis of adaptation to environmental differences stemming from changes in altitude and the maintenance of species divergence in the face of gene flow. Furthermore, owing to recombination and transgressive effects, such hybrid zones can be sources of evolutionary novelty. We document plant hybrid zones associated with altitudinal gradients and emphasize similarities and differences in their structure. We then focus on recent studies of a hybrid zone between two Senecio species that occur at high and low altitude on Mount Etna, Sicily, showing how adaptation to local environments and intrinsic selection against hybrids act to maintain it. Finally, we consider the potential of altitudinal hybrid zones for generating evolutionary novelty through adaptive introgression and hybrid speciation. Examples of homoploid hybrid species of Senecio and Pinus that originated from altitudinal hybrid zones are discussed.

Segregation of male-sterility alleles across a species boundary

Hybrid zones may serve as bridges permitting gene flow between species, including alleles influencing the evolution of breeding systems. Using greenhouse crosses, we assessed the likelihood that a hybrid zone could serve as a conduit for transfer of nuclear male-sterility alleles between a gynodioecious species and a hermaphroditic species with very rare females in some populations. Segregation patterns in progeny of crosses between rare females of hermaphroditic Schiedea menziesii and hermaphroditic plants of gynodioecious Schiedea salicaria heterozygous at the male-sterility locus, and between female S. salicaria and hermaphroditic plants from the hybrid zone, were used to determine whether male-sterility was controlled at the same locus in the parental species and the hybrid zone. Segregations of females and hermaphrodites in approximately equal ratios from many of the crosses indicate that the same nuclear male-sterility allele occurs in the parent species and the hybrid zone. These rare male-sterility alleles in S. menziesii may result from gene flow from S. salicaria through the hybrid zone, presumably facilitated by wind pollination in S. salicaria. Alternatively, rare male-sterility alleles might result from a reversal from gynodioecy to hermaphroditism in S. menziesii, or possibly de novo evolution of male sterility. Phylogenetic analysis indicates that some species of Schiedea have probably evolved separate sexes independently, but not in the lineage containing S. salicaria and S. menziesii. High levels of selfing and expression of strong inbreeding depression in S. menziesii, which together should favour females in populations, argue against a reversal from gynodioecy to hermaphroditism in S. menziesii.

Phylogenetic trait-based analyses of ecological networks

Ecological networks of two interacting guilds of species, such as flowering plants and pollinators, are common in nature, and studying their structure can yield insights into their resilience to environmental disturbances. Here we develop analytical methods for exploring the strengths of interactions within bipartite networks consisting of two guilds of phylogenetically related species. We then apply these methods to investigate the resilience of a plant-pollinator community to anticipated climate change. The methods allow the statistical assessment of, for example, whether closely related pollinators are more likely to visit plants with similar relative frequencies, and whether closely related pollinators tend to visit closely related plants. The methods can also incorporate trait information, allowing us to identify which plant traits are likely responsible for attracting different pollinators. These questions are important for our study of 14 prairie plants and their 22 insect pollinators. Over the last 70 years, six of the plants have advanced their flowering, while eight have not. When we experimentally forced earlier flowering times, five of the six advanced-flowering species experienced higher pollinator visitation rates, whereas only one of the eight other species had more visits; this network thus appears resilient to climate change, because those species with advanced flowering have ample pollinators earlier in the season. Using the methods developed here, we show that advanced-flowering plants did not have a distinct pollinator community from the other eight species. Furthermore, pollinator phylogeny did not explain pollinator community composition; closely related pollinators were not more likely to visit the same plant species. However, differences among pollinator communities visiting different plants were explained by plant height, floral color, and symmetry. As a result, closely related plants attracted similar numbers of pollinators. By parsing out characteristics that explain why plants share pollinators, we can identify plant species that likely share a common fate in a changing climate.

Differential pollinator effectiveness and importance in a milkweed (Asclepias, Apocynaceae) hybrid zone

PREMISE OF THE STUDY

Exceptions to the ideal of complete reproductive isolation between species are commonly encountered in diverse plant, animal, and fungal groups, but often the causative ecological processes are poorly understood. In flowering plants, the outcome of hybridization depends in part on the effectiveness of pollinators in interspecific pollen transport. In the Asclepias exaltata and A. syriaca (Apocynaceae) hybrid zone in Shenandoah National Park, Virginia, extensive introgression has been documented. The objectives of this study were to (1) determine the extent of pollinator overlap among A. exaltata, A. syriaca, and their hybrids and (2) identify the insect taxa responsible for hybridization and introgression.

METHODS

We observed focal plants of parental species and hybrids to measure visitation rate, visit duration, and per-visit pollinia removal and deposition, and we calculated pollinator effectiveness and importance.

KEY RESULTS

Visitation rates varied significantly between the 2 yr of the study. Overall, Apis mellifera, Bombus sp., and Epargyreus clarus were the most important pollinators. However, Bombus sp. was the only visitor that was observed to both remove and insert pollinia for both parent species as well as hybrids.

CONCLUSIONS

We conclude that Bombus may be a key agent of hybridization and introgression in these sympatric milkweed populations, and hybrids are neither preferred nor selected against by pollinators. Thus, we have identified a potential mechanism for how hybrids act as bridges to gene flow between A. exaltata and A. syriaca. These results provide insights into the breakdown of prezygotic isolating mechanisms.

Asymmetrical gene flow in a hybrid zone of Hawaiian Schiedea (Caryophyllaceae) species with contrasting mating systems

Asymmetrical gene flow, which has frequently been documented in naturally occurring hybrid zones, can result from various genetic and demographic factors. Understanding these factors is important for determining the ecological conditions that permitted hybridization and the evolutionary potential inherent in hybrids. Here, we characterized morphological, nuclear, and chloroplast variation in a putative hybrid zone between Schiedea menziesii and S. salicaria, endemic Hawaiian species with contrasting breeding systems. Schiedea menziesii is hermaphroditic with moderate selfing; S. salicaria is gynodioecious and wind-pollinated, with partially selfing hermaphrodites and largely outcrossed females. We tested three hypotheses: 1) putative hybrids were derived from natural crosses between S. menziesii and S. salicaria, 2) gene flow via pollen is unidirectional from S. salicaria to S. menziesii and 3) in the hybrid zone, traits associated with wind pollination would be favored as a result of pollen-swamping by S. salicaria. Schiedea menziesii and S. salicaria have distinct morphologies and chloroplast genomes but are less differentiated at the nuclear loci. Hybrids are most similar to S. menziesii at chloroplast loci, exhibit nuclear allele frequencies in common with both parental species, and resemble S. salicaria in pollen production and pollen size, traits important to wind pollination. Additionally, unlike S. menziesii, the hybrid zone contains many females, suggesting that the nuclear gene responsible for male sterility in S. salicaria has been transferred to hybrid plants. Continued selection of nuclear genes in the hybrid zone may result in a population that resembles S. salicaria, but retains chloroplast lineage(s) of S. menziesii.

Interspecific gene flow in a multispecies oak hybrid zone in the Sierra Tarahumara of Mexico

BACKGROUND AND AIMS

Interspecific gene flow can occur in many combinations among species within the genus Quercus, but simultaneous hybridization among more than two species has been rarely analysed. The present study addresses the genetic structure and morphological variation in a triple hybrid zone formed by Q. hypoleucoides, Q. scytophylla and Q. sideroxyla in north-western Mexico.

METHODS

A total of 247 trees from ten reference and 13 presumed intermediate populations were characterized using leaf shape variation and geometric morphometrics, and seven nuclear microsatellites as genetic markers. Discriminant function analysis was performed for leaf shape variation, and estimates of genetic diversity and structure, and individual Bayesian genetic assignments were obtained.

KEY RESULTS

Reference populations formed three completely distinct groups according to discriminant function analysis based on the morphological data, and showed low, but significant, genetic differentiation. Populations from the zone of contact contained individuals morphologically intermediate between pairs of species in different combinations, or even among the three species. The Bayesian admixture analysis found that three main genetic clusters best fitted the data, with good correspondence of reference populations of each species to one of the genetic clusters, but various degrees of admixture evidenced in populations from the contact area.

CONCLUSIONS

The three oak species have formed a complex hybrid zone that is geographically structured as a mosaic, and comprising a wide range of genotypes, including hybrids between different species pairs, backcrosses and probable triple hybrids.

Physiological differences among two Penstemon species and their hybrids in field and common garden environments

Hybrids can exhibit unique combinations of the physiological traits of their parents. These particular combinations may influence hybrid fitness and the evolutionary trajectory of a hybrid zone. Here, a hybrid zone between Penstemon newberryi and Penstemon davidsonii along an elevational gradient was examined, and physiological traits of parents and hybrids were measured in their native environment and a common garden. Gas exchange rates of nine different crosses were also measured. Alpine P. davidsonii had less negative pre-dawn water potential and lower water use efficiency (WUE) than its montane relative P. newberryi in a common garden and in field measurements. The species difference in WUE was attributable to lower conductance in P. newberryi in the field, but to a higher photosynthetic rate in this species in the common garden. The alpine species took less time to produce mature fruits and reached maximum photosynthetic rate at a lower temperature. Natural hybrids were intermediate for most characters. F(1) hybrids had lower conductance than progeny of natural hybrids. The intermediate WUE of natural hybrids may be one factor that allows them to persist in intermediate environments. Comparisons of different crosses suggest that the genotypic composition of hybrids influences their physiological performance.