Genomics advances the study of inbreeding depression in the wild

Genomic consequences of intensive inbreeding in an isolated wolf population

Inbreeding (mating between relatives) is a major concern for conservation as it decreases individual fitness and can increase the risk of population extinction. We used whole-genome resequencing of 97 grey wolves (Canis lupus) from the highly inbred Scandinavian wolf population to identify 'identical-by-descent' (IBD) chromosome segments as runs of homozygosity (ROH). This gave the high resolution required to precisely measure realized inbreeding as the IBD fraction of the genome in ROH (F _ROH). We found a striking pattern of complete or near-complete homozygosity of entire chromosomes in many individuals. The majority of individual inbreeding was due to long IBD segments (>5 cM) originating from ancestors ≤10 generations ago, with 10 genomic regions showing very few ROH and forming candidate regions for containing loci contributing strongly to inbreeding depression. Inbreeding estimated with an extensive pedigree (F _P) was strongly correlated with realized inbreeding measured with the entire genome (r ² = 0.86). However, inbreeding measured with the whole genome was more strongly correlated with multi-locus heterozygosity estimated with as few as 500 single nucleotide polymorphisms, and with F _ROH estimated with as few as 10,000 single nucleotide polymorphisms, than with F _P. These results document in fine detail the genomic consequences of intensive inbreeding in a population of conservation concern.

Population genomics analyses of European ibex species show lower diversity and higher inbreeding in reintroduced populations

Restoration of lost species ranges to their native distribution is key for the survival of endangered species. However, reintroductions often fail and long‐term genetic consequences are poorly understood. Alpine ibex (Capra ibex) are wild goats that recovered from <100 individuals to ~50,000 within a century by population reintroductions. We analyzed the population genomic consequences of the Alpine ibex reintroduction strategy. We genotyped 101,822 genomewide single nucleotide polymorphism loci in 173 Alpine ibex, the closely related Iberian ibex (Capra pyrenaica) and domestic goat (Capra hircus). The source population of all Alpine ibex maintained genetic diversity comparable to Iberian ibex, which experienced less severe bottlenecks. All reintroduced Alpine ibex populations had individually and combined lower levels of genetic diversity than the source population. The reintroduction strategy consisted of primary reintroductions from captive breeding and secondary reintroductions from established populations. This stepwise reintroduction strategy left a strong genomic footprint of population differentiation, which increased with subsequent rounds of reintroductions. Furthermore, analyses of genomewide runs of homozygosity showed recent inbreeding primarily in individuals of reintroduced populations. We showed that despite the rapid census recovery, Alpine ibex carry a persistent genomic signature of their reintroduction history. We discuss how genomic monitoring can serve as an early indicator of inbreeding.

New-age ideas about age-old sex: separating meiosis from mating could solve a century-old conundrum

Ever since Darwin first addressed it, sexual reproduction reigns as the 'queen' of evolutionary questions. Multiple theories tried to explain how this apparently costly and cumbersome method has become the universal mode of eukaryote reproduction. Most theories stress the adaptive advantages of sex by generating variation, they fail however to explain the ubiquitous persistence of sexual reproduction also where adaptation is not an issue. I argue that the obstacle for comprehending the role of sex stems from the conceptual entanglement of two distinct processes - gamete production by meiosis and gamete fusion by mating (mixis). Meiosis is an ancient, highly rigid and evolutionary conserved process identical and ubiquitous in all eukaryotes. Mating, by contrast, shows tremendous evolutionary variability even in closely related clades and exhibits wonderful ecological adaptability. To appreciate the respective roles of these two processes, which are normally linked and alternating, we require cases where one takes place without the other. Such cases are rather common. The heteromorphic sex chromosomes Y and W, that do not undergo meiotic recombination are an evolutionary test case for demonstrating the role of meiosis. Substantial recent genomic evidence highlights the accelerated rates of change and attrition these chromosomes undergo in comparison to those of recombining autosomes. I thus propose that the most basic role of meiosis is conserving integrity of the genome. A reciprocal case of meiosis without bi-parental mating, is presented by self-fertilization, which is fairly common in flowering plants, as well as most types of apomixis. I argue that deconstructing sex into these two distinct processes - meiosis and mating - will greatly facilitate their analysis and promote our understanding of sexual reproduction.

A model-based approach to characterize individual inbreeding at both global and local genomic scales

Inbreeding results from the mating of related individuals and may be associated with reduced fitness because it brings together deleterious variants in one individual. In general, inbreeding is estimated with respect to an arbitrary base population consisting of ancestors that are assumed unrelated. We herein propose a model-based approach to estimate and characterize individual inbreeding at both global and local genomic scales by assuming the individual genome is a mosaic of homozygous-by-descent (HBD) and non-HBD segments. The HBD segments may originate from ancestors tracing back to different periods in the past defining distinct age-related classes. The lengths of the HBD segments are exponentially distributed with class-specific parameters reflecting that inbreeding of older origin generates on average shorter stretches of observed homozygous markers. The model is implemented in a hidden Markov model framework that uses marker allele frequencies, genetic distances, genotyping error rates and the sequences of observed genotypes. Note that genotyping errors, low-fold sequencing or genotype-by-sequencing data are easily accommodated under this framework. Based on simulations under the inference model, we show that the genomewide inbreeding coefficients and the parameters of the model are accurately estimated. In addition, when several inbreeding classes are simulated, the model captures them if their ages are sufficiently different. Complementary analyses, either on data sets simulated under more realistic models or on human, dog and sheep real data, illustrate the range of applications of the approach and how it can reveal recent demographic histories among populations (e.g., very recent bottlenecks or founder effects). The method also allows to clearly identify individuals resulting from extreme consanguineous matings.

Whole-genome sequencing approaches for conservation biology: Advantages, limitations and practical recommendations

Whole-genome resequencing (WGR) is a powerful method for addressing fundamental evolutionary biology questions that have not been fully resolved using traditional methods. WGR includes four approaches: the sequencing of individuals to a high depth of coverage with either unresolved or resolved haplotypes, the sequencing of population genomes to a high depth by mixing equimolar amounts of unlabelled-individual DNA (Pool-seq) and the sequencing of multiple individuals from a population to a low depth (lcWGR). These techniques require the availability of a reference genome. This, along with the still high cost of shotgun sequencing and the large demand for computing resources and storage, has limited their implementation in nonmodel species with scarce genomic resources and in fields such as conservation biology. Our goal here is to describe the various WGR methods, their pros and cons and potential applications in conservation biology. WGR offers an unprecedented marker density and surveys a wide diversity of genetic variations not limited to single nucleotide polymorphisms (e.g., structural variants and mutations in regulatory elements), increasing their power for the detection of signatures of selection and local adaptation as well as for the identification of the genetic basis of phenotypic traits and diseases. Currently, though, no single WGR approach fulfils all requirements of conservation genetics, and each method has its own limitations and sources of potential bias. We discuss proposed ways to minimize such biases. We envision a not distant future where the analysis of whole genomes becomes a routine task in many nonmodel species and fields including conservation biology.

A High-Density Linkage Map Reveals Sexual Dimorphism in Recombination Landscapes in Red Deer (Cervus elaphus)

High-density linkage maps are an important tool to gain insight into the genetic architecture of traits of evolutionary and economic interest, and provide a resource to characterize variation in recombination landscapes. Here, we used information from the cattle genome and the 50 K Cervine Illumina BeadChip to inform and refine a high-density linkage map in a wild population of red deer (Cervus elaphus). We constructed a predicted linkage map of 38,038 SNPs and a skeleton map of 10,835 SNPs across 34 linkage groups. We identified several chromosomal rearrangements in the deer lineage relative to sheep and cattle, including six chromosome fissions, one fusion, and two large inversions. Otherwise, our findings showed strong concordance with map orders in the cattle genome. The sex-averaged linkage map length was 2739.7 cM and the genome-wide autosomal recombination rate was 1.04 cM/Mb. The female autosomal map length was 1.21 longer than that of males (2767.4 cM vs. 2280.8 cM, respectively). Sex differences in map length were driven by high female recombination rates in peri-centromeric regions, a pattern that is unusual relative to other mammal species. This effect was more pronounced in fission chromosomes that would have had to produce new centromeres. We propose two hypotheses to explain this effect: (1) that this mechanism may have evolved to counteract centromeric drive associated with meiotic asymmetry in oocyte production; and/or (2) that sequence and structural characteristics suppressing recombination in close proximity to the centromere may not have evolved at neo-centromeres. Our study provides insight into how recombination landscapes vary and evolve in mammals, and will provide a valuable resource for studies of evolution, genetic improvement, and population management in red deer and related species.

Genomic Variation of Inbreeding and Ancestry in the Remaining Two Isle Royale Wolves

Inbreeding, relatedness, and ancestry have traditionally been estimated with pedigree information, however, molecular genomic data can provide more detailed examination of these properties. For example, pedigree information provides estimation of the expected value of these measures but molecular genomic data can estimate the realized values of these measures in individuals. Here, we generate the theoretical distribution of inbreeding, relatedness, and ancestry for the individuals in the pedigree of the Isle Royale wolves, the first examination of such variation in a wild population with a known pedigree. We use the 38 autosomes of the dog genome and their estimated map lengths in our genomic analysis. Although it is known that the remaining wolves are highly inbred, closely related, and descend from only 3 ancestors, our analyses suggest that there is significant variation in the realized inbreeding and relatedness around pedigree expectations. For example, the expected inbreeding in a hypothetical offspring from the 2 remaining wolves is 0.438 but the realized 95% genomic confidence interval is from 0.311 to 0.565. For individual chromosomes, a substantial proportion of the whole chromosomes are completely identical by descent. This examination provides a background to use when analyzing molecular genomic data for individual levels of inbreeding, relatedness, and ancestry. The level of variation in these measures is a function of the time to the common ancestor(s), the number of chromosomes, and the rate of recombination. In the Isle Royale wolf population, the few generations to a common ancestor results in the high variance in genomic inbreeding.

No evidence of inbreeding depression in a Tasmanian devil insurance population despite significant variation in inbreeding

Inbreeding depression occurs when inbred individuals experience reduced fitness as a result of reduced genome-wide heterozygosity. The Tasmanian devil faces extinction due to a contagious cancer, devil facial tumour disease (DFTD). An insurance metapopulation was established in 2006 to ensure the survival of the species and to be used as a source population for re-wilding and genetic rescue. The emergence of DFTD and the rapid decline of wild devil populations have rendered the species at risk of inbreeding depression. We used 33 microsatellite loci to (1) reconstruct a pedigree for the insurance population and (2) estimate genome-wide heterozygosity for 200 individuals. Using heterozygosity-fitness correlations, we investigated the effect of heterozygosity on six diverse fitness measures (ulna length, asymmetry, weight-at-weaning, testes volume, reproductive success and survival). Despite statistically significant evidence of variation in individual inbreeding in this population, we found no associations between inbreeding and any of our six fitness measurements. We propose that the benign environment in captivity may decrease the intensity of inbreeding depression, relative to the stressful conditions in the wild. Future work will need to measure fitness of released animals to facilitate translation of this data to the broader conservation management of the species in its native range.

The distribution of runs of homozygosity and selection signatures in six commercial meat sheep breeds

Domestication and the subsequent selection of animals for either economic or morphological features can leave a variety of imprints on the genome of a population. Genomic regions subjected to high selective pressures often show reduced genetic diversity and frequent runs of homozygosity (ROH). Therefore, the objective of the present study was to use 42,182 autosomal SNPs to identify genomic regions in 3,191 sheep from six commercial breeds subjected to selection pressure and to quantify the genetic diversity within each breed using ROH. In addition, the historical effective population size of each breed was also estimated and, in conjunction with ROH, was used to elucidate the demographic history of the six breeds. ROH were common in the autosomes of animals in the present study, but the observed breed differences in patterns of ROH length and burden suggested differences in breed effective population size and recent management. ROH provided a sufficient predictor of the pedigree inbreeding coefficient, with an estimated correlation between both measures of 0.62. Genomic regions under putative selection were identified using two complementary algorithms; the fixation index and hapFLK. The identified regions under putative selection included candidate genes associated with skin pigmentation, body size and muscle formation; such characteristics are often sought after in modern-day breeding programs. These regions of selection frequently overlapped with high ROH regions both within and across breeds. Multiple yet uncharacterised genes also resided within putative regions of selection. This further substantiates the need for a more comprehensive annotation of the sheep genome as these uncharacterised genes may contribute to traits of interest in the animal sciences. Despite this, the regions identified as under putative selection in the current study provide an insight into the mechanisms leading to breed differentiation and genetic variation in meat production.

Inferring Individual Inbreeding and Demographic History from Segments of Identity by Descent in Ficedula Flycatcher Genome Sequences

Individual inbreeding and historical demography can be estimated by analyzing runs of homozygosity (ROH), which are indicative of chromosomal segments of identity by descent (IBD). Such analyses have so far been rare in natural populations due to limited genomic resources. We analyzed ROH in whole genome sequences from 287 Ficedula flycatchers representing four species, with the objectives of evaluating the causes of genome-wide variation in the abundance of ROH and inferring historical demography. ROH were clearly more abundant in genomic regions with low recombination rate. However, this pattern was substantially weaker when ROH were mapped using genetic rather than physical single nucleotide polymorphism (SNP) coordinates in the genome. Empirical results and simulations suggest that high ROH abundance in regions of low recombination was partly caused by increased power to detect the very long IBD segments typical of regions with a low recombination rate. Simulations also showed that hard selective sweeps (but not soft sweeps or background selection) likely contributed to variation in the abundance of ROH across the genome. Comparisons of the abundance of ROH among several study populations indicated that the Spanish pied flycatcher population had the smallest historical effective population size (N_e) for this species, and that a putatively recently founded island (Baltic) population had the smallest historical N_e among the collared flycatchers. Analysis of pairwise IBD in Baltic collared flycatchers indicated that this population was founded <60 generations ago. This study provides a rare genomic glimpse into demographic history and the mechanisms underlying the genome-wide distribution of ROH.