- 1Department of Psychiatry, Xijing 986 Hospital, Fourth Military Medical University, Xi’an, China
- 2Shaanxi Provincial Key Laboratory of Clinic Genetics, Fourth Military Medical University, Xi’an, China
- 3Honghui Hospital, Xi’an Jiaotong University, Xi’an, China
- 4Unit 95871 of Air Force, Hengyang, China
- 5Department of Respiratory and Critical Care Medicine, The First Affiliated Hospital of Chongqing Medical University, Chongqing, China
- 6Department of Respiratory Medicine, The 988th Hospital of the Joint Logistic Support Force of the People’s Liberation Army of China, Zhengzhou, China
- 7Department of Radiology, The Second Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China
- 8State Key Laboratory of Neurology and Oncology Drug Development, Xi’an, China
Aim: The adverse health outcomes of schizophrenia (SZ) are largely driven by the high prevalence of other non-neurological diseases. In addition to accelerated brain aging, patients with SZ also exhibit signs of systemic aging. However, the potential causal or biological mechanisms between multisystem aging and schizophrenia remain unknown.
Methods: We obtained SZ-associated single-nucleotide polymorphism (SNP) sets, aging gene data, and tissue-specific cis-expression quantitative trait locus (cis-eQTL) data of the cerebral cortex and other tissues from a previous two-stage genome-wide association study (GWAS), GeneCards database, and Genotype-Tissue Expression (GTEx) project. We employed tissue-specific Mendelian randomization (MR) analysis to elucidate the tissue-specific expression patterns of aging-related genes, and used the summary data-based MR (SMR) approach to obtain tissue aging-related genes associated with the risk of SZ development. We identified the potential aging-related pathways through which these tissue-specific cis-eQTLs may affect SZ using enrichment analyses. Finally, we explored the relationship between the identified crucial aging-related genes and predicted age difference (PAD) of the brain in our clinical patients.
Results: We found that the expression of tissue-specific aging genes, including synuclein alpha (SNCA), angiotensin I converting enzyme (ACE), BRCA1 DNA repair-associated (BRCA1), MutL homolog 1 (MLH1), vascular endothelial growth factor A (VEGFA), microtubule-associated protein tau (MAPT), and age-related maculopathy susceptibility 2 (ARMS2), may affect SZ. The tissue-specific cis-eQTL may influence SZ through aging pathways. The brain PAD was significantly higher in the high-expression group of BRCA1 than in the low-expression group.
Conclusions: This study provides valuable clues to understand the link between SZ and multiorgan system aging and improves the current understanding of multiple tissue-specific aging-related genes with SZ.
Introduction
Schizophrenia (SZ) is a highly heterogeneous disorder (1, 39). People with SZ show abnormalities in several organ systems in addition to the central nervous system (CNS) (2). To a large extent, the negative health outcomes for SZ are driven by the high rates of comorbid metabolic syndrome and related diseases (3). Previous studies have uncovered a shared genetic etiology among cardiovascular disease, frailty, and SZ, as well as altered oral microbiota and systemic immune dysfunction in patients with SZ (4–6). Whether SZ is a multisystem disorder or whether different mechanisms trigger its high rates of comorbidity but result from common risk factors still remains unknown.
Evidence suggests that aging plays an important role in SZ (7–10). Five case–control studies found that SZ is accompanied by accelerated biological aging by midlife (11). Previous studies have found that there is a common biological basis between SZ patients and normal elderly individuals with brain aging. In SZ and aging, astrocytes, glutamatergic, and GABAergic neurons show low synaptic neuron–astrocyte program expression, which is associated with cognitive flexibility and plasticity (12, 13). The cognitive impairment symptoms of SZ patients resemble those of the elderly, mainly involving decreased ability to process high-load information, episodic non-verbal memory impairment, slowed processing speed, and weakened motor coordination. These symptoms suggest that the pathological state of SZ patients is associated with accelerated brain aging (14–16). Compared with the general population, early SZ is associated not only with alterations in brain structure and function but also with multiple changes in the body (2). However, the association between the brain and body health, as well as the associated disease risk and physical multimorbidity across body systems, remains poorly characterized.
Mendelian randomization (MR) analysis is an emerging method that uses genetic variants as instrumental variables (IVs) to infer the causal effect of an exposure on an outcome (17). In order to be able to locate causality more precisely at the molecular level, we utilize the summary data-based Mendelian randomization (SMR), which could effectively integrate multisource data (18). Due to the specificity of IVs, the MR estimates are not commonly subject to confounding bias and reverse causation. MR has also been applied to detect putative causal effects of tissue-specific gene expression and a wide range of diseases using expression quantitative trait loci (eQTLs) as instruments (19, 20). Generally, comorbidity-related studies use SMR, which is essential for studying the causal relationships between different organ systems and diseases, as it helps to avoid confounding factors and establish more reliable causal links (21, 22).
The aim of this study was to investigate the causal effect of aging on SZ by using the MR method. To identify the potential target gene, the tissue-type-specific causal effects of aging on cognitive function were evaluated using cis-eQTL-based MR.
Materials and methods
Data acquisition
SZ-associated SNP sets were derived from a previous two-stage genome-wide association study (GWAS) (23). This is one of the largest available GWASs of SZ that report common variant associations at 287 distinct genomic loci, including up to 76,755 individuals with SZ and 243,649 control individuals. In the primary GWAS, they have analyzed up to 7,585,078 SNPs with MAF ≥1% in 175,799 individuals of whom 74.3% were European, 17.5% East Asian, 5.7% African-American, and 2.5% Latino. In the extended GWAS, they have meta-analyzed the primary GWAS results with summary statistics from deCODE Genetics (1,979 cases, 142,626 controls) for index SNPs with P <1 × 10−5 and identified 342 LD-independent significant SNPs located in 287 loci.
The tissue-specific cis-eQTL data of the brain cortex, brain hippocampus, brain hypothalamus, heart, liver, lung, kidney, pancreas, muscle, and adipose were obtained from the Genotype-Tissue Expression (GTEx) project (v8; https://gtexportal.org/home/). The GTEx Portal is a comprehensive public resource for researchers studying tissue- and cell-specific gene expression and regulation across individuals, development, and species, with data from three NIH projects. Ethical approval of all data was obtained in the original studies.
The aging genes were obtained from the GeneCards database (https://www.genecards.org/), which is a comprehensive and authoritative compendium of human gene information. We selected 50 genes to represent aging-related biology. These genes, which are often located at key nodes of aging molecular networks or broadly expressed in multiple aging-related tissues, correspond to the top 50 genes with the highest aging-related scores in the GeneCards database.
We recruited 43 patients with SZ from Xijing Hospital for brain MRI scanning to calculate brain age and collected peripheral whole blood samples from the patients to measure gene expression using RNA sequencing (RNA-seq) technology. The diagnosis of SZ was determined according to the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) and confirmed by two experienced clinical psychiatrists after a comprehensive assessment of all available information. Detailed inclusion and exclusion criteria were previously documented. RNA-seq data derived from peripheral blood samples were utilized in this study, as previously described in detail. Specifically, 2.5 mL of whole blood was collected into PAXgene Blood RNA Tubes and immediately stored at −80°C. Subsequent RNA sequencing was performed on the Illumina NovaSeq 6000 platform. Raw sequencing data underwent quality control using Fastp (v.0.18.0; fastp: an ultra-fast all-in-one FASTQ preprocessor; https://github.com/OpenGene/fastp), with low-quality reads being excluded. The cleaned reads were then aligned to the human reference genome (hg19) using HISAT2.2.4 (HISAT: a fast spliced aligner with low memory requirements; https://daehwankimlab.github.io/hisat2/). Finally, the raw count data were normalized using DESeq2 (moderated estimation of fold change and dispersion for RNA-seq data with DESeq2; http://www.bioconductor.org/packages/release/bioc/html/DESeq2.html). The study protocol was reviewed and approved by the Xijing Hospital Institutional Ethics Committee and conformed to the ethical standards for medical research involving human subjects, as laid out in the 1964 Declaration of Helsinki and its later amendments. Participants provided written informed consent prior to taking part in the study.
Brain age calculation
Imaging data acquisition and preprocessing
Patients and healthy controls underwent 3D high-resolution structural MRI scans, with raw images stored in DICOM format. Imaging data were acquired using a GE Discovery MR750 3.0T scanner and an 8-channel standard phased-array head coil (40).
Brain age model
Leveraging a 3D-convolutional neural network (3D-CNN) algorithm, this model accurately predicts brain age. Based on the classic VGGNet architecture, it is optimized into a simple fully convolutional neural network (SFCN) for efficient brain imaging data processing and analysis. Predicted age difference (PAD)=calculated age − chronological age. The specific establishment and training process of the brain age prediction model is detailed in the Supplementary Material.
MR analysis
We first assessed different tissue-dependent effects of aging gene expression on SZ through tissue-specific MR analysis and estimated the putative causal effects of the expression of 50 aging genes in 10 tissues based on the GTEx database. MR analysis of the expression of aging genes in eight tissues on SZ was then conducted. FDR correction for SMR P-value was applied using the Benjamini–Hochberg method. Tissue-dependent effects of the expression of aging genes on SZ are significant when FDR_SMR <0.05 and P_HEIDI >0.05.
Then, we applied a tissue MR analysis to determine all tissue-specific eQTLs that have causal effects on SZ. The SMR selected tissue-specific cis-eQTL significantly associated with SZ with a genome-wide threshold of 5 × 10−8. All specific MR analysis details and codes are provided in the Supplementary Material.
Enrichment analysis
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were used to detect whether the tissue-specific cis-eQTLs through the aging pathways influence SZ. The “clusterProfiler” package in R software (https://guangchuangyu.github.io/software/clusterProfiler/) was utilized for the GO and KEGG enrichment analyses of cis-eQTL, with a selection criterion of q-value <0.05.
Comparison of brain PAD between the high- and low-expression groups of aging genes
We divided the data of SZ patients into two groups, high- and low-expression, based on the top 25% and bottom 25% of tissue-specific senescence gene expression. Then, we analyzed their data on brain age to verify whether the screened genes significantly affect the brain PAD of SZ patients. A substantial brain age gap between the high- and low-expression groups points to this gene being a significant contributor to brain aging in schizophrenia patients.
Statistical analyses
Data are presented as mean ± standard deviation (SD). Comparison between the two groups was performed using the Student’s t-test and the Wilcoxon test. All statistical analyses were conducted using SPSS software (version 30.0.0). A P-value of less than 0.05 was considered statistically significant.
Results
Tissue-dependent effects of aging gene expression on SZ
The MR analyses suggested putative causal effects of seven aging gene expressions in seven tissues on SZ. Among these were ACE in the lung, VEGFA in the pancreas, MAPT in the spleen, and SNCA in the heart’s left ventricle (FDR_SMR < 0.05, P_HEIDI > 0.05) (Table 1). Additionally, we detected several genes with suggestive associations, such as MLH1 in the muscle and BRCA1 in the liver.
The influence of tissue-specific aging eQTLs on SZ
Among the seven tissues analyzed, we identified that the ACE gene in lung tissue (Figure 1A), the VEGFA gene in pancreatic tissue (Figure 1B), the MAPT gene in the spleen (Figure 1C), and the SNCA gene in the heart (Figure 1D) exhibited causal relationships with SZ.
Figure 1. The relationship between the effect size of eQTLs and GWAS of tissue-specific genes. eQTL, expression quantitative trait locus; GWAS, genome-wide association study. (A) The relationship between the effect size of eQTLs and GWAS of tissue-specific genes in the lung. (B) The relationship between the effect size of eQTLs and GWAS of tissue-specific genes in the pancreas. (C) The relationship between the effect size of eQTLs and GWAS of tissue-specific genes in the spleen. (D) The relationship between the effect size of eQTLs and GWAS of tissue-specific genes in the left ventricle.
Enrichment analysis of tissue-specific eQTLs
To further elucidate the connections between peripheral organs and SZ, we conducted enrichment analyses of tissue-specific genes associated with SZ (Supplementary File 1). As shown by the GO enrichment analysis, we found that organs like the liver, lungs, and heart may influence SZ through aging pathways. In the legend of Figure 2, GeneRatio was used to represent the proportion of tissue-specific genes annotated to different GO/KEGG pathways in an organ, while Padj was used to represent the enrichment significance adjusted by the Benjamini–Hochberg procedure. For instance, liver-related pathways involve antigen processing and presentation (e.g., peptide antigen via MHC class I), DNA strand elongation, and oligopeptide transport. In the lungs, the cAMP-responsive element binding protein (CREB) pathway, the isomerase activity pathway, and the four-way junction DNA binding pathway are relevant. The pancreas-specific genes were associated with ATPase complexes and SWI/SNF superfamily type complex pathways. KEGG pathway analysis revealed that liver-specific genes associated with SZ were mainly enriched in lysosomes and ABC transporter proteins, which have essential roles in aging-related pathways. The details are shown in Figure 2.
Figure 2. The GO and KEGG enrichment analyses of tissue-specific aging genes. GO, Gene Ontology. KEGG, Kyoto Encyclopedia of Genes and Genomes.
Validation of the effects of selected genes on brain PAD in SZ patients
The brain PAD was significantly higher in the high-expression group of BRCA1 than in the low-expression group (Figures 3A–G). Although statistically significant results were obtained only for the BRCA1 gene, the VEGFA and SNCA still showed a favorable trend, probably due to the limitation of a small sample size.
Figure 3. Analysis of the differences in brain PAD between the high- and low-expression groups of tissue-specific aging genes in patients with SZ. PAD, predicted age difference; SZ, schizophrenia. (A) Analysis of the differences in brain PAD between the high- and low-expression groups of ACE in patients with SZ. (B) Analysis of the differences in brain PAD between the high- and low-expression groups of VEGFA in patients with SZ. (C) Analysis of the differences in brain PAD between the high- and low-expression groups of MAPT in patients with SZ. (D) Analysis of the differences in brain PAD between the high- and low-expression groups of SNCA in patients with SZ. (E) Analysis of the differences in brain PAD between the high- and low-expression groups of MLH1 in patients with SZ. (F) Analysis of the differences in brain PAD between the high- and low-expression groups of BRCA1 in patients with SZ. (G) Analysis of the differences in brain PAD between the high- and low-expression groups of ARMS2 in patients with SZ.
Discussion
Previous studies have mostly focused on the interplay between SZ and brain aging or the impact of a single organ or system on SZ. No research has yet examined how the body’s multiple organ systems may influence the brain via aging pathways to affect SZ. In this study, we used genetic tools and found that organ-specific aging genes in peripheral organs—ACE in the lung, VEGFA in the pancreas, MAPT in the spleen, SNCA in the heart’s left ventricle, MLH1 in muscle, BRCA1 in the liver, and ARMS2 in adipose tissue—might influence SZ through aging-related pathways. Expression data of these genes in the whole blood and enrichment analysis results further confirmed this finding.
The SNCA gene encodes α-synuclein, which is involved in synaptic transmission and neurotransmitter release. α-Synuclein is highly expressed at the presynaptic terminal and is involved in a variety of cellular functions, including synaptic vesicle transport, membrane binding, and signal transduction, as well as affecting the homeostasis of the mitochondria and lysosomes (24). Joung et al. found that although the role of the SNCA gene is more pronounced in the brain, it also affects cardiac aging. The heart–brain axis is a bidirectional communication system that enables interaction between the heart and the brain through various pathways, including the autonomic nervous system and the immune system. Aging can lead to functional and structural changes in the heart, which in turn can affect brain function. A previous study found that cardiovascular diseases can result in reduced cerebral blood flow, potentially leading to cognitive decline symptoms observed in patients with SZ (25). Additionally, the SNCA gene may influence monocyte metabolism by participating in the regulation of the unfolded protein response and endoplasmic reticulum stress. Abnormal aggregation of α-synuclein may also activate monocytes and macrophages, prompting them to secrete pro-inflammatory cytokines and exacerbate systemic inflammation. These factors may all contribute to the development and progression of SZ. Our research also found that monocyte count genes are associated with cortical thickness (24, 26, 27).
The liver–brain axis and the muscle–brain axis are also bidirectional communication systems. Aging can lead to changes in liver function and structure, which may promote the development of liver and muscle diseases and affect brain function through the liver–brain and muscle–brain axes (28, 29). Although the primary functions of BRCA1 and MLH1 are related to DNA repair (30), they can still influence the function and structure of the liver and muscles, which in turn can affect brain function (31, 32). More importantly, in our clinical patients with SZ, the BRCA1 gene expression in the blood is associated with brain PAD, which further suggests the potential role of the aging-related gene BRCA1 in the liver–brain axis and SZ.
In lung–brain interaction, changes in ACE gene expression impact pulmonary blood flow and oxygenation, which in turn affect brain function. ACE inhibitors are widely used to treat hypertension and heart failure. Cao et al. found that the effects of ACE inhibitors on blood pressure can indirectly influence the lung–brain axis, thereby impacting brain function. Overexpression of ACE is associated with increased oxidative metabolism and heightened immune responses. This can affect the energy production and immune function of monocytes, leading to inflammation and subsequent brain changes, which is consistent with our study (27, 33).
ARMS2 is a gene linked to age-related macular degeneration. It encodes a secreted protein functioning in the extracellular matrix, likely maintaining matrix homeostasis through interactions with various matrix proteins (34). Adipose tissue undergoes significant changes during aging, including adipocyte volume reduction, decreased lipid storage capacity, and alterations in adipokine secretion levels (35). Adipose tissue may upregulate ARMS2 to regulate and maintain its structure and function. ARMS2 may interact with matrix proteins to influence monocyte migration, phagocytosis, and immune responses. Its mitochondria-associated functions could indirectly affect monocyte energy metabolism and activity (36), potentially causing chronic inflammation that impacts brain function (27).
The VEGFA gene is crucial for angiogenesis, inflammation, and oxidative stress. As pancreatic tissue ages, vascular regression and reduced blood supply occur. By promoting angiogenesis, VEGFA may improve pancreatic blood supply and slow its aging. VEGFA might also alter chromatin accessibility via the SWI/SNF complex, regulating genes linked to neuronal survival, differentiation, and synaptic plasticity. Moreover, VEGFA overexpression may activate pathways that influence monocyte and macrophage function, potentially leading to a pro-inflammatory microenvironment (27, 37, 38).
However, one limitation of the present study is the relatively modest sample size (n=43) of the clinical cohort employed for the validation analyses examining the relationship between peripheral aging-related gene expression and brain PAD in schizophrenia patients. This limited sample size may have constrained statistical power, and future investigations incorporating multicenter collaborations and expanded sample sizes are essential to validate our preliminary observations.
In conclusion, our study offers novel insights into the relationships between SZ and the aging-related genes SNCA, ACE, BRCA1, MLH1, VEGFA, MAPT, and ARMS2 in multiple organs. Notably, the BRCA1 gene may be associated with accelerated brain aging in individuals with SZ. Our findings provide valuable clues for understanding the link between peripheral organ aging and SZ.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by Xijing Hospital Institutional Ethics Committee. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
Y-KH: Data curation, Investigation, Writing – original draft. M-YL: Data curation, Resources, Writing – original draft. D-LY: Software, Visualization, Writing – review & editing. J-XX: Data curation, Investigation, Writing – original draft. X-HW: Data curation, Investigation, Writing – original draft. D-BW: Data curation, Writing – original draft. Y-LL: Investigation, Writing – original draft. C-CW: Conceptualization, Methodology, Resources, Supervision, Writing – review & editing. L-BC: Writing – review & editing, Methodology. Y-JC: Conceptualization, Writing – review & editing. H-JZ: Conceptualization, Resources, Supervision, Writing – review & editing, Software.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This study was supported by the “Sanxin” talent projects Xijing 986 Hospital Department, Air Force Medical University (986SX24XMPY03), the State Key Laboratory of Neurology and Oncology Drug Development Open Project (SKLSIM-F-2024-66), New Aviator Program of Air Force Medical University, National Natural Science Foundation of China (82572176, 82271949), Natural Science Basic Research 14 Program of Shaanxi (2025SYS-SYSZD-061) and Inner Mongolia Natural Science Foundation (2025MS08146).
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
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.2025.1730143/full#supplementary-material
References
1. Smeland OB, Frei O, Dale AM, and Andreassen OA. The polygenic architecture of schizophrenia - rethinking pathogenesis and nosology. Nat Rev Neurol. (2020) 16:366–79. doi: 10.1038/s41582-020-0364-0
2. Pillinger T, D’Ambrosio E, McCutcheon R, and Howes OD. Is psychosis a multisystem disorder? A meta-review of central nervous system, immune, cardiometabolic, and endocrine alterations in first-episode psychosis and perspective on potential models. Mol Psychiatry. (2019) 24:776–94. doi: 10.1038/s41380-018-0058-9
3. Taube C, Mentzel C, Glue P, and Barak Y. Aging of persons with schizophrenia: analysis of a national dataset. Int psychogeriatrics. (2024) 36:43–50. doi: 10.1017/S1041610223000145
4. Rodevand L, Rahman Z, Hindley GFL, Smeland OB, Frei O, Tekin TF, et al. Characterizing the shared genetic underpinnings of schizophrenia and cardiovascular disease risk factors. Am J Psychiatry. (2023) 180:815–26. doi: 10.1176/appi.ajp.20220660
5. Liu X, Ling Z, Cheng Y, Wu L, Shao L, Gao J, et al. Oral fungal dysbiosis and systemic immune dysfunction in Chinese patients with schizophrenia. Transl Psychiatry. (2024) 14:475. doi: 10.1038/s41398-024-03183-5
6. Deng MG, Wang K, Liu F, Zhou X, Nie JQ, Zhao ZH, et al. Shared genetic architecture and causal relationship between frailty and schizophrenia. Schizophr (Heidelb). (2025) 11:24. doi: 10.1038/s41537-024-00550-5
7. Gerretsen P, Plitman E, Rajji TK, and Graff-Guerrero A. The effects of aging on insight into illness in schizophrenia: a review. Int J geriatric Psychiatry. (2014) 29:1145–61. doi: 10.1002/gps.4154
8. Parlikar R and Shivakumar V. Accelerated aging and frailty in schizophrenia. Int psychogeriatrics. (2024) 36:18–20. doi: 10.1017/S1041610223000479
9. Okusaga OO. Accelerated aging in schizophrenia patients: the potential role of oxidative stress. Aging Dis. (2014) 5:256–62. doi: 10.14336/AD.2014.0500256
10. Seeman MV. Subjective overview of accelerated aging in schizophrenia. Int J Environ Res Public Health. (2022) 20. doi: 10.3390/ijerph20010737
11. Caspi A, Shireby G, Mill J, Moffitt TE, Sugden K, and Hannon E. Accelerated pace of aging in schizophrenia: five case-control studies. Biol Psychiatry. (2024) 95:1038–47. doi: 10.1016/j.biopsych.2023.10.023
12. Ling E, Nemesh J, Goldman M, Kamitaki N, Reed N, Handsaker RE, et al. A concerted neuron-astrocyte program declines in ageing and schizophrenia. Nature. (2024) 2024:627. doi: 10.1038/s41586-024-07109-5
13. Hua JPY, Abram SV, Loewy RL, Stuart B, Fryer SL, Vinogradov S, et al. Brain age gap in early illness schizophrenia and the clinical high-risk syndrome: associations with experiential negative symptoms and conversion to psychosis. Schizophr Bull. (2024) 50:1159–70. doi: 10.1093/schbul/sbae074
14. Constantinides C, Han LKM, Alloza C, Antonucci LA, Arango C, Ayesa-Arriola R, et al. Brain ageing in schizophrenia: evidence from 26 international cohorts via the ENIGMA Schizophrenia consortium. Mol Psychiatry. (2023) 28:1201–9. doi: 10.1038/s41380-022-01897-w
15. Kirkpatrick B, Messias E, Harvey PD, Fernandez-Egea E, and Bowie CR. Is schizophrenia a syndrome of accelerated aging? Schizophr Bull. (2008) 34:1024–32. doi: 10.1093/schbul/sbm140
16. Nguyen TT, Eyler LT, and Jeste DV. Systemic biomarkers of accelerated aging in schizophrenia: A critical review and future directions. Schizophr Bull. (2018) 44:398–408. doi: 10.1093/schbul/sbx069
17. Evans DM and Davey Smith G. Mendelian randomization: new applications in the coming age of hypothesis-free causality. Annu Rev Genomics Hum Genet. (2015) 16:327–50. doi: 10.1146/annurev-genom-090314-050016
18. Dong XX, Chen DL, Chen HM, Li DL, Hu DN, Lanca C, et al. DNA methylation biomarkers and myopia: a multi-omics study integrating GWAS, mQTL and eQTL data. Clin Epigenet. (2024) 16:157. doi: 10.1186/s13148-024-01772-1
19. Qiu X, Guo R, Wang Y, Zheng S, Wang B, and Gong Y. Mendelian randomization reveals potential causal relationships between cellular senescence-related genes and multiple cancer risks. Commun Biol. (2024) 7:1069. doi: 10.1038/s42003-024-06755-9
20. Bryois J, Calini D, Macnair W, Foo L, Urich E, Ortmann W, et al. Cell-type-specific cis-eQTLs in eight human brain cell types identify novel risk genes for psychiatric and neurological disorders. Nat Neurosci. (2022) 25:1104–12. doi: 10.1038/s41593-022-01128-z
21. Shi X, Li M, Yao J, Li MD, and Yang Z. Alcohol drinking, DNA methylation and psychiatric disorders: A multi-omics Mendelian randomization study to investigate causal pathways. Addict (Abingdon England). (2024) 119:1226–37. doi: 10.1111/add.16465
22. Wang D, Li D, Dang X, Mu C, Liu C, Zeng Y, et al. Mendelian randomization reveals causalities between DNA methylation and schizophrenia. Biol Psychiatry. (2025). doi: 10.1016/j.biopsych.2025.03.012
23. Trubetskoy V, Pardinas AF, Qi T, Panagiotaropoulou G, Awasthi S, Bigdeli TB, et al. Mapping genomic loci implicates genes and synaptic biology in schizophrenia. Nature. (2022) 604:502–8. doi: 10.1038/s41586-022-04434-5
24. Joung J, Heo Y, Kim Y, Kim J, Choi H, Jeon T, et al. Cell enlargement modulated by GATA4 and YAP instructs the senescence-associated secretory phenotype. Nat Commun. (2025) 16:1696. doi: 10.1038/s41467-025-56929-0
25. Shen J and Jiang C. Unraveling the heart-brain axis: shared genetic mechanisms in cardiovascular diseases and Schizophrenia. Schizophr (Heidelberg Germany). (2024) 10:113. doi: 10.1038/s41537-024-00533-6
26. Soda T, Pasqua T, De Sarro G, and Moccia F. Cognitive impairment and synaptic dysfunction in cardiovascular disorders: the new frontiers of the heart-brain axis. Biomedicines. (2024) 12. doi: 10.3390/biomedicines12102387
27. Cui LB, Wang XY, Fu YF, Liu XF, Wei Y, Zhao SW, et al. Transcriptional level of inflammation markers associates with short-term brain structural changes in first-episode schizophrenia. BMC Med. (2023) 21:250. doi: 10.1186/s12916-023-02963-y
28. Duan L, Li H, Li S, Shi Y, and Feng Y. Causal association between sarcopenia and cognitive impairment contributes to the muscle-brain axis: A bidirectional Mendelian randomization study. Geriatrics gerontol Int. (2025) 25:116–22. doi: 10.1111/ggi.15045
29. Cutuli D, Decandia D, Giacovazzo G, and Coccurello R. Physical exercise as disease-modifying alternative against alzheimer’s disease: A gut-muscle-brain partnership. Int J Mol Sci. (2023) 24. doi: 10.3390/ijms241914686
30. Tomer G, Buermeyer AB, Nguyen MM, and Liskay RM. Contribution of human mlh1 and pms2 ATPase activities to DNA mismatch repair. J Biol Chem. (2002) 277:21801–9. doi: 10.1074/jbc.M111342200
31. Wang Z, Yin Z, Sun G, Zhang D, and Zhang J. Genetic evidence for the liver-brain axis: lipid metabolism and neurodegenerative disease risk. Lipids Health Dis. (2025) 24:41. doi: 10.1186/s12944-025-02455-3
32. Wei C, Li X, Jin Y, Zhang Y, and Yuan Q. Does the liver facilitate aging-related cognitive impairment: Conversation between liver and brain during exercise? J Cell Physiol. (2024) 239:e31287. doi: 10.1002/jcp.31287
33. Cao DY, Spivia WR, Veiras LC, Khan Z, Peng Z, Jones AE, et al. ACE overexpression in myeloid cells increases oxidative metabolism and cellular ATP. J Biol Chem. (2020) 295:1369–84. doi: 10.1016/S0021-9258(17)49895-4
34. Kortvely E, Hauck SM, Duetsch G, Gloeckner CJ, Kremmer E, Alge-Priglinger CS, et al. ARMS2 is a constituent of the extracellular matrix providing a link between familial and sporadic age-related macular degenerations. Invest Ophthalmol Visual Sci. (2010) 51:79–88. doi: 10.1167/iovs.09-3850
35. Nguyen TT and Corvera S. Adipose tissue as a linchpin of organismal ageing. Nat Metab. (2024) 6:793–807. doi: 10.1038/s42255-024-01046-3
36. Fritsche LG, Loenhardt T, Janssen A, Fisher SA, Rivera A, Keilhauer CN, et al. Age-related macular degeneration is associated with an unstable ARMS2 (LOC387715) mRNA. Nat Genet. (2008) 40:892–6. doi: 10.1038/ng.170
37. Failla CM, Carbone ML, Ramondino C, Bruni E, and Orecchia A. Vascular endothelial growth factor (VEGF) family and the immune system: activators or inhibitors? Biomedicines. (2024) 13. doi: 10.3390/biomedicines13010006
38. Torres-Espin A, Radabaugh HL, Treiman S, Fitzsimons SS, Harvey D, Chou A, et al. Sexually dimorphic differences in angiogenesis markers are associated with brain aging trajectories in humans. Sci Trans Med. (2024) 16:eadk3118. doi: 10.1126/scitranslmed.adk3118
39. Zhang X, Hatoum L, Ying J, and Huang C. Assessing glymphatic-associated fluid dynamics in psychiatric disorders: evidence from neuroimaging - a review. Psychoradiology. (2025) 5. doi: 10.1093/psyrad/kkaf031
Keywords: aging, brain neuroimaging, brain-age, genome-wide association study, schizophrenia, transcriptome
Citation: Han Y-K, Liu M-Y, Yang D-l, Xie J-X, Wang X-H, Wang D-B, Liang Y-L, Wang C-C, Cui L-B, Chen Y-J and Zhang H-J (2026) Brain imaging data and summary data-based Mendelian randomization analysis reveal the impact of multiorgan aging on schizophrenia. Front. Psychiatry 16:1730143. doi: 10.3389/fpsyt.2025.1730143
Received: 22 October 2025; Accepted: 31 December 2025; Revised: 30 December 2025;
Published: 02 February 2026.
Edited by:
Guglielmo Lucchese, Universitätsmedizin Greifswald, GermanyReviewed by:
Shiwu Li, Chinese Academy of Sciences (CAS), ChinaXingzhong Zhao, Fudan University, China
Copyright © 2026 Han, Liu, Yang, Xie, Wang, Wang, Liang, Wang, Cui, Chen and Zhang. 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: Hai-Jun Zhang, emhqYmVpamluZzIwMDhAMTYzLmNvbQ==; Yu-Jing Chen, YzE4MDAzNDA5NDAyQDE2My5jb20=
†These authors have contributed equally to this work
Cui-Cui Wang1