Microsatellite DNA Variation and the Evolution, Domestication and Phylogeography of Taurine and Zebu Cattle (Bos Taurus and Bos Indicus)

Genetic variation at 20 microsatellite loci was surveyed to determine the evolutionary relationships and molecular biogeography of 20 different cattle populations from Africa, Europe and Asia. Phylogenetic reconstruction and multivariate analysis highlighted a marked distinction between humpless (taurine) and humped (zebu) cattle, providing strong support for a separate origin for domesticated zebu cattle. A molecular clock calculation using bison (Bison sp.) as an outgroup gave an estimated divergence time between the two subspecies of 610,000-850,000 years. Substantial differences in the distribution of alleles at 10 of these loci were observed between zebu and taurine cattle. These markers subsequently proved very useful for investigations of gene flow and admixture in African populations. When these data were considered in conjunction with previous mitochondrial and Y chromosomal studies, a distinctive male-mediated pattern of zebu genetic introgression was revealed. The introgression of zebu-specific alleles in African cattle afforded a high resolution perspective on the hybrid nature of African cattle populations and also suggested that certain West African populations of valuable disease-tolerant taurine cattle are under threat of genetic absorption by migrating zebu herds.     

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.     

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.     

Paternally inherited markers in bovine hybrid populations

The genetic integrity of crossfertile bovine- or cattle-like species may be endangered by species hybridization. Previously, amplified fragment length polymorphism, satellite fragment length polymorphism and microsatellite assays have been used to analyze the species composition of nuclear DNA in taurine cattle, zebu, banteng and bison populations, while mitochondrial DNA reveals the origin of the maternal lineages. Here, we describe species-specific markers of the paternally transmitted Y-chromosome for the direct detection of male-mediated introgression. Convenient PCR-restriction fragment length polymorphism and competitive PCR assays are shown to differentiate the Y-chromosomes of taurine cattle, American bison and European bison, and to detect the banteng origin of Indonesian Madura and Bali cattle bulls.     

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.     

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.     

Extensive female-mediated gene flow from sub-Saharan Africa into near eastern Arab populations

We have analyzed and compared mitochondrial DNA variation of populations from the Near East and Africa and found a very high frequency of African lineages present in the Yemen Hadramawt: more than a third were of clear sub-Saharan origin. Other Arab populations carried approximately 10% lineages of sub-Saharan origin, whereas non-Arab Near Eastern populations, by contrast, carried few or no such lineages, suggesting that gene flow has been preferentially into Arab populations. Several lines of evidence suggest that most of this gene flow probably occurred within the past approximately 2,500 years. In contrast, there is little evidence for male-mediated gene flow from sub-Saharan Africa in Y-chromosome haplotypes in Arab populations, including the Hadramawt. Taken together, these results are consistent with substantial migration from eastern Africa into Arabia, at least in part as a result of the Arab slave trade, and mainly female assimilation into the Arabian population as a result of miscegenation and manumission.     

Y-specific microsatellite polymorphisms in a range of bovid species

At least five dinucleotide (CA)n microsatellite repeat arrays have been assigned to the bovine Y-chromosome, with one marker (INRA124) shown to be polymorphic. We describe here the assessment of a panel of four Y-specific microsatellite markers for polymorphism in a range of cattle and related species. It was possible to amplify all the markers in the animals sampled and all showed variation. Three of the microsatellite loci (INRA124, INRA189 and BM861) displayed putative taurine- and zebu-specific alleles which can be useful indicators of male-mediated gene flow in hybrid populations. In the future these microsatellites, in combination with other Y-specific markers should provide a high-resolution Y haplotype system for evolutionary studies in both domesticated cattle and other related species.     

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.     

A microsatellite survey of cattle from a centre of origin: the Near East

Eight humpless cattle breeds from the Near East, three from Europe, one from West Africa and two zebu breeds from India were screened with 20 microsatellite loci. Breeds from the Near East revealed considerable levels of introgression from zebu cattle, which was apparent most in populations from the East and which declined in populations further West. This nonrandom pattern is suggestive of the introduction of zebu cattle from the East. Notwithstanding the overlay of zebu alleles, it was possible to demonstrate that Near Eastern cattle exhibited significantly higher levels of allelic diversity than breeds from other regions, which is consistent with the view that this region represents a primary domestication centre for Bos taurus cattle. The hypothesis that B. taurus and B. indicus cattle have separate domestic origins is also supported by the survey, a large genetic divergence being apparent between the nonhybrid taurine and zebu groups.     

Assessing introgression of Sahelian zebu genes into native Bos taurus breeds in Burkina Faso

A total of 350 samples were analyzed to estimate zebu gene proportions into two different taurine cattle breeds of Burkina Faso (Lobi and N’Dama) using 38 microsatellites and various statistical methodologies. West African and East African zebu samples were sequentially used as reference parental populations. Furthermore, N’Dama cattle from Congo, the composite South African Bonsmara cattle breed and a pool of European cattle were used successively as second parental populations. Independently of the methodology applied: (a) the use of West African zebu samples gave higher admixture coefficients than the East African zebu; (b) the higher zebu proportions were estimated when the European cattle was used as parental population 2; and (c) the use of the N’Dama population from Congo as parental population 2 gave the more consistent zebu proportion estimates for both the Lobi and the N’Dama breeds. In any case, the zebu admixture proportions estimated were not negligible and were always higher in the N’Dama cattle than in the Lobi cattle of Burkina Faso. This suggested that the introgression of Sahelian zebu genes into the taurine cattle of Southern West Africa can follow a complex pattern that can depend on local agro-ecological features. The current research pointed out that the estimation of admixture coefficients is highly dependent on both the assumptions underlying the methodologies applied and the selection of parental populations. Our analyses suggest that either too high or nil genetic identity between the parental and the expectedly derived populations must be avoided.     

High-resolution analysis of human Y-chromosome variation shows a sharp discontinuity and limited gene flow between northwestern Africa and the Iberian Peninsula

In the present study we have analyzed 44 Y-chromosome biallelic polymorphisms in population samples from northwestern (NW) Africa and the Iberian Peninsula, which allowed us to place each chromosome unequivocally in a phylogenetic tree based on >150 polymorphisms. The most striking results are that contemporary NW African and Iberian populations were found to have originated from distinctly different patrilineages and that the Strait of Gibraltar seems to have acted as a strong (although not complete) barrier to gene flow. In NW African populations, an Upper Paleolithic colonization that probably had its origin in eastern Africa contributed 75% of the current gene pool. In comparison, approximately 78% of contemporary Iberian Y chromosomes originated in an Upper Paleolithic expansion from western Asia, along the northern rim of the Mediterranean basin. Smaller contributions to these gene pools (constituting 13% of Y chromosomes in NW Africa and 10% of Y chromosomes in Iberia) came from the Middle East during the Neolithic and, during subsequent gene flow, from Sub-Saharan to NW Africa. Finally, bidirectional gene flow across the Strait of Gibraltar has been detected: the genetic contribution of European Y chromosomes to the NW African gene pool is estimated at 4%, and NW African populations may have contributed 7% of Iberian Y chromosomes. The Islamic rule of Spain, which began in a.d. 711 and lasted almost 8 centuries, left only a minor contribution to the current Iberian Y-chromosome pool. The high-resolution analysis of the Y chromosome allows us to separate successive migratory components and to precisely quantify each historical layer.     

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.     

Y-chromosome lineages in Cabo Verde Islands witness the diverse geographic origin of its first male settlers

The Y-chromosome haplogroup composition of the population of the Cabo Verde Archipelago was profiled by using 32 single-nucleotide polymorphism markers and compared with potential source populations from Iberia, west Africa, and the Middle East. According to the traditional view, the major proportion of the founding population of Cabo Verde was of west African ancestry with the addition of a minor fraction of male colonizers from Europe. Unexpectedly, more than half of the paternal lineages (53.5%) of Cabo Verdeans clustered in haplogroups I, J, K, and R1, which are characteristic of populations of Europe and the Middle East, while being absent in the probable west African source population of Guiné-Bissau. Moreover, a high frequency of J* lineages in Cabo Verdeans relates them more closely to populations of the Middle East and probably provides the first genetic evidence of the legacy of the Jews. In addition, the considerable proportion (20.5%) of E3b(xM81) lineages indicates a possible gene flow from the Middle East or northeast Africa, which, at least partly, could be ascribed to the Sephardic Jews. In contrast to the predominance of west African mitochondrial DNA haplotypes in their maternal gene pool, the major west African Y-chromosome lineage E3a was observed only at a frequency of 15.9%. Overall, these results indicate that gene flow from multiple sources and various sex-specific patterns have been important in the formation of the genomic diversity in the Cabo Verde islands.An erratum to this article can be found at     

Microsatellite analysis reveals high genetic diversity but low genetic structure in Ethiopian indigenous cattle populations

Ethiopian cattle are under threat from uncontrolled mating practices and are at high risk of becoming genetically homogeneous. Therefore, to evaluate genetic diversity, population structure and degree of admixture, 30 microsatellite markers were genotyped using 351 DNA samples from 10 Ethiopian cattle populations and the Holstein breed. The mean number of alleles per cattle population ranged from 6.93 ± 2.12 in Sheko to 7.50 ± 2.35 in Adwa. The mean observed and expected heterozygosities were 0.674 ± 0.015 and 0.726 ± 0.019 respectively. Ethiopian cattle populations have maintained a high level of within-population genetic differentiation (98.7%), the remainder being accounted for by differentiation among populations (1.3%). A highly significant deficiency in heterozygotes was detected within populations ( F _IS = 0.071; P  <   0.001) and total inbreeding ( F _IT = 0.083; P  <   0.001). The study populations were highly admixed but distinct from pure Bos taurus and Bos indicus breeds. The various levels of admixture and high genetic diversity make Ethiopian cattle populations suitable for future genetic improvement and utilization under a wide range of agro-ecologies in Ethiopia.     

Allozyme variation in 16 natural populations of Faidherbia albida (Del.) A. Chev

The objectives of this study were (1) to investigate the genetic variability; and (2) to determine the phylogenetic relationship in the leguminous tree Faidherbia albida through out its range of distribution in East, West and Southern Africa. A total of 16 populations were subjected to enzyme electrophoresis and 6 enzyme systems (Adh, Mdh, G6pdh, Idh, 6Pgdh, and alpha-Est) encoded by 12 loci were scored. A total of fifty one different alleles were detected, with an average of 2.5 alleles per locus. Forty three percent of the loci were polymorphic at a 95% criterion. The average expected heterozygosity (gene diversity index H(e)) was 0.141. All provenances showed significant deviation from Hardy-Weinberg expectation. The UPGMA cluster analysis, based on Modified Rogers distance revealed close similarities between Eastern and Southern African provenances, except Debre zeit from Ethiopia, which was closest to West African populations than to the East African populations. Also, Bignona from West Africa was peculiarly closer to East African provenances. Differentiation of West African populations from Eastern and Southern African provenances was quite evident, emphasising the clinical pattern of variation in these two major geographical races of F. albida in Africa. Wright's F-statistics showed an overall significant deficit of heterozygotes, a common feature in mixed mating, entomophilous, widespread species such as F. albida. The dendrogram analysis showed wide separation among the three Ethiopian provenances indicating a high level of genetic differentiation and diversity among them.     

High-resolution predictive mapping for Rhipicephalus appendiculatus (Acari: Ixodidae) in the Horn of Africa

The brown ear tick Rhipicephalus appendiculatus, vector of East Coast fever (ECF) and related cattle diseases caused by Theileria parva has never been reported from the Horn of Africa. Habitat suitability for this tick species was predicted using Maxent modelling technique based on R. appendiculatus records in Sub-Saharan Africa. Two models were developed: the first is based on the tropical R. appendiculatus distribution and the one is based on the distribution records in the temperate region of Sub-Saharan Africa. The tropical model shows favourable habitat in much of the Ethiopian highlands. The whole Djibouti, the south eastern Ethiopian lowlands, majority of Somalia and Eritrea were found to be not suitable for the survival and development of this tick species. Highly suitable areas occur in areas which have moderate temperature and high precipitation. Introductions of R. appendiculatus into the Horn of Africa probably have been prevented by the natural barrier between the known R. appendiculatus distribution range in East Africa and the Horn of Africa. The effect of an introduction of R. appendiculatus and thereby ECF into the Horn of Africa could be catastrophic since the cattle in this area have no immunity against ECF, and mortality might be considerable in all age groups of cattle.     

Genetic Structure and Differentiation of Three Chinese Indigenous Cattle Populations

Levels of genetic differentiation, gene flow, and genetic structure of three indigenous cattle populations (Luxi, Bohai, and Minnan) and two reference cattle populations (Chinese Holstein and Qinhai yak) in China were estimated using the information from 12 microsatellites, and 141 microsatellite alleles were identified. The mean number of alleles per locus ranged from 2.9005 in yak to 4.9722 in Holstein. The observed heterozygosity ranged from 0.5325 (yak) to 0.7719 (Holstein); 29 private alleles were detected. The global heterozygote deficit across all populations amounted to 58.5% (p < 0.001). The overall significant (p < 0.001) deficit of heterozygotes because of inbreeding within breeds amounted to 43.2%. The five cattle populations were highly differentiated (F _st = 26.9%, p < 0.001) at all loci. The heterozygote deficit within the population was highest in Luxi cattle and lowest in yak. The average number of effective migrants exchanged per generation was highest (1.149) between Luxi and Holstein, and lowest (0.509) between Luxi and yak. With the application of prior population information, cluster analysis achieved posterior probabilities from 91% to 98% of correctly assigning individuals to populations. Combining the information of cluster analysis, gene flow, and Structure analysis, the five cattle populations belong to three genetic clusters, a taurine (Luxi and Chinese Holstein), a zebu (Bohai and Minnan), and a yak cluster. This indicates that Bohai black is closer to Bos indicus than Luxi cattle. The evolution and development of three indigenous cattle populations are discussed.     

The distribution of human genetic diversity: a comparison of mitochondrial, autosomal, and Y-chromosome data

We report a comparison of worldwide genetic variation among 255 individuals by using autosomal, mitochondrial, and Y-chromosome polymorphisms. Variation is assessed by use of 30 autosomal restriction-site polymorphisms (RSPs), 60 autosomal short-tandem-repeat polymorphisms (STRPs), 13 Alu-insertion polymorphisms and one LINE-1 element, 611 bp of mitochondrial control-region sequence, and 10 Y-chromosome polymorphisms. Analysis of these data reveals substantial congruity among this diverse array of genetic systems. With the exception of the autosomal RSPs, in which an ascertainment bias exists, all systems show greater gene diversity in Africans than in either Europeans or Asians. Africans also have the largest total number of alleles, as well as the largest number of unique alleles, for most systems. GST values are 11%-18% for the autosomal systems and are two to three times higher for the mtDNA sequence and Y-chromosome RSPs. This difference is expected because of the lower effective population size of mtDNA and Y chromosomes. A lower value is seen for Y-chromosome STRs, reflecting a relative lack of continental population structure, as a result of rapid mutation and genetic drift. Africa has higher GST values than does either Europe or Asia for all systems except the Y-chromosome STRs and Alus. All systems except the Y-chromosome STRs show less variation between populations within continents than between continents. These results are reassuring in their consistency and offer broad support for an African origin of modern human populations.     

Phylogeography of mitochondrial DNA and Y-chromosome haplogroups reveal asymmetric gene flow in populations of Eastern India

Polymorphisms in mitochondrial (mt) DNA and Y-chromosomes of seven socially and linguistically diverse castes and tribes of Eastern India were examined to determine their genetic relationships, their origin, and the influence of demographic factors on population structure. Samples from the Orissa Brahmin, Karan, Khandayat, Gope, Juang, Saora, and Paroja were analyzed for mtDNA hypervariable sequence (HVS) I and II, eight Y-chromosome short tandem repeats (Y-STRs), and lineage-defining mutations diagnostic for Indian- and Eurasian-specific haplogroups. Our results reveal that haplotype diversity and mean pairwise differences (MPD) was higher in caste groups of the region (>0.998, for both systems) compared to tribes (0.917-0.996 for Y-STRs, and 0.958-0.988 for mtDNA haplotypes). The majority of paternal lineages belong to the R1a1, O2a, and H haplogroups (62.7%), while 73.2% of maternal lineages comprise the Indian-specific M*, M5, M30, and R* mtDNA haplogroups, with a sporadic occurrence of West Eurasian lineages. Our study reveals that Orissa Brahmins (a higher caste population) have a genetic affinity with Indo-European speakers of Eastern Europe, although the Y-chromosome data show that the genetic distances of populations are not correlated to their position in the caste hierarchy. The high frequency of the O2a haplogroup and absence of East Asian-specific mtDNA lineages in the Juang and Saora suggest that a migration of Austro-Asiatic tribes to mainland India was exclusively male-mediated which occurred during the demographic expansion of Neolithic farmers in southern China. The phylogeographic analysis of mtDNA and Y-chromosomes revealed varied ancestral sources for the diverse genetic components of the populations of Eastern India.