Investigation of X- and Y-specific single nucleotide polymorphisms in taurine (Bos taurus) and indicine (Bos indicus) cattle

Initially, domesticated African cattle were of taurine type. Today, we find both African Bos taurus and Bos indicus cattle, as well as their crossbreeds, on the continent of Africa and they all share the same set of African taurine mitochondrial DNA haplogroups. In this study, we report genetic variation as substitutions and insertions/deletions (indels) on both the X and Y chromosomes, and use the variation to assess hybridization between taurine and indicine cattle. Six African cattle breeds (four Sanga breeds, including Raya Azebu, Danakil, Caprivi, Nguni; and two Zebu breeds, including Kilimanjaro Zebu and South Kavirondo Zebu) were screened for six new X-chromosomal markers, specifically three single nucleotide polymorphisms and three indels in the DDX3X (previously DBX ) and ZFX genes, and five previously identified Y-chromosomal markers in the DDX3Y (previously DBY ) and ZFY genes. In total, 90 (57 bulls and 33 cows) samples from the African breeds were analysed. We identify five diagnostic haplotypes of indicine and taurine origins on both the X and Y chromosomes. For each breed, the level of indicine introgression varies; in addition to pure taurine, indicine and hybrid X-chromosome individuals, recombinant X-chromosome variants were also detected. These markers are useful molecular tools for assessing the level of indicine admixture in African cattle breeds.     

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

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

Classification of the southern African sanga and East African shorthorned zebu

Humped African cattle, which are differentiated into zebu and sanga types, have traditionally been classified as Bos indicus . This paper discusses existing evidence and presents new evidence supporting the classification of southern African sangas as Bos taurus and East African zebus as ' taurindicus '. Classification is based on karyotype, frequencies of DNA markers and protein polymorphisms. The Boran, an East African zebu, has an acrocentric Y chromosome typical of Bos indicus . The southern African sanga breeds have a submetacentric Y chromosome typical of Bos taurus . Frequencies of four DNA markers support the hypothesis that the Tuli, a southern African sanga, had taurine ancestors and the Boran had both taurine and indicine ancestors. Frequencies for several protein polymorphisms strongly suggest that southern African sangas have more in common with taurine than with indicine breeds, while East African zebus are an admixture of African taurine and Asian indicine breeds.     

Y-SNPs Haplotype Diversity in Four Chinese Cattle Breeds

To investigate the genetic diversity of Chinese cattle, 96 male samples of 4 Chinese native cattle breeds were investigated using 5 single nucleotide polymorphisms specific to the bovine Y chromosome. Two previously described haplotypes (taurine Y2 and indicine Y3) were detected in 74 and 22 animals, respectively. The haplotype frequencies varied amongst the four native breeds. The taurine Y2 haplotype dominated in the Qinchuan, Dabieshan, and Yunba breeds. However, the indicine Y3 haplotype occurred in high frequency in the Enshi breed. Among the four native breeds, Yunba had the highest haplotype diversity (0.4330 ± 0.0750), followed by Qinchuan (0.2899 ± 0.1028) and Enshi (0.2222 ± 0.1662), Dabieshan was the least differentiated (0.1079 ± 0.0680). Compared with some foreign cattle breeds, the low level of haplotype diversity was detected in our breeds (0.2633 ± 0.1030).     

The extent of linkage disequilibrium in beef cattle breeds using high-density SNP genotypes

Background

The extent of linkage disequilibrium (LD) between molecular markers impacts genome-wide association studies and implementation of genomic selection. The availability of high-density single nucleotide polymorphism (SNP) genotyping platforms makes it possible to investigate LD at an unprecedented resolution. In this work, we characterised LD decay in breeds of beef cattle of taurine, indicine and composite origins and explored its variation across autosomes and the X chromosome.

Findings

In each breed, LD decayed rapidly and r² was less than 0.2 for marker pairs separated by 50 kb. The LD decay curves clustered into three groups of similar LD decay that distinguished the three main cattle types. At short distances between markers (< 10 kb), taurine breeds showed higher LD (r² = 0.45) than their indicine (r² = 0.25) and composite (r² = 0.32) counterparts. This higher LD in taurine breeds was attributed to a smaller effective population size and a stronger bottleneck during breed formation. Using all SNPs on only the X chromosome, the three cattle types could still be distinguished. However for taurine breeds, the LD decay on the X chromosome was much faster and the background level much lower than for indicine breeds and composite populations. When using only SNPs that were polymorphic in all breeds, the analysis of the X chromosome mimicked that of the autosomes.

Conclusions

The pattern of LD mirrored some aspects of the history of breed populations and showed a sharp decay with increasing physical distance between markers. We conclude that the availability of the HD chip can be used to detect association signals that remained hidden when using lower density genotyping platforms, since LD dropped below 0.2 at distances of 50 kb.     

Polymorphisms of two Y chromosome microsatellites in Chinese cattle

Two Y chromosome specific microsatellites UMN2404 and UMN0103 were genotyped and assessed for polymorphisms in a total of 423 unrelated males from 25 indigenous Chinese cattle breeds. Consistently, both microsatellites displayed specific indicine and taurine alleles in each bull examined. The indicine and taurine alleles were detected in 248 males (58.6%), and 175 males (41.4%), respectively, although these frequencies varied amongst different breeds examined. The indicine alleles dominated in the southern group (92.4%), while the taurine alleles dominated in the northern group (95.5%). Hainan Island was possibly the site for the origin of Chinese zebu, and Tibetan cattle were probably independently domesticated from another strain of Bos primigenius. The geographical distribution of these frequencies reveals a pattern of male indicine introgression and a hybrid zone of indicine and taurine cattle in China. The declining south-to-north and east-to-west gradient of male indicine introgression in China could be explained by historical data, geographical segregation and temperature and weather conditions.     

Y chromosome diversity and paternal origin of Chinese cattle

To determine the Y chromosome genetic diversity and paternal origin of Chinese cattle, 369 bulls from 17 Chinese native cattle breeds and 30 bulls from Holstein and four bulls from Burma were analyzed using a recently discovered USP9Y marker that could distinguish between taurine and indicine cattle more efficiently. In total, the taurine Y1, Y2 haplogroup and indicine Y3 haplogroup were detected in 7 (1.9 %), 193 (52.3 %) and 169 (45.8 %) individuals of 17 Chinese native breeds, respectively, although these frequencies varied amongst the Chinese native cattle breeds examined. Y2 dominates in northern China (91.4 %), while Y3 dominates in southern China (81.2 %). Central China is an admixture zone with Y2 predominating overall (72.0 %). Our results demonstrate that Chinese cattle have two paternal origins, one from B. taurus (Y2) and the other from B. indicus (Y3). The Y1 haplogroup may originate from the imported beef cattle breeds in western countries. The geographical distributions of the Y2 and Y3 haplogroup frequencies reveal a pattern of male indicine introgression from south to north China, and male taurine introgression from north to south China.     

Dynamics of genomic architecture during composite breed development in cattle

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

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.     

Paternal origins of Chinese cattle

To determine the genetic diversity and paternal origin of Chinese cattle, 302 males from 16 Chinese native cattle breeds as well as 30 Holstein males and four Burma males as controls were analysed using four Y‐SNPs and two Y‐STRs. In Chinese bulls, the taurine Y1 and Y2 haplogroups and indicine Y3 haplogroup were detected in seven, 172 and 123 individuals respectively, and these frequencies varied among the Chinese cattle breeds examined. Y2 dominates in northern China (91.4%), and Y3 dominates in southern China (90.8%). Central China is an admixture zone, although Y2 predominates overall (72.0%). The geographical distributions of the Y2 and Y3 haplogroup frequencies revealed a pattern of male indicine introgression from south to north China. The three Y haplogroups were further classified into one Y1 haplotype, five Y2 haplotypes and one Y3 haplotype in Chinese native bulls. Due to the interplay between taurine and indicine types, Chinese cattle represent an extensive reservoir of genetic diversity. The Y haplotype distribution of Chinese cattle exhibited a clear geographical structure, which is consistent with mtDNA, historical and geographical information.     

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

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

Male-mediated introgression of Bos indicus genes into Argentine and Bolivian Creole cattle breeds

The geographic distribution and frequency of Bos taurus and Bos indicus Y chromosome haplotypes amongst Argentine and Bolivian Creole cattle breeds were studied, using cytogenetic and molecular genetic techniques. A complete correspondence between Y chromosome morphology and the haplotype of the Y-linked microsatellite marker INRA 124 was found in all males examined. The taurine and indicine haplotypes were detected in 85.7 and 14.3% of the males studied, respectively, although these frequencies varied amongst the different breeds examined. The geographic distribution of this polymorphism suggests a pattern of zebu introgression in South America. The highest frequencies of the Zebu Y-chromosome are found in Brazilian populations (43-90%), in the eastern part of the continent, while it is absent in the southernmost breeds from Uruguay and Argentina. Bolivian breeds, at the centre of the continent, exhibit intermediate values (17-41%). This east/west and north/south gradient of male Zebu introgression could be explained by historical events and environmental factors.     

A genomic map of climate adaptation in Mediterranean cattle breeds

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

Inviability of hybrids between D. melanogaster and D. simulans results from the absence of simulans X not the presence of simulans Y chromosome.

Interspecific crosses between D. melanogaster and D. simulans or its sibling species result in unisexual inviability of the hybrids. Mostly, crosses of D. melanogaster females x D. simulans males produce hybrid females. On the other hand, only hybrid males are viable in the reciprocal crosses. A classical question is the cause of the unisexual hybrid inviability on the chromosomal level. Is it due to the absence of a D. simulans X chromosome or is it due to the presence of a D. simulans Y chromosome? A lack of adequate chromosomal rearrangements available in D. simulans has made it difficult to answer this question. However, it has been assumed that the lethality results from the absence of the D. simulans X rather than the presence of the D. simulans Y. Recently I synthesized the first D. simulans compound-XY chromosome that consists of almost the entire X and Y chromosomes. Males carrying the compound-XY and no free Y chromosome are fertile. By utilizing the compound-XY chromosome, the viability of hybrids with various constitutions of cytoplasm and sex chromosomes has been examined. The results consistently demonstrate that the absence of a D. simulans X chromosome in hybrid genome, and not the presence of the Y chromosome, is a determinant of the hybrid inviability.     

The Y chromosome as a target for acquired and amplified genetic material in evolution

The special properties of the Y chromosome stem form the fact that it is a non-recombining degenerate derivative of the X chromosome. The absence of homologous recombination between the X and the Y chromosome leads to gradual degeneration of various Y chromosome genes on an evolutionary timescale. The absence of recombination, however, also favors the accumulation of transposable elements on the Y chromosome during its evolution, as seen with both Drosophila and mammalian Y chromosomes. Alongside these processes, the acquisition and amplification of autosomal male benefit genes occur. This review will focus on recent studies that reveal the autosome-acquired genes on the Y chromosome of both Drosophila and humans. The evolution of the acquired and amplified genes on the Y chromosome is also discussed. Molecular and comparative analyses of Y-linked repeats in the Drosophila melanogaster genome demonstrate that there was a period of their degeneration followed by a period of their integration into RNAi silencing, which was beneficial for male fertility. Finally, the function of non-coding RNA produced by amplified Y chromosome genetic elements will be discussed.     

Interrogating cell division errors using random and chromosome-specific missegregation approaches

Accurate segregation of the duplicated genome in mitosis is essential for maintaining genetic stability. Errors in this process can cause numerical and/or structural chromosome abnormalities – hallmark genomic features commonly associated with both tumorigenesis and developmental disorders. A cell-based approach was recently developed permitting inducible missegregation of the human Y chromosome by selectively disrupting kinetochore assembly onto the Y centromere. Although this strategy initially requires several steps of genetic manipulation, it is easy to use, highly efficient and specific for the Y without affecting the autosomes or the X, and does not require cell cycle synchronization or mitotic perturbation. Here we describe currently available tools for studying chromosome segregation errors, aneuploidy, and micronuclei, as well as discuss how the Y-specific missegregation system has been used to elucidate how chromosomal micronucleation can trigger a class of extensive rearrangements termed chromothripsis. The combinatorial use of these different tools will allow unresolved aspects of cell division defects and chromosomal instability to be experimentally explored.     

The Genetic Architecture of Climatic Adaptation of Tropical Cattle

Adaptation of global food systems to climate change is essential to feed the world. Tropical cattle production, a mainstay of profitability for farmers in the developing world, is dominated by heat, lack of water, poor quality feedstuffs, parasites, and tropical diseases. In these systems European cattle suffer significant stock loss, and the cross breeding of taurine x indicine cattle is unpredictable due to the dilution of adaptation to heat and tropical diseases. We explored the genetic architecture of ten traits of tropical cattle production using genome wide association studies of 4,662 animals varying from 0% to 100% indicine. We show that nine of the ten have genetic architectures that include genes of major effect, and in one case, a single location that accounted for more than 71% of the genetic variation. One genetic region in particular had effects on parasite resistance, yearling weight, body condition score, coat colour and penile sheath score. This region, extending 20 Mb on BTA5, appeared to be under genetic selection possibly through maintenance of haplotypes by breeders. We found that the amount of genetic variation and the genetic correlations between traits did not depend upon the degree of indicine content in the animals. Climate change is expected to expand some conditions of the tropics to more temperate environments, which may impact negatively on global livestock health and production. Our results point to several important genes that have large effects on adaptation that could be introduced into more temperate cattle without detrimental effects on productivity.     

Fine‐mapping the POLL locus in Brahman cattle yields the diagnostic marker CSAFG29

The POLL locus has been mapped to the centromeric region of bovine chromosome 1 (BTA1) in both taurine breeds and taurine–indicine crosses in an interval of approximately 1 Mb. It has not yet been mapped in pure‐bred zebu cattle. Despite several efforts, neither causative mutations in candidate genes nor a singular diagnostic DNA marker has been identified. In this study, we genotyped a total of 68 Brahman cattle and 20 Hereford cattle informative for the POLL locus for 33 DNA microsatellites, 16 of which we identified de novo from the bovine genome sequence, mapping the POLL locus to the region of the genes IFNAR2 and SYNJ1. The 303‐bp allele of the new microsatellite, CSAFG29, showed strong association with the POLL allele. We then genotyped 855 Brahman cattle for CSAFG29 and confirmed the association between the 303‐bp allele and POLL. To determine whether the same association was found in taurine breeds, we genotyped 334 animals of the Angus, Hereford and Limousin breeds and 376 animals of the Brangus, Droughtmaster and Santa Gertrudis composite taurine–zebu breeds. The association between the 303‐bp allele and POLL was confirmed in these breeds; however, an additional allele (305 bp) was also associated but not fully predictive of POLL. Across the data, CSAFG29 was in sufficient linkage disequilibrium to the POLL allele in Australian Brahman cattle that it could potentially be used as a diagnostic marker in that breed, but this may not be the case in other breeds. Further, we provide confirmatory evidence that the scur phenotype generally occurs in animals that are heterozygous for the POLL allele.