Abstract

Clinical conditions commonly associated with mitochondrial disorders (CAMDs) are often present in autism spectrum disorders (ASD) and intellectual disability (ID). Therefore, the mitochondrial dysfunction hypothesis has been proposed as a transversal mechanism that may function in both disorders. Here, we investigated the presence of conditions associated with mitochondrial disorders and mitochondrial DNA (mtDNA) alterations in 122 subjects who presented ASD with ID (ASD group), 115 subjects who presented ID but not ASD (ID group) and 112 healthy controls (HC group). We assessed in the three study groups the presence of the clinical conditions through a questionnaire and the mtDNA content of two mitochondrial genes, MT-ND1 and MT-ND4, by qPCR. The mtDNA sequences of 98 ASD and 95 ID subjects were obtained by mtDNA-targeted next generation sequencing and analysed through the MToolBox pipeline to identify mtDNA mutations. Subjects with ASD and ID showed higher frequencies of constipation, edema, seizures, vision alterations, strabismus and sphincter incontinence than HCs subjects. ASD and ID subjects showed significantly lower mtDNA content than HCs in both MT-ND1 and MT-ND4 genes. In addition, we identified 49 putative pathogenic variants with a heteroplasmy level higher than 60%: 8 missense, 29 rRNA and 12 tRNA variants. A total of 28.6% of ASD and 30.5% of ID subjects carried at least one putative pathogenic mtDNA mutation. The high frequency of CAMDs, the low mtDNA content and the presence of putative pathogenic mtDNA mutations observed in both ASD and ID subjects are evidence of mitochondrial dysfunction in ASD and ID.

Introduction

Autism spectrum disorders (ASDs) are a group of severe neurodevelopmental disorders with an onset in the early stages of development. People affected by ASD tend to suffer persistent deficits in social communication, social interaction and restricted or repetitive patterns of behavior, interests, or activities (1). The median worldwide prevalence of ASD is 1–2%, and a similar prevalence is reported in the dataset of adults only (2). Intellectual disability (ID) is characterized by both intellectual and adaptive functional deficits in conceptual, social and practical domains, affecting 2–3% of the population (3). The onset of ID occurs during the developmental period. Both ASD and ID are strongly associated and show a high comorbidity (2,4). In children, 18 to 70% of patients with ASD also present ID and up to 40% of subjects with ID also present ASD (5,6). In adulthood, approximately one-fifth of subjects with ASD are diagnosed with ID (7).

ASD also occurs with other frequent medical comorbidities, such as epilepsy, sensory abnormalities, motor deficits, sleep abnormalities, gastrointestinal disturbances, feeding and nutritional problems, metabolic syndrome and immune conditions (7–9). In fact, more than 70% of subjects with ASD have concurrent medical, developmental or psychiatric disorders (2). Interestingly, many of the above mentioned clinical characteristics are conditions that are commonly associated with mitochondrial disorders (CAMDs). Mitochondrial disorders comprise a collection of individual rare syndromes produced by defects in either the mitochondrial DNA (mtDNA) or nuclear DNA that result in mitochondrial dysfunction. Impaired mitochondria affect most organ systems to different extents, producing a wide spectrum of phenotypes; however, little is known about the cascade of events involved in such diverse, rare syndromes (10,11). Interestingly, mitochondrial dysfunction has been described in both subjects with ASD (12,13) and ID (14) and has been suggested as a possible neurobiological subtype of ASD (15).

Mitochondria are key organelles required for energy production, and their proper function is therefore crucial for tissues with a high metabolic demand, such as the nervous system. Mitochondria have their own genome, known as mtDNA, which comprises a 16 569-bp, double-stranded, maternally inherited, circular molecule that contains the genetic information necessary for the synthesis of 13 essential polypeptides in the mitochondrial respiratory chain (16). Unlike nuclear DNA that is only present in two copies per cell, the mtDNA content varies, depending on the energy requirements. Each mitochondrion may present 2 to 10 mtDNA molecules, and each cell can contain several mitochondria (17). Interestingly, the mtDNA content is related to anthropometric indices (18) and has been proposed as a biomarker of mitochondrial dysfunction involved in conditions such as diabetes, obesity, cancer, HIV complications, aging (19) and autism (20).

The genetic architecture of ASD and ID is complex, with several types of genetic variants known to be associated with both conditions, although many of the genetic factors are unknown (3,4). Furthermore, despite the evidence of mitochondrial dysfunction in subjects with ASD and ID, the study of mtDNA has been neglected even it may contribute to the genetic bases of both conditions. Therefore, we tested the hypothesis that, in adulthood, CAMDs comorbidity is frequently present in ASD and ID compared with a control group. In addition, we investigated whether alterations in the mtDNA sequence were present in ASD and ID and whether the mtDNA content in ASD and ID cases differed from controls.

Results

CAMDs

We compared the frequency of CAMDs among the three groups (Table 1). Constipation, edema (referred to as venous insufficiency and other venous alterations), seizures, vision alterations, strabismus and sphincter incontinence were more frequent in the ASD group than in the HC group. All these conditions were also more frequent in the ID group than in the HC group. Other factors, such as temperature changes, diabetes, hypertension and cancer, were also more frequent in the ID group than in the HC group, but did not reach statistical significance. Conditions present in less than 5% of subjects that were not reported in the ASD and ID groups are: headache, migraine, diarrhea, abdominal pain, nausea, kinetosis, severe fatigue, fibromyalgia, dysautonomia, arthritis, muscular weakness, deafness, stroke, heart disease, hypercholesterolemia, kidney disease, hypoglycemia, and hypothyroidism.

Table 1.

Percentage of subjects in the three study groups presenting conditions commonly associated with mitochondrial disorders (CAMDs)

ConditionsHCn = 112ASDn = 122IDn = 115Statistics
Compared groupsχ2P
Specific conditions with frequency >5% in the study groups
 Constipation17.073.857.4HC vs. ASD75.403<0.001
HC vs. ID38.496<0.001
ASD vs. ID7.0600.008
 Edema3.615.627.8HC vs. ASD9.3570.002
HC vs. ID24.478<0.001
ASD vs. ID5.2620.021
 Seizures063.942.6HC vs. ASD107.059<0.001
HC vs. ID59.918<0.001
ASD vs. ID10.8200.001
 Vision alterations5.418.033.9HC vs. ASD8.7620.003
HC vs. ID28.458<0.001
ASD vs. ID7.8100.005
 Strabismus1.810.713.0HC vs. ASD7.6470.006
HC vs. ID10.2680.001
ASD vs. ID0.3280.567
 Sphincters incontinence066.431.3HC vs. ASD113.400<0.001
HC vs. ID40.994<0.001
ASD vs. ID29.160<0.001
Other conditions grouped into organ systems
 Nervous system15.2100100HC vs. ASD174.207<0.001
HC vs. ID167.747<0.001
ASD vs. ID
 Gastro-intestinal4.55.71.7HC vs. ASD0.1950.659
HC vs. ID0.6460.422
ASD vs. ID1.6120.204
 Cardio-vascular29.513.932.2HC vs. ASD5.2870.021
HC vs. ID0.1950.659
ASD vs. ID11.194<0.001
 Musculo-skeletal2.700HC vs. ASD
HC vs. ID
ASD vs. ID
 Endocrine12.58.217.4HC vs. ASD1.1750.278
HC vs. ID0.0710.790
ASD vs. ID4.527<0.033
ConditionsHCn = 112ASDn = 122IDn = 115Statistics
Compared groupsχ2P
Specific conditions with frequency >5% in the study groups
 Constipation17.073.857.4HC vs. ASD75.403<0.001
HC vs. ID38.496<0.001
ASD vs. ID7.0600.008
 Edema3.615.627.8HC vs. ASD9.3570.002
HC vs. ID24.478<0.001
ASD vs. ID5.2620.021
 Seizures063.942.6HC vs. ASD107.059<0.001
HC vs. ID59.918<0.001
ASD vs. ID10.8200.001
 Vision alterations5.418.033.9HC vs. ASD8.7620.003
HC vs. ID28.458<0.001
ASD vs. ID7.8100.005
 Strabismus1.810.713.0HC vs. ASD7.6470.006
HC vs. ID10.2680.001
ASD vs. ID0.3280.567
 Sphincters incontinence066.431.3HC vs. ASD113.400<0.001
HC vs. ID40.994<0.001
ASD vs. ID29.160<0.001
Other conditions grouped into organ systems
 Nervous system15.2100100HC vs. ASD174.207<0.001
HC vs. ID167.747<0.001
ASD vs. ID
 Gastro-intestinal4.55.71.7HC vs. ASD0.1950.659
HC vs. ID0.6460.422
ASD vs. ID1.6120.204
 Cardio-vascular29.513.932.2HC vs. ASD5.2870.021
HC vs. ID0.1950.659
ASD vs. ID11.194<0.001
 Musculo-skeletal2.700HC vs. ASD
HC vs. ID
ASD vs. ID
 Endocrine12.58.217.4HC vs. ASD1.1750.278
HC vs. ID0.0710.790
ASD vs. ID4.527<0.033

HC: healthy controls; ASD: autism spectrum disorder; ID: intellectual disability.

Significant differences are indicated in boldface.

Table 1.

Percentage of subjects in the three study groups presenting conditions commonly associated with mitochondrial disorders (CAMDs)

ConditionsHCn = 112ASDn = 122IDn = 115Statistics
Compared groupsχ2P
Specific conditions with frequency >5% in the study groups
 Constipation17.073.857.4HC vs. ASD75.403<0.001
HC vs. ID38.496<0.001
ASD vs. ID7.0600.008
 Edema3.615.627.8HC vs. ASD9.3570.002
HC vs. ID24.478<0.001
ASD vs. ID5.2620.021
 Seizures063.942.6HC vs. ASD107.059<0.001
HC vs. ID59.918<0.001
ASD vs. ID10.8200.001
 Vision alterations5.418.033.9HC vs. ASD8.7620.003
HC vs. ID28.458<0.001
ASD vs. ID7.8100.005
 Strabismus1.810.713.0HC vs. ASD7.6470.006
HC vs. ID10.2680.001
ASD vs. ID0.3280.567
 Sphincters incontinence066.431.3HC vs. ASD113.400<0.001
HC vs. ID40.994<0.001
ASD vs. ID29.160<0.001
Other conditions grouped into organ systems
 Nervous system15.2100100HC vs. ASD174.207<0.001
HC vs. ID167.747<0.001
ASD vs. ID
 Gastro-intestinal4.55.71.7HC vs. ASD0.1950.659
HC vs. ID0.6460.422
ASD vs. ID1.6120.204
 Cardio-vascular29.513.932.2HC vs. ASD5.2870.021
HC vs. ID0.1950.659
ASD vs. ID11.194<0.001
 Musculo-skeletal2.700HC vs. ASD
HC vs. ID
ASD vs. ID
 Endocrine12.58.217.4HC vs. ASD1.1750.278
HC vs. ID0.0710.790
ASD vs. ID4.527<0.033
ConditionsHCn = 112ASDn = 122IDn = 115Statistics
Compared groupsχ2P
Specific conditions with frequency >5% in the study groups
 Constipation17.073.857.4HC vs. ASD75.403<0.001
HC vs. ID38.496<0.001
ASD vs. ID7.0600.008
 Edema3.615.627.8HC vs. ASD9.3570.002
HC vs. ID24.478<0.001
ASD vs. ID5.2620.021
 Seizures063.942.6HC vs. ASD107.059<0.001
HC vs. ID59.918<0.001
ASD vs. ID10.8200.001
 Vision alterations5.418.033.9HC vs. ASD8.7620.003
HC vs. ID28.458<0.001
ASD vs. ID7.8100.005
 Strabismus1.810.713.0HC vs. ASD7.6470.006
HC vs. ID10.2680.001
ASD vs. ID0.3280.567
 Sphincters incontinence066.431.3HC vs. ASD113.400<0.001
HC vs. ID40.994<0.001
ASD vs. ID29.160<0.001
Other conditions grouped into organ systems
 Nervous system15.2100100HC vs. ASD174.207<0.001
HC vs. ID167.747<0.001
ASD vs. ID
 Gastro-intestinal4.55.71.7HC vs. ASD0.1950.659
HC vs. ID0.6460.422
ASD vs. ID1.6120.204
 Cardio-vascular29.513.932.2HC vs. ASD5.2870.021
HC vs. ID0.1950.659
ASD vs. ID11.194<0.001
 Musculo-skeletal2.700HC vs. ASD
HC vs. ID
ASD vs. ID
 Endocrine12.58.217.4HC vs. ASD1.1750.278
HC vs. ID0.0710.790
ASD vs. ID4.527<0.033

HC: healthy controls; ASD: autism spectrum disorder; ID: intellectual disability.

Significant differences are indicated in boldface.

mtDNA copy number

The mtDNA copy number, namely, the number of mtDNA molecules per white blood cell, differed among the three study groups (ASD, ID, and HC). Age, BMI and the number of cigarettes smoked per day were correlated with the mtDNA copy number; therefore, these factors were included as covariates in the general linear analysis comparing the mtDNA content between groups. The mtDNA copy number measured using either the MT-ND1 or the MT-ND4 genes was lower in the ASD and ID groups than in the HC group and was also significantly lower in the ID group than in the ASD group (Fig. 1). The deletion ratio was higher in both the ASD and ID groups than in the HC group, although significant differences were only present when the HC and ASD groups were compared.

Estimated marginal means of the MT-ND1 and MT-ND4 genes and the deletion ratio in the three study groups. HC: healthy controls; ASD: autism spectrum disorder; ID: intellectual disability; n.s.: not significant; *P < 0.05; **P < 0.01; ***P < 0.001.
Figure 1.

Estimated marginal means of the MT-ND1 and MT-ND4 genes and the deletion ratio in the three study groups. HC: healthy controls; ASD: autism spectrum disorder; ID: intellectual disability; n.s.: not significant; *P < 0.05; **P < 0.01; ***P < 0.001.

Putative pathogenic mtDNA variants

Among the 193 sequenced subjects, MToolBox (21) identified 664 variants of clinical interest. After applying the selection criteria for pathogenicity prediction, we identified 57 (29.5%) subjects who carried 49 putative pathogenic variants, with a heteroplasmy level ≥ 60%. In the ASD group, 28 subjects carried 37 variants, whereas 29 subjects in the ID group carried 28 variants (Table 2). Notably, 11 variants were present in more than one subject and 12 subjects carried more than one putative pathogenic variant.

Table 2.

Putatively pathogenic variants present in the study sample

NVariantAA changeFreq.dbSNP IDSamplesLocusHomo/ heteroplasmyDisease score
Non-synonymous variants
 1m.8663A>CQ46P0ID074MT-ATP60.790.889
 2m.8936T>AL137H0A003, ID009, ID036, ID055, ID069, ID073, ID074, ID104, ID1120.60–0.970.855
 3m.11039C>TL94F0ID071MT-ND40.990.702
 4m.11268C>TT170I0.0010A095, ID01510.739
 5m.13855C>AL507M0ID070MT-ND510.686
 6m.15132T>CM129T0.0002A117MT-CYB10.821
 7m.15482T>CS246P0ID03810.559
 8m.15597T>CV284A0A05910.581
rRNA variants
 1m.680T>C0.0006A099MT-RNR11
 2m.720T>C0A0700.99
 3m.736C>T0ID0571
 4m.851A>G0rs28502491A092, ID0110.96–1
 5m.870C>T0.0005ID1161
 6m.1007G>A0.0010rs111033213A0510.99
 7m.1120C>T0.0005rs727505171ID0071
 8m.1211G>A0.0010rs397515725A063, A095, ID015, ID0421
 9m.1303G>A0.0023A0211
 10m.1537C>T0.0003ID1081
 11m.1555A>Ga0.0032rs267606617A043, A054, ID0180.99–1
 12m.1717T>Ca0.0007ID076MT-RNR21
 13m.1821A>G0.0008ID1051
 14m.1834T>C0.0007A078, A089, A1561
 15m.2141T>C0.0031ID1121
 16m.2222T>C0ID061, ID0871
 17m.2223A>G0ID0381
 18m.2234C>T0ID0171
 19m.2244T>G0ID0111
 20m.2398A>G0A0681
 21m.2407T>C0ID0241
 22m.2524A>C0.0001A0961
 23m.2702G>A0.0020A0671
 24m.2746T>C0.0005rs28663331A0901
 25m.2775A>G0ID061, ID0871
 26m.2792A>C0.0001A0571
 27m.2969A>G0.0002ID0651
 28m.3203A>G0.0009ID0551
 29m.3213A>Ga0.0003A0981
tRNA variants
 1m.4310A>G0.0008ID085MT-TI1
 2m.4336T>C0.0220rs41456348A059, ID004MT-TQ0.98–10.4
 3m.4388A>G0.0002rs375986475ID036MT-TQ1
 4m.4435A>Ga0.0010A079, ID037MT-TM0.95–10.35
 5m.5605A>G0.0003A152MT-TA1
 6m.10007T>C0.0010rs201906571A009MT-TG1
 7m.12142A>Ga0.0001A125MT-TH1
 8m.12193A>G0A0841
 9m.14684C>T0.0009A003MT-TE1
 10m.14687A>G0.0098rs200189658ID11210.65
 11m.15913C>T0.0003ID065MT-TT0.99
 12m.15946C>T0.0032rs202014122A014, A0450.99–1
NVariantAA changeFreq.dbSNP IDSamplesLocusHomo/ heteroplasmyDisease score
Non-synonymous variants
 1m.8663A>CQ46P0ID074MT-ATP60.790.889
 2m.8936T>AL137H0A003, ID009, ID036, ID055, ID069, ID073, ID074, ID104, ID1120.60–0.970.855
 3m.11039C>TL94F0ID071MT-ND40.990.702
 4m.11268C>TT170I0.0010A095, ID01510.739
 5m.13855C>AL507M0ID070MT-ND510.686
 6m.15132T>CM129T0.0002A117MT-CYB10.821
 7m.15482T>CS246P0ID03810.559
 8m.15597T>CV284A0A05910.581
rRNA variants
 1m.680T>C0.0006A099MT-RNR11
 2m.720T>C0A0700.99
 3m.736C>T0ID0571
 4m.851A>G0rs28502491A092, ID0110.96–1
 5m.870C>T0.0005ID1161
 6m.1007G>A0.0010rs111033213A0510.99
 7m.1120C>T0.0005rs727505171ID0071
 8m.1211G>A0.0010rs397515725A063, A095, ID015, ID0421
 9m.1303G>A0.0023A0211
 10m.1537C>T0.0003ID1081
 11m.1555A>Ga0.0032rs267606617A043, A054, ID0180.99–1
 12m.1717T>Ca0.0007ID076MT-RNR21
 13m.1821A>G0.0008ID1051
 14m.1834T>C0.0007A078, A089, A1561
 15m.2141T>C0.0031ID1121
 16m.2222T>C0ID061, ID0871
 17m.2223A>G0ID0381
 18m.2234C>T0ID0171
 19m.2244T>G0ID0111
 20m.2398A>G0A0681
 21m.2407T>C0ID0241
 22m.2524A>C0.0001A0961
 23m.2702G>A0.0020A0671
 24m.2746T>C0.0005rs28663331A0901
 25m.2775A>G0ID061, ID0871
 26m.2792A>C0.0001A0571
 27m.2969A>G0.0002ID0651
 28m.3203A>G0.0009ID0551
 29m.3213A>Ga0.0003A0981
tRNA variants
 1m.4310A>G0.0008ID085MT-TI1
 2m.4336T>C0.0220rs41456348A059, ID004MT-TQ0.98–10.4
 3m.4388A>G0.0002rs375986475ID036MT-TQ1
 4m.4435A>Ga0.0010A079, ID037MT-TM0.95–10.35
 5m.5605A>G0.0003A152MT-TA1
 6m.10007T>C0.0010rs201906571A009MT-TG1
 7m.12142A>Ga0.0001A125MT-TH1
 8m.12193A>G0A0841
 9m.14684C>T0.0009A003MT-TE1
 10m.14687A>G0.0098rs200189658ID11210.65
 11m.15913C>T0.0003ID065MT-TT0.99
 12m.15946C>T0.0032rs202014122A014, A0450.99–1

N: nucleotide; AA: amino acid; Freq. frequency; dbSNP: Single Nucleotide Polymorphism data base; A: autism spectrum disorder; ID: intellectual disability; MT: mitochondrial; CO1: cytochrome c oxidase subunit 1; ATP6: ATP synthase subunit 6; ND3: NADH dehydrogenase subunit 3; CYB: cytochrome b; RNR1: 12S ribosomal RNA; RNR2: 16S ribosomal RNA; TF = tRNA-Phenylalanine; TI = tRNA-Isoleucine; TQ = tRNA-Glutamine; TM = tRNA-Metionine; TA = tRNA-Alanine; TD: tRNA-Aspartic acid; TG = tRNA-Glycine; TH = tRNA-Histidine; TE = tRNA-Glutamic; TT = tRNA-Threonine.

a

Reported by Wang et al. (36).

Frequency of the variants was retrieved from 8686 complete European human genomes present in HmtDB in July 2017 (44).

Table 2.

Putatively pathogenic variants present in the study sample

NVariantAA changeFreq.dbSNP IDSamplesLocusHomo/ heteroplasmyDisease score
Non-synonymous variants
 1m.8663A>CQ46P0ID074MT-ATP60.790.889
 2m.8936T>AL137H0A003, ID009, ID036, ID055, ID069, ID073, ID074, ID104, ID1120.60–0.970.855
 3m.11039C>TL94F0ID071MT-ND40.990.702
 4m.11268C>TT170I0.0010A095, ID01510.739
 5m.13855C>AL507M0ID070MT-ND510.686
 6m.15132T>CM129T0.0002A117MT-CYB10.821
 7m.15482T>CS246P0ID03810.559
 8m.15597T>CV284A0A05910.581
rRNA variants
 1m.680T>C0.0006A099MT-RNR11
 2m.720T>C0A0700.99
 3m.736C>T0ID0571
 4m.851A>G0rs28502491A092, ID0110.96–1
 5m.870C>T0.0005ID1161
 6m.1007G>A0.0010rs111033213A0510.99
 7m.1120C>T0.0005rs727505171ID0071
 8m.1211G>A0.0010rs397515725A063, A095, ID015, ID0421
 9m.1303G>A0.0023A0211
 10m.1537C>T0.0003ID1081
 11m.1555A>Ga0.0032rs267606617A043, A054, ID0180.99–1
 12m.1717T>Ca0.0007ID076MT-RNR21
 13m.1821A>G0.0008ID1051
 14m.1834T>C0.0007A078, A089, A1561
 15m.2141T>C0.0031ID1121
 16m.2222T>C0ID061, ID0871
 17m.2223A>G0ID0381
 18m.2234C>T0ID0171
 19m.2244T>G0ID0111
 20m.2398A>G0A0681
 21m.2407T>C0ID0241
 22m.2524A>C0.0001A0961
 23m.2702G>A0.0020A0671
 24m.2746T>C0.0005rs28663331A0901
 25m.2775A>G0ID061, ID0871
 26m.2792A>C0.0001A0571
 27m.2969A>G0.0002ID0651
 28m.3203A>G0.0009ID0551
 29m.3213A>Ga0.0003A0981
tRNA variants
 1m.4310A>G0.0008ID085MT-TI1
 2m.4336T>C0.0220rs41456348A059, ID004MT-TQ0.98–10.4
 3m.4388A>G0.0002rs375986475ID036MT-TQ1
 4m.4435A>Ga0.0010A079, ID037MT-TM0.95–10.35
 5m.5605A>G0.0003A152MT-TA1
 6m.10007T>C0.0010rs201906571A009MT-TG1
 7m.12142A>Ga0.0001A125MT-TH1
 8m.12193A>G0A0841
 9m.14684C>T0.0009A003MT-TE1
 10m.14687A>G0.0098rs200189658ID11210.65
 11m.15913C>T0.0003ID065MT-TT0.99
 12m.15946C>T0.0032rs202014122A014, A0450.99–1
NVariantAA changeFreq.dbSNP IDSamplesLocusHomo/ heteroplasmyDisease score
Non-synonymous variants
 1m.8663A>CQ46P0ID074MT-ATP60.790.889
 2m.8936T>AL137H0A003, ID009, ID036, ID055, ID069, ID073, ID074, ID104, ID1120.60–0.970.855
 3m.11039C>TL94F0ID071MT-ND40.990.702
 4m.11268C>TT170I0.0010A095, ID01510.739
 5m.13855C>AL507M0ID070MT-ND510.686
 6m.15132T>CM129T0.0002A117MT-CYB10.821
 7m.15482T>CS246P0ID03810.559
 8m.15597T>CV284A0A05910.581
rRNA variants
 1m.680T>C0.0006A099MT-RNR11
 2m.720T>C0A0700.99
 3m.736C>T0ID0571
 4m.851A>G0rs28502491A092, ID0110.96–1
 5m.870C>T0.0005ID1161
 6m.1007G>A0.0010rs111033213A0510.99
 7m.1120C>T0.0005rs727505171ID0071
 8m.1211G>A0.0010rs397515725A063, A095, ID015, ID0421
 9m.1303G>A0.0023A0211
 10m.1537C>T0.0003ID1081
 11m.1555A>Ga0.0032rs267606617A043, A054, ID0180.99–1
 12m.1717T>Ca0.0007ID076MT-RNR21
 13m.1821A>G0.0008ID1051
 14m.1834T>C0.0007A078, A089, A1561
 15m.2141T>C0.0031ID1121
 16m.2222T>C0ID061, ID0871
 17m.2223A>G0ID0381
 18m.2234C>T0ID0171
 19m.2244T>G0ID0111
 20m.2398A>G0A0681
 21m.2407T>C0ID0241
 22m.2524A>C0.0001A0961
 23m.2702G>A0.0020A0671
 24m.2746T>C0.0005rs28663331A0901
 25m.2775A>G0ID061, ID0871
 26m.2792A>C0.0001A0571
 27m.2969A>G0.0002ID0651
 28m.3203A>G0.0009ID0551
 29m.3213A>Ga0.0003A0981
tRNA variants
 1m.4310A>G0.0008ID085MT-TI1
 2m.4336T>C0.0220rs41456348A059, ID004MT-TQ0.98–10.4
 3m.4388A>G0.0002rs375986475ID036MT-TQ1
 4m.4435A>Ga0.0010A079, ID037MT-TM0.95–10.35
 5m.5605A>G0.0003A152MT-TA1
 6m.10007T>C0.0010rs201906571A009MT-TG1
 7m.12142A>Ga0.0001A125MT-TH1
 8m.12193A>G0A0841
 9m.14684C>T0.0009A003MT-TE1
 10m.14687A>G0.0098rs200189658ID11210.65
 11m.15913C>T0.0003ID065MT-TT0.99
 12m.15946C>T0.0032rs202014122A014, A0450.99–1

N: nucleotide; AA: amino acid; Freq. frequency; dbSNP: Single Nucleotide Polymorphism data base; A: autism spectrum disorder; ID: intellectual disability; MT: mitochondrial; CO1: cytochrome c oxidase subunit 1; ATP6: ATP synthase subunit 6; ND3: NADH dehydrogenase subunit 3; CYB: cytochrome b; RNR1: 12S ribosomal RNA; RNR2: 16S ribosomal RNA; TF = tRNA-Phenylalanine; TI = tRNA-Isoleucine; TQ = tRNA-Glutamine; TM = tRNA-Metionine; TA = tRNA-Alanine; TD: tRNA-Aspartic acid; TG = tRNA-Glycine; TH = tRNA-Histidine; TE = tRNA-Glutamic; TT = tRNA-Threonine.

a

Reported by Wang et al. (36).

Frequency of the variants was retrieved from 8686 complete European human genomes present in HmtDB in July 2017 (44).

Non-synonymous variants

In the present study, eight non-synonymous variants fulfilled the cut-off prioritization criteria previously described as pathogenic. All were missense variants, none of them have been previously associated with a disease or health condition and, notably, three of them were not previously reported: m.8663A > C and m.8936T > A in the MT-ATP6 gene and m.15597T > C in the MT-CYB gene.

rRNA variants

We identified 29 putative pathogenic rRNA variants based on the criterion of nucleotide variability ≤ 0.0097 and a heteroplasmy level higher than 60%. All variants, with the exception of m.2398A > G, have been previously reported. Interestingly, two of the 29 variants were reported to be associated with disease in MITOMAP: m.1537C > T was reported to be associated with maternally inherited or aminoglycoside-induced deafness, and the subject bearing this mutation actually suffered from this condition; and m.1555A > G was confirmed to be associated with deafness in MITOMAP, ClinVar and OMIM (http://omim.org/entry/561000#0001). However, only the oldest of the three subjects presenting the variant showed deafness (Table 3).

Table 3.

Clinical characteristics reported to be associated with mtDNA variants present in the study sample

VariantLocusHeteroplasmySubjectsAgeGenderDeafnessMigraineAD/PDHypertensionVisual lossMMRF
m.1537C>TMT-RNR11ID10860M+NE
m.1555A>GMT-RNR11ID01833MNE+
1A04343FNE
1A05454F+NE
m.4336T>CMT-TQ0.98A05938MNE+
1ID00463FNE+
m.4388A>GMT-TQ1ID03636FNE+
m.4435A>GMT-TM1ID03740FNE+
0.96A07939MNE+
m.14687A>GMT-TE1ID11221MNE++
VariantLocusHeteroplasmySubjectsAgeGenderDeafnessMigraineAD/PDHypertensionVisual lossMMRF
m.1537C>TMT-RNR11ID10860M+NE
m.1555A>GMT-RNR11ID01833MNE+
1A04343FNE
1A05454F+NE
m.4336T>CMT-TQ0.98A05938MNE+
1ID00463FNE+
m.4388A>GMT-TQ1ID03636FNE+
m.4435A>GMT-TM1ID03740FNE+
0.96A07939MNE+
m.14687A>GMT-TE1ID11221MNE++

ID: intellectual disability; A: autism spectrum disorder; M: male; F: female; NE: not evaluable; AD: Alzheimer’s disease; PD: Parkinson’s disease; MMRF: Mitochondrial Myopathy and Respiratory Failure; MT: mitochondrial; RNR1: 12S ribosomal RNA; TQ: tRNA-Glutamine; TM: tRNA-Metionine; TE = tRNA-Glutamic acid.

Plus (+) or minus (−) symbol indicates the presence or absence of the clinical characteristic in the study subject.

Reported clinical conditions with each variant are shaded: 1537 C > T, deafness (MITOMAP); 1555 A > G, maternally inherited or aminoglycoside-induced deafness (MITOMAP); 4336 T > C, AD and PD, hearing loss and migraine (MITOMAP), sensorineural deafness and migraine (ClinVar and OMIM); 4388 A > G, possible hypertension factor (MITOMAP); 4435 A > G, hypertension and modifying factor of Leber’s hereditary optic neuropathy phenotype (MITOMAP); 14687 A > G, MMRF (MITOMAP).

Other conditions not reported to be associated with the variants but present in the study subjects: epilepsy (ID108, A043, A054, A059, ID004, ID112); constipation (ID018, A043, A054, A059, ID036, ID037, ID112); cardiac disease (A059, ID004) and diabetes (ID112).

Table 3.

Clinical characteristics reported to be associated with mtDNA variants present in the study sample

VariantLocusHeteroplasmySubjectsAgeGenderDeafnessMigraineAD/PDHypertensionVisual lossMMRF
m.1537C>TMT-RNR11ID10860M+NE
m.1555A>GMT-RNR11ID01833MNE+
1A04343FNE
1A05454F+NE
m.4336T>CMT-TQ0.98A05938MNE+
1ID00463FNE+
m.4388A>GMT-TQ1ID03636FNE+
m.4435A>GMT-TM1ID03740FNE+
0.96A07939MNE+
m.14687A>GMT-TE1ID11221MNE++
VariantLocusHeteroplasmySubjectsAgeGenderDeafnessMigraineAD/PDHypertensionVisual lossMMRF
m.1537C>TMT-RNR11ID10860M+NE
m.1555A>GMT-RNR11ID01833MNE+
1A04343FNE
1A05454F+NE
m.4336T>CMT-TQ0.98A05938MNE+
1ID00463FNE+
m.4388A>GMT-TQ1ID03636FNE+
m.4435A>GMT-TM1ID03740FNE+
0.96A07939MNE+
m.14687A>GMT-TE1ID11221MNE++

ID: intellectual disability; A: autism spectrum disorder; M: male; F: female; NE: not evaluable; AD: Alzheimer’s disease; PD: Parkinson’s disease; MMRF: Mitochondrial Myopathy and Respiratory Failure; MT: mitochondrial; RNR1: 12S ribosomal RNA; TQ: tRNA-Glutamine; TM: tRNA-Metionine; TE = tRNA-Glutamic acid.

Plus (+) or minus (−) symbol indicates the presence or absence of the clinical characteristic in the study subject.

Reported clinical conditions with each variant are shaded: 1537 C > T, deafness (MITOMAP); 1555 A > G, maternally inherited or aminoglycoside-induced deafness (MITOMAP); 4336 T > C, AD and PD, hearing loss and migraine (MITOMAP), sensorineural deafness and migraine (ClinVar and OMIM); 4388 A > G, possible hypertension factor (MITOMAP); 4435 A > G, hypertension and modifying factor of Leber’s hereditary optic neuropathy phenotype (MITOMAP); 14687 A > G, MMRF (MITOMAP).

Other conditions not reported to be associated with the variants but present in the study subjects: epilepsy (ID108, A043, A054, A059, ID004, ID112); constipation (ID018, A043, A054, A059, ID036, ID037, ID112); cardiac disease (A059, ID004) and diabetes (ID112).

tRNA variants

We identified 12 putative pathogenic tRNA variants, four of them have been previously associated with several conditions; deafness, Alzheimer’s disease, hypertension and visual loss were present in some study subjects, but mitochondrial myopathy and respiratory failure was not (Table 3).

mtDNA variants of uncertain significance

In addition, we have identified other putative pathogenic variants with a heteroplasmy percentage ≤ 60% (Supplementary Material, Table S1). Furthermore, several variants were also present in subjects with ASD and ID. Among them, it is worth mentioning variants located in the MT-D-LOOP and MT-ORIL non-previously reported or with low nucleotide variability (Supplementary Material, Table S2). Some variants have previously been reported to be associated with various conditions in MITOMAP (Supplementary Material, Table S3), although they did not fulfill the selection criteria for identification as putatively pathogenic using the MToolBox priorization criteria.

Each ASD and ID subject was assigned to a European mitochondrial haplogroup (Supplementary Material, Table S4). The frequencies of the most common haplogroups in the two study groups, ASD and ID, did not differ from the reported frequencies in healthy subjects (data not shown) selected from a population-based sample in the same geographic region (22). Furthermore, putative pathogenic variants were not specific for a single haplotype, as we observed the following ratios of putative pathogenic variants within a haplogroup/number of total subjects within the same haplogroup: 23/91H, 3/18J, 2/12K, 2/5L, 2/4M, 2/3N, 7/16T, 13/32U, 2/10V and 1/6W.

Datasets of mtDNA sequences and clinical data from participants are available at the European Genome-phenome Archive (EGA, https://ega-archive.org) with the study reference number EGAS00001002750.

Discussion

CAMDs

Few studies have been performed on the physical health of adults with ASD, although some medical comorbidities that are present in infancy persist into adulthood (9,23). Subjects with ID experience more extensive psychiatric and behavioral disorders and other medical conditions in addition to visual, dental and hearing problems (24). However, the percentages of comorbid conditions in ASD and ID studies vary, depending on the methodology of the study, and are often difficult to compare. Even among ASD studies, the results are difficult to compare due to different percentages of ID among study subjects (7,25). In the present study, we identified constipation, edema, seizures, vision alterations, strabismus and sphincter incontinence as the most prevalent comorbid conditions present in both subjects with ASD and ID compared with HCs. Constipation, seizures and sphincter incontinence were also more prevalent in subjects with ASD than in subjects with ID, whereas edema and vision alterations were more prevalent in subjects with ID than in subjects with ASD. Notably, in our study subjects with ID were on average 12 years older than subjects with ASD, and, therefore, we cannot exclude the possibility that edema and vision alterations were more prevalent in subjects with ID because they were older. Moreover, ID was more severe in patients who fulfilled the criteria for ASD than in subjects who did not, and, therefore, the high percentage of subjects presenting constipation, seizures and sphincter incontinence may be associated with the severe ID rather than ASD.

mtDNA copy number

We analysed MT-ND1 and MT-ND4, two of the seven mtDNA-encoded genes of complex I. Both subjects with ASD and ID presented lower MT-ND1 and MT-ND4 contents than HC subjects, indicating that the number of mtDNA copies was reduced in subjects with ASD and ID compared with that in HCs. Previous studies conducted in peripheral blood cells from children with ASD have reported an elevated mtDNA copy number compared with that in healthy children (26) and unaffected siblings (27). Notably, all these studies were conducted in children, whereas our study was conducted in adult subjects. The peripheral blood mtDNA content is associated with sex and age and is influenced by several variables, such as platelet and white blood cell counts and estroprogesterone intake (28). Regarding age, the mtDNA content increases until the fifth decade of life and then declines in older subjects (28). In the present study, we compared the mtDNA content after correcting for covariates that correlate with the mtDNA content, and we conclusively identified lower mtDNA content in subjects with ASD and ID than in HCs. The low peripheral blood mtDNA content in ASD and ID suggests a transversal mitochondrial dysfunction that can affect several tissues.

Regarding the deletion ratio, which is directly related to the presence of the common deletion, we observed that patients with ASD exhibited significantly higher ratios than HCs. Therefore, the common deletion was more frequent in the ASD group than in the HC group. We included all the confounding factors in the analysis, and significant differences were only observed for the ASD group. In a study that included 10 children with ASD, deletions of the MT-ND4 gene were reported in five subjects, and the deletion ratio was higher in 2 of 10 children (29).

mtDNA variants

Mitochondrial function, structure and cell size are altered by mtDNA heteroplasmy (30). Since the disease manifests at the heteroplasmic threshold in the range of 0.60 to 0.95 for most pathogenic mutations (31), we selected putative pathogenic variants exhibiting heteroplasmy levels higher than 60%. The present study identified that 26.6% of ASD patients and 30.5% of ID subjects had a putative pathogenic mutation suggesting that part of the genetic risk factors involved in ASD and ID can be located in mtDNA. Mitochondrial dysfunction and mtDNA implications have been broadly proposed in the etiology of ASD, however, only a few studies have analysed the mtDNA sequence in subjects with ASD (32–35). Recently, Wang et al. primarily compared the mtDNA sequence between children with ASD and their non-autistic siblings and identified that predicted pathogenic mutations conferred an increased risk of ASD. In addition, this risk was most pronounced in families with probands who exhibited a diminished intelligence quotient and/or impaired social behavior (36). The authors used whole-exome sequencing, unlike our study that performed mtDNA-targeted deep sequencing, and therefore, we obtained a higher mean deep coverage of 512X than the 141X coverage obtained in the previous study. Additionally, we analyse the mtDNA data using the MToolBox phylogeny-based prioritization workflow which uses the mtDNA sequences from 14 144 healthy individuals to select the variants that are more likely to be pathogenic (37). Most of the putative pathogenic variants we identified have not been previously been associated with ASD, with the exception of m.1555A > G, MT-RNR1; m.1717T > C, MT-RNR2; m.3213A > G, MT-RNR2; m.4435A > G, MT-TM and m.12142A > G, MT-TH that were identified by Wang et al. Each variant was present in a unique family; however, high-confidence heteroplasmic data were not reported for all three members in four out the five variants, and only 1717T > C, 4435A > G and 12142A > G were reported in ASD patients while no data were available for the other two variants. Regarding ID, our study is the first to analyse the entire mtDNA of subjects with ID using next generation sequencing via mtDNA-targeted deep sequencing. We reported a large number of putative pathogenic mtDNA variants; however, we could not confirm that all variants were implicated in ID or ASD because pathogenicity must be demonstrated by functional studies. In addition, high percentages of pathogenic mtDNA mutations have been described in the general population. Approximately one in 200 healthy subjects carries a pathogenic mtDNA variant; however, the heteroplasmy levels were less than 0.60 (38). Recently, the analyses of the 1000 Genomes Project have identified that 20% of subjects harbor heteroplasmies reported to be implicated in disease although positions with heteroplasmy levels larger than 0.60 showed a reduction in pathogenicity (39). The present study identified that one in four subjects with ASD or ID carry a pathogenic mtDNA variant with heteroplasmy > 0.60. Therefore, in addition to the known nuclear genetic alterations involved in ASD and ID, we identified mtDNA variants that may contribute to the genetic architecture of both clinical conditions.

It is worth mentioning that six variants present in four subjects with ASD (4%) and six subjects with ID (6%) have previously been associated with clinical conditions; all variants were present in rRNA and tRNA genes and had heteroplasmy levels higher than 0.95. Some of the clinical conditions were present in the subjects, whereas others were not. For example, the m.4336T > C has been reported to be associated with Alzheimer’s disease, and one of the two carriers, who was 63 years old, suffered from this condition. We cannot exclude the possibility that the other subject may develop the disease in the future, since at the time of the study this subject was 38 years old. Similarly, the m.1555A > G variant is associated with deafness, and the oldest subject of the three carriers was deaf, whereas the other two younger subjects were not. These two patients are predisposed to aminoglycoside-induced ototoxicity, and, therefore, aminoglycoside treatments must be avoided. Regarding tRNA variants, they are reported in MITOMAP as ‘possibly benign’ or ‘likely benign’ based on database frequencies, nature of the nucleotide change and conservation score. However, clinical history, heteroplasmy data or functional studies are not included in retrieving tRNA pathogenicity scores and, therefore, they should be further evaluated to discard pathogenicity.

Much additional work is needed to elucidate the role of putative pathogenic mtDNA mutations in mitochondrial function and to determine whether these variants are related to ID, ASD, and/or CAMDs. A limitation of the present study is that some of the rare rRNA and tRNA variants for which RNA disease scores are not yet available may be specific to our population and not related with clinical conditions. Finally, we studied ASD subjects with concurrent ID because our sample was obtained from institutionalized subjects who presented with low-functioning ASD. Therefore, our results should not be generalized to all subjects with ASD, although they may be generalized to subjects with ID.

In summary, our study identified 1) a high frequency of conditions commonly associated with mitochondrial disorders that occurred concomitantly with 2) low mtDNA content and 3) the presence putative pathogenic mtDNA mutations in subjects with ASD and ID. Based on these findings, subjects with ASD and ID may present mitochondrial dysfunction. Future research will elucidate whether the screening for mtDNA alterations in ASD and ID can be as appropriate as the screening for nuclear DNA alterations.

Materials and Methods

Study design

In this cross-sectional study, DNA samples and clinical data from adult Caucasian white subjects recruited at the Intellectual Disability and Developmental Disorders Research Unit of the Group Pere Mata in Reus (www.peremata.com, Catalonia, Spain) and from adult healthy controls were analysed. All participants were over 18 years old. Informed written consent was obtained from the relatives or other legal guardians of all institutionalized subjects and from healthy individuals. The Ethical Committee of Hospital Sant Joan de Reus approved the study.

Participants

The sample consisted of three groups: 122 subjects with a severe or profound intellectual disability who fulfilled the diagnostic criteria for autism spectrum disorders (ASD), 115 subjects with a severe or profound ID who did not fulfill the diagnostic criteria for ASD, and 112 healthy controls (HCs) (Table 4). ASD and ID were diagnosed according to the Diagnostic and Statistical Manual of Mental Disorders, 4th Edition, Text Revision (DSM-IV-TR) criteria. The ASD diagnosis was confirmed using the Autism Diagnostic Interview-Revised (ADI-R) (40). The Childhood Autism Rating Scale (CARS) (41) was used to measure autism severity in the ASD group and to explore autism features in the ID group. Given the difficulty in diagnosing ASD in persons with a severe and profound intellectual disability, only patients showing both evaluations (DSM-IV-TR and ADI-R) positive were included in the ASD group.

Table 4.

Characteristics of the sample

HCN = 112ASDN = 122IDN = 115Compared groupsStatistics
Gender
 Male, N (%)50 (44.6)80 (65.5)74 (64.3)HC vs. ASDχ2=10.4; P=0.001
HC vs. IDχ2=13.0; P=0.003
 Female, N (%)62 (55.4)42 (34.4)41 (35.7)
ASD vs. IDχ2=0.1; P=0.843
Age in years (mean, SD)42.4 (11.4)40.7 (8.3)52.1 (12.4)HC vs. ASDt= 1.3; P=0.194
HC vs. IDt = −6.1; P<0.001
ASD vs. IDt = −8.3; P<0.001
BMI in kg/m2 (mean, SD)24.1 (3.1)25.1 (4.9)27.3 (5.5)HC vs. ASDt = −1.7; P=0.086
HC vs. IDt = −5.4; P<0.001
ASD vs. IDt = −3.3; p<0.001
Tobacco consumption in c/d, (mean, SD)3.2 (6.5)0.2 (1.2)0.8 (3.4)HC vs. ASDt = 4.9; P<0.001
HC vs. IDt = 3.4; P=0.001
ASD vs. IDt = −2.1; P=0.041
Chromosomal alterations, N (%)37 (30.3)10 (8.7)ASD vs. IDχ2=17.7; P<0.001
ID Severity
 Severe, N (%)16 (13.1)69 (60%)ASD vs. IDχ2=55.3; P<0.001
 Profound, N (%)106 (86.9)46 (40%)
CARS score (mean, SD)42.6 (5.9)21.8 (4.4)ASD vs. IDt=28.9; P<0.001
Treatment, N (%)
 Antipsychotic66 (54.1)64 (55.6)ASD vs. IDχ2=0.1; P=0.810
 Anticholinergic16 (13.1)17 (14.8)χ2=0.1; P=0.711
 Benzodiazepines72 (59.0)61 (53.0)χ2=0.9; P=0.354
 Antidepressant11 (9.0)18 (15.7)χ2=2.4; P=0.119
 Antiepileptic50 (41.0)33 (28.7)χ2=3.9; P=0.047
 Mood stabilizers57 (46.7)49 (42.6)χ2=0.4; P=0.525
 Antihypertensive2 (1.6)13 (12.7)χ2=9.3; P=0.002
 Analgesics15 (12.3)21 (21.3)χ2=1.6; P=0.201
 Laxative51 (41.8)61 (53.0)χ2=3.0; P=0.083
 Antithyroid2 (1.6)4 (3.5)χ2=0.8; P=0.368
HCN = 112ASDN = 122IDN = 115Compared groupsStatistics
Gender
 Male, N (%)50 (44.6)80 (65.5)74 (64.3)HC vs. ASDχ2=10.4; P=0.001
HC vs. IDχ2=13.0; P=0.003
 Female, N (%)62 (55.4)42 (34.4)41 (35.7)
ASD vs. IDχ2=0.1; P=0.843
Age in years (mean, SD)42.4 (11.4)40.7 (8.3)52.1 (12.4)HC vs. ASDt= 1.3; P=0.194
HC vs. IDt = −6.1; P<0.001
ASD vs. IDt = −8.3; P<0.001
BMI in kg/m2 (mean, SD)24.1 (3.1)25.1 (4.9)27.3 (5.5)HC vs. ASDt = −1.7; P=0.086
HC vs. IDt = −5.4; P<0.001
ASD vs. IDt = −3.3; p<0.001
Tobacco consumption in c/d, (mean, SD)3.2 (6.5)0.2 (1.2)0.8 (3.4)HC vs. ASDt = 4.9; P<0.001
HC vs. IDt = 3.4; P=0.001
ASD vs. IDt = −2.1; P=0.041
Chromosomal alterations, N (%)37 (30.3)10 (8.7)ASD vs. IDχ2=17.7; P<0.001
ID Severity
 Severe, N (%)16 (13.1)69 (60%)ASD vs. IDχ2=55.3; P<0.001
 Profound, N (%)106 (86.9)46 (40%)
CARS score (mean, SD)42.6 (5.9)21.8 (4.4)ASD vs. IDt=28.9; P<0.001
Treatment, N (%)
 Antipsychotic66 (54.1)64 (55.6)ASD vs. IDχ2=0.1; P=0.810
 Anticholinergic16 (13.1)17 (14.8)χ2=0.1; P=0.711
 Benzodiazepines72 (59.0)61 (53.0)χ2=0.9; P=0.354
 Antidepressant11 (9.0)18 (15.7)χ2=2.4; P=0.119
 Antiepileptic50 (41.0)33 (28.7)χ2=3.9; P=0.047
 Mood stabilizers57 (46.7)49 (42.6)χ2=0.4; P=0.525
 Antihypertensive2 (1.6)13 (12.7)χ2=9.3; P=0.002
 Analgesics15 (12.3)21 (21.3)χ2=1.6; P=0.201
 Laxative51 (41.8)61 (53.0)χ2=3.0; P=0.083
 Antithyroid2 (1.6)4 (3.5)χ2=0.8; P=0.368

HC: healthy controls; ASD: autism spectrum disorder; ID: intellectual disability; N: number of cases; BMI: body mass index; c/d: cigarettes smoked per day; CARS: Childhood Autism Rating Scale.

Significant differences are indicated in boldface.

Table 4.

Characteristics of the sample

HCN = 112ASDN = 122IDN = 115Compared groupsStatistics
Gender
 Male, N (%)50 (44.6)80 (65.5)74 (64.3)HC vs. ASDχ2=10.4; P=0.001
HC vs. IDχ2=13.0; P=0.003
 Female, N (%)62 (55.4)42 (34.4)41 (35.7)
ASD vs. IDχ2=0.1; P=0.843
Age in years (mean, SD)42.4 (11.4)40.7 (8.3)52.1 (12.4)HC vs. ASDt= 1.3; P=0.194
HC vs. IDt = −6.1; P<0.001
ASD vs. IDt = −8.3; P<0.001
BMI in kg/m2 (mean, SD)24.1 (3.1)25.1 (4.9)27.3 (5.5)HC vs. ASDt = −1.7; P=0.086
HC vs. IDt = −5.4; P<0.001
ASD vs. IDt = −3.3; p<0.001
Tobacco consumption in c/d, (mean, SD)3.2 (6.5)0.2 (1.2)0.8 (3.4)HC vs. ASDt = 4.9; P<0.001
HC vs. IDt = 3.4; P=0.001
ASD vs. IDt = −2.1; P=0.041
Chromosomal alterations, N (%)37 (30.3)10 (8.7)ASD vs. IDχ2=17.7; P<0.001
ID Severity
 Severe, N (%)16 (13.1)69 (60%)ASD vs. IDχ2=55.3; P<0.001
 Profound, N (%)106 (86.9)46 (40%)
CARS score (mean, SD)42.6 (5.9)21.8 (4.4)ASD vs. IDt=28.9; P<0.001
Treatment, N (%)
 Antipsychotic66 (54.1)64 (55.6)ASD vs. IDχ2=0.1; P=0.810
 Anticholinergic16 (13.1)17 (14.8)χ2=0.1; P=0.711
 Benzodiazepines72 (59.0)61 (53.0)χ2=0.9; P=0.354
 Antidepressant11 (9.0)18 (15.7)χ2=2.4; P=0.119
 Antiepileptic50 (41.0)33 (28.7)χ2=3.9; P=0.047
 Mood stabilizers57 (46.7)49 (42.6)χ2=0.4; P=0.525
 Antihypertensive2 (1.6)13 (12.7)χ2=9.3; P=0.002
 Analgesics15 (12.3)21 (21.3)χ2=1.6; P=0.201
 Laxative51 (41.8)61 (53.0)χ2=3.0; P=0.083
 Antithyroid2 (1.6)4 (3.5)χ2=0.8; P=0.368
HCN = 112ASDN = 122IDN = 115Compared groupsStatistics
Gender
 Male, N (%)50 (44.6)80 (65.5)74 (64.3)HC vs. ASDχ2=10.4; P=0.001
HC vs. IDχ2=13.0; P=0.003
 Female, N (%)62 (55.4)42 (34.4)41 (35.7)
ASD vs. IDχ2=0.1; P=0.843
Age in years (mean, SD)42.4 (11.4)40.7 (8.3)52.1 (12.4)HC vs. ASDt= 1.3; P=0.194
HC vs. IDt = −6.1; P<0.001
ASD vs. IDt = −8.3; P<0.001
BMI in kg/m2 (mean, SD)24.1 (3.1)25.1 (4.9)27.3 (5.5)HC vs. ASDt = −1.7; P=0.086
HC vs. IDt = −5.4; P<0.001
ASD vs. IDt = −3.3; p<0.001
Tobacco consumption in c/d, (mean, SD)3.2 (6.5)0.2 (1.2)0.8 (3.4)HC vs. ASDt = 4.9; P<0.001
HC vs. IDt = 3.4; P=0.001
ASD vs. IDt = −2.1; P=0.041
Chromosomal alterations, N (%)37 (30.3)10 (8.7)ASD vs. IDχ2=17.7; P<0.001
ID Severity
 Severe, N (%)16 (13.1)69 (60%)ASD vs. IDχ2=55.3; P<0.001
 Profound, N (%)106 (86.9)46 (40%)
CARS score (mean, SD)42.6 (5.9)21.8 (4.4)ASD vs. IDt=28.9; P<0.001
Treatment, N (%)
 Antipsychotic66 (54.1)64 (55.6)ASD vs. IDχ2=0.1; P=0.810
 Anticholinergic16 (13.1)17 (14.8)χ2=0.1; P=0.711
 Benzodiazepines72 (59.0)61 (53.0)χ2=0.9; P=0.354
 Antidepressant11 (9.0)18 (15.7)χ2=2.4; P=0.119
 Antiepileptic50 (41.0)33 (28.7)χ2=3.9; P=0.047
 Mood stabilizers57 (46.7)49 (42.6)χ2=0.4; P=0.525
 Antihypertensive2 (1.6)13 (12.7)χ2=9.3; P=0.002
 Analgesics15 (12.3)21 (21.3)χ2=1.6; P=0.201
 Laxative51 (41.8)61 (53.0)χ2=3.0; P=0.083
 Antithyroid2 (1.6)4 (3.5)χ2=0.8; P=0.368

HC: healthy controls; ASD: autism spectrum disorder; ID: intellectual disability; N: number of cases; BMI: body mass index; c/d: cigarettes smoked per day; CARS: Childhood Autism Rating Scale.

Significant differences are indicated in boldface.

Clinical data

ASD and ID

A trained psychiatrist reviewed the clinical record of each participant to collect sociodemographic, anthropometric, medication use, family history and neurodevelopment features. In conjunction with the referring physician of each participant, the psychiatrist completed a questionnaire focused on CAMDs, which are divided into the following categories: migraine headaches, peripheral neurovascular disorders, gastrointestinal dysmotility, neurological disorders, cardiac abnormalities, skeletal muscle disorders, endocrine disorders, and constitutional disorders. The same questionnaire has been used in previous studies to identify maternal inheritance in disorders with presumed mitochondrial dysfunction (42–44) and the presence of clinical features associated with mitochondrial disorders among relatives of patients with schizophrenia (45). Some conditions present in the questionnaire, such as questions regarding headache, migraine or fatigue, could not be properly assessed because of the severity of the ID of the subjects.

HCs

The same psychiatrist who obtained data from subjects in the ASD and ID groups collected sociodemographic, anthropometric and family history data, and completed the CAMDs questionnaire for HCs. Participants who had no personal, first- or second-degree relatives with antecedents of severe mental disorders were selected for inclusion in the control group.

DNA extraction

Total DNA was obtained from peripheral blood mononuclear cells from fasting subjects using the Gentra® PureGene reagent (Qiagen, Barcelona, Spain), according to the manufacturer’s instructions. DNA concentrations were quantified using a NanoDrop spectrophotometer (Thermo Fisher Scientific, Madrid, Spain).

mtDNA-targeted next generation sequencing

The entire mtDNA was amplified using three pairs of mtDNA-specific primers that produced three overlapping PCR fragments; A: 6928 base pairs (bp), positions 569–7497; B: 7050 bp, positions 5061–12 111; and C: 6867 bp, position 11 107–1405. Long-range PCR was performed using 10 ng of DNA with Expand Long Range dNTPack (Roche, Barcelona, Spain), and the PCR products were purified using the QIAquick PCR Purification kit (Qiagen, Barcelona, Spain), according to the manufacturer's instructions. Nineteen blood samples from the ASD group and 15 blood samples from the ID group were discarded, as they did not amplify any of the mtDNA fragments. Finally, the three PCR fragments from each individual were mixed in equimolar ratios, and each sample was prepared for sequencing in an Ion Torrent Personal Genome Machine (PGM, Fisher Scientific, Madrid, Spain), according to the manufacturer’s user guide. Seven multiplexed sequencing runs (32 samples per run) were performed on Ion 316 chips (Fisher Scientific, Madrid, Spain). After filtering, sequences from five subjects with ASD and five subjects with ID were discarded because of low-quality reads. We obtained an average of 60 000 high-quality reads from 98 subjects with ASD and 95 subjects with ID that were analysed using MToolBox (21) that assembles mtDNA starting from NGS data, assigns haplogroup and prioritizes mtDNA variants of clinical interest according to Santorsola et al. (37) for non synonymous variants and Diroma et al. (46) for rRNA and tRNA variants. Briefly, using 14 144 complete sequences from healthy individuals available in HmtDB (47), the process reports rare or private variants, the nucleotide variability, a disease score based on several predictors of pathogenicity for non-synonymous variants and an RNA pathogenicity prediction for rRNA and tRNA variants based on published data (37,46). Finally, information about the variants associated with clinical outcomes present in MITOMAP (www.mitomap.org), ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/), OMIM (https://omim.org) and dbSNP (https://www.ncbi.nlm.nih.gov/SNP/) was obtained through MToolBox, together with data on homoplasmy and heteroplasmy (homoplasmy refers to a cell that has a uniform collection of mtDNA—either completely normal mtDNA or completely mutant mtDNA—and heteroplasmy refers to a cell with different proportions of normal and mutant mtDNA). Regarding heteroplasmy, we set the cut-off to 0.1, and therefore, variants with low levels of heteroplasmy (< 10%) were not reported. The mean coverage obtained was 512. Variants are reported according to the revised Cambridge Reference Sequence (rCRS) of the human mtDNA, GenBank sequence NC_012920 gi: 251831106.

Selection criteria for the pathogenicity prediction

Non-synonymous variants

We selected putative pathogenic non-synonymous variants that fulfilled the two cut-off prioritization criteria described by Santorsola et al., (37) diseases score > 0.4311 and nucleotide variability < 0.0026. The disease score is based on the weighted mean of the probabilities that an amino acid substitution may affect gene/protein function provided by the pathogenicity predictors SIFT, Polyphen-2, MutPred, SNPs&GO, PhD-SNP and PANTHER.

rRNA and tRNA variants

The RNA prediction score retrieved by MToolBox was based on published data, and the pathogenic range was established between 0.35 and 1 (46). Thus, variants that have not previously been associated with diseases have no RNA prediction score. Based on our results, we selected putative pathogenic rRNA or tRNA variants presenting an RNA prediction score ≥ 0.35 or a nucleotide variability ≤ 0.0097.

Heteroplasmy level

For most of the described pathogenic mutations, the disease manifests when the heteroplasmic threshold ranges from 0.60 to 0.95, depending on the mutation and cell type (31). Therefore, we only considered variants that presented heteroplasmy levels higher than 60% as putatively pathogenic.

Quantification of the mtDNA copy number

We measured the mtDNA content in total DNA using quantitative real-time PCR (qPCR). Two mtDNA regions were selected as target regions: MT-ND1 and MT-ND4 genes. Nuclear RPPH1 (corresponding to the RNase P enzyme), a single-copy nuclear gene, was selected as the reference gene. The mtDNA content was calculated for each sample and region using the 2-ΔCq method (ΔCq = Cq of the mitochondrial gene—Cq of the reference gene). The mtDNA content was calculated as the ratio between the number of copies of the mitochondrial genome and the number of copies of the nuclear genome (mtDNA/nDNA). MT-ND1 is a rarely deleted region and MT-ND4 is located in the major arc of mtDNA and deleted in 97% of all common deletions (17). Therefore, the MT-ND1/MT-ND4 ratio indicates the deletion ratio. Details regarding qPCR reactions and quality control have been published previously (17,48).

Statistical analyses

The normality of the distributions of continuous variables was explored using the Kolmogorov-Smirnov normality test. Chi-square tests were used to compare the presence of CAMDs among the study groups. Spearman’s correlation analysis was used to explore the correlation between the mtDNA content of the two studied mtDNA regions and age, body mass index (BMI) and smoking habits. A general linear model was used to identify whether the mtDNA content in the three study groups differed after adjusting for age, BMI and number of cigarettes smoked per day as covariates. Data were processed using IBM SPSS Statistics for Windows, Version 23.0 (IBM Corp., Armonk, NY) and Prism, Version 5 (Graphpad, La Jolla, CA).

Supplementary Material

Supplementary Material is available at HMG online.

Acknowledgements

We are grateful to the subjects who participated in the study, and we also acknowledge the technicians from the Biobanc-IISPV in Reus (http://www.iispv.cat) for managing the samples and the physicians who helped collect the data related to the CAMDs conditions.

Conflict of Interest statement. None declared.

Funding

Instituto de Salud Carlos III of the Spanish Ministry of Science and Innovation in Spain [grant numbers PS09/01052 and PI12/01885 to L.M.] and FEDER. H.T. was the recipient of a FI-DGR, and G.M. was the recipient of a BP-DGR scholarship, both of which were from the Generalitat de Catalunya.

References

1

American Psychiatric Association
(
2013
) Diagnostic and Statistical Manual of Mental Disorders.

2

Lai
M.-C.
,
Lombardo
M.V.
,
Baron-Cohen
S.
(
2014
)
Autism
.
Lancet (London, England)
,
383
,
896
910
.

3

Srour
M.
,
Shevell
M.
(
2014
)
Genetics and the investigation of developmental delay/intellectual disability
.
Arch. Dis. Child
.,
99
,
386
389
.

4

de la Torre-Ubieta
L.
,
Won
H.
,
Stein
J.L.
,
Geschwind
D.H.
(
2016
)
Advancing the understanding of autism disease mechanisms through genetics
.
Nat. Med
.,
22
,
345
361
.

5

La Malfa
G.
,
Lassi
S.
,
Bertelli
M.
,
Salvini
R.
,
Placidi
G.F.
(
2004
)
Autism and intellectual disability: a study of prevalence on a sample of the Italian population
.
J. Intellect. Disabil. Res
.,
48
,
262
267
.

6

Levy
S.E.
,
Giarelli
E.
,
Lee
L.-C.
,
Schieve
L.A.
,
Kirby
R.S.
,
Cunniff
C.
,
Nicholas
J.
,
Reaven
J.
,
Rice
C.E.
(
2010
)
Autism spectrum disorder and co-occurring developmental, psychiatric, and medical conditions among children in multiple populations of the United States
.
J. Dev. Behav. Pediatr
.,
31
,
267
275
.

7

Croen
L.A.
,
Zerbo
O.
,
Qian
Y.
,
Massolo
M.L.
,
Rich
S.
,
Sidney
S.
,
Kripke
C.
(
2015
)
The health status of adults on the autism spectrum
.
Autism
,
19
,
814
823
.

8

Geschwind
D.H.
(
2009
)
Advances in autism
.
Ann. Rev. Med
.,
60
,
367
380
.

9

Nicolaidis
C.
,
Kripke
C.C.
,
Raymaker
D.
(
2014
)
Primary care for adults on the autism spectrum
.
Med. Clin. North Am
.,
98
,
1169
1191
.

10

Vafai
S.B.
,
Mootha
V.K.
(
2012
)
Mitochondrial disorders as windows into an ancient organelle
.
Nature
,
491
,
374
383
.

11

El-Hattab
A.W.
,
Scaglia
F.
(
2016
)
Mitochondrial cytopathies
.
Cell Calcium
,
60
,
199
206
.

12

Frye
R.E.
,
Rossignol
D.A.
(
2011
)
Mitochondrial dysfunction can connect the diverse medical symptoms associated with autism spectrum disorders
.
Pediatr. Res
.,
69
,
41R
7R
.

13

Rossignol
D.A.
,
Frye
R.E.
(
2012
)
Mitochondrial dysfunction in autism spectrum disorders: a systematic review and meta-analysis
.
Mol. Psychiatry
,
17
,
290
314
.

14

Valenti
D.
,
de Bari
L.
,
De Filippis
B.
,
Henrion-Caude
A.
,
Vacca
R.A.
(
2014
)
Mitochondrial dysfunction as a central actor in intellectual disability-related diseases: An overview of Down syndrome, autism, Fragile X and Rett syndrome
.
Neurosci. Biobehav. Rev
.,
46
,
202
217
.

15

Goh
S.
,
Dong
Z.
,
Zhang
Y.
,
DiMauro
S.
,
Peterson
B.S.
(
2014
)
Mitochondrial dysfunction as a neurobiological subtype of autism spectrum disorder: evidence from brain imaging
.
JAMA Psychiatry
,
71
,
665
671
.

16

Verge
B.
,
Alonso
Y.
,
Valero
J.
,
Miralles
C.
,
Vilella
E.
,
Martorell
L.
(
2011
)
Mitochondrial DNA (mtDNA) and schizophrenia
.
Eur. Psychiatry
,
26
,
45
56
.

17

Torrell
H.
,
Montaña
E.
,
Abasolo
N.
,
Roig
B.
,
Gaviria
A.M.
,
Vilella
E.
,
Martorell
L.
(
2013
)
Mitochondrial DNA (mtDNA) in brain samples from patients with major psychiatric disorders: gene expression profiles, mtDNA content and presence of the mtDNA common deletion
.
Am. J. Med. Genet. B, Neuropsychiatr. Genet
.,
162
,
213
223
.

18

Meng
S.
,
Wu
S.
,
Liang
L.
,
Liang
G.
,
Giovannucci
E.
,
De Vivo
I.
,
Nan
H.
(
2016
)
Leukocyte mitochondrial DNA copy number, anthropometric indices, and weight change in US women
.
Oncotarget
,
7
,
60676
60686
.

19

Malik
A.N.
,
Czajka
A.
(
2013
)
Is mitochondrial DNA content a potential biomarker of mitochondrial dysfunction?
Mitochondrion
,
13
,
481
492
.

20

Golzio
C.
,
Katsanis
N.
(
2016
)
Mitochondrial Copy Number as a Biomarker for Autism?
Pediatrics
,
137
,
e20160049
.

21

Calabrese
C.
,
Simone
D.
,
Diroma
M.A.
,
Santorsola
M.
,
Guttà
C.
,
Gasparre
G.
,
Picardi
E.
,
Pesole
G.
,
Attimonelli
M.
(
2014
)
MToolBox: a highly automated pipeline for heteroplasmy annotation and prioritization analysis of human mitochondrial variants in high-throughput sequencing
.
Bioinformatics
,
30
,
3115
3117
.

22

Mosquera-Miguel
A.
,
Torrell
H.
,
Abasolo
N.
,
Arrojo
M.
,
Paz
E.
,
Ramos-Ríos
R.
,
Agra
S.
,
Páramo
M.
,
Brenlla
J.
,
Martínez
S.
et al. (
2012
)
No evidence that major mtDNA European haplogroups confer risk to schizophrenia
.
Am. J. Med. Genet. B, Neuropsychiatr. Genet
.,
159B
,
414
421
.

23

Cashin
A.
,
Buckley
T.
,
Trollor
J.N.
,
Lennox
N.
(
2016
)
A scoping review of what is known of the physical health of adults with autism spectrum disorder
.
J. Intellect. Disabil
., 10.1177/1744629516665242.

24

Brown
M.
,
Jacobstein
D.
,
Yoon
I.S.
,
Anthony
B.
,
Bullock
K.
(
2016
)
Systemwide initiative documents robust health screening for adults with intellectual disability
.
Intellect. Develop. Disabil
.,
54
,
354
365
.

25

Jones
K.B.
,
Cottle
K.
,
Bakian
A.
,
Farley
M.
,
Bilder
D.
,
Coon
H.
,
McMahon
W.M.
(
2016
)
A description of medical conditions in adults with autism spectrum disorder: A follow-up of the 1980s Utah/UCLA Autism Epidemiologic Study
.
Autism
,
20
,
551
561
.

26

Chen
S.
,
Li
Z.
,
He
Y.
,
Zhang
F.
,
Li
H.
,
Liao
Y.
,
Wei
Z.
,
Wan
G.
,
Xiang
X.
,
Hu
M.
et al. (
2015
)
Elevated mitochondrial DNA copy number in peripheral blood cells is associated with childhood autism
.
BMC Psychiatry
,
15
,
50.

27

Yoo
H.J.
,
Park
M.
,
Kim
S.A.
(
2017
)
Difference in mitochondrial DNA copy number in peripheral blood cells between probands with autism spectrum disorders and their unaffected siblings
.
World J. Biol. Psychiatry 
,
18
,
151
156
.

28

Knez
J.
,
Winckelmans
E.
,
Plusquin
M.
,
Thijs
L.
,
Cauwenberghs
N.
,
Gu
Y.
,
Staessen
J.A.
,
Nawrot
T.S.
,
Kuznetsova
T.
(
2015
)
Correlates of peripheral blood mitochondrial DNA content in a general population
.
Am. J. Epidemiol
.,
183
,
kwv175
.

29

Giulivi
C.
,
Zhang
Y.
,
Omanska-Klusek
A.
,
Ross-Inta
C.
,
Wong
S.
,
Hertz-Picciotto
I.
,
Tassone
F.
,
Pessah
I.N.
(
2010
)
Mitochondrial dysfunction in autism
.
JAMA
,
304
,
2389
2396
.

30

Picard
M.
,
Zhang
J.
,
Hancock
S.
,
Derbeneva
O.
,
Golhar
R.
,
Golik
P.
,
O’Hearn
S.
,
Levy
S.
,
Potluri
P.
,
Lvova
M.
et al. (
2014
) Progressive increase in mtDNA 3243A> G heteroplasmy causes abrupt transcriptional reprogramming. Proc. Natl Acad. Sci. USA, 111, E4033–4042.

31

Smith
P.M.
,
Lightowlers
R.N.
(
2011
)
Altering the balance between healthy and mutated mitochondrial DNA
.
J. Inherit. Metabol. Dis
.,
34
,
309
313
.

32

Napoli
E.
,
Wong
S.
,
Giulivi
C.
(
2013
)
Evidence of reactive oxygen species-mediated damage to mitochondrial DNA in children with typical autism
.
Mol. Autism
,
4
,
2.

33

Avdjieva-Tzavella
D.
,
Mihailova
S.
,
Lukanov
C.
,
Naumova
E.
,
Simeonov
E.
,
Tincheva
R.
,
Toncheva
D.
(
2012
)
Mitochondrial DNA mutations in two bulgarian children with autistic spectrum disorders
.
Balkan J. Med. Genet
.,
15
,
47
54
.

34

Pons
R.
,
Andreu
A.L.
,
Checcarelli
N.
,
Vilà
M.R.
,
Engelstad
K.
,
Sue
C.M.
,
Shungu
D.
,
Haggerty
R.
,
De Vivo
D.C.
,
DiMauro
S.
(
2004
)
Mitochondrial DNA abnormalities and autistic spectrum disorders
.
J. Pediatr
.,
144
,
81
85
.

35

Patowary
A.
,
Nesbitt
R.
,
Archer
M.
,
Bernier
R.
,
Brkanac
Z.
(
2017
)
Next generation sequencing mitochondrial DNA analysis in autism spectrum disorder
.
Autism Res
.,
10
,
1338
1343
.

36

Wang
Y.
,
Picard
M.
,
Gu
Z.
,
Girirajan
S.
(
2016
)
Genetic evidence for elevated pathogenicity of mitochondrial DNA heteroplasmy in autism spectrum disorder
.
PLoS Genet
.,
12
,
e1006391.

37

Santorsola
M.
,
Calabrese
C.
,
Girolimetti
G.
,
Diroma
M.A.
,
Gasparre
G.
,
Attimonelli
M.
(
2016
)
A multi-parametric workflow for the prioritization of mitochondrial DNA variants of clinical interest
.
Hum. Genet
.,
135
,
121
136
.

38

Elliott
H.R.
,
Samuels
D.C.
,
Eden
J.A.
,
Relton
C.L.
,
Chinnery
P.F.
(
2008
)
Pathogenic mitochondrial DNA mutations are common in the general population
.
Am. J. Hum. Genet
.,
83
,
254
260
.

39

Ye
K.
,
Lu
J.
,
Ma
F.
,
Keinan
A.
,
Gu
Z.
(
2014
)
Extensive pathogenicity of mitochondrial heteroplasmy in healthy human individuals
.
Proc. Natl Acad. Sci. USA
,
111
,
10654
10659
.

40

Lord
C.
,
Rutter
M.
,
Le Couteur
A.
(
1994
)
Autism diagnostic interview-revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders
.
J. Autism Dev. Disord
.,
24
,
659
685
.

41

Schopler
E.
,
Reichler
R.J.
,
DeVellis
R.F.
,
Daly
K.
(
1980
)
Toward objective classification of childhood autism: Childhood Autism Rating Scale (CARS)
.
J. Autism Develop. Disord
.,
10
,
91
103
.

42

Boles
R.G.
,
Adams
K.
,
Li
B.U.K.
(
2005
)
Maternal inheritance in cyclic vomiting syndrome
.
Am. J. Med. Genet. Part A
,
133A
,
71
77
.

43

Boles
R.G.
,
Burnett
B.B.
,
Gleditsch
K.
,
Wong
S.
,
Guedalia
A.
,
Kaariainen
A.
,
Eloed
J.
,
Stern
A.
,
Brumm
V.
(
2005
)
A high predisposition to depression and anxiety in mothers and other matrilineal relatives of children with presumed maternally inherited mitochondrial disorders
.
Am. J. Med. Genet. Part B, Neuropsychiatr. Genet
.,
137B
,
20
24
.

44

Burnett
B.B.
,
Gardner
A.
,
Boles
R.G.
(
2005
)
Mitochondrial inheritance in depression, dysmotility and migraine?
J. Affect. Disord
.,
88
,
109
116
.

45

Verge
B.
,
Alonso
Y.
,
Miralles
C.
,
Valero
J.
,
Vilella
E.
,
Boles
R.G.
,
Martorell
L.
(
2012
)
New evidence for the involvement of mitochondrial inheritance in schizophrenia: results from a cross-sectional study evaluating the risk of illness in relatives of schizophrenia patients
.
J. Clin. Psychiatr
.,
73
,
684
690
.

46

Diroma
M.A.
,
Lubisco
P.
,
Attimonelli
M.
(
2016
)
A comprehensive collection of annotations to interpret sequence variation in human mitochondrial transfer RNAs
.
BMC Bioinformatics
,
17
,
338.

47

Rubino
F.
,
Piredda
R.
,
Calabrese
F.M.
,
Simone
D.
,
Lang
M.
,
Calabrese
C.
,
Petruzzella
V.
,
Tommaseo-Ponzetta
M.
,
Gasparre
G.
,
Attimonelli
M.
(
2012
)
HmtDB, a genomic resource for mitochondrion-based human variability studies
.
Nucleic Acids Res
.,
40
,
D1150
D1159
.

48

Torrell
H.
,
Alonso
Y.
,
Garrabou
G.
,
Mulet
D.
,
Catalán
M.
,
Valiente-Pallejà
A.
,
Carreño-Gago
L.
,
García-Arumí
E.
,
Montaña
E.
,
Vilella
E.
,
Martorell
L.
(
2017
)
Mitochondrial dysfunction in a family with psychosis and chronic fatigue syndrome
.
Mitochondrion
,
34
,
1
8
.

Author notes

Alba Valiente-Pallejà, Helena Torrell and Gerard Muntané have equally contributed to the work.

Supplementary data