Stephen R. Keller

History, chance and adaptation during biological invasion: separating stochastic phenotypic evolution from response to selection

Introduced species often exhibit changes in genetic variation, population structure, selection regime and phenotypic traits as they colonize and expand into new ranges. For these reasons, species invasions are increasingly recognized as promising systems for studying adaptive evolution over contemporary time scales. However, changes in phenotypic traits during invasion occur under non-equilibrium demographic conditions and may reflect the influences of prior evolutionary history and chance events, as well as selection. We briefly review the evidence for phenotypic evolution and the role of selection during invasion. While there is ample evidence for evolutionary change, it is less clear if selection is the primary mechanism. We then discuss the likelihood that stochastic events shift phenotypic distributions during invasion, and argue that hypotheses of adaptation should be tested against appropriate null models. We suggest two experimental frameworks for separating stochastic evolution from adaptation: statistically accounting for phenotypic variation among putative invasion sources identified by using phylogenetic or assignment methods and by comparing estimates of differentiation within and among ranges for both traits and neutral markers (Q(ST) vs. F(ST)). Designs that incorporate a null expectation can reveal the role of history and chance in the evolutionary process, and provide greater insights into evolution during species invasions.

Genomic admixture increases fitness during a biological invasion

Abstract During biological invasions, multiple introductions can provide opportunities for admixture among genetically distinct lineages. Admixture is predicted to contribute to invasion success by directly increasing fitness through hybrid vigour or by enhancing evolutionary potential within populations. Here, we demonstrate genome‐wide admixture during an invasion that substantially boosted fitness in the cosmopolitan weed, Silene vulgaris. We identified three divergent demes in the native European range that expanded from glacial refugia and experienced historical admixture in a well‐known suture zone. During recent invasion of North America, multiple introductions created additional opportunities for admixture. In common garden experiments, recombinant genotypes from North America experienced a two‐fold increase in fitness relative to nonrecombinants, whereas recombinant genotypes from Europe showed no lasting fitness benefits. This contrast implicates hybrid vigour behind the boost in fitness and supports the hypothesis that admixture can lead to fitness increases that may catapult invasion into a new range.

HISTORICAL RANGE EXPANSION DETERMINES THE PHYLOGENETIC DIVERSITY INTRODUCED DURING CONTEMPORARY SPECIES INVASION

Abstract For a species rapidly expanding its geographic range, such as during biological invasion, most alleles in the introduced range will have their evolutionary origins in the native range. Yet, the way in which historical processes occurring over evolutionary time in the native range contribute to the diversity sampled during contemporary invasion is largely unknown. We used chloroplast DNA (cpDNA) gene genealogies and coalescent methods to study two congeneric plants, Silene latifolia and S. vulgaris. We examined how phylogenetic diversity was shaped by demographic growth and historical range expansions in the native European range, and how this history affected the diversity sampled during their recent invasion of North America. Genealogies from both species depart from neutrality, likely as a result of demographic expansion in the ancestral range, the timing of which corresponds to shortly after each species originated. However, the species differ in the spatial distribution of cpDNA lineages across the native range. Silene latifolia shows a highly significant phylogeographic structure that most likely reflects different avenues of the post-glacial expansion into northern Europe from Mediterranean refugia. By contrast, cpDNA lineages in S. vulgaris have been widely scattered across Europe during, or since, the most recent post-glacial expansion. These different evolutionary histories resulted in dramatic differences in how phylogenetic diversity was sampled during invasion of North America. In S. latifolia, relatively few, discrete invasion events from a structured native range resulted in a rather severe genetic bottleneck, but also opportunities for admixture among previously isolated lineages. In S. vulgaris, lack of genetic structure was accompanied by more representative sampling of phylogenetic diversity during invasion, and reduced potential for admixture. Our results provide clear insights into how historical processes may feed forward to influence the phylogenetic diversity of species invading new geographic ranges.

Adaptation and colonization history affect the evolution of clines in two introduced species

Phenotypic and genetic clines have long been synonymous with adaptive evolution. However, other processes (for example, migration, range expansion, invasion) may generate clines in traits or loci across geographical and environmental gradients. It is therefore important to distinguish between clines that represent adaptive evolution and those that result from selectively neutral demographic or genetic processes. We tested for the differentiation of phenotypic traits along environmental gradients using two species in the genus Silene, whilst statistically controlling for colonization history and founder effects. We sampled seed families from across the native and introduced ranges, genotyped individuals and estimated phenotypic differentiation in replicated common gardens. The results suggest that post-glacial expansion of S. vulgaris and S. latifolia involved both neutral and adaptive genetic differentiation (clines) of life history traits along major axes of environmental variation in Europe and North America. Phenotypic clines generally persisted when tested against the neutral expectation, although some clines disappeared (and one cline emerged) when the effects of genetic ancestry were statistically removed. Colonization history, estimated using genetic markers, is a useful null model for tests of adaptive trait divergence, especially during range expansion and invasion when selection and gene flow may not have reached equilibrium.

Ecological genomics meets community-level modelling of biodiversity: mapping the genomic landscape of current and future environmental adaptation

Local adaptation is a central feature of most species occupying spatially heterogeneous environments, and may factor critically in responses to environmental change. However, most efforts to model the response of species to climate change ignore intraspecific variation due to local adaptation. Here, we present a new perspective on spatial modelling of organism-environment relationships that combines genomic data and community-level modelling to develop scenarios regarding the geographic distribution of genomic variation in response to environmental change. Rather than modelling species within communities, we use these techniques to model large numbers of loci across genomes. Using balsam poplar (Populus balsamifera) as a case study, we demonstrate how our framework can accommodate nonlinear responses of loci to environmental gradients. We identify a threshold response to temperature in the circadian clock gene GIGANTEA-5 (GI5), suggesting that this gene has experienced strong local adaptation to temperature. We also demonstrate how these methods can map ecological adaptation from genomic data, including the identification of predicted differences in the genetic composition of populations under current and future climates. Community-level modelling of genomic variation represents an important advance in landscape genomics and spatial modelling of biodiversity that moves beyond species-level assessments of climate change vulnerability.

Genomic diversity, population structure, and migration following rapid range expansion in the Balsam poplar, Populus balsamifera.

Rapid range expansions can cause pervasive changes in the genetic diversity and structure of populations. The postglacial history of the Balsam Poplar, Populus balsamifera, involved the colonization of most of northern North America, an area largely covered by continental ice sheets during the last glacial maximum. To characterize how this expansion shaped genomic diversity within and among populations, we developed 412 SNP markers that we assayed for a range‐wide sample of 474 individuals sampled from 34 populations. We complemented the SNP data set with DNA sequence data from 11 nuclear loci from 94 individuals, and used coalescent analyses to estimate historical population size, demographic growth, and patterns of migration. Bayesian clustering identified three geographically separated demes found in the Northern, Central, and Eastern portions of the species’ range. These demes varied significantly in nucleotide diversity, the abundance of private polymorphisms, and population substructure. Most measures supported the Central deme as descended from the primary refuge of diversity. Both SNPs and sequence data suggested recent population growth, and coalescent analyses of historical migration suggested a massive expansion from the Centre to the North and East. Collectively, these data demonstrate the strong influence that range expansions exert on genomic diversity, both within local populations and across the range. Our results suggest that an in‐depth knowledge of nucleotide diversity following expansion requires sampling within multiple populations, and highlight the utility of combining insights from different data types in population genomic studies.

Climate-driven local adaptation of ecophysiology and phenology in balsam poplar, Populus balsamifera L. (Salicaceae)

PREMISE OF THE STUDY During past episodes of climate change, many plant species experienced large-scale range expansions. Expanding populations likely encountered strong selection as they colonized new environments. In this study we examine the extent to which populations of the widespread forest tree Populus balsamifera L. have become locally adapted as the species expanded into its current range since the last glaciation. METHODS We tested for adaptive variation in 13 ecophysiology and phenology traits on clonally propagated genotypes originating from a range-wide sample of 20 subpopulations. The hypothesis of local adaption was tested by comparing among-population variation at ecologically important traits (Q(ST)) to expected variation based on demographic history (F(ST)) estimated from a large set of nuclear single nucleotide polymorphism loci. KEY RESULTS Evidence for divergence in excess of neutral expectations was present for eight of 13 traits. Bud phenology, petiole length, and leaf nitrogen showed the greatest divergence (all Q(ST) > 0.6), whereas traits related to leaf water usage showed the least (all Q(ST) ≤ 0.30) and were not different from neutrality. Strong correlations were present between traits, geography, and climate, and they revealed a general pattern of northern subpopulations adapted to shorter, drier growing seasons compared with populations in the center or eastern regions of the range. CONCLUSIONS Our study demonstrates pronounced adaptive variation in ecophysiology and phenology among balsam poplar populations. These results suggest that as this widespread forest tree species expanded its range since the end of the last glacial maximum, it evolved rapidly in response to geographically variable selection.

De novo transcriptome assembly and polymorphism detection in the flowering plant Silene vulgaris (Caryophyllaceae)

Members of the angiosperm genus Silene are widely used in studies of ecology and evolution, but available genomic and population genetic resources within Silene remain limited. Deep transcriptome (i.e. expressed sequence tag or EST) sequencing has proven to be a rapid and cost‐effective means to characterize gene content and identify polymorphic markers in non‐model organisms. In this study, we report the results of 454 GS‐FLX Titanium sequencing of a polyA‐selected and normalized cDNA library from Silene vulgaris. The library was generated from a single pool of transcripts, combining RNA from leaf, root and floral tissue from three genetically divergent European subpopulations of S. vulgaris. A single full‐plate 454 run produced 959 520 reads totalling 363.6 Mb of sequence data with an average read length of 379.0 bp after quality trimming and removal of custom library adaptors. We assembled 832 251 (86.7%) of these reads into 40 964 contigs, which have a total length of 25.4 Mb and can be organized into 18 178 graph‐based clusters or ‘isogroups’. Assembled sequences were annotated based on homology to genes in multiple public databases. Analysis of sequence variants identified 13 432 putative single‐nucleotide polymorphisms (SNPs) and 1320 simple sequence repeats (SSRs) that are candidates for microsatellite analysis. Estimates of nucleotide diversity from 1577 contigs were used to generate genome‐wide distributions that revealed several outliers with high diversity. All of these resources are publicly available through NCBI and/or our website (http://silenegenomics.biology.virginia.edu) and should provide valuable genomic and population genetic tools for the Silene research community.

The adaptive potential of Populus balsamifera L. to phenology requirements in a warmer global climate

The manner in which organisms adapt to climate change informs a broader understanding of the evolution of biodiversity as well as conservation and mitigation plans. We apply common garden and association mapping approaches to quantify genetic variance and identify loci affecting bud flush and bud set, traits that define a trees season for height growth, in the boreal forest tree Populus balsamifera L. (balsam poplar). Using data from 478 genotypes grown in each of two common gardens, one near the southern edge and another near the northern edge of P. balsamiferas range, we found that broad‐sense heritability for bud flush and bud set was generally high (H2 > 0.5 in most cases), suggesting that abundant genetic variation exists for phenological response to changes in the length of the growing season. To identify the molecular genetic basis of this variation, we genotyped trees for 346 candidate single nucleotide polymorphisms (SNPs) from 27 candidate genes for the CO/FT pathway in poplar. Mixed‐model analyses of variance identified SNPs in 10 genes to be associated with variation in either bud flush or bud set. Multiple SNPs within FRIGIDA were associated with bud flush, whereas multiple SNPs in LEAFY and GIGANTEA 5 were associated with bud set. Although there was strong population structure in stem phenology, the geographic distribution of multilocus association SNP genotypes was widespread except at the most northern populations, indicating that geographic regions may harbour sufficient diversity in functional genes to facilitate adaption to future climatic conditions in many sites.

EPISTATIC AND CYTONUCLEAR INTERACTIONS GOVERN OUTBREEDING DEPRESSION IN THE AUTOTETRAPLOID CAMPANULASTRUM AMERICANUM

Abstract The consequences of combining divergent genomes among populations of a diploid species often involve F1 hybrid vigor followed by hybrid breakdown in later recombinant generations. As many as 70% of plant species are thought to have polyploid origins; yet little is known about the genetic architecture of divergence in polyploids and how it may differ from diploid species. We investigated the genetic architecture of population divergence using controlled crosses among five populations of the autotetraploid herb, Campanulastrum americanum. Plants were reciprocally hybridized to produce F1, F2, and F1-backcross generations that were grown with parental types in a greenhouse and measured for performance. In contrast to diploid expectations, most F1 hybrids lacked heterosis and instead showed strong outbreeding depression for early life traits. Recombinant hybrid generations often showed a recovery of performance to levels approximating, or at times even exceeding, the parental values. This pattern was also evident for an index of cumulative fitness. Analyses of line means indicated nonadditive gene action, especially forms of digenic epistasis, often influenced hybrid performance. However, standard diploid genetic models were not adequate for describing the underlying genetic architecture in a number of cases. Differences between reciprocal hybrids indicated that cytoplasmic and/or cytonuclear interactions also contributed to divergence. An enhanced role of epistasis in population differentiation may be the norm in polyploids, which have more gene copies. This study, the first of its kind on a natural autotetraploid, suggests that gene duplication may cause polyploid populations to diverge in a fundamentally different way than diploids.