Environment predicts the maintenance of reproductive isolation in a mosaic hybrid zone of rubber rabbitbrush.

Widely distributed plants of western North America experience divergent selection across environmental gradients, have complex histories shaped by biogeographic barriers and distributional shifts and often illustrate continuums of reproductive isolation. Rubber rabbitbrush (Ericameria nauseosa) is a foundational shrub species that occurs across diverse environments of western North America. Its remarkable phenotypic diversity is currently ascribed to two subspecies-Ericameria nauseosa nauseosa and Ericameria nauseosa consimilis-and 22 named varieties. To understand how genetic variation is partitioned across subspecies, varieties, and environments, we used high throughput sequencing of reduced representation libraries. We found clear evidence for divergence between the two subspecies, despite largely sympatric distributions. Numerous locations exhibiting admixed ancestry were not geographically localized but were widely distributed across a mosaic hybrid zone. The occurrence of hybrid and subspecific ancestries was strongly predicted by environmental variables as well as the proximity to major ecotones between ecoregions. Although this repeatability illustrates the importance of environmental factors in shaping reproductive isolation, variability in the prevalence of hybridization also indicates these factors likely differ across ecological contexts. There was mixed evidence for the evolutionary cohesiveness of varieties, but several genetically distinct and narrow endemic varieties exhibited admixed subspecific ancestries, hinting at the possibility for transgressive hybridization to contribute to phenotypic novelty and the colonization of new environments in E. nauseosa.

Does a history of population co-occurrence predict plant performance, community productivity, or invasion resistance?

A history of species co-occurrence in plant communities is hypothesized to lead to greater niche differentiation, more efficient resource partitioning, and more productive, resistant communities as a result of evolution in response to biotic interactions. A similar question can be asked of co-occurring populations: do individual species or community responses differ when communities are founded with plants sharing a history of population co-occurrence (sympatric) or originating from different locations (allopatric)? Using shrub, grass, and forb species from six locations in the western Great Basin, North America, we compared establishment, productivity, reproduction, phenology, and resistance to invaders for experimental communities with either sympatric or allopatric population associations. Each community type was planted with six taxa in outdoor mesocosms, measured over three growing seasons, and invaded with the annual grass Bromus tectorum in the final season. For most populations, the allopatric or sympatric status of neighbors was not important. However, in some cases, it was beneficial for some species from some locations to be planted with allopatric neighbors, while others benefited from sympatric neighbors, and some of these responses had large effects. For instance, the Elymus population that benefited the most from allopatry grew 50% larger with allopatric neighbors than in single origin mesocosms. This response affected invasion resistance, as B. tectorum biomass was strongly affected by productivity and phenology of Elymus spp., as well as Poa secunda. Our results demonstrate that, while community composition can affect plant performance in semi-arid plant communities, assembling communities from sympatric populations is not sufficient to ensure high productivity and invasion resistance. Instead, we observed an idiosyncratic interaction between sampling effects and evolutionary history, with the potential for seed source of individual populations to have community-level effects.

Hybridization in the absence of an ecotone favors hybrid success in woodrats (Neotoma spp.).

Hybridization is a common process that has broadly impacted the evolution of multicellular eukaryotes; however, how ecological factors influence this process remains poorly understood. Here, we report the findings of a 3-year recapture study of the Bryant's woodrat (Neotoma bryanti) and desert woodrat (Neotoma lepida), two species that hybridize within a creosote bush (Larrea tridentata) shrubland in Whitewater, CA, USA. We used a genotype-by-sequencing approach to characterize the ancestry distribution of individuals across this hybrid zone coupled with Cormack-Jolly-Seber modeling to describe demography. We identified a high frequency of hybridization at this site with ~40% of individuals possessing admixed ancestry, which is the result of multigenerational backcrossing and advanced hybrid-hybrid crossing. F1, F2, and advanced generation hybrids had apparent survival rates similar to parental N. bryanti, while parental and backcross N. lepida had lower apparent survival rates and were far less abundant. Compared to bimodal hybrid zones where hybrids are often rare and selected against, we find that hybrids at Whitewater are common and have comparable survival to the dominant parental species, N. bryanti. The frequency of hybridization at Whitewater is therefore likely limited by the abundance of the less common parental species, N. lepida, rather than selection against hybrids.

The gut microbiome reflects ancestry despite dietary shifts across a hybrid zone.

The microbiome is critical to an organism's phenotype, and its composition is shaped by, and a driver of, eco-evolutionary interactions. We investigated how host ancestry, habitat and diet shape gut microbial composition in a mammalian hybrid zone between Neotoma lepida and N. bryanti that occurs across an ecotone between distinct vegetation communities. We found that habitat is the primary determinant of diet, while host genotype is the primary determinant of the gut microbiome-a finding further supported by intermediate microbiome composition in first-generation hybrids. Despite these distinct primary drivers, microbial richness was correlated with diet richness, and individuals that maintained higher dietary richness had greater gut microbial community stability. Both relationships were stronger in the relative dietary generalist of the two parental species. Our findings show that host ancestry interacts with dietary habits to shape the microbiome, ultimately resulting in the phenotypic plasticity that host-microbial interactions allow.

Reduced representation sequencing to understand the evolutionary history of Torrey pine (Pinus torreyana parry) with implications for rare species conservation.

Understanding the contribution of neutral and adaptive evolutionary processes to population differentiation is often necessary for better informed management and conservation of rare species. In this study, we focused on Pinus torreyana Parry (Torrey pine), one of the world's rarest pines, endemic to one island and one mainland population in California. Small population size, low genetic diversity, and susceptibility to abiotic and biotic stresses suggest Torrey pine may benefit from interpopulation genetic rescue to preserve the species' evolutionary potential. We leveraged reduced representation sequencing to tease apart the respective contributions of stochastic and deterministic evolutionary processes to population differentiation. We applied these data to model spatial and temporal demographic changes in effective population sizes and genetic connectivity, to identify loci possibly under selection, and evaluate genetic rescue as a potential conservation strategy. Overall, we observed exceedingly low standing variation within both Torrey pine populations, reflecting consistently low effective population sizes across time, and limited genetic differentiation, suggesting maintenance of gene flow between populations following divergence. However, genome scans identified more than 2000 candidate SNPs potentially under divergent selection. Combined with previous observations indicating population phenotypic differentiation, this indicates natural selection has probably contributed to the evolution of population genetic differences. Thus, while reduced genetic diversity, small effective population size, and genetic connectivity between populations suggest genetic rescue could mitigate the adverse effects of rarity, evidence for adaptive differentiation suggests genetic mixing could disrupt adaptation. Further work evaluating the fitness consequences of inter-population admixture is necessary to empirically evaluate the trade-offs associated with genetic rescue in Torrey pine.

Testing for fitness epistasis in a transplant experiment identifies a candidate adaptive locus in Timema stick insects.

Identifying the genetic basis of adaptation is a central goal of evolutionary biology. However, identifying genes and mutations affecting fitness remains challenging because a large number of traits and variants can influence fitness. Selected phenotypes can also be difficult to know a priori, complicating top-down genetic approaches for trait mapping that involve crosses or genome-wide association studies. In such cases, experimental genetic approaches, where one maps fitness directly and attempts to infer the traits involved afterwards, can be valuable. Here, we re-analyse data from a transplant experiment involving Timema stick insects, where five physically clustered single-nucleotide polymorphisms associated with cryptic body coloration were shown to interact to affect survival. Our analysis covers a larger genomic region than past work and revealed a locus previously not identified as associated with survival. This locus resides near a gene, Punch (Pu), involved in pteridine pigments production, implying that it could be associated with an unmeasured coloration trait. However, by combining previous and newly obtained phenotypic data, we show that this trait is not eye or body coloration. We discuss the implications of our results for the discovery of traits, genes and mutations associated with fitness in other systems, as well as for supergene evolution. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Phenotypes and environment predict seedling survival for seven co-occurring Great Basin plant taxa growing with invasive grass.

Trait-environment correlations can arise from local adaptation and can identify genetically and environmentally appropriate seeds for restoration projects. However, anthropogenic changes can disrupt the relationships between traits and fitness. Finding the best seed sources for restoration may rely on describing plant traits adaptive in disturbed and invaded environments, recognizing that while traits may differ among species and functional groups, there may be similarities in the strategies that increase seedling establishment. Focusing on three grass genera, two shrub species, and two forb genera, we collected seeds of all taxa from 16 common sites in the sagebrush steppe of the western United States. We measured seed and seedling characteristics, including seed size, emergence timing, and root and shoot traits, and compiled a suite of environmental variables for each collection site. We described trait-environment associations and asked how traits or environment of origin were associated with seedling survival in invaded gardens. Sampling seven taxa from the same sites allowed us to ask how trait-environment-performance associations differ among taxa and whether natural selection favors similar traits across multiple taxa and functional groups. All taxa showed trait-environment associations consistent with local adaptation, and both environment of origin and phenotypes predicted survival in competitive restoration settings, with some commonalities among taxa. Notably, rapid emergence and larger seeds increased survival for multiple taxa. Environmental factors at collection sites, including lower slopes (especially for grasses), greater mean annual temperatures (especially for shrubs and forbs), and greater precipitation seasonality were frequently associated with increased survival. We noted one collection site with high seedling survival across all seven taxa, suggesting that conditions within some sites may result in selection for traits that increase establishment for multiple species. Thus, choosing native plant sources with the most adaptive traits, along with matching climates, will likely improve the restoration of invaded communities.

Trio-binned genomes of the woodrats Neotoma bryanti and Neotoma lepida reveal novel gene islands and rapid copy number evolution of xenobiotic metabolizing genes.

The genomic architecture underlying the origins and maintenance of biodiversity is an increasingly accessible feature of species, due in large part to third-generation sequencing and novel analytical toolsets. Applying these techniques to woodrats (Neotoma spp.) provides a unique opportunity to study how herbivores respond to environmental change. Neotoma bryanti and N. lepida independently achieved a major dietary feat in the aftermath of a natural climate change event: switching to the novel, toxic food source creosote bush (Larrea tridentata). To better understand the genetic mechanisms underlying this ability, we employed a trio binning sequencing approach with a N. bryanti × N. lepida F1 hybrid, allowing the simultaneous assembly of genomes representing each parental species. The resulting phased, chromosome-level, highly complete haploid references enabled us to explore the genomic architecture of several gene families-cytochromes P450, UDP-glucuronosyltransferases (UGTs), and ATP-binding cassette (ABC) transporters-known to play key roles in the metabolism of naturally occurring toxic dietary compounds. In addition to duplication events in the ABCG and UGT2B subfamilies, we found expansions in three P450 gene families (2A, 2B, 3A), including the evolution of multiple novel gene islands within the 2B and 3A subfamilies, which may have provided the crucial substrate for dietary adaptation. Our assemblies demonstrate that trio binning from an F1 hybrid rodent effectively recovers parental genomes from species that diverged more than a million years ago.

Phylogenomic analyses resolve relationships among garter snakes (Thamnophis: Natricinae: Colubridae) and elucidate biogeographic history and morphological evolution.

Garter snakes (Thamnophis) are a successful group of natricines endemic to North America. They have become important natural models for ecological and evolutionary research, yet prior efforts to resolve phylogenetic relationships have resulted in conflicting topologies and weak support for certain relationships. Here, we use genomic data generated with a reduced representation double-digest RADseq approach to reassess evolutionary relationships across Thamnophis. We then use the resulting phylogeny to better understand how biogeography and feeding ecology have influenced lineage diversification and morphological evolution. We recovered highly congruent and strongly supported topologies from maximum likelihood and Bayesian analyses, but some discordance with a multispecies coalescent approach. All phylogenomic estimates split Thamnophis into two clades largely defined by northern and southern North American species. Divergence time estimates and biogeographic analyses indicate a mid-Miocene origin of Thamnophis in Mexico. In addition, historic vicariant events thought to explain biogeographic patterns in other lineages (e.g., Isthmus of Tehuantepec, Rocky Mountain Range, and Trans-Mexican Volcanic Belt) appear to have influenced patterns of diversification in Thamnophis as well. Analyses of morphological traits associated with feeding ecology showed moderate to strong phylogenetic signal. Nevertheless, phylogenetic ANOVA suggested significant differences in certain cranial morphologies between aquatic specialists and garter snakes that are terrestrial-aquatic generalists, independent of evolutionary history. Our new estimate of Thamnophis phylogeny yields an improved understanding of the biogeographic history and morphological evolution of garter snakes, and provides a robust framework for future research on these snakes.

Genomic and common garden approaches yield complementary results for quantifying environmental drivers of local adaptation in rubber rabbitbrush, a foundational Great Basin shrub.

The spatial structure of genomic and phenotypic variation across populations reflects historical and demographic processes as well as evolution via natural selection. Characterizing such variation can provide an important perspective for understanding the evolutionary consequences of changing climate and for guiding ecological restoration. While evidence for local adaptation has been traditionally evaluated using phenotypic data, modern methods for generating and analyzing landscape genomic data can directly quantify local adaptation by associating allelic variation with environmental variation. Here, we analyze both genomic and phenotypic variation of rubber rabbitbrush (Ericameria nauseosa), a foundational shrub species of western North America. To quantify landscape genomic structure and provide perspective on patterns of local adaptation, we generated reduced representation sequencing data for 17 wild populations (222 individuals; 38,615 loci) spanning a range of environmental conditions. Population genetic analyses illustrated pronounced landscape genomic structure jointly shaped by geography and environment. Genetic-environment association (GEA) analyses using both redundancy analysis (RDA) and a machine-learning approach (Gradient Forest) indicated environmental variables (precipitation seasonality, slope, aspect, elevation, and annual precipitation) influenced spatial genomic structure and were correlated with allele frequency shifts indicative of local adaptation at a consistent set of genomic regions. We compared our GEA-based inference of local adaptation with phenotypic data collected by growing seeds from each population in a greenhouse common garden. Population differentiation in seed weight, emergence, and seedling traits was associated with environmental variables (e.g., precipitation seasonality) that were also implicated in GEA analyses, suggesting complementary conclusions about the drivers of local adaptation across different methods and data sources. Our results provide a baseline understanding of spatial genomic structure for E. nauseosa across the western Great Basin and illustrate the utility of GEA analyses for detecting the environmental causes and genetic signatures of local adaptation in a widely distributed plant species of restoration significance.

Hierarchical genetic structure and implications for conservation of the world's largest salmonid, Hucho taimen.

Population genetic analyses can evaluate how evolutionary processes shape diversity and inform conservation and management of imperiled species. Taimen (Hucho taimen), the world's largest freshwater salmonid, is threatened, endangered, or extirpated across much of its range due to anthropogenic activity including overfishing and habitat degradation. We generated genetic data using high throughput sequencing of reduced representation libraries for taimen from multiple drainages in Mongolia and Russia. Nucleotide diversity estimates were within the range documented in other salmonids, suggesting moderate diversity despite widespread population declines. Similar to other recent studies, our analyses revealed pronounced differentiation among the Arctic (Selenge) and Pacific (Amur and Tugur) drainages, suggesting historical isolation among these systems. However, we found evidence for finer-scale structure within the Pacific drainages, including unexpected differentiation between tributaries and the mainstem of the Tugur River. Differentiation across the Amur and Tugur basins together with coalescent-based demographic modeling suggests the ancestors of Tugur tributary taimen likely diverged in the eastern Amur basin, prior to eventual colonization of the Tugur basin. Our results suggest the potential for differentiation of taimen at different geographic scales, and suggest more thorough geographic and genomic sampling may be needed to inform conservation and management of this iconic salmonid.

Phytochemistry reflects different evolutionary history in traditional classes versus specialized structural motifs.

Foundational hypotheses addressing plant-insect codiversification and plant defense theory typically assume a macroevolutionary pattern whereby closely related plants have similar chemical profiles. However, numerous studies have documented variation in the degree of phytochemical trait lability, raising the possibility that phytochemical evolution is more nuanced than initially assumed. We utilize proton nuclear magnetic resonance (1H NMR) data, chemical classification, and double digest restriction-site associated DNA sequencing (ddRADseq) to resolve evolutionary relationships and characterize the evolution of secondary chemistry in the Neotropical plant clade Radula (Piper; Piperaceae). Sequencing data substantially improved phylogenetic resolution relative to past studies, and spectroscopic characterization revealed the presence of 35 metabolite classes. Metabolite classes displayed phylogenetic signal, whereas the crude 1H NMR spectra featured little evidence of phylogenetic signal in multivariate tests of chemical resonances. Evolutionary correlations were detected in two pairs of compound classes (flavonoids with chalcones; p-alkenyl phenols with kavalactones), where the gain or loss of a class was dependent on the other's state. Overall, the evolution of secondary chemistry in Radula is characterized by strong phylogenetic signal of traditional compound classes and weak phylogenetic signal of specialized chemical motifs, consistent with both classic evolutionary hypotheses and recent examinations of phytochemical evolution in young lineages.

Multigenerational backcrossing and introgression between two woodrat species at an abrupt ecological transition.

When organisms experience secondary contact after allopatric divergence, genomic regions can introgress differentially depending on their relationships with adaptation, reproductive isolation, recombination, and drift. Analyses of genome-wide patterns of divergence and introgression could provide insight into the outcomes of hybridization and the potential relationship between allopatric divergence and reproductive isolation. Here, we generate population genetic data (26,262 SNPs; 353 individuals) using a reduced-representation sequencing approach to quantify patterns of ancestry, differentiation, and introgression between a pair of ecologically distinct mammals-the desert woodrat (N. lepida) and Bryant's woodrat (N. bryanti)-that hybridize at a sharp ecotone in southern California. Individual ancestry estimates confirmed that hybrids were rare in this bimodal hybrid zone, and entirely consisted of a few F1 individuals and a broad range of multigenerational backcrosses. Genomic cline analyses indicated more than half of loci had elevated introgression from one genomic background into the other. However, introgression was not associated with relative or absolute measures of divergence, and loci with extreme values for both were not typically found near detoxification enzymes previously implicated in dietary specialization for woodrats. The decoupling of differentiation and introgression suggests that processes other than adaptation, such as drift, may underlie the extreme clines at this contact zone.

A suite of rare microbes interacts with a dominant, heritable, fungal endophyte to influence plant trait expression.

Endophytes are microbes that live, for at least a portion of their life history, within plant tissues. Endophyte assemblages are often composed of a few abundant taxa and many infrequently observed, low-biomass taxa that are, in a word, rare. The ways in which most endophytes affect host phenotype are unknown; however, certain dominant endophytes can influence plants in ecologically meaningful ways-including by affecting growth and immune system functioning. In contrast, the effects of rare endophytes on their hosts have been unexplored, including how rare endophytes might interact with abundant endophytes to shape plant phenotype. Here, we manipulate both the suite of rare foliar endophytes (including both fungi and bacteria) and Alternaria fulva-a vertically transmitted and usually abundant fungus-within the fabaceous forb Astragalus lentiginosus. We report that rare, low-biomass endophytes affected host size and foliar %N, but only when the heritable fungal endophyte (A. fulva) was not present. A. fulva also reduced plant size and %N, but these deleterious effects on the host could be offset by a negative association we observed between this heritable fungus and a foliar pathogen. These results demonstrate how interactions among endophytic taxa determine the net effects on host plants and suggest that the myriad rare endophytes within plant leaves may be more than a collection of uninfluential, commensal organisms, but instead have meaningful ecological roles.

Genomic variation in the American pika: signatures of geographic isolation and implications for conservation.

BACKGROUND: Distributional responses by alpine taxa to repeated, glacial-interglacial cycles throughout the last two million years have significantly influenced the spatial genetic structure of populations. These effects have been exacerbated for the American pika (Ochotona princeps), a small alpine lagomorph constrained by thermal sensitivity and a limited dispersal capacity. As a species of conservation concern, long-term lack of gene flow has important consequences for landscape genetic structure and levels of diversity within populations. Here, we use reduced representation sequencing (ddRADseq) to provide a genome-wide perspective on patterns of genetic variation across pika populations representing distinct subspecies. To investigate how landscape and environmental features shape genetic variation, we collected genetic samples from distinct geographic regions as well as across finer spatial scales in two geographically proximate mountain ranges of eastern Nevada. RESULTS: Our genome-wide analyses corroborate range-wide, mitochondrial subspecific designations and reveal pronounced fine-scale population structure between the Ruby Mountains and East Humboldt Range of eastern Nevada. Populations in Nevada were characterized by low genetic diversity (π = 0.0006-0.0009; θW = 0.0005-0.0007) relative to populations in California (π = 0.0014-0.0019; θW = 0.0011-0.0017) and the Rocky Mountains (π = 0.0025-0.0027; θW = 0.0021-0.0024), indicating substantial genetic drift in these isolated populations. Tajima's D was positive for all sites (D = 0.240-0.811), consistent with recent contraction in population sizes range-wide. CONCLUSIONS: Substantial influences of geography, elevation and climate variables on genetic differentiation were also detected and may interact with the regional effects of anthropogenic climate change to force the loss of unique genetic lineages through continued population extirpations in the Great Basin and Sierra Nevada.

Model-based genotype and ancestry estimation for potential hybrids with mixed-ploidy.

Non-random mating among individuals can lead to spatial clustering of genetically similar individuals and population stratification. This deviation from panmixia is commonly observed in natural populations. Consequently, individuals can have parentage in single populations or involving hybridization between differentiated populations. Accounting for this mixture and structure is important when mapping the genetics of traits and learning about the formative evolutionary processes that shape genetic variation among individuals and populations. Stratified genetic relatedness among individuals is commonly quantified using estimates of ancestry that are derived from a statistical model. Development of these models for polyploid and mixed-ploidy individuals and populations has lagged behind those for diploids. Here, we extend and test a hierarchical Bayesian model, called entropy, which can use low-depth sequence data to estimate genotype and ancestry parameters in autopolyploid and mixed-ploidy individuals (including sex chromosomes and autosomes within individuals). Our analysis of simulated data illustrated the trade-off between sequencing depth and genome coverage and found lower error associated with low-depth sequencing across a larger fraction of the genome than with high-depth sequencing across a smaller fraction of the genome. The model has high accuracy and sensitivity as verified with simulated data and through analysis of admixture among populations of diploid and tetraploid Arabidopsis arenosa.

Genome-wide RAD sequencing resolves the evolutionary history of serrate leaf Juniperus and reveals discordance with chloroplast phylogeny.

Juniper (Juniperus) is an ecologically important conifer genus of the Northern Hemisphere, the members of which are often foundational tree species of arid regions. The serrate leaf margin clade is native to topologically variable regions in North America, where hybridization has likely played a prominent role in their diversification. Here we use a reduced-representation sequencing approach (ddRADseq) to generate a phylogenomic data set for 68 accessions representing all 22 species in the serrate leaf margin clade, as well as a number of close and distant relatives, to improve understanding of diversification in this group. Phylogenetic analyses using three methods (SVDquartets, maximum likelihood, and Bayesian) yielded highly congruent and well-resolved topologies. These phylogenies provided improved resolution relative to past analyses based on Sanger sequencing of nuclear and chloroplast DNA, and were largely consistent with taxonomic expectations based on geography and morphology. Calibration of a Bayesian phylogeny with fossil evidence produced divergence time estimates for the clade consistent with a late Oligocene origin in North America, followed by a period of elevated diversification between 12 and 5 Mya. Comparison of the ddRADseq phylogenies with a phylogeny based on Sanger-sequenced chloroplast DNA revealed five instances of pronounced discordance, illustrating the potential for chloroplast introgression, chloroplast transfer, or incomplete lineage sorting to influence organellar phylogeny. Our results improve understanding of the pattern and tempo of diversification in Juniperus, and highlight the utility of reduced-representation sequencing for resolving phylogenetic relationships in non-model organisms with reticulation and recent divergence.

The cost of travel: How dispersal ability limits local adaptation in host-parasite interactions.

Classical theory suggests that parasites will exhibit higher fitness in sympatric relative to allopatric host populations (local adaptation). However, evidence for local adaptation in natural host-parasite systems is often equivocal, emphasizing the need for infection experiments conducted over realistic geographic scales and comparisons among species with varied life history traits. Here, we used infection experiments to test how two trematode (flatworm) species (Paralechriorchis syntomentera and Ribeiroia ondatrae) with differing dispersal abilities varied in the strength of local adaptation to their amphibian hosts. Both parasites have complex life cycles involving sequential transmission among aquatic snails, larval amphibians and vertebrate definitive hosts that control dispersal across the landscape. By experimentally pairing 26 host-by-parasite population infection combinations from across the western USA with analyses of host and parasite spatial genetic structure, we found that increasing geographic distance-and corresponding increases in host population genetic distance-reduced infection success for P. syntomentera, which is dispersed by snake definitive hosts. For the avian-dispersed R. ondatrae, in contrast, the geographic distance between the parasite and host populations had no influence on infection success. Differences in local adaptation corresponded to parasite genetic structure; although populations of P. syntomentera exhibited ~10% mtDNA sequence divergence, those of R. ondatrae were nearly identical (<0.5%), even across a 900 km range. Taken together, these results offer empirical evidence that high levels of dispersal can limit opportunities for parasites to adapt to local host populations.

How specialized is a soil specialist? Early life history responses of a rare Eriogonum to site-level variation in volcanic soils.

PREMISE: Understanding edaphic specialization is crucial for conserving rare plants that may need relocation due to habitat loss. Focusing on Eriogonum crosbyae, a rare soil specialist in the Great Basin of the United States, we asked how site-level variation among volcanic soil outcrops affected plant growth and population distribution. METHODS: We measured emergence, survival, size, and biomass allocation of E. crosbyae seedlings planted in soils collected from 42 outcrops of actual and potential habitat. We also measured phenotypic variation in the wild, documented abiotic and biotic components of E. crosbyae habitat, re-surveyed Nevada populations, and evaluated occupancy changes over time. RESULTS: Plants responded plastically to edaphic variation, growing larger and allocating relatively more to aboveground tissues in soils with greater nutrient availability and growing smaller in soils higher in copper in the field and the greenhouse. However, the chemical and physical soil properties we measured did not predict site occupancy, nor was plant phenotype in the greenhouse different when plants were grown in soils from sites with different occupation status. We observed occupation status reversals at five locations. CONCLUSIONS: Eriogonum crosbyae performed well in soils formed on hydrothermally altered rocks that are inhospitable to many other plants. Extirpation/colonization events observed were consistent with metapopulation dynamics, which may partially explain the patchy distribution of E. crosbyae among outcrops of potential habitat. While soil properties did not predict site occupancy, early life stages showed sensitivity to soil variation, indicating that seedling dynamics may be important to consider for the conservation of this soil specialist.

Ecology shapes epistasis in a genotype-phenotype-fitness map for stick insect colour.

Genetic interactions such as epistasis are widespread in nature and can shape evolutionary dynamics. Epistasis occurs due to nonlinearity in biological systems, which can arise via cellular processes that convert genotype to phenotype and via selective processes that connect phenotype to fitness. Few studies in nature have connected genotype to phenotype to fitness for multiple potentially interacting genetic variants. Thus, the causes of epistasis in the wild remain poorly understood. Here, we show that epistasis for fitness is an emergent and predictable property of nonlinear selective processes. We do so by measuring the genetic basis of cryptic colouration and survival in a field experiment with stick insects. We find that colouration shows a largely additive genetic basis but with some effects of epistasis that enhance differentiation between colour morphs. In terms of fitness, different combinations of loci affecting colouration confer high survival in one host-plant treatment. Specifically, nonlinear correlational selection for specific combinations of colour traits in this treatment drives the emergence of pairwise and higher-order epistasis for fitness at loci underlying colour. In turn, this results in a rugged fitness landscape for genotypes. In contrast, fitness epistasis was dampened in another treatment, where selection was weaker. Patterns of epistasis that are shaped by ecologically based selection could be common and central to understanding fitness landscapes, the dynamics of evolution and potentially other complex systems.