Mating system and early viability resistance to habitat fragmentation in a bird-pollinated eucalypt

Habitat fragmentation has been shown to disrupt ecosystem processes such as plant-pollinator mutualisms. Consequently, mating patterns in remnant tree populations are expected to shift towards increased inbreeding and reduced pollen diversity, with fitness consequences for future generations. However, mating patterns and phenotypic assessments of open-pollinated progeny have rarely been combined in a single study. Here, we collected seeds from 37 Eucalyptus incrassata trees from contrasting stand densities following recent clearance in a single South Australian population (intact woodland=12.6 trees ha⁻¹; isolated pasture=1.7 trees ha⁻¹; population area=10 km²). 649 progeny from these trees were genotyped at eight microsatellite loci. We estimated genetic diversity, spatial genetic structure, indirect contemporary pollen flow and mating patterns for adults older than the clearance events and open-pollinated progeny sired post-clearance. A proxy of early stage progeny viability was assessed in a common garden experiment. Density had no impact on mating patterns, adult and progeny genetic diversity or progeny growth, but was associated with increased mean pollen dispersal. Weak spatial genetic structure among adults suggests high historical gene flow. We observed preliminary evidence for inbreeding depression related to stress caused by fungal infection, but which was not associated with density. Higher observed heterozygosities in adults compared with progeny may relate to weak selection on progeny and lifetime-accumulated mortality of inbred adults. E. incrassata appears to be resistant to the negative mating pattern and fitness changes expected within fragmented landscapes. This pattern is likely explained by strong outcrossing and regular long-distance pollen flow.     

Linking extinction-colonization dynamics to genetic structure in a salamander metapopulation

Theory predicts that founder effects have a primary role in determining metapopulation genetic structure. However, ecological factors that affect extinction-colonization dynamics may also create spatial variation in the strength of genetic drift and migration. We tested the hypothesis that ecological factors underlying extinction-colonization dynamics influenced the genetic structure of a tiger salamander (Ambystoma tigrinum) metapopulation. We used empirical data on metapopulation dynamics to make a priori predictions about the effects of population age and ecological factors on genetic diversity and divergence among 41 populations. Metapopulation dynamics of A. tigrinum depended on wetland area, connectivity and presence of predatory fish. We found that newly colonized populations were more genetically differentiated than established populations, suggesting that founder effects influenced genetic structure. However, ecological drivers of metapopulation dynamics were more important than age in predicting genetic structure. Consistent with demographic predictions from metapopulation theory, genetic diversity and divergence depended on wetland area and connectivity. Divergence was greatest in small, isolated wetlands where genetic diversity was low. Our results show that ecological factors underlying metapopulation dynamics can be key determinants of spatial genetic structure, and that habitat area and isolation may mediate the contributions of drift and migration to divergence and evolution in local populations.     

Functional characterization of gynodioecy in Fragaria vesca ssp. bracteata (Rosaceae)

Background and Aims

Gynodioecy is a phylogenetically widespread and important sexual system where females coexist with hermaphrodites. Because dioecy can arise from gynodioecy, characterization of gynodioecy in close relatives of dioecious and sub-dioecious species can provide insight into this transition. Thus, we sought to determine whether Fragaria vesca ssp. bracteata, a close relative to F. chiloensis and F. virginiana, exhibits the functional and population genetic hallmarks of a gynodioecious species.

Methods

We compared reproductive allocation of females and hermaphrodites grown in the greenhouse and estimated genetic diversity (allelic diversity, heterozygosity) and inbreeding coefficients for field-collected adults of both sexes using simple sequence repeat (SSR) markers. We estimated mating system and early seed fitness from open-pollinated families of both sex morphs.

Key Results

Under greenhouse conditions, females and hermaphrodites allocated similarly to all reproductive traits except flower number, and, as a consequence, females produced 30 % fewer seeds per plant than hermaphrodites. Under natural conditions, hermaphrodites produce seeds by self-fertilization approx. 75 % of the time, and females produced outcrossed seeds with very little biparental inbreeding. Consistent with inbreeding depression, seeds from open-pollinated hermaphrodites were less likely to germinate than those from females, and family-level estimates of hermaphrodite selfing rates were negatively correlated with germination success and speed. Furthermore, estimates of inbreeding depression based on genetic markers and population genetic theory indicate that inbreeding depression in the field could be high.

Conclusions

The joint consideration of allocation and mating system suggests that compensation may be sufficient to maintain females given the current understanding of sex determination. Fragaria vesca ssp. bracteata exhibited similar sex morph-dependent patterns of mating system and genetic diversity, but less reproductive trait dimorphism, than its sub-dioecious and dioecious congeners.     

Increased inbreeding and strong kinship structure in Taxus baccata estimated from both AFLP and SSR data

Habitat fragmentation can have severe genetic consequences for trees, such as increased inbreeding and decreased effective population size. In effect, local populations suffer from reduction of genetic variation, and thus loss of adaptive capacity, which consequently increases their risk of extinction. In Europe, Taxus baccata is among a number of tree species experiencing strong habitat fragmentation. However, there is little empirical data on the population genetic consequences of fragmentation for this species. This study aimed to characterize local genetic structure in two natural remnants of English yew in Poland based on both amplified fragment length polymorphism (AFLP) and microsatellite (SSR) markers. We introduced a Bayesian approach that estimates the average inbreeding coefficient using AFLP (dominant) markers. Results showed that, in spite of high dispersal potential (bird-mediated seed dispersal and wind-mediated pollen dispersal), English yew populations show strong kinship structure, with a spatial extent of 50–100 m, depending on the population. The estimated inbreeding levels ranged from 0.016 to 0.063, depending on the population and marker used. Several patterns were evident: (1) AFLP markers showed stronger kinship structure than SSRs; (2) AFLP markers provided higher inbreeding estimates than SSRs; and (3) kinship structure and inbreeding were more pronounced in denser populations regardless of the marker used. Our results suggest that, because both kinship structure and (bi-parental) inbreeding exist in populations of English yew, gene dispersal can be fairly limited in this species. Furthermore, at a local scale, gene dispersal intensity can be more limited in a dense population.     

High selfing and high inbreeding depression in peripheral populations of Juncus atratus

The mating system of a plant is the prime determinant of its population genetic structure. However, mating system effects may be modified by postzygotic mechanisms like inbreeding depression. Furthermore, historical as well as contemporary ecological factors and population characteristics, like the location within the species range can contribute to genetic variability. Using microsatellite markers we assessed the population genetic structure of the wind-pollinated Juncus atratus in 16 populations from peripheral and nearly central areas of the distribution range and studied the mating system of the species. In three peripheral populations, outcrossing rates at seeds stage were low (mean t(m) = 5.6%), suggesting a highly autogamous mating system. Despite this fact, on adult stage both individual heterozygosity (mean H(O) = 0.48) and gene diversity (mean H(E) = 0.58) were high even in small populations. Inbreeding coefficients were consistently low among all populations (mean F(IS) = 0.15). Within the three peripheral populations indirect estimates of lifetime inbreeding depression were surprisingly high (delta(eq) = 0.96) and inbreeding depression could be shown to act mostly on early seedling establishment. Similar conditions of autogamy combined with high inbreeding depression are typical for plants with a large lifetime genomic mutation rate that cannot avoid selfing by geitonogamy. However, the results presented here are unexpected for small-statured, herbaceous plants. Substantial genetic differentiation among all populations was found (mean F(ST) = 0.24). An isolation-by-distance pattern was apparent on large scale but not on local scale suggesting that the overall pattern was largely influenced by historical factors, e.g. colonization, whereas locally genetic drift was of greater importance than gene flow. Peripheral populations exhibited lower genetic diversity and higher inbreeding coefficients when compared with subcentral populations.     

Microevolutionary Patterns and Molecular Markers: The Genetics of Geographic Variation in Ascaris suum

Molecular markers have been used only rarely to characterize the population genetic structure of nematodes. Published studies have suggested that different taxa may show distinct genetic architectures. Isoenzyme and RAPD markers have been used to investigate geographic variation of Ascaris suum at the level of infrapopulations (nematodes within individual hosts), within localities, and among geographic regions. Independent estimates of genetic differentiation among population samples based on isoenzyme and RAPD data showed similar patterns and substantial correlation. Heterozygote deficiencies within infrapopulations and large values for inbreeding coefficients among infrapopulations suggested that the composition of these populations was not consistent with a model of random recruitment from a large panmictic pool of life-cycle stages. Both isoenzyme and RAPD markers revealed moderate levels of genetic differentiation among samples representing infrapopulations and localities. Of total gene diversity, 9.4% (isoenzyme) and 9.2% (RAPD) was partitioned among infrapopulations. Geographic localities accounted for 7.8% (isoenzyme) and 6.2% (RAPD) of total diversity. Only infrapopulations from the same farm had low levels of differentiation.     

Genetic diversity at neutral and adaptive loci determines individual fitness in a long-lived territorial bird

There is compelling evidence about the manifest effects of inbreeding depression on individual fitness and populations' risk of extinction. The majority of studies addressing inbreeding depression on wild populations are generally based on indirect measures of inbreeding using neutral markers. However, the study of functional loci, such as genes of the major histocompatibility complex (MHC), is highly recommended. MHC genes constitute an essential component of the immune system of individuals, which is directly related to individual fitness and survival. In this study, we analyse heterozygosity fitness correlations of neutral and adaptive genetic variation (22 microsatellite loci and two loci of the MHC class II, respectively) with the age of recruitment and breeding success of a decimated and geographically isolated population of a long-lived territorial vulture. Our results indicate a negative correlation between neutral genetic diversity and age of recruitment, suggesting that inbreeding may be delaying reproduction. We also found a positive correlation between functional (MHC) genetic diversity and breeding success, together with a specific positive effect of the most frequent pair of cosegregating MHC alleles in the population. Globally, our findings demonstrate that genetic depauperation in small populations has a negative impact on the individual fitness, thus increasing the populations' extinction risk.     

Temporal genetic homogeneity among shore crab (Carcinus maenas) larval events supplied to an estuarine system on the Portuguese northwest coast

Despite the importance of larval biology in the life histories of many marine animals, relatively little information exists on the dynamics and genetic composition of larval cohorts. The supply of megalopae larvae of the shore crab, Carcinus maenas, was measured on a daily basis during 8 months spread along two larval periods (2006 and 2007) at the Ria de Aveiro estuary, on the Portuguese northwest coast. A total of 10 microsatellite DNA loci were employed to explore the genetic structure, variability and relatedness of temporally distinct megalopal events, selected from the major pulses of supply. Larval variation was also compared genetically with that of a previously studied adult crabs sample, at the same loci, collected in 2006 and 2007 along the Iberian Peninsula. Results revealed a lack of genetic differentiation and identical diversity levels among larval events over time. No evidence of reduced genetic diversity between megalopae relative to the diversity assessed from the pooled sample of adults was found. Moreover, there was no evidence of any family relatedness among larvae from temporal events. The results obtained for C. maenas contradict predictions made by the sweepstakes reproduction hypothesis, in which large variance in reproductive success is expected, which is presumably detectable as sharp genetic discontinuities among separate larval events. Data here indicate conversely a high degree of temporal genetic stability among larval supply to a given estuary under variable oceanographic conditions, consistent with the hypothesis that sampled larvae were drawn from a large number of adults that do not differ in reproductive success.     

Unravelling genetics at the top: mountain islands or isolated belts?

Background and Aims

In mountain plant populations, local adaptation has been described as one of the main responses to climate warming, allowing plants to persist under stressful conditions. This is especially the case for marginal populations at their lowest elevation, as they are highly vulnerable. Adequate levels of genetic diversity are required for selection to take place, while high levels of altitudinal gene flow are seen as a major limiting factor potentially precluding local adaptation processes. Thus, a compromise between genetic diversity and gene flow seems necessary to guarantee persistence under oncoming conditions. It is therefore critical to determine if gene flow occurs preferentially between mountains at similar altitudinal belts, promoting local adaptation at the lowest populations, or conversely along altitude within each mountain.

Methods

Microsatellite markers were used to unravel genetic diversity and population structure, inbreeding and gene flow of populations at two nearby altitudinal gradients of Silene ciliata, a Mediterranean high-mountain cushion plant.

Key Results

Genetic diversity and inbreeding coefficients were similar in all populations. Substantial gene flow was found both along altitudinal gradients and horizontally within each elevation belt, although greater values were obtained along altitudinal gradients. Gene flow may be responsible for the homogeneous levels of genetic diversity found among populations. Bayesian cluster analyses also suggested that shifts along altitudinal gradients are the most plausible scenario.

Conclusions

Past population shifts associated with glaciations and interglacial periods in temperate mountains may partially explain current distributions of genetic diversity and population structure. In spite of the predominance of gene flow along the altitudinal gradients, local genetic differentiation of one of the lower populations together with the detection of one outlier locus might support the existence of different selection forces at low altitudes.     

Adult survival and microsatellite diversity in possums: effects of major histocompatibility complex-linked microsatellite diversity but not multilocus inbreeding estimators

Adult survival is perhaps the fitness parameter most important to population growth in long-lived species. Intrinsic and extrinsic covariates of survival are therefore likely to be important drivers of population dynamics. We used long-term mark-recapture data to identify genetic, individual and environmental covariates of local survival in a natural population of mountain brushtail possums (Trichosurus cunninghami). Rainfall and intra-individual diversity at microsatellite DNA markers were associated with increased local survival of adults and juveniles. We contrasted the performance of several microsatellite heterozygosity measures, including internal relatedness (IR), homozygosity by loci (HL) and the mean multilocus estimate of the squared difference in microsatellite allele sizes within an individual (mean d ²). However, the strongest effect on survival was not associated with multilocus microsatellite diversity (which would indicate a genome-wide inbreeding effect), but a subset of two loci. This included a major histocompatibility complex (MHC)-linked marker and a putatively neutral microsatellite locus. For both loci, diversity measures incorporating allele size information had stronger associations with survival than measures based on heterozygosity, whether or not allele frequency information was included (such as IR). Increased survival was apparent among heterozygotes at the MHC-linked locus, but the benefits of heterozygosity to survival were reduced in heterozygotes with larger differences in allele size. The effect of heterozygosity on fitness-related traits was supported by data on endoparasites in a subset of the individuals studied in this population. There was no apparent density dependence in survival, nor an effect of sex, age or immigrant status. Our findings suggest that in the apparent absence of inbreeding, variation at specific loci can generate strong associations between fitness and diversity at linked markers.     

Complexity and context of MHC-correlated mating preferences in wild populations

There is now substantial and growing evidence for a role of the major histocompatibility complex (MHC) in shaping individual mate preferences. In view of both its codominant expression and its function in immune response, it is often expected that females aim to avoid inbreeding or maximize offspring MHC-heterozygosity by selecting as mates those males which share fewest or no MHC alleles with themselves. However, it is becoming increasingly clear that this view is over-simplistic: not only is MHC dissimilarity just one of several (perhaps many) criteria important in mate choice decision-making, extremely MHC-dissimilar males may be avoided, and furthermore, specific alleles or combinations might be preferred if they bestow particular advantages. These points are raised in two papers in this Molecular Ecology issue, in which patterns of reproductive success in tiger salamanders ( Ambystoma tigrinum ) and three-spined sticklebacks ( Gasterosteus aculeatus ) are each inconsistent with a generalized preference for MHC dissimilarity. Together, these studies demonstrate several adaptive reasons for decision rules that do not necessarily result in maximizing mate dissimilarity.     

Strong spatial genetic structure in peripheral but not core populations of Sitka spruce [Picea sitchensis (Bong.) Carr.

We examined spatial genetic structure within eight populations of Sitka spruce classified as core or peripheral based on ecological niche, and continuous or disjunct based on species distribution. In each population, 200 trees were spatially mapped and genotyped for eight cDNA-based sequence tagged site (STS) codominant markers. Spatial autocorrelation was assessed by estimating p(ij), the average co-ancestry coefficient, between individuals within distance intervals. The distribution of alleles and genotypes within core populations was almost random, with nonsignificant co-ancestry values among trees as close as 50 m in core populations. In contrast, the distribution of alleles and genotypes within peripheral populations revealed an aggregation of similar multilocus genotypes, with co-ancestry values greater than 0.20 among trees up to 50 m apart and significant, positive values between trees up to 500 m. The relatively high density of reproductive adults in core populations may lead to highly overlapping seed shadows that limit development of spatial genetic structure. However, in peripheral populations with a lower density of adults, the distribution of alleles and genotypes was highly structured, likely due to offspring establishment near maternal trees and subsequent biparental inbreeding, as well as more recent population establishment at the leading edge of post-Pleistocene range expansion. Conserving genetic diversity in peripheral populations may require larger reserves for in situ conservation than required in core populations. These data on spatial genetic structure can be used to provide guidance for sampling strategies for both ex situ conservation and research collections.     

Quantification and reduction of bias from sampling larvae to infer population and landscape genetic structure

A well-designed sampling scheme is critical for obtaining accurate results from population genetic studies. Larval samples contain only the genetic material of successful breeders, often of a single year, and may be biased towards particular families. To quantify the bias of using larval samples to infer population and landscape genetic structure and explore how this bias may be reduced using sibship analysis, we analysed eight microsatellite loci from 484 tissue samples of larvae and adults of Columbia spotted frogs (Rana luteiventris) and long-toed salamanders (Ambystoma macrodactylum) at nine breeding sites in north Idaho. Differences in allele frequencies between adult and larval samples were not detected after full-siblings were removed from the larval data set for Columbia spotted frogs; for long-toed salamanders, these differences remained at two out of four ponds. Data from Columbia spotted frog larvae indicated higher levels of differentiation among populations (median difference in F_ST = 0.020, P < 0.01), as predicted by population genetic theory, whereas data from larval samples of long-toed salamanders showed some evidence of lower levels of differentiation among populations (median difference in F_ST = 0.012, P = 0.06). For both species, removing all but one individual from each full-sibling family led to parameter estimates that were closer to those calculated from adult samples for both population and landscape genetic measures. Removal of full-siblings is likely to improve estimates of population genetic parameters; however, knowledge of the species’ breeding system is essential for understanding additional sources of bias when inferring population genetic structure from larval samples.     

Sex ratio rather than population size affects genetic diversity in Antennaria dioica

• Habitat fragmentation and small population size can lead to genetic erosion in threatened plant populations. Classical theory implies that dioecy can counteract genetic erosion as it decreases the magnitude of inbreeding and genetic drift due to obligate outcrossing. However, in small populations, sex ratios may be strongly male‐ or female‐biased, leading to substantial reductions in effective population size. This may theoretically result in a unimodal relationship between sex ratios and genetic diversity; yet, empirical studies on this relationship are scarce.
• Using AFLP markers, we studied genetic diversity, structure and differentiation in 14 highly fragmented Antennaria dioica populations from the Central European lowlands. Our analyses focused on the relationship between sex ratio, population size and genetic diversity.
• Although most populations were small (mean: 35.5 patches), genetic diversity was moderately high. We found evidence for isolation‐by‐distance, but overall differentiation of the populations was rather weak. Females dominated 11 populations, which overall resulted in a slightly female‐biased sex ratio (61.5%). There was no significant relationship between population size and genetic diversity. The proportion of females was not unimodally but positively linearly related to genetic diversity.
• The high genetic diversity and low genetic differentiation suggest that A. dioica has been widely distributed in the Central European lowlands in the past, while fragmentation occurred only in the last decades. Sex ratio has more immediate consequences on genetic diversity than population size. An increasing proportion of females can increase genetic diversity in dioecious plants, probably due to a higher amount of sexual reproduction.

     

Thyroid hormone responsive QTL and the evolution of paedomorphic salamanders

The transformation of ancestral phenotypes into novel traits is poorly understood for many examples of evolutionary novelty. Ancestrally, salamanders have a biphasic life cycle with an aquatic larval stage, a brief and pronounced metamorphosis, followed by a terrestrial adult stage. Repeatedly during evolution, metamorphic timing has been delayed to exploit growth-permissive environments, resulting in paedomorphic salamanders that retain larval traits as adults. We used thyroid hormone (TH) to rescue metamorphic phenotypes in paedomorphic salamanders and then identified quantitative trait loci (QTL) for life history traits that are associated with amphibian life cycle evolution: metamorphic timing and adult body size. We demonstrate that paedomorphic tiger salamanders (Ambystoma tigrinum complex) carry alleles at three moderate effect QTL (met1–3) that vary in responsiveness to TH and additively affect metamorphic timing. Salamanders that delay metamorphosis attain significantly larger body sizes as adults and met2 explains a significant portion of this variation. Thus, substitution of alleles at TH-responsive loci suggests an adaptive pleiotropic basis for two key life-history traits in amphibians: body size and metamorphic timing. Our study demonstrates a likely pathway for the evolution of novel paedomorphic species from metamorphic ancestors via selection of TH-response alleles that delay metamorphic timing and increase adult body size.     

Genetic diversity and structure in two species of Leavenworthia with self-incompatible and self-compatible populations

Self-fertilization is a common mating system in plants and is known to reduce genetic diversity, increase genetic structure and potentially put populations at greater risk of extinction. In this study, we measured the genetic diversity and structure of two cedar glade endemic species, Leavenworthia alabamica and L. crassa. These species have self-incompatible (SI) and self-compatible (SC) populations and are therefore ideal for understanding how the mating system affects genetic diversity and structure. We found that L. alabamica and L. crassa had high species-level genetic diversity (H_e=0.229 and 0.183, respectively) and high genetic structure among their populations (F_ST=0.45 and 0.36, respectively), but that mean genetic diversity was significantly lower in SC compared with SI populations (SC vs SI, H_e for L. alabamica was 0.065 vs 0.206 and for L. crassa was 0.084 vs 0.189). We also found significant genetic structure using maximum-likelihood clustering methods. These data indicate that the loss of SI leads to the loss of genetic diversity within populations. In addition, we examined genetic distance relationships between SI and SC populations to analyze possible population history and origins of self-compatibility. We find there may have been multiple origins of self-compatibility in L. alabamica and L. crassa. However, further work is required to test this hypothesis. Finally, given their high genetic structure and that individual populations harbor unique alleles, conservation strategies seeking to maximize species-level genetic diversity for these or similar species should protect multiple populations.     

Multilocus heterozygosity, parental relatedness and individual fitness components in a wild mountain goat, Oreamnos americanus population

Matings between relatives lead to a decrease in offspring genetic diversity which can reduce fitness, a phenomenon known as inbreeding depression. Because alpine ungulates generally live in small structured populations and often exhibit a polygynous mating system, they are susceptible to inbreeding. Here, we used marker-based measures of pairwise genetic relatedness and inbreeding to investigate the fitness consequences of matings between relatives in a long-term study population of mountain goats ( Oreamnos americanus ) at Caw Ridge, Alberta, Canada. We first assessed whether individuals avoided mating with kin by comparing actual and random mating pairs according to their estimated genetic relatedness, which was derived from 25 unlinked polymorphic microsatellite markers and reflected pedigree relatedness. We then examined whether individual multilocus heterozygosity H , used as a measure of inbreeding, was predicted by parental relatedness and associated with yearling survival and the annual probability of giving birth to a kid in adult females. Breeding pairs identified by genetic parentage analyses of offspring that survived to 1 year of age were less genetically related than expected under random matings. Parental relatedness was negatively correlated with offspring H , and more heterozygous yearlings had higher survival to 2 years of age. The probability of giving birth was not affected by H in adult females. Because kids that survived to yearling age were mainly produced by less genetically related parents, our results suggest that some individuals experienced inbreeding depression in early life. Future research will be required to quantify the levels of gene flow between different herds, and evaluate their effects on population genetic diversity and dynamics.     

Diversity array technology markers: genetic diversity analyses and linkage map construction in rapeseed (Brassica napus L.)

We developed Diversity Array Technology (DArT) markers for application in genetic studies of Brassica napus and other Brassica species with A or C genomes. Genomic representation from 107 diverse genotypes of B. napus L. var. oleifera (rapeseed, AACC genomes) and B. rapa (AA genome) was used to develop a DArT array comprising 11 520 clones generated using PstI/BanII and PstI/BstN1 complexity reduction methods. In total, 1547 polymorphic DArT markers of high technical quality were identified and used to assess molecular diversity among 89 accessions of B. napus, B. rapa, B. juncea, and B. carinata collected from different parts of the world. Hierarchical cluster and principal component analyses based on genetic distance matrices identified distinct populations clustering mainly according to their origin/pedigrees. DArT markers were also mapped in a new doubled haploid population comprising 131 lines from a cross between spring rapeseed lines 'Lynx-037DH' and 'Monty-028DH'. Linkage groups were assigned on the basis of previously mapped simple sequence repeat (SSRs), intron polymorphism (IP), and gene-based markers. The map consisted of 437 DArT, 135 SSR, 6 IP, and 6 gene-based markers and spanned 2288 cM. Our results demonstrate that DArT markers are suitable for genetic diversity analysis and linkage map construction in rapeseed.     

Genetic diversity and structure of an endemic and critically endangered stream river salamander (Caudata: Ambystoma leorae) in Mexico

Small or isolated populations are highly susceptible to stochastic events. They are prone and vulnerable to random demographic or environmental fluctuations that could lead to extinction due to the loss of alleles through genetic drift and increased inbreeding. We studied Ambystoma leorae an endemic and critically threatened species. We analyzed the genetic diversity and structure, effective population size, presence of bottlenecks and inbreeding coefficient of 96 individuals based on nine microsatellite loci. We found high levels of genetic diversity expressed as heterozygosity (H_o = 0.804, H_e = 0.613, H_e* = 0.626 and H_Nei = 0.622). The population presents few alleles (4–9 per locus) and genotypes (3–14 per locus) compared with other mole salamanders species. We identified three genetically differentiated subpopulations with a significant level of genetic structure (F_ST = 0.021, R_ST = 0.044 y D_est = 0.010, 95 % CI). We also detected a reduction signal in population size and evidence of a genetic bottleneck (M = 0.367). The effective population size is small (Ne = 45.2), but similar to another mole salamanders with restricted distributions or with recently fragmented habitat. The inbreeding coefficient levels detected are low (F_IS = ?0.619–0.102) as is gene flow. Despite, high levels of genetic diversity A. leorae is critically endangered because it is a small isolated population.