Bovine polledness--an autosomal dominant trait with allelic heterogeneity

The persistent horns are an important trait of speciation for the family Bovidae with complex morphogenesis taking place briefly after birth. The polledness is highly favourable in modern cattle breeding systems but serious animal welfare issues urge for a solution in the production of hornless cattle other than dehorning. Although the dominant inhibition of horn morphogenesis was discovered more than 70 years ago, and the causative mutation was mapped almost 20 years ago, its molecular nature remained unknown. Here, we report allelic heterogeneity of the POLLED locus. First, we mapped the POLLED locus to a ～381-kb interval in a multi-breed case-control design. Targeted re-sequencing of an enlarged candidate interval (547 kb) in 16 sires with known POLLED genotype did not detect a common allele associated with polled status. In eight sires of Alpine and Scottish origin (four polled versus four horned), we identified a single candidate mutation, a complex 202 bp insertion-deletion event that showed perfect association to the polled phenotype in various European cattle breeds, except Holstein-Friesian. The analysis of the same candidate interval in eight Holsteins identified five candidate variants which segregate as a 260 kb haplotype also perfectly associated with the POLLED gene without recombination or interference with the 202 bp insertion-deletion. We further identified bulls which are progeny tested as homozygous polled but bearing both, 202 bp insertion-deletion and Friesian haplotype. The distribution of genotypes of the two putative POLLED alleles in large semi-random sample (1,261 animals) supports the hypothesis of two independent mutations.     

Novel Features of the Prenatal Horn Bud Development in Cattle (Bos taurus)

Whereas the genetic background of horn growth in cattle has been studied extensively, little is known about the morphological changes in the developing fetal horn bud. In this study we histologically analyzed the development of horn buds of bovine fetuses between ~70 and ~268 days of pregnancy and compared them with biopsies taken from the frontal skin of the same fetuses. In addition we compared the samples from the wild type (horned) fetuses with samples taken from the horn bud region of age-matched genetically hornless (polled) fetuses. In summary, the horn bud with multiple layers of vacuolated keratinocytes is histologically visible early in fetal life already at around day 70 of gestation and can be easily differentiated from the much thinner epidermis of the frontal skin. However, at the gestation day (gd) 212 the epidermis above the horn bud shows a similar morphology to the epidermis of the frontal skin and the outstanding layers of vacuolated keratinocytes have disappeared. Immature hair follicles are seen in the frontal skin at gd 115 whereas hair follicles below the horn bud are not present until gd 155. Interestingly, thick nerve bundles appear in the dermis below the horn bud at gd 115. These nerve fibers grow in size over time and are prominent shortly before birth. Prominent nerve bundles are not present in the frontal skin of wild type or in polled fetuses at any time, indicating that the horn bud is a very sensitive area. The samples from the horn bud region from polled fetuses are histologically equivalent to samples taken from the frontal skin in horned species. This is the first study that presents unique histological data on bovine prenatal horn bud differentiation at different developmental stages which creates knowledge for a better understanding of recent molecular findings.     

The 80-kb DNA duplication on BTA1 is the only remaining candidate mutation for the polled phenotype of Friesian origin

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.     

Associations of single nucleotide polymorphisms in candidate genes with the polled trait in Datong domestic yaks

The domestic yak (Bos grunniens) is an iconic symbol of animal husbandry at high altitudes. Yaks exhibit unique external characteristics including long hair and large horns. However, hornless yaks can be found in different breeds and different populations. The hornless trait is also known as polled, and the POLL locus has been fine‐mapped to chromosome 1 in cattle (Bos taurus), although the underlying genetic basis of the polled trait is still unclear in the yak. Thus, we performed an association study to identify the genetic polymorphisms responsible for the polled trait in the yak. Fifty polled Datong domestic yaks and 51 horned individuals were selected randomly from a huge herd and were used as the case and control groups respectively for the association analysis. Twelve genes located in the candidate region of the POLL locus in cattle were used as references to detect DNA polymorphisms related to yak polledness, which were analyzed by sequencing and a high‐resolution melting test. We applied Fisher's exact test and haplotype analysis to show that a 147‐kb segment that included three protein‐coding genes C1H21orf62, GCFC1 and SYNJ1 was the most likely location of the POLL mutation in domestic yaks.     

Genetic diversity and investigation of polledness in divergent goat populations using 52 088 SNPs

The recent availability of a genome‐wide SNP array for the goat genome dramatically increases the power to investigate aspects of genetic diversity and to conduct genome‐wide association studies in this important domestic species. We collected and analysed genotypes from 52 088 SNPs in Boer, Cashmere and Rangeland goats that had both polled and horned individuals. Principal components analysis revealed a clear genetic division between animals for each population, and model‐based clustering successfully detected evidence of admixture that matched aspects of their recorded history. For example, shared co‐ancestry was detected, suggesting Boer goats have been introgressed into the Rangeland population. Further, allele frequency data successfully tracked the altered genetic profile that has taken place after 40 years of breeding Australian Cashmere goats using the Rangeland animals as the founding population. Genome‐wide association mapping of the POLL locus revealed a strong signal on goat chromosome 1. The 769‐kb critical interval contained the polled intersex syndrome locus, confirming the genetic basis in non‐European animals is the same as identified previously in Saanen goats. Interestingly, analysis of the haplotypes carried by a small set of sex‐reversed animals, known to be associated with polledness, revealed some animals carried the wild‐type chromosome associated with the presence of horns. This suggests a more complex basis for the relationship between polledness and the intersex condition than initially thought while validating the application of the goat SNP50 BeadChip for fine‐mapping traits in goat.     

Polledness in Argentinean Creole cattle,five centuries surviving

Polledness has been shown to have autosomal Mendelian inheritance, with the polled locus being dominant to the horned locus. This trait was mapped to the BTA1 centromeric end in several breeds. One of the distinctive attributes of Creole cattle, such as the Argentinean Creole, is the presence of long, lyre‐shaped horns. However, polled native animals were reported before the introduction of modern selected European breeds. Here, we studied the origin of the polled mutation, either independent or introgressed, in a Creole line from the Creole cattle founder group at the IIACS‐INTA Leales Experimental Station (northwest Argentina). The study sample (65 animals: 26 horned and 39 polled) was genotyped using high‐density SNP microarrays and three previously reported genetic markers (P_202ID, P_80kbID and P_G). A genome‐wide association study, selection signatures, linkage disequilibrium analysis and copy number variations were used to detect the responsible region and the segregating haplotypes/alleles. The interval mapped in the Leales herd (1.23–2.13 Mb) overlapped with the region previously reported in several European cattle breeds, suggesting that the same locus could be segregating in this population. The previously reported variants P_F and P_G were not detected, thus dismissing the Holstein‐Friesian and Nellore origins of the polled phenotype in this native breed. Conversely, the presence of the Celtic variant P_C suggests an almost complete co‐segregation. The cluster analysis rejected the hypothesis of recent introgression, which is compatible with the historical record of polled Creole cattle in northwest Argentina.     

A Single Nucleotide Polymorphism within the Interferon Gamma Receptor 2 Gene Perfectly Coincides with Polledness in Holstein Cattle

Polledness is a high impact trait in modern milk and beef production to meet the demands of animal welfare and work safety. Previous studies have mapped the polled-locus to the proximal region of the bovine chromosome 1 (BTA1) and narrowed it down to approximately 1 Mb. Sequencing of the positional candidate genes within the 1 Mb polled region and whole genome sequencing of Holsteins revealed a single nucleotide polymorphism (SNP) AC000158: g.1390292G>A within intron 3 of the interferon gamma receptor 2 gene (IFNGR2) in perfect co-segregation with polledness in Holsteins. This complete association was validated in 443 animals of the same breed. This SNP allows reliable genotyping of horned, heterozygous and homozygous polled Holsteins, even in animals that could not be resolved using the previously published haplotype for Holstein.     

Comparative iTRAQ proteomics revealed proteins associated with horn development in yak

Background

The practice of dehorning yak raises animal safety concerns, which have been addressed by selective breeding to obtain genetically hornless yak. The POLLED locus in yak has been studied extensively; however, little is known regarding the proteins that regulate horn bud development.

Methods

A differential proteomic analysis was performed to compare the skin from the horn bud region of polled yak fetuses and the horn bud tissue of horned yak fetuses using isobaric tags for relative and absolute quantitation (iTRAQ) technology coupled with 2D LC-MS/MS.

Results

One hundred differentially abundant proteins (DAPs) were identified. Of these, 29 were up-regulated and 71 were down-regulated in skin from the horn bud region of polled fetuses when compared to the horn bud tissue of horned fetuses. Bioinformatics analyses showed that the up-regulated DAPs were mainly associated with metabolic activities, while the down-regulated DAPs were significantly enriched in cell adhesion and cell movement activities.

Conclusions

We concluded that some important proteins were associated with cell adhesion, cell motility, keratinocyte differentiation, cytoskeleton organization, osteoblast differentiation, and fatty acid metabolism during horn bud development. These results advance our understanding of the molecular mechanisms underlying horn development.

     

Genetic association of marbling score with intragenic nucleotide variants at selection signals of the bovine genome

Selection signals of Korean cattle might be attributed largely to artificial selection for meat quality. Rapidly increased intragenic markers of newly annotated genes in the bovine genome would help overcome limited findings of genetic markers associated with meat quality at the selection signals in a previous study. The present study examined genetic associations of marbling score (MS) with intragenic nucleotide variants at selection signals of Korean cattle. A total of 39 092 nucleotide variants of 407 Korean cattle were utilized in the association analysis. A total of 129 variants were selected within newly annotated genes in the bovine genome. Their genetic associations were analyzed using the mixed model with random polygenic effects based on identical-by-state genetic relationships among animals in order to control for spurious associations produced by population structure. Genetic associations of MS were found (P<3.88×10⁻⁴) with six intragenic nucleotide variants on bovine autosomes 3 (cache domain containing 1, CACHD1), 5 (like-glycosyltransferase, LARGE), 16 (cell division cycle 42 binding protein kinase alpha, CDC42BPA) and 21 (snurportin 1, SNUPN; protein tyrosine phosphatase, non-receptor type 9, PTPN9; chondroitin sulfate proteoglycan 4, CSPG4). In particular, the genetic associations with CDC42BPA and LARGE were confirmed using an independent data set of Korean cattle. The results implied that allele frequencies of functional variants and their proximity variants have been augmented by directional selection for greater MS and remain selection signals in the bovine genome. Further studies of fine mapping would be useful to incorporate favorable alleles in marker-assisted selection for MS of Korean cattle.     

Mapping the four‐horned locus and testing the polled locus in three Chinese sheep breeds

Four‐horned sheep are an ideal animal model for illuminating the genetic basis of horn development. The objective of this study was to locate the genetic region responsible for the four‐horned phenotype and to verify a previously reported polled locus in three Chinese breeds. A genome‐wide association study (GWAS) was performed using 34 two‐horned and 32 four‐horned sheep from three Chinese indigenous breeds: Altay, Mongolian and Sishui Fur sheep. The top two significant single nucleotide polymorphisms (SNPs) associated with the four‐horned phenotype were both located in a region spanning positions 132.6 to 132.7 Mb on sheep chromosome 2. Similar locations for the four‐horned trait were previously identified in Jacob, Navajo‐Churro, Damara and Sishui Fur sheep, suggesting a common genetic component underlying the four‐horned phenotype. The two identified SNPs were both downstream of the metaxin 2 (MTX2) gene and the HOXD gene cluster. For the top SNP—OAR2:g.132619300G>A—the strong associations of the AA and AG genotypes with the four‐horned phenotype and the GG genotype with the two‐horned phenotype indicated the dominant inheritance of the four‐horned trait. No significant SNPs for the polled phenotype were identified in the GWAS analysis, and a PCR analysis for the detection of the 1.8‐kb insertion associated with polled sheep in other breeds failed to verify the association with polledness in the three Chinese breeds. This study supports the hypothesis that two different loci are responsible for horn existence and number. This study contributes to the understanding of the molecular regulation of horn development and enriches the knowledge of qualitative traits in domestic animals.     

牛、绵羊角的遗传定位及遗传机制研究进展

角属于动物颅骨附属物,为反刍动物所特有。牛(Bos taurus)、绵羊(Ovis aries)角的表型包括野生型两角表型、人工驯化的无角表型、畸形角等多种。牛和绵羊是阐明角的质量性状和数量性状之间的关系以及质量性状的多基因调控机制等方面的理想动物模型。近年来,对角性状研究不断深入,在阐明新器官起源进化、自然选择、性别选择和人工选择对角表型的影响等方面取得了一系列进展。本文详细介绍了角的研究概况、多角表型遗传定位、无角位点基因遗传定位和畸形角等,并对目前牛和绵羊角的遗传机制及存在的问题进行了分析,以期为反刍动物角性状和其他特异性性状遗传机制研究提供参考。     

A single nucleotide polymorphism on chromosome 10 is highly predictive for the polled phenotype in Australian Merino sheep

The aim of this study was to fine map the genomic location of the Horns locus in the Australian Merino sheep population and to identify markers that can be used to predict the horn phenotype. A linkage disequilibrium analysis of horn data from Australian Merino sheep mapped the Horns locus to a small region on chromosome 10. A single nucleotide polymorphism in the region was found to be highly predictive for the polled phenotype in an experimental population of Merino sheep. This was owing to a dominance effect of one of the alleles when inherited maternally. It was suggested that a genetic test would provide a good predictor of the polled phenotype. Finally, an evaluation of industry data showed that the SNP is at very different frequencies in Poll Merino sheep that have been bred for polledness (based on phenotype alone) compared with the Merino sheep breed.     

SNP-based association mapping of the polled gene in divergent cattle breeds

Naturally, hornless cattle are called polled. Although the POLL locus could be assigned to a c. 1.36‐Mb interval in the centromeric region of BTA1, the underlying genetic basis for the polled trait is still unknown. Here, an association mapping design was set up to refine the candidate region of the polled trait for subsequent high‐throughput sequencing. The case group comprised 101 homozygous polled animals from nine divergent cattle breeds, the majority represented by Galloway, Angus, Fleckvieh and Holstein Friesian. Additionally, this group included some polled individuals of Blonde d’Aquitaine, Charolais, Hereford, Jersey and Limousin breeds. The control group comprised horned Belgian Blue, Fleckvieh, Holstein Friesian and Illyrian Bu?a cattle. A genome‐wide scan using 49 163 SNPs was performed, which revealed one shared homozygous haplotype block consisting of nine neighbouring SNPs in all polled animals. This segment defines a 381‐kb interval on BTA1 that we consider to be the most likely location of the POLL mutation. Our results further demonstrate that the polled‐associated haplotype is also frequent in horned animals included in this study, and thus the haplotype as such cannot be used for population‐wide genetic testing. The actual trait‐associated haplotype may be revealed by using higher‐density SNP arrays. For the final identification of the causal mutation, we suggest high‐throughput sequencing of the entire candidate region, because the identification of functional candidate genes is difficult owing to the lack of a comparable model.     

Editor's choice: Functional classification of 15 million SNPs detected from diverse chicken populations

Next-generation sequencing has prompted a surge of discovery of millions of genetic variants from vertebrate genomes. Besides applications in genetic association and linkage studies, a fraction of these variants will have functional consequences. This study describes detection and characterization of 15 million SNPs from chicken genome with the goal to predict variants with potential functional implications (pfVars) from both coding and non-coding regions. The study reports: 183K amino acid-altering SNPs of which 48% predicted as evolutionary intolerant, 13K splicing variants, 51K likely to alter RNA secondary structures, 500K within most conserved elements and 3K from non-coding RNAs. Regions of local fixation within commercial broiler and layer lines were investigated as potential selective sweeps using genome-wide SNP data. Relationships with phenotypes, if any, of the pfVars were explored by overlaying the sweep regions with known QTLs. Based on this, the candidate genes and/or causal mutations for a number of important traits are discussed. Although the fixed variants within sweep regions were enriched with non-coding SNPs, some non-synonymous-intolerant mutations reached fixation, suggesting their possible adaptive advantage. The results presented in this study are expected to have important implications for future genomic research to identify candidate causal mutations and in poultry breeding.     

Identification of small and large genomic candidate variants in bovine pulmonary hypoplasia and anasarca syndrome

The pulmonary hypoplasia and anasarca syndrome (PHA) is a congenital lethal disorder, which until now has been reported in cattle and sheep. PHA is characterized by extensive subcutaneous fetal edema combined with hypoplasia or aplasia of the lungs and dysplasia of the lymphatic system. PHA is assumed to be of genetic etiology. This study presents the occurrence of PHA in two different cattle breeds and their genetic causation. Two PHA cases from one sire were observed in Slovenian Cika cattle. Under the assumption of monogenic inheritance, genome-wide homozygosity mapping scaled down the critical regions to 3% of the bovine genome including a 43.6 Mb-sized segment on chromosome 6. Whole-genome sequencing of one case, variant filtering against controls and genotyping of a larger cohort of Cika cattle led to the detection of a likely pathogenic protein-changing variant perfectly associated with the disease: a missense variant on chromosome 6 in ADAMTS3 (NM_001192797.1: c.1222C>T), which affects an evolutionary conserved residue (NP_001179726.1: p.(His408Tyr)). A single PHA case was found in Danish Holstein cattle and was whole-genome sequenced along with its parents. However, as there was no plausible private protein-changing variant, mining for structural variation revealed a likely pathogenic trisomy of the entire chromosome 20. The identified ADAMTS3 associated missense variant and the trisomy 20 are two different genetic causes, which shows a compelling genetic heterogeneity for bovine PHA.