RAD-sequencing for estimating genomic relatedness matrix-based heritability in the wild: A case study in roe deer

Next-generation data filtering in the genomics era

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

Dissecting the genetic architecture of quantitative traits using genome-wide identity-by-descent sharing

Additive and dominance genetic variances underlying the expression of quantitative traits are important quantities for predicting short-term responses to selection, but they are notoriously challenging to estimate in most non-model wild populations. Specifically, large-sized or panmictic populations may be characterized by low variance in genetic relatedness among individuals which, in turn, can prevent accurate estimation of quantitative genetic parameters. We used estimates of genome-wide identity-by-descent (IBD) sharing from autosomal SNP loci to estimate quantitative genetic parameters for ecologically important traits in nine-spined sticklebacks (Pungitius pungitius) from a large, outbred population. Using empirical and simulated datasets, with varying sample sizes and pedigree complexity, we assessed the performance of different crossing schemes in estimating additive genetic variance and heritability for all traits. We found that low variance in relatedness characteristic of wild outbred populations with high migration rate can impair the estimation of quantitative genetic parameters and bias heritability estimates downwards. On the other hand, the use of a half-sib/full-sib design allowed precise estimation of genetic variance components and revealed significant additive variance and heritability for all measured traits, with negligible dominance contributions. Genome-partitioning and QTL mapping analyses revealed that most traits had a polygenic basis and were controlled by genes at multiple chromosomes. Furthermore, different QTL contributed to variation in the same traits in different populations suggesting heterogeneous underpinnings of parallel evolution at the phenotypic level. Our results provide important guidelines for future studies aimed at estimating adaptive potential in the wild, particularly for those conducted in outbred large-sized populations.

How Can Genomics Help or Hinder Wildlife Conservation?

Genomic data are becoming increasingly affordable and easy to collect, and new tools for their analysis are appearing rapidly. Conservation biologists are interested in using this information to assist in management and planning but are typically limited financially and by the lack of genomic resources available for non-model taxa. It is therefore important to be aware of the pitfalls as well as the benefits of applying genomic approaches. Here, we highlight recent methods aimed at standardizing population assessments of genetic variation, inbreeding, and forms of genetic load and methods that help identify past and ongoing patterns of genetic interchange between populations, including those subjected to recent disturbance. We emphasize challenges in applying some of these methods and the need for adequate bioinformatic support. We also consider the promises and challenges of applying genomic approaches to understand adaptive changes in natural populations to predict their future adaptive capacity.

Contrasting whole-genome and reduced representation sequencing for population demographic and adaptive inference: an alpine mammal case study

Genomes capture the adaptive and demographic history of a species, but the choice of sequencing strategy and sample size can impact such inferences. We compared whole genome and reduced representation sequencing approaches to study the population demographic and adaptive signals of the North American mountain goat (Oreamnos americanus). We applied the restriction site-associated DNA sequencing (RADseq) approach to 254 individuals and whole genome resequencing (WGS) approach to 35 individuals across the species range at mid-level coverage (9X) and to 5 individuals at high coverage (30X). We used ANGSD to estimate the genotype likelihoods and estimated the effective population size (Ne), population structure, and explicitly modelled the demographic history with δaδi and MSMC2. The data sets were overall concordant in supporting a glacial induced vicariance and extremely low Ne in mountain goats. We evaluated a set of climatic variables and geographic location as predictors of genetic diversity using redundancy analysis. A moderate proportion of total variance (36% for WGS and 21% for RADseq data sets) was explained by geography and climate variables; both data sets support a large impact of drift and some degree of local adaptation. The empirical similarities of WGS and RADseq presented herein reassuringly suggest that both approaches will recover large demographic and adaptive signals in a population; however, WGS offers several advantages over RADseq, such as inferring adaptive processes and calculating runs-of-homozygosity estimates. Considering the predicted climate-induced changes in alpine environments and the genetically depauperate mountain goat, the long-term adaptive capabilities of this enigmatic species are questionable.

Genomic signatures of inbreeding in a critically endangered parrot, the kākāpō

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kākāpō (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kākāpō, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

Pathogen-mediated selection favours the maintenance of innate immunity gene polymorphism in a widespread wild ungulate

Toll-like receptors (TLR) play a central role in recognition and host frontline defence against a wide range of pathogens. A number of recent studies have shown that TLR genes (Tlrs) often exhibit large polymorphism in natural populations. Yet, there is little knowledge on how this polymorphism is maintained and how it influences disease susceptibility in the wild. In previous work, we showed that some Tlrs exhibit similarly high levels of genetic diversity as genes of the Major Histocompatibility Complex (MHC), and signatures of contemporary balancing selection in roe deer (Capreolus capreolus), the most abundant cervid species in Europe. Here, we investigated the evolutionary mechanisms by which pathogen-mediated selection could shape this innate immunity genetic diversity by examining the relationships between Tlr (Tlr2, Tlr4 and Tlr5) genotypes (heterozygosity status and presence of specific alleles) and infections with Toxoplasma and Chlamydia, two widespread intracellular pathogens known to cause reproductive failure in ungulates. We showed that Toxoplasma and Chlamydia exposures vary significantly across years and landscape features with few co-infection events detected and that the two pathogens exert antagonistic selection on Tlr2 polymorphism. By contrast, we found limited support for Tlr heterozygote advantage. Our study confirmed the importance of looking beyond Mhc genes in wildlife immunogenetic studies. It also emphasized the necessity to consider multiple pathogen challenges and their spatiotemporal variation to improve our understanding of vertebrate defence evolution against pathogens.

Examining adaptive evolution of immune activity: opportunities provided by gastropods in the age of 'omics'

Parasites threaten all free-living organisms, including molluscs. Understanding the evolution of immune defence traits in natural host populations is crucial for predicting their long-term performance under continuous infection risk. Adaptive trait evolution requires that traits are subject to selection (i.e. contribute to organismal fitness) and that they are heritable. Despite broad interest in the evolutionary ecology of immune activity in animals, the understanding of selection on and evolutionary potential of immune defence traits is far from comprehensive. For instance, empirical observations are only rarely in line with theoretical predictions of immune activity being subject to stabilizing selection. This discrepancy may be because ecoimmunological studies can typically cover only a fraction of the complexity of an animal immune system. Similarly, molecular immunology/immunogenetics studies provide a mechanistic understanding of immunity, but neglect variation that arises from natural genetic differences among individuals and from environmental conditions. Here, we review the current literature on natural selection on and evolutionary potential of immune traits in animals, signal how merging ecological immunology and genomics will strengthen evolutionary ecological research on immunity, and indicate research opportunities for molluscan gastropods for which well-established ecological understanding and/or 'immune-omics' resources are already available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Heritability Estimates of Antler and Body Traits in White-Tailed Deer (Odocoileus virginianus) From Genomic-Relatedness Matrices

Estimating heritability (h2) is required to predict the response to selection and is useful in species that are managed or farmed using trait information. Estimating h2 in free-ranging populations is challenging due to the need for pedigrees; genomic-relatedness matrices (GRMs) circumvent this need and can be implemented in nearly any system where phenotypic and genome-wide single-nucleotide polymorphism (SNP) data are available. We estimated the heritability of 5 body and 3 antler traits in a free-ranging population of white-tailed deer (Odocoileus virginianus) on Anticosti Island, Quebec, Canada. We generated classic and robust GRMs from >10,000 SNPs: hind foot length, dressed body mass, and peroneus muscle mass had high h2 values of 0.62, 0.44, and 0.55, respectively. Heritability in male-only antler features ranged from 0.07 to 0.33. We explored the influence of filtering by minor allele frequency and data completion on h2: GRMs derived from fewer SNPs had reduced h2 estimates and the relatedness coefficients significantly deviated from those generated with more SNPs. As a corollary, we discussed limitations to the application of GRMs in the wild, notably how skewed GRMs, specifically many unrelated individuals, can increase variance around h2 estimates. This is the first study to estimate h2 on a free-ranging population of white-tailed deer and should be informative for breeding designs and management as these traits could respond to selection.

Demography and adaptation promoting evolutionary transitions in a mammalian genus that diversified during the Pleistocene

Species that evolved in temperate regions during the Pleistocene experienced periods of extreme climatic transitions. Consequent population fragmentation and dynamics had the potential to generate small, isolated populations where the influence of genetic drift would be expected to be strong. We use comparative genomics to assess the evolutionary influence of historical demographics and natural selection through a series of transitions associated with the formation of the genus Capreolus, speciation within this genus during the Quaternary and during divergence among European roe deer (C. capreolus) populations. Our analyses were facilitated by the generation of a new high-coverage reference genome for the Siberian roe deer (C. pygargus). We find progressive reductions in effective population size (N_e ), despite very large census sizes in modern C. capreolus populations and show that low N_e has impacted the C. capreolus genome, reducing diversity and increasing linkage disequilibrium. Even so, we find evidence for natural selection shared among C. capreolus populations, including a historically documented founder population that has been through a severe bottleneck. During each phylogenetic transition there is evidence for selection (from dN/dS and nucleotide diversity tests), including at loci associated with diapause (delayed embryonic development), a phenotype restricted to this genus among the even-toed ungulates. Together these data allow us to assess expectations for the origin and diversification of a mammalian genus during a period of extreme environmental change.

Pedigree-free quantitative genetic approach provides evidence for heritability of movement tactics in wild roe deer

Assessing the evolutionary potential of animal populations in the wild is crucial to understanding how they may respond to selection mediated by rapid environmental change (e.g. habitat loss and fragmentation). A growing number of studies have investigated the adaptive role of behaviour, but assessments of its genetic basis in a natural setting remain scarce. We combined intensive biologging technology with genome-wide data and a pedigree-free quantitative genetic approach to quantify repeatability, heritability and evolvability for a suite of behaviours related to the risk avoidance-resource acquisition trade-off in a wild roe deer (Capreolus capreolus) population inhabiting a heterogeneous, human-dominated landscape. These traits, linked to the stress response, movement and space-use behaviour, were all moderately to highly repeatable. Furthermore, the repeatable among-individual component of variation in these traits was partly due to additive genetic variance, with heritability estimates ranging from 0.21 ± 0.08 to 0.70 ± 0.11 and evolvability ranging from 1.1% to 4.3%. Changes in the trait mean can therefore occur under hypothetical directional selection over just a few generations. To the best of our knowledge, this is the first empirical demonstration of additive genetic variation in space-use behaviour in a free-ranging population based on genomic relatedness data. We conclude that wild animal populations may have the potential to adjust their spatial behaviour to human-driven environmental modifications through microevolutionary change.