Genome-wide association analysis identifies 13 new risk loci for schizophrenia

Novel genetic loci of inhibitory control in ADHD and healthy children and genetic correlations with ADHD

Cumulative evidence has showed the deficits of inhibitory control in patients with attention deficit hyperactivity disorder (ADHD), which is considered as an endophenotype of ADHD. Genetic study of inhibitory control could advance gene discovery and further facilitate the understanding of ADHD genetic basis, but the studies were limited in both the general population and ADHD patients. To reveal genetic risk variants of inhibitory control and its potential genetic relationship with ADHD, we conducted genome-wide association studies (GWAS) on inhibitory control using three datasets, which included 783 and 957 ADHD patients and 1350 healthy children. Subsequently, we employed polygenic risk scores (PRS) to explore the association of inhibitory control with ADHD and related psychiatric disorders. Firstly, we identified three significant loci for inhibitory control in the healthy dataset, two loci in the case dataset, and one locus in the meta-analysis of three datasets. Besides, we found more risk genes and variants by applying transcriptome-wide association study (TWAS) and conditional FDR method. Then, we constructed a network by connecting the genes identified in our study, leading to the identification of several vital genes. Lastly, we identified a potential relationship between inhibitory control and ADHD and autism by PRS analysis and found the direct and mediated contribution of the identified genetic loci on ADHD symptoms by mediation analysis. In conclusion, we revealed some genetic risk variants associated with inhibitory control and elucidated the benefit of inhibitory control as an endophenotype, providing valuable insights into the mechanisms underlying ADHD.

Monozygotic twins discordant for schizophrenia differ in maturation and synaptic transmission

Schizophrenia affects approximately 1% of the world population. Genetics, epigenetics, and environmental factors are known to play a role in this psychiatric disorder. While there is a high concordance in monozygotic twins, about half of twin pairs are discordant for schizophrenia. To address the question of how and when concordance in monozygotic twins occur, we have obtained fibroblasts from two pairs of schizophrenia discordant twins (one sibling with schizophrenia while the second one is unaffected by schizophrenia) and three pairs of healthy twins (both of the siblings are healthy). We have prepared iPSC models for these 3 groups of patients with schizophrenia, unaffected co-twins, and the healthy twins. When the study started the co-twins were considered healthy and unaffected but both the co-twins were later diagnosed with a depressive disorder. The reprogrammed iPSCs were differentiated into hippocampal neurons to measure the neurophysiological abnormalities in the patients. We found that the neurons derived from the schizophrenia patients were less arborized, were hypoexcitable with immature spike features, and exhibited a significant reduction in synaptic activity with dysregulation in synapse-related genes. Interestingly, the neurons derived from the co-twin siblings who did not have schizophrenia formed another distinct group that was different from the neurons in the group of the affected twin siblings but also different from the neurons in the group of the control twins. Importantly, their synaptic activity was not affected. Our measurements that were obtained from schizophrenia patients and their monozygotic twin and compared also to control healthy twins point to hippocampal synaptic deficits as a central mechanism in schizophrenia.

Cellular adhesion molecules in drug-naïve and previously medicated patients with schizophrenia-spectrum disorders

BACKGROUND

Endothelial inflammation may be involved in the pathogenesis of schizophrenia, and cellular adhesion molecules (CAMs) on endothelial cells may facilitate leukocyte binding and transendothelial migration of cells and inflammatory factors. The aim of the present study was to assess levels of soluble cellular adhesion molecules, including intercellular adhesion molecule (ICAM)-1, vascular adhesion molecule (VCAM)-1, mucosal addressin cell adhesion molecule (MADCAM), junctional adhesion molecule (JAM-A) and neural cadherin (N-CAD) in patients with schizophrenia compared to healthy controls.

METHODS

The study population consists of 138 patients with schizophrenia-spectrum disorder, of whom 54 were drug-naïve, compared to 317 general population controls. The potential confounders age, gender, smoking and body mass index (BMI) were adjusted for in linear regression models.

RESULTS

The total patient group showed significantly higher levels of ICAM-1 (p < 0.001) and VCAM-1 (p < 0.001) compared to controls. Previously medicated patients showed higher ICAM-1 levels compared to drug-naïve patients (p = 0.042) and controls (p < 0.001), and elevated VCAM-1 levels compared to controls (p < 0.001). Drug-naive patients had elevated levels of VCAM-1 (p = 0.031) compared to controls.

CONCLUSIONS

In our study, patients with schizophrenia - including the drug-naïve - have higher levels of soluble CAMs compared to healthy controls. These findings suggest activation of the endothelial system as in inflammation.

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.

Unravelling the genetic basis of Schizophrenia

Neuronal development is a highly regulated mechanism that is central to organismal function in animals. In humans, disruptions to this process can lead to a range of neurodevelopmental phenotypes, including Schizophrenia (SCZ). SCZ has a significant genetic component, whereby an individual with an SCZ affected family member is eight times more likely to develop the disease than someone with no family history of SCZ. By examining a combination of genomic, transcriptomic and epigenomic datasets, large-scale 'omics' studies aim to delineate the relationship between genetic variation and abnormal cellular activity in the SCZ brain. Herein, we provide a brief overview of some of the key omics methods currently being used in SCZ research, including RNA-seq, the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and high-throughput chromosome conformation capture (3C) approaches (e.g., Hi-C), as well as single-cell/nuclei iterations of these methods. We also discuss how these techniques are being employed to further our understanding of the genetic basis of SCZ, and to identify associated molecular pathways, biomarkers, and candidate drug targets.

Role of NT5DC2 in tyrosine hydroxylase phosphorylation based on the analysis of NT5DC2-binding proteins

The gene encoding 5'-nucleotidase domain-containing protein 2 (NT5DC2) has been associated with neuropsychiatric disorders related to the abnormality of dopamine activity in the brain. However, its physiological functions remain unclear. In this study, we analyzed the features of NT5DC2 that influence its binding with tyrosine hydroxylase (TH) and its effects on dihydroxyphenylalanine (DOPA) synthesis, using NT5DC2 overexpressed in PC12D cells by the pCMV vector. Western blot analysis revealed that the purified NT5DC2-DYKDDDDK-tag (NT5DC2-tag) protein can bind with the phosphorylated form of recombinant human TH type 1 (rhTH1), apart from the endogenous TH in PC12D cells. Proteomic analysis by mass spectrometry revealed that the purified NT5DC2-tag protein has the potential to bind to 41 proteins with multiple phosphorylation sites in PC12D cells (NT5DC2 binding proteins: positive, 391 sites/41 proteins; and negative, 85 sites/27 proteins). Overexpression of NT5DC2 in PC12D cells decreased DOPA levels in the medium. When the lysate of PC12D cells overexpressing NT5DC2 was incubated at 37 °C, the phosphorylated form of endogenous TH in PC12D cells decreased. This decrease was also detected when phosphorylated rhTH1 was incubated with purified NT5DC2-tag. Overall, our results suggest that NT5DC2 regulates DOPA synthesis by promoting the dephosphorylation of TH, similar to a phosphatase. Therefore, our study provides useful information for understanding various disorders associated with abnormalities in dopamine levels in the brain.

Comparative genomics of macaques and integrated insights into genetic variation and population history

The crab-eating macaques ( Macaca fascicularis ) and rhesus macaques ( M. mulatta ) are widely studied nonhuman primates in biomedical and evolutionary research. Despite their significance, the current understanding of the complex genomic structure in macaques and the differences between species requires substantial improvement. Here, we present a complete genome assembly of a crab-eating macaque and 20 haplotype-resolved macaque assemblies to investigate the complex regions and major genomic differences between species. Segmental duplication in macaques is ∼42% lower, while centromeres are ∼3.7 times longer than those in humans. The characterization of ∼2 Mbp fixed genetic variants and ∼240 Mbp complex loci highlights potential associations with metabolic differences between the two macaque species (e.g., CYP2C76 and EHBP1L1 ). Additionally, hundreds of alternative splicing differences show post-transcriptional regulation divergence between these two species (e.g., PNPO ). We also characterize 91 large-scale genomic differences between macaques and humans at a single-base-pair resolution and highlight their impact on gene regulation in primate evolution (e.g., FOLH1 and PIEZO2 ). Finally, population genetics recapitulates macaque speciation and selective sweeps, highlighting potential genetic basis of reproduction and tail phenotype differences (e.g., STAB1 , SEMA3F , and HOXD13 ). In summary, the integrated analysis of genetic variation and population genetics in macaques greatly enhances our comprehension of lineage-specific phenotypes, adaptation, and primate evolution, thereby improving their biomedical applications in human diseases.

Transcriptional control of CCAAT/enhancer binding protein zeta gene in chicken adipose tissue

CCAAT/enhancer binding protein zeta (C/EBPZ) was differentially expressed in abdominal adipose tissues of fat and lean broilers and regulated adipogenesis in chicken. The objective of this study was to elucidate the transcriptional regulation of C/EBPZ gene in chicken adipose tissue. A 2,031-base pair (bp) chicken C/EBPZ sequence (2,025 nucleotides upstream to 6 nucleotides downstream from the initiator codon, -2,025/+6) was studied. The sequence exhibited a significant promoter activity (P < 0.05) and had some cis-acting elements, notably, a core promoter was identified in nucleotides -94 to +6. Additionally, DNA pull-down assay showed that proteins interacted with chicken C/EBPZ promoter (-173/+6) in preadipocytes were implicated in transcription, post-transcriptional regulation and translation. In addition, KLF2 facilitated the activities of chicken C/EBPZ promoter (-2,025/+6, -1,409/+6, -793/+6, -485/+6, -173/+6, and -94/+6) in preadipocytes (P < 0.05). The expression levels of KLF2 and C/EBPZ in chicken abdominal adipose tissue were substantially associated (r = 0.5978278, P < 0.0001), and KLF2 increased C/EBPZ expression in vitro (P < 0.05). Additionally, chromatin immunoprecipitation (ChIP)-PCR analysis revealed that KLF2 has the ability to interact with the chicken C/EBPZ promoter regions at least at the positions -1,245/-1,048 and -571/-397. Mutation analysis showed that the CGCAGCGCCCG motif located in the chicken C/EBPZ promoter at positions -45 to -35 is involved in regulating transcription and facilitates trans activation by KLF2. These results provided some information of transcription control of C/EBPZ in chicken adipose tissue.

Nuclear microRNA-mediated transcriptional control determines adult microglial homeostasis and brain function

Microglia are versatile regulators in brain development and disorders. Emerging evidence links microRNA (miRNA)-mediated regulation to microglial function; however, the exact underlying mechanism remains largely unknown. Here, we uncover the enrichment of miR-137, a neuropsychiatric-disorder-associated miRNA, in the microglial nucleus, and reveal its unexpected nuclear functions in maintaining the microglial global transcriptomic state, phagocytosis, and inflammatory response. Mechanistically, microglial Mir137 deletion increases chromatin accessibility, which contains binding motifs for the microglial master transcription factor Pu.1. Through biochemical and bioinformatics analyses, we propose that miR-137 modulates Pu.1-mediated gene expression by suppressing Pu.1 binding to chromatin. Importantly, we find that increased Pu.1 binding upregulates the target gene Jdp2 (Jun dimerization protein 2) and that knockdown of Jdp2 significantly suppresses the impaired phagocytosis and pro-inflammatory response in Mir137 knockout microglia. Collectively, our study provides evidence supporting the notion that nuclear miR-137 acts as a transcriptional modulator and that this microglia-specific function is essential for maintaining normal adult brain function.

Scaled and Efficient Derivation of Loss of Function Alleles in Risk Genes for Neurodevelopmental and Psychiatric Disorders in Human iPSC

Translating genetic findings for neurodevelopmental and psychiatric disorders (NPD) into actionable disease biology would benefit from large-scale and unbiased functional studies of NPD genes. Leveraging the cytosine base editing (CBE) system, here we developed a pipeline for clonal loss-of-function (LoF) allele mutagenesis in human induced pluripotent stem cells (hiPSCs) by introducing premature stop-codons (iSTOP) that lead to mRNA nonsense-mediated-decay (NMD) or protein truncation. We tested the pipeline for 23 NPD genes on 3 hiPSC lines and achieved highly reproducible, efficient iSTOP editing in 22 NPD genes. Using RNAseq, we confirmed their pluripotency, absence of chromosomal abnormalities, and NMD. Interestingly, for three schizophrenia risk genes (SETD1A, TRIO, CUL1), despite the high efficiency of base editing, we only obtained heterozygous LoF alleles, suggesting their essential roles for cell growth. We replicated the reported neural phenotypes of SHANK3-haploinsufficiency and found CUL1-LoF reduced neurite branches and synaptic puncta density. This iSTOP pipeline enables a scaled and efficient LoF mutagenesis of NPD genes, yielding an invaluable shareable resource.

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.

Developmental psychopathology as a meta-paradigm: From zero-sum science to epistemological pluralism in theory and research

In a thoughtful commentary in this journal a decade ago, Michael Rutter reviewed 25 years of progress in the field before concluding that developmental psychopathology (DP) initiated a paradigm shift in clinical science. This deduction requires that DP itself be a paradigm. According to Thomas Kuhn, canonical paradigms in the physical sciences serve unifying functions by consolidating scientists' thinking and scholarship around single, closed sets of discipline-defining epistemological assumptions and methods. Paradigm shifts replace these assumptions and methods with a new field-defining framework. In contrast, the social sciences are multiparadigmatic, with thinking and scholarship unified locally around open sets of epistemological assumptions and methods with varying degrees of inter-, intra-, and subdisciplinary reach. DP challenges few if any of these local paradigms. Instead, DP serves an essential pluralizing function, and is therefore better construed as a metaparadigm. Seen in this way, DP holds tremendous untapped potential to move the field from zero-sum thinking and scholarship to positive-sum science and epistemological pluralism. This integrative vision, which furthers Dante Cicchetti's legacy of interdisciplinarity, requires broad commitment among scientists to reject zero-sum scholarship in which portending theories, useful principles, and effective interventions are jettisoned based on confirmation bias, errors in logic, and ideology.

Causal Relationship Between Ferritin and Neuropsychiatric Disorders: A Two-Sample Mendelian Randomization Study

Background

Previous observational research has indicated a correlation between ferritin levels and neuropsychiatric disorders, although the causal relationship remains uncertain.

Objective

The objective of this study was to investigate the potential causal link between plasma ferritin levels and neuropsychiatric disorders.

Methods

A two-sample Mendelian randomization (MR) study was conducted, wherein genetic instruments associated with ferritin were obtained from a previously published genome-wide association study (GWAS). Summary statistics pertaining to neuropsychiatric disorders were derived from five distinct GWAS datasets. The primary MR analysis employed the inverse variance weighted (IVW) method and was corroborated by additional methods including MR-Egger, weighted median, simple mode, and weighted mode. Sensitivity analyses were employed to identify potential pleiotropy and heterogeneity in the results.

Results

The fixed effects IVW method revealed a statistically significant causal relationship between plasma ferritin level and the occurrence of Alzheimer's disease (odds ratio [OR] = 1.06, 95% confidence interval [CI]: 1.00-1.12, p = 0.037), as well as Parkinson's disease (OR = 1.06, 95% CI: 1.00-1.13, p = 0.041). Various sensitivity analyses were conducted, which demonstrated no substantial heterogeneity or pleiotropy. Conversely, no compelling evidence was found to support a causal association between ferritin and amyotrophic lateral sclerosis, schizophrenia, or major depressive disorder.

Conclusions

This MR study provides evidence at the genetic level for a causal relationship between plasma ferritin and an increased risk of Alzheimer's disease and Parkinson's disease. The exact genetic mechanisms underlying this connection necessitate further investigation.

Polygenic effects on brain functional endophenotype for deficit and non-deficit schizophrenia

Deficit schizophrenia (DS) is a subtype of schizophrenia (SCZ). The polygenic effects on the neuroimaging alterations in DS still remain unknown. This study aims to calculate the polygenic risk scores for schizophrenia (PRS-SCZ) in DS, and further explores the potential associations with functional features of brain. PRS-SCZ was calculated according to the Whole Exome sequencing and Genome-wide association studies (GWAS). Resting-state fMRI, as well as biochemical features and neurocognitive data were obtained from 33 DS, 47 NDS and 41 HCs, and association studies of genetic risk with neuroimaging were performed in this sample. The analyses of amplitude of low-frequency fluctuation (ALFF), regional homogeneity (ReHo) and functional connectivity (FC) were performed to detect the functional alterations between DS and NDS. In addition, correlation analysis was used to investigate the relationships between functional features (ALFF, ReHo, FC) and PRS-SCZ. The PRS-SCZ of DS was significantly lower than that in NDS and HC. Compared to NDS, there was a significant increase in the ALFF of left inferior temporal gyrus (ITG.L) and left inferior frontal gyrus (IFG.L) and a significant decrease in the ALFF of right precuneus (PCUN.R) and ReHo of right middle frontal gyrus (MFG.R) in DS. FCs were widely changed between DS and NDS, mainly concentrated in default mode network, including ITG, PCUN and angular gyrus (ANG). Correlation analysis revealed that the ALFF of left ITG, the ReHo of right middle frontal gyrus, the FC value between insula and ANG, left ITG and right corpus callosum, left ITG and right PCUN, as well as the scores of Trail Making Test-B, were associated with PRS-SCZ in DS. The present study demonstrated the differential polygenic effects on functional changes of brain in DS and NDS, providing a potential neuroimaging-genetic perspective for the pathogenesis of schizophrenia.

rvTWAS: identifying gene-trait association using sequences by utilizing transcriptome-directed feature selection

Toward the identification of genetic basis of complex traits, transcriptome-wide association study (TWAS) is successful in integrating transcriptome data. However, TWAS is only applicable for common variants, excluding rare variants in exome or whole-genome sequences. This is partly because of the inherent limitation of TWAS protocols that rely on predicting gene expressions. Our previous research has revealed the insight into TWAS: the 2 steps in TWAS, building and applying the expression prediction models, are essentially genetic feature selection and aggregations that do not have to involve predictions. Based on this insight disentangling TWAS, rare variants' inability of predicting expression traits is no longer an obstacle. Herein, we developed "rare variant TWAS," or rvTWAS, that first uses a Bayesian model to conduct expression-directed feature selection and then uses a kernel machine to carry out feature aggregation, forming a model leveraging expressions for association mapping including rare variants. We demonstrated the performance of rvTWAS by thorough simulations and real data analysis in 3 psychiatric disorders, namely schizophrenia, bipolar disorder, and autism spectrum disorder. We confirmed that rvTWAS outperforms existing TWAS protocols and revealed additional genes underlying psychiatric disorders. Particularly, we formed a hypothetical mechanism in which zinc finger genes impact all 3 disorders through transcriptional regulations. rvTWAS will open a door for sequence-based association mappings integrating gene expressions.

miR-137 regulates PTP61F, affecting insulin signaling, metabolic homeostasis, and starvation resistance in Drosophila

miR-137 is a highly conserved brain-enriched microRNA (miRNA) that has been associated with neuronal function and proliferation. Here, we show that Drosophila miR-137 null mutants display increased body weight with enhanced triglyceride content and decreased locomotor activity. In addition, when challenged by nutrient deprivation, miR-137 mutants exhibit reduced motivation to feed and prolonged survival. We show through genetic epistasis and rescue experiments that this starvation resistance is due to a disruption in insulin signaling. Our studies further show that miR-137 null mutants exhibit a drastic reduction in levels of the phosphorylated/activated insulin receptor, InR (InR-P). We investigated if this is due to the predicted miR-137 target, Protein Tyrosine Phosphatase 61F (PTP61F), ortholog of mammalian TC-PTP/PTP1B, which are known to dephosphorylate InR-P. Indeed, levels of an endogenously tagged GFP-PTP61F are significantly elevated in miR-137 null mutants, and we show that overexpression of PTP61F alone is sufficient to mimic many of the metabolic phenotypes of miR-137 mutants. Finally, we knocked-down elevated levels of PTP61F in the miR-137 null mutant background and show that this rescues levels of InR-P, restores normal body weight and triglyceride content, starvation sensitivity, as well as attenuates locomotor and starvation-induced feeding defects. Our study supports a model in which miR-137 is critical for dampening levels of PTP61F, thereby maintaining normal insulin signaling and energy homeostasis.

Linking Polygenic Risk of Schizophrenia to Variation in Magnetic Resonance Imaging Brain Measures: A Comprehensive Systematic Review

BACKGROUND AND HYPOTHESIS

Schizophrenia is highly heritable, with a polygenic effect of many genes conferring risk. Evidence on whether cumulative risk also predicts alterations in brain morphology and function is inconsistent. This systematic review examined evidence for schizophrenia polygenic risk score (sczPRS) associations with commonly used magnetic resonance imaging (MRI) measures. We expected consistent evidence to emerge for significant sczPRS associations with variation in structure and function, specifically in frontal, temporal, and insula cortices that are commonly implicated in schizophrenia pathophysiology.

STUDY DESIGN

In accordance with Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines, we searched MEDLINE, Embase, and PsycINFO for peer-reviewed studies published between January 2013 and March 2022. Studies were screened against predetermined criteria and National Institutes of Health (NIH) quality assessment tools.

STUDY RESULTS

In total, 57 studies of T1-weighted structural, diffusion, and functional MRI were included (age range = 9-80 years, Nrange = 64-76 644). We observed moderate, albeit preliminary, evidence for higher sczPRS predicting global reductions in cortical thickness and widespread variation in functional connectivity, and to a lesser extent, region-specific reductions in frontal and temporal volume and thickness. Conversely, sczPRS does not predict whole-brain surface area or gray/white matter volume. Limited evidence emerged for sczPRS associations with diffusion tensor measures of white matter microstructure in a large community sample and smaller cohorts of children and young adults. These findings were broadly consistent across community and clinical populations.

CONCLUSIONS

Our review supports the hypothesis that schizophrenia is a disorder of disrupted within and between-region brain connectivity, and points to specific whole-brain and regional MRI metrics that may provide useful intermediate phenotypes.

Identifying diseases associated with Post-COVID syndrome through an integrated network biology approach

A growing body of research shows that COVID-19 is now recognized as a multi-organ disease with a wide range of manifestations that can have long-lasting repercussions, referred to as post-COVID-19 syndrome. It is unknown why the vast majority of COVID-19 patients develop post-COVID-19 syndrome, or why patients with pre-existing disorders are more likely to experience severe COVID-19. This study used an integrated network biology approach to obtain a comprehensive understanding of the relationship between COVID-19 and other disorders. The approach involved building a PPI network with COVID-19 genes and identifying highly interconnected regions. The molecular information contained within these subnetworks, as well as the pathway annotations, were used to reveal the link between COVID-19 and other disorders. Using Fisher's exact test and disease-specific gene information, significant COVID-19-disease associations were discovered. The study discovered diseases that affect multiple organs and organ systems, thus proving the theory of multiple organ damage caused by COVID-19. Cancers, neurological disorders, hepatic diseases, cardiac disorders, pulmonary diseases, and hypertensive diseases are just a few of the conditions linked to COVID-19. Pathway enrichment analysis of shared proteins revealed the shared molecular mechanism of COVID-19 and these diseases. The findings of the study shed new light on the major COVID-19-associated disease conditions and how their molecular mechanisms interact with COVID-19. The novelty of studying disease associations in the context of COVID-19 provides new insights into the management of rapidly evolving long-COVID and post-COVID syndromes, which have significant global implications.Communicated by Ramaswamy H. Sarma.

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

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

Therapeutic potential for KCC2-targeted neurological diseases

Patients with neurological diseases, such as schizophrenia, tend to show low K+-Cl- co-transporter 2 (KCC2) levels in the brain. The cause of these diseases has been associated with stress and neuroinflammation. However, since the pathogenesis of these diseases is not yet fully investigated, drug therapy is still limited to symptomatic therapy. Targeting KCC2, which is mainly expressed in the brain, seems to be an appropriate approach in the treatment of these diseases. In this review, we aimed to discuss about stress and inflammation, KCC2 and Gamma-aminobutyric acid (GABA) function, diseases which decrease the KCC2 levels in the brain, factors that regulate KCC2 activity, and the possibility to overcome neuronal dysfunction targeting KCC2. We also aimed to discuss the relationships between neurological diseases and LPS caused by Porphyromonas gingivalis (P. g), which is a type of oral bacterium. Clinical trials on oxytocin, sirtuin 1 (SIRT1) activator, and transient receptor potential cation channel subfamily V Member 1 activator have been conducted to develop effective treatment methods. We believe that KCC2 modulators that regulate mitochondria, such as oxytocin, glycogen synthase kinase 3β (GSK3β), and SIRT1, can be potential targets for neurological diseases.

Effect of rs1108580 of DBH and rs1006737 of CACNA1C on Cognition and Tardive Dyskinesia in a North Indian Schizophrenia Cohort

Genetic perturbations in dopamine neurotransmission and calcium signaling pathways are implicated in the etiology of schizophrenia. We aimed to test the association of a functional splice variant each in Dopamine β-Hydroxylase (DBH; rs1108580) and Calcium voltage-gated channel subunit alpha1 C (CACNA1C; rs1006737) genes in these pathways with schizophrenia (506 cases, 443 controls); Abnormal Involuntary Movement Scale (AIMS) scores in subjects assessed for tardive dyskinesia (76 TD-positive, 95 TD-negative) and Penn Computerized Neurocognitive Battery (PennCNB) scores (334 cases, 234 controls). The effect of smoking status and SNP genotypes on AIMS scores were assessed using ANOVA; health status and SNP genotypes on three performance functions of PennCNB cognitive domains were assessed by ANCOVA with age and sex as covariates. Association with Positive and Negative Syndrome Scale (PANSS) scores in the TD cohort and cognitive scores in healthy controls of the cognition cohort were tested by linear regression. None of the markers were associated with schizophrenia. Smoking status [F(2, 139) = 10.6; p = 5 × 10-5], rs1006737 [F(2, 139) = 7.1; p = 0.001], TD status*smoking [F(2, 139) = 8.0; p = 5.0 × 10-4] and smoking status*rs1006737 [F(4, 139) = 2.7; p = 0.03] had an effect on AIMS score. Furthermore, rs1006737 was associated with orofacial [F(2, 139) = 4.6; p = 0.01] and limb-truncal TD [(F(2, 139) = 3.8; p = 0.02]. Main effect of rs1108580 on working memoryprocessing speed [F(2, 544) = 3.8; p = 0.03] and rs1006737 on spatial abilityefficiency [F(1, 550) = 9.4; p = 0.02] was identified. Health status*rs1006737 interaction had an effect on spatial memoryprocessing speed [F(1, 550) = 6.9; p = 0.01]. Allelic/genotypic association (p = 0.01/0.03) of rs1006737 with disorganized/concrete factor and allelic association of rs1108580 (p = 0.04) with a depressive factor of PANSS was observed in the TD-negative subcohort. Allelic association of rs1006737 with sensorimotor dexterityaccuracy (p = 0.03), attentionefficiency (p = 0.05), and spatial abilityefficiency (p = 0.02); allelic association of rs1108580 with face memoryaccuracy (p = 0.05) and emotionefficiency (p = 0.05); and allelic/genotypic association with emotionaccuracy (p = 0.003/0.009) were observed in healthy controls of the cognition cohort. These association findings may have direct implications for personalized medicine and cognitive remediation.

New insights into the structure and function of CNNM proteins

Magnesium (Mg2+ ) is the most abundant divalent cation in cells and plays key roles in almost all biological processes. CBS-pair domain divalent metal cation transport mediators (CNNMs) are a newly characterized class of Mg2+ transporters present throughout biology. Originally discovered in bacteria, there are four CNNM proteins in humans, which are involved in divalent cation transport, genetic diseases, and cancer. Eukaryotic CNNMs are composed of four domains: an extracellular domain, a transmembrane domain, a cystathionine-β-synthase (CBS)-pair domain, and a cyclic nucleotide-binding homology domain. The transmembrane and CBS-pair core are the defining features of CNNM proteins with over 20 000 protein sequences known from over 8000 species. Here, we review the structural and functional studies of eukaryotic and prokaryotic CNNMs that underlie our understanding of their regulation and mechanism of ion transport. Recent structures of prokaryotic CNNMs confirm the transmembrane domain mediates ion transport with the CBS-pair domain likely playing a regulatory role through binding divalent cations. Studies of mammalian CNNMs have identified new binding partners. These advances are driving progress in understanding this deeply conserved and widespread family of ion transporters.

Evaluating 17 methods incorporating biological function with GWAS summary statistics to accelerate discovery demonstrates a tradeoff between high sensitivity and high positive predictive value

Where sufficiently large genome-wide association study (GWAS) samples are not currently available or feasible, methods that leverage increasing knowledge of the biological function of variants may illuminate discoveries without increasing sample size. We comprehensively evaluated 17 functional weighting methods for identifying novel associations. We assessed the performance of these methods using published results from multiple GWAS waves across each of five complex traits. Although no method achieved both high sensitivity and positive predictive value (PPV) for any trait, a subset of methods utilizing pleiotropy and expression quantitative trait loci nominated variants with high PPV (>75%) for multiple traits. Application of functionally weighting methods to enhance GWAS power for locus discovery is unlikely to circumvent the need for larger sample sizes in truly underpowered GWAS, but these results suggest that applying functional weighting to GWAS can accurately nominate additional novel loci from available samples for follow-up studies.

One versus two biological parents with mental disorders: Relationship to educational attainment in the next generation

BACKGROUND

Both maternal and, separately, paternal mental illness are associated with diminished academic attainment among children. However, the differential impacts of diagnostic type and degree of parental burden (e.g. one v. both parents affected) on these functional outcomes are unknown.

METHODS

Using the Swedish national patient (NPR) and multi-generation (MGR) registers, 2 226 451 children (1 290 157 parental pairs), born 1 January 1973-31 December 1997, were followed through 31 December 2013. Diagnostic status of all cohort members was defined for eleven psychiatric disorders, and families classed by exposure: (1) parents affected with any disorder, (2) parents affected with a disorder group (e.g. neuropsychiatric disorders), and (3) parents affected with a specific disorder (e.g. ADHD). Pairs were further defined as 'unaffected,' 'single-affected,', or 'dual-affected.' Among offspring, the study evaluated fulfillment of four academic milestones, from compulsory (primary) school through University (college). Sensitivity analyses considered the impact of child's own mental health, as well as parental education, on main effects.

RESULTS

Marked reductions in the odds of achievement were observed, emerging at the earliest levels of schooling for both single-affected [adjusted odds ratio (aOR), 0.50; 95% CI 0.49-0.51] and dual-affected (aOR 0.29, 95% CI 0.28-0.30) pairs and persisting thereafter [aOR range (single), 0.52-0.65; aOR range (dual), 0.30-0.40]. This pattern was repeated for analyses within diagnosis/diagnostic group. Main results were robust to adjustment for offspring mental health and parent education level.

CONCLUSIONS

Parental mental illness is associated with profound reductions in educational attainment in the subsequent generation, with children from dual-affected families at uniquely high risk.

A higher dysregulation burden of brain DNA methylation in female patients implicated in the sex bias of Schizophrenia

Sex differences are pervasive in schizophrenia (SCZ), but the extent and magnitude of DNA methylation (DNAm) changes underlying these differences remain uncharacterized. In this study, sex-stratified differential DNAm analysis was performed in postmortem brain samples from 117 SCZ and 137 controls, partitioned into discovery and replication datasets. Three differentially methylated positions (DMPs) were identified (adj.p < 0.05) in females and 29 DMPs in males without overlap between them. Over 81% of these sex-stratified DMPs were directionally consistent between sexes but with different effect sizes. Females experienced larger magnitude of DNAm changes and more DMPs (based on data of equal sample size) than males, contributing to a higher dysregulation burden of DNAm in females SCZ. Additionally, despite similar proportions of female-related DMPs (fDMPs, 8%) being under genetic control compared with males (10%), significant enrichment of DMP-related single nucleotide polymorphisms (SNPs) in signals of genome-wide association studies was identified only in fDMPs. One DMP in each sex connected the SNPs and gene expression of CALHM1 in females and CCDC149 in males. PPI subnetworks revealed that both female- and male-related differential DNAm interacted with synapse-related dysregulation. Immune-related pathways were unique for females and neuron-related pathways were associated with males. This study reveals remarkable quantitative differences in DNAm-related sexual dimorphism in SCZ and that females have a higher dysregulation burden of SCZ-associated DNAm than males.

MicroRNA137-loaded lipid nanoparticles regulate synaptic proteins in the prefrontal cortex

Genome-wide association studies indicate that allele variants in MIR137, the host gene of microRNA137 (miR137), confer an increased risk of schizophrenia (SCZ). Aberrant expression of miR137 and its targets, many of which regulate synaptic functioning, are also associated with an increased risk of SCZ. Thus, miR137 represents an attractive target aimed at correcting the molecular basis for synaptic dysfunction in individuals with high genetic risk for SCZ. Advancements in nanotechnology utilize lipid nanoparticles (LNPs) to transport and deliver therapeutic RNA. However, there remains a gap in using LNPs to regulate gene and protein expression in the brain. To study the delivery of nucleic acids by LNPs to the brain, we found that LNPs released miR137 cargo and inhibited target transcripts of interest in neuroblastoma cells. Biodistribution of LNPs loaded with firefly luciferase mRNA remained localized to the mouse prefrontal cortex (PFC) injection site without circulating to off-target organs. LNPs encapsulating Cre mRNA preferentially co-expressed in neuronal over microglial or astrocytic cells. Using quantitative proteomics, we found miR137 modulated glutamatergic synaptic protein networks that are commonly dysregulated in SCZ. These studies support engineering the next generation of brain-specific LNPs to deliver RNA therapeutics and improve symptoms of central nervous system disorders.

A novel risk variant block across introns 36-45 of CACNA1C for schizophrenia: a cohort-wise replication and cerebral region-wide validation study

OBJECTIVES

Numerous genome-wide association studies have identified CACNA1C as one of the top risk genes for schizophrenia. As a necessary post-genome-wide association study (GWAS) follow-up, here, we focused on this risk gene, carefully investigated its novel risk variants for schizophrenia, and explored their potential functions.

METHODS

We analyzed four independent samples (including three European and one African-American) comprising 5648 cases and 6936 healthy subjects to identify replicable single nucleotide polymorphism-schizophrenia associations. The potential regulatory effects of schizophrenia-risk alleles on CACNA1C mRNA expression in 16 brain regions (n = 348), gray matter volumes (GMVs) of five subcortical structures (n = 34 431), and surface areas and thickness of 34 cortical regions (n = 36 936) were also examined.

RESULTS

A novel 17-variant block across introns 36-45 of CACNA1C was significantly associated with schizophrenia in the same effect direction across at least two independent samples (1.8 × 10-4 ≤ P ≤ 0.049). Most risk variants within this block showed significant associations with CACNA1C mRNA expression (1.6 × 10-3 ≤ P ≤ 0.050), GMVs of subcortical structures (0.016 ≤ P ≤ 0.048), cortical surface areas (0.010 ≤ P ≤ 0.050), and thickness (0.004 ≤ P ≤ 0.050) in multiple brain regions.

CONCLUSION

We have identified a novel and functional risk variant block at CACNA1C for schizophrenia, providing further evidence for the important role of this gene in the pathogenesis of schizophrenia.

Polygenic liabilities underlying job stress and exhaustion over a 10-year follow-up: A general population study

We investigated whether individuals, who have a high polygenic loading for schizophrenia and major depression (PGL) but have not developed the respective disorders, are still susceptible to experience milder forms of ill-being in terms of job strain or exhaustion. We used the population-based Young Finns Study data (n = 928). PGL was assessed with a cumulative score of the polygenic risk scores for schizophrenia and depression. Participants (24-49-year-olds) evaluated their exhaustion levels and perceived job characteristics over a 10-year follow-up (2001, 2007, 2011). Participants with diagnosed psychotic or affective disorders were excluded. We found that high PGL did not predict less favorable perceptions of job environment (job strain, demands, control, satisfaction, social support at work) but high PGL predicted a higher trajectory of exhaustion in early adulthood and middle age. Additionally, high (vs. low) PGL predicted a stronger increase in exhaustion at increased levels of job strain. These findings remained after controlling for sex, socioeconomic factors, health behaviors, and cognitive performance. In conclusion, individuals with high PGL may have an elevated liability to experience exhaustion especially in early adulthood and middle age (despite they perceive their job environment similarly than others), and especially and at high levels of job strain.

Genetic Influences on Outcomes of Psychotherapy in Borderline Personality Disorder: A Narrative Review of Implications for Personalized Treatment

Borderline personality disorder (BPD) manifests as instability in mood, relationships, self-image, and behavior, representing a challenging mental health issue. This review scrutinizes genetic factors influencing BPD and the corresponding treatment outcomes. The primary objective of this narrative review is to illuminate the association between genetic factors and BPD treatment outcomes, discussing the potential of genetic testing for personalized therapy. The review is derived from observational and experimental studies on BPD, genetic factors, and psychotherapy from 2000 to 2023, sourced primarily through PubMed. Reviews and meta-analyses were excluded. Our review suggests that genetic factors account for 40-60% of BPD variation, with significant roles played by epigenetic alterations like DNA methylation and microRNAs, particularly in the context of childhood trauma. Gene-environment interactions are also vital for BPD's development. Treatments such as dialectical behavior therapy, mentalization-based therapy, and schema therapy have shown efficacy, with success variability possibly linked to genetic factors. However, existing research is constrained by recall bias, diverse methodologies, and limited sample sizes. Future research necessitates long-term follow-up, diverse populations, and controlled variables to enhance our comprehension of BPD treatment outcomes' genetic foundations. The review underlines the promise of personalized medicine in BPD treatment, driven by genetic insights.

The prediction of Alzheimer's disease through multi-trait genetic modeling

To better capture the polygenic architecture of Alzheimer's disease (AD), we developed a joint genetic score, MetaGRS. We incorporated genetic variants for AD and 24 other traits from two independent cohorts, NACC (n = 3,174, training set) and UPitt (n = 2,053, validation set). One standard deviation increase in the MetaGRS is associated with about 57% increase in the AD risk [hazard ratio (HR) = 1.577, p = 7.17 E-56], showing little difference from the HR for AD GRS alone (HR = 1.579, p = 1.20E-56), suggesting similar utility of both models. We also conducted APOE-stratified analyses to assess the role of the e4 allele on risk prediction. Similar to that of the combined model, our stratified results did not show a considerable improvement of the MetaGRS. Our study showed that the prediction power of the MetaGRS significantly outperformed that of the reference model without any genetic information, but was effectively equivalent to the prediction power of the AD GRS.

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.

A pangenome reference of 36 Chinese populations

Human genomics is witnessing an ongoing paradigm shift from a single reference sequence to a pangenome form, but populations of Asian ancestry are underrepresented. Here we present data from the first phase of the Chinese Pangenome Consortium, including a collection of 116 high-quality and haplotype-phased de novo assemblies based on 58 core samples representing 36 minority Chinese ethnic groups. With an average 30.65× high-fidelity long-read sequence coverage, an average contiguity N50 of more than 35.63 megabases and an average total size of 3.01 gigabases, the CPC core assemblies add 189 million base pairs of euchromatic polymorphic sequences and 1,367 protein-coding gene duplications to GRCh38. We identified 15.9 million small variants and 78,072 structural variants, of which 5.9 million small variants and 34,223 structural variants were not reported in a recently released pangenome reference¹. The Chinese Pangenome Consortium data demonstrate a remarkable increase in the discovery of novel and missing sequences when individuals are included from underrepresented minority ethnic groups. The missing reference sequences were enriched with archaic-derived alleles and genes that confer essential functions related to keratinization, response to ultraviolet radiation, DNA repair, immunological responses and lifespan, implying great potential for shedding new light on human evolution and recovering missing heritability in complex disease mapping.

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

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

Predicting eating disorder and anxiety symptoms using disorder-specific and transdiagnostic polygenic scores for anorexia nervosa and obsessive-compulsive disorder

BACKGROUND

Clinical, epidemiological, and genetic findings support an overlap between eating disorders, obsessive-compulsive disorder (OCD), and anxiety symptoms. However, little research has examined the role of genetics in the expression of underlying phenotypes. We investigated whether the anorexia nervosa (AN), OCD, or AN/OCD transdiagnostic polygenic scores (PGS) predict eating disorder, OCD, and anxiety symptoms in a large developmental cohort in a sex-specific manner.

METHODS

Using summary statistics from Psychiatric Genomics Consortium AN and OCD genome-wide association studies, we conducted an AN/OCD transdiagnostic genome-wide association meta-analysis. We then calculated AN, OCD, and AN/OCD PGS in participants from the Avon Longitudinal Study of Parents and Children to predict eating disorder, OCD, and anxiety symptoms, stratified by sex (combined N = 3212-5369 per phenotype).

RESULTS

The PGS prediction of eating disorder, OCD, and anxiety phenotypes differed between sexes, although effect sizes were small. AN and AN/OCD PGS played a more prominent role in predicting eating disorder and anxiety risk than OCD PGS, especially in girls. AN/OCD PGS provided a small boost over AN PGS in the prediction of some anxiety symptoms. All three PGS predicted higher compulsive exercise across different developmental timepoints [β = 0.03 (s.e. = 0.01) for AN and AN/OCD PGS at age 14; β = 0.05 (s.e. = 0.02) for OCD PGS at age 16] in girls.

CONCLUSIONS

Compulsive exercise may have a transdiagnostic genetic etiology, and AN genetic risk may play a role in the presence of anxiety symptoms. Converging with prior twin literature, our results also suggest that some of the contribution of genetic risk may be sex-specific.

Genes implicated by a methylome-wide schizophrenia study in neonatal blood show differential expression in adult brain samples

Schizophrenia is a disabling disorder involving genetic predisposition in combination with environmental influences that likely act via dynamic alterations of the epigenome and the transcriptome but its detailed pathophysiology is largely unknown. We performed cell-type specific methylome-wide association study of neonatal blood (N = 333) from individuals who later in life developed schizophrenia and controls. Suggestively significant associations (P < 1.0 × 10-6) were detected in all cell-types and in whole blood with methylome-wide significant associations in monocytes (P = 2.85 × 10-9-4.87 × 10-9), natural killer cells (P = 1.72 × 10-9-7.82 × 10-9) and B cells (P = 3.8 × 10-9). Validation of methylation findings in post-mortem brains (N = 596) from independent schizophrenia cases and controls showed significant enrichment of transcriptional differences (enrichment ratio = 1.98-3.23, P = 2.3 × 10-3-1.0 × 10-5), with specific highly significant differential expression for, for example, BDNF (t = -6.11, P = 1.90 × 10-9). In addition, expression difference in brain significantly predicted schizophrenia (multiple correlation = 0.15-0.22, P = 3.6 × 10-4-4.5 × 10-8). In summary, using a unique design combining pre-disease onset (neonatal) blood methylomic data and post-disease onset (post-mortem) brain transcriptional data, we have identified genes of likely functional relevance that are associated with schizophrenia susceptibility, rather than confounding disease associated artifacts. The identified loci may be of clinical value as a methylation-based biomarker for early detection of increased schizophrenia susceptibility.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Brain transcriptome-wide association study implicates novel risk genes underlying schizophrenia risk

BACKGROUND

To identify risk genes whose expression are regulated by the reported risk variants and to explore the potential regulatory mechanism in schizophrenia (SCZ).

METHODS

We systematically integrated three independent brain expression quantitative traits (eQTLs) (CommonMind, GTEx, and BrainSeq Phase 2, a total of 1039 individuals) and GWAS data (56 418 cases and 78 818 controls), with the use of transcriptome-wide association study (TWAS). Diffusion magnetic resonance imaging was utilized to quantify the integrity of white matter bundles and determine whether polygenic risk of novel genes linked to brain structure was present in patients with first-episode antipsychotic SCZ.

RESULTS

TWAS showed that eight risk genes (CORO7, DDAH2, DDHD2, ELAC2, GLT8D1, PCDHA8, THOC7, and TYW5) reached transcriptome-wide significance (TWS) level. These findings were confirmed by an independent integrative approach (i.e. Sherlock). We further conducted conditional analyses and identified the potential risk genes that driven the TWAS association signal in each locus. Gene expression analysis showed that several TWS genes (including CORO7, DDAH2, DDHD2, ELAC2, GLT8D1, THOC7 and TYW5) were dysregulated in the dorsolateral prefrontal cortex of SCZ cases compared with controls. TWS genes were mainly expressed on the surface of glutamatergic neurons, GABAergic neurons, and microglia. Finally, SCZ cases had a substantially greater TWS genes-based polygenic risk (PRS) compared to controls, and we showed that fractional anisotropy of the cingulum-hippocampus mediates the influence of TWS genes PRS on SCZ.

CONCLUSIONS

Our findings identified novel SCZ risk genes and highlighted the importance of the TWS genes in frontal-limbic dysfunctions in SCZ, indicating possible therapeutic targets.

Genetic variations in DOCK4 contribute to schizophrenia susceptibility in a Chinese cohort: A genetic neuroimaging study

BACKGROUND

Emerging evidence suggests that the DOCK4 gene increases susceptibility to schizophrenia. However, no study has hitherto repeated this association in Chinese, and further investigated the relationship between DOCK4 and clinical symptoms in schizophrenic patients using clinical scales and functional magnetic resonance imaging (fMRI).

METHODS

In this study, we genotyped three single nucleotide polymorphisms (SNPs) (rs2074127, rs2217262, and rs2074130) within the DOCK4 gene using a case-control design (including 1289 healthy controls and 1351 patients with schizophrenia). 55 first-episode schizophrenia (FES) patients and 59 healthy participants were divided by the genotypes of rs2074130 into CC and CT+TT groups. We further investigated the association with clinical symptoms and neural characteristics (brain activation/connectivity and nodal network metrics).

RESULTS

Our results showed significant associations between all selected SNPs and schizophrenia (all P < 0.05). In patients, letter fluency and motor speed scores of T allele carriers were significantly higher than the CC group (all P < 0.05). Interestingly, greater brain activity, functional connectivity, and betweenness centrality (BC) in language processing and motor coordination were also observed in the corresponding brain zones in patients with the T allele based on a two-way ANCOVA model. Moreover, a potential positive correlation was found between brain activity/connectivity of these brain regions and verbal fluency and motor speed.

CONCLUSION

Our findings suggest that the DOCK4 gene may contribute to the onset of schizophrenia and lead to language processing and motor coordination dysfunction in this patient population from China.

The contribution of copy number variants to psychiatric symptoms and cognitive ability

Copy number variants (CNVs) are deletions and duplications of DNA sequence. The most frequently studied CNVs, which are described in this review, are recurrent CNVs that occur in the same locations on the genome. These CNVs have been strongly implicated in neurodevelopmental disorders, namely autism spectrum disorder (ASD), intellectual disability (ID), and developmental delay (DD), but also in schizophrenia. More recent work has also shown that CNVs increase risk for other psychiatric disorders, namely, depression, bipolar disorder, and post-traumatic stress disorder. Many of the same CNVs are implicated across all of these disorders, and these neuropsychiatric CNVs are also associated with cognitive ability in the general population, as well as with structural and functional brain alterations. Neuropsychiatric CNVs also show incomplete penetrance, such that carriers do not always develop any psychiatric disorder, and may show only mild symptoms, if any. Variable expressivity, whereby the same CNVs are associated with many different phenotypes of varied severity, also points to highly complex mechanisms underlying disease risk in CNV carriers. Comprehensive and longitudinal phenotyping studies of individual CNVs have provided initial insights into these mechanisms. However, more work is needed to estimate and predict the effect of non-recurrent, ultra-rare CNVs, which also contribute to psychiatric and cognitive outcomes. Moreover, delineating the broader phenotypic landscape of neuropsychiatric CNVs in both clinical and general population cohorts may also offer important mechanistic insights.

Gene and schizophrenia in the pregenome and postgenome-wide association studies era: a bibliometric analysis and network visualization

This study aimed to perform a bibliometric analysis on genetic studies in schizophrenia in the pregenome-wide association studies (GWAS) and post-GWAS era. We searched the literature on genes and schizophrenia using the Scopus database. The documents increased with time, especially after the human genome project and International HapMap Project, with the highest citation in 2008. The top occurrence author keywords were discovered to be different in the pre-GWAS and post-GWAS eras, reflecting the progress of genetic studies connected to schizophrenia. Emerging keywords highlighted a trend towards an application of precision medicine, showing an interplay of environmental exposures as well as genetic factors in schizophrenia pathogenesis, progression, and response to therapy. In conclusion, the gene and schizophrenia literature has grown rapidly after the human genome project, and the temporal variation in the author keywords pattern reflects the trend of genetic studies related to schizophrenia in the pre-GWAS and post-GWAS era.

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.

Genomic analysis of firework fear and noise reactivity in standard poodles

BACKGROUND

Fear of firework noises and other loud, sudden noises (noise reactivity) is a significant problem for many dogs and may have a negative effect on both welfare and, in severe cases, the life expectancy of dogs. A wide range of behavior traits, including fear-related behaviors, have high heritability estimates in dogs. The aim of this study was to estimate genomic heritability for fear of fireworks and loud noises in dogs.

RESULTS

A genomic heritability estimate was performed based on genome-wide SNPs from standard poodles with records of fear of fireworks and noise reactivity. The study was based on questionnaires answered by owners, who also volunteered to return a cheek swab from their dog for DNA analyses. SNP-based heritability was estimated to be 0.28 for firework fear and 0.16 for noise reactivity. We also identified an interesting region on chromosome 17 that was weakly associated with both traits.

CONCLUSIONS

We have estimated low to medium genomic heritabilities for fear of fireworks and noise reactivity in standard poodles. We have also identified an interesting region on chromosome 17, which harbors genes that have been shown to be involved in different psychiatric traits with anxiety components in humans. The region was associated with both traits; however, the association was weak and need further verification from other studies.

Inflammation and cognition in severe mental illness: patterns of covariation and subgroups

A potential relationship between dysregulation of immune/inflammatory pathways and cognitive impairment has been suggested in severe mental illnesses (SMI), such as schizophrenia (SZ) and bipolar (BD) spectrum disorders. However, multivariate relationships between peripheral inflammatory/immune-related markers and cognitive domains are unclear, and many studies do not account for inter-individual variance in both cognitive functioning and inflammatory/immune status. This study aimed to investigate covariance patterns between inflammatory/immune-related markers and cognitive domains and further elucidate heterogeneity in a large SMI and healthy control (HC) cohort (SZ = 343, BD = 289, HC = 770). We applied canonical correlation analysis (CCA) to identify modes of maximum covariation between a comprehensive selection of cognitive domains and inflammatory/immune markers. We found that poor verbal learning and psychomotor processing speed was associated with higher levels of interleukin-18 system cytokines and beta defensin 2, reflecting enhanced activation of innate immunity, a pattern augmented in SMI compared to HC. Applying hierarchical clustering on covariance patterns identified by the CCA revealed a high cognition-low immune dysregulation subgroup with predominantly HC (24% SZ, 45% BD, 74% HC) and a low cognition-high immune dysregulation subgroup predominantly consisting of SMI patients (76% SZ, 55% BD, 26% HC). These subgroups differed in IQ, years of education, age, CRP, BMI (all groups), level of functioning, symptoms and defined daily dose (DDD) of antipsychotics (SMI cohort). Our findings suggest a link between cognitive impairment and innate immune dysregulation in a subset of individuals with severe mental illness.

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.

Genetics of bipolar disorder: insights into its complex architecture and biology from common and rare variants

Bipolar disorder (BD) is a common mental disorder characterized by recurrent mood episodes, which causes major socioeconomic burdens globally. Though its disease pathogenesis is largely unknown, the high heritability of BD indicates strong contributions from genetic factors. In this review, we summarize the recent achievements in the genetics of BD, particularly those from genome-wide association study (GWAS) of common variants and next-generation sequencing analysis of rare variants. These include the identification of dozens of robust disease-associated loci, deepening of our understanding of the biology of BD, objective description of correlations with other psychiatric disorders and behavioral traits, formulation of methods for predicting disease risk and drug response, and the discovery of a single gene associated with bipolar disorder and schizophrenia spectrum with a large effect size. On the other hand, the findings to date have not yet made a clear contribution to the improvement of clinical psychiatry of BD. We overview the remaining challenges as well as possible paths to resolve them, referring to studies of other major neuropsychiatric disorders.

Independent Associated SNPs at SORCS3 and Its Protein Interactors for Multiple Brain-Related Disorders and Traits

Sortilin-related vacuolar protein sorting 10 (VPS10) domain containing receptor 3 (SORCS3) is a neuron-specific transmembrane protein involved in the trafficking of proteins between intracellular vesicles and the plasma membrane. Genetic variation at SORCS3 is associated with multiple neuropsychiatric disorders and behavioural phenotypes. Here, we undertake a systematic search of published genome-wide association studies to identify and catalogue associations between SORCS3 and brain-related disorders and traits. We also generate a SORCS3 gene-set based on protein-protein interactions and investigate the contribution of this gene-set to the heritability of these phenotypes and its overlap with synaptic biology. Analysis of association signals at SORSC3 showed individual SNPs to be associated with multiple neuropsychiatric and neurodevelopmental brain-related disorders and traits that have an impact on the experience of feeling, emotion or mood or cognitive function, while multiple LD-independent SNPs were associated with the same phenotypes. Across these SNPs, alleles associated with the more favourable outcomes for each phenotype (e.g., decreased risk of neuropsychiatric illness) were associated with increased expression of the SORCS3 gene. The SORCS3 gene-set was enriched for heritability contributing to schizophrenia (SCZ), bipolar disorder (BPD), intelligence (IQ) and education attainment (EA). Eleven genes from the SORCS3 gene-set were associated with more than one of these phenotypes at the genome-wide level, with RBFOX1 associated with SCZ, IQ and EA. Functional annotation revealed that the SORCS3 gene-set is enriched for multiple ontologies related to the structure and function of synapses. Overall, we find many independent association signals at SORCS3 with brain-related disorders and traits, with the effect possibly mediated by reduced gene expression, resulting in a negative impact on synaptic function.

Impact of traumatic life events and polygenic risk scores for major depression and posttraumatic stress disorder on Iraq/Afghanistan Veterans

Traumatic experiences and genetic heritability are among the most widely acknowledged risk factors leading to the development of psychopathology; including posttraumatic stress disorder (PTSD) and major depressive disorder (MDD). The purpose of this study was to investigate if polygenic risk scores (PRS) among Veterans interacted with traumatic stress to predict PTSD and MDD. 1,389 Iraq-Afghanistan military service Veterans from the Mental Illness Research Education and Clinical Center dataset were analyzed. Genome-wide association study (GWAS) statistics were utilized to generate PRS for PTSD (PRS_PTSD) and PRS for MDD (PRS_MDD) in order to analyze PRS-by-environment (PRSxE) with trauma exposure to predict PTSD and MDD diagnoses. Trauma exposure and PRS_PTSD, were independently associated with a current PTSD diagnosis (p < 0.001 and p < 0.001, respectively). The interaction between trauma exposure and PRS_MDD to predict a current diagnosis of PTSD trended towards significance (p = 0.053). Stratifying by trauma thresholds, among those within the lowest trauma load, the association of PRS_MDD with PTSD was found to be nominally significant (p = 0.03). For a MDD diagnosis, there was a significant association with trauma exposure (p < 0.001); and the association with PRS_MDD was found to be nominally significant (p = 0.03). No significant PRSxE effects were found with MDD. Our findings corroborate previous research highlighting trauma exposure, and genetic heritability, as risk factors for the development of PTSD and MDD in a Veteran population. Additionally, findings suggest that genetic vulnerability may be less important as trauma exposure increases, with high levels of trauma likely to result in PTSD and MDD, regardless of genetic vulnerability.

Depression trajectories and cytokines in schizophrenia spectrum disorders - A longitudinal observational study

Depression occurs frequently in all phases of schizophrenia spectrum disorders. Altered activity in the immune system is seen in both depression and schizophrenia. We aimed to uncover depressive trajectories in a sample of 144 adult individuals with schizophrenia spectrum disorders followed for one year, in order to identify possible cytokine profile differences. Patients were assessed longitudinally with the Positive and Negative Syndrome Scale (PANSS) and the Calgary Depression Scale for Schizophrenia (CDSS), where a score above 6 predicts depression. The serum cytokine concentrations for tumor necrosis factor (TNF)-alpha, interferon (IFN)-gamma, interleukin (IL)-1beta, IL-2, IL-4, IL-6, IL-10, IL-12p70 and IL-17A were measured using immunoassays. Latent growth curve models, multilevel models and latent class growth analysis (LCGA) were applied. The LCGA model supported three latent classes (trajectories) with differing CDSS profiles during the one-year follow-up: a high CDSS group (40.8 % of participants), a moderate CDSS group (43.9 %) and a low CDSS group (15.3 %). Five single PANSS items predicted affiliation to depressive trajectory: hallucinations, difficulty in abstract thinking, anxiety, guilt feelings and tension. In the high CDSS group, despite diminishing psychotic symptoms, depressive symptoms persisted throughout one year. The pro-inflammatory cytokines IFN-γ, IL-1β and TNF-α were differentially distributed between the depressive trajectories, although levels remained remarkably stable throughout 12 months. Significant changes were found for the anti-inflammatory cytokine IL-10 at baseline with an accompanying difference in change over time. More research is required to optimize future treatment stratification and investigate the contribution of inflammation in depressed patients with schizophrenia spectrum disorders.

Genomic insights into schizophrenia

Schizophrenia is a common, complex, heterogeneous psychiatric syndrome which can have profound impacts on affected individuals and imposes significant burdens on society. Despite intensive research, it has been challenging to understand basic mechanisms and to identify novel therapeutic targets. Given its high heritability and the complexity and inaccessibility of the human brain, much hope has been invested in the application of genomics as a route to better understanding. This work has identified many common and rare risk alleles and laid the foundations for a new generation of mechanistic studies. Genomics has also thrown new light on the relationship between schizophrenia and other psychiatric disorders and revealed its previously unappreciated aetiological relationship with childhood neurodevelopmental disorders, providing further evidence that it has its origins in disturbances of brain development. In addition, genomic findings suggest that the condition reflects fundamental disturbances in neuronal, and particularly synaptic, function that impact broadly on brain function, rather than being a disorder of specific brain regions and circuits. Finally, genomics has provided a plausible solution to the evolutionary paradox of how the condition persists in the face of high heritability and reduced fecundity.