Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection

Genetic correlation and causal associations between psychiatric disorders and lung cancer risk

BACKGROUND

Patients with certain psychiatric disorders have increased lung cancer incidence. However, establishing a causal relationship through traditional epidemiological methods poses challenges.

METHODS

Available summary statistics of genome-wide association studies of cigarette smoking, lung cancer, and eight psychiatric disorders, including attention deficit/hyperactivity disorder (ADHD), autism, depression, major depressive disorder, bipolar disorder, insomnia, neuroticism, and schizophrenia (range N: 46,350-1,331,010) were leveraged to estimate genetic correlations using Linkage Disequilibrium Score Regression and assess causal effect of each psychiatric disorder on lung cancer using two-sample Mendelian randomization (MR) models, comprising inverse-variance weighted (IVW), weighted median, MR-Egger, pleiotropy residual sum and outlier testing (MR-PRESSO), and a constrained maximum likelihood approach (cML-MR).

RESULTS

Significant positive correlations were observed between each psychiatric disorder and both smoking and lung cancer (all FDR < 0.05), except for the correlation between autism and lung cancer. Both univariable and the cML-MA MR analyses demonstrated that liability to schizophrenia, depression, ADHD, or insomnia was associated with an increased risk of overall lung cancer. Genetic liability to insomnia was linked specifically to squamous cell carcinoma (SCC), while genetic liability to ADHD was associated with an elevated risk of both SCC and small cell lung cancer (all P < 0.05). The later was further supported by multivariable MR analyses, which accounted for smoking.

LIMITATIONS

Participants were constrained to European ancestry populations. Causal estimates from binary psychiatric disorders may be biased.

CONCLUSION

Our findings suggest appropriate management of several psychiatric disorders, particularly ADHD, may potentially reduce the risk of developing lung cancer.

The role and molecular mechanisms of the early growth response 3 gene in schizophrenia

Schizophrenia is a chronic, debilitating mental illness caused by both genetic and environmental factors. Genetic factors play a major role in schizophrenia development. Early growth response 3 (EGR3) is a member of the EGR family, which is associated with schizophrenia. Accumulating studies have investigated the relationship between EGR3 and schizophrenia. However, the role of EGR3 in schizophrenia pathogenesis remains unclear. In the present review, we focus on the progress of research related to the role of EGR3 in schizophrenia, including association studies between EGR3 and schizophrenia, abnormal gene expressional analysis of EGR3 in schizophrenia, biological function studies of EGR3 in schizophrenia, the molecular regulatory mechanism of EGR3 and schizophrenia susceptibility candidate genes, and possible role of EGR3 in the immune system function in schizophrenia. In summary, EGR3 is a schizophrenia risk candidate factor and has comprehensive regulatory roles in schizophrenia pathogenesis. Further studies investigating the molecular mechanisms of EGR3 in schizophrenia are warranted for understanding the pathophysiology of this disorder as well as the development of new therapeutic strategies for the treatment and control of this disorder.

Discovery of 3-((4-Benzylpyridin-2-yl)amino)benzamides as Potent GPR52 G Protein-Biased Agonists

Orphan GPR52 is emerging as a promising neurotherapeutic target. Optimization of previously reported lead 4a employing an iterative drug design strategy led to the identification of a series of unique GPR52 agonists, such as 10a (PW0677), 15b (PW0729), and 24f (PW0866), with improved potency and efficacy. Intriguingly, compounds 10a and 24f showed greater bias for G protein/cAMP signaling and induced significantly less in vitro desensitization than parent compound 4a, indicating that reducing GPR52 β-arrestin activity with biased agonism results in sustained GPR52 activation. Further exploration of compounds 15b and 24f indicated improved potency and efficacy, and excellent target selectivity, but limited brain exposure warranting further optimization. These balanced and biased GPR52 agonists provide important pharmacological tools to study GPR52 activation, signaling bias, and therapeutic potential for neuropsychiatric and neurological diseases.

Unraveling NEK4 as a Potential Drug Target in Schizophrenia and Bipolar I Disorder: A Proteomic and Genomic Approach

BACKGROUND AND HYPOTHESIS

Investigating the shared brain protein and genetic components of schizophrenia (SCZ) and bipolar I disorder (BD-I) presents a unique opportunity to understand the underlying pathophysiological processes and pinpoint potential drug targets.

STUDY DESIGN

To identify overlapping susceptibility brain proteins in SCZ and BD-I, we carried out proteome-wide association studies (PWAS) and Mendelian Randomization (MR) by integrating human brain protein quantitative trait loci with large-scale genome-wide association studies for both disorders. We utilized transcriptome-wide association studies (TWAS) to determine the consistency of mRNA-protein dysregulation in both disorders. We applied pleiotropy-informed conditional false discovery rate (pleioFDR) analysis to identify common risk genetic loci for SCZ and BD-I. Additionally, we performed a cell-type-specific analysis in the human brain to detect risk genes notably enriched in distinct brain cell types. The impact of risk gene overexpression on dendritic arborization and axon length in neurons was also examined.

STUDY RESULTS

Our PWAS identified 42 proteins associated with SCZ and 14 with BD-I, among which NEK4, HARS2, SUGP1, and DUS2 were common to both conditions. TWAS and MR analysis verified the significant risk gene NEK4 for both SCZ and BD-I. PleioFDR analysis further supported genetic risk loci associated with NEK4 for both conditions. The cell-type specificity analysis revealed that NEK4 is expressed on the surface of glutamatergic neurons, and its overexpression enhances dendritic arborization and axon length in cultured primary neurons.

CONCLUSIONS

These findings underscore a shared genetic origin for SCZ and BD-I, offering novel insights for potential therapeutic target identification.

Influence of antipsychotic drugs on microRNA expression in schizophrenia patients - A systematic review

Schizophrenia (SCZ) is a severe psychiatric disorder with unclear pathophysiology. Moreover, there is no specific biological marker to help clinicians to define a diagnosis, and medication is decided according to the psychiatrist's experience. In this scenario, microRNAs (miRNAs), which are small noncoding RNA molecules that regulate several genes, emerge as potential peripheral biomarkers to help not only the evaluation of the disease state but also the treatment response. Here, we systematically reviewed indexed literature and evaluated follow-up studies investigating the changes in miRNA expression due to antipsychotic treatment. We also assessed target genes and performed pathway enrichment analysis of miRNAs listed in this systematic review. A total of 11 studies were selected according to research criteria, and we observed that 28 miRNAs play a relevant role in schizophrenia pathogenesis or response to antipsychotic treatment, seven of those of extreme interest as possible biomarkers either for condition or treatment. Predicted targets of the miRNAs reviewed here were previously associated with schizophrenia in genome-wide studies, and pathway analysis showed enrichment for genes related to neural processes. With this review, we expect to highlight the importance of miRNAs in schizophrenia pathogenesis and its treatment and point out interesting miRNAs to future studies.

Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

Identification of important gene signatures in schizophrenia through feature fusion and genetic algorithm

Schizophrenia is a debilitating psychiatric disorder that can significantly affect a patient's quality of life and lead to permanent brain damage. Although medical research has identified certain genetic risk factors, the specific pathogenesis of the disorder remains unclear. Despite the prevalence of research employing magnetic resonance imaging, few studies have focused on the gene level and gene expression profile involving a large number of screened genes. However, the high dimensionality of genetic data presents a great challenge to accurately modeling the data. To tackle the current challenges, this study presents a novel feature selection strategy that utilizes heuristic feature fusion and a multi-objective optimization genetic algorithm. The goal is to improve classification performance and identify the key gene subset for schizophrenia diagnostics. Traditional gene screening techniques are inadequate for accurately determining the precise number of key genes associated with schizophrenia. Our innovative approach integrates a filter-based feature selection method to reduce data dimensionality and a multi-objective optimization genetic algorithm for improved classification tasks. By combining the filtering and wrapper methods, our strategy leverages their respective strengths in a deliberate manner, leading to superior classification accuracy and a more efficient selection of relevant genes. This approach has demonstrated significant improvements in classification results across 11 out of 14 relevant datasets. The performance on the remaining three datasets is comparable to the existing methods. Furthermore, visual and enrichment analyses have confirmed the practicality of our proposed method as a promising tool for the early detection of schizophrenia.

NRN1 epistasis with BDNF and CACNA1C: mediation effects on symptom severity through neuroanatomical changes in schizophrenia

The expression of Neuritin-1 (NRN1), a neurotrophic factor crucial for neurodevelopment and synaptic plasticity, is enhanced by the Brain Derived Neurotrophic Factor (BDNF). Although the receptor of NRN1 remains unclear, it is suggested that NRN1's activation of the insulin receptor (IR) pathway promotes the transcription of the calcium voltage-gated channel subunit alpha1 C (CACNA1C). These three genes have been independently associated with schizophrenia (SZ) risk, symptomatology, and brain differences. However, research on how they synergistically modulate these phenotypes is scarce. We aimed to study whether the genetic epistasis between these genes affects the risk and clinical presentation of the disorder via its effect on brain structure. First, we tested the epistatic effect of NRN1 and BDNF or CACNA1C on (i) the risk for SZ, (ii) clinical symptoms severity and functionality (onset, PANSS, CGI and GAF), and (iii) brain cortical structure (thickness, surface area and volume measures estimated using FreeSurfer) in a sample of 86 SZ patients and 89 healthy subjects. Second, we explored whether those brain clusters influenced by epistatic effects mediate the clinical profiles. Although we did not find a direct epistatic impact on the risk, our data unveiled significant effects on the disorder's clinical presentation. Specifically, the NRN1-rs10484320 x BDNF-rs6265 interplay influenced PANSS general psychopathology, and the NRN1-rs4960155 x CACNA1C-rs1006737 interaction affected GAF scores. Moreover, several interactions between NRN1 SNPs and BDNF-rs6265 significantly influenced the surface area and cortical volume of the frontal, parietal, and temporal brain regions within patients. The NRN1-rs10484320 x BDNF-rs6265 epistasis in the left lateral orbitofrontal cortex fully mediated the effect on PANSS general psychopathology. Our study not only adds clinical significance to the well-described molecular relationship between NRN1 and BDNF but also underscores the utility of deconstructing SZ into biologically validated brain-imaging markers to explore their mediation role in the path from genetics to complex clinical manifestation.

The FXR1 network acts as signaling scaffold for actomyosin remodeling

It is currently not known whether mRNAs fulfill structural roles in the cytoplasm. Here, we report the FXR1 network, an mRNA-protein (mRNP) network present throughout the cytoplasm, formed by FXR1-mediated packaging of exceptionally long mRNAs. These mRNAs serve as underlying condensate scaffold and concentrate FXR1 molecules. The FXR1 network contains multiple protein binding sites and functions as a signaling scaffold for interacting proteins. We show that it is necessary for RhoA signaling-induced actomyosin reorganization to provide spatial proximity between kinases and their substrates. Point mutations in FXR1, found in its homolog FMR1, where they cause Fragile X syndrome, disrupt the network. FXR1 network disruption prevents actomyosin remodeling-an essential and ubiquitous process for the regulation of cell shape, migration, and synaptic function. These findings uncover a structural role for cytoplasmic mRNA and show how the FXR1 RNA-binding protein as part of the FXR1 network acts as organizer of signaling reactions.

Brain cell-type shifts in Alzheimer's disease, autism, and schizophrenia interrogated using methylomics and genetics

Few neuropsychiatric disorders have replicable biomarkers, prompting high-resolution and large-scale molecular studies. However, we still lack consensus on a more foundational question: whether quantitative shifts in cell types-the functional unit of life-contribute to neuropsychiatric disorders. Leveraging advances in human brain single-cell methylomics, we deconvolve seven major cell types using bulk DNA methylation profiling across 1270 postmortem brains, including from individuals diagnosed with Alzheimer's disease, schizophrenia, and autism. We observe and replicate cell-type compositional shifts for Alzheimer's disease (endothelial cell loss), autism (increased microglia), and schizophrenia (decreased oligodendrocytes), and find age- and sex-related changes. Multiple layers of evidence indicate that endothelial cell loss contributes to Alzheimer's disease, with comparable effect size to APOE genotype among older people. Genome-wide association identified five genetic loci related to cell-type composition, involving plausible genes for the neurovascular unit (P2RX5 and TRPV3) and excitatory neurons (DPY30 and MEMO1). These results implicate specific cell-type shifts in the pathophysiology of neuropsychiatric disorders.

Evaluating performance and applications of sample-wise cell deconvolution methods on human brain transcriptomic data

Sample-wise deconvolution methods estimate cell-type proportions and gene expressions in bulk tissue samples, yet their performance and biological applications remain unexplored, particularly in human brain transcriptomic data. Here, nine deconvolution methods were evaluated with sample-matched data from bulk tissue RNA sequencing (RNA-seq), single-cell/nuclei (sc/sn) RNA-seq, and immunohistochemistry. A total of 1,130,767 nuclei per cells from 149 adult postmortem brains and 72 organoid samples were used. The results showed the best performance of dtangle for estimating cell proportions and bMIND for estimating sample-wise cell-type gene expressions. For eight brain cell types, 25,273 cell-type eQTLs were identified with deconvoluted expressions (decon-eQTLs). The results showed that decon-eQTLs explained more schizophrenia GWAS heritability than bulk tissue or single-cell eQTLs did alone. Differential gene expressions associated with Alzheimer's disease, schizophrenia, and brain development were also examined using the deconvoluted data. Our findings, which were replicated in bulk tissue and single-cell data, provided insights into the biological applications of deconvoluted data in multiple brain disorders.

The pivotal role of the X-chromosome in the genetic architecture of the human brain

Genes on the X-chromosome are extensively expressed in the human brain. However, little is known for the X-chromosome's impact on the brain anatomy, microstructure, and functional network. We examined 1,045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X-chromosome interactions, while proposing an atlas outlining dosage compensation (DC) for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we discovered unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X-chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Asian-European differentiation of schizophrenia-associated genes driven by admixture and natural selection

The European-centered genome-wide association studies of schizophrenia (SCZ) may not be well applied to non-European populations. We analyzed 1,592 reported SCZ-associated genes using the public genome data and found an overall higher Asian-European differentiation on the SCZ-associated variants than at the genome-wide level. Notable examples included 15 missense variants, a regulatory variant SLC5A10-rs1624825, and a damaging variant TSPAN18-rs1001292. Independent local adaptations in recent 25,000 years, after the Asian-European divergence, could have contributed to such genetic differentiation, as were identified at a missense mutation LTN1-rs57646126-A in Asians, and a non-risk allele ZSWIM6-rs72761442-G in Europeans. Altai-Neanderthal-derived alleles may have opposite effects on SCZ susceptibility between ancestries. Furthermore, adaptive introgression was detected on the non-risk haplotype at 1q21.2 in Europeans, while in Asians it was observed on the SCZ risk haplotype at 3p21.31 which is also potentially ultra-violet protective. This study emphasizes the importance of including more representative Asian samples in future SCZ studies.

Breaking down causes, consequences, and mediating effects of telomere length variation on human health

BACKGROUND

Telomeres form repeated DNA sequences at the ends of chromosomes, which shorten with each cell division. Yet, factors modulating telomere attrition and the health consequences thereof are not fully understood. To address this, we leveraged data from 326,363 unrelated UK Biobank participants of European ancestry.

RESULTS

Using linear regression and bidirectional univariable and multivariable Mendelian randomization (MR), we elucidate the relationships between leukocyte telomere length (LTL) and 142 complex traits, including diseases, biomarkers, and lifestyle factors. We confirm that telomeres shorten with age and show a stronger decline in males than in females, with these factors contributing to the majority of the 5.4% of LTL variance explained by the phenome. MR reveals 23 traits modulating LTL. Smoking cessation and high educational attainment associate with longer LTL, while weekly alcohol intake, body mass index, urate levels, and female reproductive events, such as childbirth, associate with shorter LTL. We also identify 24 traits affected by LTL, with risk for cardiovascular, pulmonary, and some autoimmune diseases being increased by short LTL, while longer LTL increased risk for other autoimmune conditions and cancers. Through multivariable MR, we show that LTL may partially mediate the impact of educational attainment, body mass index, and female age at childbirth on proxied lifespan.

CONCLUSIONS

Our study sheds light on the modulators, consequences, and the mediatory role of telomeres, portraying an intricate relationship between LTL, diseases, lifestyle, and socio-economic factors.

Enrichment of a subset of Neanderthal polymorphisms in autistic probands and siblings

Homo sapiens and Neanderthals underwent hybridization during the Middle/Upper Paleolithic age, culminating in retention of small amounts of Neanderthal-derived DNA in the modern human genome. In the current study, we address the potential roles Neanderthal single nucleotide polymorphisms (SNP) may be playing in autism susceptibility in samples of black non-Hispanic, white Hispanic, and white non-Hispanic people using data from the Simons Foundation Powering Autism Research (SPARK), Genotype-Tissue Expression (GTEx), and 1000 Genomes (1000G) databases. We have discovered that rare variants are significantly enriched in autistic probands compared to race-matched controls. In addition, we have identified 25 rare and common SNPs that are significantly enriched in autism on different ethnic backgrounds, some of which show significant clinical associations. We have also identified other SNPs that share more specific genotype-phenotype correlations but which are not necessarily enriched in autism and yet may nevertheless play roles in comorbid phenotype expression (e.g., intellectual disability, epilepsy, and language regression). These results strongly suggest Neanderthal-derived DNA is playing a significant role in autism susceptibility across major populations in the United States.

Structural characterization of a polymorphic repeat at the CACNA1C schizophrenia locus

Genetic variation within intron 3 of the CACNA1C calcium channel gene is associated with schizophrenia and bipolar disorder, but analysis of the causal variants and their effect is complicated by a nearby variable-number tandem repeat (VNTR). Here, we used 155 long-read genome assemblies from 78 diverse individuals to delineate the structure and population variability of the CACNA1C intron 3 VNTR. We categorized VNTR sequences into 7 Types of structural alleles using sequence differences among repeat units. Only 12 repeat units at the 5' end of the VNTR were shared across most Types, but several Types were related through a series of large and small duplications. The most diverged Types were rare and present only in individuals with African ancestry, but the multiallelic structural polymorphism Variable Region 2 was present across populations at different frequencies, consistent with expansion of the VNTR preceding the emergence of early hominins. VR2 was in complete linkage disequilibrium with fine-mapped schizophrenia variants (SNPs) from genome-wide association studies (GWAS). This risk haplotype was associated with decreased CACNA1C gene expression in brain tissues profiled by the GTEx project. Our work suggests that sequence variation within a human-specific VNTR affects gene expression, and provides a detailed characterization of new alleles at a flagship neuropsychiatric locus.

Functional classes of SNPs related to psychiatric disorders and behavioral traits contrast with those related to neurological disorders

We investigated the functional classes of genomic regions containing SNPS contributing most to the SNP-heritability of important psychiatric and neurological disorders and behavioral traits, as determined from recent genome-wide association studies. We employed linkage-disequilibrium score regression with several brain-specific genomic annotations not previously utilized. The classes of genomic annotations conferring substantial SNP-heritability for the psychiatric disorders and behavioral traits differed systematically from the classes associated with neurological disorders, and both differed from the classes enriched for height, a biometric trait used here as a control outgroup. The SNPs implicated in these psychiatric disorders and behavioral traits were highly enriched in CTCF binding sites, in conserved regions likely to be enhancers, and in brain-specific promoters, regulatory sites likely to affect responses to experience. The SNPs relevant for neurological disorders were highly enriched in constitutive coding regions and splice regulatory sites.

Sex affects transcriptional associations with schizophrenia across the dorsolateral prefrontal cortex, hippocampus, and caudate nucleus

Schizophrenia is a complex neuropsychiatric disorder with sexually dimorphic features, including differential symptomatology, drug responsiveness, and male incidence rate. Prior large-scale transcriptome analyses for sex differences in schizophrenia have focused on the prefrontal cortex. Analyzing BrainSeq Consortium data (caudate nucleus: n = 399, dorsolateral prefrontal cortex: n = 377, and hippocampus: n = 394), we identified 831 unique genes that exhibit sex differences across brain regions, enriched for immune-related pathways. We observed X-chromosome dosage reduction in the hippocampus of male individuals with schizophrenia. Our sex interaction model revealed 148 junctions dysregulated in a sex-specific manner in schizophrenia. Sex-specific schizophrenia analysis identified dozens of differentially expressed genes, notably enriched in immune-related pathways. Finally, our sex-interacting expression quantitative trait loci analysis revealed 704 unique genes, nine associated with schizophrenia risk. These findings emphasize the importance of sex-informed analysis of sexually dimorphic traits, inform personalized therapeutic strategies in schizophrenia, and highlight the need for increased female samples for schizophrenia analyses.

Cell-type and dynamic state govern genetic regulation of gene expression in heterogeneous differentiating cultures

Identifying the molecular effects of human genetic variation across cellular contexts is crucial for understanding the mechanisms underlying disease-associated loci, yet many cell-types and developmental stages remain underexplored. Here we harnessed the potential of heterogeneous differentiating cultures ( HDCs ), an in vitro system in which pluripotent cells asynchronously differentiate into a broad spectrum of cell-types. We generated HDCs for 53 human donors and collected single-cell RNA-sequencing data from over 900,000 cells. We identified expression quantitative trait loci in 29 cell-types and characterized regulatory dynamics across diverse differentiation trajectories. This revealed novel regulatory variants for genes involved in key developmental and disease-related processes while replicating known effects from primary tissues, and dynamic regulatory effects associated with a range of complex traits.

The cliff edge model of the evolution of schizophrenia: Mathematical, epidemiological, and genetic evidence

How has schizophrenia, a condition that significantly reduces an individual's evolutionary fitness, remained common across generations and cultures? Numerous theories about the evolution of schizophrenia have been proposed, most of which are not consistent with modern epidemiological and genetic evidence. Here, we briefly review this evidence and explore the cliff edge model of schizophrenia. It suggests that schizophrenia is the extreme manifestation of a polygenic trait or a combination of traits that, within a normal range of variation, confer cognitive, linguistic, and/or social advantages. Only beyond a certain threshold, these traits precipitate the onset of schizophrenia and reduce fitness. We provide the first mathematical model of this qualitative concept and show that it requires only very weak positive selection of the underlying trait(s) to explain today's schizophrenia prevalence. This prediction, along with expectations about the effect size of schizophrenia risk alleles, are surprisingly well matched by empirical evidence. The cliff edge model predicts a dynamic change of selection of risk alleles, which explains the contradictory findings of evolutionary genetic studies.

Phenotypic and Genetic Associations Between Preschool Fine Motor Skills and Later Neurodevelopment, Psychopathology, and Educational Achievement

BACKGROUND

Fine motor skills are heritable and comprise important milestones in development, and some evidence suggests that impairments in fine motor skills are associated with neurodevelopmental conditions, psychiatric disorders, and poor educational outcomes.

METHODS

In a preregistered study of 9625 preschool children from TEDS (Twins Early Development Study), fine motor assessments (drawing, block building, folding, and questionnaires) were conducted at 2, 3, and 4 years of age. A cross-age fine motor score was derived using principal component analysis. Multivariate regression analysis was used to examine the relationships between the fine motor score and neurodevelopmental traits, psychopathology, and educational outcomes at 3 later ages (7-8, 12, and 16 years) and cross-age psychopathology composite scores. Polygenic scores (PGSs) were created for attention-deficit/hyperactivity disorder (ADHD), autism, schizophrenia, anxiety, major depressive disorder, obsessive-compulsive disorder, and years of education. We ran single-PGS models and a multi-PGS model.

RESULTS

Fine motor skills were negatively associated with neurodevelopmental traits and psychopathology across childhood and adolescence and positively associated with educational achievement in adolescence (β = 0.25, p < .001). Superior fine motor skills were associated with a higher years-of-education PGS (β = 0.07, p < .001), a lower ADHD PGS (β = -0.04, p = .011), and a higher anxiety PGS (β = 0.03, p = .040). Similarly, the multi-PGS model retained the PGSs for years of education (β = 0.07), ADHD (β = -0.03), and anxiety (β = 0.01). A non-preregistered analysis in an independent preschool sample replicated the ADHD PGS association, but not the years of education or anxiety PGS associations.

CONCLUSIONS

Fine motor skills are linked genetically and phenotypically to later neurodevelopment, psychopathology, and educational outcomes. Future work should investigate the mechanisms that underlie the role of fine motor development in later outcomes.

Hippocampal Pyk2 regulates specific social skills: Implications for schizophrenia

Pyk2 has been shown previously to be involved in several psychological and cognitive alterations related to stress, Huntington's disease, and Alzheimer's disease. All these disorders are accompanied by different types of impairments in sociability, which has recently been linked to improper mitochondrial function. We hypothesize that Pyk2, which regulates mitochondria, could be associated with the regulation of mitochondrial dynamics and social skills. In the present manuscript, we report that a reduction of Pyk2 levels in mouse pyramidal neurons of the hippocampus decreased social dominance and aggressivity. Furthermore, social interactions induced robust Pyk2-dependent hippocampal changes in several oxidative phosphorylation complexes. We also observed that Pyk2 levels were increased in the CA1 pyramidal neurons of schizophrenic subjects, occurring alongside changes in different direct and indirect regulators of mitochondrial function including DISC1 and Grp75. Accordingly, overexpressing Pyk2 in hippocampal CA1 pyramidal cells mimicked some specific schizophrenia-like social behaviors in mice. In summary, our results indicate that Pyk2 might play a role in regulating specific social skills likely via mitochondrial dynamics and that there might be a link between Pyk2 levels in hippocampal neurons and social disturbances in schizophrenia.

Parvalbumin interneuron deficits in schizophrenia

Parvalbumin-expressing (PV+) interneurons represent one of the most abundant subclasses of cortical interneurons. Owing to their specific electrophysiological and synaptic properties, PV+ interneurons are essential for gating and pacing the activity of excitatory neurons. In particular, PV+ interneurons are critically involved in generating and maintaining cortical rhythms in the gamma frequency, which are essential for complex cognitive functions. Deficits in PV+ interneurons have been frequently reported in postmortem studies of schizophrenia patients, and alterations in gamma oscillations are a prominent electrophysiological feature of the disease. Here, I summarise the main features of PV+ interneurons and review clinical and preclinical studies linking the developmental dysfunction of cortical PV+ interneurons with the pathophysiology of schizophrenia.

Enrichment of the Local Synaptic Translatome for Genetic Risk Associated With Schizophrenia and Autism Spectrum Disorder

BACKGROUND

Genes that encode synaptic proteins or messenger RNA targets of the RNA-binding protein FMRP (fragile X messenger ribonucleoprotein) have been linked to schizophrenia and autism spectrum disorder (ASD) through the enrichment of genetic variants that confer risk for these disorders. FMRP binds many transcripts with synaptic functions and is thought to regulate their local translation, a process that enables rapid and compartmentalized protein synthesis required for development and plasticity.

METHODS

We used summary statistics from large-scale genome-wide association studies of schizophrenia (74,776 cases, 101,023 controls) and ASD (18,381 cases, 27,969 controls) to test the hypothesis that the subset of synaptic genes that encode localized transcripts is more strongly associated with each disorder than nonlocalized transcripts. We also postulated that this subset of synaptic genes is responsible for associations attributed to FMRP targets.

RESULTS

Schizophrenia associations were enriched in genes encoding localized synaptic transcripts compared to the remaining synaptic genes or to the remaining localized transcripts; this also applied to ASD associations, although only for transcripts observed after stimulation by fear conditioning. The genetic associations with either disorder captured by these gene sets were independent of those derived from FMRP targets. Schizophrenia association was related to FMRP interactions with messenger RNAs in somata, but not in dendrites, while ASD association was related to FMRP binding in either compartment.

CONCLUSIONS

Our data suggest that synaptic transcripts capable of local translation are particularly relevant to the pathogenesis of schizophrenia and ASD, but they do not characterize the associations attributed to current sets of FMRP targets.

xWAS analysis in neuropsychiatric disorders by integrating multi-molecular phenotype quantitative trait loci and GWAS summary data

BACKGROUND

Integrating quantitative trait loci (QTL) data related to molecular phenotypes with genome-wide association study (GWAS) data is an important post-GWAS strategic approach employed to identify disease-associated molecular features. Various types of molecular phenotypes have been investigated in neuropsychiatric disorders. However, these findings pertaining to distinct molecular features are often independent of each other, posing challenges for having an overview of the mapped genes.

METHODS

In this study, we comprehensively summarized published analyses focusing on four types of risk-related molecular features (gene expression, splicing transcriptome, protein abundance, and DNA methylation) across five common neuropsychiatric disorders. Subsequently, we conducted supplementary analyses with the latest GWAS dataset and corresponding deficient molecular phenotypes using Functional Summary-based Imputation (FUSION) and summary data-based Mendelian randomization (SMR). Based on the curated and supplemented results, novel reliable genes and their functions were explored.

RESULTS

Our findings revealed that eQTL exhibited superior ability in prioritizing risk genes compared to the other QTL, followed by sQTL. Approximately half of the genes associated with splicing transcriptome, protein abundance, and DNA methylation were successfully replicated by eQTL-associated genes across all five disorders. Furthermore, we identified 436 novel reliable genes, which enriched in pathways related with neurotransmitter transportation such as synaptic, dendrite, vesicles, axon along with correlations with other neuropsychiatric disorders. Finally, we identified ten multiple molecular involved regulation patterns (MMRP), which may provide valuable insights into understanding the contribution of molecular regulation network targeting these disease-associated genes.

CONCLUSIONS

The analyses prioritized novel and reliable gene sets related with five molecular features based on published and supplementary results for five common neuropsychiatric disorders, which were missed in the original GWAS analysis. Besides, the involved MMRP behind these genes could be given priority for further investigation to elucidate the pathogenic molecular mechanisms underlying neuropsychiatric disorders in future studies.

Massively parallel reporter assays and mouse transgenic assays provide complementary information about neuronal enhancer activity

Genetic studies find hundreds of thousands of noncoding variants associated with psychiatric disorders. Massively parallel reporter assays (MPRAs) and in vivo transgenic mouse assays can be used to assay the impact of these variants. However, the relevance of MPRAs to in vivo function is unknown and transgenic assays suffer from low throughput. Here, we studied the utility of combining the two assays to study the impact of non-coding variants. We carried out an MPRA on over 50,000 sequences derived from enhancers validated in transgenic mouse assays and from multiple fetal neuronal ATAC-seq datasets. We also tested over 20,000 variants, including synthetic mutations in highly active neuronal enhancers and 177 common variants associated with psychiatric disorders. Variants with a high impact on MPRA activity were further tested in mice. We found a strong and specific correlation between MPRA and mouse neuronal enhancer activity including changes in neuronal enhancer activity in mouse embryos for variants with strong MPRA effects. Mouse assays also revealed pleiotropic variant effects that could not be observed in MPRA. Our work provides a large catalog of functional neuronal enhancers and variant effects and highlights the effectiveness of combining MPRAs and mouse transgenic assays.

Dissecting Schizophrenia Biology Using Pleiotropy with Cognitive Genomics

Given the increasingly large number of loci discovered by psychiatric GWAS, specification of the key biological pathways underlying these loci has become a priority for the field. We have previously leveraged the pleiotropic genetic relationships between schizophrenia and two cognitive phenotypes (educational attainment and cognitive task performance) to differentiate two subsets of illness-relevant SNPs: (1) those with "concordant" alleles, which are associated with reduced cognitive ability/education and increased schizophrenia risk; and (2) those with "discordant" alleles linked to reduced educational and/or cognitive levels but lower schizophrenia susceptibility. In the present study, we extend our prior work, utilizing larger input GWAS datasets and a more powerful statistical approach to pleiotropic meta-analysis, the Pleiotropic Locus Exploration and Interpretation using Optimal test (PLEIO). Our pleiotropic meta-analysis of schizophrenia and the two cognitive phenotypes revealed 768 significant loci (159 novel). Among these, 347 loci harbored concordant SNPs, 270 encompassed discordant SNPs, and 151 "dual" loci contained concordant and discordant SNPs. Competitive gene-set analysis using MAGMA related concordant SNP loci with neurodevelopmental pathways (e.g., neurogenesis), whereas discordant loci were associated with mature neuronal synaptic functions. These distinctions were also observed in BrainSpan analysis of temporal enrichment patterns across developmental periods, with concordant loci containing more prenatally expressed genes than discordant loci. Dual loci were enriched for genes related to mRNA translation initiation, representing a novel finding in the schizophrenia literature.

A preliminary metabolomics study of the database for biological samples of schizophrenia among Chinese ethnic minorities

BACKGROUND

Schizophrenia (SCZ) is a profound mental disorder with a multifactorial etiology, including genetics, environmental factors, and demographic influences such as ethnicity and geography. Among these, the studies of SCZ also shows racial and regional differences.

METHODS

We first established a database of biological samples for SCZ in China's ethnic minorities, followed by a serum metabolomic analysis of SCZ patients from various ethnic groups within the same region using the LC-HRMS platform.

RESULTS

Analysis identified 47 metabolites associated with SCZ, with 46 showing significant differences between Miao and Han SCZ patients. These metabolites, primarily fatty acids, amino acids, benzene, and derivatives, are involved in fatty acid metabolism pathways. Notably, L-Carnitine, L-Cystine, Aspartylphenylalanine, and Methionine sulfoxide demonstrated greater diagnostic efficacy in Miao SCZ patients compared to Han SCZ patients.

CONCLUSION

Preliminary findings suggest that there are differences in metabolic levels among SCZ patients of different ethnicities in the same region, offering insights for developing objective diagnostic or therapeutic monitoring strategies that incorporate ethnic considerations of SCZ.

Discrimination between healthy participants and people with panic disorder based on polygenic scores for psychiatric disorders and for intermediate phenotypes using machine learning

OBJECTIVE

Panic disorder is a modestly heritable condition. Currently, diagnosis is based only on clinical symptoms; identifying objective biomarkers and a more reliable diagnostic procedure is desirable. We investigated whether people with panic disorder can be reliably diagnosed utilizing combinations of multiple polygenic scores for psychiatric disorders and their intermediate phenotypes, compared with single polygenic score approaches, by applying specific machine learning techniques.

METHODS

Polygenic scores for 48 psychiatric disorders and intermediate phenotypes based on large-scale genome-wide association studies (n = 7556-1,131,881) were calculated for people with panic disorder (n = 718) and healthy controls (n = 1717). Discrimination between people with panic disorder and healthy controls was based on the 48 polygenic scores using five methods for classification: logistic regression, neural networks, quadratic discriminant analysis, random forests and a support vector machine. Differences in discrimination accuracy (area under the curve) due to an increased number of polygenic score combinations and differences in the accuracy across five classifiers were investigated.

RESULTS

All five classifiers performed relatively well for distinguishing people with panic disorder from healthy controls by increasing the number of polygenic scores. Of the 48 polygenic scores, the polygenic score for anxiety UK Biobank was the most useful for discrimination by the classifiers. In combinations of two or three polygenic scores, the polygenic score for anxiety UK Biobank was included as one of polygenic scores in all classifiers. When all 48 polygenic scores were used in combination, the greatest areas under the curve significantly differed among the five classifiers. Support vector machine and logistic regression had higher accuracy than quadratic discriminant analysis and random forests. For each classifier, the greatest area under the curve was 0.600 ± 0.030 for logistic regression (polygenic score combinations N = 14), 0.591 ± 0.039 for neural networks (N = 9), 0.603 ± 0.033 for quadratic discriminant analysis (N = 10), 0.572 ± 0.039 for random forests (N = 25) and 0.617 ± 0.041 for support vector machine (N = 11). The greatest areas under the curve at the best polygenic score combination significantly differed among the five classifiers. Random forests had the lowest accuracy among classifiers. Support vector machine had higher accuracy than neural networks.

CONCLUSIONS

These findings suggest that increasing the number of polygenic score combinations up to approximately 10 effectively improved the discrimination accuracy and that support vector machine exhibited greater accuracy among classifiers. However, the discrimination accuracy for panic disorder, when based solely on polygenic score combinations, was found to be modest.

Patterns of stressful life events and polygenic scores for five mental disorders and neuroticism among adults with depression

The dominant ('general') version of the diathesis-stress theory of depression views stressors and genetic vulnerability as independent risks. In the Australian Genetics of Depression Study (N = 14,146; 75% female), we tested whether polygenic scores (PGS) for major depression, bipolar disorder, schizophrenia, anxiety, ADHD, and neuroticism were associated with reported exposure to 32 childhood, past-year, lifetime, and accumulated stressful life events (SLEs). In false discovery rate-corrected models, the clearest PGS-SLE relationships were for the ADHD- and depression-PGSs, and to a lesser extent, the anxiety- and schizophrenia-PGSs. We describe the associations for childhood and accumulated SLEs, and the 2-3 strongest past-year/lifetime SLE associations. Higher ADHD-PGS was associated with all childhood SLEs (emotional abuse, emotional neglect, physical neglect; ORs = 1.09-1.14; p's < 1.3 × 10-5), more accumulated SLEs, and reported exposure to sudden violent death (OR = 1.23; p = 3.6 × 10-5), legal troubles (OR = 1.15; p = 0.003), and sudden accidental death (OR = 1.14; p = 0.006). Higher depression-PGS was associated with all childhood SLEs (ORs = 1.07-1.12; p's < 0.013), more accumulated SLEs, and severe human suffering (OR = 1.17; p = 0.003), assault with a weapon (OR = 1.12; p = 0.003), and living in unpleasant surroundings (OR = 1.11; p = 0.001). Higher anxiety-PGS was associated with childhood emotional abuse (OR = 1.08; p = 1.6 × 10-4), more accumulated SLEs, and serious accident (OR = 1.23; p = 0.004), physical assault (OR = 1.08; p = 2.2 × 10-4), and transportation accident (OR = 1.07; p = 0.001). Higher schizophrenia-PGS was associated with all childhood SLEs (ORs = 1.12-1.19; p's < 9.3-8), more accumulated SLEs, and severe human suffering (OR = 1.16; p = 0.003). Higher neuroticism-PGS was associated with living in unpleasant surroundings (OR = 1.09; p = 0.007) and major financial troubles (OR = 1.06; p = 0.014). A reversed pattern was seen for the bipolar-PGS, with lower odds of reported physical assault (OR = 0.95; p = 0.014), major financial troubles (OR = 0.93; p = 0.004), and living in unpleasant surroundings (OR = 0.92; p = 0.007). Genetic risk for several mental disorders influences reported exposure to SLEs among adults with moderately severe, recurrent depression. Our findings emphasise that stressors and diatheses are inter-dependent and challenge diagnosis and subtyping (e.g., reactive/endogenous) based on life events.

Morphogenetic theory of mental and cognitive disorders: the role of neurotrophic and guidance molecules

Mental illness and cognitive disorders represent a serious problem for the modern society. Many studies indicate that mental disorders are polygenic and that impaired brain development may lay the ground for their manifestation. Neural tissue development is a complex and multistage process that involves a large number of distant and contact molecules. In this review, we have considered the key steps of brain morphogenesis, and the major molecule families involved in these process. The review provides many indications of the important contribution of the brain development process and correct functioning of certain genes to human mental health. To our knowledge, this comprehensive review is one of the first in this field. We suppose that this review may be useful to novice researchers and clinicians wishing to navigate the field.

Genetic contribution to disease-course severity and progression in the SUPER-Finland study, a cohort of 10,403 individuals with psychotic disorders

Genetic factors contribute to the susceptibility of psychotic disorders, but less is known how they affect psychotic disease-course development. Utilizing polygenic scores (PGSs) in combination with longitudinal healthcare data with decades of follow-up we investigated the contributing genetics to psychotic disease-course severity and diagnostic shifts in the SUPER-Finland study, encompassing 10 403 genotyped individuals with a psychotic disorder. To longitudinally track the study participants' past disease-course severity, we created a psychiatric hospitalization burden metric using the full-coverage and nation-wide Finnish in-hospital registry (data from 1969 and onwards). Using a hierarchical model, ranking the psychotic diagnoses according to clinical severity, we show that high schizophrenia PGS (SZ-PGS) was associated with progression from lower ranked psychotic disorders to schizophrenia (OR = 1.32 [1.23-1.43], p = 1.26e-12). This development manifested already at psychotic illness onset as a higher psychiatric hospitalization burden, the proxy for disease-course severity. In schizophrenia (n = 5 479), both a high SZ-PGS and a low educational attainment PGS (EA-PGS) were associated with increased psychiatric hospitalization burden (p = 1.00e-04 and p = 4.53e-10). The SZ-PGS and the EA-PGS associated with distinct patterns of hospital usage. In individuals with high SZ-PGS, the increased hospitalization burden was composed of longer individual hospital stays, while low EA-PGS associated with shorter but more frequent hospital visits. The negative effect of a low EA-PGS was found to be partly mediated via substance use disorder, a major risk factor for hospitalizations. In conclusion, we show that high SZ-PGS and low EA-PGS both impacted psychotic disease-course development negatively but resulted in different disease-course trajectories.

Copy number variant risk loci for schizophrenia converge on the BDNF pathway

OBJECTIVES

Schizophrenia genetics is intricate, with common and rare variants' contributions not fully understood. Certain copy number variations (CNVs) elevate risk, pivotal for understanding mental disorder models. Despite CNVs' genome-wide distribution and variable gene and protein effects, we must explore beyond affected genes to interaction partners and molecular pathways.

METHODS

In this study, we developed machine-readable interactive pathways to enable analysis of functional effects of genes within CNV loci and identify ten common pathways across CNVs with high schizophrenia risk using the WikiPathways database, schizophrenia risk gene collections from GWAS studies, and a gene-disease association database.

RESULTS

For CNVs that are pathogenic for schizophrenia, we found overlapping pathways, including BDNF signalling, cytoskeleton, and inflammation. Common schizophrenia risk genes identified by different studies are found in all CNV pathways, but not enriched.

CONCLUSIONS

Our findings suggest that specific pathways - BDNF signalling - are critical contributors to schizophrenia risk conferred by rare CNVs. Our approach highlights the importance of not only investigating deleted or duplicated genes within pathogenic CNV loci, but also study their direct interaction partners, which may explain pleiotropic effects of CNVs on schizophrenia risk and offer a broader field for interventions.

Orphan GPR52 as an emerging neurotherapeutic target

GPR52 is a highly conserved, brain-enriched, Gs/olf-coupled orphan G protein-coupled receptor (GPCR) that controls various cyclic AMP (cAMP)-dependent physiological and pathological processes. Stimulation of GPR52 activity might be beneficial for the treatment of schizophrenia, psychiatric disorders and other human neurological diseases, whereas inhibition of its activity might provide a potential therapeutic approach for Huntington's disease. Excitingly, HTL0048149 (HTL'149), an orally available GPR52 agonist, has been advanced into phase I human clinical trials for the treatment of schizophrenia. In this concise review, we summarize the current understanding of GPR52 receptor distribution as well as its structure and functions, highlighting the recent advances in drug discovery efforts towards small-molecule GPR52 ligands. The opportunities and challenges presented by targeting GPR52 for novel therapeutics are also briefly discussed.

Genetic and Phenotypic Features of Schizophrenia in the UK Biobank

Importance

Large-scale biobanks provide important opportunities for mental health research, but selection biases raise questions regarding the comparability of individuals with those in clinical research settings.

Objective

To compare the genetic liability to psychiatric disorders in individuals with schizophrenia in the UK Biobank with individuals in the Psychiatric Genomics Consortium (PGC) and to compare genetic liability and phenotypic features with participants recruited from clinical settings.

Design, Setting, and Participants

This cross-sectional study included participants from the population-based UK Biobank and schizophrenia samples recruited from clinical settings (CLOZUK, CardiffCOGS, Cardiff F-Series, and Cardiff Affected Sib-Pairs). Data were collected between January 1993 and July 2021. Data analysis was conducted between July 2021 and June 2023.

Main Outcomes and Measures

A genome-wide association study of UK Biobank schizophrenia case-control status was conducted, and the results were compared with those from the PGC via genetic correlations. To test for differences with the clinical samples, polygenic risk scores (PRS) were calculated for schizophrenia, bipolar disorder, depression, and intelligence using PRS-CS. PRS and phenotypic comparisons were conducted using pairwise logistic regressions. The proportions of individuals with copy number variants associated with schizophrenia were compared using Firth logistic regression.

Results

The sample of 517 375 participants included 1438 UK Biobank participants with schizophrenia (550 [38.2%] female; mean [SD] age, 54.7 [8.3] years), 499 475 UK Biobank controls (271 884 [54.4%] female; mean [SD] age, 56.5 [8.1] years), and 4 schizophrenia research samples (4758 [28.9%] female; mean [SD] age, 38.2 [21.0] years). Liability to schizophrenia in UK Biobank was highly correlated with the latest genome-wide association study from the PGC (genetic correlation, 0.98; SE, 0.18) and showed the expected patterns of correlations with other psychiatric disorders. The schizophrenia PRS explained 6.8% of the variance in liability for schizophrenia case status in UK Biobank. UK Biobank participants with schizophrenia had significantly lower schizophrenia PRS than 3 of the clinically ascertained samples and significantly lower rates of schizophrenia-associated copy number variants than the CLOZUK sample. UK Biobank participants with schizophrenia had higher educational attainment and employment rates than the clinically ascertained schizophrenia samples, lower rates of smoking, and a later age of onset of psychosis.

Conclusions and Relevance

Individuals with schizophrenia in the UK Biobank, and likely other volunteer-based biobanks, represent those less severely affected. Their inclusion in wider studies should enhance the representation of the full spectrum of illness severity.

Accelerated Cortical Thinning in Schizophrenia Is Associated With Rare and Common Predisposing Variation to Schizophrenia and Neurodevelopmental Disorders

BACKGROUND

Schizophrenia is a highly heritable disorder characterized by increased cortical thinning throughout the life span. Studies have reported a shared genetic basis between schizophrenia and cortical thickness. However, no genes whose expression is related to abnormal cortical thinning in schizophrenia have been identified.

METHODS

We conducted linear mixed models to estimate the rates of accelerated cortical thinning across 68 regions from the Desikan-Killiany atlas in individuals with schizophrenia compared with healthy control participants from a large longitudinal sample (ncases = 169 and ncontrols = 298, ages 16-70 years). We studied the correlation between gene expression data from the Allen Human Brain Atlas and accelerated thinning estimates across cortical regions. Finally, we explored the functional and genetic underpinnings of the genes that contribute most to accelerated thinning.

RESULTS

We found a global pattern of accelerated cortical thinning in individuals with schizophrenia compared with healthy control participants. Genes underexpressed in cortical regions that exhibit this accelerated thinning were downregulated in several psychiatric disorders and were enriched for both common and rare disrupting variation for schizophrenia and neurodevelopmental disorders. In contrast, none of these enrichments were observed for baseline cross-sectional cortical thickness differences.

CONCLUSIONS

Our findings suggest that accelerated cortical thinning, rather than cortical thickness alone, serves as an informative phenotype for neurodevelopmental disruptions in schizophrenia. We highlight the genetic and transcriptomic correlates of this accelerated cortical thinning, emphasizing the need for future longitudinal studies to elucidate the role of genetic variation and the temporal-spatial dynamics of gene expression in brain development and aging in schizophrenia.

Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder

Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed health problems, prior pharmacological treatments, and polygenic scores (PGS) has potential to inform risk stratification. We examined self-reported SB and ideation using the Columbia Suicide Severity Rating Scale (C-SSRS) among 3,942 SCZ and 5,414 BPI patients receiving care within the Veterans Health Administration (VHA). These cross-sectional data were integrated with electronic health records (EHRs), and compared across lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. PGS were constructed using available genomic data for related traits. Genome-wide association studies were performed to identify and prioritize specific loci. Only 20% of the veterans who reported SB had a corroborating ICD-9/10 EHR code. Among those without prior SB, more than 20% reported new-onset SB at follow-up. SB were associated with a range of additional clinical diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking initiation, suicide attempt, and major depressive disorder were associated with SB. The GWAS for SB yielded no significant loci. Among individuals with a diagnosed mental illness, self-reported SB were strongly associated with clinical variables across several EHR domains. Analyses point to sequelae of substance-related and psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in health records, underscoring the value of regular screening with direct, in-person assessments, especially among high-risk individuals.

The Risk Genes for Neuropsychiatric Disorders negr1 and opcml Are Expressed throughout Zebrafish Brain Development

The immunoglobulin LAMP/OBCAM/NTM (IgLON) family of cell adhesion molecules comprises five members known for their involvement in establishing neural circuit connectivity, fine-tuning, and maintenance. Mutations in IgLON genes result in alterations in these processes and can lead to neuropsychiatric disorders. The two IgLON family members NEGR1 and OPCML share common links with several of them, such as schizophrenia, autism, and major depressive disorder. However, the onset and the underlying molecular mechanisms have remained largely unresolved, hampering progress in developing therapies. NEGR1 and OPCML are evolutionarily conserved in teleosts like the zebrafish (Danio rerio), which is excellently suited for disease modelling and large-scale screening for disease-ameliorating compounds. To explore the potential applicability of zebrafish for extending our knowledge on NEGR1- and OPCML-linked disorders and to develop new therapeutic strategies, we investigated the spatio-temporal expression of the two genes during early stages of development. negr1 and opcml are expressed maternally and subsequently in partially distinct domains of conserved brain regions. Other areas of expression in zebrafish have not been reported in mammals to date. Our results indicate that NEGR1 and OPCML may play roles in neural circuit development and function at stages earlier than previously anticipated. A detailed functional analysis of the two genes based on our findings could contribute to understanding the mechanistic basis of related psychiatric disorders.

Genome-wide association study identifies 30 obsessive-compulsive disorder associated loci

Obsessive-compulsive disorder (OCD) affects ~1% of the population and exhibits a high SNP-heritability, yet previous genome-wide association studies (GWAS) have provided limited information on the genetic etiology and underlying biological mechanisms of the disorder. We conducted a GWAS meta-analysis combining 53,660 OCD cases and 2,044,417 controls from 28 European-ancestry cohorts revealing 30 independent genome-wide significant SNPs and a SNP-based heritability of 6.7%. Separate GWAS for clinical, biobank, comorbid, and self-report sub-groups found no evidence of sample ascertainment impacting our results. Functional and positional QTL gene-based approaches identified 249 significant candidate risk genes for OCD, of which 25 were identified as putatively causal, highlighting WDR6, DALRD3, CTNND1 and genes in the MHC region. Tissue and single-cell enrichment analyses highlighted hippocampal and cortical excitatory neurons, along with D1- and D2-type dopamine receptor-containing medium spiny neurons, as playing a role in OCD risk. OCD displayed significant genetic correlations with 65 out of 112 examined phenotypes. Notably, it showed positive genetic correlations with all included psychiatric phenotypes, in particular anxiety, depression, anorexia nervosa, and Tourette syndrome, and negative correlations with a subset of the included autoimmune disorders, educational attainment, and body mass index.. This study marks a significant step toward unraveling its genetic landscape and advances understanding of OCD genetics, providing a foundation for future interventions to address this debilitating disorder.

Genome-wide association study identifies new loci associated with OCD

To date, four genome-wide association studies (GWAS) of obsessive-compulsive disorder (OCD) have been published, reporting a high single-nucleotide polymorphism (SNP)-heritability of 28% but finding only one significant SNP. A substantial increase in sample size will likely lead to further identification of SNPs, genes, and biological pathways mediating the susceptibility to OCD. We conducted a GWAS meta-analysis with a 2-3-fold increase in case sample size (OCD cases: N = 37,015, controls: N = 948,616) compared to the last OCD GWAS, including six previously published cohorts (OCGAS, IOCDF-GC, IOCDF-GC-trio, NORDiC-nor, NORDiC-swe, and iPSYCH) and unpublished self-report data from 23andMe Inc. We explored the genetic architecture of OCD by conducting gene-based tests, tissue and celltype enrichment analyses, and estimating heritability and genetic correlations with 74 phenotypes. To examine a potential heterogeneity in our data, we conducted multivariable GWASs with MTAG. We found support for 15 independent genome-wide significant loci (14 new) and 79 protein-coding genes. Tissue enrichment analyses implicate multiple cortical regions, the amygdala, and hypothalamus, while cell type analyses yielded 12 cell types linked to OCD (all neurons). The SNP-based heritability of OCD was estimated to be 0.08. Using MTAG we found evidence for specific genetic underpinnings characteristic of different cohort-ascertainment and identified additional significant SNPs. OCD was genetically correlated with 40 disorders or traits-positively with all psychiatric disorders and negatively with BMI, age at first birth and multiple autoimmune diseases. The GWAS meta-analysis identified several biologically informative genes as important contributors to the aetiology of OCD. Overall, we have begun laying the groundwork through which the biology of OCD will be understood and described.

Temporal changes in brain morphology related to inflammation and schizophrenia: an omnigenic Mendelian randomization study

BACKGROUND

Over the past several decades, more research focuses have been made on the inflammation/immune hypothesis of schizophrenia. Building upon synaptic plasticity hypothesis, inflammation may contribute the underlying pathophysiology of schizophrenia. Yet, pinpointing the specific inflammatory agents responsible for schizophrenia remains a complex challenge, mainly due to medication and metabolic status. Multiple lines of evidence point to a wide-spread genetic association across genome underlying the phenotypic variations of schizophrenia.

METHOD

We collected the latest genome-wide association analysis (GWAS) summary data of schizophrenia, cytokines, and longitudinal change of brain. We utilized the omnigenic model which takes into account all genomic SNPs included in the GWAS of trait, instead of traditional Mendelian randomization (MR) methods. We conducted two round MR to investigate the inflammatory triggers of schizophrenia and the resulting longitudinal changes in the brain.

RESULTS

We identified seven inflammation markers linked to schizophrenia onset, which all passed the Bonferroni correction for multiple comparisons (bNGF, GROA(CXCL1), IL-8, M-CSF, MCP-3 (CCL7), TNF-β, CRP). Moreover, CRP were found to significantly influence the linear rate of brain morphology changes, predominantly in the white matter of the cerebrum and cerebellum.

CONCLUSION

With an omnigenic approach, our study sheds light on the immune pathology of schizophrenia. Although these findings need confirmation from future studies employing different methodologies, our work provides substantial evidence that pervasive, low-level neuroinflammation may play a pivotal role in schizophrenia, potentially leading to notable longitudinal changes in brain morphology.

Predicting cell population-specific gene expression from genomic sequence

Most regulatory elements, especially enhancer sequences, are cell population-specific. One could even argue that a distinct set of regulatory elements is what defines a cell population. However, discovering which non-coding regions of the DNA are essential in which context, and as a result, which genes are expressed, is a difficult task. Some computational models tackle this problem by predicting gene expression directly from the genomic sequence. These models are currently limited to predicting bulk measurements and mainly make tissue-specific predictions. Here, we present a model that leverages single-cell RNA-sequencing data to predict gene expression. We show that cell population-specific models outperform tissue-specific models, especially when the expression profile of a cell population and the corresponding tissue are dissimilar. Further, we show that our model can prioritize GWAS variants and learn motifs of transcription factor binding sites. We envision that our model can be useful for delineating cell population-specific regulatory elements.

Examination of a novel expression-based gene-SNP annotation strategy to identify tissue-specific contributions to heritability in multiple traits

Complex traits show clear patterns of tissue-specific expression influenced by single nucleotide polymorphisms (SNPs), yet current strategies aggregate SNP effects to genes by employing simple physical proximity-based windows. Here, we examined whether incorporating SNPs with effects on tissue-specific cis-expression would improve our ability to detect trait-relevant tissues across 31 complex traits using stratified linkage disequilibrium score regression (S-LDSC). We found that a physical proximity annotation produced more significant tissue enrichments and larger S-LDSC regression coefficients, as compared to an expression-based annotation. Furthermore, we showed that our expression-based annotation did not outperform an annotation strategy in which an equal number of randomly chosen SNPs were annotated to genes within the same genomic window, suggesting extensive redundancy among SNP effect estimates due to linkage disequilibrium. That said, current sample sizes limit estimation of cis-genetic SNP effects; therefore, we recommend reexamination of the expression-based annotation when larger tissue-specific expression datasets become available. To examine the influence of sample size, we used a large whole blood eQTL reference panel (N = 31,684) applying a similar expression-based annotation strategy. We found that significant cis-expression QTLs in whole blood did not outperform the physical proximity annotation when estimating tissue-specific SNP heritability enrichment for either high- or low-density lipoprotein phenotypes but performed similarly for inflammatory bowel disease. Finally, we report new and updated tissue enrichment estimates across 31 complex traits, such as significant heritability enrichment of the frontal cortex for cognitive performance, educational attainment, and intelligence, providing further evidence of this structure's importance in higher cognitive function.

Rare variants in pharmacogenes influence clozapine metabolism in individuals with schizophrenia

Clozapine is the only licensed medication for treatment-resistant schizophrenia (TRS). Few predictors for variation in response to clozapine have been identified, but clozapine metabolism is known to influence therapeutic response and adverse side effects. Here, we expand on genome-wide studies of clozapine metabolism, previously focused on common genetic variation, by analysing whole-exome sequencing data from 2062 individuals with schizophrenia taking clozapine in the UK. We investigated whether rare genomic variation in genes and gene sets involved in the clozapine metabolism pathway influences plasma concentrations of clozapine metabolites, assessed through the longitudinal analysis of 6585 pharmacokinetic assays. We observed a statistically significant association between the burden of rare damaging coding variants (MAF ≤ 1 %) in gene sets broadly related to drug pharmacokinetics and lower clozapine (β = -0.054, SE = 0.019, P-value = 0.005) concentrations in plasma. We estimate that the effects in clozapine plasma concentrations of a single damaging allele in this gene set are akin to reducing the clozapine dose by about 35 mg/day. The gene-based analysis identified rare variants in CYP1A2, which encodes the enzyme responsible for converting clozapine to norclozapine, as having the strongest effects of any gene on clozapine metabolism (β = 0.324, SE = 0.124, P = 0.009). Our findings support the hypothesis that rare genetic variants in known drug-metabolising enzymes and transporters can markedly influence clozapine plasma concentrations; these results suggest that pharmacogenomic efforts trying to predict clozapine metabolism and personalise drug therapy could benefit from the inclusion of rare damaging variants in pharmacogenes beyond those already identified and catalogued as PGx star alleles.

A polygenic resilience score moderates the genetic risk for schizophrenia: Replication in 18,090 cases and 28,114 controls from the Psychiatric Genomics Consortium

Identifying heritable factors that moderate the genetic risk for schizophrenia (SCZ) could help clarify why some individuals remain unaffected despite having relatively high genetic liability. Previously, we developed a framework to mine genome-wide association (GWAS) data for common genetic variants that protect high-risk unaffected individuals from SCZ, leading to derivation of the first-ever "polygenic resilience score" for SCZ (resilient controls n = 3786; polygenic risk score-matched SCZ cases n = 18,619). Here, we performed a replication study to verify the moderating effect of our polygenic resilience score on SCZ risk (OR = 1.09, p = 4.03 × 10-5 ) using newly released GWAS data from 23 independent case-control studies collated by the Psychiatric Genomics Consortium (PGC) (resilient controls n = 2821; polygenic risk score-matched SCZ cases n = 5150). Additionally, we sought to optimize our polygenic resilience-scoring formula to improve subsequent modeling of resilience to SCZ and other complex disorders. We found significant replication of the polygenic resilience score, and found that strict pruning of SNPs based on linkage disequilibrium to known risk SNPs and their linked loci optimizes the performance of the polygenic resilience score.

Parkinson's disease and schizophrenia interactomes contain temporally distinct gene clusters underlying comorbid mechanisms and unique disease processes

Genome-wide association studies suggest significant overlaps in Parkinson's disease (PD) and schizophrenia (SZ) risks, but the underlying mechanisms remain elusive. The protein-protein interaction network ('interactome') plays a crucial role in PD and SZ and can incorporate their spatiotemporal specificities. Therefore, to study the linked biology of PD and SZ, we compiled PD- and SZ-associated genes from the DisGeNET database, and constructed their interactomes using BioGRID and HPRD. We examined the interactomes using clustering and enrichment analyses, in conjunction with the transcriptomic data of 26 brain regions spanning foetal stages to adulthood available in the BrainSpan Atlas. PD and SZ interactomes formed four gene clusters with distinct temporal identities (Disease Gene Networks or 'DGNs'1-4). DGN1 had unique SZ interactome genes highly expressed across developmental stages, corresponding to a neurodevelopmental SZ subtype. DGN2, containing unique SZ interactome genes expressed from early infancy to adulthood, correlated with an inflammation-driven SZ subtype and adult SZ risk. DGN3 contained unique PD interactome genes expressed in late infancy, early and late childhood, and adulthood, and involved in mitochondrial pathways. DGN4, containing prenatally-expressed genes common to both the interactomes, involved in stem cell pluripotency and overlapping with the interactome of 22q11 deletion syndrome (comorbid psychosis and Parkinsonism), potentially regulates neurodevelopmental mechanisms in PD-SZ comorbidity. Our findings suggest that disrupted neurodevelopment (regulated by DGN4) could expose risk windows in PD and SZ, later elevating disease risk through inflammation (DGN2). Alternatively, variant clustering in DGNs may produce disease subtypes, e.g., PD-SZ comorbidity with DGN4, and early/late-onset SZ with DGN1/DGN2.

The structure of psychiatric comorbidity without selection and assortative mating

The widespread comorbidity observed across psychiatric disorders may be the result of processes such as assortative mating, gene-environment correlation, or selection into population studies. Between-family analyses of comorbidity are subject to these sources of bias, whereas within-family analyses are not. Because of Mendelian inheritance, alleles are randomly assigned within families, conditional on parental alleles. We exploit this variation to compare the structure of comorbidity across broad psychiatric polygenic scores when calculated either between-family (child polygenic scores) or within-family (child polygenic scores regressed on parental polygenic scores) in over 25,000 genotyped parent-offspring trios from the Norwegian Mother Father and Child Cohort study (MoBa). We fitted a series of factor models to the between- and within-family data, which consisted of a single genetic p-factor and a varying number of uncorrelated subfactors. The best-fitting model was identical for between- and within-family analyses and included three subfactors capturing variants associated with neurodevelopment, psychosis, and constraint, in addition to the genetic p-factor. Partner genetic correlations, indicating assortative mating, were not present for the genetic p-factor, but were substantial for the psychosis (b = 0.081;95% CI [0.038,0.124]) and constraint (b = 0.257;95% CI [0.075,0.439]) subfactors. When average factor levels for MoBa mothers and fathers were compared to a population mean of zero we found evidence of sex-specific participation bias, which has implications for the generalizability of findings from cohort studies. Our results demonstrate the power of the within-family design for better understanding the mechanisms driving psychiatric comorbidity and their consequences on population health.

Coexpression network analysis of the adult brain sheds light on the pathogenic mechanism of DDR1 in schizophrenia and bipolar disorder

DDR1 has been linked to schizophrenia (SCZ) and bipolar disorder (BD) in association studies. DDR1 encodes 58 distinct transcripts, which can be translated into five isoforms (DDR1a-e) and are expressed in the brain. However, the transcripts expressed in each brain cell type, their functions and their involvement in SCZ and BD remain unknown. Here, to infer the processes in which DDR1 transcripts are involved, we used transcriptomic data from the human brain dorsolateral prefrontal cortex of healthy controls (N = 936) and performed weighted gene coexpression network analysis followed by enrichment analyses. Then, to explore the involvement of DDR1 transcripts in SCZ (N = 563) and BD (N = 222), we studied the association of coexpression modules with disease and performed differential expression and transcript significance analyses. Some DDR1 transcripts were distributed across five coexpression modules identified in healthy controls (MHC). MHC1 and MHC2 were enriched in the cell cycle and proliferation of astrocytes and OPCs; MHC3 and MHC4 were enriched in oligodendrocyte differentiation and myelination; and MHC5 was enriched in neurons and synaptic transmission. Most of the DDR1 transcripts associated with SCZ and BD pertained to MHC1 and MHC2. Altogether, our results suggest that DDR1 expression might be altered in SCZ and BD via the proliferation of astrocytes and OPCs, suggesting that these processes are relevant in psychiatric disorders.

Processing speed mediates the relationship between DDR1 and psychosocial functioning in euthymic patients with bipolar disorder presenting psychotic symptoms

The DDR1 locus is associated with the diagnosis of schizophrenia and with processing speed in patients with schizophrenia and first-episode psychosis. Here, we investigated whether DDR1 variants are associated with bipolar disorder (BD) features. First, we performed a case‒control association study comparing DDR1 variants between patients with BD and healthy controls. Second, we performed linear regression analyses to assess the associations of DDR1 variants with neurocognitive domains and psychosocial functioning. Third, we conducted a mediation analysis to explore whether neurocognitive impairment mediated the association between DDR1 variants and psychosocial functioning in patients with BD. Finally, we studied the association between DDR1 variants and white matter microstructure. We did not find any statistically significant associations in the case‒control association study; however, we found that the combined genotypes rs1264323AA-rs2267641AC/CC were associated with worse neurocognitive performance in patients with BD with psychotic symptoms. In addition, the combined genotypes rs1264323AA-rs2267641AC/CC were associated with worse psychosocial functioning through processing speed. We did not find correlations between white matter microstructure abnormalities and the neurocognitive domains associated with the combined genotypes rs1264323AA-rs2267641AC/CC. Overall, the results suggest that DDR1 may be a marker of worse neurocognitive performance and psychosocial functioning in patients with BD, specifically those with psychotic symptoms.