Population structure and conservation genetics of the Oregon spotted frog, Rana pretiosa

Importance of Genetic-Fitness Correlations for the Conservation of Amphibians

Endangered animals suffer from isolation of their habitats. Isolation leads to a reduction in population size as well as a decrease in genetic diversity and a concomitant increase in the risk of extinction. Amphibians are the most endangered vertebrate class. Besides habitat loss, fragmentation and isolation, amphibians are threatened by emerging diseases e.g., chytrid fungus or Ranavirus. By employing experiments, researchers investigate whether changes in genetic diversity within or among isolated populations affect amphibian fitness. While genetic diversity estimates are based on molecular markers, typically microsatellites, fitness is mostly measured as tadpole performance in rearing experiments often under varying environmental conditions. Tadpole performances (e.g., body mass, growth rate and survival) have been found to be negatively affected by low genetic diversity, as several studies have found a positive association between genetic diversity and these fitness traits. Moreover, infection with pathogens also seems to be more likely in individuals or populations with lower genetic diversity. Overall, these genetic-fitness correlations seem to be more pronounced or detectable in smaller, declining populations but not in larger populations. Genomic studies, which sample a larger fraction of the genome, are still scarce in the conservation genetic literature on amphibians. These are likely to increase in upcoming years and may reveal adaptive variants that protect against dangerous pathogens or environmental changes. Altogether, genetic-fitness correlation studies should be a priority in order to develop effective management plans for the genetic rescue of isolated, imperilled amphibian populations.

Extremely low mtDNA diversity and high genetic differentiation reveal the precarious genetic status of dugongs in New Caledonia, South Pacific

New Caledonia is home to one of the largest remaining populations of dugongs (Dugong dugon) and is located at the southeastern limit of the species' range. Local knowledge suggests that current levels of removal due to anthropogenic pressures are unsustainable, while trends suggest an ongoing decline in the population. Considering this unfavorable conservation context, this study aimed to assess the New Caledonian dugong population's resilience by determining its level of genetic diversity and degree of isolation relative to other populations. Mitochondrial DNA (mtDNA) control region sequences (n = 55) collected from live and dead dugongs in New Caledonia were compared to a global dataset of previously published sequences (n = 631) representing dugong populations throughout the species range. The New Caledonian dugong population displayed the lowest level of mtDNA diversity documented worldwide (3 haplotypes with 1 base pair difference), suggesting a recent origin of the current population through limited colonization events. Population structure analyses indicate a strong genetic differentiation with all the putative populations represented in the global dataset, including large neighbouring Australian populations. These results show that the dugong population in New Caledonia is particularly isolated, fragile, and vulnerable to anthropogenic threats and diseases with low potential for resilience through incoming gene flow. Our findings call for an instant conservation response and consideration for IUCN population assessment to support the long-term survival of the New Caledonian dugong population.

Conservation genetics of relict tropical species of Magnolia (section Macrophylla)

Special conservation efforts should be made for relict species, as they usually have small population sizes and restricted distributions, placing them in critical extinction risk. To achieve conservation, information about genetic diversity distribution is needed. Here, using nine nuclear microsatellites, we analyzed 23 populations of five recently described species of Magnolia distributed in Mexico, which were previously assigned to Magnolia dealbata. We aimed to determine the level of genetic diversity and the distribution of genetic variation and proposed conservation measures. Compared to other endemic and relict species, we found a moderate level of genetic diversity in most populations; however, we identified two populations with no genetic variation. Additionally, we found evidence of positive values of inbreeding likely due to geitonogamy. We found a strong population structure, low effective population size, and no evidence of bottlenecks. Patterns of genetic differentiation did not support the morphological distinction of five species, so we hypothesized that the gene pools may instead represent well-differentiated populations of a single species. We argue that the pattern of genetic differentiation is explained by the natural fragmentation of the cloud forests after glaciation events, and the effects of genetic drift in small populations poorly connected by gene flow. Despite the moderate levels of genetic diversity, special attention is needed to guarantee conservation, with emphasis on the populations in the central region of the country as well as the valuable populations identified in the southwestern region.

Conservation genetics of two threatened frogs from the Mambilla highlands, Nigeria

Amphibians are the vertebrate group with the highest number of species threatened with extinction, and habitat loss and fragmentation are considered to be among the leading causes of their declines and extinctions. Little is known of the population biology of amphibian species inhabiting montane forests in Central and West Africa, where anthropogenic activities such as farming and cattle raising are major threats to native biodiversity. We used Amplified Fragment Length Polymorphisms (AFLPs) to assess the population genetic structure of two poorly known species, Cardioglossa schioetzi and Leptodactylodon bicolor (both in the Arthroleptidae), in and around Ngel Nyaki Forest Reserve on the Mambilla Plateau in eastern Nigeria. The landscape comprises continuous forest on steep slopes and small riparian forest fragments in a grassland matrix. While increased fragmentation is well documented for these and other forests in the mountains of Cameroon and Nigeria over the past century, there are no previous assessments of the impact of forest fragmentation on montane amphibian populations in this region. Our estimates of genetic diversity are similar across populations within each species with levels of heterozygosity values consistent with local population declines. Except for a pair of populations (C. schioetzi) we did not observe genetic differentiation between forest and riparian forest fragment populations, nor across sites within continuous forest (L. bicolor). Our results demonstrate recent gene flow between forest fragments and the adjacent protected forests and suggest that small forest corridors connecting these may lessen the genetic consequences of at least 30 years of intense and severe fragmentation in Ngel Nyaki.

Regional variation in drivers of connectivity for two frog species (Rana pretiosa and R. luteiventris) from the U.S. Pacific Northwest

Comparative landscape genetics has uncovered high levels of variability in which landscape factors affect connectivity among species and regions. However, the relative importance of species traits versus environmental variation for predicting landscape patterns of connectivity is unresolved. We provide evidence from a landscape genetics study of two sister taxa of frogs, the Oregon spotted frog (Rana pretiosa) and the Columbia spotted frog (Rana luteiventris) in Oregon and Idaho, USA. Rana pretiosa is relatively more dependent on moisture for dispersal than R. luteiventris, so if species traits influence connectivity, we predicted that connectivity among R. pretiosa populations would be more positively associated with moisture than R. luteiventris. However, if environmental differences are important drivers of gene flow, we predicted that connectivity would be more positively related to moisture in arid regions. We tested these predictions using eight microsatellite loci and gravity models in two R. pretiosa regions and four R. luteiventris regions (n = 1,168 frogs). In R. pretiosa, but not R. luteiventris, connectivity was positively related to mean annual precipitation, supporting our first prediction. In contrast, connectivity was not more positively related to moisture in more arid regions. Various temperature metrics were important predictors for both species and in all regions, but the directionality of their effects varied. Therefore, the pattern of variation in drivers of connectivity was consistent with predictions based on species traits rather than on environmental variation.

No signs of inbreeding despite long-term isolation and habitat fragmentation in the critically endangered Montseny brook newt (Calotriton arnoldi)

Endemic species with restricted geographic ranges potentially suffer the highest risk of extinction. If these species are further fragmented into genetically isolated subpopulations, the risk of extinction is elevated. Habitat fragmentation is generally considered to have negative effects on species survival, despite some evidence for neutral or even positive effects. Typically, non-negative effects are ignored by conservation biology. The Montseny brook newt (Calotriton arnoldi) has one of the smallest distribution ranges of any European amphibian (8 km²) and is considered critically endangered by the International Union for Conservation of Nature. Here we apply molecular markers to analyze its population structure and find that habitat fragmentation owing to a natural barrier has resulted in strong genetic division of populations into two sectors, with no detectable migration between sites. Although effective population size estimates suggest low values for all populations, we found low levels of inbreeding and relatedness between individuals within populations. Moreover, C. arnoldi displays similar levels of genetic diversity to its sister species Calotriton asper, from which it separated around 1.5 million years ago and which has a much larger distribution range. Our extensive study shows that natural habitat fragmentation does not result in negative genetic effects, such as the loss of genetic diversity and inbreeding on an evolutionary timescale. We hypothesize that species in such conditions may evolve strategies (for example, special mating preferences) to mitigate the effects of small population sizes. However, it should be stressed that the influence of natural habitat fragmentation on an evolutionary timescale should not be conflated with anthropogenic habitat loss or degradation when considering conservation strategies.

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.

Significant population genetic structuring but a lack of phylogeographic structuring in the endemic Tasmanian tree frog (Litoria burrowsae)

Conservation of frogs is of global concern, owing to declines resulting from habitat destruction, global climate change, and disease. Knowledge of genetic variation in frog species is therefore desirable for the identification of management units. Here we surveyed mitochondrial DNA sequence variation in the Tasmanian endemic hylid frog Litoria burrowsae, which is infected by chytrid fungus, Batrachochytrium dendrobatidis, and may be declining. Neither phylogeographic structure nor deep phylogenetic divergence was detected in the species, although its populations were highly differentiated with respect to haplotype frequencies. The low-haplotype diversity in L. burrowsae suggests a recent bottleneck in the species, and population genetic structuring may reflect isolation by distance as well as founder effects associated with range expansion. Three putative management units were identified that require verification based on nuclear DNA variation and adaptation to local environments.

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.

Extensive genetic diversity and substructuring among zebrafish strains revealed through copy number variant analysis

Copy number variants (CNVs) represent a substantial source of genomic variation in vertebrates and have been associated with numerous human diseases. Despite this, the extent of CNVs in the zebrafish, an important model for human disease, remains unknown. Using 80 zebrafish genomes, representing three commonly used laboratory strains and one native population, we constructed a genome-wide, high-resolution CNV map for the zebrafish comprising 6,080 CNV elements and encompassing 14.6% of the zebrafish reference genome. This amount of copy number variation is four times that previously observed in other vertebrates, including humans. Moreover, 69% of the CNV elements exhibited strain specificity, with the highest number observed for Tubingen. This variation likely arose, in part, from Tubingen's large founding size and composite population origin. Additional population genetic studies also provided important insight into the origins and substructure of these commonly used laboratory strains. This extensive variation among and within zebrafish strains may have functional effects that impact phenotype and, if not properly addressed, such extensive levels of germ-line variation and population substructure in this commonly used model organism can potentially confound studies intended for translation to human diseases.

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.

More than meets the eye: detecting cryptic microgeographic population structure in a parasite with a complex life cycle

Nonrandom recruitment of parasites among hosts can lead to genetic differentiation among hosts and mating dynamics that promote inbreeding. It has been hypothesized that strictly aquatic parasites with intermediate hosts will behave as panmictic populations among hosts because ample opportunity exists for random mixing of unrelated individuals during transmission to the definitive host. A previous allozyme study on the marine trematode Lecithochirium fusiforme did not support this hypothesis; in that, there was genetic differentiation among, and significant heterozygote deficiencies within, definitive hosts. We revisit this system and use microsatellites to obtain multilocus genotypes. Our goal was to determine whether cryptic subgroups and/or the presence of clones could account for the apparent deviation from 'panmixia'. We find strong evidence for cryptic subdivision (three genetic clusters) that causes the Wahlund effect and differentiation among definitive hosts. After accounting for these cryptic groups, we see panmictic genetic structure among definitive hosts that is consistent with the 'high mixing in aquatic habitats' hypothesis. We see evidence for cotransmission of clones in all three clusters, but this level of clonal structure did not have a major impact in causing deviations from Hardy-Weinberg equilibrium, and only affected genetic differentiation among hosts in one cluster. A cursory examination of the data may have led to incorrect conclusions about nonrandom transmission. However, it is obvious in this system that there is more than meets the eye in relation to the actual make-up of parasite populations. In general, the methods we employ will be useful for elucidating hidden patterns in other organisms where cryptic structure may be common (e.g. those with limited morphology or complex life histories).