Evolution repeats itself in replicate long-term studies in the wild

The extent to which evolution is repeatable remains debated. Here, we study changes over time in the frequency of cryptic color-pattern morphs in 10 replicate long-term field studies of a stick insect, each spanning at least a decade (across 30 years of total data). We find predictable "up-and-down" fluctuations in stripe frequency in all populations, representing repeatable evolutionary dynamics based on standing genetic variation. A field experiment demonstrates that these fluctuations involve negative frequency-dependent natural selection (NFDS). These fluctuations rely on demographic and selective variability that pushes populations away from equilibrium, such that they can reliably move back toward it via NFDS. Last, we show that the origin of new cryptic forms is associated with multiple structural genomic variants such that which mutations arise affects evolution at larger temporal scales. Thus, evolution from existing variation is predictable and repeatable, but mutation adds complexity even for traits evolving deterministically under natural selection.

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, USA, 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. This article is protected by copyright. All rights reserved.

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 (N. 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.

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'.

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 inter-population 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 likely 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.

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.

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.

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.

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 F₁ 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.

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.

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.

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.

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.

Large-scale mutation in the evolution of a gene complex for cryptic coloration

The types of mutations affecting adaptation in the wild are only beginning to be understood. In particular, whether structural changes shape adaptation by suppressing recombination or by creating new mutations is unresolved. Here, we show that multiple linked but recombining loci underlie cryptic color morphs of Timema chumash stick insects. In a related species, these loci are found in a region of suppressed recombination, forming a supergene. However, in seven species of Timema, we found that a megabase-size "supermutation" has deleted color loci in green morphs. Moreover, we found that balancing selection likely contributes more to maintaining this mutation than does introgression. Our results show how suppressed recombination and large-scale mutation can help to package gene complexes into discrete units of diversity such as morphs, ecotypes, or species.

Rarity does not limit genetic variation or preclude subpopulation structure in the geographically restricted desert forb Astragalus lentiginosus var. piscinensis

PREMISE OF THE STUDY

Characteristics of rare taxa include small population sizes and limited geographical ranges. The genetic consequences of rarity are poorly understood for most taxa. A small geographical range could result in reduced opportunity for isolation by distance or environment, thereby limiting genetic structure and variation, but few studies explore genetic structure at small spatial scales with sufficient resolution to test this hypothesis. Moreover, few comparative genetic studies exist among infrataxa differing in rarity. Here, we compare genetic variation among varieties of Astragalus lentiginosus differing in range size. Additionally, we ask if genetic structure exists in A. lentiginosus var. piscinensis, a rare taxon consisting of several thousand individuals that persist on ~8 km² of alkaline soil.

METHODS

We compared genetic variation among 11 varieties of A. lentiginosus differing in range size using a genotyping by sequencing (GBS) approach, which generated 11,475 single nucleotide polymorphisms (SNPs). We characterized genetic structure among subpopulations of A. lentiginosus var. piscinensis using a second GBS data set of 7274 SNPs and explored associations between genetic structure and environmental variation.

KEY RESULTS

We found no association between genetic variation and range size among varieties of A. lentiginosus. Additionally, despite the extremely small range of A. lentiginosus var. piscinensis, we report a well-defined genetic structure among subpopulations associated with microhabitat variation in soil composition.

CONCLUSIONS

Our results suggest that fine scale genetic structure may exist within other rare Astragalus taxa and that rarity does not preclude the maintenance of genetic diversity in this genus.

The genetic legacy of 50 years of desert bighorn sheep translocations

Conservation biologists have increasingly used translocations to mitigate population declines and restore locally extirpated populations. Genetic data can guide the selection of source populations for translocations and help evaluate restoration success. Bighorn sheep (Ovis canadensis) are a managed big game species that suffered widespread population extirpations across western North America throughout the early 1900s. Subsequent translocation programs have successfully re-established many formally extirpated bighorn herds, but most of these programs pre-date genetically informed management practices. The state of Nevada presents a particularly well-documented case of decline followed by restoration of extirpated herds. Desert bighorn sheep (O. c. nelsoni) populations declined to less than 3,000 individuals restricted to remnant herds in the Mojave Desert and a few locations in the Great Basin Desert. Beginning in 1968, the Nevada Department of Wildlife translocated ~2,000 individuals from remnant populations to restore previously extirpated areas, possibly establishing herds with mixed ancestries. Here, we examined genetic diversity and structure among remnant herds and the genetic consequences of translocation from these herds using a genotyping-by-sequencing approach to genotype 17,095 loci in 303 desert bighorn sheep. We found a signal of population genetic structure among remnant Mojave Desert populations, even across geographically proximate mountain ranges. Further, we found evidence of a genetically distinct, potential relict herd from a previously hypothesized Great Basin lineage of desert bighorn sheep. The genetic structure of source herds was clearly reflected in translocated populations. In most cases, herds retained genetic evidence of multiple translocation events and subsequent admixture when founded from multiple remnant source herds. Our results add to a growing literature on how population genomic data can be used to guide and monitor restoration programs.

Resource stability and geographic isolation are associated with genome divergence in western Palearctic crossbills

While many conifers produce annually variable seed crops, serotinous species (which hold seeds in cones for multiple years) represent unusually stable food resources for seed predators. Such stability is conducive to residency and potentially population divergence of consumers as exemplified by the Cassia crossbill (Loxia sinesciuris) in North America. We used genotyping by sequencing (GBS) to test whether three Mediterranean subspecies of common crossbills (L. curvirostra) associated with the serotinous Aleppo pine (Pinus halepensis) were more genetically distinct than European crossbills associated with nonserotinous conifers. We assembled a Cassia crossbill draft genome as a reference for mapping GBS reads and as a first step towards a more contiguous genome assembly. We found clear patterns of genetic divergence for each of the P. halepensis-associated subspecies. Geographic isolation, as promoted by resource stability and residency, is associated with genetic divergence of two of these subspecies. However, geographic isolation cannot account for divergence of L. c. hispana. Instead, resource stability likely contributed to divergence by reducing dispersal and increasing resource competition that may limit breeding by immigrants. In contrast, we found no differentiation among common crossbills associated with less stable resources, and only slight differentiation between common crossbills and parrot crossbills (L. pytyopsittacus). The substantial morphological divergence between common and parrot crossbills has likely originated or been maintained by selection despite gene flow generated by spatiotemporal resource fluctuation. Our results indicate that phenological as well as morphological characteristics of conifers have influenced crossbill diversification, and suggest a possible link between resource stability and population divergence.

Inconsistent reproductive isolation revealed by interactions between Catostomus fish species

Interactions between species are central to evolution and ecology, but we do not know enough about how outcomes of interactions between species vary across geographic locations, in heterogeneous environments, or over time. Ecological dimensions of interactions between species are known to vary, but evolutionary interactions such as the establishment and maintenance of reproductive isolation are often assumed to be consistent across instances of an interaction between species. Hybridization among Catostomus fish species occurs over a large and heterogeneous geographic area and across taxa with distinct evolutionary histories, which allows us to assess consistency in species interactions. We analyzed hybridization among six Catostomus species across the Upper Colorado River basin (US mountain west) and found extreme variation in hybridization across locations. Different hybrid crosses were present in different locations, despite similar species assemblages. Within hybrid crosses, hybridization varied from only first generation hybrids to extensive hybridization with backcrossing. Variation in hybridization outcomes might result from uneven fitness of hybrids across locations, polymorphism in genetic incompatibilities, chance, unidentified historical contingencies, or some combination thereof. Our results suggest caution in assuming that one or a few instances of hybridization represent all interactions between the focal species, as species interactions vary substantially across locations.

Absence of population structure across elevational gradients despite large phenotypic variation in mountain chickadees (Poecile gambeli)

Montane habitats are characterized by predictably rapid heterogeneity along elevational gradients and are useful for investigating the consequences of environmental heterogeneity for local adaptation and population genetic structure. Food-caching mountain chickadees inhabit a continuous elevation gradient in the Sierra Nevada, and birds living at harsher, high elevations have better spatial memory ability and exhibit differences in male song structure and female mate preference compared to birds inhabiting milder, low elevations. While high elevation birds breed, on average, two weeks later than low elevation birds, the extent of gene flow between elevations is unknown. Despite phenotypic variation and indirect evidence for local adaptation, population genetic analyses based on 18 073 single nucleotide polymorphisms across three transects of high and low elevation populations provided no evidence for genetic differentiation. Analyses based on individual genotypes revealed no patterns of clustering, pairwise estimates of genetic differentiation (F_ST, Nei's D) were very low, and AMOVA revealed no evidence for genetic variation structured by transect or by low and high elevation sites within transects. In addition, we found no consistent evidence for strong parallel allele frequency divergence between low and high elevation sites within the three transects. Large elevation-related phenotypic variation may be maintained by strong selection despite gene flow and future work should focus on the mechanisms underlying such variation.

Vertical stratification of the foliar fungal community in the world's tallest trees

PREMISE OF THE STUDY

The aboveground tissues of plants host numerous, ecologically important fungi, yet patterns in the spatial distribution of these fungi remain little known. Forest canopies in particular are vast reservoirs of fungal diversity, but intracrown variation in fungal communities has rarely been explored. Knowledge of how fungi are distributed throughout tree crowns will contribute to our understanding of interactions between fungi and their host trees and is a first step toward investigating drivers of community assembly for plant-associated fungi. Here we describe spatial patterns in fungal diversity within crowns of the world's tallest trees, coast redwoods (Sequoia sempervirens).

METHODS

We took a culture-independent approach, using the Illumina MiSeq platform, to characterize the fungal assemblage at multiple heights within the crown across the geographical range of the coast redwood.

KEY RESULTS

Within each tree surveyed, we uncovered evidence for vertical stratification in the fungal community; different portions of the tree crown harbored different assemblages of fungi. We also report between-tree variation in the fungal community within redwoods.

CONCLUSIONS

Our results suggest the potential for vertical stratification of fungal communities in the crowns of other tall tree species and should prompt future study of the factors giving rise to this stratification.

Intraspecific phytochemical variation shapes community and population structure for specialist caterpillars

Chemically mediated plant-herbivore interactions contribute to the diversity of terrestrial communities and the diversification of plants and insects. While our understanding of the processes affecting community structure and evolutionary diversification has grown, few studies have investigated how trait variation shapes genetic and species diversity simultaneously in a tropical ecosystem. We investigated secondary metabolite variation among subpopulations of a single plant species, Piper kelleyi (Piperaceae), using high-performance liquid chromatography (HPLC), to understand associations between plant phytochemistry and host-specialized caterpillars in the genus Eois (Geometridae: Larentiinae) and associated parasitoid wasps and flies. In addition, we used a genotyping-by-sequencing approach to examine the genetic structure of one abundant caterpillar species, Eois encina, in relation to host phytochemical variation. We found substantive concentration differences among three major secondary metabolites, and these differences in chemistry predicted caterpillar and parasitoid community structure among host plant populations. Furthermore, E. encina populations located at high elevations were genetically different from other populations. They fed on plants containing high concentrations of prenylated benzoic acid. Thus, phytochemistry potentially shapes caterpillar and wasp community composition and geographic variation in species interactions, both of which can contribute to diversification of plants and insects.

Fine-scale genetic structure among greater sage-grouse leks in central Nevada

Background

Mating systems that reduce dispersal and lead to non-random mating might increase the potential for genetic structure to arise at fine geographic scales. Greater sage-grouse (Centrocercus urophasianus) have a lek-based mating system and exhibit high site fidelity and skewed mating ratios. We quantified population structure by analyzing variation at 27,866 single-nucleotide polymorphisms in 140 males from ten leks (within five lek complexes) occurring in a small geographic region in central Nevada.

Results

Lek complexes, and to a lesser extent individual leks, formed statistically identifiable clusters in ordination analyses, providing evidence for fine-scale geographic genetic differentiation. Lek geography predicted genetic differentiation even at a small geographic scale, which could be sharpened by strong site fidelity. Relatedness was also higher among individuals within lek complexes (and leks), suggesting that reproductive skew, where few males participate in most of the successful matings, could also potentially contribute to genetic differentiation. Models incorporating a habitat resistance surface as a proxy for potentially reduced movement due to landscape features indicated that both geographic distance and habitat suitability (i.e. preferred habitat) predicted genetic structure, with no significant effect of man-made barriers to movement (i.e. power lines and roads). Finally, we illustrate how data sets containing fewer loci (<4000) had less statistical precision and failed to detect the full degree of genetic structure.

Conclusion

Our results suggest that habitat features and lek site geography of sage-grouse shape fine scale genetic structure, and highlight how larger data sets can have increased precision and accuracy for quantifying ecologically relevant genetic structure over small geographic scales.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0702-4) contains supplementary material, which is available to authorized users.

Selection on a genetic polymorphism counteracts ecological speciation in a stick insect

The interplay between selection and aspects of the genetic architecture of traits (such as linkage, dominance, and epistasis) can either drive or constrain speciation [1-3]. Despite accumulating evidence that speciation can progress to "intermediate" stages-with populations evolving only partial reproductive isolation-studies describing selective mechanisms that impose constraints on speciation are more rare than those describing drivers. The stick insect Timema cristinae provides an example of a system in which partial reproductive isolation has evolved between populations adapted to different host plant environments, in part due to divergent selection acting on a pattern polymorphism [4, 5]. Here, we demonstrate how selection on a green/melanistic color polymorphism counteracts speciation in this system. Specifically, divergent selection between hosts does not occur on color phenotypes because melanistic T. cristinae are cryptic on the stems of both host species, are resistant to a fungal pathogen, and have a mating advantage. Using genetic crosses and genome-wide association mapping, we quantify the genetic architecture of both the pattern and color polymorphism, illustrating their simple genetic control. We use these empirical results to develop an individual-based model that shows how the melanistic phenotype acts as a "genetic bridge" that increases gene flow between populations living on different hosts. Our results demonstrate how variation in the nature of selection acting on traits, and aspects of trait genetic architecture, can impose constraints on both local adaptation and speciation.

Highly variable reproductive isolation among pairs of Catostomus species

Hybridization between diverged taxa tests the strength of reproductive isolation and can therefore reveal mechanisms of reproductive isolation. However, it remains unclear how consistent reproductive isolation is across species' ranges and to what extent reproductive isolation might remain polymorphic as species diverge. To address these questions, we compared outcomes of hybridization across species pairs of Catostomus fishes in three rivers in the Upper Colorado River basin, where an introduced species, C. commersoni, hybridizes with at least two native species, C. discobolus and C. latipinnis. We observed substantial heterogeneity in outcomes of hybridization, both between species pairs and across geographically separate rivers within each species pair. We also observed hybridization of additional related species with our focal species, suggesting that reproductive isolation in this group involves interactions of multiple evolutionary and ecological factors. These findings suggest that a better understanding of the determinants of variation in reproductive isolation is needed and that studies of reproductive isolation in hybrids should consider how the dynamics and mechanisms of reproductive isolation vary over ecological space and over evolutionary time. Our results also have implications for the conservation and management of native catostomids in the Colorado River basin. Heterogeneity in outcomes of hybridization suggests that the threat posed by hybridization and genetic introgression to the persistence of native species probably varies with extent of reproductive isolation, both across rivers and across species pairs.

New dimensions of tropical diversity: an inordinate fondness for insect molecules, taxa, and trophic interactions

Most known insect species are involved in chemically mediated plant-insect multi-trophic interactions, and recent syntheses point to a substantial gap in our understanding of trophic interaction diversity, especially in the tropics. One approach to filling this gap is to examine relationships between genomic, metabolomic, taxonomic, and trophic interaction diversity via quantifying and comparing these dimensions of biodiversity at multiple scales. Innovative approaches to research on the origins and maintenance of tropical insect diversity should merge traditional approaches to natural history and taxonomy with modern measures of interaction diversity, genetic variation, and phytochemical diversity. These approaches will elucidate relationships between plant chemistry, specialization, climate, and different dimensions of biodiversity.

Stick insect genomes reveal natural selection's role in parallel speciation

Natural selection can drive the repeated evolution of reproductive isolation, but the genomic basis of parallel speciation remains poorly understood. We analyzed whole-genome divergence between replicate pairs of stick insect populations that are adapted to different host plants and undergoing parallel speciation. We found thousands of modest-sized genomic regions of accentuated divergence between populations, most of which are unique to individual population pairs. We also detected parallel genomic divergence across population pairs involving an excess of coding genes with specific molecular functions. Regions of parallel genomic divergence in nature exhibited exceptional allele frequency changes between hosts in a field transplant experiment. The results advance understanding of biological diversification by providing convergent observational and experimental evidence for selection's role in driving repeatable genomic divergence.

Experimental evidence for ecological selection on genome variation in the wild

Understanding natural selection's effect on genetic variation is a major goal in biology, but the genome-scale consequences of contemporary selection are not well known. In a release and recapture field experiment we transplanted stick insects to native and novel host plants and directly measured allele frequency changes within a generation at 186 576 genetic loci. We observed substantial, genome-wide allele frequency changes during the experiment, most of which could be attributed to random mortality (genetic drift). However, we also documented that selection affected multiple genetic loci distributed across the genome, particularly in transplants to the novel host. Host-associated selection affecting the genome acted on both a known colour-pattern trait as well as other (unmeasured) phenotypes. We also found evidence that selection associated with elevation affected genome variation, although our experiment was not designed to test this. Our results illustrate how genomic data can identify previously underappreciated ecological sources and phenotypic targets of selection.

When directional selection reduces geographic variation in traits mediating species interactions

Although we often focus on the causes of geographic variation, understanding processes that act to reduce geographic variation is also important. Here, we consider a process whereby adaptive foraging across the landscape and directional selection exerted by a conifer seed predator, the common crossbill (Loxia curvirostra), potentially act to homogenize geographic variation in the defensive traits of its prey. We measured seed predation and phenotypic selection exerted by crossbills on black pine (Pinus nigra) at two sites in the Pindos Mountains, Greece. Seed predation by crossbills was over an order of magnitude higher at the site where cone scale thickness was significantly thinner, which was also the cone trait that was the target of selection at the high predation site. Additional comparisons of selection differentials demonstrate that crossbills exert selection on black pine that is consistent in form across space and time, and increases in strength with increasing seed predation. If predators distribute themselves in relation to the defensive traits of their prey and the strength of selection predators exert is proportional to the amount of predation, then predators may act to homogenize trait variation among populations of their prey in a process analogous to coevolutionary alternation with escalation.

Development of genetic diversity, differentiation and structure over 500 years in four ponderosa pine populations

Population history plays an important role in shaping contemporary levels of genetic variation and geographic structure. This is especially true in small, isolated range-margin populations, where effects of inbreeding, genetic drift and gene flow may be more pronounced than in large continuous populations. Effects of landscape fragmentation and isolation distance may have implications for persistence of range-margin populations if they are demographic sinks. We studied four small, disjunct populations of ponderosa pine over a 500-year period. We coupled demographic data obtained through dendroecological methods with microsatellite data to discern how and when contemporary levels of allelic diversity, among and within-population levels of differentiation, and geographic structure, arose. Alleles accumulated rapidly following initial colonization, demonstrating proportionally high levels of gene flow into the populations. At population sizes of approximately 100 individuals, allele accumulation saturated. Levels of genetic differentiation among populations (F(ST) and Jost's D(est)) and diversity within populations (F(IS)) remained stable through time. There was no evidence of geographic genetic structure at any time in the populations' history. Proportionally, high gene flow in the early stages of population growth resulted in rapid accumulation of alleles and quickly created relatively homogenous genetic patterns among populations. Our study demonstrates that contemporary levels of genetic diversity were formed quickly and early in population development. How contemporary genetic diversity accumulates over time is a key facet of understanding population growth and development. This is especially relevant given the extent and speed at which species ranges are predicted to shift in the coming century.

Do highly divergent loci reside in genomic regions affecting reproductive isolation? A test using next-generation sequence data in Timema stick insects

Background

Genetic divergence during speciation with gene flow is heterogeneous across the genome, with some regions exhibiting stronger differentiation than others. Exceptionally differentiated regions are often assumed to experience reduced introgression, i.e., reduced flow of alleles from one population into another because such regions are affected by divergent selection or cause reproductive isolation. In contrast, the remainder of the genome can be homogenized by high introgression. Although many studies have documented variation across the genome in genetic differentiation, there are few tests of this hypothesis that explicitly quantify introgression. Here, we provide such a test using 38,304 SNPs in populations of Timema cristinae stick insects. We quantify whether loci that are highly divergent between geographically separated (‘allopatric’) populations exhibit unusual patterns of introgression in admixed populations. To the extent this is true, highly divergent loci between allopatric populations contribute to reproductive isolation in admixed populations.

Results

As predicted, we find a substantial association between locus-specific divergence between allopatric populations and locus-specific introgression in admixed populations. However, many loci depart from this relationship, sometimes strongly so. We also report evidence for selection against foreign alleles due to local adaptation.

Conclusions

Loci that are strongly differentiated between allopatric populations sometimes contribute to reproductive isolation in admixed populations. However, geographic variation in selection and local adaptation, in aspects of genetic architecture (such as organization of genes, recombination rate variation, number and effect size of variants contributing to adaptation, etc.), and in stochastic evolutionary processes such as drift can cause strong differentiation of loci that do not always contribute to reproductive isolation. The results have implications for the theory of ‘genomic islands of speciation’.

Genomic consequences of multiple speciation processes in a stick insect

Diverse geographical modes and mechanisms of speciation are known, and individual speciation genes have now been identified. Despite this progress, genome-wide outcomes of different evolutionary processes during speciation are less understood. Here, we integrate ecological and spatial information, mating trials, transplantation data and analysis of 86 130 single nucleotide polymorphisms (SNPs) in eight populations (28 pairwise comparisons) of Timema cristinae stick insects to test the effects of different factors on genomic divergence in a system undergoing ecological speciation. We find patterns consistent with effects of numerous factors, including geographical distance, gene flow, divergence in host plant use and climate, and selection against maladaptive hybridization (i.e. reinforcement). For example, the number of highly differentiated ‘outlier loci’, allele-frequency clines and the overall distribution of genomic differentiation were recognizably affected by these factors. Although host use has strong effects on phenotypic divergence and reproductive isolation, its effects on genomic divergence were subtler and other factors had pronounced effects. The results demonstrate how genomic data can provide new insights into speciation and how genomic divergence can be complex, yet predictable. Future work could adopt experimental, mapping and functional approaches to directly test which genetic regions are affected by selection and determine their physical location in the genome.

Genomics of isolation in hybrids

Hybrid zones are common in nature and can offer critical insights into the dynamics and components of reproductive isolation. Hybrids between diverged lineages are particularly informative about the genetic architecture of reproductive isolation, because introgression in an admixed population is a direct measure of isolation. In this paper, we combine simulations and a new statistical model to determine the extent to which different genetic architectures of isolation leave different signatures on genome-level patterns of introgression. We found that reproductive isolation caused by one or several loci of large effect caused greater heterogeneity in patterns of introgression than architectures involving many loci with small fitness effects, particularly when isolating factors were closely linked. The same conditions that led to heterogeneous introgression often resulted in a reasonable correspondence between outlier loci and the genetic loci that contributed to isolation. However, demographic conditions affected both of these results, highlighting potential limitations to the study of the speciation genomics. Further progress in understanding the genomics of speciation will require large-scale empirical studies of introgression in hybrid zones and model-based analyses, as well as more comprehensive modelling of the expected levels of isolation with different demographies and genetic architectures of isolation.

De novo characterization of the Timema cristinae transcriptome facilitates marker discovery and inference of genetic divergence

Adaptation to different ecological environments can promote speciation. Although numerous examples of such 'ecological speciation' now exist, the genomic basis of the process, and the role of gene flow in it, remains less understood. This is, at least in part, because systems that are well characterized in terms of their ecology often lack genomic resources. In this study, we characterize the transcriptome of Timema cristinae stick insects, a system that has been researched intensively in terms of ecological speciation, but for which genomic resources have not been previously developed. Specifically, we obtained >1 million 454 sequencing reads that assembled into 84,937 contigs representing approximately 18,282 unique genes and tens of thousands of potential molecular markers. Second, as an illustration of their utility, we used these genomic resources to assess multilocus genetic divergence within both an ecotype pair and a species pair of Timema stick insects. The results suggest variable levels of genetic divergence and gene flow among taxon pairs and genes and illustrate a first step towards future genomic work in Timema.

Cross-species transferability of SSR loci developed from transciptome sequencing in lodgepole pine

With the advent of next generation sequencing technologies, transcriptome level sequence collections are arising as prominent resources for the discovery of gene-based molecular markers. In a previous study more than 15,000 simple sequence repeats (SSRs) in expressed sequence tag (EST) sequences resulting from 454 pyrosequencing of Pinus contorta cDNA were identified. From these we developed PCR primers for approximately 4000 candidate SSRs. Here, we tested 184 of these SSRs for successful amplification across P. contorta and eight other pine species and examined patterns of polymorphism and allelic variability for a subset of these SSRs. Cross-species transferability was high, with high percentages of loci producing PCR products in all species tested. In addition, 50% of the loci we screened across panels of individuals from three of these species were polymorphic and allelically diverse. We examined levels of diversity in a subset of these SSRs by collecting genotypic data across several populations of Pinus ponderosa in northern Wyoming. Our results indicate the utility of mining pyrosequenced EST collections for gene-based SSRs and provide a source of molecular markers that should bolster evolutionary genetic investigations across the genus Pinus.

Assembly, gene annotation and marker development using 454 floral transcriptome sequences in Ziziphus celata (Rhamnaceae), a highly endangered, Florida endemic plant

Large-scale DNA sequence data may enable development of genetic resources in endangered species, thereby facilitating conservation efforts. Ziziphus celata, a federally endangered, self-incompatible plant species occurring in Florida, USA, is one species for which genetic resources are necessary to facilitate new introductions and augmentations essential for recovery of the species. We used 454 pyrosequencing of a Z. celata normalized floral cDNA library to create a genomic resource for gene and marker discovery. A half-plate GS-FLX Titanium run yielded 655 337 reads averaging 250 bp. A total of 474 025 reads were assembled de novo into 84 645 contigs averaging 408 bp, while 181 312 reads remained unassembled. Forty-seven and 43% of contig consensus sequences had BLAST matches to known proteins in the Uniref50 and TAIR9 annotated protein databases, respectively; many contigs fully represented orthologous proteins in TAIR9. A total of 22 707 unique genes were sequenced, indicating substantial coverage of the Z. celata transcriptome. We detected single-nucleotide polymorphisms and simple sequence repeats (SSRs) and developed thousands of SSR primers for use in future genetic studies. As a first step towards understanding self-incompatibility in Z. celata, we identified sequences belonging to the gene family encoding self-incompatibility. This study demonstrates the efficacy of 454 transcriptome sequencing for rapid gene and marker discovery in an endangered plant.

Low levels of population genetic structure in Pinus contorta (Pinaceae) across a geographic mosaic of co-evolution

PREMISE OF THE STUDY

Population genetic analyses provide information on the population context in which evolutionary processes operate and are important for understanding the evolution of geographically variable traits. Earlier studies showed that cone structure of lodgepole pine in the Rocky Mountains diverged among populations because of geographic variation in coevolutionary interactions involving mammalian and avian seed predators. Analyses of population genetic variation are needed to determine whether this divergence has arisen despite extensive gene flow and whether populations to the east and west of the Rocky Mountains have evolved convergent phenotypes independently.

METHODS

We investigated genetic structuring across 22 stands of lodgepole pine in the central Rocky Mountains and in isolated peripheral populations that experience different seed predators and exhibit parallel divergence in cone traits using a set of nine simple sequence repeats and 235 AFLP loci.

KEY RESULTS

Our analyses reveal high levels of genetic diversity within and low genetic differentiation among populations. Nonetheless, geographic and genetic distances were correlated, and isolated populations to the east and west of the Rocky Mountains had higher levels of differentiation than did populations in the central part of the range.

CONCLUSIONS

These data indicate not only that adaptive divergence of cone traits across a geographic mosaic of coevolution has occurred despite minimal genetic differentiation, but also that isolated populations to the east and west of the Rocky Mountains have evolved distinctive cones independently and in parallel. The population structure quantified here will inform future research aimed at detecting genetic variants associated with divergent adaptive traits.

Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery

Background

Massively parallel sequencing of cDNA is now an efficient route for generating enormous sequence collections that represent expressed genes. This approach provides a valuable starting point for characterizing functional genetic variation in non-model organisms, especially where whole genome sequencing efforts are currently cost and time prohibitive. The large and complex genomes of pines (Pinus spp.) have hindered the development of genomic resources, despite the ecological and economical importance of the group. While most genomic studies have focused on a single species (P. taeda), genomic level resources for other pines are insufficiently developed to facilitate ecological genomic research. Lodgepole pine (P. contorta) is an ecologically important foundation species of montane forest ecosystems and exhibits substantial adaptive variation across its range in western North America. Here we describe a sequencing study of expressed genes from P. contorta, including their assembly and annotation, and their potential for molecular marker development to support population and association genetic studies.

Results

We obtained 586,732 sequencing reads from a 454 GS XLR70 Titanium pyrosequencer (mean length: 306 base pairs). A combination of reference-based and de novo assemblies yielded 63,657 contigs, with 239,793 reads remaining as singletons. Based on sequence similarity with known proteins, these sequences represent approximately 17,000 unique genes, many of which are well covered by contig sequences. This sequence collection also included a surprisingly large number of retrotransposon sequences, suggesting that they are highly transcriptionally active in the tissues we sampled. We located and characterized thousands of simple sequence repeats and single nucleotide polymorphisms as potential molecular markers in our assembled and annotated sequences. High quality PCR primers were designed for a substantial number of the SSR loci, and a large number of these were amplified successfully in initial screening.

Conclusions

This sequence collection represents a major genomic resource for P. contorta, and the large number of genetic markers characterized should contribute to future research in this and other pines. Our results illustrate the utility of next generation sequencing as a basis for marker development and population genomics in non-model species.

Coevolution between crossbills and black pine: the importance of competitors, forest area and resource stability

Studies of predator-prey interactions have found that geographically structured coevolution has played an important role in the adaptive diversification of crossbills (Loxia spp.). We extend those studies by considering common crossbills (L. curvirostra) in the Mediterranean where they rely on seeds in the cones of black pine (Pinus nigra). On the continent, where tree squirrels (Sciurus vulgaris) are present, enhanced defenses against crossbills were most evident in larger areas of pine forest. On islands in the absence of tree squirrels, crossbills and black pine have coevolved in a predator-prey arms race on Cyprus but not Corsica. In contrast to other conifers that island endemic crossbills rely upon, black pine does not hold seeds in its cones year round. Consequently, key to the strong crossbill-pine interaction on Cyprus is likely the presence of an alternative conifer that provides seeds during early summer when black pine seeds are scarce.