Meta-analysis of six dairy cattle breeds reveals biologically relevant candidate genes for mastitis resistance.

BACKGROUND: Mastitis is a disease that incurs significant costs in the dairy industry. A promising approach to mitigate its negative effects is to genetically improve the resistance of dairy cattle to mastitis. A meta-analysis of genome-wide association studies (GWAS) across multiple breeds for clinical mastitis (CM) and its indicator trait, somatic cell score (SCS), is a powerful method to identify functional genetic variants that impact mastitis resistance. RESULTS: We conducted meta-analyses of eight and fourteen GWAS on CM and SCS, respectively, using 30,689 and 119,438 animals from six dairy cattle breeds. Methods for the meta-analyses were selected to properly account for the multi-breed structure of the GWAS data. Our study revealed 58 lead markers that were associated with mastitis incidence, including 16 loci that did not overlap with previously identified quantitative trait loci (QTL), as curated at the Animal QTLdb. Post-GWAS analysis techniques such as gene-based analysis and genomic feature enrichment analysis enabled prioritization of 31 candidate genes and 14 credible candidate causal variants that affect mastitis. CONCLUSIONS: Our list of candidate genes can help to elucidate the genetic architecture underlying mastitis resistance and provide better tools for the prevention or treatment of mastitis, ultimately contributing to more sustainable animal production.

Sequence-based GWAS meta-analyses for beef production traits.

BACKGROUND: Combining the results of within-population genome-wide association studies (GWAS) based on whole-genome sequences into a single meta-analysis (MA) is an accurate and powerful method for identifying variants associated with complex traits. As part of the H2020 BovReg project, we performed sequence-level MA for beef production traits. Five partners from France, Switzerland, Germany, and Canada contributed summary statistics from sequence-based GWAS conducted with 54,782 animals from 15 purebred or crossbred populations. We combined the summary statistics for four growth, nine morphology, and 15 carcass traits into 16 MA, using both fixed effects and z-score methods. RESULTS: The fixed-effects method was generally more informative to provide indication on potentially causal variants, although we combined substantially different traits in each MA. In comparison with within-population GWAS, this approach highlighted (i) a larger number of quantitative trait loci (QTL), (ii) QTL more frequently located in genomic regions known for their effects on growth and meat/carcass traits, (iii) a smaller number of genomic variants within the QTL, and (iv) candidate variants that were more frequently located in genes. MA pinpointed variants in genes, including MSTN, LCORL, and PLAG1 that have been previously associated with morphology and carcass traits. We also identified dozens of other variants located in genes associated with growth and carcass traits, or with a function that may be related to meat production (e.g., HS6ST1, HERC2, WDR75, COL3A1, SLIT2, MED28, and ANKAR). Some of these variants overlapped with expression or splicing QTL reported in the cattle Genotype-Tissue Expression atlas (CattleGTEx) and could therefore regulate gene expression. CONCLUSIONS: By identifying candidate genes and potential causal variants associated with beef production traits in cattle, MA demonstrates great potential for investigating the biological mechanisms underlying these traits. As a complement to within-population GWAS, this approach can provide deeper insights into the genetic architecture of complex traits in beef cattle.

Activation of cryptic splicing in bovine WDR19 is associated with reduced semen quality and male fertility.

Cattle are ideally suited to investigate the genetics of male reproduction, because semen quality and fertility are recorded for all ejaculates of artificial insemination bulls. We analysed 26,090 ejaculates of 794 Brown Swiss bulls to assess ejaculate volume, sperm concentration, sperm motility, sperm head and tail anomalies and insemination success. The heritability of the six semen traits was between 0 and 0.26. Genome-wide association testing on 607,511 SNPs revealed a QTL on bovine chromosome 6 that was associated with sperm motility (P = 2.5 x 10-27), head (P = 2.0 x 10-44) and tail anomalies (P = 7.2 x 10-49) and insemination success (P = 9.9 x 10-13). The QTL harbors a recessive allele that compromises semen quality and male fertility. We replicated the effect of the QTL on fertility (P = 7.1 x 10-32) in an independent cohort of 2481 Brown Swiss bulls. The analysis of whole-genome sequencing data revealed that a synonymous variant (BTA6:58373887C>T, rs474302732) in WDR19 encoding WD repeat-containing protein 19 was in linkage disequilibrium with the fertility-associated haplotype. WD repeat-containing protein 19 is a constituent of the intraflagellar transport complex that is essential for the physiological function of motile cilia and flagella. Bioinformatic and transcription analyses revealed that the BTA6:58373887 T-allele activates a cryptic exonic splice site that eliminates three evolutionarily conserved amino acids from WDR19. Western blot analysis demonstrated that the BTA6:58373887 T-allele decreases protein expression. We make the remarkable observation that, in spite of negative effects on semen quality and bull fertility, the BTA6:58373887 T-allele has a frequency of 24% in the Brown Swiss population. Our findings are the first to uncover a variant that is associated with quantitative variation in semen quality and male fertility in cattle.

A 1-bp deletion in bovine QRICH2 causes low sperm count and immotile sperm with multiple morphological abnormalities.

BACKGROUND: Semen quality and insemination success are monitored in artificial insemination bulls to ensure high male fertility rates. Only ejaculates that fulfill minimum quality requirements are processed and eventually used for artificial inseminations. We examined 70,990 ejaculates from 1343 Brown Swiss bulls to identify bulls from which all ejaculates were rejected due to low semen quality. This procedure identified a bull that produced 12 ejaculates with an aberrantly small number of sperm (0.2 ± 0.2 × 109 sperm per mL) which were mostly immotile due to multiple morphological abnormalities. RESULTS: The genome of this bull was sequenced at a 12× coverage to investigate a possible genetic cause. Comparing the sequence variant genotypes of this bull with those from 397 fertile bulls revealed a 1-bp deletion in the coding sequence of the QRICH2 gene which encodes the glutamine rich 2 protein, as a compelling candidate causal variant. This 1-bp deletion causes a frameshift in translation and a premature termination codon (ENSBTAP00000018337.1:p.Cys1644AlafsTer52). The analysis of testis transcriptomes from 76 bulls showed that the transcript with the premature termination codon is subject to nonsense-mediated mRNA decay. The 1-bp deletion resides in a 675-kb haplotype that includes 181 single nucleotide polymorphisms (SNPs) from the Illumina BovineHD Bead chip. This haplotype segregates at a frequency of 5% in the Brown Swiss cattle population. Our analysis also identified another bull that carried the 1-bp deletion in the homozygous state. Semen analyses from the second bull confirmed low sperm concentration and immotile sperm with multiple morphological abnormalities that primarily affect the sperm flagellum and, to a lesser extent, the sperm head. CONCLUSIONS: A recessive loss-of-function allele of the bovine QRICH2 gene likely causes low sperm concentration and immotile sperm with multiple morphological abnormalities. Routine sperm analyses unambiguously identify homozygous bulls for this allele. A direct gene test can be implemented to monitor the frequency of the undesired allele in cattle populations.

Associated regions for multiple birth in Brown Swiss and Original Braunvieh cattle on chromosomes 15 and 11.

Twin and multiple births have negative effects on the performance and health of cows and calves. To decipher the genetic architecture of this trait in the two Swiss Brown Swiss cattle populations, we performed various association analyses based on de-regressed breeding values. Genome-wide association analyses were executed using ~600 K imputed SNPs for the maternal multiple birth trait in ~3500 Original Braunvieh and ~7800 Brown Swiss animals. Significantly associated QTL were observed on different chromosomes for both breeds. We have identified on chromosome 11 a QTL that explains ~6% of the total genetic variance of the maternal multiple birth trait in Original Braunvieh. For the Brown Swiss breed, we have discovered a QTL on chromosome 15 that accounts for ~4% of the total genetic variance. For Original Braunvieh, subsequent haplotype analysis revealed a 90-kb window on chromosome 11 at 88 Mb, where a likely regulatory region is located close to the ID2 gene. In Brown Swiss, a 130-kb window at 75 Mb on chromosome 15 was identified. Analysis of whole-genome sequence data using linkage-disequilibrium estimation revealed possible causal variants for the identified QTL. A presumably regulatory variant in the non-coding 5' region of the ID2 gene was strongly associated with the haplotype for Original Braunvieh. In Brown Swiss, an intron variant in PRDM11, one 3' UTR variant in SYT13 and three intergenic variants 5' upstream of SYT13 were identified as candidate variants for the trait multiple birth maternal. In this study, we report for the first time QTL for the trait of multiple births in Original Braunvieh and Brown Swiss cattle. Moreover, our findings are another step towards a better understanding of the complex genetic architecture of this polygenic trait.

Mining massive genomic data of two Swiss Braunvieh cattle populations reveals six novel candidate variants that impair reproductive success.

BACKGROUND: This study was carried out on the two Braunvieh populations reared in Switzerland, the dairy Brown Swiss (BS) and the dual-purpose Original Braunvieh (OB). We performed a genome-wide analysis of array data of trios (sire, dam, and offspring) from the routine genomic selection to identify candidate regions showing missing homozygosity and phenotypic associations with five fertility, ten birth, and nine growth-related traits. In addition, genome-wide single SNP regression studies based on 114,890 single nucleotide polymorphisms (SNPs) for each of the two populations were performed. Furthermore, whole-genome sequencing data of 430 cattle including 70 putative haplotype carriers were mined to identify potential candidate variants that were validated by genotyping the current population using a custom array. RESULTS: Using a trio-based approach, we identified 38 haplotype regions for BS and five for OB that segregated at low to moderate frequencies. For the BS population, we confirmed two known haplotypes, BH1 and BH2. Twenty-four variants that potentially explained the missing homozygosity and associated traits were detected, in addition to the previously reported TUBD1:p.His210Arg variant associated with BH2. For example, for BS we identified a stop-gain variant (p.Arg57*) in the MRPL55 gene in the haplotype region on chromosome 7. This region is associated with the 'interval between first and last insemination' trait in our data, and the MRPL55 gene is known to be associated with early pregnancy loss in mice. In addition, we discuss candidate missense variants in the CPT1C, MARS2, and ACSL5 genes for haplotypes mapped in BS. In OB, we highlight a haplotype region on chromosome 19, which is potentially caused by a frameshift variant (p.Lys828fs) in the LIG3 gene, which is reported to be associated with early embryonic lethality in mice. Furthermore, we propose another potential causal missense variant in the TUBGCP5 gene for a haplotype mapped in OB. CONCLUSIONS: We describe, for the first time, several haplotype regions that segregate at low to moderate frequencies and provide evidence of causality by trait associations in the two populations of Swiss Braunvieh. We propose a list of six protein-changing variants as potentially causing missing homozygosity. These variants need to be functionally validated and incorporated in the breeding program.

A major QTL at the LHCGR/FSHR locus for multiple birth in Holstein cattle.

BACKGROUND: Twin and multiple births are rare in cattle and have a negative impact on the performance and health of cows and calves. Therefore, selection against multiple birth would be desirable in dairy cattle breeds such as Holstein. We applied different methods to decipher the genetic architecture of this trait using de-regressed breeding values for maternal multiple birth of ~ 2500 Holstein individuals to perform genome-wide association analyses using ~ 600 K imputed single nucleotide polymorphisms (SNPs). RESULTS: In the population studied, we found no significant genetic trend over time of the estimated breeding values for multiple birth, which indicates that this trait has not been selected for in the past. In addition to several suggestive non-significant quantitative trait loci (QTL) on different chromosomes, we identified a major QTL on chromosome 11 for maternal multiple birth that explains ~ 16% of the total genetic variance. Using a haplotype-based approach, this QTL was fine-mapped to a 70-kb window on chromosome 11 between 31.00 and 31.07 Mb that harbors two functional candidate genes (LHCGR and FSHR). Analysis of whole-genome sequence data by linkage-disequilibrium estimation revealed a regulatory variant in the 5'-region of LHCGR as a possible candidate causal variant for the identified major QTL. Furthermore, the identified haplotype showed significant effects on stillbirth and days to first service. CONCLUSIONS: QTL detection and subsequent identification of causal variants in livestock species remain challenging in spite of the availability of large-scale genotype and phenotype data. Here, we report for the first time a major QTL for multiple birth in Holstein cattle and provide evidence for a linked variant in the non-coding region of a functional candidate gene. This discovery, which is a first step towards the understanding of the genetic architecture of this polygenic trait, opens the path for future selection against this undesirable trait, and thus contributes to increased animal health and welfare.

CNGB3 Missense Variant Causes Recessive Achromatopsia in Original Braunvieh Cattle.

Sporadic occurrence of inherited eye disorders has been reported in cattle but so far pathogenic variants were found only for rare forms of cataract but not for retinopathies. The aim of this study was to characterize the phenotype and the genetic aetiology of a recessive form of congenital day-blindness observed in several cases of purebred Original Braunvieh cattle. Electroretinography in an affected calf revealed absent cone-mediated function, whereas the rods continue to function normally. Brain areas involved in vision were morphologically normal. When targeting cones by immunofluorescence, a decrease in cone number and an accumulation of beta subunits of cone cyclic-nucleotide gated channel (CNGB3) in the outer plexiform layer of affected animals was obvious. Achromatopsia is a monogenic Mendelian disease characterized by the loss of cone photoreceptor function resulting in day-blindness, total color-blindness, and decreased central visual acuity. After SNP genotyping and subsequent homozygosity mapping with twelve affected cattle, we performed whole-genome sequencing and variant calling of three cases. We identified a single missense variant in the bovine CNGB3 gene situated in a ~2.5 Mb homozygous genome region on chromosome 14 shared between all cases. All affected cattle were homozygous carriers of the p.Asp251Asn mutation that was predicted to be deleterious, affecting an evolutionary conserved residue. In conclusion, we have evidence for the occurrence of a breed-specific novel CNGB3-related form of recessively inherited achromatopsia in Original Braunvieh cattle which we have designated OH1 showing an allele frequency of the deleterious allele of ~8%. The identification of carriers will enable selection against this inherited disorder. The studied cattle might serve as an animal model to further elucidate the function of CNGB3 in mammals.

Christmas disease in a Hovawart family resembling human hemophilia B Leyden is caused by a single nucleotide deletion in a highly conserved transcription factor binding site of the F9 gene promoter.

Hemophilia B is a classical monogenic, X-chromosomal, recessively transmitted bleeding disorder caused by genetic variants within the coagulation factor IX gene (F9). Although hemophilia B has been described in dogs, it has not yet been reported in the Hovawart breed. Here we describe the identification of a Hovawart family transmitting typical signs of an X-linked bleeding disorder. Five males were reported to suffer from recurrent hemorrhagic episodes. A blood sample from one of these males with only 2% of the normal concentration of plasma factor IX together with samples from seven relatives were provided. Next-generation sequencing of the mother and grandmother revealed a single nucleotide deletion in the F9 promoter. Genotyping of the deletion in 1,298 dog specimens including 720 Hovawarts revealed that the mutant allele was only present in the aforementioned Hovawart family. The deletion is located 73 bp upstream of the F9 start codon in the conserved overlapping DNA binding sites of hepatocyte nuclear factor 4α (HNF-4α) and androgen receptor (AR). The deletion only abolished binding of HNF-4α, while AR binding was unaffected as demonstrated by electrophoretic mobility shift assay using human HNF-4α and AR with double-stranded DNA probes encompassing the mutant promoter region. Luciferase reporter assays using wildtype and mutated promoter fragment constructs transfected into Hep G2 cells showed a significant reduction in expression from the mutant promoter. The data provide evidence that the deletion in the Hovawart family caused a rare type of hemophilia B resembling human hemophilia B Leyden.

Four novel candidate causal variants for deficient homozygous haplotypes in Holstein cattle.

Mendelian variants can determine both insemination success and neonatal survival and thus influence fertility and rearing success of cattle. We present 24 deficient homozygous haplotype regions in the Holstein population of Switzerland and provide an overview of the previously identified haplotypes in the global Holstein breed. This study encompasses massive genotyping, whole-genome sequencing (WGS) and phenotype association analyses. We performed haplotype screenings on almost 53 thousand genotyped animals including 114 k SNP data with two different approaches. We revealed significant haplotype associations to several survival, birth and fertility traits. Within haplotype regions, we mined WGS data of hundreds of bovine genomes for candidate causal variants, which were subsequently evaluated by using a custom genotyping array in several thousand breeding animals. With this approach, we confirmed the known deleterious SMC2:p.Phe1135Ser missense variant associated with Holstein haplotype (HH) 3. For two previously reported deficient homozygous haplotypes that show negative associations to female fertility traits, we propose candidate causative loss-of-function variants: the HH13-related KIR2DS1:p.Gln159* nonsense variant and the HH21-related NOTCH3:p.Cys44del deletion. In addition, we propose the RIOX1:p.Ala133_Glu142del deletion as well as the PCDH15:p.Leu867Val missense variant to explain the unexpected low number of homozygous haplotype carriers for HH25 and HH35, respectively. In conclusion, we demonstrate that with mining massive SNP data in combination with WGS data, we can map several haplotype regions and unravel novel recessive protein-changing variants segregating at frequencies of 1 to 5%. Our findings both confirm previously identified loci and expand the spectrum of undesired alleles impairing reproduction success in Holstein cattle, the world's most important dairy breed.