Search

Background: The indigenous cattle populations from Greece and Cyprus have decreased to small numbers and are currently at risk of extinction due to socio-economic reasons, geographic isolation and crossbreeding with commercial breeds. This study represents the first comprehensive genome-wide analysis of 10 indigenous cattle populations from continental Greece and the Greek islands, and one from Cyprus, and compares them with 104 international breeds using more than 46,000 single nucleotide polymorphisms (SNPs)., Results: We estimated several parameters of genetic diversity (e.g. heterozygosity and allelic diversity) that indicated a severe loss of genetic diversity for the island populations compared to the mainland populations, which is mainly due to the declining size of their population in recent years and subsequent inbreeding. This high inbreeding status also resulted in higher genetic differentiation within the Greek and Cyprus cattle group compared to the remaining geographical breed groups. Supervised and unsupervised cluster analyses revealed that the phylogenetic patterns in the indigenous Greek breeds were consistent with their geographical origin and historical information regarding crosses with breeds of Anatolian or Balkan origin. Cyprus cattle showed a relatively high indicine ancestry. Greek island populations are placed close to the root of the tree as defined by Gir and the outgroup Yak, whereas the mainland breeds share a common historical origin with Buša. Unsupervised clustering and D-statistics analyses provided strong support for Bos indicus introgression in almost all the investigated local cattle breeds along the route from Anatolia up to the southern foothills of the Alps, as well as in most cattle breeds along the Apennine peninsula to the southern foothills of the Alps., Conclusions: All investigated Cyprus and Greek breeds present complex mosaic genomes as a result of historical and recent admixture events between neighbor and well-separated breeds. While the contribution of some mainland breeds to the genetic diversity pool seems important, some island and fragmented mainland breeds suffer from a severe decline of population size and loss of alleles due to genetic drift. Conservation programs that are a compromise between what is feasible and what is desirable should focus not only on the still highly diverse mainland breeds but also promote and explore the conservation possibilities for island breeds.

Background: Breeding genetically hornless, i.e. polled, cattle provides an animal welfare-friendly and non-invasive alternative to the dehorning of calves. However, the molecular regulation of the development of horns in cattle is still poorly understood. Studying genetic characters such as polledness and scurs, can provide valuable insights into this process. Scurs are hornlike formations that occur occasionally in a wide variety of sizes and forms as an unexpected phenotype when breeding polled cattle., Methods: We present a unique dataset of 885 Holstein-Friesian cattle with polled parentage. The horn phenotype was carefully examined, and the phenotypic heterogeneity of the trait is described. Using a direct gene test for polledness, the polled genotype of the animals was determined. Subsequently, the existence of a putative scurs locus was investigated using high-density genotype data of a selected subset of 232 animals and two mapping approaches: mixed linear model-based association analyses and combined linkage disequilibrium and linkage analysis., Results: The results of an exploratory data analysis indicated that the expression of scurs depends on age at phenotyping, sex and polled genotype. Scurs were more prevalent in males than in females. Moreover, homozygous polled animals did not express any pronounced scurs and we found that the Friesian polled allele suppresses the development of scurs more efficiently than the Celtic polled allele. Combined linkage and linkage disequilibrium mapping revealed four genome-wide significant loci that affect the development of scurs, one on BTA5 and three on BTA12. Moreover, suggestive associations were detected on BTA16, 18 and 23. The mixed linear model-based association analysis supports the results of the combined linkage and linkage disequilibrium analysis. None of the mapping approaches provided convincing evidence for a monogenic inheritance of scurs., Conclusions: Our results contradict the initial and still broadly accepted model for the inheritance of horns and scurs. We hypothesise an oligogenetic model to explain the development of scurs and polledness.

Background: The mode of inheritance of horn status in sheep is far more complex than a superficial analysis might suggest. Observations, which were mostly based on crossbreeding experiments, indicated that the allele that results in horns is dominant in males and recessive in females, and some authors even speculated about the involvement of more than two alleles. However, all recent genome-wide association analyses point towards a very strong effect of a single autosomal locus on ovine chromosome 10, which was narrowed down to a putatively causal insertion polymorphism in the 3'-untranslated region of the relaxin/insulin-like family peptide receptor 2 gene (RXFP2). The main objective of this study was to test this insertion polymorphism as the causal mutation in diverse sheep breeds, including breeds with a variable and/or sex-dependent horn status., Results: After re-sequencing a region of about 246 kb that covered the RFXP2 gene and its flanking regions for 24 sheep from six completely horned and six completely polled breeds, we identified the same insertion polymorphism that was previously published as segregating with horn status in these breeds. Multiplex PCR genotyping of 489 sheep from 34 breeds and some crosses between sheep breeds showed a nearly perfect segregation of the insertion polymorphism with horn status in sheep breeds of Central and Western European origin. In these breeds and their crossings, heterozygous males were horned and heterozygous females were polled. However, this segregation pattern was not, or at least not completely, reproducible in breeds with sex-dependent and/or variable horn status, especially in sheep that originated from even more southern European regions and from Africa. In such breeds, we observed almost all possible combinations of genotype, sex and horn status phenotype., Conclusions: The 1.78-kb insertion polymorphism in the 3'-untranslated region of RXFP2 and SNPs in the 3'-UTR, exon 14 and intron 11 of this gene that we analyzed in this study cannot be considered as the only cause of polledness in sheep and are not useful as a universal marker to define the genetic horn status in sheep.

Background: Bovine progressive degenerative myeloencephalopathy (Weaver syndrome) is a neurodegenerative disorder in Brown Swiss cattle that is characterized by progressive hind leg weakness and ataxia, while sensorium and spinal reflexes remain unaffected. Although the causal mutation has not been identified yet, an indirect genetic test based on six microsatellite markers and consequent exclusion of Weaver carriers from breeding have led to the complete absence of new cases for over two decades. Evaluation of disease status by imputation of 41 diagnostic single nucleotide polymorphisms (SNPs) and a common haplotype published in 2013 identified several suspected carriers in the current breeding population, which suggests a higher frequency of the Weaver allele than anticipated. In order to prevent the reemergence of the disease, this study aimed at mapping the gene that underlies Weaver syndrome and thus at providing the basis for direct genetic testing and monitoring of today's Braunvieh/Brown Swiss herds., Results: Combined linkage/linkage disequilibrium mapping on Bos taurus chromosome (BTA) 4 based on Illumina Bovine SNP50 genotypes of 43 Weaver-affected, 31 Weaver carrier and 86 Weaver-free animals resulted in a maximum likelihood ratio test statistic value at position 49,812,384 bp. The confidence interval (0.853 Mb) determined by the 2-LOD drop-off method was contained within a 1.72-Mb segment of extended homozygosity. Exploitation of whole-genome sequence data from two official Weaver carriers and 1145 other bulls that were sequenced in Run4 of the 1000 bull genomes project showed that only a non-synonymous SNP (rs800397662) within the PNPLA8 gene at position 49,878,773 bp was concordant with the Weaver carrier status. Targeted SNP genotyping confirmed this SNP as a candidate causal mutation for Weaver syndrome. Genotyping for the candidate causal mutation in a random sample of 2334 current Braunvieh animals suggested a frequency of the Weaver allele of 0.26 %., Conclusions: Through combined use of exhaustive sequencing data and SNP genotyping results, we were able to provide evidence that supports the non-synonymous mutation at position 49,878,773 bp as the most likely causal mutation for Weaver syndrome. Further studies are needed to uncover the exact mechanisms that underlie this syndrome.

Background: The absence of horns, called polled phenotype, is the favored trait in modern cattle husbandry. To date, polled cattle are obtained primarily by dehorning calves. Dehorning is a practice that raises animal welfare issues, which can be addressed by selecting for genetically hornless cattle. In the past 20 years, there have been many studies worldwide to identify unique genetic markers in complete association with the polled trait in cattle and recently, two different alleles at the POLLED locus, both resulting in the absence of horns, were reported: (1) the Celtic allele, which is responsible for the polled phenotype in most breeds and for which a single candidate mutation was detected and (2) the Friesian allele, which is responsible for the polled phenotype predominantly in the Holstein-Friesian breed and in a few other breeds, but for which five candidate mutations were identified in a 260-kb haplotype. Further studies based on genome-wide sequencing and high-density SNP (single nucleotide polymorphism) genotyping confirmed the existence of the Celtic and Friesian variants and narrowed down the causal Friesian haplotype to an interval of 145 kb., Results: Almost 6000 animals were genetically tested for the polled trait and we detected a recombinant animal which enabled us to reduce the Friesian POLLED haplotype to a single causal mutation, namely a 80-kb duplication. Moreover, our results clearly disagree with the recently reported perfect co-segregation of the POLLED mutation and a SNP at position 1 390 292 bp on bovine chromosome 1 in the Holstein-Friesian population., Conclusion: We conclude that the 80-kb duplication, as the only remaining variant within the shortened Friesian haplotype, represents the most likely causal mutation for the polled phenotype of Friesian origin.

Background: In a previous study in the Fleckvieh dual purpose cattle breed, we mapped a quantitative trait locus (QTL) affecting milk yield (MY1), milk protein yield (PY1) and milk fat yield (FY1) during first lactation to the distal part of bovine chromosome 5 (BTA5), but the confidence interval was too large for positional cloning of the causal gene. Our objective here was to refine the position of this QTL and to define the candidate region for high-throughput sequencing., Methods: In addition to those previously studied, new Fleckvieh families were genotyped, in order to increase the number of recombination events. Twelve new microsatellites and 240 SNP markers covering the most likely QTL region on BTA5 were analysed. Based on haplotype analysis performed in this complex pedigree, families segregating for the low frequency allele of this QTL (minor allele) were selected. Single- and multiple-QTL analyses using combined linkage and linkage disequilibrium methods were performed., Results: Single nucleotide polymorphism haplotype analyses on representative family sires and their ancestors revealed that the haplotype carrying the minor QTL allele is rare and most probably originates from a unique ancestor in the mapping population. Analyses of different subsets of families, created according to the results of haplotype analysis and availability of SNP and microsatellite data, refined the previously detected QTL affecting MY1 and PY1 to a region ranging from 117.962 Mb to 119.018 Mb (1.056 Mb) on BTA5. However, the possibility of a second QTL affecting only PY1 at 122.115 Mb was not ruled out., Conclusion: This study demonstrates that targeting families segregating for a less frequent QTL allele is a useful method. It improves the mapping resolution of the QTL, which is due to the division of the mapping population based on the results of the haplotype analysis and to the increased frequency of the minor allele in the families. Consequently, we succeeded in refining the region containing the previously detected QTL to 1 Mb on BTA5. This candidate region contains 27 genes with unknown or partially known function(s) and is small enough for high-throughput sequencing, which will allow future detailed analyses of candidate genes.

A joint analysis of five paternal half-sib Holstein families that were part of two different granddaughter designs (ADR- or Inra-design) was carried out for five milk production traits and somatic cell score in order to conduct a QTL confirmation study and to increase the experimental power. Data were exchanged in a coded and standardised form. The combined data set (JOINT-design) consisted of on average 231 sires per grandsire. Genetic maps were calculated for 133 markers distributed over nine chromosomes. QTL analyses were performed separately for each design and each trait. The results revealed QTL for milk production on chromosome 14, for milk yield on chromosome 5, and for fat content on chromosome 19 in both the ADR- and the Inra-design (confirmed within this study). Some QTL could only be mapped in either the ADR- or in the Inra-design (not confirmed within this study). Additional QTL previously undetected in the single designs were mapped in the JOINT-design for fat yield (chromosome 19 and 26), protein yield (chromosome 26), protein content (chromosome 5), and somatic cell score (chromosome 2 and 19) with genomewide significance. This study demonstrated the potential benefits of a combined analysis of data from different granddaughter designs.