Kazusa Miyahara1
Mizuki Hino1,2
Zhiqian Yu3
Chiaki Ono3
Atsuko Nagaoka2Masataka Hatano2
Risa Shishido2
Hirooki Yabe2
Hiroaki Tomita3,4
Yasuto Kunii1,2*- 1Department of Disaster Psychiatry, International Research Institute of Disaster Science, Tohoku University, Sendai, Japan
- 2Department of Neuropsychiatry, School of Medicine, Fukushima Medical University, Fukushima, Japan
- 3Department of Psychiatry, Graduate School of Medicine, Tohoku University, Sendai, Miyagi, Japan
- 4Department of Psychiatry, Tohoku University Hospital, Sendai, Miyagi, Japan
Background: Evaluating and controlling confounders are necessary when investigating molecular pathogenesis using human postmortem brain tissue. Particularly, tissue pH and RNA integrity number (RIN) are valuable indicators for controlling confounders. However, the influences of these indicators on the expression of each gene in postmortem brain have not been fully investigated. Therefore, we aimed to assess these effects on gene expressions of human brain samples.
Methods: We isolated total RNA from occipital lobes of 13 patients with schizophrenia and measured the RIN and tissue pH. Gene expression was analyzed and gene sets affected by tissue pH and RIN were identified. Moreover, we examined the functions of these genes by enrichment analysis and upstream regulator analysis.
Results: We identified 2,043 genes (24.7%) whose expressions were highly correlated with pH; 3,004 genes (36.3%) whose expressions were highly correlated with RIN; and 1,293 genes (15.6%) whose expressions were highly correlated with both pH and RIN. Genes commonly affected by tissue pH and RIN were highly associated with energy production and the immune system. In addition, genes uniquely affected by tissue pH were highly associated with the cell cycle, whereas those uniquely affected by RIN were highly associated with RNA processing.
Conclusion: The current study elucidated the influence of pH and RIN on gene expression profiling and identified gene sets whose expressions were affected by tissue pH or RIN. These findings would be helpful in the control of confounders for future postmortem brain studies.
1. Introduction
There has been increasing interest in studies on neuropsychiatric disorders conducted on the postmortem brain over the years. Schizophrenia was once called “the graveyard of neuropathologists” in 1970s, because only a few neuropathological characteristics were observed in the postmortem brain of diagnosed patients (1). However, genetic and molecular analyses have been well developed since the beginning of the 21st century (2) and postmortem brains of patients with schizophrenia were re-examined to validate genetic and molecular findings. Now, investigating the pathology of the human postmortem brain at the molecular level is critical for elucidating the pathophysiology of schizophrenia. For example, abnormal expression level of transcripts related to various biological components including γ-aminobutyric acid (GABA) signaling (3, 4), glutamatergic signaling (5), inflammatory system (6), astrocyte (7), oligodendrocyte (8), microglia (9) have been repeatedly reported and it is certain that these system are profoundly included in pathophysiology of schizophrenia. Moreover, effects of schizophrenia-associated single nucleotide polymorphisms (SNPs) on gene expression in the brain have been investigated and the increase of SNPs has been associated with aberrant gene expression and the abnormal brain structure of patients with schizophrenia (10–12). Furthermore, comprehensive analyses on gene expression have clarified the pathways that are closely associated with schizophrenia and identified the potential therapeutic agents (13, 14). There is a growing trend toward open accessibility to a large dataset of human brain such as CommonMind Consortium (15) and BrainCloud (16), lessening the hurdles for conducting postmortem brain studies. It is predicted that studies on postmortem brain would increase and would contribute to clarifying the pathophysiology of schizophrenia and other neuropsychiatric disorders.
Despite the benefits of postmortem brain studies, several inevitable issues may arise, such as heterogeneity among subjects. Various conditions other than pathophysiology of diseases, such as antemortem history, agonal state, and postmortem interval (PMI) are reflected on molecules in the postmortem brain. Therefore, it is essential to adequately evaluate and control for confounders when conducting postmortem brain studies.
Tissue pH and RNA integrity number (RIN) (17) are the most important confounders among the others, such as age of death, sex, PMI, and antemortem drugs, since these two reflects the severity of the tissue degradation (18). Tissue pH mostly depends on the agonal state (19) but not on the PMI (20). The agonal factor score, developed by Tomita et al., quantifies the influence of the agonal factors and is associated with tissue pH, demonstrating that the manner of death affects tissue pH (21, 22). Additionally, tissue pH is known to be influenced by smoking history (23) and intake of antipsychotics (24). Correspondingly, RIN has been associated with the agonal state (21) and the condition of the postmortem brain tissue after death (25). RIN and tissue pH have been known to affect each other. In previous studies, poor RNA quality was closely associated with a drop in pH (20, 26). The association between low tissue pH and RIN may be partially explained by the effect of agonal state, including hypoxia, coma, and dehydration. Severe agonal states would expose brain cells to hypoxia, causing changes of gene expression level and cell death. Hypoxia results in lower pH due to increase of lactate production (27, 28), while cell death also causes acidosis via release of intercellular vesicle (18). In addition, destruction of RNA along with cell damage itself may result in low tissue pH, since RNA consists of phosphoric acid. Therefore, both tissue pH and RIN reflects the severity of agonal state, which may be indicators of severity of agonal states and according changes of gene expressions. Previous studies have proven that various confounders are reflected on tissue pH and RIN. Thus, these two are essential confounders for heterogeneity adjustment.
Most postmortem brain studies include tissue pH and/or RIN as confounders; however, the effects of these two factors on gene expressions have not been fully investigated. Several studies have assessed their effects on mRNA by reverse transcription polymerase chain reaction (29–31), but the number of target genes is limited and usually under 100. A few studies (32–35) had comprehensively evaluated the effect of tissue pH on gene expression using microarray; Iwamoto et al. (33) concluded that tissue pH was associated with the altered expressions of genes involved in energy metabolism, mitochondrial function, proteolytic activities, stress-response proteins, and transcription factors, although they did not include RIN as a cofounder (34). Therefore, further studies are needed to precisely evaluate the influence of tissue pH and RIN on global gene expressions provided that control and experimental groups are examined separately.
This study aimed to identify the gene sets that were affected by tissue pH and RIN using the postmortem brain of patients with schizophrenia. We conducted microarray analysis, which was able to measure the expression levels of low-expressed genes unlike RNA-seq (36), and identified genes that were significantly correlated with tissue pH, RIN, or both. Furthermore, we assessed the function of each gene set by pathway analysis and upstream regulator analysis.
2. Subjects and methods
2.1. Human postmortem brain tissue
Postmortem brain tissue samples from patients with schizophrenia were obtained from the Fukushima Brain Bank, Department of Neuropsychiatry, School of Medicine, Fukushima Medical University, Fukushima, Japan (37). The use of postmortem human brain tissue in this study was approved by the Ethics Committee of Fukushima Medical University and Tohoku University, and complied with the principles of the Declaration of Helsinki and its later amendments. The ethics approval numbers for this study were 2021–1-127 (Tohoku University) and 1,685 (Fukushima Medical University). All procedures were performed with the written consent of close relatives. Detailed demographic information on brain tissue from the 13 subjects with schizophrenia used in this study is summarized in Table 1. The patients with schizophrenia fulfilled the diagnostic criteria established by the American Psychiatric Association (Diagnostic and Statistical Manual of Mental Disorders: DSM-IV).
Table 1. Demographic data of the samples.
2.2. Tissue pH measurement of postmortem brain
The pH of postmortem brain tissue was measured as described in a previous study (38) with slight modifications. In brief, approximately 20 mg of frozen tissue was dissected and homogenized in a 5-times volume (1 ml/100 mg tissue) of nuclease-free water (Ambion, Austin, TX, USA). The pH of the homogenate was measured using Twin pH-B212 (Horiba, Kyoto, Japan).
2.3. RNA isolation and RNA quality measurements
Total RNA was isolated by using an AllPrep DNA/RNA/Protein Mini Kit (Qiagen, Valencia, CA, USA), and genomic DNA in the RNA samples was digested by RNase-free DNase I. The RIN was measured using the 2,100 Bioanalyzer (Agilent, Santa Clara, CA, USA) as described in a previous study (39).
2.4. Gene expression microarray analysis
Gene expression analysis was conducted as described in a previous study (40). In brief, biotinylated cRNA was prepared using the Ambion Illumina RNA Amplification Kit (Ambion, Austin, TX, USA). The biotinylated cRNA were hybridized to Human 6v2 Expression BeadChips (Illumina, San Diego, CA, USA). The BeadChips were then labeled with streptavidin–Cy3 (GE Healthcare, Chicago, IL, USA), washed, and scanned with Bead Station 500X (Illumina). Signal intensity (SI) of each bead was measured and normalized using average normalization by BeadStudio software Ver. 3.1 (Illumina) as previous studies (40, 41). Among the 48,702 probes, 8,286 probes had an average SI ≥ 100.
2.5. Statistical analysis
Correlation analysis of gene expression levels between samples was performed without applying a cutoff based on SI. For all other analyses, beads with an average SI less than 100 were excluded. Correlation analysis was performed by calculating Pearson’s correlation coefficients. For the multiple comparison tests, Benjamini–Hochberg corrections were performed. Statistical significance was defined as q value <0.05. As to gene sets which were correlated with tissue pH and RIN, we separately conducted pathway analysis using Reactome Pathway Database1 (42) and upstream regulator analysis using Ingenuity Pathway Analysis (IPA). We constructed a network by using Cytoscape (43) and visualized the relation between genes and reactome pathways. R programming and Excel were used for data processing.
3. Results
3.1. Influence of tissue pH and RIN on correlation of gene expression profiling
First, we investigated whether tissue pH and RIN similarly affected gene expression profiling on each subject. Figure 1 shows a correlation of gene expression profiling along with tissue pH and RIN of each subject. Correlation coefficients were high between the pair with similar tissue pH and RIN.
Figure 1. Heatmap of the correlation of gene expression profiling along with tissue pH and RIN of each subject. The numbers represent the correlation coefficients.
3.2. Correlation analysis between gene expression and covariates
To identify gene sets susceptible to tissue pH and RIN, we conducted correlation analysis and elucidated genes whose expression was significantly correlated with tissue pH and/or RIN (Supplementary Tables S1–S3). The expression of 2,043 genes; 3,004 genes; and 1,293 genes showed a significant correlation with tissue pH, RIN, and both, respectively (Figures 2A–C). Conversely, 4,532 genes were not significantly correlated with pH or RIN. Among the genes susceptible to tissue pH and RIN, LMBR1 and PPP3CC were most highly correlated with tissue pH and RIN, respectively (Figure 3A: LMBR1; r = 0.91, Figure 3B: PPP3CC; r = 0.95).
Figure 2. (A) Venn diagram of the number of genes susceptible to tissue pH, RIN, and both. (B,C) Volcano plots of correlation coefficient and value of p evaluated between gene expressions and (B) pH and (C) RIN. Significant genes are highlighted.
Figure 3. Scatter plots showing the correlation between (A) tissue pH and gene expression of LMBR1 (r = 0.91), and (B) RIN and PPP3CC (r = 0.95). SI, signal intensity.
3.3. Pathway analysis on covariates-susceptible genes
To examine the function of identified gene sets susceptible to tissue pH and RIN, we carried out pathway analysis using Reactome Pathway Database. Genes susceptible to tissue pH were significantly enriched in pathways that were related to energy production, immune system, and cell cycle (Figure 4A). Genes susceptible to RIN were significantly enriched in pathways that were related to energy production, immune system, and RNA processing (Figure 4B). Genes susceptible to both tissue pH and RIN were significantly enriched in pathways that were related to energy production and immune system, including nuclear factor kappa B (NF-kappaB) signaling (Figure 4C). Supplementary Figures S1, S2 show concept networks of the gene set susceptible to tissue pH and RIN, which reflect correlation coefficient of each gene with top 6 pathways. The full results of enrichment analysis are shown in Supplementary Tables S4–S6. As a result of upstream regulator analysis, hypoxia inducible factors (HIF) 1a was found to be a significant regulator of genes susceptible to tissue pH (p = 9.1 × 10−4), although regulation of its downstream molecules was varied (Supplementary Table S7). Significantly-correlated genes belonging to the hypoxia-or NF-kappaB-related pathway are visualized as a network with coefficient correlation and value of p (pH: Figures 5A,B; RIN: Figures 5C,D). Genes belonging to pathway related to HIFs included those coding proteins associated with ubiquitin–proteasome pathway, such as proteasome 26S subunit ubiquitin receptor, non-ATPase 4 (PSMD4) and proteasome 20S subunit beta 6 (PSMB6). On the other hand, genes belonging to pathway related to NF-kappaB include those coding tumor necrosis factor receptor superfamily (TNFRSF), such as TNFRSF1B and TNFRSF6B.
Figure 4. Results of pathway analyses of the gene sets susceptible to (A) tissue pH, (B) RIN, and (C) both.
Figure 5. Network of genes significantly correlated with pH belonging to the (A) hypoxia-or (B) NF-kappaB-related pathway and network of genes significantly correlated with RIN belonging to the (C) hypoxia-or (D) NF-kappaB-related pathway. Node color indicates the coefficient correlation. Node size indicates the significance of the correlation; larger sizes indicate smaller p-values.
4. Discussion
In this study, we used the postmortem brains of 13 patients with schizophrenia and examined the influence of tissue pH and RIN on global gene expression profiling. We identified 2,043 genes (24.7%) susceptible to tissue pH; 3,004 genes (36.3%) susceptible to RIN; and 1,293 genes (15.6%) susceptible to both, and the remaining genes were not susceptible to pH or RIN. Additionally, we conducted pathway analysis and clarified that genes susceptible to pH were enriched in pathways that were related to energy production, immune system, and cell cycle, whereas genes susceptible to RIN were significantly enriched in pathways that were related to energy production, immune system, and RNA processing. To the best of our knowledge, this study is the first to report on these findings and would have an impact on the future postmortem studies for controlling confounders.
The characteristics of the gene sets susceptible to tissue pH and RIN had commonalities. Both gene sets were enriched in pathways associated with energy production and immune system. This may be attributed to the hypoxia exposure during the agonal state affecting tissue pH, RIN, and specific gene expression profiling. Hypoxia exposure is known to induce abnormal energy production and immune system functioning (44, 45).
Hypoxia-dependent change of energy production is mainly due to HIFs, a transcription factor activated during conditions of hypoxia (46). Although HIFs are usually hydroxylated by prolyl hydroxylases (PHD) and degradated by ubiquitin-proteasome system, HIFs are activated in hypoxic conditions because exposure decreases PHD enzymatic activity. HIFs are known to impact energy metabolism. For example, HIFs upregulates lactate dehydrogenase and downregulates pyruvate dehydrogenase, resulting in decreased energy production in the tricarboxylic acid cycle and electron transport system (44). Because tissue pH and RIN were significantly correlated with gene expression related to energy production, it is suggested that these confounders have reflected metabolic change during the agonal state. Since variable proteasome proteins were positively correlated with tissue pH and RIN as shown in Figures 5A,C, activation of proteasome may be crucial for adaptation to hypoxia, which is consistent with previous studies reporting that PSMD4 and PSMB6 are activated in hypoxia environment (47, 48). In addition to regulation of metabolic enzyme due to HIFs, miRNAs located downstream of HIFs (called hypoxamiRs) are also acting in adaptation to hypoxia (49, 50). hypoxamiRs are small non-coding RNAs that regulate gene expression by forming complexes and degrading mRNAs with complementary sequences. Such hypoxamiRNAs effect on metabolic and immune pathways by downregulation of specific mRNA (51, 52). The above adaptive responses to hypoxia exposure occur when downstream genes of HIFs, including hypoximiRNAs, alter their gene expression patterns. In vivo experiments have shown that HIF1 and HIF2 play a role in the production of energy sufficient for survival in hypoxic environments (53) and in recovery from hypoxic exposure in the brain (54).
Meanwhile, the link between schizophrenia and aberrant energy production has been repeatedly reported (55, 56). In previous postmortem brain studies, the hypoxia-related pathway was activated and mitochondrial dysfunction was seen in the postmortem brains of patients with schizophrenia, despite controlling multiple confounders, including tissue pH and RIN (57, 58). Thus, it is likely that schizophrenia, to some extent, impacts the metabolism process in the human brain. Nevertheless, tissue pH and RIN may be critical markers for correcting the agonal metabolic change to precisely distinguish pathophysiology of schizophrenia and the artifact of the postmortem brain study.
Hypoxia exposure also induces immune response. Inhibition of PHD activates NF-kappaB, an essential transcription factor for inducing immune response (45). Since HIFs and NF-kappaB promote expression of the each other (59), hypoxia response and immune response are considered closely related (60). Among the genes belonging to pathways of NF-kappaB, TNFRSF were negatively correlated with pH and RIN, as shown in Figures 5B,D. TNFRSF has been known to include essential molecules for immune activity (61), and the current results suggest TNFRSF also plays a crucial role for reaction to hypoxia. NF-kappaB is a molecule often associated with schizophrenia. NF-kappaB was observed to be over activated while its inhibitor, human immunodeficiency virus enhancer binding protein 2 (HIVEP2), was inactivated in the postmortem brain of patients with schizophrenia (62, 63). Additionally, Murphy CE et al. reported that low tissue pH was associated with the high expression of NF-kappaB in patients with schizophrenia (64), which is consistent with our findings that tissue pH affects genes related to the immune system including NF-kappaB signaling. Therefore, activation of NF-kappaB is important for the pathophysiology of schizophrenia and its expression level should be controlled for tissue pH and RIN.
The pathways unique to pH-susceptible genes were related to cell cycle. Investigations on cells outside the brain have shown that hypoxia exposure affects cell cycle transition (65–67). For example, NF-kappa B pathway activated by hypoxia exposure promotes cell cycle to enhance the cell proliferation in pulmonary arterial smooth muscle cells. It is likely that similar cell cycle promotion had occurred during the agonal state as a response to hypoxia exposure. Conversely, RNA processing was the pathway unique to RIN-susceptible genes. This may be attributed to the disruption of RNA dynamics during the agonal state due to the reduction of RIN induced abnormal expression of the gene set related to RNA processing.
Our finding, that nearly 5,000 genes were affected by tissue pH or RIN, underlines the importance of controlling the confounders. These gene sets include genes commonly used as reference genes for postmortem brain studies (68), such as ACTB, B2M, POLR2A, and RPLP0CTB. In these for genes, expression levels of ACTB and RPLP0 were correlated with both tissue pH (ACTB: r = 0.75; RPLP0: r = −0.69) and RIN (ACTB: r = −0.90; RPLP0: r = 0.70). Expression level of POLR2A were correlated with only tissue pH (r = −0.68), while that of B2M were correlated with only RIN (r = 0.72). On the other hand, other reference genes like GAPDH and TFRC were not included. These results suggest that some reference genes were susceptible to tissue pH, RIN or both. Therefore, when conducting normalization of gene expression level using reference genes, it should be noted that reference genes may be affected by agonal states. In addition, the current identified gene sets include genes that have been linked to schizophrenia. For example, we demonstrated that GAD, protein phosphatase 3 catalytic subunit gamma (PPP3CC) and glutamate ionotropic receptor NMDA type subunit (GRIN) were susceptible to tissue pH and RIN, all of which were associated with schizophrenia in previous postmortem brain studies (3, 69, 70). When examining the expression level of these genes in postmortem brain, it is essential to control tissue pH and RIN as confounders. Cyclic AMP-responsive element binding protein (CREBBP) was identified as RIN-susceptible gene; however, its expression level has not been investigated while controlling RIN as a continuous confounder (71, 72). Therefore, the current insight that CREBBP is upregulated in patients with schizophrenia (73) may not be a solid theory and further studies are required.
Moreover, our findings may be the key for integrating inconsistent insights about specific gene expression level. For example, synaptophysin (SYP), an essential molecule for neurotransmission, has been reported to be reduced in patients with schizophrenia (74–76); however, a contradicting report (77) did not include tissue pH as a confounder. Because this study demonstrated that SYP was susceptible to tissue pH, considering adequate confounders may consolidate ambiguous findings.
The gene set identified susceptible to tissue pH in this study may reflect not only agonal state but also pathophysiological change due to schizophrenia. Although we considered tissue pH and RIN as confounders, several previous studies have reported that low pH may be a part of pathophysiology of psychosis (66, 78, 79). A study (78) conducted a meta-analysis evaluating tissue pH of postmortem brain from patients with schizophrenia and bipolar disorder and demonstrated that tissue pH of the patients was significantly lower than that of controls even when several confounders such as PMI and antemortem intake of antipsychotics were considered. Lower tissue pH in brain of patients with schizophrenia has been associated with aberrant metabolism and hypoxia-related changes in gene expression (58), which are overlapped with the transcriptome profile associated with pH and RIN in our current study. Therefore, caution is needed to interpret gene expression data from the postmortem brain, especially when the genes associated with pH and RIN would be detected as differentially expressed between tissues from control and patients with schizophrenia.
Our study has several limitations. First, effects of other confounders such as PMI and tissue processing could not be excluded. Second, we were unable to investigate the effects of tissue pH and RIN on all transcripts. Future studies are needed to fully investigate the effect of tissue pH and RIN using RNA-seq. Third, we evaluated gene expression profile in a single region. Influence of tissue pH and RIN in multiple regions remains unknown. Furthermore, as we used postmortem brain of patients with schizophrenia, there might be a potential effect of pathophysiology of the disease. Thus, further investigations are required to confirm our findings.
It should be noted that the exact magnitude of effect was not examined for individual genes. Ideally, it would be necessary to measure individual genes expression levels accurately by quantitative polymerase chain reaction and quantify protein expression levels by immunohistochemistry or western blot. Such additional experiments would clarify the details of the changes in gene expression patterns and are necessary for future postmortem brain research.
In conclusion, our study provides new insight into the effect of tissue pH and RIN on gene expression profile in postmortem brain. Our findings indicate that tissue pH and RIN oscillate the expression level of specific gene sets and that these two are critical confounders in postmortem brain analysis. These findings promote the need to adequately control confounders, subsequently paving the way to upgraded postmortem brain studies for elucidating pathophysiology of schizophrenia.
Data availability statement
The data presented in the study are deposited in the Gene Expression Omnibus database, accession number GSE235055 and are available at the following URL: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE235055.
Ethics statement
The studies involving human participants were reviewed and approved by ethics committees of Fukushima Medical University and the Tohoku University Graduate School of Medicine. The patients/participants provided their written informed consent to participate in this study.
Author contributions
KM, MHa, MHi, CO, ZY, HT, and YK designed the study. MHi, ZY, and YK performed the experiments. MHa, MHi, RS, AN, HY, HT, and YK collected postmortem brain samples and clinical information. KM, MHi, ZY, and YK undertook the statistical analysis. KM wrote the first draft. All authors contributed to and have approved the final manuscript.
Funding
This work was supported by Japan Agency for Medical Research and Development (AMED) – JP22dm0207074 (YK) and AMED – JP22wm0425019 (HY); the Grant-in-Aid for Scientific Research on Innovative Areas from the Ministry of Education, Culture, Sports, Science, and Technology of Japan under Grant Number JP221S0003, JP24116007 (HT), and JP21H00180 (YK), the Promotion of Science Grant-in-Aid for Scientific Research on Priority Areas (20023005); and the Collaborative Research Project of Brain Research Institute, Niigata University under Grant Number 2017–2806 (HT), and 22002 (YK).
Acknowledgments
We thank H. Onuma for her contribution in coordinating for donations. We also thank the families of the deceased for the donation of brain tissues and their time and effort devoted to the consent process and interviews.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2023.1156524/full#supplementary-material
Abbreviations
CREBBP, cyclic AMP-responsive element binding protein; GAD, glutamic acid decarboxylase; HIF, hypoxia inducible factors; IPA, ingenuity Pathway Analysis; PMI, postmortem interval; PSMB6, proteasome 20S subunit beta 6; PSMD4, proteasome 26S subunit ubiquitin receptor, non-ATPase 4; RIN, RNA integrity number; SI, signal intensity; TNFRSF, tumor necrosis factor receptor superfamily.
Footnotes
References
1. Plum, F. Prospects for research on schizophrenia. 3. Neurophysiology. Neuropathological findings. Neurosci Res Program Bul. (1972) 10:384–8.
2. Harrison, PJ, and Weinberger, DR. Schizophrenia genes, gene expression, and neuropathology: on the matter of their convergence. Mol Psychiatry. (2005) 10:40–68. doi: 10.1038/sj.mp.4001558
3. de Jonge, JC, Vinkers, CH, Hulshoff Pol, HE, and Marsman, A. GABAergic mechanisms in schizophrenia: linking postmortem and in vivo studies. Front Psych. (2017) 8:118. doi: 10.3389/fpsyt.2017.00118
4. Yu, Z, Ueno, K, Funayama, R, Sakai, M, Nariai, N, Kojima, K, et al. Sex-specific differences in the transcriptome of the human dorsolateral prefrontal cortex in schizophrenia. Mol Neurobiol. (2023) 60:1083–98. doi: 10.1007/s12035-022-03109-6
5. De Rosa, A, Fontana, A, Nuzzo, T, Garofalo, M, Di Maio, A, Punzo, D, et al. Machine learning algorithm unveils glutamatergic alterations in the post-mortem schizophrenia brain. Schizophrenia (Heidelb). (2022) 8:8. doi: 10.1038/s41537-022-00231-1
6. Trepanier, MO, Hopperton, KE, Mizrahi, R, Mechawar, N, and Bazinet, RP. Postmortem evidence of cerebral inflammation in schizophrenia: a systematic review. Mol Psychiatry. (2016) 21:1009–26. doi: 10.1038/mp.2016.90
7. Toker, L, Mancarci, BO, Tripathy, S, and Pavlidis, P. Transcriptomic evidence for alterations in astrocytes and Parvalbumin interneurons in subjects with bipolar disorder and schizophrenia. Biol Psychiatry. (2018) 84:787–96. doi: 10.1016/j.biopsych.2018.07.010
8. Iwamoto, K, Bundo, M, Yamada, K, Takao, H, Iwayama-Shigeno, Y, Yoshikawa, T, et al. DNA methylation status of SOX10 correlates with its downregulation and oligodendrocyte dysfunction in schizophrenia. J Neurosci. (2005) 25:5376–81. doi: 10.1523/JNEUROSCI.0766-05.2005
9. Sakai, M, Takahashi, Y, Yu, Z, and Tomita, H. Microglial gene expression alterations in the brains of patients with psychiatric disorders. Adv Neuroimmune Biol. (2016) 6:83–93. doi: 10.3233/NIB-160110
10. Komatsu, H, Takeuchi, H, Kikuchi, Y, Ono, C, Yu, Z, Iizuka, K, et al. Ethnicity-dependent effects of schizophrenia risk variants of the OLIG2 gene on OLIG2 transcription and white matter integrity. Schizophr Bull. (2020) 46:1619–28. doi: 10.1093/schbul/sbaa049
11. Chen, J, Calhoun, VD, Lin, D, Perrone-Bizzozero, NI, Bustillo, JR, Pearlson, GD, et al. Shared genetic risk of schizophrenia and gray matter reduction in 6p22.1. Schizophr Bull. (2019) 45:222–32. doi: 10.1093/schbul/sby010
12. Dobbyn, A, Huckins, LM, Boocock, J, Sloofman, LG, Glicksberg, BS, Giambartolomei, C, et al. Landscape of conditional eQTL in dorsolateral prefrontal cortex and co-localization with schizophrenia GWAS. Am J Hum Genet. (2018) 102:1169–84. doi: 10.1016/j.ajhg.2018.04.011
13. Radulescu, E, Jaffe, AE, Straub, RE, Chen, Q, Shin, JH, Hyde, TM, et al. Identification and prioritization of gene sets associated with schizophrenia risk by co-expression network analysis in human brain. Mol Psychiatry. (2020) 25:791–804. doi: 10.1038/s41380-018-0304-1
14. Yu, H, Guo, Y, Chen, J, Chen, X, Jia, P, and Zhao, Z. Rewired pathways and disrupted pathway crosstalk in schizophrenia transcriptomes by multiple differential coexpression methods. Genes. (2021) 12:665. doi: 10.3390/genes12050665
15. Hoffman, GE, Bendl, J, Voloudakis, G, Montgomery, KS, Sloofman, L, Wang, Y-C, et al. Common mind consortium provides transcriptomic and epigenomic data for schizophrenia and bipolar disorder. Scientific Data. (2019) 6
16. Colantuoni, C, Lipska, BK, Ye, T, Hyde, TM, Tao, R, Leek, JT, et al. Temporal dynamics and genetic control of transcription in the human prefrontal cortex. Nature. (2011) 478:519–23. doi: 10.1038/nature10524
17. Schroeder, A, Mueller, O, Stocker, S, Salowsky, R, Leiber, M, Gassmann, M, et al. The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. (2006) 7:3. doi: 10.1186/1471-2199-7-3
18. Sheedy, D, Harding, A, Say, M, Stevens, J, and Kril, JJ. Histological assessment of cerebellar granule cell layer in postmortem brain; a useful marker of tissue integrity? Cell Tissue Bank. (2012) 13:521–7. doi: 10.1007/s10561-011-9265-1
19. Elaine, K, Perry, RHP, and Tomlinson, BE. The influence of agonal status on some neurochemical activities of postmortem human brain tissue. Neurosci Lett. (1982) 29:303–7. doi: 10.1016/0304-3940(82)90334-2
20. Kingsbury, AE, Foster, OJ, Nisbet, AP, Cairns, N, Bray, L, Eve, DJ, et al. Tissue pH as an indicator of mRNA preservation in human post-mortem brain. Brain Res Mol Brain Res. (1995) 28:311–8. doi: 10.1016/0169-328X(94)00219-5
21. Tomita, H, Vawter, MP, Walsh, DM, Evans, SJ, Choudary, PV, Li, J, et al. Effect of agonal and postmortem factors on gene expression profile: quality control in microarray analyses of postmortem human brain. Biol Psychiatry. (2004) 55:346–52. doi: 10.1016/j.biopsych.2003.10.013
22. Stan, AD, Ghose, S, Gao, XM, Roberts, RC, Lewis-Amezcua, K, Hatanpaa, KJ, et al. Human postmortem tissue: what quality markers matter? Brain Res. (2006) 1123:1–11. doi: 10.1016/j.brainres.2006.09.025
23. Lipska, BK, Deep-Soboslay, A, Weickert, CS, Hyde, TM, Martin, CE, Herman, MM, et al. Critical factors in gene expression in postmortem human brain: focus on studies in schizophrenia. Biol Psychiatry. (2006) 60:650–8. doi: 10.1016/j.biopsych.2006.06.019
24. Halim, ND, Lipska, BK, Hyde, TM, Deep-Soboslay, A, Saylor, EM, Herman, MM, et al. Increased lactate levels and reduced pH in postmortem brains of schizophrenics: medication confounds. Neurosci Methods. (2008) 169:208–13. doi: 10.1016/j.jneumeth.2007.11.017
25. Barton, AJ, Pearson, RC, Najlerahim, A, and Harrison, PJ. Pre-and postmortem influences on brain RNA. J Neurochem. (1993) 61:1–11. doi: 10.1111/j.1471-4159.1993.tb03532.x
26. Chevyreva, I, Faull, RL, Green, CR, and Nicholson, LF. Assessing RNA quality in postmortem human brain tissue. Exp Mol Pathol. (2008) 84:71–7. doi: 10.1016/j.yexmp.2007.08.019
27. Mutch, WA, and Hansen, AJ. Extracellular pH changes during spreading depression and cerebral ischemia: mechanisms of brain pH regulation. J Cereb Blood Flow Metab. (1984) 4:17–27. doi: 10.1038/jcbfm.1984.3
28. Schurr, A, Payne, RS, Miller, JJ, and Rigor, BM. Glia are the main source of lactate utilized by neurons for recovery of function posthypoxia. Brain Res. (1997) 774:221–4. doi: 10.1016/S0006-8993(97)81708-8
29. Weickert, CS, Sheedy, D, Rothmond, DA, Dedova, I, Fung, S, Garrick, T, et al. Selection of reference gene expression in a schizophrenia brain cohort. Aust N Z J Psychiatry. (2010) 44:59–70. doi: 10.3109/00048670903393662
30. Preece, P, and Cairns, NJ. Quantifying mRNA in postmortem human brain: influence of gender, age at death, postmortem interval, brain pH, agonal state and inter-lobe mRNA variance. Brain Res Mol Brain Res. (2003) 118:60–71. doi: 10.1016/S0169-328X(03)00337-1
31. Birdsill, AC, Walker, DG, Lue, L, Sue, LI, and Beach, TG. Postmortem interval effect on RNA and gene expression in human brain tissue. Cell Tissue Bank. (2011) 12:311–8. doi: 10.1007/s10561-010-9210-8
32. Li, JZ, Vawter, MP, Walsh, DM, Tomita, H, Evans, SJ, Choudary, PV, et al. Systematic changes in gene expression in postmortem human brains associated with tissue pH and terminal medical conditions. Hum Mol Genet. (2004) 13:609–16. doi: 10.1093/hmg/ddh065
33. Iwamoto, K, Bundo, M, and Kato, T. Altered expression of mitochondria-related genes in postmortem brains of patients with bipolar disorder or schizophrenia, as revealed by large-scale DNA microarray analysis. Hum Mol Genet. (2005) 14:241–53. doi: 10.1093/hmg/ddi022
34. Atz, M, Walsh, D, Cartagena, P, Li, J, Evans, S, Choudary, P, et al. Methodological considerations for gene expression profiling of human brain. J Neurosci Methods. (2007) 163:295–309. doi: 10.1016/j.jneumeth.2007.03.022
35. Hoffman, GE, Jaffe, AE, Gandal, MJ, Collado-Torres, L, Sieberts, SK, Devlin, B, et al. Comment on: what genes are differentially expressed in individuals with schizophrenia? A Syst Rev Mol Psychiatry. (2023) 28:523–5. doi: 10.1038/s41380-022-01781-7
36. Xu, W, Seok, J, Mindrinos, MN, Schweitzer, AC, Jiang, H, Wilhelmy, J, et al. Human transcriptome array for high-throughput clinical studies. Proc Natl Acad Sci U S A. (2011) 108:3707–12. doi: 10.1073/pnas.1019753108
37. Kunii, Y, Matsumoto, J, Izumi, R, Nagaoka, A, Hino, M, Shishido, R, et al. Evidence for altered phosphoinositide signaling-associated molecules in the postmortem prefrontal cortex of patients with schizophrenia. Int J Mol Sci. (2021) 22:8280. doi: 10.3390/ijms22158280
38. Ono, CT, Yu, Z, Kikuchi, Y, Kunii, Y, Hino, M, Matsumoto, J, et al. Minimal amount of tissue-based pH measurement to improve quality control in neuropsychiatric post-mortem brain studies. Psychiatry Clin Neurosci. (2019) 73:566–73. doi: 10.1111/pcn.12863
39. Ono, C, Yu, Z, Kasahara, Y, Kikuchi, Y, Ishii, N, and Tomita, H. Fluorescently activated cell sorting followed by microarray profiling of helper T cell subtypes from human peripheral blood. PLoS One. (2014) 9:e111405. doi: 10.1371/journal.pone.0111405
40. Yu, Z, Ono, C, Kim, HB, Komatsu, H, Tanabe, Y, Sakae, N, et al. Four mood stabilizers commonly induce FEZ1 expression in human astrocytes. Bipolar Disord. (2011) 13:486–99. doi: 10.1111/j.1399-5618.2011.00946.x
41. Yu, Z, Sakai, M, Fukushima, H, Ono, C, Kikuchi, Y, Koyama, R, et al. Contextual fear conditioning regulates synapse-related gene transcription in mouse microglia. Brain Res Bull. (2022) 189:57–68. doi: 10.1016/j.brainresbull.2022.08.017
42. Jassal, B, Matthews, L, Viteri, G, Gong, C, Lorente, P, Fabregat, A, et al. The reactome pathway knowledgebase. Nucleic Acids Res. (2020) 48:D498–503. doi: 10.1093/nar/gkz1031
43. Shannon, PT, Grimes, M, Kutlu, B, Bot, JJ, and Galas, DJ. RCytoscape: tools for exploratory network analysis. BMC Bioinform. (2013) 14:217. doi: 10.1186/1471-2105-14-217
44. Marchiq, I, and Pouyssegur, J. Hypoxia, cancer metabolism and the therapeutic benefit of targeting lactate/H(+) symporters. J Mol Med (Berl). (2016) 94:155–71. doi: 10.1007/s00109-015-1307-x
45. Cummins, EP, Berra, E, Comerford, KM, Ginouves, A, Fitzgerald, KT, Seeballuck, F, et al. Prolyl hydroxylase-1 negatively regulates Ikappa B kinase-beta, giving insight into hypoxia-induced NFkappaB activity. Proc Natl Acad Sci U S A. (2006) 103:18154–9. doi: 10.1073/pnas.0602235103
46. Hayashi, Y, Yokota, A, Harada, H, and Huang, G. Hypoxia/pseudohypoxia-mediated activation of hypoxia-inducible factor-1alpha in cancer. Cancer Sci. (2019) 110:1510–7. doi: 10.1111/cas.13990
47. Wang, J, Xu, L, Yun, X, Yang, K, Liao, D, Tian, L, et al. Proteomic analysis reveals that proteasome subunit beta 6 is involved in hypoxia-induced pulmonary vascular remodeling in rats. PLoS One. (2013) 8:e67942. doi: 10.1371/journal.pone.0067942
48. Aydo Gan, TS, Dayi, G, and Kockar, F. Upregulation of PSMD4 gene by hypoxia in prostate cancer cells. Turk J Biol. (2020) 44:275–83. doi: 10.3906/biy-2002-71
49. Azzouzi, HE, Leptidis, S, Doevendans, PA, and De Windt, LJ. Hypoxami Rs: regulators of cardiac hypoxia and energy metabolism. Trends Endocrinol Metab. (2015) 26:502–8. doi: 10.1016/j.tem.2015.06.008
50. Nallamshetty, S, Chan, SY, and Loscalzo, J. Hypoxia: a master regulator of micro RNA biogenesis and activity. Free Radic Biol Med. (2013) 64:20–30. doi: 10.1016/j.freeradbiomed.2013.05.022
51. Aschrafi, A, Kar, AN, Natera-Naranjo, O, Mac Gibeny, MA, Gioio, AE, and Kaplan, BB. Micro RNA-338 regulates the axonal expression of multiple nuclear-encoded mitochondrial mRNAs encoding subunits of the oxidative phosphorylation machinery. Cell Mol Life Sci. (2012) 69:4017–27. doi: 10.1007/s00018-012-1064-8
52. Garchow, B, Maque Acosta, Y, and Kiriakidou, M. HIF-1alpha and mi R-210 differential and lineage-specific expression in systemic lupus erythematosus. Mol Immunol. (2021) 133:128–34. doi: 10.1016/j.molimm.2021.02.019
53. Kondo, A, Yamamoto, S, Nakaki, R, Shimamura, T, Hamakubo, T, Sakai, J, et al. Extracellular acidic pH activates the sterol regulatory element-binding protein 2 to promote tumor progression. Cell Rep. (2017) 18:2228–42. doi: 10.1016/j.celrep.2017.02.006
54. Adzigbli, L, Sokolov, EP, Wimmers, K, Sokolova, IM, and Ponsuksili, S. Effects of hypoxia and reoxygenation on mitochondrial functions and transcriptional profiles of isolated brain and muscle porcine cells. Sci Rep. (2022) 12:19881. doi: 10.1038/s41598-022-24386-0
55. Hirai, S, Miwa, H, Tanaka, T, Toriumi, K, Kunii, Y, Shimbo, H, et al. High-sucrose diets contribute to brain angiopathy with impaired glucose uptake and psychosis-related higher brain dysfunctions in mice. Sci Adv. (2021) 7:eabl6077. doi: 10.1126/sciadv.abl6077
56. Ben-Shachar, D. Mitochondrial multifaceted dysfunction in schizophrenia; complex I as a possible pathological target. Schizophr Res. (2017) 187:3–10. doi: 10.1016/j.schres.2016.10.022
57. Prabakaran, S, Swatton, JE, Ryan, MM, Huffaker, SJ, Huang, JT, Griffin, JL, et al. Mitochondrial dysfunction in schizophrenia: evidence for compromised brain metabolism and oxidative stress. Mol Psychiatry. (2004) 9:684–97. doi: 10.1038/sj.mp.4001511
58. Pruett, BS, and Meador-Woodruff, JH. Evidence for altered energy metabolism, increased lactate, and decreased pH in schizophrenia brain: a focused review and meta-analysis of human postmortem and magnetic resonance spectroscopy studies. Schizophr Res. (2020) 223:29–42. doi: 10.1016/j.schres.2020.09.003
59. Rius, J, Guma, M, Schachtrup, C, Akassoglou, K, Zinkernagel, AS, Nizet, V, et al. NF-kappa B links innate immunity to the hypoxic response through transcriptional regulation of HIF-1alpha. Nature. (2008) 453:807–11. doi: 10.1038/nature06905
60. Winning, S, and Fandrey, J. Oxygen sensing in innate immune cells: how inflammation broadens classical hypoxia-inducible factor regulation in myeloid cells. Antioxid Redox Signal. (2022) 37:956–71. doi: 10.1089/ars.2022.0004
61. Zhou, T, Mountz, JD, and Kimberly, RP. Immunobiology of tumor necrosis factor receptor superfamily. Immunol Res. (2002) 26:323–36. doi: 10.1385/IR:26:1-3:323
62. Roussos, P, Katsel, P, Davis, KL, Giakoumaki, SG, Lencz, T, Malhotra, AK, et al. Convergent findings for abnormalities of the NF-kappa B signaling pathway in schizophrenia. Neuropsychopharmacology. (2013) 38:533–9. doi: 10.1038/npp.2012.215
63. Murphy, CE, Kondo, Y, Walker, AK, Rothmond, DA, Matsumoto, M, and Shannon, WC. Regional, cellular and species difference of two key neuroinflammatory genes implicated in schizophrenia. Brain Behav Immun. (2020) 88:826–39. doi: 10.1016/j.bbi.2020.05.055
64. Murphy, CE, Lawther, AJ, Webster, MJ, Asai, M, Kondo, Y, Matsumoto, M, et al. Nuclear factor kappa B activation appears weaker in schizophrenia patients with high brain cytokines than in non-schizophrenic controls with high brain cytokines. J Neuroinflammation. (2020) 17:215. doi: 10.1186/s12974-020-01890-6
65. Jing, X, Jiang, T, Dai, L, Wang, X, Jia, L, Wang, H, et al. Hypoxia-induced autophagy activation through NF-kappa B pathway regulates cell proliferation and migration to induce pulmonary vascular remodeling. Exp Cell Res. (2018) 368:174–83. doi: 10.1016/j.yexcr.2018.04.026
66. Jeong, W, Jung, S, Bazer, FW, and Kim, J. Hypoxia-dependent accumulation of hypoxia-inducible factor-1 alpha induces transient cell cycle arrest in porcine trophectoderm cells. Theriogenology. (2018) 115:9–15. doi: 10.1016/j.theriogenology.2018.04.016
67. Cao, L, Wang, J, Bo, L, Li, DW, Yin, N, Zhou, AW, et al. Effects of hypoxia on the growth and development of the fetal ovine hepatocytes in primary culture. Biomed Environ Sci. (2019) 32:592–601. doi: 10.3967/bes2019.077
68. Silberberg, G, Baruch, K, and Navon, R. Detection of stable reference genes for real-time PCR analysis in schizophrenia and bipolar disorder. Anal Biochem. (2009) 391:91–7. doi: 10.1016/j.ab.2009.05.026
69. Wockner, LF, Noble, EP, Lawford, BR, Young, RM, Morris, CP, Whitehall, VL, et al. Genome-wide DNA methylation analysis of human brain tissue from schizophrenia patients. Transl Psychiatry. (2014) 4:e339. doi: 10.1038/tp.2013.111
70. Catts, VS, Lai, YL, Weickert, CS, Weickert, TW, and Catts, SV. A quantitative review of the postmortem evidence for decreased cortical N-methyl-D-aspartate receptor expression levels in schizophrenia: how can we link molecular abnormalities to mismatch negativity deficits? Biol Psychol. (2016) 116:57–67. doi: 10.1016/j.biopsycho.2015.10.013
71. Maycox, PR, Kelly, F, Taylor, A, Bates, S, Reid, J, Logendra, R, et al. Analysis of gene expression in two large schizophrenia cohorts identifies multiple changes associated with nerve terminal function. Mol Psychiatry. (2009) 14:1083–94. doi: 10.1038/mp.2009.18
72. Iwamoto, K, Kakiuchi, C, Bundo, M, Ikeda, K, and Kato, T. Molecular characterization of bipolar disorder by comparing gene expression profiles of postmortem brains of major mental disorders. Mol Psychiatry. (2004) 9:406–16. doi: 10.1038/sj.mp.4001437
73. Ohayon, S, Yitzhaky, A, and Hertzberg, L. Gene expression meta-analysis reveals the up-regulation of CREB1 and CREBBP in Brodmann area 10 of patients with schizophrenia. Psychiatry Res. (2020) 292:113311. doi: 10.1016/j.psychres.2020.113311
74. Egbujo, CN, Sinclair, D, and Hahn, CG. Dysregulations of synaptic vesicle trafficking in schizophrenia. Curr Psychiatry Rep. (2016) 18:77. doi: 10.1007/s11920-016-0710-5
75. Varea, E, Guirado, R, Gilabert-Juan, J, Marti, U, Castillo-Gomez, E, Blasco-Ibanez, JM, et al. Expression of PSA-NCAM and synaptic proteins in the amygdala of psychiatric disorder patients. J Psychiatr Res. (2012) 46:189–97. doi: 10.1016/j.jpsychires.2011.10.011
76. Shen, YC, Tsai, HM, Ruan, JW, Liao, YC, Chen, SF, and Chen, CH. Genetic and functional analyses of the gene encoding synaptophysin in schizophrenia. Schizophr Res. (2012) 137:14–9. doi: 10.1016/j.schres.2012.01.028
77. Honer, WGFP, Chen, C, Arango, V, Mann, JJ, and Dwork, AJ. Synaptic and plasticity-associated proteins in anterior frontal cortex in severe mental illness. Neuroscience. (1999) 91:1247–55. doi: 10.1016/S0306-4522(98)00679-4
78. Hagihara, H, Catts, VS, Katayama, Y, Shoji, H, Takagi, T, Huang, FL, et al. Decreased brain pH as a shared endophenotype of psychiatric disorders. Neuropsychopharmacology. (2018) 43:459–68. doi: 10.1038/npp.2017.167
Keywords: schizophrenia, postmortem brain, microarray, pH, RNA integrity number
Citation: Miyahara K, Hino M, Yu Z, Ono C, Nagaoka A, Hatano M, Shishido R, Yabe H, Tomita H and Kunii Y (2023) The influence of tissue pH and RNA integrity number on gene expression of human postmortem brain. Front. Psychiatry. 14:1156524. doi: 10.3389/fpsyt.2023.1156524
Edited by:
Ju Wang, Tianjin Medical University, ChinaReviewed by:
Yulin Kang, Chinese Research Academy of Environmental Sciences, ChinaQi Wang, Southern Medical University, China
Copyright © 2023 Miyahara, Hino, Yu, Ono, Nagaoka, Hatano, Shishido, Yabe, Tomita and Kunii. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yasuto Kunii, a3VuaWlAbWVkLnRvaG9rdS5hYy5qcA==