Genome wide association study identifies four loci for early onset schizophrenia

Whole-exome sequencing in a Chinese sample provides preliminary evidence for the link between rare/low-frequency immune-related variants and early-onset schizophrenia

Rare and low-frequency variants contribute to schizophrenia (SCZ), and may influence its age-at-onset (AAO). We examined the association of rare or low-frequency deleterious coding variants in Chinese patients with SCZ. We collected DNA samples in 197 patients with SCZ spectrum disorder and 82 healthy controls (HC), and performed exome sequencing. The AAO variable was ascertained in the majority of SCZ participants for identify the early-onset (EOS, AAO<=18) and adult-onset (AOS, AAO>18) subgroups. We examined the overall association of rare/low-frequency, damaging variants in SCZ versus HC, EOS versus HC, and AOS versus HC at the gene and gene-set levels using Sequence Kernel Association Test. The quantitative rare-variant association test of AAO was conducted. Resampling was used to obtain empirical p-values and to control for family-wise error rate (FWER). In binary-trait association tests, we identified 5 potential candidate risk genes and 10 gene ontology biological processes (GOBP) terms, among which PADI2 reached FWER-adjusted significance. In quantitative rare-variant association tests, we found marginally significant correlations of AAO with alterations in 4 candidate risk genes, and 5 GOBP pathways. Together, the biological and functional profiles of these genes and gene sets supported the involvement of perturbations of neural systems in SCZ, and altered immune functions in EOS.

Identifying frequency-dependent imaging genetic associations via hypergraph-structured multi-task sparse canonical correlation analysis

Identifying complex associations between genetic variations and imaging phenotypes is a challenging task in the research of brain imaging genetics. The previous study has proved that neuronal oscillations within distinct frequency bands are derived from frequency-dependent genetic modulation. Thus it is meaningful to explore frequency-dependent imaging genetic associations, which may give important insights into the pathogenesis of brain disorders. In this work, the hypergraph-structured multi-task sparse canonical correlation analysis (HS-MTSCCA) was developed to explore the associations between multi-frequency imaging phenotypes and single-nucleotide polymorphisms (SNPs). Specifically, we first created a hypergraph for the imaging phenotypes of each frequency and the SNPs, respectively. Then, a new hypergraph-structured constraint was proposed to learn high-order relationships among features in each hypergraph, which can introduce biologically meaningful information into the model. The frequency-shared and frequency-specific imaging phenotypes and SNPs could be identified using the multi-task learning framework. We also proposed a useful strategy to tackle this algorithm and then demonstrated its convergence. The proposed method was evaluated on four simulation datasets and a real schizophrenia dataset. The experimental results on synthetic data showed that HS-MTSCCA outperforms the other competing methods according to canonical correlation coefficients, canonical weights, and cosine similarity. And the results on real data showed that HS-MTSCCA could obtain superior canonical coefficients and canonical weights. Furthermore, the identified frequency-shared and frequency-specific biomarkers could provide more interesting and meaningful information, demonstrating that HS-MTSCCA is a powerful method for brain imaging genetics.

miRNA cargo in circulating vesicles from neurons is altered in individuals with schizophrenia and associated with severe disease

While RNA expression appears to be altered in several brain disorders, the constraints of postmortem analysis make it impractical for well-powered population studies and biomarker development. Given that the unique molecular composition of neurons are reflected in their extracellular vesicles (EVs), we hypothesized that the fractionation of neuron derived EVs provides an opportunity to specifically profile their encapsulated contents noninvasively from blood. To investigate this hypothesis, we determined miRNA expression in microtubule associated protein 1B (MAP1B)-enriched serum EVs derived from neurons from a large cohort of individuals with schizophrenia and nonpsychiatric comparison participants. We observed dysregulation of miRNA in schizophrenia subjects, in particular those with treatment-resistance and severe cognitive deficits. These data support the hypothesis that schizophrenia is associated with alterations in posttranscriptional regulation of synaptic gene expression and provides an example of the potential utility of tissue-specific EV analysis in brain disorders.

Early-onset schizophrenia is associated with immune-related rare variants in a Chinese sample

Background

Rare variants are likely to contribute to schizophrenia (SCZ), given the large discrepancy between the heritability estimated from twin and GWAS studies. Furthermore, the nature of the rare-variant contribution to SCZ may vary with the "age-at-onset" (AAO), since early-onset has been suggested as being indicative of neurodevelopment deviance.

Objective

To examine the association of rare deleterious coding variants in early- and adult-onset SCZ in a Chinese sample.

Method

Exome sequencing was performed on DNA from 197 patients with SCZ spectrum disorder and 82 healthy controls (HC) of Chinese ancestry recruited in Hong Kong. We also gathered AAO information in the majority of SCZ samples. Patients were classified into early-onset (EOS, AAO<18) and adult-onset (AOS, AAO>18). We collapsed the rare variants to improve statistical power and examined the overall association of rare variants in SCZ versus HC, EOS versus HC, and AOS versus HC at the gene and gene-set levels by Sequence Kernel Association Test. The quantitative rare-variant association test of AAO was also conducted. We focused on variants which were predicted to have a medium or high impact on the protein-encoding process as defined by Ensembl. We applied a 100000-time permutation test to obtain empirical p-values, with significance threshold set at p < 1e -3 to control family-wise error rates. Moreover, we compared the burden of targeted rare variants in significant risk genes and gene sets in cases and controls.

Results

Based on several binary-trait association tests (i.e., SCZ vs HC, EOS vs HC and AOS vs HC), we identified 7 candidate risk genes and 20 gene ontology biological processes (GOBP) terms, which exhibited higher burdens in SCZ than in controls. Based on quantitative rare-variant association tests, we found that alterations in 5 candidate risk genes and 7 GOBP pathways were significantly correlated with AAO. Based on biological and functional profiles of the candidate risk genes and gene sets, our findings suggested that, in addition to the involvement of perturbations in neural systems in SCZ in general, altered immune responses may be specifically implicated in EOS.

Conclusion

Disrupted immune responses may exacerbate abnormal perturbations during neurodevelopment and trigger the early onset of SCZ. We provided evidence of rare variants increasing SCZ risk in the Chinese population.

Common genetic variants contribute to heritability of age at onset of schizophrenia

Schizophrenia (SCZ) is a complex disorder that typically arises in late adolescence or early adulthood. Age at onset (AAO) of SCZ is associated with long-term outcomes of the disease. We explored the genetic architecture of AAO with a genome-wide association study (GWAS), heritability, polygenic risk score (PRS), and copy number variant (CNV) analyses in 4 740 subjects of European ancestry. Although no genome-wide significant locus was identified, SNP-based heritability of AAO was estimated to be between 17 and 21%, indicating a moderate contribution of common variants. We also performed cross-trait PRS analyses with a set of mental disorders and identified a negative association between AAO and common variants for SCZ, childhood maltreatment and attention-deficit/hyperactivity disorder. We also investigated the role of copy number variants (CNVs) in AAO and found an association with the length and number of deletions (P-value = 0.03), whereas the presence of CNVs previously reported in SCZ was not associated with earlier onset. To our knowledge, this is the largest GWAS of AAO of SCZ to date in individuals from European ancestry, and the first study to determine the involvement of common variants in the heritability of AAO. Finally, we evidenced the role played by higher SCZ load in determining AAO but discarded the role of pathogenic CNVs. Altogether, these results shed light on the genetic architecture of AAO, which needs to be confirmed with larger studies.

The anatomical networks based on probabilistic structurally connectivity in bipolar disorder across mania, depression, and euthymic states

BACKGROUNDS

There have pieces of evidence of the distinct aberrant functional network topology profile in bipolar disorder (BD) across mania, depression, and euthymic episodes. However, the underlying anatomical network topology pattern in BD across different episodes is unclear.

METHODS

We calculated the whole-brain probabilistic structurally connectivity across 143 subjects (72 with BD [34 depression; 13 mania; 25 euthymic] and 53 healthy controls), and used graph theory to examine the trait- and state-related topology alterations of the structural connectome in BD. The correlation analysis was further conducted to explore the relationship between detected network measures and clinical symptoms.

RESULTS

There no omnibus alteration of any global network metrics were observed across all diagnostic groups. In the regional network metrics level, bipolar depression showed increased clustering coefficient in the right lingual gyrus compared with all other groups, and the increased clustering coefficient in the right lingual gyrus positively correlated with depression, anxiety, and illness burden symptoms but negatively correlated with mania symptoms; manic and euthymic patients showed decreased clustering coefficient in the left inferior occipital gyrus compared with HCs.

LIMITATIONS

The moderate sample size of all patient groups (especially for subjects with mania) might have contributed to the negative findings of the trait feature in this study.

CONCLUSIONS

We demonstrated the altered regional connectivity pattern in the occipital lobe of the bipolar depression and mania episode, especially the lingual gyrus. The association of the clustering coefficient in the lingual gyrus with clinical symptoms helps monitor the state of BD.

Emerging epigenetic dynamics in gut-microglia brain axis: experimental and clinical implications for accelerated brain aging in schizophrenia

Brain aging, which involves a progressive loss of neuronal functions, has been reported to be premature in probands affected by schizophrenia (SCZ). Evidence shows that SCZ and accelerated aging are linked to changes in epigenetic clocks. Recent cross-sectional magnetic resonance imaging analyses have uncovered reduced brain reserves and connectivity in patients with SCZ compared to typically aging individuals. These data may indicate early abnormalities of neuronal function following cyto-architectural alterations in SCZ. The current mechanistic knowledge on brain aging, epigenetic changes, and their neuropsychiatric disease association remains incomplete. With this review, we explore and summarize evidence that the dynamics of gut-resident bacteria can modulate molecular brain function and contribute to age-related neurodegenerative disorders. It is known that environmental factors such as mode of birth, dietary habits, stress, pollution, and infections can modulate the microbiota system to regulate intrinsic neuronal activity and brain reserves through the vagus nerve and enteric nervous system. Microbiota-derived molecules can trigger continuous activation of the microglial sensome, groups of receptors and proteins that permit microglia to remodel the brain neurochemistry based on complex environmental activities. This remodeling causes aberrant brain plasticity as early as fetal developmental stages, and after the onset of first-episode psychosis. In the central nervous system, microglia, the resident immune surveillance cells, are involved in neurogenesis, phagocytosis of synapses and neurological dysfunction. Here, we review recent emerging experimental and clinical evidence regarding the gut-brain microglia axis involvement in SCZ pathology and etiology, the hypothesis of brain reserve and accelerated aging induced by dietary habits, stress, pollution, infections, and other factors. We also include in our review the possibilities and consequences of gut dysbiosis activities on microglial function and dysfunction, together with the effects of antipsychotics on the gut microbiome: therapeutic and adverse effects, role of fecal microbiota transplant and psychobiotics on microglial sensomes, brain reserves and SCZ-derived accelerated aging. We end the review with suggestions that may be applicable to the clinical setting. For example, we propose that psychobiotics might contribute to antipsychotic-induced therapeutic benefits or adverse effects, as well as reduce the aging process through the gut-brain microglia axis. Overall, we hope that this review will help increase the understanding of SCZ pathogenesis as related to chronobiology and the gut microbiome, as well as reveal new concepts that will serve as novel treatment targets for SCZ.

Consolidation of metabolomic, proteomic, and GWAS data in connective model of schizophrenia

Despite of multiple systematic studies of schizophrenia based on proteomics, metabolomics, and genome-wide significant loci, reconstruction of underlying mechanism is still a challenging task. Combination of the advanced data for quantitative proteomics, metabolomics, and genome-wide association study (GWAS) can enhance the current fundamental knowledge about molecular pathogenesis of schizophrenia. In this study, we utilized quantitative proteomic and metabolomic assay, and high throughput genotyping for the GWAS study. We identified 20 differently expressed proteins that were validated on an independent cohort of patients with schizophrenia, including ALS, A1AG1, PEDF, VTDB, CERU, APOB, APOH, FASN, GPX3, etc. and almost half of them are new for schizophrenia. The metabolomic survey revealed 18 group-specific compounds, most of which were the part of transformation of tyrosine and steroids with the prevalence to androgens (androsterone sulfate, thyroliberin, thyroxine, dihydrotestosterone, androstenedione, cholesterol sulfate, metanephrine, dopaquinone, etc.). The GWAS assay mostly failed to reveal significantly associated loci therefore 52 loci with the smoothened p < 10^-5 were fractionally integrated into proteome-metabolome data. We integrated three omics layers and powered them by the quantitative analysis to propose a map of molecular events associated with schizophrenia psychopathology. The resulting interplay between different molecular layers emphasizes a strict implication of lipids transport, oxidative stress, imbalance in steroidogenesis and associated impartments of thyroid hormones as key interconnected nodes essential for understanding of how the regulation of distinct metabolic axis is achieved and what happens in the conditioned proteome and metabolome to produce a schizophrenia-specific pattern.

Genetic insights into childhood-onset schizophrenia: The yield of clinical exome sequencing

Childhood-onset schizophrenia (COS) is a rare form of schizophrenia with an onset prior to 13 years of age. Although genetic factors play a role in COS etiology, only a few causal variants have been reported to date. This study presents a diagnostic exome sequencing (ES) in 37 Israeli Jewish families with a proband diagnosed with COS. By implementing a trio/duo ES approach and applying a well-established diagnostic pipeline, we detected clinically significant variants in 7 probands (19 %). These single nucleotide variants and indels were mostly inherited. The implicated genes were ANKRD11, GRIA2, CHD2, CLCN3, CLTC, IGF1R and MICU1. In a secondary analysis that compared COS patients to 4721 healthy controls, we observed that patients had a significant enrichment of rare loss of function (LoF) variants in LoF intolerant genes associated with developmental diseases. Taken together, ES could be considered as a valuable tool in the genetic workup for COS patients.

Biological hypotheses, risk factors, and biomarkers of schizophrenia

Both the discovery of biomarkers of schizophrenia and the verification of biological hypotheses of schizophrenia are an essential part of the process of understanding the etiology of this mental disorder. Schizophrenia has long been considered a neurodevelopmental disease whose symptoms are caused by impaired synaptic signal transduction and brain neuroplasticity. Both the onset and chronic course of schizophrenia are associated with risk factors-induced disruption of brain function and the establishment of a new homeostatic setpoint characterized by biomarkers. Different risk factors and biomarkers can converge to the same symptoms of schizophrenia, suggesting that the primary cause of the disease can be highly individual. Schizophrenia-related biomarkers include measurable biochemical changes induced by stress (elevated allostatic load), mitochondrial dysfunction, neuroinflammation, oxidative and nitrosative stress, and circadian rhythm disturbances. Here is a summary of selected valid biological hypotheses of schizophrenia formulated based on risk factors and biomarkers, neurodevelopment, neuroplasticity, brain chemistry, and antipsychotic medication. The integrative neurodevelopmental-vulnerability-neurochemical model is based on current knowledge of the neurobiology of the onset and progression of the disease and the effects of antipsychotics and psychotomimetics and reflects the complex and multifactorial nature of schizophrenia.

NRN1 Gene as a Potential Marker of Early-Onset Schizophrenia: Evidence from Genetic and Neuroimaging Approaches

Included in the neurotrophins family, the Neuritin 1 gene (NRN1) has emerged as an attractive candidate gene for schizophrenia (SZ) since it has been associated with the risk for the disorder and general cognitive performance. In this work, we aimed to further investigate the association of NRN1 with SZ by exploring its role on age at onset and its brain activity correlates. First, we developed two genetic association analyses using a family-based sample (80 early-onset (EO) trios (offspring onset ≤ 18 years) and 71 adult-onset (AO) trios) and an independent case–control sample (120 healthy subjects (HS), 87 EO and 138 AO patients). Second, we explored the effect of NRN1 on brain activity during a working memory task (N-back task; 39 HS, 39 EO and 39 AO; matched by age, sex and estimated IQ). Different haplotypes encompassing the same three Single Nucleotide Polymorphisms(SNPs, rs3763180–rs10484320–rs4960155) were associated with EO in the two samples (GCT, TCC and GTT). Besides, the GTT haplotype was associated with worse N-back task performance in EO and was linked to an inefficient dorsolateral prefrontal cortex activity in subjects with EO compared to HS. Our results show convergent evidence on the NRN1 association with EO both from genetic and neuroimaging approaches, highlighting the role of neurotrophins in the pathophysiology of SZ.

Association of NKAPL rs1635 With Cognitive Function in Early-Onset Schizophrenia

BACKGROUND: Schizophrenia is a severe mental disorder with high heritability, and cognitive dysfunction is one of the core features. Growing evidence suggests the genetic risk of schizophrenia may contribute to cognitive impairments. The variant rs1635 (nucleotide sequence: c.455C>A; amino acid sequence: T152N) located on the (NFKB activating protein like) NKAPL gene confers risk for schizophrenia and might play a role in the neurodevelopmental process, which is particularly relevant to cognitive function. However, the relationship between rs1635 and cognitive function remains unclear.

METHODS: A total of 130 patients with early-onset schizophrenia (EOS) and 300 patients with adult-onset schizophrenia (AOS) of Han Chinese were recruited and underwent neurocognitive tests by using the MATRICS Consensus Cognitive Battery (MCCB). The NKAPL rs1635 was genotyped by using DNA sequencing. The peripheral blood NKAPL mRNA expression level was examined in 152T or 152N carriers (n = 20) in EOS patients, by using the qRT-PCR. The phosphorylation level of NAKPL T152N polymorphism was detected by cell experiments. In utero electroporation of mouse embryos was examined to explore the effect of Nkapl on neuronal migration.

RESULTS: Compared with rs1635 AA and AC carriers, CC (the CC genotype encodes the protein NKAPL-152T) carriers of EOS patients performed better in cognitive domain of speed of processing (t = 2.644, p = 0.009), trail making test (t = 2.221, p = 0.028) and category fluency (t = 2.578, p = 0.011). However, patients with AOS exhibited no significant differences in seven domains among the three genotype groups. There were no significant differences in cognitive performance between EOS and AOS. In EOS patients, NKAPL mRNA level in NKAPL-152N carriers is significantly lower than that of NKAPL-152T carriers. The phosphorylation level of NKAPL-152N is significantly decreased compared to NKAPL-152T. In utero electroporation showed that Nkapl deletion impairs the embryonic radial migration process.

CONCLUSION: The present study found that NKAPL rs1635 was associated with cognitive impairments and peripheral blood mRNA expression level in EOS patients. The NKAPL full-length protein is required for embryonic cortical neuronal migration. The phosphorylation level of NKAPL-152N is significantly decreased. The NKAPL T152N may affect the NAKPL mRNA expression level and embryonic cortical neuronal migration by regulating the NAKPL protein phosphorylation. These data suggest that NKAPL rs1635 affects cognitive function by regulating early brain development in early-onset schizophrenia.

Association between MAP3K4 gene polymorphisms and the risk of schizophrenia susceptibility in a Northeast Chinese Han population

Schizophrenia stands out as one of the most devastating psychiatric disorders. Previous findings have shown that schizophrenia is a polygenic genetic disorder. Thus, abnormal neurodevelopment and neurogenesis may be associated with the etiology of schizophrenia, so genes which affect these processes may be potential candidate genes of schizophrenia. Mitogen-activated protein kinase kinase kinase 4 (MAP3K4) gene is a member of the mitogen-activated protein kinase family. Taking into account previous findings, MAP3K4 plays a crucial role in the fundamental pathology of various nervous system diseases. In the present study, we aim to explore the association of MAP3K4 and schizophrenia in an independent case-control sample including 627 schizophrenic patients and 1175 healthy controls from a Northeast Chinese Han population. Both the allelic and genotypic association analyses showed that 6 SNPs in MAP3K4 were significantly associated with schizophrenia (rs590988, rs625977, rs9295134, rs12110787, rs1001808 and rs9355870). After rigorous Bonferroni correction, 4 SNPs (rs9295134, rs12110787, rs1001808 and rs9355870) were still significantly associated with the disease. The haplotype composed of these four SNPs also showed significantly global and individual association with schizophrenia. These results suggest that MAP3K4 is a susceptibility gene for schizophrenia in the Northeast Chinese Han population.

Identification of a Risk Locus at 7p22.3 for Schizophrenia and Bipolar Disorder in East Asian Populations

Background: Shared psychopathological features and mechanisms have been observed between schizophrenia (SZ) and bipolar disorder (BD), but their common risk genes and full genetic architectures remain to be fully characterized. The genome-wide association study (GWAS) datasets offer the opportunity to explore this scientific question using combined genetic data from enormous samples, ultimately allowing a better understanding of the onset and development of these illnesses.

Methods: We have herein performed a genome-wide meta-analysis in two GWAS datasets of SZ and BD respectively (24,600 cases and 40,012 controls in total, discovery sample), followed by replication analyses in an independent sample of 4,918 SZ cases and 5,506 controls of Han Chinese origin (replication sample). The risk SNPs were then explored for their correlations with mRNA expression of nearby genes in multiple expression quantitative trait loci (eQTL) datasets.

Results: The single nucleotide polymorphisms (SNPs) rs1637749 and rs3800908 at 7p22.3 region were significant in both discovery and replication samples, and exhibited genome-wide significant associations when combining all East Asian SZ and BD samples (29,518 cases and 45,518 controls). The risk SNPs were also significant in GWAS of SZ and BD among Europeans. Both risk SNPs significantly predicted lower expression of MRM2 in the whole blood and brain samples in multiple datasets, which was consistent with its reduced mRNA level in the brains of SZ patients compared with normal controls. The risk SNPs were also associated with MAD1L1 expression in the whole blood sample.

Discussion: We have identified a novel genome-wide risk locus associated with SZ and BD in East Asians, adding further support for the putative common genetic risk of the two illnesses. Our study also highlights the necessity and importance of mining public datasets to explore risk genes for complex psychiatric diseases.

Weak Association Between the Glutamate Decarboxylase 1 Gene (GAD1) and Schizophrenia in Han Chinese Population

Objectives

Schizophrenia (SZ) is a complex psychiatric disorder with high heritability, and genetic components are thought to be pivotal risk factors for this illness. The glutamate decarboxylase 1 gene (GAD1) was hypothesized to be a candidate risk locus for SZ given its crucial role in the GABAergic neurotransmission system, and previous studies have examined the associations of single nucleotide polymorphisms (SNPs) spanning the GAD1 gene with SZ. However, inconsistent results were obtained. We hence examined the associations between GAD1 SNPs and SZ in two independent case-control samples of Han Chinese ancestry.

Materials and Methods

Two Han Chinese SZ case-control samples, referred as the discovery sample and the replication sample, respectively, were recruited for the current study. The discovery sample comprised of 528 paranoid SZ cases (with age of first onset ≥ 18) and 528 healthy controls; the independent replication sample contained 1,256 early onset SZ cases (with age of first onset < 18) and 2,661 healthy controls. Logistic regression analysis was performed to examine the associations between GAD1 SNPs and SZ.

Results

Ten SNPs covering GAD1 gene were analyzed in the discovery sample, and two SNPs showed nominal associations with SZ (rs2241165, P = 0.0181, OR = 1.261; rs2241164, P = 0.0225, OR = 1.219). SNP rs2241164 was also nominally significant in the independent replication sample (P = 0.0462, OR = 1.110), and the significance became stronger in a subsequent meta-analysis combining both discovery and replication samples (P = 0.00398, OR = 1.138). Nevertheless, such association could not survive multiple corrections, although the effect size of rs2241164 was comparable with other SZ risk loci identified in genome-wide association studies (GWAS) in Han Chinese population. We also examined the associations between GAD1 SNPs and SZ in published datasets of SZ GWAS in East Asians and Europeans, and no significant associations were observed.

Conclusion

We observed weak associations between GAD1 SNPs and risk of SZ in Han Chinese populations. Further analyses in larger Han Chinese samples with more detailed phenotyping are necessary to elucidate the genetic correlation between GAD1 SNPs and SZ.