Gene expression imputation across multiple brain regions provides insights into schizophrenia risk

Aberrant patterns of spontaneous brain activity in schizophrenia: A resting-state fMRI study and classification analysis

BACKGROUND

Schizophrenia is a prevalent mental disorder, leading to severe disability. Currently, the absence of objective biomarkers hinders effective diagnosis. This study was conducted to explore the aberrant spontaneous brain activity and investigate the potential of abnormal brain indices as diagnostic biomarkers employing machine learning methods.

METHODS

A total of sixty-one schizophrenia patients and seventy demographically matched healthy controls were enrolled in this study. The static indices of resting-state functional magnetic resonance imaging (rs-fMRI) including amplitude of low frequency fluctuations (ALFF), fractional ALFF (fALFF), regional homogeneity (ReHo), and degree centrality (DC) were calculated to evaluate spontaneous brain activity. Subsequently, a sliding-window method was then used to conduct temporal dynamic analysis. The comparison of static and dynamic rs-fMRI indices between the patient and control groups was conducted using a two-sample t-test. Finally, the machine learning analysis was applied to estimate the diagnostic value of abnormal indices of brain activity.

RESULTS

Schizophrenia patients exhibited a significant increase ALFF value in inferior frontal gyrus, alongside significant decreases in fALFF values observed in left postcentral gyrus and right cerebellum posterior lobe. Pervasive aberrations in ReHo indices were observed among schizophrenia patients, particularly in frontal lobe and cerebellum. A noteworthy reduction in voxel-wise concordance of dynamic indices was observed across gray matter regions encompassing the bilateral frontal, parietal, occipital, temporal, and insular cortices. The classification analysis achieved the highest values for area under curve at 0.87 and accuracy at 81.28% when applying linear support vector machine and leveraging a combination of abnormal static and dynamic indices in the specified brain regions as features.

CONCLUSIONS

The static and dynamic indices of brain activity exhibited as potential neuroimaging biomarkers for the diagnosis of schizophrenia.

Antidepressant Exposure and DNA Methylation: Insights from a Methylome-Wide Association Study

Importance

Understanding antidepressant mechanisms could help design more effective and tolerated treatments.

Objective

Identify DNA methylation (DNAm) changes associated with antidepressant exposure.

Design

Case-control methylome-wide association studies (MWAS) of antidepressant exposure were performed from blood samples collected between 2006-2011 in Generation Scotland (GS). The summary statistics were tested for enrichment in specific tissues, gene ontologies and an independent MWAS in the Netherlands Study of Depression and Anxiety (NESDA). A methylation profile score (MPS) was derived and tested for its association with antidepressant exposure in eight independent cohorts, alongside prospective data from GS.

Setting

Cohorts; GS, NESDA, FTC, SHIP-Trend, FOR2107, LBC1936, MARS-UniDep, ALSPAC, E-Risk, and NTR.

Participants

Participants with DNAm data and self-report/prescription derived antidepressant exposure.

Main Outcomes and Measures

Whole-blood DNAm levels were assayed by the EPIC/450K Illumina array (9 studies, N exposed = 661, N unexposed = 9,575) alongside MBD-Seq in NESDA (N exposed = 398, N unexposed = 414). Antidepressant exposure was measured by self- report and/or antidepressant prescriptions.

Results

The self-report MWAS (N = 16,536, N exposed = 1,508, mean age = 48, 59% female) and the prescription-derived MWAS (N = 7,951, N exposed = 861, mean age = 47, 59% female), found hypermethylation at seven and four DNAm sites (p < 9.42x10 -8 ), respectively. The top locus was cg26277237 ( KANK1, p self-report = 9.3x10 -13 , p prescription = 6.1x10 -3 ). The self-report MWAS found a differentially methylated region, mapping to DGUOK-AS1 ( p adj = 5.0x10 -3 ) alongside significant enrichment for genes expressed in the amygdala, the "synaptic vesicle membrane" gene ontology and the top 1% of CpGs from the NESDA MWAS (OR = 1.39, p < 0.042). The MPS was associated with antidepressant exposure in meta-analysed data from external cohorts (N studies = 9, N = 10,236, N exposed = 661, f3 = 0.196, p < 1x10 -4 ).

Conclusions and Relevance

Antidepressant exposure is associated with changes in DNAm across different cohorts. Further investigation into these changes could inform on new targets for antidepressant treatments.

3 Key Points

Question: Is antidepressant exposure associated with differential whole blood DNA methylation?Findings: In this methylome-wide association study of 16,536 adults across Scotland, antidepressant exposure was significantly associated with hypermethylation at CpGs mapping to KANK1 and DGUOK-AS1. A methylation profile score trained on this sample was significantly associated with antidepressant exposure (pooled f3 [95%CI]=0.196 [0.105, 0.288], p < 1x10 -4 ) in a meta-analysis of external datasets. Meaning: Antidepressant exposure is associated with hypermethylation at KANK1 and DGUOK-AS1 , which have roles in mitochondrial metabolism and neurite outgrowth. If replicated in future studies, targeting these genes could inform the design of more effective and better tolerated treatments for depression.

Prioritizing susceptibility genes for the prognosis of male-pattern baldness with transcriptome-wide association study

BACKGROUND

Male-pattern baldness (MPB) is the most common cause of hair loss in men. It can be categorized into three types: type 2 (T2), type 3 (T3), and type 4 (T4), with type 1 (T1) being considered normal. Although various MPB-associated genetic variants have been suggested, a comprehensive study for linking these variants to gene expression regulation has not been performed to the best of our knowledge.

RESULTS

In this study, we prioritized MPB-related tissue panels using tissue-specific enrichment analysis and utilized single-tissue panels from genotype-tissue expression version 8, as well as cross-tissue panels from context-specific genetics. Through a transcriptome-wide association study and colocalization analysis, we identified 52, 75, and 144 MPB associations for T2, T3, and T4, respectively. To assess the causality of MPB genes, we performed a conditional and joint analysis, which revealed 10, 11, and 54 putative causality genes for T2, T3, and T4, respectively. Finally, we conducted drug repositioning and identified potential drug candidates that are connected to MPB-associated genes.

CONCLUSIONS

Overall, through an integrative analysis of gene expression and genotype data, we have identified robust MPB susceptibility genes that may help uncover the underlying molecular mechanisms and the novel drug candidates that may alleviate MPB.

Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.

A miR-137-Related Biological Pathway of Risk for Schizophrenia Is Associated With Human Brain Emotion Processing

BACKGROUND

miR-137 is a microRNA involved in brain development, regulating neurogenesis and neuronal maturation. Genome-wide association studies have implicated miR-137 in schizophrenia risk but do not explain its involvement in brain function and underlying biology. Polygenic risk for schizophrenia mediated by miR-137 targets is associated with working memory, although other evidence points to emotion processing. We characterized the functional brain correlates of miR-137 target genes associated with schizophrenia while disentangling previously reported associations of miR-137 targets with working memory and emotion processing.

METHODS

Using RNA sequencing data from postmortem prefrontal cortex (N = 522), we identified a coexpression gene set enriched for miR-137 targets and schizophrenia risk genes. We validated the relationship of this set to miR-137 in vitro by manipulating miR-137 expression in neuroblastoma cells. We translated this gene set into polygenic scores of coexpression prediction and associated them with functional magnetic resonance imaging activation in healthy volunteers (n1 = 214; n2 = 136; n3 = 2075; n4 = 1800) and with short-term treatment response in patients with schizophrenia (N = 427).

RESULTS

In 4652 human participants, we found that 1) schizophrenia risk genes were coexpressed in a biologically validated set enriched for miR-137 targets; 2) increased expression of miR-137 target risk genes was mediated by low prefrontal miR-137 expression; 3) alleles that predict greater gene set coexpression were associated with greater prefrontal activation during emotion processing in 3 independent healthy cohorts (n1, n2, n3) in interaction with age (n4); and 4) these alleles predicted less improvement in negative symptoms following antipsychotic treatment in patients with schizophrenia.

CONCLUSIONS

The functional translation of miR-137 target gene expression linked with schizophrenia involves the neural substrates of emotion processing.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia

Multimorbidity represents an increasingly important public health challenge with far-reaching implications for health management and policy. Mental health and metabolic diseases have a well-established epidemiological association. In this study, we investigate the genetic intersection between type 2 diabetes and schizophrenia. We use Mendelian randomization to examine potential causal relationships between the two conditions and related endophenotypes. We report no compelling evidence that type 2 diabetes genetic liability potentially causally influences schizophrenia risk and vice versa. Our findings show that increased body mass index (BMI) has a protective effect against schizophrenia, in contrast to the well-known risk-increasing effect of BMI on type 2 diabetes risk. We identify evidence of colocalization of association signals for these two conditions at 11 genomic loci, six of which have opposing directions of effect for type 2 diabetes and schizophrenia. To elucidate these colocalizing signals, we integrate multi-omics data from bulk and single-cell gene expression studies, along with functional information. We identify putative effector genes and find that they are enriched for homeostasis and lipid-related pathways. We also highlight drug repurposing opportunities including N-methyl-D-aspartate (NMDA) receptor antagonists. Our findings provide insights into shared biological mechanisms for type 2 diabetes and schizophrenia, highlighting common factors that influence the risk of the two conditions in opposite directions and shedding light on the complex nature of this comorbidity.

Gene expression imputation provides clinical and biological insights into treatment-resistant schizophrenia polygenic risk

Genome-wide association studies (GWAS) have revealed the polygenic nature of treatment-resistant schizophrenia TRS. Gene expression imputation allowed the translation of GWAS results into regulatory mechanisms and the construction of gene expression (GReX) risk scores (GReX-RS).  In the present study we computed GReX-RS from the largest GWAS of TRS to assess its association with clinical features. We perform transcriptome imputation in the largest GWAS of TRS to find GReX associated with TRS using brain tissues. Then, for each tissue, we constructed a GReX-RS of the identified genes in a sample of 254 genotyped first episode of psychosis (FEP) patients to test its association with clinical phenotypes, including clinical symptomatology, global functioning and cognitive performance. Our analysis provides evidence that the polygenic basis of TRS includes genetic variants that modulate the expression of certain genes in certain brain areas (substantia nigra, hippocampus, amygdala and frontal cortex), which at the same time are related to clinical features in FEP patients, mainly persistence of negative symptoms and cognitive alterations in sustained attention, which have also been suggested as clinical predictors of TRS. Our results provide a clinical explanation of the polygenic architecture of TRS and give more insight into the biological mechanisms underlying TRS.

The schizophrenia syndrome, circa 2024: What we know and how that informs its nature

With new data about different aspects of schizophrenia being continually generated, it becomes necessary to periodically revisit exactly what we know. Along with a need to review what we currently know about schizophrenia, there is an equal imperative to evaluate the construct itself. With these objectives, we undertook an iterative, multi-phase process involving fifty international experts in the field, with each step building on learnings from the prior one. This review assembles currently established findings about schizophrenia (construct, etiology, pathophysiology, clinical expression, treatment) and posits what they reveal about its nature. Schizophrenia is a heritable, complex, multi-dimensional syndrome with varying degrees of psychotic, negative, cognitive, mood, and motor manifestations. The illness exhibits a remitting and relapsing course, with varying degrees of recovery among affected individuals with most experiencing significant social and functional impairment. Genetic risk factors likely include thousands of common genetic variants that each have a small impact on an individual's risk and a plethora of rare gene variants that have a larger individual impact on risk. Their biological effects are concentrated in the brain and many of the same variants also increase the risk of other psychiatric disorders such as bipolar disorder, autism, and other neurodevelopmental conditions. Environmental risk factors include but are not limited to urban residence in childhood, migration, older paternal age at birth, cannabis use, childhood trauma, antenatal maternal infection, and perinatal hypoxia. Structural, functional, and neurochemical brain alterations implicate multiple regions and functional circuits. Dopamine D-2 receptor antagonists and partial agonists improve psychotic symptoms and reduce risk of relapse. Certain psychological and psychosocial interventions are beneficial. Early intervention can reduce treatment delay and improve outcomes. Schizophrenia is increasingly considered to be a heterogeneous syndrome and not a singular disease entity. There is no necessary or sufficient etiology, pathology, set of clinical features, or treatment that fully circumscribes this syndrome. A single, common pathophysiological pathway appears unlikely. The boundaries of schizophrenia remain fuzzy, suggesting the absence of a categorical fit and need to reconceptualize it as a broader, multi-dimensional and/or spectrum construct.

Decreased CNNM2 expression in prefrontal cortex affects sensorimotor gating function, cognition, dendritic spine morphogenesis and risk of schizophrenia

Genome-wide association studies (GWASs) have reported multiple single nucleotide polymorphisms (SNPs) associated with schizophrenia, yet the underlying molecular mechanisms are largely unknown. In this study, we aimed to identify schizophrenia relevant genes showing alterations in mRNA and protein expression associated with risk SNPs at the 10q24.32-33 GWAS locus. We carried out the quantitative trait loci (QTL) and summary data-based Mendelian randomization (SMR) analyses, using the PsychENCODE dorsolateral prefrontal cortex (DLPFC) expression QTL (eQTL) database, as well as the ROSMAP and Banner DLPFC protein QTL (pQTL) datasets. The gene CNNM2 (encoding a magnesium transporter) at 10q24.32-33 was identified to be a robust schizophrenia risk gene, and was highly expressed in human neurons according to single cell RNA-seq (scRNA-seq) data. We further revealed that reduced Cnnm2 in the mPFC of mice led to impaired cognition and compromised sensorimotor gating function, and decreased Cnnm2 in primary cortical neurons altered dendritic spine morphogenesis, confirming the link between CNNM2 and endophenotypes of schizophrenia. Proteomics analyses showed that reduced Cnnm2 level changed expression of proteins associated with neuronal structure and function. Together, these results identify a robust gene in the pathogenesis of schizophrenia.

Transcriptional risk scores in Alzheimer's disease: From pathology to cognition

INTRODUCTION

Our previously developed blood-based transcriptional risk scores (TRS) showed associations with diagnosis and neuroimaging biomarkers for Alzheimer's disease (AD). Here, we developed brain-based TRS.

METHODS

We integrated AD genome-wide association study summary and expression quantitative trait locus data to prioritize target genes using Mendelian randomization. We calculated TRS using brain transcriptome data of two independent cohorts (N = 878) and performed association analysis of TRS with diagnosis, amyloidopathy, tauopathy, and cognition. We compared AD classification performance of TRS with polygenic risk scores (PRS).

RESULTS

Higher TRS values were significantly associated with AD, amyloidopathy, tauopathy, worse cognition, and faster cognitive decline, which were replicated in an independent cohort. The AD classification performance of PRS was increased with the inclusion of TRS up to 16% with the area under the curve value of 0.850.

DISCUSSION

Our results suggest brain-based TRS improves the AD classification of PRS and may be a potential AD biomarker.

HIGHLIGHTS

Transcriptional risk score (TRS) is developed using brain RNA-Seq data. Higher TRS values are shown in Alzheimer's disease (AD). TRS improves the AD classification power of PRS up to 16%. TRS is associated with AD pathology presence. TRS is associated with worse cognitive performance and faster cognitive decline.

Rediscovering tandem repeat variation in schizophrenia: challenges and opportunities

Tandem repeats (TRs) are prevalent throughout the genome, constituting at least 3% of the genome, and often highly polymorphic. The high mutation rate of TRs, which can be orders of magnitude higher than single-nucleotide polymorphisms and indels, indicates that they are likely to make significant contributions to phenotypic variation, yet their contribution to schizophrenia has been largely ignored by recent genome-wide association studies (GWAS). Tandem repeat expansions are already known causative factors for over 50 disorders, while common tandem repeat variation is increasingly being identified as significantly associated with complex disease and gene regulation. The current review summarizes key background concepts of tandem repeat variation as pertains to disease risk, elucidating their potential for schizophrenia association. An overview of next-generation sequencing-based methods that may be applied for TR genome-wide identification is provided, and some key methodological challenges in TR analyses are delineated.

Mapping the landscape of lineage-specific dynamic regulation of gene expression using single-cell transcriptomics and application to genetics of complex disease

Single-cell transcriptome data can provide insights into how genetic variation influences biological processes involved in human biology and disease. However, the identification of gene-level associations in distinct cell types faces several challenges, including the limited reference resource from population scale studies, data sparsity in single-cell RNA sequencing, and the complex cell-state pattern of expression within individual cell types. Here we develop genetic models of cell type specific and cell state adjusted gene expression in mid-brain neurons in the process of specializing from induced pluripotent stem cells. The resulting framework quantifies the dynamics of the genetic regulation of gene expression and estimates its cell type specificity. As an application, we show that the approach detects known and new genes associated with schizophrenia and enables insights into context-dependent disease mechanisms. We provide a genomic resource from a phenome-wide application of our models to more than 1500 phenotypes from the UK Biobank. Using longitudinal genetically determined expression, we implement a predictive causality framework, evaluating the prediction of future values of a target gene expression using prior values of a putative regulatory gene. Collectively, this work demonstrates the insights that can be gained into the molecular underpinnings of diseases by quantifying the genetic control of gene expression at single-cell resolution.

Spatiotemporal expression patterns of anxiety disorder-associated genes

Anxiety disorders (ADs) are the most common form of mental disorder that affects millions of individuals worldwide. Although physiological studies have revealed the neural circuits related to AD symptoms, how AD-associated genes are spatiotemporally expressed in the human brain still remains unclear. In this study, we integrated genome-wide association studies of four human AD subtypes-generalized anxiety disorder, social anxiety disorder, panic disorder, and obsessive-compulsive disorder-with spatial gene expression patterns. Our investigation uncovered a novel division among AD-associated genes, marked by significant and distinct expression enrichments in the cerebral nuclei, limbic, and midbrain regions. Each gene cluster was associated with specific anxiety-related behaviors, signaling pathways, region-specific gene networks, and cell types. Notably, we observed a significant negative correlation in the temporal expression patterns of these gene clusters during various developmental stages. Moreover, the specific brain regions enriched in each gene group aligned with neural circuits previously associated with negative decision-making and anxious temperament. These results suggest that the two distinct gene clusters may underlie separate neural systems involved in anxiety. As a result, our findings bridge the gap between genes and neural circuitry, shedding light on the mechanisms underlying AD-associated behaviors.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia

Multimorbidity represents an increasingly important public health challenge with far-reaching implications for health management and policy. Mental health and metabolic diseases have a well-established epidemiological association. In this study, we investigate the genetic intersection between type 2 diabetes and schizophrenia. We use Mendelian randomization to examine potential causal relationships between the two conditions and related endophenotypes. We report no compelling evidence that type 2 diabetes genetic liability potentially causally influences schizophrenia risk and vice versa. Our findings show that increased body mass index (BMI) has a protective effect against schizophrenia, in contrast to the well-known risk-increasing effect of BMI on type 2 diabetes risk. We identify evidence of colocalization of association signals for these two conditions at 11 genomic loci, six of which have opposing directions of effect for type 2 diabetes and schizophrenia. To elucidate these colocalizing signals, we integrate multi-omics data from bulk and single-cell gene expression studies, along with functional information. We identify high-confidence effector genes and find that they are enriched for homeostasis and lipid-related pathways. We also highlight drug repurposing opportunities including N-methyl-D-aspartate (NMDA) receptor antagonists. Our findings provide insights into shared biological mechanisms for type 2 diabetes and schizophrenia, highlighting common factors that influence the risk of the two conditions in opposite directions and shedding light on the complex nature of this comorbidity.

Human forebrain organoid-based multi-omics analyses of PCCB as a schizophrenia associated gene linked to GABAergic pathways

Identifying genes whose expression is associated with schizophrenia (SCZ) risk by transcriptome-wide association studies (TWAS) facilitates downstream experimental studies. Here, we integrated multiple published datasets of TWAS, gene coexpression, and differential gene expression analysis to prioritize SCZ candidate genes for functional study. Convergent evidence prioritized Propionyl-CoA Carboxylase Subunit Beta (PCCB), a nuclear-encoded mitochondrial gene, as an SCZ risk gene. However, the PCCB's contribution to SCZ risk has not been investigated before. Using dual luciferase reporter assay, we identified that SCZ-associated SNPs rs6791142 and rs35874192, two eQTL SNPs for PCCB, showed differential allelic effects on transcriptional activities. PCCB knockdown in human forebrain organoids (hFOs) followed by RNA sequencing analysis revealed dysregulation of genes enriched with multiple neuronal functions including gamma-aminobutyric acid (GABA)-ergic synapse. The metabolomic and mitochondrial function analyses confirmed the decreased GABA levels resulted from inhibited tricarboxylic acid cycle in PCCB knockdown hFOs. Multielectrode array recording analysis showed that PCCB knockdown in hFOs resulted into SCZ-related phenotypes including hyper-neuroactivities and decreased synchronization of neural network. In summary, this study utilized hFOs-based multi-omics analyses and revealed that PCCB downregulation may contribute to SCZ risk through regulating GABAergic pathways, highlighting the mitochondrial function in SCZ.

Dysfunction of cAMP-Protein Kinase A-Calcium Signaling Axis in Striatal Medium Spiny Neurons: A Role in Schizophrenia and Huntington's Disease Neuropathology

Background

Striatal medium spiny neurons (MSNs) are preferentially lost in Huntington's disease. Genomic studies also implicate a direct role for MSNs in schizophrenia, a psychiatric disorder known to involve cortical neuron dysfunction. It remains unknown whether the two diseases share similar MSN pathogenesis or if neuronal deficits can be attributed to cell type-dependent biological pathways. Transcription factor BCL11B, which is expressed by all MSNs and deep layer cortical neurons, was recently proposed to drive selective neurodegeneration in Huntington's disease and identified as a candidate risk gene in schizophrenia.

Methods

Using human stem cell-derived neurons lacking BCL11B as a model, we investigated cellular pathology in MSNs and cortical neurons in the context of these disorders. Integrative analyses between differentially expressed transcripts and published genome-wide association study datasets identified cell type-specific disease-related phenotypes.

Results

We uncover a role for BCL11B in calcium homeostasis in both neuronal types, while deficits in mitochondrial function and PKA (protein kinase A)-dependent calcium transients are detected only in MSNs. Moreover, BCL11B-deficient MSNs display abnormal responses to glutamate and fail to integrate dopaminergic and glutamatergic stimulation, a key feature of striatal neurons in vivo. Gene enrichment analysis reveals overrepresentation of disorder risk genes among BCL11B-regulated pathways, primarily relating to cAMP-PKA-calcium signaling axis and synaptic signaling.

Conclusions

Our study indicates that Huntington's disease and schizophrenia are likely to share neuronal pathophysiology where dysregulation of intracellular calcium homeostasis is found in both striatal and cortical neurons. In contrast, reduction in PKA signaling and abnormal dopamine/glutamate receptor signaling is largely specific to MSNs.

Prioritization of potential causative genes for schizophrenia in placenta

Our earlier work has shown that genomic risk for schizophrenia converges with early life complications in affecting risk for the disorder and sex-biased neurodevelopmental trajectories. Here, we identify specific genes and potential mechanisms that, in placenta, may mediate such outcomes. We performed TWAS in healthy term placentae (N = 147) to derive candidate placental causal genes that we confirmed with SMR; to search for placenta and schizophrenia-specific associations, we performed an analogous analysis in fetal brain (N = 166) and additional placenta TWAS for other disorders/traits. The analyses in the whole sample and stratifying by sex ultimately highlight 139 placenta and schizophrenia-specific risk genes, many being sex-biased; the candidate molecular mechanisms converge on the nutrient-sensing capabilities of placenta and trophoblast invasiveness. These genes also implicate the Coronavirus-pathogenesis pathway and showed increased expression in placentae from a small sample of SARS-CoV-2-positive pregnancies. Investigating placental risk genes for schizophrenia and candidate mechanisms may lead to opportunities for prevention that would not be suggested by study of the brain alone.

Induction of dopaminergic neurons for neuronal subtype-specific modeling of psychiatric disease risk

Dopaminergic neurons are critical to movement, mood, addiction, and stress. Current techniques for generating dopaminergic neurons from human induced pluripotent stem cells (hiPSCs) yield heterogenous cell populations with variable purity and inconsistent reproducibility between donors, hiPSC clones, and experiments. Here, we report the rapid (5 weeks) and efficient (~90%) induction of induced dopaminergic neurons (iDANs) through transient overexpression of lineage-promoting transcription factors combined with stringent selection across five donors. We observe maturation-dependent increase in dopamine synthesis and electrophysiological properties consistent with midbrain dopaminergic neuron identity, such as slow-rising after- hyperpolarization potentials, an action potential duration of ~3 ms, tonic sub-threshold oscillatory activity, and spontaneous burst firing at a frequency of ~1.0-1.75 Hz. Transcriptome analysis reveals robust expression of genes involved in fetal midbrain dopaminergic neuron identity. Specifically expressed genes in iDANs, as well as those from isogenic induced GABAergic and glutamatergic neurons, were enriched in loci conferring heritability for cannabis use disorder, schizophrenia, and bipolar disorder; however, each neuronal subtype demonstrated subtype-specific heritability enrichments in biologically relevant pathways, and iDANs alone were uniquely enriched in autism spectrum disorder risk loci. Therefore, iDANs provide a critical tool for modeling midbrain dopaminergic neuron development and dysfunction in psychiatric disease.

Integrating genetics and transcriptomics to study major depressive disorder: a conceptual framework, bioinformatic approaches, and recent findings

Major depressive disorder (MDD) is a complex and heterogeneous psychiatric syndrome with genetic and environmental influences. In addition to neuroanatomical and circuit-level disturbances, dysregulation of the brain transcriptome is a key phenotypic signature of MDD. Postmortem brain gene expression data are uniquely valuable resources for identifying this signature and key genomic drivers in human depression; however, the scarcity of brain tissue limits our capacity to observe the dynamic transcriptional landscape of MDD. It is therefore crucial to explore and integrate depression and stress transcriptomic data from numerous, complementary perspectives to construct a richer understanding of the pathophysiology of depression. In this review, we discuss multiple approaches for exploring the brain transcriptome reflecting dynamic stages of MDD: predisposition, onset, and illness. We next highlight bioinformatic approaches for hypothesis-free, genome-wide analyses of genomic and transcriptomic data and their integration. Last, we summarize the findings of recent genetic and transcriptomic studies within this conceptual framework.

Consensus molecular environment of schizophrenia risk genes in coexpression networks shifting across age and brain regions

Schizophrenia is a neurodevelopmental brain disorder whose genetic risk is associated with shifting clinical phenomena across the life span. We investigated the convergence of putative schizophrenia risk genes in brain coexpression networks in postmortem human prefrontal cortex (DLPFC), hippocampus, caudate nucleus, and dentate gyrus granule cells, parsed by specific age periods (total N = 833). The results support an early prefrontal involvement in the biology underlying schizophrenia and reveal a dynamic interplay of regions in which age parsing explains more variance in schizophrenia risk compared to lumping all age periods together. Across multiple data sources and publications, we identify 28 genes that are the most consistently found partners in modules enriched for schizophrenia risk genes in DLPFC; twenty-three are previously unidentified associations with schizophrenia. In iPSC-derived neurons, the relationship of these genes with schizophrenia risk genes is maintained. The genetic architecture of schizophrenia is embedded in shifting coexpression patterns across brain regions and time, potentially underwriting its shifting clinical presentation.

Mapping anorexia nervosa genes to clinical phenotypes

BACKGROUND

Anorexia nervosa (AN) is a psychiatric disorder with complex etiology, with a significant portion of disease risk imparted by genetics. Traditional genome-wide association studies (GWAS) produce principal evidence for the association of genetic variants with disease. Transcriptomic imputation (TI) allows for the translation of those variants into regulatory mechanisms, which can then be used to assess the functional outcome of genetically regulated gene expression (GReX) in a broader setting through the use of phenome-wide association studies (pheWASs) in large and diverse clinical biobank populations with electronic health record phenotypes.

METHODS

Here, we applied TI using S-PrediXcan to translate the most recent PGC-ED AN GWAS findings into AN-GReX. For significant genes, we imputed AN-GReX in the Mount Sinai BioMe™ Biobank and performed pheWASs on over 2000 outcomes to test the clinical consequences of aberrant expression of these genes. We performed a secondary analysis to assess the impact of body mass index (BMI) and sex on AN-GReX clinical associations.

RESULTS

Our S-PrediXcan analysis identified 53 genes associated with AN, including what is, to our knowledge, the first-genetic association of AN with the major histocompatibility complex. AN-GReX was associated with autoimmune, metabolic, and gastrointestinal diagnoses in our biobank cohort, as well as measures of cholesterol, medications, substance use, and pain. Additionally, our analyses showed moderation of AN-GReX associations with measures of cholesterol and substance use by BMI, and moderation of AN-GReX associations with celiac disease by sex.

CONCLUSIONS

Our BMI-stratified results provide potential avenues of functional mechanism for AN-genes to investigate further.

Stem Cell Models for Context-Specific Modeling in Psychiatric Disorders

Genome-wide association studies reveal the complex polygenic architecture underlying psychiatric disorder risk, but there is an unmet need to validate causal variants, resolve their target genes(s), and explore their functional impacts on disorder-related mechanisms. Disorder-associated loci regulate transcription of target genes in a cell type- and context-specific manner, which can be measured through expression quantitative trait loci. In this review, we discuss methods and insights from context-specific modeling of genetically and environmentally regulated expression. Human induced pluripotent stem cell-derived cell type and organoid models have uncovered context-specific psychiatric disorder associations by investigating tissue-, cell type-, sex-, age-, and stressor-specific genetic regulation of expression. Techniques such as massively parallel reporter assays and pooled CRISPR (clustered regularly interspaced short palindromic repeats) screens make it possible to functionally fine-map genome-wide association study loci and validate their target genes at scale. Integration of disorder-associated contexts with these patient-specific human induced pluripotent stem cell models makes it possible to uncover gene by environment interactions that mediate disorder risk, which will ultimately improve our ability to diagnose and treat psychiatric disorders.

Human forebrain organoids-based multi-omics analyses reveal PCCB's regulation on GABAergic system contributing to schizophrenia

Identifying genes whose expression is associated with schizophrenia (SCZ) risk by transcriptome-wide association studies (TWAS) facilitates downstream experimental studies. Here, we integrated multiple published datasets of TWAS (including FUSION, PrediXcan, summary-data-based Mendelian randomization (SMR), joint-tissue imputation approach with Mendelian randomization (MR-JTI)), gene coexpression, and differential gene expression analysis to prioritize SCZ candidate genes for functional study. Convergent evidence prioritized Propionyl-CoA Carboxylase Subunit Beta ( PCCB ), a nuclear-encoded mitochondrial gene, as an SCZ risk gene. However, the PCCB ’s contribution to SCZ risk has not been investigated before. Using dual luciferase reporter assay, we identified that SCZ-associated SNP rs35874192, an eQTL SNP for PCCB , showed differential allelic effects on transcriptional activities. PCCB knockdown in human forebrain organoids (hFOs) followed by RNA-seq revealed dysregulation of genes enriched with multiple neuronal functions including gamma-aminobutyric acid (GABA)-ergic synapse, as well as genes dysregulated in postmortem brains of SCZ patients or in cerebral organoids derived from SCZ patients. The metabolomic and mitochondrial function analyses confirmed the deceased GABA levels resulted from reduced tricarboxylic acid cycle in PCCB knockdown hFOs. Multielectrode array recording analysis showed that PCCB knockdown in hFOs resulted into SCZ-related phenotypes including hyper-neuroactivities and decreased synchronization of neural network. In summary, this study utilized hFOs-based multi-omics data and revealed that PCCB downregulation may contribute to SCZ risk through regulating GABAergic system, highlighting the mitochondrial function in SCZ.

Non-additive effects of schizophrenia risk genes reflect convergent downstream function

Genetic studies of schizophrenia (SCZ) reveal a complex polygenic risk architecture comprised of hundreds of risk variants, the majority of which are common in the population at-large and confer only modest increases in disorder risk. Precisely how genetic variants with individually small predicted effects on gene expression combine to yield substantial clinical impacts in aggregate is unclear. Towards this, we previously reported that the combinatorial perturbation of four SCZ risk genes ("eGenes", whose expression is regulated by common variants) resulted in gene expression changes that were not predicted by individual perturbations, being most non-additive among genes associated with synaptic function and SCZ risk. Now, across fifteen SCZ eGenes, we demonstrate that non-additive effects are greatest within groups of functionally similar eGenes. Individual eGene perturbations reveal common downstream transcriptomic effects ("convergence"), while combinatorial eGene perturbations result in changes that are smaller than predicted by summing individual eGene effects ("sub-additive effects"). Unexpectedly, these convergent and sub-additive downstream transcriptomic effects overlap and constitute a large proportion of the genome-wide polygenic risk score, suggesting that functional redundancy of eGenes may be a major mechanism underlying non-additivity. Single eGene perturbations likewise fail to predict the magnitude or directionality of cellular phenotypes resulting from combinatorial perturbations. Overall, our results indicate that polygenic risk cannot be extrapolated from experiments testing one risk gene at a time and must instead be empirically measured. By unravelling the interactions between complex risk variants, it may be possible to improve the clinical utility of polygenic risk scores through more powerful prediction of symptom onset, clinical trajectory, and treatment response, or to identify novel targets for therapeutic intervention.

Cerebellar Functional Dysconnectivity in Drug-Naïve Patients With First-Episode Schizophrenia

BACKGROUND

Cerebellar functional dysconnectivity has long been implicated in schizophrenia. However, the detailed dysconnectivity pattern and its underlying biological mechanisms have not been well-charted. This study aimed to conduct an in-depth characterization of cerebellar dysconnectivity maps in early schizophrenia.

STUDY DESIGN

Resting-state fMRI data were processed from 196 drug-naïve patients with first-episode schizophrenia and 167 demographically matched healthy controls. The cerebellum was parcellated into nine functional systems based on a state-of-the-art atlas, and seed-based connectivity for each cerebellar system was examined. The observed connectivity alterations were further associated with schizophrenia risk gene expressions using data from the Allen Human Brain Atlas.

STUDY RESULTS

Overall, we observed significantly increased cerebellar connectivity with the sensorimotor cortex, default-mode regions, ventral part of visual cortex, insula, and striatum. In contrast, decreased connectivity was shown chiefly within the cerebellum, and between the cerebellum and the lateral prefrontal cortex, temporal lobe, and dorsal visual areas. Such dysconnectivity pattern was statistically similar across seeds, with no significant group by seed interactions identified. Moreover, connectivity strengths of hypoconnected but not hyperconnected regions were significantly correlated with schizophrenia risk gene expressions, suggesting potential genetic underpinnings for the observed hypoconnectivity.

CONCLUSIONS

These findings suggest a common bidirectional dysconnectivity pattern across different cerebellar subsystems, and imply that such bidirectional alterations may relate to different biological mechanisms.

Mendelian Randomization Study Using Dopaminergic Neuron-Specific eQTL Identifies Novel Risk Genes for Schizophrenia

Multiple integrative studies have been performed to identify the potential target genes of the non-coding schizophrenia (SCZ) risk variants. However, all the integrative studies used expression quantitative trait loci (eQTL) data from bulk tissues. Considering the cell type-specific regulatory effect of many genetic variants, it is important to conduct integrative studies using cell type-specific eQTL data. Here, we conduct a Mendelian randomization (MR) study by integrating genome-wide associations of SCZ (74,776 cases and 101,023 controls) and eQTL data (N = 215) from dopaminergic neurons, which were differentiated from human-induced pluripotent stem cell (iPSC) lines. For eQTL from young post-mitotic dopaminergic neurons (differentiation of iPSC for 30 days, D30), we identified 34 genes whose genetically regulated expression in dopaminergic neurons may have a causal role in SCZ. Among which, ARL3 showed the most significant associations with SCZ. For eQTL from more mature dopaminergic neurons (D52), we identified 37 potential SCZ causal genes, and ARL3 and GNL3 showed the most significant associations. Only 12 genes showed significant associations with SCZ in both D30 and D52 eQTL datasets, indicating the time point-specific genetic regulatory effects in young post-mitotic dopaminergic neurons and more mature dopaminergic neurons. Comparing the results from dopaminergic neurons with bulk brain tissues prioritized 2 high-confidence risk genes, including DDHD2 and GALNT10. Our study identifies multiple risk genes whose genetically regulated expression in dopaminergic neurons may have a causal role in SCZ. Further mechanistic investigation will provide pivotal insights into SCZ pathophysiology.

The genetic architecture of schizophrenia: review of large-scale genetic studies

Schizophrenia is a complex and often chronic psychiatric disorder with high heritability. Diagnosis of schizophrenia is still made clinically based on psychiatric symptoms; no diagnostic tests or biomarkers are available. Pathophysiology-based diagnostic scheme and treatments are also not available. Elucidation of the pathogenesis is needed for development of pathology-based diagnostics and treatments. In the past few decades, genetic research has made substantial advances in our understanding of the genetic architecture of schizophrenia. Rare copy number variations (CNVs) and rare single-nucleotide variants (SNVs) detected by whole-genome CNV analysis and whole-genome/-exome sequencing analysis have provided the great advances. Common single-nucleotide polymorphisms (SNPs) detected by large-scale genome-wide association studies have also provided important information. Large-scale genetic studies have been revealed that both rare and common genetic variants play crucial roles in this disorder. In this review, we focused on CNVs, SNVs, and SNPs, and discuss the latest research findings on the pathogenesis of schizophrenia based on these genetic variants. Rare variants with large effect sizes can provide mechanistic hypotheses. CRISPR-based genetics approaches and induced pluripotent stem cell technology can facilitate the functional analysis of these variants detected in patients with schizophrenia. Recent advances in long-read sequence technology are expected to detect variants that cannot be detected by short-read sequence technology. Various studies that bring together data from common variant and transcriptomic datasets provide biological insight. These new approaches will provide additional insight into the pathophysiology of schizophrenia and facilitate the development of pathology-based therapeutics.

Harnessing Transcriptomic Signals for Amyotrophic Lateral Sclerosis to Identify Novel Drugs and Enhance Risk Prediction

Introduction

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease. This study integrates the latest ALS genome-wide association study (GWAS) summary statistics with functional genomic annotations with the aim of providing mechanistic insights into ALS risk loci, inferring drug repurposing opportunities, and enhancing prediction of ALS risk and clinical characteristics.

Methods

Genes associated with ALS were identified using GWAS summary statistic methodology including SuSiE SNP-based fine-mapping, and transcriptome- and proteome-wide association study (TWAS/PWAS) analyses. Using several approaches, gene associations were integrated with the DrugTargetor drug-gene interaction database to identify drugs that could be repurposed for the treatment of ALS. Furthermore, ALS gene associations from TWAS were combined with observed blood expression in two external ALS case-control datasets to calculate polytranscriptomic scores and evaluate their utility for prediction of ALS risk and clinical characteristics, including site of onset, age at onset, and survival.

Results

SNP-based fine-mapping, TWAS and PWAS identified 117 genes associated with ALS, with TWAS and PWAS providing novel mechanistic insights. Drug repurposing analyses identified five drugs significantly enriched for interactions with ALS associated genes, with directional analyses highlighting α-glucosidase inhibitors may exacerbate ALS pathology. Additionally, drug class enrichment analysis showed calcium channel blockers may reduce ALS risk. Across the two observed expression target samples, ALS polytranscriptomic scores significantly predicted ALS risk (R2 = 4%; p-value = 2.1×10-21).

Conclusions

Functionally-informed analyses of ALS GWAS summary statistics identified novel mechanistic insights into ALS aetiology, highlighted several therapeutic research avenues, and enabled statistically significant prediction of ALS risk.

Gene co-expression architecture in peripheral blood in a cohort of remitted first-episode schizophrenia patients

A better understanding of schizophrenia subtypes is necessary to stratify the patients according to clinical attributes. To explore the genomic architecture of schizophrenia symptomatology, we analyzed blood co-expression modules and their association with clinical data from patients in remission after a first episode of schizophrenia. In total, 91 participants of the 2EPS project were included. Gene expression was assessed using the Clariom S Human Array. Weighted-gene co-expression network analysis (WGCNA) was applied to identify modules of co-expressed genes and to test its correlation with global functioning, clinical symptomatology, and premorbid adjustment. Among the 25 modules identified, six modules were significantly correlated with clinical data. These modules could be clustered in two groups according to their correlation with clinical data. Hub genes in each group showing overlap with risk genes for schizophrenia were enriched in biological processes related to metabolic processes, regulation of gene expression, cellular localization and protein transport, immune processes, and neurotrophin pathways. Our results indicate that modules with significant associations with clinical data showed overlap with gene sets previously identified in differential gene-expression analysis in brain, indicating that peripheral tissues could reveal pathogenic mechanisms. Hub genes involved in these modules revealed multiple signaling pathways previously related to schizophrenia, which may represent the complex interplay in the pathological mechanisms behind the disease. These genes could represent potential targets for the development of peripheral biomarkers underlying illness traits in clinical remission stages after a first episode of schizophrenia.

Identification of Novel Metabolic Subtypes Using Multi-Trait Limited Mixed Regression in the Chinese Population

The aggregation and interaction of metabolic risk factors leads to highly heterogeneous pathogeneses, manifestations, and outcomes, hindering risk stratification and targeted management. To deconstruct the heterogeneity, we used baseline data from phase II of the Fangshan Family-Based Ischemic Stroke Study (FISSIC), and a total of 4632 participants were included. A total of 732 individuals who did not have any component of metabolic syndrome (MetS) were set as a reference group, while 3900 individuals with metabolic abnormalities were clustered into subtypes using multi-trait limited mixed regression (MFMR). Four metabolic subtypes were identified with the dominant characteristics of abdominal obesity, hypertension, hyperglycemia, and dyslipidemia. Multivariate logistic regression showed that the hyperglycemia-dominant subtype had the highest coronary heart disease (CHD) risk (OR: 6.440, 95% CI: 3.177-13.977) and that the dyslipidemia-dominant subtype had the highest stroke risk (OR: 2.450, 95% CI: 1.250-5.265). Exome-wide association studies (EWASs) identified eight SNPs related to the dyslipidemia-dominant subtype with genome-wide significance, which were located in the genes APOA5, BUD13, ZNF259, and WNT4. Functional analysis revealed an enrichment of top genes in metabolism-related biological pathways and expression in the heart, brain, arteries, and kidneys. Our findings provide directions for future attempts at risk stratification and evidence-based management in populations with metabolic abnormalities from a systematic perspective.

Identification of genetic loci that overlap between schizophrenia and metabolic syndrome

Patients with schizophrenia (SCZ) frequently exhibit an elevated risk of metabolic syndrome (MetS), which may lead to a worse clinical outcome. Even though these two phenotypes are genetically linked, little is known about their shared genetic determinants. Here, we investigated whether SCZ and MetS share genetic risk factors. To examine the genetic overlap between the two disorders, we applied a comprehensive genetic overlap analysis by integrating GWAS data for SCZ (n = 320,404) and MetS (n = 291,107) at the genome, genetic variants, and gene levels. At the genome level, we observed polygenic overlap between SCZ and MetS by utilizing LDSC (r_g=-0.09, P = 1 × 10^-4) and GNOVA (rho=-0.04, P = 1.39 × 10^-8) analysis. At the SNP level, we performed conjunctional FDR (conjFDR) analysis to identify genetic variants simultaneously associated with two phenotypes. Based on conjFDR < 0.05, we identified 22 loci shared between SCZ and MetS. At the gene level, we further demonstrated that SCZ- and MetS-inferred gene expression overlapped across 49 GTEx tissues and highlighted the PCCB and KCTD13 genes as putative mediators of the genetic association. Overall, these findings shed novel light on the association between SCZ and MetS, and potentially enhance our knowledge of the high comorbidity and genetic processes that overlap between the two disorders.

Bioinformatics detection of modulators controlling splicing factor-dependent intron retention in the human brain

Alternative RNA splicing is an important means of genetic control and transcriptome diversity. However, when alternative splicing events are studied independently, coordinated splicing modulated by common factors is often not recognized. As a result, the molecular mechanisms of how splicing regulators promote or repress splice site recognition in a context-dependent manner are not well understood. The functional coupling between multiple gene regulatory layers suggests that splicing is modulated by additional genetic or epigenetic components. Here, we developed a bioinformatics approach to identify causal modulators of splicing activity based on the variation of gene expression in large RNA sequencing datasets. We applied this approach in a neurological context with hundreds of dorsolateral prefrontal cortex samples. Our model is strengthened with the incorporation of genetic variants to impute gene expression in a Mendelian randomization-based approach. We identified novel modulators of the splicing factor SRSF1, including UIMC1 and the long noncoding RNA CBR3-AS1, that function over dozens of SRSF1 intron retention splicing targets. This strategy can be widely used to identify modulators of RNA-binding proteins involved in tissue-specific alternative splicing.

Lessons Learned From Parsing Genetic Risk for Schizophrenia Into Biological Pathways

The clinically heterogeneous presentation of schizophrenia is compounded by the heterogeneity of risk factors and neurobiological correlates of the disorder. Genome-wide association studies in schizophrenia have uncovered a remarkably high number of genetic variants, but the biological pathways they impact upon remain largely unidentified. Among the diverse methodological approaches employed to provide a more granular understanding of genetic risk for schizophrenia, the use of biological labels, such as gene ontologies, regulome approaches, and gene coexpression have all provided novel perspectives into how genetic risk translates into the neurobiology of schizophrenia. Here, we review the salient aspects of parsing polygenic risk for schizophrenia into biological pathways. We argue that parsed scores, compared to standard polygenic risk scores, may afford a more biologically plausible and accurate physiological modeling of the different dimensions involved in translating genetic risk into brain mechanisms, including multiple brain regions, cell types, and maturation stages. We discuss caveats, opportunities, and pitfalls inherent in the parsed risk approach.

Upper cortical layer-driven network impairment in schizophrenia

Schizophrenia is one of the most widespread and complex mental disorders. To characterize the impact of schizophrenia, we performed single-nucleus RNA sequencing (snRNA-seq) of >220,000 neurons from the dorsolateral prefrontal cortex of patients with schizophrenia and matched controls. In addition, >115,000 neurons were analyzed topographically by immunohistochemistry. Compositional analysis of snRNA-seq data revealed a reduction in abundance of GABAergic neurons and a concomitant increase in principal neurons, most pronounced for upper cortical layer subtypes, which was substantiated by histological analysis. Many neuronal subtypes showed extensive transcriptomic changes, the most marked in upper-layer GABAergic neurons, including down-regulation in energy metabolism and up-regulation in neurotransmission. Transcription factor network analysis demonstrated a developmental origin of transcriptomic changes. Last, Visium spatial transcriptomics further corroborated upper-layer neuron vulnerability in schizophrenia. Overall, our results point toward general network impairment within upper cortical layers as a core substrate associated with schizophrenia symptomatology.

A novel genetic variant potentially altering the expression of MANBA in the cerebellum associated with attention deficit hyperactivity disorder in Han Chinese children

OBJECTIVES

To obtain additional insight into the genetic factors of attention deficit hyperactivity disorder (ADHD).

METHODS

First, we performed a transcriptome-wide association study (TWAS) integrating human cerebellum-specific variant-expression/splicing correlations to identify ADHD susceptibility genes. Then, the associations between expression/splicing quantitative trait loci (eQTLs/sQTLs) of the transcriptome-wide significant genes and ADHD were observed in a case-control study of Han Chinese children. Furthermore, dual luciferase reporter gene assays were performed to validate the regulatory function of ADHD risk variants. Additionally, the transcription level of target genes in blood was detected by real-time quantitative polymerase chain reaction (RT-qPCR) assay.

RESULTS

TWAS identified that the genetically regulated expression of MANBA in the cerebellum was significantly associated with ADHD risk. Furthermore, we observed a higher risk of ADHD and more severe clinical symptoms in subjects harbouring heterozygous (TC) or mutant homozygous (TT) genotypes of MANBA rs1054037 than CC carriers. The dual luciferase reporter gene assay revealed that the mutation of rs1054037(C > T) potentially upregulated MANBA expression by eliminating the binding site for hsa-miR-5591-3P. Finally, RT-qPCR showed that MANBA expression in blood samples of patients was significantly higher than that of controls.

CONCLUSIONS

Taken together, these results suggest a role of MANBA in the development of ADHD.

The shared genetic basis of mood instability and psychiatric disorders: A cross-trait genome-wide association analysis

Recent genome-wide association studies of mood instability (MOOD) have found significant positive genetic correlation with major depression (DEP) and weak correlations with other psychiatric disorders. We investigated the polygenic overlap between MOOD and psychiatric disorders beyond genetic correlation to better characterize putative shared genetic determinants. GWAS summary statistics for schizophrenia (SCZ, n = 105,318), bipolar disorder (BIP, n = 413,466), DEP (n = 450,619), attention-deficit hyperactivity disorder (ADHD, n = 53,293), and MOOD (n = 363,705) were analyzed using the bivariate causal mixture model and conjunctional false discovery rate methods. MOOD correlated positively with all psychiatric disorders, but with wide variation in strength (r_g  = 0.10-0.62). Of 10.4 K genomic variants influencing MOOD, 4 K-9.4 K influenced psychiatric disorders. Furthermore, MOOD was jointly associated with DEP at 163 loci, SCZ at 110, BIP at 60 and ADHD at 25. Fifty-three jointly associated loci were overlapping across two or more disorders, seven of which had discordant effect directions on psychiatric disorders. Genes mapped to loci associated with MOOD and all four disorders were enriched in a single gene-set, "synapse organization." The extensive polygenic overlap indicates shared molecular underpinnings across MOOD and psychiatric disorders. However, distinct patterns of genetic correlation and effect directions may relate to differences in the core clinical features of each disorder.

Genetic control of RNA splicing and its distinct role in complex trait variation

Most genetic variants identified from genome-wide association studies (GWAS) in humans are noncoding, indicating their role in gene regulation. Previous studies have shown considerable links of GWAS signals to expression quantitative trait loci (eQTLs) but the links to other genetic regulatory mechanisms, such as splicing QTLs (sQTLs), are underexplored. Here, we introduce an sQTL mapping method, testing for heterogeneity between isoform-eQTLeffects (THISTLE), with improved power over competing methods. Applying THISTLE together with a complementary sQTL mapping strategy to brain transcriptomic (n = 2,865) and genotype data, we identified 12,794 genes with cis-sQTLs at P < 5 × 10⁻⁸, approximately 61% of which were distinct from eQTLs. Integrating the sQTL data into GWAS for 12 brain-related complex traits (including diseases), we identified 244 genes associated with the traits through cis-sQTLs, approximately 61% of which could not be discovered using the corresponding eQTL data. Our study demonstrates the distinct role of most sQTLs in the genetic regulation of transcription and complex trait variation.

A powerful method for splicing quantitative trait loci (sQTL) mapping, THISTLE, is presented and applied to a collection of 2,865 brain samples. Integration with GWAS identifies 244 genes associated via cis-sQTLs, of which 61% were not identified using expression QTLs.

Shared mechanisms across the major psychiatric and neurodegenerative diseases

Several common psychiatric and neurodegenerative diseases share epidemiologic risk; however, whether they share pathophysiology is unclear and is the focus of our investigation. Using 25 GWAS results and LD score regression, we find eight significant genetic correlations between psychiatric and neurodegenerative diseases. We integrate the GWAS results with human brain transcriptomes (n = 888) and proteomes (n = 722) to identify cis- and trans- transcripts and proteins that are consistent with a pleiotropic or causal role in each disease, referred to as causal proteins for brevity. Within each disease group, we find many distinct and shared causal proteins. Remarkably, 30% (13 of 42) of the neurodegenerative disease causal proteins are shared with psychiatric disorders. Furthermore, we find 2.6-fold more protein-protein interactions among the psychiatric and neurodegenerative causal proteins than expected by chance. Together, our findings suggest these psychiatric and neurodegenerative diseases have shared genetic and molecular pathophysiology, which has important ramifications for early treatment and therapeutic development.

Integrating human brain proteomes with genome-wide association data implicates novel proteins in post-traumatic stress disorder

Genome-wide association studies (GWAS) have identified several risk loci for post-traumatic stress disorder (PTSD); however, how they confer PTSD risk remains unclear. We aimed to identify genes that confer PTSD risk through their effects on brain protein abundance to provide new insights into PTSD pathogenesis. To that end, we integrated human brain proteomes with PTSD GWAS results to perform a proteome-wide association study (PWAS) of PTSD, followed by Mendelian randomization, using a discovery and confirmatory study design. Brain proteomes (N = 525) were profiled from the dorsolateral prefrontal cortex using mass spectrometry. The Million Veteran Program (MVP) PTSD GWAS (n = 186,689) was used for the discovery PWAS, and the Psychiatric Genomics Consortium PTSD GWAS (n = 174,659) was used for the confirmatory PWAS. To understand whether genes identified at the protein-level were also evident at the transcript-level, we performed a transcriptome-wide association study (TWAS) using human brain transcriptomes (N = 888) and the MVP PTSD GWAS results. We identified 11 genes that contribute to PTSD pathogenesis via their respective cis-regulated brain protein abundance. Seven of 11 genes (64%) replicated in the confirmatory PWAS and 4 of 11 also had their cis-regulated brain mRNA levels associated with PTSD. High confidence level was assigned to 9 of 11 genes after considering evidence from the confirmatory PWAS and TWAS. Most of the identified genes are expressed in other PTSD-relevant brain regions and several are preferentially expressed in excitatory neurons, astrocytes, and oligodendrocyte precursor cells. These genes are novel, promising targets for mechanistic and therapeutic studies to find new treatments for PTSD.

Powerful and robust inference of complex phenotypes' causal genes with dependent expression quantitative loci by a median-based Mendelian randomization

Isolating the causal genes from numerous genetic association signals in genome-wide association studies (GWASs) of complex phenotypes remains an open and challenging question. In the present study, we proposed a statistical approach, the effective-median-based Mendelian randomization (MR) framework, for inferring the causal genes of complex phenotypes with the GWAS summary statistics (named EMIC). The effective-median method solved the high false-positive issue in the existing MR methods due to either correlation among instrumental variables or noises in approximated linkage disequilibrium (LD). EMIC can further perform a pleiotropy fine-mapping analysis to remove possible false-positive estimates. With the usage of multiple cis-expression quantitative trait loci (eQTLs), EMIC was also more powerful than the alternative methods for the causal gene inference in the simulated datasets. Furthermore, EMIC rediscovered many known causal genes of complex phenotypes (schizophrenia, bipolar disorder, and total cholesterol) and reported many new and promising candidate causal genes. In sum, this study provided an efficient solution to discriminate the candidate causal genes from vast amounts of GWAS signals with eQTLs. EMIC has been implemented in our integrative software platform KGGSEE.

An epigenetic association analysis of childhood trauma in psychosis reveals possible overlap with methylation changes associated with PTSD

Patients with a severe mental disorder report significantly higher levels of childhood trauma (CT) than healthy individuals. Studies have suggested that CT may affect brain plasticity through epigenetic mechanisms and contribute to developing various psychiatric disorders. We performed a blood-based epigenome-wide association study using the Childhood Trauma Questionnaire-short form in 602 patients with a current severe mental illness, investigating DNA methylation association separately for five trauma subtypes and the total trauma score. The median trauma score was set as the predefined cutoff for determining whether the trauma was present or not. Additionally, we compared our genome-wide results with methylation probes annotated to candidate genes previously associated with CT. Of the patients, 83.2% reported CT above the cutoff in one or more trauma subtypes, and emotional neglect was the trauma subtype most frequently reported. We identified one significant differently methylated position associated with the gene TANGO6 for physical neglect. Seventeen differentially methylated regions (DMRs) were associated with different trauma categories. Several of these DMRs were annotated to genes previously associated with neuropsychiatric disorders such as post-traumatic stress disorder and cognitive impairments. Our results support a biomolecular association between CT and severe mental disorders. Genes that were previously identified as differentially methylated in CT-exposed subjects with and without psychosis did not show methylation differences in our analysis. We discuss this inconsistency, the relevance of our findings, and the limitations of our study.

A conditional gene-based association framework integrating isoform-level eQTL data reveals new susceptibility genes for schizophrenia

Linkage disequilibrium and disease-associated variants in the non-coding regions make it difficult to distinguish the truly associated genes from the redundantly associated genes for complex diseases. In this study, we proposed a new conditional gene-based framework called eDESE that leveraged an improved effective chi-squared statistic to control the type I error rates and remove the redundant associations. eDESE initially performed the association analysis by mapping variants to genes according to their physical distance. We further demonstrated that the isoform-level eQTLs could be more powerful than the gene-level eQTLs in the association analysis using a simulation study. Then the eQTL-guided strategies, that is, mapping variants to genes according to their gene/isoform-level variant-gene cis-eQTLs associations, were also integrated with eDESE. We then applied eDESE to predict the potential susceptibility genes of schizophrenia and found that the potential susceptibility genes were enriched with many neuronal or synaptic signaling-related terms in the Gene Ontology knowledgebase and antipsychotics-gene interaction terms in the drug-gene interaction database (DGIdb). More importantly, seven potential susceptibility genes identified by eDESE were the target genes of multiple antipsychotics in DrugBank. Comparing the potential susceptibility genes identified by eDESE and other benchmark approaches (i.e., MAGMA and S-PrediXcan) implied that strategy based on the isoform-level eQTLs could be an important supplement for the other two strategies (physical distance and gene-level eQTLs). We have implemented eDESE in our integrative platform KGGSEE (http://pmglab.top/kggsee/#/) and hope that eDESE can facilitate the prediction of candidate susceptibility genes and isoforms for complex diseases in a multi-tissue context.

Predicted gene expression in ancestrally diverse populations leads to discovery of susceptibility loci for lifestyle and cardiometabolic traits

One mechanism by which genetic factors influence complex traits and diseases is altering gene expression. Direct measurement of gene expression in relevant tissues is rarely tenable; however, genetically regulated gene expression (GReX) can be estimated using prediction models derived from large multi-omic datasets. These approaches have led to the discovery of many gene-trait associations, but whether models derived from predominantly European ancestry (EA) reference panels can map novel associations in ancestrally diverse populations remains unclear. We applied PrediXcan to impute GReX in 51,520 ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) participants (35% African American, 45% Hispanic/Latino, 10% Asian, and 7% Hawaiian) across 25 key cardiometabolic traits and relevant tissues to identify 102 novel associations. We then compared associations in PAGE to those in a random subset of 50,000 White British participants from UK Biobank (UKBB50k) for height and body mass index (BMI). We identified 517 associations across 47 tissues in PAGE but not UKBB50k, demonstrating the importance of diverse samples in identifying trait-associated GReX. We observed that variants used in PrediXcan models were either more or less differentiated across continental-level populations than matched-control variants depending on the specific population reflecting sampling bias. Additionally, variants from identified genes specific to either PAGE or UKBB50k analyses were more ancestrally differentiated than those in genes detected in both analyses, underlining the value of population-specific discoveries. This suggests that while EA-derived transcriptome imputation models can identify new associations in non-EA populations, models derived from closely matched reference panels may yield further insights. Our findings call for more diversity in reference datasets of tissue-specific gene expression.

Transcriptome-wide association study reveals increased neuronal FLT3 expression is associated with Tourette's syndrome

Tourette's Syndrome (TS) is a neurodevelopmental disorder that is characterized by motor and phonic tics. A recent TS genome-wide association study (GWAS) identified a genome-wide significant locus. However, determining the biological mechanism of GWAS signals remains difficult. To characterize effects of expression quantitative trait loci (eQTLs) in TS and understand biological underpinnings of the disease. Here, we conduct a TS transcriptome-wide association study (TWAS) consisting of 4819 cases and 9488 controls. We demonstrate that increased expression of FLT3 in the dorsolateral prefrontal cortex (DLPFC) is associated with TS. We further show that there is global dysregulation of FLT3 across several brain regions and probabilistic causal fine-mapping of the TWAS signal prioritizes FLT3 with a posterior inclusion probability of 0.849. After, we proxy the expression with 100 lymphoblastoid cell lines, and demonstrate that TS cells has a 1.72 increased fold change compared to controls. A phenome-wide association study also points toward FLT3 having links with immune-related pathways such as monocyte count. We further identify several splicing events in MPHOSPH9, CSGALNACT2 and FIP1L1 associated with TS, which are also implicated in immune function. This analysis of expression and splicing begins to explore the biology of TS GWAS signals.

Bioinformatics and Network-based Approaches for Determining Pathways, Signature Molecules, and Drug Substances connected to Genetic Basis of Schizophrenia etiology

Knowledge of heterogeneous etiology and pathophysiology of schizophrenia (SZP) is reasonably inadequate and non-deterministic due to its inherent complexity and underlying vast dynamics related to genetic mechanisms. The evolution of large-scale transcriptome-wide datasets and subsequent development of relevant, robust technologies for their analyses show promises toward elucidating the genetic basis of disease pathogenesis, its early risk prediction, and predicting drug molecule targets for therapeutic intervention. In this research, we have scrutinized the genetic basis of SZP through functional annotation and network-based system biology approaches. We have determined 96 overlapping differentially expressed genes (DEGs) from 2 microarray datasets and subsequently identified their interconnecting networks to reveal transcriptome signatures like hub proteins (FYN, RAD51, SOCS3, XIAP, AKAP13, PIK3C2A, CBX5, GATA3, EIF3K, and CDKN2B), transcription factors and miRNAs. In addition, we have employed gene set enrichment to highlight significant gene ontology (e.g., positive regulation of microglial cell activation) and relevant pathways (such as axon guidance and focal adhesion) interconnected to the genes associated with SZP. Finally, we have suggested candidate drug substances like Luteolin HL60 UP as a possible therapeutic target based on these key molecular signatures.

Anti-MOG autoantibody-associated schizophreniform psychosis

OBJECTIVES

Autoimmune mechanisms are related to disease development in a subgroup of patients with psychosis. The contribution of immunoglobulin G (IgG) antibodies against myelin oligodendrocyte glycoprotein (MOG) is mainly unclear in this context.

METHODS

Therefore, two patients with psychosis and anti-MOG antibodies - detected in fixed cell-based and live cell-based assays - are presented.

RESULTS

Patient 1 suffered from late-onset psychosis with singular white matter lesions in magnetic resonance imaging (MRI) and intermittent electroencephalography (EEG) slowing. Patient 2 suffered from a chronic paranoid-hallucinatory disorder with intermittent confusional states, non-specific white matter alterations on MRI, a disorganised alpha rhythm on EEG, and elevated cerebrospinal fluid protein. Both patients had anti-MOG antibody titres of 1 : 320 in serum (reference < 1 : 20).

CONCLUSIONS

The arguments for and against a causal role for anti-MOG antibodies are discussed. The antibodies could be relevant, but due to moderate titres, they may have caused a rather 'subtle clinical picture' consisting of psychosis instead of 'classical' MOG encephalomyelitis.

What genes are differentially expressed in individuals with schizophrenia? A systematic review

Schizophrenia is a severe, complex mental disorder characterized by a combination of positive symptoms, negative symptoms, and impaired cognitive function. Schizophrenia is highly heritable (~80%) with multifactorial etiology and complex polygenic genetic architecture. Despite the large number of genetic variants associated with schizophrenia, few causal variants have been established. Gaining insight into the mechanistic influences of these genetic variants may facilitate our ability to apply these findings to prevention and treatment. Though there have been more than 300 studies of gene expression in schizophrenia over the past 15 years, none of the studies have yielded consistent evidence for specific genes that contribute to schizophrenia risk. The aim of this work is to conduct a systematic review and synthesis of case-control studies of genome-wide gene expression in schizophrenia. Comprehensive literature searches were completed in PubMed, EmBase, and Web of Science, and after a systematic review of the studies, data were extracted from those that met the following inclusion criteria: human case-control studies comparing the genome-wide transcriptome of individuals diagnosed with schizophrenia to healthy controls published between January 1, 2000 and June 30, 2020 in the English language. Genes differentially expressed in cases were extracted from these studies, and overlapping genes were compared to previous research findings from the genome-wide association, structural variation, and tissue-expression studies. The transcriptome-wide analysis identified different genes than those previously reported in genome-wide association, exome sequencing, and structural variation studies of schizophrenia. Only one gene, GBP2, was replicated in five studies. Previous work has shown that this gene may play a role in immune function in the etiology of schizophrenia, which in turn could have implications for risk profiling, prevention, and treatment. This review highlights the methodological inconsistencies that impede valid meta-analyses and synthesis across studies. Standardization of the use of covariates, gene nomenclature, and methods for reporting results could enhance our understanding of the potential mechanisms through which genes exert their influence on the etiology of schizophrenia. Although these results are promising, collaborative efforts with harmonization of methodology will facilitate the identification of the role of genes underlying schizophrenia.