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

Ancient DNA confirms diverse origins of early post-Columbian cattle in the Americas

Before the arrival of Europeans, domestic cattle (Bos taurus) did not exist in the Americas, and most of our knowledge about how domestic bovines first arrived in the Western Hemisphere is based on historical documents. Sixteenth-century colonial accounts suggest that the first cattle were brought in small numbers from the southern Iberian Peninsula via the Canary archipelago to the Caribbean islands where they were bred locally and imported to other circum-Caribbean regions. Modern American heritage cattle genetics and limited ancient mtDNA data from archaeological colonial cattle suggest a more complex story of mixed ancestries from Europe and Africa. So far little information exists to understand the nature and timing of the arrival of these mixed-ancestry populations. In this study we combine ancient mitochondrial and nuclear DNA from a robust sample of some of the earliest archaeological specimens from Caribbean and Mesoamerican sites to clarify the origins and the dynamics of bovine introduction into the Americas. Our analyses support first arrival of cattle from diverse locales and potentially confirm the early arrival of African-sourced cattle in the Americas, followed by waves of later introductions from various sources over several centuries.

Genetic Structure Analysis of 155 Transboundary and Local Populations of Cattle (Bos taurus, Bos indicus and Bos grunniens) Based on STR Markers

Every week, 1-2 breeds of farm animals, including local cattle, disappear in the world. As the keepers of rare allelic variants, native breeds potentially expand the range of genetic solutions to possible problems of the future, which means that the study of the genetic structure of these breeds is an urgent task. Providing nomadic herders with valuable resources necessary for life, domestic yaks have also become an important object of study. In order to determine the population genetic characteristics, and clarify the phylogenetic relationships of modern representatives of 155 cattle populations from different regions of the world, we collected a large set of STR data (10,250 individuals), including unique native cattle, 12 yak populations from Russia, Mongolia and Kyrgyzstan, as well as zebu breeds. Estimation of main population genetic parameters, phylogenetic analysis, principal component analysis and Bayesian cluster analysis allowed us to refine genetic structure and provided insights in relationships of native populations, transboundary breeds and populations of domestic yak. Our results can find practical application in conservation programs of endangered breeds, as well as become the basis for future fundamental research.

Mitochondrial DNA diversity of D-loop region in three native Turkish cattle breeds

This study aimed to reveal the genetic variability of mitochondrial DNA (mtDNA) displacement-loop (D-loop) region in 62 animals belonging to three native Turkish cattle breeds, namely Anatolian Black (AB), East Anatolian Red (EAR) and Zavot (ZAV), and to conduct phylogenetic relationship analyses to obtain deeper information on their genetic origin and breeding history by comparison of 6 taurine and 11 indicine breeds, together with 66 polymorphic sites, a total of 31 haplotypes, of which 15, 10 and 6 were detected in AB, EAR and ZAV, respectively. Mean nucleotide and haplotype diversity were 0.01 and 0.891, respectively, whereas the genetic differentiation derived from Wright's F ST index was 0.174 across the breeds. A significant level of total variation (17.42 %) was observed among breeds in molecular variance analysis. Six main haplogroups (T, T1, T2, T3, Q and I2) were detected in Anatolian cattle populations, where T3 was the most frequent among breeds (43.55 %), whereas I2, an indicine specific haplogroup, was observed only in ZAV. At the breed level, phylogenetic analyses supported by 198 sequences of 17 cattle breeds and 3 outgroup species retrieved from the GenBank clustered native Turkish cattle breeds with the taurine group rather than the indicine one, as expected. However, indicine admixture at low frequency (8.89 %) was detected in the ZAV breed for the first time due to more likely gene flow from indicine cattle breeds raised in neighbour countries, particularly Iran. This finding should be further investigated in all native Turkish and indicine cattle breeds from nearby countries to clarify gene flow and indicine admixture in Anatolian cattle.

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.

Mitochondrial DNA diversity of the Sardinian local cattle stock

The aim of this research was to characterize the genetic diversity of the Sarda (Sa, n = 131), Sardo Bruna (SB, n = 44) and Sardo Modicana (SM, n = 26) cattle breeds, reared in the island of Sardinia (Italy). A portion of the mitochondrial DNA hypervariable region was sequenced, in order to identify a potential signature of African introgression. The F_ST coefficients among populations ranged between 0.056 for Sa vs SB and 0.167 for SB vs SM. AMOVA analysis indicated there was a significant differentiation of the three breeds, although most of diversity was gathered at the within-breed level. The Median Joining Network of the Sardinian sequences showed a potential founder effect signature. A MJ network including Sardinian cattle plus African, Italian, Iberian and Asian sequences, revealed the presence of haplogroup T3, already detected in Sa cattle, and the presence of Hg T1 and Hg T1′2′3, in Sa and SB. The presence of a private haplotype belonging to haplogroup T1, which is characteristic of African taurine breeds, may be due to the introgression of Sardinian breeds with African cattle, either directly (most probable source: North African cattle) or indirectly (through a Mediterranean intermediary already introgressed with African blood).

Genetic variation of Nigerian cattle inferred from maternal and paternal genetic markers

The African cattle provide unique genetic resources shaped up by both diverse tropical environmental conditions and human activities, the assessment of their genetic diversity will shade light on the mechanism of their remarkable adaptive capacities. We therefore analyzed the genetic diversity of cattle samples from Nigeria using both maternal and paternal DNA markers. Nigerian cattle can be assigned to 80 haplotypes based on the mitochondrial DNA (mtDNA) D-loop sequences and haplotype diversity was 0.985 + 0.005. The network showed two major matrilineal clustering: the dominant cluster constituting the Nigerian cattle together with other African cattle while the other clustered Eurasian cattle. Paternal analysis indicates only zebu haplogroup in Nigerian cattle with high genetic diversity 1.000 ± 0.016 compared to other cattle. There was no signal of maternal genetic structure in Nigerian cattle population, which may suggest an extensive genetic intermixing within the country. The absence of Bos indicus maternal signal in Nigerian cattle is attributable to vulnerability bottleneck of mtDNA lineages and concordance with the view of male zebu genetic introgression in African cattle. Our study shades light on the current genetic diversity in Nigerian cattle and population history in West Africa.

Population genetic structure of two herds of Aberdeen Angus cattle breed in Colombia

Background: Two biotypes of Aberdeen Angus cattle breed, known as Old Type and New Type, that differ in their origin and beef production are formally recognized. In Colombia, this breed has been commercialized for approximately 80 years. Studies on the origin, kinship and levels of genetic diversity of this breed in Colombian herds are scarce, yet important for planning crossing and management strategies. Objective: To measure the genetic diversity and structure of two Colombian herds of Old Type and New Type biotypes of Aberdeen Angus from Huila and Cundinamarca provinces and assess mitochondrial introgression with other breeds. Methods: A set of ten microsatellites and sequences of the Mitochondrial Control Region were characterized. Estimators of genetic diversity and population differentiation along with tests of population assignment were applied. Results: Nuclear loci were highly polymorphic as shown by the Polymorphic Information Content (0.599) and the Probability of Identity (1.896 10-08). Both populations were highly diverse and clearly differentiated into two groups corresponding to the Old Type and New Type phenotypes. In contrast, mitochondrial data failed to distinguish these two groups and showed extensive admixture. Conclusions: This study optimized a set of ten highly polymorphic nuclear markers that may be used for parentage and population genetic studies of Aberdeen Angus. Genetic differentiation in these loci agreed with phenotypic differences of the Old and New Types. However, mitochondrial data indicated ancestry of multiple European breeds in the origin of Colombian Aberdeen Angus.

Mitochondrial DNA Variations in Colombian Creole Sheep Confirm an Iberian Origin and Shed Light on the Dynamics of Introduction Events of African Genotypes

The genetic origins and diversity of Creole sheep from five regions of Colombia were investigated based on mitochondrial DNA (mtDNA) variations across 89 sequences from five breeds: one wool Creole sheep (CL) and four hair Creole sheep, including Ethiopian (OPCE), Sudan (OPCS), Pelibuey (OPCP) and Wayúu (OPCW). A global comparison was done using 62 haplotypes from Iberian, African, Indian, Caribbean, Mexican, Caucasian and European sheep based on sequences retrieved from GenBank. This study aimed to identify the maternal origin of Colombian Creole sheep and their genetic relationships at a global level. The results showed 31 different haplotypes from Colombian Creole sheep, which can be assigned to maternal lineage B, the most common lineage found in European sheep breeds and the only one found in several Iberian breed (e.g., Churra, Spanish Merino) that most likely participated in the Creole formation. Additional analyses showed that wool and hair sheep retained a broad genetic identity despite being geographically separated. The global-level phylogenetic analysis revealed that Colombian Creole sheep belong to a distinct and defined genetic lineage that is likely the result of a founder effect with ecotypes of Iberian descent and the subsequent introduction of foreign breeds. This is consistent with historical reports on the presence of sheep in South America and, particularly, Colombia.

Molecular inference in the colonization of cattle in Ecuador

The aim of the present paper was to evaluate genetic diversity within and among Ecuadorian bovine breeds using microsatellite markers. The genetic identity and the exotic influences from taurine or zebuine cosmopolitan breeds on Macabea, Loja, Santa Elena, Manabí and Isla de Puná breeds were evaluated. Twenty-eight microsatellite markers were tested across 254 individuals belonging to the five Ecuadorian populations to investigate intra and inter population genetic diversity levels. Nineteen markers in common with a set of 1595 samples of 35 previously typed breeds were used to determine the potential origin of Ecuadorian bovine breeds and to identify and quantify their genetic relationships. The findings of F_IS > F_ST (0.0814 > 0.0499), even in the context of low but significant F_ST values, may be indicative of the intrapopulation origin of the variability of allelic frequencies across populations. Conclusively, the study of genetic distances and population structure suggested the origin of Ecuadorian populations may likely stem from Southern Spanish breeds, with no significant recent influence from cosmopolitan Taurine or Zebuine breeds. These findings may provide a solid basis for the demonstration of an Ecuadorian breeds identity in the framework of the genetic conservation of the American creole populations.

Y-chromosome genetic diversity of Bos indicus cattle in close proximity to the centre of domestication

Y-chromosome genetic diversity in and around its domestication origin and a better understanding of indicine-specific microsatellite alleles are imperative concerns but less -targeted. We analysed Y-chromosome markers in 301 bulls representing 19 native Indian cattle (Bos indicus) and identified new alleles and haplotypes. Compared to other indicine studies, the high Y-haplotype diversity found in Indian cattle supports the hypothesis of greater genetic variability across the centre of origin decreasing along migratory routes with increasing distance. Hence, a considerable paternal genetic diversity of Indian cattle appears to have been lost in transboundary commercial indicine breeds. The Khillar and Gir are the most diversified populations where the first tends to be the well-differentiated traditional breed carrying strikingly distinct Y-lineages with typical BM861-158 bp allele, characteristics of taurine cattle, while retaining standard indicine lineages for all other markers. Geographical distribution found to be an unreliable predictor of parental variation, and Y-lineages seemed closely related to Indian breed function/utility. The comprehensive Y-chromosome information will be useful to examine the demographic expansion/spread of Bos indicus lineages from close proximity to the domestication centre across different countries worldwide and such diversity should be preserved through effective management and conservation programs.

Organization and Management of Conservation Programs and Research in Domestic Animal Genetic Resources

Creating national committees for domestic animal genetic resources within genetic resource national commissions is recommended to organize in situ and ex situ conservation initiatives. In situ conservation is a high priority because it retains traditional zootechnical contexts and locations to ensure the long-term survival of breeds. In situ actions can be based on subsidies, technical support, structure creation, or trademark definition. Provisional or permanent relocation of breeds may prevent immediate extinction when catastrophes, epizootics, or social conflicts compromise in situ conservation. Ex situ in vivo (animal preservation in rescue or quarantine centers) and in vitro methods (germplasm, tissues/cells, DNA/genes storage) are also potential options. Alert systems must detect emergencies and summon the national committee to implement appropriate procedures. Ex situ coordinated centers must be prepared to permanently or provisionally receive extremely endangered collections. National germplasm banks must maintain sufficient samples of national breeds (duplicated) in their collections to restore extinct populations at levels that guarantee the survival of biodiversity. A conservation management survey, describing national and international governmental and non-governmental structures, was developed. Conservation research initiatives for international domestic animal genetic resources from consortia centralize the efforts of studies on molecular, genomic or geo-evolutionary breed characterization, breed distinction, and functional gene identification. Several consortia also consider ex situ conservation relying on socioeconomic or cultural aspects. The CONBIAND network (Conservation for the Biodiversity of Local Domestic Animals for Sustainable Rural Development) exemplifies conservation efficiency maximization in a low-funding setting, integrating several Latin American consortia with international cooperation where limited human, material, and economic resources are available.

Abundant Genetic Diversity of Yunling Cattle Based on Mitochondrial Genome

Simple Summary

Mitochondrial DNA (mtDNA) analysis is a critical tool in assessing the maternal origin, phylogeny and population structure of domestic animals. Yunling cattle are a composite breed that was created by local Yunnan cattle breeds (as a maternal line) 30 years ago. It can be said that Yunling cattle represent important reservoirs of genetic diversity of local Yunnan cattle. Yunling cattle are characterized as a tropical/subtropical breed with fast growth, heat tolerance and parasite resistance, and can survive in a harsh environment and on low-quality roughage. Assessing the genetic characteristics of Yunling cattle is of particular importance for reasonable breeding strategies for Yunling cattle and the design of a local Yunnan cattle conservation program.

Abstract

Yunling cattle are a composite beef cattle breed, combining Brahman (1/2), Murray Grey (1/4) and Local Yunnan cattle (1/4), that was developed in Yunnan, China in the 1980s. Understanding the genetic information of Yunling cattle is of great significance to the development of reasonable breeding strategies for this breed. In the present study, we assessed the current genetic status of Yunling cattle in Yunnan Province (China) by analyzing the variability of the whole mitochondrial genome of 129 individuals. Altogether, 129 sequences displayed 47 different haplotypes. The haplotype diversity and the average number of nucleotide differences were 0.964 and 128.074, respectively. Phylogenetic analyses classified Yunling cattle into seven haplogroups: T1, T2, T3, T4, T6, I1 and I2. Haplogroup I1 was found to be predominant (41.86%), followed by T3 (28.68%). Furthermore, we also identified a novel haplogroup, T6, and defined the sub-haplogroup I1a in Yunling cattle. According to the formation process of Yunling cattle (local Yunnan cattle as the maternal line), the high genetic diversities in the mitochondria of Yunling cattle could be due to the complex maternal origin of local Yunnan cattle. Further studies about local Yunnan breeds are necessary to determine the exact source of haplogroup T6 in Yunling cattle. Our results will be useful for the evaluation and effective management of Yunling cattle.

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.

Deciphering the patterns of genetic admixture and diversity in southern European cattle using genome-wide SNPs

The divergence between indicine cattle (Bos indicus) and taurine cattle (Bos taurus) is estimated to have occurred approximately 250,000 years ago, but a small number of European cattle breeds still display shared ancestry with indicine cattle. Additionally, following the divergence of African and European taurine, the gene flow between African taurine and southern European cattle has also been proposed. However, the extent to which non-European cattle ancestry is diffused across southern European cattle has not been investigated thoroughly. Also, in recent times, many local breeds have suffered severe reductions in effective population size. Therefore, in the present study, we investigated the pattern of genetic diversity in various European cattle based on single nucleotide polymorphisms (SNP) identified from whole-genome sequencing data. Additionally, we also employed unlinked and phased SNP-based approaches on high-density SNP array data to characterize non-European cattle ancestry in several southern European cattle breeds. Using heterozygosity-based parameters, we concluded that, on average, nucleotide diversity is greater in southern European cattle than western European (British and commercial) cattle. However, an abundance of long runs of homozygosity (ROH) and the pattern of Linkage disequilibrium decay suggested recent bottlenecks in Maltese and Romagnola. High nucleotide diversity outside ROH indicated a highly diverse founder population for southern European and African taurine. We also show that Iberian cattle display shared ancestry with African cattle. Furthermore, we show that Podolica is an ancient cross-bred between Indicine zebu and European taurine. Additionally, we also inferred similar ancestry profile of non-European cattle ancestry in different Balkan and Italian cattle breeds which might be an indication of the common origin of indicine ancestry in these breeds. Finally, we discuss several plausible demographic scenarios which might account for the presence of non-European cattle ancestry in these cattle breeds.

Consequences of breed formation on patterns of genomic diversity and differentiation: the case of highly diverse peripheral Iberian cattle

BACKGROUND

Iberian primitive breeds exhibit a remarkable phenotypic diversity over a very limited geographical space. While genomic data are accumulating for most commercial cattle, it is still lacking for these primitive breeds. Whole genome data is key to understand the consequences of historic breed formation and the putative role of earlier admixture events in the observed diversity patterns.

RESULTS

We sequenced 48 genomes belonging to eight Iberian native breeds and found that the individual breeds are genetically very distinct with FST values ranging from 4 to 16% and have levels of nucleotide diversity similar or larger than those of their European counterparts, namely Jersey and Holstein. All eight breeds display significant gene flow or admixture from African taurine cattle and include mtDNA and Y-chromosome haplotypes from multiple origins. Furthermore, we detected a very low differentiation of chromosome X relative to autosomes within all analyzed taurine breeds, potentially reflecting male-biased gene flow.

CONCLUSIONS

Our results show that an overall complex history of admixture resulted in unexpectedly high levels of genomic diversity for breeds with seemingly limited geographic ranges that are distantly located from the main domestication center for taurine cattle in the Near East. This is likely to result from a combination of trading traditions and breeding practices in Mediterranean countries. We also found that the levels of differentiation of autosomes vs sex chromosomes across all studied taurine and indicine breeds are likely to have been affected by widespread breeding practices associated with male-biased gene flow.

Development of amplified fragment length polymorphism (AFLP) markers for the identification of Cholistani cattle

The identification issue of livestock can be resolved by using molecular identification tools that are acceptable to preserve and maintain pure breeds worldwide. The application of a molecular identification methodology is more important for developing nations, e.g., Pakistan, where uncontrolled crossbreeding has become a common practice and the import of exotic animals and germplasm is ever increasing. This presents a risk to local breeds as also stated by the FAO. Therefore, the current study was designed to develop standard molecular markers for Cholistani cattle to ascertain their purity for breeding purpose. In this study 50 and 48 unrelated males were sampled for Cholistani and each crossbred cattle, respectively. Candidate molecular markers present in Cholistani but absent in crossbred cattle and vice versa were detected using the amplified fragment length polymorphism (AFLP) method. Eleven markers were developed and were converted to single nucleotide polymorphism (SNP) markers for genotyping. The allele frequencies in both breeds were determined for discrimination ability using polymerase-chain-reaction-restriction-fragment-polymorphism (PCR-AFLP). The probability of identifying the Cholistani breed was 0.905 and the probability of misjudgment was 0.073 using a panel of markers. The identified markers can ascertain the breed purity and are likely to extend the facility for breed purity testing before entering into a genetic improvement program in the country.

New insights into the genetic composition and phylogenetic relationship of wolves and dogs in the Iberian Peninsula

This study investigates the gene pool of Portuguese autochthonous dog breeds and their wild counterpart, the Iberian wolf subspecies (Canis lupus signatus), using standard molecular markers. A combination of paternal and maternal molecular markers was used to investigate the genetic composition, genetic differentiation and genetic relationship of native Portuguese dogs and the Iberian wolf. A total of 196 unrelated dogs, including breed and village dogs from Portugal, and other dogs from Spain and North Africa, and 56 Iberian wolves (wild and captive) were analyzed for nuclear markers, namely Y chromosome SNPs, Y chromosome STR loci, autosomal STR loci, and a mitochondrial fragment of the control region I. Our data reveal new variants for the molecular markers and confirm significant genetic differentiation between Iberian wolf and native domestic dogs from Portugal. Based on our sampling, no signs of recent introgression between the two subspecies were detected. Y chromosome data do not reveal genetic differentiation among the analyzed dog breeds, suggesting they share the same patrilineal origin. Moreover, the genetic distinctiveness of the Iberian wolf from other wolf populations is further confirmed with the description of new mtDNA variants for this endemism. Our research also discloses new molecular markers for wolf and dog subspecies assignment, which might become particularly relevant in the case of forensic or noninvasive genetic studies. The Iberian wolf represents a relic of the once widespread wolf population in Europe and our study reveals that it is a reservoir of unique genetic diversity of the grey wolf, Canis lupus. These results stress the need for conservation plans that will guarantee the sustainability of this threatened top predator in Iberia.

Identification of a new Y chromosome haplogroup in Spanish native cattle

The aim of this work was to perform a thorough analysis of the diversity of Y-haplotypes in Spanish cattle. A total of 207 Bos taurus males were sampled across 25 European breeds, with a special focus on rare, local Spanish populations. Animals were genotyped with five Y-specific microsatellites (INRA189, UMN0103, UMN0307, BM861 and BYM1), two indels (ZFY10 and USP9Y) and one SNP (UTY19). A new haplogroup, distinct from those described by Götherström et al. (2005), was identified and named Y1.2. Samples representing the three B. taurus Y-haplogroups were genotyped for four additional Y chromosome SNPs (rs121919254, rs121919281, rs121919323 and rs137049553). Among these SNPs, only rs121919281 was informative in B. taurus and helped to confirm the new Y1.2 haplogroup. Analysis of a larger dataset of standardized haplotypes for 1507 individuals from 57 populations from Spain, other European countries and Africa showed the new Y1.2 haplogroup to be found exclusively in Spanish breeds. This finding reinforces the importance of local Spanish cattle as reservoirs of genetic diversity as well as the importance of the Iberian Peninsula in the history of cattle.

Genetic diversity and patterns of population structure in Creole goats from the Americas

Biodiversity studies are more efficient when large numbers of breeds belonging to several countries are involved, as they allow for an in-depth analysis of the within- and between-breed components of genetic diversity. A set of 21 microsatellites was used to investigate the genetic composition of 24 Creole goat breeds (910 animals) from 10 countries to estimate levels of genetic variability, infer population structure and understand genetic relationships among populations across the American continent. Three commercial transboundary breeds were included in the analyses to investigate admixture with Creole goats. Overall, the genetic diversity of Creole populations (mean number of alleles = 5.82 ± 1.14, observed heterozygosity = 0.585 ± 0.074) was moderate and slightly lower than what was detected in other studies with breeds from other regions. The Bayesian clustering analysis without prior information on source populations identified 22 breed clusters. Three groups comprised more than one population, namely from Brazil (Azul and Graúna; Moxotó and Repartida) and Argentina (Long and shorthair Chilluda, Pampeana Colorada and Angora-type goat). Substructure was found in Criolla Paraguaya. When prior information on sample origin was considered, 92% of the individuals were assigned to the source population (threshold q ≥ 0.700). Creole breeds are well-differentiated entities (mean coefficient of genetic differentiation = 0.111 ± 0.048, with the exception of isolated island populations). Dilution from admixture with commercial transboundary breeds appears to be negligible. Significant levels of inbreeding were detected (inbreeding coefficient > 0 in most Creole goat populations, P < 0.05). Our results provide a broad perspective on the extant genetic diversity of Creole goats, however further studies are needed to understand whether the observed geographical patterns of population structure may reflect the mode of goat colonization in the Americas.

Genetic diversity of the Mexican Lidia bovine breed and its divergence from the Spanish population

Lidia bovine breed exists since the XIV century in the Iberian Peninsula. These animals were initially produced for meat but some, showing an aggressive behaviour which difficulted their management, were used to participate in popular traditional and social events. A specialization of the breed giving rise to the original Lidia population is documented in Spain since mid-XVIII century. Following the same tradition than in the Spanish population, Mexico used aggressive animals at the beginning of the XX century until two families of breeders started importing Lidia breed bovines from Spain with the aim of specializing their production. Each family (Llaguno and González) followed different breeding managements, and currently, most of the Lidia Mexican population derives from the Llaguno line. Although genetic structure and diversity of the Spanish population have been studied (using autosomal microsatellite markers, Y chromosome DNA markers and mitochondrial DNA sequences), the Mexican population is not analysed. The aim of the study was to assess both the genetic structure and diversity of the Mexican Lidia breed and its relationship with the original Spanish population using the same molecular tools. A total of 306 animals belonging to 20 breeders issued from both existing Mexican families were genotyped, and the genetic information was compared to the previously existing Spanish information. Slightly higher levels of genetic diversity in Mexican population were found when comparing to the Spanish population, and the variability among populations accounted for differences within them showing mean values of 0.18 and 0.12, respectively. Animals from the Mexican breeders, belonging to each of the two families, clustered together, and there was little evidence of admixture with the Spanish population. The analysis of Y chromosome diversity showed a high frequency of the H6 haplotype in the Mexican population, whereas this haplotype is rare in the Spanish, which is only found in the Miura (100%) and Casta Navarra (38%) lineages. Mitochondrial DNA revealed similar haplotypic pattern in both Spanish and Mexican populations, which is in accordance with most of the Mediterranean bovine breeds. In conclusion, as the Mexican Lidia population had initially a small number of founders and its current population has been reared isolated from their Spanish ancestors since a long time, these bottleneck effects and a combination of mixed cattle origin are the factors that might erase any trace of the Spanish origin of this population.

An independent event of Neolithic cattle domestication on the South-eastern Balkans: evidence from prehistoric aurochs and cattle populations

Neolithic/Chalcolithic livestock domestication is an important issue for understanding the mode of life and economics of ancient human communities. The Balkans appears to be a crucial point for clarifying the socio-economical interrelations between the Oldest Middle Eastern/Anatolian and newly formed cultures in Europe. Two main hypotheses regarding the early history of cattle domestication, from their ancestor - the aurochs, have been discussed: multipoint domestication centers or single point origin and subsequent worldwide dissemination. In this study, we provide molecular data about the Balkan aurochs for the first time as well as additional information for the Neolithic/Chalcolithic cattle populations in this geographic location. A total of seventeen samples from different ancient settlements were analyzed according to D-loop control region. The results did not show different genetic profile of wild and domestic populations. All haplotypes were found to belong to the basic macro-haplogroup T. The majority of specimens (n = 14) were defined to form a new Balkan-specific T6 haplogroup. Only two of the ancient samples analyzed were assigned to the T3 haplotype predominating in Europe. We attempt to throw new light on the earliest cattle domestication events in Europe, thus, the results presented are discussed in two directions: (a) The possibility of local independent domestication processes in Neolithic South-Eastern Europe; (b) The single point domestication in the Middle East and subsequent cattle dissemination in Europe. Our data does not exclude the possibility for independent domestication events followed by a second wave of parallel dissemination of cattle herds via the Mediterranean route.

New Q lineage found in bovine (Bos taurus) of Iberian Peninsula

The northern Iberian Peninsula is home to a variety of autochthonous cattle breeds, such as the Terreña and Pirenaica. With the objective of characterizing the matrilineal lineages of these breeds, a study of mitochondrial DNA was performed. The D-loop of 155 individuals was analyzed and most of the individuals were carriers of the T3 haplogroup, while haplogroups T and T1 were much less frequent. A Pirenaica individual belonging to the Q haplogroup was found. To verify the presence of the Q haplogroup individual, the entire mitochondrial DNA was sequenced and compared with two descendants. The individuals were assigned to the Q1 sub-haplogroup. These findings extend the geographic distribution of the Q haplogroup to the south west of the European continent. This new Q1 lineage has seven polymorphisms in the coding region, so this lineage is probably as old as the Q lineages described to date.

Prospects and challenges for the conservation of farm animal genomic resources, 2015-2025

Livestock conservation practice is changing rapidly in light of policy developments, climate change and diversifying market demands. The last decade has seen a step change in technology and analytical approaches available to define, manage and conserve Farm Animal Genomic Resources (FAnGR). However, these rapid changes pose challenges for FAnGR conservation in terms of technological continuity, analytical capacity and integrative methodologies needed to fully exploit new, multidimensional data. The final conference of the ESF Genomic Resources program aimed to address these interdisciplinary problems in an attempt to contribute to the agenda for research and policy development directions during the coming decade. By 2020, according to the Convention on Biodiversity's Aichi Target 13, signatories should ensure that “…the genetic diversity of …farmed and domesticated animals and of wild relatives …is maintained, and strategies have been developed and implemented for minimizing genetic erosion and safeguarding their genetic diversity.” However, the real extent of genetic erosion is very difficult to measure using current data. Therefore, this challenging target demands better coverage, understanding and utilization of genomic and environmental data, the development of optimized ways to integrate these data with social and other sciences and policy analysis to enable more flexible, evidence-based models to underpin FAnGR conservation. At the conference, we attempted to identify the most important problems for effective livestock genomic resource conservation during the next decade. Twenty priority questions were identified that could be broadly categorized into challenges related to methodology, analytical approaches, data management and conservation. It should be acknowledged here that while the focus of our meeting was predominantly around genetics, genomics and animal science, many of the practical challenges facing conservation of genomic resources are societal in origin and are predicated on the value (e.g., socio-economic and cultural) of these resources to farmers, rural communities and society as a whole. The overall conclusion is that despite the fact that the livestock sector has been relatively well-organized in the application of genetic methodologies to date, there is still a large gap between the current state-of-the-art in the use of tools to characterize genomic resources and its application to many non-commercial and local breeds, hampering the consistent utilization of genetic and genomic data as indicators of genetic erosion and diversity. The livestock genomic sector therefore needs to make a concerted effort in the coming decade to enable to the democratization of the powerful tools that are now at its disposal, and to ensure that they are applied in the context of breed conservation as well as development.

Microsatellite and Mitochondrial DNA Study of Native Eastern European Cattle Populations: The Case of the Romanian Grey

The Eastern European Grey cattle are regarded as the direct descendants of the aurochs (Bos taurus primigenius). Nowadays in Romania, less than 100 Grey animals are being reared and included in the national gene reserve. We examined the genetic diversity among Romanian Grey, Brown, Spotted and Black and White cattle breeds, with a particular focus on Romanian Grey through the use of (i) 11 bovine specific microsatellite markers on 83 animals and (ii) 638 bp length of mitochondrial DNA (mtDNA) D-loop region sequence data from a total of 81 animals. Both microsatellite and mtDNA analysis revealed a high level of genetic variation in the studied breeds. In Romanian Grey a total of 100 alleles were found, the mean number of observed alleles per locus was 9.091; the average observed heterozygosity was 0.940; the Wright's fixation index (FIS) was negative (-0.189) and indicates that there is no inbreeding and no selection pressure. MtDNA analysis revealed 52 haplotypes with 67 variable sites among the Romanian cattle breeds without any insertion or deletion. Haplotype diversity was 0.980 ± 0.007 and ranged from 0.883 ± 0.056 (Brown) to 0.990 ± 0.028 (Spotted and Black and White). The highest genetic variability of the mtDNA was recorded in the Grey breed, where 18 haplotypes were identified. The most frequent mtDNA D-loop region belonged to T3 haplogroup (80.247%), which was found across all studied breeds, while T2 haplotypes (16.049%) was only found in Grey, Spotted and Black and White genotypes. The T1 haplotypes (3.704%) were found in the Grey and Spotted. The current results contribute to the general knowledge on genetic diversity found in Eastern European cattle breeds and could prove a valuable tool for the conservation efforts of animal genetic resources (FAnGR).

Impact of Parental Bos taurus and Bos indicus Origins on Copy Number Variation in Traditional Chinese Cattle Breeds

Copy number variation (CNV) is an important component of genomic structural variation and plays a role not only in evolutionary diversification but also in domestication. Chinese cattle were derived from Bos taurus and Bos indicus, and several breeds presumably are of hybrid origin, but the evolution of CNV regions (CNVRs) has not yet been examined in this context. Here, we of CNVRs, mtDNA D-loop sequence variation, and Y-chromosomal single nucleotide polymorphisms to assess the impact of maternal and paternal B. taurus and B. indicus origins on the distribution of CNVRs in 24 Chinese domesticated bulls. We discovered 470 genome-wide CNVRs, only 72 of which were shared by all three Y-lineages (B. taurus: Y₁, Y₂; B. indicus: Y₃), whereas 265 were shared by inferred taurine or indicine paternal lineages, and 228 when considering their maternal taurine or indicine origins. Phylogenetic analysis uncovered eight taurine/indicine hybrids, and principal component analysis on CNVs corroborated genomic exchange during hybridization. The distribution patterns of CNVRs tended to be lineage-specific, and correlation analysis revealed significant positive or negative co-occurrences of CNVRs across lineages. Our study suggests that CNVs in Chinese cattle partly result from selective breeding during domestication, but also from hybridization and introgression.

Temporal fluctuation in North East Baltic Sea region cattle population revealed by mitochondrial and Y-chromosomal DNA analyses

BACKGROUND

Ancient DNA analysis offers a way to detect changes in populations over time. To date, most studies of ancient cattle have focused on their domestication in prehistory, while only a limited number of studies have analysed later periods. Conversely, the genetic structure of modern cattle populations is well known given the undertaking of several molecular and population genetic studies.

RESULTS

Bones and teeth from ancient cattle populations from the North-East Baltic Sea region dated to the Prehistoric (Late Bronze and Iron Age, 5 samples), Medieval (14), and Post-Medieval (26) periods were investigated by sequencing 667 base pairs (bp) from the mitochondrial DNA (mtDNA) and 155 bp of intron 19 in the Y-chromosomal UTY gene. Comparison of maternal (mtDNA haplotypes) genetic diversity in ancient cattle (45 samples) with modern cattle populations in Europe and Asia (2094 samples) revealed 30 ancient mtDNA haplotypes, 24 of which were shared with modern breeds, while 6 were unique to the ancient samples. Of seven Y-chromosomal sequences determined from ancient samples, six were Y2 and one Y1 haplotype. Combined data including Swedish samples from the same periods (64 samples) was compared with the occurrence of Y-chromosomal haplotypes in modern cattle (1614 samples).

CONCLUSIONS

The diversity of haplogroups was highest in the Prehistoric samples, where many haplotypes were unique. The Medieval and Post-Medieval samples also show a high diversity with new haplotypes. Some of these haplotypes have become frequent in modern breeds in the Nordic Countries and North-Western Russia while other haplotypes have remained in only a few local breeds or seem to have been lost. A temporal shift in Y-chromosomal haplotypes from Y2 to Y1 was detected that corresponds with the appearance of new mtDNA haplotypes in the Medieval and Post-Medieval period. This suggests a replacement of the Prehistoric mtDNA and Y chromosomal haplotypes by new types of cattle.

Low levels of taurine introgression in the current Brazilian Nelore and Gir indicine cattle populations

Background

Nelore and Gir are the two most important indicine cattle breeds for production of beef and milk in Brazil. Historical records state that these breeds were introduced in Brazil from the Indian subcontinent, crossed to local taurine cattle in order to quickly increase the population size, and then backcrossed to the original breeds to recover indicine adaptive and productive traits. Previous investigations based on sparse DNA markers detected taurine admixture in these breeds. High-density genome-wide analyses can provide high-resolution information on the genetic composition of current Nelore and Gir populations, estimate more precisely the levels and nature of taurine introgression, and shed light on their history and the strategies that were used to expand these breeds.

Results

We used the high-density Illumina BovineHD BeadChip with more than 777 K single nucleotide polymorphisms (SNPs) that were reduced to 697 115 after quality control filtering to investigate the structure of Nelore and Gir populations and seven other worldwide populations for comparison. Multidimensional scaling and model-based ancestry estimation clearly separated the indicine, European taurine and African taurine ancestries. The average level of taurine introgression in the autosomal genome of Nelore and Gir breeds was less than 1% but was 9% for the Brahman breed. Analyses based on the mitochondrial SNPs present in the Illumina BovineHD BeadChip did not clearly differentiate taurine and indicine haplotype groupings.

Conclusions

The low level of taurine ancestry observed for both Nelore and Gir breeds confirms the historical records of crossbreeding and supports a strong directional selection against taurine haplotypes via backcrossing. Random sampling in production herds across the country and subsequent genotyping would be useful for a more complete view of the admixture levels in the commercial Nelore and Gir populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0109-5) contains supplementary material, which is available to authorized users.

On the History of Cattle Genetic Resources

Cattle are our most important livestock species because of their production and role in human culture. Many breeds that differ in appearance, performance and environmental adaptation are kept on all inhabited continents, but the historic origin of the diverse phenotypes is not always clear. We give an account of the history of cattle by integrating archaeological record and pictorial or written sources, scarce until 300 years ago, with the recent contributions of DNA analysis. We describe the domestication of their wild ancestor, migrations to eventually all inhabited continents, the developments during prehistory, the antiquity and the Middle Ages, the relatively recent breed formation, the industrial cattle husbandry in the Old and New World and the current efforts to preserve the cattle genetic resources. Surveying the available information, we propose three main and overlapping phases during the development of the present genetic diversity: (i) domestication and subsequent wild introgression; (ii) natural adaptation to a diverse agricultural habitat; and (iii) breed development.

Microsatellite analysis supports clonal propagation and reduced divergence of Trypanosoma vivax from asymptomatic to fatally infected livestock in South America compared to West Africa

BACKGROUND

Mechanical transmission of the major livestock pathogen Trypanosoma vivax by other biting flies than tsetse allows its spread from Africa to the New World. Genetic studies are restricted to a small number of isolates and based on molecular markers that evolve too slowly to resolve the relationships between American and West African populations and, thus, unable us to uncover the recent history of T. vivax in the New World.

METHODS

T. vivax genetic diversity, population structure and the source of outbreaks was investigated through the microsatellite multiloci (7 loci) genotype (MLGs) analysis in South America (47isolates from Brazil, Venezuela and French Guiana) and West Africa (12 isolates from The Gambia, Burkina Faso, Ghana, Benin and Nigeria). Relationships among MLGs were explored using phylogenetic, principal component and STRUCTURE analyses.

RESULTS

Although closely phylogenetically related, for the first time, genetic differences were detected between T. vivax isolates from South America (11 genotypes/47 isolates) and West Africa (12 genotypes/12 isolates) with no MLGs in common. Diversity was far greater across West Africa than in South America, where genotypes from Brazil (MLG1-6), Venezuela (MLG7-10) and French Guiana (MLG11) shared similar but not identical allele composition. No MLG was exclusive to asymptomatic (endemic areas) or sick (outbreaks in non-endemic areas) animals, but only MLGs1, 2 and 3 were responsible for severe haematological and neurological disorders.

CONCLUSIONS

Our results revealed closely related genotypes of T. vivax in Brazil and Venezuela, regardless of endemicity and clinical conditions of the infected livestock. The MLGs analysis from T. vivax across SA and WA support clonal propagation, and is consistent with the hypothesis that the SA populations examined here derived from common ancestors recently introduced from West Africa. The molecular markers defined here are valuable to assess the genetic diversity, to track the source and dispersion of outbreaks, and to explore the epidemiological and pathological significance of T. vivax genotypes.

Analysis of conservation priorities of Iberoamerican cattle based on autosomal microsatellite markers

Background

Determining the value of livestock breeds is essential to define conservation priorities, manage genetic diversity and allocate funds. Within- and between-breed genetic diversity need to be assessed to preserve the highest intra-specific variability. Information on genetic diversity and risk status is still lacking for many Creole cattle breeds from the Americas, despite their distinct evolutionary trajectories and adaptation to extreme environmental conditions.

Methods

A comprehensive genetic analysis of 67 Iberoamerican cattle breeds was carried out with 19 FAO-recommended microsatellites to assess conservation priorities. Contributions to global diversity were investigated using alternative methods, with different weights given to the within- and between-breed components of genetic diversity. Information on Iberoamerican plus 15 worldwide cattle breeds was used to investigate the contribution of geographical breed groups to global genetic diversity.

Results

Overall, Creole cattle breeds showed a high level of genetic diversity with the highest level found in breeds admixed with zebu cattle, which were clearly differentiated from all other breeds. Within-breed kinships revealed seven highly inbred Creole breeds for which measures are needed to avoid further genetic erosion. However, if contribution to heterozygosity was the only criterion considered, some of these breeds had the lowest priority for conservation decisions. The Weitzman approach prioritized highly differentiated breeds, such as Guabalá, Romosinuano, Cr. Patagonico, Siboney and Caracú, while kinship-based methods prioritized mainly zebu-related breeds. With the combined approaches, breed ranking depended on the weights given to the within- and between-breed components of diversity. Overall, the Creole groups of breeds were generally assigned a higher priority for conservation than the European groups of breeds.

Conclusions

Conservation priorities differed significantly according to the weight given to within- and between-breed genetic diversity. Thus, when establishing conservation programs, it is necessary to also take into account other features. Creole cattle and local isolated breeds retain a high level of genetic diversity. The development of sustainable breeding and crossbreeding programs for Creole breeds, and the added value resulting from their products should be taken into consideration to ensure their long-term survival.

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 diversity analysis of the Uruguayan Creole cattle breed using microsatellites and mtDNA markers

The Uruguayan Creole cattle population (N = 600) is located in a native habitat in south-east Uruguay. We analyzed its genetic diversity and compared it to other populations of American Creole cattle. A random sample of 64 animals was genotyped for a set of 17 microsatellite loci, and the D-loop hyper-variable region of mtDNA was sequenced for 28 calves of the same generation. We identified an average of 5.59 alleles per locus, with expected heterozygosities between 0.466 and 0.850 and an expected mean heterozygosity of 0.664. The polymorphic information content ranged from 0.360 to 0.820, and the global FIS index was 0.037. The D-loop analysis revealed three haplotypes (UY1, UY2 and UY3), belonging to the European matriline group, with a haplotype diversity of 0.532. The history of the population, changes in the effective population size, bottlenecks, and genetic drift are possible causes of the genetic variability patterns that we detected.

New World cattle show ancestry from multiple independent domestication events

Previous archeological and genetic research has shown that modern cattle breeds are descended from multiple independent domestication events of the wild aurochs (Bos primigenius) ∼10,000 y ago. Two primary areas of domestication in the Middle East/Europe and the Indian subcontinent resulted in taurine and indicine lines of cattle, respectively. American descendants of cattle brought by European explorers to the New World beginning in 1493 generally have been considered to belong to the taurine lineage. Our analyses of 47,506 single nucleotide polymorphisms show that these New World cattle breeds, as well as many related breeds of cattle in southern Europe, actually exhibit ancestry from both the taurine and indicine lineages. In this study, we show that, although European cattle are largely descended from the taurine lineage, gene flow from African cattle (partially of indicine origin) contributed substantial genomic components to both southern European cattle breeds and their New World descendants. New World cattle breeds, such as Texas Longhorns, provide an opportunity to study global population structure and domestication in cattle. Following their introduction into the Americas in the late 1400s, semiferal herds of cattle underwent between 80 and 200 generations of predominantly natural selection, as opposed to the human-mediated artificial selection of Old World breeding programs. Our analyses of global cattle breed population history show that the hybrid ancestry of New World breeds contributed genetic variation that likely facilitated the adaptation of these breeds to a novel environment.

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.

Mitochondrial analysis sheds light on the origin of hair sheep

A total of 180 mtDNA sequences from hair Caribbean (93), West African (73) and Canarian-wooled (14) sheep were analysed to shed light on the origin of hair sheep. A comparison of 360 Iberian sheep sequences retrieved from GenBank was performed to assess a possible European origin of the Caribbean hair sheep. These 180 sequences gave 48 different haplotypes (16 in Caribbean sheep). All Caribbean and Canarian-wooled sequences and 91.8% of the West African samples belonged to haplogroup B. The sheep analysed showed wide haplotypic identity. Caribbean sheep shared roughly two-thirds of their samples with Canarian-wooled and West African samples, respectively. Principal component analysis showed that the Caribbean and the Canarian-wooled sheep clustered together. Additional analyses showed that hair and Iberian sheep had wide genetic identity. It was not possible to ascertain a single Canarian, African or European origin of the Caribbean hair sheep using mtDNA markers only. European, African and Caribbean hair sheep maternal genetic backgrounds likely result from related domestication events.

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.

A novel USP9Y polymorphism allowing a rapid and unambiguous classification of Bos taurus Y chromosomes into haplogroups

A new sequence-tagged site (STS) was identified within intron 26 of the bovine USP9Y gene, showing an 81-base pair insertion (g.76439_76440ins81 in sequence with GenBank accession FJ195366) able to distinguish Y2 and Y3 Bos Y haplogroups from Y1. Moreover, four Y3-specific sequence variants allow a distinction from haplogroup Y2. The typing of a Bison bison Y chromosome indicates that the ancestral allele for the USP9Y 81-bp insertion is the short Y1 version. The results from typing the new STS in 1230 cattle Y chromosomes are fully consistent with their classification through standard methods. Thanks to the newly identified STS, it is now possible to assign cattle Y chromosomes to the currently known haplogroups using a single marker.

Genetic diversity, structure, and breed relationships in Iberian cattle

In Iberia there are 51 officially recognized cattle breeds of which 15 are found in Portugal and 38 in Spain. We present here a comprehensive analysis of the genetic diversity and structure of Iberian cattle. Forty of these breeds were genotyped with 19 highly polymorphic microsatellite markers. Asturiana de los Valles displayed the greatest allelic diversity and Mallorquina the least. Unbiased heterozygosity values ranged from 0.596 to 0.787. The network based on Reynolds distances was star-shaped with few pairs of interrelated breeds and a clear cluster of 4 breeds (Alistana/Arouquesa/Marinhoa/Mirandesa). The analysis of the genetic structure of Iberian cattle indicated that the most probable number of population clusters included in the study would be 36. Distance results were supported by the STRUCTURE software indicating a relatively recent origin or possible crossbreeding or both between pairs or small groups of breeds. Five clusters included 2 different breeds (Betizu/Pirenaica, Morucha/Avileña, Parda de Montaña/Bruna de los Pirineos, Barrosã/Cachena, and Toro de Lidia/Brava de Lide), 3 breeds (Berrenda en Negro, Negra Andaluza, and Mertolenga) were divided in 2 independent clusters each, and 2 breeds were considered admixed (Asturiana de los Valles and Berrenda en Colorado). Individual assignation to breeds was not possible in the 2 admixed breeds and the pair Parda de Montaña/Bruna de los Pirineos. The relationship between Iberian cattle reflects their geographical origin rather than their morphotypes. Exceptions to this geographic clustering are most probably a consequence of crossbreeding with foreign breeds. The relative genetic isolation within their geographical origin, the consequent genetic drift, the adaptation to specific environment and production systems, and the influence of African and European cattle have contributed to the current genetic status of Iberian cattle, which are grouped according to their geographical origin. The greater degree of admixture observed in some breeds should be taken into account before using molecular markers for genetic assignment of individuals to breeds.

Drift across the Atlantic: genetic differentiation and population structure in Brazilian and Portuguese native goat breeds

Brazilian goat breeds are believed to derive mainly from animals brought by Portuguese settlers since the 16th century. We used microsatellite markers in a sample of 436 animals to study genetic variability and differentiation of the six Portuguese (PT) and six Brazilian (BR) goat breeds currently recognized in the two countries. These breeds were also compared with an outgroup represented by a sample of Alpine (ALP) goats. The effective number of alleles and allelic richness were slightly higher in PT than in BR breeds. The global F(ST) was nearly 0.11 when PT and BR breeds were considered, with a mean pairwise F(ST) of about 0.03 among PT breeds, 0.07 among BR breeds and 0.15 between PT and BR breeds. The dendrogram illustrating relationships between populations and the correspondence analysis indicate the existence of two very distinct clusters, corresponding to the countries of origin of the breeds studied, which are nearly equidistant from the Alpine outgroup. The analysis with structure confirmed the separation between PT and BR breeds but suggests that some BR breeds, especially Graúna and Canindé, may share a common ancestry with PT breeds. The divergence observed between PT and BR breeds may result from founder effects and genetic drift but could also reflect the introduction in Brazil of goats originating from other regions, e.g., West Africa.