A nonsense mutation in PLD4 is associated with a zinc deficiency-like syndrome in Fleckvieh cattle

Detecting and characterizing copy number variation in a large commercial U.S. Holstein cattle population

BACKGROUND

Copy number variations (CNVs) are an important source of genomic variation that play an active role in modulating biological processes by altering gene expression and dosage. These structural variants involve duplications and deletions of segments usually exceeding 1 kilobase in size, dispersed throughout the genome of humans and livestock individuals. We mapped CNVs from high-density single-nucleotide polymorphism (SNP) genotyping array data on 3,601 Holsteins. Following, we explored their association with reported quantitative trait loci (QTLs), genes, and biological pathways, unveiling the potential biological contributions of CNVs to economically important traits in the dairy industry and breeding programs.

RESULTS

We identified 4,113 non-redundant high-confidence CNVs, of which 78% were deletions and 22% duplications, distributed across all bovine autosomal chromosomes (BTA). Out of the 1,184 compiled CNV regions (CNVRs) covering 3.02% of the autosomal genome, 199 novel CNVRs were mapped. QTLs overlapping with CNVRs detected in this study were enriched for 140 economically important traits, such as milk yield, conception and pregnancy rates, susceptibility to diseases and length of productive life, indicating that CNVs likely underlie productive, reproductive and health performance in Holstein dairy cattle. Moreover, detected CNVRs overlapped with 2,788 annotated genes, including those linked to milk production, fertility, and immune response in cattle, such as DGAT1, AFF1, and ADAMTS13 genes. Furthermore, the gene set analysis revealed GO terms related to metabolic processes, immune system, response to stimulus, and cellular binding activities. Notably, enriched GO terms highlighted relevant genes to cattle health and reproduction overlapping CNVRs, such as DEFB4A, GATA3, GNB1, and PPP1R1B.

CONCLUSIONS

We mapped and demonstrated the characteristics of genome-wide distributed CNVs in a large commercial Holstein population genotyped with a high-density SNP array. Collectively, the results emphasize the role of CNVs as a valuable resource of genetic variation and contribute to better understand the genetic architecture of economic complex traits in dairy cattle. Furthermore, these findings may provide opportunities for the development of novel and enhanced genomic selection strategies in Holstein cattle.

Structural and mechanistic insights into disease-associated endolysosomal exonucleases PLD3 and PLD4

Endolysosomal exonucleases PLD3 and PLD4 (phospholipases D3 and D4) are associated with autoinflammatory and autoimmune diseases. We report structures of these enzymes, and the molecular basis of their catalysis. The structures reveal an intra-chain dimer topology forming a basic active site at the interface. Like other PLD superfamily members, PLD3 and PLD4 carry HxKxxxxD/E motifs and participate in phosphodiester-bond cleavage. The enzymes digest ssDNA and ssRNA in a 5'-to-3' manner and are blocked by 5'-phosphorylation. We captured structures in apo, intermediate, and product states and revealed a "link-and-release" two-step catalysis. We also unexpectedly demonstrated phosphatase activity via a covalent 3-phosphohistidine intermediate. PLD4 contains an extra hydrophobic clamp that stabilizes substrate and could affect oligonucleotide substrate preference and product release. Biochemical and structural analysis of disease-associated mutants of PLD3/4 demonstrated reduced enzyme activity or thermostability and the possible basis for disease association. Furthermore, these findings provide insight into therapeutic design.

Structural and mechanistic insights into disease-associated endolysosomal exonucleases PLD3 and PLD4

Endolysosomal exonucleases PLD3 and PLD4 (phospholipases D3 and D4) are associated with autoinflammatory and autoimmune diseases. We report structures of these enzymes, and the molecular basis of their catalysis. The structures reveal an intra-chain dimer topology forming a basic active site at the interface. Like other PLD superfamily members, PLD3 and PLD4 carry HxKxxxxD/E motifs and participate in phosphodiester-bond cleavage. The enzymes digest ssDNA and ssRNA in a 5'-to-3' manner and are blocked by 5'-phosphorylation. We captured structures in apo, intermediate, and product states and revealed a 'link-and-release' two-step catalysis. We also unexpectedly demonstrated phosphatase activity via a covalent 3' phosphistidine intermediate. PLD4 contains an extra hydrophobic clamp that stabilizes substrate and could affect oligonucleotide substrate preference and product release. Biochemical and structural analysis of disease-associated mutants of PLD3/4 demonstrated reduced enzyme activity or thermostability and the possible basis for disease association. Furthermore, these findings provide insight into therapeutic design.

Chromosome-length genome assembly and linkage map of a critically endangered Australian bird: the helmeted honeyeater

BACKGROUND

The helmeted honeyeater (Lichenostomus melanops cassidix) is a Critically Endangered bird endemic to Victoria, Australia. To aid its conservation, the population is the subject of genetic rescue. To understand, monitor, and modulate the effects of genetic rescue on the helmeted honeyeater genome, a chromosome-length genome and a high-density linkage map are required.

RESULTS

We used a combination of Illumina, Oxford Nanopore, and Hi-C sequencing technologies to assemble a chromosome-length genome of the helmeted honeyeater, comprising 906 scaffolds, with length of 1.1 Gb and scaffold N50 of 63.8 Mb. Annotation comprised 57,181 gene models. Using a pedigree of 257 birds and 53,111 single-nucleotide polymorphisms, we obtained high-density linkage and recombination maps for 25 autosomes and Z chromosome. The total sex-averaged linkage map was 1,347 cM long, with the male map being 6.7% longer than the female map. Recombination maps revealed sexually dimorphic recombination rates (overall higher in males), with average recombination rate of 1.8 cM/Mb. Comparative analyses revealed high synteny of the helmeted honeyeater genome with that of 3 passerine species (e.g., 32 Hi-C scaffolds mapped to 30 zebra finch autosomes and Z chromosome). The genome assembly and linkage map suggest that the helmeted honeyeater exhibits a fission of chromosome 1A into 2 chromosomes relative to zebra finch. PSMC analysis showed a ∼15-fold decline in effective population size to ∼60,000 from mid- to late Pleistocene.

CONCLUSIONS

The annotated chromosome-length genome and high-density linkage map provide rich resources for evolutionary studies and will be fundamental in guiding conservation efforts for the helmeted honeyeater.

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.

Macrophages in dermal disease progression of phospholipase D4-deficient Fleckvieh calves

A new gene defect in Fleckvieh calves leads to a syndrome with partial phenotype overlap with bovine hereditary zinc deficiency. A mutation in a gene encoding phospholipase D4 (PLD4), an endosomal exonuclease, causes the disorder. In mice, PLD4 activity indirectly regulates the Toll-like receptor 9 (TLR9) pathway via degradation of microbial DNA. PLD4 absence thus results in visceral macrophage activation comparable to human macrophage activation syndrome. In this study, disease progression and the role of macrophages in affected calves were monitored clinically, clinicopathologically, and histologically over time. Breeding data identified 73 risk matings of heterozygous carriers resulting in 54 potentially PLD4-deficient calves born on farms. PLD4 status was examined via 5'-exonuclease assay, detecting 6 calves carrying the defect. These were purchased and monitored daily until final necropsy. The calves developed progressive skin lesions starting with small scaling areas terminating in severe crusting dermatitis, especially in areas with mechanical exposure. Histological and immunohistochemical analyses indicated that macrophages with cytoplasmic vacuolation increased considerably in skin sections obtained weekly during the disease course. Macrophage increase correlated with increased dermal lesion severity. Macrophage activation was confirmed by prominent phagocytic activity in the superficial dermis using electron microscopy. Dermal mRNA abundance of CCL2 and CCL3 measured by quantitative polymerase chain reaction verified macrophage activation. Further increase in mRNA of downstream molecule MyD88 and cytokine IL12b connected bovine PLD4 deficiency to increased TLR9 pathway activation. In contrast to human macrophage activation syndrome, the main feature of bovine PLD4 deficiency was local disease in organs with contact to microbial DNA (skin, intestine, lungs).

Reverse Genetic Screen for Deleterious Recessive Variants in the Local Simmental Cattle Population of Switzerland

We herein report the result of a large-scale reverse genetic screen in the Swiss Simmental population, a local dual-purpose cattle breed. We aimed to detect possible recessively inherited variants affecting protein-coding genes, as such deleterious variants can impair fertility and rearing success significantly. We used 115,000 phased SNP data of almost 10 thousand cattle with pedigree data. This revealed evidence for 11 genomic regions of 1.17 Mb on average, with haplotypes (SH1 to SH11) showing a significant depletion in homozygosity and an allele frequency between 3.2 and 10.6%. For the proposed haplotypes, it was unfortunately not possible to evaluate associations with fertility traits as no corresponding data were available. For each haplotype region, possible candidate genes were listed based on their known function in development and disease. Subsequent mining of single-nucleotide variants and short indels in the genomes of 23 sequenced haplotype carriers allowed us to identify three perfectly linked candidate causative protein-changing variants: a SH5-related DIS3:p.Ile678fs loss-of-function variant, a SH8-related CYP2B6:p.Ile313Asn missense variant, and a SH9-related NUBPL:p.Ser143Tyr missense variant. None of these variants occurred in homozygous state in any of more than 5200 sequenced cattle of various breeds. Selection against these alleles in order to reduce reproductive failure and animal loss is recommended.

A 44-kb deleted-type copy number variation is associated with decreasing complement component activity and calf mortality in Japanese Black cattle

BACKGROUND

Calf mortality generally occurs in calves prior to weaning, which is a serious problem in cattle breeding. Several causative variants of monogenic Mendelian disorders in calf mortality have been identified, whereas genetic factors affecting the susceptibility of calves to death are not well known. To identify variants associated with calf mortality in Japanese Black cattle, we evaluated calf mortality as a categorical trait with a threshold model and performed a genome-wide copy number variation (CNV) association study on calf mortality.

RESULTS

We identified a 44-kb deleted-type CNV ranging from 103,317,687 to 103,361,802 bp on chromosome 5, which was associated with the mortality of 1-180-day-old calves. The CNV harbored C1RL, a pseudogene, and an IncRNA localized in the C1R and C1S gene cluster, which is a component of the classical complement activation pathway for immune complexes for infectious pathogens. The average complement activity in CNVR_221 homozygotes at postnatal day 7 was significantly lower than that of wild-type animals and heterozygotes. The frequency of the risk allele in dead calves suffering from diarrhea and pneumonia and in healthy cows was 0.35 and 0.28, respectively (odds ratio = 2.2, P = 0.016), suggesting that CNVR_221 was associated with the mortality of Japanese Black calves suffering from an infectious disease.

CONCLUSIONS

This study identified a deleted-type CNV associated with the mortality of 1-180-day-old calves. The complement activity in CNVR_221 homozygotes was significantly lower than that in heterozygotes and wild type animals. The frequency of the risk allele was higher in dead calves suffering from an infectious disease than in healthy cows. These results suggest that the existence of CNVR_221 in calves could be attributed to a reduction in complement activity, which in turn leads to susceptibility to infections. Thus, the risk allele could serve as a useful marker to reduce the mortality of infected Japanese Black calves.

Quantification and characterization of the 5' exonuclease activity of the lysosomal nuclease PLD3 by a novel cell-based assay

Phospholipase D3 (PLD3) and phospholipase D4 (PLD4), the most recently described lysosomal nucleases, are associated with Alzheimer’s disease, spinocerebellar ataxia, and systemic lupus erythematosus. They exhibit 5′ exonuclease activity on single-stranded DNA, hydrolyzing it at the acidic pH associated with the lysosome. However, their full cellular function is inadequately understood. To examine these enzymes, we developed a robust and automatable cell-based assay based on fluorophore- and fluorescence-quencher-coupled oligonucleotides for the quantitative determination of acidic 5′ exonuclease activity. We validated the assay under knockout and PLD-overexpression conditions and then applied it to characterize PLD3 and PLD4 biochemically. Our experiments revealed PLD3 as the principal acid 5′ exonuclease in HeLa cells, where it showed a markedly higher specific activity compared with PLD4. We further used our newly developed assay to determine the substrate specificity and inhibitory profile of PLD3 and found that proteolytic processing of PLD3 is dispensable for its hydrolytic activity. We followed the expression, proteolytic processing, and intracellular distribution of genetic PLD3 variants previously associated with Alzheimer’s disease and investigated each variant's effect on the 5′ nuclease activity of PLD3, finding that some variants lead to reduced activity, but others not. The development of a PLD3/4-specific biochemical assay will be instrumental in understanding better both nucleases and their incompletely understood roles in vitro and in vivo.

A 63-bp insertion in exon 2 of the porcine KIF21A gene is associated with arthrogryposis multiplex congenita

A recessive form of arthrogryposis multiplex congenita (AMC) was detected 20 years ago in the Swiss Large White (SLW) pig population. A diagnostic marker test enabled the identification of carrier animals, but the underlying causal mutation remains unknown. To identify the mutation underlying AMC, we collected SNP chip genotyping data for 11 affected piglets and 23 healthy pigs. Association testing using 47 829 SNPs confirmed that AMC maps to SSC5 (P = 9.4 × 10^-13 ). Subsequent autozygosity mapping revealed a common 6.06 Mb region (from 66 757 970 to 72 815 151 bp) of extended homozygosity in 11 piglets affected by AMC. Using WGS data, we detected a 63-bp insertion compatible with the recessive inheritance of AMC in the second exon of KIF21A gene encoding Kinesin Family Member 21A. The 63-bp insertion is predicted to introduce a premature stop codon in KIF21A gene (p.Val41_Phe42insTer) that truncates 1614 amino acids (~97%) from the protein. We found that this deleterious allele still segregates at a frequency of 0.1% in the SLW pig population. Carrier animals can now be detected unambiguously and excluded from breeding.

An IL17RA frameshift variant in a Holstein cattle family with psoriasis-like skin alterations and immunodeficiency

BACKGROUND

Skin lesions and dermatoses in cattle are often associated with infections due to bacteria, fungi or environmental risk factors. Dermatoses with genetic etiology have been described in cattle. Among these rare disorders, there are primary congenital dermatoses that are associated with inherited nutritional deficiencies, such as bovine hereditary zinc deficiency or zinc deficiency-like syndrome. This study presents three cases of Holstein cattle with congenital skin lesions observed on a single farm that resemble zinc deficiency-like syndrome. Close clinical and pathological examinations took place in two cases. Pedigree analysis indicated autosomal recessive inheritance and whole-genome sequencing of both affected calves was performed.

RESULTS

The two calves showed retarded growth and suffered from severe ulcerative dermatitis with hyperkeratosis, alopecia furunculosis and subcutaneous abscess formation. Blood analysis showed correspondent leukocytosis with neutrophilia whereas minerals, macro- and micronutrients were within the reference ranges. Variant calling and filtering against the 1000 Bull Genomes variant catalogue resulted in the detection of a single homozygous protein-changing variant exclusively present in both sequenced genomes. This single-nucleotide deletion in exon 3 of IL17RA on bovine chromosome 5 was predicted to have a deleterious impact on the encoded protein due to a frameshift leading to a truncated gene product. Genotyping of the affected cattle family confirmed recessive inheritance.

CONCLUSIONS

A loss-of-function mutation of the IL17RA transmembrane protein could be identified as most likely pathogenic variant for the psoriasis-like skin alterations observed in the two affected Holstein calves. In man, rare recessive diseases associated with IL17RA include immunodeficiency 51 and chronic mucocutaneous candidiasis. This supports the observed immunodeficiency of the presented cases. This study reports the first naturally occurring IL17RA-associated animal model.

A highly contiguous genome assembly of the bat hawkmoth Hyles vespertilio (Lepidoptera: Sphingidae)

BACKGROUND

Adapted to different ecological niches, moth species belonging to the Hyles genus exhibit a spectacular diversity of larval color patterns. These species diverged ∼7.5 million years ago, making this rather young genus an interesting system to study a wide range of questions including the process of speciation, ecological adaptation, and adaptive radiation.

RESULTS

Here we present a high-quality genome assembly of the bat hawkmoth Hyles vespertilio, the first reference genome of a member of the Hyles genus. We generated 51× Pacific Biosciences long reads with an average read length of 8.9 kb. Pacific Biosciences reads longer than 4 kb were assembled into contigs, resulting in a 651.4-Mb assembly consisting of 530 contigs with an N50 value of 7.5 Mb. The circular mitochondrial contig has a length of 15,303 bp. The H. vespertilio genome is very repeat-rich and exhibits a higher repeat content (50.3%) than other Bombycoidea species such as Bombyx mori (45.7%) and Manduca sexta (27.5%). We developed a comprehensive gene annotation workflow to obtain consensus gene models from different evidence including gene projections, protein homology, transcriptome data, and ab initio predictions. The resulting gene annotation is highly complete with 94.5% of BUSCO genes being completely present, which is higher than the BUSCO completeness of the B. mori (92.2%) and M. sexta (90%) annotations.

CONCLUSIONS

Our gene annotation strategy has general applicability to other genomes, and the H. vespertilio genome provides a valuable molecular resource to study a range of questions in this genus, including phylogeny, incomplete lineage sorting, speciation, and hybridization. A genome browser displaying the genome, alignments, and annotations is available at https://genome-public.pks.mpg.de/cgi-bin/hgTracks?db=HLhylVes1.

Phospholipase D and the Mitogen Phosphatidic Acid in Human Disease: Inhibitors of PLD at the Crossroads of Phospholipid Biology and Cancer

Lipids are key building blocks of biological membranes and are involved in complex signaling processes such as metabolism, proliferation, migration, and apoptosis. Extracellular signaling by growth factors, stress, and nutrients is transmitted through receptors that activate lipid-modifying enzymes such as the phospholipases, sphingosine kinase, or phosphoinositide 3-kinase, which then modify phospholipids, sphingolipids, and phosphoinositides. One such important enzyme is phospholipase D (PLD), which cleaves phosphatidylcholine to yield phosphatidic acid and choline. PLD isoforms have dual role in cells. The first involves maintaining cell membrane integrity and cell signaling, including cell proliferation, migration, cytoskeletal alterations, and invasion through the PLD product PA, and the second involves protein-protein interactions with a variety of binding partners. Increased evidence of elevated PLD expression and activity linked to many pathological conditions, including cancer, neurological and inflammatory diseases, and infection, has motivated the development of dual- and isoform-specific PLD inhibitors. Many of these inhibitors are reported to be efficacious and safe in cells and mouse disease models, suggesting the potential for PLD inhibitors as therapeutics for cancer and other diseases. Current knowledge and ongoing research of PLD signaling networks will help to evolve inhibitors with increased efficacy and safety for clinical studies.

Zinc Deficiency-Like Syndrome in Fleckvieh Calves: Clinical and Pathological Findings and Differentiation from Bovine Hereditary Zinc Deficiency

Background

Zinc deficiency‐like (ZDL) syndrome is an inherited defect of Fleckvieh calves, with striking similarity to bovine hereditary zinc deficiency (BHZD). However, the causative mutation in a phospholipase D4 encoding gene (PLD4) shows no connection to zinc metabolism.

Objectives

To describe clinical signs, laboratory variables, and pathological findings of ZDL syndrome and their utility to differentiate ZDL from BHZD and infectious diseases with similar phenotype.

Animals

Nine hospitalized calves with crusting dermatitis and confirmed mutation in PLD4 and medical records from 25 calves with crusting dermatitis or suspected zinc deficiency.

Methods

Prospective and retrospective case series.

Results

The 9 calves (age: 5–53 weeks) displayed a moderate to severe crusting dermatitis mainly on the head, ventrum, and joints. Respiratory and digestive tract inflammations were frequently observed. Zinc supplementation did not lead to remission of clinical signs in 4 calves. Laboratory variables revealed slight anemia in 8 calves, hypoalbuminemia in 6 calves, but reduced serum zinc concentrations in only 3 calves. Mucosal erosions/ulcerations were present in 7 calves and thymus atrophy or reduced thymic weights in 8 calves. Histologically, skin lesions were indistinguishable from BHZD. Retrospective analysis of medical records revealed the presence of this phenotype since 1988 and pedigree analysis revealed a common ancestor of several affected calves.

Conclusions and Clinical Importance

ZDL syndrome should be suspected in Fleckvieh calves with crusting dermatitis together with diarrhea or respiratory tract inflammations without response to oral zinc supplementation. Definite diagnosis requires molecular genetic confirmation of the PLD4 mutation.

Identification of a haplotype associated with cholesterol deficiency and increased juvenile mortality in Holstein cattle

Over the last decades, several genetic disorders have been discovered in cattle. However, the genetic background of disorders in calves is less reported. Recently, German cattle farmers reported on calves from specific matings with chronic diarrhea and retarded growth of unknown etiology. Affected calves did not respond to any medical treatment and died within the first months of life. These calves were underdeveloped in weight and showed progressive and severe emaciation despite of normal feed intake. Hallmark findings of the blood biochemical analysis were pronounced hypocholesterolemia and deficiency of fat-soluble vitamins. Results of the clinical and blood biochemical examination had striking similarities with findings reported in human hypobetalipoproteinemia. Postmortem examination revealed near-complete atrophy of the body fat reserves including the spinal canal and bone marrow. To identify the causal region, we performed a genome-wide association study with 9 affected and 21,077 control animals genotyped with the Illumina BovineSNP50 BeadChip (Illumina Inc., San Diego, CA), revealing a strong association signal on BTA 11. Subsequent autozygosity mapping identified a disease-associated haplotype encompassing 1.01 Mb. The segment of extended homozygosity contains 6 transcripts, among them the gene APOB, which is causal for cholesterol disorders in humans. However, results from multi-sample variant calling of 1 affected and 47 unaffected animals did not detect any putative causal mutation. The disease-associated haplotype has an important adverse effect on calf mortality in the homozygous state when comparing survival rates of risk matings vs. non-risk matings. Blood cholesterol values of animals are significantly associated with the carrier status indicating a codominant inheritance. The frequency of the haplotype in the current Holstein population was estimated to be 4.2%. This study describes the identification and phenotypic manifestation of a new Holstein haplotype characterized by pronounced hypocholesterolemia, chronic emaciation, growth retardation, and increased mortality in young cattle, denominated as cholesterol deficiency haplotype. Our genomic investigations and phenotypic examinations provide additional evidence for a mutation within the APOB gene causing cholesterol deficiency in Holstein cattle.

A missense mutation in TUBD1 is associated with high juvenile mortality in Braunvieh and Fleckvieh cattle

Background

Haplotypes with reduced or missing homozygosity may harbor deleterious alleles that compromise juvenile survival. A scan for homozygous haplotype deficiency revealed a short segment on bovine chromosome 19 (Braunvieh haplotype 2, BH2) that was associated with high juvenile mortality in Braunvieh cattle. However, the molecular genetic underpinnings and the pathophysiology of BH2 remain to be elucidated.

Results

The frequency of BH2 was 6.5 % in 8,446 Braunvieh animals from the national bovine genome databases. Both perinatal and juvenile mortality of BH2 homozygous calves were higher than the average in Braunvieh cattle resulting in a depletion of BH2 homozygous adult animals (P = 9.3x10⁻¹²). The analysis of whole-genome sequence data from 54 Braunvieh animals uncovered a missense mutation in TUBD1 (rs383232842, p.H210R) that was compatible with recessive inheritance of BH2. The availability of sequence data of 236 animals from diverse bovine populations revealed that the missense mutation also segregated at a low frequency (1.7 %) in the Fleckvieh breed. A validation study in 37,314 Fleckvieh animals confirmed high juvenile mortality of homozygous calves (P = 2.2x10⁻¹⁵). Our findings show that the putative disease allele is located on an ancestral haplotype that segregates in Braunvieh and Fleckvieh cattle. To unravel the pathophysiology of BH2, six homozygous animals were examined at the animal clinic. Clinical and pathological findings revealed that homozygous calves suffered from chronic airway disease possibly resulting from defective cilia in the respiratory tract.

Conclusions

A missense mutation in TUBD1 is associated with high perinatal and juvenile mortality in Braunvieh and Fleckvieh cattle. The mutation is located on a common haplotype likely originating from an ancient ancestor of Braunvieh and Fleckvieh cattle. Our findings demonstrate for the first time that deleterious alleles may segregate across closed cattle breeds without recent admixture. Homozygous calves suffer from chronic airway disease resulting in poor growth performance and high juvenile mortality. The respiratory manifestations resemble key features of diseases resulting from impaired function of airway cilia.

Electronic supplementary material

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

A frameshift mutation in GON4L is associated with proportionate dwarfism in Fleckvieh cattle

Background

Low birth weight and postnatal growth restriction are the most evident symptoms of dwarfism. Accompanying skeletal aberrations may compromise the general condition and locomotion of affected individuals. Several paternal half-sibs with a low birth weight and a small size were born in 2013 in the Fleckvieh cattle population.

Results

Affected calves were strikingly underweight at birth in spite of a normal gestation length and had craniofacial abnormalities such as elongated narrow heads and brachygnathia inferior. In spite of a normal general condition, their growth remained restricted during rearing. We genotyped 27 affected and 10,454 unaffected animals at 44,672 single nucleotide polymorphisms and performed association tests followed by homozygosity mapping, which allowed us to map the locus responsible for growth failure to a 1.85-Mb segment on bovine chromosome 3. Analysis of whole-genome re-sequencing data from one affected and 289 unaffected animals revealed a 1-bp deletion (g.15079217delC, rs723240647) in the coding region of the GON4L gene that segregated with the dwarfism-associated haplotype. We showed that the deletion induces intron retention and premature termination of translation, which can lead to a severely truncated protein that lacks domains that are likely essential to normal protein function. The widespread use of an undetected carrier bull for artificial insemination has resulted in a tenfold increase in the frequency of the deleterious allele in the female population.

Conclusions

A frameshift mutation in GON4L is associated with autosomal recessive proportionate dwarfism in Fleckvieh cattle. The mutation has segregated in the population for more than 50 years without being recognized as a genetic disorder. However, the widespread use of an undetected carrier bull for artificial insemination caused a sudden accumulation of homozygous calves with dwarfism. Our findings provide the basis for genome-based mating strategies to avoid the inadvertent mating of carrier animals and thereby prevent the birth of homozygous calves with impaired growth.

Electronic supplementary material

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

A frameshift mutation in ARMC3 is associated with a tail stump sperm defect in Swedish Red (Bos taurus) cattle

BACKGROUND

Artificial insemination is widely used in many cattle breeding programs. Semen samples of breeding bulls are collected and closely examined immediately after collection at artificial insemination centers. Only ejaculates without anomalous findings are retained for artificial insemination. Although morphological aberrations of the spermatozoa are a frequent reason for discarding ejaculates, the genetic determinants underlying poor semen quality are scarcely understood.

RESULTS

A tail stump sperm defect was observed in three bulls of the Swedish Red cattle breed. The spermatozoa of affected bulls were immotile because of severely disorganized tails indicating disturbed spermatogenesis. We genotyped three affected bulls and 18 unaffected male half-sibs at 46,035 SNPs and performed homozygosity mapping to map the fertility disorder to an 8.42 Mb interval on bovine chromosome 13. The analysis of whole-genome re-sequencing data of an affected bull and 300 unaffected animals from eleven cattle breeds other than Swedish Red revealed a 1 bp deletion (Chr13: 24,301,425 bp, ss1815612719) in the eleventh exon of the armadillo repeat containing 3-encoding gene (ARMC3) that was compatible with the supposed recessive mode of inheritance. The deletion is expected to alter the reading frame and to induce premature translation termination (p.A451fs26). The mutated protein is shortened by 401 amino acids (46 %) and lacks domains that are likely essential for normal protein function.

CONCLUSIONS

We report the phenotypic and genetic characterization of a sterilizing tail stump sperm defect in the Swedish Red cattle breed. Exploiting high-density genotypes and massive re-sequencing data enabled us to identify the most likely causal mutation for the fertility disorder in bovine ARMC3. Our results provide the basis for monitoring the mutated variant in the Swedish Red cattle population and for the early identification of infertile animals.

Microglial phospholipase D4 deficiency influences myelination during brain development

Phospholipase D4 (PLD4) is expressed in activated microglia that transiently appear in white matter during postnatal brain development. Previous knockdown experiments using cultured microglia showed PLD4 involvement in phagocytosis and proliferation. To elucidate the role of PLD4 in vivo, PLD4-deficient mice were generated and the cerebella were examined at postnatal day 5 (P5) and P7, when PLD4 expression is highest in microglia. Wild type microglia showed strong immunoreactivity for microglial marker CD68 at P5, whereas CD68 signals were weak in PLD4-deficient microglia, suggesting that loss of PLD4 affects microglial activation. At P5 and P7, immunostaining for anti-myelin basic protein (MBP) antibody indicated a mild but significant delay in myelination in PLD4-deficient cerebellum. Similar change was also observed in the corpus callosum at P7. However, this difference was not apparent at P10, suggesting that microglial PLD4-deficiency primarily influences the early myelination stage. Thus, microglia may have a transient role in myelination via a PLD4-related mechanism during development.

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.

Physiological and pathophysiological roles for phospholipase D

Individual members of the mammalian phospholipase D (PLD) superfamily undertake roles that extend from generating the second messenger signaling lipid, phosphatidic acid, through hydrolysis of the membrane phospholipid, phosphatidylcholine, to functioning as an endonuclease to generate small RNAs and facilitating membrane vesicle trafficking through seemingly nonenzymatic mechanisms. With recent advances in genome-wide association studies, RNA interference screens, next-generation sequencing approaches, and phenotypic analyses of knockout mice, roles for PLD family members are being uncovered in autoimmune, infectious neurodegenerative, and cardiovascular disease, as well as in cancer. Some of these disease settings pose opportunities for small molecule inhibitory therapeutics, which are currently in development.

Homozygous haplotype deficiency reveals deleterious mutations compromising reproductive and rearing success in cattle

Background

Cattle breeding populations are susceptible to the propagation of recessive diseases. Individual sires generate tens of thousands of progeny via artificial insemination. The frequency of deleterious alleles carried by such sires may increase considerably within few generations. Deleterious alleles manifest themselves often by missing homozygosity resulting from embryonic/fetal, perinatal or juvenile lethality of homozygotes.

Results

A scan for homozygous haplotype deficiency in 25,544 Fleckvieh cattle uncovered four haplotypes affecting reproductive and rearing success. Exploiting whole-genome resequencing data from 263 animals facilitated to pinpoint putatively causal mutations in two of these haplotypes. A mutation causing an evolutionarily unlikely substitution in SUGT1 was perfectly associated with a haplotype compromising insemination success. The mutation was not found in homozygous state in 10,363 animals (P = 1.79 × 10⁻⁵) and is thus likely to cause lethality of homozygous embryos. A frameshift mutation in SLC2A2 encoding glucose transporter 2 (GLUT2) compromises calf survival. The mutation leads to premature termination of translation and activates cryptic splice sites resulting in multiple exon variants also with premature translation termination. The affected calves exhibit stunted growth, resembling the phenotypic appearance of Fanconi-Bickel syndrome in humans (OMIM 227810), which is also caused by mutations in SLC2A2.

Conclusions

Exploiting comprehensive genotype and sequence data enabled us to reveal two deleterious alleles in SLC2A2 and SUGT1 that compromise pre- and postnatal survival in homozygous state. Our results provide the basis for genome-assisted approaches to avoiding inadvertent carrier matings and to improving reproductive and rearing success in Fleckvieh cattle.

Electronic supplementary material

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

The phospholipase D superfamily as therapeutic targets

The phospholipase D (PLD) lipid-signaling enzyme superfamily has long been studied for its roles in cell communication and a wide range of cell biological processes. With the advent of loss-of-function genetic mouse models that have revealed that PLD1 and PLD2 ablation is overtly tolerable, small-molecule PLD1/2 inhibitors that do not cause unacceptable clinical toxicity, a PLD2 polymorphism that has been linked to altered physiology, and growing delineation of processes that are subtly altered in mice lacking PLD1/2 activity, the stage is being set for assessment of PLD1/2 inhibition for therapeutic purposes. Based on findings to date, PLD1/2 inhibition may be of more utility in acute rather than chronic settings, although this generalization will depend on the specific risks and benefits in each disease setting.

In frame exon skipping in UBE3B is associated with developmental disorders and increased mortality in cattle

Background

Inherited developmental diseases can cause severe animal welfare and economic problems in dairy cattle. The use of a small number of bulls for artificial insemination (AI) carries a risk that recessive defects rapidly enrich in the population. In recent years, an increasing number of Finnish Ayrshire calves have been identified with signs of ptosis, intellectual disability, retarded growth and mortality, which constitute an inherited disorder classified as PIRM syndrome.

Results

We established a cohort of nine PIRM-affected calves and 38 unaffected half-siblings and performed a genome-wide association study (GWAS) to map the disease to a 700-kb region on bovine chromosome 17 (p = 1.55 × 10^-9). Whole genome re-sequencing of an unaffected carrier, its affected progeny and 43 other unaffected animals from another breed identified a G > A substitution mutation at the last nucleotide of exon 23 in the ubiquitin protein ligase E3B encoding gene (UBE3B). UBE3B transcript analysis revealed in-frame exon skipping in the affected animals resulting in an altered protein lacking 40 amino acids, of which 20 are located in the conserved HECT-domain, the catalytic site of the UBE3B protein. Mutation screening in 129 Ayrshire AI bulls currently used in Finland indicated a high carrier frequency (17.1%). We also found that PIRM syndrome might be connected to the recently identified AH1 haplotype, which has a frequency of 26.1% in the United States Ayrshire population.

Conclusion

We describe PIRM syndrome in cattle, which is associated with the mutated UBE3B gene. The bovine phenotype resembles human Kaufman oculocerebrofacial syndrome, which is also caused by mutations in UBE3B. PIRM syndrome might be connected with the recently identified AH1 haplotype, which is associated with reduced fertility in the US Ayrshire population. This study enables the development of a genetic test to efficiently reduce the high frequency of mutant UBE3B in Ayrshires, significantly improving animal health and reducing economic loss.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-890) contains supplementary material, which is available to authorized users.