W. Chris Funk

Harnessing genomics for delineating conservation units

Genomic data have the potential to revolutionize the delineation of conservation units (CUs) by allowing the detection of adaptive genetic variation, which is otherwise difficult for rare, endangered species. In contrast to previous recommendations, we propose that the use of neutral versus adaptive markers should not be viewed as alternatives. Rather, neutral and adaptive markers provide different types of information that should be combined to make optimal management decisions. Genetic patterns at neutral markers reflect the interaction of gene flow and genetic drift that affects genome-wide variation within and among populations. This population genetic structure is what natural selection operates on to cause adaptive divergence. Here, we provide a new framework to integrate data on neutral and adaptive markers to protect biodiversity.

Genetic rescue to the rescue

Genetic rescue can increase the fitness of small, imperiled populations via immigration. A suite of studies from the past decade highlights the value of genetic rescue in increasing population fitness. Nonetheless, genetic rescue has not been widely applied to conserve many of the threatened populations that it could benefit. In this review, we highlight recent studies of genetic rescue and place it in the larger context of theoretical and empirical developments in evolutionary and conservation biology. We also propose directions to help shape future research on genetic rescue. Genetic rescue is a tool that can stem biodiversity loss more than has been appreciated, provides population resilience, and will become increasingly useful if integrated with molecular advances in population genomics.

High levels of cryptic species diversity uncovered in Amazonian frogs

One of the greatest challenges for biodiversity conservation is the poor understanding of species diversity. Molecular methods have dramatically improved our ability to uncover cryptic species, but the magnitude of cryptic diversity remains unknown, particularly in diverse tropical regions such as the Amazon Basin. Uncovering cryptic diversity in amphibians is particularly pressing because amphibians are going extinct globally at an alarming rate. Here, we use an integrative analysis of two independent Amazonian frog clades, Engystomops toadlets and Hypsiboas treefrogs, to test whether species richness is underestimated and, if so, by how much. We sampled intensively in six countries with a focus in Ecuador (Engystomops: 252 individuals from 36 localities; Hypsiboas: 208 individuals from 65 localities) and combined mitochondrial DNA, nuclear DNA, morphological, and bioacoustic data to detect cryptic species. We found that in both clades, species richness was severely underestimated, with more undescribed species than described species. In Engystomops, the two currently recognized species are actually five to seven species (a 150–250% increase in species richness); in Hypsiboas, two recognized species represent six to nine species (a 200–350% increase). Our results suggest that Amazonian frog biodiversity is much more severely underestimated than previously thought.

High dispersal in a frog species suggests that it is vulnerable to habitat fragmentation

Global losses of amphibian populations are a major conservation concern and have generated substantial debate over their causes. Habitat fragmentation is considered one important cause of amphibian decline. However, if fragmentation is to be invoked as a mechanism of amphibian decline, it must first be established that dispersal is prevalent among contiguous amphibian populations using formal movement estimators. In contrast, if dispersal is naturally low in amphibians, fragmentation can be disregarded as a cause of amphibian declines and conservation efforts can be focused elsewhere. We examined dispersal rates in Columbia spotted frogs (Rana luteiventris) using capture–recapture analysis of over 10 000 frogs in combination with genetic analysis of microsatellite loci in replicate basins. We found that frogs had exceptionally high juvenile dispersal rates (up to 62% annually) over long distances (>5 km), large elevation gains (>750 m), and steep inclines (36° incline over 2 km) that were corroborated by genetic data showing high gene flow. These findings show that dispersal is an important life-history feature of some amphibians and suggest that habitat fragmentation is a serious threat to amphibian persistence.

Small effective population size in the long-toed salamander

The effective population sizes (Ne) of six populations of the long‐toed salamander (Ambystoma macrodactylum) from Montana and Idaho, USA were estimated from allozyme data from samples collected in 1978, 1996 and 1997 using the temporal allele frequency method. Five of the six estimates ranged from 23 to 207 (mean = 123 ± 79); one estimate was indistinguishable from infinity. In order to infer the actual Ne of salamander populations, we compared the frequency distribution of our observed Ne estimates with distributions obtained from simulated populations of known Ne. Our observed Ne estimate distribution was consistent with distributions from simulated populations with Ne values of 10, 25, and 50, suggesting an actual Ne for each of the six salamander populations of less than 100. This Ne estimate agrees with most other Ne estimates for amphibians. We conclude by discussing the conservation implications of small Ne values in amphibians in the context of increasing isolation of populations due to habitat fragmentation.

Genetic Differentiation among Long-Toed Salamander (Ambystoma macrodactylum) Populations

Abstract We examined the genetic population structure of Long-Toed Salamanders (Ambystoma macrodactylum) from the Bitterroot Mountains of Idaho and Montana to better understand their evolutionary history and genetic population structure. Populations show high levels of within-population genetic variation at six polymorphic allozyme loci (H̄s = 0.09 for all 18 loci examined; range 0.04–0.14). There is very little divergence among populations within basins, suggesting panmixia within basins. In contrast, genetic differentiation among all populations is high (Gst = 0.30). We used computer simulations to examine population structures that could have led to the observed distribution of genetic variation, assuming selective neutrality of the allozymes. To test the assumption of selective neutrality of the markers used in this study, we compared the observed divergence among the allozymes to that expected from simulations of independently segregating and selectively neutral markers. The observed genetic divergence among populations is compatible with that expected for neutral genetic markers sampled from panmictic populations within basins that exchange less than one migrant among basins each generation.

Unbroken: RADseq remains a powerful tool for understanding the genetics of adaptation in natural populations

Recently, Lowry et al. addressed the ability of RADseq approaches to detect loci under selection in genome scans. While the authors raise important considerations, such as accounting for the extent of linkage disequilibrium in a study system, we strongly disagree with their overall view of the ability of RADseq to inform our understanding of the genetic basis of adaptation. The family of RADseq protocols has radically improved the field of population genomics, expanding by several orders of magnitude the number of markers available while substantially reducing the cost per marker. Researchers whose goal is to identify regions of the genome under selection must consider the LD of the experimental system; however, there is no magical LD cutoff below which researchers should refuse to use RADseq. Lowry et al. further made two major arguments: a theoretical argument that modeled the likelihood of detecting selective sweeps with RAD markers, and gross summaries based on an anecdotal collection of RAD studies. Unfortunately, their simulations were off by two orders of magnitude in the worst case, while their anecdotes merely showed that it is possible to get widely divergent densities of RAD tags for any particular experiment, either by design or due to experimental efficacy. We strongly argue that RADseq remains a powerful and efficient approach that provides sufficient marker density for studying selection in many natural populations. Given limited resources, we argue that researchers should consider a wide range of trade‐offs among genomic techniques, in light of their study question and the power of different techniques to answer it.

COMPARATIVE ANALYSES OF EFFECTIVE POPULATION SIZE WITHIN AND AMONG SPECIES: RANID FROGS AS A CASE STUDY

It has recently become practicable to estimate the effective sizes (Ne) of multiple populations within species. Such efforts are valuable for estimating Ne in evolutionary modeling and conservation planning. We used microsatellite loci to estimate Ne of 90 populations of four ranid frog species (20–26 populations per species, mean n per population = 29). Our objectives were to determine typical values of Ne for populations of each species, compare Ne estimates among the species, and test for correlations between several geographic variables and Ne within species. We used single‐sample linkage disequilibrium (LD), approximate Bayesian computation (ABC), and sibship assignment (SA) methods to estimate contemporary Ne for each population. Three of the species—Rana pretiosa, R. luteiventris, and R. cascadae— have consistently small effective population sizes (<50). Ne in Lithobates pipiens spans a wider range, with some values in the hundreds or thousands. There is a strong east‐to‐west trend of decreasing Ne in L. pipiens. The smaller effective sizes of western populations of this species may be related to habitat fragmentation and population bottlenecking.

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small Ne (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high FST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high FST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small Ne in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.

Range-wide phylogeographic analysis of the spotted frog complex (Rana luteiventris and Rana pretiosa) in northwestern North America

The dynamic geological and climatic history of northwestern North America has made it a focal region for phylogeography. We conducted a range-wide phylogeographic analysis of the spotted frog complex (Rana luteiventris and Rana pretiosa) across its range in northwestern North America to understand its evolutionary history and the distribution of clades to inform conservation of R. pretiosa and Great Basin R. luteiventris, candidates for listing under the US Endangered Species Act. Mitochondrial DNA sequence data from a segment of the cytochrome b gene were obtained from 308 R. luteiventris and R. pretiosa from 96 sites. Phylogenetic analysis revealed one main R. pretiosa clade and three main R. luteiventris clades, two of which overlapped in southeastern Oregon. The three R. luteiventris clades were separated from each other by high levels of sequence divergence (average of 4.75-4.97%). Two divergent clades were also uncovered within the Great Basin. Low genetic variation in R. pretiosa and the southeastern Oregon clade of R. luteiventris suggests concern about their vulnerability to extinction.