Rapid selection against inbreeding in a wild population of a rare frog

Variations in Genetic Diversity of Invasive Species Lithobates catesbeianus in China

The introduction and subsequent range expansion of the American bullfrog (Lithobates catesbeianus) is part of a rising trend of troublesome biological invasions happening in China. This detrimental amphibious invasive species has strong adaptability. After its introduction and spread, it established its own ecological niche in many provinces of China, and its range has continued to expand to more areas. Previous studies recorded the introduction time of bullfrogs and calculated the changes in their genetic diversity in China using mitochondria, but the specific introduction route in China is still unknown. Expanding upon previous research, we employed whole-genome scans (utilizing 2b-RAD genomic sequencing) to examine single nucleotide polymorphisms (SNPs) and microsatellites within Lithobates catesbeianus to screen the genomes of these invasive amphibian species from eight Chinese provinces and two U.S. states, including Kansas, where bullfrogs originate. A total of 1,336,475 single nucleotide polymorphic loci and 17 microsatellite loci were used to calculate the genetic diversity of bullfrogs and their migration pathways. Our results suggest that the population in Hunan was the first to be introduced and to spread, and there may have been multiple introductions of subpopulations. Additionally, the genetic diversity of both the SNP and microsatellite loci in the Chinese bullfrog population was lower than that of the US population due to bottleneck effects, but the bullfrogs can adapt and spread rapidly. This study will offer crucial insights for preventing and controlling future introductions into the natural habitats in China. Additionally, it will assist in devising more precise strategies to manage the existing populations and curtail their continued expansion, as well as aim to improve clarity and originality while mitigating plagiarism risk.

A centenary tale: population genetic insights into the introduction history of the oriental fire-bellied toad (Bombina orientalis) in Beijing

Background

The successful establishment of a species population following a single introduction of a few individuals to a non-native area has been limited. Nevertheless, the oriental fire-bellied toad (Bombina orientalis) population in Beijing is purportedly descended from a single introduction of about 200 individuals translocated from Yantai, Shandong Province, China, in 1927.

Results

To resolve the introduction process and to understand the genetic consequences since that introduction approximately 90 years ago, we investigated the population’s genetic diversity and structure using 261 toads from Beijing and two native Shandong populations and inferred the species’ introduction history using simulation-based approaches. Analysis of mitochondrial DNA (mtDNA) sequences showed the two haplotypes found in Beijing nested within Yantai haplotypes, thus corroborating the historical record of the translocation source. The mtDNA and 11 nuclear microsatellite markers revealed both considerably lower genetic diversity in Beijing than in the source population and strong genetic differentiation between them. Although the current census population in Beijing may be in the range of a few thousand, the effective population size was estimated at only 20–57. Simulations also suggest that this population may have descended from 40–60 founders.

Conclusions

The Beijing population’s genetic patterns were consistent with the consequences of a severe bottleneck during introduction followed by genetic drift. The introduction trajectory constructed for this B. orientalis population reveals the genetic footprints of a small population sustained in isolation for nearly a century. Our results provide an intriguing example of establishment success from limited founders and may inform ex situ conservation efforts as well as the management of biological invasions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02072-z.

Annual climatic fluctuations and short-term genetic variation in the eastern spadefoot toad

In addition to variations on the spatial scale, short- and long-term temporal variations, too, can impose intense selection on the overall genetic diversity and composition of a population. We hypothesized that the allelic composition in populations of the eastern spadefoot toad (Pelobates syriacus) would change among successive years in accordance with the short-term changes in environmental conditions. Surprisingly, the effect of short-term climate fluctuations on genetic composition have rarely been addressed in the literature, and to our knowledge the effect of annual climatic fluctuations have not been considered meaningful. Our findings show that climatic variation among successive years, primarily the amount of rainfall and rainy days, can significantly alter both microsatellite allelic composition and diversity. We suggest that environmental (i.e. fluctuating) selection is differential across the globe, and that its intensity is expected to be greatest in regions where short-term climatic conditions are least stable.

Evolutionary principles guiding amphibian conservation

The Anthropocene has witnessed catastrophic amphibian declines across the globe. A multitude of new, primarily human-induced drivers of decline may lead to extinction, but can also push species onto novel evolutionary trajectories. If these are recognized by amphibian biologists, they can be engaged in conservation actions. Here, we summarize how principles stemming from evolutionary concepts have been applied for conservation purposes, and address emerging ideas at the vanguard of amphibian conservation science. In particular, we examine the consequences of increased drift and inbreeding in small populations and their implications for practical conservation. We then review studies of connectivity between populations at the landscape level, which have emphasized the limiting influence of anthropogenic structures and degraded habitat on genetic cohesion. The rapid pace of environmental changes leads to the central question of whether amphibian populations can cope either by adapting to new conditions or by shifting their ranges. We gloomily conclude that extinction seems far more likely than adaptation or range shifts for most species. That said, conservation strategies employing evolutionary principles, such as selective breeding, introduction of adaptive variants through translocations, ecosystem interventions aimed at decreasing phenotype-environment mismatch, or genetic engineering, may effectively counter amphibian decline in some areas or for some species. The spread of invasive species and infectious diseases has often had disastrous consequences, but has also provided some premier examples of rapid evolution with conservation implications. Much can be done in terms of setting aside valuable amphibian habitat that should encompass both natural and agricultural areas, as well as designing protected areas to maximize the phylogenetic and functional diversity of the amphibian community. We conclude that an explicit consideration and application of evolutionary principles, although certainly not a silver bullet, should increase effectiveness of amphibian conservation in both the short and long term.

Invasive crayfish does not influence spawning microhabitat selection of brown frogs

Microhabitat selection is a key component of amphibian breeding biology and can be modulated in response to the features of breeding sites and the presence of predators. Despite invasive alien species being among the major threats to amphibians, there is limited information on the role of invasive species in shaping amphibians’ breeding microhabitat choice. The invasive red swamp crayfish (Procambarus clarkii) is a major predator of amphibians’ larvae, including those of the brown frogs Rana dalmatina and Rana latastei. Although qualitative information about the spawning site preferences and breeding microhabitat choice of brown frogs is available in the literature, only a few studies performed quantitative analyses, and the relationship between microhabitat choice and the presence of alien predators has not been investigated yet. The aims of this study were: (1) to characterize the microhabitats selected for clutch deposition by R. dalmatina and R. latastei and (2) to test if the position and the aggregation of egg clutches differ in sites invaded or not invaded by P. clarkii. During spring 2017, we surveyed multiple times 15 breeding sites of both brown frogs in Northern Italy; in each site we assessed the features of the microhabitat where each egg clutch was laid, considering its position (distance from the shore, depth of the water column) and the degree of aggregation of clutches. In each site we also assessed the presence/absence of the invasive crayfish and the relative abundance in the breeding period. We detected egg clutches in all sites; the crayfish occurred in eight ponds. Our results showed substantial differences between the spawning microhabitat features of the two brown frogs: Rana latastei clutches showed a higher degree of aggregation and were associated with deeper areas of the ponds , while Rana dalmatina deposited more spaced out clutches in areas of the ponds that were less deep. For both species, spawning microhabitat features were not significantly different between sites with and without P. clarkii. Although we did not detect behavioural responses to P. clarkii in the choice of spawning microhabitat , additional studies are required to assess whether these frogs modulate other behavioural traits (e.g. during larval development) in response to the invasive predator.

Temporal trends in genetic data and effective population size support efficacy of management practices in critically endangered dusky gopher frogs (Lithobates sevosus)

Monitoring temporal changes in population genetic diversity and effective population size can provide vital information on future viability. The dusky gopher frog, Lithobates sevosus, is a critically endangered species found only in coastal Mississippi, with low genetic variability as a consequence of isolation and population size reduction. Conservation management practices have been implemented, but their efficacy has not been addressed. We genotyped individuals collected 1997-2014 to determine temporal trends in population genetic variation, structure, and effective size. Observed and expected heterozygosity and allelic richness revealed temporally stable, but low, levels of genetic variation. Positive levels of inbreeding were found in each year. There was weak genetic structure among years, which can be attributed to increased effects of genetic drift and inbreeding in small populations. L. sevosus exhibited an increase in effective population size, and currently has an estimated effective size of 33.0-58.6 individuals, which is approximately half the census size. This large ratio could possibly be explained by genetic compensation. We found that management practices have been effective at maintaining and improving effective size and genetic diversity, but that additional strategies need to be implemented to enhance viability of the species.

Reducing bias in population and landscape genetic inferences: the effects of sampling related individuals and multiple life stages

In population or landscape genetics studies, an unbiased sampling scheme is essential for generating accurate results, but logistics may lead to deviations from the sample design. Such deviations may come in the form of sampling multiple life stages. Presently, it is largely unknown what effect sampling different life stages can have on population or landscape genetic inference, or how mixing life stages can affect the parameters being measured. Additionally, the removal of siblings from a data set is considered best-practice, but direct comparisons of inferences made with and without siblings are limited. In this study, we sampled embryos, larvae, and adult Ambystoma maculatum from five ponds in Missouri, and analyzed them at 15 microsatellite loci. We calculated allelic richness, heterozygosity and effective population sizes for each life stage at each pond and tested for genetic differentiation (F ST and D C ) and isolation-by-distance (IBD) among ponds. We tested for differences in each of these measures between life stages, and in a pooled population of all life stages. All calculations were done with and without sibling pairs to assess the effect of sibling removal. We also assessed the effect of reducing the number of microsatellites used to make inference. No statistically significant differences were found among ponds or life stages for any of the population genetic measures, but patterns of IBD differed among life stages. There was significant IBD when using adult samples, but tests using embryos, larvae, or a combination of the three life stages were not significant. We found that increasing the ratio of larval or embryo samples in the analysis of genetic distance weakened the IBD relationship, and when using D C , the IBD was no longer significant when larvae and embryos exceeded 60% of the population sample. Further, power to detect an IBD relationship was reduced when fewer microsatellites were used in the analysis.

High genetic diversity of common toad (Bufo bufo) populations under strong natural fragmentation on a Northern archipelago

The last decades have shown a surge in studies focusing on the interplay between fragmented habitats, genetic variation, and conservation. In the present study, we consider the case of a temperate pond-breeding anuran (the common toad Bufo bufo) inhabiting a naturally strongly fragmented habitat at the Northern fringe of the species' range: islands offshore the Norwegian coast. A total of 475 individuals from 19 populations (three mainland populations and 16 populations on seven adjacent islands) were genetically characterized using nine microsatellite markers. As expected for a highly fragmented habitat, genetic distances between populations were high (pairwise F st values ranging between 0.06 and 0.33), with however little differences between populations separated by ocean and populations separated by terrestrial habitat (mainland and on islands). Despite a distinct cline in genetic variation from mainland populations to peripheral islands, the study populations were characterized by overall high genetic variation, in line with effective population sizes derived from single-sample estimators which were on average about 20 individuals. Taken together, our results reinforce the notion that spatial and temporal scales of fragmentation need to be considered when studying the interplay between landscape fragmentation and genetic erosion.

Amphibian molecular ecology and how it has informed conservation

Molecular ecology has become one of the key tools in the modern conservationist's kit. Here we review three areas where molecular ecology has been applied to amphibian conservation: genes on landscapes, within-population processes, and genes that matter. We summarize relevant analytical methods, recent important studies from the amphibian literature, and conservation implications for each section. Finally, we include five in-depth examples of how molecular ecology has been successfully applied to specific amphibian systems.

Within- and among-population impact of genetic erosion on adult fitness-related traits in the European tree frog Hyla arborea

Assessing in wild populations how fitness is impacted by inbreeding and genetic drift is a major goal for conservation biology. An approach to measure the detrimental effects of inbreeding on fitness is to estimate correlations between molecular variation and phenotypic performances within and among populations. Our study investigated the effect of individual multilocus heterozygosity on body size, body condition and reproductive investment of males (that is, chorus attendance) and females (that is, clutch mass and egg size) in both small fragmented and large non-fragmented populations of European tree frog (Hyla arborea). Because adult size and/or condition and reproductive investment are usually related, genetic erosion may have detrimental effects directly on reproductive investment, and also on individual body size and condition that in turn may affect reproductive investment. We confirmed that the reproductive investment was highly size-dependent for both sexes. Larger females invested more in offspring production, and larger males attended the chorus in the pond more often. Our results did not provide evidence for a decline in body size, condition and reproductive effort with decreased multilocus heterozygosity both within and among populations. We showed that the lack of heterozygosity-fitness correlations within populations probably resulted from low inbreeding levels (inferior to ca. 20% full-sib mating rate), even in the small fragmented populations. The detrimental effects of fixation load were either low in adults or hidden by environmental variation among populations. These findings will be useful to design specific management actions to improve population persistence.

The female perspective of mating in A. femoralis, a territorial frog with paternal care--a spatial and genetic analysis

The adaptive significance of sequential polyandry is a challenging question in evolutionary and behavioral biology. Costs and benefits of different mating patterns are shaped by the spatial distribution of individuals and by genetic parameters such as the pairwise relatedness between potential mating partners. Thus, females should become less choosy as costs of mating and searching for mates increase. We used parentage assignments to investigate spatial and genetic patterns of mating across a natural population of the Neotropical frog Allobates femoralis, a species characterized by male territoriality and care and female iteroparity. There was no correlation between genetic and spatial distances between adult individuals across the population. In 72% of cases, females mated with males available within a radius of 20 m. Mean pairwise relatedness coefficients of successful reproducers did not differ from random mating but had a lower variance than expected by chance, suggesting maximal reproductive output at intermediate genetic divergence. We also found evidence for selection in favor of more heterozygous individuals between the embryo and adult stage. The level of sequential polyandry significantly increased with the number of spatially available males. Females that had more candidate males also produced more adult progeny. We hypothesize that the benefits associated with female multiple mating outweigh the costs of in- and outbreeding depression, and consequently precluded the evolution of 'choosy' mate selection in this species.

Applications and implications of neutral versus non-neutral markers in molecular ecology

The field of molecular ecology has expanded enormously in the past two decades, largely because of the growing ease with which neutral molecular genetic data can be obtained from virtually any taxonomic group. However, there is also a growing awareness that neutral molecular data can provide only partial insight into parameters such as genetic diversity, local adaptation, evolutionary potential, effective population size, and taxonomic designations. Here we review some of the applications of neutral versus adaptive markers in molecular ecology, discuss some of the advantages that can be obtained by supplementing studies of molecular ecology with data from non-neutral molecular markers, and summarize new methods that are enabling researchers to generate data from genes that are under selection.

Variation in repeat length and heteroplasmy of the mitochondrial DNA control region along a core-edge gradient in the eastern spadefoot toad (Pelobates syriacus)

Peripheral populations are those situated at the distribution margins of a species and are often subjected to more extreme abiotic and biotic conditions than those at the core. Here, we hypothesized that shorter repeat length and fewer heteroplasmic mitochondrial DNA (mtDNA) copies, which are associated with more efficient mitochondrial function, may be related to improved survival under extreme environmental conditions. We sampled eastern spadefoot toads (mostly as tadpoles) from 43 rain pools distributed along a 300-km gradient from core to edge of the species' distribution. We show that mean pool tandem repeat length and heteroplasmy increase from edge to core, even after controlling for body size. We evaluate several alternative hypotheses and propose the Fisher hypothesis as the most likely explanation. However, additional sequential sampling and experimental studies are required to determine whether selection under extreme conditions, or alternative mechanisms, could account for the gradient in heteroplasmy and repeat length in the mtDNA control region.