Demographic and genetic estimates of effective population and breeding size in the amphibian Rana temporaria

Agricultural intensification alters marbled newt genetic diversity and gene flow through density and dispersal reduction

Recent agricultural intensification threatens global biodiversity with amphibians being one of the most impacted groups. Because of their biphasic life cycle, amphibians are particularly vulnerable to habitat loss and fragmentation that often result in small, isolated populations and loss of genetic diversity. Here, we studied how landscape heterogeneity affects genetic diversity, gene flow and demographic parameters in the marbled newt, Triturus marmoratus, over a hedgerow network landscape in Western France. While the northern part of the study area consists of preserved hedged farmland, the southern part was more profoundly converted for intensive arable crops production after WWII. Based on 67 sampled ponds and 10 microsatellite loci, we characterized regional population genetic structure and evaluated the correlation between landscape variables and (i) local genetic diversity using mixed models and (ii) genetic distance using multiple regression methods and commonality analysis. We identified a single genetic population characterized by a spatially heterogeneous isolation-by-distance pattern. Pond density in the surrounding landscape positively affected local genetic diversity while arable crop land cover negatively affected gene flow and connectivity. We used demographic inferences to quantitatively assess differences in effective population density and dispersal between the contrasted landscapes characterizing the northern and southern parts of the study area. Altogether, results suggest recent land conversion affected T. marmoratus through reduction in both effective population density and dispersal due to habitat loss and reduced connectivity.

Integrative demographic study of the Iberian painted frog (Discoglossus galganoi): inter-annual variation in the effective to census population size ratio, with insights on mating system and breeding success

In the face of worldwide amphibian declines, integrative studies combining individual-based information and genetic data represent a powerful approach to produce robust, reliable, and comparable assessments of demographic dynamics. The Iberian painted frog (Discoglossus galganoi) is endemic to Spain and Portugal and shows decreasing population trends across its range, but few studies have attempted to estimate census sizes or assess genetic diversity in wild populations, and little is known about their reproductive biology. We applied an integrative approach based on the combination of capture-mark-recapture data and multilocus genotypes to monitor a breeding population of D. galganoi in central Spain during two consecutive breeding seasons, focusing on the estimation of demographic parameters and their temporal variation. Specifically, we estimated the number of adults (N_a ), the effective population size (N_e ), and the effective number of breeders (N_b ), as well as survival and migration rates. We documented a >50% decrease in the estimated number of adults of both sexes between the breeding seasons of 2018 and 2019, probably associated with reduced rainfall in the latter. Estimates of N_b and the N_b /N_a ratio were low in both seasons, with a 20-30% decrease in N_b and a 47% increase in the N_b /N_a ratio in 2019. Based on the reconstruction of pedigrees from larval and adult genotypes, we provide the first genetic evidence of polygamy in males and females of D. galganoi and the first estimates of breeding success in the species.

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.

Reducing inbreeding rates with a breeding circle: Theory and practice in Veluws Heideschaap

Breeding circles allow genetic management in closed populations without pedigrees. In a breeding circle, breeding is split over sub‐populations. Each sub‐population receives breeding males from a single sub‐population and supplies breeding males to one other sub‐population. Donor‐recipient combinations of sub‐populations remain the same over time. Here, we derive inbreeding levels both mathematically and by computer simulation and compare them to actual inbreeding rates derived from DNA information in a real sheep population. In Veluws Heideschaap, a breeding circle has been in operation for over 30 years. Mathematically, starting with inbreeding levels and kinships set to zero, inbreeding rates per generation (ΔF) initially were 0.29%–0.47% within flocks but later converged to 0.18% in all flocks. When, more realistically, inbreeding levels at the start were high and kinship between flocks low, inbreeding levels immediately dropped to the kinship levels between flocks and rates more gradually converged to 0.18%. In computer simulations with overlapping generations, inbreeding levels and rates followed the same pattern, but converged to a lower ΔF of 0.12%. ΔF was determined in the real population with a 12 K SNP chip in recent generations. ΔF in the real population was 0.29%, based on markers to 0.41% per generation based on heterozygosity levels. This is two to three times the theoretically derived values. These increased rates in the real population are probably due to selection and/or the presence of dominant rams siring a disproportionate number of offspring. When these were simulated, ΔF agreed better: 0.35% for selection, 0.38% for dominant rams and 0.67% for both together. The realized inbreeding rates are a warning that in a real population inbreeding rates in a breeding circle can be higher than theoretically expected due to selection and dominant rams. Without a breeding circle, however, inbreeding rates would have been even higher.

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.

Genomic data recover previously undetectable fragmentation effects in an endangered amphibian

A critical consideration when using molecular ecological methods to detect trends and parameterize models at very fine spatial and temporal scales has always been the technical limits of resolution. Key landscape features, including most anthropogenic modifications, can cause biologically important, but very recent changes in gene flow that require substantial statistical power to detect. The problem is one of temporal scale: Human change is rapid and recent, while genetic changes accumulate slowly. We generated SNPs from thousands of nuclear loci to characterize the population structure of New York-endangered eastern tiger salamanders (Ambystoma tigrinum) on Long Island and quantify the impacts of roads on population fragmentation. In stark contrast to a recent microsatellite study, we uncovered highly structured populations over an extremely small spatial scale (approximately 40 km² ) in an increasingly human-modified landscape. Geographic distance and the presence of roads between ponds were both strong predictors of genetic divergence, suggesting that both natural and anthropogenic factors contribute to the observed patterns of genetic variation. All ponds supported small to modest effective breeding populations, and pond surface area showed a strong positive correlation with population size. None of these patterns emerged in an earlier study of the same system using microsatellite loci, and we determined that at least 300-400 SNPs were needed to recover the fine-scale population structure present in this system. Conservation assessments using earlier genetic techniques in other species may similarly lack the statistical power for small-scale inferences and benefit from reassessments using genomic tools.

Jailed in the mountains: Genetic diversity and structure of an endemic newt species across the Pyrenees

The Pyrenees represent a natural laboratory for biogeographic, evolutionary and ecological research of mountain fauna as a result of the high variety of habitats and the profound effect of the glacial and interglacial periods. There is a paucity of studies providing a detailed insight into genetic processes and better knowledge on the patterns of genetic diversity and how they are maintained under high altitude conditions. This is of particular interest when considering the course of past climate conditions and glaciations in a species which is considered site tenacious, with long generation times. Here we analyzed the genetic patterns of diversity and structure of the endemic Pyrenean brook newt (Calotriton asper) along its distribution range, with special emphasis on the distinct habitat types (caves, streams, and lakes), and the altitudinal and geographical ranges, using a total set of 900 individuals from 44 different localities across the Pyrenean mountain range genotyped for 19 microsatellite loci. We found evidence for a negative longitudinal and positive altitudinal gradient of genetic diversity in C. asper populations. The fact that genetic diversity was markedly higher westwards is in accordance with other Pyrenean species. However, the impact of altitudinal gradient on the genetic diversity seems to differ from other species, and mostly from other amphibians. We found that lower altitudes can act as a barrier probably because the lowlands do not provide a suitable habitat for C. asper. Regarding the distinct habitat types, caves had significantly lower values of genetic diversity compared to streams or lakes. The mean F_ST value was relatively high (0.304) with maximum values as high as 0.771, suggesting a highly structured total population. Indeed, populations were grouped into five subclusters, the eastern populations (cluster 1) remained grouped into two subclusters and the central-western Pyrenees (cluster 2) into three subclusters. The increase of isolation with geographical distance is consistent with the population structure detected. In conclusion, C. asper seems to be adapted to high altitude mountain habitats, and its genetic diversity is higher in the western Pyrenees. In terms of conservation priority, we consider more relevant the populations that represent a reservoir of genetic diversity.

Population genetic analysis of the recently rediscovered Hula painted frog (Latonia nigriventer) reveals high genetic diversity and low inbreeding

After its recent rediscovery, the Hula painted frog (Latonia nigriventer) has remained one of the world's rarest and least understood amphibian species. Together with its apparently low dispersal capability and highly disturbed niche, the low abundance of this living fossil calls for urgent conservation measures. We used 18 newly developed microsatellite loci and four different models to calculate the effective population size (N_e) of a total of 125 Hula painted frog individuals sampled at a single location. We compare the N_e estimates to the estimates of potentially reproducing adults in this population (N_ad) determined through a capture-recapture study on 118 adult Hula painted frogs captured at the same site. Surprisingly, our data suggests that, despite N_ad estimates of only ~234-244 and N_e estimates of ~16.6-35.8, the species appears to maintain a very high genetic diversity (H_O = 0.771) and low inbreeding coefficient (F_IS = -0.018). This puzzling outcome could perhaps be explained by the hypotheses of either genetic rescue from one or more unknown Hula painted frog populations nearby or by recent admixture of genetically divergent subpopulations. Independent of which scenario is correct, the original locations of these populations still remain to be determined.

Reliable effective number of breeders/adult census size ratios in seasonal-breeding species: Opportunity for integrative demographic inferences based on capture-mark-recapture data and multilocus genotypes

The ratio of the effective number of breeders (Nb) to the adult census size (Na), Nb/Na, approximates the departure from the standard capacity of a population to maintain genetic diversity in one reproductive season. This information is relevant for assessing population status, understanding evolutionary processes operating at local scales, and unraveling how life-history traits affect these processes. However, our knowledge on Nb/Na ratios in nature is limited because estimation of both parameters is challenging. The sibship frequency (SF) method is adequate for reliable Nb estimation because it is based on sibship and parentage reconstruction from genetic marker data, thereby providing demographic inferences that can be compared with field-based information. In addition, capture-mark-recapture (CMR) robust design methods are well suited for Na estimation in seasonal-breeding species. We used tadpole genotypes of three pond-breeding amphibian species (Epidalea calamita, Hyla molleri, and Pelophylax perezi, n = 73-96 single-cohort tadpoles/species genotyped at 15-17 microsatellite loci) and candidate parental genotypes (n = 94-300 adults/species) to estimate Nb by the SF method. To assess the reliability of Nb estimates, we compared sibship and parentage inferences with field-based information and checked for the convergence of results in replicated subsampled analyses. Finally, we used CMR data from a 6-year monitoring program to estimate annual Na in the three species and calculate the Nb/Na ratio. Reliable ratios were obtained for E. calamita (Nb/Na = 0.18-0.28) and P. perezi (0.5), but in H. molleri, Na could not be estimated and genetic information proved insufficient for reliable Nb estimation. Integrative demographic studies taking full advantage of SF and CMR methods can provide accurate estimates of the Nb/Na ratio in seasonal-breeding species. Importantly, the SF method provides results that can be readily evaluated for reliability. This represents a good opportunity for obtaining robust demographic inferences with wide applications for evolutionary and conservation research.

Population genetic and field-ecological analyses return similar estimates of dispersal over space and time in an endangered amphibian

The explosive growth of empirical population genetics has seen a proliferation of analytical methods leading to a steady increase in our ability to accurately measure key population parameters, including genetic isolation, effective population size, and gene flow, in natural systems. Assuming they yield similar results, population genetic methods offer an attractive complement to, or replacement of, traditional field‐ecological studies. However, empirical assessments of the concordance between direct field‐ecological and indirect population genetic studies of the same populations are uncommon in the literature. In this study, we investigate genetic isolation, rates of dispersal, and population sizes for the endangered California tiger salamander, Ambystoma californiense, across multiple breeding seasons in an intact vernal pool network. We then compare our molecular results to a previously published study based on multiyear, mark–recapture data from the same breeding sites. We found that field and genetic estimates of population size were only weakly correlated, but dispersal rates were remarkably congruent across studies and methods. In fact, dispersal probability functions derived from genetic data and traditional field‐ecological data were a significant match, suggesting that either method can be used effectively to assess population connectivity. These results provide one of the first explicit tests of the correspondence between landscape genetic and field‐ecological approaches to measuring functional population connectivity and suggest that even single‐year genetic samples can return biologically meaningful estimates of natural dispersal and gene flow.

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.

Low reproductive skew despite high male-biased operational sex ratio in a glass frog with paternal care

BACKGROUND

Reproductive skew, the uneven distribution of reproductive success among individuals, is a common feature of many animal populations. Several scenarios have been proposed to favour either high or low levels of reproductive skew. Particularly a male-biased operational sex ratio and the asynchronous arrival of females is expected to cause high variation in reproductive success among males. Recently it has been suggested that the type of benefits provided by males (fixed vs. dilutable) could also strongly impact individual mating patterns, and thereby affecting reproductive skew. We tested this hypothesis in Hyalinobatrachium valerioi, a Neotropical glass frog with prolonged breeding and paternal care.

RESULTS

We monitored and genetically sampled a natural population in southwestern Costa Rica during the breeding season in 2012 and performed parentage analysis of adult frogs and tadpoles to investigate individual mating frequencies, possible mating preferences, and estimate reproductive skew in males and females. We identified a polygamous mating system, where high proportions of males (69 %) and females (94 %) reproduced successfully. The variance in male mating success could largely be attributed to differences in time spent calling at the reproductive site, but not to body size or relatedness. Female H. valerioi were not choosy and mated indiscriminately with available males.

CONCLUSIONS

Our findings support the hypothesis that dilutable male benefits - such as parental care - can favour female polyandry and maintain low levels of reproductive skew among males within a population, even in the presence of direct male-male competition and a highly male-biased operational sex ratio. We hypothesize that low male reproductive skew might be a general characteristic in prolonged breeders with paternal care.

Simple life-history traits explain key effective population size ratios across diverse taxa

Effective population size (Ne) controls both the rate of random genetic drift and the effectiveness of selection and migration, but it is difficult to estimate in nature. In particular, for species with overlapping generations, it is easier to estimate the effective number of breeders in one reproductive cycle (Nb) than Ne per generation. We empirically evaluated the relationship between life history and ratios of Ne, Nb and adult census size (N) using a recently developed model (agene) and published vital rates for 63 iteroparous animals and plants. Nb/Ne varied a surprising sixfold across species and, contrary to expectations, Nb was larger than Ne in over half the species. Up to two-thirds of the variance in Nb/Ne and up to half the variance in Ne/N was explained by just two life-history traits (age at maturity and adult lifespan) that have long interested both ecologists and evolutionary biologists. These results provide novel insights into, and demonstrate a close general linkage between, demographic and evolutionary processes across diverse taxa. For the first time, our results also make it possible to interpret rapidly accumulating estimates of Nb in the context of the rich body of evolutionary theory based on Ne per generation.

What remains from a 454 run: estimation of success rates of microsatellite loci development in selected newt species (Calotriton asper, Lissotriton helveticus, and Triturus cristatus) and comparison with Illumina-based approaches

The development of microsatellite loci has become more efficient using next-generation sequencing (NGS) approaches, and many studies imply that the amount of applicable loci is large. However, few studies have sought to quantify the number of loci that are retained for use out of the thousands of sequence reads initially obtained. We analyzed the success rate of microsatellite loci development for three amphibian species using a 454 NGS approach on tetra-nucleotide motif-enriched species-specific libraries. The number of sequence reads obtained differed strongly between species and ranged from 19,562 for Triturus cristatus to 55,626 for Lissotriton helveticus, with 52,075 reads obtained for Calotriton asper. PHOBOS was used to identify sequences with tetra-nucleotide repeat motifs with a minimum repeat number of ten and high quality primer binding sites. Of 107 sequences for T. cristatus, 316 for C. asper and 319 for L. helveticus, we tested the amplification success, polymorphism, and degree of heterozygosity for 41 primer combinations each for C. asper and T. cristatus, and 22 for L. helveticus. We found 11 polymorphic loci for T. cristatus, 20 loci for C. asper, and 15 loci for L. helveticus. Extrapolated, the number of potentially amplifiable loci (PALs) resulted in estimated species-specific success rates of 0.15% (T. cristatus), 0.30% (C. asper), and 0.39% (L. helveticus). Compared with representative Illumina NGS approaches, our applied 454-sequencing approach on specifically enriched sublibraries proved to be quite competitive in terms of success rates and number of finally applicable loci.

Radically different phylogeographies and patterns of genetic variation in two European brown frogs, genus Rana

We reconstruct range-wide phylogeographies of two widespread and largely co-occurring Western Palearctic frogs, Rana temporaria and R. dalmatina. Based on tissue or saliva samples of over 1000 individuals, we compare a variety of genetic marker systems, including mitochondrial DNA, single-copy protein-coding nuclear genes, microsatellite loci, and single nucleotide polymorphisms (SNPs) of transcriptomes of both species. The two focal species differ radically in their phylogeographic structure, with R. temporaria being strongly variable among and within populations, and R. dalmatina homogeneous across Europe with a single strongly differentiated population in southern Italy. These differences were observed across the various markers studied, including microsatellites and SNP density, but especially in protein-coding nuclear genes where R. dalmatina had extremely low heterozygosity values across its range, including potential refugial areas. On the contrary, R. temporaria had comparably high range-wide values, including many areas of probable postglacial colonization. A phylogeny of R. temporaria based on various concatenated mtDNA genes revealed that two haplotype clades endemic to Iberia form a paraphyletic group at the base of the cladogram, and all other haplotypes form a monophyletic group, in agreement with an Iberian origin of the species. Demographic analysis suggests that R. temporaria and R. dalmatina have genealogies of roughly the same time to coalescence (TMRCA ~3.5 mya for both species), but R. temporaria might have been characterized by larger ancestral and current effective population sizes than R. dalmatina. The high genetic variation in R. temporaria can therefore be explained by its early range expansion out of Iberia, with subsequent cycles of differentiation in cryptic glacial refugial areas followed by admixture, while the range expansion of R. dalmatina into central Europe is a probably more recent event.

Evolutionary dynamics of a rapidly receding southern range boundary in the threatened California Red-Legged Frog (Rana draytonii)

Populations forming the edge of a species range are often imperiled by isolation and low genetic diversity, with proximity to human population centers being a major determinant of edge stability in modern landscapes. Since the 1960s, the California red‐legged frog (Rana draytonii) has undergone extensive declines in heavily urbanized southern California, where the range edge has rapidly contracted northward while shifting its cardinal orientation to an east‐west trending axis. We studied the genetic structure and diversity of these frontline populations, tested for signatures of contemporary disturbance, specifically fire, and attempted to disentangle these signals from demographic events extending deeper into the past. Consistent with the genetic expectations of the ‘abundant‐center’ model, we found that diversity, admixture, and opportunity for random mating increases in populations sampled successively further away from the range boundary. Demographic simulations indicate that bottlenecks in peripheral isolates are associated with processes extending tens to a few hundred generations in the past, despite the demographic collapse of some due to recent fire‐flood events. While the effects of recent disturbance have left little genetic imprint on these populations, they likely contribute to an extinction debt that will lead to continued range contraction unless management intervenes to stall or reverse the process.

Can genetic estimators provide robust estimates of the effective number of breeders in small populations?

The effective population size (N(e)) is proportional to the loss of genetic diversity and the rate of inbreeding, and its accurate estimation is crucial for the monitoring of small populations. Here, we integrate temporal studies of the gecko Oedura reticulata, to compare genetic and demographic estimators of N(e). Because geckos have overlapping generations, our goal was to demographically estimate N(bI), the inbreeding effective number of breeders and to calculate the N(bI)/N(a) ratio (N(a) =number of adults) for four populations. Demographically estimated N(bI) ranged from 1 to 65 individuals. The mean reduction in the effective number of breeders relative to census size (N(bI)/N(a)) was 0.1 to 1.1. We identified the variance in reproductive success as the most important variable contributing to reduction of this ratio. We used four methods to estimate the genetic based inbreeding effective number of breeders N(bI(gen)) and the variance effective populations size N(eV(gen)) estimates from the genotype data. Two of these methods - a temporal moment-based (MBT) and a likelihood-based approach (TM3) require at least two samples in time, while the other two were single-sample estimators - the linkage disequilibrium method with bias correction LDNe and the program ONeSAMP. The genetic based estimates were fairly similar across methods and also similar to the demographic estimates excluding those estimates, in which upper confidence interval boundaries were uninformative. For example, LDNe and ONeSAMP estimates ranged from 14-55 and 24-48 individuals, respectively. However, temporal methods suffered from a large variation in confidence intervals and concerns about the prior information. We conclude that the single-sample estimators are an acceptable short-cut to estimate N(bI) for species such as geckos and will be of great importance for the monitoring of species in fragmented landscapes.

Toward responsible stock enhancement: broadcast spawning dynamics and adaptive genetic management in white seabass aquaculture

The evolutionary effects captive-bred individuals that can have on wild conspecifics are necessary considerations for stock enhancement programs, but breeding protocols are often developed without the knowledge of realized reproductive behavior. To help fill that gap, parentage was assigned to offspring produced by a freely mating group of 50 white seabass (Atractoscion nobilis), a representative broadcast spawning marine finfish cultured for conservation. Similar to the well-known and closely related red drum (Sciaenops ocellatus), A. nobilis exhibited large variation in reproductive success. More males contributed and contributed more equally than females within and among spawns in a mating system best described as lottery polygyny. Two females produced 27% of the seasonal offspring pool and female breeding effective size averaged 1.85 per spawn and 12.38 seasonally, whereas male breeding effective size was higher (6.42 and 20.87, respectively), with every male contributing 1-7% of offspring. Further, females batch spawned every 1-5 weeks, while males displayed continuous reproductive readiness. Sex-specific mating strategies resulted in multiple successful mate pairings and a breeding effective to census size ratio of ≥0.62. Understanding a depleted species' mating system allowed management to more effectively utilize parental genetic variability for culture, but the fitness consequences of long-term stocking can be difficult to address.

Discordant patterns of nuclear and mitochondrial introgression in Iberian populations of the European common frog (Rana temporaria)

Amphibians often show complex histories of intraspecific and interspecific genetic introgression, which might differ in mitochondrial and nuclear genes. In our study of the genetic differentiation of the European common frog, Rana temporaria (159 specimens from 23 populations were analyzed for 24 presumptive allozyme loci; 82 specimens were sequenced for a 540-bp fragment of the mitochondrial 16S rRNA gene), multilocus correspondence analysis (CA) and Bayesian assignment tests of the nuclear data were concordant in identifying 2 population groups corresponding to 1) the Pyrenees in the east and 2) the Galicia and Asturias regions in the west, the latter corresponding to the subspecies R. temporaria parvipalmata. Geographically intermediate populations were genetically intermediate in the allozyme CA and, less clearly in the Bayesian assignment, with mitochondrial haplotypes exclusively belonging to the parvipalmata group. This indicates different degrees of introgression in the mitochondrial and nuclear genomes. Although Pyrenean high-altitude populations are morphologically distinct from low-altitude populations, these 2 groups were not separate clusters in any analysis. This suggests that the morphological differences may be due to fast adaptation to elevational gradients, likely under maintenance of gene flow, and that the underlying genetic changes are not detectable by the analyzed markers. We argue that a parsimonious explanation for the observed pattern along the east-west axis in northern Spain may be competition between invading and resident populations, with no need to invoke selection. However, in order to conclusively rule out selective processes, additional and finer scale data are required to test for asymmetric mating preference/behaviour, sex-biased gene flow, or sex-biased survival of potential hybrids.

Study of closely related species within the Physalaemus cuvieri group (Anura): contribution of microsatellite markers

Various species of the Physalaemus cuvieri group of frogs are difficult to distinguish morphologically, making molecular analysis an attractive alternative for indentifying members of this group, which is considered to be at risk because of loss of habitat. The genetic structure of natural populations of P. ephippifer and P. albonotatus species was investigated and analyzed, together with that of five previously studied populations of P. cuvieri. Nine microsatellite loci were used in the analyses. The overall G(ST) value (0.46) revealed high genetic variation among the populations, as expected for different species. Bayesian analysis implemented by the STRUCTURE software clustered the seven populations into seven groups (K = 7). All the P. albonotatus and P. ephippifer specimens were grouped into a single cluster, both species showing clear differentiation from P. cuvieri. The different grouping based on these microsatellites of some P. cuvieri individuals from Porto Nacional and from Passo Fundo suggests that they could be a new species, indicating a necessity for taxonomic reevaluation. Despite the intrinsic difficulties in analyzing closely related species, the nine microsatellite loci were found to be adequate for distinguishing these three species of the P. cuvieri group and their populations.

Phenotypic plasticity allows the Mediterranean parsley frog Pelodytes punctatus to exploit two temporal niches under continuous gene flow

Environmental changes, such as climate change, lead to the opening of new niches. In such situations, species that adapt to new niches can survive and/or expand their ranges. However, gene flow can hamper genetic adaptation to new environments. Alternatively, recent models have highlighted the importance of phenotypic plasticity in tracking environmental change. Here, we investigate whether phenotypic plasticity or genetic evolution (or both) allows an amphibian species to exploit two divergent climatic niches. In the Mediterranean region, the parsley frog Pelodytes punctatus breeds both in spring, as do most other species, and in autumn, a temporal niche not exploited by most other species, but which may become increasingly important with global warming. Conditions of development are dramatically different between the two seasons and deeply impact tadpole life-history traits. To determine whether these temporal niches are exploited by two genetically differentiated subpopulations, or whether the bimodal phenology arises in a panmictic population displaying plastic life-history traits, we use two complementary approaches. We measure both molecular genetic differentiation and quantitative-trait differentiation between spring and autumn cohorts, using microsatellites and common garden experiments, respectively. Seasonal cohorts were not genetically differentiated and differences in tadpole life history between cohorts were not maintained in laboratory conditions. We conclude that phenotypic plasticity, rather than genetic adaptation, allows Parsley frog to exploit two contrasting temporal niches.

Comparative analyses of effective population size within and among species: ranid frogs as a case study

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

Genetic diversity and molecular differentiation of Chinese toad based on microsatellite markers

Genetic diversity and population structure of 9 populations of Bufo gargarizans with total 111 samples in China were assessed using seven microsatellite loci. The analysed microsatellite markers produced 161 alleles, varied from 9 to 38 alleles each locus. The number of alleles per population per locus ranged from 4.43 to 10.29. Polymorphic information content showed that all seven loci were highly informative (mean = 0.810 + or - 0.071). The average observed heterozygosity was less than the expected (0.353 + or - 0.051 and 0.828 + or - 0.067, respectively). All tested populations gave significant departures from Hardy-Weinberg equilibrium. Genetic differentiation among the populations was considerably high with the overall and pairwise F (ST) values (mean = 0.160 + or - 0.039), and showed fairly high level of inbreeding (indicated by a mean F (IS) value of 0.504 + or - 0.051) and global heterozygote deficit. In comparison to other amphibian studies; however, our results suggested that the level of genetic structuring in B. gargarizans was relatively low in the geographical scale of the study area. Interestingly, the speculated population bottleneck was found to be absent and the analyses provide only weak evidence for a recent contraction in size even though there was severe inbreeding (indicated by the F (IS) value) in the Chinese toad populations.

Opportunity for sexual selection and effective population size in the lek-breeding European treefrog (Hyla arborea)

Sexual selection in lek-breeding species might drastically lower male effective population size, with potentially important consequences for evolutionary and conservation biology. Using field-monitoring and parental-assignment methods, we analyzed sex-specific variances in breeding success in a population of European treefrogs, to (1) help understanding the dynamics of genetic variance at sex-specific loci, and (2) better quantify the risk posed by genetic drift in this species locally endangered by habitat fragmentation. The variance in male mating success turned out to be markedly lower than values obtained from other amphibian species with polygamous mating systems. The ratio of effective breeding size to census breeding size was only slightly lower in males (0.44) than in females (0.57), in line with the patterns of genetic diversity previously reported from H. arborea sex chromosomes. Combining our results with data on age at maturity and adult survival, we show that the negative effect of the mating system is furthermore compensated by the effect of delayed maturity, so that the estimated instantaneous effective size broadly corresponded to census breeding size. We conclude that the lek-breeding system of treefrogs impacts only weakly the patterns of genetic diversity on sex-linked genes and the ability of natural populations to resist genetic drift.