Introduction

Mutations in mitochondrial DNA (mtDNA) or nuclear genes encoding vital mitochondrial proteins result in a clinically heterogeneous group of mitochondrial disorders, which occur at an incidence of 1 in 5000 neonates.1 Importantly, the genetic basis of mitochondrial disorders is still unexplained in many patients, despite encouraging advances in genetic research.2

Although the existence of large deletions in mtDNA is well known, deletions affecting the nuclear genes are not commonly described in patients with mitochondrial disorders. A homozygous deletion spanning three genes including NDUFAF2 was found in a patient with fatal multisystem disorder.3 Leary et al.4 reported heterozygous de novo deletion of the SCO2 allele in combination with point mutation in a patient with neonatal hypertrophic cardiomyopathy (HCMP). Homozygous or heterozygous deletions of several exons of other OXPHOS-related disease genes were identified by application of targeted array CGH in several patients with mitochondrial diseases.5, 6

Structural variants of human DNA at least ≥1 kb in length are called copy number variations (CNVs). According to the latest studies, CNVs are significant genetic factors that may influence the apparent phenotype of every individual.7 CNVs can expose sequence variants with an unclear effect in genes of unknown function, change epigenetic factors or alter genomic neighbourhood.8, 9, 10 These facts together with both the nuclear and mitochondrial genetic backgrounds8, 11 may further contribute to the complex clinical phenotype of mitochondrial disorders.

In this study, we report the findings of two large heterozygous CNVs on 22q13.33 in combination with TYMP or SCO2 point mutations in two patients with mitochondrial disorders. Moreover, we discuss the phenotype–genotype correlation of our patients considering recent data regarding the heterogeneity of the human genome.

Patients and methods

Patient 1

The boy was born at term after an uncomplicated delivery as a second child of non-consanguineous parents of Caucasian origin. The birth weight was 3380 g and the length was 50 cm. The postnatal adaptation was uneventful, and until the age of 7.5 months, he developed normally. Thereafter, myoclonic movements on extremities were noted, and severe hypotonic syndrome developed. Failure to thrive, bilateral ptosis, divergent strabismus and ataxia were observed from the age of 9 months, accompanied by the regression of his psychomotor development. At the age of 14 months, the boy was microcephalic, with a head circumference of 44.5 cm (<3rd percentile), weight 7.9 kg (<3rd percentile) and length 79 cm (50th percentile), and his psychomotor development was on the level of a 6–9-month-old infant. Brain MRI revealed delayed myelination. EMG and echocardiography were normal. Neither hepatopathy nor elevation of creatine kinase was noted. Metabolic investigation revealed mild hyperlactaciduria (81 mmol/l creatinine, controls <30) accompanied by increased urinary excretion of the Krebs-cycle intermediates (fumarate 21 mg/g creatinine, controls <15; aconitate 150 mg/g creatinine, controls <100; 2-oxoglutarate 216 mg/g creatinine, controls <200). Lactate concentration was elevated in CSF (3 mmol/l, controls <2). The patient died due to further deterioration of the disease at the age of 20 months.

Patient 2

The boy was born at term after an uncomplicated delivery as a first child of non-consanguineous parents of Caucasian origin. The birth weight was 3550 g and the length was 52 cm. The postnatal adaptation was uneventful. Severe failure to thrive developed at the age of 4–5 years. At the age of 10 years, his weight was only 16 kg. He had frequent attacks of abdominal pain and alternation of diarrhoea and constipation due to gastrointestinal dysmotility, leading to extreme nausea and intermittent vomiting. X-ray investigation of the gastrointestinal tract found jejunal diverticulosis. Physical endurance was diminished: he was not able to climb stairs and his ability to ride a bicycle was very limited. At the age of 12 years, brain MRI showed white matter abnormalities compatible with the diagnosis of leukodystrophy, and an EMG study revealed polyneuropathy. At the age of 13 years, his weight was 17.5 kg (−4.67 SD), height 135 cm (−3.78 SD) and head circumference 52 cm (3rd percentile). Neurological investigation documented severe hypotonic syndrome, general hyporeflexia and muscle hypotrophy. Psychological investigation showed normal mental capabilities. Retinitis pigmentosa developed at the age of 16 years. At present, the 16-year-old boy is on total parenteral nutrition, and his weight is 24 kg. The severe gastrointestinal dysmotility precludes any enteral feeding. He is wheelchair-bound. He attends school on an intermittent basis. There was neither hepatopathy nor elevation of creatine kinase. During acute gastroenteritis, protein-losing enteropathy developed, leading to hypoalbuminemia (19 g/l, controls 35–53) with oedema and ascites. Metabolic investigation revealed mild hyperlactaciduria (81 mmol/l creatinine, controls <30) accompanied by increased urinary excretion of the Krebs-cycle intermediates (fumarate 18 mg/g creatinine, controls <15; aconitate 140 mg/g creatinine, controls <100) and 3-methylglutaconic aciduria (45 mg/g creatinine, controls <15). The lactate concentration in blood was not elevated. We found increased levels of uracil (93 mmol/mol creatinine, controls <15), thymine (15.4 mmol/mol creatinine, controls <15), deoxyuridine (39 mmol/mol creatinine, controls not detectable) and thymidine (12.6 mmol/mol creatinine, controls not detectable).

Molecular genetic analysis

Genomic DNA was extracted from peripheral blood lymphocytes using standard procedures. The nuclear DNA of these patients was analysed using a Genome-Wide Human SNP 6.0 microarray chip (Affymetrix, Santa Clara, CA, USA). The presence of just one copy of the SCO2 and TYMP genes in both affected patients was verified by a real-time PCR CNV assay (Hs00093549_cn, Hs00001601_cn, Hs00137275_cn, Hs00574610_cn; Applied Biosystems, Foster City, CA, USA). Sequencing of coding exons and their flanking intronic regions of SCO2 (NG_016235.1) and TYMP (NG_011860.1) PCR products was performed in both directions on an ABI PRISM 3100/3100-Avant Genetic Analyser (Applied Biosystems). The identified missense mutations were verified by PCR-RFLP analysis, and their impact was evaluated by in silico analysis methods using the PolyPhen-2 (http://genetics.bwh.harvard.edu/pph2), MutPred (http://mutpred.mutdb.org) and SNPs&GO (http://snps-and-go.biocomp.unibo.it/snps-and-go) databases.

Other methods used are described in Supplementary file 1 Methods.

Ethics

This study was approved by the Ethics Committee of the General University Hospital in Prague. All blood and tissue samples were analysed with the informed consent of the parents of the patients.

Results

In patient 1, a paternally inherited 87-kb deletion on Chr. 22: g.[(49275958_49362964)del (NCBI Build 36.1)] spanning eight genes (LMF2, NCAPH2, SCO2, TYMP, ODF3B, KLHDC7B, c22orf41 and CPT1B) was identified by microarrays (see Supplementary Figure 1) and confirmed by a real-time PCR CNV assay (see Supplementary Figure 2). This deletion corresponds to a known CNV-variation_4139, which occurs with a frequency of 0.4% in control samples.12 In combination with this deletion, a novel point mutation c.667G>A (p.Asp223Asn) in the SCO2 gene was found. The alignment of SCO2 proteins from multiple species showed that the affected Asp223 is completely conserved. The novel missense SCO2 mutation was not present in 100 Czech healthy control samples. No pathogenic mutation was found in the TYMP gene. In muscle and heart mitochondria, a markedly decreased amount of SCO2 protein was found by immunoblot analysis (Figure 1). Lower levels of copper content (Table 1) were detected in autoptic liver (10% of control value), brain, heart and muscle (60–70% of control value), further suggesting cytochrome c oxidase (COX) deficiency due to SCO2 dysfunction13 in this patient. In addition, the pathogenicity of this mutation was predicted by in silico analysis with PolyPhen-2, MutPred and SNPs&GO.

Figure 1
figure 1

Steady-state levels of SCO2 protein in muscle (a) and heart (b) mitochondria of patient 1. A total of 10 μg of mitochondrial protein from muscle and heart mitochondria was separated using SDS-PAGE (12% polyacrylamide), electroblotted onto PVDF membranes and detected with polyclonal antiserum raised against human SCO2 or with monoclonal antibody raised against the COX1, COX2 subunits and porin. Two aliquots of control mitochondria corresponding to indicated dilutions of control samples were loaded on the same gels.

Table 1 Total copper content in autoptic tissues of patient 1 with SCO2 deficiency (expressed in μg in 1 g of dry tissue)

In patient 2, a maternally inherited 175-kb deletion on Chr. 22: g.[(49275958_49451008)del (NCBI Build 36.1)] spanning 12 genes (LMF2, NCAPH2, SCO2, TYMP, ODF3B, KLHDC7B, c22orf41, CPT1B, CHKB, LOC100144603, MAPKIP2 and ARSA) along with the paternally inherited point mutation c.261G>T (p.Glu87Asp) in the TYMP gene were identified (see Supplementary Figure 1). Although the nucleotide substitution c.261G>T has not been described previously, the resulting amino-acid replacement p.Glu87Asp has already been characterised.14 The 175-kb deletion spans 12 genes and corresponds to known CNV-variation_5192, which occurs with almost 7% frequency in control samples.12 In this patient, the mitochondrial neurogastrointestinal encephalopathy (MNGIE) diagnosis was biochemically confirmed by the diminished thymidine phosphorylase activity in patient lymphocytes. The levels of thymidine in plasma and urine were increased, and the activity of COX and its ratio to citrate synthase (CS) activity in lymphocytes were decreased (Table 2). No pathogenic mutation was found in the SCO2 gene.

Table 2 Biochemical markers assessed in patient 2 with MNGIE

To evaluate the occurrence of CNVs spanning the SCO2 gene in the Czech population, a real-time PCR CNV assay was carried out in 50 control samples. However, no CNV encompassing SCO2 was found in the 50 tested Czech control samples (see Supplementary Figure 3).

Discussion

We described two large CNVs on chromosome 22q13.33 that span up to 12 genes found in two patients with mitochondrial disorders. According to the Database of Genomic Variants, only three documented CNVs of 87–175 kb on chromosome 22q13.33 and overlapping the TYMP and SCO2 genes, including CNV-variation_4139 and CNV-variation_5192 found in the presented patients, have been identified so far with total frequency of 2.2% in control samples.12 Only 2 out of the 12 genes, apart from SCO2 and TYMP, affected by the deletions are associated with human disease (Supplementary Table 1). Although mutations in both ARSA and CHKB result in autosomal recessive disorders, their clinical features, neonatal manifestation of cardiomyopathy with increased creatine kinase activity found in CHKB-related muscular dystrophy15 or mental decline specific for juvenile form of metachromatic leukodystrophy, are distinct from symptoms observed in our patients.

In patient 1, we diagnosed SCO2 deficiency based on genetic analyses. With regard to knowledge of the human SCO2 protein structure, the newly identified hemizygous mutation p.D223N is likely to affect the coordination of copper (I) by the two cysteines of the CXXXC conserved motif and by the conserved His224.16 Furthermore, markedly decreased levels of SCO2 protein and border-low copper content found in tissues of patient 1 are in accordance with previously reported findings in SCO2 patients.13, 17 To date, we have diagnosed the late non-neonatal onset of SCO2 deficiency in eight patients,18 but ataxia was not noted in any of these or other reported SCO2 patients (Supplementary Table 2). All of the observed symptoms are typical of SCO2 deficiency. Interestingly, based on the review of published cases, HCMP, which is the key finding in SCO2-deficient patients with neonatal onset, is present only in about 50% of the SCO2 patients with non-neonatal presentation (Supplementary Table 2).

With regard to the biochemical and genetic findings, MNGIE was diagnosed in patient 2. The manifested symptoms of patient 2 were similar to those described previously.19 Whereas no pathogenic mutation was found in the SCO2 gene of patient 2, we assume that mildly decreased COX activity in his lymphocytes might be a consequence of point mtDNA-mutation accumulation due to diminished thymidine phosphorylase activity as reported previously in MNGIE patients.20

The number of mitochondrial proteins in mammals has been estimated to be between 1130 and 1500, but pathogenic defects have been found in only a fraction of these proteins.2 Additionally, there is no satisfactory explanation for the phenotype–genotype variability in patients suffering from mitochondrial diseases. However, we believe that the use of advanced whole-genome analysis techniques will shed light on the genetic causes of mitochondrial defects, as well as on the regulation of tissue-specific energy demands.9, 10 Further improvement in our knowledge of mechanisms based on the presence of variable gene dosage, epigenetic modifications, repetitive sequences, non-coding RNAs, cis- and trans-splicing, diverse penetrance of diseases, the complex effects of the hundreds of loss-of-function variants present in every individual, and general agreement on the classification of particular genetic variants will clarify the phenotypic diversity in humans.21, 22