Batten disease: biochemical and molecular characterization revealing novel <Emphasis Type="Italic">PPT1</Emphasis> and <Emphasis Type="Italic">TPP1</Emphasis> gene mutations in Indian patients

Background

Neuronal ceroid lipofuscinoses type I and type II (NCL1 and NCL2) also known as Batten disease are the commonly observed neurodegenerative lysosomal storage disorder caused by mutations in the PPT1 and TPP1 genes respectively. Till date, nearly 76 mutations in PPT1 and approximately 140 mutations, including large deletion/duplications, in TPP1 genes have been reported in the literature. The present study includes 34 unrelated Indian patients (12 females and 22 males) having epilepsy, visual impairment, cerebral atrophy, and cerebellar atrophy.

Methods

The biochemical investigation involved measuring the palmitoyl protein thioesterase 1 and tripeptidy peptidase l enzyme activity from the leukocytes. Based on the biochemical analysis all patients were screened for variations in either PPT1 gene or TPP1 gene using bidirectional Sanger sequencing. In cases where Sanger sequencing results was uninformative Multiplex Ligation-dependent Probe Amplification technique was employed. The online tools performed the protein homology modeling and orthologous conservation of the novel variants.

Results

Out of 34 patients analyzed, the biochemical assay confirmed 12 patients with NCL1 and 22 patients with NCL2. Molecular analysis of PPT1 gene in NCL1 patients revealed three known mutations (p.Val181Met, p.Asn110Ser, and p.Trp186Ter) and four novel variants (p.Glu178Asnfs*13, p.Pro238Leu, p.Cys45Arg, and p.Val236Gly). In the case of NCL2 patients, the TPP1 gene analysis identified seven known mutations and eight novel variants. Overall these 15 variants comprised seven missense variants (p.Met345Leu, p.Arg339Trp, p.Arg339Gln, p.Arg206Cys, p.Asn286Ser, p.Arg152Ser, p.Tyr459Ser), four frameshift variants (p.Ser62Argfs*19, p.Ser153Profs*19, p.Phe230Serfs*28, p.Ile484Aspfs*7), three nonsense variants (p.Phe516*, p.Arg208*, p.Tyr157*) and one intronic variant (g.2023_2024insT). No large deletion/duplication was identified in three NCL1 patients where Sanger sequencing study was normal.

Conclusion

The given study reports 34 patients with Batten disease. In addition, the study contributes four novel variants to the spectrum of PPT1 gene mutations and eight novel variants to the TPP1 gene mutation data. The novel pathogenic variant p.Pro238Leu occurred most commonly in the NCL1 cohort while the occurrence of a known pathogenic mutation p.Arg206Cys dominated in the NCL2 cohort. This study provides an insight into the molecular pathology of NCL1 and NCL2 disease for Indian origin patients.

     

Complete chloroplast genome sequence of <Emphasis Type="Italic">Capsicum</Emphasis><Emphasis Type="Italic">baccatum</Emphasis> var. <Emphasis Type="Italic">baccatum</Emphasis>

Molecular markers derived from the complete chloroplast genome can provide effective tools for species identification and phylogenetic resolution. Complete chloroplast (cp) genome sequences of Capsicum species have been reported. We herein report the complete chloroplast genome sequence of Capsicum baccatum var. baccatum, a wild Capsicum species. The total length of the chloroplast genome is 157,145 bp with 37.7 % overall GC content. One pair of inverted repeats, 25,910 bp in length, was separated by a small single-copy region (17,974 bp) and large single-copy region (87,351 bp). This region contains 86 protein-coding genes, 30 tRNA genes, 4 rRNA genes, and 11 genes contain one or two introns. Pair-wise alignments of chloroplast genome were performed for genome-wide comparison. Analysis revealed a total of 134 simple sequence repeat (SSR) motifs and 282 insertions or deletions variants in the C. baccatum var. baccatum cp genome. The types and abundances of repeat units in Capsicum species were relatively conserved, and these loci could be used in future studies to investigate and conserve the genetic diversity of the Capsicum species.     

High-throughput sequencing of <Emphasis Type="Italic">Prunus ferganensis</Emphasis> indicates that it is a geographical population of <Emphasis Type="Italic">P. persica</Emphasis>

Peach belongs to the genus Prunus, which includes Prunus persica and its relative species, P. mira, P. davidiana, P. kansuensis, and P. ferganensis. Of these, P. ferganensis have been classified as a species, subspecies, or geographical population of P. persica. To explore the genetic difference between P. ferganensis and P. persica, high-throughput sequencing was used in different peach accessions belonging to different species. First, low-depth sequencing data of peach accessions belonging to four categories revealed that similarity between P. ferganensis and P. persica was similar to that between P. persica accessions from different geographical populations. Then, to further detect the genomic variation in P. ferganensis, the P. ferganensis accession “Xinjiang Pan Tao 1” and the P. persica accession “Xia Miao 1” were sequenced with high depth, and sequence reads were assembled. The results showed that the collinearity of “Xinjiang Pan Tao 1” with the reference genome “Lovell” was higher than that of “Xia Miao 1” and “Lovell” peach. Additionally, the number of genetic variants, including single nucleotide polymorphisms (SNPs), structural variations (SVs), and the specific genes annotated from unmapped sequence in “Xia Miao 1” was higher than that in “Xinjiang Pan Tao 1” peach. The data showed that there was a close distance between “Xinjiang Pan Tao 1” (P. ferganensis) and reference genome which belong to P. persica, comparing “Xia Miao 1” (P. persica) and reference ones. The results accompany with phylogenetic tree and structure analysis confirmed that P. ferganensis should be considered as a geographic population of P. persica rather than a subspecies or a distinct species. Furthermore, gene ontology analysis was performed using the gene comprising large-effect variation to understand the phenotypic difference between two accessions. The result revealed that the pathways of gene function affected by SVs but SNPs and insertion-deletions markedly differed between the two peach accessions.     

Exome sequencing identifies nonsegregating nonsense <Emphasis Type="Italic">ATM</Emphasis> and <Emphasis Type="Italic">PALB2</Emphasis>variants in familial pancreatic cancer

We sequenced 11 germline exomes from five families with familial pancreatic cancer (FPC). One proband had a germline nonsense variant in ATM with somatic loss of the variant allele. Another proband had a nonsense variant in PALB2 with somatic loss of the variant allele. Both variants were absent in a relative with FPC. These findings question the causal mechanisms of ATM and PALB2 in these families and highlight challenges in identifying the causes of familial cancer syndromes using exome sequencing.     

Exome sequencing identifies novel and recurrent mutations in <Emphasis Type="Italic">GJA8</Emphasis> and <Emphasis Type="Italic">CRYGD</Emphasis>associated with inherited cataract

Background

Inherited cataract is a clinically important and genetically heterogeneous cause of visual impairment. Typically, it presents at an early age with or without other ocular/systemic signs and lacks clear phenotype-genotype correlation rendering both clinical classification and molecular diagnosis challenging. Here we have utilized trio-based whole exome sequencing to discover mutations in candidate genes underlying autosomal dominant cataract segregating in three nuclear families.

Results

In family A, we identified a recurrent heterozygous mutation in exon-2 of the gene encoding γD-crystallin (CRYGD; c.70C > A, p.Pro24Thr) that co-segregated with `coralliform' lens opacities. Families B and C were found to harbor different novel variants in exon-2 of the gene coding for gap-junction protein α8 (GJA8; c.20T > C, p.Leu7Pro and c.293A > C, p.His98Pro). Each novel variant co-segregated with disease and was predicted in silico to have damaging effects on protein function.

Conclusions

Exome sequencing facilitates concurrent mutation-profiling of the burgeoning list of candidate genes for inherited cataract, and the results can provide enhanced clinical diagnosis and genetic counseling for affected families.

     

Filtering genetic variants and placing informative <Emphasis Type="Italic">priors</Emphasis> based on putative biological function

High-density genetic marker data, especially sequence data, imply an immense multiple testing burden. This can be ameliorated by filtering genetic variants, exploiting or accounting for correlations between variants, jointly testing variants, and by incorporating informative priors. Priors can be based on biological knowledge or predicted variant function, or even be used to integrate gene expression or other omics data. Based on Genetic Analysis Workshop (GAW) 19 data, this article discusses diversity and usefulness of functional variant scores provided, for example, by PolyPhen2, SIFT, or RegulomeDB annotations. Incorporating functional scores into variant filters or weights and adjusting the significance level for correlations between variants yielded significant associations with blood pressure traits in a large family study of Mexican Americans (GAW19 data set). Marker rs218966 in gene PHF14 and rs9836027 in MAP4 significantly associated with hypertension; additionally, rare variants in SNUPN significantly associated with systolic blood pressure. Variant weights strongly influenced the power of kernel methods and burden tests. Apart from variant weights in test statistics, prior weights may also be used when combining test statistics or to informatively weight p values while controlling false discovery rate (FDR). Indeed, power improved when gene expression data for FDR-controlled informative weighting of association test p values of genes was used. Finally, approaches exploiting variant correlations included identity-by-descent mapping and the optimal strategy for joint testing rare and common variants, which was observed to depend on linkage disequilibrium structure.     

<Emphasis Type="Italic">Agrobacterium tumefaciens</Emphasis>-mediated transformation of <Emphasis Type="Italic">Botryosphaeria dothidea</Emphasis>

Botryosphaeria dothidea is a severe causal agent of die-back and cankers of many woody plants and causes great losses in many regions. The pathogenic mechanism of this pathogen has not been well explored due to lack of mutants and genetic information. In this study, we developed an Agrobacterium tumefaciens-mediated transformation (ATMT) protocol for B. dothidea protoplasts using vector pBHt2 containing the hph gene as a selection marker under the control of trp C promoter. Using this protocol we successfully generated the B. dothidea transformants with efficiency about 23 transformants per 10⁵ protoplasts. This is the first report of genetic transformation of B. dothidea via ATMT and this protocol provides an effective tool for B. dothidea genome manipulation, gene identification and functional analysis.     

A new source of resistance to 2-furaldehyde from <Emphasis Type="Italic">Scheffersomyces</Emphasis> (<Emphasis Type="Italic">Pichia</Emphasis>) <Emphasis Type="Italic">stipitis</Emphasis> for sustainable lignocellulose-to-biofuel conversion

Aldehyde inhibitory compounds derived from lignocellulosic biomass pretreatment have been identified as a major class of toxic chemicals that interfere with microbial growth and subsequent fermentation for advanced biofuel production. Development of robust next-generation biocatalyst is a key for a low-cost biofuel production industry. Scheffersomyces (Pichia) stipitis is a naturally occurring C-5 sugar utilization yeast; however, little is known about the genetic background underlying its potential tolerance to biomass conversion inhibitors. We investigated and identified five uncharacterized putative aryl-alcohol dehydrogenase genes (SsAADs) from this yeast as a new source of resistance against biomass fermentation inhibitor 2-furaldehyde (furfural) by gene expression, gene cloning, and direct enzyme assay analysis using partially purified proteins. All five proteins from S. stipitis showed furfural reduction using cofactor NADH. An optimum active temperature was observed at 40 °C for SsAad1p; 30 °C for SsAad3p, SsAad4p, and SsAad5p; and 20 °C for SsAad2p. SsAad2p, SsAad3p, and SsAad4p showed tolerance to a wide range of pH from 4.5 to 8, but SsAad1p and SsAad5p were sensitive to pH changes beyond 7. Genes SsAAD2, SsAAD3, and SsAAD4 displayed significantly enhanced higher levels of expression in response to the challenge of furfural. Their encoding proteins also showed higher levels of specific activity toward furfural and were suggested as core functional enzymes contributing aldehyde resistance in S. stipitis.     

Semi-dwarf barley (<Emphasis Type="Italic">Hordeum vulgare</Emphasis> L.) <Emphasis Type="Italic">brh2</Emphasis> and <Emphasis Type="Italic">ari-l</Emphasis> mutants are deficient in a U-box E3 ubiquitin ligase

Lodging is the process where crop plants fall over and lie on the ground due to strong winds and heavy precipitation. This problem reduces yield and increases the risk of fungal infections and pre-harvest germination. In order to avoid lodging, plant breeders utilize short-culm mutants, which often have a robust culm that can support the weight of a heavy spike. In barley (Hordeum vulgare L.), thousands of short-culm mutants have been isolated in breeding programs around the world. Our long-term goal is to reveal the genetic network underlying culm length, with the objective to provide an enlarged repertoire of genes and alleles suitable for future breeding of lodging resistant barley. In the present work we studied a group of allelic brh2 and ari-l mutants, which have a relatively strong semi-dwarf phenotype and are phenotypically similar to previously identified mutants deficient in brassinosteroid signalling or metabolism. The Brh2 gene is located in the centromeric region of chromosome 4H and we applied a candidate gene approach to identify the gene. Brh2 is orthologous to TUD1 in rice (Orysa sativa L.), which encodes a U-box E3 ubiquitin ligase. We identified one missense mutation, one nonsense mutation and four deletions of the complete Brh2 gene. The mutants could respond to exogenously applied brassinolide, which suggests that the apparent brassinosteroid deficient phenotype of barley brh2 and ari-l mutants is related to brassinosteroid metabolism rather than signalling.     

Adaptive mutation related to cellulose producibility in <Emphasis Type="Italic">Komagataeibacter medellinensis</Emphasis> (<Emphasis Type="Italic">Gluconacetobacter xylinus</Emphasis>) NBRC 3288

Gluconacetobacter xylinus (formerly Acetobacter xylinum and presently Komagataeibacter medellinensis) is known to produce cellulose as a stable pellicle. However, it is also well known to lose this ability very easily. We investigated the on and off mechanisms of cellulose producibility in two independent cellulose-producing strains, R1 and R2. Both these strains were isolated through a repetitive static culture of a non-cellulose-producing K. medellinensis NBRC 3288 parental strain. Two cellulose synthase operons, types I and II, of this strain are truncated by the frameshift mutation in the bcsBI gene and transposon insertion in the bcsCII gene, respectively. The draft genome sequencing of R1 and R2 strains revealed that in both strains the bcsBI gene was restored by deletion of a nucleotide in its C-rich region. This result suggests that the mutations in the bcsBI gene are responsible for the on and off mechanism of cellulose producibility. When we looked at the genomic DNA sequences of other Komagataeibacter species, several non-cellulose-producing strains were found to contain similar defects in the type I and/or type II cellulose synthase operons. Furthermore, the phylogenetic relationship among cellulose synthase genes conserved in other bacterial species was analyzed. We observed that the cellulose genes in the Komagataeibacter shared sequence similarities with the γ-proteobacterial species but not with the α-proteobacteria and that the type I and type II operons could be diverged from a same ancestor in Komagataeibacter.     

Analysis of <Emphasis Type="Italic">COQ2</Emphasis>gene in multiple system atrophy

Background

Loss of function COQ2 mutations results in primary CoQ10 deficiency. Recently, recessive mutations of the COQ2 gene have been identified in two unrelated Japanese families with multiple system atrophy (MSA). It has also been proposed that specific heterozygous variants in the COQ2 gene may confer susceptibility to sporadic MSA. To assess the frequency of COQ2 variants in patients with MSA, we sequenced the entire coding region and investigated all exonic copy number variants of the COQ2 gene in 97 pathologically-confirmed and 58 clinically-diagnosed MSA patients from the United States.

Results

We did not find any homozygous or compound heterozygous pathogenic COQ2 mutations including deletion or multiplication within our series of MSA patients. In two patients, we identified two heterozygous COQ2 variants (p.S54W and c.403?+?10G?>?T) of unknown significance, which were not observed in 360 control subjects. We also identified one heterozygous carrier of a known loss of function p.S146N substitution in a severe MSA-C pathologically-confirmed patient.

Conclusions

The COQ2 p.S146N substitution has been previously reported as a pathogenic mutation in primary CoQ10 deficiency (including infantile multisystem disorder) in a recessive manner. This variant is the third primary CoQ10 deficiency mutation observed in an MSA case (p.R387X and p.R197H). Therefore it is possible that in the heterozygous state it may increase susceptibility to MSA. Further studies, including reassessing family history in patients of primary CoQ10 deficiency for the possible occurrence of MSA, are now warranted to resolve the role of COQ2 variation in MSA.

     

Comparative analysis of diploid species of <Emphasis Type="Italic">Avena</Emphasis> L. using cytogenetic and biochemical markers: <Emphasis Type="Italic">Avena canariensis</Emphasis> Baum et Fedak and <Emphasis Type="Italic">A. longiglumis</Emphasis> Dur.

The diploid oat species containing the A genome of two types (Al and Ac) were studied by electrophoresis of grain storage proteins (avenins), chromosome C-banding, and in situ hybridization with probes pTa71 and pTa794. The karyotypes of the studied species displayed similar C-banding patterns but differed in size and morphology of several chromosomes, presumably, resulting from structural rearrangements that took place during the divergence of A genomes from a common ancestor. In situ hybridization demonstrated an identical location of the 45S and 5S rRNA gene loci in Avena canariensis and A. longiglumis similar to that in the A. strigosa genome. However, the 5S rDNA locus in A. longiglumis (5S rDNA1) was considerably decreased in the chromosome 3Al long arm. The analysis demonstrated that these oat species were similar in the avenin component composition, although individual accessions differed in the electrophoretic mobilities of certain components. A considerable similarity of A. canariensis and A. longiglumis to the Avena diploid species carrying the As genome variant was demonstrated.     

<Emphasis Type="Italic">BcMF11</Emphasis> and its homologous sequences may form a lncRNA family in <Emphasis Type="Italic">Brassica</Emphasis> diploids

BcMF11 is a long non-coding RNA that has been identified in Brassica rapa and shown to be involved in pollen development. Here, when re-cloned the gene sequence, multiple paralogous copies of BcMF11 were identified in B. rapa (A genome). Multiple paralogous copies of BcMF11 were also found in B. nigra (B genome) and Brassica oleracea (C genome), the other two primary diploids of Brassica U triangle. While in the early diverging Brassicaceae lineage including Arabidopsis thaliana, no BcMF11 homolog was found. Phylogenetic analysis showed that the BcMF11 homologous sequences cloned from A genome or C genome could be clustered into a separate branch, respectively. However, there was no distinct cluster defined for BcMF11 homologous sequences cloned from B genome. The expression of BcMF11 in B. rapa was investigated and revealed a different result in the previous study. In addition, 12 expressed sequence tags from B. napus and B. rapa showing high similarities with BcMF11 were identified in the NCBI database, which further verified that rather than the useless repeat fragments in the genome, the BcMF11 homologous genes could transcribe. It is possible that BcMF11 and its homologous sequences may form a large gene family which might be originated in the recent ancestral lineage of Brassica.     

A novel homozygous mutation in <Emphasis Type="Italic">SZT2</Emphasis> gene in Saudi family with developmental delay,macrocephaly and epilepsy

Epileptic encephalopathies are genetically heterogeneous disorders which leads to epilepsy and cause neurological disorders. Seizure threshold 2 (SZT2) gene located on chromosome 1p34.2 encodes protein mainly expressed predominantly in the parietal and frontal cortex and dorsal root ganglia in the brain. Previous studies in mice showed that mutation in this gene can confers low seizure threshold, enhance epileptogenesis and in human may leads to facial dysmorphism, intellectual disability, seizure and macrocephaly. Objective of this study was to find out novel gene or novel mutation related to the gene phenotype. We have identified a large consanguineous Saudi family segregating developmental delay, intellectual disability, epilepsy, high forehead and macrocephaly. Exome sequencing was performed in affected siblings of the family to study the novel mutation. Whole exome sequencing data analysis, confirmed by subsequent Sanger sequencing validation study. Our results showed a novel homozygous mutation (c.9368G>A) in a substitution of a conserved glycine residue into a glutamic acid in the exon 67 of SZT2 gene. The mutation was ruled out in 100 unrelated healthy controls. The missense variant has not yet been reported as pathogenic in literature or variant databases. In conclusion, the here detected homozygous SZT2 variant might be the causative mutation that further explain epilepsy and developmental delay in this Saudi family.     

Genetic Testing of Non-familial Deaf Patients for <Emphasis Type="Italic">CIB2</Emphasis> and <Emphasis Type="Italic">GJB2</Emphasis> Mutations: Phenotype and Genetic Counselling

CIB2 and GJB2 genes variants contribute significantly in familial cases of prelingual recessive hearing loss (HL). This study was aimed to determine the CIB2 and GJB2 variants and associated phenotype in 150 non-familial individuals with HL. After getting informed consent, 150 non-familial deaf patients were enrolled and blood samples were obtained for DNA extraction. Pure tone air conduction audiometry was performed. Coding exons of CIB2 and GJB2 genes were Sanger sequenced. A tetra primer ARMS assay was developed for recurrent CIB2 variant. Four bi-allelic GJB2 variants, c.71G>A p.(Trp24*), c.231G>A p.(Trp77*), c.235delC p.(Leu79Cysfs3*) and c.35delG p.(Gly11Leufs24*), were found in nine hearing impaired individuals. We also found four homozygotes and five carriers of c.380G>A p. (Arg127His) variant of controversial clinical significance. CIB2 sequencing revealed single recurrent variant c.272T>C p. (Phe91Ser) segregating with HL in ten individuals. Among our patients, c.71G>A (p.Trp24*) was the most common variant, accounted for 45% of GJB2 variants. Two known GJB2 variants, c.235delC p. (Leu79Cysfs3*) and c.310del14 p. (Lys105Argfs2*), are reported here for the first time in Pakistani population. Our data further support the benign nature of c.380G>A p. (Arg127His) variant. For CIB2, c.272T>C p. (Phe91Ser) is the second common cause of HL among our sporadic cases. Phenotypically, in our patients, individuals homozygous for GJB2 variants had profound HL, whereas CIB2 homozygotes had severe to profound prelingual HL. Our results suggest that GJB2 and CIB2 are common cause of HL in different Pakistani ethnicities.     

Development of an efficient genetic manipulation strategy for sequential gene disruption and expression of different heterologous <Emphasis Type="Italic">GFP</Emphasis> genes in <Emphasis Type="Italic">Candida tropicalis</Emphasis>

The diploid yeast Candida tropicalis, which can utilize n-alkane as a carbon and energy source, is an attractive strain for both physiological studies and practical applications. However, it presents some characteristics, such as rare codon usage, difficulty in sequential gene disruption, and inefficiency in foreign gene expression, that hamper strain improvement through genetic engineering. In this work, we present a simple and effective method for sequential gene disruption in C. tropicalis based on the use of an auxotrophic mutant host defective in orotidine monophosphate decarboxylase (URA3). The disruption cassette, which consists of a functional yeast URA3 gene flanked by a 0.3 kb gene disruption auxiliary sequence (gda) direct repeat derived from downstream or upstream of the URA3 gene and of homologous arms of the target gene, was constructed and introduced into the yeast genome by integrative transformation. Stable integrants were isolated by selection for Ura⁺ and identified by PCR and sequencing. The important feature of this construct, which makes it very attractive, is that recombination between the flanking direct gda repeats occurs at a high frequency (10^?8) during mitosis. After excision of the URA3 marker, only one copy of the gda sequence remains at the recombinant locus. Thus, the resulting ura3 strain can be used again to disrupt a second allelic gene in a similar manner. In addition to this effective sequential gene disruption method, a codon-optimized green fluorescent protein-encoding gene (GFP) was functionally expressed in C. tropicalis. Thus, we propose a simple and reliable method to improve C. tropicalis by genetic manipulation.     

A <Emphasis Type="Italic">de novo</Emphasis> germline mutation in <Emphasis Type="Italic">MYH7</Emphasis> causes a progressive dominant myopathy in pigs

Background

About 9% of the offspring of a clinically healthy Piétrain boar named ‘Campus’ showed a progressive postural tremor called Campus syndrome (CPS). Extensive backcross experiments suggested a dominant mode of inheritance, and the founder boar was believed to be a gonadal mosaic. A genome-scan mapped the disease-causing mutation to an 8 cM region of porcine chromosome 7 containing the MHY7 gene. Human distal myopathy type 1 (MPD1), a disease partially resembling CPS in pigs, has been associated with mutations in the MYH7 gene.

Results

The porcine MYH7 gene structure was predicted based on porcine reference genome sequence, porcine mRNA, and in comparison to the human ortholog. The gene structure was highly conserved with the exception of the first exon. Mutation analysis of a contiguous genomic interval of more than 22 kb spanning the complete MYH7 gene revealed an in-frame insertion within exon 30 of MYH7 (c.4320_4321insCCCGCC) which was perfectly associated with the disease phenotype and confirmed the dominant inheritance. The mutation is predicted to insert two amino acids (p.Ala1440_Ala1441insProAla) in a very highly conserved region of the myosin tail. The boar ‘Campus’ was shown to be a germline and somatic mosaic as assessed by the presence of the mutant allele in seven different organs.

Conclusion

This study illustrates the usefulness of recently established genomic resources in pigs. We have identified a spontaneous mutation in MYH7 as the causative mutation for CPS. This paper describes the first case of a disorder caused by a naturally occurring mutation in the MYH7 gene of a non-human mammalian species. Our study confirms the previous classification as a primary myopathy and provides a defined large animal model for human MPD1. We provide evidence that the CPS mutation occurred during the early development of the boar ‘Campus’. Therefore, this study provides an example of germline mosaicism with an asymptomatic founder.