Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

Using landscape genomics to delineate future adaptive potential for climate change in the Yosemite toad (Anaxyrus canorus)

An essential goal in conservation biology is delineating population units that maximize the probability of species persisting into the future and adapting to future environmental change. However, future-facing conservation concerns are often addressed using retrospective patterns that could be irrelevant. We recommend a novel landscape genomics framework for delineating future "Geminate Evolutionary Units" (GEUs) in a focal species: (1) identify loci under environmental selection, (2) model and map adaptive conservation units that may spawn future lineages, (3) forecast relative selection pressures on each future lineage, and (4) estimate their fitness and likelihood of persistence using geo-genomic simulations. Using this process, we delineated conservation units for the Yosemite toad (Anaxyrus canorus), a U.S. federally threatened species that is highly vulnerable to climate change. We used a genome-wide dataset, redundancy analysis, and Bayesian association methods to identify 24 candidate loci responding to climatic selection (R 2 ranging from 0.09 to 0.52), after controlling for demographic structure. Candidate loci included genes such as MAP3K5, involved in cellular response to environmental change. We then forecasted future genomic response to climate change using the multivariate machine learning algorithm Gradient Forests. Based on all available evidence, we found three GEUs in Yosemite National Park, reflecting contrasting adaptive optima: YF-North (high winter snowpack with moderate summer rainfall), YF-East (low to moderate snowpack with high summer rainfall), and YF-Low-Elevation (low snowpack and rainfall). Simulations under the RCP 8.5 climate change scenario suggest that the species will decline by 29% over 90 years, but the highly diverse YF-East lineage will be least impacted for two reasons: (1) geographically it will be sheltered from the largest climatic selection pressures, and (2) its standing genetic diversity will promote a faster adaptive response. Our approach provides a comprehensive strategy for protecting imperiled non-model species with genomic data alone and has wide applicability to other declining species.

Landscape genetics of a sub-alpine toad: climate change predicted to induce upward range shifts via asymmetrical migration corridors

Climate change is expected to have a major hydrological impact on the core breeding habitat and migration corridors of many amphibians in the twenty-first century. The Yosemite toad (Anaxyrus canorus) is a species of meadow-specializing amphibian endemic to the high-elevation Sierra Nevada Mountains of California. Despite living entirely on federal lands, it has recently faced severe extirpations, yet our understanding of climatic influences on population connectivity is limited. In this study, we used a previously published double-digest RADseq dataset along with numerous remotely sensed habitat features in a landscape genetics framework to answer two primary questions in Yosemite National Park: (1) Which fine-scale climate, topographic, soil, and vegetation features most facilitate meadow connectivity? (2) How is climate change predicted to influence both the magnitude and net asymmetry of genetic migration? We developed an approach for simultaneously modeling multiple toad migration paths, akin to circuit theory, except raw environmental features can be separately considered. Our workflow identified the most likely migration corridors between meadows and used the unique cubist machine learning approach to fit and forecast environmental models of connectivity. We identified the permuted modeling importance of numerous snowpack-related features, such as runoff and groundwater recharge. Our results highlight the importance of considering phylogeographic structure, and asymmetrical migration in landscape genetics. We predict an upward elevational shift for this already high-elevation species, as measured by the net vector of anticipated genetic movement, and a north-eastward shift in species distribution via the network of genetic migration corridors across the park.

Gene Pool Boundaries for the Yosemite Toad (Anaxyrus canorus) Reveal Asymmetrical Migration Within Meadow Neighborhoods

The Yosemite toad (Anaxyrus [Bufo] canorus) is a federally threatened species of meadow-specializing amphibian endemic to the high-elevation Sierra Nevada Mountains of California. The species is one of the first amphibians to undergo a large demographic collapse that was well-documented, and is reputed to remain in low abundance throughout its range. Recent phylogeographic work has demonstrated that Pleistocene toad lineages diverged and then admixed to differing extents across an elevational gradient. Although lineage divisions may have significant effects on evolutionary trajectories over large spatial and temporal scales, present-day population dynamics must be delineated in order to manage and conserve the species effectively. In this study, we used a double-digest RADseq dataset to address three primary questions: (1) Are single meadows or neighborhoods of nearby meadows most correlated with population boundaries? (2) Does asymmetrical migration occur among neighborhoods of nearby meadows? (3) What topographic or hydrological variables predict such asymmetrical migration in these meadow neighborhoods? Hierarchical STRUCTURE and AMOVA analyses suggested that populations are typically circumscribed by a single meadow, although 84% of meadows exist in neighborhoods of at least two meadows connected by low levels of migration, and over half (53%) of neighborhoods examined display strong asymmetrical migration. Meadow neighborhoods often contain one or more large and flat “hub” meadows that experience net immigration, surrounded by smaller and topographically rugged “satellite” meadows with net emigration. Hubs tend to contain more genetic diversity and could be prioritized for conservation and habitat management and as potential sources for reestablishment efforts.

Phylogenetics and biogeography of the two-wing flyingfish (Exocoetidae: Exocoetus)

Two-wing flyingfish (Exocoetus spp.) are widely distributed, epipelagic, mid-trophic organisms that feed on zooplankton and are preyed upon by numerous predators (e.g., tunas, dolphinfish, tropical seabirds), yet an understanding of their speciation and systematics is lacking. As a model of epipelagic fish speciation and to investigate mechanisms that increase biodiversity, we studied the phylogeny and biogeography of Exocoetus, a highly abundant holoepipelagic fish taxon of the tropical open ocean. Morphological and molecular data were used to evaluate the phylogenetic relationships, species boundaries, and biogeographic patterns of the five putative Exocoetus species. We show that the most widespread species (E. volitans) is sister to all other species, and we find no evidence for cryptic species in this taxon. Sister relationship between E. monocirrhus (Indo-Pacific) and E. obtusirostris (Atlantic) indicates the Isthmus of Panama and/or Benguela Barrier may have played a role in their divergence via allopatric speciation. The sister species E. peruvianus and E. gibbosus are found in different regions of the Pacific Ocean; however, our molecular results do not show a clear distinction between these species, indicating recent divergence or ongoing gene flow. Overall, our phylogeny reveals that the most spatially restricted species are more recently derived, suggesting that allopatric barriers may drive speciation, but subsequent dispersal and range expansion may affect the distributions of species.

Population Genetic Structure of the Tropical Two-Wing Flyingfish (Exocoetus volitans)

Delineating populations of pantropical marine fish is a difficult process, due to widespread geographic ranges and complex life history traits in most species. Exocoetus volitans, a species of two-winged flyingfish, is a good model for understanding large-scale patterns of epipelagic fish population structure because it has a circumtropical geographic range and completes its entire life cycle in the epipelagic zone. Buoyant pelagic eggs should dictate high local dispersal capacity in this species, although a brief larval phase, small body size, and short lifespan may limit the dispersal of individuals over large spatial scales. Based on these biological features, we hypothesized that E. volitans would exhibit statistically and biologically significant population structure defined by recognized oceanographic barriers. We tested this hypothesis by analyzing cytochrome b mtDNA sequence data (1106 bps) from specimens collected in the Pacific, Atlantic and Indian oceans (n = 266). AMOVA, Bayesian, and coalescent analytical approaches were used to assess and interpret population-level genetic variability. A parsimony-based haplotype network did not reveal population subdivision among ocean basins, but AMOVA revealed limited, statistically significant population structure between the Pacific and Atlantic Oceans (ΦST = 0.035, p<0.001). A spatially-unbiased Bayesian approach identified two circumtropical population clusters north and south of the Equator (ΦST = 0.026, p<0.001), a previously unknown dispersal barrier for an epipelagic fish. Bayesian demographic modeling suggested the effective population size of this species increased by at least an order of magnitude ~150,000 years ago, to more than 1 billion individuals currently. Thus, high levels of genetic similarity observed in E. volitans can be explained by high rates of gene flow, a dramatic and recent population expansion, as well as extensive and consistent dispersal throughout the geographic range of the species.

Genetic and experimental evidence for a mixed-age, mixed-origin bank of kelp microscopic stages in southern California

Laboratory studies have demonstrated that the microscopic stages of kelps can rapidly resume development from a delayed state. Like terrestrial seeds or aquatic resting eggs, banks of delayed kelp stages may supplement population recovery after periods of stress, playing an important role for kelp populations that experience adult sporophyte absences due to seasonal or interannual disturbances. We found that removing the microscopic stages from natural rock substratum could prevent the appearance of juvenile kelp sporophytes for three months and the establishment of a diverse kelp assemblage for over four months within a southern California kelp forest. Juveniles were observed within one month in plots where microscopic stages were left intact, which may confer an advantage for the resulting sporophytes as they attain larger sizes before later recruiting neighbors. Microsatellite diversity was high (expected heterozygosity HE approximately 0.9) for juveniles and adults within our sites. Using a microsatellite-based parentage analysis for the dominant kelp, Macrocystis pyrifera, we estimated that a portion of the new M. pyrifera sporophyte recruits had originated from their parents at least seven months after their parents had disappeared. Similar delay durations have been demonstrated in recent laboratory studies. Additionally, our results suggest that zoospore dispersal distances > 50 m may be supported by including additional microsatellite loci in the analysis. We propose a mixed-age and, potentially, a mixed-origin bank of M. pyrifera gametophytes promotes maximal genetic diversity in recovering populations and reduces population genetic subdivision and self-fertilization rates for intact populations by promoting the survival of zoospores dispersed > 10 m and during inhospitable environmental conditions.

Genetic signature of rapid IHHNV (infectious hypodermal and hematopoietic necrosis virus) expansion in wild Penaeus shrimp populations

Infectious hypodermal and hematopoietic necrosis virus (IHHNV) is a widely distributed single-stranded DNA parvovirus that has been responsible for major losses in wild and farmed penaeid shrimp populations on the northwestern Pacific coast of Mexico since the early 1990's. IHHNV has been considered a slow-evolving, stable virus because shrimp populations in this region have recovered to pre-epizootic levels, and limited nucleotide variation has been found in a small number of IHHNV isolates studied from this region. To gain insight into IHHNV evolutionary and population dynamics, we analyzed IHHNV capsid protein gene sequences from 89 Penaeus shrimp, along with 14 previously published sequences. Using Bayesian coalescent approaches, we calculated a mean rate of nucleotide substitution for IHHNV that was unexpectedly high (1.39×10⁻⁴ substitutions/site/year) and comparable to that reported for RNA viruses. We found more genetic diversity than previously reported for IHHNV isolates and highly significant subdivision among the viral populations in Mexican waters. Past changes in effective number of infections that we infer from Bayesian skyline plots closely correspond to IHHNV epizootiological historical records. Given the high evolutionary rate and the observed regional isolation of IHHNV in shrimp populations in the Gulf of California, we suggest regular monitoring of wild and farmed shrimp and restriction of shrimp movement as preventative measures for future viral outbreaks.

Parallelization and optimization of genetic analyses in isolation by distance web service

Background

The Isolation by Distance Web Service (IBDWS) is a user-friendly web interface for analyzing patterns of isolation by distance in population genetic data. IBDWS enables researchers to perform a variety of statistical tests such as Mantel tests and reduced major axis regression (RMA), and returns vector based graphs. The more than 60 citations since 2005 confirm the popularity and utility of this website. Despite its usefulness, the data sets with over 65 populations can take hours or days to complete due to the computational intensity of the statistical tests. This is especially troublesome for web-based software analysis, since users tend to expect real-time results on the order of seconds, or at most, minutes. Moreover, as genetic data continue to increase and diversify, so does the demand for more processing power. In order to increase the speed and efficiency of IBDWS, we first determined which aspects of the code were most time consuming and whether they might be amenable to improvements by parallelization or algorithmic optimization.

Results

Runtime tests uncovered two areas of IBDWS that consumed significant amounts of time: randomizations within the Mantel test and the RMA calculations. We found that these sections of code could be restructured and parallelized to improve efficiency. The code was first optimized by combining two similar randomization routines, implementing a Fisher-Yates shuffling algorithm, and then parallelizing those routines. Tests of the parallelization and Fisher-Yates algorithmic improvements were performed on a variety of data sets ranging from 10 to 150 populations. All tested algorithms showed runtime reductions and a very close fit to the predicted speedups based on time-complexity calculations. In the case of 150 populations with 10,000 randomizations, data were analyzed 23 times faster.

Conclusion

Since the implementation of the new algorithms in late 2007, datasets have continued to increase substantially in size and many exceed the largest population sizes we used in our test sets. The fact that the website has continued to work well in "real-world" tests, and receives a considerable number of new citations provides the strongest testimony to the effectiveness of our improvements. However, we soon expect the need to upgrade the number of nodes in our cluster significantly as dataset sizes continue to expand. The parallel implementation can be found at .

Identification of co-occurring Branchinecta fairy shrimp species from encysted embryos using multiplex polymerase chain reaction

Morphological identification of many fairy shrimp species is difficult because distinguishing characters are restricted to adults. We developed two multiplex polymerase chain reaction assays that differentiate among three Branchinecta fairy shrimp with distributional overlap in southern California vernal pools. Two of the species are federally listed as threatened. Molecular identification of Branchinecta from cysts allows for species surveys to be conducted during the dry season, expanding the timeframe for population assessment and providing a less intrusive method of sampling sensitive vernal pool habitats.

Understanding the genetic effects of recent habitat fragmentation in the context of evolutionary history: phylogeography and landscape genetics of a southern California endemic Jerusalem cricket (Orthoptera: Stenopelmatidae: Stenopelmatus)

Habitat loss and fragmentation due to urbanization are the most pervasive threats to biodiversity in southern California. Loss of habitat and fragmentation can lower migration rates and genetic connectivity among remaining populations of native species, reducing genetic variability and increasing extinction risk. However, it may be difficult to separate the effects of recent anthropogenic fragmentation from the genetic signature of prehistoric fragmentation due to previous natural geological and climatic changes. To address these challenges, we examined the phylogenetic and population genetic structure of a flightless insect endemic to cismontane southern California, Stenopelmatus'mahogani' (Orthoptera: Stenopelmatidae). Analyses of mitochondrial DNA sequence data suggest that diversification across southern California began during the Pleistocene, with most haplotypes currently restricted to a single population. Patterns of genetic divergence correlate with contemporary urbanization, even after correcting for (geographical information system) GIS-based reconstructions of fragmentation during the Pleistocene. Theoretical simulations confirm that contemporary patterns of genetic structure could be produced by recent urban fragmentation using biologically reasonable assumptions about model parameters. Diversity within populations was positively correlated with current fragment size, but not prehistoric fragment size, suggesting that the effects of increased drift following anthropogenic fragmentation are already being seen. Loss of genetic connectivity and diversity can hinder a population's ability to adapt to ecological perturbations commonly associated with urbanization, such as habitat degradation, climatic changes and introduced species. Consequently, our results underscore the importance of preserving and restoring landscape connectivity for long-term persistence of low vagility native species.

Population structure of the pumpkin fruit flyBactrocera depressa(Tephritidae) in Korea and Japan: Pliocene allopatry or recent invasion?

Because of their widespread agricultural impact and rapid range expansions, true fruit flies (Tephritidae) are the subject of quarantine and control efforts worldwide. Among these flies, the pumpkin fruit fly Bactrocera depressa, which infests squash and other cucurbitaceous plants in Korea, Japan and Taiwan, was recently isolated from produce shipments entering Japan and identified as a regulatory target. This species was described in 1933 from collections in Japan and discovered in 1974 in Korea, suggesting that it may have recently invaded mainland Asia. We analysed the genetic structure of Asian populations of B. depressa using sequence variation for mitochondrial gene cytochrome-oxidase I and three nuclear loci: elongation factor 1alpha, tubulinbeta1 and tubulinbeta3, using frequency-based approaches, nested clade analysis and assignment tests. Contrary to the hypothesis of recent invasion, high levels of genetic subdivision were found among five Korean and three Japanese populations. Nested clade analysis suggested a variety of processes operating over different time scales, including ancient isolation between Korea and Japan and more recent range expansions within each country. Contrary to a priori expectations, the results also suggested the recent introduction of a mitochondrial haplotype into Yokohama, Japan that is related closely to a widespread haplotype found throughout Korea. Assignment tests also supported these conclusions. The combination of a genealogical approach and probabilistic assignments of individuals to populations of origin was able to provide statistical support for the identification of cryptic introductions within an otherwise widespread indigenous species.

Dispersal, gene flow, and population structure

The accuracy of gene flow estimates is unknown in most natural populations because direct estimates of dispersal are often not possible. These estimates can be highly imprecise or even biased because population genetic structure reflects more than a simple balance between genetic drift and gene flow. Most of the models used to estimate gene flow also assume very simple patterns of movement. As a result, multiple interpretations of population structure involving contemporary gene flow, departures from equilibrium, and other factors are almost always possible. One way to isolate the relative contribution of gene flow to population genetic differentiation is to utilize comparative methods. Population genetic statistics such as FST, heterozygosity and Nei's D can be compared between species with differing dispersal abilities if these species are otherwise phylogenetically, geographically and demographically comparable. Accordingly, the available literature was searched for all groups that meet these criteria to determine whether broad conclusions regarding the relationships between dispersal, population genetic structure, and gene flow estimates are possible. Allozyme and mtDNA data were summarized for 27 animal groups in which dispersal differences can be characterized. In total, genetic data were obtained for 333 species of vertebrates and invertebrates from terrestrial, freshwater and marine habitats. Across these groups, dispersal ability was consistently related to population structure, with a mean rank correlation of -0.72 between ranked dispersal ability and FST. Gene flow estimates derived from private alleles were also correlated with dispersal ability, but were less widely available. Direct-count heterozygosity and average values of Nei's D showed moderate degrees of correlation with dispersal ability. Thus, despite regional, taxonomic and methodological differences among the groups of species surveyed, available data demonstrate that dispersal makes a measurable contribution to population genetic differentiation in the majority of animal species in nature, and that gene flow estimates are rarely so overwhelmed by population history, departures from equilibrium, or other microevolutionary forces as to be uninformative.

Genetic population structure of the fairy shrimp Branchinecta coloradensis (Anostraca) in the Rocky Mountains of Colorado

Dispersal rates for freshwater invertebrates are often inferred from population genetic data. Although genetic approaches can indicate the amount of isolation in natural populations, departures from an equilibrium between drift and gene flow often lead to biased gene flow estimates. I investigated the genetic population structure of the pond-dwelling fairy shrimp Branchinecta coloradensis in the Rocky Mountains of Colorado, U.S.A., using allozymes. Glaciation in this area and the availability of direct dispersal estimates from previous work permit inferences regarding the relative impacts of history and contemporary gene flow on population structure. Hierarchical F statistics were used to quantify differentiation within and between valleys (thetaSV and thetaVT, respectively). Between valleys separated by 5-10 km, a high degree of differentiation (thetaVT = 0.77) corresponds to biologically reasonable gene flow estimates of 0.07 individuals per generation, although it is possible that this value represents founder effects and nonequilibrium conditions. On a local scale (<=110 m), populations are genetically similar (thetaSV = 0.13) and gene flow is estimated to be 1.7 individuals exchanged between ponds each generation. This is very close to an ecological estimate of dispersal for B. coloradensis via salamanders. Gene flow estimates from previous studies on other Anostraca are also similar on comparable geographic scales. Thus, population structure in B. coloradensis appears to be at or near equilibrium on a local scale, and possibly on a regional scale as well.

Seasonal movement patterns in a subalpine population of the tiger salamander, Ambystoma tigrinum nebulosum

Seasonal movements of a subalpine population of metamorphic tiger salamanders, Ambystoma tigrinum nebulosum, were explored from 1990 to 1992. Metamorphic adults bred in permanent and semipermanent habitats during June of each year. After breeding, soma individuals returned to the terrestrial environment, but many remained and congregated in nonpermanent ponds. Metamorphic adults in nonpermanent ponds included those that moved from permanent ponds, those that bred and remained in semipermanent ponds, and those that did not breed but migrated from the terrestrial environment after the breeding season. Dietary analyses indicated that metamorphic adults in nonpermanent ponds have significantly greater numbers, biomass, and calories of prey in their gut than the few metamorphic adults remaining in permanent ponds. This difference was due to the presence of fairy shrimp (Branchinecta coloradensis), which composed 91 % of metamorphic adult diets in nonpermanent ponds and accounted for three times the average caloric value of all prey in gut samples from metamorphic adults in permanent ponds. Intraspecific competition also may have contributed to movement into nonpermanent ponds: post-breeding densities of metamorphic adults in permanent ponds were inversely related to the densities of other morphs (paedo-morphic adults and large larvae). We conclude that metamorphic A. t. nebulosum in this population utilize nonpermanent ponds after the breeding season because of the abundance of high-quality aquatic prey and reduced competition from con-specific morphs.