A linkage map for the hybridising toads Bombina bombina and B. variegata (Anura: Discoglossidae)

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.

Population Genetic Structure Analysis Reveals Decreased but Moderate Diversity for the Oriental Fire-Bellied Toad Introduced to Beijing after 90 Years of Independent Evolution

Simple Summary

Habitat isolation and loss are significant factors that lead to the decline of wildlife populations worldwide, and habitat loss further leads to the shrinkage of populations, which increases the risk of inbreeding and the genetic decline of the populations. To explore the independent evolutionary characteristics of different populations, this study analyzed the genetic disparity of the introduced oriental fire-bellied toad in Beijing from a source population in Shandong Province. The results show that, despite originating from a small artificially introduced population, the toads in the Beijing region have maintained a moderate genetic diversity after 90 years of independent evolution, indicating that this species has a high capacity for survival and adaptation.

Abstract

Detailed molecular genetic research on amphibian populations has a significant role in understanding the genetic adaptability to local environments. The oriental fire-bellied toads (Bombina orientalis) were artificially introduced to Beijing from Shandong Province in 1927, and since then, this separated population went through an independent evolution. To explore the differentiation of the introduced population with its original population, this study analyzed the genetic structure of the oriental fire-bellied toad, based on the mitochondrial genome control region and six microsatellite sites. The results showed that the haplotype diversity and nucleotide diversity of the mitochondrial D-loop region partial sequences of the Beijing Botanical Garden population and the Baiwangshan population were lower than those of the Shangdong Kunyushan population. Microsatellite marker analysis also showed that the observed heterozygosity and expected heterozygosity of the Beijing populations were lower than those of the Kunyushan population. The phylogenetic trees and network diagrams of haplotypes indicated that the three populations were not genetically separated. However, the structure analysis showed a genetic differentiation and categorized the sampling individuals into Beijing and Shandong genetic clusters, which indicated a tendency for isolated evolution in the Beijing population. Although the Beijing populations showed a decline in genetic diversity, it was still at a moderate level, which could maintain the survival of the population. Thus, there is no need to reintroduce new individuals from the Kunyushan source population.

Para-allopatry in hybridizing fire-bellied toads (Bombina bombina and B. variegata): Inference from transcriptome-wide coalescence analyses

Ancient origins, profound ecological divergence, and extensive hybridization make the fire‐bellied toads Bombina bombina and B. variegata (Anura: Bombinatoridae) an intriguing test case of ecological speciation. Previous modeling has proposed that the narrow Bombina hybrid zones represent strong barriers to neutral introgression. We test this prediction by inferring the rate of gene exchange between pure populations on either side of the intensively studied Kraków transect. We developed a method to extract high confidence sets of orthologous genes from de novo transcriptome assemblies, fitted a range of divergence models to these data and assessed their relative support with analytic likelihood calculations. There was clear evidence for postdivergence gene flow, but, as expected, no perceptible signal of recent introgression via the nearby hybrid zone. The analysis of two additional Bombina taxa (B. v. scabra and B. orientalis) validated our parameter estimates against a larger set of prior expectations. Despite substantial cumulative introgression over millions of years, adaptive divergence of the hybridizing taxa is essentially unaffected by their lack of reproductive isolation. Extended distribution ranges also buffer them against small‐scale environmental perturbations that have been shown to reverse the speciation process in other, more recent ecotypes.

All or nothing: Survival, reproduction and oxidative balance in Spotted Wing Drosophila (Drosophila suzukii) in response to cold

Winter severity and overwintering capacity are key ecological factors in successful invasions, especially in ectotherms. The integration of physiological approaches into the study of invasion processes is emerging and promising. Physiological information describes the mechanisms underlying observed survival and reproductive capacities, and it can be used to predict an organism's response to environmental perturbations such as cold temperatures. We investigated the effects of various cold treatments on life history and physiological traits of an invasive pest species, Drosophila suzukii, such as survival, fertility and oxidative balance. This species, a native of temperate Asian areas, is known to survive where cold temperatures are particularly harsh and has been recently introduced into Europe and North America. We found that cold treatments had a strong impact on adult survival but no effect on female's fertility. Although only minor changes were observed after cold treatment on studied physiological traits, a strong sex-based difference was observed in both survival and physiological markers (antioxidant defences and oxidative markers). Females exhibited higher survival, reduced oxidative defences, less damage to nucleic acids, and more damage to lipids. These results suggest that D. suzukii relies on a pathway other than oxidative balance to resist cold injury. Altogether, our results provide information concerning the mechanisms of successful invasion by D. suzukii. These findings may assist in the development of population models that predict the current and future geographic ranges of this species.

First-generation linkage map for the common frog Rana temporaria reveals sex-linkage group

The common frog (Rana temporaria) has become a model species in the fields of ecology and evolutionary biology. However, lack of genomic resources has been limiting utility of this species for detailed evolutionary genetic studies. Using a set of 107 informative microsatellite markers genotyped in a large full-sib family (800 F1 offspring), we created the first linkage map for this species. This partial map-distributed over 15 linkage groups-has a total length of 1698.8 cM. In line with the fact that males are the heterogametic sex in this species and a reduction of recombination is expected, we observed a lower recombination rate in the males (map length: 1371.5 cM) as compared with females (2089.8 cM). Furthermore, three loci previously documented to be sex-linked (that is, carrying male-specific alleles) in adults from the wild mapped to the same linkage group. The linkage map described in this study is one of the densest ones available for amphibians. The discovery of a sex linkage group in Rana temporaria, as well as other regions with strongly reduced male recombination rates, should help to uncover the genetic underpinnings of the sex-determination system in this species. As the number of linkage groups found (n=15) is quite close to the actual number of chromosomes (n=13), the map should provide a useful resource for further evolutionary, ecological and conservation genetic work in this and other closely related species.

Isolation and characterization of the first eight microsatellite loci in Gammarus fossarum (Crustacea, Amphipoda) and cross-amplification in Gammarus pulex and Gammarus orinos

We characterized the first microsatellite markers for Gammarus fossarum. Eight loci gave satisfactory amplification patterns in two stream populations (Southern France) with number of alleles ranging from 2 to 10 and expected heterozygosity from 0.076 to 0.857. We performed cross-amplification in two closely related gammarid species, Gammarus pulex and Gammarus orinos. Among the eight tested microsatellite loci, four correctly amplified in G. pulex and three in G. orinos.

Comparative analysis of mitochondrial genomes in Bombina (Anura; Bombinatoridae)

The complete mitochondrial genomes of two basal anurans, Bombina bombina and B. variegata (Anura; Bombinatoridae), were sequenced. The gene order of their mitochondrial DNA (mtDNA) is identical to that of canonical vertebrate mtDNA. In contrast, we show that there are structural differences in regulatory regions and protein coding genes between the mtDNA of these two closely related species. Corrected sequence divergence between the mtDNA of B. bombina and B. variegata amounts to 8.7% (2.3% divergence in amino acids). Comparisons with two East Asian congeners show that the control region contains two repeat regions, LV1 and LV2, present in all species except for B. bombina, in which LV2 has been secondarily lost. The rRNAs and tRNAs are characterized by low nucleotide divergence. The protein coding genes are considerably more disparate, although functional constraint is high but variable among genes, as evidenced by dN/dS ratios. A mtDNA phylogeny established the distribution of autapomorphic nonsynonomous substitutions in the mitogenomes of B. bombina and B. variegata. Nine of 98 nonsynonomous substitutions led to radical amino acid replacements that may alter mitochondrial protein function. Most radical substitutions were found in ND2, ND4, or ND5, encoding mitochondrial subunits of complex I of the electron transport system. The extensive divergence between the mitogenomes of B. bombina and B. variegata is discussed in terms of its possible role in impeding gene flow in natural hybrid zones between these two species.

Cytonuclear discordance across a leopard frog contact zone

We conducted a genetic analysis of the extent of hybridization within and outside a contact zone between two North American leopard frogs, Rana blairi and Rana pipiens by comparing distribution patterns of mitochondrial and nuclear haplotypes. The contact zone, located between South Dakota, Nebraska, and Iowa, USA, was previously examined using morphological data in the 1970s, leading to the conclusion that hybridization was rare between R. pipiens and R. blairi. Our genetic analysis of 51 populations (611 samples) shows strong cytonuclear discordance. Mitochondrial-haplotype distribution matches the same pattern as the documented species spatial distributions based on morphology. However, the geographic distribution of the nuclear haplotypes reveals asymmetrical swamping of the R. pipiens nuclear haplotypes by R. blairi haplotypes. Phylogenetic analyses of both mitochondrial and nuclear markers provide strong evidence for the presence of R. blairi-R. pipiens introgression for the nuclear marker. A pattern of mitochondrial isolation with nuclear introgression is extremely unusual, and predicted to occur much less often than the reverse.

Genetic architecture of conspecific sperm precedence in Allonemobius fasciatus and A. socius

The evolution of barriers to gene exchange is centrally important to speciation. We used the crickets Allonemobius fasciatus and A. socius to investigate the genetic architecture of conspecific sperm precedence (CSP), a postinsemination prezygotic reproductive barrier. With amplified fragment-length polymorphism (AFLP) markers and controlled crosses we constructed linkage maps and estimated positions of QTL associated with CSP. The majority of QTL have low to moderate effects, although a few QTL exist in A. socius with large effects, and the numbers of QTL are comparable to numbers of genes accounting for species differences in other studies. The QTL are spread across many unlinked markers, yet QTL placed with linked markers are on a small number of linkage groups that could reflect the role of the large Allonemobius sex chromosome in prezygotic isolation. Although many QTL had positive effects on conspecific sperm utilization several QTL also exerted negative effects, which could be explained by intraspecific sexual conflict, sperm competition, or epistasis of introgressed genes on novel backgrounds. One unexpected outcome was that A. socius CSP alleles have a stronger effect than those from A. fasciatus in hybrid females, causing hybrids to behave like A. socius with regard to sperm utilization. Implications of this asymmetry in the Allonemobius hybrid zone are discussed.

Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits

A central challenge in evolutionary biology is to identify genes underlying ecologically important traits and describe the fitness consequences of naturally occurring variation at these loci. To address this goal, several novel approaches have been developed, including 'population genomics,' where a large number of molecular markers are scored in individuals from different environments with the goal of identifying markers showing unusual patterns of variation, potentially due to selection at linked sites. Such approaches are appealing because of (1) the increasing ease of generating large numbers of genetic markers, (2) the ability to scan the genome without measuring phenotypes and (3) the simplicity of sampling individuals without knowledge of their breeding history. Although such approaches are inherently applicable to non-model systems, to date these studies have been limited in their ability to uncover functionally relevant genes. By contrast, quantitative genetics has a rich history, and more recently, quantitative trait locus (QTL) mapping has had some success in identifying genes underlying ecologically relevant variation even in novel systems. QTL mapping, however, requires (1) genetic markers that specifically differentiate parental forms, (2) a focus on a particular measurable phenotype and (3) controlled breeding and maintenance of large numbers of progeny. Here we present current advances and suggest future directions that take advantage of population genomics and quantitative genetic approaches - in both model and non-model systems. Specifically, we discuss advantages and limitations of each method and argue that a combination of the two provides a powerful approach to uncovering the molecular mechanisms responsible for adaptation.

Mating patterns in a hybrid zone of fire-bellied toads (Bombina): inferences from adult and full-sib genotypes

We present two novel methods to infer mating patterns from genetic data. They differ from existing statistical methods of parentage inference in that they apply to populations that deviate from Hardy-Weinberg and linkage equilibrium, and so are suited for the study of assortative mating in hybrid zones. The core data set consists of genotypes at several loci for a number of full-sib clutches of unknown parentage. Our inference is based throughout on estimates of allelic associations within and across loci, such as heterozygote deficit and pairwise linkage disequilibrium. In the first method, the most likely parents of a given clutch are determined from the genotypic distribution of the associated adult population, given an explicit model of nonrandom mating. This leads to estimates of the strength of assortment. The second approach is based solely on the offspring genotypes and relies on the fact that a linear relation exists between associations among the offspring and those in the population of breeding pairs. We apply both methods to a sample from the hybrid zone between the fire-bellied toads Bombina bombina and B. variegata (Anura: Disco glossidae) in Croatia. Consistently, both approaches provide no evidence for a departure from random mating, despite adequate statistical power. Instead, B. variegata-like individuals among the adults contributed disproportionately to the offspring cohort, consistent with their preference for the type of breeding habitat in which this study was conducted.