A decade of amphibian population genetic studies: synthesis and recommendations

Genomic insights into the Montseny brook newt (Calotriton arnoldi), a Critically Endangered glacial relict

The Montseny brook newt (Calotriton arnoldi), considered the most endangered amphibian in Europe, is a relict salamandrid species endemic to a small massif located in northeastern Spain. Although conservation efforts should always be guided by genomic studies, those are yet scarce among urodeles, hampered by the extreme sizes of their genomes. Here, we present the third available genome assembly for the order Caudata, and the first genomic study of the species and its sister taxon, the Pyrenean brook newt (Calotriton asper), combining whole-genome and ddRADseq data. Our results reveal significant demographic oscillations which accurately mirrored Europe's climatic history. Although severe bottlenecks have led to depauperate genomic diversity and long runs of homozygosity along a gigantic genome, inbreeding might have been avoided by assortative mating strategies. Other life history traits, however, seem to have been less advantageous, and the lack of land dispersal has driven to exceptional levels of population fragmentation.

Fine-scale spatial genetic structure and dispersal among Italian smooth newt populations in a rural landscape

Amphibians are particularly sensitive to habitat loss and fragmentation caused by the intensification and modernization of farming occurring in the second half of the twentieth century in the Mediterranean basin. However, artificial water bodies, associated with traditional husbandry, proved to be important surrogate for amphibian feeding and reproduction. Here, multilocus genotypes were used to investigate the spatial population structure of Lissotriton vulgaris meridionalis and the role of drinking troughs in supporting viable breeding populations within a rural landscape interested by traditional husbandry and agriculture. Our genetic analysis highlighted the conservation value and the potential stepping-stone function of artificial aquatic sites in the dispersal of the species and for the gene flow maintenance. Indeed, populations of drinking troughs show allelic richness and heterozygosity levels comparable to those from natural ponds and there is no great evidence of genetic bottlenecks. A complex system of artificial aquatic sites and few natural wetlands was identified sustaining a well-structured network of demes highly interconnected with themselves and natural aquatic sites. The conservation of the identified genetic clusters may be useful to prevent further population declines and future loss of genetic diversity within the study area characterized by scarce natural wetlands that frequently dried because of agricultural practices and strong seasonality. Site-specific protection measures are needed to contrast the progressive disappearance of drinking troughs observed in the last years in Italy because of the abandonment of traditional farming practices in favour of modern agriculture and intensive farming.

Population genetic structure of the Natterjack toad (Epidalea calamita) in Ireland: implications for conservation management

Molecular methods can play a crucial role in species management and conservation. Despite the usefulness of genetic approaches, they are often not explicitly included as part of species recovery plans and conservation practises. The Natterjack toad (Epidalea calamita) is regionally Red-Listed as Endangered in Ireland. The species is declining and is now present at just seven sites within a highly restricted range. This study used 13 highly polymorphic microsatellite markers to analyse the population genetic diversity and structure. Genetic diversity was high with expected heterozygosity between 0.55 and 0.61 and allelic richness between 4.77 and 5.92. Effective population sizes were small (Ne < 100 individuals), but not abnormal for pond breeding amphibians. However, there was no evidence of historical or contemporary genetic bottlenecks or high levels of inbreeding. We identified a positive relationship between Ne and breeding pond surface area, suggesting that environmental factors are a key determinant of population size. Significant genetic structuring was detected throughout the species’ range, and we identified four genetic entities that should be considered in the species’ conservation strategies. Management should focus on preventing further population declines and future loss of genetic diversity overall and within genetic entities while maintaining adequate local effective population size through site-specific protection, human-mediated translocations and head-start programs. The apparent high levels of genetic variation give hope for the conservation of Ireland’s rarest amphibian if appropriately protected and managed.

Predicting amphibian intraspecific diversity with machine learning: Challenges and prospects for integrating traits, geography, and genetic data

The growing availability of genetic data sets, in combination with machine learning frameworks, offers great potential to answer long-standing questions in ecology and evolution. One such question has intrigued population geneticists, biogeographers, and conservation biologists: What factors determine intraspecific genetic diversity? This question is challenging to answer because many factors may influence genetic variation, including life history traits, historical influences, and geography, and the relative importance of these factors varies across taxonomic and geographic scales. Furthermore, interpreting the influence of numerous, potentially correlated variables is difficult with traditional statistical approaches. To address these challenges, we analysed repurposed data using machine learning and investigated predictors of genetic diversity, focusing on Nearctic amphibians as a case study. We aggregated species traits, range characteristics, and >42,000 genetic sequences for 299 species using open-access scripts and various databases. After identifying important predictors of nucleotide diversity with random forest regression, we conducted follow-up analyses to examine the roles of phylogenetic history, geography, and demographic processes on intraspecific diversity. Although life history traits were not important predictors for this data set, we found significant phylogenetic signal in genetic diversity within amphibians. We also found that salamander species at northern latitudes contained low genetic diversity. Data repurposing and machine learning provide valuable tools for detecting patterns with relevance for conservation, but concerted efforts are needed to compile meaningful data sets with greater utility for understanding global biodiversity.

Hybridizing salamanders experience accelerated diversification

Whether hybridization generates or erodes species diversity has long been debated, but to date most studies have been conducted at small taxonomic scales. Salamanders (order Caudata) represent a taxonomic order in which hybridization plays a prevalent ecological and evolutionary role. We employed a recently developed model of trait-dependent diversification to test the hypothesis that hybridization impacts the diversification dynamics of species that are currently hybridizing. We find strong evidence supporting this hypothesis, showing that hybridizing salamander lineages have significantly greater net-diversification rates than non-hybridizing lineages. This pattern is driven by concurrently increased speciation rates and decreased extinction rates in hybridizing lineages. Our results support the hypothesis that hybridization can act as a generative force in macroevolutionary diversification.

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.

Systematic review on the conservation genetics of African savannah elephants

Background

In this paper we review the conservation genetics of African savannah elephants, aiming to understand the spatio-temporal research trends and their underlying factors. As such, we explore three questions associated to the conservation genetics and molecular ecology of these elephants: (1) what are the research trends concerning the conservation genetics of Loxodonta africana? (2) Do richer countries conduct more research on the genetics of African elephants? (3) Which attributes influence where scholars conduct their research?

Materials and Methods

We examined available peer-reviewed publications from 1993 to 2014 in complementary online databases, including the ISI/Web of Science (WoS), Scopus and Google Scholar (GS), and searched for publications in scientific journals as well as in the reference section of these publications. We analyzed the annual trend of publications in this field of research, including the number of authors, levels of collaboration among authors, year of publication, publishing journal and the countries from where genetic samples were collected. Additionally, we identified main research clusters, authors, and institutional collaborations, based on co-citation and co-occurrence networks.

Results

We found that during the study period there was a positive trend in the number of publications and a reduction in the number of authors per paper. Twenty-five countries contributed, with the majority of publications authored by researchers in the USA, Kenya and South Africa. The majority of samples were collected in Kenya, Tanzania and South Africa. Research outputs are associated with the existence of long-term conservation/research projects and research potential as measured by the literacy rate and the number of higher education institutions in a country. Five research clusters were identified, focusing on the origin and evolution of the species, methodological issues and the relatedness among elephant species.

Conclusions

Research in this field should be expanded to additional countries harboring elephant populations to enable a more comprehensive understanding of the population structure and genetic differentiation of the species, and to cope with challenges associated with the conservation of the species such as illegal hunting, habitat fragmentation, species reintroduction and climate change.

Patterns of Genetic Variability in Island Populations of the Cane Toad (Rhinella marina) from the Mouth of the Amazon

The Amazonian coast has several unique geological characteristics resulting from the interaction between drainage pattern of the Amazon River and the Atlantic Ocean. It is one of the most extensive and sedimentologically dynamic regions of the world, with a large number of continental islands mostly formed less than 10,000 years ago. The natural distribution of the cane toad (Rhinella marina), one of the world's most successful invasive species, in this complex Amazonian system provides an intriguing model for the investigation of the effects of isolation or the combined effects of isolation and habitat dynamic changes on patterns of genetic variability and population differentiation. We used nine fast-evolving microsatellite loci to contrast patterns of genetic variability in six coastal (three mainlands and three islands) populations of the cane toad near the mouth of the Amazon River. Results from Bayesian multilocus clustering approach and Discriminant Analyses of Principal Component were congruent in showing that each island population was genetically differentiated from the mainland populations. All FST values obtained from all pairwise comparisons were significant, ranging from 0.048 to 0.186. Estimates of both recent and historical gene flow were not significantly different from zero across all population pairs, except the two mainland populations inhabiting continuous habitats. Patterns of population differentiation, with a high level of population substructure and absence/restricted gene flow, suggested that island populations of R. marina are likely isolated since the Holocene sea-level rise. However, considering the similar levels of genetic variability found in both island and mainland populations, it is reliable to assume that they were also isolated for longer periods. Given the genetic uniqueness of each cane toad population, together with the high natural vulnerability of the coastal regions and intense human pressures, we suggest that these populations should be treated as discrete units for conservation management purposes.

A European Concern? Genetic Structure and Expansion of Golden Jackals (Canis aureus) in Europe and the Caucasus

In the first continent-wide study of the golden jackal (Canis aureus), we characterised its population genetic structure and attempted to identify the origin of European populations. This provided a unique insight into genetic characteristics of a native carnivore population with rapid large-scale expansion. We analysed 15 microsatellite markers and a 406 base-pair fragment of the mitochondrial control region. Bayesian-based and principal components methods were applied to evaluate whether the geographical grouping of samples corresponded with genetic groups. Our analysis revealed low levels of genetic diversity, reflecting the unique history of the golden jackal among Europe’s native carnivores. The results suggest ongoing gene flow between south-eastern Europe and the Caucasus, with both contributing to the Baltic population, which appeared only recently. The population from the Peloponnese Peninsula in southern Greece forms a common genetic cluster with samples from south-eastern Europe (ΔK approach in STRUCTURE, Principal Components Analysis [PCA]), although the results based on BAPS and the estimated likelihood in STRUCTURE indicate that Peloponnesian jackals may represent a distinct population. Moreover, analyses of population structure also suggest either genetic distinctiveness of the island population from Samos near the coast of Asia Minor (BAPS, most STRUCTURE, PCA), or possibly its connection with the Caucasus population (one analysis in STRUCTURE). We speculate from our results that ancient Mediterranean jackal populations have persisted to the present day, and have merged with jackals colonising from Asia. These data also suggest that new populations of the golden jackal may be founded by long-distance dispersal, and thus should not be treated as an invasive alien species, i.e. an organism that is “non-native to an ecosystem, and which may cause economic or environmental harm or adversely affect human health”. These insights into the genetic structure and ancestry of Baltic jackals have important implications for management and conservation of jackals in Europe. The golden jackal is listed as an Annex V species in the EU Habitats Directive and as such, considering also the results presented here, should be legally protected in all EU member states.

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.

The Phylogeography and Population Demography of the Yunnan Caecilian (Ichthyophis bannanicus): Massive Rivers as Barriers to Gene Flow

Ichthyophis bannanicus is the only caecilian species in China. In this study, the phylogeography and population demography of I. bannanicus were explored, based on the mitochondrial DNA genes (cyt b and ND2) and 15 polymorphic microsatellite loci. Altogether 158 individuals were collected from five populations in Yunnan province, Guangxi province, Guangdong province, and Northern Vietnam. Phylogeographical and population structure analysis identified either two groups (Xishuangbanna, Northern Vietnam-Yulin-Yangchun-Deqing) or three groups (Xishuangbanna, Northern Vietnam-Yulin-Yangchun, and Deqing), indicating that the Red River and Pearl River systems may have acted as gene-flow barriers for I. bannanicus. Historical population expansion that happened 15–17 Ka ago was detected for mtDNA data and was possibly triggered by warmer weather after the Last Glacial Maximum. However, the Bayesian simulations of population history based on microsatellite data pinpointed population decline in all populations since 19,123 to 1,029 years ago, demonstrating a significant influence of anthropogenic habitat alteration on I. bannanicus.

Conflicting effects of microhabitats on Long-toed Salamander (Ambystoma macrodactylum) movement: implications for landscape connectivity

Understanding dispersal requires multiple lines of investigation, from the study of broad patterns of population connectivity to the identification of factors impacting movement at local scales. To determine the potential effects of different microhabitats on dispersal in the Long-toed Salamander (Ambystoma macrodactylum Baird, 1850), we experimentally evaluated mobility, moisture loss, and habitat choice in response to five common substrates (deciduous and coniferous leaf litter, grass, moss, and sand). Specifically, we examined differences in the efficiency with which salamanders moved across substrates when motivated to move. We then quantified moisture loss in each substrate and evaluated habitat preference. Our results point to a trade-off between substrates that are easily traversed and those that offer high protection against desiccation. Habitat choice appeared to balance these two aspects of performance, with salamanders favouring a substrate that offered both low resistance to movement and high protection against desiccation. This result was context-dependent, as preferences shifted towards wetter but less easily traversed substrates when supplemental cover objects were made available. Overall, our study highlights the potential for individuals to respond to a given substrate in ways that can both facilitate and limit dispersal and thus underscores the need to consider different aspects of individual performance and behaviour when studying population connectivity.

Urban park characteristics, genetic variation, and historical demography of white-footed mouse (Peromyscus leucopus) populations in New York City

Severe fragmentation is a typical fate of native remnant habitats in cities, and urban wildlife with limited dispersal ability are predicted to lose genetic variation in isolated urban patches. However, little information exists on the characteristics of urban green spaces required to conserve genetic variation. In this study, we examine whether isolation in New York City (NYC) parks results in genetic bottlenecks in white-footed mice (Peromyscus leucopus), and test the hypotheses that park size and time since isolation are associated with genetic variability using nonlinear regression and information-theoretic model selection. White-footed mice have previously been documented to exhibit male-biased dispersal, which may create disparities in genetic variation between males and females in urban parks. We use genotypes of 18 neutral microsatellite data and four different statistical tests to assess this prediction. Given that sex-biased dispersal may create disparities between population genetic patterns inferred from bi- vs. uni-parentally inherited markers, we also sequenced a 324 bp segment of the mitochondrial D-loop for independent inferences of historical demography in urban P. leucopus. We report that isolation in urban parks does not necessarily result in genetic bottlenecks; only three out of 14 populations in NYC parks exhibited a signature of a recent bottleneck at 18 neutral microsatellite loci. Mouse populations in larger urban parks, or parks that have been isolated for shorter periods of time, also do not generally contain greater genetic variation than populations in smaller parks. These results suggest that even small networks of green spaces may be sufficient to maintain the evolutionary potential of native species with certain characteristics. We also found that isolation in urban parks results in weak to nonexistent sex-biased dispersal in a species known to exhibit male-biased dispersal in less fragmented environments. In contrast to nuclear loci, mitochondrial D-loop haplotypes exhibited a mutational pattern of demographic expansion after a recent bottleneck or selective sweep. Estimates of the timing of this expansion suggest that it occurred concurrent with urbanization of NYC over the last few dozens to hundreds of years. Given the general non-neutrality of mtDNA in many systems and evidence of selection on related coding sequences in urban P. leucopus, we argue that the P. leucopus mitochondrial genome experienced recent negative selection against haplotypes not favored in isolated urban parks. In general, rapid adaptive evolution driven by urbanization, global climate change, and other human-caused factors is underappreciated by evolutionary biologists, but many more cases will likely be documented in the near future.

First report of Hepatozoon (Apicomplexa: Adeleorina) in caecilians, with description of a new species

Hepatozoon spp. are identified for the first time in the amphibian order Gymnophiona, or caecilians, from the Seychelles island of Silhouette. Estimate of relationships derived from partial 18S rRNA gene sequences indicate these are not related to Hepatozoon spp. from frogs or to other Hepatozoon spp. from reptiles in the Seychelles. Assessment of mature gamonts from blood smears indicate that these can be recognized as a new species, Hepatozoon seychellensis n. sp.

Conservation genetics of extremely isolated urban populations of the northern dusky salamander (Desmognathus fuscus) in New York City

Urbanization is a major cause of amphibian decline. Stream-dwelling plethodontid salamanders are particularly susceptible to urbanization due to declining water quality and hydrological changes, but few studies have examined these taxa in cities. The northern dusky salamander (Desmognathus fuscus) was once common in the New York City metropolitan area, but has substantially declined throughout the region in recent decades. We used five tetranucleotide microsatellite loci to examine population differentiation, genetic variation, and bottlenecks among five remnant urban populations of dusky salamanders in NYC. These genetic measures provide information on isolation, prevalence of inbreeding, long-term prospects for population persistence, and potential for evolutionary responses to future environmental change. All populations were genetically differentiated from each other, and the most isolated populations in Manhattan have maintained very little genetic variation (i.e. <20% heterozygosity). A majority of the populations also exhibited evidence of genetic bottlenecks. These findings contrast with published estimates of high genetic variation within and lack of structure between populations of other desmognathine salamanders sampled over similar or larger spatial scales. Declines in genetic variation likely resulted from population extirpations and the degradation of stream and terrestrial paths for dispersal in NYC. Loss of genetic variability in populations isolated by human development may be an underappreciated cause and/or consequence of the decline of this species in urbanized areas of the northeast USA.