Genetic diversity and admixture of indigenous cattle from North Ethiopia: implications of historical introgressions in the gateway region to Africa

Microsatellite variation was surveyed to determine the genetic diversity, population structure and admixture of seven North Ethiopian cattle breeds by combining multiple microsatellite data sets of Indian and West African zebu, and European, African and Near-Eastern taurine in genetic analyses. Based on allelic distribution, we identified four diagnostic alleles (HEL1-123 bp, CSSM66-201 bp, BM2113-150 bp and ILSTS6-285 bp) specific to the Near-Eastern taurine. Results of genetic relationship and population structure analyses confirmed the previously established marked genetic distinction between taurine and zebu, and indicated further divergence among the bio-geographical groupings of breeds such as North Ethiopian, Indian and West African zebu, and African, European and Near-Eastern taurine. Using the diagnostic alleles for bio-geographical groupings and a Bayesian method for population structure inference, we estimated the genetic influences of major historical introgressions in North Ethiopian cattle. The breeds have been heavily (>90%) influenced by zebu, followed by African, European and the Near-Eastern taurine. Overall, North Ethiopian cattle show a high level of within-population genetic variation (e.g. observed heterozygosity = 0.659-0.687), which is in the upper range of that reported for domestic cattle and indicates their potential for future breeding applications, even in a global context. Rather low but significant population differentiation (F(ST) = 1.1%, P < 0.05) was recorded as a result of multiple introgression events and strong genetic exchanges among the North Ethiopian breeds.     

Detection of selection signatures within candidate regions underlying trypanotolerance in outbred cattle populations

Breeding indigenous African taurine cattle tolerant to trypanosomosis is a straightforward approach to control costs generated by this disease. A recent study identified quantitative trait loci (QTL) underlying trypanotolerance traits in experimental crosses between tolerant N'Dama and susceptible Boran zebu cattle. As trypanotolerance is thought to result from local adaptation of indigenous cattle breeds, we propose an alternative and complementary approach to study the genetic architecture of this trait, based on the identification of selection signatures within QTL or candidate genes. A panel of 92 microsatellite markers was genotyped on 509 cattle belonging to four West African trypanotolerant taurine breeds and 10 trypanosusceptible European or African cattle breeds. Some of these markers were located within previously identified QTL regions or candidate genes, while others were chosen in regions assumed to be neutral. A detailed analysis of the genetic structure of these different breeds was carried out to confirm a priori grouping of populations based on previous data. Tests based on the comparison of the observed heterozygosities and variances in microsatellite allelic size among trypanotolerant and trypanosusceptible breeds led to the identification of two significantly less variable microsatellite markers. BM4440, one of these two outlier loci, is located within the confidence interval of a previously described QTL underlying a trypanotolerance-related trait.
Detection of selection signatures appears to be a straightforward approach for unravelling the molecular determinism of trypanosomosis pathogenesis. We expect that a whole genome approach will help confirm these results and achieve a higher resolving power.     

Comprehensive analysis of the mitochondrial DNA diversity in Chinese cattle

Complete mitochondrial DNA D‐loop sequences of 1105 individuals were used to assess the diversity of maternal lineages of cattle populations in China. In total, 250 taurine and 88 zebu haplotypes were identified. Five main haplogroups—T1a, T2, T3, T4 and T5—were identified in Bos taurus, whereas Bos indicus harbored two haplogroups—I1 and I2. Our results suggest that the distribution of T1a in Asia was concentrated mainly in the northeast region (northeast China, Korea and Japan); haplogroups T2, T3 and T4 were predominant in Chinese cattle; and T5 was sporadically detected in Mongolian and Pingwu cattle. In contrast to the widespread presence of I1, I2 was distributed only in southwestern China (Yunnan‐Guizhou Plateau and the Tibet Autonomous Region) and Xinjiang Uygur Autonomous Region. This is the first time that all five taurine haplogroups and two zebu haplogroups have been found in Mongolian cattle. In addition, eight individuals in Tibetan cattle carried the Bos grunniens mtDNA type. The high mtDNA diversity (H = 0.904 ± 0.008) and the weak genetic structure among the 57 Chinese cattle breeds/populations are consistent with their complex historical background, migration route and ecological environment.     

Admixture and diversity in West African cattle populations

We present a population genetic analysis of microsatellite variation in 16 West African cattle populations. West Africa represents a unique juxtaposition of different climatic and ecological zones in a relatively small geographical area. While more humid coastal regions are inhabited by the tsetse fly, a vector which spreads trypanosomiasis among cattle, the disease is not transmitted in the drier areas outside this zone. This is the most thorough study of genetic diversity in cattle within this area, which contains genetically important trypanotolerant Bos taurus breeds. Genetic relationships among the many breeds are examined and levels of diversity are assessed. Admixture levels were determined using a variety of methods. Ancestry informative or population-associated alleles (PAAs) were selected using populations from India, the Near East and Europe. Multivariate analysis, the admix program and model-based Bayesian admixture analysis approaches were also employed. These analyses reveal the direct impact of ecological factors and the profound effect of admixture on the cattle of this region. They also highlight the importance of efforts to prevent further dilution of African taurine breeds by B. indicus cattle.     

Genome-wide genetic diversity,population structure and admixture analysis in African and Asian cattle breeds

Knowledge about genetic diversity and population structure is useful for designing effective strategies to improve the production, management and conservation of farm animal genetic resources. Here, we present a comprehensive genome-wide analysis of genetic diversity, population structure and admixture based on 244 animals sampled from 10 cattle populations in Asia and Africa and genotyped for 69 903 autosomal single-nucleotide polymorphisms (SNPs) mainly derived from the indicine breed. Principal component analysis, STRUCTURE and distance analysis from high-density SNP data clearly revealed that the largest genetic difference occurred between the two domestic lineages (taurine and indicine), whereas Ethiopian cattle populations represent a mosaic of the humped zebu and taurine. Estimation of the genetic influence of zebu and taurine revealed that Ethiopian cattle were characterized by considerable levels of introgression from South Asian zebu, whereas Bangladeshi populations shared very low taurine ancestry. The relationships among Ethiopian cattle populations reflect their history of origin and admixture rather than phenotype-based distinctions. The high within-individual genetic variability observed in Ethiopian cattle represents an untapped opportunity for adaptation to changing environments and for implementation of within-breed genetic improvement schemes. Our results provide a basis for future applications of genome-wide SNP data to exploit the unique genetic makeup of indigenous cattle breeds and to facilitate their improvement and conservation.     

Historical demographic profiles and genetic variation of the East African Butana and Kenana indigenous dairy zebu cattle

Butana and Kenana breeds from Sudan are part of the East African zebu Bos indicus type of cattle. Unlike other indigenous zebu cattle in Africa, they are unique due to their reputation for high milk production and are regarded as dairy cattle, the only ones of their kind on the African continent. In this study, we sequenced the complete mitochondrial DNA (mtDNA) D‐loop of 70 animals to understand the maternal genetic variation, demographic profiles and history of the two breeds in relation to the history of cattle pastoralism on the African continent. Only taurine mtDNA sequences were identified. We found very high mtDNA diversity but low level of maternal genetic structure within and between the two breeds. Bayesian coalescent‐based analysis revealed different historical and demographic profiles for the two breeds, with an earlier population expansion in the Butana vis a vis the Kenana. The maternal ancestral populations of the two breeds may have diverged prior to their introduction into the African continent, with first the arrival of the ancestral Butana population. We also reveal distinct demographic history between the two breeds with the Butana showing a decline in its effective population size (N_e) in the recent past ~590 years. Our results provide new insights on the early history of cattle pastoralism in Sudan indicative of a large ancient effective population size.     

Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

Phylogenies using mtDNA and SRY provide evidence for male-mediated introgression in Asian domestic cattle

Using nucleotide sequences of the mitochondrial DNA (mtDNA) cytochrome b and SRY genes, we examined the genetic status of two major groups of domestic cattle, the humpless taurine (Bos taurus) and humped zebu (B. indicus), using 10 cattle populations in Asia. Several sequence polymorphisms specific for each major group were found, although the frequency of these polymorphisms varied in each population. Six major mtDNA-SRY composite types were observed. The Mishima, Mongolian, Korean, Chinese Yellow and Sri Lanka cattle populations had a full match between the mtDNA and SRY sequences, specifically the taurine/taurine type or zebu/zebu type. A non-match type (zebu/taurine type) was found at a high frequency in the Bangladesh (83.4%) and Nepal populations (83.3%). Our results suggest that these non-match type populations developed from genetic hybridization of different strains. Also, the domestication history of modern Asian domestic cattle could be explained by male-mediated introgression. Additionally, our results suggest the occurrence of introgression of mtDNA from other Bibos or Poephagus species into native cattle populations. The existence of other mtDNA-SRY composite types, such as the Bali-zebu and yak-zebu types in Indonesia (85.7%) and Nepal (16.7%), respectively, suggests that genetic introgression also occurred from other genera into domestic cattle during the process of domestication.     

The Origins of Domesticated Cattle

The origin of domestic cattle has perplexed archaeologists for more than a century. Researchers have proposed various theories, which offer alternative spatial and chronological models for the origin and spread of domesticated cattle. One point of discussion is whether domestic cattle had a single or multiple origins; however, most authorities considered that the first steps towards cattle domestication were taken in southwest Asia and that domesticated cattle entered Europe with pastoralists migrating from this region. Domesticated taurine cattle were thought to have entered Africa in successive waves from southwest Asia, while zebu cattle migrated into Africa at a later date from Arabia and the Indian subcontinent. Analysis of mitochondrial DNA (mtDNA) shows that taurine and zebu cattle divergence before the Holocene and were probably domesticated independently. Recent mtDNA sequence data shows that African and European taurine cattle were probably domesticated independently, but that there was a process of genetic introgression between taurine and zebu cattle in Africa. Ancient DNA studies over the last 10 years suggest that Northern European aurochsens apparently contributed little or nothing to domestic cattle while Southern European aurochsens apparently made a significant input. However, Middle Eastern aurochsen, unfortunately not typed yet, are expected to be to be very similar to European breeds as well, both because archeological data suggest that the major center of domestication for European Bos taurus breeds was the Fertile Crescent (9), and also because a mtDNA sequence from a Syrian specimen dated at 8,000–9,000 years ago shows a typical European haplotype found both in modern breeds and the Italian aurochsen. Evidence seems to suggest that small to moderate levels of local gene flow from wild Bos primigenius females in selected breeds were either accepted or may be reinforced by Neolithic breeders.     

Combination of multiple microsatellite data sets to investigate genetic diversity and admixture of domestic cattle

Microsatellite markers are commonly used for population genetic analyses of livestock. However, up to now, combinations of microsatellite data sets or comparison of population genetic parameters from different studies and breeds has proven difficult. Often different genotyping methods have been employed, preventing standardization of microsatellite allele calling. In other cases different sets of markers have been genotyped, providing differing estimates of population genetic parameters. Here, we address these issues and illustrate a general two-step regression approach in cattle using three different sets of microsatellite data, to combine population genetics estimates of diversity and admixture. This regression-based method is independent of the loci genotyped but requires common breeds in the data sets. We show that combining microsatellite data sets can provide new insights on the origin and geographical distribution of genetic diversity and admixture in cattle, which will facilitate global management of this livestock species.     

Geographic distribution of haplotype diversity at the bovine casein locus

The genetic diversity of the casein locus in cattle was studied on the basis of haplotype analysis. Consideration of recently described genetic variants of the casein genes which to date have not been the subject of diversity studies, allowed the identification of new haplotypes. Genotyping of 30 cattle breeds from four continents revealed a geographically associated distribution of haplotypes, mainly defined by frequencies of alleles at CSN1S1 and CSN3. The genetic diversity within taurine breeds in Europe was found to decrease significantly from the south to the north and from the east to the west. Such geographic patterns of cattle genetic variation at the casein locus may be a result of the domestication process of modern cattle as well as geographically differentiated natural or artificial selection. The comparison of African Bos taurus and Bos indicus breeds allowed the identification of several Bos indicus specific haplotypes (CSN1S1*C-CSN2*A²-CSN3*A_I/CSN3*H) that are not found in pure taurine breeds. The occurrence of such haplotypes in southern European breeds also suggests that an introgression of indicine genes into taurine breeds could have contributed to the distribution of the genetic variation observed.     

凉地方牛群体母系遗传背景研究

黄牛（Bos taurus）的驯养作为农耕文化重要内容之一,在我国历史悠久,长久以来形成的役用型逐步被肉用型所替代。为了揭示平凉地方牛群体的遗传背景,分析其是否具有生产优质牛肉的遗传基础,本研究测定了88头平凉地方牛群体mtDNA D-loop HVS序列,并对包括平凉地方牛群体在内的我国23个地方牛群体单倍型分布及系统发生关系进行了分析。结果表明,mtDNA D-loop 高变区,在平凉地方牛群体共有52个单倍型,而23个地方牛群体共有95个单倍型,这些单倍型在系统发生树和中介网络关系中分布于两个分枝,即瘤牛型和普通牛进化枝。因此,我们认为平凉地方牛群体和我国其他牛种一样,存在普通牛和瘤牛两个母系起源的遗传背景。     

Geographic distribution and frequency of a taurine Bos taurus and an indicine Bos indicus Y specific allele amongst sub-saharan African cattle breeds

We report for the first time, and for the whole of sub-Saharan Africa, the geographical distribution and the frequency of an indicine and a taurine Y specific allele amongst African cattle breeds. A total of 984 males from 69 indigenous African populations from 22 countries were analysed at the microsatellite locus INRA 124. The taurine allele is probably the oldest one on the continent. However, the taurine and the indicine alleles were present in 291 males (30%), and 693 males (70%), respectively. More particularly, 96% of zebu males (n = 470), 50% of taurine males (n = 263), 29% of sanga males (crossbreed Bos taurus x Bos indicus, n = 263) and 95% of zebu x sanga crossbred males (n = 56) had the indicine allele. The Borgou, a breed classified as zebu x taurine cross showed only the zebu allele (n = 12). The indicine allele dominates today in the Abyssinian region, a large part of the Lake Victoria region and the sahelian belt of West Africa. All the sanga males (n = 64) but only one from the Abyssinian region had the indicine allele. The taurine allele is the commonest only among the sanga breeds of the southern African region and the trypanotolerant taurine breeds of West Africa. In West Africa and in the southern Africa regions, zones of introgression were detected with breeds showing both Y chromosome alleles. Our data also reveal a pattern of male zebu introgression in Mozambique and Zimbabwe, probably originating from the Mozambique coast. The sanga cattle from the Lake Victoria region and the Kuri cattle of Lake Chad, cattle populations surrounded by zebu breeds were, surprisingly, completely devoid of the indicine allele. Human migration, phenotypic preferences by the pastoralists, adaptation to specific habitats and to specific diseases are the main factors explaining the present-day distribution of the alleles in sub-Saharan Africa.     

Genetic characterization of Argentine and Bolivian Creole cattle breeds assessed through microsatellites

In the present report, the polymorphisms from 9 microsatellites were used to assess genetic diversity and relationships in 4 Creole cattle breeds from Argentina and Bolivia, 4 European taurine breeds, and 2 American zebu populations. The Creole populations display a relatively high level of genetic variation as estimated by allelic diversity and heterozygosity, whereas the British breeds displayed reduced levels of genetic diversity. The analysis of molecular variance indicated that 7.8% of variance can be explained by differences among taurine and zebu breeds. Consistent with these results, the first principal component (PC), which comprised the 40% of the total variance, clearly distinguishes these 2 groups. In addition, all constructed phylogenetic trees cluster together Nelore and Brahman breeds with robust bootstrap values. Only 1% of variance was due to difference between American Creole and European taurine cattle. Although this secondary split was supported by the classical genetic distance and the second PC (15%), the topology of trees is not particularly robust. The presence of zebu-specific alleles in Creole cattle allowed estimating a moderate degree of zebu admixture. When these data were compared with mitochondrial and Y chromosomal studies, a clear pattern of male-mediated introgression was revealed. The results presented here contribute to the understanding of origin and history of the American Creole cattle.     

Dynamics of genomic architecture during composite breed development in cattle

Some livestock breeds face the challenge of reduced genetic variation, increased inbreeding depression owing to genetic drift and selection. Hybridization can reverse these processes and increase levels of productivity and adaptation to various environmental stressors. Samples from American Brangus were used to evaluate the indicine/taurine composition through nine generations (~45 years) after the hybridization process was completed. The purpose was to determine how hybridization alters allelic combinations of a breed over time when genetic factors such as selection and drift are operating. Furthermore, we explored genomic regions with deviations from the expected composition from the progenitor breeds and related these regions to traits under selection. The Brangus composition deviated from the theoretical expectation, defined by the breed association, of 62.5% taurine, showing taurine composition to be 70.4 ± 0.6%. Taurine and indicine proportion were not consistent across chromosomes. Furthermore, these non-uniform areas were found to be associated with traits that were probably under selection such as intermuscular fat and average daily gain. Interestingly, the sex chromosomes were predominantly taurine, which could be due to the composite being formed particularly in the final cross that resulted in progeny designated as purebred Brangus. This work demonstrated the process of new breed formation on a genomic level. It suggests that factors like genetic drift, selection and complementarity shift the genetic architecture into a uniquely different population. These findings are important to better understand how hybridization and crossbreeding systems shape the genetic architecture of composite populations.     

A genomic map of climate adaptation in Mediterranean cattle breeds

Domestic species such as cattle (Bos taurus taurus and B. t. indicus) represent attractive biological models to characterize the genetic basis of short‐term evolutionary response to climate pressure induced by their post‐domestication history. Here, using newly generated dense SNP genotyping data, we assessed the structuring of genetic diversity of 21 autochtonous cattle breeds from the whole Mediterranean basin and performed genome‐wide association analyses with covariables discriminating the different Mediterranean climate subtypes. This provided insights into both the demographic and adaptive histories of Mediterranean cattle. In particular, a detailed functional annotation of genes surrounding variants associated with climate variations highlighted several biological functions involved in Mediterranean climate adaptation such as thermotolerance, UV protection, pathogen resistance or metabolism with strong candidate genes identified (e.g., NDUFB3, FBN1, METTL3, LEF1, ANTXR2 and TCF7). Accordingly, our results suggest that main selective pressures affecting cattle in Mediterranean area may have been related to variation in heat and UV exposure, in food resources availability and in exposure to pathogens, such as anthrax bacteria (Bacillus anthracis). Furthermore, the observed contribution of the three main bovine ancestries (indicine, European and African taurine) in these different populations suggested that adaptation to local climate conditions may have either relied on standing genomic variation of taurine origin, or adaptive introgression from indicine origin, depending on the local breed origins. Taken together, our results highlight the genetic uniqueness of local Mediterranean cattle breeds and strongly support conservation of these populations.     

Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping

Background

Natural selection has molded evolution across all taxa. At an arguable date of around 330,000 years ago there were already at least two different types of cattle that became ancestors of nearly all modern cattle, the Bos taurus taurus more adapted to temperate climates and the tropically adapted Bos taurus indicus. After domestication, human selection exponentially intensified these differences. To better understand the genetic differences between these subspecies and detect genomic regions potentially under divergent selection, animals from the International Bovine HapMap Experiment were genotyped for over 770,000 SNP across the genome and compared using smoothed F_ST. The taurine sample was represented by ten breeds and the contrasting zebu cohort by three breeds.

Results

Each cattle group evidenced similar numbers of polymorphic markers well distributed across the genome. Principal components analyses and unsupervised clustering confirmed the well-characterized main division of domestic cattle. The top 1% smoothed F_ST, potentially associated to positive selection, contained 48 genomic regions across 17 chromosomes. Nearly half of the top F_ST signals (n = 22) were previously detected using a lower density SNP assay. Amongst the strongest signals were the BTA7:~50 Mb and BTA14:~25 Mb; both regions harboring candidate genes and different patterns of linkage disequilibrium that potentially represent intrinsic differences between cattle types. The bottom 1% of the smoothed F_ST values, potentially associated to balancing selection, included 24 regions across 13 chromosomes. These regions often overlap with copy number variants, including the highly variable region at BTA23:~24 Mb that harbors a large number of MHC genes. Under these regions, 318 unique Ensembl genes are annotated with a significant overrepresentation of immune related pathways.

Conclusions

Genomic regions that are potentially linked to purifying or balancing selection processes in domestic cattle were identified. These regions are of particular interest to understand the natural and human selective pressures to which these subspecies were exposed to and how the genetic background of these populations evolved in response to environmental challenges and human manipulation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-876) contains supplementary material, which is available to authorized users.     

平凉地方牛群体母系遗传背景

牛(Bos)的驯养作为农耕文化重要内容之一,在我国历史悠久,长久以来形成的役用型逐步被肉用型所替代.为了揭示平凉地方牛群体的遗传背景,分析其是否具有生产优质牛肉的遗传基础,本研究测定了88头平凉地方牛群体mtDNA D-loop HVS序列,对包括平凉地方牛群体在内的我国23个地方牛群体单倍型分布及系统发生关系进行了分析.结果表明,mtDNA D-loop高变区,在平凉地方牛群体共有52个单倍型,23个地方牛群体共有95个单倍型,这些单倍型在系统发生树和中介网络关系中分布于两个分枝,即瘤牛型和普通牛进化枝.因此,我们认为平凉地方牛群体和我国其他牛种一样,存在普通牛和瘤牛两个母系起源的遗传背景.     

Genetic diversity and genomic signatures of selection among cattle breeds from Siberia,eastern and northern Europe

Domestication in the near eastern region had a major impact on the gene pool of humpless taurine cattle (Bos taurus). As a result of subsequent natural and artificial selection, hundreds of different breeds have evolved, displaying a broad range of phenotypic traits. Here, 10 Eurasian B. taurus breeds from different biogeographic and production conditions, which exhibit different demographic histories and have been under artificial selection at various intensities, were investigated using the Illumina BovineSNP50 panel to understand their genetic diversity and population structure. In addition, we scanned genomes from eight breeds for signatures of diversifying selection. Our population structure analysis indicated six distinct breed groups, the most divergent being the Yakutian cattle from Siberia. Selection signals were shared (experimental P‐value < 0.01) with more than four breeds on chromosomes 6, 7, 13, 16 and 22. The strongest selection signals in the Yakutian cattle were found on chromosomes 7 and 21, where a miRNA gene and genes related to immune system processes are respectively located. In general, genomic regions indicating selection overlapped with known QTL associated with milk production (e.g. on chromosome 19), reproduction (e.g. on chromosome 24) and meat quality (e.g. on chromosome 7). The selection map created in this study shows that native cattle breeds and their genetic resources represent unique material for future breeding.     

Genetics and domestic cattle origins

Genetics has the potential to provide a novel layer of information pertaining to the origins and relationships of domestic cattle. While it is important not to overstate the power of archeological inference from genetic data, some previously widespread conjectures are inevitably contradicted with the addition of new information. Conjectures regarding domesticated cattle that fall into this category include a single domestication event with the development of Bos indicus breeds from earlier Bos taurus domesticates; the domestication of a third type of cattle in Africa having an intermediate morphology between the two taxa; and the special status of the Jersey breed as a European type with some exotic influences. In reality, a wideranging survey of the genetic variation of modern cattle reveals that they all derive from either zebu or taurine progenitors or are hybrids of the two. The quantitative divergence between Bos indicus and Bos taurus strongly supports a predomestic separation; that between African and European taurines also suggests genetic input from native aurochsen populations on each continent. Patterns of genetic variants assayed from paternally, maternally, and biparentally inherited genetic systems reveal that extensive hybridization of the two subspecies is part of the ancestry of Northern Indian, peripheral European, and almost all African cattle breeds. In Africa, which is the most extensive hybrid zone, the sexual asymmetry of the process of zebu introgression into native taurine breeds is strikingly evident. © 1998 Wiley-Liss, Inc.