Document Detail

Optimum contribution selection for conserved populations with historic migration.
Jump to Full Text
MedLine Citation:
PMID:  23153196     Owner:  NLM     Status:  MEDLINE    
Abstract/OtherAbstract:
BACKGROUND: In recent decades, local varieties of domesticated animal species have been frequently crossed with economically superior breeds which has resulted in considerable genetic contributions from migrants. Optimum contribution selection by maximizing gene diversity while constraining breeding values of the offspring or vice versa could eventually lead to the extinction of local breeds with historic migration because maximization of gene diversity or breeding values would be achieved by maximization of migrant contributions. Therefore, other objective functions are needed for these breeds.
RESULTS: Different objective functions and side constraints were compared with respect to their ability to reduce migrant contributions, to increase the genome equivalents originating from native founders, and to conserve gene diversity. Additionally, a new method for monitoring the development of effective size for breeds with incomplete pedigree records was applied. Approaches were compared for Vorderwald cattle, Hinterwald cattle, and Limpurg cattle. Migrant contributions could be substantially decreased for these three breeds, but the potential to increase the native genome equivalents is limited.
CONCLUSIONS: The most promising approach was constraining migrant contributions while maximizing the conditional probability that two alleles randomly chosen from the offspring population are not identical by descent, given that both descend from native founders.
Authors:
Robin Wellmann; Sonja Hartwig; Jörn Bennewitz
Related Documents :
16202086 - Chloroplast and microsatellite dna diversities reveal the introduction history of brazi...
23572906 - Genetic diversity in indian common bean (phaseolus vulgaris l.) using random amplified ...
12180766 - A cautionary tale on ancient migration detection: mitochondrial dna variation in santa ...
17530696 - Morphological differentiation of aboriginal human populations from tierra del fuego (pa...
6334906 - Hla-dr2 and dr4 further defined by two new hla-d specificities (htc) derived from israe...
1634996 - Additive gene actions on the fiber number in the anterior optic tract of drosophila mel...
Publication Detail:
Type:  Journal Article     Date:  2012-11-15
Journal Detail:
Title:  Genetics, selection, evolution : GSE     Volume:  44     ISSN:  1297-9686     ISO Abbreviation:  Genet. Sel. Evol.     Publication Date:  2012  
Date Detail:
Created Date:  2012-12-12     Completed Date:  2013-04-11     Revised Date:  2013-11-04    
Medline Journal Info:
Nlm Unique ID:  9114088     Medline TA:  Genet Sel Evol     Country:  France    
Other Details:
Languages:  eng     Pagination:  34     Citation Subset:  IM    
Affiliation:
Institute of Animal Husbandry and Animal Breeding, University of Hohenheim, D-70599 Stuttgart, Germany. r.wellmann@uni-hohenheim.de
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:
Alleles
Animal Migration*
Animals
Animals, Domestic / genetics
Breeding
Founder Effect
Genetic Variation / genetics
Models, Genetic*
Pedigree
Population / genetics
Selection, Genetic*
Comments/Corrections

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): Genet Sel Evol
Journal ID (iso-abbrev): Genet. Sel. Evol
ISSN: 0999-193X
ISSN: 1297-9686
Publisher: BioMed Central
Article Information
Download PDF
Copyright ©2012 Wellmann et al.; licensee BioMed Central Ltd.
open-access:
Received Day: 11 Month: 6 Year: 2012
Accepted Day: 17 Month: 10 Year: 2012
collection publication date: Year: 2012
Electronic publication date: Day: 15 Month: 11 Year: 2012
Volume: 44 Issue: 1
First Page: 34 Last Page: 34
PubMed Id: 23153196
ID: 3807754
Publisher Id: 1297-9686-44-34
DOI: 10.1186/1297-9686-44-34

Optimum contribution selection for conserved populations with historic migration
Robin Wellmann1 Email: r.wellmann@uni-hohenheim.de
Sonja Hartwig1 Email: Sonja.Hartwig@uni-hohenheim.de
Jörn Bennewitz1 Email: j.bennewitz@uni-hohenheim.de
1Institute of Animal Husbandry and Animal Breeding, University of Hohenheim, D-70599 Stuttgart, Germany

Background

Many local varieties of domesticated animal species have been established in the last centuries. However, due to agricultural innovations since the beginning of the 19th century and subsequent intensification of production, many landraces are no longer adapted to their changing environments [1,2]. They have been crossed with superior breeds in order to improve the economic value of the breeding stock. Gene flow usually occured from the economically most important breeds to the landraces, but not backwards. Consequently, most historic breeds are now extinct and the remaining ones have considerable genetic contributions from a small number of economically superior breeds. Efforts are needed to prevent the remaining historic breeds and their gene pools to become extinct. Conservation efforts can have different objectives. Objectives of breeding programs can be to breed back the historic breeds by removing genetic contributions of migrants, to conserve the breeds in their present appearance, or to increase the economic values. In any case, genetic contributions arising from more frequent breeds are not subject to conservation efforts since their genes are widespread.

Meuwissen [3] proposed to maximize the expected mean breeding value of the offspring while constraining its gene diversity to a predefined value. A related but not equivalent approach is to maximize the gene diversity in the offspring with or without constraining its expected mean breeding value to a predefined value. In this paper, the latter approach is applied and generalized. This approach seems more appropriate for conserved populations because for these populations the focus is on conservation. In general, the method consists of calculating an optimum contribution ca (or the desired number of offspring) for each breeding individual a such that the offspring population maximizes an appropriate objective function ϕunder some side conditions.

In the classical approach [4] (Approach A) the gene diversity GD in the offspring O(c) is maximized, where the vector c contains the genetic contribution of each breeding individual to the offspring population. Thus, ϕA(c) = GD(O(c)). Gene diversity of a population is the probability that two alleles randomly chosen from the population are not identical by descent (IBD). However, this objective function may be not appropriate for conserved populations because maximization of gene diversity could be achieved by maximization of genetic contributions of migrants. Thus, this approach could eventually lead to extinction of the native breeds. Gene diversity should not fall below a certain level in order to avoid inbreeding depression. Gene diversity is, however, not the parameter that should be maximized in conserved populations. In conserved populations, we are interested in the conservation of alleles that come from native founders, as migrant alleles usually originate from non-endangered breeds. That is, we want to maximize the probability ϕB that both alleles are not IBD and descended from native founders (Approach B), or the probability ϕC that both alleles are not IBD and at least one of them descended from native founders (Approach C). We also considered the possibility of maximizing the conditional probability ϕD that both alleles are not IBD, given that both descended from native founders (Approach D). For Approach D, we constrained the mean migrant contribution in the offspring population.

Lacy [5] introduced the concept of founder genome equivalents (FGE). The FGE of a population is the minimum number of unrelated founders that would be needed to establish a population that has the same gene diversity as the population under study. Recall that gene diversity is the probability that two alleles chosen at random are not IBD. However, a more important parameter to characterize the value of a breed for conservation purposes is the conditional probability that two randomly chosen alleles are not IBD, given that both descended from native founders. We call it the conditional gene diversity of the population. Large conditional gene diversity indicates that many native founder alleles have been retained in the population even though they may be at low frequencies. This has led to the following definition of the native genome equivalents (NGE) of a population as the minimum number of unrelated founders that would be needed to establish a population that has the same conditional gene diversity as the population under study. It can be interpreted as the FGE that originate from native founders and that are still present in the population. Besides maintaining the economic value of the breed, the main objective of a conservation program for a population with historic migration is to maximize the NGE and to minimize the genetic contributions of migrants simultaneously.

In this paper, we compare objective functions ϕA,ϕB, ϕC and ϕD with respect to their ability to conserve the gene diversity, to increase the FGE originating from native founders (i.e. the NGE), and to decrease the genetic contributions of migrants. Algorithms for solving these optimization problems are also derived and implemented in the R package PedAnalysis. Methods were applied and effective sizes were calculated for three German cattle breeds: Vorderwald, Hinterwald and Limpurg.


Methods
Definitions

Since the methods were applied to populations with overlapping generations, all definitions are based on birth cohorts rather than generations. A birth cohort J is a set of individuals born in a particular time interval, e.g. the individuals Bt born in year t, or the population Pt at time t[6]. Since the date of death is unknown in most cases, the population Pt consists of all individuals up to a particular age T. This age T could be the average age of individuals when their last offspring was born, or, for simplicity, it could be the generation interval I. Thus, population Pt consists of all individuals born in the time interval t-T,t.

The gene diversity GD(J) of birth cohort J is the probability that two alleles chosen at random from the birth cohort are not IBD. We can write

[Formula ID: bmcM1]
(1) 
GD(J)=P(XJ≠YJ),

where alleles XJ and YJ are randomly chosen with replacement from birth cohort J, and founder alleles are assumed to be pairwise different. An equivalent representation is GD(J)=1−f¯J, where f¯J is the average coancestry in birth cohort J.

Each allele descends from a particular founder. Take to be the set of founder alleles. We distinguish between native founders and migrants, whereby a native founder is a founder that is not a migrant. A native founder is typically an individual with unknown pedigree that belongs to the population and was born before a certain date ts. A migrant is typically an individual that either comes from an other population (breed), or an individual with unknown pedigree that was born after the date. The date tscould be chosen shortly after establishment of the stud book when a sufficient portion of the population was recorded. We can write

[Formula ID: bmcM2]
(2) 
A=F∪ℳ,

where is the set of alleles that come from native founders and > is the set of alleles that come from migrants.

We define the conditional gene diversitycondGD(J) of birth cohort J as the conditional probability that two alleles randomly chosen from the birth cohort are not IBD, given that both descend from native founders. That is,

[Formula ID: bmcM3]
(3) 
condGD(J)=P(XJ≠YJ|XJ∈F,YJ∈F).

The founder genome equivalentsFGE(J) of birth cohort J is defined as the minimum number of founders that would be needed to establish a population that has the same gene diversity as the individuals in birth cohort J. It can be computed as

[Formula ID: bmcM4]
(4) 
FGE(J)=12f¯J=12(1−GD(J)),

see [4]. Analogously, we define the native genome equivalentsNGE(J) of birth cohort J as the minimum number of founders that would be needed to create a population that has the same conditional gene diversity as the individuals in birth cohort J. We have

[Formula ID: bmcM5]
(5) 
NGE(J)=12(1−condGD(J)).

However, a problem with this definition is that native founders of the population are assumed to be unrelated, which is not true. As a consequence, in the first generation the NGE would be almost as large as the total population size. However, due to the invalid assumption of unrelated founders, the limited effective size causes the NGE to decrease tremendously shortly after the last native founders have entered the population. In order to avoid this artifact, we extrapolate the history of the breed back in time and use as the reference population not the founders listed in the stud book, but the population at an earlier time t0. That is, all individuals are assumed to be unrelated in year t0. In the applications, the base year was t0 = 1800. We define the conditional gene diversity of an age cohort Jtat time t ≥ ts with respect to base year t0 as

[Formula ID: bmcM6]
(6) 
condGDt0(Jt)=1−12histNets−t0IcondGD(Jt)condGD(Pts),

where Pts is the population at time ts, I is the generation interval, and histNe is the historic effective size of the population. The historic effective size can be estimated from marker data [7]. The term that defines the conditional gene diversity is the product of two factors. The first is the estimated gene diversity in the population at time ts, and the second is the factor by which the conditional gene diversity decreased between ts and t. Consequently, the NGE with respect to base yeart0 can be calculated as

[Formula ID: bmcM7]
(7) 
NGEt0(Jt)=12(1−condGDt0(Jt)).

A further parameter that can be of interest is the effective size of the population. The effective size Ne(t1t2) of a population within a time interval t1t2 is the size of an idealized random mating population of constant size that causes the same decrease of gene diversity as the true population within t2−t1I generations. However, in breeds with steady gene flow from other populations, the gene diversity does not decrease below a certain level, so this definition of the effective size does not make much sense for populations with migration. Therefore, we use a slightly different definition. We define the native effective sizeNeN(t1t2) as the size of an idealized random mating population of constant size that causes the same decrease of the conditional gene diversity condGD(Pt) as the true population within t2−t1I generations. The effective population size at time t, defined as NeN(t) = limε→0NeN([t − ε, t + ε]), was calculated as described in [8], except that it was calculated from the conditional gene diversity. The native effective size quantifies the decrease of genome equivalents originating from native founders because the NGE depend only on the conditional gene diversity, as can be seen from the previous two equations. In a population without migration, Ne and NeNare equal. However, in a population with steady gene flow from other populations, NeN is smaller than Ne because the gene diversity approaches a plateau level, so Ne(t) goes to infinity.

The population Ptat time t, which consists of all individuals up to an age of T years, has gene diversity GD(Pt), native genome equivalents NGEt0(Pt), and genetic contribution CF(Pt) from native founders. Note that CF(J)=P(XJ∈F), so CF(J) is the probability that a randomly chosen allele from age cohort J descends from a native founder. Besides monitoring of these quantities, a major task for a conservation program is the calculation of optimal genetic contributions for the breeding individuals that maximize the conditional gene diversity in the offspring and simultaneously maximize the genetic contribution from native founders in the offspring. Moreover, a sufficient level of gene diversity must be maintained in order to avoid inbreeding depression. In general, however, the quantities NGEt0(J) and CF(J) cannot be maximized simultaneously, so an objective function is needed that considers each appropriately.

The usual approach (Approach A) for populations without migration is the calculation of genetic contributions ctA for the breeding individuals of population Ptsuch that the gene diversity

[Formula ID: bmcM8]
(8) 
ϕA(J)=GD(J)=P(XJ≠YJ)

is maximized by a hypothetical (infinitely large) offspring population Ot(ctA). This approach is called minimum kinship selection [9]. Note that the gene diversity GD(Ot(ctA))=ϕA(Ot(ctA)) of the hypothetical offspring is known as the potential diversity of the population at time t[6]. A more appealing approach for populations with migration is to use genetic contributions ctB for the breeding individuals such that the probability

[Formula ID: bmcM9]
(9) 
ϕB(J)=P(XJ≠YJandXJ∈FandYJ∈F)

is maximized by the resulting offspring population Ot(ctB). This is the probability that two randomly chosen alleles from the offspring are not IBD and are from native founders (Approach B). As a third approach, we consider maximization of the probability that two randomly chosen alleles from the offspring are not IBD and at least one of them descends from a native founder (Approach C). In this case, genetic contributions ctC for the breeding individuals are calculated such that the offspring population Ot(ctC) maximizes

[Formula ID: bmcM10]
(10) 
ϕC(J)=P(XJ≠YJand(XJ∈ForYJ∈F)).

Finally, we consider maximizing the conditional gene diversity in the offspring population. That is, genetic contributions ctD for the breeding individuals were calculated such that the conditional probability

[Formula ID: bmcM11]
(11) 
ϕD(J)=P(XJ≠YJ|XJ∈FandYJ∈F)

is maximized. This approach is intuitively appealing because it maximizes NGE. It has, however, the disadvantage that the conditional gene diversity can be large even for offspring populations with very large migrant contributions. This is due to conditioning on the event that the randomly chosen alleles XJ and YJoriginate from native founders. This can be seen as follows. Take a solution ctD of the optimization problem and suppose that at least one migrant is a potential breeding individual. Then it can be shown mathematically that the genetic contribution of this migrant to the offspring population can be arbitrarily increased without changing the value of the objective function. Thus, the solution of the optimization problem may be not unique, and one solution maximizes migrant contributions. In order to avoid this, we put an additional constraint on the maximum permissible value for the genetic contribution from migrants to the offspring population.

Computations

To calculate the parameters defined in the previous section, the following quantities are needed. First, the coancestry fi,j is needed for each pair of individuals i,j. It is the probability that two alleles randomly chosen from the individuals are IBD. That is,

[Formula ID: bmcM12]
(12) 
fi,j=P(Xi=Xj),

where allele Xi is randomly chosen from the two alleles of individual i at a particular locus.

Now we define an equivalence relation on the set of founder alleles. Two alleles xi,xj are equivalent (xiMxj) if they are IBD or if both are migrant alleles. For two alleles randomly chosen from individuals i, j, the probability for this to occur is

[Formula ID: bmcM13]
(13) 
fi,jM=P(Xi≡MXj)=PXi=Xjor(Xi∈ℳ,Xj∈ℳ).

A second equivalence relation is defined as follows. Two alleles xi,xj are equivalent (xiFMxj) if both are native founder alleles or if both are migrant alleles. For two alleles randomly chosen from individuals i, j, the probability for this to occur is

[Formula ID: bmcM14]
(14) 
fi,jFM=P(Xi≡FMXj)=P(Xi∈F,Xj∈F)or(Xi∈ℳ,Xj∈ℳ).

These probabilities have the advantage that they can easily be computed with existing software, e.g. with function kinship() from the R-package kinship. For calculation of fi,jM, the parents of all migrants were identified with the same dummy individual and for this individual a pedigree with several generations of selfing was added. The coancestry of individuals i, j, computed from this extended pedigree is equal to fi,jM. Equality holds only approximately because only a finite number of generations of selfing was added. For calculation of fi,jFM, the parents of all migrants were identified with one single dummy individual, the parents of all native founders were identified with another single dummy individual, and for both individuals pedigrees with several generations of selfing were added. The coancestry of individuals i, j, computed from this extended pedigree, is equal to fi,jFM. For example, consider two full sibs i, j whose sire is a migrant and whose dam is a native founder. Their coancestry is fij=14, but fijM=38, and fijFM=12.

Let fPt be the Nt×Nt coancestry submatrix for the Nt individuals from population Pt that is obtained from the true pedigree (i.e., fPt=(fij)i,j∈Pt). The Nt×Nt matrix that contains the probabilities fi,jM for each pair of individuals i, j from population Ptis denoted as fPtM=(fijM)i,j∈Pt, and the Nt×Nt matrix that contains the probabilities fi,jFM is denoted as fPtFM=(fijFM)i,j∈Pt. That is, rows and columns that correspond to individuals not born in time interval [t-T,t] and dummy individuals were excluded from the matrix.

Additionally, the Nt-dimensional vector Ct=(Ct1,…,CtNt)T is needed and contains the genetic contribution of native founders for each individual of population Pt. Note that CF(Pt) is the mean of vector Ct. Let f¯Pt,f¯PtM, and f¯PtM be the means of the respective matrices. It is well known that the gene diversity can be computed as [4]

[Formula ID: bmcM15]
(15) 
GD(Pt)=1−f¯Pt.

Proofs of all numbered equations are presented in Additional file 1, in which it is shown that the conditional gene diversity satisfies

[Formula ID: bmcM16]
(16) 
condGD(Pt)=f¯PtFM−f¯PtMCF(Pt)−1−f¯PtFM2.

Let OtN(c) be an arbitrary (hypothetical) offspring population of size N that is obtained from population Pt such that each breeding individual a ∈  Pt has genetic contribution cato the offspring population. The probability that an allele randomly chosen from the offspring population descends from a native founder is

[Formula ID: bmcM17]
(17) 
CF(Ot(c))=CF(OtN(c))=cTCt,

and the conditional gene diversity in the offspring population is

[Formula ID: bmcM18]
(18) 
condGD(Ot(c))=limN→∞condGD(OtN(c))=cT(fPtFM−fPtM)ccTCt−1−cTfPtFMc2.

It is well known that

[Formula ID: bmcM19]
(19) 
limN→∞ϕA(OtN(c))=1−cTfPtc,

so the optimum contributions ctA for the breeding individuals with respect to objective function ϕAminimize cTfPtc under side conditions ca ≥ 0 and ∑aca=1. Additional side conditions can be added to fulfil biological and practical requirements. Moreover, we have

[Formula ID: bmcM20]
(20) 
limN→∞ϕB(OtN(c))=cT(fPtFM−fPtM)c,

so the optimum contributions ctB for the breeding individuals with respect to objective function ϕB minimize cT(11T−(fPtFM−fPtM))c under the side conditions described above, where 1 is a vector with ones. Since

[Formula ID: bmcM21]
(21) 
limN→∞ϕC(OtN(c))=1−cTfPtMc,

the optimum contributions ctC for the breeding individuals with respect to objective function ϕC minimize cTfPtMc under the side conditions. Finally, we have

[Formula ID: bmcM22]
(22) 
limN→∞ϕD(OtN(c))=cT(fPtFM−fPtM)ccTQtc,

where Qt=12Ct1T+1CtT−11T+fPtFM is a Nt × Nt matrix. This function was maximized under the side conditions described above. Moreover, the additional side constraint cTCt≥cF was applied, where cF is the minimum permissible contribution of native founders to the offspring population. This is a quadratic fractional programming problem with linear constraints, so the objective function could have multiple local maxima. As mentioned in the previous section, one solution of the optimization problem maximizes migrant contributions, so the inequality constraint could be replaced by the equality constraint cTCt=cF. For each offspring population J that satisfies this equality constraint, the objective function (i.e. the conditional gene diversity) satisfies

[Formula ID: bmcM23]
(23) 
ϕD(J)≈ϕB(J)P(XJ∈F)2=ϕB(J)cF2∝ϕB(J),

where the approximation is exact if the events XJ∈F and YJ∈F are independent. Therefore, an approximate solution was obtained by maximizing objective function ϕB under the additional constraint cTCt=cF. The resulting contributions for the breeding individuals were used as starting values for general nonlinear optimization in order to obtain the exact solution. In the applications, the threshold value cF was quite arbitrarily chosen as the 75% quantile of the genetic contributions from native founders to individuals in the population. The same quantile was used for all breeds and years in order to make the results comparable. Results could be improved by choosing breed dependent threshold values.

We used the interior point method ipop in R-package kernlab (see [10]) for objective functions ϕB and ϕD, whereas for objective functions ϕAand ϕCwith positive definite matrices we used solve.QP from R-package quadprog. It implements the dual method of Goldfarb and Idnani [11,12].

Materials

Only three local cattle varieties of Baden and Württemberg in the south-west of Germany have been preserved from extinction. These are the Vorderwald cattle, Hinterwald cattle, and Limpurg cattle. Other local breeds were replaced by Simmentaler Fleckvieh after their introduction at the beginning of the 19th century because the small landraces were not suitable for tillage [1].

The small Hinterwald cattle could be preserved as an almost pure breed until the beginning of the 20th century [13,14] because the poor soil quality in its region of origin was not suitable for larger breeds. Nevertheless, this breed adopted the colour of the Simmentaler Fleckvieh during the 19th century [15]. The Hinterwald cattle were occasionally crossed with the Vorderwald cattle [16] and with Fleckvieh.

The red-and-white marked, colour-sided [17]Vorderwald cattle were frequently crossed with Simmentaler cattle. Consequently, the white stripe along the back became rare already around 1900 [16]. After the Second World War, Vorderwald cattle were also crossed with Ayrshire, Red Holstein and Montbéliard cattle in order to improve milk yield. These crosses were registered as Vorderwald cattle. Extinction probabilities for Vorderwald and Hinterwald cattle were estimated by [18].

The yellow coloured Limpurg cattle were not only frequently crossed with Simmentaler cattle [19], but also occasionally with Braunvieh and Gelbvieh cattle [15] in order to increase body size. Nevertheless, the population size decreased dramatically. Only 17 Limpurg cows were registered in 1967, so the breeding association was dissolved. Several Limpurg cattle, however, were rediscovered in 1986 and a new stud book was established. Not only Limpurg cattle were registered, but also Fleckvieh crosses, and some Gelbvieh and Glan-Donnersberger bulls [16].

The data consisted of the pedigrees and additional information on 25 412 Hinterwald cattle, 185 315 Vorderwald cattle, and 4 150 Limpurg cattle. Vorderwald cattle without offspring were removed from the data in order to reduce the data set. Pedigrees of Hinterwald and Vorderwald cattle trace back only to 1948 because the stud books were renewed after the Second World War. Pedigrees of Limpurg cattle trace back only to 1970. Cattle from other breeds were considered to be migrants. Additionally, Hinterwald and Vorderwald cattle with unknown pedigree born after ts = 1970 were also considered migrants, although some may have purebred ancestors. Limpurg cattle with unknown pedigree were considered to be migrants if they were born after ts = 1988. The generation intervals were similar for the three breeds (unpublished results). Here, we assumed a generation interval of I = 5.3 years for all breeds.


Results

The left hand side of Figure 1 shows the development of the native effective size NeN for the three breeds. Around 1990, the effective size of Limpurg cattle was only about 20, which was due to the small population size. However in most cases, the effective size was above 50 for all three breeds. In 2011, 7952 Vorderwald cows, 2328 Hinterwald cows, and 471 Limpurg cows were registered. Interestingly, there appeared to be no relationship between the effective size and the total population size when the number of individuals exceeds the minimum number required to reach an Neof approximately 50.

The right hand side of Figure 1 shows for each breed how the genetic contributions of migrants changed over time. Migrant contributions are shown for the true population Pt and for the hypothetical offspring populations that would be obtained if optimum contribution selection were applied to population Pt. The solid lines show that migrant contributions increased steadily for all three breeds. The dashed line for offspring A shows that all three breeds would become extinct if optimum contribution selection were used to maximize the gene diversity in the offspring. In contrast, objective functions ϕB and ϕC would reduce migrant contributions substantially by more than 50% in all three breeds. According to the constraint applied for objective function ϕD, the corresponding line shows the 25% quantile of the migrant contributions in the population.

The left hand side of Figure 2 shows the development of NGE for the true population and for the hypothetical offspring populations. We used the year t0 = 1800 as the base year. The historic Ne is not known for these breeds. In the figure, we assumed a historic Neof 150 for each breed, which is in good accordance with the results obtained by [20] for various cattle breeds during this period of time. For the population in 2005, the computed NGE with respect to base year t0 was 3.1 for Limpurg cattle, 3.3 for Hinterwald cattle, and 3.2 for Vorderwald cattle. For comparison, the NGE computed under the assumption of unrelated founders was 7.3 for Limpurg cattle, 8.3 for Hinterwald cattle, and 8.0 for Vorderwald cattle. Figure 2 shows that the NGE would decrease by using objective function ϕB. This suggests that the individuals with the smallest migrant contributions are closely related, so they share the same founder alleles. The use of objective function ϕCwould cause a small increase of the NGE for all three breeds. If the constraint on migrant contributions is not too serious, then objective function ϕDwould cause the largest increase in NGE. However, the potential to increase NGE is limited.

The right hand side of Figure 2 shows the changes in gene diversity. It can be seen that the gene diversity is high for all three breeds. This is caused by migration. Note that the native effective population size quantifies the decrease of genome equivalents arising from native founders, so the gene diversity can be constant (or increase due to migration) even if the native effective population size is small. As expected, optimum contribution selection with objective functions ϕB, ϕC, or ϕD would cause a moderate but an acceptable loss of gene diversity.


Discussion

Most of the time, the native effective size NeNwas above 50 for the three breeds and due to migration, Ne was larger than NeN. An effective size of at least 50 is considered acceptable, although an Neof 100 is recommended to be on the safe side [21]. Many cattle breeds have effective sizes between 50 and 100 regardless of the total population size. Therefore, in order to conserve the overall gene diversity, it is generally recommended to conserve a large number of breeds with small population sizes rather than a small number of breeds with large population sizes. In this case, different alleles would be preserved in different subpopulations. These populations can be used as resources to identify advantageous genes that can be introgressed into commercial populations. Conserved populations must be sufficiently large to allow for this. However, breeds that are close to the economic viability threshold and populations that are expected to occupy niches that are different from that of established commercial breeds, should have larger population sizes in order to enable a sufficient selection response. Examples of the importance of farm animal genetic resources are the introgression of the polled gene into economically important cattle breeds, the introduction of indicine cattle breeds to South America because of their adaption to extreme environments, and introgression of genes for disease resistance into highly productive susceptible breeds [22].

The current NeNof the Vorderwald cattle was smaller than the estimates of the effective size obtained by [23] with other methods. The reason is probably that other methods do not distinguish between migrants and native founders. Genome equivalents arising from native founders are likely to decline faster than those arising from migrants because migrants are usually from economically superior breeds. The sufficiently large NeN show that for all three breeds, migration from other breeds was much larger than it was needed to avoid unacceptably high inbreeding depression. As a consequence, these breeds share only a small portion of their genes with the corresponding historic breeds of the same name. We showed that it is still possible to substantially increase the genetic contribution from the historic breeds by optimum contribution selection.

For optimum contribution selection, the choice of the objective function was crucial. Maximization of gene diversity (Approach A) turned out to substantially increase the migrant contributions and thus would lead to the extinction of these breeds. Approach B has the desired effect to substantially decrease the migrant contributions but does not put enough weight on the conservation of gene diversity. It is not recommended because it would reduce NGE and cause the largest loss of gene diversity. Approach C is recommended for conserved populations because for all three breeds the use of this objective function substantially decreased the migrant contributions, increased the NGE, and caused only a moderate decrease of gene diversity. Approach D can also be recommended, although it requires choice of a threshold for the migrant contributions. If the threshold is chosen appropriately, then this approach causes the largest increase in NGE. However, the potential to increase the NGE was small for the breeds considered. Interestingly, for the current populations, optimum contributions for Approach A were slightly negatively correlated with the optimum contributions obtained for the other approaches, whereas the optimum contributions for the approaches B, C, and D were pairwise positively correlated (not shown).

Amador et al.[24] proposed two other approaches to reduce migrant contributions. Their first approach was to minimize migrant contributions in the offspring population. Their second approach was to minimize the probability that two alleles randomly chosen from the offspring population are IBD and descend from migrants. This objective function was computed from partial coancestry coefficients [25], but could also be computed by the methodology introduced in this paper. For both approaches, the maximum rate of inbreeding was restricted. However, provided that an acceptable rate of inbreeding can be achieved, it is not obvious why it is desirable that alleles originating from migrants should be not IBD in the offspring population. In contrast, all approaches proposed in this paper aim at increasing the probability that alleles originating from native founders are not IBD. Amador et al.[24] concluded that even with only a few generations without management, a small amount of introgression can spread into the population and it may be almost impossible to recover. This was not observed in our study. The reason is probably that the total population sizes of the cattle breeds were much larger than their effective sizes, which increased the probability to find individuals with small migrant contributions. Moreover, the cattle populations may deviate from random mating populations because some breeders avoid the use of bulls with high migrant contributions.

Another approach could be to minimize the effective number of non-founders Nenf, as defined by Caballero and Toro [4], in the offspring population. This approach would be equivalent to maximization of f¯O−12Nef(O) which would be achieved by increasing the average relationship f¯O in the offspring population O and by increasing the effective number Nef of founders in the offspring generation. Thus, the rate of inbreeding would have to be restricted by this alternative approach. This approach, however, would by definition not be optimal with respect to the objective functions introduced in this paper.

Our results show that migrant contributions can be substantially decreased for all three breeds, but the potential to increase the NGE is limited. The reduction of migrant contributions would be largely achieved in the first generation of management. In subsequent generations, some further improvement would be possible due to biological restrictions in previous generations. However, thereafter the management method becomes equivalent to an equalization of family sizes and no further reduction of migrant contributions could be achieved. Moreover, pedigree-based optimum contribution selection cannot remove genetic contributions of migrants that arose before recording of pedigrees started. However, removal of migrant contributions that arose earlier can be done subsequent to pedigree-based optimum contribution selection by identification of chromosome segments that are also present in the migrant breeds and by removal of those individuals with large migrant contributions from the breeding pool. Since migrants are usually males, haplotype variants of the Y-chromosome can be used as markers for paternal lineage [26] to identify the migrant breeds. For individuals that are not removed from the breeding pool (i.e individuals with small migrant contributions), optimum contributions can be calculated based on genomic relationships. In order to avoid that this approach causes the frequencies of migrant alleles to increase, the set of breeding individuals could be enlarged with individuals of the migrant breeds. After the optimum contributions have been computed, the contributions of these additional migrant individuals are set to zero, and the optimum contributions for individuals of the breed of interest are rescaled, so that they add up to one. Thereafter, it would be beneficial to combine closely related breeds with low gene diversity in order to reduce extinction probabilities [27], and to split breeds with a high gene diversity into several subpopulations in order to reduce the decrease of overall gene diversity [28]. Breeds with highest value for conservation should be given priority [29]. These breeds are likely found near the domestication center (since genetic diversity declines with increasing distance from the domestication centre [30]), far from the native areas of economically superior breeds, or live in harsh environmental conditions. Candidates are also breeds that are used for uncommon purposes (e.g. fighting cattle, cattle breeds used for cow racing).


Conclusions

The usual recommendation to optimize contributions for breeding individuals by maximizing gene diversity in the offspring is not suitable for populations with historic migration because maximization of gene diversity would be achieved by maximization of migrant contributions. Thus, this approach, applied to populations with migration, would rapidely lead to their extinction. Two approaches can be recommended. The first is to maximize the probability that two alleles randomly chosen from the offspring population are not IBD and that at least one of them descended from a native founder (Approach C). The other approach is to constrain migrant contributions while maximizing the conditional probability that two alleles randomly chosen from the offspring population are not IBD, given that both descended from native founders (Approach D). Migrant contributions could be substantially decreased for the three breeds investigated here, but the potential to increase the NGE is limited.

Programs for pedigree-based optimum contribution selection and for the analyses presented in this paper are available in R package PedAnalysis from the first author. Since migrants are usually from genetically superior breeds, optimum contribution selection is likely to reduce breeding values if there is no constraint on the expected breeding value of the offspring. The program for optimum contribution selection allows adding the constraint that the expected mean breeding value of the offspring does not fall below a certain value. Moreover it is possible to put a constraint on the maximum number of offspring per male and female.


Competing interests

The authors declare that they have no competing interests.


Authors’ contributions

RW developed and implemented the methods and drafted the manuscript. SH prepared the data files and applied the methods. JB conceived the project and helped in drafting the manuscript. All authors have read and approved the final manuscript.


Supplementary Material Additional file 1

Proofs. The file provides the proofs of all numbered equations.


Click here for additional data file (1297-9686-44-34-S1.pdf)


Acknowledgements

The data were kindly provided by Henning Hamann, Landesamt für Geoinformation und Landentwicklung Baden-Württemberg, Germany. The authors thank the referees for pointing us to another paper dealing with this subject.


References
Brüggemann H,Viehhaltung im BauernbetriebYear: 1954DLG-Verlag: Frankfurt am Main
Oldenbroek K,Oldenbroek KIntroductionUtilisation and conservation of farm animal genetic resourcesYear: 2007Wageningen: Wageningen Academic Publishers1327
Meuwissen THE,Maximising the response of selection with a pre-defined rate of inbreedingJ Anim SciYear: 1997449349409110204
Caballero A,Toro MA,Interrelations between effective population size and other pedigree tools for the management of conserved populationsGenet ResYear: 20004433134310.1017/S001667239900444910893869
Lacy RC,Analysis of founder representation in pedigrees: founder equivalents and founder genome equivalentsZoo BiolYear: 19894411112310.1002/zoo.1430080203
Wellmann R,Pfeiffer I,Pedigree analysis for conservation of genetic diversity and purgingGenet ResYear: 20094420921910.1017/S0016672309000202
Meuwissen THE,Goddard ME,Multipoint Identity-by-Descent prediction using dense markers to map QTL and estimate effective population sizeGeneticsYear: 2007442551256010.1534/genetics.107.07095317565953
Wellmann R,Bennewitz J,Identification and characterization of hierarchical structures in dog breeding schemes, a novel method applied to the Norfolk TerrierJ Anim SciYear: 2011443846385810.2527/jas.2010-370921821815
Meuwissen THE,Oldenbroek KOperation of conservation schemesUtilisation and conservation of farm animal genetic resourcesYear: 2007Wageningen: Wageningen Academic Publishers167193
Vanderbei RJ,LOQO: An interior point code for quadratic programmingOptim Method SoftwYear: 19994445148410.1080/10556789908805759
Goldfarb D,Idnani A,Hennart JPDual and primal-dual methods for solving strictly convex quadratic programsNumerical AnalysisYear: 1982Berlin: Springer-Verlag226239
Goldfarb D,Idnani A,A numerically stable dual method for solving strictly convex quadratic programsMath ProgramYear: 19834413310.1007/BF02591962
Martin W,Zeeb R,Martin-Zeeb Handbuch der LandwirtschaftYear: 1922Stuttgart: Verlag Eugen Ulmer
Schlipf JA,Schlipf’s Handbuch der LandwirtschaftYear: 1914Berlin: Verlag Paul Parey
Hansen J,Lehrbuch der RinderzuchtYear: 1921Berlin: Verlag Paul Parey
Sambraus HH,Gefährdete NutztierrassenYear: 2010Stuttgart: Ulmer Eugen Verlag
Mangold P,Reicherter E,Neues illustriertes HaustierarzneibuchYear: 1925Reutlingen: Verlag Enklin und Laib
Bennewitz J,Meuwissen THE,Estimation of extinction probabilities of five German cattle breeds by population viability analysisJ Dairy SciYear: 2005442949296110.3168/jds.S0022-0302(05)72975-116027209
Wilckens M,Die Rinderrassen Mittel-EuropasYear: 1885Berlin: Verlag von Paul Parey
Villa-Angulo R,Matukumalli LK,Gill CA,Choi J,Van Tassell CP,Grefenstette JJ,High-resolution haplotype block structure in the cattle genomeBMC GenetYear: 2009441919393054
Meuwissen THE,Genetic management of small populations: A reviewActa Agr Scand AYear: 2009447179
Toro M,Mäki-Tanila A,Oldenbroek KGenomics reveals domestication history and facilitates breed developmentUtilisation and conservation of farm animal genetic resourcesYear: 2007Wageningen: Wageningen Academic Publishers1327
vit LfL ADR,Erfassungsprojekt Erhebung von Populationsdaten tiergenetischer Ressourcen in Deutschland: Tierart RindYear: 2010Verden/Grub/Bonn: Bundesministerium für Landwirtschaft und Ernährung
Amador C,Toro MA,Fernández J,Removing exogeneous information using pedigree dataConserv GenetYear: 2011441565157310.1007/s10592-011-0255-4
García-Cortés LA,Martínez-Ávila JC,Toro MA,Fine decomposition of the inbreeding and the coancestry coefficients by using the tabular methodConserv GenetYear: 2010441945195210.1007/s10592-010-0084-x
Svensson E,Götherström A,Temporal fluctuations of Y-chromosomal variation in Bos taurusBiol LettYear: 20084475275410.1098/rsbl.2008.034218713712
Bennewitz J,Simianer H,Meuwissen THE,Investigations on merging breeds in genetic conservation schemesJ Dairy SciYear: 2008442512251910.3168/jds.2007-092418487675
Kimura M,Crow JF,On the maximum avoidance of inbreedingGenet ResYear: 19634439941510.1017/S0016672300003797
Bennewitz J,Kantanen J,Tapio I,Li MH,Kalm E,Vilkki J,Ammosov I,Ivanova Z,Kiselyova T,Popov R,Meuwissen TH,Estimation of breed contributions to present and future genetic diversity of 44 North Eurasian cattle breeds using core set diversity measuresGenet Sel EvolYear: 20064420122010.1186/1297-9686-38-2-20116492375
Groeneveld LF,Lenstra JA,Eding H,Toro MA,Scherf B,Pilling D,Negrini R,Finlay EK,Jianlin H,Groeneveld E,Weigend S,GLOBALDIV Consortium: Genetic diversity in farm animals - a reviewAnim GenetYear: 20104463120500753

Figures

[Figure ID: F1]
Figure 1 

Native effective size and migrant contributions. Native effective size NeN(left) and genetic contributions from migrants (right) in the real population Ptand in the hypothetical offspring populations for selection strategies A, B, C, and D (right) for (a): Limpurg cattle, (b): Hinterwald cattle, (c): Vorderwald cattle.



[Figure ID: F2]
Figure 2 

Native genome equivalents and gene diversities. Genome equivalents originating from native founders NGE (left) and gene diversity (right) in the real population Ptand in the hypothetical offspring populations for selection strategies A, B, C, and D for (a): Limpurg cattle, (b): Hinterwald cattle, (c): Vorderwald cattle.



Article Categories:
  • Research


Previous Document:  Pro-death and pro-survival properties of ouabain in U937 lymphoma derived cells.
Next Document:  Detecting Supercritical CO2 in Brine at Sequestration Pressure with an Optical Fiber Sensor.