A missense variant in NDUFA6 confers schizophrenia risk by affecting YY1 binding and NAGA expression

rs2736098, a synonymous polymorphism, is associated with carcinogenesis and cell count in multiple tissue types by regulating TERT expression

rs2736098 is a synonymous polymorphism in TERT (telomerase reverse transcriptase), an enzyme involved in tumor onset of multiple tissues, and should play no roles in carcinogenesis. However, a search in cancer somatic mutation database indicated that the mutation frequency at rs2736098 is much higher than the average one for TERT. Moreover, there are significant H3K4me1 and H3K27Ac signals, two universal histone modifications for active enhancers, surrounding rs2736098. Therefore, we hypothesized that rs2736098 might be within an enhancer region, regulate TERT expression and influence cancer risk. Through luciferase assay, it was verified that the enhancer activity of rs2736098C allele is significantly higher than that of T in multiple tissues. Transfection of plasmids containing TERT coding region with two different alleles indicated that rs2736098C allele can induce a significantly higher TERT expression than T. By chromatin immunoprecipitation, it was observed that the fragment spanning rs2736098 can interact with USF1 (upstream transcription factor 1). The two alleles of rs2736098 present evidently different binding affinity with nuclear proteins. Database and literature search indicated that rs2736098 is significantly associated with carcinogenesis in multiple tissues and count of multiple cell types. All these facts indicated that rs2736098 is also an oncogenic polymorphism and plays important role in cell proliferation.

Validation of Enhancer Regions in Primary Human Neural Progenitor Cells using Capture STARR-seq

Genome-wide association studies (GWAS) and expression analyses implicate noncoding regulatory regions as harboring risk factors for psychiatric disease, but functional characterization of these regions remains limited. We performed capture STARR-sequencing of over 78,000 candidate regions to identify active enhancers in primary human neural progenitor cells (phNPCs). We selected candidate regions by integrating data from NPCs, prefrontal cortex, developmental timepoints, and GWAS. Over 8,000 regions demonstrated enhancer activity in the phNPCs, and we linked these regions to over 2,200 predicted target genes. These genes are involved in neuronal and psychiatric disease-associated pathways, including dopaminergic synapse, axon guidance, and schizophrenia. We functionally validated a subset of these enhancers using mutation STARR-sequencing and CRISPR deletions, demonstrating the effects of genetic variation on enhancer activity and enhancer deletion on gene expression. Overall, we identified thousands of highly active enhancers and functionally validated a subset of these enhancers, improving our understanding of regulatory networks underlying brain function and disease.

Functional genomics identify causal variant underlying the protective CTSH locus for Alzheimer's disease

Alzheimer's disease (AD) is the most prevalent age-related neurodegenerative disease, which has a high heritability of up to 79%. Exploring the genetic basis is essential for understanding the pathogenic mechanisms underlying AD development. Recent genome-wide association studies (GWASs) reported an AD-associated signal in the Cathepsin H (CTSH) gene in European populations. However, the exact functional/causal variant(s), and the genetic regulating mechanism of CTSH in AD remain to be determined. In this study, we carried out a comprehensive study to characterize the role of CTSH variants in the pathogenesis of AD. We identified rs2289702 in CTSH as the most significant functional variant that is associated with a protective effect against AD. The genetic association between rs2289702 and AD was validated in independent cohorts of the Han Chinese population. The CTSH mRNA expression level was significantly increased in AD patients and AD animal models, and the protective allele T of rs2289702 was associated with a decreased expression level of CTSH through the disruption of the binding affinity of transcription factors. Human microglia cells with CTSH knockout showed a significantly increased phagocytosis of Aβ peptides. Our study identified CTSH as being involved in AD genetic susceptibility and uncovered the genetic regulating mechanism of CTSH in pathogenesis of AD.

Identifying Novel Risk Conferring Genes Involved in Glycosylation Processes with Familial Schizophrenia in an Indian Cohort: Prediction of ADAMTS9 gene Variant for Structural Stability

Schizophrenia is a complex neuropsychiatric disorder and heritability is as high as 80% making it the most heritable mental disorder. Although GWAS has identified numerous variants, the pathophysiology is still elusive. Here, an attempt was made to identify genetic risk factors in familial cases of schizophrenia that are associated with a common causative pathway. To achieve this objective, exome sequencing was done in 4 familial cases and identified six unique coding variants in five genes. Among these genes, PIGQ gene has two pathogenic variants, one nonsense and in-frame deletion. One missense variant in GALNT16 and one in GALNT5 have variable damaging score, however, the other variants, in ADAMTS9 and in LTBP4 have the highest damaging score. Further analysis showed that the variant of LTBP4 was not present in the functional domain. The other missense variant in the ADAMTS9 gene was found to be significant and was present in the thrombospondin repeat motif, one of the important motifs. Detailed molecular dynamics simulation study on this variant showed a damaging effect on structural stability. Since, all these genes culminated into the glycosylation process, it was evident that an aberrant glycosylation process may be one of the risk factors. Although, extracellular matrix formation through glycosylation have been shown to be associated, the involvement of ADAMTS9 and PIGQ gene mediated glycosylation has not been reported. In this paper, a novel link between ADAMTS9 and PIGQ gene with schizophrenia have been reported. Therefore, this novel observation has contributed immensely to the existing knowledge on risk factor of Schizophrenia.

Increased parietal and occipital lobe gyrification predicts conversion to syndromal psychosis in a clinical high-risk cohort

BACKGROUND

Local gyrification index (lGI), indicative of the degree of cortical folding is a proxy marker for early cortical neurodevelopmental abnormalities. We studied the difference in lGI between those who do and do not convert to psychosis (non-converters) in a clinical high-risk (CHR) cohort, and whether lGI predicts conversion to psychosis.

METHODS

Seventy-two CHR participants with attenuated positive symptom syndrome were followed up for two years. The difference in baseline whole-brain lGI was examined on the T1-weighted MRIs between, i)CHR (N = 72) and healthy controls (N = 19), ii)Converters to psychosis (N = 24) and non-converters (N = 48), adjusting for age and sex, on Freesurfer-6.0. The significant cluster obtained in the converters versus non-converters comparison was registered as a region of interest to individual images of all 72 participants and lGI values were extracted from this region. A cox proportional hazards model was applied with these values to study whether lGI predicts conversion to psychosis.

RESULTS

lGI was not different between CHR and healthy controls. lGI was increased in converters in the right-sided inferior parietal and lateral occipital areas (corrected cluster-wise-p-value = 0.009, cohen's f = 0.42) compared to non-converters, which significantly increased the risk of onset of psychosis (p = 0.029, hazard ratio = 1.471).

CONCLUSIONS

Increased gyrification in the right-sided inferior parietal and lateral occipital area differentiates converters to psychosis in CHR, significantly increasing the risk of conversion to psychosis. This measure may reflect underlying traits in parts of the brain that develop earliest in-utero (parietal and occipital), conferring a heightened vulnerability to convert to syndromal psychosis subsequently.

ERVWE1 Reduces Hippocampal Neuron Density and Impairs Dendritic Spine Morphology through Inhibiting Wnt/JNK Non-Canonical Pathway via miR-141-3p in Schizophrenia

Human endogenous retroviruses (HERVs) are remnants of ancestral germline infections by exogenous retroviruses. Human endogenous retroviruses W family envelope gene (HERV-W env, also called ERVWE1), located on chromosome 7q21-22, encodes an envelope glycoprotein from the HERV-W family. Mounting evidence suggests that aberrant expression of ERVWE1 involves the etiology of schizophrenia. Moreover, the genetic and morphological studies indicate that dendritic spine deficits may contribute to the onset of schizophrenia. Here, we reported that ERVWE1 changed the density and morphology of the dendritic spine through inhibiting Wingless-type (Wnt)/c-Jun N-terminal kinases (JNK) non-canonical pathway via miR-141-3p in schizophrenia. In this paper, we found elevated levels of miR-141-3p and a significant positive correlation with ERVWE1 in schizophrenia. Moreover, serum Wnt5a and actin-related protein 2 (Arp2) levels decreased and demonstrated a significant negative correlation with ERVWE1 in schizophrenia. In vitro experiments disclosed that ERVWE1 up-regulated miR-141-3p expression by interacting with transcription factor (TF) Yin Yang 1 (YY1). YY1 modulated miR-141-3p expression by binding to its promoter. The luciferase assay revealed that YY1 enhanced the promoter activity of miR-141-3p. Using the miRNA target prediction databases and luciferase reporter assays, we demonstrated that miR-141-3p targeted Wnt5a at its 3' untranslated region (3' UTR). Furthermore, ERVWE1 suppressed the expression of Arp2 through non-canonical pathway, Wnt5a/JNK signaling pathway. In addition, ERVWE1 inhibited Wnt5a/JNK/Arp2 signal pathway through miR-141-3p. Finally, functional assays showed that ERVWE1 induced the abnormalities in hippocampal neuron morphology and spine density through inhibiting Wnt/JNK non-canonical pathway via miR-141-3p in schizophrenia. Our findings indicated that miR-141-3p, Wnt5a, and Arp2 might be potential clinical blood-based biomarkers or therapeutic targets for schizophrenia. Our work also provided new insight into the role of ERVWE1 in schizophrenia pathogenesis.

The schizophrenia-associated missense variant rs13107325 regulates dendritic spine density

The missense variant rs13107325 (C/T, p.Ala391Thr) in SLC39A8 consistently showed robust association with schizophrenia in recent genome-wide association studies (GWASs), suggesting the potential pathogenicity of this non-synonymous risk variant. Nevertheless, how this missense variant confers schizophrenia risk remains unknown. Here we constructed a knock-in mouse model (by introducing a threonine at the 393th amino acid of mouse SLC39A8 (SLC39A8-p.393T), which corresponds to rs13107325 (p.Ala391Thr) of human SLC39A8) to explore the potential roles and biological effects of this missense variant in schizophrenia pathogenesis. We assessed multiple phenotypes and traits (associated with rs13107325) of the knock-in mice, including body and brain weight, concentrations of metal ions (including cadmium, zinc, manganese, and iron) transported by SLC39A8, blood lipids, proliferation and migration of neural stem cells (NSCs), cortical development, behaviors and cognition, transcriptome, dendritic spine density, and synaptic transmission. Many of the tested phenotypes did not show differences in SLC39A8-p.393T knock-in and wild-type mice. However, we found that zinc concentration in brain and blood of SLC39A8-p.393T knock-in mice was dysregulated compared with wild-types, validating the functionality of rs13107325. Further analysis indicated that cortical dendritic spine density of the SLC39A8-p.393T knock-in mice was significantly decreased compared with wild-types, indicating the important role of SLC39A8-p.393T in dendritic spine morphogenesis. These results indicated that SLC39A8-p.393T knock-in resulted in decreased dendritic spine density, thus mimicking the dendritic spine pathology observed in schizophrenia. Our study indicates that rs13107325 might confer schizophrenia risk by regulating zinc concentration and dendritic spine density, a featured characteristic that was frequently reported to be decreased in schizophrenia.

Functional variant rs2270363 on 16p13.3 confers schizophrenia risk by regulating NMRAL1

Recent genome-wide association studies have reported multiple schizophrenia risk loci, yet the functional variants and their roles in schizophrenia remain to be characterized. Here we identify a functional single nucleotide polymorphism (rs2270363: G>A) at the schizophrenia risk locus 16p13.3. rs2270363 lies in the E-box element of the promoter of NMRAL1 and disrupts binding of the basic helix-loop-helix leucine zipper family proteins, including USF1, MAX and MXI1. We validated the regulatory effects of rs2270363 using reporter gene assays and electrophoretic mobility shift assay. Besides, expression quantitative trait loci analysis showed that the risk allele (A) of rs2270363 was significantly associated with elevated NMRAL1 expression in the human brain. Transcription factors knockdown and CRISPR-Cas9-mediated editing further confirmed the regulatory effects of the genomic region containing rs2270363 on NMRAL1. Intriguingly, NMRAL1 was significantly downregulated in the brain of schizophrenia patients compared with healthy subjects, and knockdown of Nmral1 expression affected proliferation and differentiation of mouse neural stem cells, as well as genes and pathways associated with brain development and synaptic transmission. Of note, Nmral1 knockdown resulted in significant decrease of dendritic spine density, revealing the potential pathophysiological mechanisms of NMRAL1 in schizophrenia. Finally, we independently confirmed the association between rs2270363 and schizophrenia in the Chinese population and found that the risk allele of rs2270363 was the same in European and Chinese populations. These lines of evidence suggest that rs2270363 may confer schizophrenia risk by regulating NMRAL1, a gene whose expression dysregulation might be involved in the pathogenesis of schizophrenia by affecting neurodevelopment and synaptic plasticity.

Functional genomic analysis delineates regulatory mechanisms of GWAS-identified bipolar disorder risk variants

Background

Genome-wide association studies (GWASs) have identified multiple risk loci for bipolar disorder (BD). However, pinpointing functional (or causal) variants in the reported risk loci and elucidating their regulatory mechanisms remain challenging.

Methods

We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) data from human brain tissues (or neuronal cell lines) and position weight matrix (PWM) data to identify functional single-nucleotide polymorphisms (SNPs). Then, we verified the regulatory effects of these transcription factor (TF) binding–disrupting SNPs (hereafter referred to as “functional SNPs”) through a series of experiments, including reporter gene assays, allele-specific expression (ASE) analysis, TF knockdown, CRISPR/Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis. Finally, we overexpressed PACS1 (whose expression was most significantly associated with the identified functional SNPs rs10896081 and rs3862386) in mouse primary cortical neurons to investigate if PACS1 affects dendritic spine density.

Results

We identified 16 functional SNPs (in 9 risk loci); these functional SNPs disrupted the binding of 7 TFs, for example, CTCF and REST binding was frequently disrupted. We then identified the potential target genes whose expression in the human brain was regulated by these functional SNPs through eQTL analysis. Of note, we showed dysregulation of some target genes of the identified TF binding–disrupting SNPs in BD patients compared with controls, and overexpression of PACS1 reduced the density of dendritic spines, revealing the possible biological mechanisms of these functional SNPs in BD.

Conclusions

Our study identifies functional SNPs in some reported risk loci and sheds light on the regulatory mechanisms of BD risk variants. Further functional characterization and mechanistic studies of these functional SNPs and candidate genes will help to elucidate BD pathogenesis and develop new therapeutic approaches and drugs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01057-3.

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

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

The schizophrenia-associated variant in SLC39A8 alters protein glycosylation in the mouse brain

A missense mutation (A391T) in SLC39A8 is strongly associated with schizophrenia in genomic studies, though the molecular connection to the brain is unknown. Human carriers of A391T have reduced serum manganese, altered plasma glycosylation, and brain MRI changes consistent with altered metal transport. Here, using a knock-in mouse model homozygous for A391T, we show that the schizophrenia-associated variant changes protein glycosylation in the brain. Glycosylation of Asn residues in glycoproteins (N-glycosylation) was most significantly impaired, with effects differing between regions. RNAseq analysis showed negligible regional variation, consistent with changes in the activity of glycosylation enzymes rather than gene expression. Finally, nearly one-third of detected glycoproteins were differentially N-glycosylated in the cortex, including members of several pathways previously implicated in schizophrenia, such as cell adhesion molecules and neurotransmitter receptors that are expressed across all cell types. These findings provide a mechanistic link between a risk allele and potentially reversible biochemical changes in the brain, furthering our molecular understanding of the pathophysiology of schizophrenia and a novel opportunity for therapeutic development.

Functional genomics elucidates regulatory mechanisms of Parkinson's disease-associated variants

BACKGROUND

Genome-wide association studies (GWASs) have identified multiple risk loci for Parkinson's disease (PD). However, identifying the functional (or potential causal) variants in the reported risk loci and elucidating their roles in PD pathogenesis remain major challenges. To identify the potential causal (or functional) variants in the reported PD risk loci and to elucidate their regulatory mechanisms, we report a functional genomics study of PD.

METHODS

We first integrated chromatin immunoprecipitation sequencing (ChIP-Seq) (from neuronal cells and human brain tissues) data and GWAS-identified single-nucleotide polymorphisms (SNPs) in PD risk loci. We then conducted a series of experiments and analyses to validate the regulatory effects of these (i.e., functional) SNPs, including reporter gene assays, allele-specific expression (ASE), transcription factor (TF) knockdown, CRISPR-Cas9-mediated genome editing, and expression quantitative trait loci (eQTL) analysis.

RESULTS

We identified 44 SNPs (from 11 risk loci) affecting the binding of 12 TFs and we validated the regulatory effects of 15 TF binding-disrupting SNPs. In addition, we also identified the potential target genes regulated by these TF binding-disrupting SNPs through eQTL analysis. Finally, we showed that 4 eQTL genes of these TF binding-disrupting SNPs were dysregulated in PD cases compared with controls.

CONCLUSION

Our study systematically reveals the gene regulatory mechanisms of PD risk variants (including widespread disruption of CTCF binding), generates the landscape of potential PD causal variants, and pinpoints promising candidate genes for further functional characterization and drug development.

Transcriptomics Reveals the Effect of Thymol on the Growth and Toxin Production of Fusarium graminearum

Fusarium graminearum is a harmful pathogen causing head blight in cereals such as wheat and barley, and thymol has been proven to inhibit the growth of many pathogens. This study aims to explore the fungistatic effect of thymol on F. graminearum and its mechanism. Different concentrations of thymol were used to treat F. graminearum. The results showed that the EC50 concentration of thymol against F. graminearum was 40 μg/mL. Compared with the control group, 40 μg/mL of thymol reduced the production of Deoxynivalenol (DON) and 3-Ac-DON by 70.1% and 78.2%, respectively. Our results indicate that thymol can effectively inhibit the growth and toxin production of F. graminearum and cause an extensive transcriptome response. Transcriptome identified 16,727 non-redundant unigenes and 1653 unigenes that COG did not annotate. The correlation coefficients between samples were all >0.941. When FC was 2.0 times, a total of 3230 differential unigenes were identified, of which 1223 were up-regulated, and 2007 were down-regulated. Through the transcriptome, we confirmed that the expression of many genes involved in F. graminearum growth and synthesis of DON and other secondary metabolites were also changed. The gluconeogenesis/glycolysis pathway may be a potential and important way for thymol to affect the growth of F. graminearum hyphae and the production of DON simultaneously.

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

Genome-wide association studies have identified 3p21.1 as a robust risk locus for schizophrenia. However, the underlying molecular mechanisms remain elusive. Here a functional regulatory variant (rs2535629) is identified that disrupts CTCF binding at 3p21.1. It is confirmed that rs2535629 is also significantly associated with schizophrenia in Chinese population and the regulatory effect of rs2535629 is validated. Expression quantitative trait loci analysis indicates that rs2535629 is associated with the expression of three distal genes (GLT8D1, SFMBT1, and NEK4) in the human brain, and CRISPR-Cas9-mediated genome editing confirmed the regulatory effect of rs2535629 on GLT8D1, SFMBT1, and NEK4. Interestingly, differential expression analysis of GLT8D1, SFMBT1, and NEK4 suggested that rs2535629 may confer schizophrenia risk by regulating SFMBT1 expression. It is further demonstrated that Sfmbt1 regulates neurodevelopment and dendritic spine density, two key pathological characteristics of schizophrenia. Transcriptome analysis also support the potential role of Sfmbt1 in schizophrenia pathogenesis. The study identifies rs2535629 as a plausibly causal regulatory variant at the 3p21.1 risk locus and demonstrates the regulatory mechanism and biological effect of this functional variant, indicating that this functional variant confers schizophrenia risk by altering CTCF binding and regulating expression of SFMBT1, a distal gene which plays important roles in neurodevelopment and synaptic morphogenesis.

Proteome-wide Association Study Provides Insights Into the Genetic Component of Protein Abundance in Psychiatric Disorders

BACKGROUND

Genome-wide association studies have identified multiple risk variants for psychiatric disorders. Nevertheless, how the risk variants confer risk of psychiatric disorders remains largely unknown.

METHODS

We performed proteome-wide association studies to identify genes whose cis-regulated protein abundance change in the human brain were associated with psychiatric disorders.

RESULTS

By integrating genome-wide associations of four common psychiatric disorders and two independent brain proteomes (n = 376 and n = 152, respectively) from the dorsolateral prefrontal cortex, we identified 61 genes (including 48 genes for schizophrenia, 12 genes for bipolar disorder, 5 genes for depression, and 2 genes for attention-deficit/hyperactivity disorder) whose genetically regulated protein abundance levels were associated with risk of psychiatric disorders. Comparison with transcriptome-wide association studies identified 18 overlapping genes that showed significant associations with psychiatric disorders at both proteome-wide and transcriptome-wide levels, suggesting that genetic risk variants likely confer risk of psychiatric disorders by regulating messenger RNA expression and protein abundance of these genes.

CONCLUSIONS

Our study not only provides new insights into the genetic component of protein abundance in psychiatric disorders but also highlights several high-confidence risk proteins (including CNNM2 and CTNND1) for schizophrenia and depression. These high-confidence risk proteins represent promising therapeutic targets for future drug development.

Chromatin Alterations in Neurological Disorders and Strategies of (Epi)Genome Rescue

Neurological disorders (NDs) comprise a heterogeneous group of conditions that affect the function of the nervous system. Often incurable, NDs have profound and detrimental consequences on the affected individuals' lives. NDs have complex etiologies but commonly feature altered gene expression and dysfunctions of the essential chromatin-modifying factors. Hence, compounds that target DNA and histone modification pathways, the so-called epidrugs, constitute promising tools to treat NDs. Yet, targeting the entire epigenome might reveal insufficient to modify a chosen gene expression or even unnecessary and detrimental to the patients' health. New technologies hold a promise to expand the clinical toolkit in the fight against NDs. (Epi)genome engineering using designer nucleases, including CRISPR-Cas9 and TALENs, can potentially help restore the correct gene expression patterns by targeting a defined gene or pathway, both genetically and epigenetically, with minimal off-target activity. Here, we review the implication of epigenetic machinery in NDs. We outline syndromes caused by mutations in chromatin-modifying enzymes and discuss the functional consequences of mutations in regulatory DNA in NDs. We review the approaches that allow modifying the (epi)genome, including tools based on TALENs and CRISPR-Cas9 technologies, and we highlight how these new strategies could potentially change clinical practices in the treatment of NDs.