Genetic Control of Expression and Splicing in Developing Human Brain Informs Disease Mechanisms

Integrative Multi-Omics Approach for Improving Causal Gene Identification

Transcriptome-wide association studies (TWAS) have been widely used to identify thousands of likely causal genes for diseases and complex traits using predicted expression models. However, most existing TWAS methods rely on gene expression alone and overlook other regulatory mechanisms of gene expression, including DNA methylation and splicing, that contribute to the genetic basis of these complex traits and diseases. Here we introduce a multi-omics method that integrates gene expression, DNA methylation, and splicing data to improve the identification of associated genes with our traits of interest. Through simulations and by analyzing genome-wide association study (GWAS) summary statistics for 24 complex traits, we show that our integrated method, which leverages these complementary omics biomarkers, achieves higher statistical power, and improves the accuracy of likely causal gene identification in blood tissues over individual omics methods. Finally, we apply our integrated model to a lung cancer GWAS data set, demonstrating the integrated models improved identification of prioritized genes for lung cancer risk.

Primary cartilage transcriptional signatures reflect cell-type-specific molecular pathways underpinning osteoarthritis

Translational efforts in osteoarthritis are hampered by a gap in our understanding of disease processes at the molecular level. Here, we present evidence of pronounced transcriptional changes in high- and low-disease-grade cartilage tissue, pointing to embryonic processes involved in disease progression. We identify shared transcriptional programs between osteoarthritis cartilage and cell populations in the human embryonic and fetal limb, pointing to increases in pre-hypertrophic chondrocytes' transcriptional programs in low-grade cartilage and increases in osteoblastic signatures in high-grade disease tissue. We find that osteoarthritis genetic risk signals are enriched in six gene co-expression modules and show that these transcriptional signatures reflect cell-type-specific expression along the endochondral ossification developmental trajectory. Using this network approach in combination with causal inference analysis, we present evidence of a causal effect on osteoarthritis risk for variants associated with the expression of ten genes that have not been previously reported as effector genes in genome-wide association studies in osteoarthritis. Our findings point to key molecular pathways as drivers of cartilage degeneration and identify high-value drug targets and repurposing opportunities.

Genetic architecture of RNA editing, splicing and gene expression in schizophrenia

Genome wide association studies (GWAS) have been conducted over the past decades to investigate the underlying genetic origin of neuropsychiatric diseases, such as schizophrenia (SCZ). While these studies demonstrated the significance of disease-phenotype associations, there is a pressing need to fully characterize the functional relevance of disease-associated genetic variants. Functional genetic loci can affect transcriptional and post-transcriptional phenotypes that may contribute to disease pathology. Here, we investigate the associations between genetic variation and RNA editing, splicing, and overall gene expression through identification of quantitative trait loci (QTL) in the CommonMind Consortium SCZ cohort. We find that editing QTL (edQTL), splicing QTL (sQTL) and expression QTL (eQTL) possess both unique and common gene targets, which are involved in many disease-relevant pathways, including brain function and immune response. We identified two QTL that fall into all three QTL categories (seedQTL), one of which, rs146498205, targets the lincRNA gene, RP11-156P1.3. In addition, we observe that the RNA binding protein AKAP1, with known roles in neuronal regulation and mitochondrial function, had enriched binding sites among edQTL, including the seedQTL, rs146498205. We conduct colocalization with various brain disorders and find that all QTL have top colocalizations with SCZ and related neuropsychiatric diseases. Furthermore, we identify QTL within biologically relevant GWAS loci, such as in ELA2, an important tRNA processing gene associated with SCZ risk. This work presents the investigation of multiple QTL types in parallel and demonstrates how they target both distinct and overlapping SCZ-relevant genes and pathways.

Brain multi-omic Mendelian randomisation to identify novel drug targets for gliomagenesis

BACKGROUND

Genetic variants associated with molecular traits that are also associated with liability to glioma can provide causal evidence for the identification and prioritisation of drug targets.

METHODS

We performed comprehensive two-sample Mendelian randomisation (Wald ratio and/or IVW) and colocalisation analyses of molecular traits on glioma. Instrumentable traits (QTLs P < 5 × 10-8) were identified amongst 11 985 gene expression measures, 13 285 splicing isoforms and 10 198 protein abundance measures, derived from 15 brain regions. Glioma summary-level data was extracted from a genome-wide association meta-analysis of 12 496 cases and 18 190 controls.

RESULTS

We found evidence for causal effect of 22 molecular traits (across 18 genes/proteins) on glioma risk. Thirteen molecular traits have been previously linked with glioma risk and five were novel; HBEGF (5q31.3) expression and all glioma [OR 1.36 (95%CI 1.19-1.55); P = 4.41 × 10-6]; a CEP192 (18p11.21) splice isoform and glioblastoma [OR 4.40 (95%CI 2.28-8.48); P = 9.78 × 10-4]; a FAIM (3q22.3) splice isoform and all glioma [OR 2.72-3.43; P = 1.03 × 10-5 to 1.09 × 10-5]; a SLC8A1 (2p22.1) splice isoform and all glioma [OR 0.37 (95%CI 0.24-0.56; P = 5.72 × 10-6]; D2HGDH (2q37.3) protein and all glioma [OR 0.86 (95%CI 0.80-0.92); P = 5.94 × 10-6)].

CONCLUSIONS

We provide robust causal evidence for prioritising genes and their protein products in glioma research. Our results highlight the importance of alternative splicing as a mechanism in gliomagenesis and as an avenue for exploration of drug targets.

Response splicing QTLs in primary human chondrocytes identifies putative osteoarthritis risk genes

Osteoarthritis affects millions worldwide, yet effective treatments remain elusive due to poorly understood molecular mechanisms. While genome-wide association studies (GWAS) have identified over 100 OA-associated loci, identifying the genes impacted at each locus remains challenging. Several studies have mapped expression quantitative trait loci (eQTL) in chondrocytes and colocalized them with OA GWAS variants to identify putative OA risk genes; however, the degree to which genetic variants influence OA risk via alternative splicing has not been explored. We investigated the role of alternative splicing in OA pathogenesis using RNA-seq data from 101 human chondrocyte samples treated with PBS (control) or fibronectin fragment (FN-f), an OA trigger. We identified 590 differentially spliced genes between conditions, with FN-f inducing splicing events similar to those in primary OA tissue. We used CRISPR/Cas9 to mimic an SNRNP70 splicing event observed in OA and FN-f-treated chondrocytes and found that it induced an OA-like expression pattern. Integration with genotyping data revealed 7,188 splicing quantitative trait loci (sQTL) affecting 3,056 genes. While many sQTLs were shared, we identified 738 and 343 condition-specific sQTLs for control and FN-f, respectively. We identified 15 RNA binding proteins whose binding sites were enriched at sQTL splice junctions and found that expression of those RNA binding proteins correlated with exon inclusion. Colocalization with OA GWAS identified 6 putative risk genes, including a novel candidate, PBRM1. Our study highlights the significant impact of alternative splicing in OA and provides potential therapeutic targets for future research.

A multi-omics Mendelian randomization study identifies new therapeutic targets for alcohol use disorder and problem drinking

Integrating proteomic and transcriptomic data with genetic architectures of problematic alcohol use and alcohol consumption behaviours can advance our understanding and help identify therapeutic targets. We conducted systematic screens using genome-wise association study data from ~3,500 cortical proteins (N = 722) and ~6,100 genes in 8 canonical brain cell types (N = 192) with 4 alcohol-related outcomes (N ≤ 537,349), identifying 217 cortical proteins and 255 cell-type genes associated with these behaviours, with 36 proteins and 37 cell-type genes being new. Although there was limited overlap between proteome and transcriptome targets, downstream neuroimaging revealed shared neurophysiological pathways. Colocalization with independent genome-wise association study data further prioritized 16 proteins, including CAB39L and NRBP1, and 12 cell-type genes, implicating mechanisms such as mTOR signalling. In addition, genes such as SAMHD1, VIPAS39, NUP160 and INO80E were identified as having favourable neuropsychiatric profiles. These findings provide insights into the genetic landscapes governing problematic alcohol use and alcohol consumption behaviours, highlighting promising therapeutic targets for future research.

Network-wide risk convergence in gene co-expression identifies reproducible genetic hubs of schizophrenia risk

The omnigenic model posits that genetic risk for traits with complex heritability involves cumulative effects of peripheral genes on mechanistic "core genes," suggesting that in a network of genes, those closer to clusters including core genes should have higher GWAS signals. In gene co-expression networks, we confirmed that GWAS signals accumulate in genes more connected to risk-enriched gene clusters, highlighting across-network risk convergence. This was strongest in adult psychiatric disorders, especially schizophrenia (SCZ), spanning 70% of network genes, suggestive of super-polygenic architecture. In snRNA-seq cell type networks, SCZ risk convergence was strongest in L2/L3 excitatory neurons. We prioritized genes most connected to SCZ-GWAS genes, which showed robust association to a CRISPRa measure of PGC3 regulation and were consistently identified across several brain regions. Several genes, including dopamine-associated ones, were prioritized specifically in the striatum. This strategy thus retrieves current drug targets and can be used to prioritize other potential drug targets.

Human brain proteome-wide association study provides insights into the genetic components of protein abundance in obesity

BACKGROUNDS

Genome-wide association studies have identified multiple genetic variants associated with obesity. However, most obesity-associated loci were waiting to be translated into new biological insights. Given the critical role of brain in obesity development, we sought to explore whether obesity-associated genetic variants could be mapped to brain protein abundances.

METHODS

We performed proteome-wide association studies (PWAS) and colocalization analyses to identify genes whose cis-regulated brain protein abundances were associated with obesity-related traits, including body fat percentage, trunk fat percentage, body mass index, visceral adipose tissue, waist circumference, and waist-to-hip ratio. We then assessed the druggability of the identified genes and conducted pathway enrichment analysis to explore their functional relevance. Finally, we evaluated the effects of the significant PWAS genes at the brain transcriptional level.

RESULTS

By integrating human brain proteomes from discovery (ROSMAP, N = 376) and validation datasets (BANNER, N = 198) with genome-wide summary statistics of obesity-related phenotypes (N ranged from 325,153 to 806,834), we identified 51 genes whose cis-regulated brain protein abundance was associated with obesity. These 51 genes were enriched in 11 metabolic processes, e.g., small molecule metabolic process and metabolic pathways. Fourteen of the 51 genes had high drug repurposing value. Ten of the 51 genes were also associated with obesity at the transcriptome level, suggesting that genetic variants likely confer risk of obesity by regulating mRNA expression and protein abundance of these genes.

CONCLUSIONS

Our study provides new insights into the genetic component of human brain protein abundance in obesity. The identified proteins represent promising therapeutic targets for future drug development.

Oxygen-induced stress reveals context-specific gene regulatory effects in human brain organoids

The interaction between genetic variants and environmental stressors is key to understanding the mechanisms underlying neurological diseases. In this study, we used human brain organoids to explore how varying oxygen levels expose context-dependent gene regulatory effects. By subjecting a genetically diverse panel of 21 brain organoids to hypoxic and hyperoxic conditions, we identified thousands of gene regulatory changes that are undetectable under baseline conditions, with 1,745 trait-associated genes showing regulatory effects only in response to oxygen stress. To capture more nuanced transcriptional patterns, we employed topic modeling, which revealed context-specific gene regulation linked to dynamic cellular processes and environmental responses, offering a deeper understanding of how gene regulation is modulated in the brain. These findings underscore the importance of genotype-environment interactions in genetic studies of neurological disorders and provide new insights into the hidden regulatory mechanisms influenced by environmental factors in the brain.

The identification of biomarkers for Alzheimer's disease using a systems biology approach based on lncRNA-circRNA-miRNA-mRNA ceRNA networks

In addition to being the most prevalent form of neurodegeneration among the elderly, AD is a devastating multifactorial disease. Currently, treatments address only its symptoms. Several clinical studies have shown that the disease begins to manifest decades before the first symptoms appear, indicating that studying early changes is crucial to improving early diagnosis and discovering novel treatments. Our study used bioinformatics and systems biology to identify biomarkers in AD that could be used for diagnosis and prognosis. The procedure was performed on data from the GEO database, and GO and KEGG enrichment analysis were performed. Then, we set up a network of interactions between proteins. Several miRNA prediction tools including miRDB, miRWalk, and TargetScan were used. The ceRNA network led to the identification of eight mRNAs, four circRNAs, seven miRNAs, and seven lncRNAs. Multiple mechanisms, including the cell cycle and DNA replication, have been linked to the promotion of AD development by the ceRNA network. By using the ceRNA network, it should be possible to extract prospective biomarkers and therapeutic targets for the treatment of AD. It is possible that the processes involved in DNA cell cycle and the replication of DNA contribute to the development of Alzheimer's disease.

Schizophrenia genomics: genetic complexity and functional insights

Determining the causes of schizophrenia has been a notoriously intractable problem, resistant to a multitude of investigative approaches over centuries. In recent decades, genomic studies have delivered hundreds of robust findings that implicate nearly 300 common genetic variants (via genome-wide association studies) and more than 20 rare variants (via whole-exome sequencing and copy number variant studies) as risk factors for schizophrenia. In parallel, functional genomic and neurobiological studies have provided exceptionally detailed information about the cellular composition of the brain and its interconnections in neurotypical individuals and, increasingly, in those with schizophrenia. Taken together, these results suggest unexpected complexity in the mechanisms that drive schizophrenia, pointing to the involvement of ensembles of genes (polygenicity) rather than single-gene causation. In this Review, we describe what we now know about the genetics of schizophrenia and consider the neurobiological implications of this information.

Colocalised Genetic Associations Reveal Alternative Splicing Variants as Candidate Causal Links for Breast Cancer Risk in 10 Loci

Genome-wide association studies (GWASs) have revealed numerous loci associated with breast cancer risk, yet the precise causal variants, their impact on molecular mechanisms, and the affected genes often remain elusive. We hypothesised that specific variants exert their influence by affecting cis-regulatory alternative splice elements. An analysis of splicing quantitative trait loci (sQTL) in healthy breast tissue from female individuals identified multiple variants linked to alterations in splicing ratios. Through colocalisation analysis, we pinpointed 43 variants within twelve genes that serve as candidate causal links between sQTL and GWAS findings. In silico splice analysis highlighted a potential mechanism for three genes-FDPS, SGCE, and MRPL11-where variants in proximity to or on the splice site modulate usage, resulting in alternative splice transcripts. Further in vitro/vivo studies are imperative to fully understand how these identified changes contribute to breast oncogenesis. Moreover, investigating their potential as biomarkers for breast cancer risk could enhance screening strategies and early detection methods for breast cancer.

Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing

Quantitative trait loci analyses have revealed an important role for genetic variants in regulating alternative splicing (AS) and alternative cleavage and polyadenylation (APA) in humans. Yet, these studies are generally performed with mature mRNA, so they report on the outcome rather than the processes of RNA maturation and thus may overlook how variants directly modulate pre-mRNA processing. The order in which the many introns of a human gene are removed can substantially influence AS, while nascent RNA polyadenylation can affect RNA stability and decay. However, how splicing order and poly(A) tail length are regulated by genetic variation has never been explored. Here, we used direct RNA nanopore sequencing to investigate allele-specific pre-mRNA maturation in 12 human lymphoblastoid cell lines. We found frequent splicing order differences between alleles and uncovered significant single nucleotide polymorphism (SNP)-splicing order associations in 17 genes. This included SNPs located in or near splice sites as well as more distal intronic and exonic SNPs. Moreover, several genes showed allele-specific poly(A) tail lengths, many of which also had a skewed allelic abundance ratio. HLA class I transcripts, which encode proteins that play an essential role in antigen presentation, showed the most allele-specific splicing orders, which frequently co-occurred with allele-specific AS, APA or poly(A) tail length differences. Together, our results expose new layers of genetic regulation of pre-mRNA maturation and highlight the power of long-read RNA sequencing for allele-specific analyses.

Genetic Implication of Prenatal GABAergic and Cholinergic Neuron Development in Susceptibility to Schizophrenia

BACKGROUND

The ganglionic eminences (GE) are fetal-specific structures that give rise to gamma-aminobutyric acid (GABA)- and acetylcholine-releasing neurons of the forebrain. Given the evidence for GABAergic, cholinergic, and neurodevelopmental disturbances in schizophrenia, we tested the potential involvement of GE neuron development 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 GE to test the enrichment of schizophrenia risk variation in genes with high expression specificity for developing GE cell populations. 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 GE, using these to 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 concentrated in predicted regulatory genomic sequence mapped in developing neuronal populations of the GE.

CONCLUSIONS

Our study implicates prenatal development of specific populations of GABAergic and cholinergic neurons in later susceptibility to schizophrenia, and provides a map of predicted regulatory genomic elements operating in cells of the GE.

The association between genetic variations and morphology-based brain networks changes in Alzheimer's disease

Alzheimer's disease (AD) is a highly heritable disease. The morphological changes of cortical cortex (such as, cortical thickness and surface area) in AD always accompany by the change of the functional connectivity to other brain regions and influence the short- and long-range brain network connections, causing functional deficits of AD. In this study, the first hypothesis is that genetic variations might affect morphology-based brain networks, leading to functional deficits; the second hypothesis is that protein-protein interaction (PPI) between the candidate proteins and known interacting proteins to AD might exist and influence AD. 600 470 variants and structural magnetic resonance imaging scans from 175 AD patients and 214 healthy controls were obtained from the Alzheimer's Disease Neuroimaging Initiative-1 database. A co-sparse reduced-rank regression model was fit to study the relationship between non-synonymous mutations and morphology-based brain networks. After that, PPIs between selected genes and BACE1, an enzyme that was known to be related to AD, are explored by using molecular dynamics (MD) simulation and co-immunoprecipitation (Co-IP) experiments. Eight genes affecting morphology-based brain networks were identified. The results of MD simulation showed that the PPI between TGM4 and BACE1 was the strongest among them and its interaction was verified by Co-IP. Hence, gene variations influence morphology-based brain networks in AD, leading to functional deficits. This finding, validated by MD simulation and Co-IP, suggests that the effect is robust.

Critical reasoning on the co-expression module QTL in the dorsolateral prefrontal cortex

Expression quantitative trait locus (eQTL) analysis is a popular method of gaining insight into the function of regulatory variation. While cis-eQTL resources have been instrumental in linking genome-wide association study variants to gene function, complex trait heritability may be additionally mediated by other forms of gene regulation. Toward this end, novel eQTL methods leverage gene co-expression (module-QTL) to investigate joint regulation of gene modules by single genetic variants. Here we broadly define a "module-QTL" as the association of a genetic variant with a summary measure of gene co-expression. This approach aims to reduce the multiple testing burden of a trans-eQTL search through the consolidation of gene-based testing and provide biological context to eQTLs shared between genes. In this article we provide an in-depth examination of the co-expression module eQTL (module-QTL) through literature review, theoretical investigation, and real-data application of the module-QTL to three large prefrontal cortex genotype-RNA sequencing datasets. We find module-QTLs in our study that are disease associated and reproducible are not additionally informative beyond cis- or trans-eQTLs for module genes. Through comparison to prior studies, we highlight promises and limitations of the module-QTL across study designs and provide recommendations for further investigation of the module-QTL framework.

Examining the associations between microglia genetic capacity, early life exposures and white matter development at the level of the individual

There are inter-individual differences in susceptibility to the influence of early life experiences for which the underlying neurobiological mechanisms are poorly understood. Microglia play a role in environmental surveillance and may influence individual susceptibility to environmental factors. As an index of neurodevelopment, we estimated individual slopes of mean white matter fractional anisotropy (WM-FA) across three time-points (age 4.5, 6.0, and 7.5 years) for 351 participants. Individual variation in microglia reactivity was derived from an expression-based polygenic score(ePGS) comprised of Single Nucleotide Polymorphisms (SNPs) functionally related to the expression of microglia-enriched genes.A higher ePGS denotes an increased genetic capacity for the expression of microglia-related genes, and thus may confer a greater capacity to respond to the early environment and to influence brain development. We hypothesized that this ePGS would associate with the WM-FA index of neurodevelopment and moderate the influence of early environmental factors.Our findings show sex dependency, where a significant association between WM-FA and microglia ePGS was only obtained for females.We then examined associations with perinatal factors known to decrease (optimal birth outcomes and familial conditions) or increase (systemic inflammation) the risk for later mental health problems.In females, individuals with high microglia ePGS showed a negative association between systemic inflammation and WM-FA and a positive association between more advantageous environmental conditions and WM-FA. The microglia ePGS in females thus accounted for variations in the influence of the quality of the early environment on WM-FA.Finally, WM-FA slopes mediated the association of microglia ePGS with interpersonal problems and social hostility in females. Our findings suggest the genetic capacity for microglia function as a potential factor underlying differential susceptibility to early life exposuresthrough influences on neurodevelopment.

Microglial function interacts with the environment to affect sex-specific depression risk

There is a two-fold higher incidence of depression in females compared to men with recent studies suggesting a role for microglia in conferring this sex-dependent depression risk. In this study we investigated the nature of this relation. Using GWAS enrichment, gene-set enrichment analysis and Mendelian randomization, we found minimal evidence for a direct relation between genes functionally related to microglia and sex-dependent genetic risk for depression. We then used expression quantitative trait loci and single nucleus RNA-sequencing resources to generate polygenic scores (PGS) representative of individual variation in microglial function in the adult (UK Biobank; N = 54753-72682) and fetal (ALSPAC; N = 1452) periods. The adult microglial PGS moderated the association between BMI (UK Biobank; beta = 0.001, 95 %CI 0.0009 to 0.003, P = 7.74E-6) and financial insecurity (UK Biobank; beta = 0.001, 95 %CI 0.005 to 0.015, P = 2E-4) with depressive symptoms in females. The fetal microglia PGS moderated the association between maternal prenatal depressive symptoms and offspring depressive symptoms at 24 years in females (ALSPAC; beta = 0.04, 95 %CI 0.004 to 0.07, P = 0.03). We found no evidence for an interaction between the microglial PGS and depression risk factors in males. Our results illustrate a role for microglial function in the conferral of sex-dependent depression risk following exposure to a depression risk factor.

Sex-Differential Gene Expression in Developing Human Cortex and Its Intersection With Autism Risk Pathways

Background

Sex-differential biology may contribute to the consistently male-biased prevalence of autism spectrum disorder (ASD). Gene expression differences between males and females in the brain can indicate possible molecular and cellular mechanisms involved, although transcriptomic sex differences during human prenatal cortical development have been incompletely characterized, primarily due to small sample sizes.

Methods

We performed a meta-analysis of sex-differential expression and co-expression network analysis in 2 independent bulk RNA sequencing datasets generated from cortex of 273 prenatal donors without known neuropsychiatric disorders. To assess the intersection between neurotypical sex differences and neuropsychiatric disorder biology, we tested for enrichment of ASD-associated risk genes and expression changes, neuropsychiatric disorder risk genes, and cell type markers within identified sex-differentially expressed genes (sex-DEGs) and sex-differential co-expression modules.

Results

We identified 101 significant sex-DEGs, including Y-chromosome genes, genes impacted by X-chromosome inactivation, and autosomal genes. Known ASD risk genes, implicated by either common or rare variants, did not preferentially overlap with sex-DEGs. We identified 1 male-specific co-expression module enriched for immune signaling that is unique to 1 input dataset.

Conclusions

Sex-differential gene expression is limited in prenatal human cortex tissue, although meta-analysis of large datasets allows for the identification of sex-DEGs, including autosomal genes that encode proteins involved in neural development. Lack of sex-DEG overlap with ASD risk genes in the prenatal cortex suggests that sex-differential modulation of ASD symptoms may occur in other brain regions, at other developmental stages, or in specific cell types, or may involve mechanisms that act downstream from mutation-carrying genes.

A PLURIPOTENT STEM CELL PLATFORM FOR IN VITRO SYSTEMS GENETICS STUDIES OF MOUSE DEVELOPMENT

The directed differentiation of pluripotent stem cells (PSCs) from panels of genetically diverse individuals is emerging as a powerful experimental system for characterizing the impact of natural genetic variation on developing cell types and tissues. Here, we establish new PSC lines and experimental approaches for modeling embryonic development in a genetically diverse, outbred mouse stock (Diversity Outbred mice). We show that a range of inbred and outbred PSC lines can be stably maintained in the primed pluripotent state (epiblast stem cells -- EpiSCs) and establish the contribution of genetic variation to phenotypic differences in gene regulation and directed differentiation. Using pooled in vitro fertilization, we generate and characterize a genetic reference panel of Diversity Outbred PSCs (n = 230). Finally, we demonstrate the feasibility of pooled culture of Diversity Outbred EpiSCs as "cell villages", which can facilitate the differentiation of large numbers of EpiSC lines for forward genetic screens. These data can complement and inform similar efforts within the stem cell biology and human genetics communities to model the impact of natural genetic variation on phenotypic variation and disease-risk.

The Genesis of Schizophrenia: An Origin Story

Schizophrenia is routinely referred to as a neurodevelopmental disorder, but the role of brain development in a disorder typically diagnosed during early adult life is enigmatic. The authors revisit the neurodevelopmental model of schizophrenia with genomic insights from the most recent schizophrenia clinical genetic association studies, transcriptomic and epigenomic analyses from human postmortem brain studies, and analyses from cellular models that recapitulate neurodevelopment. Emerging insights into schizophrenia genetic risk continue to converge on brain development, particularly stages of early brain development, that may be perturbed to deviate from a typical, normative course, resulting in schizophrenia clinical symptomatology. As the authors explicate, schizophrenia genetic risk is likely dynamic and context dependent, with effects of genetic risk varying spatiotemporally, across the neurodevelopmental continuum. Optimizing therapeutic strategies for the heterogeneous collective of individuals with schizophrenia may likely be guided by leveraging markers of genetic risk and derivative functional insights, well before the emergence of psychosis. Ultimately, rather than a focus on therapeutic intervention during adolescence or adulthood, principles of prediction and prophylaxis in the pre- and perinatal and neonatal stages may best comport with the biology of schizophrenia to address the early-stage perturbations that alter the normative neurodevelopmental trajectory.

Functional neurogenomics in autism spectrum disorders: A decade of progress

Advances in autism spectrum disorder (ASD) genetics have identified many genetic causes, reflecting remarkable progress while at the same time identifying challenges such as heterogeneity and pleiotropy, which complicate attempts to connect genetic risk to mechanisms. High-throughput functional genomic approaches have yielded progress by defining a molecular pathology in the brain of individuals with ASD and in identifying convergent biological pathways through which risk genes are predicted to act. These studies indicate that ASD genetic risk converges in early brain development, primarily during the period of cortical neurogenesis. Over development, genetic perturbations in turn lead to broad neuronal signaling dysregulation, most prominent in glutamatergic cortical-cortical projecting neurons and somatostatin positive interneurons, which is accompanied by glial dyshomeostasis throughout the cerebral cortex. Connecting these developmental perturbations to disrupted neuronal and glial function in the postnatal brain is an important direction in current research. Coupling functional genomic approaches with advances in induced pluripotent stem cell-derived neural organoid development provides a promising avenue for connecting brain pathology to developmental mechanisms.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

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.

Molecular cascades and cell type-specific signatures in ASD revealed by single-cell genomics

Genomic profiling in postmortem brain from autistic individuals has consistently revealed convergent molecular changes. What drives these changes and how they relate to genetic susceptibility in this complex condition are not well understood. We performed deep single-nucleus RNA sequencing (snRNA-seq) to examine cell composition and transcriptomics, identifying dysregulation of cell type-specific gene regulatory networks (GRNs) in autism spectrum disorder (ASD), which we corroborated using single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) and spatial transcriptomics. Transcriptomic changes were primarily cell type specific, involving multiple cell types, most prominently interhemispheric and callosal-projecting neurons, interneurons within superficial laminae, and distinct glial reactive states involving oligodendrocytes, microglia, and astrocytes. Autism-associated GRN drivers and their targets were enriched in rare and common genetic risk variants, connecting autism genetic susceptibility and cellular and circuit alterations in the human brain.

Identification and characterization of whole blood gene expression and splicing quantitative trait loci during early to mid-lactation of dairy cattle

BACKGROUND

Characterization of regulatory variants (e.g., gene expression quantitative trait loci, eQTL; gene splicing QTL, sQTL) is crucial for biologically interpreting molecular mechanisms underlying loci associated with complex traits. However, regulatory variants in dairy cattle, particularly in specific biological contexts (e.g., distinct lactation stages), remain largely unknown. In this study, we explored regulatory variants in whole blood samples collected during early to mid-lactation (22-150 days after calving) of 101 Holstein cows and analyzed them to decipher the regulatory mechanisms underlying complex traits in dairy cattle.

RESULTS

We identified 14,303 genes and 227,705 intron clusters expressed in the white blood cells of 101 cattle. The average heritability of gene expression and intron excision ratio explained by cis-SNPs is 0.28 ± 0.13 and 0.25 ± 0.13, respectively. We identified 23,485 SNP-gene expression pairs and 18,166 SNP-intron cluster pairs in dairy cattle during early to mid-lactation. Compared with the 2,380,457 cis-eQTLs reported to be present in blood in the Cattle Genotype-Tissue Expression atlas (CattleGTEx), only 6,114 cis-eQTLs (P < 0.05) were detected in the present study. By conducting colocalization analysis between cis-e/sQTL and the results of genome-wide association studies (GWAS) from four traits, we identified a cis-e/sQTL (rs109421300) of the DGAT1 gene that might be a key marker in early to mid-lactation for milk yield, fat yield, protein yield, and somatic cell score (PP4 > 0.6). Finally, transcriptome-wide association studies (TWAS) revealed certain genes (e.g., FAM83H and TBC1D17) whose expression in white blood cells was significantly (P < 0.05) associated with complex traits.

CONCLUSIONS

This study investigated the genetic regulation of gene expression and alternative splicing in dairy cows during early to mid-lactation and provided new insights into the regulatory mechanisms underlying complex traits of economic importance.

xWAS analysis in neuropsychiatric disorders by integrating multi-molecular phenotype quantitative trait loci and GWAS summary data

BACKGROUND

Integrating quantitative trait loci (QTL) data related to molecular phenotypes with genome-wide association study (GWAS) data is an important post-GWAS strategic approach employed to identify disease-associated molecular features. Various types of molecular phenotypes have been investigated in neuropsychiatric disorders. However, these findings pertaining to distinct molecular features are often independent of each other, posing challenges for having an overview of the mapped genes.

METHODS

In this study, we comprehensively summarized published analyses focusing on four types of risk-related molecular features (gene expression, splicing transcriptome, protein abundance, and DNA methylation) across five common neuropsychiatric disorders. Subsequently, we conducted supplementary analyses with the latest GWAS dataset and corresponding deficient molecular phenotypes using Functional Summary-based Imputation (FUSION) and summary data-based Mendelian randomization (SMR). Based on the curated and supplemented results, novel reliable genes and their functions were explored.

RESULTS

Our findings revealed that eQTL exhibited superior ability in prioritizing risk genes compared to the other QTL, followed by sQTL. Approximately half of the genes associated with splicing transcriptome, protein abundance, and DNA methylation were successfully replicated by eQTL-associated genes across all five disorders. Furthermore, we identified 436 novel reliable genes, which enriched in pathways related with neurotransmitter transportation such as synaptic, dendrite, vesicles, axon along with correlations with other neuropsychiatric disorders. Finally, we identified ten multiple molecular involved regulation patterns (MMRP), which may provide valuable insights into understanding the contribution of molecular regulation network targeting these disease-associated genes.

CONCLUSIONS

The analyses prioritized novel and reliable gene sets related with five molecular features based on published and supplementary results for five common neuropsychiatric disorders, which were missed in the original GWAS analysis. Besides, the involved MMRP behind these genes could be given priority for further investigation to elucidate the pathogenic molecular mechanisms underlying neuropsychiatric disorders in future studies.

From Vessels to Neurons-The Role of Hypoxia Pathway Proteins in Embryonic Neurogenesis

Embryonic neurogenesis can be defined as a period of prenatal development during which divisions of neural stem and progenitor cells give rise to neurons. In the central nervous system of most mammals, including humans, the majority of neocortical neurogenesis occurs before birth. It is a highly spatiotemporally organized process whose perturbations lead to cortical malformations and dysfunctions underlying neurological and psychiatric pathologies, and in which oxygen availability plays a critical role. In case of deprived oxygen conditions, known as hypoxia, the hypoxia-inducible factor (HIF) signaling pathway is activated, resulting in the selective expression of a group of genes that regulate homeostatic adaptations, including cell differentiation and survival, metabolism and angiogenesis. While a physiological degree of hypoxia is essential for proper brain development, imbalanced oxygen levels can adversely affect this process, as observed in common obstetrical pathologies such as prematurity. This review comprehensively explores and discusses the current body of knowledge regarding the role of hypoxia and the HIF pathway in embryonic neurogenesis of the mammalian cortex. Additionally, it highlights existing gaps in our understanding, presents unanswered questions, and provides avenues for future research.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

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.

Enhancer-targeted CRISPR-Activation Rescues Haploinsufficient Autism Susceptibility Genes

Autism Spectrum Disorder (ASD) is a highly heritable condition with diverse clinical presentations. Approximately 20% of ASD's genetic susceptibility is imparted by de novo mutations of major effect, most of which cause haploinsufficiency. We mapped enhancers of two high confidence autism genes - CHD8 and SCN2A and used CRISPR-based gene activation (CRISPR-A) in hPSC-derived excitatory neurons and cerebral forebrain organoids to correct the effects of haploinsufficiency, taking advantage of the presence of a wildtype allele of each gene and endogenous gene regulation. We found that CRISPR-A induced a sustained increase in CHD8 and SCN2A expression in treated neurons and organoids, with rescue of gene expression levels and mutation-associated phenotypes, including gene expression and physiology. These data support gene activation via targeting enhancers of haploinsufficient genes, as a therapeutic intervention in ASD and other neurodevelopmental disorders.

The spectrum of pre-mRNA splicing in autism

Disruptions in spatiotemporal gene expression can result in atypical brain function. Specifically, autism spectrum disorder (ASD) is characterized by abnormalities in pre-mRNA splicing. Abnormal splicing patterns have been identified in the brains of individuals with ASD, and mutations in splicing factors have been found to contribute to neurodevelopmental delays associated with ASD. Here we review studies that shed light on the importance of splicing observed in ASD and that explored the intricate relationship between splicing factors and ASD, revealing how disruptions in pre-mRNA splicing may underlie ASD pathogenesis. We provide an overview of the research regarding all splicing factors associated with ASD and place a special emphasis on five specific splicing factors-HNRNPH2, NOVA2, WBP4, SRRM2, and RBFOX1-known to impact the splicing of ASD-related genes. In the discussion of the molecular mechanisms influenced by these splicing factors, we lay the groundwork for a deeper understanding of ASD's complex etiology. Finally, we discuss the potential benefit of unraveling the connection between splicing and ASD for the development of more precise diagnostic tools and targeted therapeutic interventions. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Evolution and Genomics > Computational Analyses of RNA RNA-Based Catalysis > RNA Catalysis in Splicing and Translation.

A distinct class of pan-cancer susceptibility genes revealed by an alternative polyadenylation transcriptome-wide association study

Alternative polyadenylation plays an important role in cancer initiation and progression; however, current transcriptome-wide association studies mostly ignore alternative polyadenylation when identifying putative cancer susceptibility genes. Here, we perform a pan-cancer 3' untranslated region alternative polyadenylation transcriptome-wide association analysis by integrating 55 well-powered (n > 50,000) genome-wide association studies datasets across 22 major cancer types with alternative polyadenylation quantification from 23,955 RNA sequencing samples across 7,574 individuals. We find that genetic variants associated with alternative polyadenylation are co-localized with 28.57% of cancer loci and contribute a significant portion of cancer heritability. We further identify 642 significant cancer susceptibility genes predicted to modulate cancer risk via alternative polyadenylation, 62.46% of which have been overlooked by traditional expression- and splicing- studies. As proof of principle validation, we show that alternative alleles facilitate 3' untranslated region lengthening of CRLS1 gene leading to increased protein abundance and promoted proliferation of breast cancer cells. Together, our study highlights the significant role of alternative polyadenylation in discovering new cancer susceptibility genes and provides a strong foundational framework for enhancing our understanding of the etiology underlying human cancers.

Polygenic profiles define aspects of clinical heterogeneity in attention deficit hyperactivity disorder

Attention deficit hyperactivity disorder (ADHD) is a complex disorder that manifests variability in long-term outcomes and clinical presentations. The genetic contributions to such heterogeneity are not well understood. Here we show several genetic links to clinical heterogeneity in ADHD in a case-only study of 14,084 diagnosed individuals. First, we identify one genome-wide significant locus by comparing cases with ADHD and autism spectrum disorder (ASD) to cases with ADHD but not ASD. Second, we show that cases with ASD and ADHD, substance use disorder and ADHD, or first diagnosed with ADHD in adulthood have unique polygenic score (PGS) profiles that distinguish them from complementary case subgroups and controls. Finally, a PGS for an ASD diagnosis in ADHD cases predicted cognitive performance in an independent developmental cohort. Our approach uncovered evidence of genetic heterogeneity in ADHD, helping us to understand its etiology and providing a model for studies of other disorders.

Exploring noncoding variants in genetic diseases: from detection to functional insights

Previous studies on genetic diseases predominantly focused on protein-coding variations, overlooking the vast noncoding regions in the human genome. The development of high-throughput sequencing technologies and functional genomics tools has enabled the systematic identification of functional noncoding variants. These variants can impact gene expression, regulation, and chromatin conformation, thereby contributing to disease pathogenesis. Understanding the mechanisms that underlie the impact of noncoding variants on genetic diseases is indispensable for the development of precisely targeted therapies and the implementation of personalized medicine strategies. The intricacies of noncoding regions introduce a multitude of challenges and research opportunities. In this review, we introduce a spectrum of noncoding variants involved in genetic diseases, along with research strategies and advanced technologies for their precise identification and in-depth understanding of the complexity of the noncoding genome. We will delve into the research challenges and propose potential solutions for unraveling the genetic basis of rare and complex diseases.

Molecular quantitative trait loci in reproductive tissues impact male fertility in cattle

Breeding bulls are well suited to investigate inherited variation in male fertility because they are genotyped and their reproductive success is monitored through semen analyses and thousands of artificial inseminations. However, functional data from relevant tissues are lacking in cattle, which prevents fine-mapping fertility-associated genomic regions. Here, we characterize gene expression and splicing variation in testis, epididymis, and vas deferens transcriptomes of 118 mature bulls and conduct association tests between 414,667 molecular phenotypes and 21,501,032 genome-wide variants to identify 41,156 regulatory loci. We show broad consensus in tissue-specific and tissue-enriched gene expression between the three bovine tissues and their human and murine counterparts. Expression- and splicing-mediating variants are more than three times as frequent in testis than epididymis and vas deferens, highlighting the transcriptional complexity of testis. Finally, we identify genes (WDR19, SPATA16, KCTD19, ZDHHC1) and molecular phenotypes that are associated with quantitative variation in male fertility through transcriptome-wide association and colocalization analyses.

Population Neuroscience: Principles and Advances

In population neuroscience, three disciplines come together to advance our knowledge of factors that shape the human brain: neuroscience, genetics, and epidemiology (Paus, Human Brain Mapping 31:891-903, 2010). Here, I will come back to some of the background material reviewed in more detail in our previous book (Paus, Population Neuroscience, 2013), followed by a brief overview of current advances and challenges faced by this integrative approach.

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.

Transcriptomics : Approaches to Quantifying Gene Expression and Their Application to Studying the Human Brain

To date, the field of transcriptomics has been characterized by rapid methods development and technological advancement, with new technologies continuously rendering older ones obsolete.This chapter traces the evolution of approaches to quantifying gene expression and provides an overall view of the current state of the field of transcriptomics, its applications to the study of the human brain, and its place in the broader emerging multiomics landscape.

Concurrent impact of de novo mutations on cranial and cortical development in nonsyndromic craniosynostosis

OBJECTIVE

Nonsyndromic craniosynostosis (nsCS), characterized by premature cranial suture fusion, is considered a primary skull disorder in which impact on neurodevelopment, if present, results from the mechanical hindrance of brain growth. Despite surgical repair of the cranial defect, neurocognitive deficits persist in nearly half of affected children. Therefore, the authors performed a functional genomics analysis of nsCS to determine when, where, and in what cell types nsCS-associated genes converge during development.

METHODS

The authors integrated whole-exome sequencing data from 291 nsCS proband-parent trios with 29,803 single-cell transcriptomes of the prenatal and postnatal neurocranial complex to inform when, where, and in what cell types nsCS-mutated genes might exert their pathophysiological effects.

RESULTS

The authors found that nsCS-mutated genes converged in cranial osteoprogenitors and pial fibroblasts and their transcriptional networks that regulate both skull ossification and cerebral neurogenesis. Nonsyndromic CS-mutated genes also converged in inhibitory neurons and gene coexpression modules that overlapped with autism and other developmental disorders. Ligand-receptor cell-cell communication analysis uncovered crosstalk between suture osteoblasts and neurons via the nsCS-associated BMP, FGF, and noncanonical WNT signaling pathways.

CONCLUSIONS

These data implicate a concurrent impact of nsCS-associated de novo mutations on cranial morphogenesis and cortical development via cell- and non-cell-autonomous mechanisms in a developmental nexus of fetal osteoblasts, pial fibroblasts, and neurons. These results suggest that neurodevelopmental outcomes in nsCS patients may be driven more by mutational status than surgical technique.

FEZ1 participates in human embryonic brain development by modulating neuronal progenitor subpopulation specification and migrations

Mutations in the human fasciculation and elongation protein zeta 1 (FEZ1) gene are found in schizophrenia and Jacobsen syndrome patients. Here, using human cerebral organoids (hCOs), we show that FEZ1 expression is turned on early during brain development and is detectable in both neuroprogenitor subtypes and immature neurons. FEZ1 deletion disrupts expression of neuronal and synaptic development genes. Using single-cell RNA sequencing, we detected abnormal expansion of homeodomain-only protein homeobox (HOPX)- outer radial glia (oRG), concurrent with a reduction of HOPX+ oRG, in FEZ1-null hCOs. HOPX- oRGs show higher cell mobility as compared to HOPX+ oRGs. Ectopic localization of neuroprogenitors to the outer layer is seen in FEZ1-null hCOs. Anomalous encroachment of TBR2+ intermediate progenitors into CTIP2+ deep layer neurons further indicated abnormalities in cortical layer formation these hCOs. Collectively, our findings highlight the involvement of FEZ1 in early cortical brain development and how it contributes to neurodevelopmental disorders.

Multi-ancestry genome-wide association study of cannabis use disorder yields insight into disease biology and public health implications

As recreational use of cannabis is being decriminalized in many places and medical use widely sanctioned, there are growing concerns about increases in cannabis use disorder (CanUD), which is associated with numerous medical comorbidities. Here we performed a genome-wide association study of CanUD in the Million Veteran Program (MVP), followed by meta-analysis in 1,054,365 individuals (ncases = 64,314) from four broad ancestries designated by the reference panel used for assignment (European n = 886,025, African n = 123,208, admixed American n = 38,289 and East Asian n = 6,843). Population-specific methods were applied to calculate single nucleotide polymorphism-based heritability within each ancestry. Statistically significant single nucleotide polymorphism-based heritability for CanUD was observed in all but the smallest population (East Asian). We discovered genome-wide significant loci unique to each ancestry: 22 in European, 2 each in African and East Asian, and 1 in admixed American ancestries. A genetically informed causal relationship analysis indicated a possible effect of genetic liability for CanUD on lung cancer risk, suggesting potential unanticipated future medical and psychiatric public health consequences that require further study to disentangle from other known risk factors such as cigarette smoking.

Isoform-level transcriptome-wide association uncovers genetic risk mechanisms for neuropsychiatric disorders in the human brain

Methods integrating genetics with transcriptomic reference panels prioritize risk genes and mechanisms at only a fraction of trait-associated genetic loci, due in part to an overreliance on total gene expression as a molecular outcome measure. This challenge is particularly relevant for the brain, in which extensive splicing generates multiple distinct transcript-isoforms per gene. Due to complex correlation structures, isoform-level modeling from cis-window variants requires methodological innovation. Here we introduce isoTWAS, a multivariate, stepwise framework integrating genetics, isoform-level expression and phenotypic associations. Compared to gene-level methods, isoTWAS improves both isoform and gene expression prediction, yielding more testable genes, and increased power for discovery of trait associations within genome-wide association study loci across 15 neuropsychiatric traits. We illustrate multiple isoTWAS associations undetectable at the gene-level, prioritizing isoforms of AKT3, CUL3 and HSPD1 in schizophrenia and PCLO with multiple disorders. Results highlight the importance of incorporating isoform-level resolution within integrative approaches to increase discovery of trait associations, especially for brain-relevant traits.

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.

The molecular genetic landscape of human brain size variation

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

Hofbauer cell function in the term placenta associates with adult cardiovascular and depressive outcomes

Pathological placental inflammation increases the risk for several adult disorders, but these mediators are also expressed under homeostatic conditions, where their contribution to adult health outcomes is unknown. Here we define an inflammation-related expression signature, primarily expressed in Hofbauer cells of the term placenta and use expression quantitative trait loci to create a polygenic score (PGS) predictive of its expression. Using this PGS in the UK Biobank we conduct a phenome-wide association study, followed by Mendelian randomization and identify protective, sex-dependent effects of the placental module on cardiovascular and depressive outcomes. Genes differentially regulated by intra-amniotic infection and preterm birth are over-represented within the module. We also identify aspirin as a putative modulator of this inflammation-related signature. Our data support a model where disruption of placental Hofbauer cell function, due to preterm birth or prenatal infection, contributes to the increased risk of depression and cardiovascular disease observed in these individuals.

An integrative pipeline for circular RNA quantitative trait locus discovery with application in human T cells

MOTIVATION

Molecular quantitative trait locus (QTL) mapping has proven to be a powerful approach for prioritizing genetic regulatory variants and causal genes identified by genome-wide association studies. Recently, this success has been extended to circular RNA (circRNA), a potential group of RNAs that can serve as markers for the diagnosis, prognosis, or therapeutic targets of various human diseases. However, a well-developed computational pipeline for circRNA QTL (circQTL) discovery is still lacking.

RESULTS

We introduce an integrative method for circQTL mapping and implement it as an automated pipeline based on Nextflow, named cscQTL. The proposed method has two main advantages. Firstly, cscQTL improves the specificity by systematically combining outputs of multiple circRNA calling algorithms to obtain highly confident circRNA annotations. Secondly, cscQTL improves the sensitivity by accurately quantifying circRNA expression with the help of pseudo references. Compared to the single method approach, cscQTL effectively identifies circQTLs with an increase of 20%-100% circQTLs detected and recovered all circQTLs that are highly supported by the single method approach. We apply cscQTL to a dataset of human T cells and discover genetic variants that control the expression of 55 circRNAs. By colocalization tests, we further identify circBACH2 and circYY1AP1 as potential candidates for immune disease regulation.

AVAILABILITY AND IMPLEMENTATION

cscQTL is freely available at: https://github.com/datngu/cscQTL and https://doi.org/10.5281/zenodo.7851982.

Genetic Control of Alternative Splicing and its Distinct Role in Colorectal Cancer Mechanisms

BACKGROUND & AIMS

Dysregulation of alternative splicing is implicated in many human diseases, and understanding the genetic variation underlying transcript splicing is essential to dissect the molecular mechanisms of cancers. We aimed to provide a comprehensive functional dissection of splicing quantitative trait loci (sQTLs) in cancer and focus on elucidating its distinct role in colorectal cancer (CRC) mechanisms.

METHODS

We performed a comprehensive sQTL analysis to identify genetic variants that control messenger RNA splicing across 33 cancer types from The Cancer Genome Atlas and independently validated in our 154 CRC tissues. Then, large-scale, multicenter, multi-ethnic case-control studies (34,585 cases and 76,023 controls) were conducted to examine the association of these sQTLs with CRC risk. A series of biological experiments in vitro and in vivo were performed to investigate the potential mechanisms of the candidate sQTLs and target genes.

RESULTS

The molecular characterization of sQTL revealed its distinct role in cancer susceptibility. Tumor-specific sQTL further showed better response to cancer development. In addition, functionally informed polygenic risk score highlighted the potentiality of sQTLs in the CRC prediction. Complemented by large-scale population studies, we identified that the risk allele (T) of a multi-ancestry-associated sQTL rs61746794 significantly increased the risk of CRC in Chinese (odds ratio, 1.20; 95% CI, 1.12-1.29; P = 8.82 × 10-7) and European (odds ratio, 1.11; 95% CI, 1.07-1.16; P = 1.13 × 10-7) populations. rs61746794-T facilitated PRMT7 exon 16 splicing mediated by the RNA-binding protein PRPF8, thus increasing the level of canonical PRMT7 isoform (PRMT7-V2). Overexpression of PRMT7-V2 significantly enhanced the growth of CRC cells and xenograft tumors compared with PRMT7-V1. Mechanistically, PRMT7-V2 functions as an epigenetic writer that catalyzes the arginine methylation of H4R3 and H3R2, subsequently regulating diverse biological processes, including YAP, AKT, and KRAS pathway. A selective PRMT7 inhibitor, SGC3027, exhibited antitumor effects on human CRC cells.

CONCLUSIONS

Our study provides an informative sQTLs resource and insights into the regulatory mechanisms linking splicing variants to cancer risk and serving as biomarkers and therapeutic targets.