African matrilineages in American Creole cattle: evidence of two independent continental sources

Phenotypic Characterization of Creole Cattle in the Andean Highlands Using Bio-Morphometric Measures and Zoometric Indices

Several Creole cattle biotypes can be found in the Andean highlands, and most of them are considered as being in risk of extinction. The main aim of the present study was to perform a phenotypic characterization of the Creole cattle in the Andean highlands using bio-morphometric measures and zoometric indices. Individuals from three different biotypes (Black 'Negro' (n = 57), Colour-Sided 'Callejón' (n = 20), and Brindle 'Atigrado' (n = 18)) from an experimental research center located in the Peruvian highlands were enrolled in the study. In total, seventeen morphometric parameters were evaluated and ten zoometric indices were calculated in each biotype. To test the relationship between biometric traits, correlation analyses were carried out between morphometric parameters. Differences were observed regarding different morphometric variables such as head length (HL) and rump length (RL) among cattle biotypes (p ≤ 0.05). The coefficient of variation (CV; %) regarding different morphometric parameters ranged between 11.32 for neck length (NL) and 3.63 for height at withers (HaW), which indicated low-moderate variability among morphometric variables. Differences were observed in the longitudinal pelvic index (LPI) when different zoometric indices were compared among biotypes (p ≤ 0.05). The CV regarding different zoometric indices, which ranged between 10.78 for the cephalic index (CEI) and 5.05 for LPI, indicated low variability among indices. No differences were observed in any other morphometric parameter or zoometric index among cattle biotypes or genders (p > 0.05). Finally, multiple correlations were observed between morphometric variables (p ≤ 0.05). In conclusion, it was determined that Peruvian Andean Creole cattle can be considered as a dairy-related biotype with a slight tendency for beef production (dual-purpose). The great homogeneity regarding zoometric characteristics among biotypes and genders may indicate that the Andean Creole cattle have been maintained quite isolated, avoiding the genetic influence of other foreign breeds. Finally, the phenotypic characterization including bio-morphometric measurements and zoometric indices obtained from the different Creole bovine biotypes is crucial in order to begin different conservation programs to preserve cattle breeds in the Peruvian Andean highlands.

The Complete Mitochondrial Genome of a Neglected Breed, the Peruvian Creole Cattle (Bos taurus), and Its Phylogenetic Analysis

Cattle spread throughout the American continent during the colonization years, originating creole breeds that adapted to a wide range of climate conditions. The population of creole cattle in Peru is decreasing mainly due to the introduction of more productive breeds in recent years. During the last 15 years, there has been significant progress in cattle genomics. However, little is known about the genetics of the Peruvian creole cattle (PCC) despite its importance to (i) improving productivity in the Andean region, (ii) agricultural labor, and (iii) cultural traditions. In addition, the origin and phylogenetic relationship of the PCC are still unclear. In order to promote the conservation of the PCC, we sequenced the mitochondrial genome of a creole bull, which also possessed exceptional fighting skills and was employed for agricultural tasks, from the highlands of Arequipa for the first time. The total mitochondrial genome sequence is 16,339 bp in length with the base composition of 31.43% A, 28.64% T, 26.81% C, and 13.12% G. It contains 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. Among the 37 genes, 28 were positioned on the H-strand and 9 were positioned on the L-strand. The most frequently used codons were CUA (leucine), AUA (isoleucine), AUU (isoleucine), AUC (isoleucine), and ACA (threonine). Maximum likelihood reconstruction using complete mitochondrial genome sequences showed that the PCC is related to native African breeds. The annotated mitochondrial genome of PCC will serve as an important genetic data set for further breeding work and conservation strategies.

Recovery of mitogenomes from whole genome sequences to infer maternal diversity in 1883 modern taurine and indicine cattle

Maternal diversity based on a sub-region of mitochondrial genome or variants were commonly used to understand past demographic events in livestock. Additionally, there is growing evidence of direct association of mitochondrial genetic variants with a range of phenotypes. Therefore, this study used complete bovine mitogenomes from a large sequence database to explore the full spectrum of maternal diversity. Mitogenome diversity was evaluated among 1883 animals representing 156 globally important cattle breeds. Overall, the mitogenomes were diverse: presenting 11 major haplogroups, expanding to 1309 unique haplotypes, with nucleotide diversity 0.011 and haplotype diversity 0.999. A small proportion of African taurine (3.5%) and indicine (1.3%) haplogroups were found among the European taurine breeds and composites. The haplogrouping was largely consistent with the population structure derived from alternate clustering methods (e.g. PCA and hierarchical clustering). Further, we present evidence confirming a new indicine subgroup (I1a, 64 animals) mainly consisting of breeds originating from China and characterised by two private mutations within the I1 haplogroup. The total genetic variation was attributed mainly to within-breed variance (96.9%). The accuracy of the imputation of missing genotypes was high (99.8%) except for the relatively rare heteroplasmic genotypes, suggesting the potential for trait association studies within a breed.

The genetic ancestry of American Creole cattle inferred from uniparental and autosomal genetic markers

Cattle imported from the Iberian Peninsula spread throughout America in the early years of discovery and colonization to originate Creole breeds, which adapted to a wide diversity of environments and later received influences from other origins, including zebu cattle in more recent years. We analyzed uniparental genetic markers and autosomal microsatellites in DNA samples from 114 cattle breeds distributed worldwide, including 40 Creole breeds representing the whole American continent, and samples from the Iberian Peninsula, British islands, Continental Europe, Africa and American zebu. We show that Creole breeds differ considerably from each other, and most have their own identity or group with others from neighboring regions. Results with mtDNA indicate that T1c-lineages are rare in Iberia but common in Africa and are well represented in Creoles from Brazil and Colombia, lending support to a direct African influence on Creoles. This is reinforced by the sharing of a unique Y-haplotype between cattle from Mozambique and Creoles from Argentina. Autosomal microsatellites indicate that Creoles occupy an intermediate position between African and European breeds, and some Creoles show a clear Iberian signature. Our results confirm the mixed ancestry of American Creole cattle and the role that African cattle have played in their development.

Genetic diversity of BoLA-DRB3 in South American Zebu cattle populations

Background

Bovine leukocyte antigens (BoLAs) are used extensively as markers of disease and immunological traits in cattle. However, until now, characterization of BoLA gene polymorphisms in Zebu breeds using high resolution typing methods has been poor. Here, we used a polymerase chain reaction sequence-based typing (PCR-SBT) method to sequence exon 2 of the BoLA class II DRB3 gene from 421 cattle (116 Bolivian Nellore, 110 Bolivian Gir, and 195 Peruvian Nellore-Brahman). Data from 1416 Taurine and Zebu samples were also included in the analysis.

Results

We identified 46 previously reported alleles and no novel variants. Of note, 1/3 of the alleles were detected only in Zebu cattle. Comparison of the degree of genetic variability at the population and sequence levels with genetic distance in the three above mentioned breeds and nine previously reported breeds revealed that Zebu breeds had a gene diversity score higher than 0.86, a nucleotide diversity score higher than 0.06, and a mean number of pairwise differences greater than 16, being similar to those estimated for other cattle breeds. A neutrality test revealed that only Nellore-Brahman cattle showed the even gene frequency distribution expected under a balanced selection scenario. The F_ST index and the exact G test showed significant differences across all cattle populations (F_ST = 0.057; p <  0.001). Neighbor-joining trees and principal component analysis identified two major clusters: one comprising mainly European Taurine breeds and a second comprising Zebu breeds. This is consistent with the historical and geographical origin of these breeds. Some of these differences may be explained by variation of amino acid motifs at antigen-binding sites.

Conclusions

The results presented herein show that the historical divergence between Taurine and Zebu cattle breeds is a result of origin, selection, and adaptation events, which would explain the observed differences in BoLA-DRB3 gene diversity between the two major bovine types. This allelic information will be important for investigating the relationship between the major histocompatibility complex and disease, and contribute to an ongoing effort to catalog bovine MHC allele frequencies according to breed and location.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0618-7) contains supplementary material, which is available to authorized users.

Genetic diversity, population structure, and correlations between locally adapted zebu and taurine breeds in Brazil using SNP markers

Genetic diversity is one of the most important issues in studies on conservation of cattle breeds and endangered species. The objective of this study was to estimate the levels of genetic differentiation between locally adapted taurine (Bos taurus taurus) and zebu (Bos taurus indicus) breeds in Brazil, which were genotyped for more than 777,000 SNPs. The fixation index (F _ST), principal component analysis (PCA), and Bayesian clustering were estimated. The F _ST highlighted genetic differentiation between taurine and zebu breeds. The taurine lines, Caracu and Caracu Caldeano, had significant genetic differentiation (F _ST close to 5%) despite their recent selection for different uses (meat and milk). This genetic variability can be used for conservation of locally adapted animals, as well as for breeding programs on zebu breeds. Introgression of zebu in locally adapted breeds was identified, especially in Curraleiro Pé-Duro breed. The Gyr breed, however, had low breed purity at genomic level due to its very heterogeneous mixing pattern.

Assessment of biodiversity in Chilean cattle using the distribution of major histocompatibility complex class II BoLA-DRB3 allele

Bovine leukocyte antigens (BoLAs) are used extensively as markers for bovine disease and immunological traits. In this study, we estimated BoLA-DRB3 allele frequencies using 888 cattle from 10 groups, including seven cattle breeds and three crossbreeds: 99 Red Angus, 100 Black Angus, 81 Chilean Wagyu, 49 Hereford, 95 Hereford × Angus, 71 Hereford × Jersey, 20 Hereford × Overo Colorado, 113 Holstein, 136 Overo Colorado, and 124 Overo Negro cattle. Forty-six BoLA-DRB3 alleles were identified, and each group had between 12 and 29 different BoLA-DRB3 alleles. Overo Negro had the highest number of alleles (29); this breed is considered in Chile to be an 'Old type' European Holstein Friesian descendant. By contrast, we detected 21 alleles in Holstein cattle, which are considered to be a 'Present type' Holstein Friesian cattle. Chilean cattle groups and four Japanese breeds were compared by neighbor-joining trees and a principal component analysis (PCA). The phylogenetic tree showed that Red Angus and Black Angus cattle were in the same clade, crossbreeds were closely related to their parent breeds, and Holstein cattle from Chile were closely related to Holstein cattle in Japan. Overall, the tree provided a thorough description of breed history. It also showed that the Overo Negro breed was closely related to the Holstein breed, consistent with historical data indicating that Overo Negro is an 'Old type' Holstein Friesian cattle. This allelic information will be important for investigating the relationship between major histocompatibility complex (MHC) and disease.

Worldwide patterns of ancestry, divergence, and admixture in domesticated cattle

The domestication and development of cattle has considerably impacted human societies, but the histories of cattle breeds and populations have been poorly understood especially for African, Asian, and American breeds. Using genotypes from 43,043 autosomal single nucleotide polymorphism markers scored in 1,543 animals, we evaluate the population structure of 134 domesticated bovid breeds. Regardless of the analytical method or sample subset, the three major groups of Asian indicine, Eurasian taurine, and African taurine were consistently observed. Patterns of geographic dispersal resulting from co-migration with humans and exportation are recognizable in phylogenetic networks. All analytical methods reveal patterns of hybridization which occurred after divergence. Using 19 breeds, we map the cline of indicine introgression into Africa. We infer that African taurine possess a large portion of wild African auroch ancestry, causing their divergence from Eurasian taurine. We detect exportation patterns in Asia and identify a cline of Eurasian taurine/indicine hybridization in Asia. We also identify the influence of species other than Bos taurus taurus and B. t. indicus in the formation of Asian breeds. We detect the pronounced influence of Shorthorn cattle in the formation of European breeds. Iberian and Italian cattle possess introgression from African taurine. American Criollo cattle originate from Iberia, and not directly from Africa with African ancestry inherited via Iberian ancestors. Indicine introgression into American cattle occurred in the Americas, and not Europe. We argue that cattle migration, movement and trading followed by admixture have been important forces in shaping modern bovine genomic variation.

The DNA of domesticated plants and animals contains information about how species were domesticated, exported, and bred by early farmers. Modern breeds were developed by lengthy and complex processes; however, our use of 134 breeds and new analytical models enabled us to reveal some of the processes that created modern cattle diversity. In Asia, Africa, North and South America, humpless (Bos t. taurus or taurine) and humped (Bos t. indicus or indicine) cattle were crossbred to produce hybrids adapted to the environment and local production systems. The history of Asian cattle involves the domestication and admixture of several species whereas African taurines arose through the introduction of domesticated Fertile Crescent taurines and their hybridization with wild African aurochs. African taurine genetic background is commonly observed among European Mediterranean breeds. The absence of indicine introgression within most European taurine breeds, but presence within three Italian breeds is consistent with at least two separate migration waves of cattle to Europe, one from the Middle East which captured taurines in which indicine introgression had already occurred and the second from western Africa into Spain with no indicine introgression. This second group seems to have radiated from Spain into the Mediterranean resulting in a cline of African taurine introgression into European taurines.

Ancient mtDNA analysis of early 16(th) century Caribbean cattle provides insight into founding populations of New World creole cattle breeds

The Columbian Exchange resulted in a widespread movement of humans, plants and animals between the Old and New Worlds. The late 15^th to early 16^th century transfer of cattle from the Iberian Peninsula and Canary Islands to the Caribbean laid the foundation for the development of American creole cattle (Bos taurus) breeds. Genetic analyses of modern cattle from the Americas reveal a mixed ancestry of European, African and Indian origins. Recent debate in the genetic literature centers on the ‘African’ haplogroup T1 and its subhaplogroups, alternatively tying their origins to the initial Spanish herds, and/or from subsequent movements of taurine cattle through the African slave trade. We examine this problem through ancient DNA analysis of early 16^th century cattle bone from Sevilla la Nueva, the first Spanish colony in Jamaica. In spite of poor DNA preservation, both T3 and T1 haplogroups were identified in the cattle remains, confirming the presence of T1 in the earliest Spanish herds. The absence, however, of “African-derived American” haplotypes (AA/T1c1a1) in the Sevilla la Nueva sample, leaves open the origins of this sub-haplogroup in contemporary Caribbean cattle.

Genetic footprints of Iberian cattle in America 500 years after the arrival of Columbus

Background

American Creole cattle presumably descend from animals imported from the Iberian Peninsula during the period of colonization and settlement, through different migration routes, and may have also suffered the influence of cattle directly imported from Africa. The introduction of European cattle, which began in the 18th century, and later of Zebu from India, has threatened the survival of Creole populations, some of which have nearly disappeared or were admixed with exotic breeds. Assessment of the genetic status of Creole cattle is essential for the establishment of conservation programs of these historical resources.

Methodology/Principal Findings

We sampled 27 Creole populations, 39 Iberian, 9 European and 6 Zebu breeds. We used microsatellite markers to assess the origins of Creole cattle, and to investigate the influence of different breeds on their genetic make-up. The major ancestral contributions are from breeds of southern Spain and Portugal, in agreement with the historical ports of departure of ships sailing towards the Western Hemisphere. This Iberian contribution to Creoles may also include some African influence, given the influential role that African cattle have had in the development of Iberian breeds, but the possibility of a direct influence on Creoles of African cattle imported to America can not be discarded. In addition to the Iberian influence, the admixture with other European breeds was minor. The Creoles from tropical areas, especially those from the Caribbean, show clear signs of admixture with Zebu.

Conclusions/Significance

Nearly five centuries since cattle were first brought to the Americas, Creoles still show a strong and predominant signature of their Iberian ancestors. Creole breeds differ widely from each other, both in genetic structure and influences from other breeds. Efforts are needed to avoid their extinction or further genetic erosion, which would compromise centuries of selective adaptation to a wide range of environmental conditions.

Origin and spread of Bos taurus: new clues from mitochondrial genomes belonging to haplogroup T1

Background

Most genetic studies on modern cattle have established a common origin for all taurine breeds in the Near East, during the Neolithic transition about 10 thousand years (ka) ago. Yet, the possibility of independent and/or secondary domestication events is still debated and is fostered by the finding of rare mitochondrial DNA (mtDNA) haplogroups like P, Q and R. Haplogroup T1, because of its geographic distribution, has been the subject of several investigations pointing to a possible independent domestication event in Africa and suggesting a genetic contribution of African cattle to the formation of Iberian and Creole cattle. Whole mitochondrial genome sequence analysis, with its proven effectiveness in improving the resolution of phylogeographic studies, is the most appropriate tool to investigate the origin and structure of haplogroup T1.

Methodology

A survey of >2200 bovine mtDNA control regions representing 28 breeds (15 European, 10 African, 3 American) identified 281 subjects belonging to haplogroup T1. Fifty-four were selected for whole mtDNA genome sequencing, and combined with ten T1 complete sequences from previous studies into the most detailed T1 phylogenetic tree available to date.

Conclusions

Phylogenetic analysis of the 64 T1 mitochondrial complete genomes revealed six distinct sub-haplogroups (T1a–T1f). Our data support the overall scenario of a Near Eastern origin of the T1 sub-haplogroups from as much as eight founding T1 haplotypes. However, the possibility that one sub-haplogroup (T1d) arose in North Africa, in domesticated stocks, shortly after their arrival from the Near East, can not be ruled out. Finally, the previously identified “African-derived American" (AA) haplotype turned out to be a sub-clade of T1c (T1c1a1). This haplotype was found here for the first time in Africa (Egypt), indicating that it probably originated in North Africa, reached the Iberian Peninsula and sailed to America, with the first European settlers.

Assessing priorities for conservation in Tuscan cattle breeds using microsatellites

Preservation of rare genetic stocks requires assessment of within-population genetic diversity and between-population differentiation to make inferences on their degree of uniqueness. A total of 194 Tuscan cattle (44 Calvana, 35 Chianina, 25 Garfagnina, 31 Maremmana, 31 Mucca Pisana and 28 Pontremolese) individuals were genotyped for 34 microsatellite markers. Moreover, 56 samples belonging to Argentinean Creole and Asturiana de la Montaña cattle breeds were used as an outgroup. Genetic diversity was quantified in terms of molecular coancestry and allelic richness. STRUCTURE analyses showed that the Tuscan breeds have well-differentiated genetic backgrounds, except for the Calvana and Chianina breeds, which share the same genetic ancestry. The between-breed Nei's minimum distance (Dm) matrices showed that the pair Calvana-Chianina was less differentiated (0.049 ± 0.006). The endangered Tuscan breeds (Calvana, Garfagnina, Mucca Pisana and Pontremolese) made null or negative contributions to diversity, except for the Mucca Pisana contribution to allelic richness (CT = 1.8%). The Calvana breed made null or negative within-breed contributions (W = 0.0%; CW = -0.4%). The Garfagnina and Pontremolese breeds made positive contributions to between-breed diversity but negative and high within-breed contributions, thus suggesting population bottleneck with allelic losses and increase of homozygosity in the population. Exclusion of the four endangered Tuscan cattle breeds did not result in losses in genetic diversity (T = -0.7%; CT = -1.2%), whereas exclusion of the non-endangered breeds (Chianina and Maremmana) did (T = 2.1%; CT = 3.9%); the simple exclusion of the Calvana breed from the former group led to losses in genetic diversity (T = 0.47%; CT = 2.34%), indicating a diverse significance for this breed. We showed how quantifying both within-population diversity and between-population differentiation in terms of allelic frequencies and allelic richness provides different and complementary information on the genetic backgrounds assessed and may help to implement priorities and strategies for conservation in livestock.

Genetic characterization of Latin-American Creole cattle using microsatellite markers

Genetic diversity in and relationships among 26 Creole cattle breeds from 10 American countries were assessed using 19 microsatellites. Heterozygosities, F-statistics estimates, genetic distances, multivariate analyses and assignment tests were performed. The levels of within-breed diversity detected in Creole cattle were considerable and higher than those previously reported for European breeds, but similar to those found in other Latin American breeds. Differences among breeds accounted for 8.4% of the total genetic variability. Most breeds clustered separately when the number of pre-defined populations was 21 (the most probable K value), with the exception of some closely related breeds that shared the same cluster and others that were admixed. Despite the high genetic diversity detected, significant inbreeding was also observed within some breeds, and heterozygote excess was detected in others. These results indicate that Creoles represent important reservoirs of cattle genetic diversity and that appropriate conservation measures should be implemented for these native breeds in order to minimize inbreeding and uncontrolled crossbreeding.

A new insight into cattle's maternal origin in six Asian countries

The domestication of cattle fuelled the development of agricultural society in the history of human being. The evolution and genetic relationship of cattle can be elucidated by investigating the variation of mitochondrial DNA (mtDNA) D-loop sequence. In this study, we built a cattle phylogeny with a pool of 856 individual D-loop sequences, of which 264 Chinese cattle D-loop sequences were obtained in this study (141 ones were first analyzed, and 123 were first submitted) and the rest sequences of cattle from six Asian countries (Japan, Korea, Mongolia, Nepal, India and China) were retrieved from GenBank. Our results indicated that cattle from six Asian countries fell into three clades, Bos taurus (taurine), Bos indicus (zebu) and yak. Four main haplogroups T1A, T2, T3 (including T3A and T3B) and T5 were found in taurine, and two haplogroups I1 and I2 in zebu. Furthermore, we found that I1 and I2 haplogroups were separated by four variable sites rather than five ones and four haplogroups or sub-haplogroups of T1A, T3A, T3B and T5 were found for the first time in these Asian cattle. These data brought us a new insight into cattle's genetic structure in these six Asian countries. The geographical distribution of haplogroups was also outlined to provide systematic information on cattle genetic resources.

Genetic diversity in farm animals--a review

Domestication of livestock species and a long history of migrations, selection and adaptation have created an enormous variety of breeds. Conservation of these genetic resources relies on demographic characterization, recording of production environments and effective data management. In addition, molecular genetic studies allow a comparison of genetic diversity within and across breeds and a reconstruction of the history of breeds and ancestral populations. This has been summarized for cattle, yak, water buffalo, sheep, goats, camelids, pigs, horses, and chickens. Further progress is expected to benefit from advances in molecular technology.

Origins and genetic diversity of New World Creole cattle: inferences from mitochondrial and Y chromosome polymorphisms

The ancestry of New World cattle was investigated through the analysis of mitochondrial and Y chromosome variation in Creoles from Argentina, Brazil, Mexico, Paraguay and the United States of America. Breeds that influenced the Creoles, such as Iberian native, British and Zebu, were also studied. Creoles showed high mtDNA diversity (H = 0.984 +/- 0.003) with a total of 78 haplotypes, and the European T3 matriline was the most common (72.1%). The African T1a haplogroup was detected (14.6%), as well as the ancestral African-derived AA matriline (11.9%), which was absent in the Iberian breeds. Genetic proximity among Creoles, Iberian and Atlantic Islands breeds was inferred through their sharing of mtDNA haplotypes. Y-haplotype diversity in Creoles was high (H = 0.779 +/- 0.019), with several Y1, Y2 and Y3 haplotypes represented. Iberian patrilines in Creoles were more difficult to infer and were reflected by the presence of H3Y1 and H6Y2. Y-haplotypes confirmed crossbreeding with British cattle, mainly of Hereford with Pampa Chaqueño and Texas Longhorn. Male-mediated Bos indicus introgression into Creoles was found in all populations, except Argentino1 (herd book registered) and Pampa Chaqueño. The detection of the distinct H22Y3 patriline with the INRA189-90 allele in Caracú suggests introduction of bulls directly from West Africa. Further studies of Spanish and African breeds are necessary to elucidate the origins of Creole cattle, and determine the exact source of their African lineages.

Ancestral matrilineages and mitochondrial DNA diversity of the Lidia cattle breed

To clarify the genetic ancestry and the mitochondrial DNA (mtDNA) diversity of the Lidia cattle breed, a 521-bp D-loop fragment was sequenced in 527 animals belonging to 70 herds distributed across 29 lineages. The mtDNA diversity recorded was similar to that seen for Middle Eastern breeds and greater than that recorded for the majority of European breeds. Haplotype T3 was the most common (81%), followed by the African T1 haplotype (17%); very low frequencies were recorded for haplotypes T and T2. The results agree with there being two major ancestral lines for the Lidia breed, European and African, similar to that seen for other Mediterranean breeds. A wide range of variation in haplotype frequencies was seen between the examined lineages. Haplotype T3 was present in all those analysed; in five it was the only one present, and in only one lineage (Miura) was its frequency lower than that of T1. T1*, a haplotype reported in Criollo breeds and to date in only a single European breed (the Retinta breed from Spain), was found in a single animal belonging to the Concha y Sierra lineage. Network analysis of the Lidia breed revealed the presence of two major haplotypes: T3 and T1. The Lidia breed appears to be more closely related to prehistoric Iberian and Italian than to British aurochs.