Temperature and sex shape reproductive barriers in a climate change hotspot

Climate change is altering species ranges and reproductive interactions in existing ranges, offering species new scope to mate and hybridize. The outcomes will depend on how environmental factors shape reproductive barriers across life stages, yet this is rarely assessed across the environments that species encounter in nature. We assess prezygotic and postzygotic barriers, and their dependence on temperature and parental sex, in species of a reef-building tubeworm (Galeolaria) from a fast-warming biodiversity hotspot in southern Australia. By replicating pure and reciprocal hybrid crosses across 5 temperatures spanning species' thermal ranges, we estimate thermal tolerance curves (defining niches) for crosses and reproductive isolation at each temperature. By also replicating crosses at 3 life stages, we partition the contributions of prezygotic barriers at fertilization, postzygotic barriers at embryogenesis, and postzygotic barriers at larval development to reproductive isolation. We show that barriers are weaker at fertilization and embryogenesis, but stronger and more temperature sensitive at larval development, as species diverge in thermal niche. Asymmetry of barriers between parental sexes, moreover, suggests a complex interplay between niche differentiation and maternal inheritance. Our findings point to a key role for temperature in reproductive isolation, but also challenges for predicting the fate of isolation in future climates.

Does urbanisation lead to parallel demographic shifts across the world in a cosmopolitan plant?

Urbanisation is occurring globally, leading to dramatic environmental changes that are altering the ecology and evolution of species. In particular, the expansion of human infrastructure and the loss and fragmentation of natural habitats in cities is predicted to increase genetic drift and reduce gene flow by reducing the size and connectivity of populations. Alternatively, the 'urban facilitation model' suggests that some species will have greater gene flow into and within cities leading to higher diversity and lower differentiation in urban populations. These alternative hypotheses have not been contrasted across multiple cities. Here, we used the genomic data from the GLobal Urban Evolution project (GLUE), to study the effects of urbanisation on non-adaptive evolutionary processes of white clover (Trifolium repens) at a global scale. We found that white clover populations presented high genetic diversity and no evidence of reduced Ne linked to urbanisation. On the contrary, we found that urban populations were less likely to experience a recent decrease in effective population size than rural ones. In addition, we found little genetic structure among populations both globally and between urban and rural populations, which showed extensive gene flow between habitats. Interestingly, white clover displayed overall higher gene flow within urban areas than within rural habitats. Our study provides the largest comprehensive test of the demographic effects of urbanisation. Our results contrast with the common perception that heavily altered and fragmented urban environments will reduce the effective population size and genetic diversity of populations and contribute to their isolation.

Genomic Tools in Biological Invasions: Current State and Future Frontiers

Human activities are accelerating rates of biological invasions and climate-driven range expansions globally, yet we understand little of how genomic processes facilitate the invasion process. Although most of the literature has focused on underlying phenotypic correlates of invasiveness, advances in genomic technologies are showing a strong link between genomic variation and invasion success. Here, we consider the ability of genomic tools and technologies to (i) inform mechanistic understanding of biological invasions and (ii) solve real-world issues in predicting and managing biological invasions. For both, we examine the current state of the field and discuss how genomics can be leveraged in the future. In addition, we make recommendations pertinent to broader research issues, such as data sovereignty, metadata standards, collaboration, and science communication best practices that will require concerted efforts from the global invasion genomics community.

Assessing the invasive potential of different source populations of ragweed (Ambrosia artemisiifolia L.) through genomically informed species distribution modelling

The genetic composition of founding populations is likely to play a key role in determining invasion success. Individual genotypes may differ in habitat preference and environmental tolerance, so their ability to colonize novel environments can be highly variable. Despite the importance of genetic variation on invasion success, its influence on the potential distribution of invaders is rarely investigated. Here, we integrate population genomics and ecological niche models (ENMs) into a single framework to predict the distribution of globally invasive common ragweed (Ambrosia artemisiifolia) in Australia. We identified three genetic clusters for ragweed and used these to construct cluster-specific ENMs and characterize within-species niche differentiation. The potential range of ragweed in Australia depended on the genetic composition and continent of origin of the introduced population. Invaders originating from warmer, wetter climates had a broader potential distribution than those from cooler, drier ones. By quantifying this change, we identified source populations most likely to expand the ragweed distribution. As prevention remains the most effective method of invasive species management, our work provides a valuable way of ranking the threat posed by different populations to better inform management decisions.

Temporal collections to study invasion biology

Biological invasions represent an extraordinary opportunity to study evolution. This is because accidental or deliberate species introductions have taken place for centuries across large geographical scales, frequently prompting rapid evolutionary transitions in invasive populations. Until recently, however, the utility of invasions as evolutionary experiments has been hampered by limited information on the makeup of populations that were part of earlier invasion stages. Now, developments in ancient and historical DNA technologies, as well as the quickening pace of digitization for millions of specimens that are housed in herbaria and museums globally, promise to help overcome this obstacle. In this review, we first introduce the types of temporal data that can be used to study invasions, highlighting the timescale captured by each approach and their respective limitations. We then discuss how ancient and historical specimens as well as data available from prior invasion studies can be used to answer questions on mechanisms of (mal)adaptation, rates of evolution, or community-level changes during invasions. By bridging the gap between contemporary and historical invasive populations, temporal data can help us connect pattern to process in invasion science. These data will become increasingly important if invasions are to achieve their full potential as experiments of evolution in nature.

Haplotype-Resolved, Chromosome-Level Assembly of White Clover (Trifolium repens L., Fabaceae)

White clover (Trifolium repens L.; Fabaceae) is an important forage and cover crop in agricultural pastures around the world and is increasingly used in evolutionary ecology and genetics to understand the genetic basis of adaptation. Historically, improvements in white clover breeding practices and assessments of genetic variation in nature have been hampered by a lack of high-quality genomic resources for this species, owing in part to its high heterozygosity and allotetraploid hybrid origin. Here, we use PacBio HiFi and chromosome conformation capture (Omni-C) technologies to generate a chromosome-level, haplotype-resolved genome assembly for white clover totaling 998 Mbp (scaffold N50 = 59.3 Mbp) and 1 Gbp (scaffold N50 = 58.6 Mbp) for haplotypes 1 and 2, respectively, with each haplotype arranged into 16 chromosomes (8 per subgenome). We additionally provide a functionally annotated haploid mapping assembly (968 Mbp, scaffold N50 = 59.9 Mbp), which drastically improves on the existing reference assembly in both contiguity and assembly accuracy. We annotated 78,174 protein-coding genes, resulting in protein BUSCO completeness scores of 99.6% and 99.3% against the embryophyta_odb10 and fabales_odb10 lineage datasets, respectively.

Large haploblocks underlie rapid adaptation in the invasive weed Ambrosia artemisiifolia

Adaptation is the central feature and leading explanation for the evolutionary diversification of life. Adaptation is also notoriously difficult to study in nature, owing to its complexity and logistically prohibitive timescale. Here, we leverage extensive contemporary and historical collections of Ambrosia artemisiifolia-an aggressively invasive weed and primary cause of pollen-induced hayfever-to track the phenotypic and genetic causes of recent local adaptation across its native and invasive ranges in North America and Europe, respectively. Large haploblocks-indicative of chromosomal inversions-contain a disproportionate share (26%) of genomic regions conferring parallel adaptation to local climates between ranges, are associated with rapidly adapting traits, and exhibit dramatic frequency shifts over space and time. These results highlight the importance of large-effect standing variants in rapid adaptation, which have been critical to A. artemisiifolia's global spread across vast climatic gradients.

Climate adaptation and vulnerability of foundation species in a global change hotspot

Climate change is altering species ranges, and relative abundances within ranges, as populations become differentially adapted and vulnerable to the climates they face. Understanding present species ranges, whether species harbour and exchange adaptive variants, and how variants are distributed across landscapes undergoing rapid change, is therefore crucial to predicting responses to future climates and informing conservation strategies. Such insights are nonetheless lacking for most species of conservation concern. We assess genomic patterns of neutral variation, climate adaptation and climate vulnerability (offsets in predicted distributions of putatively adaptive variants across present and future landscapes) for sister foundation species, the marine tubeworms Galeolaria caespitosa and Galeolaria gemineoa, in a sentinel region for climate change impacts. We find that species are genetically isolated despite uncovering sympatry in their ranges, show parallel and nonparallel signals of thermal adaptation on spatial scales smaller than gene flow across their ranges, and are predicted to face different risks of maladaptation under future temperatures across their ranges. Our findings have implications for understanding local adaptation in the face of gene flow, and generate spatially explicit predictions for climatic disruption of adaptation and species distributions in coastal ecosystems that could guide experimental validation and conservation planning.

Genome assembly of an Australian native grass species reveals a recent whole-genome duplication and biased gene retention of genes involved in stress response

BACKGROUND

The adaptive significance of polyploidy has been extensively debated, and chromosome-level genome assemblies of polyploids can provide insight into this. The Australian grass Bothriochloa decipiens belongs to the BCD clade, a group with a complex history of hybridization and polyploid. This is the first genome assembly and annotation of a species that belongs to this fascinating yet complex group.

FINDINGS

Using Illumina short reads, 10X Genomics linked reads, and Hi-C sequencing data, we assembled a highly contiguous genome of B. decipiens, with a total length of 1,218.22 Mb and scaffold N50 of 42.637 Mb. Comparative analysis revealed that the species experienced a relatively recent whole-genome duplication. We clustered the 20 major scaffolds, representing the 20 chromosomes, into the 2 subgenomes of the parental species using unique repeat signatures. We found evidence of biased fractionation and differences in the activity of transposable elements between the subgenomes prior to hybridization. Duplicates were enriched for genes involved in transcription and response to external stimuli, supporting a biased retention of duplicated genes following whole-genome duplication.

CONCLUSIONS

Our results support the hypotheses of a biased retention of duplicated genes following polyploidy and point to differences in repeat activity associated with subgenome dominance. B. decipiens is a widespread species with the ability to establish across many soil types, making it a prime candidate for climate change- resilient ecological restoration of Australian grasslands. This reference genome is a valuable resource for future population genomic research on Australian grasses.

Uncovering the genomic basis of an extraordinary plant invasion

Invasive species are a key driver of the global biodiversity crisis, but the drivers of invasiveness, including the role of pathogens, remain debated. We investigated the genomic basis of invasiveness in Ambrosia artemisiifolia (common ragweed), introduced to Europe in the late 19th century, by resequencing 655 ragweed genomes, including 308 herbarium specimens collected up to 190 years ago. In invasive European populations, we found selection signatures in defense genes and lower prevalence of disease-inducing plant pathogens. Together with temporal changes in population structure associated with introgression from closely related Ambrosia species, escape from specific microbial enemies likely favored the plant's remarkable success as an invasive species.

Quantitative trait loci mapping reveals an oligogenic architecture of a rapidly adapting trait during the European invasion of common ragweed

Biological invasions offer a unique opportunity to investigate evolution over contemporary timescales. Rapid adaptation to local climates during range expansion can be a major determinant of invasion success, yet fundamental questions remain about its genetic basis. This study sought to investigate the genetic basis of climate adaptation in invasive common ragweed (Ambrosia artemisiifolia). Flowering time adaptation is key to this annual species' invasion success, so much so that it has evolved repeated latitudinal clines in size and phenology across its native and introduced ranges despite high gene flow among populations. Here, we produced a high-density linkage map (4493 SNPs) and paired this with phenotypic data from an F2 mapping population (n = 336) to identify one major and two minor quantitative trait loci (QTL) underlying flowering time and height differentiation in this species. Within each QTL interval, several candidate flowering time genes were also identified. Notably, the major flowering time QTL detected in this study was found to overlap with a previously identified haploblock (putative inversion). Multiple genetic maps of this region identified evidence of suppressed recombination in specific genotypes, consistent with inversions. These discoveries support the expectation that a concentrated genetic architecture with fewer, larger, and more tightly linked alleles should underlie rapid local adaptation during invasion, particularly when divergently adapting populations experience high levels of gene flow.

Pollen competition in hybridizing Cakile species: How does a latecomer win the race?

PREMISE

Hybridization between cross-compatible species depends on the extent of competition between alternative mates. Even if stigmatic compatibility allows for hybridization, hybridization requires the heterospecific pollen to be competitive. Here, we determined whether conspecific pollen has an advantage in the race to fertilize ovules and the potential handicap to be overcome by heterospecific pollen in invasive Cakile species.

METHODS

We used fluorescence microscopy to measure pollen tube growth after conspecific and heterospecific hand-pollination treatments. We then determined siring success in the progeny relative to the timing of heterospecific pollen arrival on the stigma using CAPS markers.

RESULTS

In the absence of pollen competition, pollination time and pollen recipient species had a significant effect on the ratio of pollen tube growth. In long-styled C. maritima (outcrosser), pollen tubes grew similarly in both directions. In short-styled C. edentula (selfer), conspecific and heterospecific pollen tubes grew differently. Cakile edentula pollen produced more pollen tubes, revealing the potential for a mating asymmetry whereby C. edentula pollen had an advantage relative to C. maritima. In the presence of pollen competition, siring success was equivalent when pollen deposition was synchronous. However, a moderate 1-h advantage in the timing of conspecific pollination resulted in almost complete assortative mating, while an equivalent delay in conspecific pollination resulted in substantial hybrid formation.

CONCLUSIONS

Hybridization can aid the establishment of invasive species through the transfer of adaptive alleles from cross-compatible species, but also lead to extinction through demographic or genetic swamping. Time of pollen arrival on the stigma substantially affected hybridization rate, pointing to the importance of pollination timing in driving introgression and genetic swamping.

Global urban environmental change drives adaptation in white clover

Urbanization transforms environments in ways that alter biological evolution. We examined whether urban environmental change drives parallel evolution by sampling 110,019 white clover plants from 6169 populations in 160 cities globally. Plants were assayed for a Mendelian antiherbivore defense that also affects tolerance to abiotic stressors. Urban-rural gradients were associated with the evolution of clines in defense in 47% of cities throughout the world. Variation in the strength of clines was explained by environmental changes in drought stress and vegetation cover that varied among cities. Sequencing 2074 genomes from 26 cities revealed that the evolution of urban-rural clines was best explained by adaptive evolution, but the degree of parallel adaptation varied among cities. Our results demonstrate that urbanization leads to adaptation at a global scale.

Local adaptation and spatiotemporal patterns of genetic diversity revealed by repeated sampling of Caenorhabditis elegans across the Hawaiian Islands

The nematode Caenorhabditis elegans is among the most widely studied organisms, but relatively little is known about its natural ecology. Genetic diversity is low across much of the globe but high in the Hawaiian Islands and across the Pacific Rim. To characterize the niche and genetic diversity of C. elegans on the Hawaiian Islands and to explore how genetic diversity might be influenced by local adaptation, we repeatedly sampled nematodes over a three‐year period, measured various environmental parameters at each sampling site, and whole‐genome sequenced the C. elegans isolates that we identified. We found that the typical Hawaiian C. elegans niche comprises moderately moist native forests at high elevations (500–1,500 m) where ambient air temperatures are cool (15–20°C). Compared to other Caenorhabditis species found on the Hawaiian Islands (e.g., Caenorhabditis briggsae and Caenorhabditis tropicalis), we found that C. elegans were enriched in native habitats. We measured levels of genetic diversity and differentiation among Hawaiian C. elegans and found evidence of seven genetically distinct groups distributed across the islands. Then, we scanned these genomes for signatures of local adaptation and identified 18 distinct regions that overlap with hyper‐divergent regions, which may be maintained by balancing selection and are enriched for genes related to environmental sensing, xenobiotic detoxification, and pathogen resistance. These results provide strong evidence of local adaptation among Hawaiian C. elegans and contribute to our understanding of the forces that shape genetic diversity on the most remote volcanic archipelago in the world.

Allen Orr and the genetics of adaptation

Over most of the 20th century, evolutionary biologists predominantly subscribed to a strong form of "micro-mutationism," in which adaptive phenotypic divergence arises from allele frequency changes at many loci, each with a small effect on the phenotype. To be sure, there were well-known examples of large-effect alleles contributing to adaptation, yet such cases were generally regarded as atypical and unrepresentative of evolutionary change in general. In 1998, Allen Orr published a landmark theoretical paper in Evolution, which showed that both small- and large-effect mutations are likely to contribute to "adaptive walks" of a population to an optimum. Coupled with a growing set of empirical examples of large-effect alleles contributing to divergence (e.g., from QTL studies), Orr's paper provided a mathematical formalism that converted many evolutionary biologists from micro-mutationism to a more pluralistic perspective on the genetic basis of evolutionary change. We revisit the theoretical insights emerging from Orr's paper within the historical context leading up to 1998, and track the influence of this paper on the field of evolutionary biology through an examination of its citations over the last two decades and an analysis of the extensive body of theoretical and empirical research that Orr's pioneering paper inspired.

The tip of the iceberg: Genome wide marker analysis reveals hidden hybridization during invasion

Biological invasions are accelerating, and invasive species can have large economic impacts as well as severe consequences for biodiversity. During invasions, species can interact, potentially resulting in hybridization. Here, we examined two Cakile species, C. edentula and C. maritima (Brassicaceae), that co-occur and may hybridize during range expansion in separate regions of the globe. Cakile edentula invaded each location first, while C. maritima established later, apparently replacing the former. We assessed the evidence for hybridization in western North America and Australia, where both species have been introduced, and identified source populations with 4561 SNPs using Genotype-by-Sequencing. Our results indicate that C. edentula in Australia originated from one region of eastern North America while in western North America it is probably from multiple sources. Cakile maritima in Australia is derived from at least two different parts of Europe while the introduction in western North America is from one. Although morphological evidence of hybridization is generally limited to mixed species populations in Australia and virtually absent elsewhere, our genetic analysis revealed relatively high levels of hybridization in Australia (58% hybrids using Admixture) and supported the presence of hybrids in western North America (16% hybrids using Admixture) and New Zealand. Hybrids might be commonly overlooked in invaders, as identification based solely on morphological traits may represent only the tip of the iceberg. Our study reveals a repeated pattern of invasion, hybridization and apparent replacement of one species by another, which offers an opportunity to investigate the role of hybridization and introgression during invasion.

Genetically Based Trait Differentiation but Lack of Trade-offs between Stress Tolerance and Performance in Introduced Canada Thistle

Trade-offs between performance and tolerance of abiotic and biotic stress have been proposed to explain both the success of invasive species and frequently observed size differences between native and introduced populations. Canada thistle seeds collected from across the introduced North American and the native European range were grown in benign and stressful conditions (nutrient stress, shading, simulated herbivory, drought, and mowing), to evaluate whether native and introduced individuals differ in performance or stress tolerance. An additional experiment assessed the strength of maternal effects by comparing plants derived from field-collected seeds with those derived from clones grown in the glasshouse. Introduced populations tended to be larger in size, but no trade-off of stress tolerance with performance was detected; introduced populations had either superior performance or equivalent trait values and survivorship in the treatment common gardens. We also detected evidence of parallel latitudinal clines of some traits in both the native and introduced ranges and associations with climate variables in some treatments, consistent with recent climate adaptation within the introduced range. Our results are consistent with rapid adaptation of introduced populations, but, contrary to predictions, the evolution of invasive traits did not come at the cost of reduced stress tolerance.

Parallel flowering time clines in native and introduced ragweed populations are likely due to adaptation

As introduced species expand their ranges, they often encounter differences in climate which are often correlated with geography. For introduced species, encountering a geographically variable climate sometimes leads to the re-establishment of clines seen in the native range. However, clines can also be caused by neutral processes, and so it is important to gather additional evidence that population differentiation is the result of selection as opposed to nonadaptive processes. Here, we examine phenotypic and genetic differences in ragweed from the native (North America) and introduced (European) ranges. We used a common garden to assess phenotypic differentiation in size and flowering time in ragweed populations. We found significant parallel clines in flowering time in both North America and Europe. Height and branch number had significant clines in North America, and, while not statistically significant, the patterns in Europe were the same. We used SNP data to assess population structure in both ranges and to compare phenotypic differentiation to neutral genetic variation. We failed to detect significant patterns of isolation by distance, geographic patterns in population structure, or correlations between the major axes of SNP variation and phenotypes or latitude of origin. We conclude that the North American clines in size and the parallel clines seen for flowering time are most likely the result of adaptation.

A lack of genetically compatible mates constrains the spread of an invasive weed

Introduced populations often experience lag times before invasion, but the mechanisms constraining rapid expansions of introduced populations are unclear. Solidago altissima is a North American native plant with highly invasive Japanese populations and introduced Australian populations that are not invasive despite the climatic and ecological suitability of the region. By contrasting Australian with Japanese populations, we tested the hypothesis that Australian population growth is limited by a lack of long-distance dispersal via seeds owing to a limited number of compatible mates. In the field, Australian populations rarely produced viable seeds. A cross-pollination experiment found that Australian plants are fertile, yet lack compatible mates within Australia. Genetic analysis revealed that Australian individuals descend from a small set of self-incompatible genetic clones, which explains the negligible seed set within Australia. Our results show that low genetic diversity, leading to mate incompatibility, inhibits invasiveness of Australian S.  altissima, and provides compelling evidence for genetic, rather than ecological, factors constraining invasion in Australia.

Rapid repeatable phenotypic and genomic adaptation following multiple introductions

Uncovering the genomic basis of repeated adaption can provide important insights into the constraints and biases that limit the diversity of genetic responses. Demographic processes such as admixture or bottlenecks affect genetic variation underlying traits experiencing selection. The impact of these processes on the genetic basis of adaptation remains, however, largely unexamined empirically. We here test repeatability in phenotypes and genotypes along parallel climatic clines within the native North American and introduced European and Australian Ambrosia artemisiifolia ranges. To do this, we combined multiple lines of evidence from phenotype-environment associations, F_ST -like outlier tests, genotype-environment associations and genotype-phenotype associations. We used 853 individuals grown in common garden from 84 sampling locations, targeting 19 phenotypes, >83 k SNPs and 22 environmental variables. We found that 17%-26% of loci with adaptive signatures were repeated among ranges, despite alternative demographic histories shaping genetic variation and genetic associations. Our results suggest major adaptive changes can occur on short timescales, with seemingly minimum impacts due to demographic changes linked to introduction. These patterns reveal some predictability of evolutionary change during range expansion, key in a world facing ongoing climate change, and rapid invasive spread.

Mating system impacts the genetic architecture of adaptation to heterogeneous environments

Self-fertilisation has consequences for variation across the genome as it reduces effective population size, effect recombination rates and pollen flow, with implications for local adaptation. We conducted simulations of divergent stabilising selection on a quantitative trait with drift, pollen flow, mutation, recombination and different outcrossing rates. We quantified trait divergence and the genetic architecture of adaptation. We conducted an F_ST outlier analysis to identify candidate loci and quantified the impact of mating system on detectability. Selfing promoted trait divergence mainly through reductions in pollen flow. Moreover, trait architecture became more diffuse with selfing. Average effect size of trait loci was lower, while the number of loci, and their clustering distance increased. The genetic architecture of selfers was also more diffuse than outcrossers for equivalent migration rates. However, when deleterious alleles were included, architectures became more concentrated in selfers, likely to be because of reductions in population size caused by mutational meltdown and impacts of background selection on N_e . Our simulations demonstrate that mating system has important impacts on adaptive divergence of traits and the genetic landscape underlying that divergence. Selfing has a significant effect on detectability of regions of the genome important for adaptation because of neutral divergence and diffuse trait architecture.

Rapid growth and defence evolution following multiple introductions

Rapid adaptation can aid invasive populations in their competitive success. Resource allocation trade-off hypotheses predict higher resource availability or the lack of natural enemies in introduced ranges allow for increased growth and reproduction, thus contributing to invasive success. Evidence for such hypotheses is however equivocal and tests among multiple ranges over productivity gradients are required to provide a better understanding of the general applicability of these theories.Using common gardens, we investigated the adaptive divergence of various constitutive and inducible defence-related traits between the native North American and introduced European and Australian ranges, while controlling for divergence due to latitudinal trait clines, individual resource budgets, and population differentiation, using >11,000 SNPs.Rapid, repeated clinal adaptation in defence-related traits was apparent despite distinct demographic histories. We also identified divergence among ranges in some defence-related traits, although differences in energy budgets among ranges may explain some, but not all, defence-related trait divergence. We do not identify a general reduction in defence in concert with an increase in growth among the multiple introduced ranges as predicted trade-off hypotheses. Synthesis: The rapid spread of invasive species is affected by a multitude of factors, likely including adaptation to climate and escape from natural enemies. Unravelling the mechanisms underlying invasives' success enhances understanding of eco-evolutionary theory and is essential to inform management strategies in the face of ongoing climate change.

OPEN RESEARCH BADGES

This article has been awarded Open Materials, Open Data, Preregistered Research Designs Badges. All materials and data are publicly accessible via the Open Science Framework at https://doi.org/10.6084/m9.figshare.8028875.v1, https://github.com/lotteanna/defence_adaptation,https://doi.org/10.1101/435271.

Purifying selection does not drive signatures of convergent local adaptation of lodgepole pine and interior spruce

BACKGROUND

Lodgepole pine (Pinus contorta) and interior spruce (Picea glauca, Picea engelmannii, and their hybrids) are distantly related conifer species. Previous studies identified 47 genes containing variants associated with environmental variables in both species, providing evidence of convergent local adaptation. However, if the intensity of purifying selection varies with the environment, clines in nucleotide diversity could evolve through linked (background) selection that would yield allele frequency-environment signatures resembling local adaptation. If similar geographic patterns in the strength of purifying selection occur in these species, this could result in the convergent signatures of local adaptation, especially if the landscape of recombination is conserved. In the present study, we investigated whether spatially/environmentally varying purifying selection could give rise to the convergent signatures of local adaptation that had previously reported.

RESULTS

We analyzed 86 lodgepole pine and 50 interior spruce natural populations spanning heterogeneous environments in western Canada where previous analyses had found signatures of convergent local adaptation. We estimated nucleotide diversity and Tajima's D for each gene within each population and calculated the strength of correlations between nucleotide diversity and environmental variables. Overall, these estimates in the genes with previously identified convergent local adaptation signatures had no similar pattern between pine and spruce. Clines in nucleotide diversity along environmental variables were found for interior spruce, but not for lodgepole pine. In spruce, genes with convergent adaption signatures showed a higher strength of correlations than genes without convergent adaption signatures, but there was no such disparity in pine, which suggests the pattern in spruce may have arisen due to a combination of selection and hybridization.

CONCLUSIONS

The results rule out purifying/background selection as a driver of convergent local adaption signatures in lodgepole pine and interior spruce.

Rapid and repeated local adaptation to climate in an invasive plant

Biological invasions provide opportunities to study evolutionary processes occurring over contemporary timescales. To explore the speed and repeatability of adaptation, we examined the divergence of life-history traits to climate, using latitude as a proxy, in the native North American and introduced European and Australian ranges of the annual plant Ambrosia artemisiifolia. We explored niche changes following introductions using climate niche dynamic models. In a common garden, we examined trait divergence by growing seeds collected across three ranges with highly distinct demographic histories. Heterozygosity-fitness associations were used to explore the effect of invasion history on potential success. We accounted for nonadaptive population differentiation using 11 598 single nucleotide polymorphisms. We revealed a centroid shift to warmer, wetter climates in the introduced ranges. We identified repeated latitudinal divergence in life-history traits, with European and Australian populations positioned at either end of the native clines. Our data indicate rapid and repeated adaptation to local climates despite the recent introductions and a bottleneck limiting genetic variation in Australia. Centroid shifts in the introduced ranges suggest adaptation to more productive environments, potentially contributing to trait divergence between the ranges.

Quantifying how constraints limit the diversity of viable routes to adaptation

Convergent adaptation occurs at the genome scale when independently evolving lineages use the same genes to respond to similar selection pressures. These patterns of genetic repeatability provide insights into the factors that facilitate or constrain the diversity of genetic responses that contribute to adaptive evolution. A first step in studying such factors is to quantify the observed amount of repeatability relative to expectations under a null hypothesis. Here, we formulate a novel index to quantify the constraints driving the observed amount of repeated adaptation in pairwise contrasts based on the hypergeometric distribution, and then generalize this for simultaneous analysis of multiple lineages. This index is explicitly based on the probability of observing a given amount of repeatability by chance under a given null hypothesis and is readily compared among different species and types of trait. We also formulate an index to quantify the effective proportion of genes in the genome that have the potential to contribute to adaptation. As an example of how these indices can be used to draw inferences, we assess the amount of repeatability observed in existing datasets on adaptation to stress in yeast and climate in conifers. This approach provides a method to test a wide range of hypotheses about how different kinds of factors can facilitate or constrain the diversity of genetic responses observed during adaptive evolution.

How many ways can evolution solve the same adaptive problem? While convergent adaptation is evident in many organisms at the phenotypic level, we are only beginning to understand how commonly this convergence extends to the genome scale. Quantifying the repeatability of adaptation at the genome scale is therefore critical for assessing how constraints affect the diversity of viable genetic responses. Here, we develop probability-based indices to quantify the deviation between observed repeatability and expectations under a range of null hypotheses, and an estimator of the proportion of loci in the genome that can contribute to adaptation. We demonstrate the usage of these indices with individual-based simulations and example datasets from yeast and conifers and discuss how they differ from previously developed approaches to studying repeatability. Because these indices are unitless, they provide a general approach to quantifying and comparing how constraints drive convergence at the genome scale across a wide range of traits and taxa.

Modularity of genes involved in local adaptation to climate despite physical linkage

BACKGROUND

Linkage among genes experiencing different selection pressures can make natural selection less efficient. Theory predicts that when local adaptation is driven by complex and non-covarying stresses, increased linkage is favored for alleles with similar pleiotropic effects, with increased recombination favored among alleles with contrasting pleiotropic effects. Here, we introduce a framework to test these predictions with a co-association network analysis, which clusters loci based on differing associations. We use this framework to study the genetic architecture of local adaptation to climate in lodgepole pine, Pinus contorta, based on associations with environments.

RESULTS

We identify many clusters of candidate genes and SNPs associated with distinct environments, including aspects of aridity and freezing, and discover low recombination rates among some candidate genes in different clusters. Only a few genes contain SNPs with effects on more than one distinct aspect of climate. There is limited correspondence between co-association networks and gene regulatory networks. We further show how associations with environmental principal components can lead to misinterpretation. Finally, simulations illustrate both benefits and caveats of co-association networks.

CONCLUSIONS

Our results support the prediction that different selection pressures favor the evolution of distinct groups of genes, each associating with a different aspect of climate. But our results went against the prediction that loci experiencing different sources of selection would have high recombination among them. These results give new insight into evolutionary debates about the extent of modularity, pleiotropy, and linkage in the evolution of genetic architectures.

Bioinformatically predicted deleterious mutations reveal complementation in the interior spruce hybrid complex

BACKGROUND

Mutation load is expected to be reduced in hybrids via complementation of deleterious alleles. While local adaptation of hybrids confounds phenotypic tests for reduced mutation load, it may be possible to assess variation in load by analyzing the distribution of putatively deleterious alleles. Here, we use this approach in the interior spruce (Picea glauca x P. engelmannii) hybrid complex, a group likely to suffer from high mutation load and in which hybrids exhibit local adaptation to intermediate conditions. We used PROVEAN to bioinformatically predict whether non-synonymous alleles are deleterious, based on conservation of the position and abnormality of the amino acid change.

RESULTS

As expected, we found that predicted deleterious alleles were at lower average allele frequencies than alleles not predicted to be deleterious. We were unable to detect a phenotypic effect on juvenile growth rate of the many rare alleles predicted to be deleterious. Both the proportion of alleles predicted to be deleterious and the proportion of loci homozygous for predicted deleterious alleles were higher in P. engelmannii (Engelmann spruce) than in P. glauca (white spruce), due to higher diversity and frequencies of rare alleles in Engelmann. Relative to parental species, the proportion of alleles predicted to be deleterious was intermediate in hybrids, and the proportion of loci homozygous for predicted deleterious alleles was lowest.

CONCLUSION

Given that most deleterious alleles are recessive, this suggests that mutation load is reduced in hybrids due to complementation of deleterious alleles. This effect may enhance the fitness of hybrids.

Expression Divergence Is Correlated with Sequence Evolution but Not Positive Selection in Conifers

The evolutionary and genomic determinants of sequence evolution in conifers are poorly understood, and previous studies have found only limited evidence for positive selection. Using RNAseq data, we compared gene expression profiles to patterns of divergence and polymorphism in 44 seedlings of lodgepole pine (Pinus contorta) and 39 seedlings of interior spruce (Picea glauca × engelmannii) to elucidate the evolutionary forces that shape their genomes and their plastic responses to abiotic stress. We found that rapidly diverging genes tend to have greater expression divergence, lower expression levels, reduced levels of synonymous site diversity, and longer proteins than slowly diverging genes. Similar patterns were identified for the untranslated regions, but with some exceptions. We found evidence that genes with low expression levels had a larger fraction of nearly neutral sites, suggesting a primary role for negative selection in determining the association between evolutionary rate and expression level. There was limited evidence for differences in the rate of positive selection among genes with divergent versus conserved expression profiles and some evidence supporting relaxed selection in genes diverging in expression between the species. Finally, we identified a small number of genes that showed evidence of site-specific positive selection using divergence data alone. However, estimates of the proportion of sites fixed by positive selection (α) were in the range of other plant species with large effective population sizes suggesting relatively high rates of adaptive divergence among conifers.

When can stress facilitate divergence by altering time to flowering?

Stressors and heterogeneity are ubiquitous features of natural environments, and theory suggests that when environmental qualities alter flowering schedules through phenotypic plasticity, assortative mating can result that promotes evolutionary divergence. Therefore, it is important to determine whether common ecological stressors induce similar changes in flowering time. We review previous studies to determine whether two important stressors, water restriction and herbivory, induce consistent flowering time responses among species; for example, how often do water restriction and herbivory both delay flowering? We focus on the direction of change in flowering time, which affects the potential for divergence in heterogeneous environments. We also tested whether these stressors influenced time to flowering and nonphenology traits using Mimulus guttatus. The literature review suggests that water restriction has variable effects on flowering time, whereas herbivory delays flowering with exceptional consistency. In the Mimulus experiment, low water and herbivory advanced and delayed flowering, respectively. Overall, our results temper theoretical predictions for evolutionary divergence due to habitat‐induced changes in flowering time; in particular, we discuss how accounting for variation in the direction of change in flowering time can either increase or decrease the potential for divergence. In addition, we caution against adaptive interpretations of stress‐induced phenology shifts.

Comparative genomics in the Asteraceae reveals little evidence for parallel evolutionary change in invasive taxa

Asteraceae, the largest family of flowering plants, has given rise to many notorious invasive species. Using publicly available transcriptome assemblies from 35 Asteraceae, including six major invasive species, we examined evidence for micro- and macro-evolutionary genomic changes associated with invasion. To detect episodes of positive selection repeated across multiple introductions, we conducted comparisons between native and introduced genotypes from six focal species and identified genes with elevated rates of amino acid change (dN/dS). We then looked for evidence of positive selection at a broader phylogenetic scale across all taxa. As invasive species may experience founder events during colonization and spread, we also looked for evidence of increased genetic load in introduced genotypes. We rarely found evidence for parallel changes in orthologous genes in the intraspecific comparisons, but in some cases we identified changes in members of the same gene family. Using among-species comparisons, we detected positive selection in 0.003-0.69% and 2.4-7.8% of the genes using site and stochastic branch-site models, respectively. These genes had diverse putative functions, including defence response, stress response and herbicide resistance, although there was no clear pattern in the GO terms. There was no indication that introduced genotypes have a higher proportion of deleterious alleles than native genotypes in the six focal species, suggesting multiple introductions and admixture mitigated the impact of drift. Our findings provide little evidence for common genomic responses in invasive taxa of the Asteraceae and hence suggest that multiple evolutionary pathways may lead to adaptation during introduction and spread in these species.

Biological invasions, climate change and genomics

The rate of biological invasions is expected to increase as the effects of climate change on biological communities become widespread. Climate change enhances habitat disturbance which facilitates the establishment of invasive species, which in turn provides opportunities for hybridization and introgression. These effects influence local biodiversity that can be tracked through genetic and genomic approaches. Metabarcoding and metagenomic approaches provide a way of monitoring some types of communities under climate change for the appearance of invasives. Introgression and hybridization can be followed by the analysis of entire genomes so that rapidly changing areas of the genome are identified and instances of genetic pollution monitored. Genomic markers enable accurate tracking of invasive species' geographic origin well beyond what was previously possible. New genomic tools are promoting fresh insights into classic questions about invading organisms under climate change, such as the role of genetic variation, local adaptation and climate pre-adaptation in successful invasions. These tools are providing managers with often more effective means to identify potential threats, improve surveillance and assess impacts on communities. We provide a framework for the application of genomic techniques within a management context and also indicate some important limitations in what can be achieved.

Conservation and divergence of gene expression plasticity following c. 140 million years of evolution in lodgepole pine (Pinus contorta) and interior spruce (Picea glauca×Picea engelmannii)

Species respond to environmental stress through a combination of genetic adaptation and phenotypic plasticity, both of which may be important for survival in the face of climatic change. By characterizing the molecular basis of plastic responses and comparing patterns among species, it is possible to identify how such traits evolve. Here, we used de novo transcriptome assembly and RNAseq to explore how patterns of gene expression differ in response to temperature, moisture, and light regime treatments in lodgepole pine (Pinus contorta) and interior spruce (a natural hybrid population of Picea glauca and Picea engelmannii). We found wide evidence for an effect of treatment on expression within each species, with 6413 and 11,658 differentially expressed genes identified in spruce and pine, respectively. Comparing patterns of expression among these species, we found that 74% of all orthologs with differential expression had a pattern that was conserved in both species, despite 140 million yr of evolution. We also found that the specific treatments driving expression patterns differed between genes with conserved versus diverged patterns of expression. We conclude that natural selection has probably played a role in shaping plastic responses to environment in these species.

The molecular basis of invasiveness: differences in gene expression of native and introduced common ragweed (Ambrosia artemisiifolia) in stressful and benign environments

Although the evolutionary and ecological processes that contribute to plant invasion have been the focus of much research, investigation into the molecular basis of invasion is just beginning. Common ragweed (Ambrosia artemisiifolia) is an annual weed native to North America and has been introduced to Europe where it has become invasive. Using a custom-designed NimbleGen oligoarray, we examined differences in gene expression between five native and six introduced populations of common ragweed in three different environments (control, light stress and nutrient stress), as well as two different time points. We identified candidate genes that may contribute to invasiveness in common ragweed based on differences in expression between native and introduced populations from Europe. Specifically, we found 180 genes where range explained a significant proportion of the variation in gene expression and a further 103 genes with a significant range by treatment interaction. Several of these genes are potentially involved in the metabolism of secondary compounds, stress response and the detoxification of xenobiotics. Previously, we found more rapid growth and greater reproductive success in introduced populations, particularly in benign and competitive (light stress) environments, and many of these candidate genes potentially underlie these growth differences. We also found expression differences among populations within each range, reflecting either local adaptation or neutral processes, although no associations with climate or latitude were identified. These data provide a first step in identifying genes that are involved with introduction success in an aggressive annual weed.

Reproductive isolation during domestication

It has been hypothesized that reproductive isolation should facilitate evolution under domestication. However, a systematic comparison of reproductive barrier strength between crops and their progenitors has not been conducted to test this hypothesis. Here, we present a systematic survey of reproductive barriers between 32 economically important crop species and their progenitors to better understand the role of reproductive isolation during the domestication process. We took a conservative approach, avoiding those types of reproductive isolation that are poorly known for these taxa (e.g., differences in flowering time). We show that the majority of crops surveyed are isolated from their progenitors by one or more reproductive barriers, despite the fact that the most important reproductive barrier in natural systems, geographical isolation, was absent, at least in the initial stages of domestication for most species. Thus, barriers to reproduction between crops and wild relatives are closely associated with domestication and may facilitate it, thereby raising the question whether reproductive isolation could be viewed as a long-overlooked "domestication trait." Some of the reproductive barriers observed (e.g., polyploidy and uniparental reproduction), however, may have been favored for reasons other than, or in addition to, their effects on gene flow.

Genomics of Compositae weeds: EST libraries, microarrays, and evidence of introgression

PREMISE OF STUDY

Weeds cause considerable environmental and economic damage. However, genomic characterization of weeds has lagged behind that of model plants and crop species. Here we describe the development of genomic tools and resources for 11 weeds from the Compositae family that will serve as a basis for subsequent population and comparative genomic analyses. Because hybridization has been suggested as a stimulus for the evolution of invasiveness, we also analyze these genomic data for evidence of hybridization.

METHODS

We generated 22 expressed sequence tag (EST) libraries for the 11 targeted weeds using Sanger, 454, and Illumina sequencing, compared the coverage and quality of sequence assemblies, and developed NimbleGen microarrays for expression analyses in five taxa. When possible, we also compared the distributions of Ks values between orthologs of congeneric taxa to detect and quantify hybridization and introgression.

RESULTS

Gene discovery was enhanced by sequencing from multiple tissues, normalization of cDNA libraries, and especially greater sequencing depth. However, assemblies from short sequence reads sometimes failed to resolve close paralogs. Substantial introgression was detected in Centaurea and Helianthus, but not in Ambrosia and Lactuca.

CONCLUSIONS

Transcriptome sequencing using next-generation platforms has greatly reduced the cost of genomic studies of nonmodel organisms, and the ESTs and microarrays reported here will accelerate evolutionary and molecular investigations of Compositae weeds. Our study also shows how ortholog comparisons can be used to approximately estimate the genome-wide extent of introgression and to identify genes that have been exchanged between hybridizing taxa.

Genetic control of invasive plants species using selfish genetic elements

Invasive plants cause substantial environmental damage and economic loss. Here, we explore the possibility that a selfish genetic element found in plants called cytoplasmic male sterility (CMS) could be exploited for weed control. CMS is caused by mutations in the mitochondrial genome that sterilize male reproductive organs. We developed an analytical model and a spatial simulation to assess the use of CMS alleles to manage weed populations. Specifically, we examined how fertility, selfing, pollen limitation and dispersal influenced extinction rate and time until extinction in populations where CMS arises. We found that the introduction of a CMS allele can cause rapid population extinction, but only under a restricted set of conditions. Both models suggest that the CMS strategy will be appropriate for species where pollen limitation is negligible, inbreeding depression is high and the fertility advantage of females over hermaphrodites is substantial. In general, spatial structure did not have a strong influence on the simulation outcome, although low pollen dispersal and intermediate levels of seed dispersal tended to reduce population extinction rates. Given these results, the introduction of CMS alleles into a population of invasive plants probably represents an effective control method for only a select number of species.

Natural selection on floral traits through male and female function in wild populations of the heterostylous daffodil Narcissus triandrus

Measurements of natural selection in hermaphrodite populations require the analysis of performance through both female and male sex functions. Here, we investigate selection on three floral traits: flower number, flower length, and corona width through both sex functions in natural populations of the tristylous daffodil Narcissus triandrus. Selection through female function was examined in six populations, and in two of these we also estimated male selection gradients using multilocus microsatellite genotyping of parents and offspring. We detected significant directional selection for flower number through female function, and significant stabilizing selection for corona width and flower length through male function. Variation in male reproductive success was strongly influenced by the distance between mates and was significantly higher than variation in female reproductive success in one population, a result consistent with Bateman's principle. However, variation through both sex functions was similar in the other population and there was a significant negative correlation between female and male fitness indicating sex-specific trade-offs in reproductive success. Selection on floral design in N. triandrus was stronger through male than female function probably because floral morphology plays an important role in promoting effective cross-pollen transfer in populations of this heterostylous species.

Geographic variation in floral morphology and style-morph ratios in a sexually polymorphic daffodil

Narcissus triandrus is a tristylous daffodil from the Iberian Peninsula that has striking geographical variation in floral morphology and style-morph ratios. Here, we investigate the relation between this variation and ecological factors to understand mechanisms governing morph ratios. We estimated morph ratios in 124 populations throughout the range of N. triandrus and measured 13 morphological traits in 35-78 populations. Sampling of morph ratios conducted over 2-10-yr intervals demonstrated strong temporal stability. Variation in floral and vegetative traits enabled statistical prediction of morph-frequency variation among populations. Latitudinal gradients in precipitation and temperature were correlated with plant and flower size, with larger flowers in northern populations associated with bumblebee visitation and stylar dimorphism. Flowers of the L- and M-morphs differed significantly from the S-morph in several size-related characters, unlike other tristylous species. This pattern and the similarity in anther positions of the L- and M-morphs suggest that the M-morph of N. triandrus originated through genetic modifiers that shortened styles of the L-morph. Our findings support the hypothesis that geographical variation in style-morph ratios is largely governed by climatic gradients in the Iberian Peninsula, which influence the floral morphology and pollination biology of N. triandrus populations.

Population structure and genetic diversity in tristylous Narcissus triandrus: insights from microsatellite and chloroplast DNA variation

We investigated cpDNA sequence and nuclear microsatellite variation among populations of the wild daffodil Narcissus triandrus to examine the role of historical vs. contemporary forces in shaping population structure, morphological differentiation and sexual-system evolution. This wide-ranging heterostylous species of the Iberian Peninsula is largely composed of two allopatric varieties (vars. cernuus and triandrus), and populations with either stylar trimorphism or dimorphism. Dimorphic populations only occur in var. triandrus, are mainly restricted to the northwestern portion of the species range, and uniformly lack the mid-styled morph (M-morph). Chloroplast DNA (cpDNA) sequence variation revealed strong geographical structuring and evidence for a fragmentation event associated with differentiation of the two varieties. In var. triandrus, population fragmentation, restricted gene flow and isolation-by-distance were also inferred. Significant differences in genetic diversity and population structure between the two varieties likely reflect historical and contemporary differences in demography and gene flow among populations. Discordance between cpDNA markers and both microsatellite and morphological variation indicate that hybridization has occurred between the two varieties at contact zones. There were no differences in genetic diversity or population structure between dimorphic and trimorphic populations, and chloroplast haplotypes were not associated with either sexual system, indicating transitions in morph structure within each maternal lineage. M-morph frequencies were positively correlated with differentiation at microsatellite loci, indicating that the evolutionary processes influencing these neutral markers also influence alleles controlling the style morphs.

Mating patterns and demography in the tristylous daffodil Narcissus triandrus

Mating patterns in plant populations are influenced by interactions between reproductive traits and ecological conditions, both factors that are likely to vary geographically. Narcissus triandrus, a wide-ranging heterostylous herb, exhibits populations with either two (dimorphic) or three (trimorphic) style morphs and displays substantial geographical variation in demographic attributes and floral morphology. Here, we investigate this variation to determine if demography, morphology, and mating system differ between the two sexual systems. Our surveys in Portugal and NW Spain indicated that dimorphic populations were less dense, of smaller size, and had larger plants and flowers compared to trimorphic populations. Outcrossing rates estimated using allozyme markers revealed similar outcrossing rates in dimorphic and trimorphic populations (t(m) dimorphic=0.759; t(m) trimorphic=0.710). All populations experienced significant inbreeding in progeny (mean F=0.143). In contrast, parental estimates of inbreeding were not significantly different from zero (mean F=0.062), implying that few inbred offspring survive to reproductive maturity due to inbreeding depression. Although the majority of inbreeding results from selfing, significant levels of biparental inbreeding were also detected in eight of the nine populations (mean s(s)-s(m)=0.081). Density was negatively associated with levels of selfing but positively associated with biparental inbreeding. Population size was positively associated with outcrossing but not biparental inbreeding. There were no consistent differences among the style morphs in outcrossing or biparental inbreeding indicating that the maintenance of trimorphism vs dimorphism is unlikely to be associated with inbreeding of maternal parents.

Female reproductive success and the evolution of mating-type frequencies in tristylous populations

In tristylous populations, mating-type frequencies are governed by negative frequency-dependent selection typically resulting in equal morph ratios at equilibrium. However, Narcissus triandrus generally exhibits long-styled (L)-biased populations with a deficiency of the mid-styled (M)-morph. Here we used a pollen-transfer model and measurements of female fertility in natural populations to investigate whether these uneven morph ratios were associated with variation in female reproductive success. Our theoretical analysis demonstrated that morph ratio bias can result from maternal fitness differences among the morphs, and that these effects were magnified by asymmetrical mating. In nine out of 15 populations of N. triandrus, seed set differed significantly among the morphs, but pollen limitation occurred in only two of 11 populations investigated. Average seed set of the M-morph was positively associated with its frequency in populations. Flower size was negatively correlated with the seed set of the M-morph. Our results suggest that interactions between mating patterns and female fertility are responsible for variation in morph frequencies and loss of the M-morph from tristylous populations of N. triandrus.