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.     

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.     

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.     

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.     

Effect of non-random mating on genomic and BLUP selection schemes

Background

The risk of long-term unequal contribution of mating pairs to the gene pool is that deleterious recessive genes can be expressed. Such consequences could be alleviated by appropriately designing and optimizing breeding schemes i.e. by improving selection and mating procedures.

Methods

We studied the effect of mating designs, random, minimum coancestry and minimum covariance of ancestral contributions on rate of inbreeding and genetic gain for schemes with different information sources, i.e. sib test or own performance records, different genetic evaluation methods, i.e. BLUP or genomic selection, and different family structures, i.e. factorial or pair-wise.

Results

Results showed that substantial differences in rates of inbreeding due to mating design were present under schemes with a pair-wise family structure, for which minimum coancestry turned out to be more effective to generate lower rates of inbreeding. Specifically, substantial reductions in rates of inbreeding were observed in schemes using sib test records and BLUP evaluation. However, with a factorial family structure, differences in rates of inbreeding due mating designs were minor. Moreover, non-random mating had only a small effect in breeding schemes that used genomic evaluation, regardless of the information source.

Conclusions

It was concluded that minimum coancestry remains an efficient mating design when BLUP is used for genetic evaluation or when the size of the population is small, whereas the effect of non-random mating is smaller in schemes using genomic evaluation.     

Optimization of selection contribution and mate allocations in monoecious tree breeding populations

Background

The combination of optimized contribution dynamic selection and various mating schemes was investigated over seven generations for a typical tree breeding scenario. The allocation of mates was optimized using a simulated annealing algorithm for various object functions including random mating (RM), positive assortative mating (PAM) and minimization of pair-wise coancestry between mates (MCM) all combined with minimization of variance in family size and coancestry. The present study considered two levels of heritability (0.05 and 0.25), two restrictions on relatedness (group coancestry; 1 and 2%) and two maximum permissible numbers of crosses in each generation (100 and 400). The infinitesimal genetic model was used to simulate the genetic architecture of the trait that was the subject of selection. A framework of the long term genetic contribution of ancestors was used to examine the impacts of the mating schemes on population parameters.

Results

MCM schemes produced on average, an increased rate of genetic gain in the breeding population, although the difference between schemes was small but significant after seven generations (up to 7.1% more than obtained with RM). In addition, MCM reduced the level of inbreeding by as much as 37% compared with RM, although the rate of inbreeding was similar after three generations of selection. PAM schemes yielded levels of genetic gain similar to those produced by RM, but the increase in the level of inbreeding was substantial (up to 43%).

Conclusion

The main reason why MCM schemes yielded higher genetic gains was the improvement in managing the long term genetic contribution of founders in the population; this was achieved by connecting unrelated families. In addition, the accumulation of inbreeding was reduced by MCM schemes since the variance in long term genetic contributions of founders was smaller than in the other schemes. Consequently, by combining an MCM scheme with an algorithm that optimizes contributions of the selected individuals, a higher long term response is obtained while reducing the risk within the breeding program.     

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.     

Using genomic tools to maintain diversity and fitness in conservation programmes

Conservation programmes aim at maximizing the survival probability of populations, by minimizing the loss of genetic diversity, which allows populations to adapt to changes, and controlling inbreeding increases. The best known strategy to achieve these goals is optimizing the contributions of the parents to minimize global coancestry in their offspring. Results on neutral scenarios showed that management based on molecular coancestry could maintain more diversity than management based on genealogical coancestry when a large number of markers were available. However, if the population has deleterious mutations, managing using optimal contributions can lead to a decrease in fitness, especially using molecular coancestry, because both beneficial and harmful alleles are maintained, compromising the long‐term viability of the population. We introduce here two strategies to avoid this problem: The first one uses molecular coancestry calculated removing markers with low minor allele frequencies, as they could be linked to selected loci. The second one uses a coancestry based on segments of identity by descent, which measures the proportion of genome segments shared by two individuals because of a common ancestor. We compare these strategies under two contrasting mutational models of fitness effects, one assuming many mutations of small effect and another with few mutations of large effect. Using markers at intermediate frequencies maintains a larger fitness than using all markers, but leads to maintaining less diversity. Using the segment‐based coancestry provides a compromise solution between maintaining diversity and fitness, especially when the population has some inbreeding load.     

Optimisation of contribution of candidate parents to maximise genetic gain and restricting inbreeding using semidefinite programming (Open Access publication)

An approach for optimising genetic contributions of candidates to control inbreeding in the offspring generation using semidefinite programming (SDP) was proposed. Formulations were done for maximising genetic gain while restricting inbreeding to a preset value and for minimising inbreeding without regard of gain. Adaptations to account for candidates with fixed contributions were also shown. Using small but traceable numerical examples, the SDP method was compared with an alternative based upon Lagrangian multipliers (RSRO). The SDP method always found the optimum solution that maximises genetic gain at any level of restriction imposed on inbreeding, unlike RSRO which failed to do so in several situations. For these situations, the expected gains from the solution obtained with RSRO were between 1.5–9% lower than those expected from the optimum solution found with SDP with assigned contributions varying widely. In conclusion SDP is a reliable and flexible method for solving contribution problems.     

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.     

Minimization of rate of inbreeding for small populations with overlapping generations

We propose a method that minimizes the rate of inbreeding (delta F) for small unselected populations with overlapping generations and several reproductive age classes. It minimizes the increase in coancestry of parents and optimizes the contribution of each selection candidate. The carrying capacity of the population is limited to a fixed number of animals per year. When survival rate equalled 100%, only animals from the oldest age class were selected, which maximized the number of parents per generation, slowed down the turnover of generations and minimized the increase of coancestry across sublines. However, the population became split into sublines separated by age classes, which substantially increased inbreeding within sublines. Sublines were prevented by a restriction of selecting at least one sire and one dam from the second-oldest age class, which resulted in an L times lower delta F, where L equals the average generation interval of sires and dams. Minimum coancestry mating resulted in lower levels of inbreeding than random mating, but delta F was approximately the same. For schemes where the oldest animals were selected, delta F increased by 18-52% compared with the proposed method.     

Effect of linkage on the control of inbreeding in selection programmes

Selection and mating methods for controlling inbreeding in selection programmes are based on relationships obtained from pedigrees. The efficiency of these methods has always been tested by studies using genetic models of independent loci. However, under linkage the rate of inbreeding obtained from pedigrees can be different from the probability of identity by descent of genes. We simulated a quantitative trait under artificial selection controlled by a large number of genes spread on genome regions of different sizes. A method to control inbreeding based on minimising the average coancestry of selected individuals with a restriction in the loss of selection response, and a mating procedure to control inbreeding were applied. These methods, that use coancestry relationships, were not effective in controlling inbreeding when the genome sizes were smaller than five morgans or so. However, for larger genome sizes the methods were sufficiently efficient. For very tight linkage, methods that utilise molecular information from markers should be used. We finally discuss the effects of the selection of individual major genes on the neutral variability of adjacent genome regions.     

Optimum contribution selection in large general tree breeding populations with an application to Scots pine

Development of selection methods that optimises selection differential subject to a constraint on the increase of inbreeding (or coancestry) in a population is an important part of breeding programmes. One such method that has received much attention in animal breeding is the optimum contribution (OC) dynamic selection method. We implemented the OC algorithm and applied it to a diallel progeny trial of Pinus sylvestris L. (Scots pine) focussing on two traits (total tree height and stem diameter). The OC method resulted in a higher increase in genetic gain (8–30%) compared to the genetic gain achieved using standard restricted selection method at the same level of coancestry constraint. Genetic merit obtained at two different levels of restriction on coancestry showed that the benefit of OC was highest when restriction was strict. At the same level of genetic merit, OC decreased coancestry with 56 and 39% for diameter and height, respectively, compared to the level of coancestry obtained using unrestricted truncation selection. Inclusion of a dominance term in the statistical model resulted in changes in contribution rank of trees with 7 and 13% for diameter and height, respectively, compared to results achieved by using a pure additive model. However, the genetic gain was higher for the pure additive model than for the model including dominance for both traits.     

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.     

Population structure of the Trakehner Horse breed

The objective of this study was to examine the population structure of the Trakehner Horse breed. A total of 13 793 pedigree records were used for analysing the active breeding population and their ancestors dating back to 1950. Ancestors that were born before 1950 were called as base animals. The average generation interval was calculated as 10.2 years. The effective population size (Ne) was estimated by the increase in average year-wise inbreeding coefficient and average coancestry, respectively. Two methods were applied to estimate the effective population size: 1. Numerator-relationship-matrix (NRM), which did not consider missing ancestries. 2. Uncertain-parentage-matrix (UPM), which considered a probabilistic correction for unknown ancestors. There were no major differences between these two methods with respect to the rate of increase in inbreeding although the global levels using the UPM method were observed to be higher. Estimates for the inbreeding coefficients and the average coancestries varied little between both methods. The estimates of the effective population size per generation based on the rate of inbreeding ranged from 169 (NRM) to 150 (UPM) and 158 (NRM) to 144 (UPM) calculated by the average coancestry. From the early 1990s onwards, a strong increase in the rate of inbreeding was observed. This may be due to an increasing variance of the family size of sires and may be interpreted as a consequence of the growing use of artificial insemination. Analysing coancestries within and between the centrally managed regional breeding societies in Germany further revealed the Trakehner horse breed to be a genetically fragmented population with a main partition corresponding to formerly divided East and West Germany. The average rate of gene contributions (Thoroughbred (xx), Arab Horse breed (ox)) to the defined actual breeding population was calculated to be 22.3% xx-genes and 11.7% ox-genes.     

sofsog: a suite of programs to avoid inbreeding in plantation designs

Cost-effective ways of controlling inbreeding in conservation or productive plantations imply the allocation of individuals reducing the possibility of close relatives' mating and, consequently, limiting inbreeding. sofsog is a suite of programs, which helps to design plantation sites. First, if the plantation scheme involves several plots, it allows distribution of individuals available among different sites minimizing within-site global coancestry. Then, it yields a plantation design for each site, either following the classical permutated neighbourhood strategy or the recently developed method by Fernández and González-Martínez. This new method allows the implementation of different pollen dispersion kernels, and to include in the designing strategy any available information on individual relationships, reproductive success, differences in phenology, etc., via weighting or penalization matrices. Additionally, the package includes a tool for calculating the molecular coancestry (Identity By State) from codominant marker data.     

Controlling genetic variability by mathematical programming in a selection scheme on an open-pollinated population in Eucalyptus globulus

A model using integer quadratic mathematical programming has been developed to control the inbreeding level (or genetic diversity) through group coancestry in a selection programme for a forestry population structured in terms of maternal families coming from different locations. A method to calculate the average group coancestry between- and within-families for these open-pollinated populations is also proposed. This model has been applied to data from a breeding programme of Australian Eucalyptus globulus. The strategy proved to be effective as reductions of up to 50% for the group coancestry of the selected individuals were reached with a loss of only 5% of the maximum attainable selection differential (corresponding to truncation selection). Received: 14 October 1999 / Accepted: 26 July 2000     

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.     

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.     

Methods for minimizing the loss of genetic diversity in conserved populations with overlapping generations

Minimization of the average coancestry in a population has been theoretically proven to be the most efficient method to preserve genetic diversity. In the present study, based on a population genetic model, two methods to minimize the average coancestry in populations with overlapping generations were developed. For a given parental coancestry structure, the first method (OG) minimizes the average coancestry in the next generation, and the second method (LT) is designed to minimize the long-term accumulation of coancestry. The efficiencies of the two methods were examined by stochastic simulation. Compared to random choice of parents, the annual effective population sizes under the two proposed methods increased 2–3 folds. The difference among the two methods was small in a population with short generation interval. For populations with long generation intervals, the OG method showed a slightly larger annual effective size in an initial few years. However, in the subsequent years, the LT method gave a 5–15% larger annual effective size than the OG method. From these results, it is suggested that the LT method would be preferred to the OG method in most practical situations.