Maintaining genetic diversity using molecular coancestry: the effect of marker density and effective population size

Background

The most efficient method to maintain genetic diversity in populations under conservation programmes is to optimize, for each potential parent, the number of offspring left to the next generation by minimizing the global coancestry. Coancestry is usually calculated from genealogical data but molecular markers can be used to replace genealogical coancestry with molecular coancestry. Recent studies showed that optimizing contributions based on coancestry calculated from a large number of SNP markers can maintain higher levels of diversity than optimizing contributions based on genealogical data. In this study, we investigated how SNP density and effective population size impact the use of molecular coancestry to maintain diversity.

Results

At low SNP densities, the genetic diversity maintained using genealogical coancestry for optimization was higher than that maintained using molecular coancestry. The performance of molecular coancestry improved with increasing marker density, and, for the scenarios evaluated, it was as efficient as genealogical coancestry if SNP density reached at least 3 times the effective population size.However, increasing SNP density resulted in reduced returns in terms of maintained diversity. While a benefit of 12% was achieved when marker density increased from 10 to 100 SNP/Morgan, the benefit was only 2% when it increased from 100 to 500 SNP/Morgan.

Conclusions

The marker density of most SNP chips already available for farm animals is sufficient for molecular coancestry to outperform genealogical coancestry in conservation programmes aimed at maintaining genetic diversity. For the purpose of effectively maintaining genetic diversity, a marker density of around 500 SNPs/Morgan can be considered as the most cost effective density when developing SNP chips for new species. Since the costs to develop SNP chips are decreasing, chips with 500 SNPs/Morgan should become available in a short-term horizon for non domestic species.     

Genome-Wide Estimates of Coancestry and Inbreeding in a Closed Herd of Ancient Iberian Pigs

Maintaining genetic variation and controlling the increase in inbreeding are crucial requirements in animal conservation programs. The most widely accepted strategy for achieving these objectives is to maximize the effective population size by minimizing the global coancestry obtained from a particular pedigree. However, for most natural or captive populations genealogical information is absent. In this situation, microsatellites have been traditionally the markers of choice to characterize genetic variation, and several estimators of genealogical coefficients have been developed using marker data, with unsatisfactory results. The development of high-throughput genotyping techniques states the necessity of reviewing the paradigm that genealogical coancestry is the best parameter for measuring genetic diversity. In this study, the Illumina PorcineSNP60 BeadChip was used to obtain genome-wide estimates of rates of coancestry and inbreeding and effective population size for an ancient strain of Iberian pigs that is now in serious danger of extinction and for which very accurate genealogical information is available (the Guadyerbas strain). Genome-wide estimates were compared with those obtained from microsatellite and from pedigree data. Estimates of coancestry and inbreeding computed from the SNP chip were strongly correlated with genealogical estimates and these correlations were substantially higher than those between microsatellite and genealogical coefficients. Also, molecular coancestry computed from SNP information was a better predictor of genealogical coancestry than coancestry computed from microsatellites. Rates of change in coancestry and inbreeding and effective population size estimated from molecular data were very similar to those estimated from genealogical data. However, estimates of effective population size obtained from changes in coancestry or inbreeding differed. Our results indicate that genome-wide information represents a useful alternative to genealogical information for measuring and maintaining genetic diversity.     

A comparison of management strategies for conservation with regard to population fitness

Computer simulations have been carried out tocompare, under realistic genetic models, twomethods proposed in the literature to retaingenetic diversity in conservation programmes.In a two-step method, contributions of parentsare set up to produce minimum coancestry(kinship) among the offspring, and this isindependent of the mating system subsequentlyapplied. In a single-step method,contributions and matings are decidedsimultaneously in order to minimise coancestry.The comparison is made in terms of maintainedgenetic diversity and in terms of populationfitness. We conclude that the two methodsmaintain approximately the same geneticdiversity but the latter induces higher levelsof inbreeding, reducing the fitness of thepopulation. Avoidance of close relatives'matings improves this latter method, but thefitness levels do not reach those of thetwo-step scheme. We also investigate theperformances of different mating strategies incombination with minimum coancestry (two-stepmethod), concluding that these mating systemsdo not substantially affect the effectivenessof the management. Finally, we illustrate howminimum group coancestry can be restrictedto a minimum loss of fitness, if a measure ofthis is available for the individuals.     

Genetic management of captive populations: the advantages of circular mating

We performed computer simulations to evaluate the effectiveness of circular mating as a genetic management option for captive populations. As a benchmark, we used the method proposed by Fernández and Caballero according to which parental contributions are set to produce minimum coancestry among the offspring and matings are performed so as to minimize mean pairwise coancestry (referred to as the Gc/mc method). In contrast to other methods, fitness does not vary with population size in the case of circular mating, and can be higher than under random mating. Whether circular mating is an effective method in conserving captive populations depends on the trade-off between different considerations. On the one hand, circular mating shows the highest allelic diversity and the lowest mean pairwise coancestry for all population sizes. It also shows a relatively higher efficiency of purging deleterious alleles. More importantly, circular mating can significantly increase the success probability of populations released to the wild relative to the Gc/mc method. On the other hand, circular mating has the drawback of showing high inbreeding rates and low fitness in early generations, which can result to an increase in the extinction probability of the captive populations. However, this increase is slight unless population size and litter size are both very low. Overall, if the slight increase in extinction probability can be tolerated then circular mating fulfils the primary goals of a captive breeding program, i.e., it maintains high levels of genetic diversity and increases the success probability of reintroduced populations.     

An experimental evaluation with Drosophila melanogaster of a novel dynamic system for the management of subdivided populations in conservation programs

A dynamic method (DM) recently proposed for the management of captive subdivided populations was evaluated using the pilot species Drosophila melanogaster. By accounting for the particular genetic population structure, the DM determines the optimal mating pairs, their contributions to progeny and the migration pattern that minimize the overall coancestry in the population with a control of inbreeding levels. After a pre-management period such that one of the four subpopulations had higher inbreeding and differentiation than the others, three management methods were compared for 10 generations over three replicates: (1) isolated subpopulations (IS), (2) one-migrant-per-generation rule (OMPG), (3) DM aimed to produce the same or lower inbreeding coefficient than OMPG. The DM produced the lowest coancestry and equal or lower inbreeding than the OMPG method throughout the experiment. The initially lower fitness and lower variation for nine microsatellite loci of the highly inbred subpopulation were restored more quickly with the DM than with the OMPG method. We provide, therefore, an empirical illustration of the usefulness of the DM as a conservation protocol for captive subdivided populations when pedigree information is available (or can be deduced) and manipulation of breeding pairs is possible.     

Management of genetic diversity in small farm animal populations

Many local breeds of farm animals have small populations and, consequently, are highly endangered. The correct genetic management of such populations is crucial for their survival. Managing an animal population involves two steps: first, the individuals who will be permitted to leave descendants are to be chosen and the number offspring they will be permitted to produce has to be determined; second, the mating scheme has to be identified. Strategies dealing with the first step are directed towards the maximisation of effective population size and, therefore, act jointly on the reduction in the loss of genetic variation and in the increase of inbreeding. In this paper, the most relevant methods are summarised, including the so-called 'Optimum Contribution' methodology (contributions are proportional to the coancestry of each individual with the rest), which has been shown to be the best. Typically, this method is applied to pedigree information, but molecular marker data can be used to complete or replace the genealogy. When the population is subjected to explicit selection on any trait, the above methodology can be used by balancing the response to selection and the increase in coancestry/inbreeding. Different mating strategies also exist. Some of the mating schemes try to reduce the level of inbreeding in the short term by preventing mating between relatives. Others involve regular (circular) schemes that imply higher levels of inbreeding within populations in the short term, but demonstrate better performance in the long term. In addition, other tools such as cryopreservation and reproductive techniques aid in the management of small populations. In the future, genomic marker panels may replace the pedigree information in measuring the coancestry. The paper also includes the results of several experiments and field studies on the effectiveness and on the consequences of the use of the different strategies.     

Management of subdivided populations in conservation programs: development of a novel dynamic system

Within the context of a conservation program the management of subdivided populations implies a compromise between the control of the global genetic diversity, the avoidance of high inbreeding levels, and, sometimes, the maintenance of a certain degree of differentiation between subpopulations. We present a dynamic and flexible methodology, based on genealogical information, for the maximization of the genetic diversity (measured through the global population coancestry) in captive subdivided populations while controlling/restricting the levels of inbreeding. The method is able to implement specific restrictions on the desired relative levels of coancestry between and within subpopulations. By accounting for the particular genetic population structure, the method determines the optimal contributions (i.e., number of offspring) of each individual, the number of migrants, and the particular subpopulations involved in the exchange of individuals. Computer simulations are used to illustrate the procedure and its performance in a range of reasonable scenarios. The method performs well in most situations and is shown to be more efficient than the commonly accepted one-migrant-per-generation strategy.     

Invasion success and genetic diversity of introduced populations of guppies Poecilia reticulata in Australia

High genetic diversity is thought to characterize successful invasive species, as the potential to adapt to new environments is enhanced and inbreeding is reduced. In the last century, guppies, Poecilia reticulata, repeatedly invaded streams in Australia and elsewhere. Quantitative genetic studies of one Australian guppy population have demonstrated high additive genetic variation for autosomal and Y-linked morphological traits. The combination of colonization success, high heritability of morphological traits, and the possibility of multiple introductions to Australia raised the prediction that neutral genetic diversity is high in introduced populations of guppies. In this study we examine genetic diversity at nine microsatellite and one mitochondrial locus for seven Australian populations. We used mtDNA haplotypes from the natural range of guppies and from domesticated varieties to identify source populations. There were a minimum of two introductions, but there was no haplotype diversity within Australian populations, suggesting a founder effect. This was supported by microsatellite markers, as allelic diversity and heterozygosity were severely reduced compared to one wild source population, and evidence of recent bottlenecks was found. Between Australian populations little differentiation of microsatellite allele frequencies was detected, suggesting that population admixture has occurred historically, perhaps due to male-biased gene flow followed by bottlenecks. Thus success of invasion of Australia and high additive genetic variance in Australian guppies are not associated with high levels of diversity at molecular loci. This finding is consistent with the release of additive genetic variation by dominance and epistasis following inbreeding, and with disruptive and negative frequency-dependent selection on fitness traits.     

Efficiency of the use of pedigree and molecular marker information in conservation programs

The value of molecular markers and pedigree records, separately or in combination, to assist in the management of conserved populations has been tested. The general strategy for managing the population was to optimize contributions of parents to the next generation for minimizing the global weighted coancestry. Strategies differed in the type of information used to compute global coancestries, the number and type of evaluated individuals, and the system of mating. Genealogical information proved to be very useful (at least for 10 generations of management) to arrange individuals' contributions via the minimization of global coancestry. In fact, the level of expected heterozygosity after 10 generations yielded by this strategy was 88-100% of the maximum possible improvement obtained if the genotype for all loci was known. Marker information was of very limited value if used alone. The amount and degree of polymorphism of markers to be used to compute molecular coancestry had to be high to mimic the performance of the strategy relying on pedigree, especially in the short term (for example, >10 markers per chromosome with 10 alleles each were needed if only the parents' genotype was available). When both sources of information are combined to calculate the coancestry conditional on markers, clear increases in effective population size (Ne) were found, but observed diversity levels (either gene or allelic diversity) in the early generations were quite similar to the ones obtained with pedigree alone. The advantage of including molecular information is greater when information is available on a greater number of individuals (offspring and parents vs. parents only). However, for realistic situations (i.e., large genomes) the benefits of using information on offspring are small. The same conclusions were reached when comparing the use of the different types of information (genealogical or/and molecular) to perform minimum coancestry matings.     

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.     

Genetic diversity at neutral and adaptive loci determines individual fitness in a long-lived territorial bird

There is compelling evidence about the manifest effects of inbreeding depression on individual fitness and populations' risk of extinction. The majority of studies addressing inbreeding depression on wild populations are generally based on indirect measures of inbreeding using neutral markers. However, the study of functional loci, such as genes of the major histocompatibility complex (MHC), is highly recommended. MHC genes constitute an essential component of the immune system of individuals, which is directly related to individual fitness and survival. In this study, we analyse heterozygosity fitness correlations of neutral and adaptive genetic variation (22 microsatellite loci and two loci of the MHC class II, respectively) with the age of recruitment and breeding success of a decimated and geographically isolated population of a long-lived territorial vulture. Our results indicate a negative correlation between neutral genetic diversity and age of recruitment, suggesting that inbreeding may be delaying reproduction. We also found a positive correlation between functional (MHC) genetic diversity and breeding success, together with a specific positive effect of the most frequent pair of cosegregating MHC alleles in the population. Globally, our findings demonstrate that genetic depauperation in small populations has a negative impact on the individual fitness, thus increasing the populations' extinction risk.     

Removing exogenous information using pedigree data

Management of certain populations requires the preservation of its pure genetic background. When, for different reasons, undesired alleles are introduced, the original genetic conformation must be recovered. The present study tested, through computer simulations, the power of recovery (the ability for removing the foreign information) from genealogical data. Simulated scenarios comprised different numbers of exogenous individuals taking part of the founder population and different numbers of unmanaged generations before the removal program started. Strategies were based on variables arising from classical pedigree analyses such as founders’ contribution and partial coancestry. The efficiency of the different strategies was measured as the proportion of native genetic information remaining in the population. Consequences on the inbreeding and coancestry levels of the population were also evaluated. Minimisation of the exogenous founders’ contributions was the most powerful method, removing the largest amount of genetic information in just one generation. However, as a side effect, it led to the highest values of inbreeding. Scenarios with a large amount of initial exogenous alleles (i.e. high percentage of non native founders), or many generations of mixing became very difficult to recover, pointing out the importance of being careful about introgression events in populations where these are undesired.     

Native breeds demonstrate high contributions to the molecular variation in northern European sheep

Population contribution to genetic diversity can be estimated using neutral variation. However, population expansion or hybridization of diverged ancestries may weaken correlation between neutral and non-neutral variation. Microsatellite variation was studied at 25 loci in 20 native and 12 modern or imported northern European sheep breeds. Breed contributions to total gene diversity, allelic richness and mean allele-sharing distance between individuals were measured. Indications of changes in population size and admixtures of divergent ancestries were investigated and the extent of inbreeding was estimated. The northern European sheep demonstrated signs of reduction in effective population size. Many old, small populations made a substantial positive contribution to total molecular variation, but populations with several divergent major ancestries did not contribute substantially to molecular variation, with the exception of the Norwegian Rygja sheep. However, several diverged major ancestries may cause it to contribute less to non-neutral variation than expected from the microsatellite data. Breed uniqueness and within-breed variability generally had opposite effects on breed contributions to molecular diversity. The degree of inbreeding did not reflect the breed contribution to total gene diversity or allelic richness, but inbred populations increased the mean allele-sharing distance between individuals. Our study indicates breed conservation to be especially important in maintaining allelic variation in northern European sheep and supports the evolutionary importance of peripheral populations.     

Genome-Wide Estimates of Coancestry,Inbreeding and Effective Population Size in the Spanish Holstein Population

Estimates of effective population size in the Holstein cattle breed have usually been low despite the large number of animals that constitute this breed. Effective population size is inversely related to the rates at which coancestry and inbreeding increase and these rates have been high as a consequence of intense and accurate selection. Traditionally, coancestry and inbreeding coefficients have been calculated from pedigree data. However, the development of genome-wide single nucleotide polymorphisms has increased the interest of calculating these coefficients from molecular data in order to improve their accuracy. In this study, genomic estimates of coancestry, inbreeding and effective population size were obtained in the Spanish Holstein population and then compared with pedigree-based estimates. A total of 11,135 animals genotyped with the Illumina BovineSNP50 BeadChip were available for the study. After applying filtering criteria, the final genomic dataset included 36,693 autosomal SNPs and 10,569 animals. Pedigree data from those genotyped animals included 31,203 animals. These individuals represented only the last five generations in order to homogenise the amount of pedigree information across animals. Genomic estimates of coancestry and inbreeding were obtained from identity by descent segments (coancestry) or runs of homozygosity (inbreeding). The results indicate that the percentage of variance of pedigree-based coancestry estimates explained by genomic coancestry estimates was higher than that for inbreeding. Estimates of effective population size obtained from genome-wide and pedigree information were consistent and ranged from about 66 to 79. These low values emphasize the need of controlling the rate of increase of coancestry and inbreeding in Holstein selection programmes.     

Molecular and pedigree measures of relatedness provide similar estimates of inbreeding depression in a bottlenecked population

Individual‐based estimates of the degree of inbreeding or parental relatedness from pedigrees provide a critical starting point for studies of inbreeding depression, but in practice wild pedigrees are difficult to obtain. Because inbreeding increases the proportion of genomewide loci that are identical by descent, inbreeding variation within populations has the potential to generate observable correlations between heterozygosity measured using molecular markers and a variety of fitness related traits. Termed heterozygosity‐fitness correlations (HFCs), these correlations have been observed in a wide variety of taxa. The difficulty of obtaining wild pedigree data, however, means that empirical investigations of how pedigree inbreeding influences HFCs are rare. Here, we assess evidence for inbreeding depression in three life‐history traits (hatching and fledging success and juvenile survival) in an isolated population of Stewart Island robins using both pedigree‐ and molecular‐derived measures of relatedness. We found results from the two measures were highly correlated and supported evidence for significant but weak inbreeding depression. However, standardized effect sizes for inbreeding depression based on the pedigree‐based kin coefficients (k) were greater and had smaller standard errors than those based on molecular genetic measures of relatedness (RI), particularly for hatching and fledging success. Nevertheless, the results presented here support the use of molecular‐based measures of relatedness in bottlenecked populations when information regarding inbreeding depression is desired but pedigree data on relatedness are unavailable.     

A note on the rationale for estimating genealogical coancestry from molecular markers

Background

Genetic relatedness or similarity between individuals is a key concept in population, quantitative and conservation genetics. When the pedigree of a population is available and assuming a founder population from which the genealogical records start, genetic relatedness between individuals can be estimated by the coancestry coefficient. If pedigree data is lacking or incomplete, estimation of the genetic similarity between individuals relies on molecular markers, using either molecular coancestry or molecular covariance. Some relationships between genealogical and molecular coancestries and covariances have already been described in the literature.

Methods

We show how the expected values of the empirical measures of similarity based on molecular marker data are functions of the genealogical coancestry. From these formulas, it is easy to derive estimators of genealogical coancestry from molecular data. We include variation of allelic frequencies in the estimators.

Results

The estimators are illustrated with simulated examples and with a real dataset from dairy cattle. In general, estimators are accurate and only slightly biased. From the real data set, estimators based on covariances are more compatible with genealogical coancestries than those based on molecular coancestries. A frequently used estimator based on the average of estimated coancestries produced inflated coancestries and numerical instability. The consequences of unknown gene frequencies in the founder population are briefly discussed, along with alternatives to overcome this limitation.

Conclusions

Estimators of genealogical coancestry based on molecular data are easy to derive. Estimators based on molecular covariance are more accurate than those based on identity by state. A correction considering the random distribution of allelic frequencies improves accuracy of these estimators, especially for populations with very strong drift.     

On the expected relationship between inbreeding,fitness, and extinction

We assessed the expected relationship between the level and the cost of inbreeding, measured either in terms of fitness, inbreeding depression or probability of extinction. First, we show that the assumption of frequent, slightly deleterious mutations do agree with observations and experiments, on the contrary to the assumption of few, moderately deleterious mutations. For the same inbreeding coefficient, populations can greatly differ in fitness according to the following: (i) population size; larger populations show higher fitness (ii) the history of population size; in a population that recovers after a bottleneck, higher inbreeding can lead to higher fitness and (iii) population demography; population growth rate and carrying capacity determine the relationship between inbreeding and extinction. With regards to the relationship between inbreeding depression and inbreeding coefficient, the population size that minimizes inbreeding depression depends on the level of inbreeding: inbreeding depression can even decrease when population size increases. It is therefore clear that to infer the costs of inbreeding, one must know both the history of inbreeding (e.g. past bottlenecks) and population demography.     

Positive assortative mating with selection restrictions on group coancestry enhances gain while conserving genetic diversity in long-term forest tree breeding

Selection and mating principles in a closed breeding population (BP) were studied by computer simulation. The BP was advanced, either by random assortment of mates (RAM), or by positive assortative mating (PAM). Selection was done with high precision using clonal testing. Selection considered both genetic gain and gene diversity by "group-merit selection", i.e. selection for breeding value weighted by group coancestry of the selected individuals. A range of weights on group coancestry was applied during selection to vary parent contributions and thereby adjust the balance between gain and diversity. This resulted in a series of scenarios with low to high effective population sizes measured by status effective number. Production populations (PP) were selected only for gain, as a subset of the BP. PAM improved gain in the PP substantially, by increasing the additive variance (i.e. the gain potential) of the BP. This effect was more pronounced under restricted selection when parent contributions to the next generation were more balanced with within-family selection as the extreme, i.e. when a higher status effective number was maintained in the BP. In that case, the additional gain over the BP mean for the clone PP and seed PPs was 32 and 84% higher, respectively, for PAM than for RAM in generation 5. PAM did not reduce gene diversity of the BP but increased inbreeding, and in that way caused a departure from Hardy-Weinberg equilibrium. The effect of inbreeding was eliminated by recombination during the production of seed orchard progeny. Also, for a given level of inbreeding in the seed orchard progeny or in a mixture of genotypes selected for clonal deployment, gain was higher for PAM than for RAM. After including inbreeding depression in the simulation, inbreeding was counteracted by selection, and the enhancement of PAM on production population gain was slightly reduced. In the presence of inbreeding depression the greatest PP gain was achieved at still higher levels of status effective number, i.e. when more gene diversity was conserved in the BP. Thus, the combination of precise selection and PAM resulted in close to maximal short-term PP gain, while conserving maximal gene diversity in the BP.Communicated by O. Savolainen     

Estimating levels of inbreeding using AFLP markers

In the absence of detailed pedigree records, researchers have attempted to estimate individuals' levels of inbreeding using molecular markers, generally making use of heterozygosity measures based on microsatellite markers. Here we report and validate a method for estimating an individual's inbreeding coefficient, f, using amplified fragment length polymorphism (AFLP) markers. We use simulations to confirm that our measure scales appropriately with f when allele frequencies can be estimated from a subset of outbred individuals. We also present an approach for obtaining satisfactory estimates even in the absence of an independent set of known outbred individuals from which to estimate allele frequencies. We then test our method against empirical data from 179 wild and captive-bred old-field mice, Peromyscus polionotus subgriseus, comprising pedigree-based estimates of f, along with genetic data from 94 AFLP markers and 12 microsatellites. Inbreeding estimates based on both AFLP and microsatellite markers were found to correlate strongly with pedigree-based inbreeding coefficients. Owing to their ease of amplification in any species, AFLP markers may prove to be a valuable new tool for estimating f in natural populations and for examining correlations between heterozygosity and fitness.