Regulatory Variant rs2535629 in ITIH3 Intron Confers Schizophrenia Risk By Regulating CTCF Binding and SFMBT1 Expression

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.

The contribution of silencer variants to human diseases

BACKGROUND

Although disease-causal genetic variants have been found within silencer sequences, we still lack a comprehensive analysis of the association of silencers with diseases. Here, we profiled GWAS variants in 2.8 million candidate silencers across 97 human samples derived from a diverse panel of tissues and developmental time points, using deep learning models.

RESULTS

We show that candidate silencers exhibit strong enrichment in disease-associated variants, and several diseases display a much stronger association with silencer variants than enhancer variants. Close to 52% of candidate silencers cluster, forming silencer-rich loci, and, in the loci of Parkinson's-disease-hallmark genes TRIM31 and MAL, the associated SNPs densely populate clustered candidate silencers rather than enhancers displaying an overall twofold enrichment in silencers versus enhancers. The disruption of apoptosis in neuronal cells is associated with both schizophrenia and bipolar disorder and can largely be attributed to variants within candidate silencers. Our model permits a mechanistic explanation of causative SNP effects by identifying altered binding of tissue-specific repressors and activators, validated with a 70% of directional concordance using SNP-SELEX. Narrowing the focus of the analysis to individual silencer variants, experimental data confirms the role of the rs62055708 SNP in Parkinson's disease, rs2535629 in schizophrenia, and rs6207121 in type 1 diabetes.

CONCLUSIONS

In summary, our results indicate that advances in deep learning models for the discovery of disease-causal variants within candidate silencers effectively "double" the number of functionally characterized GWAS variants. This provides a basis for explaining mechanisms of action and designing novel diagnostics and therapeutics.

Characteristics of changes in plasma proteome profiling after sleeve gastrectomy

Bariatric surgery (BS), recognized as the most effective intervention for morbid obesity and associated metabolic comorbidities, encompasses both weight loss-dependent and weight loss-independent mechanisms to exert its metabolic benefits. In this study, we employed plasma proteomics technology, a recently developed mass spectrometric approach, to quantitatively assess 632 circulating proteins in a longitudinal cohort of 9 individuals who underwent sleeve gastrectomy (SG). Through time series clustering and Gene Ontology (GO) enrichment analysis, we observed that complement activation, proteolysis, and negative regulation of triglyceride catabolic process were the primary biological processes enriched in down-regulated proteins. Conversely, up-regulated differentially expressed proteins (DEPs) were significantly associated with negative regulation of peptidase activity, fibrinolysis, keratinocyte migration, and acute-phase response. Notably, we identified seven proteins (ApoD, BCHE, CNDP1, AFM, ITIH3, SERPINF1, FCN3) that demonstrated significant alterations at 1-, 3-, and 6-month intervals post SG, compared to baseline. These proteins play essential roles in metabolism, immune and inflammatory responses, as well as oxidative stress. Consequently, they hold promising potential as therapeutic targets for combating obesity and its associated comorbidities.

Genome-wide meta-analysis, functional genomics and integrative analyses implicate new risk genes and therapeutic targets for anxiety disorders

Anxiety disorders are the most prevalent mental disorders. However, the genetic etiology of anxiety disorders remains largely unknown. Here we conducted a genome-wide meta-analysis on anxiety disorders by including 74,973 (28,392 proxy) cases and 400,243 (146,771 proxy) controls. We identified 14 risk loci, including 10 new associations near CNTNAP5, MAP2, RAB9BP1, BTN1A1, PRR16, PCLO, PTPRD, FARP1, CDH2 and RAB27B. Functional genomics and fine-mapping pinpointed the potential causal variants, and expression quantitative trait loci analysis revealed the potential target genes regulated by the risk variants. Integrative analyses, including transcriptome-wide association study, proteome-wide association study and colocalization analyses, prioritized potential causal genes (including CTNND1 and RAB27B). Evidence from multiple analyses revealed possibly causal genes, including RAB27B, BTN3A2, PCLO and CTNND1. Finally, we showed that Ctnnd1 knockdown affected dendritic spine density and resulted in anxiety-like behaviours in mice, revealing the potential role of CTNND1 in anxiety disorders. Our study identified new risk loci, potential causal variants and genes for anxiety disorders, providing insights into the genetic architecture of anxiety disorders and potential therapeutic targets.

A Functional Schizophrenia-associated genetic variant near the TSNARE1 and ADGRB1 genes

Recent collaborative genome wide association studies (GWAS) have identified >200 independent loci contributing to risk for schizophrenia (SCZ). The genes closest to these loci have diverse functions, supporting the potential involvement of multiple relevant biological processes; yet there is no direct evidence that individual variants are functional or directly linked to specific genes. Nevertheless, overlap with certain epigenetic marks suggest that most GWAS-implicated variants are regulatory. Based on the strength of association with SCZ and the presence of regulatory epigenetic marks, we chose one such variant near TSNARE1 and ADGRB1, rs4129585, to test for functional potential and assay differences that may drive the pathogenicity of the risk allele. We observed that the variant-containing sequence drives reporter expression in relevant neuronal populations in zebrafish. Next, we introduced each allele into human induced pluripotent cells and differentiated 4 isogenic clones homozygous for the risk allele and 5 clones homozygous for the non-risk allele into neural precursor cells. Employing RNA-seq, we found that the two alleles yield significant transcriptional differences in the expression of 109 genes at FDR <0.05 and 259 genes at FDR <0.1. We demonstrate that these genes are highly interconnected in pathways enriched for synaptic proteins, axon guidance, and regulation of synapse assembly. Exploration of genes near rs4129585 suggests that this variant does not regulate TSNARE1 transcripts, as previously thought, but may regulate the neighboring ADGRB1, a regulator of synaptogenesis. Our results suggest that rs4129585 is a functional common variant that functions in specific pathways likely involved in SCZ risk.

An updated catalog of CTCF variants associated with neurodevelopmental disorder phenotypes

Introduction

CTCF-related disorder (CRD) is a neurodevelopmental disorder (NDD) caused by monoallelic pathogenic variants in CTCF. The first CTCF variants in CRD cases were documented in 2013. To date, 76 CTCF variants have been further described in the literature. In recent years, due to the increased application of next-generation sequencing (NGS), growing numbers of CTCF variants are being identified, and multiple genotype-phenotype databases cataloging such variants are emerging.

Methods

In this study, we aimed to expand the genotypic spectrum of CRD, by cataloging NDD phenotypes associated with reported CTCF variants. Here, we systematically reviewed all known CTCF variants reported in case studies and large-scale exome sequencing cohorts. We also conducted a meta-analysis using public variant data from genotype-phenotype databases to identify additional CTCF variants, which we then curated and annotated.

Results

From this combined approach, we report an additional 86 CTCF variants associated with NDD phenotypes that have not yet been described in the literature. Furthermore, we describe and explain inconsistencies in the quality of reported variants, which impairs the reuse of data for research of NDDs and other pathologies.

Discussion

From this integrated analysis, we provide a comprehensive and annotated catalog of all currently known CTCF mutations associated with NDD phenotypes, to aid diagnostic applications, as well as translational and basic research.

3D genome alterations and editing in pathology

The human genome is folded into a multi-level 3D structure that controls many nuclear functions including gene expression. Recently, alterations in 3D genome organization were associated with several genetic diseases and cancer. As a consequence, experimental approaches are now being developed to modify the global 3D genome organization and that of specific loci. Here, we discuss emerging experimental approaches of 3D genome editing that may prove useful in biomedicine.