Conservation genetics and genomics of amphibians and reptiles

Evolution, Diversity, and Conservation of Herpetofauna

Amphibians and reptiles play a critical role in the evolution of Tetrapoda, showcasing significant diversity in terms of their genetics, species, morphology, life history traits, and evolutionary functions [...].

Losing the forest for the tree? On the wisdom of subpopulation management

Animal habitats are changing around the world in many ways, presenting challenges to the survival of species. Zoo animal populations are also challenged by small population sizes and limited genetic diversity. Some ex situ populations are managed as subpopulations based on presumed subspecies or geographic locality and related concerns over genetic purity or taxonomic integrity. However, these decisions can accelerate the loss of genetic diversity and increase the likelihood of population extinction. Here I challenge the wisdom of subpopulation management, pointing out significant concerns in the literature with delineation of species, subspecies, and evolutionarily significant units. I also review literature demonstrating the value of gene flow for preserving adaptive potential, the often-misunderstood role of hybridization in evolution, and the likely overstated concerns about outbreeding depression, and preservation of local adaptations. I argue that the most effective way to manage animal populations for the long term be they in human care, in the wild, or if a captive population is being managed for reintroduction, is to manage for maximum genetic diversity rather than managing subpopulations focusing on taxonomic integrity, genetic purity, or geographic locale because selection in the future, rather than the past, will determine what genotypes and phenotypes are the most fit. Several case studies are presented to challenge the wisdom of subpopulation management and stimulate thinking about the preservation of genomes rather than species, subspecies, or lineages because those units evolved in habitats that are likely very different from those habitats today and in the future.

Involvement of CD27 in innate and adaptive immunities of Nile tilapia (Oreochromis niloticus)

CD27 is a member of the TNF-receptor superfamily and plays various roles in immunities. However, the detailed information and mechanism of CD27 in bony fish immunity remain unclear. Therefore, in this research, certain interesting roles of CD27 in Nile tilapia (On-CD27) were determined. On-CD27 was largely expressed in the immune organs, head kidney, and spleen, and was sharply induced during bacterial infection. The in vitro tests suggested On-CD27 was involved in regulating inflammatory responses, activating immune-related signal pathways, and inducing apoptosis and pyroptosis progress. The scRNA data and in vivo experiments indicated that On-CD27 is mainly expressed in CD4⁺ T cells and involved in both innate and adaptive immunities. The present data provide a theoretical principle for further research on the mechanisms of CD27 in the innate and adaptive immunities of fish.

Development and characterization of twenty microsatellite markers for Phelsuma inexpectata (Squamata: Gekkonidae), a critically endangered gecko endemic to Reunion Island

BACKGROUND

The Manapany day gecko (Phelsuma inexpectata) is endemic to the south of Reunion Island. Threatened by habitat fragmentation and loss, human activities and invasive species, P. inexpectata is considered as critically endangered. Conservation measures are required but data on the species are missing, notably on its genetic diversity and population structure for which no specific markers are available to date. Here, we aimed to develop molecular markers to allow genetic studies of P. inexpectata.

METHODS AND RESULTS

We developed and characterized 20 polymorphic microsatellite markers based on 23 P. inexpectata individuals sampled from 10 sites. Then, the markers were tested on a total of 101 individuals, 30 from a natural site and 71 from an anthropized site. The mean values of N_a, H_o and H_e were 2.3 (± 0.2), 0.353 (± 0.053) and 0.345 (± 0.046) in the natural site and 2.8 (± 0.3), 0.345 (± 0.051) and 0.338 (± 0.048) in the anthropized site, respectively. Based on the combined loci, the probability of identity (PID) for unrelated specimens were 2.7 × 10^-7 and 2.6 × 10^-7 in the natural and anthropized site, respectively.

CONCLUSIONS

This work provides the first set of microsatellite markers for P. inexpectata, constituting a valuable tool to conduct classical genetic studies on the species, such as estimating genetic diversity, population structure and kinship relationships among individuals. Such studies will provide relevant information on P. inexpectata and will therefore be helpful in the implementation of conservation measures for this threatened species.

Population Subdivision in the Gopher Frog (Rana capito) across the Fragmented Longleaf Pine-Wiregrass Savanna of the Southeastern USA

Delineating genetically distinct population segments of threatened species and quantifying population connectivity are important steps in developing effective conservation and management strategies aimed at preventing extinction. The gopher frog (Rana capito) is a xeric-adapted, pond-breeding species endemic to the Gulf and Atlantic coastal plains of the southeastern United States. This species has experienced extensive habitat loss and fragmentation in the formerly widespread longleaf pine-wiregrass savanna where it lives, resulting in individual abundance declines and population extinctions throughout its range. We used individual-based clustering methods along with Bayesian inference of historical migration based on almost 1500 multilocus microsatellite genotypes to examine genetic structure in this taxon. Clustering analyses identified panhandle and peninsular populations in Florida as distinct genetic clusters separated by the Aucilla River, consistent with the division between the Coastal Plain and peninsular mitochondrial lineages, respectively. Analysis of historical migration indicated an east–west population divergence event followed by immigration to the east. Together, our results indicate that the genetically distinct Coastal Plain and peninsular Florida lineages should be considered separately for conservation and management purposes.

Dominance can increase genetic variance after a population bottleneck: a synthesis of the theoretical and empirical evidence

Drastic reductions in population size, or population bottlenecks, can lead to a reduction in additive genetic variance and adaptive potential. Genetic variance for some quantitative genetic traits, however, can increase after a population reduction. Empirical evaluations of quantitative traits following experimental bottlenecks indicate that non-additive genetic effects, including both allelic dominance at a given locus and epistatic interactions among loci, may impact the additive variance contributed by alleles that ultimately influences phenotypic expression and fitness. The dramatic effects of bottlenecks on overall genetic diversity have been well studied, but relatively little is known about how dominance and demographic events like bottlenecks can impact additive genetic variance. Herein, we critically examine how the degree of dominance among alleles affects additive genetic variance after a bottleneck. We first review and synthesize studies that document the impact of empirical bottlenecks on dominance variance. We then extend earlier work by elaborating on two theoretical models that illustrate the relationship between dominance and the potential increase in additive genetic variance immediately following a bottleneck. Furthermore, we investigate the parameters that influence the maximum level of genetic variation (associated with adaptive potential) after a bottleneck, including the number of founding individuals. Finally, we validated our methods using forward-time population genetic simulations of loci with varying dominance and selection levels. The fate of non-additive genetic variation following bottlenecks could have important implications for conservation and management efforts in a wide variety of taxa, and our work should help contextualize future studies (e.g., epistatic variance) in population genomics.

An integrative taxonomic revision of slug-eating snakes (Squamata: Pareidae: Pareineae) reveals unprecedented diversity in Indochina

Slug-eating snakes of the subfamily Pareinae are an insufficiently studied group of snakes specialized in feeding on terrestrial mollusks. Currently Pareinae encompass three genera with 34 species distributed across the Oriental biogeographic region. Despite the recent significant progress in understanding of Pareinae diversity, the subfamily remains taxonomically challenging. Here we present an updated phylogeny of the subfamily with a comprehensive taxon sampling including 30 currently recognized Pareinae species and several previously unknown candidate species and lineages. Phylogenetic analyses of mtDNA and nuDNA data supported the monophyly of the three genera Asthenodipsas, Aplopeltura, and Pareas. Within both Asthenodipsas and Pareas our analyses recovered deep differentiation with each genus being represented by two morphologically diagnosable clades, which we treat as subgenera. We further apply an integrative taxonomic approach, including analyses of molecular and morphological data, along with examination of available type materials, to address the longstanding taxonomic questions of the subgenus Pareas, and reveal the high level of hidden diversity of these snakes in Indochina. We restrict the distribution of P. carinatus to southern Southeast Asia, and recognize two subspecies within it, including one new subspecies proposed for the populations from Thailand and Myanmar. We further revalidate P. berdmorei, synonymize P. menglaensis with P. berdmorei, and recognize three subspecies within this taxon, including the new subspecies erected for the populations from Laos and Vietnam. Furthermore, we describe two new species of Pareas from Vietnam: one belonging to the P. carinatus group from southern Vietnam, and a new member of the P. nuchalis group from the central Vietnam. We provide new data on P. temporalis, and report on a significant range extension for P. nuchalis. Our phylogeny, along with molecular clock and ancestral area analyses, reveal a complex diversification pattern of Pareinae involving a high degree of sympatry of widespread and endemic species. Our analyses support the "upstream" colonization hypothesis and, thus, the Pareinae appears to have originated in Sundaland during the middle Eocene and then colonized mainland Asia in early Oligocene. Sundaland and Eastern Indochina appear to have played the key roles as the centers of Pareinae diversification. Our results reveal that both vicariance and dispersal are responsible for current distribution patterns of Pareinae, with tectonic movements, orogeny and paleoclimatic shifts being the probable drivers of diversification. Our study brings the total number of Pareidae species to 41 and further highlights the importance of comprehensive taxonomic revisions not only for the better understanding of biodiversity and its evolution, but also for the elaboration of adequate conservation actions.

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.

Land use and life history constrain adaptive genetic variation and reduce the capacity for climate change adaptation in turtles

BACKGROUND

Rapid anthropogenic climate change will require species to adapt to shifting environmental conditions, with successful adaptation dependent upon current patterns of genetic variation. While landscape genomic approaches allow for exploration of local adaptation in non-model systems, most landscape genomics studies of adaptive capacity are limited to exploratory identification of potentially important functional genes, often without a priori expectations as to the gene functions that may be most important for climate change responses. In this study, we integrated targeted sequencing of genes of known function and genotyping of single-nucleotide polymorphisms to examine spatial, environmental, and species-specific patterns of potential local adaptation in two co-occuring turtle species: the Blanding's turtle (Emydoidea blandingii) and the snapping turtle (Chelydra serpentina).

RESULTS

We documented divergent patterns of spatial clustering between neutral and putatively adaptive genetic variation in both species. Environmental associations varied among gene regions and between species, with stronger environmental associations detected for genes involved in stress response and for the more specialized Blanding's turtle. Land cover appeared to be more important than climate in shaping spatial variation in functional genes, indicating that human landscape alterations may affect adaptive capacity important for climate change responses.

CONCLUSIONS

Our study provides evidence that responses to climate change will be contingent on species-specific adaptive capacity and past history of exposure to human land cover change.

Are genomic updates of well-studied species worth the investment for conservation? A case study of the Critically Endangered Magdalena River turtle

As genomic-scale data sets become economically feasible for most organisms, a key question for conservation biology is whether the increased resolution offered by new genomic approaches justifies repeating earlier studies based on traditional markers, rather than investing those same time and monetary resources in less-known species. Genomic studies offer clear advantages when the objective is to identify adaptive loci that may be critical to conservation policy-makers. However, the answer is far less certain for the population and landscape studies based on neutral loci that dominate the conservation genetics research agenda. We used RADseq to revisit earlier molecular studies of the IUCN Critically Endangered Magdalena River turtle (Podocnemis lewyana), documenting the conservation insights gained by increasing the number of neutral markers by several orders of magnitude. Earlier research indicated that P. lewyana has the lowest genetic diversity known for any chelonian, and little or no population differentiation among independent rivers. In contrast, the RADseq data revealed discrete population structure with isolation-by-distance within river segments and identified precise population breaks clearly delineating management units. It also confirmed that the species does not have extremely low heterozygosity and that effective population sizes are probably sufficient to maintain long-term evolutionary potential. Contrary to earlier inferences from more limited population genetic markers, our genomic data suggest that management strategies should shift from active genetic rescue to more passive protection without extreme interventions. We conclude with a list of examples of conservation studies in other vertebrates indicating that for many systems a genomic update is worth the investment.

Herpetological phylogeographic analyses support a Miocene focal point of Himalayan uplift and biological diversification

The Himalaya are among the youngest and highest mountains in the world, but the exact timing of their uplift and origins of their biodiversity are still in debate. The Himalayan region is a relatively small area but with exceptional diversity and endemism. One common hypothesis to explain the rich montane diversity is uplift-driven diversification-that orogeny creates conditions favoring rapid in situ speciation of resident lineages. We test this hypothesis in the Himalayan region using amphibians and reptiles, two environmentally sensitive vertebrate groups. In addition, analysis of diversification of the herpetofauna provides an independent source of information to test competing geological hypotheses of Himalayan orogenesis. We conclude that the origins of the Himalayan herpetofauna date to the early Paleocene, but that diversification of most groups was concentrated in the Miocene. There was an increase in both rates and modes of diversification during the early to middle Miocene, together with regional interchange (dispersal) between the Himalaya and adjacent regions. Our analyses support a recently proposed stepwise geological model of Himalayan uplift beginning in the Paleocene, with a subsequent rapid increase of uplifting during the Miocene, finally giving rise to the intensification of the modern South Asian Monsoon.

Conservation genetics of the yellow-bellied toad (Bombina variegata): population structure, genetic diversity and landscape effects in an endangered amphibian

Revealing patterns of genetic diversity and barriers for gene flow are key points for successful conservation in endangered species. Methods based on molecular markers are also often used to delineate conservation units such as evolutionary significant units and management units. Here we combine phylo-geographic analyses (based on mtDNA) with population and landscape genetic analyses (based on microsatellites) for the endangered yellow-bellied toad Bombina variegata over a wide distribution range in Germany. Our analyses show that two genetic clusters are present in the study area, a northern and a southern/central one, but that these clusters are not deeply divergent. The genetic data suggest high fragmentation among toad occurrences and consequently low genetic diversity. Genetic diversity and genetic connectivity showed a negative relationship with road densities and urban areas surrounding toad occurrences, indicating that these landscape features act as barriers to gene flow. To preserve a maximum of genetic diversity, we recommend considering both genetic clusters as management units, and to increase gene flow among toad occurrences with the aim of restoring and protecting functional meta-populations within each of the clusters. Several isolated populations with especially low genetic diversity and signs of inbreeding need particular short-term conservation attention to avoid extinction. We also recommend to allow natural gene flow between both clusters but not to use individuals from one cluster for translocation or reintroduction into the other. Our results underscore the utility of molecular tools for species conservation, highlight outcomes of habitat fragmentation onto the genetic structure of an endangered amphibian and reveal particularly threatened populations in need for urgent conservation efforts.

Conservation Genomics of the Threatened Western Spadefoot, Spea hammondii, in Urbanized Southern California

Populations of the western spadefoot (Spea hammondii) in southern California occur in one of the most urbanized and fragmented landscapes on the planet and have lost up to 80% of their native habitat. Orange County is one of the last strongholds for this pond-breeding amphibian in the region, and ongoing restoration efforts targeting S. hammondii have involved habitat protection and the construction of artificial breeding ponds. These efforts have successfully increased breeding activity, but genetic characterization of the populations, including estimates of effective population size and admixture between the gene pools of constructed artificial and natural ponds, has never been undertaken. Using thousands of genome-wide single-nucleotide polymorphisms, we characterized the population structure, genetic diversity, and genetic connectivity of spadefoots in Orange County to guide ongoing and future management efforts. We identified at least 2, and possibly 3 major genetic clusters, with additional substructure within clusters indicating that individual ponds are often genetically distinct. Estimates of landscape resistance suggest that ponds on either side of the Los Angeles Basin were likely interconnected historically, but intense urban development has rendered them essentially isolated, and the resulting risk of interruption to natural metapopulation dynamics appears to be high. Resistance surfaces show that the existing artificial ponds were well-placed and connected to natural populations by low-resistance corridors. Toad samples from all ponds (natural and artificial) returned extremely low estimates of effective population size, possibly due to a bottleneck caused by a recent multi-year drought. Management efforts should focus on maintaining gene flow among natural and artificial ponds by both assisted migration and construction of new ponds to bolster the existing pond network in the region.

Perspectives on studying molecular adaptations of amphibians in the genomic era

Understanding the genetic mechanisms underlying particular adaptations/phenotypes of organisms is one of the core issues of evolutionary biology. The use of genomic data has greatly advanced our understandings on this issue, as well as other aspects of evolutionary biology, including molecular adaptation, speciation, and even conservation of endangered species. Despite the well-recognized advantages, usages of genomic data are still limited to non-mammal vertebrate groups, partly due to the difficulties in assembling large or highly heterozygous genomes. Although this is particularly the case for amphibians, nonetheless, several comparative and population genomic analyses have shed lights into the speciation and adaptation processes of amphibians in a complex landscape, giving a promising hope for a wider application of genomics in the previously believed challenging groups of organisms. At the same time, these pioneer studies also allow us to realize numerous challenges in studying the molecular adaptations and/or phenotypic evolutionary mechanisms of amphibians. In this review, we first summarize the recent progresses in the study of adaptive evolution of amphibians based on genomic data, and then we give perspectives regarding how to effectively identify key pathways underlying the evolution of complex traits in the genomic era, as well as directions for future research.

Population demographic history and adaptability of the vulnerable Lolokou Sucker Frog

Amphibians are experiencing worldwide declines due to increasing anthropogenetic disturbances. However, the genetic variability and hence adaptability are still unknown for most frogs. We integrated the mitochondrial (ND2 gene), nuclear (TYR gene) and major histocompatibility complex (MHC) loci, to clarify the demographic patterns and immune-gene diversity of the Lolokou Sucker Frog (Amolops loloensis). Demographic analysis of the ND2 and TYR genes suggested that the Lolokou Sucker Frog experienced a population expansion within the last 10,000 years. High MHC diversity was detected, which has likely resulted from positive selection, indicating the current diversity bodes well for the species' adaptive potential to pathogenic challenges. These findings broaden our knowledge on the population history and evolution adaptation of the reclusive torrent frog, and conservation implications are provided.

Phylogenomics and molecular species delimitation reveals great cryptic diversity of leaf-toed geckos (Phyllodactylidae: Phyllodactylus), ancient origins, and diversification in Mexico

We utilize the efficient GBS technique to obtain thousands of nuclear loci and SNPs to reconstruct the evolutionary history of Mexican leaf-toed geckos (Phyllodactylus). Through the incorporation of unprecedented sampling for this group of geckos, in combination with genomic data analysis, we generate mostly consistent phylogenetic hypotheses using two approaches: supermatrix and coalescent-based inference. All topologies depict three, mutually exclusive major clades. Clade I comprises P. bordai and all species closer to P. bordai than to any other Phyllodactylus. Clade II comprises P. nocticolus and all species closer to P. nocticolus than to any other Phyllodactylus. Clade III comprises P. tuberculosus and all species closer to P. tuberculosus than to any other Phyllodactylus. Analyses estimate the age for the most recent common ancestor of Phyllodactylus in the Eocene (~43 mya), and the ancestors of each major clade date to the Eocene-Oligocene transition (32-36 mya). This group includes one late-Eocene lineage (P. bordai), Oligocene lineages (P. paucituberculatus, P. delcampi), but also topological patterns that indicate a recent radiation occurred during the Pleistocene on islands in the Gulf of California. The wide spatial and temporal scale indicates a complex and unique biogeographic history for each major clade. The 33 species delimited by BPP and stepping-stone BFD*coalescent based genomic approaches reflect this history. This diversity delimited for Mexican leaf-toed geckos demonstrates a vast underestimation in the number of species based on morphological data alone.

Sequencing smart: De novo sequencing and assembly approaches for a non-model mammal

Background

Whilst much sequencing effort has focused on key mammalian model organisms such as mouse and human, little is known about the relationship between genome sequencing techniques for non-model mammals and genome assembly quality. This is especially relevant to non-model mammals, where the samples to be sequenced are often degraded and of low quality. A key aspect when planning a genome project is the choice of sequencing data to generate. This decision is driven by several factors, including the biological questions being asked, the quality of DNA available, and the availability of funds. Cutting-edge sequencing technologies now make it possible to achieve highly contiguous, chromosome-level genome assemblies, but rely on high-quality high molecular weight DNA. However, funding is often insufficient for many independent research groups to use these techniques. Here we use a range of different genomic technologies generated from a roadkill European polecat (Mustela putorius) to assess various assembly techniques on this low-quality sample. We evaluated different approaches for de novo assemblies and discuss their value in relation to biological analyses.

Results

Generally, assemblies containing more data types achieved better scores in our ranking system. However, when accounting for misassemblies, this was not always the case for Bionano and low-coverage 10x Genomics (for scaffolding only). We also find that the extra cost associated with combining multiple data types is not necessarily associated with better genome assemblies.

Conclusions

The high degree of variability between each de novo assembly method (assessed from the 7 key metrics) highlights the importance of carefully devising the sequencing strategy to be able to carry out the desired analysis. Adding more data to genome assemblies does not always result in better assemblies, so it is important to understand the nuances of genomic data integration explained here, in order to obtain cost-effective value for money when sequencing genomes.

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.

Patterns, Mechanisms and Genetics of Speciation in Reptiles and Amphibians

In this contribution, the aspects of reptile and amphibian speciation that emerged from research performed over the past decade are reviewed. First, this study assesses how patterns and processes of speciation depend on knowing the taxonomy of the group in question, and discuss how integrative taxonomy has contributed to speciation research in these groups. This study then reviews the research on different aspects of speciation in reptiles and amphibians, including biogeography and climatic niches, ecological speciation, the relationship between speciation rates and phenotypic traits, and genetics and genomics. Further, several case studies of speciation in reptiles and amphibians that exemplify many of these themes are discussed. These include studies of integrative taxonomy and biogeography in South American lizards, ecological speciation in European salamanders, speciation and phenotypic evolution in frogs and lizards. The final case study combines genomics and biogeography in tortoises. The field of amphibian and reptile speciation research has steadily moved forward from the assessment of geographic and ecological aspects, to incorporating other dimensions of speciation, such as genetic mechanisms and evolutionary forces. A higher degree of integration among all these dimensions emerges as a goal for future research.

Multi-tissue transcriptomes of caecilian amphibians highlight incomplete knowledge of vertebrate gene families

RNA sequencing (RNA-seq) has become one of the most powerful tools to unravel the genomic basis of biological adaptation and diversity. Although challenging, RNA-seq is particularly promising for research on non-model, secretive species that cannot be observed in nature easily and therefore remain comparatively understudied. Among such animals, the caecilians (order Gymnophiona) likely constitute the least known group of vertebrates, despite being an old and remarkably distinct lineage of amphibians. Here, we characterize multi-tissue transcriptomes for five species of caecilians that represent a broad level of diversity across the order. We identified vertebrate homologous elements of caecilian functional genes of varying tissue specificity that reveal a great number of unclassified gene families, especially for the skin. We annotated several protein domains for those unknown candidate gene families to investigate their function. We also conducted supertree analyses of a phylogenomic dataset of 1,955 candidate orthologous genes among five caecilian species and other major lineages of vertebrates, with the inferred tree being in agreement with current views of vertebrate evolution and systematics. Our study provides insights into the evolution of vertebrate protein-coding genes, and a basis for future research on the molecular elements underlying the particular biology and adaptations of caecilian amphibians.

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.

Ecological disturbance influences adaptive divergence despite high gene flow in golden perch (Macquaria ambigua): Implications for management and resilience to climate change

Populations that are adaptively divergent but maintain high gene flow may have greater resilience to environmental change as gene flow allows the spread of alleles that have already been tested elsewhere. In addition, populations naturally subjected to ecological disturbance may already hold resilience to future environmental change. Confirming this necessitates ecological genomic studies of high dispersal, generalist species. Here we perform one such study on golden perch (Macquaria ambigua) in the Murray-Darling Basin (MDB), Australia, using a genome-wide SNP data set. The MDB spans across arid to wet and temperate to subtropical environments, with low to high ecological disturbance in the form of low to high hydrological variability. We found high gene flow across the basin and three populations with low neutral differentiation. Genotype-environment association analyses detected adaptive divergence predominantly linked to an arid region with highly variable riverine flow, and candidate loci included functions related to fat storage, stress and molecular or tissue repair. The high connectivity of golden perch in the MDB will likely allow locally adaptive traits in its most arid and hydrologically variable environment to spread and be selected in localities that are predicted to become arid and hydrologically variable in future climates. High connectivity in golden perch is likely due to their generalist life history and efforts of fisheries management. Our study adds to growing evidence of adaptation in the face of gene flow and highlights the importance of considering ecological disturbance and adaptive divergence in biodiversity management.

Conservation genetics of the eastern yellow-bellied racer (Coluber constrictor flaviventris) and bullsnake (Pituophis catenifer sayi): River valleys are critical features for snakes at northern range limits

On the North American Great Plains, several snake species reach their northern range limit where they rely on sparsely distributed hibernacula located in major river valleys. Independent colonization histories for the river valleys and barriers to gene flow caused by the lack of suitable habitat between them may have produced genetically differentiated snake populations. To test this hypothesis, we used 10 microsatellite loci to examine the population structure of two species of conservation concern in Canada: the eastern yellow-bellied racer (Coluber constrictor flaviventris) and bullsnake (Pituophis catenifer sayi) in 3 major river valleys in southern Saskatchewan. Fixation indices (F_ST) showed that populations in river valleys were significantly differentiated for both species (racers, F_ST = 0.096, P = 0.001; bullsnakes F_ST = 0.045–0.157, P = 0.001). Bayesian assignment (STRUCTURE) and ordination (DAPC) strongly supported genetically differentiated groups in the geographically distinct river valleys. Finer-scale subdivision of populations within river valleys was not apparent based on our data, but is a topic that should be investigated further. Our findings highlight the importance of major river valleys for snakes at the northern extent of their ranges, and raise the possibility that populations in each river valley may warrant separate management strategies.

Genetic Structure and Eco-Geographical Differentiation of Wild Sheep Fescue (Festuca ovina L.) in Xinjiang, Northwest China

Glaciation and mountain orogeny have generated new ecologic opportunities for plants, favoring an increase in the speciation rate. Moreover, they also act as corridors or barriers for plant lineages and populations. High genetic diversity ensures that species are able to survive and adapt. Gene flow is one of the most important determinants of the genetic diversity and structure of out-crossed species, and it is easily affected by biotic and abiotic factors. The aim of this study was to characterize the genetic diversity and structure of an alpine species, Festuca ovina L., in Xinjiang, China. A total of 100 individuals from 10 populations were analyzed using six amplified fragment length polymorphism (AFLP) primer pairs. A total of 583 clear bands were generated, of which 392 were polymorphic; thus, the percentage of polymorphic bands (PPB) was 67.24%. The total and average genetic diversities were 0.2722 and 0.2006 (0.1686-0.2225), respectively. The unweighted group method with arithmetic mean (UPGMA) tree, principal coordinates analysis (PCoA) and Structure analyses revealed that these populations or individuals could be clustered into two groups. The analysis of molecular variance analysis (AMOVA) suggested that most of the genetic variance existed within a population, and the genetic differentiation (Fst) among populations was 20.71%. The Shannon differentiation coefficient (G'st) among populations was 0.2350. Limited gene flow (Nm = 0.9571) was detected across all sampling sites. The Fst and Nm presented at different levels under the genetic barriers due to fragmentation. The population genetic diversity was significant relative to environmental factors such as temperature, altitude and precipitation.

Multiple uprising invasions of Pelophylax water frogs, potentially inducing a new hybridogenetic complex

The genetic era has revolutionized our perception of biological invasions. Yet, it is usually too late to understand their genesis for efficient management. Here, we take the rare opportunity to reconstruct the scenario of an uprising invasion of the famous water frogs (Pelophylax) in southern France, through a fine-scale genetic survey. We identified three different taxa over less than 200 km²: the autochthonous P. perezi, along with the alien P. ridibundus and P. kurtmuelleri, which have suddenly become invasive. As a consequence, the latter hybridizes and may now form a novel hybridogenetic complex with P. perezi, which could actively promote its replacement. This exceptional situation makes a textbook application of genetics to early-detect, monitor and understand the onset of biological invasions before they pose a continental-wide threat. It further emphasizes the alarming rate of amphibian translocations, both at global and local scales, as well as the outstanding invasive potential of Pelophylax aliens.

Integrating Genomic Data Sets for Knowledge Discovery: An Informed Approach to Management of Captive Endangered Species

Many endangered captive populations exhibit reduced genetic diversity resulting in health issues that impact reproductive fitness and quality of life. Numerous cost effective genomic sequencing and genotyping technologies provide unparalleled opportunity for incorporating genomics knowledge in management of endangered species. Genomic data, such as sequence data, transcriptome data, and genotyping data, provide critical information about a captive population that, when leveraged correctly, can be utilized to maximize population genetic variation while simultaneously reducing unintended introduction or propagation of undesirable phenotypes. Current approaches aimed at managing endangered captive populations utilize species survival plans (SSPs) that rely upon mean kinship estimates to maximize genetic diversity while simultaneously avoiding artificial selection in the breeding program. However, as genomic resources increase for each endangered species, the potential knowledge available for management also increases. Unlike model organisms in which considerable scientific resources are used to experimentally validate genotype-phenotype relationships, endangered species typically lack the necessary sample sizes and economic resources required for such studies. Even so, in the absence of experimentally verified genetic discoveries, genomics data still provides value. In fact, bioinformatics and comparative genomics approaches offer mechanisms for translating these raw genomics data sets into integrated knowledge that enable an informed approach to endangered species management.

Sorting duplicated loci disentangles complexities of polyploid genomes masked by genotyping by sequencing

Many plants and animals of polyploid origin are currently enjoying a genomics explosion enabled by modern sequencing and genotyping technologies. However, routine filtering of duplicated loci in most studies using genotyping by sequencing introduces an unacceptable, but often overlooked, bias when detecting selection. Retained duplicates from ancient whole-genome duplications (WGDs) may be found throughout genomes, whereas retained duplicates from recent WGDs are concentrated at distal ends of some chromosome arms. Additionally, segmental duplicates can be found at distal ends or nearly anywhere in a genome. Evidence shows that these duplications facilitate adaptation through one of two pathways: neo-functionalization or increased gene expression. Filtering duplicates removes distal ends of some chromosomes, and distal ends are especially known to harbour adaptively important genes. Thus, filtering of duplicated loci impoverishes the interpretation of genomic data as signals from contiguous duplicated genes are ignored. We review existing strategies to genotype and map duplicated loci; we focus in detail on an overlooked strategy of using gynogenetic haploids (1N) as a part of new genotyping by sequencing studies. We provide guidelines on how to use this haploid strategy for studies on polyploid-origin vertebrates including how it can be used to screen duplicated loci in natural populations. We conclude by discussing areas of research that will benefit from better inclusion of polyploid loci; we particularly stress the sometimes overlooked fact that basing genomic studies on dense maps provides value added in the form of locating and annotating outlier loci or colocating outliers into islands of divergence.

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.