Genome-scale data reveal that endemic Poecilia populations from small sulphidic springs display no evidence of inbreeding

Populations with limited ranges can be highly vulnerable to changes in their environment and are, thus, of high conservation concern. Populations that experience human-induced range reductions are often highly inbred and lack genetic diversity, but it is unknown whether this is also the case for populations with naturally small ranges. The fishes Poecilia sulphuraria (listed as critically endangered) and Poecilia thermalis, which are endemic to small hydrogen sulphide-rich springs in southern Mexico, are examples of such populations with inherently small habitats. We used geometric morphometrics and population genetics to quantify phenotypic and genetic variation within and among two populations of P. sulphuraria and one population of P. thermalis. Principal component analyses revealed phenotypic and genetic differences among the populations. Evidence for inbreeding was low compared to populations that have undergone habitat reduction. The genetic data were also used to infer the demographic history of these populations to obtain estimates for effective population sizes and migration rates. Effective population sizes were large given the small habitats of these populations. Our results imply that these three endemic extremophile populations should each be considered separately for conservation purposes. Additionally, this study suggests that populations in naturally small habitats may have lower rates of inbreeding and higher genetic diversity than expected, and therefore may be better equipped to handle environmental perturbations than anticipated. We caution, however, that the inferred lack of inbreeding and the large effective population sizes could potentially be a result of colonization by genetically diverse ancestors.

The other white-nose syndrome transcriptome: Tolerant and susceptible hosts respond differently to the pathogen Pseudogymnoascus destructans

Mitigation of emerging infectious diseases that threaten global biodiversity requires an understanding of critical host and pathogen responses to infection. For multihost pathogens where pathogen virulence or host susceptibility is variable, host–pathogen interactions in tolerant species may identify potential avenues for adaptive evolution in recently exposed, susceptible hosts. For example, the fungus Pseudogymnoascus destructans causes white‐nose syndrome (WNS) in hibernating bats and is responsible for catastrophic declines in some species in North America, where it was recently introduced. Bats in Europe and Asia, where the pathogen is endemic, are only mildly affected. Different environmental conditions among Nearctic and Palearctic hibernacula have been proposed as an explanation for variable disease outcomes, but this hypothesis has not been experimentally tested. We report the first controlled, experimental investigation of response to P. destructans in a tolerant, European species of bat (the greater mouse‐eared bat, Myotis myotis). We compared body condition, disease outcomes and gene expression in control (sham‐exposed) and exposed M. myotis that hibernated under controlled environmental conditions following treatment. Tolerant M. myotis experienced extremely limited fungal growth and did not exhibit symptoms of WNS. However, we detected no differential expression of genes associated with immune response in exposed bats, indicating that immune response does not drive tolerance of P. destructans in late hibernation. Variable responses to P. destructans among bat species cannot be attributed solely to environmental or ecological factors. Instead, our results implicate coevolution with the pathogen, and highlight the dynamic nature of the “white‐nose syndrome transcriptome.” Interspecific variation in response to exposure by the host (and possibly pathogen) emphasizes the importance of context in studies of the bat‐WNS system, and robust characterization of genetic responses to exposure in various hosts and the pathogen should precede any attempts to use particular bat species as generalizable “model hosts.”

Conservation implications of limited genetic diversity and population structure in Tasmanian devils (Sarcophilus harrisii)

Tasmanian devils face a combination of threats to persistence, including Devil Facial Tumor Disease (DFTD), an epidemic transmissible cancer. We used RAD sequencing to investigate genome-wide patterns of genetic diversity and geographic population structure. Consistent with previous results, we found very low genetic diversity in the species as a whole, and we detected two broad genetic clusters occupying the northwestern portion of the range, and the central and eastern portions. However, these two groups overlap across a broad geographic area, and differentiation between them is modest (F_ST = 0.1081). Our results refine the geographic extent of the zone of mixed ancestry and substructure within it, potentially informing management of genetic variation that existed in pre-diseased populations of the species. DFTD has spread across both genetic clusters, but recent evidence points to a genomic response to selection imposed by DFTD. Any allelic variation for resistance to DFTD may be able to spread across the devil population under selection by DFTD, and/or be present as standing variation in both genetic regions.

Genome-wide nucleotide diversity of hatchery-reared Atlantic and Mediterranean strains of brown trout Salmo trutta compared to wild Mediterranean populations

A genome-wide assessment of diversity is provided for wild Mediterranean brown trout Salmo trutta populations from headwater tributaries of the Orb River and from Atlantic and Mediterranean hatchery-reared strains that have been used for stocking. Double-digest restriction-site-associated DNA sequencing (dd-RADseq) was performed and the efficiency of de novo and reference-mapping approaches to obtain individual genotypes was compared. Large numbers of single nucleotide polymorphism (SNP) markers with similar genome-wide distributions were discovered using both approaches (196 639 v. 121 016 SNPs, respectively), with c. 80% of the loci detected de novo being also found with reference mapping, using the Atlantic salmon Salmo salar genome as a reference. Lower mapping density but larger nucleotide diversity (π) was generally observed near extremities of linkage groups, consistent with regions of residual tetrasomic inheritance observed in salmonids. Genome-wide diversity estimates revealed reduced polymorphism in hatchery strains (π = 0·0040 and π = 0·0029 in Atlantic and Mediterranean strains, respectively) compared to wild populations (π = 0·0049), a pattern that was congruent with allelic richness estimated from microsatellite markers. Finally, pronounced heterozygote deficiency was found in hatchery strains (Atlantic F_IS = 0·18; Mediterranean F_IS = 0·42), indicating that stocking practices may affect the genetic diversity in wild populations. These new genomic resources will provide important tools to define better conservation strategies in S. trutta.

swinger: a user-friendly computer program to establish captive breeding groups that minimize relatedness without pedigree information

Captive breeding programmes are often a necessity for the continued persistence of a population or species. They typically have the goal of maintaining genetic diversity and minimizing inbreeding. However, most captive breeding programmes have been based on the assumption that the founding breeders are unrelated and outbred, even though in situ anthropogenic impacts often mean these founders may have high relatedness and substantial inbreeding. In addition, polygamous group-breeding species in captivity often have uncertain pedigrees, making it difficult to select the group composition for subsequent breeding. Molecular-based estimates of relatedness and inbreeding may instead be used to select breeding groups (≥two individuals) that minimize relatedness and filter out inbred individuals. swinger constructs breeding groups based on molecular estimates of relatedness and inbreeding. The number of possible combinations of breeding groups quickly becomes intractable by hand. swinger was designed to overcome this major issue in ex situ conservation biology. The user can specify parameters within swinger to reach breeding solutions that suit the mating system of the target species and available resources. We provide evidence of the efficiency of the software with an empirical example and using simulations. The only data required are a typical molecular marker data set, such as a microsatellite or SNP data set, from which estimates of inbreeding and pairwise relatedness may be obtained. Such molecular data sets are becoming easier to gather from non-model organisms with next-generation sequencing technology. swinger is an open-source software with a user-friendly interface and is available at http://www.molecularecology.flinders.edu.au/molecular-ecology-lab/software/swinger/swinger/ and https://github.com/Yuma248/Swinger.

Genomic Consequences of Population Decline in the Endangered Florida Scrub-Jay

Understanding the population genetic consequences of declining population size is important for conserving the many species worldwide facing severe decline [1]. Thorough empirical studies on the impacts of population reduction at a genome-wide scale in the wild are scarce because they demand huge field and laboratory investments [1, 2]. Previous studies have demonstrated the importance of gene flow in introducing genetic variation to small populations [3], but few have documented both genetic and fitness consequences of decreased immigration through time in a natural population [4-6]. Here we assess temporal variation in gene flow, inbreeding, and fitness using longitudinal genomic, demographic, and phenotypic data from a long-studied population of federally Threatened Florida scrub-jays (Aphelocoma coerulescens). We exhaustively sampled and genotyped the study population over two decades, providing one of the most detailed longitudinal investigations of genetics in a wild animal population to date. Immigrants were less heterozygous than residents but still introduced genetic variation into our study population. Owing to regional population declines, immigration into the study population declined from 1995-2013, resulting in increased levels of inbreeding and reduced fitness via inbreeding depression, even as the population remained demographically stable. Our results show that, contrary to conventional wisdom, small peripheral populations that already have undergone a genetic bottleneck may play a vital role in preserving genetic diversity of larger and seemingly stable populations. These findings underscore the importance of investing in the persistence of small populations and maintaining population connectivity in conservation of fragmented species.

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

Role of genomics and transcriptomics in selection of reintroduction source populations

The use and importance of reintroduction as a conservation tool to return a species to its historical range from which it has been extirpated will increase as climate change and human development accelerate habitat loss and population extinctions. Although the number of reintroduction attempts has increased rapidly over the past 2 decades, the success rate is generally low. As a result of population differences in fitness-related traits and divergent responses to environmental stresses, population performance upon reintroduction is highly variable, and it is generally agreed that selecting an appropriate source population is a critical component of a successful reintroduction. Conservation genomics is an emerging field that addresses long-standing challenges in conservation, and the potential for using novel molecular genetic approaches to inform and improve conservation efforts is high. Because the successful establishment and persistence of reintroduced populations is highly dependent on the functional genetic variation and environmental stress tolerance of the source population, we propose the application of conservation genomics and transcriptomics to guide reintroduction practices. Specifically, we propose using genome-wide functional loci to estimate genetic variation of source populations. This estimate can then be used to predict the potential for adaptation. We also propose using transcriptional profiling to measure the expression response of fitness-related genes to environmental stresses as a proxy for acclimation (tolerance) capacity. Appropriate application of conservation genomics and transcriptomics has the potential to dramatically enhance reintroduction success in a time of rapidly declining biodiversity and accelerating environmental change.

Wolf population genetics in Europe: a systematic review, meta-analysis and suggestions for conservation and management

The grey wolf (Canis lupus) is an iconic large carnivore that has increasingly been recognized as an apex predator with intrinsic value and a keystone species. However, wolves have also long represented a primary source of human-carnivore conflict, which has led to long-term persecution of wolves, resulting in a significant decrease in their numbers, genetic diversity and gene flow between populations. For more effective protection and management of wolf populations in Europe, robust scientific evidence is crucial. This review serves as an analytical summary of the main findings from wolf population genetic studies in Europe, covering major studies from the 'pre-genomic era' and the first insights of the 'genomics era'. We analyse, summarize and discuss findings derived from analyses of three compartments of the mammalian genome with different inheritance modes: maternal (mitochondrial DNA), paternal (Y chromosome) and biparental [autosomal microsatellites and single nucleotide polymorphisms (SNPs)]. To describe large-scale trends and patterns of genetic variation in European wolf populations, we conducted a meta-analysis based on the results of previous microsatellite studies and also included new data, covering all 19 European countries for which wolf genetic information is available: Norway, Sweden, Finland, Estonia, Latvia, Lithuania, Poland, Czech Republic, Slovakia, Germany, Belarus, Russia, Italy, Croatia, Bulgaria, Bosnia and Herzegovina, Greece, Spain and Portugal. We compared different indices of genetic diversity in wolf populations and found a significant spatial trend in heterozygosity across Europe from south-west (lowest genetic diversity) to north-east (highest). The range of spatial autocorrelation calculated on the basis of three characteristics of genetic diversity was 650-850 km, suggesting that the genetic diversity of a given wolf population can be influenced by populations up to 850 km away. As an important outcome of this synthesis, we discuss the most pressing issues threatening wolf populations in Europe, highlight important gaps in current knowledge, suggest solutions to overcome these limitations, and provide recommendations for science-based wolf conservation and management at regional and Europe-wide scales.

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction-site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small Ne (range = 2.1-89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome-wide divergence. Nonetheless, outlier tests identified 3.6-6.6% of loci as high FST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high FST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small Ne in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.

Whole-genome sequencing reveals small genomic regions of introgression in an introduced crater lake population of threespine stickleback

Invasive species pose a major threat to biological diversity. Although introduced populations often experience population bottlenecks, some invasive species are thought to be originated from hybridization between multiple populations or species, which can contribute to the maintenance of high genetic diversity. Recent advances in genome sequencing enable us to trace the evolutionary history of invasive species even at whole‐genome level and may help to identify the history of past hybridization that may be overlooked by traditional marker‐based analysis. Here, we conducted whole‐genome sequencing of eight threespine stickleback (Gasterosteus aculeatus) individuals, four from a recently introduced crater lake population and four of the putative source population. We found that both populations have several small genomic regions with high genetic diversity, which resulted from introgression from a closely related species (Gasterosteus nipponicus). The sizes of the regions were too small to be detected with traditional marker‐based analysis or even some reduced‐representation sequencing methods. Further amplicon sequencing revealed linkage disequilibrium around an introgression site, which suggests the possibility of selective sweep at the introgression site. Thus, interspecies introgression might predate introduction and increase genetic variation in the source population. Whole‐genome sequencing of even a small number of individuals can therefore provide higher resolution inference of history of introduced populations.

A nonlethal sampling method to obtain, generate and assemble whole blood transcriptomes from small, wild mammals

The acquisition of tissue samples from wild populations is a constant challenge in conservation biology, especially for endangered species and protected species where nonlethal sampling is the only option. Whole blood has been suggested as a nonlethal sample type that contains a high percentage of bodywide and genomewide transcripts and therefore can be used to assess the transcriptional status of an individual, and to infer a high percentage of the genome. However, only limited quantities of blood can be nonlethally sampled from small species and it is not known if enough genetic material is contained in only a few drops of blood, which represents the upper limit of sample collection for some small species. In this study, we developed a nonlethal sampling method, the laboratory protocols and a bioinformatic pipeline to sequence and assemble the whole blood transcriptome, using Illumina RNA-Seq, from wild greater mouse-eared bats (Myotis myotis). For optimal results, both ribosomal and globin RNAs must be removed before library construction. Treatment of DNase is recommended but not required enabling the use of smaller amounts of starting RNA. A large proportion of protein-coding genes (61%) in the genome were expressed in the blood transcriptome, comparable to brain (65%), kidney (63%) and liver (58%) transcriptomes, and up to 99% of the mitogenome (excluding D-loop) was recovered in the RNA-Seq data. In conclusion, this nonlethal blood sampling method provides an opportunity for a genomewide transcriptomic study of small, endangered or critically protected species, without sacrificing any individuals.

Morphological and genomic comparisons of Hawaiian and Japanese Black-footed Albatrosses (Phoebastria nigripes) using double digest RADseq: implications for conservation

Evaluating the genetic and demographic independence of populations of threatened species is important for determining appropriate conservation measures, but different technologies can yield different conclusions. Despite multiple studies, the taxonomic status and extent of gene flow between the main breeding populations of Black-footed Albatross (Phoebastria nigripes), a Near-Threatened philopatric seabird, are still controversial. Here, we employ double digest RADseq to quantify the extent of genomewide divergence and gene flow in this species. Our genomewide data set of 9760 loci containing 3455 single nucleotide polymorphisms yielded estimates of genetic diversity and gene flow that were generally robust across seven different filtering and sampling protocols and suggest a low level of genomic variation (θ per site = ∼0.00002-0.00028), with estimates of effective population size (N e = ∼500-15 881) falling far below current census size. Genetic differentiation was small but detectable between Japan and Hawaii (F ST ≈ 0.038-0.049), with no F ST outliers. Additionally, using museum specimens, we found that effect sizes of morphological differences by sex or population rarely exceeded 4%. These patterns suggest that the Hawaiian and Japanese populations exhibit small but significant differences and should be considered separate management units, although the evolutionary and adaptive consequences of this differentiation remain to be identified.

A field guide to whole-genome sequencing, assembly and annotation

Genome sequencing projects were long confined to biomedical model organisms and required the concerted effort of large consortia. Rapid progress in high-throughput sequencing technology and the simultaneous development of bioinformatic tools have democratized the field. It is now within reach for individual research groups in the eco-evolutionary and conservation community to generate de novo draft genome sequences for any organism of choice. Because of the cost and considerable effort involved in such an endeavour, the important first step is to thoroughly consider whether a genome sequence is necessary for addressing the biological question at hand. Once this decision is taken, a genome project requires careful planning with respect to the organism involved and the intended quality of the genome draft. Here, we briefly review the state of the art within this field and provide a step-by-step introduction to the workflow involved in genome sequencing, assembly and annotation with particular reference to large and complex genomes. This tutorial is targeted at scientists with a background in conservation genetics, but more generally, provides useful practical guidance for researchers engaging in whole-genome sequencing projects.

Nuclear and mitochondrial genetic structure in the Eurasian beaver (Castor fiber) - implications for future reintroductions

Many reintroduction projects for conservation fail, and there are a large number of factors that may contribute to failure. Genetic analysis can be used to help stack the odds of a reintroduction in favour of success, by conducting assessment of source populations to evaluate the possibility of inbreeding and outbreeding depression and by conducting postrelease monitoring. In this study, we use a panel of 306 SNP (single nucleotide polymorphism) markers and 487-489 base pairs of mitochondrial DNA control region sequence data to examine 321 individuals from possible source populations of the Eurasian beaver for a reintroduction to Scotland. We use this information to reassess the phylogenetic history of the Eurasian beavers, to examine the genetic legacy of past reintroductions on the Eurasian landmass and to assess the future power of the genetic markers to conduct ongoing monitoring via parentage analysis and individual identification. We demonstrate the capacity of medium density genetic data (hundreds of SNPs) to provide information suitable for applied conservation and discuss the difficulty of balancing the need for high genetic diversity against phylogenetic best fit when choosing source population(s) for reintroduction.