skip to main content


Title: Drift, selection and adaptive variation in small populations of a threatened rattlesnake
Abstract

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation ofS. catenatusand its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long‐ and short‐term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre‐date the effects of recent drift, and that functional variation in these loci persists despite small short‐term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of “drift debt,” a nonequilibrium state where present‐day levels of variation overestimate the amount of functional genetic diversity present in future populations.

 
more » « less
Award ID(s):
1638872 1145978 1638902
NSF-PAR ID:
10456889
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Molecular Ecology
Volume:
29
Issue:
14
ISSN:
0962-1083
Page Range / eLocation ID:
p. 2612-2625
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Effective population size affects the efficacy of selection, rate of evolution by drift and neutral diversity levels. When species are subdivided into multiple populations connected by gene flow, evolutionary processes can depend on global or local effective population sizes. Theory predicts that high levels of diversity might be maintained by gene flow, even very low levels of gene flow, consistent with species long‐term effective population size, but tests of this idea are mostly lacking. Here, we show thatLycaeidesbutterfly populations maintain low contemporary (variance) effective population sizes (e.g. ~200 individuals) and thus evolve rapidly by genetic drift. However, populations harboured high levels of genetic diversity consistent with an effective population size several orders of magnitude larger. We hypothesized that the differences in the magnitude and variability of contemporary versus long‐term effective population sizes were caused by gene flow of sufficient magnitude to maintain diversity but only subtly affect evolution on generational timescales. Consistent with this hypothesis, we detected low but nontrivial gene flow among populations. Furthermore, using short‐term population‐genomic time‐series data, we documented patterns consistent with predictions from this hypothesis, including a weak but detectable excess of evolutionary change in the direction of the mean (migrant gene pool) allele frequencies across populations and consistency in the direction of allele frequency change over time. The documented decoupling of diversity levels and short‐term change by drift inLycaeideshas implications for our understanding of contemporary evolution and the maintenance of genetic variation in the wild.

     
    more » « less
  2. Abstract

    Disentangling the effects of neutral and adaptive processes in maintaining phenotypic variation across environmental gradients is challenging in natural populations. Song sparrows (Melospiza melodia) on the California Channel Islands occupy a pronounced east‐west climate gradient within a small spatial scale, providing a unique opportunity to examine the interaction of genetic isolation (reduced gene flow) and the environment (selection) in driving variation. We used reduced representation genomic libraries to infer the role of neutral processes (drift and restricted gene flow) and divergent selection in driving variation in thermoregulatory traits with an emphasis on the mechanisms that maintain bill divergence among islands. Analyses of 22,029 neutral SNPs confirm distinct population structure by island with restricted gene flow and relatively large effective population sizes, suggesting bill differences are probably not a product of genetic drift. Instead, we found strong support for local adaptation using 3294 SNPs in differentiation‐based and environmental association analyses coupled with genome‐wide association tests. Specifically, we identified several putatively adaptive and candidate loci in or near genes involved in bill development pathways (e.g.,BMP,CaM,Wnt), confirming the highly complex and polygenic architecture underlying bill morphology. Furthermore, we found divergence in genes associated with other thermoregulatory traits (i.e., feather structure, plumage colour, and physiology). Collectively, these results suggest strong divergent selection across an island archipelago results in genomic changes in a suite of traits associated with climate adaptation over small spatial scales. Future research should move beyond studying univariate traits to better understand multidimensional responses to complex environmental conditions.

     
    more » « less
  3. Abstract Context

    Processes that shape genomic and ecological divergence can reveal important evolutionary dynamics to inform the conservation of threatened species.Fontaineais a genus of rainforest shrubs and small trees including critically endangered and threatened species restricted to narrow, but complex geographic and ecological regions. Several species ofFontaineaare subject to spatially explicit conditions and experience limited intra-specific gene flow, likely generating genetic differentiation and local adaptation.

    Objectives

    Here, we explored the genetic and ecological mechanisms underlying patterns of diversification in two, closely related threatenedFontaineaspecies. Our aim was to compare spatial patterns of genetic variation between the vulnerableFontainea australis(Southern Fontainea) and critically endangeredF. oraria(Coastal Fontainea), endemic to the heterogeneous subtropical region of central, eastern Australia, where large-scale clearing has severely reduced rainforest habitat to a fraction (< 1%) of its pre-European settlement extent.

    Methods

    We used a set of 10,000 reduced-representation markers to infer genetic relationships and the drivers of spatial genetic variation across the two species. In addition, we employed a combination of univariate and multivariate genome-environment association analysis using a set of topo-climatic variables to explore potential patterns of local adaptation as a factor impacting genomic divergence.

    Results

    Our study revealed that Coastal Fontainea have a close genetic relationship with Southern Fontainea. We showed that isolation by distance has played a key role in their genetic variation, indicating that vicariance can explain the spatial genetic distribution of the two species. Genotype-environment analyses showed a strong association with temperature and topographic features, suggesting adaptation to localised thermal environments. We used a multivariate redundancy analysis to identify a range of putatively adapted loci associated with local environmental conditions.

    Conclusions

    Divergent selection at the local-habitat scale as a result of dispersal limitations and environmental heterogeneity (including physical barriers) are likely contributors to adaptive divergence between the twoFontaineaspecies. Our findings have presented evidence to indicate that Southern and Coastal Fontainea were comprised of distinct genetic groups and ecotypes, that together may form a single species continuum, with further phenotype research suggested to confirm the current species boundaries. Proactive conservation actions, including assisted migration to enhance the resilience of populations lacking stress-tolerant single nucleotide polymorphisms (SNPs) may be required to secure the long-term future of both taxa. This is especially vital for the critically endangered Coastal Fontainea given projections of habitat decline for the species under future climate scenarios.

     
    more » « less
  4. Abstract

    Multiple highly divergent lineages have been identified withinLigia occidentalis sensu lato, a rocky supralittoral isopod distributed along a ~3000 km latitudinal gradient that encompasses several proposed marine biogeographic provinces and ecoregions in the eastern Pacific. Highly divergent lineages have nonoverlapping geographic distributions, with distributional limits that generally correspond with sharp environmental changes. Crossbreeding experiments suggest postmating reproductive barriers exist among some of them, and surveys of mitochondrial and nuclear gene markers do not show evidence of hybridization. Populations are highly isolated, some of which appear to be very small; thus, the effects of drift are expected to reduce the efficiency of selection. Large genetic divergences among lineages, marked environmental differences in their ranges, reproductive isolation, and/or high isolation of populations may have resulted in morphological differences inL. occidentalis, not detected yet by traditional taxonomy. We used landmark‐based geometric morphometric analyses to test for differences in body shape among highly divergent lineages ofL. occidentalis, and among populations within these lineages. We analyzed a total of 492 individuals from 53 coastal localities from the southern California Bight to Central Mexico, including the Gulf of California. We conducted discriminant function analyses (DFAs) on body shape morphometrics to assess morphological variation among genetically differentiated lineages and their populations. We also tested for associations between phylogeny and morphological variation, and whether genetic divergence is correlated to multivariate morphological divergence. We detected significant differences in body shape among highly divergent lineages, and among populations within these lineages. Nonetheless, neither lineages nor populations can be discriminated on the basis of body shape, because correct classification rates of cross‐validatedDFAs were low. Genetic distance and phylogeny had weak to no effect on body shape variation. The supralittoral environment appears to exert strong stabilizing selection and/or strong functional constraints on body shape inL. occidentalis, thereby leading to morphological stasis in this isopod.

     
    more » « less
  5. Abstract

    Theory predicts that threatened species living in small populations will experience high levels of inbreeding that will increase their genetic load, but recent work suggests that the impact of load may be minimized by purging resulting from long‐term population bottlenecks. Empirical studies that examine this idea using genome‐wide estimates of inbreeding and genetic load in threatened species are limited. Here we use individual genome resequencing data to compare levels of inbreeding, levels of genetic load (estimated as mutation load) and population history in threatened Eastern massasauga rattlesnakes (Sistrurus catenatus), which exist in small isolated populations, and closely related yet outbred Western massasauga rattlesnakes (Sistrurus tergeminus). In terms of inbreeding,Scatenatusgenomes had a greater number of runs of homozygosity of varying sizes, indicating sustained inbreeding through repeated bottlenecks when compared toStergeminus. At the species level, outbredStergeminushad higher genome‐wide levels of mutation load in the form of greater numbers of derived deleterious mutations compared toScatenatus, presumably due to long‐term purging of deleterious mutations inScatenatus. In contrast, mutations that escaped species‐level drift effects withinScatenatuspopulations were in general more frequent and more often found in homozygous genotypes than inStergeminus, suggesting a reduced efficiency of purifying selection in smallerScatenatuspopulations for most mutations. Our results support an emerging idea that the historical demography of a threatened species has a significant impact on the type of genetic load present, which impacts implementation of conservation actions such as genetic rescue.

     
    more » « less