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.     

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.     

The genetics of brown coat color and white spotting in domestic yaks (Bos grunniens)

Domestic yaks (Bos grunniens) exhibit two major coat color variations: a brown vs. wild‐type black pigmentation and a white spotting vs. wild‐type solid color pattern. The genetic basis for these variations in color and distribution remains largely unknown and may be complicated by a breeding history involving hybridization between yaks and cattle. Here, we investigated 92 domestic yaks from China using a candidate gene approach. Sequence variations in MC1R, PMEL and TYRP1 were surveyed in brown yaks; TYRP1 was unassociated with the coloration and excluded. Recessive mutations from MC1R, or p.Gln34*, p.Met73Leu and possibly p.Arg142Pro, are reported in bovids for the first time and accounted for approximately 40% of the brown yaks in this study. The remaining 60% of brown individuals correlated with a cattle‐derived deletion mutation from PMEL (p.Leu18del) in a dominant manner. Degrees of white spotting found in yaks vary from color sidedness and white face, to completely white. After examining the candidate gene KIT, we suggest that color‐sided and all‐white yaks are caused by the serial translations of KIT (Cs₆ or Cs₂₉) as reported for cattle. The white‐faced phenotype in yaks is associated with the KIT haplotype S^wf. All KIT mutations underlying the serial phenotypes of white spotting in yaks are identical to those in cattle, indicating that cattle are the likely source of white spotting in yaks. Our results reveal the complex genetic origins of domestic yak coat color as either native in yaks through evolution and domestication or as introduced from cattle through interspecific hybridization.     

Fine mapping of the polled locus to a 1-Mb region on bovine Chromosome 1q12

The absence of horns in Bos taurus is under genetic control of the autosomal dominant polled locus which has been genetically mapped to the centromeric region of cattle Chromosome 1. Recently a 4-Mb BAC contig of this chromosomal region has been constructed. Toward positional cloning of the bovine polled locus, we identified 20 additional microsatellite markers spread over the contig map by random sequencing of bacterial artificial chromosome (BAC) subclones. A total of 26 markers were genotyped in 30 two-generation half-sib families of six different German cattle breeds segregating for the hornless phenotype including 336 informative meioses for the polled character. Our fine-mapping study involving 19 recombinant haplotypes allowed us to narrow the critical region for the bovine polled locus to a 1-Mb segment with a centromeric boundary at RP42-218J17_MS1 and a telomeric boundary at BM6438. For marker-assisted selection purposes, the first evidence of informative flanking markers helps to predict polled genotypes with a higher degree of accuracy within families until testing of the causative mutation is available.     

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.     

Assessment of cattle genetic introgression into domestic yak populations using mitochondrial and microsatellite DNA markers

Hybridization between yak Poephagus grunniens and taurine Bos taurus or indicine B. indicus cattle has been widely practiced throughout the yak geographical range, and gene flow is expected to have occurred between these species. To assess the impact of cattle admixture on domestic yak, we examined 1076 domestic yak from 29 populations collected in China, Bhutan, Nepal, India, Pakistan, Kyrgyzstan, Mongolia and Russia using mitochondrial DNA and 17 autosomal microsatellite loci. A cattle diagnostic marker‐based analysis reveals cattle‐specific mtDNA and/or autosomal microsatellite allele introgression in 127 yak individuals from 22 populations. The mean level of cattle admixture across the populations, calculated using allelic information at 17 autosomal microsatellite loci, remains relatively low (mY_cattle = 2.66 ± 0.53% and Q_cattle = 0.69 ± 2.58%), although it varies a lot across populations as well as among individuals within population. Although the level of cattle admixture shows a clear geographical structure, with higher levels of admixture in the Qinghai‐Tibetan Plateau and Mongolian and Russian regions, and lower levels in the Himalayan and Pamir Plateau region, our results indicate that the level of cattle admixture is not significantly correlated with the altitude across geographical regions as well as within geographical region. Although yak‐cattle hybridization is primarily driven to produce F₁ hybrids, our results show that the subsequent gene flow between yak and cattle took place and has affected contemporary genetic make‐up of domestic yak. To protect yak genetic integrity, hybridization between yak and cattle should be tightly controlled.     

Comprehensive investigation of nucleotide diverdity in yaks

To understand the maternal genetic diversity of Tianzhu white yak better, we analyzed mtDNA D‐loop sequences of 209 Tianzhu white yaks, which are from the central region of Tianzhu white yak habitat. Accordingly, a total of 45 haplotypes were identified in Tianzhu white yaks in this study, and 18 of them were unique. The nucleotide diversity and haplotype diversity of population studied were 0.024 ± 0.003 and 0.946 ± 0.007 respectively, revealing that Tianzhu white yak possess a relatively high genetic diversity. The phylogenetic analysis, combining D‐loop sequences in this study with 533 previous published D‐loop sequences of 13 yak breeds, indicated that Tianzhu white yaks fell mainly into haplogroup A and that a small portion belonged to haplogroups B, C, D and E. Moreover, six haplotypes of 20 individuals identified in Tianzhu white yak were in the taurine haplogroup, indicating hybridization between Bos taurus and Tianzhu white yaks. In summary, this study supplies a comprehensive maternal genetic pattern for Tianzhu white yak and provides a basic reference for future breeding programs to conserve the purebred Tianzhu white yak.     

四个牦牛品种MC1R基因部分序列的多态性研究

为探索4个牦牛品种MC1R基因多态性的相关信息,选取甘南牦牛、天祝白牦牛、青海高原牦牛、大通牦牛4个品种共408头个体为研究对象,采用PCR-SSCP方法分析牦牛MC1R基因部分序列的基因多态性。结果表明,与GenBank中牛MCIR基因序列(登录号:AF445641.1)比对发现,该扩增片段在3 891 bp处发生C→G的突变,在3 912 bp处发生T→C的突变,共发现CC、DD、EE、CD、CE和DE 6种基因型。4个牦牛品种中CD、CE和DE 3种基因型在青海高原牦牛和大通牦牛中占主要优势,这3种基因型频率总和在青海高原牦牛和大通牦牛群体中分别是0.778和0.781。DD和CD两基因型是甘南牦牛群里中的优势基因型,其基因型频率分别是0.351和0.328。天祝白牦牛中优势基因型是DD,其基因型频率是0.500。D等位基因是4个地方品种牦牛中的优势等位基因。4个地方品种在该基因座上都处于Hardy-Weinberg平衡状态(P>0.05)。青海高原牦牛和大通牦牛两个群体处于高度多态(PIC>0.5),甘南牦牛和天祝白牦牛处于中度多态(0.25相似文献   

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.     

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.     

Structural Variants Selected during Yak Domestication Inferred from Long-Read Whole-Genome Sequencing

Structural variants (SVs) represent an important genetic resource for both natural and artificial selection. Here we present a chromosome-scale reference genome for domestic yak (Bos grunniens) that has longer contigs and scaffolds (N50 44.72 and 114.39 Mb, respectively) than reported for any other ruminant genome. We further obtained long-read resequencing data for 6 wild and 23 domestic yaks and constructed a genetic SV map of 372,220 SVs that covers the geographic range of the yaks. The majority of the SVs contains repetitive sequences and several are in or near genes. By comparing SVs in domestic and wild yaks, we identified genes that are predominantly related to the nervous system, behavior, immunity, and reproduction and may have been targeted by artificial selection during yak domestication. These findings provide new insights in the domestication of animals living at high altitude and highlight the importance of SVs in animal domestication.     

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.     

Screening for signatures of selection of Tianzhu white yak using genome‐wide re‐sequencing

The Tianzhu white yak, a domestic yak indigenous to the Qilian Mountains, migrated inland from the Qinghai‐Tibet Plateau. Specific ecological and long‐term artificial selection influenced the evolution of its pure white coat and physiological characteristics. Therefore, it is not only a natural population that represents a genomic selective region of environmental adaptability but is also an animal model for studying the pigmentation of the yak coat. A total of 24 261 829 variants, including 22 445 252 SNPs, were obtained from 29 yaks by genome‐wide re‐sequencing. According to the results of a selective sweep analysis of Tianzhu white yak in comparison to Tibetan yaks, nine candidate genes under selection in Tianzhu white yak were identified by combining π, Tajima's D, πA/πB and F_ST statistics, with threshold standards of 5%. These genes include PDCD1, NUP210, ABCG8, NEU4, LOC102287650, D2HGDH, COL4A1, RTP5 and HDAC11. Five of the nine genes were classified into 12 molecular signaling pathways, and most of these signaling pathways are involved in environmental information processing, organismal systems and metabolism. A majority of these genes has not been implicated in previous studies of yak coat color and high‐altitude animals. Our findings are helpful not only for explaining the molecular mechanism of yak coat pigmentation but also for exploring the genetic changes in Tianzhu white yak due to environmental adaptation.     

Phylogeographical analysis shows the need to protect the wild yaks' last refuge in Nepal

The wild yak Bos mutus was believed to be regionally extinct in Nepal for decades until our team documented two individuals from Upper Humla, north‐western Nepal, in 2014. The International Union for Conservation of Nature (IUCN) seeks further evidence for the conclusive confirmation of that sighting. We conducted line transects and opportunistic sign surveys in the potential wild yak habitats of Humla, Dolpa, and Mustang districts between 2015 and 2017 and collected genetic samples (present and historic) of wild and domestic yaks Bos grunniens. We also sighted another wild yak in Upper Humla in 2015. Phylogenetic and haplotype network analyses based on mitochondrial D‐loop sequences (~450 bp) revealed that wild yaks in Humla share the haplotype with wild yaks from the north‐western region of the Qinghai‐Tibetan Plateau in China. While hybridization with domestic yaks is a major long‐term threat, illegal hunting for meat and trophy put the very small populations of wild yaks in Nepal at risk. Our study indicates that the unprotected habitat of Upper Humla is the last refuge for wild yaks in Nepal. We recommend wild yak conservation efforts in the country to focus on Upper Humla by (i) assigning a formal status of protected area to the region, (ii) raising awareness in the local communities for wild yak conservation, and (iii) providing support for adaptation of herding practice and pastureland use to ensure the viability of the population.     

Novel Y‐chromosome polymorphisms in Chinese domestic yak

Y‐chromosome‐specific haplotypes (Y‐haplotypes) constructed using single nucleotide polymorphisms (Y‐SNPs) in the MSY (male‐specific region of the Y‐chromosome) are valuable in population genetic studies. But sequence variants in the yak MSY region have been poorly characterized so far. In this study, we screened a total of 16 Y‐chromosome‐specific gene segments from the ZFY, SRY, UTY, USP9Y, AMELY and OFD1Y genes to identify Y‐SNPs in domestic yaks. Six novel Y‐SNPs distributed in the USP9Y (g.223C>T), UTY19 (g.158A>C and g.169C>T), AMELY2 (g.261C>T), OFD1Y9 (g.165A>G) and SRY4 (g.104G>A) loci, which can define three Y‐haplotypes (YH1, YH2 and YH3) in yaks, were discovered. YH1 was the dominant and presumably most ancient haplotype based on the comparison of UTY19 locus with other bovid species. Interestingly, we found informative UTY19 markers (g.158A>C and g.169C>T) that can effectively distinguish the three yak Y‐haplotypes. The nucleotide diversity was 1.7 × 10^?4 ± 0.3 × 10^?4, indicating rich Y‐chromosome diversity in yaks. We identified two highly divergent lineages (YH1 and YH2 vs. YH3) that share similar frequencies (YH1 +  YH2: 0.82–0.89, YH3: 0.11–0.18) among all three populations. In agreement with previous mtDNA studies, we supported the hypothesis that the two highly divergent lineages (YH1 and YH2 vs. YH3) derived from a single gene pool, which can be explained by the reunion of at least two paternal populations with the divergent lineages already accumulated before domestication. We estimated a divergence time of 408 110 years between the two divergent lineages, which is consistent with the data from mitochondrial DNA in yaks.     

Accumulation of deleterious mutations in the domestic yak genome

Deleterious mutations play an important functional role, affecting trait phenotypes in ways that decrease the fitness of organisms. Estimating the frequency of occurrence and abundance has been a topic of much interest, especially in crops and livestock. The processes of domestication and breeding allow deleterious mutations to persist at high frequency, and identifying such deleterious mutations is particularly important for breed improvement. Here, we assessed genome‐wide patterns of deleterious variation in 59 domestic and 13 wild yaks using genome resequencing data. Based on the intersection of results given by three methods (provean , polyphen 2 and sift 4g ), we identified 3187 putative deleterious mutation sites affecting 2586 genes in domestic yaks and 2067 affecting 1701 genes in wild yaks. Multiple lines of evidence indicate a significant increase in the load of deleterious mutations in domesticated yaks compared to wild yaks. Private deleterious genes were found to be associated with the perception of smell and detection of chemical stimulus. We also identified 36 genes related to Mendelian genetic diseases involved in sensory perception, skeletal development and the nervous and immune systems. This study not only adds to the understanding of the genetic basis of yak domestication but also provides a rich catalog of variants that will facilitate future breeding‐related research on the yak genome and on other bovid species.     

Haematobia irritans parasitism of F1 yak × beef cattle (Bos grunniens × Bos taurus) hybrids

The horn fly Haematobia irritans (Diptera: Muscidae) is a blood obligate ectoparasite of bovids that causes annual losses to the U.S. beef cattle industry of over US$1.75 billion. Climate warming, the anthropogenic dispersion of bovids and the cross‐breeding of beef cattle with other bovid species may facilitate novel horn fly–host interactions. In particular, hybridizing yaks [Bos grunniens (Artiodactyla: Bovidae)] with beef cows (Bos taurus) for heterosis and carcass improvements may increase the exposure of yak × beef hybrids to horn flies. The present paper reports on the collection of digital images of commingled beef heifers (n = 12) and F1 yak × beef hybrid bovids (heifers, n = 7; steers, n = 5) near Laramie, Wyoming (～ 2200 m a.s.l.) in 2018. The total numbers of horn flies on beef heifers and F1 yak × beef heifers [mean ± standard error (SE): 88 ± 13 and 70 ± 17, respectively] did not differ significantly; however, F1 yak × beef steers had greater total horn fly abundance (mean ± SE: 159 ± 39) than female bovids. The present report of this experiment is the first such report in the literature and suggests that F1 yak × beef bovids are as susceptible as cattle to horn fly parasitism. Therefore, similar monitoring and treatment practices should be adopted by veterinarians, entomologists and producers.     

Genetic Structure and Differentiation of Three Chinese Indigenous Cattle Populations

Levels of genetic differentiation, gene flow, and genetic structure of three indigenous cattle populations (Luxi, Bohai, and Minnan) and two reference cattle populations (Chinese Holstein and Qinhai yak) in China were estimated using the information from 12 microsatellites, and 141 microsatellite alleles were identified. The mean number of alleles per locus ranged from 2.9005 in yak to 4.9722 in Holstein. The observed heterozygosity ranged from 0.5325 (yak) to 0.7719 (Holstein); 29 private alleles were detected. The global heterozygote deficit across all populations amounted to 58.5% (p < 0.001). The overall significant (p < 0.001) deficit of heterozygotes because of inbreeding within breeds amounted to 43.2%. The five cattle populations were highly differentiated (F _st = 26.9%, p < 0.001) at all loci. The heterozygote deficit within the population was highest in Luxi cattle and lowest in yak. The average number of effective migrants exchanged per generation was highest (1.149) between Luxi and Holstein, and lowest (0.509) between Luxi and yak. With the application of prior population information, cluster analysis achieved posterior probabilities from 91% to 98% of correctly assigning individuals to populations. Combining the information of cluster analysis, gene flow, and Structure analysis, the five cattle populations belong to three genetic clusters, a taurine (Luxi and Chinese Holstein), a zebu (Bohai and Minnan), and a yak cluster. This indicates that Bohai black is closer to Bos indicus than Luxi cattle. The evolution and development of three indigenous cattle populations are discussed.