In frame exon skipping in UBE3B is associated with developmental disorders and increased mortality in cattle
BMC Genomics volume 15, Article number: 890 (2014)
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.
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.
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.
Inherited developmental diseases can cause severe animal welfare and economic problems in dairy cattle breeding. Where a limited number of sires are used for artificial insemination (AI), recessive genetic defects can rapidly enrich in cattle populations. Examples include degenerative axonopathy in Tyrolean Grey cattle  and compromised growth and regulation of the inflammatory response in Belgian Blue cattle . We recently reported foetal growth retardation and stillbirth in half the progeny of an Ayrshire bull carrying a genetic deletion in the MIMT1 gene . Identification of the causative mutation and development of a gene test enable such carrier animals to be excluded from breeding.
In the past few years, an increasing number of Finnish Ayrshire calves have been identified with a combination of severe symptoms including ptosis, intellectual disability, retarded growth and mortality, a disorder classified as PIRM syndrome. Bovine PIRM resembles the human autosomal recessive neurodevelopmental disorder Kaufman oculocerebrofacial syndrome, also known as blepharophimosis-ptosis-intellectual disability syndrome (MIM 615057, MIM 244450), caused by ubiquitin protein ligase E3B (UBE3B) mutations [4–6].
UBE3B belongs to the family of ubiquitin E3 ligases involved in protein ubiquitination, a post-translational protein regulation pathway that plays a key role in several biological processes during organogenesis and neurodevelopment. Mutations of other E3 ligases are associated with a variety of human developmental diseases. Increased copy number of the HUWE1 gene (HECT, UBA and WWE domain containing 1, E3 ubiquitin protein ligase) causes cognitive impairment in males (MIM 300706) . Missense mutations in CBL (Cbl proto-oncogene, E3 ubiquitin protein ligase) cause impaired growth, developmental delay, cryptorchidism and predisposition to juvenile myelomonocytic leukemia [8, 9]. The best known example of these defects is Angelman syndrome characterised by intellectual disability, absence of speech, motor dysfunction and seizures (MIM 105830) caused by loss of function of the imprinted gene UBE3A (ubiquitin protein ligase E3A) [10, 11].
Here we report that PIRM syndrome (ptosis, intellectual disability, retarded growth and mortality) in cattle is associated with an exon skipping mutation in UBE3B and this mutation is present at high frequency in the sample of AI bulls tested. Moreover, our data suggest an association between the recently identified AH1 haplotype  and PIRM. Our findings have practical implications for cattle breeding and provide a new model for human Kaufman oculocerebrofacial syndrome.
PIRM syndrome in the Ayrshire population
Farmers and breeding counsellors have reported an increasing number of calves with developmental defects including ptosis, post-natal growth retardation and increased juvenile mortality in the Finnish Ayrshire population between 2011 and 2014. Some affected calves also suffered from feeding problems, minor structural changes of the head and muscular hypotonia (Figure 1, Additional file 1: Table S1). Many affected calves failed to thrive and died at a very young age if not euthanized before. Breeders also reported learning difficulties indicating intellectual disability. For example calves had difficulties learning how to use feeding buckets. Surviving calves required special care during the neonatal period and later showed growth retardation. Usually, farmers culled affected animals before breeding. Both sexes were equally affected. The phenotype has been defined as PIRM syndrome according to its typical features (ptosis, intellectual disability, retarded growth, mortality). The analysis of pedigree records of affected animals and their close relatives indicated an autosomal recessive mode of inheritance.
PIRM syndrome maps to a 713 kb segment on bovine chromosome 17
To identify the genomic region associated with PIRM syndrome, nine affected and 37 unaffected calves descended from one AI bull were genotyped with a bovine high-density genotyping array. After quality control, genotypes for 623,881 SNPs were phased using Beagle's hidden Markov model based algorithm. The haplotypes obtained were then used in a genome wide association study. A sliding window-based approach was used to compare haplotype frequency in cases and controls, which revealed a strong association on bovine chromosome 17 (Figure 2A). The most significant association (P = 1.55 × 10-9) resulted from four adjacent haplotype windows located between 65,659,074 bp and 65,981,740 bp. To narrow down the associated region, the genotypes of affected animals were screened for segments of homozygosity. A common 713 kb region (65,645,831 bp - 66,358,629 bp) with extended homozygosity was present in all affected animals while none of the unaffected animals showed homozygosity, suggesting a recessive pattern of inheritance. The risk haplotype encompasses 14 genes (Figure 2B-C, Additional file 2: Table S2).
A synonymous mutation in UBE3B coincides exactly with PIRM syndrome
To identify the underlying mutation responsible for PIRM, the whole genomes of an obligate carrier and one of its affected progeny were sequenced to average read depths of 10.4 and 10.1. Sequence depth analysis across the 713 kb haplotype did not reveal large structural elements in the PIRM associated region (Additional file 3: Figure S1). To help identify candidate causal mutations, we used re-sequencing data from 43 previously sequenced Fleckvieh animals . There is no evidence that the PIRM syndrome occurs in breeds other than Ayrshire. Thus, the mutation causing the PIRM syndrome should not segregate among the sequenced Fleckvieh animals. Multi-sample variant calling in the 713 kb segment of extended homozygosity yielded genotypes for 1825 polymorphic sites (1684 SNPs, 141 Indels, Additional file 4: File S1). Among these only four SNPs were compatible with recessive inheritance that is homozygous for the reference allele in 43 FV animals, heterozygous in the carrier bull and homozygous for the non-reference allele in the affected animal. The functional effects of those variants were predicted based on annotation of the UMD3.1 bovine genome assembly . One of the four compatible variants was intergenic, two were located in intronic regions of TRPV4 and UBE3B, and one variant was located in the coding region of UBE3B (rs475678587, Chr17:65,921,497 bp) (Table 1). However, two of the four variants (rs440561578, rs467377722) segregated also among 191 non-Fleckvieh animals that were sequenced in the context of the 1000 bull genomes project  and can be therefore excluded as being causative. In conclusion, only a coding and an intronic variant of the UBE3B gene segregate with the PIRM syndrome.
We annotated the bovine UBE3B gene and found that it consists of 32 exons, of which exons 1 to 6 are non-coding (Additional file 5: Table S3). Variant rs475678587 is a G > A polymorphism in the third base of codon 692 that would appear to cause a synonymous substitution, p.E692E in exon 23 of UBE3B. The rs475678587 polymorphism was validated by Sanger sequencing in the carrier bull, its nine affected and 37 unaffected descendants. In addition, nine new cases and two control animals were also genotyped by a KASP genotyping assay. All 18 affected animals were homozygous for the rs475678587 A variant, whereas the unaffected animals were either heterozygous or wild type (reference allele). The suspected carrier bull and 23 of its descendants were heterozygous.
Incidence of the rs475678587 A variant in AI bulls
To determine the frequency of the rs475678587 A variant in the Finnish Ayrshire population we analysed 129 AI bulls and found a 17.1% carrier frequency. Assuming random mating, this would generate one affected calf per 138 offspring. Of those tested, 29 AI bulls had a known haplotype status for AH1. The AH1 haplotype was perfectly associated with the rs475678587 A mutation in this cohort: All 11 bulls that were carriers of the rs475678587 A mutation carry also the AH1 and 18 wild type bulls did not have the AH1.
The rs475678587 A mutation affects splicing of UBE3B
The rs475678587 G > A polymorphism is located at the junction of exon 23 and intron 23 (Additional file 5: Table S3, Figure 3B) and could therefore affect RNA splicing. The effect of the mutation on UBE3B splicing was investigated by RT-PCR in samples of cerebral cortex, tectum, hippocampus, cerebellum, lung, liver, heart, kidney, spleen, and ovary from three affected and two unaffected animals, using two primer pairs. Primer pair 2 (Additional file 6: Table S4), which flanked exon 23, resulted in the amplification of two fragments from the affected animals and only one fragment from the unaffected animals (Figure 4A). Sequencing of the smaller RT-PCR product revealed in-frame exon 23 skipping. The other fragment of expected size corresponded to the reference sequence (the University of Maryland reference sequence UMD3.1, ) with the exception of the rs475678587 A variant that was detected in affected animals.
The UBE3B gene was expressed in all examined tissues and the altered splicing pattern was observed in all tissues of affected animals, excluding tissue specificity. RT-PCR data obtained with primer pair 2 suggested lower relative expression levels of UBE3B in affected animals, but RT-PCR results obtained with primer pair 3 (Additional file 6: Table S4) that flanked exons 27–29 showed no difference between unaffected and affected animals (Figure 4B).
In silico analysis predicts in-frame deletion of 40 amino acids
Amino acid alignment of the normal and mutated proteins showed that loss of exon 23 results in an in-frame deletion of 40 amino acids between residues 652–692, deleting 20 amino acids from the E2 binding subdomain of HECT in UBE3B (Figure 3). Three-dimensional modelling of the mutated HECT domain revealed a lack of one α-helix structure and considerable structural differences compared to the normal HECT domain (Additional file 7: Figure S2).
We report a new congenital developmental syndrome in the Ayrshire breed and indicate its association with an exon skipping mutation in the UBE3B gene. The complex disease phenotype is characterized by ptosis, intellectual disability, retarded growth and mortality, and named PIRM accordingly. Genetic analyses based on across-breed comparisons efficiently shortlisted the number of candidate causative variants and ultimately identified an in-frame exon skipping mutation in UBE3B. Similar genetic comparisons approaches recently enabled the rapid identification of disease-causing alleles in cattle [15, 17]. Exon skipping results in a partial truncation of the HECT-domain in the UBE3B protein, likely compromising its function.
Loss-of-function mutations in UBE3B cause severe neurodevelopmental disorders in humans such as developmental delay, intellectual disability and characteristic facial dysmorphisms, e.g., ptosis, blepharophimosis and telecanthus. Affected individuals also suffer from severe growth retardation, hypotonia, microcephaly, neonatal respiratory distress, difficulties with feeding and gastrointestinal tract and hypocholesterolemia (MIM 615057, MIM 244450, [4–6]). Mice engineered to lack UBE3B expression (UBE3B-/-) display severe growth retardation and significant reduction of total cholesterol and lathosterol. Additionally, increased embryonal and perinatal lethality was reported for UBE3B-/- mice . Strikingly similar pathological signs were observed in PIRM Ayrshire cattle. UBE3B mutations in heterozygous form do not cause any clinically detectable symptoms in cattle, which accords with studies in human and mouse.
The functional role of UBE3B in different biochemical pathways is still unknown. Our RT-PCR results show that, as in human [18, 19] and mouse [5, 20], bovine UBE3B is expressed in several tissues. Earlier results show that UBE3B is involved in protein degradation under stress conditions. UBE3B mRNA expression was found to be up-regulated in damaged areas of chicken inner ear after noise trauma  and the UBE3B orthologue OXI-1 in C. elegans is associated with oxidative stress-response [5, 21]. Increased sensitivity toward oxidative stress might be associated with neuronal dysfunction . In addition UBE3B is connected with total cholesterol in plasma  and importantly UBE3B lesions in human and mouse are associated with hypocholesterolemia [4, 5]. Cholesterol is essential for normal development. Decreased cholesterol synthesis might be particularly harmful to the brain, which cannot utilize circulating cholesterol because of the blood–brain barrier [5, 23].
Both wild type and mutated transcripts are expressed in PIRM-affected animals. This is most likely due to ineffective splicing. Similar in-frame exon skipping mutations have been found in other developmental disorders, such as anophthalmia/microphthalmia , cholesteryl ester storage disease  and dystrophic epidermolysis bullosa . Bioinformatic analysis predicted that the mutated UBE3B protein lacks 40 amino acids, including 20 in the N-terminal E2 subdomain of the highly conserved HECT-domain. HECT domains consist of two subdomains, the large N-terminal subdomain which contains the E2 binding site, and the small C-terminal subdomain that harbours the catalytic Cys residue required for ubiquitin transfer to the substrate. The E2 binding site mediates the interaction with the cognate E2 ubiquitin-conjugating enzyme.
Cooper et al.  recently reported a segment with homozygous haplotype deficiency on bovine chromosome 17 (AH1) in Ayrshire cattle, which encompasses the rs475678587 polymorphism. In our study, the rs475678587 A mutation was perfectly associated with AH1. Our findings suggest that in-frame exon skipping in UBE3B might be connected to the lack of AH1 homozygous animals. Missing homozygosity for AH1 could be associated with juvenile mortality, growth retardation and developmental disorders of PIRM-affected calves. Reduced fertility in AH1 carrier matings  additionally could imply embryonic losses that agree with findings of UBE3B deficiency in mice . The estimated frequency of the AH1 haplotype was 26% in the US Ayrshire herd  even higher than the cohort we examined, and thus a matter of considerable concern to cattle breeders and ranchers. Our study does however provide an avenue for further investigation.
We showed here that bovine PIRM syndrome, which resembles human Kaufman oculocerebrofacial syndrome, is associated with a UBE3B mutation. The ubiquitin/proteosome system (UPS) regulates many cellular signalling pathways such as the Notch and Hedgehog that play a key role in neurogenesis . Therefore, deregulation of UPS can affect neuronal function and lead to neurological disorders . Our data support the importance of UBE3B protein for normal development in mammals. Moreover, PIRM syndrome is the first inherited disease of the ubiquitin-dependent pathway identified in cattle. The bovine PIRM model could provide the basis for comparative studies with mouse and human.
The rs475678587 A mutation was found in 17.1% of analysed Ayrshire AI bulls and is probably associated with the prevalent AH1 haplotype and therefore can cause serious economic problems for breeders and animal welfare problems if AI carrier bulls are used for breeding.
Blood sampling and clinical studies were carried out according to standard Finnish veterinary protocols. Tissue samples were collected after slaughter or euthanasia. All animal experiments were approved by the Animal Ethics Committee of the State Provincial Office of Southern Finland (ESAVI/641/04.10.07/2014).
Samples of venous blood were collected (with EDTA) from 18 affected and 41 unaffected animals for mutation mapping and frequency analysis. Most of the clinical examinations and symptom observations were done in farms by local veterinarians, farmers and breeding counsellors. Semen samples of 129 randomly selected bulls that are currently or have recently been used for AI in Finland were collected to gain an estimate of the carrier frequency in AI bulls. Tissue samples were collected from the cerebral cortex, tectum, hippocampus, cerebellum, lung, liver, heart, kidney, spleen, and ovary of three PIRM affected and two control animals for RNA expression analysis.
High- density genotyping, quality control and haplotype inference
Genomic DNA of forty-seven half sibs (9 affected, 38 unaffected) were isolated using a semi-automated Chemagen extraction robot (Chemagen Biopolymer-Technologie AG) and genotyped with the Illumina BovineHD Bead chip comprising 777,962 SNPs. Genotypes were called using Illumina BeadStudio data analysis software and default parameters. Quality control was carried out using the PLINK v1.07 whole genome association analysis toolset . Chromosomal positions of SNPs were determined on the basis of the UMD3.1 bovine genome assembly . 1224 Y-chromosomal, 343 mitochondrial and 1735 SNPs with unknown chromosomal position were excluded from further analysis. The genotypes of one unaffected animal were omitted because genotyping failed in more than 10% of SNPs. We further excluded 7235 SNPs because genotyping failed in more than 10% of individuals, and 149,129 SNPs with minor allele frequency less than 0.5%. The final dataset contained 46 animals (9 affected, 37 unaffected) and 623,881 SNPs with an average per individual call-rate of 99.38%. Beagle genetic analysis software (version 3.2.1)  was used with its default settings for imputation of sporadically missing genotypes and for haplotype inference.
Haplotype-based association study
A sliding window of 80 adjacent SNPs (corresponding to an average haplotype length of ~0.33 Mb) was shifted along the entire genome in steps of 15 SNPs. Within each window, haplotypes with frequency >5% were tested for association with PIRM syndrome using Fisher's exact test of allelic association.
Generation of sequence data
Genomic DNA was prepared from semen of the supposed carrier bull and blood from one of his affected progeny by standard protocols using proteinase K digestion and phenol-chloroform extraction. Paired-end libraries were prepared using the paired-end TruSeq DNA sample preparation kit (Illumina inc., San Diego, CA, USA) and sequenced using the HiSeq 2000 system (Illumina, San Diego, CA, USA). The read length was 101 bp. Resulting reads were processed during the sequencing step with the Illumina BaseCaller. Reads were aligned to the University of Maryland reference sequence (UMD3.1)  using the Burrows Wheeler Aligner (version 0.6.1-r104)  with default parameters. Individual files in SAM (Sequence Alignment/Map) format were converted into BAM (Binary Alignment/Map) format using SAMtools (version 0.1.18) . Duplicate reads were identified and marked with the MarkDuplicates command of Picard. Net fold coverage was 10.4 for the bull and 10.1 for its affected descendant.
Variant calling and imputation
Polymorphic sites in the region of interest (Chr17:60,000,000 bp – 70,000,000 bp) including short insertions and deletions were called for the two animals sequenced together with 43 previously sequenced Fleckvieh (FV) key ancestors  using the multi-sample approach implemented in mpileup of SAMtools along with BCFtools. Read duplicates (see above) and positions with coverage of more than 720 reads (corresponding to 2 × N samples x average coverage) were not considered in variant calling. Beagle (version 3.2.1) phasing and imputation (see above) was used to improve the primary genotype calling by SAMtools.
Identification of a candidate causal variant
Multi-sample variant calling yielded genotypes for 54,240 variants (49,245 SNPs and 4995 Indels) within the 10 Mb region of bovine chromosome 17. Among those, 1825 variants were located in the 713 kb segment (65,645,831 bp - 66,358,629 bp) of extended homozygosity. To identify mutations compatible with the presumed recessive mode of inheritance, these sites were filtered for variants that met three conditions: (i) the affected animal was homozygous for the reference allele, (ii) the presumed heterozygous carrier bull was heterozygous and (iii) all 43 sequenced FV animals were homozygous for the reference allele, as there are no indications that a similar phenotype is segregating in the FV population. Functional effects of the candidate causal variants were predicted based on gene annotation of the UMD3.1 bovine genome assembly .
Manual re-annotation of bovine UBE3B gene
The genomic structure of UBE3B was predicted based on the University of Maryland UMD3.1 bovine genome sequence assembly  and the Dana-Farber Cancer Institute bovine gene index release 12.0  using GENOMETHREADER software tool . The GENOMETHREADER output was viewed and edited using the Apollo sequence annotation editor .
Validation of rs475678587 polymorphism
PCR primer pair 1 was designed with Primer 3  for exon 23 of bovine UBE3B to scrutinise the rs475678587 polymorphism by classical Sanger sequencing in the carrier bull and its nine affected and 37 unaffected descendants. Genomic PCR products were sequenced using a 3730 × l DNA Analyser (Applied Biosystems) and data analysed with the Variant Reporter v1.0 program (Applied Biosystems). All primer pairs are listed in Additional file 6: Table S4.
Analysis of p.E692E-polymorphism in a larger cohort
DNA from semen samples of 129 AI bulls was extracted with a Chemagen extraction robot (Chemagen Biopolymer-Technologie AG). In brief, 200 μl diluted and frozen semen was washed twice with 1000 μl PBS (centrifugation for 5 min at 10,000 g) and resuspended in 500 μl lysis buffer (Chemagic DNA Blood Kit special, article No. CMG-703-1) containing 2 μl proteinase K (20 mg/ml) and 20 μl DTT (1 M). Samples were incubated overnight at 55°C and extraction continued according to the manufacturer’s instructions with 1 ml isolation buffer and 150 μl elution volume. DNA from blood samples of nine new affected calves and two control animals for the expression analysis was extracted as described earlier. Polymorphism frequency analysis was studied using KASP (Kompetitive Allele Specific PCR) reagents (LGC) and a 7500 Fast Real-Time PCR instrument (Applied Biosystems) according to the manufacturer's instructions. For quality control, every run included two samples of each rs475678587 polymorphism group determined earlier by Sanger-sequencing. From the tested 129 AI bulls 29 had genotype tested AH1 record provided individually in the Canadian Dairy Network website .
Total RNA of tissue samples was isolated with RNeasy Mini Kit (Qiagen) and converted to cDNA with High Capacity RNA-to-cDNA Kit (Applied Biosystems) according to the manufacturer's instructions. UBE3B gene expression was studied by RT-PCR analysis with two primer pairs designed with Primer 3 . PCR products of cerebral cortex and liver obtained with primer pair 2 were extracted from gel with GenElute™ Gel Extraction Kit (Sigma-Aldrich) and the DNA sequence determined. The identity of PCR products obtained with primer pair 3 was confirmed by sequence analysis from cerebral cortex and liver samples.
The effect of the p.E692E-polymorphism on mRNA splicing was predicted using the web based tool ESEfinder 3.0 [34, 35]. UBE3B protein alignment was carried out using the ClusterW2 tool  and the effect of the absence of exon 23 on protein structure was investigated using the protein homology recognition engine V2.0 - PHYRE2 .
Availability of supporting data
The sequencing data of 43 Fleckvieh animals are publically available in the Sequence Read Archive of NCBI (http://www.ncbi.nlm.nih.gov/sra) under accession numbers SRX527690-SRX527732. Four compatible variants are submitted to the The Single Nucleotide Polymorphism database of NCBI (http://www.ncbi.nlm.nih.gov/snp) under accession numbers rs440561578, rs467377722, rs463975690 and rs475678587.
Drogemuller C, Reichart U, Seuberlich T, Oevermann A, Baumgartner M, Kuhni Boghenbor K, Stoffel MH, Syring C, Meylan M, Muller S, Muller M, Gredler B, Solkner J, Leeb T: An unusual splice defect in the mitofusin 2 gene (MFN2) is associated with degenerative axonopathy in Tyrolean Grey cattle. PLoS One. 2011, 6 (4): e18931-10.1371/journal.pone.0018931.
Sartelet A, Druet T, Michaux C, Fasquelle C, Geron S, Tamma N, Zhang Z, Coppieters W, Georges M, Charlier C: A splice site variant in the bovine RNF11 gene compromises growth and regulation of the inflammatory response. PLoS Genet. 2012, 8 (3): e1002581-10.1371/journal.pgen.1002581.
Flisikowski K, Venhoranta H, Nowacka-Woszuk J, McKay SD, Flyckt A, Taponen J, Schnabel R, Schwarzenbacher H, Szczerbal I, Lohi H, Fries R, Taylor JF, Switonski M, Andersson M: A novel mutation in the maternally imprinted PEG3 domain results in a loss of MIMT1 expression and causes abortions and stillbirths in cattle (Bos taurus). PLoS One. 2010, 5 (11): e15116-10.1371/journal.pone.0015116.
Flex E: Loss of function of the E3 ubiquitin-protein ligase UBE3B causes Kaufman oculocerebrofacial syndrome. J Med Genet. 2013, 17: 2013-
Basel-Vanagaite L, Dallapiccola B, Ramirez-Solis R, Segref A, Thiele H, Edwards A, Arends MJ, Miro X, White JK, Desir J, Abramowicz M, Dentici ML, Lepri F, Hofmann K, Har-Zahav A, Ryder E, Karp NA, Estabel J, Gerdin AK, Podrini C, Ingham NJ, Altmuller J, Nurnberg G, Frommolt P, Abdelhak S, Pasmanik-Chor M, Konen O, Kelley RI, Shohat M, Nurnberg P, et al: Deficiency for the ubiquitin ligase UBE3B in a blepharophimosis-ptosis-intellectual-disability syndrome. Am J Hum Genet. 2012, 91 (6): 998-1010. 10.1016/j.ajhg.2012.10.011.
Basel-Vanagaite L, Yilmaz R, Tang S, Reuter MS, Rahner N, Grange DK, Mortenson M, Koty P, Feenstra H, Farwell Gonzalez KD, Sticht H, Boddaert N, Desir J, Anyane-Yeboa K, Zweier C, Reis A, Kubisch C, Jewett T, Zeng W, Borck G: Expanding the clinical and mutational spectrum of Kaufman oculocerebrofacial syndrome with biallelic UBE3B mutations. Hum Genet. 2014, 133 (7): 939-949. 10.1007/s00439-014-1436-2.
Froyen G, Belet S, Martinez F, Santos-Reboucas CB, Declercq M, Verbeeck J, Donckers L, Berland S, Mayo S, Rosello M, Pimentel MM, Fintelman-Rodrigues N, Hovland R, Rodrigues dos Santos S, Raymond FL, Bose T, Corbett MA, Sheffield L, van Ravenswaaij-Arts CM, Dijkhuizen T, Coutton C, Satre V, Siu V, Marynen P: Copy-number gains of HUWE1 due to replication- and recombination-based rearrangements. Am J Hum Genet. 2012, 91 (2): 252-264. 10.1016/j.ajhg.2012.06.010.
Niemeyer CM, Kang MW, Shin DH, Furlan I, Erlacher M, Bunin NJ, Bunda S, Finklestein JZ, Sakamoto KM, Gorr TA, Mehta P, Schmid I, Kropshofer G, Corbacioglu S, Lang PJ, Klein C, Schlegel PG, Heinzmann A, Schneider M, Stary J, van den Heuvel-Eibrink MM, Hasle H, Locatelli F, Sakai D, Archambeault S, Chen L, Russell RC, Sybingco SS, Ohh M, Braun BS, et al: Germline CBL mutations cause developmental abnormalities and predispose to juvenile myelomonocytic leukemia. Nat Genet. 2010, 42 (9): 794-800. 10.1038/ng.641.
Martinelli S, De Luca A, Stellacci E, Rossi C, Checquolo S, Lepri F, Caputo V, Silvano M, Buscherini F, Consoli F, Ferrara G, Digilio MC, Cavaliere ML, van Hagen JM, Zampino G, van der Burgt I, Ferrero GB, Mazzanti L, Screpanti I, Yntema HG, Nillesen WM, Savarirayan R, Zenker M, Dallapiccola B, Gelb BD, Tartaglia M: Heterozygous germline mutations in the CBL tumor-suppressor gene cause a Noonan syndrome-like phenotype. Am J Hum Genet. 2010, 87 (2): 250-257. 10.1016/j.ajhg.2010.06.015.
Fang P, Lev-Lehman E, Tsai TF, Matsuura T, Benton CS, Sutcliffe JS, Christian SL, Kubota T, Halley DJ, Meijers-Heijboer H, Langlois S, Graham JM, Beuten J, Willems PJ, Ledbetter DH, Beaudet AL: The spectrum of mutations in UBE3A causing Angelman syndrome. Hum Mol Genet. 1999, 8 (1): 129-135. 10.1093/hmg/8.1.129.
Matsuura T, Sutcliffe JS, Fang P, Galjaard RJ, Jiang YH, Benton CS, Rommens JM, Beaudet AL: De novo truncating mutations in E6-AP ubiquitin-protein ligase gene (UBE3A) in Angelman syndrome. Nat Genet. 1997, 15 (1): 74-77. 10.1038/ng0197-74.
Cooper TA, Wiggans GR, Null DJ, Hutchison JL, Cole JB: Genomic evaluation, breed identification, and discovery of a haplotype affecting fertility for Ayrshire dairy cattle. J Dairy Sci. 2014, 97 (6): 3878-3882. 10.3168/jds.2013-7427.
Jansen S, Aigner B, Pausch H, Wysocki M, Eck S, Benet-Pages A, Graf E, Wieland T, Strom TM, Meitinger T, Fries R: Assessment of the genomic variation in a cattle population by re-sequencing of key animals at low to medium coverage. BMC Genomics. 2013, 14: 446-2164. 10.1186/1471-2164-14-446. 14-446
Florea L, Souvorov A, Kalbfleisch TS, Salzberg SL: Genome assembly has a major impact on gene content: a comparison of annotation in two Bos taurus assemblies. PLoS One. 2011, 6 (6): e21400-10.1371/journal.pone.0021400.
Daetwyler HD, Capitan A, Pausch H, Stothard P, van Binsbergen R, Brondum RF, Liao X, Djari A, Rodriguez SC, Grohs C, Esquerre D, Bouchez O, Rossignol MN, Klopp C, Rocha D, Fritz S, Eggen A, Bowman PJ, Coote D, Chamberlain AJ, Anderson C, VanTassell CP, Hulsegge I, Goddard ME, Guldbrandtsen B, Lund MS, Veerkamp RF, Boichard DA, Fries R, Hayes BJ: Whole-genome sequencing of 234 bulls facilitates mapping of monogenic and complex traits in cattle. Nat Genet. 2014, 46 (8): 858-865. 10.1038/ng.3034.
Zimin AV, Delcher AL, Florea L, Kelley DR, Schatz MC, Puiu D, Hanrahan F, Pertea G, Van Tassell CP, Sonstegard TS, Marcais G, Roberts M, Subramanian P, Yorke JA, Salzberg SL: A whole-genome assembly of the domestic cow, Bos taurus. Genome Biol. 2009, 10 (4): R42-10.1186/gb-2009-10-4-r42.
Jung S, Pausch H, Langenmayer MC, Schwarzenbacher H, Majzoub-Altweck M, Gollnick NS, Fries R: A nonsense mutation in PLD4 is associated with a zinc deficiency-like syndrome in Fleckvieh cattle. BMC Genomics. 2014, 15 (1): 623-10.1186/1471-2164-15-623. -2164-15-623
Gong TW, Huang L, Warner SJ, Lomax MI: Characterization of the human UBE3B gene: structure, expression, evolution, and alternative splicing. Genomics. 2003, 82 (2): 143-152. 10.1016/S0888-7543(03)00111-3.
Lomax MI, Huang L, Cho Y, Gong TL, Altschuler RA: Differential display and gene arrays to examine auditory plasticity. Hear Res. 2000, 147 (1–2): 293-302.
Mileski A: Expression Pattern of Ubiquitin Ligase E3 (UBE3B) in Mouse. 35th Annual Meeting and Exhibition of the American Association for Dental Research, 30th Annual Meeting of the Canadian Association for Dental Research and 83rd Annual Meeting and Exhibition of American Dental Education Association. 2006
Yanase S, Ishi N: Cloning of the oxidative stress-responsive genes in Caenorhabditis elegans. J Radiat Res. 1999, 40 (1): 39-47. 10.1269/jrr.40.39.
Asselbergs FW, Guo Y, van Iperen EP, Sivapalaratnam S, Tragante V, Lanktree MB, Lange LA, Almoguera B, Appelman YE, Barnard J, Baumert J, Beitelshees AL, Bhangale TR, Chen YD, Gaunt TR, Gong Y, Hopewell JC, Johnson T, Kleber ME, Langaee TY, Li M, Li YR, Liu K, McDonough CW, Meijs MF, Middelberg RP, Musunuru K, Nelson CP, O'Connell JR, Padmanabhan S, et al: Large-scale gene-centric meta-analysis across 32 studies identifies multiple lipid loci. Am J Hum Genet. 2012, 91 (5): 823-838. 10.1016/j.ajhg.2012.08.032.
Porter FD, Herman GE: Malformation syndromes caused by disorders of cholesterol synthesis. J Lipid Res. 2011, 52 (1): 6-34. 10.1194/jlr.R009548.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81 (3): 559-575. 10.1086/519795.
Browning BL, Browning SR: A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet. 2009, 84 (2): 210-223. 10.1016/j.ajhg.2009.01.005.
Li H, Durbin R: Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010, 26 (5): 589-595. 10.1093/bioinformatics/btp698.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genome Project Data Processing Subgroup: The sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.
Quackenbush J, Cho J, Lee D, Liang F, Holt I, Karamycheva S, Parvizi B, Pertea G, Sultana R, White J: The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species. Nucleic Acids Res. 2001, 29 (1): 159-164. 10.1093/nar/29.1.159.
Gremme G: Engineering a software tool for gene structure prediction in higher organisms. Infor Softw Technol. 2005, 47 (15): 965-978. 10.1016/j.infsof.2005.09.005.
Lewis SE, Searle SM, Harris N, Gibson M, Lyer V, Richter J, Wiel C, Bayraktaroglir L, Birney E, Crosby MA, Kaminker JS, Matthews BB, Prochnik SE, Smithy CD, Tupy JL, Rubin GM, Misra S, Mungall CJ, Clamp ME: Apollo: a sequence annotation editor. Genome Biol. 2002, 3 (12): RESEARCH0082-
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.
Canadian Dairy Network. http://www.cdn.ca/home.php,
Smith PJ, Zhang C, Wang J, Chew SL, Zhang MQ, Krainer AR: An increased specificity score matrix for the prediction of SF2/ASF-specific exonic splicing enhancers. Hum Mol Genet. 2006, 15 (16): 2490-2508. 10.1093/hmg/ddl171.
Cartegni L, Wang J, Zhu Z, Zhang MQ, Krainer AR: ESEfinder: a web resource to identify exonic splicing enhancers. Nucleic Acids Res. 2003, 31 (13): 3568-3571. 10.1093/nar/gkg616.
Kelley LA, Sternberg MJ: Protein structure prediction on the Web: a case study using the Phyre server. Nat Protoc. 2009, 4 (3): 363-371. 10.1038/nprot.2009.2.
The authors would like to thank Sari Alhainen whose input in sample collection was indispensable. Thanks to veterinarians Mari Niemi, Kalle Hakala, Leena Pohjola and Sami Savolainen for collection of blood samples and professor Antti Iivanainen who also participated in sample collection. We would like to express our gratitude to Semex Finland Oy and VikingGenetics for providing the AI bull semen samples. Thanks to Minna Toivonen and Seppo Niskanen for pedigree information of the bulls. The efforts and cooperation of all farmers and the assistance of Sini Karjalainen and Lotta Yli-Kohtamäki in DNA extraction is highly appreciated. This work was supported by the Finnish Veterinary Foundation.
The authors declare that they have no competing interests.
HV performed molecular studies and participated in the study design, sample collection and extraction procedures and prepared the manuscript. HP analysed the SNP genotyping and NGS data, participated in the study design and in drafting the manuscript. KF performed the in silico analysis, participated in the study design and drafting of the manuscript. CW performed NGS. JT and HR participated in the sample collection and preparation of the manuscript. AK and AS helped write the manuscript. RF carried out the re-annotation. HL participated supervision, study design, data analysis and preparation of the manuscript. MA conceived the study, participated in the study design, supervision, sample collection, data analysis and preparation of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: Table S2: Gene content within the segment of extended homozygosity on bovine chromosome 17. The gene content was assessed based on the UMD3.1-assembly of the bovine genome sequence. A total of 14 transcripts/genes were identified within the segment of extended homozygosity. (DOCX 21 KB)
Additional file 3: Figure S1: Average read depth in the PIRM-associated region. Each dot represents the average read depth of 5 adjacent variants displayed as deviation from the average sequence coverage for the affected animal (A), the supposed carrier of the mutation (B) and an unaffected Fleckvieh animal (C). The grey shaded area represents the segment of extended homozygosity including UBE3B.(TIFF 2 MB)
Additional file 4: File S1: Genotypes for 1825 variants located in the disease-associated region in 43 Fleckvieh (FV1-FV43) and 2 Ayrshire (AYR_case, AYR_heterozygous) animals. Details of the variant calling pipeline are described in the material & methods section. (ZIP 478 KB)
Additional file 7: Figure S2: Three dimensional modelling of the UBE3B HECT domain. The photos show a PHYRE2 analysis of the normal UBE3B protein (A) and lacking 40 amino-acid of which 20 belong to E2 subdomain of HECT (B). E2 subdomain is shown in red. (TIFF 907 KB)
About this article
Cite this article
Venhoranta, H., Pausch, H., Flisikowski, K. et al. In frame exon skipping in UBE3B is associated with developmental disorders and increased mortality in cattle. BMC Genomics 15, 890 (2014). https://doi.org/10.1186/1471-2164-15-890