Expression quantitative trait loci in the developing human brain and their enrichment in neuropsychiatric disorders

Towards understanding sex differences in autism spectrum disorders

Autism Spectrum Disorder (ASD) is a neurodevelopmental condition characterized by social deficits, repetitive behaviours and lack of empathy. Its significant genetic heritability and potential comorbidities often lead to diagnostic and therapeutic challenges. This review addresses the biological basis of ASD, focusing on the sex differences in gene expression and hormonal influences. ASD is more commonly diagnosed in males at a ratio of 4:1, indicating a potential oversight in female-specific ASD research and a risk of underdiagnosis in females. We consider how ASD manifests differently across sexes by exploring differential gene expression in female and male brains and consider how variations in steroid hormones influence ASD characteristics. Synaptic function, including excitation/inhibition ratio imbalance, is influenced by gene mutations and this is explored as a key factor in the cognitive and behavioural manifestations of ASD. We also discuss the role of micro RNAs (miRNAs) and highlight a novel mutation in miRNA-873, which affects a suite of key synaptic genes, neurexin, neuroligin, SHANK and post-synaptic density proteins, implicated in the pathology of ASD. Our review suggests that genetic predisposition, sex differences in brain gene expression, and hormonal factors significantly contribute to the presentation, identification and severity of ASD, necessitating sex-specific considerations in diagnosis and treatments. These findings advocate for personalized interventions to improve the outcomes for individuals with ASD.

Cross-ancestry atlas of gene, isoform, and splicing regulation in the developing human brain

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

Deciphering the role of a SINE-VNTR-Alu retrotransposon polymorphism as a biomarker of Parkinson's disease progression

SINE-VNTR-Alu (SVA) retrotransposons are transposable elements which represent a source of genetic variation. We previously demonstrated that the presence/absence of a human-specific SVA, termed SVA_67, correlated with the progression of Parkinson's disease (PD). In the present study, we demonstrate that SVA_67 acts as expression quantitative trait loci, thereby exhibiting a strong regulatory effect across the genome using whole genome and transcriptomic data from the Parkinson's progression markers initiative cohort. We further show that SVA_67 is polymorphic for its variable number tandem repeat domain which correlates with both regulatory properties in a luciferase reporter gene assay in vitro and differential expression of multiple genes in vivo. Additionally, this variation's utility as a biomarker is reflected in a correlation with a number of PD progression markers. These experiments highlight the plethora of transcriptomic and phenotypic changes associated with SVA_67 polymorphism which should be considered when investigating the missing heritability of neurodegenerative diseases.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia

Multimorbidity represents an increasingly important public health challenge with far-reaching implications for health management and policy. Mental health and metabolic diseases have a well-established epidemiological association. In this study, we investigate the genetic intersection between type 2 diabetes and schizophrenia. We use Mendelian randomization to examine potential causal relationships between the two conditions and related endophenotypes. We report no compelling evidence that type 2 diabetes genetic liability potentially causally influences schizophrenia risk and vice versa. Our findings show that increased body mass index (BMI) has a protective effect against schizophrenia, in contrast to the well-known risk-increasing effect of BMI on type 2 diabetes risk. We identify evidence of colocalization of association signals for these two conditions at 11 genomic loci, six of which have opposing directions of effect for type 2 diabetes and schizophrenia. To elucidate these colocalizing signals, we integrate multi-omics data from bulk and single-cell gene expression studies, along with functional information. We identify putative effector genes and find that they are enriched for homeostasis and lipid-related pathways. We also highlight drug repurposing opportunities including N-methyl-D-aspartate (NMDA) receptor antagonists. Our findings provide insights into shared biological mechanisms for type 2 diabetes and schizophrenia, highlighting common factors that influence the risk of the two conditions in opposite directions and shedding light on the complex nature of this comorbidity.

Decreased CNNM2 expression in prefrontal cortex affects sensorimotor gating function, cognition, dendritic spine morphogenesis and risk of schizophrenia

Genome-wide association studies (GWASs) have reported multiple single nucleotide polymorphisms (SNPs) associated with schizophrenia, yet the underlying molecular mechanisms are largely unknown. In this study, we aimed to identify schizophrenia relevant genes showing alterations in mRNA and protein expression associated with risk SNPs at the 10q24.32-33 GWAS locus. We carried out the quantitative trait loci (QTL) and summary data-based Mendelian randomization (SMR) analyses, using the PsychENCODE dorsolateral prefrontal cortex (DLPFC) expression QTL (eQTL) database, as well as the ROSMAP and Banner DLPFC protein QTL (pQTL) datasets. The gene CNNM2 (encoding a magnesium transporter) at 10q24.32-33 was identified to be a robust schizophrenia risk gene, and was highly expressed in human neurons according to single cell RNA-seq (scRNA-seq) data. We further revealed that reduced Cnnm2 in the mPFC of mice led to impaired cognition and compromised sensorimotor gating function, and decreased Cnnm2 in primary cortical neurons altered dendritic spine morphogenesis, confirming the link between CNNM2 and endophenotypes of schizophrenia. Proteomics analyses showed that reduced Cnnm2 level changed expression of proteins associated with neuronal structure and function. Together, these results identify a robust gene in the pathogenesis of schizophrenia.

Overview: Research on the Genetic Architecture of the Developing Cerebral Cortex in Norms and Diseases

The human brain is characterized by high cell numbers, diverse cell types with diverse functions, and intricate connectivity with an exceedingly broad surface of the cortex. Human-specific brain development was accomplished by a long timeline for maturation from the prenatal period to the third decade of life. The long timeline makes complicated architecture and circuits of human cerebral cortex possible, and it makes human brain vulnerable to intrinsic and extrinsic insults resulting in the development of variety of neuropsychiatric disorders. Unraveling the molecular and cellular processes underlying human brain development under the elaborate regulation of gene expression in a spatiotemporally specific manner, especially that of the cortex will provide a biological understanding of human cognition and behavior in health and diseases. Global research consortia and the advancing technologies in brain science including functional genomics equipped with emergent neuroinformatics such as single-cell multiomics, novel human models, and high-volume databases with high-throughput computation facilitate the biological understanding of the development of the human brain cortex. Knowing the process of interplay of the genome and the environment in cortex development will lead us to understand the human-specific cognitive function and its individual diversity. Thus, it is worthwhile to overview the recent progress in neurotechnology to foresee further understanding of the human brain and norms and diseases.

Genetic implication of prenatal GABAergic and cholinergic neuron development in susceptibility to schizophrenia

Background

The ganglionic eminences are fetal-specific structures that give rise to gamma-aminobutyric acid (GABA)- and acetylcholine- releasing neurons of the forebrain. Given evidence for GABAergic and cholinergic disturbances in schizophrenia, as well as an early neurodevelopmental component to the disorder, we tested the potential involvement of developing cells of the ganglionic eminences in mediating genetic risk for the condition.

Study Design

We combined data from a recent large-scale genome-wide association study of schizophrenia with single cell RNA sequencing data from the human ganglionic eminences to test enrichment of schizophrenia risk variation in genes with high expression specificity for particular developing cell populations within these structures. We additionally performed the single nuclei Assay for Transposase-Accessible Chromatin with Sequencing (snATAC-Seq) to map potential regulatory genomic regions operating in individual cell populations of the human ganglionic eminences, using these to additionally test for enrichment of schizophrenia common genetic variant liability and to functionally annotate non-coding variants associated with the disorder.

Study Results

Schizophrenia common variant liability was enriched in genes with high expression specificity for developing neuron populations that are predicted to form dopamine D1 and D2 receptor expressing GABAergic medium spiny neurons of the striatum, cortical somatostatin-positive GABAergic interneurons, calretinin-positive GABAergic neurons and cholinergic neurons. Consistent with these findings, schizophrenia genetic risk was also concentrated in predicted regulatory genomic sequence mapped in developing neuronal populations of the ganglionic eminences.

Conclusions

Our study provides evidence for a role of prenatal GABAergic and cholinergic neuron development in later susceptibility to schizophrenia.

The genetic relationships between brain structure and schizophrenia

Genetic risks for schizophrenia are theoretically mediated by genetic effects on brain structure but it has been unclear which genes are associated with both schizophrenia and cortical phenotypes. We accessed genome-wide association studies (GWAS) of schizophrenia (N = 69,369 cases; 236,642 controls), and of three magnetic resonance imaging (MRI) metrics (surface area, cortical thickness, neurite density index) measured at 180 cortical areas (N = 36,843, UK Biobank). Using Hi-C-coupled MAGMA, 61 genes were significantly associated with both schizophrenia and one or more MRI metrics. Whole genome analysis with partial least squares demonstrated significant genetic covariation between schizophrenia and area or thickness of most cortical regions. Genetic similarity between cortical areas was strongly coupled to their phenotypic covariance, and genetic covariation between schizophrenia and brain phenotypes was strongest in the hubs of structural covariance networks. Pleiotropically associated genes were enriched for neurodevelopmental processes and positionally concentrated in chromosomes 3p21, 17q21 and 11p11. Mendelian randomization analysis indicated that genetically determined variation in a posterior cingulate cortical area could be causal for schizophrenia. Parallel analyses of GWAS on bipolar disorder, Alzheimer's disease and height showed that pleiotropic association with MRI metrics was stronger for schizophrenia compared to other disorders.

Trans -eQTL hotspots shape complex traits by modulating cellular states

Regulatory genetic variation shapes gene expression, providing an important mechanism connecting DNA variation and complex traits. The causal relationships between gene expression and complex traits remain poorly understood. Here, we integrated transcriptomes and 46 genetically complex growth traits in a large cross between two strains of the yeast Saccharomyces cerevisiae . We discovered thousands of genetic correlations between gene expression and growth, suggesting functional connections. Local regulatory variation was a minor source of these genetic correlations. Instead, genetic correlations tended to arise from multiple independent trans -acting regulatory loci. Trans -acting hotspots that affect the expression of numerous genes accounted for particularly large fractions of genetic growth variation and of genetic correlations between gene expression and growth. Genes with genetic correlations were enriched for similar biological processes across traits, but with heterogeneous direction of effect. Our results reveal how trans -acting regulatory hotspots shape complex traits by altering cellular states.

RET enhancer haplotype-dependent remodeling of the human fetal gut development program

Hirschsprung disease (HSCR) is associated with deficiency of the receptor tyrosine kinase RET, resulting in loss of cells of the enteric nervous system (ENS) during fetal gut development. The major contribution to HSCR risk is from common sequence variants in RET enhancers with additional risk from rare coding variants in many genes. Here, we demonstrate that these RET enhancer variants specifically alter the human fetal gut development program through significant decreases in gene expression of RET, members of the RET-EDNRB gene regulatory network (GRN), other HSCR genes, with an altered transcriptome of 2,382 differentially expressed genes across diverse neuronal and mesenchymal functions. A parsimonious hypothesis for these results is that beyond RET's direct effect on its GRN, it also has a major role in enteric neural crest-derived cell (ENCDC) precursor proliferation, its deficiency reducing ENCDCs with relative expansion of non-ENCDC cells. Thus, genes reducing RET proliferative activity can potentially cause HSCR. One such class is the 23 RET-dependent transcription factors enriched in early gut development. We show that their knockdown in human neuroblastoma SK-N-SH cells reduces RET and/or EDNRB gene expression, expanding the RET-EDNRB GRN. The human embryos we studied had major remodeling of the gut transcriptome but were unlikely to have had HSCR: thus, genetic or epigenetic changes in addition to those in RET are required for aganglionosis.

Genomic insights into the comorbidity between type 2 diabetes and schizophrenia

Multimorbidity represents an increasingly important public health challenge with far-reaching implications for health management and policy. Mental health and metabolic diseases have a well-established epidemiological association. In this study, we investigate the genetic intersection between type 2 diabetes and schizophrenia. We use Mendelian randomization to examine potential causal relationships between the two conditions and related endophenotypes. We report no compelling evidence that type 2 diabetes genetic liability potentially causally influences schizophrenia risk and vice versa. Our findings show that increased body mass index (BMI) has a protective effect against schizophrenia, in contrast to the well-known risk-increasing effect of BMI on type 2 diabetes risk. We identify evidence of colocalization of association signals for these two conditions at 11 genomic loci, six of which have opposing directions of effect for type 2 diabetes and schizophrenia. To elucidate these colocalizing signals, we integrate multi-omics data from bulk and single-cell gene expression studies, along with functional information. We identify high-confidence effector genes and find that they are enriched for homeostasis and lipid-related pathways. We also highlight drug repurposing opportunities including N-methyl-D-aspartate (NMDA) receptor antagonists. Our findings provide insights into shared biological mechanisms for type 2 diabetes and schizophrenia, highlighting common factors that influence the risk of the two conditions in opposite directions and shedding light on the complex nature of this comorbidity.

Mapping microRNA expression quantitative trait loci in the prenatal human brain implicates miR-1908-5p expression in bipolar disorder and other brain-related traits

MicroRNA (miRNA) are small non-coding RNA involved in post-transcriptional gene regulation. Given their known involvement in early neurodevelopment processes, we here sought to identify common genetic variants associated with altered miRNA expression in the prenatal human brain. We performed small RNA sequencing on brain tissue from 112 genome-wide genotyped fetuses from the second trimester of gestation, identifying high-confidence (false discovery rate < 0.05) expression quantitative trait loci for 30 mature miRNA. Integrating our findings with genome-wide association study data for brain-related disorders, we implicate increased prenatal expression of miR-1908-5p as a risk mechanism for bipolar disorder and find that predicted mRNA targets of miR-1908-5p that are expressed in the fetal brain are enriched for common variant genetic association with the condition. Extending these analyses to other brain-related traits, we find that common genetic variation associated with increased miR-1908-5p expression in fetal brain is additionally associated with depressive symptoms, irritability, increased right cerebellum exterior volume and increased sleep duration in the general population. Our findings provide support to the view that altered miRNA expression can influence susceptibility to neuropsychiatric illness and suggest an early neurodevelopmental risk mechanism for bipolar disorder.

SLCO5A1 and synaptic assembly genes contribute to impulsivity in juvenile myoclonic epilepsy

Elevated impulsivity is a key component of attention-deficit hyperactivity disorder (ADHD), bipolar disorder and juvenile myoclonic epilepsy (JME). We performed a genome-wide association, colocalization, polygenic risk score, and pathway analysis of impulsivity in JME (n = 381). Results were followed up with functional characterisation using a drosophila model. We identified genome-wide associated SNPs at 8q13.3 (P = 7.5 × 10-9) and 10p11.21 (P = 3.6 × 10-8). The 8q13.3 locus colocalizes with SLCO5A1 expression quantitative trait loci in cerebral cortex (P = 9.5 × 10-3). SLCO5A1 codes for an organic anion transporter and upregulates synapse assembly/organisation genes. Pathway analysis demonstrates 12.7-fold enrichment for presynaptic membrane assembly genes (P = 0.0005) and 14.3-fold enrichment for presynaptic organisation genes (P = 0.0005) including NLGN1 and PTPRD. RNAi knockdown of Oatp30B, the Drosophila polypeptide with the highest homology to SLCO5A1, causes over-reactive startling behaviour (P = 8.7 × 10-3) and increased seizure-like events (P = 6.8 × 10-7). Polygenic risk score for ADHD genetically correlates with impulsivity scores in JME (P = 1.60 × 10-3). SLCO5A1 loss-of-function represents an impulsivity and seizure mechanism. Synaptic assembly genes may inform the aetiology of impulsivity in health and disease.

Sex Significantly Impacts the Function of Major Depression-Linked Variants In Vivo

BACKGROUND

Genome-wide association studies have discovered blocks of common variants-likely transcriptional-regulatory-associated with major depressive disorder (MDD), though the functional subset and their biological impacts remain unknown. Likewise, why depression occurs in females more frequently than males is unclear. We therefore tested the hypothesis that risk-associated functional variants interact with sex and produce greater impact in female brains.

METHODS

We developed techniques to directly measure regulatory variant activity and sex interactions using massively parallel reporter assays in the mouse brain in vivo, in a cell type-specific manner, and applied these approaches to measure activity of >1000 variants from >30 MDD loci.

RESULTS

We identified extensive sex-by-allele effects in mature hippocampal neurons, suggesting that sex-differentiated impacts of genetic risk may underlie sex bias in disease. Unbiased informatics approaches indicated that functional MDD variants recurrently disrupt a number of transcription factor binding motifs, including those of sex hormone receptors. We confirmed a role for the latter by performing massively parallel reporter assays in neonatal mice on the day of birth (during a sex-differentiating hormone surge) and hormonally quiescent juveniles.

CONCLUSIONS

Our study provides novel insights into the influence of age, biological sex, and cell type on regulatory variant function and provides a framework for in vivo parallel assays to functionally define interactions between organismal variables such as sex and regulatory variation. Moreover, we experimentally demonstrate that a portion of the sex differences seen in MDD occurrence may be a product of sex-differentiated effects at associated regulatory variants.

Multiple genes in a single GWAS risk locus synergistically mediate aberrant synaptic development and function in human neurons

The mechanistic tie between genome-wide association study (GWAS)-implicated risk variants and disease-relevant cellular phenotypes remains largely unknown. Here, using human induced pluripotent stem cell (hiPSC)-derived neurons as a neurodevelopmental model, we identify multiple schizophrenia (SZ) risk variants that display allele-specific open chromatin (ASoC) and are likely to be functional. Editing the strongest ASoC SNP, rs2027349, near vacuolar protein sorting 45 homolog (VPS45) alters the expression of VPS45, lncRNA AC244033.2, and a distal gene, C1orf54. Notably, the transcriptomic changes in neurons are associated with SZ and other neuropsychiatric disorders. Neurons carrying the risk allele exhibit increased dendritic complexity and hyperactivity. Interestingly, individual/combinatorial gene knockdown shows that these genes alter cellular phenotypes in a non-additive synergistic manner. Our study reveals that multiple genes at a single GWAS risk locus mediate a compound effect on neural function, providing a mechanistic link between a non-coding risk variant and disease-related cellular phenotypes.

Genetic insights into human cortical organization and development through genome-wide analyses of 2,347 neuroimaging phenotypes

Our understanding of the genetics of the human cerebral cortex is limited both in terms of the diversity and the anatomical granularity of brain structural phenotypes. Here we conducted a genome-wide association meta-analysis of 13 structural and diffusion magnetic resonance imaging-derived cortical phenotypes, measured globally and at 180 bilaterally averaged regions in 36,663 individuals and identified 4,349 experiment-wide significant loci. These phenotypes include cortical thickness, surface area, gray matter volume, measures of folding, neurite density and water diffusion. We identified four genetic latent structures and causal relationships between surface area and some measures of cortical folding. These latent structures partly relate to different underlying gene expression trajectories during development and are enriched for different cell types. We also identified differential enrichment for neurodevelopmental and constrained genes and demonstrate that common genetic variants associated with cortical expansion are associated with cephalic disorders. Finally, we identified complex interphenotype and inter-regional genetic relationships among the 13 phenotypes, reflecting the developmental differences among them. Together, these analyses identify distinct genetic organizational principles of the cortex and their correlates with neurodevelopment.

Relationship between sex biases in gene expression and sex biases in autism and Alzheimer's disease

Sex differences in the brain may play an important role in sex-differential prevalence of neuropsychiatric conditions. In order to understand the transcriptional basis of sex differences, we analyzed multiple, large-scale, human postmortem brain RNA-seq datasets using both within-region and pan-regional frameworks. We find evidence of sex-biased transcription in many autosomal genes, some of which provide evidence for pathways and cell population differences between chromosomally male and female individuals. These analyses also highlight regional differences in the extent of sex-differential gene expression. We observe an increase in specific neuronal transcripts in male brains and an increase in immune and glial function-related transcripts in female brains. Integration with single-cell data suggests this corresponds to sex differences in cellular states rather than cell abundance. Integration with case-control gene expression studies suggests a female molecular predisposition towards Alzheimer's disease, a female-biased disease. Autism, a male-biased diagnosis, does not exhibit a male predisposition pattern in our analysis. Finally, we provide region specific analyses of sex differences in brain gene expression to enable additional studies at the interface of gene expression and diagnostic differences.

Chromosomal inversion polymorphisms shape human brain morphology

The impact of chromosomal inversions on human brain morphology remains underexplored. We studied 35 common inversions classified from genotypes of 33,018 adults with European ancestry. The inversions at 2p22.3, 16p11.2, and 17q21.31 reach genome-wide significance, followed by 8p23.1 and 6p21.33, in their association with cortical and subcortical morphology. The 17q21.31, 8p23.1, and 16p11.2 regions comprise the LRRC37, OR7E, and NPIP duplicated gene families. We find the 17q21.31 MAPT inversion region, known for harboring neurological risk, to be the most salient locus among common variants for shaping and patterning the cortex. Overall, we observe the inverted orientations decreasing brain size, with the exception that the 2p22.3 inversion is associated with increased subcortical volume and the 8p23.1 inversion is associated with increased motor cortex. These significant inversions are in the genomic hotspots of neuropsychiatric loci. Our findings are generalizable to 3,472 children and demonstrate inversions as essential genetic variation to understand human brain phenotypes.

Analysis of subcellular RNA fractions demonstrates significant genetic regulation of gene expression in human brain post-transcriptionally

Gaining insight into the genetic regulation of gene expression in human brain is key to the interpretation of genome-wide association studies for major neurological and neuropsychiatric diseases. Expression quantitative trait loci (eQTL) analyses have largely been used to achieve this, providing valuable insights into the genetic regulation of steady-state RNA in human brain, but not distinguishing between molecular processes regulating transcription and stability. RNA quantification within cellular fractions can disentangle these processes in cell types and tissues which are challenging to model in vitro. We investigated the underlying molecular processes driving the genetic regulation of gene expression specific to a cellular fraction using allele-specific expression (ASE). Applying ASE analysis to genomic and transcriptomic data from paired nuclear and cytoplasmic fractions of anterior prefrontal cortex, cerebellar cortex and putamen tissues from 4 post-mortem neuropathologically-confirmed control human brains, we demonstrate that a significant proportion of genetic regulation of gene expression occurs post-transcriptionally in the cytoplasm, with genes undergoing this form of regulation more likely to be synaptic. These findings have implications for understanding the structure of gene expression regulation in human brain, and importantly the interpretation of rapidly growing single-nucleus brain RNA-sequencing and eQTL datasets, where cytoplasm-specific regulatory events could be missed.

The relationship between case-control differential gene expression from brain tissue and genetic associations in schizophrenia

Large numbers of genetic loci have been identified that are known to contain common risk alleles for schizophrenia, but linking associated alleles to specific risk genes remains challenging. Given that most alleles that influence liability to schizophrenia are thought to do so by altered gene expression, intuitively, case-control differential gene expression studies should highlight genes with a higher probability of being associated with schizophrenia and could help identify the most likely causal genes within associated loci. Here, we test this hypothesis by comparing transcriptome analysis of the dorsolateral prefrontal cortex from 563 schizophrenia cases and 802 controls with genome-wide association study (GWAS) data from the third wave study of the Psychiatric Genomics Consortium. Genes differentially expressed in schizophrenia were not enriched for common allelic association statistics compared with other brain-expressed genes, nor were they enriched for genes within associated loci previously reported to be prioritized by genetic fine-mapping. Genes prioritized by Summary-based Mendelian Randomization were underexpressed in cases compared to other genes in the same GWAS loci. However, the overall strength and direction of expression change predicted by SMR were not related to that observed in the differential expression data. Overall, this study does not support the hypothesis that genes identified as differentially expressed from RNA sequencing of bulk brain tissue are enriched for those that show evidence for genetic associations. Such data have limited utility for prioritizing genes in currently associated loci in schizophrenia.

Cross-ancestry genome-wide association meta-analyses of hippocampal and subfield volumes

The hippocampus is critical for memory and cognition and neuropsychiatric disorders, and its subfields differ in architecture and function. Genome-wide association studies on hippocampal and subfield volumes are mainly conducted in European populations; however, other ancestral populations are under-represented. Here we conduct cross-ancestry genome-wide association meta-analyses in 65,791 individuals for hippocampal volume and 38,977 for subfield volumes, including 7,009 individuals of East Asian ancestry. We identify 339 variant-trait associations at P < 1.13 × 10^-9 for 44 hippocampal traits, including 23 new associations. Common genetic variants have similar effects on hippocampal traits across ancestries, although ancestry-specific associations exist. Cross-ancestry analysis improves the fine-mapping precision and the prediction performance of polygenic scores in under-represented populations. These genetic variants are enriched for Wnt signaling and neuron differentiation and affect cognition, emotion and neuropsychiatric disorders. These findings may provide insight into the genetic architectures of hippocampal and subfield volumes.

Hidden Genetic Regulation of Human Complex Traits via Brain Isoforms

Alternative splicing exists in most multi-exonic genes, and exploring these complex alternative splicing events and their resultant isoform expressions is essential. However, it has become conventional that RNA sequencing results have often been summarized into gene-level expression counts mainly due to the multiple ambiguous mapping of reads at highly similar regions. Transcript-level quantification and interpretation are often overlooked, and biological interpretations are often deduced based on combined transcript information at the gene level. Here, for the most variable tissue of alternative splicing, the brain, we estimate isoform expressions in 1,191 samples collected by the Genotype-Tissue Expression (GTEx) Consortium using a powerful method that we previously developed. We perform genome-wide association scans on the isoform ratios per gene and identify isoform-ratio quantitative trait loci (irQTL), which could not be detected by studying gene-level expressions alone. By analyzing the genetic architecture of the irQTL, we show that isoform ratios regulate educational attainment via multiple tissues including the frontal cortex (BA9), cortex, cervical spinal cord, and hippocampus. These tissues are also associated with different neuro-related traits, including Alzheimer's or dementia, mood swings, sleep duration, alcohol intake, intelligence, anxiety or depression, etc. Mendelian randomization (MR) analysis revealed 1,139 pairs of isoforms and neuro-related traits with plausible causal relationships, showing much stronger causal effects than on general diseases measured in the UK Biobank (UKB). Our results highlight essential transcript-level biomarkers in the human brain for neuro-related complex traits and diseases, which could be missed by merely investigating overall gene expressions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-023-00100-6.

Leisure Activities, Genetic Risk, and Frailty: Evidence from the Chinese Adults Aged 80 Years or Older

INTRODUCTION

About half of adults aged ≥80 years suffer from frailty. Exercise is considered effective in preventing frailty but may be inapplicable to adults aged ≥80 years due to physical limitations. As an alternative, we aimed to explore the association of leisure activities with frailty and identify potential interaction with established polygenic risk score (PRS) among adults aged ≥80 years.

METHODS

Analyses were performed in a prospective cohort study of 7,471 community-living older adults aged ≥80 years who were recruited between 2002 and 2014 from 23 provinces in China. Leisure activity was assessed using a seven-question leisure activity index and frailty was defined as a frailty index ≥0.25 using a validated 39-item health-related scale. The PRS was constructed using 59 single-nucleotide polymorphisms associated with frailty in a subsample of 2,541 older adults. Cox proportional hazards models were used to explore the associations of leisure activities, PRS with frailty.

RESULTS

The mean age of participants was 89.4 ± 6.6 years (range: 80-116). In total, 2,930 cases of frailty were identified during 42,216 person-years of follow-up. Each 1 unit increase in the leisure activity index was associated with 12% lower risk of frailty (hazard ratio: 0.88 [95% confidence interval, 0.85-0.91]). Participants with high genetic risk (PRS >2.47 × 10-4) suffered from 26% higher risk of frailty. Interaction between leisure activity and genetic risk was not observed.

CONCLUSION

Evidence is presented for the independent association of leisure activities and genetic risk with frailty. Engagement in leisure activities is suggested to be associated with lower risk of frailty across all levels of genetic risk among adults aged ≥80 years.

Multi-ancestry phenome-wide association of complement component 4 variation with psychiatric and brain phenotypes in youth

BACKGROUND

Increased expression of the complement component 4A (C4A) gene is associated with a greater lifetime risk of schizophrenia. In the brain, C4A is involved in synaptic pruning; yet, it remains unclear the extent to which upregulation of C4A alters brain development or is associated with the risk for psychotic symptoms in childhood. Here, we perform a multi-ancestry phenome-wide association study in 7789 children aged 9-12 years to examine the relationship between genetically regulated expression (GREx) of C4A, childhood brain structure, cognition, and psychiatric symptoms.

RESULTS

While C4A GREx is not related to childhood psychotic experiences, cognition, or global measures of brain structure, it is associated with a localized reduction in regional surface area (SA) of the entorhinal cortex. Furthermore, we show that reduced entorhinal cortex SA at 9-10 years predicts a greater number and severity of psychosis-like events at 1-year and 2-year follow-up time points. We also demonstrate that the effects of C4A on the entorhinal cortex are independent of genome-wide polygenic risk for schizophrenia.

CONCLUSIONS

Our results suggest neurodevelopmental effects of C4A on childhood medial temporal lobe structure, which may serve as a biomarker for schizophrenia risk prior to symptom onset.

Cross-ancestry, cell-type-informed atlas of gene, isoform, and splicing regulation in the developing human brain

Genomic regulatory elements active in the developing human brain are notably enriched in genetic risk for neuropsychiatric disorders, including autism spectrum disorder (ASD), schizophrenia, and bipolar disorder. However, prioritizing the specific risk genes and candidate molecular mechanisms underlying these genetic enrichments has been hindered by the lack of a single unified large-scale gene regulatory atlas of human brain development. Here, we uniformly process and systematically characterize gene, isoform, and splicing quantitative trait loci (xQTLs) in 672 fetal brain samples from unique subjects across multiple ancestral populations. We identify 15,752 genes harboring a significant xQTL and map 3,739 eQTLs to a specific cellular context. We observe a striking drop in gene expression and splicing heritability as the human brain develops. Isoform-level regulation, particularly in the second trimester, mediates the greatest proportion of heritability across multiple psychiatric GWAS, compared with eQTLs. Via colocalization and TWAS, we prioritize biological mechanisms for ~60% of GWAS loci across five neuropsychiatric disorders, nearly two-fold that observed in the adult brain. Finally, we build a comprehensive set of developmentally regulated gene and isoform co-expression networks capturing unique genetic enrichments across disorders. Together, this work provides a comprehensive view of genetic regulation across human brain development as well as the stage-and cell type-informed mechanistic underpinnings of neuropsychiatric disorders.

Investigating the shared genetic architecture and causal relationship between pain and neuropsychiatric disorders

Pain often occurs in parallel with neuropsychiatric disorders. However, the underlying mechanisms and potential causality have not been well studied. We collected the genome-wide association study (GWAS) summary statistics of 26 common pain and neuropsychiatric disorders with sample size ranging from 17,310 to 482,730 in European population. The genetic correlation between pair of pain and neuropsychiatric disorders, as well as the relevant cell types were investigated by linkage disequilibrium (LD) score regression analyses. Then, transcriptome-wide association study (TWAS) was applied to identify the potential shared genes by integrating the gene expression information and GWAS. In addition, Mendelian randomization (MR) analyses were conducted to infer the potential causality between pain and neuropsychiatric disorders. Among the 169 pairwise pain and neuropsychiatric disorders, 55 pairs showed positive correlations (median r_g = 0.43) and 9 pairs showed negative correlations (median r_g =  -0.31). Using MR analyses, 26 likely causal associations were identified, including that neuroticism and insomnia were risk factors for most of short-term pain, and multisite chronic pain was risk factor for neuroticism, insomnia, major depressive disorder and attention deficit/hyperactivity disorder, and vice versa. The signals of pain and neuropsychiatric disorders tended to be enriched in the functional regions of cell types from central nervous system (CNS). A total of 19 genes shared in at least one pain and neuropsychiatric disorder pair were identified by TWAS, including AMT, NCOA6, and UNC45A, which involved in glycine degradation, insulin secretion, and cell proliferation, respectively. Our findings provided the evidence of shared genetic structure, causality and potential shared pathogenic mechanisms between pain and neuropsychiatric disorders, and enhanced our understanding of the comorbidities of pain and neuropsychiatric disorders.

Dissecting the 22q13 region to explore the genetic and phenotypic diversity of patients with Phelan-McDermid syndrome

SHANK3-related Phelan-McDermid syndrome (PMS) is caused by a loss of the distal part of chromosome 22, including SHANK3, or by a pathological SHANK3 variant. There is an important genetic and phenotypic diversity among patients who can present with developmental delay, language impairments, autism, epilepsy, and other symptoms. SHANK3, encoding a synaptic scaffolding protein, is deleted in the majority of patients with PMS and is considered a major gene involved in the neurological impairments of the patients. However, differences in deletion size can influence clinical features, and in some rare cases, deletions at the 22q13 locus in individuals with SHANK3-unrelated PMS do not encompass SHANK3. These individuals with SHANK3-unrelated PMS still display a PMS-like phenotype. This suggests the participation of other 22q13 genes in the pathogenesis of PMS. Here, we review the biological function and potential implication in PMS symptoms of 110 genes located in the 22q13 region, focusing on 35 genes with evidence for association with neurodevelopmental disorders, including 13 genes for epilepsy and 11 genes for microcephaly and/or macrocephaly. Our review is restricted to the 22q13 region, but future large-scale studies using whole genome sequencing and deep-phenotyping are warranted to develop predictive models of clinical trajectories and to target specific medical and educational care for each individual with PMS.

De novo human brain enhancers created by single-nucleotide mutations

Advanced human cognition is attributed to increased neocortex size and complexity, but the underlying evolutionary and regulatory mechanisms are largely unknown. Using human and macaque embryonic neocortical H3K27ac data coupled with a deep learning model of enhancers, we identified ~4000 enhancer gains in humans, which, per our model, can often be attributed to single-nucleotide essential mutations. Our analyses suggest that functional gains in embryonic brain development are associated with de novo enhancers whose putative target genes exhibit increased expression in progenitor cells and interneurons and partake in critical neural developmental processes. Essential mutations alter enhancer activity through altered binding of key transcription factors (TFs) of embryonic neocortex, including ISL1, POU3F2, PITX1/2, and several SOX TFs, and are associated with central nervous system disorders. Overall, our results suggest that essential mutations lead to gain of embryonic neocortex enhancers, which orchestrate expression of genes involved in critical developmental processes associated with human cognition.

Single-Nuclei RNA Sequencing of 5 Regions of the Human Prenatal Brain Implicates Developing Neuron Populations in Genetic Risk for Schizophrenia

BACKGROUND

While a variety of evidence supports a prenatal component in schizophrenia, there are few data regarding the cell populations involved. We sought to identify cells of the human prenatal brain mediating genetic risk for schizophrenia by integrating cell-specific gene expression measures generated through single-nuclei RNA sequencing with recent large-scale genome-wide association study (GWAS) and exome sequencing data for the condition.

METHODS

Single-nuclei RNA sequencing was performed on 5 brain regions (frontal cortex, ganglionic eminence, hippocampus, thalamus, and cerebellum) from 3 fetuses from the second trimester of gestation. Enrichment of schizophrenia common variant genetic liability and rare damaging coding variation was assessed in relation to gene expression specificity within each identified cell population.

RESULTS

Common risk variants were prominently enriched within genes with high expression specificity for developing neuron populations within the frontal cortex, ganglionic eminence, and hippocampus. Enrichments were largely independent of genes expressed in neuronal populations of the adult brain that have been implicated in schizophrenia through the same methods. Genes containing an excess of rare damaging variants in schizophrenia had higher expression specificity for developing glutamatergic neurons of the frontal cortex and hippocampus that were also enriched for common variant liability.

CONCLUSIONS

We found evidence for a distinct contribution of prenatal neuronal development to genetic risk for schizophrenia, involving specific populations of developing neurons within the second-trimester fetal brain. Our study significantly advances the understanding of the neurodevelopmental origins of schizophrenia and provides a resource with which to investigate the prenatal antecedents of other psychiatric and neurologic disorders.

MicroRNA-eQTLs in the developing human neocortex link miR-4707-3p expression to brain size

Expression quantitative trait loci (eQTL) data have proven important for linking non-coding loci to protein-coding genes. But eQTL studies rarely measure microRNAs (miRNAs), small non-coding RNAs known to play a role in human brain development and neurogenesis. Here, we performed small-RNA sequencing across 212 mid-gestation human neocortical tissue samples, measured 907 expressed miRNAs, discovering 111 of which were novel, and identified 85 local-miRNA-eQTLs. Colocalization of miRNA-eQTLs with GWAS summary statistics yielded one robust colocalization of miR-4707-3p expression with educational attainment and brain size phenotypes, where the miRNA expression increasing allele was associated with decreased brain size. Exogenous expression of miR-4707-3p in primary human neural progenitor cells decreased expression of predicted targets and increased cell proliferation, indicating miR-4707-3p modulates progenitor gene regulation and cell fate decisions. Integrating miRNA-eQTLs with existing GWAS yielded evidence of a miRNA that may influence human brain size and function via modulation of neocortical brain development.

Bidirectional genetic overlap between bipolar disorder and intelligence

BACKGROUND

Bipolar disorder (BD) is a highly heritable psychiatric illness exhibiting substantial correlation with intelligence.

METHODS

To investigate the shared genetic signatures between BD and intelligence, we utilized the summary statistics from genome-wide association studies (GWAS) to conduct the bivariate causal mixture model (MiXeR) and conjunctional false discovery rate (conjFDR) analyses. Subsequent expression quantitative trait loci (eQTL) mapping in human brain and enrichment analyses were also performed.

RESULTS

Analysis with MiXeR suggested that approximately 10.3K variants could influence intelligence, among which 7.6K variants were correlated with the risk of BD (Dice: 0.80), and 47% of these variants predicted BD risk and intelligence in consistent allelic directions. The conjFDR analysis identified 37 distinct genomic loci that were jointly associated with BD and intelligence with a conjFDR < 0.01, and 16 loci (43%) had the same directions of allelic effects in both phenotypes. Brain eQTL analyses found that genes affected by the "concordant loci" were distinct from those modulated by the "discordant loci". Enrichment analyses suggested that genes related to the "concordant loci" were significantly enriched in pathways/phenotypes related with synapses and sleep quality, whereas genes associated with the "discordant loci" were enriched in pathways related to cell adhesion, calcium ion binding, and abnormal emotional phenotypes.

CONCLUSIONS

We confirmed the polygenic overlap with mixed directions of allelic effects between BD and intelligence and identified multiple genomic loci and risk genes. This study provides hints for the mesoscopic phenotypes of BD and relevant biological mechanisms, promoting the knowledge of the genetic and phenotypic heterogeneity of BD. The essential value of leveraging intelligence in BD investigations is also highlighted.

Keeping the balance: Trade-offs between human brain evolution, autism, and schizophrenia

The unique qualities of the human brain are a product of a complex evolutionary process. Evolution, famously described by François Jacob as a “tinkerer,” builds upon existing genetic elements by modifying and repurposing them for new functions. Genetic changes in DNA may lead to the emergence of new genes or cause altered gene expression patterns. Both gene and regulatory element mutations may lead to new functions. Yet, this process may lead to side-effects. An evolutionary trade-off occurs when an otherwise beneficial change, which is important for evolutionary success and is under strong positive selection, concurrently results in a detrimental change in another trait. Pleiotropy occurs when a gene affects multiple traits. Antagonistic pleiotropy is a phenomenon whereby a genetic variant leads to an increase in fitness at one life-stage or in a specific environment, but simultaneously decreases fitness in another respect. Therefore, it is conceivable that the molecular underpinnings of evolution of highly complex traits, including brain size or cognitive ability, under certain conditions could result in deleterious effects, which would increase the susceptibility to psychiatric or neurodevelopmental diseases. Here, we discuss possible trade-offs and antagonistic pleiotropies between evolutionary change in a gene sequence, dosage or activity and the susceptibility of individuals to autism spectrum disorders and schizophrenia. We present current knowledge about genes and alterations in gene regulatory landscapes, which have likely played a role in establishing human-specific traits and have been implicated in those diseases.

Impact of schizophrenia GWAS loci converge onto distinct pathways in cortical interneurons vs glutamatergic neurons during development

Remarkable advances have been made in schizophrenia (SCZ) GWAS, but gleaning biological insight from these loci is challenging. Genetic influences on gene expression (e.g., eQTLs) are cell type-specific, but most studies that attempt to clarify GWAS loci's influence on gene expression have employed tissues with mixed cell compositions that can obscure cell-specific effects. Furthermore, enriched SCZ heritability in the fetal brain underscores the need to study the impact of SCZ risk loci in specific developing neurons. MGE-derived cortical interneurons (cINs) are consistently affected in SCZ brains and show enriched SCZ heritability in human fetal brains. We identified SCZ GWAS risk genes that are dysregulated in iPSC-derived homogeneous populations of developing SCZ cINs. These SCZ GWAS loci differential expression (DE) genes converge on the PKC pathway. Their disruption results in PKC hyperactivity in developing cINs, leading to arborization deficits. We show that the fine-mapped GWAS locus in the ATP2A2 gene of the PKC pathway harbors enhancer marks by ATACseq and ChIPseq, and regulates ATP2A2 expression. We also generated developing glutamatergic neurons (GNs), another population with enriched SCZ heritability, and confirmed their functionality after transplantation into the mouse brain. Then, we identified SCZ GWAS risk genes that are dysregulated in developing SCZ GNs. GN-specific SCZ GWAS loci DE genes converge on the ion transporter pathway, distinct from those for cINs. Disruption of the pathway gene CACNA1D resulted in deficits of Ca2+ currents in developing GNs, suggesting compromised neuronal function by GWAS loci pathway deficits during development. This study allows us to identify cell type-specific and developmental stage-specific mechanisms of SCZ risk gene function, and may aid in identifying mechanism-based novel therapeutic targets.

A novel genetic variant potentially altering the expression of MANBA in the cerebellum associated with attention deficit hyperactivity disorder in Han Chinese children

OBJECTIVES

To obtain additional insight into the genetic factors of attention deficit hyperactivity disorder (ADHD).

METHODS

First, we performed a transcriptome-wide association study (TWAS) integrating human cerebellum-specific variant-expression/splicing correlations to identify ADHD susceptibility genes. Then, the associations between expression/splicing quantitative trait loci (eQTLs/sQTLs) of the transcriptome-wide significant genes and ADHD were observed in a case-control study of Han Chinese children. Furthermore, dual luciferase reporter gene assays were performed to validate the regulatory function of ADHD risk variants. Additionally, the transcription level of target genes in blood was detected by real-time quantitative polymerase chain reaction (RT-qPCR) assay.

RESULTS

TWAS identified that the genetically regulated expression of MANBA in the cerebellum was significantly associated with ADHD risk. Furthermore, we observed a higher risk of ADHD and more severe clinical symptoms in subjects harbouring heterozygous (TC) or mutant homozygous (TT) genotypes of MANBA rs1054037 than CC carriers. The dual luciferase reporter gene assay revealed that the mutation of rs1054037(C > T) potentially upregulated MANBA expression by eliminating the binding site for hsa-miR-5591-3P. Finally, RT-qPCR showed that MANBA expression in blood samples of patients was significantly higher than that of controls.

CONCLUSIONS

Taken together, these results suggest a role of MANBA in the development of ADHD.

The neurodevelopmental gene MSANTD2 belongs to a gene family formed by recurrent molecular domestication of Harbinger transposons at the base of vertebrates

The formation of new genes is a major source of organism evolutionary innovation. Beyond their mutational effects, transposable elements can be co-opted by host genomes to form different types of sequences including novel genes, through a mechanism named molecular domestication.We report the formation of four genes through molecular domestication of Harbinger transposons, three in a common ancestor of jawed vertebrates about 500 million years ago and one in sarcopterygians approx. 430 million years ago. Additionally, one processed pseudogene arose approx. 60 million years ago in simians. In zebrafish, Harbinger-derived genes are expressed during early development but also in adult tissues, and predominantly co-expressed in male brain. In human, expression was detected in multiple organs, with major expression in the brain particularly during fetal development. We used CRISPR/Cas9 with direct gene knock-out in the F0 generation and the morpholino antisense oligonucleotide knock-down technique to study in zebrafish the function of one of these genes called MSANTD2, which has been suggested to be associated to neuro-developmental diseases such as autism spectrum disorders and schizophrenia in human. MSANTD2 inactivation led to developmental delays including tail and nervous system malformation at one day post fertilization. Affected embryos showed dead cell accumulation, major anatomical defects characterized by impaired brain ventricle formation and alterations in expression of some characteristic genes involved in vertebrate nervous system development. Hence, the characterization of MSANTD2 and other Harbinger-derived genes might contribute to a better understanding of the genetic innovations having driven the early evolution of the vertebrate nervous system.

Locus-level antagonistic selection shaped the polygenic architecture of human complex diseases

BACKGROUND

We aimed to evaluate the potential role of antagonistic selection in polygenic diseases: if one variant increases the risk of one disease and decreases the risk of another disease, the signals of genetic risk elimination by natural selection will be distorted, which leads to a higher frequency of risk alleles.

METHODS

We applied local genetic correlations and transcriptome-wide association studies to identify genomic loci and genes adversely associated with at least two diseases. Then, we used different population genetic metrics to measure the signals of natural selection for these loci and genes.

RESULTS

First, we identified 2120 cases of antagonistic pleiotropy (negative local genetic correlation) among 87 diseases in 716 genomic loci (antagonistic loci). Next, by comparing with non-antagonistic loci, we observed that antagonistic loci explained an excess proportion of disease heritability (median 6%), showed enhanced signals of balancing selection, and reduced signals of directional polygenic adaptation. Then, at the gene expression level, we identified 31,991 cases of antagonistic pleiotropy among 98 diseases at 4368 genes. However, evidence of altered signals of selection pressure and heritability distribution at the gene expression level is limited.

CONCLUSION

We conclude that antagonistic pleiotropy is widespread among human polygenic diseases, and it has distorted the evolutionary signal and genetic architecture of diseases at the locus level.

Characterizing the polygenic overlaps of bipolar disorder subtypes with schizophrenia and major depressive disorder

BACKGROUND

Large-scale studies have shown that bipolar I disorder (BD-I) and bipolar II disorder (BD-II) have differences in genetic association with schizophrenia (SCZ) and major depressive disorder (MDD). However, the underlying shared genetic architectures between BD subtypes and both SCZ and MDD remain largely unknown.

METHODS

We applied univariate and bivariate causal mixture models (MiXeR) to estimate the polygenicity and polygenic overlaps on large GWASs summary statistics of BD-I (n = 25,060), BD-II (n = 6781), SCZ (n = 69,369) and MDD (n = 170,756). Then, conjunctional false discovery rate approach was used to identify specific shared genetic loci between BD subtypes and both SCZ and MDD.

RESULTS

Univariate MiXeR revealed that BD-II was substantially more polygenic (22.37 K causal variants) as compared to BD-I, SCZ and MDD (7.87-12.43 K causal variants). Bivariate MiXeR revealed substantial polygenic overlaps between BD-I and SCZ (Dice-coefficient = 0.83) and between BD-I and MDD (Dice-coefficient = 0.76), which are beyond the genetic correlation (rg = 0.71 and 0.36). Conjunctional FDR analysis identified 236 distinct shared loci between BD-I and BD-II (2 loci), BD-I and SCZ (227 loci), BD-I and MDD (19 loci), BD-II and SCZ (1 locus), and BD-II and MDD (3 loci). Most of these shared loci have concordant effect directions among BD subtypes, SCZ and MDD.

LIMITATIONS

The bivariate MiXeR model was not applied for the BD-II because of insufficient power and inadequate model fit.

CONCLUSIONS

These findings provide evidence for extensive polygenic effects across BD subtypes, SCZ and MDD, which further our understanding of the potential genetic basis for the comorbid symptoms across these disorders.

Functional Genomics Analysis to Disentangle the Role of Genetic Variants in Major Depression

Understanding the molecular basis of major depression is critical for identifying new potential biomarkers and drug targets to alleviate its burden on society. Leveraging available GWAS data and functional genomic tools to assess regulatory variation could help explain the role of major depression-associated genetic variants in disease pathogenesis. We have conducted a fine-mapping analysis of genetic variants associated with major depression and applied a pipeline focused on gene expression regulation by using two complementary approaches: cis-eQTL colocalization analysis and alteration of transcription factor binding sites. The fine-mapping process uncovered putative causally associated variants whose proximal genes were linked with major depression pathophysiology. Four colocalizing genetic variants altered the expression of five genes, highlighting the role of SLC12A5 in neuronal chlorine homeostasis and MYRF in nervous system myelination and oligodendrocyte differentiation. The transcription factor binding analysis revealed the potential role of rs62259947 in modulating P4HTM expression by altering the YY1 binding site, altogether regulating hypoxia response. Overall, our pipeline could prioritize putative causal genetic variants in major depression. More importantly, it can be applied when only index genetic variants are available. Finally, the presented approach enabled the proposal of mechanistic hypotheses of these genetic variants and their role in disease pathogenesis.

dbBIP: a comprehensive bipolar disorder database for genetic research

Bipolar disorder (BIP) is one of the most common hereditary psychiatric disorders worldwide. Elucidating the genetic basis of BIP will play a pivotal role in mechanistic delineation. Genome-wide association studies (GWAS) have successfully reported multiple susceptibility loci conferring BIP risk, thus providing insight into the effects of its underlying pathobiology. However, difficulties remain in the extrication of important and biologically relevant data from genetic discoveries related to psychiatric disorders such as BIP. There is an urgent need for an integrated and comprehensive online database with unified access to genetic and multi-omics data for in-depth data mining. Here, we developed the dbBIP, a database for BIP genetic research based on published data. The dbBIP consists of several modules, i.e.: (i) single nucleotide polymorphism (SNP) module, containing large-scale GWAS genetic summary statistics and functional annotation information relevant to risk variants; (ii) gene module, containing BIP-related candidate risk genes from various sources and (iii) analysis module, providing a simple and user-friendly interface to analyze one’s own data. We also conducted extensive analyses, including functional SNP annotation, integration (including summary-data-based Mendelian randomization and transcriptome-wide association studies), co-expression, gene expression, tissue expression, protein–protein interaction and brain expression quantitative trait loci analyses, thus shedding light on the genetic causes of BIP. Finally, we developed a graphical browser with powerful search tools to facilitate data navigation and access. The dbBIP provides a comprehensive resource for BIP genetic research as well as an integrated analysis platform for researchers and can be accessed online at http://dbbip.xialab.info.

Database URL:  http://dbbip.xialab.info

Genomics, convergent neuroscience and progress in understanding autism spectrum disorder

More than a hundred genes have been identified that, when disrupted, impart large risk for autism spectrum disorder (ASD). Current knowledge about the encoded proteins - although incomplete - points to a very wide range of developmentally dynamic and diverse biological processes. Moreover, the core symptoms of ASD involve distinctly human characteristics, presenting challenges to interpreting evolutionarily distant model systems. Indeed, despite a decade of striking progress in gene discovery, an actionable understanding of pathobiology remains elusive. Increasingly, convergent neuroscience approaches have been recognized as an important complement to traditional uses of genetics to illuminate the biology of human disorders. These methods seek to identify intersection among molecular-level, cellular-level and circuit-level functions across multiple risk genes and have highlighted developing excitatory neurons in the human mid-gestational prefrontal cortex as an important pathobiological nexus in ASD. In addition, neurogenesis, chromatin modification and synaptic function have emerged as key potential mediators of genetic vulnerability. The continued expansion of foundational 'omics' data sets, the application of higher-throughput model systems and incorporating developmental trajectories and sex differences into future analyses will refine and extend these results. Ultimately, a systems-level understanding of ASD genetic risk holds promise for clarifying pathobiology and advancing therapeutics.

Schizophrenia Genomics: Convergence on Synaptic Development, Adult Synaptic Plasticity, or Both?

Large-scale genomic studies of schizophrenia have identified hundreds of genetic loci conferring risk to the disorder. This progress offers an important route toward defining the biological basis of the condition and potentially developing new treatments. In this review, we discuss insights from recent genome-wide association study, copy number variant, and exome sequencing analyses of schizophrenia, together with functional genomics data from the pre- and postnatal brain, in relation to synaptic development and function. These data provide strong support for the view that synaptic dysfunction within glutamatergic and GABAergic (gamma-aminobutyric acidergic) neurons of the cerebral cortex, hippocampus, and other limbic structures is a central component of schizophrenia pathophysiology. Implicated genes and functional genomic data suggest that disturbances in synaptic connectivity associated with susceptibility to schizophrenia begin in utero but continue throughout development, with some alleles conferring risk to the disorder through direct effects on synaptic function in adulthood. This model implies that novel interventions for schizophrenia could include broad preventive approaches aimed at enhancing synaptic health during development as well as more targeted treatments aimed at correcting synaptic function in affected adults.

Treatment Response and GWAS Risk Allele rs2514218 (C) of the Dopamine D2 Receptor Gene in Inpatients with Schizophrenia

INTRODUCTION

The pathophysiological mechanisms of acute schizophrenia are largely unknown, but it is widely accepted that dopamine D2 receptors (DRD2s) are involved in psychosis treatments for schizophrenic patients. We suggest that genetic variation in these receptors may play a role in patients' responses to commonly used antipsychotics, particularly D2-blockers.

METHODS

This study included adult patients with ICD-10 diagnoses of schizophrenia and current acute psychosis who were treated with antipsychotics. All patients underwent genotyping for DRD2 rs2514218 polymorphism. The definition of overall treatment response was based on changes in treatment scheme: no changes indicated a good response, and changes indicated a limited response.

RESULTS

There were 275 inpatients (38.1% of whom were female; mean age = 32.7 years, SD = 11.1 years) who met the inclusion criteria. Of the participants, 99 were good responders (34% of whom were female), and 176 were limited responders (40% of whom were female). No differences in demographic, premorbid, or disease characteristics were found. The number of patients that were homozygous for the risk allele was significantly greater in the limited response group than in the good response group.

CONCLUSION

Our findings suggest that the risk variant at the DRD2 locus can be used as an indicator for patients' responses to antipsychotics without direct DRD2-blocking, thereby shortening the time needed for drug selection.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Polymorphisms in the hypoxia inducible factor binding site of the macrophage migration inhibitory factor gene promoter in schizophrenia

Background

Macrophage migration inhibitory factor (MIF) is a multifunctional cytokine that promotes neurogenesis and neuroprotection. MIF is predominantly expressed in astrocytes in the brain. The serum MIF level and microsatellites/single nucleotide polymorphisms (SNPs) in the MIF gene promoter region are known to be associated with schizophrenia (SCZ). Interestingly, previous studies reported that hypoxia, an environmental risk factor for SCZ, induced MIF expression through binding of the hypoxia inducible factor (HIF)-1 to the hypoxia response element (HRE) in the MIF promoter.

Methods

We investigated the involvement of MIF in SCZ while focusing on the HIF pathway. First, we conducted an association study of the SNP rs17004038 (C>A) in the HRE of the MIF promoter between 1758 patients with SCZ and 1507 controls. Next, we investigated the effect of hypoxia on MIF expression in primary cultured astrocytes derived from neonatal mice forebrain.

Results

SNP rs17004038 was significantly associated with SCZ (p = 0.0424, odds ratio = 1.445), indicating that this SNP in the HRE of the MIF promoter was a genetic risk factor for SCZ. Hypoxia induced MIF mRNA expression and MIF protein production and increased HIF-1 binding to the MIF promoter, while the activity of the MIF promoter was suppressed by mutations in the HRE and by deletion of the HRE in astrocytes.

Conclusion

These results suggest that SNP rs17004038 in the HRE of the MIF promoter was significantly associated with SCZ and may be involved in the pathophysiology of SCZ via suppression of hypoxia and HIF pathway-induced MIF expression.

Interpreting Functional Impact of Genetic Variations by Network QTL for Genotype–Phenotype Association Study

An enormous challenge in the post-genome era is to annotate and resolve the consequences of genetic variation on diverse phenotypes. The genome-wide association study (GWAS) is a well-known method to identify potential genetic loci for complex traits from huge genetic variations, following which it is crucial to identify expression quantitative trait loci (eQTL). However, the conventional eQTL methods usually disregard the systematical role of single-nucleotide polymorphisms (SNPs) or genes, thereby overlooking many network-associated phenotypic determinates. Such a problem motivates us to recognize the network-based quantitative trait loci (QTL), i.e., network QTL (nQTL), which is to detect the cascade association as genotype → network → phenotype rather than conventional genotype → expression → phenotype in eQTL. Specifically, we develop the nQTL framework on the theory and approach of single-sample networks, which can identify not only network traits (e.g., the gene subnetwork associated with genotype) for analyzing complex biological processes but also network signatures (e.g., the interactive gene biomarker candidates screened from network traits) for characterizing targeted phenotype and corresponding subtypes. Our results show that the nQTL framework can efficiently capture associations between SNPs and network traits (i.e., edge traits) in various simulated data scenarios, compared with traditional eQTL methods. Furthermore, we have carried out nQTL analysis on diverse biological and biomedical datasets. Our analysis is effective in detecting network traits for various biological problems and can discover many network signatures for discriminating phenotypes, which can help interpret the influence of nQTL on disease subtyping, disease prognosis, drug response, and pathogen factor association. Particularly, in contrast to the conventional approaches, the nQTL framework could also identify many network traits from human bulk expression data, validated by matched single-cell RNA-seq data in an independent or unsupervised manner. All these results strongly support that nQTL and its detection framework can simultaneously explore the global genotype–network–phenotype associations and the underlying network traits or network signatures with functional impact and importance.

Transcriptional programs regulating neuronal differentiation are disrupted in DLG2 knockout human embryonic stem cells and enriched for schizophrenia and related disorders risk variants

Coordinated programs of gene expression drive brain development. It is unclear which transcriptional programs, in which cell-types, are affected in neuropsychiatric disorders such as schizophrenia. Here we integrate human genetics with transcriptomic data from differentiation of human embryonic stem cells into cortical excitatory neurons. We identify transcriptional programs expressed during early neurogenesis in vitro and in human foetal cortex that are down-regulated in DLG2^−/− lines. Down-regulation impacted neuronal differentiation and maturation, impairing migration, morphology and action potential generation. Genetic variation in these programs is associated with neuropsychiatric disorders and cognitive function, with associated variants predominantly concentrated in loss-of-function intolerant genes. Neurogenic programs also overlap schizophrenia GWAS enrichment previously identified in mature excitatory neurons, suggesting that pathways active during prenatal cortical development may also be associated with mature neuronal dysfunction. Our data from human embryonic stem cells, when combined with analysis of available foetal cortical gene expression data, de novo rare variants and GWAS statistics for neuropsychiatric disorders and cognition, reveal a convergence on transcriptional programs regulating excitatory cortical neurogenesis.

Coordinated programs of gene expression drive brain development. Here, the authors use human embryonic stem cells and foetal cortical tissue as well as available GWAS statistics and analysis of genetic variants associated with neuropsychiatric disorders and cognition revealing a convergence on transcriptional programs regulating excitatory cortical neurogenesis.

A comprehensive gene-centric pleiotropic association analysis for 14 psychiatric disorders with GWAS summary statistics

BACKGROUND

Recent genome-wide association studies (GWASs) have revealed the polygenic nature of psychiatric disorders and discovered a few of single-nucleotide polymorphisms (SNPs) associated with multiple psychiatric disorders. However, the extent and pattern of pleiotropy among distinct psychiatric disorders remain not completely clear.

METHODS

We analyzed 14 psychiatric disorders using summary statistics available from the largest GWASs by far. We first applied the cross-trait linkage disequilibrium score regression (LDSC) to estimate genetic correlation between disorders. Then, we performed a gene-based pleiotropy analysis by first aggregating a set of SNP-level associations into a single gene-level association signal using MAGMA. From a methodological perspective, we viewed the identification of pleiotropic associations across the entire genome as a high-dimensional problem of composite null hypothesis testing and utilized a novel method called PLACO for pleiotropy mapping. We ultimately implemented functional analysis for identified pleiotropic genes and used Mendelian randomization for detecting causal association between these disorders.

RESULTS

We confirmed extensive genetic correlation among psychiatric disorders, based on which these disorders can be grouped into three diverse categories. We detected a large number of pleiotropic genes including 5884 associations and 2424 unique genes and found that differentially expressed pleiotropic genes were significantly enriched in pancreas, liver, heart, and brain, and that the biological process of these genes was remarkably enriched in regulating neurodevelopment, neurogenesis, and neuron differentiation, offering substantial evidence supporting the validity of identified pleiotropic loci. We further demonstrated that among all the identified pleiotropic genes there were 342 unique ones linked with 6353 drugs with drug-gene interaction which can be classified into distinct types including inhibitor, agonist, blocker, antagonist, and modulator. We also revealed causal associations among psychiatric disorders, indicating that genetic overlap and causality commonly drove the observed co-existence of these disorders.

CONCLUSIONS

Our study is among the first large-scale effort to characterize gene-level pleiotropy among a greatly expanded set of psychiatric disorders and provides important insight into shared genetic etiology underlying these disorders. The findings would inform psychiatric nosology, identify potential neurobiological mechanisms predisposing to specific clinical presentations, and pave the way to effective drug targets for clinical treatment.

Impact of fetal expression quantitative trait loci on transcriptome-wide association study of childhood leukemia

Transcriptome-wide association studies increase the yield of loci associated with disease phenotypes by focusing on expression quantitative trait loci (eQTL). The major source of eQTL data for is the Gene and Tissue Expression (GTEx) project, which is comprised entirely of adults, mainly those > 50 years of age at death. Since gene expression levels differ by developmental stage, it is not clear whether eQTLs derived from adult data sources are best suited for use in young-onset diseases such as pediatric cancers. To fill in this knowledge gap, we performed a large-scale eQTL mapping analysis in the GenCord study with newborn samples and compared it with GTEx. Under matched conditions, we found around 80% of the eQTLs in one study can be replicated in the other. However, among all eQTLs identified in GenCord (GTEx), 584 (1045) showed statistically significant differences in effect sizes in GTEx (GenCord). We further investigated how using fetal eQTL data can facilitate the genetic association study of acute lymphoblastic leukemia. GenCord and GTEx identified the same genetic loci with statistical significance; however, the overall association pattern was only weakly correlated. Our paper demonstrates age-differential eQTLs and shows their potential influence on childhood leukemia research.

Induced Pluripotent Stem Cells in Psychiatry: An Overview and Critical Perspective

A key challenge in psychiatry research is the development of high-fidelity model systems that can be experimentally manipulated to explore and test pathophysiological mechanisms of illness. In this respect, the emerging capacity to derive neural cells and circuits from human induced pluripotent stem cells (iPSCs) has generated significant excitement. This review aims to provide a critical appraisal of the potential for iPSCs in illuminating pathophysiological mechanisms in the context of other available technical approaches. We discuss the selection of iPSC phenotypes relevant to psychiatry, the information that researchers can draw on to help guide these decisions, and how researchers choose between the use of 2-dimensional cultures and the use of more complex 3-dimensional model systems. We discuss the strengths and limitations of current models and the challenges and opportunities that they present. Finally, we discuss the potential of iPSC-based model systems for clarifying the mechanisms underlying genetic risk for psychiatry and the steps that will be needed to ensure that robust and reliable conclusions can be drawn. We argue that while iPSC-based models are ideally placed to study fundamental processes occurring within and between neural cells, they are often less well suited for case-control studies, given issues relating to statistical power and the challenges in identifying which cellular phenotypes are meaningful at the level of the whole individual. Our aim is to highlight the importance of considering the hypotheses of a given study to guide decisions about which, if any, iPSC-based system is most appropriate to address it.