|
|
||||||||
1 Department of Genetics and Biotechnology, Danish Institute of Agricultural Sciences DK-8830 Tjele, Denmark
2 Department of Large Animal Sciences, Royal Veterinary and Agricultural University DK-1870 Frederiksberg C, Denmark
Corresponding author: Anders Christian Sørensen; e-mail: AndersC.Sorensen{at}agrsci.dk.
| ABSTRACT |
|---|
|
|
|---|
Key Words: inbreeding effective population size genetic diversity pedigree analysis
| INTRODUCTION |
|---|
|
|
|---|
The population structure determines the development in inbreeding. Danish Holstein has undergone repeated backcrossing to North American Holstein since 1965, so that now more than 93% of the Danish Holstein genes are of North American origin. Danish Holstein is the largest of the breeds with 300 young bulls progeny tested every year. Danish Jersey and Danish Red are smaller breeds with 60 young bulls each starting for progeny testing every year. Danish Jersey has imported US Jersey breeding material and has now around 35% US genes. Robertson and Mason (1954) predicted that Danish Red would meet inbreeding problems, as a few ancestors contributed very heavily (~18%) to the gene pool of breeding bulls, and Danish Red suffered from inbreeding depression in the 1970s. The breeders, therefore, chose to import genes from US Brown Swiss and later from Red Holstein. In the 1990s, an import from Swedish Red and White started. Due to this ongoing introduction of new genes, Danish Red is expected to show a different trend in inbreeding than the other 2 breeds. For all breeds, the importation of semen is ongoing, with approximately 7% of inseminations being performed with imported semen. The objective of this paper was to monitor current trends and predict future trends in inbreeding and to assess the genetic diversity in Danish Holstein, Danish Jersey, and Danish Red.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
![]() |
where Csire and Cdam are contributions from the paternal and maternal lines respectively:
![]() |
where ai is the proportion of known ancestors in generation i; and d is the number of generations that is taken into account. In this study, 5 generations are considered (d = 5). This index is ad hoc in the sense that a specific value cannot be translated into an expected bias in the calculated coefficient of inbreeding. However, being a harmonic mean, the index has a value of zero if one parent is unknown no matter how much pedigree is known for the other parent. There were no substantial differences between breeds in the completeness of pedigrees (Table 1
and Figure 2
). Hence, the results are comparable across breeds. More than half of the reference population has pedigree completeness for 5 generations above 0.9.
|
![]() |
where nj is the number of ancestors of individual j, and gij is the number of generations between individual j and its ancestor i (Sölkner et al., 1998). In that way, one-half is added for each known parent, one-fourth is added for each known grandparent, and so on. The complete generation equivalent quantifies how many generations have been traced. This number is around 7 for all 3 breeds (Table 1
).
The correlation between the pedigree completeness index and complete generation equivalents was between 0.94 and 0.96 for the 3 breeds, indicating that they assess more or less the same.
Methodology
Effective population size.
The effective population size, Ne, was calculated from the rate of inbreeding per generation, obtained by multiplying the annual rate of inbreeding,
Fy, with the generation interval, L,
![]() | ([1]) |
The effective population size was calculated for time intervals in which the trend in inbreeding was approximately linear. This was assessed by visual inspection of the trends in Figure 3a
.
|
fy, can be interpreted as the future rate of inbreeding per year. Therefore, the expected future effective population size can be calculated by replacing
Fy in [1] with
fy. The expected future effective population size was calculated for the same time intervals as the historical effective population size above assuming constant generation interval.
Effective number of founders.
All animals with both parents unknown are regarded as founders in this analysis. In addition, if an animal has one known and one unknown parent, the unknown parent is regarded as a founder. The total number of founders contains limited information on the genetic basis for the population. Firstly, founders are assumed to be unrelated, because their parentage is unknown. However, this is most probably not the case. Secondly, some founders have been used more intensely and therefore contribute more to the current population than other founders. The effective number of founders, fe, has been designed to correct for this second shortcoming. It is defined as the number of equally contributing founders that would be expected to generate a similar amount of genetic diversity as the studied population (Lacy, 1989). It can be calculated from the genetic contributions of the Nf founders:
![]() |
where qi is the genetic contribution of the ith founder to the reference population. When founders contribute unequally, the effective number of founders is smaller than the actual number. The effective number of founders has a limited usefulness; because once the contributions of the founders have converged, the effective number of founders remains constant, no matter what happens to the population afterwards. The genetic contributions are expected to have converged after 5 to 7 generations (Bijma and Woolliams, 1999). For example, a population where animals have deep pedigrees (more than 7 generations) can be characterized with a high effective number of founders even after a severe, recent bottleneck. Hence, in isolation, the effective number of founders is not a good measure of genetic diversity. However, it provides a basis for comparison of the effective population size and the effective number of ancestors (see below). The effective number of founders is expected to be one-half the effective population size in a population with minimum inbreeding. If the effective number of founders is different from this, it shows that the breeding structure has been changed since the founder generation.
Effective number of ancestors.
The effective number of ancestors, fa, has been defined to supplement the effective number of founders (Boichard et al., 1997). Instead of using genetic contributions of founders only, the effective number of ancestors is calculated from the genetic contributions of ancestors with the largest marginal genetic contributions. The genetic contributions of founders are independent and sum to one. That is not the case for genetic contributions of ancestors. For example, the dam of a highly used sire has at least half the contribution of her son, because the same genes are represented in both generations. Boichard et al. (1997) therefore introduced the marginal contribution. The ancestors contributing most to the reference population are considered one at a time in a recursive process. For each round of the recursion the ancestor with the highest contribution is chosen, and the contributions of all others are calculated conditional on the contribution of the chosen ancestor. Then, based on these marginal contributions another ancestor is chosen and the process continues. So the marginal contribution is the genetic contribution from an individual after correcting for contributions of other ancestors already considered in the recursive process. Thus, the sum of marginal contributions of all ancestors is one. Ancestors only have a large marginal contribution to the reference population if their genes have passed through many descendants, e.g., a sire of sons with many sons selected. The effective number of ancestors is calculated as
![]() |
where pi is the marginal genetic contribution of ancestor i. Calculated this way, the effective number of ancestors is a measure of genetic diversity that accounts for recent bottlenecks, and thus account partly for the loss of allelic diversity since the foundation population. The marginal contribution is not counted for all ancestors, but for a given number of ancestors the upper and lower limits to the effective number of ancestors. In this study, we used 1000 ancestors, in which case the upper and lower limits does not show any difference for 2 decimal digits.
The effective number of ancestors is an ad hoc measure in the sense that it does not fit into the theory of long-term genetic contributions (Wray and Thompson, 1990; Woolliams and Thompson, 1994). In this theory, the rate of inbreeding and loss of genetic diversity is a function of total genetic contributions and not of marginal genetic contributions. In essence, the marginal genetic contributions ignore the generational structure of the population, and a number of highly contributing ancestors are treated as founders when calculating the effective number of ancestors. Thus, the effective number of ancestors does not tell the full story of the genetic diversity and is also somewhat dependent on the depth of pedigrees. However, it is useful in comparison with the effective number of founders. The ratio of the 2 is an indication of the importance of bottlenecks in the development of the population. If the ratio is close to unity, the population has been stable in terms of numbers of effectively contributing animals. If the effective number of founders is larger than the effective number of ancestors, bottlenecks have played a role in population formation.
Effective number of founder genomes.
The effective number of founder genomes, fg, accounts for both unequal contributions of founders, bottlenecks, and random loss of alleles due to genetic drift (Lacy, 1989, 1995). It is defined similar to the effective number of founders with the difference that the genetic contribution of the I'th founder to the reference population, qi, is modified by the proportion of the founders genes that are found in the reference population, ri. As (1 ri) quantifies the proportion of alleles from the founder that are not expected to be present in the reference population, fg takes account of random loss of alleles during bottlenecks. The effective number of founder genomes is calculated as
![]() |
With this definition, the interpretation is the number of equally contributing founders with no loss of founder alleles that would be expected to produce the same amount of diversity as in the reference population (Lacy, 1995). Even in the case of minimum inbreeding, the effective number of founder genomes is smaller than the effective number of founders and the effective number of ancestors and smaller than half the effective population size. The degree to which the effective number of founder genomes is smaller is an indication of the degree of random loss of alleles. As alleles are lost every generation, the effective number of founder genomes decreases every generation and is therefore sensitive to depth of pedigree.
In this study, the effective number of founder genomes is calculated using the software package Pedig by Boichard (2002). It follows the genedrop procedure described by Boichard et al. (1997), and calculates the average of the inverse of the summations from each genedrop simulation. However, the program was modified according to the suggestions of Caballero and Toro (2000), so that it takes the inverse of the average of summations instead. The difference between the 2 methods of calculation is small.
Software.
The software package Pedig by Boichard (2002) was used to calculate generation intervals, effective number of founders, effective number of ancestors, effective number of founder genomes, and marginal contributions of ancestors. The software package Inbred by Berg (2003) was used to calculate coefficients of inbreeding, using the algorithm of Meuwissen and Luo (1992), pedigree completeness indices, number of complete generation equivalents, and average coancestry within birth cohorts.
| RESULTS |
|---|
|
|
|---|
|
The effective numbers of founders (Table 3
) show large differences between breeds. It is largest for Danish Red and smallest for Danish Holstein. The effective numbers of ancestors are smaller and more equal for the 3 breeds. The effective numbers of founder genomes are very small for Danish Holstein and Danish Jersey. The number is somewhat larger for Danish Red.
|
|
|
| DISCUSSION |
|---|
|
|
|---|
In animal breeding, the recommendation is to maintain an effective population size of at least 50 to 100 (FAO, 1998; Bijma, 2000). All 3 breeds fall within or below this interval, and so breeders should be concerned of controlling inbreeding in their breeding programs. This recommendation is by no means a magic number, but has been derived from theoretical arguments, where natural selection counteracts inbreeding depression. An effective population size of at least 500 is needed if genetic variation in the long term should not decrease (see Franklin and Frankham, 1998). This argument is the basis for the recommendation from associations of zoos (Wheater et al., 1993).
The different measures of diversity assess different aspects of diversity. The most important is the rate of increase of coancestry as it holds all currently available information on the future rate of inbreeding. Tools designed to balance the rate of gain and the rate of inbreeding in the future (Wray and Goddard, 1994; Meuwissen, 1997) constrain the rate of increase of coancestry when selecting animals. According to the results of this study, control of future rate of inbreeding is necessary for all breeds. These methods are currently being implemented for this purpose.
The other diversity measures are historical with limited predictive value. However, together they hold information on the history of population formation. The effective population size measures the rate of loss of genetic diversity that has occurred historically. For all breeds, their rate of loss has been large in the last decade. For all breeds, the difference between half the effective population size and the effective number of founders shows that drift has been accelerating since the founder generation, primarily due to increased selection pressure. Also for all breeds, the difference between the effective number of founders and the effective number of ancestors show that bottlenecks have occurred since the foundation of the population. According to the difference of half the effective population size and the effective number of founder genomes, random loss of founder alleles due to drift has been larger in Danish Holstein and Danish Jersey than in Danish Red.
This study illustrates that current breeding strategies can result in depletion of genetic diversity even in very large populations. All 3 breeds have considerable census sizes. However, their effective population sizes are small and decreasing. Despite the difference in breeding strategy between the Danish Red and the other breeds, the loss of genetic diversity seems to occur at the same rate. For all breeds, active management of the future rate of inbreeding is necessary.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication November 9, 2004. Accepted for publication January 18, 2005.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. Mrode, J. F. Kearney, S. Biffani, M. Coffey, and F. Canavesi Short communication: Genetic relationships between the Holstein cow populations of three European dairy countries J Dairy Sci, November 1, 2009; 92(11): 5760 - 5764. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. H. Buch, A. C. Sorensen, J. Lassen, P. Berg, L. G. Christensen, and M. K. Sorensen Factors affecting the exchange of genetic material between Nordic and US Holstein populations J Dairy Sci, August 1, 2009; 92(8): 4023 - 4034. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. D. Pedersen, A. C. Sorensen, and P. Berg Marker-assisted selection can reduce true as well as pedigree-estimated inbreeding J Dairy Sci, May 1, 2009; 92(5): 2214 - 2223. [Abstract] [Full Text] [PDF] |
||||
![]() |
M.-H. Li, I. Stranden, and J. Kantanen Genetic diversity and pedigree analysis of the Finnsheep breed J Anim Sci, May 1, 2009; 87(5): 1598 - 1605. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. K. Sorensen, E. Norberg, J. Pedersen, and L. G. Christensen Invited Review: Crossbreeding in Dairy Cattle: A Danish Perspective J Dairy Sci, November 1, 2008; 91(11): 4116 - 4128. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. P. W. de Roos, B. J. Hayes, R. J. Spelman, and M. E. Goddard Linkage Disequilibrium and Persistence of Phase in Holstein-Friesian, Jersey and Angus Cattle Genetics, July 1, 2008; 179(3): 1503 - 1512. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A. Martinez, D. Garcia, J. L. Gallego, G. Onofre, J. Perez, and J. Canon Genetic variability in Colombian Creole cattle populations estimated by pedigree information J Anim Sci, March 1, 2008; 86(3): 545 - 552. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Carolino and L. T. Gama Indicators of genetic erosion in an endangered population: The Alentejana cattle breed in Portugal J Anim Sci, January 1, 2008; 86(1): 47 - 56. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Hammami, C. Croquet, J. Stoll, B. Rekik, and N. Gengler Genetic Diversity and Joint-Pedigree Analysis of Two Importing Holstein Populations J Dairy Sci, July 1, 2007; 90(7): 3530 - 3541. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Norberg and A. C. Sorensen Inbreeding trend and inbreeding depression in the Danish populations of Texel, Shropshire, and Oxford Down J Anim Sci, February 1, 2007; 85(2): 299 - 304. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Mc Parland, J. F. Kearney, M. Rath, and D. P. Berry Inbreeding trends and pedigree analysis of Irish dairy and beef cattle populations J Anim Sci, February 1, 2007; 85(2): 322 - 331. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. C. Sorensen, P. Madsen, M. K. Sorensen, and P. Berg Udder health shows inbreeding depression in danish holsteins. J Dairy Sci, October 1, 2006; 89(10): 4077 - 4082. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Sewalem, G. J. Kistemaker, F. Miglior, and B. J. Van Doormaal Analysis of inbreeding and its relationship with functional longevity in Canadian dairy cattle. J Dairy Sci, June 1, 2006; 89(6): 2210 - 2216. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |