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

Unexpected recombination at the polled locus in a horned Holstein calf from the mating of a homozygous polled sire and a heterozygous polled cow

In this study, I report an unexpected case of a Holstein calf that developed horns even though the sire was homozygous and the dam was heterozygous for polledness. After verifying and confirming the correct parentage, the parents and offspring were genotyped with the Illumina EuroG_MD BeadChip and the SNPs in the polled region on chromosome 1 were evaluated. In addition, the father was sequenced with next generation sequencing to identify possible, previously unknown variants. The deletion of two base pairs within the causative 80-kb duplication described for the Friesian polled variant was verified by melting curve analysis and the 80-kb duplication by droplet digital PCR. Analysis of all data showed that, as expected, the calf was heterozygous for all SNP positions flanking the 80-kb duplication but was homozygous wild type in the 80-kb duplication region and therefore carried horns. This is certainly a very rare case of a recombination within the highly conserved polled region, which on the one hand confirms that only the 80-kb duplication is responsible for the expression of the Friesian polled variant, but on the other hand also shows that caution is required when interpreting the usual routine genotyping of the horn status based on linked single nucleotide polymorphisms in the polled region on chromosome 1. Based on the present case, it is recommended that, in addition to the evaluation of the SNP data of the BeadChip, an extended diagnosis with direct detection of the 2-bp deletion (1:g.2629156_2629158delGT) should be carried out in any case when detecting the Friesian polled variant.

Genomic analysis of three medieval parchments from German monasteries

In the last two decades there has been growing interest in the analysis of ancient DNA obtained from the parchment used in historic documents. The genetic insight that this data provides makes collections of historic documents an invaluable source for studying the development and spread of historical livestock populations. Additionally, the biological data may provide new information for the historical analysis that could be used to determine the provenance as well as the authenticity of these documents. In this study, we extracted DNA from three medieval parchments that were written in German monasteries in the twelfth century. The source animal of the parchments could be identified as cattle and we compared their genome sequences with those of modern populations that are part of the 1000 Bull Genomes Project. The three animals were found to carry mtDNA haplogroup T3 and show a closer genetic relationship to other historic animals than to modern breeds. We further identified 39 haplotypes and 132 SNPs variants, which are rare ( < 0.1 ) or even non-existent in modern breeds. Finally, the genetic distances between the parchment samples show a putative association with the dates when the documents were written, indicating the usefulness of genetic analysis for provenance research.

A large structural variant collection in Holstein cattle and associated database for variant discovery, characterization, and application

BACKGROUND

Structural variants (SVs) such as deletions, duplications, and insertions are known to contribute to phenotypic variation but remain challenging to identify and genotype. A more complete, accessible, and assessable collection of SVs will assist efforts to study SV function in cattle and to incorporate SV genotyping into animal evaluation.

RESULTS

In this work we produced a large and deeply characterized collection of SVs in Holstein cattle using two popular SV callers (Manta and Smoove) and publicly available Illumina whole-genome sequence (WGS) read sets from 310 samples (290 male, 20 female, mean 20X coverage). Manta and Smoove identified 31 K and 68 K SVs, respectively. In total the SVs cover 5% (Manta) and 6% (Smoove) of the reference genome, in contrast to the 1% impacted by SNPs and indels. SV genotypes from each caller were confirmed to accurately recapitulate animal relationships estimated using WGS SNP genotypes from the same dataset, with Manta genotypes outperforming Smoove, and deletions outperforming duplications. To support efforts to link the SVs to phenotypic variation, overlapping and tag SNPs were identified for each SV, using genotype sets extracted from the WGS results corresponding to two bovine SNP chips (BovineSNP50 and BovineHD). 9% (Manta) and 11% (Smoove) of the SVs were found to have overlapping BovineHD panel SNPs, while 21% (Manta) and 9% (Smoove) have BovineHD panel tag SNPs. A custom interactive database ( https://svdb-dc.pslab.ca ) containing the identified sequence variants with extensive annotations, gene feature information, and BAM file content for all SVs was created to enable the evaluation and prioritization of SVs for further study. Illustrative examples involving the genes POPDC3, ORM1, G2E3, FANCI, TFB1M, FOXC2, N4BP2, GSTA3, and COPA show how this resource can be used to find well-supported genic SVs, determine SV breakpoints, design genotyping approaches, and identify processed pseudogenes masquerading as deletions.

CONCLUSIONS

The resources developed through this study can be used to explore sequence variation in Holstein cattle and to develop strategies for studying SVs of interest. The lack of overlapping and tag SNPs from commonly used SNP chips for most of the SVs suggests that other genotyping approaches will be needed (for example direct genotyping) to understand their potential contributions to phenotype. The included SV genotype assessments point to challenges in characterizing SVs, especially duplications, using short-read data and support ongoing efforts to better characterize cattle genomes through long-read sequencing. Lastly, the identification of previously known functional SVs and additional CDS-overlapping SVs supports the phenotypic relevance of this dataset.

DNA methylation analysis to differentiate reference, breed, and parent-of-origin effects in the bovine pangenome era

BACKGROUND

Most DNA methylation studies have used a single reference genome with little attention paid to the bias introduced due to the reference chosen. Reference genome artifacts and genetic variation, including single nucleotide polymorphisms (SNPs) and structural variants (SVs), can lead to differences in methylation sites (CpGs) between individuals of the same species. We analyzed whole-genome bisulfite sequencing data from the fetal liver of Angus (Bos taurus taurus), Brahman (Bos taurus indicus), and reciprocally crossed samples. Using reference genomes for each breed from the Bovine Pangenome Consortium, we investigated the influence of reference genome choice on the breed and parent-of-origin effects in methylome analyses.

RESULTS

Our findings revealed that ∼75% of CpG sites were shared between Angus and Brahman, ∼5% were breed specific, and ∼20% were unresolved. We demonstrated up to ∼2% quantification bias in global methylation when an incorrect reference genome was used. Furthermore, we found that SNPs impacted CpGs 13 times more than other autosomal sites (P < $5 \times {10}^{ - 324}$) and SVs contained 1.18 times (P < $5 \times {10}^{ - 324}$) more CpGs than non-SVs. We found a poor overlap between differentially methylated regions (DMRs) and differentially expressed genes (DEGs) and suggest that DMRs may be impacting enhancers that target these DEGs. DMRs overlapped with imprinted genes, of which 1, DGAT1, which is important for fat metabolism and weight gain, was found in the breed-specific and sire-of-origin comparisons.

CONCLUSIONS

This work demonstrates the need to consider reference genome effects to explore genetic and epigenetic differences accurately and identify DMRs involved in controlling certain genes.

In it for the long run: perspectives on exploiting long-read sequencing in livestock for population scale studies of structural variants

Studies have demonstrated that structural variants (SV) play a substantial role in the evolution of species and have an impact on Mendelian traits in the genome. However, unlike small variants (< 50 bp), it has been challenging to accurately identify and genotype SV at the population scale using short-read sequencing. Long-read sequencing technologies are becoming competitively priced and can address several of the disadvantages of short-read sequencing for the discovery and genotyping of SV. In livestock species, analysis of SV at the population scale still faces challenges due to the lack of resources, high costs, technological barriers, and computational limitations. In this review, we summarize recent progress in the characterization of SV in the major livestock species, the obstacles that still need to be overcome, as well as the future directions in this growing field. It seems timely that research communities pool resources to build global population-scale long-read sequencing consortiums for the major livestock species for which the application of genomic tools has become cost-effective.

LincRNA#1 knockout alone does not affect polled phenotype in cattle heterozygous for the celtic POLLED allele

A long intergenic non-coding RNA (lincRNA#1) is overexpressed in the horn bud region of polled (hornless) bovine fetuses, suggesting a potential role in horn bud suppression. Genome editing was used to test whether the absence of this sequence was associated with the horned phenotype. Two gRNAs with high mutation efficiencies targeting the 5' and the 3' regions flanking the lincRNA#1 sequence were co-injected with Cas9 as ribonucleoprotein complexes into bovine zygotes (n = 121) 6 h post insemination. Of the resulting blastocysts (n = 31), 84% had the expected 3.7 kb deletion; of these embryos with the 3.7 kb deletions, 88% were biallelic knockouts. Thirty-nine presumptive edited 7-day blastocysts were transferred to 13 synchronized recipient cows resulting in ten pregnancies, five with embryos heterozygous for the dominant P_C POLLED allele at the POLLED locus, and five with the recessive pp genotype. Eight (80%) of the resulting fetuses were biallelic lincRNA#1 knockouts, with the remaining two being mosaic. RT-qPCR analysis was used to confirm the absence of lincRNA#1 expression in knockout fetuses. Phenotypic and histological analysis of the genotypically (P_Cp) POLLED, lincRNA#1 knockout fetuses revealed similar morphology to non-edited, control polled fetuses, indicating the absence of lincRNA#1 alone does not result in a horned phenotype.

The Complex and Diverse Genetic Architecture of the Absence of Horns (Polledness) in Domestic Ruminants, including Goats and Sheep

Horns are the most obvious common feature of Bovidae. The naturally occurring absence of horns in these species, also known as polledness, is of surprisingly heterogeneous nature, although they are Mendelian traits. This review compares in detail the molecular differences among the causes of inherited polledness in the domestic ruminant species of cattle, yak, sheep, and goat based on the causal gene variants that have been discovered in recent years. The genetic causes for the lack of horns in small ruminants seem not only to be more complex, e.g., in sheep, breed-specific characteristics are still unexplained, but in goats, there is also the associated disorder of intersexuality—polled intersex syndrome (PIS). In connection with animal welfare and the associated discussion about a legal ban on the dehorning of all farm animals, naturally hornless animals and the causal genetic variants are of increasing research interest in the age of genome editing. However, the low acceptance of genetic engineering in livestock, especially in European societies, limits its use in food-producing animals. Therefore, genotype-based targeted selection of naturally occurring variants is still a widely used method for spreading this desired trait within and across populations, at least in cattle and sheep.

A deletion at the polled PC locus alone is not sufficient to cause a polled phenotype in cattle

Dehorning is a common practice in the dairy industry, but raises animal welfare concerns. A naturally occurring genetic mutation (PC allele) comprised of a 212 bp duplicated DNA sequence replacing a 10-bp sequence at the polled locus is associated with the hornless phenotype (polled) in cattle. To test the hypothesis that the 10 bp deletion alone is sufficient to result in polled, a CRISPR-Cas9 dual guide RNA approach was optimized to delete a 133 bp region including the 10 bp sequence. Timing of ribonucleoprotein complex injections at various hours post insemination (hpi) (6, 8, and 18 hpi) as well as in vitro transcribed (IVT) vs synthetic gRNAs were compared. Embryos injected 6 hpi had a significantly higher deletion rate (53%) compared to those injected 8 (12%) and 18 hpi (7%), and synthetic gRNAs had a significantly higher deletion rate (84%) compared to IVT gRNAs (53%). Embryo transfers were performed, and bovine fetuses were harvested between 3 and 5 months of gestation. All fetuses had mutations at the target site, with two of the seven having biallelic deletions, and yet they displayed horn bud development indicating that the 10 bp deletion alone is not sufficient to result in the polled phenotype.

A variance component estimation approach to infer associations between Mendelian polledness and quantitative production and female fertility traits in German Simmental cattle

BACKGROUND

Managing beneficial Mendelian characteristics in dairy cattle breeding programs implies that the correlated genetic effects are considered to avoid possible adverse effects in selection processes. The Mendelian trait polledness in cattle is traditionally associated with the belief that the polled locus has unfavorable effects on breeding goal traits. This may be due to the inferior breeding values of former polled bulls and cows in cattle breeds, such as German Simmental, or to pleiotropic or linkage effects of the polled locus.

METHODS

We focused on a variance component estimation approach that uses a marker-based numerator relationship matrix reflecting gametic relationships at the polled locus to test for direct pleiotropic or linked quantitative trait loci (QTL) effects of the polled locus on relevant traits. We applied the approach to performance, health, and female fertility traits in German Simmental cattle.

RESULTS

Our results showed no evidence for any pleiotropic QTL effects of the polled locus on test-day production traits milk yield and fat percentage, on the mastitis indicator 'somatic cell score', and on several female fertility traits, i.e. 56 days non return rate, days open and days to first service. We detected a significant and unfavorable QTL effect accounting for 6.6% of the genetic variance for protein percentage only.

CONCLUSIONS

Pleiotropy does not explain the lower breeding values and phenotypic inferiority of polled German Simmental sires and cows relative to the horned population in the breed. Thus, intensified selection in the polled population will contribute to increased selection response in breeding goal traits and genetic merit and will narrow the deficit in breeding values for production traits.

Investigating the Effect of Imputed Structural Variants from Whole-Genome Sequence on Genome-Wide Association and Genomic Prediction in Dairy Cattle

Simple Summary

Structural variants are large changes to the DNA sequences that differ from individual to individual. We discovered and quality-controlled a set of 24,908 structural variants and used a technique called imputation to infer them into 35,588 Holstein and Jersey cattle. We then investigated whether the structural variants affected key dairy cattle traits such as milk production, fertility and overall conformation. Structural variants explained generally less than 10 percent of the phenotypic variation in these traits. Four of the structural variants were significantly associated with dairy cattle production traits. However, the inclusion of the structural variants in the genomic prediction model did not increase genomic prediction accuracy.

Abstract

Structural variations (SVs) are large DNA segments of deletions, duplications, copy number variations, inversions and translocations in a re-sequenced genome compared to a reference genome. They have been found to be associated with several complex traits in dairy cattle and could potentially help to improve genomic prediction accuracy of dairy traits. Imputation of SVs was performed in individuals genotyped with single-nucleotide polymorphism (SNP) panels without the expense of sequencing them. In this study, we generated 24,908 high-quality SVs in a total of 478 whole-genome sequenced Holstein and Jersey cattle. We imputed 4489 SVs with R2 > 0.5 into 35,568 Holstein and Jersey dairy cattle with 578,999 SNPs with two pipelines, FImpute and Eagle2.3-Minimac3. Genome-wide association studies for production, fertility and overall type with these 4489 SVs revealed four significant SVs, of which two were highly linked to significant SNP. We also estimated the variance components for SNP and SV models for these traits using genomic best linear unbiased prediction (GBLUP). Furthermore, we assessed the effect on genomic prediction accuracy of adding SVs to GBLUP models. The estimated percentage of genetic variance captured by SVs for production traits was up to 4.57% for milk yield in bulls and 3.53% for protein yield in cows. Finally, no consistent increase in genomic prediction accuracy was observed when including SVs in GBLUP.

Comparison of Gene Editing Versus Conventional Breeding to Introgress the POLLED Allele Into the Tropically Adapted Australian Beef Cattle Population

Dehorning is the process of physically removing horns to protect animals and humans from injury, but the process is costly, unpleasant, and faces increasing public scrutiny. Genetic selection for polled (hornless), which is genetically dominant to horned, is a long-term solution to eliminate the need for dehorning. However, due to the limited number of polled Australian Brahman bulls, the northern Australian beef cattle population remains predominantly horned. The potential to use gene editing to produce high-genetic-merit polled cattle was recently demonstrated. To further explore the concept, this study simulated introgression of the POLLED allele into a tropically adapted Australian beef cattle population via conventional breeding or gene editing (top 1% or 10% of seedstock bulls/year) for 3 polled mating schemes and compared results to baseline selection on genetic merit (Japan Ox selection index, $JapOx) alone, over the course of 20 years. The baseline scenario did not significantly decrease the 20-year HORNED allele frequency (80%), but resulted in one of the fastest rates of genetic gain ($8.00/year). Compared to the baseline, the conventional breeding scenarios where polled bulls were preferentially used for breeding, regardless of their genetic merit, significantly decreased the 20-year HORNED allele frequency (30%), but resulted in a significantly slower rate of genetic gain ($6.70/year, P ≤ 0.05). The mating scheme that required the exclusive use of homozygous polled bulls, resulted in the lowest 20-year HORNED allele frequency (8%), but this conventional breeding scenario resulted in the slowest rate of genetic gain ($5.50/year). The addition of gene editing the top 1% or 10% of seedstock bull calves/year to each conventional breeding scenario resulted in significantly faster rates of genetic gain (up to $8.10/year, P ≤ 0.05). Overall, our study demonstrates that, due to the limited number of polled Australian Brahman bulls, strong selection pressure on polled will be necessary to meaningfully increase the number of polled animals in this population. Moreover, these scenarios illustrate how gene editing could be a tool for accelerating the development of high-genetic-merit homozygous polled sires to mitigate the current trade-off of slower genetic gain associated with decreasing HORNED allele frequency in the Australian Brahman population.

A de novo frameshift mutation in ZEB2 causes polledness, abnormal skull shape, small body stature and subfertility in Fleckvieh cattle

Polledness in cattle is an autosomal dominant trait. Previous studies have revealed allelic heterogeneity at the polled locus and four different variants were identified, all in intergenic regions. In this study, we report a case of polled bull (FV-Polled1) born to horned parents, indicating a de novo origin of this polled condition. Using 50K genotyping and whole genome sequencing data, we identified on chromosome 2 an 11-bp deletion (AC_000159.1:g.52364063_52364073del; Del11) in the second exon of ZEB2 gene as the causal mutation for this de novo polled condition. We predicted that the deletion would shorten the protein product of ZEB2 by almost 91%. Moreover, we showed that all animals carrying Del11 mutation displayed symptoms similar to Mowat-Wilson syndrome (MWS) in humans, which is also associated with genetic variations in ZEB2. The symptoms in cattle include delayed maturity, small body stature and abnormal shape of skull. This is the first report of a de novo dominant mutation affecting only ZEB2 and associated with a genetic absence of horns. Therefore our results demonstrate undoubtedly that ZEB2 plays an important role in the process of horn ontogenesis as well as in the regulation of overall development and growth of animals.

Heterozygosity of the Celtic polled locus in Canadian scurred beef cattle

Polled cattle are preferable to horned or scurred animals because they are safer for handling and cause less bruising. Although DNA testing can determine horned/polled genotype, scurs may appear in polled animals. The inheritance of scurs is complex because it is a sex-influenced trait that interacts with the polled locus. We demonstrate that in 685 purebred and crossbred Canadian beef cattle, all 153 scurred animals were heterozygous polled at the Celtic variant. In addition, male obligate carriers of scurs were smooth polled when homozygous for the polled mutation. Scurred and non-scurred males were sequenced for five genes (CTDNEP1, SHBG, SOX15, FGF11, and DHRS7C) within the scur candidate region on BTA19 that are functionally related to bone development and hormone regulation. Multipoint linkage analysis was conducted using 18 microsatellite markers and two informative variants (DHRS7C g.29594018G>C and CTDNEP1 c.462G>A) in the scurred families and further supported mapping on BTA19 between BMS2142 (logarithm of the odds (LOD) = 5.42) and IDVGA46 (LOD = 3.47). These data indicate epistatic interactions between the scurred and polled loci and emphasise the necessity for a scurred DNA test to assist purebred beef producers in eradicating the scur trait.

Are scurs in heterozygous polled (Pp) cattle a complex quantitative trait?

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.

Optimized Genetic Testing for Polledness in Multiple Breeds of Cattle

Many breeds of modern cattle are naturally horned, and for sound husbandry management reasons the calves frequently undergo procedures to physically remove the horns by disbudding or dehorning. These procedures are however a welfare concern. Selective breeding for polledness - absence of horns - has been effective in some cattle breeds but not in others (Bos indicus genotypes) due in part to the complex genetics of horn phenotype. To address this problem different approaches to genetic testing which provide accurate early-in-life prediction of horn phenotype have been evaluated, initially using microsatellites (MSAT) and more recently single nucleotide polymorphism (SNP). A direct gene test is not effective given the genetic heterogeneity and large-sized sequence variants associated with polledness in different breeds. The current study investigated 39,943 animals of multiple breeds to assess the accuracy of available poll testing assays. While the standard SNP-based test was an improvement on the earlier MSAT haplotyping method, 1999 (9.69%) out of 20,636 animals tested with this SNP-based assay did not predict a genotype, most commonly associated with the Indicus-influenced breeds. The current study has developed an optimized poll gene test that resolved the vast majority of these 1999 unresolved animals, while the predicted genotypes of those previously resolved remained unchanged. Hence the optimized poll test successfully predicted a genotype in 99.96% of samples assessed. We demonstrated that a robust set of 5 SNPs can effectively determine P_C and P_F alleles and eliminate the ambiguous and undetermined results of poll gene testing previously identified as an issue in cattle.

Understanding the effects of the bovine POLLED variants

Horns are paired appendages on the head of bovine species, comprising an inner bony core and outer keratin sheath. The horn bud forms during early fetal development but ossification of the developing horn does not occur until approximately 1 month after birth. Little is known about the genetic pathways that lead to horn growth. Hornless, or polled, animals are found in all domestic bovids. Histological studies of bovine fetuses have shown that the horn bud does not form in polled individuals. There are currently four known genetic variants for polledness in cattle on BTA1. All of the variants are intergenic, but probably affect regulation of nearby genes or long non-coding RNAs. Transcriptomic studies suggest that the expression of two nearby long non-coding RNAs are affected by the Celtic POLLED variant, but further studies are required to confirm these data. Candidate genes located elsewhere in the genome are involved in regulating bone formation and epithelial-to-mesenchymal transition. Expression of one of these candidate genes, RXFP2, appears to be reduced in the fetal horn bud of polled animals carrying the Celtic variant compared with horned individuals. Investigating horn ontogenesis and the genetic pathway by which the POLLED variants prevent horn development has implications for cattle breeding. If the genetic basis of horn bud formation and polledness is better understood, then new targets may be identified for precision genome editing to create polled individuals.

Genotypes and allele frequencies of buried SNPs in a bovine single-nucleotide polymorphism array in Japanese Black cattle

Single nucleotide polymorphism (SNP) arrays are widely used for genetic and genomic analyses in cattle breeding; thus, data derived from SNP arrays have accumulated on a large scale nationwide. Commercial SNP arrays contain a considerable number of unassigned SNPs on the chromosome/position on the genome; these SNPs are excluded in subsequent analyses. Notably, the position-unassigned SNPs, or "buried SNPs" include some of the markers associated with genetic disease. In this study, we identified the position of buried SNPs using the Basic Local Alignment Search Tool against the surrounding sequences and characterized the relationship between SNPs and genetic diseases in Online Mendelian Inheritance in Animals based on the genomic position. We determined the position of 285 buried SNPs on the genome and surveyed the genotype and allele frequencies of these SNPs in 5,955 individual Japanese Black cattle. Eleven SNPs associated with genetic disease, which contained five buried SNPs, were found in the population with the risk allele frequency ranging from 0.00008396 to 0.46. These results indicate that buried SNPs in the bovine SNP array can be utilized to identify associations with genetic disorders from large scale accumulated SNP genotype data in Japanese Black cattle.

Short communication: Phenotypic and genetic effects of the polled haplotype on yield, longevity, and fertility in US Brown Swiss, Holstein, and Jersey cattle

Phenotypes from the December 2018 US national genetic evaluations were used to compute effects of the polled haplotype in US Brown Swiss (BS), Holstein (HO), and Jersey (JE) cattle on milk, fat, and protein yields, somatic cell score, single-trait productive life, daughter pregnancy rate, heifer conception rate, and cow conception rate. Lactation records pre-adjusted for nongenetic factors and direct genomic values were used to estimate phenotypic and genetic effects of the polled haplotype, respectively. No phenotypic or direct genomic values effects were different from zero for any trait in any breed. Genomic PTA (gPTA) for the lifetime net merit (NM$) selection index of bulls born since January 1, 2012, that received a marketing code from the National Association of Animal Breeders (Madison, WI), and cows born on or after January 1, 2015, were compared to determine whether there was a systematic benefit to polled or horned genetics. Horned bulls had the highest average gPTA for NM$ in all 3 breeds, but that difference was significant only in HO and JE (HO: 615.4 ± 1.9, JE: 402.3 ± 3.4). Homozygous polled BS cows had significantly higher average gPTA for NM$ than their heterozygous polled or horned contemporaries (PP = 261.4 ± 43.5, Pp = 166.1 ± 13.7, pp = 174.1 ± 1.8), but the sample size was very small (n = 9). In HO and JE, horned cows had higher gPTA for NM$ (HO = 378.3 ± 0.2, JE = 283.3 ± 0.3). Selection for polled cattle should not have a detrimental effect on yield, fertility, or longevity, but these differences show that, in the short term, selection for polled over horned cattle will result in lower rates of genetic gain.

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_{202 ID} , P_{80kb ID} 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.

Proposed U.S. regulation of gene-edited food animals is not fit for purpose

Dietary DNA is generally regarded as safe to consume, and is a routine ingredient of food obtained from any living organism. Millions of naturally-occurring DNA variations are observed when comparing the genomic sequence of any two healthy individuals of a given species. Breeders routinely select desired traits resulting from this DNA variation to develop new cultivars and varieties of food plants and animals. Regulatory agencies do not evaluate these new varieties prior to commercial release. Gene editing tools now allow plant and animal breeders to precisely introduce useful genetic variation into agricultural breeding programs. The U.S. Department of Agriculture (USDA) announced that it has no plans to place additional regulations on gene-edited plants that could otherwise have been developed through traditional breeding prior to commercialization. However, the U.S. Food and Drug Administration (FDA) has proposed mandatory premarket new animal drug regulatory evaluation for all food animals whose genomes have been intentionally altered using modern molecular technologies including gene editing technologies. This runs counter to U.S. biotechnology policy that regulatory oversight should be triggered by unreasonable risk, and not by the fact that an organism has been modified by a particular process or technique. Breeder intention is not associated with product risk. Harmonizing the regulations associated with gene editing in food species is imperative to allow both plant and animal breeders access to gene editing tools to introduce useful sustainability traits like disease resistance, climate adaptability, and food quality attributes into U.S. agricultural breeding programs.

Comparison of gene editing versus conventional breeding to introgress the POLLED allele into the US dairy cattle population

Disbudding and dehorning are commonly used cattle management practices to protect animals and humans from injury. They are unpleasant, costly processes subject to increased public scrutiny as an animal welfare issue. Horns are a recessively inherited trait, so one option to eliminate dehorning is to breed for polled (hornlessness). However, due to the low genetic merit and scarcity of polled dairy sires, this approach has not been widely adopted. In March 2018, only 3 Holstein and 0 Jersey active homozygous polled sires were registered with the National Association of Animal Breeders. Alternatively, gene editing to produce high-genetic-merit polled sires has been proposed. To further explore this concept, introgression of the POLLED allele into both the US Holstein and Jersey cattle populations via conventional breeding or gene editing (top 1% of bulls/year) was simulated for 3 polled mating schemes and compared with baseline selection on lifetime net merit (NM$) alone, over the course of 20 yr. Scenarios were replicated 10 times and the changes in HORNED allele frequency, inbreeding, genetic gain (NM$), and number of unique sires used were calculated. Gene editing decreased the frequency of the HORNED allele to <0.1 after 20 yr, which was as fast or faster than conventional breeding for both breeds. In the mating scheme that required the use of only existing homozygous polled sires, inbreeding reached 17% (Holstein) and 14% (Jersey), compared with less than 7% in the baseline scenarios. However, gene editing in the same mating scheme resulted in significantly less inbreeding, 9% (Holstein) and 8% (Jersey). Also, gene editing resulted in significantly higher NM$ after 20 yr compared with conventional breeding for both breeds. Additionally, the gene editing scenarios of both breeds used a significantly greater number of unique sires compared with either the conventional breeding or baseline scenarios. Overall, our simulations show that, given the current genetic merit of horned and polled dairy sires, the use of conventional breeding methods to decrease the frequency of the HORNED allele will increase inbreeding and slow genetic improvement. Furthermore, this study demonstrates how gene editing could be used to rapidly decrease the frequency of the HORNED allele in US dairy cattle populations while maintaining the rate of genetic gain, constraining inbreeding to acceptable levels, and simultaneously addressing an emerging animal welfare concern.

Livestock 2.0 - genome editing for fitter, healthier, and more productive farmed animals

The human population is growing, and as a result we need to produce more food whilst reducing the impact of farming on the environment. Selective breeding and genomic selection have had a transformational impact on livestock productivity, and now transgenic and genome-editing technologies offer exciting opportunities for the production of fitter, healthier and more-productive livestock. Here, we review recent progress in the application of genome editing to farmed animal species and discuss the potential impact on our ability to produce food.

Genome-wide association study provides insights into genes related with horn development in Nelore beef cattle

The causal mutation for polledness in Nelore (Bos taurus indicus) breed seems to have appeared first in Brazil in 1957. The expression of the polled trait is known to be ruled by a few groups of alleles in taurine breeds; however, the genetic basis of this trait in indicine cattle is still unclear. The aim of this study was to identify genomic regions associated with the hornless trait in a commercial Nelore population. A total of 107,294 animals had phenotypes recorded and 2,238 were genotyped/imputed for 777k SNP. The weighted single-step approach for genome-wide association study (WssGWAS) was used to estimate the SNP effects and variances accounted for by 1 Mb sliding SNP windows. A centromeric region of chromosome 1 with 3.11 Mb size (BTA1: 878,631–3,987,104 bp) was found to be associated with hornless in the studied population. A total of 28 protein-coding genes are mapped in this region, including the taurine Polled locus and the IFNAR1, IFNAR2, IFNGR2, KRTAP11-1, MIS18A, OLIG1, OLIG2, and SOD1 genes, which expression can be related to the horn formation as described in literature. The functional enrichment analysis by DAVID tool revealed cytokine-cytokine receptor interaction, JAK-STAT signaling, natural killer cell mediated cytotoxicity, and osteoclast differentiation pathways as significant (P < 0.05). In addition, a runs of homozygosity (ROH) analysis identified a ROH island in polled animals with 2.47 Mb inside the region identified by WssGWAS. Polledness in Nelore cattle is associated with one region in the genome with 3.1 Mb size in chromosome 1. Several genes are harbored in this region, and they may act together in the determination of the polled/horned phenotype. Fine mapping the locus responsible for polled trait in Nelore breed and the identification of the molecular mechanisms regulating the horn growth deserve further investigation.

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.

Simultaneous introgression of three POLLED mutations into a synthetic breed of Chinese cattle

The polled phenotype of cattle is increasingly becoming favourable mainly because of the enhanced emphasis on animal welfare, for which the causative mutations have been reported during the past years. The Shuxuan cattle are a new synthetic breed by crossing the indigenous cattle with both Simmental and Holstein semen in Sichuan of Southwest China, in which about 15% of polled individuals have newly emerged. Because official record about POLLED genotypes for the historically imported sires is unavailable, we therefore genotyped the proposed POLLED variants of P_202ID, P_80kbID and P_219ID among 48 polled and 16 horned Shuxuan cattle. It was first revealed that all three candidate mutations have been simultaneously introgressed into Shuxuan cattle, whereas the P_202ID mutation is dominant. Furthermore, one polled animal still remains to carry none of the three candidate mutations, which suggests that further mutation(s) would also exist. Additionally, we sequenced mitochondrial DNA and found that Shuxuan cattle are composed of two matrilineal origins of Bos taurus (65.6%) and B. indicus (34.4%); and there is no origin-biased distribution of polled phenotype. In conclusion, our study first supports the recently reported novel candidate mutation of P_219ID and detects simultaneous presences of all three known POLLED mutations within a cattle breed.

Whole-genome analysis of introgressive hybridization and characterization of the bovine legacy of Mongolian yaks

The yak is remarkable for its adaptation to high altitude and occupies a central place in the economies of the mountainous regions of Asia. At lower elevations, it is common to hybridize yaks with cattle to combine the yak's hardiness with the productivity of cattle. Hybrid males are sterile, however, preventing the establishment of stable hybrid populations, but not a limited introgression after backcrossing several generations of female hybrids to male yaks. Here we inferred bovine haplotypes in the genomes of 76 Mongolian yaks using high-density SNP genotyping and whole-genome sequencing. These yaks inherited ∼1.3% of their genome from bovine ancestors after nearly continuous admixture over at least the last 1,500 years. The introgressed regions are enriched in genes involved in nervous system development and function, and particularly in glutamate metabolism and neurotransmission. We also identified a novel mutation associated with a polled (hornless) phenotype originating from Mongolian Turano cattle. Our results suggest that introgressive hybridization contributed to the improvement of yak management and breeding.

Detection and validation of structural variations in bovine whole-genome sequence data

BACKGROUND

Several examples of structural variation (SV) affecting phenotypic traits have been reported in cattle. Currently the identification of SV from whole-genome sequence data (WGS) suffers from a high false positive rate. Our aim was to construct a high quality set of SV calls in cattle using WGS data. First, we tested two SV detection programs, Breakdancer and Pindel, and the overlap of these methods, on simulated sequence data to determine their precision and sensitivity. We then identified population SV from WGS of 252 Holstein and 64 Jersey bulls based on the overlapping calls from the two programs. In addition, we validated an overlapped SV set in 28 twice-sequenced Holstein individuals, and in another two validated sets (one for each breed) that were transmitted from sire to son. We also tested whether highly conserved gene sets across eukaryotes and recently expanded gene families in bovine were depleted and enriched, respectively, for SV.

RESULTS

In empirical WGS data, 17,518 SV covering 27.36 Mb were found in the Holstein population and 4285 SV covering 8.74 Mb in the Jersey population, of which 4.62 Mb of SV overlapped between Holsteins and Jerseys. A total of 11,534 candidate SV covering 5.64 Mb were validated in the 28 twice-sequenced individuals, while 3.49 and 0.67 Mb of SV were validated from Holstein and Jersey sire-son transmission, respectively. Only eight of 237 core eukaryotic genes had at least a 50-bp overlap with an SV from our validated sets, suggesting that conserved genes are depleted for SV (p < 0.05). In addition, we observed that recently expanded gene families were significantly more associated with SV than other genes. Long interspersed nuclear elements-1 were enriched for deletions when compared to the rest of the genome (p = 0.0035).

CONCLUSIONS

We reported SV from 252 Holstein and 64 Jersey individuals. A considerable proportion of Jersey population SV (53.5%) were also found in Holstein. In contrast, about 76.90% sire-son transmission validated SV were present in Jerseys and Holsteins. The enrichment of SV in expanding gene families suggests that SV can be a source of genetic variation for evolution.

Genome-wide Association Study Identifies Loci for the Polled Phenotype in Yak

The absence of horns, known as the polled phenotype, is an economically important trait in modern yak husbandry, but the genomic structure and genetic basis of this phenotype have yet to be discovered. Here, we conducted a genome-wide association study with a panel of 10 horned and 10 polled yaks using whole genome sequencing. We mapped the POLLED locus to a 200-kb interval, which comprises three protein-coding genes. Further characterization of the candidate region showed recent artificial selection signals resulting from the breeding process. We suggest that expressional variations rather than structural variations in protein probably contribute to the polled phenotype. Our results not only represent the first and important step in establishing the genomic structure of the polled region in yak, but also add to our understanding of the polled trait in bovid species.

Evaluation of breeding strategies for polledness in dairy cattle using a newly developed simulation framework for quantitative and Mendelian traits

Background

Intensified selection of polled individuals has recently gained importance in predominantly horned dairy cattle breeds as an alternative to routine dehorning. The status quo of the current polled breeding pool of genetically-closely related artificial insemination sires with lower breeding values for performance traits raises questions regarding the effects of intensified selection based on this founder pool.

Methods

We developed a stochastic simulation framework that combines the stochastic simulation software QMSim and a self-designed R program named QUALsim that acts as an external extension. Two traits were simulated in a dairy cattle population for 25 generations: one quantitative (QMSim) and one qualitative trait with Mendelian inheritance (i.e. polledness, QUALsim). The assignment scheme for qualitative trait genotypes initiated realistic initial breeding situations regarding allele frequencies, true breeding values for the quantitative trait and genetic relatedness. Intensified selection for polled cattle was achieved using an approach that weights estimated breeding values in the animal best linear unbiased prediction model for the quantitative trait depending on genotypes or phenotypes for the polled trait with a user-defined weighting factor.

Results

Selection response for the polled trait was highest in the selection scheme based on genotypes. Selection based on phenotypes led to significantly lower allele frequencies for polled. The male selection path played a significantly greater role for a fast dissemination of polled alleles compared to female selection strategies. Fixation of the polled allele implies selection based on polled genotypes among males. In comparison to a base breeding scenario that does not take polledness into account, intensive selection for polled substantially reduced genetic gain for this quantitative trait after 25 generations. Reducing selection intensity for polled males while maintaining strong selection intensity among females, simultaneously decreased losses in genetic gain and achieved a final allele frequency of 0.93 for polled.

Conclusions

A fast transition to a completely polled population through intensified selection for polled was in contradiction to the preservation of high genetic gain for the quantitative trait. Selection on male polled genotypes with moderate weighting, and selection on female polled phenotypes with high weighting, could be a suitable compromise regarding all important breeding aspects.

Electronic supplementary material

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

Polyceraty (multi-horns) in Damara sheep maps to ovine chromosome 2

Polyceraty (presence of multiple horns) is rare in modern day ungulates. Although not found in wild sheep, polyceraty does occur in a small number of domestic sheep breeds covering a wide geographical region. Damara are fat-tailed hair sheep, from the south-western region of Africa, which display polyceraty, with horn number ranging from zero to four. We conducted a genome-wide association study for horn number with 43 Damara genotyped with 606 006 SNP markers. The analysis revealed a region with multiple significant SNPs on ovine chromosome 2, in a location different from the mutation for polled in sheep on chromosome 10. The causal mutation for polyceraty was not identified; however, the region associated with polyceraty spans nine HOXD genes, which are critical in embryonic development of appendages. Mutations in HOXD genes are implicated in polydactly phenotypes in mice and humans. There was no evidence for epistatic interactions contributing to polyceraty. This is the first report on the genetic mechanisms underlying polyceraty in the under-studied Damara.

Considering genetic characteristics in German Holstein breeding programs

Recently, several research groups have demonstrated that several haplotypes may cause embryonic loss in the homozygous state. Up to now, carriers of genetic disorders were often excluded from mating, resulting in a decrease of genetic gain and a reduced number of sires available for the breeding program. Ongoing research is very likely to identify additional genetic defects causing embryonic loss and calf mortality by genotyping a large proportion of the female cattle population and sequencing key ancestors. Hence, a clear demand is present to develop a method combining selection against recessive defects (e.g., Holstein haplotypes HH1-HH5) with selection for economically beneficial traits (e.g., polled) for mating decisions. Our proposed method is a genetic index that accounts for the allele frequencies in the population and the economic value of the genetic characteristic without excluding carriers from breeding schemes. Fertility phenotypes from routine genetic evaluations were used to determine the economic value per embryo lost. Previous research has shown that embryo loss caused by HH1 and HH2 occurs later than the loss for HH3, HH4, and HH5. Therefore, an economic value of € 97 was used against HH1 and HH2 and € 70 against HH3, HH4, and HH5. For polled, € 7 per polled calf was considered. Minor allele frequencies of the defects ranged between 0.8 and 3.3%. The polled allele has a frequency of 4.1% in the German Holstein population. A genomic breeding program was simulated to study the effect of changing the selection criteria from assortative mating based on breeding values to selecting the females using the genetic index. Selection for a genetic index on the female path is a useful method to control the allele frequencies by reducing undesirable alleles and simultaneously increasing economical beneficial characteristics maintaining most of the genetic gain in production and functional traits. Additionally, we applied the genetic index to real data and found a decrease of the genetic trend for the birth years 1990 to 2006. Since 2010 the genetic index has increased due to a strong increase in the polled frequency. However, further investigation is needed to better understand the biology to determine the correct time of embryo loss and the economic value of fertility disorders.

A 1.8-kb insertion in the 3'-UTR of RXFP2 is associated with polledness in sheep

Sheep breeds show a broad spectrum of different horn phenotypes. In most modern production breeds, sheep are polled (absence of horns), whereas horns occur mainly in indigenous breeds. Previous studies mapped the responsible locus to the region of the RXFP2 gene on ovine chromosome 10. A 4-kb region of the 3'-end of RXFP2 was amplified in horned and polled animals from seven Swiss sheep breeds. Sequence analysis identified a 1833-bp genomic insertion located in the 3'-UTR region of RXFP2 present in polled animals only. An efficient PCR-based genotyping method to determine the polled genotype of individual sheep is presented. Comparative sequence analyses revealed evidence that the polled-associated insertion adds a potential antisense RNA sequence of EEF1A1 to the 3'-end of RXFP2 transcripts.

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.

A genomic duplication is associated with ectopic eomesodermin expression in the embryonic chicken comb and two duplex-comb phenotypes

Duplex-comb (D) is one of three major loci affecting comb morphology in the domestic chicken. Here we show that the two Duplex-comb alleles, V-shaped (D*V) and Buttercup (D*C), are both associated with a 20 Kb tandem duplication containing several conserved putative regulatory elements located 200 Kb upstream of the eomesodermin gene (EOMES). EOMES is a T-box transcription factor that is involved in mesoderm specification during gastrulation. In D*V and D*C chicken embryos we find that EOMES is ectopically expressed in the ectoderm of the comb-developing region as compared to wild-type embryos. The confinement of the ectopic expression of EOMES to the ectoderm is in stark contrast to the causal mechanisms underlying the two other major comb loci in the chicken (Rose-comb and Pea-comb) in which the transcription factors MNR2 and SOX5 are ectopically expressed strictly in the mesenchyme. Interestingly, the causal mutations of all three major comb loci in the chicken are now known to be composed of large-scale structural genomic variants that each result in ectopic expression of transcription factors. The Duplex-comb locus also illustrates the evolution of alleles in domestic animals, which means that alleles evolve by the accumulation of two or more consecutive mutations affecting the phenotype. We do not yet know whether the V-shaped or Buttercup allele correspond to the second mutation that occurred on the haplotype of the original duplication event.