Targeted resequencing of regulatory regions at schizophrenia risk loci: Role of rare functional variants at chromatin repressive states

Sex-specific transcriptional and proteomic signatures in schizophrenia

It has remained unclear why schizophrenia typically manifests after adolescence and which neurobiological mechanisms are underlying the cascade leading to the actual onset of the illness. Here we show that the use of induced pluripotent stem cell-derived neurons of monozygotic twins from pairs discordant for schizophrenia enhances disease-specific signal by minimizing genetic heterogeneity. In proteomic and pathway analyses, clinical illness is associated especially with altered glycosaminoglycan, GABAergic synapse, sialylation, and purine metabolism pathways. Although only 12% of all 19,462 genes are expressed differentially between healthy males and females, up to 61% of the illness-related genes are sex specific. These results on sex-specific genes are replicated in another dataset. This implies that the pathophysiology differs between males and females, and may explain why symptoms appear after adolescence when the expression of many sex-specific genes change, and suggests the need for sex-specific treatments.

Identification of the primate-specific gene BTN3A2 as an additional schizophrenia risk gene in the MHC loci

BACKGROUND

Schizophrenia is a complex mental disorder resulting in poor life quality and high social and economic burden. Despite the fact that genome-wide association studies (GWASs) have successfully identified a number of risk loci for schizophrenia, identifying the causal genes at the risk loci and elucidating their roles in disease pathogenesis remain major challenges.

METHODS

The summary data-based Mendelian randomization analysis (SMR) was used to integrate a large-scale GWAS of schizophrenia with brain expression quantitative trait loci (eQTL) data and brain methylation expression quantitative trait loci (meQTL) data, to identify novel risk gene(s) for schizophrenia. We then analyzed the mRNA expression and methylation statuses of the gene hit BTN3A2 during the early brain development. Electrophysiological analyses of CA1 pyramidal neurons were performed to evaluate the excitatory and inhibitory synaptic activity after overexpression of BTN3A2 in rat hippocampal slices. Cell surface binding assay was used to test the interaction of BTN3A2 and neurexins.

FINDINGS

We identified BTN3A2 as a potential risk gene for schizophrenia. The mRNA expression and methylation data showed that BTN3A2 expression in human brain is highest post-natally. Further electrophysiological analyses of rat hippocampal slices showed that BTN3A2 overexpression specifically suppressed the excitatory synaptic activity onto CA1 pyramidal neurons, most likely through its interaction with the presynaptic adhesion molecule neurexins.

INTERPRETATION

Increased expression of BTN3A2 might confer risk for schizophrenia by altering excitatory synaptic function. Our result constitutes a paradigm for distilling risk gene using an integrative analysis and functional characterization in the post-GWAS era. FUND: This study was supported by the Strategic Priority Research Program (B) of the Chinese Academy of Sciences (XDB02020003 to Y-GY), the National Natural Science Foundation of China (31730037 to Y-GY), and the Bureau of Frontier Sciences and Education, Chinese Academy of Sciences (QYZDJ-SSW-SMC005 to Y-GY).

Integrating Imaging Genomic Data in the Quest for Biomarkers of Schizophrenia Disease

It's increasingly important but difficult to determine potential biomarkers of schizophrenia (SCZ) disease, owing to the complex pathophysiology of this disease. In this study, a network-fusion based framework was proposed to identify genetic biomarkers of the SCZ disease. A three-step feature selection was applied to single nucleotide polymorphisms (SNPs), DNA methylation, and functional magnetic resonance imaging (fMRI) data to select important features, which were then used to construct two gene networks in different states for the SNPs and DNA methylation data, respectively. Two health networks (one is for SNP data and the other is for DNA methylation data) were combined into one health network from which health minimum spanning trees (MSTs) were extracted. Two disease networks also followed the same procedures. Those genes with significant changes were determined as SCZ biomarkers by comparing MSTs in two different states and they were finally validated from five aspects. The effectiveness of the proposed discovery framework was also demonstrated by comparing with other network-based discovery methods. In summary, our approach provides a general framework for discovering gene biomarkers of the complex diseases by integrating imaging genomic data, which can be applied to the diagnosis of the complex diseases in the future.

EZH1 is an antipsychotic-sensitive epigenetic modulator of social and motivational behavior that is dysregulated in schizophrenia

BACKGROUND

With the capacity to modulate gene networks in an environmentally-sensitive manner, the role of epigenetic systems in mental disorders has come under intense investigation. Dysregulation of epigenetic effectors, including microRNAs and histone-modifying enzymes, may better explain the role of environmental risk factors and the observed heritability rate that cannot be fully attributed to known genetic risk alleles. Here, we aimed to identify novel epigenetic targets of the schizophrenia-associated microRNA 132 (miR-132).

METHODS

Histone modifications were quantified by immunodetection in response to viral-mediated overexpression of miR-132 while a luminescent reporter system was used to validate targets of miR-132 in vitro. Genome-wide profiling, quantitative PCR and NanoSting were used to quantify gene expression in post-mortem human brains, neuronal cultures and prefrontal cortex (PFC) of mice chronically exposed to antipsychotics. Following viral-mediated depletion of Enhancer of Zeste 1 (EZH1) in the murine PFC, behaviors including sociability and motivation were assessed using a 3-chambered apparatus and forced-swim test, respectively.

RESULTS

Overexpression of miR-132 decreased global histone 3 lysine 27 tri-methylation (H3K27me3), a repressive epigenetic mark. Moreover, the polycomb-associated H3K27 methyltransferase, EZH1, is regulated by miR-132 and upregulated in the PFC of schizophrenics. Unlike its homolog EZH2, expression of EZH1 in the murine PFC decreased following chronic exposure to antipsychotics. Viral-mediated depletion of EZH1 in the mouse PFC attenuated sociability, enhanced motivational behaviors, and affected gene expression pathways related to neurotransmission and behavioral phenotypes.

CONCLUSIONS

EZH1 is dysregulated in schizophrenia, sensitive to antipsychotic medications, and a brain-enriched miR-132 target that controls neurobehavioral phenotypes.

Investigation of common, low-frequency and rare genome-wide variation in anorexia nervosa

Anorexia nervosa (AN) is a complex neuropsychiatric disorder presenting with dangerously low body weight, and a deep and persistent fear of gaining weight. To date, only one genome-wide significant locus associated with AN has been identified. We performed an exome-chip based genome-wide association studies (GWAS) in 2158 cases from nine populations of European origin and 15 485 ancestrally matched controls. Unlike previous studies, this GWAS also probed association in low-frequency and rare variants. Sixteen independent variants were taken forward for in silico and de novo replication (11 common and 5 rare). No findings reached genome-wide significance. Two notable common variants were identified: rs10791286, an intronic variant in OPCML (P=9.89 × 10-6), and rs7700147, an intergenic variant (P=2.93 × 10-5). No low-frequency variant associations were identified at genome-wide significance, although the study was well-powered to detect low-frequency variants with large effect sizes, suggesting that there may be no AN loci in this genomic search space with large effect sizes.

Characterization of the human thyroid epigenome

The thyroid gland, necessary for normal human growth and development, functions as an essential regulator of metabolism by the production and secretion of appropriate levels of thyroid hormone. However, assessment of abnormal thyroid function may be challenging suggesting a more fundamental understanding of normal function is needed. One way to characterize normal gland function is to study the epigenome and resulting transcriptome within its constituent cells. This study generates the first published reference epigenomes for human thyroid from four individuals using ChIP-seq and RNA-seq. We profiled six histone modifications (H3K4me1, H3K4me3, H3K27ac, H3K36me3, H3K9me3, H3K27me3), identified chromatin states using a hidden Markov model, produced a novel quantitative metric for model selection and established epigenomic maps of 19 chromatin states. We found that epigenetic features characterizing promoters and transcription elongation tend to be more consistent than regions characterizing enhancers or Polycomb-repressed regions and that epigenetically active genes consistent across all epigenomes tend to have higher expression than those not marked as epigenetically active in all epigenomes. We also identified a set of 18 genes epigenetically active and consistently expressed in the thyroid that are likely highly relevant to thyroid function. Altogether, these epigenomes represent a powerful resource to develop a deeper understanding of the underlying molecular biology of thyroid function and provide contextual information of thyroid and human epigenomic data for comparison and integration into future studies.