Targeted capture and resequencing of 1040 genes reveal environmentally driven functional variation in grey wolves

Conservation genomics of wolves: The global impact of RK Wayne's research

RK Wayne has arguably been the most influential geneticist of canids, famously promoting the conservation of wolves in his homeland, the United States. His influence has been felt in other countries and regions outside the contiguous United States, where he inspired others, also including former graduate students and research fellows of his, to use modern molecular techniques to examine the evolutionary biology of canids to inform the conservation and management of wolves. In this review, we focus on the implications of Wayne's work on wolves outside the United States. He envisioned a clear future for wolf conservation research, involving the study of wolves' ecological and genetic diversity, and the description of ecotypes requiring conservation. He also documented widespread hybridization among canids and introgression of DNA from domestic dogs to wolves, a process that started dozens of thousands of years ago. His work therefore calls for innovative studies, such as examining the potential fitness benefits of introgression. Inspired by his results, for example, on the purging of deleterious alleles in small populations, wolf researchers should use novel molecular tools to challenge other conservation genetics paradigms. Overall, RK Wayne's work constitutes a call for answers, which as scientists or citizens concerned with conservation matters, we are obliged to address, as we contribute to monitoring and maintaining biodiversity during our period of dramatic transformations of the biosphere.

Human subsistence and signatures of selection on chemosensory genes

Chemosensation (olfaction, taste) is essential for detecting and assessing foods, such that dietary shifts elicit evolutionary changes in vertebrate chemosensory genes. The transition from hunting and gathering to agriculture dramatically altered how humans acquire food. Recent genetic and linguistic studies suggest agriculture may have precipitated olfactory degeneration. Here, we explore the effects of subsistence behaviors on olfactory (OR) and taste (TASR) receptor genes among rainforest foragers and neighboring agriculturalists in Africa and Southeast Asia. We analyze 378 functional OR and 26 functional TASR genes in 133 individuals across populations in Uganda (Twa, Sua, BaKiga) and the Philippines (Agta, Mamanwa, Manobo) with differing subsistence histories. We find no evidence of relaxed selection on chemosensory genes in agricultural populations. However, we identify subsistence-related signatures of local adaptation on chemosensory genes within each geographic region. Our results highlight the importance of culture, subsistence economy, and drift in human chemosensory perception.

Possible origins and implications of atypical morphologies and domestication-like traits in wild golden jackals (Canis aureus)

Deciphering the origins of phenotypic variations in natural animal populations is a challenging topic for evolutionary and conservation biologists. Atypical morphologies in mammals are usually attributed to interspecific hybridisation or de-novo mutations. Here we report the case of four golden jackals (Canis aureus), that were observed during a camera-trapping wildlife survey in Northern Israel, displaying anomalous morphological traits, such as white patches, an upturned tail, and long thick fur which resemble features of domesticated mammals. Another individual was culled under permit and was genetically and morphologically examined. Paternal and nuclear genetic profiles, as well as geometric morphometric data, identified this individual as a golden jackal rather than a recent dog/wolf-jackal hybrid. Its maternal haplotype suggested past introgression of African wolf (Canis lupaster) mitochondrial DNA, as previously documented in other jackals from Israel. When viewed in the context of the jackal as an overabundant species in Israel, the rural nature of the surveyed area, the abundance of anthropogenic waste, and molecular and morphological findings, the possibility of an individual presenting incipient stages of domestication should also be considered.

Genetic diversity loss in the Anthropocene

Anthropogenic habitat loss and climate change are reducing species' geographic ranges, increasing extinction risk and losses of species' genetic diversity. Although preserving genetic diversity is key to maintaining species' adaptability, we lack predictive tools and global estimates of genetic diversity loss across ecosystems. We introduce a mathematical framework that bridges biodiversity theory and population genetics to understand the loss of naturally occurring DNA mutations with decreasing habitat. By analyzing genomic variation of 10,095 georeferenced individuals from 20 plant and animal species, we show that genome-wide diversity follows a mutations-area relationship power law with geographic area, which can predict genetic diversity loss from local population extinctions. We estimate that more than 10% of genetic diversity may already be lost for many threatened and nonthreatened species, surpassing the United Nations' post-2020 targets for genetic preservation.

Genomic Assessment of Cancer Susceptibility in the Threatened Catalina Island Fox (Urocyon littoralis catalinae)

Small effective population sizes raise the probability of extinction by increasing the frequency of potentially deleterious alleles and reducing fitness. However, the extent to which cancers play a role in the fitness reduction of genetically depauperate wildlife populations is unknown. Santa Catalina island foxes (Urocyon littoralis catalinae) sampled in 2007-2008 have a high prevalence of ceruminous gland tumors, which was not detected in the population prior to a recent bottleneck caused by a canine distemper epidemic. The disease appears to be associated with inflammation from chronic ear mite (Otodectes) infections and secondary elevated levels of Staphyloccus pseudointermedius bacterial infections. However, no other environmental factors to date have been found to be associated with elevated cancer risk in this population. Here, we used whole genome sequencing of the case and control individuals from two islands to identify candidate loci associated with cancer based on genetic divergence, nucleotide diversity, allele frequency spectrum, and runs of homozygosity. We identified several candidate loci based on genomic signatures and putative gene functions, suggesting that cancer susceptibility in this population may be polygenic. Due to the efforts of a recovery program and weak fitness effects of late-onset disease, the population size has increased, which may allow selection to be more effective in removing these presumably slightly deleterious alleles. Long-term monitoring of the disease alleles, as well as overall genetic diversity, will provide crucial information for the long-term persistence of this threatened population.

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype-environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.

K Locus Effects in Gray Wolves: Experimental Assessment of TLR3 Signaling and the Gene Expression Response to Canine Distemper Virus

In North American gray wolves, black coat color is dominantly inherited via a 3 base pair coding deletion in the canine beta defensin 3 (CBD103) gene. This 3 base pair deletion, called the KB allele, was introduced through hybridization with dogs and subsequently underwent a selective sweep that increased its frequency in wild wolves. Despite apparent positive selection, KBB wolves have lower fitness than wolves with the KyB genotype, even though the 2 genotypes show no observable differences in black coat color. Thus, the KB allele is thought to have pleiotropic effects on as-yet unknown phenotypes. Given the role of skin-expressed CBD103 in innate immunity, we hypothesized that the KB allele influences the keratinocyte gene expression response to TLR3 pathway stimulation and/or infection by canine distemper virus (CDV). To test this hypothesis, we developed a panel of primary epidermal keratinocyte cell cultures from 24 wild North American gray wolves of both Kyy and KyB genotypes. In addition, we generated an immortalized Kyy line and used CRISPR/Cas9 editing to produce a KyB line on the same genetic background. We assessed the transcriptome-wide responses of wolf keratinocytes to the TLR3 agonist polyinosinic:polycytidylic acid (polyI:C), and to live CDV. K locus genotype did not predict the transcriptional response to either challenge, suggesting that variation in the gene expression response does not explain pleiotropic effects of the KB allele on fitness. This study supports the feasibility of using cell culture methods to investigate the phenotypic effects of naturally occurring genetic variation in wild mammals.

An evolutionary conserved olfactory receptor for foodborne and semiochemical alkylpyrazines

Molecular recognition is a fundamental principle in biological systems. The olfactory detection of both food and predators via ecological relevant odorant cues are abilities of eminent evolutionary significance for many species. Pyrazines are such volatile cues, some of which act as both human-centered key food odorants (KFOs) and semiochemicals. A pyrazine-selective odorant receptor has been elusive. Here we screened 2,3,5-trimethylpyrazine, a KFO and semiochemical, and 2,5-dihydro-2,4,5-trimethylthiazoline, an innate fear-associated non-KFO, against 616 human odorant receptor variants, in a cell-based luminescence assay. OR5K1 emerged as sole responding receptor. Tested against a comprehensive collection of 178 KFOs, we newly identified 18 pyrazines and (2R/2S)-4-methoxy-2,5-dimethylfuran-3(2H)-one as agonists. Notably, OR5K1 orthologs in mouse and domesticated species displayed a human-like, potency-ranked activation pattern of pyrazines, suggesting a domestication-led co-evolution of OR5K1 and its orthologs. In summary, OR5K1 is a specialized olfactory receptor across mammals for the detection of pyrazine-based key food odors and semiochemicals.

Population Genomics Training for the Next Generation of Conservation Geneticists: ConGen 2018 Workshop

The increasing availability and complexity of next-generation sequencing (NGS) data sets make ongoing training an essential component of conservation and population genetics research. A workshop entitled “ConGen 2018” was recently held to train researchers in conceptual and practical aspects of NGS data production and analysis for conservation and ecological applications. Sixteen instructors provided helpful lectures, discussions, and hands-on exercises regarding how to plan, produce, and analyze data for many important research questions. Lecture topics ranged from understanding probabilistic (e.g., Bayesian) genotype calling to the detection of local adaptation signatures from genomic, transcriptomic, and epigenomic data. We report on progress in addressing central questions of conservation genomics, advances in NGS data analysis, the potential for genomic tools to assess adaptive capacity, and strategies for training the next generation of conservation genomicists.

Big Data in Conservation Genomics: Boosting Skills, Hedging Bets, and Staying Current in the Field

A current challenge in the fields of evolutionary, ecological, and conservation genomics is balancing production of large-scale datasets with additional training often required to handle such datasets. Thus, there is an increasing need for conservation geneticists to continually learn and train to stay up-to-date through avenues such as symposia, meetings, and workshops. The ConGen meeting is a near-annual workshop that strives to guide participants in understanding population genetics principles, study design, data processing, analysis, interpretation, and applications to real-world conservation issues. Each year of ConGen gathers a diverse set of instructors, students, and resulting lectures, hands-on sessions, and discussions. Here, we summarize key lessons learned from the 2019 meeting and more recent updates to the field with a focus on big data in conservation genomics. First, we highlight classical and contemporary issues in study design that are especially relevant to working with big datasets, including the intricacies of data filtering. We next emphasize the importance of building analytical skills and simulating data, and how these skills have applications within and outside of conservation genetics careers. We also highlight recent technological advances and novel applications to conservation of wild populations. Finally, we provide data and recommendations to support ongoing efforts by ConGen organizers and instructors-and beyond-to increase participation of underrepresented minorities in conservation and eco-evolutionary sciences. The future success of conservation genetics requires both continual training in handling big data and a diverse group of people and approaches to tackle key issues, including the global biodiversity-loss crisis.

The role of a mechanistic host in maintaining arctic rabies variant distributions: Assessment of functional genetic diversity in Alaskan red fox (Vulpes vulpes)

Populations are exposed to different types and strains of pathogens across heterogeneous landscapes, where local interactions between host and pathogen may present reciprocal selective forces leading to correlated patterns of spatial genetic structure. Understanding these coevolutionary patterns provides insight into mechanisms of disease spread and maintenance. Arctic rabies (AR) is a lethal disease with viral variants that occupy distinct geographic distributions across North America and Europe. Red fox (Vulpes vulpes) are a highly susceptible AR host, whose range overlaps both geographically distinct AR strains and regions where AR is absent. It is unclear if genetic structure exists among red fox populations relative to the presence/absence of AR or the spatial distribution of AR variants. Acquiring these data may enhance our understanding of the role of red fox in AR maintenance/spread and inform disease control strategies. Using a genotyping-by-sequencing assay targeting 116 genomic regions of immunogenetic relevance, we screened for sequence variation among red fox populations from Alaska and an outgroup from Ontario, including areas with different AR variants, and regions where the disease was absent. Presumed neutral SNP data from the assay found negligible levels of neutral genetic structure among Alaskan populations. The immunogenetically-associated data identified 30 outlier SNPs supporting weak to moderate genetic structure between regions with and without AR in Alaska. The outliers included SNPs with the potential to cause missense mutations within several toll-like receptor genes that have been associated with AR outcome. In contrast, there was a lack of genetic structure between regions with different AR variants. Combined, we interpret these data to suggest red fox populations respond differently to the presence of AR, but not AR variants. This research increases our understanding of AR dynamics in the Arctic, where host/disease patterns are undergoing flux in a rapidly changing Arctic landscape, including the continued northward expansion of red fox into regions previously predominated by the arctic fox (Vulpes lagopus).

Environmental correlates of genetic variation in the invasive European starling in North America

Populations of invasive species that colonize and spread in novel environments may differentiate both through demographic processes and local selection. European starlings (Sturnus vulgaris) were introduced to New York in 1890 and subsequently spread throughout North America, becoming one of the most widespread and numerous bird species on the continent. Genome-wide comparisons across starling individuals and populations can identify demographic and/or selective factors that facilitated this rapid and successful expansion. We investigated patterns of genomic diversity and differentiation using reduced-representation genome sequencing of 17 winter-season sampling sites. Consistent with this species' high dispersal rate and rapid expansion history, we found low geographical differentiation and few F_ST outliers even at a continental scale. Despite starting from a founding population of ~180 individuals, North American starlings show only a moderate genetic bottleneck, and models suggest a dramatic increase in effective population size since introduction. In genotype-environment associations we found that ~200 single-nucleotide polymorphisms are correlated with temperature and/or precipitation against a background of negligible genome- and range-wide divergence. Given this evidence, we suggest that local adaptation in North American starlings may have evolved rapidly even in this wide-ranging and evolutionarily young system. This survey of genomic signatures of expansion in North American starlings is the most comprehensive to date and complements ongoing studies of world-wide local adaptation in these highly dispersive and invasive birds.

A Genome-Wide Perspective on the Persistence of Red Wolf Ancestry in Southeastern Canids

The red wolf (Canis rufus), a legally recognized and critically endangered wolf, is known to interbreed with coyotes (Canis latrans). Declared extirpated in the wild in 1980, red wolves were reintroduced to northeastern North Carolina nearly a decade later. Interbreeding with coyotes was thought to be restricted to a narrow geographic region adjacent to the reintroduced population and largely believed to threaten red wolf recovery. However, red wolf ancestry was recently discovered in canids along the American Gulf Coast, igniting a broader survey of ancestry in southeastern canid populations. Here, we examine geographic and temporal patterns of genome-wide red wolf ancestry in 260 canids across the southeastern United States at over 164 000 SNP loci. We found that red wolf ancestry was most prevalent in canids sampled from Texas in the mid-1970s, although non-trivial amounts of red wolf ancestry persist in this region today. Further, red wolf ancestry was also observed in a subset of coyotes inhabiting North Carolina, despite management efforts to limit the occurrence of hybridization events. Lastly, we found no evidence of substantial red wolf ancestry in southeastern canids outside of these 2 admixture zones. Overall, this study provides a genome-wide survey of red wolf ancestry in canids across the southeastern United States, which may ultimately inform future red wolf restoration efforts.

Evidence for adaptive introgression of exons across a hybrid swarm in deer

Background

Secondary contact between closely related lineages can result in a variety of outcomes, including hybridization, depending upon the strength of reproductive barriers. By examining the extent to which different parts of the genome introgress, it is possible to infer the strength of selection and gain insight into the evolutionary trajectory of lineages. Following secondary contact approximately 8000 years ago in the Pacific Northwest, mule deer (Odocoileus hemionus hemionus) and black-tailed deer (O. h. columbianus) formed a hybrid swarm along the Cascade mountain range despite substantial differences in body size (up to two times) and habitat preference. In this study, we examined genetic population structure, extent of introgression, and selection pressures in freely interbreeding populations of mule deer and black-tailed deer using mitochondrial DNA sequences, 9 microsatellite loci, and 95 SNPs from protein-coding genes.

Results

We observed bi-directional hybridization and classified approximately one third of the 172 individuals as hybrids, almost all of which were beyond the F1 generation. High genetic differentiation between black-tailed deer and mule deer at protein-coding genes suggests that there is positive divergent selection, though selection on these loci is relatively weak. Contrary to predictions, there was not greater selection on protein-coding genes thought to be associated with immune function and mate choice. Geographic cline analyses were consistent across genetic markers, suggesting long-term stability (over hundreds of generations), and indicated that the center of the hybrid swarm is 20-30 km to the east of the Cascades ridgeline, where there is a steep ecological transition from wet, forested habitat to dry, scrub habitat.

Conclusions

Our data are consistent with a genetic boundary between mule deer and black-tailed deer that is porous but maintained by many loci under weak selection having a substantial cumulative effect. The absence of clear reproductive barriers and the consistent centering of geographic clines at a sharp ecotone suggests that ecology is a driver of hybrid swarm dynamics. Adaptive introgression in this study (and others) promotes gene flow and provides valuable insight into selection strength on specific genes and the evolutionary trajectory of hybridizing taxa.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1497-x) contains supplementary material, which is available to authorized users.

The evolutionary history of grey wolf Y chromosomes

Analyses of Y chromosome haplotypes uniquely provide a paternal picture of evolutionary histories and offer a very useful contrast to studies based on maternally inherited mitochondrial DNA (mtDNA). Here we used a bioinformatic approach based on comparison of male and female sequence coverage to identify 4.7 Mb from the grey wolf (Canis lupis) Y chromosome, probably representing most of the male-specific, nonampliconic sequence from the euchromatic part of the chromosome. We characterized this sequence and then identified ≈1,500 Y-linked single nucleotide polymorphisms in a sample of 145 resequenced male wolves, including 75 Finnish wolf genomes newly sequenced in this study, and in 24 dogs and eight other canids. We found 53 Y chromosome haplotypes, of which 26 were seen in grey wolves, that clustered in four major haplogroups. All four haplogroups were represented in samples of Finnish wolves, showing that haplogroup lineages were not partitioned on a continental scale. However, regional population structure was indicated because individual haplotypes were never shared between geographically distant areas, and genetically similar haplotypes were only found within the same geographical region. The deepest split between grey wolf haplogroups was estimated to have occurred 125,000 years ago, which is considerably older than recent estimates of the time of divergence of wolf populations. The distribution of dogs in a phylogenetic tree of Y chromosome haplotypes supports multiple domestication events, or wolf paternal introgression, starting 29,000 years ago. We also addressed the disputed origin of a recently founded population of Scandinavian wolves and observed that founding as well as most recent immigrant haplotypes were present in the neighbouring Finnish population, but not in sequenced wolves from elsewhere in the world, or in dogs.

Natural Selection and Origin of a Melanistic Allele in North American Gray Wolves

Pigmentation is often used to understand how natural selection affects genetic variation in wild populations since it can have a simple genetic basis, and can affect a variety of fitness-related traits (e.g., camouflage, thermoregulation, and sexual display). In gray wolves, the K locus, a β-defensin gene, causes black coat color via a dominantly inherited KB allele. The allele is derived from dog-wolf hybridization and is at high frequency in North American wolf populations. We designed a DNA capture array to probe the geographic origin, age, and number of introgression events of the KB allele in a panel of 331 wolves and 20 dogs. We found low diversity in KB, but not ancestral ky, wolf haplotypes consistent with a selective sweep of the black haplotype across North America. Further, North American wolf KB haplotypes are monophyletic, suggesting that a single adaptive introgression from dogs to wolves most likely occurred in the Northwest Territories or Yukon. We use a new analytical approach to date the origin of the KB allele in Yukon wolves to between 1,598 and 7,248 years ago, suggesting that introgression with early Native American dogs was the source. Using population genetic simulations, we show that the K locus is undergoing natural selection in four wolf populations. We find evidence for balancing selection, specifically in Yellowstone wolves, which could be a result of selection for enhanced immunity in response to distemper. With these data, we demonstrate how the spread of an adaptive variant may have occurred across a species' geographic range.

Digging for the spiny rat and hutia phylogeny using a gene capture approach, with the description of a new mammal subfamily

Next generation sequencing (NGS) and genomic database mining allow biologists to gather and select large molecular datasets well suited to address phylogenomics and molecular evolution questions. Here we applied this approach to a mammal family, the Echimyidae, for which generic relationships have been difficult to recover and often referred to as a star phylogeny. These South-American spiny rats represent a family of caviomorph rodents exhibiting a striking diversity of species and life history traits. Using a NGS exon capture protocol, we isolated and sequenced ca. 500 nuclear DNA exons for 35 species belonging to all major echimyid and capromyid clades. Exons were carefully selected to encompass as much diversity as possible in terms of rate of evolution, heterogeneity in the distribution of site-variation and nucleotide composition. Supermatrix inferences and coalescence-based approaches were subsequently applied to infer this family's phylogeny. The inferred topologies were the same for both approaches, and support was maximal for each node, entirely resolving the ambiguous relationships of previous analyses. Fast-evolving nuclear exons tended to yield more reliable phylogenies, as slower-evolving sequences were not informative enough to disentangle the short branches of the Echimyidae radiation. Based on this resolved phylogeny and on molecular and morphological evidence, we confirm the rank of the Caribbean hutias - formerly placed in the Capromyidae family - as Capromyinae, a clade nested within Echimyidae. We also name and define Carterodontinae, a new subfamily of Echimyidae, comprising the extant monotypic genus Carterodon from Brazil, which is the closest living relative of West Indies Capromyinae.

Natural re-colonization and admixture of wolves (Canis lupus) in the US Pacific Northwest: challenges for the protection and management of rare and endangered taxa

Admixture resulting from natural dispersal processes can potentially generate novel phenotypic variation that may facilitate persistence in changing environments or result in the loss of population-specific adaptations. Yet, under the US Endangered Species Act, policy is limited for management of individuals whose ancestry includes a protected taxon; therefore, they are generally not protected under the Act. This issue is exemplified by the recently re-established grey wolves of the Pacific Northwest states of Washington and Oregon, USA. This population was likely founded by two phenotypically and genetically distinct wolf ecotypes: Northern Rocky Mountain (NRM) forest and coastal rainforest. The latter is considered potentially threatened in southeast Alaska and thus the source of migrants may affect plans for their protection. To assess the genetic source of the re-established population, we sequenced a ~ 300 bp portion of the mitochondrial control region and ~ 5 Mbp of the nuclear genome. Genetic analysis revealed that the Washington wolves share ancestry with both wolf ecotypes, whereas the Oregon population shares ancestry with NRM forest wolves only. Using ecological niche modelling, we found that the Pacific Northwest states contain environments suitable for each ecotype, with wolf packs established in both environmental types. Continued migration from coastal rainforest and NRM forest source populations may increase the genetic diversity of the Pacific Northwest population. However, this admixed population challenges traditional management regimes given that admixture occurs between an adaptively distinct ecotype and a more abundant reintroduced interior form. Our results emphasize the need for a more precise US policy to address the general problem of admixture in the management of endangered species, subspecies, and distinct population segments.

Population genomics of grey wolves and wolf-like canids in North America

North America is currently home to a number of grey wolf (Canis lupus) and wolf-like canid populations, including the coyote (Canis latrans) and the taxonomically controversial red, Eastern timber and Great Lakes wolves. We explored their population structure and regional gene flow using a dataset of 40 full genome sequences that represent the extant diversity of North American wolves and wolf-like canid populations. This included 15 new genomes (13 North American grey wolves, 1 red wolf and 1 Eastern timber/Great Lakes wolf), ranging from 0.4 to 15x coverage. In addition to providing full genome support for the previously proposed coyote-wolf admixture origin for the taxonomically controversial red, Eastern timber and Great Lakes wolves, the discriminatory power offered by our dataset suggests all North American grey wolves, including the Mexican form, are monophyletic, and thus share a common ancestor to the exclusion of all other wolves. Furthermore, we identify three distinct populations in the high arctic, one being a previously unidentified "Polar wolf" population endemic to Ellesmere Island and Greenland. Genetic diversity analyses reveal particularly high inbreeding and low heterozygosity in these Polar wolves, consistent with long-term isolation from the other North American wolves.

The International Mouse Phenotyping Consortium (IMPC): a functional catalogue of the mammalian genome that informs conservation

The International Mouse Phenotyping Consortium (IMPC) is building a catalogue of mammalian gene function by producing and phenotyping a knockout mouse line for every protein-coding gene. To date, the IMPC has generated and characterised 5186 mutant lines. One-third of the lines have been found to be non-viable and over 300 new mouse models of human disease have been identified thus far. While current bioinformatics efforts are focused on translating results to better understand human disease processes, IMPC data also aids understanding genetic function and processes in other species. Here we show, using gorilla genomic data, how genes essential to development in mice can be used to help assess the potentially deleterious impact of gene variants in other species. This type of analyses could be used to select optimal breeders in endangered species to maintain or increase fitness and avoid variants associated to impaired-health phenotypes or loss-of-function mutations in genes of critical importance. We also show, using selected examples from various mammal species, how IMPC data can aid in the identification of candidate genes for studying a condition of interest, deliver information about the mechanisms involved, or support predictions for the function of genes that may play a role in adaptation. With genotyping costs decreasing and the continued improvements of bioinformatics tools, the analyses we demonstrate can be routinely applied.

New insights into the phylogenetics and population structure of the prairie falcon (Falco mexicanus)

Background

Management requires a robust understanding of between- and within-species genetic variability, however such data are still lacking in many species. For example, although multiple population genetics studies of the peregrine falcon (Falco peregrinus) have been conducted, no similar studies have been done of the closely-related prairie falcon (F. mexicanus) and it is unclear how much genetic variation and population structure exists across the species’ range. Furthermore, the phylogenetic relationship of F. mexicanus relative to other falcon species is contested. We utilized a genomics approach (i.e., genome sequencing and assembly followed by single nucleotide polymorphism genotyping) to rapidly address these gaps in knowledge.

Results

We sequenced the genome of a single female prairie falcon and generated a 1.17 Gb (gigabases) draft genome assembly. We generated maximum likelihood phylogenetic trees using complete mitochondrial genomes as well as nuclear protein-coding genes. This process provided evidence that F. mexicanus is an outgroup to the clade that includes the peregrine falcon and members of the subgenus Hierofalco. We annotated > 16,000 genes and almost 600,000 high-quality single nucleotide polymorphisms (SNPs) in the nuclear genome, providing the raw material for a SNP assay design featuring > 140 gene-associated markers and a molecular-sexing marker. We subsequently genotyped ~ 100 individuals from California (including the San Francisco East Bay Area, Pinnacles National Park and the Mojave Desert) and Idaho (Snake River Birds of Prey National Conservation Area). We tested for population structure and found evidence that individuals sampled in California and Idaho represent a single panmictic population.

Conclusions

Our study illustrates how genomic resources can rapidly shed light on genetic variability in understudied species and resolve phylogenetic relationships. Furthermore, we found evidence of a single, randomly mating population of prairie falcons across our sampling locations. Prairie falcons are highly mobile and relatively rare long-distance dispersal events may promote gene flow throughout the range. As such, California’s prairie falcons might be managed as a single population, indicating that management actions undertaken to benefit the species at the local level have the potential to influence the species as a whole.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4615-z) contains supplementary material, which is available to authorized users.

Identifying genome-wide immune gene variation underlying infectious disease in wildlife populations - a next generation sequencing approach in the gopher tortoise

Background

Infectious disease is the single greatest threat to taxa such as amphibians (chytrid fungus), bats (white nose syndrome), Tasmanian devils (devil facial tumor disease), and black-footed ferrets (canine distemper virus, plague). Although understanding the genetic basis to disease susceptibility is important for the long-term persistence of these groups, most research has been limited to major-histocompatibility and Toll-like receptor genes. To better understand the genetic basis of infectious disease susceptibility in a species of conservation concern, we sequenced all known/predicted immune response genes (i.e., the immunomes) in 16 Florida gopher tortoises, Gopherus polyphemus. All tortoises produced antibodies against Mycoplasma agassizii (an etiologic agent of infectious upper respiratory tract disease; URTD) and, at the time of sampling, either had (n = 10) or lacked (n = 6) clinical signs.

Results

We found several variants associated with URTD clinical status in complement and lectin genes, which may play a role in Mycoplasma immunity. Thirty-five genes deviated from neutrality according to Tajima’s D. These genes were enriched in functions relating to macromolecule and protein modifications, which are vital to immune system functioning.

Conclusions

These results are suggestive of genetic differences that might contribute to disease severity, a finding that is consistent with other mycoplasmal diseases. This has implications for management because tortoises across their range may possess genetic variation associated with a more severe response to URTD. More generally: 1) this approach demonstrates that a broader consideration of immune genes is better able to identify important variants, and; 2) this data pipeline can be adopted to identify alleles associated with disease susceptibility or resistance in other taxa, and therefore provide information on a population’s risk of succumbing to disease, inform translocations to increase genetic variation for disease resistance, and help to identify potential treatments.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4452-0) contains supplementary material, which is available to authorized users.

Development of a genotype-by-sequencing immunogenetic assay as exemplified by screening for variation in red fox with and without endemic rabies exposure

Pathogens are recognized as major drivers of local adaptation in wildlife systems. By determining which gene variants are favored in local interactions among populations with and without disease, spatially explicit adaptive responses to pathogens can be elucidated. Much of our current understanding of host responses to disease comes from a small number of genes associated with an immune response. High‐throughput sequencing (HTS) technologies, such as genotype‐by‐sequencing (GBS), facilitate expanded explorations of genomic variation among populations. Hybridization‐based GBS techniques can be leveraged in systems not well characterized for specific variants associated with disease outcome to “capture” specific genes and regulatory regions known to influence expression and disease outcome. We developed a multiplexed, sequence capture assay for red foxes to simultaneously assess ~300‐kbp of genomic sequence from 116 adaptive, intrinsic, and innate immunity genes of predicted adaptive significance and their putative upstream regulatory regions along with 23 neutral microsatellite regions to control for demographic effects. The assay was applied to 45 fox DNA samples from Alaska, where three arctic rabies strains are geographically restricted and endemic to coastal tundra regions, yet absent from the boreal interior. The assay provided 61.5% on‐target enrichment with relatively even sequence coverage across all targeted loci and samples (mean = 50×), which allowed us to elucidate genetic variation across introns, exons, and potential regulatory regions (4,819 SNPs). Challenges remained in accurately describing microsatellite variation using this technique; however, longer‐read HTS technologies should overcome these issues. We used these data to conduct preliminary analyses and detected genetic structure in a subset of red fox immune‐related genes between regions with and without endemic arctic rabies. This assay provides a template to assess immunogenetic variation in wildlife disease systems.

Profiling the immunome of little brown myotis provides a yardstick for measuring the genetic response to white-nose syndrome

White‐nose syndrome (WNS) has devastated populations of hibernating bats in eastern North America, leading to emergency conservation listings for several species including the previously ubiquitous little brown myotis (Myotis lucifugus). However, some bat populations near the epicenter of the WNS panzootic appear to be stabilizing after initial precipitous declines, which could reflect a selective immunogenetic sweep. To investigate the hypothesis that WNS exerts significant selection on the immunome of affected bat populations, we developed a novel, high‐throughput sequence capture assay targeting 138 adaptive, intrinsic, and innate immunity genes of putative adaptive significance, as well as their respective regulatory regions (~370 kbp of genomic sequence/individual). We used the assay to explore baseline immunogenetic variation in M. lucifugus and to investigate whether particular immune genes/variants are associated with WNS susceptibility. We also used our assay to detect 1,038 putatively neutral single nucleotide polymorphisms and characterize contemporary population structure, providing context for the identification of local immunogenetic adaptation. Sequence capture provided a cost‐effective, “all‐in‐one” assay to test for neutral genetic and immunogenetic structure and revealed fine‐scale, baseline immunogenetic differentiation between sampling sites <600 km apart. We identified functional immunogenetic variants in M. lucifugus associated with WNS susceptibility. This study lays the foundations for future investigations of rangewide immunogenetic adaptation to WNS in M. lucifugus and provides a blueprint for studies of evolutionary rescue in other host–pathogen systems.

EPAS1 variants in high altitude Tibetan wolves were selectively introgressed into highland dogs

BACKGROUND

Admixture can facilitate adaptation. For example, black wolves have obtained the variant causing black coat color through past hybridization with domestic dogs and have higher fitness than gray colored wolves. Another recent example of the transfer of adaptive variation between the two species has been suggested by the similarity between high altitude Tibetan mastiffs and wolves at the EPAS1 gene, a transcription factor induced in low oxygen environments.

METHODS

Here, we investigate the directionality of admixture in EPAS1 between 28 reference highland gray wolves, 15 reference domestic dogs, and 21 putatively admixed highland wolves. This experimental design represents an expanded sample of Asian dogs and wolves from previous studies. Admixture was inferred using 17,709 publicly available SNP genotypes on canine chromosome 10. We additionally conducted a scan for positive selection in the highland dog genome.

RESULTS

We find an excess of highland gray wolf ancestry at the EPAS1 locus in highland domestic dogs, suggesting adaptive introgression from wolves to dogs. The signal of admixture is limited in genomic extent to a small region on chromosome 10, indicating that it is the focus of selection in an oxygen-limited environment.

DISCUSSION

Our results suggest that an adaptive variant of EPAS1 in highland wolves was transferred to highland dogs, carrying linked variants that potentially function in hypoxia response at high elevation. The intertwined history of dogs and wolves ensures a unique evolutionary dynamic where variants that have appeared in the history of either species can be tested for their effects on fitness under natural and artificial selection. Such coupled evolutionary histories may be key to the persistence of wild canines and their domesticated kin given the increasing anthropogenic modifications that characterize the future of both species.

Targeted capture to assess neutral genomic variation in the narrow-leaf hopbush across a continental biodiversity refugium

The Adelaide geosyncline, a mountainous region in central southern Australia, is purported to be an important continental refugium for Mediterranean and semi-arid Australian biota, yet few population genetic studies have been conducted to test this theory. Here, we focus on a plant species distributed widely throughout the region, the narrow-leaf hopbush, Dodonaea viscosa ssp. angustissima, and examine its genetic diversity and population structure. We used a hybrid-capture target enrichment technique to selectively sequence over 700 genes from 89 individuals across 17 sampling locations. We compared 815 single nucleotide polymorphisms among individuals and populations to investigate population genetic structure. Three distinct genetic clusters were identified; a Flinders/Gammon ranges cluster, an Eastern cluster, and a Kangaroo Island cluster. Higher genetic diversity was identified in the Flinders/Gammon Ranges cluster, indicating that this area is likely to have acted as a refugium during past climate oscillations. We discuss these findings and consider the historical range dynamics of these populations. We also provide methodological considerations for population genomics studies that aim to use novel genomic approaches (such as target capture methods) on non-model systems. The application of our findings to restoration of this species across the region are also considered.

Spatial patterns of immunogenetic and neutral variation underscore the conservation value of small, isolated American badger populations

Small and isolated populations often exhibit low genetic diversity due to drift and inbreeding, but may simultaneously harbour adaptive variation. We investigate spatial distributions of immunogenetic variation in American badger subspecies (Taxidea taxus), as a proxy for evaluating their evolutionary potential across the northern extent of the species' range. We compared genetic structure of 20 microsatellites and the major histocompatibility complex (MHC DRB exon 2) to evaluate whether small, isolated populations show low adaptive polymorphism relative to large and well-connected populations. Our results suggest that gene flow plays a prominent role in shaping MHC polymorphism across large spatial scales, while the interplay between gene flow and selection was stronger towards the northern peripheries. The similarity of MHC alleles within subspecies relative to their neutral genetic differentiation suggests that adaptive divergence among subspecies can be maintained despite ongoing gene flow along subspecies boundaries. Neutral genetic diversity was low in small relative to large populations, but MHC diversity within individuals was high in small populations. Despite reduced neutral genetic variation, small and isolated populations harbour functional variation that likely contribute to the species evolutionary potential at the northern range. Our findings suggest that conservation approaches should focus on managing adaptive variation across the species range rather than protecting subspecies per se.