Abstract
Mitochondrial dysfunction is a neurobiological phenomenon implicated in the pathophysiology of schizophrenia and bipolar disorder that can synergistically affect synaptic neurotransmission. We hypothesized that schizophrenia and bipolar disorder share molecular alterations at the mitochondrial and synaptic level: mtDNA copy number (CN), mtDNA common deletion (CD), mtDNA total deletion, complex I activity, synapse number, and synaptic mitochondria number. These mitochondrial parameters were studied in dorsolateral prefrontal cortex (DLPFC), superior temporal gyrus (STG), primary visual cortex (V1), and nucleus accumbens (NAc) of 44 controls (CON), 27 schizophrenia (SZ), and 26 bipolar disorder (BD) subjects. The results showed- (i) the mtDNA CN is significantly higher in DLPFC of both SZ and BD, decreased in the STG of BD, and unaltered in V1 and NAc of both SZ and BD; (ii) the mtDNA CD is significantly higher in DLPFC of BD while unaltered in STG, V1 and NAc of both SZ and BD; (iii) The total deletion burden is significantly higher in DLPFC in both SZ and BD while unaltered in STG, V1, and NAc of SZ and BD; (iv) Complex I activity is significantly lower in DLPFC of both SZ and BD with no alteration in STG, V1 and NAc. (v) The number of synapses is decreased in DLPFC of both SZ and BD, while the synaptic mitochondria number was significantly lower in female SZ and female BD compared to female controls. Overall, these findings will pave the way to better understand the pathophysiology of schizophrenia and bipolar disorder for therapeutic interventions.
Introduction
Mitochondria are the primary energy-producers of the cells as they generate ∼90% of cell energy in the form of adenosine triphosphate (ATP). Energy production of mitochondria in the brain is critical, as the brain uses over 20% of the energy produced in the body, yet it is only about 2% of adult body weight. Mitochondria generate ATP through a complex set of enzymes called electron transport chain (ETC). The ETC consists of five protein complexes (Complex I, II, III, IV and V). The ETC proteins are encoded by two distinct genomes: the mitochondrial genome (mtDNA) and the nuclear genome (nDNA). mtDNA are closed circle, intron-free, and polyploid double-stranded DNA molecules consisting of 16569 nucleotides that encode for 13 proteins, 22 transfer RNAs and 2 ribosomal RNAs. Conversely, nDNA encodes for approximately 1500 genes involved in mitochondrial functions and localization [1]. Within ETC, Complex I is the largest multimeric protein complex and consists of 45 subunits, from which 6 subunits are encoded by mtDNA.
Schizophrenia (SZ) and bipolar disorder (BD) are complex polygenic neuropsychiatric disorders. One pathway implicated in the pathophysiology of SZ and BD is mitochondria dysfunction [1–3]. Clinical studies have shown that patients with SZ and BD have abnormal energy metabolites and redox states implicating dysfunctional mitochondria [4–6]. Genome-wide association studies (GWASs) and gene-set based analysis have linked the pathophysiology of SZ and BD to synaptic dysfunction and nuclear-encoded mitochondria genes (NEMG) [7]. Mitochondrial dysfunction in SZ and BD appears to stem from a complex set of causalities which include but are not limited to -a) single nucleotide polymorphism (SNP) in mtDNA and NEMGs [8]; b) decreased expression of genes responsible for mitochondrial function revealed by transcriptomics analysis [9, 10]; c) altered mtDNA deletion burdens in SZ and BD brain tissues [11, 12]; d) altered mtDNA copy number [13]; e) altered enzymatic activity of ETC proteins [14]; f) an altered number of mitochondria puncta and mitochondria volume density [1]. Based on above mentioned etiology, it is well established that genetic risk component plays a major role in dysfunctional mitochondria of SZ and BD. Induced pluripotent stem cells (iPSCs) derived neurons and cerebral organoids from patients with SZ and BD showed dysfunctional mitochondria which further strengthen the role of genetic risk component in pathophysiology of SZ and BD [15, 16].
Due to the complex nature of SZ and BD, we conducted a systematic investigation of critical mitochondria parameters related to mtDNA, as well as mitochondrial distribution and function, in SZ and BD brains. Specifically, we investigated the following mitochondrial parameters in clinically well-defined SZ and BD subjects: mtDNA copy number, mtDNA common deletions, mtDNA total deletions, Complex I activity, total mitochondria number, total synapse number and synaptic mitochondria number. Models tested for diagnostic effect included potential factors such as age, sex and antipsychotic medication(s). To achieve this goal, we used postmortem brain tissues from 44 control (CON), 27 schizophrenia (SZ) and 26 bipolar disorder (BD) subjects. Four key brain regions across this cohort have been investigated to identify the regional effect and specificity of any observed mitochondrial difference: dorsolateral prefrontal cortex (DLPFC), superior temporal gyrus (STG), primary visual cortex (V1) and nucleus accumbens (NAc). First, we determined how mtDNA copy number vary in SZ and BD compared to CON. Second, we measured mtDNA common deletion in SZ and BD compared to CON. Third, we determined the pooled effect of total deletions through our high-throughput and high-resolution pipeline called “Splice-Break” in CON, SZ and BD. Fourth, we determined the alteration of Complex I activity in SZ and BD. Finally, we calculated the overall mitochondria number, synapse number and synaptic mitochondria number through confocal microscopy and IMARIS modeling in CON, SZ and BD.
Materials and Methods
Subjects
Postmortem brain collection and preservation were conducted as previously published procedures [17]. Briefly, postmortem brain tissues were obtained from the UCI-Pritzker brain bank. Written consent was obtained from the next-of-kin for each subject. The postmortem brain collection and experimental procedures were approved by UCI Institutional Review Board (IRB). Psychological autopsies were conducted based on family informant interviews, subjects’ medication history, medical and psychiatric records, and coroners’ toxicology reports. Based on psychological autopsies, three clinically well-defined experimental groups were studied: Control (CON), Schizophrenia (SZ) and Bipolar Disorder (BD). In total, 44 CON, 27 SZ and 26 BD subjects were included in this study. Results for some of these parameters were previously published for DLPFC (parameters 1A, 2A, and 4A) with the same subjects [13], and are included for comparisons to three additional brain regions. The full demographics of all the subjects has been depicted in Table ST1.
mtDNA Copy Number and Common Deletion
mtDNA copy number (mtDNA CN), mtDNA common deletion (mtDNA CD; adjusted breakpoints:8482-13460) and mtDNA amplicon resequencing was performed on extracted genomic DNA (gDNA) of samples as described in our previously published procedures [12, 13]. The in-detail procedure has been described in Supplementary Information.
Pooled Effect of Total Deletion of mtDNA by High-resolution Pipeline “Splice-Break”
With our recently developed pipeline [18], called Splice-Break, we measured total deletions of mitochondria DNA (common deletion + other deletions) which occur at different rates in different brain regions and increases with age. The full description of the “Splice-Break” method has been explained in Supplementary Information. For more details about the “Splice-Break”, please read [18].
Complex I Activity (activity/unit of Complex I protein)
The Complex I activity was measured based on our previously published two-steps procedure [13]. For more of this procedure, please read Supplementary Information.
Immunohistochemistry (IHC) for Colocalization of Mitochondria and Synaptic Marker
The gene TOMM40 (translocase of outer mitochondrial membrane 40) codes for the protein TOM40 which is required for the import of proteins into the mitochondria. The PSD95 protein
(postsynaptic density protein 95) is exclusively localized in excitatory synapses. The colocalization of TOM40 and PSD95 was evaluated by immunohistochemistry (IHC) as described previously with slight modification [19]. Please read Supplementary Information for more details.
Statistical Analyses
The raw data for each dependent variable was imported into JMP Pro 16.0 (SAS Institute, Inc., NC, USA) and checked for outliers greater than 3 standard deviations from the mean, which were removed. Additional variables beyond diagnosis and brain region (age, PMI, pH, sex) were checked for correlations with raw data and were included in models when significant. Linear regression model for single variables were used to determine group differences in a single brain region. When similar technology was used across three cortical regions in which all subjects were assayed, a repeated measure (mixed effects linear model) based on subject, region, and the diagnosis was calculated to increase the analysis robustness and compared with single region analysis. Post hoc analysis of group effects were conducted after main effects or interaction effects were determined significant (p < 0.05). Exploratory analysis of medication effects using previous grouping of subjects with neuroactive compounds in cerebellar toxicology reports from NMS Labs (PA, USA) was used as a grouping variable in pooled BD and SZ subjects, details of these subjects’ toxicology were previously published [13].
Results
mtDNA Copy Number in SZ and BD
Different tissues have various bioenergetics profiles which makes mtDNA copy number (CN) vary in a tissue-specific (and sometimes age-dependent) manner. Decreased peripheral mtDNA CN has been reported in SZ and BD [20, 21]. Peripheral mtDNA CN has also been shown to decrease with aging, although brain tissue mtDNA CN showed no age-related changes [22, 23]. Here, we sought to determine how SZ and BD affect the mtDNA CN in four critical brain regions of neocortex. A mixed model effects analysis, with adjustments of covariates (∼age+sex+region+subject), was conducted for mtDNA CN in the DLPFC and STG regions as a repeated analysis. The V1 and NAc regions were analyzed separately because V1 data were obtained through ddPCR and NAc regions had subset of samples. For the DLPFC, both SZ and BD groups displayed significant increases (p= 0.003 and p = 0.011, respectively) in the level of mtDNA copy number compared to CON (Figure 1A). The same subjects’ STG region showed a significant decrease (p = 0.036) in mtDNA CN in BD compared to controls and a non-significant decrease in mtDNA CN in SZ compared to controls (Figure 1B-D) (Table ST2). Although the interaction of brain region (DLPFC, STG) and diagnosis was significant (p = 0.0015), but the mtDNA copy number brain region correlations (DLPFC vs STG) were not significant (Figure SF1). In separate analyses of V1 cortex and nucleus accumbens regions, there were no significant differences between groups for mtDNA CN in the V1 (p = 0.72) or the NAc (p = 0.30). Next, we determined the effect of age on mtDNA CN in DLPFC, STG and V1. Age did not have any significant effect on mtDNA CN in DLPFC, STG and V1 (Table ST7). Another potential confounding factor affecting mtDNA CN in SZ and BD is the use of antipsychotic medications. Antipsychotic medications have been shown to decrease peripheral mtDNA CN [24]. Thus, we also tested whether antipsychotic medications had any effect on mtDNA CN or not in four brain regions included in this study. Interestingly, we did not find any medication effect on mtDNA CN in DLPFC, STG, V1 and NAc of SZ and BD patients (Table ST8), suggesting these brain regions may have more stable mitochondrial CN than blood when exposed to environmental factors or medication.
mtDNA copy number in postmortem DLPFC, STG, V1 and NAc of control (CON), schizophrenia (SZ) and bipolar disorders (BD) groups. mtDNA copy number (CN) for DLPFC, STG and NAc were obtained through qPCR while V1 mtDNA copy number (CN) was obtained through ddPCR. mtDNA CN of a sample was calculated as mtDNA CN/ALB CN. Mixed model analysis was performed for statistical significance with adjustment of age and sex. DLPFC and STG were analyzed together to obtain the regional correction whereas V1 and NAc data were analyzed individually. DLPFC mtDNA CN was significantly higher in both SZ and BD compared to control group [A]. In STG, BD group had significantly lower mtDNA CN compared to control [B]. In V1 and NAc, there were no significant groups differences [C, D]. All data are represented as least square mean (LSM) ± Std. error. p<0.05 and p<0.01 are denoted by * and **, respectively, otherwise p>0.05.
Common Deletion of mtDNA in SZ and BD
The mtDNA common deletion (CD) is a large (4977 bp) deletion of mtDNA, shown by us and others to increase with age in multiple brain regions [12, 25]. The mtDNA CD is an indicator of oxidative stress due to ROS and the loss of 12 essential mitochondrial genes within the deleted region. mtDNA is near the electron transport chain, a major source of reactive oxygen species that can damage mtDNA and trigger deletions. Mixed model analysis was performed for statistical significance with adjustment of age and sex with the DLPFC, STG and V1 analyzed together whereas NAc data were analyzed individually. The regional mtDNA CD was increased in BD compared to CON (p =0.014) and increased in BD compared to SZ (p =0.046) in the DLPFC (Figure 2A). The effect was mainly due to increased CD in DLPFC for BD. The %CD in STG, V1 and NAc was not significantly altered (Figure 2B-D) (Table ST3). In contrast to the mtDNA CN results, the mtDNA CD showed highly significant correlations in the levels between brain regions. Regional analyses of the CD data assayed by qPCR showed a significant age effect (%CD increases with age) in each of the three cortical regions (Figure SF3). The positive correlation of age and the CD was significant in V1 cortex (p =0.032, r = 0.22), DLPFC (p = 0.0008, r = 0.34), and STG (p =0.008, r = 0.27) (Table ST9). The positive correlations of the mtDNA CD were highly significant across regions by subject, e.g., a high inter-regional correlation (STG-DLPFC, r = 0.75, 3.25E-18), which suggests a possible global mechanism for the common deletion dependent upon age (Figure SF2). A subset of subjects was assayed for the mtDNA CD in the nucleus accumbens and the inter-regional correlations remained significant when analyzing only this subset of subjects across 4 brain regions with common subjects assayed in all four regions. Further, we addressed if correlations of the mtDNA CD were altered by antipsychotic drugs and medications, for this, patients were grouped by detectable or not detectable toxicological psychotropic drugs and medications based on postmortem toxicology analysis. The inter-regional correlations remained significant for both groups regardless of toxicology measures. Thus, presence of medications or drugs did not alter the significant inter regional correlations of mtDNA CD levels, suggesting it either had no effect, or had an equivalent and non-significant effect on all brain regions we evaluated (data not shown).
mtDNA common deletions in postmortem DLPFC, STG, V1 and NAc of control (CON), schizophrenia (SZ) and bipolar disorders (BD) groups. mtDNA common deletion (CD) for DLPFC, STG and NAc were obtained through qPCR while V1 mtDNA copy number was obtained through ddPCR. %mtDNA CD was calculated as [mtDNA CD copies/ (wild type mtDNA copies+ mtDNA CD copies)] x 100. Mixed model analysis was performed for statistical significance with adjustment of age and sex. DLPFC, STG and V1 were analyzed together to obtain the regional correction whereas NAc data were analyzed individually. DLPFC mtDNA CD was significantly higher in BD compared to control group [A]. There were no significant groups differences in STG, V1 and NAc mtDNA CD [B, C, D]. All data are represented as least square mean (LSM) ± Std. error. p<o.o5 and p<0.01 are denoted by * and **, respectively, otherwise p>0.05.
The pooled effect of total mtDNA deletions in SZ and BD
We utilized the highly sensitive “Splice-Break” pipeline to detect the total number of unique mtDNA deletions whose breakpoints fell within the position of mtDNA 357-15925. The “Splice- Break” pipeline provided us deletion read %’s for hundreds to thousands of mtDNA deletion breakpoints, and we calculated the cumulative deletion metrics as previously described [18]. Here, we focused on the cumulative analysis of “deletions per 10K coverage”, which is the estimation of the pooled effect of total deletions species (i.e., unique set of breakpoints), as the other two metrics described in our method paper were not significantly different in group comparisons. Previously, we showed that deletions per 10K coverage parameter did not significantly vary with age, brain pH and postmortem interval (PMI) [18]. Here, we computed standardized residuals for deletions per 10K coverage with benchmark coverage and age added as a covariate. The repeated measure analysis revealed that the standardized residuals of deletions per 10K was significantly higher in DLPFC of both SZ (p=0.036) and BD (p=0.000055) compared to CON (Figure 3A). There were no other significant group comparisons within STG [CON vs SZ, p=0.78; CON vs BD, p=0.68], V1 [CON vs SZ, p=0.87; CON vs BD, p=0.71] and NAc [CON vs SZ, p=0.77; CON vs BD, p=0.40] (Figure 3B-D) (Table ST4). The deletions per 10K coverage was not significantly correlated with age (Table ST10). The medications did not affect the total deletions per 10K coverage in DLPFC, STG and V1 (Figure SF4).
Standardized residuals of deletions per 10K in postmortem DLPFC, STG, V1 and NAc of control (CON), schizophrenia (SZ) and bipolar disorders (BD) groups. Mixed model analysis was performed for statistical significance with adjustment of age and sex. DLPFC, STG and V1 were analyzed together to obtain the regional correction whereas NAc data were analyzed individually. DLPFC standardized residuals of deletion per 10K was significantly higher in both SZ and BD compared to control group [A]. There were no significant groups differences in STG, V1 and NAc [B, C, D]. All data are represented as least square mean (LSM) ± Std. error. p<0.05 and p<0.0001 are denoted by * and ****, respectively, otherwise p>0.05.
Complex I activity in SZ and BD
Complex I (NADH ubiquinone oxidoreductase) is the first enzyme complex for cellular oxidative phosphorylation and a major contributor to the proton gradient across the electron transport chain. It has been shown that mitochondria quality decreases with aging which in turn decreases the Complex I activity in various tissues including frontal cortex and muscle [13, 26]. Here, we sought to determine how schizophrenia and bipolar disorder affect the Complex I activity in four brain regions with adjustment of age and sex. Repeated measure analysis was performed with DLPFC, STG and V1 regions together while NAc data was analyzed separately. The effect of medication was analyzed for each four brain regions individually. The Complex I activity is significantly decreased in DLPFC of both SZ (p=0.028) and BD (p=0.012) groups compared to CON. SZ and BD did not affect the Complex I activity in STG, V1 and NAc (Figure 4) (Table ST5). The regional correlation of Complex I activity was not significant across DLPFC, STG, V1 and NAc (Table ST11). The STG Complex I activity was dependent on age while DLPFC, V1 and NAc Complex I activity was independent of age (Table ST11). An analysis of medications/drug effects on Complex I showed the DLPFC Complex I activity was affected by medications while there was no medications effect on STG, V1 and NAc Complex I activity (Table ST12). Thus, the effects of medication/drugs may be regionally specific.
Complex I activity in postmortem DLPFC, STG, V1 and NAc of control (CON), schizophrenia (SZ) and bipolar disorders (BD) groups. Mixed model analysis was performed for statistical significance with adjustment of age and sex. DLPFC Complex I activity was significantly lower in both SZ and BD compared to control group [A]. There were no significant groups differences in STG, V1 and NAc [B, C, D]. All data are represented as least square mean (LSM) ± Std. error. p<0.05 is denoted by *, otherwise p>0.05.
Overall Mitochondria Number, Synapse Number and Synaptic Mitochondria Number in SZ and BD
The synaptic strength and synaptic activity are dependent upon ATP production by mitochondria. Thus, mitochondria play a pivotal role in synaptic neurotransmission. Among four brain regions studied here, DLPFC has been the most affected area in terms of mtDNA CN, mtDNA deletions and Complex I activity. Therefore, we investigated if there is any change in overall mitochondria number, synapse number and synaptic mitochondria number in DLPFC of SZ and BD subjects. The overall mitochondria puncta number, synapse puncta number and synaptic mitochondria puncta colocalization were determined through immunostaining of mitochondria surface marker TOM40, synaptic marker PSD95 and colocalization of TOM40+PSD95, respectively. The mixed model analysis was conducted with adjustment of subjects’ age, pH and sex. For overall synapse number (# of PSD95 puncta), the main effect of group was significant. The number of synaptic marker (number of PSD95 puncta) was significantly lower in DLPFC of both SZ and BD (p=0.0008) subjects compared to control subjects (Figure 5B). For overall mitochondria number (# of TOM40), there was not significant main effect of group. The number of TOM40 was not significantly different in SZ (p=0.388) and BD (p=0.065) groups compared to CON group. For the synaptic mitochondria number (colocalized TOM40 with PSD95), the main effect of group was not significant. But very interestingly, the sex x group interaction was statistically significant. Both the female SZ (p=0.002) and female BD (p=0.266) subject had significantly lower colocalized TOM40 with PSD95 compared to CON female subjects. The male subjects with SZ and BD did not differ from male control subjects colocalized mitochondria (p = 0.85, p = 0.09, respectively) (Figure 5C-E) (Table ST6). The age did not have any significant effect on the number of PSD95, TOM40 and colocalized PSD95 with TOM40 (Table ST13). Similarly, the number of PSD95, TOM40 and colocalized PSD95 with TOM40 was not dependent on medications (Table ST14).
Overall mitochondria number, synapse number and synaptic mitochondria number in DLPFC of CON, SZ and BD subjects. The representative images for PSD95, TOM40 and PSD95 colocalized with TOM40 are depicted in [A]. Mixed model analysis was performed for statistical significance with adjustment of age and sex. The number of PSD95 was significantly lower in both SZ and BD compared to CON group [B]. The overall number of TOM40 was not significantly different in SZ and BD compared to CON, although there was a decrease trend in BD [C]. The sex x group effect was significant for the PSD95 colocalized with TOM40. Both the female SZ and female BD subjects had significantly lower PSD95 colocalized with TOM40 compared to CON females [D]. All data are represented as least square mean (LSM) ± Std. error. p<0.05, p<0.01 and p<0.001 are denoted by *, ** and ***, respectively, otherwise p>0.05. CON=control (n=19), SZ= schizophrenia (n=17), BD=bipolar disorder (n=15), Red spots= PSD95, green spots= TOM40, yellow spots= PSD95 colocalized with TOM40, blue=DAPI nuclear stain.
Discussion
Regional Alteration of mtDNA Copy Number in SZ and BD
In this manuscript, we report that postmortem brain tissue of schizophrenia patients showed increased mtDNA CN in DLPFC with no alterations of mtDNA CN in STG, V1 and NAc. On the other hand, bipolar patients showed increased mtDNA CN in DLPFC, decreased mtDNA CN in STG and unaltered mtDNA CN in V1 and NAc (Figure 1). These mtDNA CN results suggest that although schizophrenia and bipolar disorder have overlapping effect on mtDNA CN in various brain regions, but it is not the same for all brain regions. The decreased mtDNA CN has also been reported in postmortem hippocampus of BD patients [27]. Interestingly, peripheral mtDNA CN showed marked decrease in both SZ [21,28,29] and BD [20,30–32] while our postmortem DLPFC mtDNA CN showed marked increase. This discrepancy might be because the brain tissue has higher bioenergetics profile compared to peripheral tissue. Aging is one of the confounding factors affecting peripheral mtDNA CN. Peripheral mtDNA CN has been shown to decrease with aging [33]. Here, we have not found any effect of age on mtDNA CN in all four brain regions, although a trend for decrease with age was observed and warrants further investigation in larger cohorts (Table ST7). These results suggest that aging affect the mtDNA CN parameter in tissue specific manner. In both SZ and BD, the mtDNA copy number changes across brain regions are independent by correlation agreeing other lab reports [34]. The lack of correlation of mtDNA copy number in various brain regions suggests varying energetic demands in each region and perhaps variation in cell types. The mtDNA CN changes in this study were also independent of psychotropic medications or drugs in the postmortem brain (Table ST8).
The mtDNA common deletion is increased in DLPFC of BD while total mtDNA deletions are increased in DLPFC of both SZ and BD subjects
The mtDNA common deletion (CD) is a 4977 bp deletion which occur due to various events including mtDNA replication/repair, oxidative stress, exposure to ionizing radiation, and mitochondria disease. The mtDNA common deletion can also accumulate as a somatic heteroplasmic deletion with age (positively correlates with age). The age-related accumulation rate of mtDNA CD is tissue specific and more pronounced in tissues with higher energy demands. Here we sought to determine if there is any schizophrenia- or bipolar disorder-specific alteration of mtDNA CD, and at the same time, if there was any region-specific effect. We report here that the mtDNA CD remain unaltered in all four brain regions (DLPFC, STG, V1 and NAc) of SZ subjects. On the other hand, mtDNA CD is increased in DLPFC and unaltered in STG, V1 and NAc of BD subjects (Figure. 2). Previously, we showed that there was a global decrease of mtDNA CD in SZ subjects across 10 brain regions [12]. In contrast to the global effect previously reported in 10 brain regions in SZ, the mtDNA CD was analyzed in four brain regions reducing the statistical power of SZ subjects. The unaltered cortical mtDNA CD in SZ mirrors the previously published results from our and other labs [11,25,35]. The increased mtDNA CD in DLPFC of BD subjects is in accordance with previous findings from our lab determined with independent samples [11]. Previously we showed that the mtDNA common deletion was neither the most frequent nor the most abundant deletion in human brain tissues [18]. This has prompted us to explore the burden of total somatic deletions of mtDNA in schizophrenia and bipolar disorders. Here, we report that the total deletion burden (residuals of deletion per 10K coverage) is increased in DLPFC of both SZ and BD subjects (Figure 3). The residuals of deletion per 10K was not significantly altered in STG, V1 and NAc of both SZ and BD subjects. This is also consistent with our previous observation that the cumulative deletion burden of the DLPFC, relative to the same subject’s anterior cingulate cortex (ACC), was increased in a combined group of SZ and BD subjects when compared to subjects with no psychiatric disease, depression or alcohol dependency [18]. As expected, the deletions per 10K coverage was independent of age and medications in all four brain regions (Table ST10, Figure SF4).
Decreased Complex I activity in DLPFC of both SZ and BD
We have used a dual method to measure Complex I activity that takes into account the amount of Complex I protein present, then equalizes that amount in the second step of the assay and measures the conversion of NAD+ to NADH. This measure is analogous to a specific activity per unit of enzyme. We did not specifically test for protein abundance deficits, instead choosing to carefully normalize values of Complex I protein for activity measures. The mtDNA CD removes 7 genes encoding 4 Complex I subunits. Thus, we expected decreased Complex I activity in DLPFC of BD subjects since mtDNA CD was increased in DLPFC of BD subjects (Figure. 2). As expected, the Complex I activity was significantly decreased in DLPFC with no alteration in STG, V1 and NAc of BD subjects (Figure. 4). Although the mtDNA CD was unaltered in DLPFC of SZ subjects (Fig. 2), we found that the Complex I activity was also decreased in DLPFC with no alteration in STG, V1 and NAc of SZ subjects (Figure. 4). One possible reason might be due to the increased rate of total mtDNA deletion in DLPFC of SZ subjects as evident by increased residuals of deletion per 10K in DLPFC of SZ subjects. Our decreased DLPFC Complex I activity markedly contrasts with the peripheral platelet Complex I activity of SZ (active psychosis) and BD subjects where the Complex I activity is significantly increased [36–38]. Thus, the peripheral and brain Complex I activity might get affected by different molecular mechanism in blood and brain. The recent meta- analysis [26] showed that the SZ and BD had moderate (statistically insignificant) effect on Complex I activity. This discrepancy might be attributed to the regional heterogeneity of alteration in Complex I activity by SZ and BD. Our postmortem Complex I findings are in line with previously published reports showing decreased Complex I activity in DLPFC and temporal cortex of SZ and BD subjects [13,14,39,40].
Loss of synapses and sex-specific loss of synaptic mitochondria in both SZ and BD subjects
The energy produced by the mitochondria is utilized in neuronal communication such as development of cytoskeleton for presynapse development, synaptic vesicle transport and generation of synaptic membrane potentials. Mitochondria also plays a crucial role in synaptic neurotransmission through intracellular Ca2+ buffering. Thus, the number of mitochondria remain higher in high energy demanding sites such as axon terminals and postsynaptic areas compared to other parts of the neurons [41]. Fewer mitochondria number and/or dysfunctional mitochondria would lead to fewer synapse number or aberrant synaptic communication. In this manuscript, we identified if there is any alteration of the overall mitochondria number, synapse number and synaptic mitochondria number in SZ and BD after adjustment of age and sex. We report here that the overall synapse number are significantly decreased in DLPFC of both SZ and BD subjects. This is an important cellular hallmark in understanding the pathophysiology of SZ and BD. This finding mimics the recent studies showing the decrease of excitatory synapses or postsynaptic elements in prefrontal cortex and anterior cingulate cortex of schizophrenia subjects [42–44]. In contrast to the synapse number, the overall mitochondria number (as measured by IHC that was not restricted to excitatory synapses) was not significantly altered in SZ and BD compared to controls. Decreased mitochondria number has been reported in the axon terminals and layer 3 excitatory synapses of anterior cingulate cortex [40,44,45], and in the oligodendrocytes of prefrontal cortex of SZ subjects [46, 47]. This discrepancy might be attributed to the fact that our calculated overall mitochondria belong to all types of cells and all parts of cells compared to high energy demanding sites such as axon terminals or synapses. To calculate the synaptic mitochondria number, we measured the colocalized spots for mitochondria marker TOM40 with synapse marker PSD95. Interestingly, we have found that synaptic mitochondria number are decreased in both female SZ and female BD subject compared to female controls.
Conclusion
Among the four brain regions, the DLPFC remains the most affected brain region in terms of mitochondrial parameters studied in this manuscript. The measures of mitochondria unique deletions and copy number are both increased, while Complex I activity is decreased in the DLPFC of both schizophrenia and bipolar subjects compared to controls. Decreased excitatory synapses in schizophrenia and bipolar disorder are relevant to these findings. These indicators of mitochondria dysfunction require further investigations into cell specificity of the findings and potential relationship to genetic influences to understand the pathophysiology of schizophrenia and bipolar disorders.
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Funding and Disclosures
National Institute of Mental Disorder [R01MH08580 to MPV]; Della Martin Foundation [to SCD, postdoctoral fellowship]. Postmortem brain tissues collection at the University of California-Irvine (UCI) Brain Bank was supported by the Pritzker Neuropsychiatric Disorders Research Consortium.
Conflict of Interest Statement
The authors have nothing to disclose.
Authors Contributions
Designed, supervised study, and wrote manuscript (MPV). Conducted laboratory experiments (SCD, LM, BR, BH). Data analysis (SCD, BH, MPV, AAO). Drafted and edited manuscript (all).
Supplementary Information
Demographics of subjects
mtDNA copy number in four brain regions among control, schizophrenia and bipolar disorders.
mtDNA common deletion in four brain regions among control, schizophrenia and bipolar disorders.
mtDNA deletions per 10K in four brain regions among control, schizophrenia and bipolar disorders.
Complex I activity in four brain regions among control, schizophrenia and bipolar disorders.
Number of overall mitochondria number, synapse number and synaptic mitochondria number in DLPFC of CON, SZ and BD subjects.
Correlation of age on mtDNA copy number in DLPFC, STG, V1 and NAc.
Effect of Medications on mtDNA copy number in DLPFC, STG, V1 and NAc. 1 = Pooled BD and SZ subjects positive for neuroactive medications by toxicology screen 2= Pooled BD and SZ subjects negative for neuroactive medications by toxicology screen C = controls
Inter-region correlation and age-effect on mtDNA common deletion in DLPFC, STG, V1 and NAc.
Inter-region correlation and age-effect on deletions per 10K across DLPFC, STG, V1 and NAc.
Inter-region correlation and age-effect on Complex I activity across DLPFC, STG, V1 and NAc.
Effect of Medication/Drugs on complex I activity in DLPFC, STG, V1 and NAc.
Effect of Age on number of TOM40, PSD95 and coloc in DLPFC.
Effect of Medication/Drugs on number of Tomm40, PSD95 and Coloc in DLPFC.
mtDNA copy number correlation across DLPFC, V1 and STG.
mtDNA common deletion correlation across DLPFC, STG, V1 and NAc.
Effect of age on common deletions (%) in DLPFC, STG, V1 and NAc.
Effect of Medication/Drugs on deletions per 10K in DLPFC, STG and V1.
Supplementary Materials and Methods
mtDNA Copy Number and Common Deletion:
The genomic DNA was extracted from DLPFC, STG, V1 and NAc samples using a phenol chloroform extraction method. The mtDNA CN and mtDNA CD in DLPFC, STG and NAc samples was determined with quantitative real-time PCR (qPCR), whereas mtDNA CN and mtDNA CD in V1 samples was determined by droplet digital PCR. Quantitative real-time PCR was run using the Applied Biosystems 7900HT real-time instrument in 12-µL reaction volumes using 2 ng gDNA and 2 µl of each standard with published conditions in triplicate. For mtDNA CN, a 200-bp gBlock oligonucleotide (Integrated DNA Technologies, Coralville, IA, USA) representing the human albumin gene (ALB) on chromosome 4 is the reference DNA standard. An 11- point standard curve created with this gBlock has copies ranging from 50 - 50,000 copies per µl of the albumin gene. A 326-bp gBlock Gene Fragment (Integrated DNA Technologies, Coralville, IA) for the mitochondrial wild type sequence (mt13176-13501) is used as a copy number standard for mtDNA. An 11-point standard curve is created with copies ranging from 20,000-10,000,000 per µL of mitochondria wild type product. Mitochondrial DNA CN for the samples is determined using the ratio of mtDNA CN/ALB CN. For the mtDNA CD, a 301-bp gBlock Gene Fragment (Integrated DNA Technologies, Coralville, IA, USA) encompassing the common deletion breakpoints (mt8224–8469 and mt13447–13501) is used as a standard for the deletion. An 11-point standard curve is created with copies ranging from 50-50,000-copies/µL. CD% is calculated using the following equation: CD% = (mtDNA deletion copies)/ (mtDNA wild type copies + deletion copies) × 100.
The following primers were used:
mtDNA CN/CD based on SYBR Green chemistry:
ALB_Primer_F: 5‘-TGCATGAGAAAACGCCAGTA-3’
ALB_Primer_R: 5‘-TCTGCATGGAAGGTGAATGT-3’
mtWT_Primer_F: 5‘-TTACAATCGGCATCAACCAA-3’
mtWT_Primer_R: 5‘-GCTAATGCTAGGCTGCCAAT-3’
mtDel_Primer_F: 5‘-AGGGCCCGTATTTACCCTAT-3’
mtDel_Primer_R: 5‘-GCTAATGCTAGGCTGCCAAT-3’
mtDNA CN and mtDNA CD in V1 samples were determined using the QX200 Droplet Digital PCR System (Bio-Rad, Hercules, CA, USA). For each assay a mixture of both the forward and reverse primers and probes were combined. Stock concentrations of each primer and probe were 100 µM and to make 100 reactions worth of a 20x assay we used 18 µL of forward primer stock at 100 µM, 18 µL of reverse primer stock at 100 µM, 5 ul of probe stock at 100 µM and 59 µL of water. The final 1x concentration was 900 nM for each primer and 250 nM for each probe. Each 20 µL reaction consisted of 10 µL 2x ddPCR Supermix (Bio-Rad, Hercules, CA, USA), 1 µL of the above 20x assay consisting of the primers and probe, and the appropriate amount of gDNA and water up to the final 20 µL volume. To run the mitochondria wildtype assay, 0.01 ng of gDNA were used. To run the ALB and common deletion assays 10 ng of gDNA were used. Samples were placed in the QX200 Droplet Generator to partition each sample into 20,000 nanoliter sized droplets. Droplets were then run on the thermal cycler with the following program: 95 x 10 min, followed by 40 cycles of 94 x 30 sec then 60 x 1 min. A final enzyme deactivation for 98 x 10 min was followed by a hold at 4. Following PCR amplification, the droplets were analyzed in the QX200 Droplet Reader. QuantaSoft Software (Bio-Rad, Hercules, CA, USA) was used to measure the number of positive and negative droplets for each fluorophore.
The following primers and probe Integrated DNA Technologies, Coralville, IA, USA) were used for ddPCR (all as listed 5’-3’):
Primers for ALB: F-TGCATGAGAAAACGCCAGTA, R- TCTGCATGGAAGGTGAATGT Probe for ALB: /5HEX/TGGTGAACA/ZEN/GGCGACCATGCT/3IABkFQ/.
For the mito wildtype (WT) assay, the primers and probe were
Primers for WT: F- GGAAGCCTATTCGCAGGATT, R- CTAGGAAAGTGACAGCGAGG Probe for WT: /56-FAM/CCCCTTCCA/ZEN/AACAACAATCCCC/3IABkFQ/.
For the common deletion (CD) assay the primers and probe were-
Primers for CD: F-AGGGCCCGTATTTACCCTAT, R-GCTAATGCTAGGCTGCCAAT Probe for CD: /5HEX/CCCCTCTAG/ZEN/AGCCCACTGTAAAGC/3IABkFQ/.
Pooled Effect of Total Deletion of mtDNA by High-resolution Pipeline “Splice-Break”:
In this manuscript, we present the results of the number of deletion species per 10K coverage as this showed a significant effect across diagnosis. The number of deletions per 10K was calculated by counting the total number of deletion species (unique set of breakpoints) in a particular sample, then dividing by the benchmark coverage and multiplying by 10000. The Splice-Break pipeline included long-range (LR) PCR amplification of mtDNA, next-generation sequencing (NGS), and MapSplice (v.2.1.18) algorithm (which detect putative mtDNA deletion breakpoint “junctions”). Prior to NGS library preparation, the mtDNA was enriched for each sample using a long-range (LR) PCR utilizing back-to-back primers that hybridize to the control region of the mitochondrial genome. Primer sequences used for most samples were F 5’- CCGCACAAGAGTGCTACTCTCCTC-3’ and R 5’-GATATTGATTTCACGGAGGATGGTG-3’ (Integrated DNA Technologies). In cases where back-to-back primers failed to amplify the WT mitochondrial genome, the following alternate primers were used: 5’- AAGCCTAAATAGCCCACACG-3’ and 5’-GTGGCTTTGGAGTTGCAGTT-3’. For NGS, libraries were prepared using the 384 Seq-Well kit (San Diego, CA). Libraries are sequenced as 150-mer paired-end reads on a patterned flow cell using the Illumina NovaSeq6000 at the UCI Genomics High Throughput Facility. We used our custom bash script Splice-Break.sh (https://github.com/brookehjelm/Splice-Break/) to perform automated filtering, normalization, and annotation of mtDNA deletion breakpoints, and performed manual calculations of the cumulative deletion metrics.
Complex I Activity (activity/unit of Complex I protein):
Briefly, brain tissue (60–100 mg) was homogenized in PBS at 4°C and centrifuged at 3000 rpm for 5 min. After centrifugation, the first supernatant is collected and measured for protein concentration to ensure adequate protein for both ELISA and activity assays. Equivalent amounts of supernatant protein are extracted (Part 8201086, Abcam, Cambridge, UK), the second supernatants are collected (16000xg, 4°C, 20 min) and stored at –80°C for Complex I ELISA and activity assays. The Complex I ELISA assay (AB178011, Abcam, Cambridge, UK) quantitatively measures NADH dehydrogenase protein. Equal protein amounts are used with a specific monoclonal capture and detection antibodies reactive to CI isolates. All standards, samples, and controls are run in duplicate in 96-well plates pre-coated with a capture antibody provided in the Abcam kit. A standard curve prepared from the HeLa lyophilized cell extract as well as additional standards prepared from a postmortem human PFC brain extract are used. The standard CI ELISA protein curve is calculated from the regression of the average OD 450 nm reading of each HeLa standard against their concentrations. The second part of the assay is the CI activity assay. Each brain sample is diluted to the same CI protein (amount determined above by ELISA) to determine the rate of CI activity. All samples and controls are run in duplicate according to the manufacturer’s protocol (AB109721, Abcam, Cambridge, UK). The rate of activity is expressed as the change in absorbance per minute using this equation: Complex I activity for each sample (milli-optical density units, mOD/min) = (1,000 × ΔOD450/30 min) where Δ represents the final minus initial OD450.
Immunohistochemistry (IHC) for Colocalization of Mitochondria and Synaptic Marker:
Fresh frozen DLPFC tissue blocks (∼1 cm3) were immersion-fixed in freshly made 4% PFA (in PBS) for 24 hours at 4°C. Post-fixation, tissue blocks were cryopreserved in 20% sucrose and 0.1% Na-Azide in PBS for 3 days; then placed in 30% sucrose for 3 days at 4°C. After cryopreservation, blocks are stored frozen at -80°C until cryo-sectioning. Free-floating 10 µm tissue cryostat sections are collected in 0.1% sodium azide in PBS. Antigen-retrieval was performed through incubation of free-floating tissue sections in 1 mM EDTA buffer (pH 8) at 100°C and high pressure for 5 min. After antigen retrieval, slices were incubated in PBS-Triton (0.3%) for 1 hour. Slices were blocked with 10% goat serum. Primary antibodies for excitatory postsynaptic markers (PSD95) (MA1-045, Thermofisher) and mitochondria (TOM40) (ab185543, Abcam) were added and incubated for 24 hours at 4°C. Next day, slices were washed with PBS three times followed by respective secondary antibodies incubation for 24 hours at 4°C. After washing secondary antibodies, slices were incubated with DAPI as a nuclear counterstain for 10 min. Next, slices were incubated in 70% ethanol for 5 min followed by incubation with auto- fluorescence eliminator reagent (Sigma Aldrich, Cat#2160) for 1 min. Slides were mounted in Prolong diamond antifade solution (Thermofisher, P36970). Digital z-stack images are collected using a Zeiss LSM 780 confocal microscope with oil X 63 objective (NA=1.4). Acquired image parameters: voxel size 0.07 x 0.07 x 0.32 micron3, 1024 x 1024 pixels with a depth of 12 bit. Colocalization analysis was performed through IMARIS software (Bitplane). The spots of TOM40 and PSD95 with diameter of greater than 1.2 µm were excluded from the analysis. The colocalization was defined as the number of TOM40 spots within the distance of 0.25 µm from PSD95 spots. 5-8 images from 2-3 sections per subject were captured. The total number of subjects evaluated for IHC was 51.
Acknowledgements
We would like to thank donors’ family for donating the postmortem brain tissues.