Glyoxalase 1 Confers Susceptibility to Schizophrenia: From Genetic Variants to Phenotypes of Neural Function

Abstract

This research aimed to investigate the cause of glyoxalase 1 (Glo-1) to the susceptibility of schizophrenia. Here, significant differences in Glo-1 mRNA expression (P = 3.98 E-5) and enzymatic activity (P =1.40E-6) were found in peripheral blood of first-onset antipsychotic-naïve schizophrenia patients and controls, the following receiver operating characteristic (ROC) curves analysis showed that Glo-1 could be able to predict the schizophrenia risk (P = 4.75 E-6 in mRNA, P = 1.43 E-7 in enzymatic activity, respectively). To identify the genetic source of Glo-1 risk in schizophrenia, polymorphisms of Glo-1 were genotyped by SNaPshot technology in 1069 schizophrenia patients and 1023 healthy individuals. Significant differences were found in distributions of genotype (P = 0.020) and allele (P = 0.020) in rs1781735, in which G >T mutation showed significantly reduction in promoter activity of Glo-1 (P = 0.016). In peripheral blood of first-onset antipsychotic-naïve schizophrenia patients and healthy controls, TT genotype was associated with lower mRNA expression and enzymatic activity compared to GG genotype (P < 0.05). Then basing on the expression quantitative trait locus (eQTL) and functional magnetic resonance imaging (fMRI) analysis, TT genotype of rs1781735, associated with lower RNA expression in brain, showed decreased neuronal activation in left middle frontal gyrus in schizophrenia (P < 0.001). This study for the first time investigated the Glo-1 effect in the schizophrenia risk on multiple levels, from genetic variants, transcription, expression, protein function, all the way to phenotypes.

Introduction

Schizophrenia is a devastating and complex neurodevelopmental disorder with highly heterogeneous and multifaceted clinical manifestations (Cannon, 1996; Allen et al., 2008). It has been postulated that a combination of multiple environmental and genetic predisposing factors may confer risk for schizophrenia (Yin et al., 2014; Modai and Shomron, 2016). The variants in risk genes might lead to dysregulation of molecular pathways and aberrations in the structure and function of the brain, which could be reflected by the abnormalities in cognitive function, emotions, and behaviors (Owen et al., 2016; Klingler et al., 2021). To systematically understand how genes confer risks across these levels of schizophrenia will help to build a spectrum of abnormal states, predict vulnerable processes, design personalized diagnostic, and therapeutic tools (Klingler et al., 2021).

Carbonyl stress is an abnormal metabolic state that results in the accumulation of dicarbonyl compounds, such as methylglyoxal (MG) and glyoxal (Miyata et al., 1999) leading to the eventual formation of advanced glycation end products (AGEs). The accumulation of MG and AGEs ultimately induces severe cytotoxicity in the neurons and causes dysfunction of the nervous system. The recent studies in animals and humans have demonstrated that enhanced carbonyl stress contributes to the onset of various neurological disorders, such as schizophrenia (Itokawa et al., 2018; Yoshioka et al., 2021), Parkinson's disease (Kurz et al., 2011), Alzheimer's disease (Kuhla et al., 2007; Frandsen et al., 2020; Lv et al., 2020), mood disorder (Fujimoto et al., 2008), autism (Junaid et al., 2004; Gabriele et al., 2014; Kovac et al., 2015), panic disorder (Politi et al., 2006), and anxiety behavior (Hovatta et al., 2005; Distler et al., 2012; Du et al., 2019).

Glyoxalase 1 (Glo-1), a rate-limiting enzyme in the glyoxalase system, metabolizes the MG into S-D-lactoylglutathione and effectively decreases the formation of AGEs (He et al., 2020; Sarker et al., 2020). In contrast, the impairment of Glo-1 leads to carbonyl stress and, eventually, the accumulation of AGEs in the brain (Wang et al., 2019; Kold-Christensen and Johannsen, 2020). Therefore, Glo-1 might play a key role in the balance of cellular detoxification and toxification in reactive carbonyl species (Hara et al., 2021). This process has been reported as a possible causative factor in the etiology of schizophrenia (Mizutani et al., 2019; Ohnishi et al., 2019).

The Glo-1 gene is located at human chromosome 6p21.2, which has previously been reported as an associated region for schizophrenia in multicenter genome-wide association studies (GWAS) conducted on large samples (International Schizophrenia et al., 2009; Shi et al., 2009; Stefansson et al., 2009; Schizophrenia Working Group of the Psychiatric Genomics, 2014). In particular, Glo-1 has been identified as a target gene of miR-137 (Lv et al., 2018), which is strongly associated with schizophrenia, as reported by GWAS (Schizophrenia Psychiatric Genome-Wide Association Study (GWAS) Consortium, 2011; Ma et al., 2014; Wang et al., 2020). Reduced expression of Glo-1 and, thus, increased accumulation of glycated dicarbonyl compounds have been reported in a subclass of schizophrenia (Kouidrat et al., 2013; Katsuta et al., 2014). The Glo-1 mutation caused structural alteration of the neurites in the post-mortem human cerebral tissues of patients with schizophrenia (Mizutani et al., 2019). Furthermore, a missense substitution (Glu111Ala) and heterozygous frameshift mutations were recently shown to significantly decrease the Glo-1 protein expression and enzymatic activity in the red blood cells (RBCs) from Japanese patients with schizophrenia (Arai et al., 2010; Miyashita et al., 2014a,b; Torii et al., 2020).

Given the critical role of Glo-1 in the detoxification of carbonyl stress, we sought to test the hypothesis that Glo-1 is involved in the pathogenesis of schizophrenia. The case-control studies were performed to identify Glo-1 genetic variants, mRNA expression, and enzymatic activity in the Han population of China. To examine the effects of Glo-1 on phenotypes of brain function in schizophrenia, an effective proxy of Glo-1 expression was selected, and neural imaging analysis was performed using the functional magnetic resonance imaging (fMRI) combined with the Glo-1 brain expression pattern. Our investigation improved the understanding of the effect of Glo-1 on schizophrenia and provides a referable analysis workflow in imaging genetics research.

Subjects And Methods

Subjects and Clinical Assessment

This hospital-based case-control study consecutively recruited 1,069 patients with schizophrenia (684 men and 385 women) from the Department of Psychiatry at the Affiliated Hospital of Guangdong Medical University, Guangdong Province, China. All the subjects were Han Chinese from the local area of southern China, received a structured clinical interview and were independently diagnosed by the two experienced senior psychiatrists according to the criteria of the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5). The patients were excluded from the study if they had any physical disease, especially renal dysfunction or diabetes mellitus (because these diseases may affect Glo-1 expression). The Positive and Negative Syndrome Scale (PANSS) was used to evaluate the psychotic symptoms in the patients with schizophrenia. The control group consisted of 1,023 healthy individuals (619 men and 404 women) who were recruited from the same geographical area as the patients and had no known history of major psychiatric disorders, serious physical disease, or substance abuse. For the imaging-genetics subgroup, 91 first-onset antipsychotic-naïve patients with schizophrenia who completed all the imaging and genetics protocols were recruited. All the participants were right-handed. For the biochemical analysis, first-onset antipsychotic-naïve patients with schizophrenia and age- and sex-matched controls were recruited. There were 63 patients and 53 controls for mRNA expression analysis and 74 patients and 73 controls for enzymatic activity analysis. The study protocol was reviewed and approved by the ethics committee of the Affiliated Hospital of Guangdong Medical University (IRB number, PJ2017058). All the subjects or their relatives provided informed consent to participate in the study.

Measurement of Glo-1 Enzymatic Activity and mRNA Expression

For mRNA analysis, fresh blood samples were obtained, and total RNA was extracted from the peripheral blood using the PaxGene Blood RNA Kit (Sangon, Shanghai, China). A real-time PCR was performed using the SYBR Premix Ex Taq kit (Takara, Dalian, China) according to the instructions from the manufacturer using the primers listed in Supplementary Table 1. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as the endogenous control, and the comparative cycle threshold (Ct) method was used to quantify the Glo-1 mRNA levels. Each reaction was performed three times. For enzymatic activity analysis, the RBCs were separated from the venous blood samples by centrifugation (2,000 g, 5 min). The Glo-1 enzymatic activity was determined with a spectrophotometer at 240 nm by measuring the formation of S-D-lactoylglutathione as reported in a previous study by Gabriele et al. (2014).

Single Nucleotide Polymorphism (SNP) Selection

The criteria for identifying SNPs included a linkage disequilibrium (LD) r² value >0.8 and a minor allele frequency (MAF) of more than 0.10 in the Han population of China. Five candidate SNPs (rs1781735, rs1049346, rs1130534, rs4746, and rs9470916) were selected based on the data from the 1,000 Genomes Project and the most reported SNPs (Gale et al., 2004b; Barua et al., 2011; Peculis et al., 2013; Gabriele et al., 2014; Tao et al., 2016). Therefore, these loci were used for the subsequent analysis (detailed information is shown in Figure 1A). The SNP Rs4746 (also called rs2736654), in exon 4 of the Glo-1 gene, is a C > A change that causes Ala111Glu and is associated with a decrease in Glo-1 enzymatic activity (Barua et al., 2011). In addition, T > A in rs1130534 causes a synonymous substitution at codon 124 (GGA to GGT and Gly to Gly). It was reported that the reduced MG concentrations in human whole blood cell lysates correlated with the A allele (Peculis et al., 2013), which may represent a marker for susceptibility to Glo-1-related diseases (Gabriele et al., 2014; Tao et al., 2016). The SNP rs9470916 is located in the 3′-untranslated region of the Glo-1 gene. The variant rs1781735 is located in the promoter and has complete linkage with rs1049346 (c.-7C>T) which exerts significant effects on Glo-1 transcriptional activity (Gale et al., 2004a). Therefore, we genotyped the rs1781735 SNP, which can serve as a proxy for rs1049346.

FIGURE 1

Figure 1. (A) Distribution of single nucleotide polymorphisms (SNPs) across the Glo-1 gene selected for the association analysis. Exons are shown as black boxes (exons), and intron sequences are shown as lines. The gray boxes represent the 3' untranslated region (3'-UTR) and 5' untranslated region (5'-UTR). The vertical lines indicate the locations of these polymorphisms. Five SNPs were identified in the Glo-1 gene from the HapMap/Haploview database and the literature; (B) transcripts per kilobase million (TPM) values of the Glo-1 expression pattern in the brain tissues and whole blood in healthy participants. (C) Effect of the rs1781735 T/G polymorphism on the transcriptional activity of the Glo-1 promoter. *, P < 0.05. The error bar in the bar plot represents the SD of the mean. (D) GTEx Multitissue expression quantitative trait locus (eQTL) Plot for rs1781735. NES, normalized effect size. P-value, from a t-test that compares the observed NES from single-tissue eQTL analysis to a null NES of 0; m-value, the posterior probability that an eQTL effect exists in each tissue tested in the cross-tissue meta-analysis. Small m-value (e.g., < 0.1): the study is predicted to not have an effect. Large m-value (e.g., >0.9): the study is predicted to have an effect.

Genotyping

A human genomic DNA was isolated from the EDTA-treated peripheral blood samples using a commercially available kit (Tiangen Biotech, Beijing, China) according to the instructions from the manufacturer. The selected Glo-1 SNPs were genotyped using the SNaPshot technique as described in the previous studies (Ma et al., 2014; Fu et al., 2017), and the primers are listed in Supplementary Table 2.

Luciferase Reporter Assay

The DNA fragments of the human Glo-1 promoter containing either the T allele or the G allele at rs1781735 were cloned into the pGL3 vector, and the promoter activities were determined in SH-SY5Y cells. The Glo-1 promoter fragments encompassing nucleotides from −1,200 base pairs (bp) to +22 bp (relative to transcription start site +1) were amplified, which came from the two individual homozygotes with respect to the corresponding genotypes (GG and TT) for rs1781735. Next, we cloned the amplified fragments into the pGL3-basic plasmid vectors using the primers shown in Supplementary Table 3 and verified the constructs using bidirectional DNA sequencing. The reporter constructs containing either the G allele or the T allele were transiently cotransfected into SH-SY5Y cells along with the Renilla reporter gene (pRLTK) for internal normalization. After a 48-h incubation, the luciferase and Renilla luciferase activities were measured using the Dual-Luciferase Reporter Assay System according to the instructions from the manufacturer (Promega Corporation, Madison, WI, USA).

Gene Expression Patterns and Expression Quantitative Trait Locus (eQTL) Analysis

We replicated Glo-1 expression levels and the eQTL of the G/T polymorphism at rs1781735 using data from Genotype-Tissue Expression (GTEx; https://www.gtexportal.org) Analysis Release V8 (dbGaP: phs000424.v8.p2) in 14 brain regions (Sul et al., 2013; Human Genomics, 2015; Melé et al., 2015).

fMRI Acquisition and Data Processing

The resting-state fMRI images and high-resolution T1-weighted images were acquired using a 3.0-T GE (Discovery MR750) Signa Scanner (GE Medical Systems, WI, USA) with an 8-channel phased-array head coil at the Department of Radiology of Affiliated Hospital of Guangdong Medical University. The high-resolution T1-weighted images were obtained using three-dimensional T1-weighted magnetization-prepared rapid gradient echoes (MPRAGE) with the following parameters: time of repetition (TR) = 8.16 ms, time of echo (TE) = 3.18 ms, field of view (FOV) = 512 mm × 512 mm, flip angle = 9 degrees, slice thickness/gap = 1 mm, matrix = 256 × 256, slices = 172. The fMRI images were collected using a gradient-echo planar imaging (EPI) sequence. The imaging parameters were as follows: TR = 2,000 ms, TE = 30 ms, FOV = 230 mm × 230 mm, matrix = 64 × 64, slices = 38, slice thickness = 3.6 mm, 240 dynamics, and scan time = 8 min.

The amplitude of low-frequency fluctuation (ALFF) based on the blood oxygen level-dependent (BOLD) fMRI signals was analyzed. Data pre-processing was performed using Statistical Parametric Mapping (SPM) 12 software (http://www.fil.ion.ucl.ac.uk/spm/) based on MATLAB 2019. First, the initial 10 scan volumes were discarded to allow for steady-state magnetization. The slice timing and head motion were corrected using the remaining scan volumes. The motions exceeding 2 mm displacement and 2 degrees rotation were excluded. The standard Montreal Neurological Institute (MNI) EPI template (3 mm isotropic voxels) was used for spatial normalization. Spatial smoothing was performed using a Gaussian kernel of 6 × 6 × 6 mm full width at half-maximum. Finally, temporal handpass filtering (0.01–0.08 Hz) was conducted to remove the low-frequency drifts and physiological high-frequency noise.

The Statistical Analyses

Statistical analysis was performed using SPSS 22.0 software (SPSS, Chicago, IL, USA), and statistical graphs were generated using GraphPad Prism version 8.0.2 (GraphPad Software Inc, San Diego, CA, USA). Hardy-Weinberg equilibrium (HWE) was tested using the χ² test. The statistical analyses comparing the genotypes and allele distributions were performed using the χ² test. The generalized odds ratios (ORs) and 95% confidence interval (CI) of the alleles were calculated. The multiple comparisons were performed using false discovery rate (FDR) correction. The Linkage disequilibrium (LD) status and haplotype analysis were determined using the Haploview 4.2 program. Only those haplotypes with frequencies >3% were further analyzed. The power calculations were performed using QUANTO 1.2 software. The descriptive statistics of the clinical characteristics, mRNA expression, and enzymatic activity of Glo-1 are presented as the mean ± standard deviation (SD). The statistical analyses were performed using the independent-sample t-tests. The receiver operating characteristic (ROC) curves were constructed from the mRNA expression and enzymatic activity data in the patients with schizophrenia. The value of P < 0.05 was considered statistically significant for the above tests. Transcripts per kilobase million (TPM) of RNA-sequencing and eQTL data in the GTEx database were downloaded and extracted using R version 3.6.2. (R Core Team, Vienna, Austria).

A multitissue eQTL plot for rs1781735 was calculated in the GTEx browser (Sul et al., 2013). The normalized effect size (NES) of the eQTLs is defined as the slope of the linear regression of normalized expression data vs. the three genotype categories using single tissue eQTL analysis, representing eQTL effect size. It is computed as the effect of the alternative allele (ALT) T relative to the reference allele (REF) G in rs1781735 (chr6_38704303_G_T_b38) of Glo-1 in the reference human genome GRCh38/hg38.

The second-level analysis of fMRI data was performed using the SPM 12 software. The independent-sample t-tests were used to examine the genotype effects on ALFF with sex, age, and education as covariates. The multiple comparison corrections were performed within each analytical imaging modality using familywise error (FWE) corrected at a two-tailed threshold of P < 0.05. The results were visualized by BrainNet Viewer version 1.7 (Beijing Normal University, Beijing, China) (Xia et al., 2013).

Results

Glo-1 Confers the Risk of Schizophrenia

To explore the role of Glo-1 in the risk of schizophrenia, the biochemical effect of Glo-1 was analyzed in first-onset antipsychotic-naïve patients with schizophrenia and age- and sex-matched controls (Supplementary Table 4). Both the mRNA expression and enzymatic activity were much lower in the patients with schizophrenia than in controls (t = 4.42, P =3.98 × 10⁻⁵; t = 5.02, P = 1.40 × 10⁻⁶, respectively) (Figures 2A,B and Supplementary Table 4).

FIGURE 2

Figure 2. Discrimination of glyoxalase 1 (Glo-1) mRNA levels (A) and enzymatic activity levels (B) in peripheral blood from the first-onset antipsychotic-naïve patients with schizophrenia and healthy controls. A receiver operating characteristic (ROC) curve was constructed to evaluate the role of Glo-1 mRNA expression (C) and enzymatic activity (D) in the diagnosis of schizophrenia. Discrimination of Glo-1 mRNA levels (E) and enzymatic activity levels (F) in genotypes of rs1781735 in peripheral blood of first-onset antipsychotic-naïve patients with schizophrenia and healthy controls. AUC, the area under the ROC curve. *, P < 0.05. ***, P < 0.001. The error bar represents the SD of the mean.

The ROC curves were constructed to further evaluate the correlation of Glo-1 mRNA and enzyme activity with the diagnosis of schizophrenia. For Glo-1 mRNA, the area under the curve (AUC) was 0.747 (P = 4.75 × 10⁻⁶, 95% CI: 0.65–0.85), representing a potential diagnostic capability. The calculation of the cut-off value was based on the maximum Youden Index (sensitivity + specificity −1). The cut-off value of the mRNA expression level was 4.69 (ΔΔCT), with the sensitivity and specificity values of 0.952 and 0.604, respectively (Figure 2C). For Glo-1 enzyme activity, the AUC was 0.751 (P = 1.43 × 10⁻⁷), the 95% CI was 0.671–0.832, and the cut-off value was 3.05 (unit/ml Hb), with the sensitivity and specificity values of 0.671 and 0.797, respectively (Figure 2D). These combined results indicated that Glo-1 may confer a risk for schizophrenia.

The rs1781735 T/G Polymorphism Might Be the Genetic Source of Schizophrenia Vulnerability

To further identify the Glo-1 genetic source of risk for schizophrenia, a case-control study was performed in which 1,069 patients with schizophrenia and 1,023 healthy controls were recruited to analyze the risk effects of Glo-1 variants on schizophrenia. The demographic characteristics of the sample are shown in Supplementary Table 5. There were no significant differences in the demographic characteristics between the patients with schizophrenia and healthy controls and all selected SNPs were in HWE (P > 0.05). The significant differences were observed in SNP rs1781735 between the patients with schizophrenia and healthy controls for genotype (χ² = 8.64, P = 0.020) and allele (χ² = 6.59, P = 0.020) distributions after FDR correction (Table 1). The T allele of rs1781735 was more frequent in the patients with schizophrenia than in the healthy controls. However, no statistically significant associations of allele or genotype distributions with schizophrenia were observed for rs9470916, rs1130914, or rs4746 (all P > 0.05). The haplotype analysis was conducted for these SNPs, and the corresponding results are shown in Supplementary Table 7. The frequency of the C-T-A-G haplotype, corresponding to the rs9470916-rs1130534-rs4746-rs1781735 polymorphism, was significantly lower overall in the patients with schizophrenia than in the healthy controls (χ² = 5.13, OR = 0.86, 95% CI (0.75–0.98), P = 0.024). The LD of these SNPs is shown in Supplementary Table 6. The power calculations were performed using the MAF of each allele with an OR of 1.5 at the 0.05 level (Supplementary Table 8). The estimated statistical power was from 0.993 to 1. Our findings could be statistically strong with appropriate sample size.

TABLE 1

Table 1. Genotype and allele distributions of glyoxalase 1 (Glo-1) polymorphisms in the participants.

The rs1781735 T/G Polymorphism Might Be a Functional SNP of Glo-1 in Schizophrenia

We speculated that the rs1781735 SNP, located in the promoter region of Glo-1, may affect Glo-1 expression and ultimately contribute to the pathogenesis of schizophrenia. To evaluate the effects of rs1781735 T/G on Glo-1 transcription activity, we cloned DNA fragments of the human Glo-1 promoter, containing either the T allele or the G allele at rs1781735, into the pGL3 vector, and the promoter activities were assessed in the SH-SY5Y cells. As shown in Figure 1C, the luciferase assay indicated that the Glo-1 promoter containing the T allele of rs1781735 reduced the relative promoter activity by 11.2% compared with the G allele in the SH-SY5Y (P = 0.016) cells.

Next, we suspected that rs1781735 may affect the Glo-1 mRNA expression and enzymatic activity in the patients with schizophrenia and healthy controls. The mRNA expression and enzymatic activity data were grouped by genotypes of rs1781735. The differences in mRNA expression were found in both the patients with schizophrenic (t = 3.54, P = 0.001) and controls (t = 3.34, P = 0.003) when the GG genotype was compared with the TT genotype (Figure 2E and Supplementary Table 4). For Glo-1 enzymatic activity, the differences were observed in the patients with schizophrenia and healthy controls when the GG genotype was compared to TT genotype (t = 2.54, P = 0.015 in patients, t = 2.59, P = 0.014 in controls) and when the GG + GT genotype was compared with the TT genotype (t = 2.13, P = 0.036 in patients, t = 2.16, P = 0.034 in controls) (Figure 2F and Supplementary Table 4). The T allele in rs1781735 has a decreasing tendency of mRNA expression and enzymatic activity compared with the G allele in the peripheral blood of both first-onset antipsychotic-naïve patients with schizophrenia and healthy controls.

Collectively, these data indicate that rs1781735 is likely a functional SNP, and sequence changes at this site can influence the Glo-1 transcription activation, mRNA expression, and enzymatic activity in both the patients with schizophrenia and healthy controls.

Glo-1 Is Widely Expressed in the Brain and the rs1781735 T/G Polymorphism Could Be an eQTL for Glo-1

We further evaluated the functional relevance of the rs1781735 polymorphism in Glo-1 expression in the brain. The Glo-1 expression patterns in the brain were defined and clarified using eQTL analysis conducted in the 14 brain regions and whole blood in the post-mortem healthy participant data from the GTEx database. The Glo-1 was broadly expressed in the 14 brain regions {such as the frontal cortex [Brodmann area 9 (BA9)], cerebellar hemisphere, and anterior cingulate cortex (BA24).} and whole blood (Figure 1B). As shown in Figure 1D, rs1781735 affected the Glo-1 gene expression at the mRNA level in the 14 human brain regions and whole blood, and the TT genotype was associated with the reduced Glo-1 expression compared with the GG genotypes in the 14 human brain regions and whole blood.

The rs1781735 T/G Polymorphism Affects the Brain Function in First-Onset Antipsychotic-Naïve Patients With Schizophrenia

We hypothesized that Glo-1 might affect the brain function of schizophrenia, considering its crucial role in neurodetoxification. To examine the functional effect of Glo-1 in the brain, rs1781735 was applied as a proxy of Glo-1 expression in the neuroimaging analysis for the Glo-1 functional localization in the brain of patients with schizophrenia.

The ALFF based on BOLD fMRI signals was analyzed to evaluate the effect of Glo-1 on neuronal activation in the brain. There were no significant differences in the age, sex, years of education, symptom phenotypes evaluated by PANSS between the genotype groups (Supplementary Table 9). In first-onset antipsychotic-naïve patients with schizophrenia, TT genotype exhibited significantly decreased ALFF in the left middle frontal gyrus (BA9) (Talairach: −45, 18, 51; 236 mm³, t = 4.15, P < 0.001, FWE corrected) compared with the G carriers (combined GG and GT genotypes) (Figure 3).

FIGURE 3

Figure 3. Significant regions mapped onto the standard Ch2 Template using BrainNet Viewer. The rs1781735 genotype effect on the amplitude of low-frequency fluctuation (ALFF) map in the patients (TT genotype > G carriers). Areas colored blue gradient correspond to the t-value. ***, P < 0.001, familywise error (FWE) corrected. The error bar represents the SD of the mean.

Discussion

This study investigated the effect of Glo-1 on schizophrenia risks at multiple levels, such as genetic variants, transcription, expression, protein function, and phenotypes.

The previous studies identified downregulated Glo-1 expression in the blood of Japanese patients with schizophrenia (Arai et al., 2010), suggesting that the dysregulation of Glo-1 might be related to the molecular causes or consequences of schizophrenia. To further identify the effect of Glo-1 on the risk of schizophrenia and more accurately measure the changes in Glo-1, we attempted to reduce or eliminate the influence of other factors, such as antipsychotic drugs and disease treatment period, on the results. To maintain the homogeneity in the cohort in this study, we analyzed Glo-1 expression and enzymatic activity in the peripheral blood of antipsychotic-naïve first-episode patients with schizophrenia and a sex-, age-, and education-matched healthy control group from the same community. We demonstrated that both the Glo-1 mRNA expression and activity were significantly decreased in the blood of patients compared with the healthy controls. These results indicate credible evidence of the potential biomarker effect of Glo-1 on schizophrenia.

We demonstrated that a novel functional SNP (rs1781735) in the promoter of the Glo-1 gene is associated with the risk of schizophrenia, and a large proportion of individuals with schizophrenia are homozygous for the rs1781735 T allele. Next, we identified the T allele of rs1781735, which effectively reduced the transcription activation of Glo-1 and was associated with the reduced RNA expression and enzymatic activity in the antipsychotic-naïve first-episode patients with schizophrenia. Finally, the evidence not only from the variants in rs1718735 but also from the mRNA expression and protein function suggested that Glo-1 is involved in the pathogenesis of schizophrenia in our Chinese cohorts (the hypothesized metabolic process is illustrated in Figure 4). To our knowledge, this is the first study to report a relationship between the rs1781735 and schizophrenia in a Chinese population.

FIGURE 4

Figure 4. A schematic illustrating the hypothesized regulatory process of the rs1781735 G-T mutation in the etiology of schizophrenia. Dicarbonyl compounds, such as methylglyoxal (MG), as a toxic species could lead to neural dysfunction and symptomatology in schizophrenia. Dicarbonyl compounds can be effectively detoxified via the metabolic pathway of the glyoxalase system. MG reacts with glutathione (GSH) and is converted to S-D-lactoylglutathione (SLG) by Glo-1. SLG is catalyzed into D-lactate by glyoxalase-2 (Glo-2), and GSH is recycled. As a rate-limiting enzyme in the glyoxalase system, Glo-1 plays a critical role in the detoxification of carbonyl stress in schizophrenia.

In schizophrenia, the deficits in neural function (Shenton et al., 2001) are well-established phenotypic measures, that are related to the disease severity and outcome (Lieberman et al., 2001; Wojtalik et al., 2017). Imaging genetics is an emerging field to identify the genetic sources of brain activity and functional connectivity in vivo in a non-invasive manner (Krajcovic et al., 2019). However, there is still a gap in the understanding of the underlying molecular mechanisms by which the gene acts on the neural function phenotypes due to the absence of gene expression information in the corresponding and spontaneous brain regions. The eQTL could be an effective approach to bridge this gap by offering the genetic variations as a proxy to interpret gene expression levels in the brains (Battle et al., 2017; Huckins et al., 2019). It is necessary to combine eQTL analysis with imaging genetics research to elucidate the neural function phenotypes and the underlying molecular mechanisms.

Based on the evidence that Glo-1 is broadly expressed in the brains and that rs1718735 could be an eQTL for Glo-1, a neuroimaging analysis was further performed to investigate the relationship of Glo-1 and brain function in schizophrenia using rs1781735 as a proxy of Glo-1 expression. The rs1781735 TT genotype, which is associated with a low Glo-1 expression in brains, displayed decreased ALFF in the left middle frontal gyrus compared with the G carriers in the patients with schizophrenia. We can infer that the low expression of Glo-1 implies the decreased neural activity in the left middle frontal gyrus in patients with schizophrenia. This is in accordance with the brain expression pattern results showing that the highest Glo-1 expression in the frontal cortex (BA9) in GTEx (the brain areas involved but not completely overlapping with the neuroimaging atlas). The Gray matter reduction and dysfunction of the left middle frontal gyrus have been identified in schizophrenia (Li et al., 2018) and are associated with the symptoms of schizophrenia, such as delusion and language processing (Gao et al., 2015; Cantisani et al., 2018). Toriumi et al. reported that the combination of Glo-1 dysfunction and vitamin B6 (VB6) deficiency may cause behavioral deficits via mitochondrial dysfunction and oxidative stress in the frontal cortex in the schizophrenia mouse model (Toriumi et al., 2021). In this study, we first identified evidence that Glo-1 is associated with dysfunction in the frontal cortex in patients with schizophrenia. The combination of the evidence and our results introduce eQTL data as a bridge for the genetic imaging studies and provide a glimpse of the functional localization of Glo-1 in the brain.

The ethnic heterogeneity in the genetic cohort studies should be highlighted in this study. To date, the research on the relationship between the polymorphisms in the Glo-1 gene and schizophrenia susceptibility has been reported only in a Japanese population (Arai et al., 2010; Bangel et al., 2015). Arai found that a missense mutation polymorphism at rs4746 (Ala111Glu) and two heterozygous frameshift mutations were associated with the risk of schizophrenia, as demonstrated by DNA sequencing of 1,761 patients with schizophrenia and 1,921 control subjects. In the present study, we attempted to replicate the positive association between the rs4746 of the Glo-1 gene and schizophrenia, as reported by Arai et al. (2010), in our Han population of China, we did not find an association of rs4746 with schizophrenia. The same situation was observed for rs1781735, which has been reported to have a negative association with schizophrenia in the population of Japan. This discordance may stem from the ethnic differences that exist between the population of China and Japan. The MAF of rs4746 exhibits high population diversity: 0.354–0.475 in Europeans, 0.267 in sub-Saharan Africans, 0.239–0.395 in African Americans, and 0.033–0.125 in Asian populations (Junaid et al., 2004; Sacco et al., 2007; Wu et al., 2008). We detected 50 CC genotype carriers among the 2,092 Chinese individuals. The C allele frequency in the control group of the present study was 0.134, similar to the values in the Han Chinese population of Taiwan (the frequency for the controls was 0.124) and the HapMap database for the Han population of China, but the allele frequency was much lower in the population of Japan (0.033). For rs1781735, the frequency of the G allele in our cohort from China was 0.392, but in the cohort from Japan, it was 0.443. In addition, the diagnostic heterogeneity, cohort homogeneity, environmental exposure, sample sizes, and research strategies between the two studies may also explain these differences. Hence, the results of the genetic association studies need to be determined in other cohorts from different ethnic backgrounds.

There are possible limitations to this study. For the rs1781735 effect on mRNA expression and enzymatic activity analysis, the sample size in the rs1781735 GG sub-group is limited, especially in the control group, which may induce variation and unstable results (Supplementary Figure 4). But the result was supported by further analysis. For the mRNA expression analysis, the evidence was supported by eQTL analysis with 670 whole blood samples (Figure 2D). For the enzymatic activity, a significant result was observed in the GG+GT subgroup comparing with the TT subgroup. However, this issue still needs to be identified in a larger sample and the other cohorts. In the rs1781735 effect on the phenotypes of schizophrenia, we only perform the resting-state fMRI analysis. More comprehensive analysis of the neural characteristics, such as the brain structure, function connectivity analysis, and the relationship with the clinical characteristics, need to be investigated in further studies.

Conclusions

To our knowledge, this study is the first to investigate Glo-1 risk in schizophrenia at multiple levels, such as genetic variants, transcription, expression, protein function, and phenotypes. Significant differences in variant distribution, RNA expression, and enzymatic activity were found between the patients with schizophrenia and controls. We first reported the effect of the Glo-1 gene on the neural function of schizophrenia using a functional SNP (rs1781735) as the proxy. The most prominent Glo-1 signal intensities were observed in the left middle frontal gyrus, suggesting that Glo-1 is involved in the dysfunction of the left middle frontal gyrus in schizophrenia. Glo-1 protects the nervous system against toxic carbonyl stress, and its reduced expression may play a key role in the etiology of schizophrenia. Our investigation lays the foundation for a better understanding of the influence of Glo-1 in schizophrenia and provides a referable analysis workflow in imaging genetics research.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Affiliated Hospital of Guangdong Medical University. The patients/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

GM, YW, KL, ZBL, WC, and YL conceived and designed the experiments. JY, CL, XZ, XX, QC, and XW performed the experiments. JY, SL, XL, BH, SX, DZ, ZT, JF, DL, ZD, XY, and ZXL collected the clinical samples and analyzed the data. JY and SL analyzed the data. JY and GM wrote the original draft. All authors contributed to the article and approved the submitted version.

Funding

This study was supported by the National Natural Science Foundation of China (81670252, 81770034, and 81571157) and the Guangdong Basic and Applied Basic Foundation (2019A1515011306 and 2019A1515011424).

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/fnmol.2021.739526/full#supplementary-material

References

Allen, N. C., Bagade, S., McQueen, M. B., Ioannidis, J. P., Kavvoura, F. K., Khoury, M. J., et al. (2008). Systematic meta-analyses and field synopsis of genetic association studies in schizophrenia: the SzGene database. Nat. Genet. 40, 827–834. doi: 10.1038/ng.171

Arai, M., Yuzawa, H., Nohara, I., Ohnishi, T., Obata, N., Iwayama, Y., et al. (2010). Enhanced carbonyl stress in a subpopulation of schizophrenia. Arch. Gen. Psychiatry. 67, 589–597. doi: 10.1001/archgenpsychiatry.2010.62

Bangel, F. N., Yamada, K., Arai, M., Iwayama, Y., Balan, S., Toyota, T., et al. (2015). Genetic analysis of the glyoxalase system in schizophrenia. Prog. Neuropsychopharmacol. Biol. Psychiatry. 59, 105–110. doi: 10.1016/j.pnpbp.2015.01.014

Barua, M., Jenkins, E. C., Chen, W., Kuizon, S., Pullarkat, R. K., and Junaid, M. A. (2011). Glyoxalase I polymorphism rs2736654 causing the Ala111Glu substitution modulates enzyme activity–implications for autism. Autism. Res. 4, 262–270. doi: 10.1002/aur.197

Battle, A., Brown, C. D., Engelhardt, B. E., and Montgomery, S. B. (2017). Genetic effects on gene expression across human tissues. Nature. 550, 204–213. doi: 10.1038/nature24277

Cannon, T. D. (1996). Abnormalities of brain structure and function in schizophrenia: implications for aetiology and pathophysiology. Ann. Med. 28, 533–539. doi: 10.3109/07853899608999117

Cantisani, A., Stegmayer, K., Federspiel, A., Bohlhalter, S., Wiest, R., and Walther, S. (2018). Blood perfusion in left inferior and middle frontal gyrus predicts communication skills in schizophrenia. Psychiatry Res. Neuroimaging. 274, 7–10. doi: 10.1016/j.pscychresns.2018.02.002

Distler, M. G., Plant, L. D., Sokoloff, G., Hawk, A. J., Aneas, I., Wuenschell, G. E., et al. (2012). Glyoxalase 1 increases anxiety by reducing GABAA receptor agonist methylglyoxal. J. Clin. Invest. 122, 2306–2315. doi: 10.1172/jci61319

Du, F., Li, Y., Shen, J., Zhao, Y., Kaboli, P. J., Xiang, S., et al. (2019). Glyoxalase 1 gene improves the antistress capacity and reduces the immune inflammatory response. BMC Genet. 20, 95. doi: 10.1186/s12863-019-0795-z

Frandsen, J., Choi, S. R., and Narayanasamy, P. (2020). Neural Glyoxalase Pathway Enhancement by Morin Derivatives in an Alzheimer's Disease Model. ACS Chem. Neurosci. 11, 356–366. doi: 10.1021/acschemneuro.9b00566

Fu, J., Zuo, X., Yin, J., Luo, X., Li, Z., Lin, J., et al. (2017). Association of Polymorphisms of the Receptor for Advanced Glycation Endproducts Gene with Schizophrenia in a Han Chinese Population. Biomed. Res. Int. 2017, 6379639. doi: 10.1155/2017/6379639

Fujimoto, M., Uchida, S., Watanuki, T., Wakabayashi, Y., Otsuki, K., Matsubara, T., et al. (2008). Reduced expression of glyoxalase-1 mRNA in mood disorder patients. Neurosci. Lett. 438, 196–199. doi: 10.1016/j.neulet.2008.04.024

Gabriele, S., Lombardi, F., Sacco, R., Napolioni, V., Altieri, L., Tirindelli, M. C., et al. (2014). The GLO1 C332 (Ala111) allele confers autism vulnerability: family-based genetic association and functional correlates. J. Psychiatr. Res. 59, 108–116. doi: 10.1016/j.jpsychires.2014.07.021

Gale, C. P., Futers, T. S., and Summers, L. K. (2004a). Common polymorphisms in the glyoxalase-1 gene and their association with pro-thrombotic factors. Diab. Vasc. Dis. Res. 1, 34–39. doi: 10.3132/dvdr.2004.004

Gale, C. P., Futers, T. S., and Summers, L. K. (2004b). Common polymorphisms in the glyoxalase-1 gene and their association with pro-thrombotic factors. Diabetes Vascular Dis. Res. 1, 34–39.

Gao, B., Wang, Y., Liu, W., Chen, Z., Zhou, H., Yang, J., et al. (2015). Spontaneous activity associated with delusions of schizophrenia in the left medial superior frontal gyrus: a resting-state fMRI study. PLoS ONE 10, e0133766. doi: 10.1371/journal.pone.0133766

Hara, T., Toyoshima, M., Hisano, Y., Balan, S., Iwayama, Y., Aono, H., et al. (2021). Glyoxalase I disruption and external carbonyl stress impair mitochondrial function in human induced pluripotent stem cells and derived neurons. Transl. Psychiatry. 11, 275. doi: 10.1038/s41398-021-01392-w

He, Y., Zhou, C., Huang, M., Tang, C., Liu, X., Yue, Y., et al. (2020). Glyoxalase system: A systematic review of its biological activity, related-diseases, screening methods and small molecule regulators. Biomed. Pharmacother. 131, 110663. doi: 10.1016/j.biopha.2020.110663

Hovatta, I., Tennant, R. S., Helton, R., Marr, R. A., Singer, O., Redwine, J. M., et al. (2005). Glyoxalase 1 and glutathione reductase 1 regulate anxiety in mice. Nature. 438, 662–666. doi: 10.1038/nature04250

Huckins, L. M., Dobbyn, A., Ruderfer, D. M., Hoffman, G., Wang, W., Pardinas, A. F., et al. (2019). Gene expression imputation across multiple brain regions provides insights into schizophrenia risk. Nat. Genet. 51, 659–674. doi: 10.1038/s41588-019-0364-4

Human Genomics (2015). The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. Science. 348, 648–660. doi: 10.1126/science.1262110

International Schizophrenia, C., Purcell, S. M., Wray, N. R., Stone, J. L., Visscher, P. M., O'Donovan, M. C., et al. (2009). Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature. 460, 748–752. doi: 10.1038/nature08185

Itokawa, M., Miyashita, M., Arai, M., Dan, T., Takahashi, K., Tokunaga, T., et al. (2018). Pyridoxamine: A novel treatment for schizophrenia with enhanced carbonyl stress. Psychiat. Clini. Neurosci. 72, 35–44. doi: 10.1111/pcn.12613

Junaid, M. A., Kowal, D., Barua, M., Pullarkat, P. S., Sklower Brooks, S., and Pullarkat, R. K. (2004). Proteomic studies identified a single nucleotide polymorphism in glyoxalase I as autism susceptibility factor. Am. J. Med. Genet. A. 131, 11–17. doi: 10.1002/ajmg.a.30349

Katsuta, N., Ohnuma, T., Maeshima, H., Takebayashi, Y., Higa, M., Takeda, M., et al. (2014). Significance of measurements of peripheral carbonyl stress markers in a cross-sectional and longitudinal study in patients with acute-stage schizophrenia. Schizophr. Bull. 40, 1366–1373. doi: 10.1093/schbul/sbt234

Klingler, E., Francis, F., Jabaudon, D., and Cappello, S. (2021). Mapping the molecular and cellular complexity of cortical malformations. Science. 371:eaba4517. doi: 10.1126/science.aba4517

Kold-Christensen, R., and Johannsen, M. (2020). Methylglyoxal metabolism and aging-related disease: moving from correlation toward causation. Trends. Endocrinol. Metab. 31, 81–92. doi: 10.1016/j.tem.2019.10.003

Kouidrat, Y., Amad, A., Desailloud, R., Diouf, M., Fertout, E., Scoury, D., et al. (2013). Increased advanced glycation end-products (AGEs) assessed by skin autofluorescence in schizophrenia. J. Psychiatr. Res. 47, 1044–1048. doi: 10.1016/j.jpsychires.2013.03.016

Kovac, J., Podkrajsek, K. T., Luksic, M. M., and Battelino, T. (2015). Weak association of glyoxalase 1 (GLO1) variants with autism spectrum disorder. Eur. Child Adolesc. Psychiatry. 24, 75–82. doi: 10.1007/s00787-014-0537-8

Krajcovic, B., Fajnerova, I., Horacek, J., Kelemen, E., Kubik, S., Svoboda, J., et al. (2019). Neural and neuronal discoordination in schizophrenia: From ensembles through networks to symptoms. Acta. Physiol. (Oxf). 226, e13282. doi: 10.1111/apha.13282

Kuhla, B., Boeck, K., Schmidt, A., Ogunlade, V., Arendt, T., Munch, G., et al. (2007). Age- and stage-dependent glyoxalase I expression and its activity in normal and Alzheimer's disease brains. Neurobiol. Aging. 28, 29–41. doi: 10.1016/j.neurobiolaging.2005.11.007

Kurz, A., Rabbani, N., Walter, M., Bonin, M., Thornalley, P., Auburger, G., et al. (2011). Alpha-synuclein deficiency leads to increased glyoxalase I expression and glycation stress. Cell. Mol. Life. Sci. 68, 721–733. doi: 10.1007/s00018-010-0483-7

Li, Y., Li, W.-,x., Xie, D.-,j., Wang, Y., Cheung, E. F. C., and Chan, R. C. K. (2018). Grey matter reduction in the caudate nucleus in patients with persistent negative symptoms: An ALE meta-analysis. Schizophrenia Res. 192, 9–15. doi: 10.1016/j.schres.2017.04.005

Lieberman, J., Chakos, M., Wu, H., Alvir, J., Hoffman, E., Robinson, D., et al. (2001). Longitudinal study of brain morphology in first episode schizophrenia. Biol. Psychiat. 49, 487–499.

Lv, H., Wei, G. Y., Guo, C. S., Deng, Y. F., Jiang, Y. M., Gao, C., et al. (2020). 20S proteasome and glyoxalase 1 activities decrease in erythrocytes derived from Alzheimer's disease patients. Neural. Regen. Res. 15, 178–183. doi: 10.4103/1673-5374.264473

Lv, N., Hao, S., Luo, C., Abukiwan, A., Hao, Y., Gai, F., et al. (2018). miR-137 inhibits melanoma cell proliferation through downregulation of GLO1. Sci. China Life Sci. 61, 541–549. doi: 10.1007/s11427-017-9138-9

Ma, G., Yin, J., Fu, J., Luo, X., Zhou, H., Tao, H., et al. (2014). Association of a miRNA-137 polymorphism with schizophrenia in a Southern Chinese Han population. Biomed. Res. Int. 2014, 751267. doi: 10.1155/2014/751267

Melé, M., Ferreira, P. G., Reverter, F., DeLuca, D. S., Monlong, J., Sammeth, M., et al. (2015). Human Genomics. The human transcriptome across tissues and individuals. Science. 348, 660–665. doi: 10.1126/science.aaa0355

Miyashita, M., Arai, M., Kobori, A., Ichikawa, T., Toriumi, K., Niizato, K., et al. (2014a). Clinical features of schizophrenia with enhanced carbonyl stress. Schizophr. Bull. 40, 1040–1046. doi: 10.1093/schbul/sbt129

Miyashita, M., Arai, M., Yuzawa, H., Niizato, K., Oshima, K., Kushima, I., et al. (2014b). Replication of enhanced carbonyl stress in a subpopulation of schizophrenia. Psychiatry Clin. Neurosci. 68, 83–84. doi: 10.1111/pcn.12081

Miyata, T., van Ypersele de Strihou, C., Kurokawa, K., and Baynes, J. W. (1999). Alterations in nonenzymatic biochemistry in uremia: origin and significance of “carbonyl stress” in long-term uremic complications. Kidney Int. 55, 389–399. doi: 10.1046/j.1523-1755.1999.00302.x

Mizutani, R., Saiga, R., Takeuchi, A., Uesugi, K., Terada, Y., Suzuki, Y., et al. (2019). Three-dimensional alteration of neurites in schizophrenia. Transl. Psychiatry. 9, 85. doi: 10.1038/s41398-019-0427-4

Modai, S., and Shomron, N. (2016). Molecular Risk Factors for Schizophrenia. Trends. Mol. Med. 22, 242–253. doi: 10.1016/j.molmed.2016.01.006

Ohnishi, T., Balan, S., Toyoshima, M., Maekawa, M., Ohba, H., Watanabe, A., et al. (2019). Investigation of betaine as a novel psychotherapeutic for schizophrenia. EBioMedicine. 45, 432–446. doi: 10.1016/j.ebiom.2019.05.062

Owen, M. J., Sawa, A., and Mortensen, P. B. (2016). Schizophrenia. Lancet. 388, 86–97. doi: 10.1016/S0140-6736(15)01121-6

Peculis, R., Konrade, I., Skapare, E., Fridmanis, D., Nikitina-Zake, L., Lejnieks, A., et al. (2013). Identification of glyoxalase 1 polymorphisms associated with enzyme activity. Gene. 515, 140–143. doi: 10.1016/j.gene.2012.11.009

Politi, P., Minoretti, P., Falcone, C., Martinelli, V., and Emanuele, E. (2006). Association analysis of the functional Ala111Glu polymorphism of the glyoxalase I gene in panic disorder. Neurosci. Lett. 396, 163–166. doi: 10.1016/j.neulet.2005.11.028

Sacco, R., Papaleo, V., Hager, J., Rousseau, F., Moessner, R., Militerni, R., et al. (2007). Case-control and family-based association studies of candidate genes in autistic disorder and its endophenotypes: TPH2 and GLO1. BMC Med. Genet. 8, 11. doi: 10.1186/1471-2350-8-11

Sarker, M. K., Lee, J. H., Lee, D. H., Chun, K. H., and Jun, H. S. (2020). Attenuation of diabetic kidney injury in DPP4-deficient rats; role of GLP-1 on the suppression of AGE formation by inducing glyoxalase 1. Aging (Albany NY). 12, 593–610. doi: 10.18632/aging.102643

Schizophrenia Psychiatric Genome-Wide Association Study (GWAS) Consortium (2011). Genome-wide association study identifies five new schizophrenia loci. Nature Genetics. 43, 969–976. doi: 10.1038/ng.940

Schizophrenia Working Group of the Psychiatric Genomics C. (2014). Biological insights from 108 schizophrenia-associated genetic loci. Nature. 511, 421–427. doi: 10.1038/nature13595

Shenton, M. E., Dickey, C. C., Frumin, M., and McCarley, R. W. (2001). A review of MRI findings in schizophrenia. Schizophr. Res. 49, 1–52.

Shi, J., Levinson, D. F., Duan, J., Sanders, A. R., Zheng, Y., Pe'er, I., et al. (2009). Common variants on chromosome 6p22.1 are associated with schizophrenia. Nature 460, 753–757. doi: 10.1038/nature08192

Stefansson, H., Ophoff, R. A., Steinberg, S., Andreassen, O. A., Cichon, S., Rujescu, D., et al. (2009). Common variants conferring risk of schizophrenia. Nature. 460, 744–747. doi: 10.1038/nature08186

Sul, J. H., Han, B., Ye, C., Choi, T., and Eskin, E. (2013). Effectively identifying eQTLs from multiple tissues by combining mixed model and meta-analytic approaches. PLoS Genet. 9, e1003491. doi: 10.1371/journal.pgen.1003491

Tao, H., Si, L., Zhou, X., Liu, Z., Ma, Z., Zhou, H., et al. (2016). Role of glyoxalase I gene polymorphisms in late-onset epilepsy and drug-resistant epilepsy. J. Neurol. Sci. 363, 200–206. doi: 10.1016/j.jns.2016.01.052

Torii, Y., Iritani, S., Sekiguchi, H., Habuchi, C., Fujishiro, H., Kushima, I., et al. (2020). The accumulation of advanced glycation end-products in a schizophrenic patient with a glyoxalase 1 frameshift mutation: An autopsy study. Schizophr. Res. 223, 356–358. doi: 10.1016/j.schres.2020.09.012

Toriumi, K., Berto, S., Koike, S., Usui, N., Dan, T., Suzuki, K., et al. (2021). Combined glyoxalase 1 dysfunction and vitamin B6 deficiency in a schizophrenia model system causes mitochondrial dysfunction in the prefrontal cortex. Redox. Biol. 45, 102057. doi: 10.1016/j.redox.2021.102057

Wang, Y., Chen, R., Zhou, X., Guo, R., Yin, J., Li, Y., et al. (2020). miR-137: A Novel Therapeutic Target for Human Glioma. Mol. Ther. Nucleic. Acids. 21, 614–622. doi: 10.1016/j.omtn.2020.06.028

Wang, Z., Zhang, J., Chen, L., Li, J., Zhang, H., and Guo, X. (2019). Glycine Suppresses AGE/RAGE Signaling Pathway and Subsequent Oxidative Stress by Restoring Glo1 Function in the Aorta of Diabetic Rats and in HUVECs. Oxid. Med. Cell. Longev. 2019, 4628962. doi: 10.1155/2019/4628962

Wojtalik, J. A., Smith, M. J., Keshavan, M. S., and Eack, S. M. (2017). A Systematic and Meta-analytic Review of Neural Correlates of Functional Outcome in Schizophrenia. Schizophr. Bull. 43, 1329–1347. doi: 10.1093/schbul/sbx008

Wu, Y. Y., Chien, W. H., Huang, Y. S., Gau, S. S., and Chen, C. H. (2008). Lack of evidence to support the glyoxalase 1 gene (GLO1) as a risk gene of autism in Han Chinese patients from Taiwan. Prog. Neuropsychopharmacol. Biol. Psychiatry. 32, 1740–1744. doi: 10.1016/j.pnpbp.2008.07.019

Xia, M., Wang, J., and He, Y. (2013). BrainNet Viewer: a network visualization tool for human brain connectomics. PLoS ONE. 8, e68910. doi: 10.1371/journal.pone.0068910

Yin, J., Lin, J., Luo, X., Chen, Y., Li, Z., Ma, G., et al. (2014). miR-137: a new player in schizophrenia. Int. J. Mol. Sci. 15, 3262–3271. doi: 10.3390/ijms15023262

Yoshioka, S., Odani, H., Furuhashi, T., Tanaka, T., and Ogawa, T. (2021). Increase in the peripheral blood methylglyoxal levels in 10% of hospitalized chronic schizophrenia patients. Nagoya. J. Med. Sci. 83, 51–61. doi: 10.18999/nagjms.83.1.51