| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Genetics, Vol. 173, 1503-1510, July 2006, Copyright © 2006
doi:10.1534/genetics.105.053280
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


* Smurfit Institute of Genetics, Trinity College, Dublin 2, Ireland,
Imperial College, London SW7 2AZ, United Kingdom and
International Livestock Research Institute, Nairobi, Kenya
1 Corresponding author: Smurfit Institute of Genetics, Trinity College, Dublin 2, Ireland.
E-mail: dbradley{at}tcd.ie
| ABSTRACT |
|---|
|
|
|---|
10,000 years ago and are now found throughout northern Eurasia and certain regions of Africa. B. indicus or zebu cattle were domesticated from a different variant of wild progenitor on the Indian subcontinent. Substantial genetic exchange has taken place between the two types. Zebu cattle were brought from southern Asia to the Middle East and the breeds that are found in this area today are affected by hybridization that took place perhaps >3000 years ago (LOFTUS et al. 1999). Nuclear genetic markers including Y chromosome polymorphisms (BRADLEY et al. 1994; HANOTTE et al. 2000) and autosomal microsatellites (MACHUGH et al. 1997) have been used to show that there has been a widespread introgression of zebu cattle into the African continent, which would have initially been inhabited only by African B. taurus, pure examples of which remain only in certain parts of North and West Africa.
The genetic exchange between taurine and zebu cattle represents an appropriate model to study the effects of postadmixture haplotypic decay through recombination, mutation, and genetic drift. Four of the domestic cattle hybrid zones, the Middle East plus East, West, and southern Africa, have breeds with roughly equal amounts, but different ages of admixture (LOFTUS et al. 1999; HANOTTE et al. 2002). The first-generation population created after an admixture event will contain intact chromosomes from each of the inputting parental populations. However, as mating occurs within the hybrid population these intact chromosomes are mixed by recombination and gradually become more mosaic with each generation. Therefore, the extent of genomic blocks undisturbed by hybrid recombination will contain information about the history of admixture.
The examination of closely linked markers can be used to estimate how parental haplotypes have been broken down over time in the hybrid population. Microsatellite allele frequency spectra have been shown to differ sharply between B. indicus and B. taurus, facilitating the assignment of ancestry at specific chromosome positions. Moreover, extended haplotypes may be unambiguously assigned to X chromosomes in males.
Here we investigate differences in haplotypic mosaicism among 346 X chromosomes typed in male cattle from regions of differing phylogeographic history. We estimate microsatellite allele ancestry at 10 loci within a 10-cM region to examine the breakdown of genomic blocks and assess levels of linkage disequilibrium (LD) in populations with differing histories of admixture. We relate the different levels of haplotypic mosaicism and LD to the different time depths and manners of introgression. For simplicity, we use the term LD to describe both traditional linkage disequilibrium and pairwise disequilibrium between nonlinked markers.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
|
Statistical inference on admixture and the ancestry crossover rate in the admixed cattle populations was made in the computer program ADMIXMAP by simulating posterior samples of the model parameters using Markov chain Monte Carlo simulation (HOGGART et al. 2003, 2004). ADMIXMAP implements a full Bayesian hierarchical model for admixture at the population, individual, and locus levels and also models the allele frequencies in the unadmixed subpopulations that have contributed to the admixed population under study. The stochastic variation of ancestry between loci on each chromosome is modeled by independent Poisson arrival processes, one for each unadmixed population contributing to the admixed population under study, in our case B. taurus and B. indicus. Thus inference is made on parameters representing population and individual admixture, indicator parameters for locus ancestry for all loci and all individuals, parameters of multinomial distributions representing allele and haplotype frequencies in the two populations, and
, the sum-of-intensities of the two Poisson arrival processes. While admixture was modeled separately for each individual in a population (independent and identically distributed from the population level admixture distribution),
was assumed to be the same for all individuals in the population.
For an individual with admixture proportion
from B. taurus the intensities per morgan of the arrival processes for B. taurus and B. indicus are
and
, respectively. The ancestry crossover rate
(the rate per morgan of transitions between ancestry of one subpopulation to ancestry of the other) is given by
; thus the larger
is the more frequent the crossovers. Since the X chromosome undergoes recombination only in female gametes, a different
will be applicable to the X chromosome than to the autosomes. However, in our data sets the autosomal markers are unlinked and thus
was estimated from the linked markers typed in the 10-cM region of the X chromosome alone. In contrast, individual admixture proportions (
) were estimated from both the autosomal and X chromosome markers. When introgression occurs in a single pulse the expectation of
on the autosomes is equal to the number of generations since admixture (FALUSH et al. 2003). However, admixture in the cattle populations is not believed to have occurred as a single event but rather through continuous gene flow of B. indicus into B. taurus populations. The effect of continuous gene flow on
is unclear and thus it cannot be interpreted as the number of generations since admixture but as a relative measure of the age of admixture among the populations studied. Given the lack of interpretation of
it was assigned an uninformative prior, log
-uniform. Prior distributions for population-specific allele frequencies for B. taurus and B. indicus were set from allele counts in the five B. taurus populations (both European and West African) and two B. indicus populations listed in Table 1. While it is known that the West African B. taurus populations contain a degree of indicus admixture, this is very small (mR = 3, 7, and 10% in Guinean N'Dama, Guinea Bissau N'Dama, and Somba, respectively) and thus will have minimal effect on the allele counts. Furthermore, specifying prior distributions from the allele counts rather than fixing them at their mean values allows for error in the frequencies to be accounted for. Uninformative uniform priors were assigned for the population level admixture.
The ADMIXMAP analyses assume that all LD between markers can be attributed to admixture LD and LD due to ancestral mosaicism but ignore background LD. The X chromosome markers used in our analyses are between 0.1 and 2.5 Mb apart (Table 2), and at such distances it is a reasonable approximation to ignore background LD.
Statistical testing was carried out using SPSS (v. 12.0) and values of
were plotted with sampling location using the ArcView GIS software package. Equal interval contours join areas with similar values of
(Environmental Systems Research Institute, Redlands, CA). A synthetic map was finalized using the Adobe Illustrator (v. 8.0) package.
| RESULTS |
|---|
|
|
|---|
The strength of pairwise LD between all typed markers was estimated and the results were divided into two groups according to the significance of association between pairs of markers: P < 0.05 and P < 0.01. Figure 2 shows matrices of LD significance levels for all possible pairwise comparisons of the 10 X chromosome loci in their chromosomal order and for the unlinked markers. The samples were divided into regional groupings according to Table 1, excluding the Zebu Malagasy breed, which does not fall easily into an appropriate category and is represented only by four individuals. The unlinked autosomal microsatellite markers provide a measure of the background levels of LD in each regional group. LD measures for linked markers are shown on the left-hand side of the thick solid line in Figure 2, in the order that they are located on the X chromosome.
|
Estimates of
obtained from ADMIXMAP enable comparison of time since admixture among the different sampled populations, but cannot be taken as absolute values, as a single event is unlikely to be an accurate model of past hybridization processes. Values of
are shown on a synthetic contour map (Figure 3). Despite the relatively small number of hybrid breeds (15), some geographical congruency can be observed for values of
.
|
for the admixed breeds are shown in Table 1. We see that within the African continent there are two main trends among the hybrid populations. First, values of
tend to be higher in the south of the continent, and second, there is a decreasing cline in values of
from east to west in the northern part of the continent. To test for significant difference between hybrid breeds in the north and south we pooled the estimates of
using the inverse variance method (DEEKS et al. 2001). The pooled estimate of
for hybrid breeds in the northern part of the continent (Mbororo, Gobra, Sokoto Gudali, Arashie, and Kuri) is 7.6 (standard deviation 1.66), and the pooled estimate for hybrids in the south (Ogaden, Ankole, Afrikaner, Nguni, Watusi, Zebu Malagasey, Barotse, and Kavango) is 13.9 (standard deviation 1.85). The test for difference between the means of the groups is significant, P = 0.012. The pooled estimate for
in the two Near Eastern hybrids is 15.4 (standard deviation 3.54), with the highest value being observed in the easternmost Iraqi breed (
= 17.8). | DISCUSSION |
|---|
|
|
|---|
The amount of admixture in the breeds used in this study had previously been estimated using autosomal microsatellites (MACHUGH et al. 1997; LOFTUS et al. 1999; HANOTTE et al. 2002; FREEMAN et al. 2004) (Figure 1). Here we set out to examine the pattern of haplotype decay in the bovine genome after admixture by genotyping closely linked microsatellite markers. The level of divergence between the contributing parental populations will determine whether it is possible to assign definite ancestry to an allele in a hybrid individual. In a previous study, we found that >80% of microsatellite alleles surveyed at 20 loci in European, Near Eastern, and Indian cattle populations showed >50% frequency differential (scaled by the absolute allele frequency) between B. indicus and B. taurus (KUMAR et al. 2003). All of the markers used here are located within a 10-cM region of the X chromosome and were typed only in males. This allowed the unambiguous assignment of a haplotype to an individual without pedigree information.
We found LD in this region of the X chromosome in all of the cattle populations and between linked and unlinked markers in all of the populations except for the Near Eastern hybrids. These results concur with those of FARNIR et al. (2000), who surveyed 284 autosomal markers and found that intrachromosomal LD extends over several tens of centimorgans and is common between nonsyntenic loci in dairy cattle populations. In that case, the high levels of overall LD may have been a result of a small effective population size (LONJOU et al. 1999), which could be as low as 50 for 1.2 million HolsteinFriesian females, or admixture (FARNIR et al. 2000).
Recurrent introgression has been cited previously as the cause of higher than expected levels of LD in AfricanAmerican human populations (LAUTENBERGER et al. 2000). Here, the very high level in the West African hybrid population may partially be explained by recent and continuous gene flow to this region (PFAFF et al. 2001). Pfaff and colleagues argue that maintenance of LD over relatively large (
10 cM) chromosomal segments is a not typical of a single admixture event, but is a characteristic of a continuous gene flow pattern of admixture. While the initial introgression of B. indicus chromosomes into Africa is likely ancient, zebutaurine admixture in West Africa probably began only recently in many populations where it had previously been constrained by the tsetse/trypanosomiasis infestation zone, a barrier to introgression. Understanding the patterns of LD in the hybrid populations from the surround of the West African tsetse zone will prove useful if these populations are used to map genes that influence disease traits (LAAN and PAABO 1997).
In the absence of mutation, every portion of a chromosome in a hybrid individual should be directly traceable to one of the individual's ancestors. When two distinct gene pools meet, a mosaic genome is produced (CHAPMAN and THOMPSON 2002), containing intact segments from each contributing founder population (NORDBORG and TAVARE 2002). Genomes that have undergone recent hybridization will have largely intact ancestral haplotypes, while recombination will act to disrupt these haplotypes in subsequent generations (WIEHE et al. 2000). If admixture is ancient, reciprocal recombination, genetic drift, and mutation will have eroded the association between markers from parental populations and a greater density of markers will be required to detect significant associations.
Here we employed a Bayesian method to assign probable ancestry to each allele in every X chromosome haplotype. The amount of fragmentation in haplotypes was assessed by the posterior distribution of
, and this was used as a measure for the age of admixture. The mean values of
(Table 1) for each of the cattle breeds are between 4.27 and 21.3, which is substantially less than the number of generations since admixture is believed to have first occurred. However, as discussed earlier, in the case of cattle a single admixture event may not accurately model the historical reality and
has expectation only equal to the number of generations since admixture in the case of a single admixture event. Furthermore, the prior has the effect of favoring small values of
.
The values obtained appear to show some geographical consistency, whereby most hybrid populations from southern Africa and the Near East appear to have experienced older admixture than those in the more northern part of the African continent. Furthermore there is a decline in age of admixture from east to west in the latter. The differences between the east and the west of the continent support a scenario of earlier introgression of B. indicus material into the former populations, as the degree of recombination would be expected to accumulate with time. This finding is supported by high levels of LD in West African hybrids and by previous work that has shown that the major B. indicus introductions occurred in the east, probably in the region of the Horn (HANOTTE et al. 2002).
The difference in
-values between the south and the north of Africa may be a signal of different histories of cattle in the two regions. However, it is well established that the earliest domesticated cattle of Africa were found in the north of the continent and were B. taurus in nature and that introgression of B. indicus into the continent via the Horn of Africa was secondary. Our results suggest that these earliest hybrids may have been the animals that later gave rise to southern populations. This scenario implies that B. indicus introgression occurred before pastoralism spread extensively southward and indeed, despite extensive sampling, there is no evidence for the presence of pure African B. taurus in the south of the continent (HANOTTE et al. 2002).
On the basis of archaeological evidence the origins of B. indicus introgression in Africa are controversial. Some authors have suggested that they were found as early as 35002000 years before the present (YBP) (HASSAN 2000; MARSHALL 2000; PARIS 2000), but these studies are inconclusive, and the earliest verified archaeological date for their presence is in the Horn of Africa in the second century A.D. (MARSHALL 2000). The archaeological record shows that B. taurus cattle were present in northern Africa and the Sahara from
6000 YBP (HASSAN 2000) and spread as far south as Kenya by 3000 YBP (HASSAN 2000; MITCHELL 2005). Further southward movement may have been hampered by ecological factors such as the many sub-Saharan disease challenges to which domestic cattle would be maladapted (GIFFORD-GONZALEZ 2000; MITCHELL 2005), and relatively few cattle were present in the south until 1000 A.D. (SMITH 2000). It is therefore possible that nomadic pastoralists in eastern Africa were herding animals composed of the two cattle lineages before the expansion of domestic cattle to the south.
Our sample from Madagascar suggests that admixture here was older but this sample is very small, n = 4, and should not be overinterpreted. There is limited archaeological evidence for early arrival of B. indicus cattle to Madagascar and zebu cattle occupy a place of prominence in Malagasy ritual and culture, which is compatible with an extensive history on the island (HEMMER 1990; VAN DER ZWAN and EVERS 1998).
The employment of haplotypic mosaicism as an indicator of the antiquity of admixture receives some confirmation from the Near Eastern results. The region between the primary centers of B. indicus and B. taurus domestication is an area of ancient interaction between the two lineages. As such, the sample from Iraq yields one of the highest values of
within the haplotypes.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| LITERATURE CITED |
|---|
|
|
|---|
AMARAL, M. E., S. R. KATA and J. E. WOMACK, 2002 A radiation hybrid map of bovine X chromosome (BTAX). Mamm. Genome 13: 268271.[CrossRef]
BISHOP, M. D., S. M. KAPPES, J. W. KEELE, R. T. STONE, S. L. SUNDEN et al., 1994 A genetic linkage map for cattle. Genetics 136: 619639.[Abstract]
BRADLEY, D. G., D. E. MACHUGH, R. T. LOFTUS, R. S. SOW, C. H. HOSTE et al., 1994 Zebu-taurine variation in Y chromosomal DNA: a sensitive assay for genetic introgression in west African trypanotolerant cattle populations. Anim. Genet. 25: 712.[Medline]
BRADLEY, D. G., D. E. MACHUGH, P. CUNNINGHAM and R. T. LOFTUS, 1996 Mitochondrial diversity and the origins of African and European cattle. Proc. Natl. Acad. Sci. USA 93: 51315135.
CHAPMAN, N. H., and E. A. THOMPSON, 2002 The effect of population history on the lengths of ancestral chromosome segments. Genetics 162: 449458.
DEEKS, J., D. ALTMAN and M. BRADBURN, 2001 Statistical methods for examining heterogeneity and combining results from several studies in meta-analysis, pp. 285372 in Systematic Reviews in Health Care. Meta-Analysis in Context, edited by M. EGGER, G. DAVEY SMITH and D. ALTMAN. BMJ Books, London.
FALUSH, D., M. STEPHENS and J. K. PRITCHARD, 2003 Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164: 15671587.
FARNIR, F., W. COPPIETERS, J. J. ARRANZ, P. BERZI, N. CAMBISANO et al., 2000 Extensive genome-wide linkage disequilibrium in cattle. Genome Res. 10: 220227.
FELIUS, M., 1995 Cattle BreedsAn Encyclopedia. Misset, Doetinchem, The Netherlands.
FREEMAN, A. R., C. M. MEGHAN, D. E. MACHUGH, R. T. LOFTUS, M. D. ACHUKWI et al., 2004 Admixture and diversity in West African cattle populations. Mol. Ecol. 13: 3477.
GIFFORD-GONZALEZ, D., 2000 Animal disease challenges to the emergence of pastoralism in sub-Saharan Africa. Afr. Archaeol. Rev. 17: 95139.
HANOTTE, O., C. L. TAWAH, D. G. BRADLEY, M. OKOMO, Y. VERJEE et al., 2000 Geographic distribution and frequency of a taurine Bos taurus and an indicine Bos indicus Y specific allele amongst sub-Saharan African cattle breeds. Mol. Ecol. 9: 387396.[CrossRef][Medline]
HANOTTE, O., D. G. BRADLEY, J. W. OCHIENG, Y. VERJEE, E. W. HILL et al., 2002 African pastoralism: genetic imprints of origins and migrations. Science 296: 336339.
HASSAN, F. A., 2000 Climate and cattle in North Africa: a first approximation, pp. 6186 in The Origins and Development of African Livestock: Archaeology, Genetics, Linguistics and Ethnography, edited by R. M. BLENCH and K. C. MACDONALD. UCL Press, London.
HEMMER, H. (Editor), 1990 Domestication: The Decline of Environmental Appreciation, p. 11. Cambridge University Press, Cambridge/London/New York.
HOGGART, C. J., E. J. PARRA, M. D. SHRIVER, C. BONILLA, R. A. KITTLES et al., 2003 Control of confounding of genetic associations in stratified populations. Am. J. Hum. Genet. 72: 14921504.[CrossRef][Medline]
HOGGART, C. J., M. D. SHRIVER, R. A. KITTLES, D. G. CLAYTON and P. M. MCKEIGUE, 2004 Design and analysis of admixture mapping studies. Am. J. Hum. Genet. 74: 965978.[CrossRef][Medline]
IHARA, N., A. TAKASUGA, K. MIZOSHITA, H. TAKEDA, M. SUGIMOTO et al., 2004 A comprehensive genetic map of the cattle genome based on 3802 microsatellites. Genome Res. 14: 19871998.
KAPPES, S. M., J. W. KEELE, R. T. STONE, R. A. MCGRAW, T. S. SONSTEGARD et al., 1997 A second-generation linkage map of the bovine genome. Genome Res. 7: 235249.
KUMAR, P., A. R. FREEMAN, R. T. LOFTUS, C. GAILLARD, D. Q. FULLER et al., 2003 Admixture analysis of South Asian cattle. Heredity 91: 4350.[CrossRef][Medline]
LAAN, M., and S. PAABO, 1997 Demographic history and linkage disequilibrium in human populations. Nat. Genet. 17: 435438.[CrossRef][Medline]
LAUTENBERGER, J. A., J. C. STEPHENS, S. J. O'BRIEN and M. W. SMITH, 2000 Significant admixture linkage disequilibrium across 30 cM around the FY locus in African Americans. Am. J. Hum. Genet. 66: 969978.[CrossRef][Medline]
LOFTUS, R. T., D. E. MACHUGH, D. G. BRADLEY, P. M. SHARP and P. CUNNINGHAM, 1994 Evidence for two independent domestications of cattle. Proc. Natl. Acad. Sci. USA 91: 27572761.
LOFTUS, R. T., O. ERTUGRUL, A. H. HARBA, M. A. EL-BARODY, D. E. MACHUGH et al., 1999 A microsatellite survey of cattle from a centre of origin: the Near East. Mol. Ecol. 8: 20152022.[CrossRef][Medline]
LONJOU, C., A. COLLINS and N. E. MORTON, 1999 Allelic association between marker loci. Proc. Natl. Acad. Sci. USA 96: 16211626.
MACHUGH, D. E., M. D. SHRIVER, R. T. LOFTUS, P. CUNNINGHAM and D. G. BRADLEY, 1997 Microsatellite DNA variation and the evolution, domestication and phylogeography of taurine and zebu cattle (Bos taurus and Bos indicus). Genetics 146: 10711086.[Abstract]
MARSHALL, F., 2000 The origins and spread of domestic animals in East Africa, pp. 191221 in The Origins and Development of African Livestock: Archaeology, Genetics, Linguistics and Ethnography, edited by R. M. BLENCH and K. C. MACDONALD. UCL Press, London.
MITCHELL, P., 2005 African Connections, pp. 3363. AltaMira, Walnut Creek, CA.
NORDBORG, M., and S. TAVARE, 2002 Linkage disequilibrium: what history has to tell us. Trends Genet. 18: 8390.[CrossRef][Medline]
PARIS, F., 2000 African livestock remains from Saharan mortuary context, pp. 111126 in The Origins and Development of African Livestock: Archaeology, Genetics, Linguistics and Ethnography, edited by R. M. BLENCH and K. C. MACDONALD. UCL Press, London.
PFAFF, C. L., E. J. PARRA, C. BONILLA, K. HIESTER, P. M. MCKEIGUE et al., 2001 Population structure in admixed populations: effect of admixture dynamics on the pattern of linkage disequilibrium. Am. J. Hum. Genet. 68: 198207.[CrossRef][Medline]
RAYMOND, M., and F. ROUSSET, 1995 GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J. Hered. 86: 248249.
RITZ, L. R., M. L. GLOWATZKI-MULLIS, D. E. MACHUGH and C. GAILLARD, 2000 Phylogenetic analysis of the tribe Bovini using microsatellites. Anim. Genet. 31: 178185.[CrossRef][Medline]
SMITH, A. B., 2000 The origins of the domesticated animals of southern Africa, pp. 222238 in The Origins and Development of African Livestock: Archaeology, Genetics, Linguistics and Ethnography, edited by R. M. BLENCH and K. C. MACDONALD. UCL Press, London.
SØNSTEGARD, T. S., N. L. LOPEZ-CORRALES, S. M. KAPPES, R. T. STONE, S. AMBADY et al., 1997 An integrated genetic and physical map of the bovine X chromosome. Mamm. Genome 8: 1620.[CrossRef][Medline]
SØNSTEGARD, T. S., W. BARENDSE, G. L. BENNETT, G. A. BROCKMANN, S. DAVIS et al., 2001 Consensus and comprehensive linkage maps of the bovine sex chromosomes. Anim. Genet. 32: 115117.[CrossRef][Medline]
SUN, H. S., M. R. DENTINE, W. BARENDSE and B. W. KIRKPATRICK, 1994 UWCA19 and UWCA20: polymorphic bovine microsatellites. Anim. Genet. 25: 121.[Medline]
VAN DER ZWAN, N., and S. EVERS, 1998 Madagascar: The Zebu as Guide Through Past and Present. Berg en Dal, The Netherlands.
WIEHE, T., J. MOUNTAIN, P. PARHAM and M. SLATKIN, 2000 Distinguishing recombination and intragenic gene conversion by linkage disequilibrium patterns. Genet. Res. 75: 6173.[CrossRef][Medline]
WILLIAMS, J. L., A. EGGEN, L. FERRETTI, C. J. FARR, M. GAUTIER et al., 2002 A bovine whole-genome radiation hybrid panel and outline map. Mamm. Genome 13: 469474.
Communicating editor: C. HALEYThis article has been cited by other articles:
![]() |
M. Gautier, T. Faraut, K. Moazami-Goudarzi, V. Navratil, M. Foglio, C. Grohs, A. Boland, J.-G. Garnier, D. Boichard, G. M. Lathrop, et al. Genetic and Haplotypic Structure in 14 European and African Cattle Breeds Genetics, October 1, 2007; 177(2): 1059 - 1070. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |