History vs. current demography: explaining the genetic population structure of the common frog (Rana temporaria)

Major Histocompatibility Complex Variation and Haplotype Associated Survival in Response to Experimental Infection of Two Bd-GPL Strains Along a Latitudinal Gradient

While both innate and adaptive immune system mechanisms have been implicated in resistance against the chytrid fungus Batrachochytrium dendrobatidis (Bd), studies on the role of specific MHC haplotypes on Bd infection are rare. Here, we studied variation in MHC Class IIB loci in the common toad Bufo bufo along a latitudinal gradient across Sweden. In general, Swedish toad populations had few MHC Class IIB haplotypes and MHC diversity declined from south (13 haplotypes) to the north (four haplotypes). The low diversity may compromise the ability of northern populations to fight emerging disease, such as Bd. In a laboratory experiment, we infected newly metamorphosed toads with two strains of the Global Pandemic Lineage of the fungus (Bd-GPL) and compared survival with sham controls. Bd-infected toads had lower survival compared to controls. Moreover, survival was dependent on the Bd-strain and northern toads had lower Bd-mediated survival than southern individuals. MHC diversity was lower in northern toads. All northern experimental animals were monomorphic for a single MHC haplotype, whereas we found seven different haplotypes in southern experimental animals. In southern toads, survival was dependent on both Bd-strain and MHC haplotype suggesting differential infection dynamics depending on both Bd-strain and host immune system characteristics.

Antimicrobial peptide and sequence variation along a latitudinal gradient in two anurans

BACKGROUND

While there is evidence of both purifying and balancing selection in immune defense genes, large-scale genetic diversity in antimicrobial peptides (AMPs), an important part of the innate immune system released from dermal glands in the skin, has remained uninvestigated. Here we describe genetic diversity at three AMP loci (Temporin, Brevinin and Palustrin) in two ranid frogs (Rana arvalis and R. temporaria) along a 2000 km latitudinal gradient. We amplified and sequenced part of the Acidic Propiece domain and the hypervariable Mature Peptide domain (~ 150-200 bp) in the three genes using Illumina Miseq and expected to find decreased AMP genetic variation towards the northern distribution limit of the species similarly to studies on MHC genetic patterns.

RESULTS

We found multiple loci for each AMP and relatively high gene diversity, but no clear pattern of geographic genetic structure along the latitudinal gradient. We found evidence of trans-specific polymorphism in the two species, indicating a common evolutionary origin of the alleles. Temporin and Brevinin did not form monophyletic clades suggesting that they belong to the same gene family. By implementing codon evolution models we found evidence of strong positive selection acting on the Mature Peptide. We also found evidence of diversifying selection as indicated by divergent allele frequencies among populations and high Theta k values.

CONCLUSION

Our results suggest that AMPs are an important source of adaptive diversity, minimizing the chance of microorganisms developing resistance to individual peptides.

Genetic diversity and structure of two endangered mole salamander species of the Trans-Mexican Volcanic Belt

The most important factor leading to amphibian population declines and extinctions is habitat degradation and destruction. To help prevent further extinctions, studies are needed to make appropriate conservation decisions in small and fragmented populations. The goal of this study was to provide data from the population genetics of two micro-endemic mole salamanders from the Trans-Mexican Volcanic Belt. Nine microsatellite markers were used to study the population genetics of 152 individuals from twoAmbystomaspecies. We sampled 38 individuals in two localities forA. altamiraniandA. rivualre. We found medium to high levels of genetic diversity expressed as heterozygosity in the populations. However, all the populations presented few alleles per locus and genotypes. We found strong genetic structure between populations for each species. Effective population size was small but similar to that of the studies from other mole salamanders with restricted distributions or with recently fragmented habitats. Despite the medium to high levels of genetic diversity expressed as heterozygosity, we found few alleles, evidence of a genetic bottleneck and that the effective population size is small in all populations. Therefore, this study is important to propose better management plans and conservation efforts for these species.

Genetic diversity and demography of the critically endangered Roberts' false brook salamander (Pseudoeurycea robertsi) in Central Mexico

Land use changes are threatening the maintenance of biodiversity. Genetic diversity is one of the main indicators of biological diversity and is highly important as it shapes the capability of populations to respond to environmental changes. We studied eleven populations of Pseudoeurycea robertsi, a micro-endemic and critically endangered species from the Nevado de Toluca Volcano, a mountain that is part of the Trans-Mexican Volcanic Belt, Mexico. We sequenced the mitochondrial cytochrome b gene from 71 individuals and genotyped 9 microsatellites from 150 individuals. Our results based on the cytochrome b showed two divergent lineages, with moderate levels of genetic diversity and a recently historical demographic expansion. Microsatellite-based results indicated low levels of heterozygosity for all populations and few alleles per locus, as compared with other mole salamander species. We identified two genetically differentiated subpopulations with a significant level of genetic structure. These results provide fundamental data for the development of management plans and conservation efforts for this critically endangered species.

Sampling related individuals within ponds biases estimates of population structure in a pond-breeding amphibian

Effective conservation and management of pond-breeding amphibians depends on the accurate estimation of population structure, demographic parameters, and the influence of landscape features on breeding-site connectivity. Population-level studies of pond-breeding amphibians typically sample larval life stages because they are easily captured and can be sampled nondestructively. These studies often identify high levels of relatedness between individuals from the same pond, which can be exacerbated by sampling the larval stage. Yet, the effect of these related individuals on population genetic studies using genomic data is not yet fully understood. Here, we assess the effect of within-pond relatedness on population and landscape genetic analyses by focusing on the barred tiger salamanders (Ambystoma mavortium) from the Nebraska Sandhills. Utilizing genome-wide SNPs generated using a double-digest RADseq approach, we conducted standard population and landscape genetic analyses using datasets with and without siblings. We found that reduced sample sizes influenced parameter estimates more than the inclusion of siblings, but that within-pond relatedness led to the inference of spurious population structure when analyses depended on allele frequencies. Our landscape genetic analyses also supported different models across datasets depending on the spatial resolution analyzed. We recommend that future studies not only test for relatedness among larval samples but also remove siblings before conducting population or landscape genetic analyses. We also recommend alternative sampling strategies to reduce sampling siblings before sequencing takes place. Biases introduced by unknowingly including siblings can have significant implications for population and landscape genetic analyses, and in turn, for species conservation strategies and outcomes.

Deciphering the drivers of negative species-genetic diversity correlation in Alpine amphibians

The evolutionary and ecological importance of neutral and adaptive genetic diversity is widely recognized. Nevertheless, genetic diversity is rarely assessed for conservation planning, which often implicitly assumes a positive correlation between species and genetic diversity. Multiple drivers can cause the co-variation between the genetic diversity of one species and the richness of the whole communities, and explicit tests are needed to identify the processes that can determine species-genetic diversity correlations (SGDCs). Here, we tested whether intrapopulation genetic diversity (at neutral loci) and species richness co-vary in the amphibian communities of a southern Alpine region (Trentino, Italy), using the common frog (Rana temporaria) as focal species for the study of genetic diversity. We also analysed ecological similarity, niche overlap and interspecific interactions between the species, to unravel the processes determining SGDC. The neutral genetic diversity of common frogs was negatively related to species richness. The negative SGDC was probably due to an opposite influence of environmental gradients on the two levels of biodiversity, since the focal species and the other amphibians differ in ecological preferences, particularly in terms of thermal optimum. Conversely, we did not find evidence for a role of interspecific interactions in the negative SGDC. Our findings stress that species richness cannot be used as a universal proxy for genetic diversity, and only combining SGDC with analyses on the determinants of biodiversity can allow to identify the processes determining the relationships between genetic and species diversity.

Present and future ecological niche modeling of garter snake species from the Trans-Mexican Volcanic Belt

Land use and climate change are affecting the abundance and distribution of species. The Trans-Mexican Volcanic Belt (TMVB) is a very diverse region due to geological history, geographic position, and climate. It is also one of the most disturbed regions in Mexico. Reptiles are particularly sensitive to environmental changes due to their low dispersal capacity and thermal ecology. In this study, we define the important environmental variables (considering climate, topography, and land use) and potential distribution (present and future) of the five Thamnophis species present in TMVB. To do so, we used the maximum entropy modeling software (MAXENT). First, we modeled to select the most important variables to explain the distribution of each species, then we modeled again using only the most important variables and projected these models to the future considering a middle-moderate climate change scenario (rcp45), and land use and vegetation variables for the year 2050 (generated according to land use changes that occurred between years 2002 and 2011). Arid vegetation had an important negative effect on habitat suitability for all species, and minimum temperature of the coldest month was important for four of the five species. Thamnophis cyrtopsis was the species with the lowest tolerance to minimum temperatures. The maximum temperature of the warmest month was important for T. scalaris and T. cyrtopsis. Low percentages of agriculture were positive for T. eques and T. melanogaster but, at higher values, agriculture had a negative effect on habitat suitability for both species. Elevation was the most important variable to explain T. eques and T. melanogaster potential distribution while distance to Abies forests was the most important variable for T. scalaris and T. scaliger. All species had a high proportion of their potential distribution in the TMVB. However, according to our models, all Thamnophis species will experience reductions in their potential distribution in this region. T. scalaris will suffer the biggest reduction because this species is limited by high temperatures and will not be able to shift its distribution upward, as it is already present in the highest elevations of the TMVB.

Genetic variability and structure of an isolated population of Ambystoma altamirani, a mole salamander that lives in the mountains of one of the largest urban areas in the world

Amphibians are globally threatened by habitat loss and fragmentation; species within the order Ambystoma are not the exception, as there are 18 species of mole salamanders in México, of which 16 are endemic and all species are under some national or international status of protection. The mole salamander, Ambystoma altamirani is a microendemic species, which is distributed in central México, within the trans-Mexican volcanic belt, and is one of the most threatened species due to habitat destruction and the introduction of exotic species. Nine microsatellite markers were used to determine the genetic structure, genetic variability, effective population size, presence of bottlenecks and inbreeding coefficient of one population of A. altamirani to generate information which might help to protect and conserve this threatened species. We found two genetic subpopulations with significant level of genetic structure (F_ST = 0.005) and high levels of genetic variability (H_o = 0.883; H_e = 0.621); we also found a small population size (N_e = 8.8), the presence of historical (M = 0.486) and recent bottlenecks under IAM and TPM models, with a low, but significant coefficient of inbreeding (F_IS = -0.451). This information will help us to raise conservation strategies of this microendemic mole salamander species.

Drift, selection, or migration? Processes affecting genetic differentiation and variation along a latitudinal gradient in an amphibian

BACKGROUND

Past events like fluctuations in population size and post-glacial colonization processes may influence the relative importance of genetic drift, migration and selection when determining the present day patterns of genetic variation. We disentangle how drift, selection and migration shape neutral and adaptive genetic variation in 12 moor frog populations along a 1700 km latitudinal gradient. We studied genetic differentiation and variation at a MHC exon II locus and a set of 18 microsatellites.

RESULTS

Using outlier analyses, we identified the MHC II exon 2 (corresponding to the β-2 domain) locus and one microsatellite locus (RCO8640) to be subject to diversifying selection, while five microsatellite loci showed signals of stabilizing selection among populations. STRUCTURE and DAPC analyses on the neutral microsatellites assigned populations to a northern and a southern cluster, reflecting two different post-glacial colonization routes found in previous studies. Genetic variation overall was lower in the northern cluster. The signature of selection on MHC exon II was weaker in the northern cluster, possibly as a consequence of smaller and more fragmented populations.

CONCLUSION

Our results show that historical demographic processes combined with selection and drift have led to a complex pattern of differentiation along the gradient where some loci are more divergent among populations than predicted from drift expectations due to diversifying selection, while other loci are more uniform among populations due to stabilizing selection. Importantly, both overall and MHC genetic variation are lower at northern latitudes. Due to lower evolutionary potential, the low genetic variation in northern populations may increase the risk of extinction when confronted with emerging pathogens and climate change.

Population differentiation at a regional scale in spadefoot toads: contributions of distance and divergent selective environments

The causes of population differentiation can provide insight into the origins of early barriers to gene flow. Two key drivers of population differentiation are geographic distance and local adaptation to divergent selective environments. When reproductive isolation arises because some populations of a species are under selection to avoid hybridization while others are not, population differentiation and even speciation can result. Spadefoot toad populations Spea multiplicata that are sympatric with a congener have undergone reinforcement. This reinforcement has resulted not only in increased reproductive isolation from the congener, but also in the evolution of reproductive isolation from nearby and distant conspecific allopatric populations. We used multiple approaches to evaluate the contributions of geographic distance and divergent selective environments to population structure across this regional scale in S. multiplicata, based on genotypes from six nuclear microsatellite markers. We compared groups of populations varying in both geographic location and in the presence of a congener. Hierarchical F-statistics and results from cluster analyses and discriminant analyses of principal components all indicate that geographic distance is the stronger contributor to genetic differentiation among S. multiplicata populations at a regional scale. However, we found evidence that adaptation to divergent selective environments also contributes to population structure. Our findings highlight how variation in the balance of evolutionary forces acting across a species' range can lead to variation in the relative contributions of geographic distance and local adaptation to population differentiation across different spatial scales.

Ecological opportunities and specializations shaped genetic divergence in a highly mobile marine top predator

Environmental conditions can shape genetic and morphological divergence. Release of new habitats during historical environmental changes was a major driver of evolutionary diversification. Here, forces shaping population structure and ecotype differentiation ('pelagic' and 'coastal') of bottlenose dolphins in the North-east Atlantic were investigated using complementary evolutionary and ecological approaches. Inference of population demographic history using approximate Bayesian computation indicated that coastal populations were likely founded by the Atlantic pelagic population after the Last Glacial Maxima probably as a result of newly available coastal ecological niches. Pelagic dolphins from the Atlantic and the Mediterranean Sea likely diverged during a period of high productivity in the Mediterranean Sea. Genetic differentiation between coastal and pelagic ecotypes may be maintained by niche specializations, as indicated by stable isotope and stomach content analyses, and social behaviour. The two ecotypes were only weakly morphologically segregated in contrast to other parts of the World Ocean. This may be linked to weak contrasts between coastal and pelagic habitats and/or a relatively recent divergence. We suggest that ecological opportunity to specialize is a major driver of genetic and morphological divergence. Combining genetic, ecological and morphological approaches is essential to understanding the population structure of mobile and cryptic species.

Small woods positively influence the occurrence and abundance of the common frog (Rana temporaria) in a traditional farming landscape

Traditional farming landscape in western Europe is made of a complex mosaic of pastures, cultures, ponds and hedgerows connected with woods. Previous observations in the common frog species suggest that lowland populations are closely associated to wood cover and our aim was to test the validity of this assumption. We studied common frog occurrence and abundance in western central France (Deux-Sèvres department) close to the southern margin of lowland distribution. Our results pointed out that the proportion of woods surface around sampled areas (1 ha) was a critical determinant of common frog presence and abundance. Extensive farming, which maintains a mosaic of small woods, may provide a robust conservation tool for this species.

Physiological vagility and its relationship to dispersal and neutral genetic heterogeneity in vertebrates

Vagility is the inherent power of movement by individuals. Vagility and the available duration of movement determine the dispersal distance individuals can move to interbreed, which affects the fine-scale genetic structure of vertebrate populations. Vagility and variation in population genetic structure are normally explained by geographic variation and not by the inherent power of movement by individuals. We present a new, quantitative definition for physiological vagility that incorporates aerobic capacity, body size, body temperature and the metabolic cost of transport, variables that are independent of the physical environment. Physiological vagility is the speed at which an animal can move sustainably based on these parameters. This meta-analysis tests whether this definition of physiological vagility correlates with empirical data for maximal dispersal distances and measured microsatellite genetic differentiation with distance {[F(ST)/[1-F(ST))]/ln distance} for amphibians, reptiles, birds and mammals utilizing three locomotor modes (running, flying, swimming). Maximal dispersal distance and physiological vagility increased with body mass for amphibians, reptiles and mammals utilizing terrestrial movement. The relative slopes of these relationships indicate that larger individuals require longer movement durations to achieve maximal dispersal distances. Both physiological vagility and maximal dispersal distance were independent of body mass for flying vertebrates. Genetic differentiation with distance was greatest for terrestrial locomotion, with amphibians showing the greatest mean and variance in differentiation. Flying birds, flying mammals and swimming marine mammals showed the least differentiation. Mean physiological vagility of different groups (class and locomotor mode) accounted for 98% of the mean variation in genetic differentiation with distance in each group. Genetic differentiation with distance was not related to body mass. The physiological capacity for movement (physiological vagility) quantitatively predicts genetic isolation by distance in the vertebrates examined.

Compensatory growth strategies are affected by the strength of environmental time constraints in anuran larvae

Organisms normally grow at a sub-maximal rate. After experiencing a period of arrested growth, individuals often show compensatory growth responses by modifying their life-history, behaviour and physiology. However, the strength of compensatory responses may vary across broad geographic scales as populations differ in their exposition to varying time constraints. We examined differences in compensatory growth strategies in common frog (Rana temporaria) populations from southern and northern Sweden. Tadpoles from four populations were reared in the laboratory and exposed to low temperature to evaluate the patterns and mechanisms of compensatory growth responses. We determined tadpoles' growth rate, food intake and growth efficiency during the compensation period. In the absence of arrested growth conditions, tadpoles from all the populations showed similar (size-corrected) growth rates, food intake and growth efficiency. After being exposed to low temperature for 1 week, only larvae from the northern populations increased growth rates by increasing both food intake and growth efficiency. These geographic differences in compensatory growth mechanisms suggest that the strategies for recovering after a period of growth deprivation may depend on the strength of time constraints faced by the populations. Due to the costs of fast growth, only populations exposed to the strong time constraints are prone to develop fast recovering strategies in order to metamorphose before conditions deteriorate. Understanding how organisms balance the cost and benefits of growth strategies may help in forecasting the impact of fluctuating environmental conditions on life-history strategies of populations likely to be exposed to increasing environmental variation in the future.

Rapid range expansion increases genetic differentiation while causing limited reduction in genetic diversity in a damselfly

Many ectothermic species are currently expanding their geographic range due to global warming. This can modify the population genetic diversity and structure of these species because of genetic drift during the colonization of new areas. Although the genetic signatures of historical range expansions have been investigated in an array of species, the genetic consequences of natural, contemporary range expansions have received little attention, with the only studies available focusing on range expansions along a narrow front. We investigate the genetic consequences of a natural range expansion in the Mediterranean damselfly Coenagrion scitulum, which is currently rapidly expanding along a broad front in different directions. We assessed genetic diversity and genetic structure using 12 microsatellite markers in five centrally located populations and five recently established populations at the edge of the geographic distribution. Our results suggest that, although a marginal significant decrease in the allelic richness was found in the edge populations, genetic diversity has been preserved during the range expansion of this species. Nevertheless, edge populations were genetically more differentiated compared with core populations, suggesting genetic drift during the range expansion. The smaller effective population sizes of the edge populations compared with central populations also suggest a contribution of genetic drift after colonization. We argue and document that range expansion along multiple axes of a broad expansion front generates little reduction in genetic diversity, yet stronger differentiation of the edge populations.

Using genetic variation to infer associations with climate in the common frog, Rana temporaria

Recent and historical species' associations with climate can be inferred using molecular markers. This knowledge of population and species-level responses to climatic variables can then be used to predict the potential consequences of ongoing climate change. The aim of this study was to predict responses of Rana temporaria to environmental change in Scotland by inferring historical and contemporary patterns of gene flow in relation to current variation in local thermal conditions. We first inferred colonization patterns within Europe following the last glacial maximum by combining new and previously published mitochondrial DNA sequences. We found that sequences from our Scottish samples were identical to (92%), or clustered with, the common haplotype previously identified from Western Europe. This clade showed very low mitochondrial variation, which did not allow inference of historical colonization routes but did allow interpretation of patterns of current fine-scale population structure without consideration of confounding historical variation. Second, we assessed fine-scale microsatellite-based patterns of genetic variation in relation to current altitudinal temperature gradients. No population structure was found within altitudinal gradients (average FST=0.02), despite a mean annual temperature difference of 4.5 °C between low- and high-altitude sites. Levels of genetic diversity were considerable and did not vary between sites. The panmictic population structure observed, even along temperature gradients, is a potentially positive sign for R. temporaria persistence in Scotland in the face of a changing climate. This study demonstrates that within taxonomic groups, thought to be at high risk from environmental change, levels of vulnerability can vary, even within species.

Potential effects of climate change on the distribution of the common frog Rana temporaria at its northern range margin

Climate change is projected to be particularly strong in northern latitudes, and subarctic species are thus likely to be especially susceptible to the effects of climate warming. We forecast potential effects of climate change on the extent of the suitable habitat of the common frog,Rana temporaria,at the margin of its northern range. We investigated 179 potential breeding sites in subarctic Finland and subjected the data to detailed bioclimate envelope modelling using three state-of-the-art techniques: generalized additive models, maximum entropy and generalized boosting methods. Moreover, we included local environmental factors in the models to investigate whether they improve model performance. Under all tested climate change projections and irrespective of the modelling method, the suitable habitat forR. temporariaincreased in warming climate. The inclusion of local abiotic variables significantly improved the performance of the models. However, June temperature appeared to be the most informative variable in all modelling approaches: a major increase in the extent of suitable habitat occurred when it increased by 1°C. Overall, the modelling results indicate that the distribution of northernR. temporariais likely to be very sensitive to climate warming. The results also highlight the fact that overlooking local abiotic variation can significantly bias bioclimatic modelling results.

Highest genetic diversity at the northern range limit of the rare orchid Isotria medeoloides

Populations in previously glaciated regions are often genetically depauperate in comparison with populations at lower latitudes, due either to bottlenecks experienced in post-glacial colonization or to contemporary genetic drift in small, peripheral populations. Populations of the rare self-fertilizing North American orchid Isotria medeoloides are largest in the previously glaciated region near the northern range limit, allowing us to examine the role of historical versus contemporary processes in determining population genetic diversity and structure. If contemporary processes predominate, genetic diversity should increase with increasing census size. In contrast, if sequential bottlenecks associated with colonization are paramount, diversity should decrease with latitude and be relatively insensitive to census size. We genotyped 299 individuals from 20 populations at four variable microsatellite loci to contrast genetic diversity and structure for populations in previously glaciated regions versus previously unglaciated regions. Populations were highly inbred (F=0.95) and highly differentiated (R(ST)=0.485). Across all sampled populations, genetic diversity decreased and genetic differentiation increased with declining population size. Small southern populations were especially differentiated and genetically depauperate. In the glaciated part of the range, genetic diversity increased as populations approached the northern range limit, demonstrating the centrality of contemporary processes for this post-glacial colonist.

Current and historical drivers of landscape genetic structure differ in core and peripheral salamander populations

With predicted decreases in genetic diversity and greater genetic differentiation at range peripheries relative to their cores, it can be difficult to distinguish between the roles of current disturbance versus historic processes in shaping contemporary genetic patterns. To address this problem, we test for differences in historic demography and landscape genetic structure of coastal giant salamanders (Dicamptodon tenebrosus) in two core regions (Washington State, United States) versus the species' northern peripheral region (British Columbia, Canada) where the species is listed as threatened. Coalescent-based demographic simulations were consistent with a pattern of post-glacial range expansion, with both ancestral and current estimates of effective population size being much larger within the core region relative to the periphery. However, contrary to predictions of recent human-induced population decline in the less genetically diverse peripheral region, there was no genetic signature of population size change. Effects of current demographic processes on genetic structure were evident using a resistance-based landscape genetics approach. Among core populations, genetic structure was best explained by length of the growing season and isolation by resistance (i.e. a ‘flat’ landscape), but at the periphery, topography (slope and elevation) had the greatest influence on genetic structure. Although reduced genetic variation at the range periphery of D. tenebrosus appears to be largely the result of biogeographical history rather than recent impacts, our analyses suggest that inherent landscape features act to alter dispersal pathways uniquely in different parts of the species' geographic range, with implications for habitat management.

Approximate Bayesian computation reveals the factors that influence genetic diversity and population structure of foxsnakes

Contemporary geographical range and patterns of genetic diversity within species reflect complex interactions between multiple factors acting across spatial and temporal scales, and it is notoriously difficult to disentangle causation. Here, we quantify patterns of genetic diversity and genetic population structure using mitochondrial DNA sequences (101 individuals, cytochrome b) and microsatellites (816 individuals, 12 loci) and use Approximate Bayesian computation methods to test competing models of the demographic history of eastern and western foxsnakes. Our analyses indicate that post-glacial colonization and past population declines, probably caused by the infilling of deciduous forest and cooler temperatures since the mid-Holocene, largely underpin large-scale genetic patterns for foxsnakes. At finer geographical scales, our results point to more recent anthropogenic habitat loss as having accentuated genetic population structure by causing further declines and fragmentation.

Genetic structure of muskrat (Ondatra zibethicus) and its concordance with taxonomy in North America

Extrinsic factors such as physical barriers play an important role in shaping population genetic structure. A reduction in gene flow leading to population structuring may ultimately lead to population divergence. These divergent populations are often considered subspecies. Because genetic differentiation may represent differences between subspecies, patterns of genetic structure should reflect subspecies groupings. In this study, we examine the contemporary population genetic structure of muskrat (n = 331) and assess the relevance of 4 geographically distinct subspecies designations across northern North America using 9 microsatellite loci. We predicted that patterns of gene flow and genetic structure would reflect the described subspecies. We found evidence of genetic differentiation between western and eastern regions, and muskrats from Newfoundland (NF) showed significantly lower genetic diversity than central regions. A strong isolation by distance pattern was also detected within the eastern cluster. Our results did not differentiate Ondatra zibethicus spatulus (northwest) from O. z. albus (central), but they suggest a distinction between O. z. obscurus (NF) and O. z. zibethicus (east). This study highlights the need for more phylogenetic studies in order to better understand intraspecific divergence and the genetic characterization of subspecies.

Population differentiation of temperate amphibians in unpredictable environments

Amphibians are a globally distributed and diverse lineage, but much of our current understanding of their population genetic structure comes from studies in mesic temperate habitats. We characterize the population genetic structure of two sympatric explosive breeding amphibians in the southwestern deserts of the United States: the Great Plains toad (Anaxyrus cognatus) and Couch's spadefoot toad (Scaphiopus couchii). For both species, we find limited genetic differentiation even between populations in adjacent valleys separated by dispersal barriers such as mountainous habitats. To understand how population genetic patterns in these two arid-adapted species compare to taxa in more mesic environments, we computed a standardized measure of population differentiation for A. cognatus, S. couchii, and for pond-breeding amphibians that inhabit mesic temperate environments. Our results indicate that the arid-adapted species have lower population genetic structure at fine and moderate scales than most other amphibian species we surveyed. We hypothesize that stochasticity in the availability of appropriate breeding sites as well as landscape homogeneity may result in increased population connectivity in desert-adapted frogs. Future work examining fine-scale population structure in amphibians from a diversity of habitats will test the generality of our findings. Intraspecific comparisons among localities with varied seasonality and habitats will be particularly useful for investigating the interaction between species-typical population dynamics and environmental characteristics as determinants of population connectivity in pond-breeding amphibians.

Big mountains but small barriers: population genetic structure of the Chinese wood frog (Rana chensinensis) in the Tsinling and Daba Mountain region of northern China

Background

Amphibians in general are poor dispersers and highly philopatric, and landscape features often have important impacts on their population genetic structure and dispersal patterns. Numerous studies have suggested that genetic differentiation among amphibian populations are particularly pronounced for populations separated by mountain ridges. The Tsinling Mountain range of northern China is a major mountain chain that forms the boundary between the Oriental and Palearctic zoogeographic realms. We studied the population structure of the Chinese wood frog (Rana chensinensis) to test whether the Tsinling Mountains and the nearby Daba Mountains impose major barriers to gene flow.

Results

Using 13 polymorphic microsatellite DNA loci, 523 individuals from 12 breeding sites with geographical distances ranging from 2.6 to 422.8 kilometers were examined. Substantial genetic diversity was detected at all sites with an average of 25.5 alleles per locus and an expected heterozygosity ranging from 0.504 to 0.855, and two peripheral populations revealed significantly lower genetic diversity than the central populations. In addition, the genetic differentiation among the central populations was statistically significant, with pairwise F_ST values ranging from 0.0175 to 0.1625 with an average of 0.0878. Furthermore, hierarchical AMOVA analysis attributed most genetic variation to the within-population component, and the between-population variation can largely be explained by isolation-by-distance. None of the putative barriers detected from genetic data coincided with the location of the Tsinling Mountains.

Conclusion

The Tsinling and Daba Mountains revealed no significant impact on the population genetic structure of R. chensinensis. High population connectivity and extensive juvenile dispersal may account for the significant, but moderate differentiation between populations. Chinese wood frogs are able to use streams as breeding sites at high elevations, which may significantly contribute to the diminishing barrier effect of mountain ridges. Additionally, a significant decrease in genetic diversity in the peripheral populations supports Mayr's central-peripheral population hypothesis.