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

Genomics yields biological and phenotypic insights into bipolar disorder

Bipolar disorder is a leading contributor to the global burden of disease1. Despite high heritability (60-80%), the majority of the underlying genetic determinants remain unknown2. We analysed data from participants of European, East Asian, African American and Latino ancestries (n = 158,036 cases with bipolar disorder, 2.8 million controls), combining clinical, community and self-reported samples. We identified 298 genome-wide significant loci in the multi-ancestry meta-analysis, a fourfold increase over previous findings3, and identified an ancestry-specific association in the East Asian cohort. Integrating results from fine-mapping and other variant-to-gene mapping approaches identified 36 credible genes in the aetiology of bipolar disorder. Genes prioritized through fine-mapping were enriched for ultra-rare damaging missense and protein-truncating variations in cases with bipolar disorder4, highlighting convergence of common and rare variant signals. We report differences in the genetic architecture of bipolar disorder depending on the source of patient ascertainment and on bipolar disorder subtype (type I or type II). Several analyses implicate specific cell types in the pathophysiology of bipolar disorder, including GABAergic interneurons and medium spiny neurons. Together, these analyses provide additional insights into the genetic architecture and biological underpinnings of bipolar disorder.

Quantitative trait loci mapping of circulating metabolites in cerebrospinal fluid to uncover biological mechanisms involved in brain-related phenotypes

Genomic studies of molecular traits have provided mechanistic insights into complex disease, though these lag behind for brain-related traits due to the inaccessibility of brain tissue. We leveraged cerebrospinal fluid (CSF) to study neurobiological mechanisms in vivo, measuring 5543 CSF metabolites, the largest panel in CSF to date, in 977 individuals of European ancestry. Individuals originated from two separate cohorts including cognitively healthy subjects (n = 490) and a well-characterized memory clinic sample, the Amsterdam Dementia Cohort (ADC, n = 487). We performed metabolite quantitative trait loci (mQTL) mapping on CSF metabolomics and found 126 significant mQTLs, representing 65 unique CSF metabolites across 51 independent loci. To better understand the role of CSF mQTLs in brain-related disorders we integrated our CSF mQTL results with pre-existing summary statistics on brain traits, identifying 34 genetic associations between CSF metabolites and brain traits. Over 90% of significant mQTLs demonstrated colocalized associations with brain-specific gene expression, unveiling potential neurobiological pathways.

Dissecting Schizophrenia Biology Using Pleiotropy with Cognitive Genomics

BACKGROUND

Given the increasingly large number of loci discovered by psychiatric GWAS, specification of the key biological pathways underlying these loci has become a priority for the field. We have previously leveraged the pleiotropic genetic relationships between schizophrenia and two cognitive phenotypes (educational attainment and cognitive task performance) to differentiate two subsets of illness-relevant SNPs: (1) those with "concordant" alleles, which are associated with reduced cognitive performance and educational attainment and increased schizophrenia risk; and (2) those with "discordant" alleles linked to reduced educational and/or cognitive levels but lower schizophrenia susceptibility.

METHODS

In the present study, we extend our prior work, utilizing larger input GWAS datasets and a more powerful statistical approach to pleiotropic meta-analysis, the Pleiotropic Locus Exploration and Interpretation using Optimal test (PLEIO).

RESULTS

Our pleiotropic meta-analysis of schizophrenia and the two cognitive phenotypes revealed 768 significant loci (166 novel). Among these, 347 loci harbored concordant SNPs, 270 encompassed discordant SNPs, and 151 "dual" loci contained concordant and discordant SNPs. Competitive gene-set analysis using MAGMA linked concordant SNP loci with neurodevelopmental pathways (e.g., neurogenesis), whereas discordant loci were associated with mature neuronal synaptic functions. These distinctions were also observed in BrainSpan analysis of temporal enrichment patterns across developmental periods, with concordant loci containing more prenatally expressed genes than discordant loci. Dual loci were enriched for genes related to mRNA translation initiation, representing a novel finding in the schizophrenia literature.

CONCLUSIONS

Pleiotropic analysis permits not only enhanced statistical power for locus discovery, but also the ability to parse distinct biological processes associated with endophenotypes.

webTWAS 2.0: update platform for identifying complex disease susceptibility genes through transcriptome-wide association study

Transcriptome-wide association study (TWAS) has successfully identified numerous complex disease susceptibility genes in the post-genome-wide association study (GWAS) era. Over the past 3 years, the focus of TWAS algorithms has shifted from merely identifying associations to understanding how single nucleotide polymorphisms (SNPs) regulate gene expression, with a growing emphasis on incorporating fine-mapping techniques. Additionally, the rapid increase in GWAS summary statistics, driven largely by the UK Biobank and other consortia, has made it essential to update our webTWAS resource. To address these challenges and meet the growing needs of researchers, we developed webTWAS 2.0, an updated platform for identifying susceptibility genes for human complex diseases using TWAS. Additionally, webTWAS 2.0 provides an online TWAS analysis tool that simplifies conducting TWAS analyses. The updated resource includes 7247 GWAS summary statistics covering 1588 complex human diseases from 192 publications. It also incorporates multiple TWAS methods, such as sTF-TWAS, 3'aTWAS and GIFT, along with an updated interactive visualization tool that allows users to easily explore significant associations across different methods. Other upgrades include a personalized online analysis tool for user-submitted GWAS data and a refined search function that makes it easier to identify relevant associations and meet diverse user needs more efficiently. webTWAS 2.0 is freely accessible at http://www.webtwas.net.

Isoform-level expression of the constitutive androstane receptor (CAR or NR1I3) transcription factor better predicts the mRNA expression of the cytochrome P450s in human liver samples

Many factors cause interperson variability in the activity and expression of the cytochrome P450 (CYP) drug-metabolizing enzymes in the liver, leading to variable drug exposure and treatment outcomes. Several liver-enriched transcription factors are associated with CYP expression, with estrogen receptor α (ESR1) and constitutive androstane receptor (CAR or NR1I3) being the 2 top factors. ESR1 and NR1I3 undergo extensive alternative splicing that results in numerous splice isoforms, but how these splice isoforms associate with CYP expression is unknown. Here, we quantified 18 NR1I3 splice isoforms and the 3 most abundant ESR1 isoforms in 260 liver samples derived from African Americans (n = 125) and European Americans (n = 135). Our results showed variable splice isoform populations in the liver for both NR1I3 and ESR1. Multiple linear regression analyses revealed that compared with gene-level NR1I3, isoform-level NR1I3 expression better predicted the mRNA expression of most CYPs and 3 UDP-glucuronosyltransferases (UGTs), whereas ESR1 isoforms improved predictive models for the UGTs and CYP2D6 but not for most CYPs. Also, different NR1I3 isoforms were associated with different CYPs, and the associations varied depending on sample ancestry. Surprisingly, noncanonical NR1I3 isoforms having retained introns (introns 2 or 6) were abundantly expressed and associated with the expression of most CYPs and UGTs, whereas the reference isoform (NR1I3-205) was only associated with CYP2D6. Moreover, NR1I3 isoform diversity increased during the differentiation of induced pluripotent stem cells to hepatocytes, paralleling increasing CYP expression. These results suggest that isoform-level transcription factor expression may help to explain variation in CYP or UGT expression between individuals. SIGNIFICANCE STATEMENT: We quantified 18 NR1I3 splice isoforms and the 3 most abundant ESR1 splice isoforms in 260 liver samples derived from African American and European American donors and found variable NR1I3 and ESR1 splice isoform expression in the liver. Multiple linear regression analysis showed that, compared with gene-level expression, isoform-level expression of NR1I3 and ESR1 better predicted the mRNA expression of some cytochrome P450s and UDP-glucuronosyltransferases, highlighting the importance of isoform-level analyses to enhance our understanding of gene transcriptional regulatory networks controlling the expression of drug-metabolizing enzymes.

Genetic analysis of psychosis Biotypes: shared Ancestry-adjusted polygenic risk and unique genomic associations

The Bipolar-Schizophrenia Network for Intermediate Phenotypes (B-SNIP) created psychosis Biotypes based on neurobiological measurements in a multi-ancestry sample. These Biotypes cut across DSM diagnoses of schizophrenia, schizoaffective disorder, and bipolar disorder with psychosis. Two recently developed post hoc ancestry adjustment methods of Polygenic Risk Scores (PRSs) generate Ancestry-Adjusted PRSs (AAPRSs), which allow for PRS analysis of multi-ancestry samples. Applied to schizophrenia PRS, we found the Khera AAPRS method to show superior portability and comparable prediction accuracy as compared with the Ge method. The three Biotypes of psychosis disorders had similar AAPRSs across ancestries. In genomic analysis of Biotypes, 12 genes, and isoforms showed significant genomic associations with specific Biotypes in a Transcriptome-Wide Association Study (TWAS) of genetically regulated expression (GReX) in the adult brain and fetal brain. TWAS inflation was addressed by the inclusion of genotype principal components in the association analyses. Seven of these 12 genes/isoforms satisfied Mendelian Randomization (MR) criteria for putative causality, including four genes TMEM140, ARTN, C1orf115, CYREN, and three transcripts ENSG00000272941, ENSG00000257176, ENSG00000287733. These genes are enriched in the biological pathways of Rearranged during Transfection (RET) signaling, Neural Cell Adhesion Molecule 1 (NCAM1) interactions, and NCAM signaling for neurite out-growth. The specific associations with Biotypes suggest that pharmacological clinical trials and biological investigations might benefit from analyzing Biotypes separately.

Co-expression-wide association studies link genetically regulated interactions with complex traits

Transcriptome- and proteome-wide association studies (TWAS/PWAS) have proven successful in prioritizing genes and proteins whose genetically regulated expression modulates disease risk, but they ignore potential co-expression and interaction effects. To address this limitation, we introduce the co-expression-wide association study (COWAS) method, which can identify pairs of genes or proteins whose genetically regulated co-expression is associated with complex traits. COWAS first trains models to predict expression and co-expression conditional on genetic variation, and then tests for association between imputed co-expression and the trait of interest while also accounting for direct effects from each exposure. We applied our method to plasma proteomic concentrations from the UK Biobank, identifying dozens of interacting protein pairs associated with cholesterol levels, Alzheimer's disease, and Parkinson's disease. Notably, our results demonstrate that co-expression between proteins may affect complex traits even if neither protein is detected to influence the trait when considered on its own. We also show how COWAS can help disentangle direct and interaction effects, providing a richer picture of the molecular networks that mediate genetic effects on disease outcomes.

Multimodal analysis of RNA sequencing data powers discovery of complex trait genetics

RNA sequencing has the potential to reveal many modalities of transcriptional regulation, such as various splicing phenotypes, but studies on gene regulation are often limited to gene expression due to the complexity of extracting and analyzing multiple RNA phenotypes. Here, we present Pantry, a framework to efficiently generate diverse RNA phenotypes from RNA sequencing data and perform downstream integrative analyses with genetic data. Pantry generates phenotypes from six modalities of transcriptional regulation (gene expression, isoform ratios, splice junction usage, alternative TSS/polyA usage, and RNA stability) and integrates them with genetic data via QTL mapping, TWAS, and colocalization testing. We apply Pantry to Geuvadis and GTEx data, finding that 4768 of the genes with no identified eQTL in Geuvadis have QTL in at least one other transcriptional modality, resulting in a 66% increase in genes over eQTL mapping. We further found that the QTL exhibit modality-specific functional properties that are further reinforced by joint analysis of different RNA modalities. We also show that generalizing TWAS to multiple RNA modalities approximately doubles the discovery of unique gene-trait associations, and enhances identification of regulatory mechanisms underlying GWAS signal in 42% of previously associated gene-trait pairs.

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.

Isoform-level analyses of 6 cancers uncover extensive genetic risk mechanisms undetected at the gene-level

Integrating genome-wide association study (GWAS) and transcriptomic datasets can help identify potential mediators for germline genetic risk of cancer. However, traditional methods have been largely unsuccessful because of an overreliance on total gene expression. These approaches overlook alternative splicing, which can produce multiple isoforms from the same gene, each with potentially different effects on cancer risk. Here, we integrate genetic and multi-tissue isoform-level gene expression data from the Genotype Tissue-Expression Project (GTEx, N = 108-574) with publicly available European-ancestry GWAS summary statistics (all N > 20,000 cases) to identify both isoform- and gene-level risk associations with six cancers (breast, endometrial, colorectal, lung, ovarian, prostate) and six related cancer subtype classifications (N = 12 total). Compared to traditional methods leveraging total gene expression, directly modeling isoform expression through transcriptome-wide association studies (isoTWAS) substantially increases discovery of transcriptomic mechanisms underlying genetic associations. Using the same RNA-seq datasets, isoTWAS identified 164% more significant unique gene associations compared to TWAS (6,163 and 2,336, respectively), with isoTWAS-prioritized genes enriched 4-fold for evolutionarily-constrained genes (P = 6.1 × 10-13). isoTWAS tags transcriptomic associations at 52% more independent GWAS loci compared to TWAS across the six cancers. Additionally, isoform expression mediates an estimated 63% greater proportion of cancer risk SNP heritability compared to gene expression when evaluating cis-genetic influence on isoform expression. We highlight several notable isoTWAS associations that demonstrate GWAS colocalization at the isoform level but not at the gene level, including, CLPTM1L (lung cancer), LAMC1 (colorectal), and BABAM1 (breast). These results underscore the critical importance of modeling isoform-level expression to maximize discovery of genetic risk mechanisms for cancers.

A Deep Dive into Statistical Modeling of RNA Splicing QTLs Reveals New Variants that Explain Neurodegenerative Disease

Genome-wide association studies (GWAS) have identified thousands of putative disease causing variants with unknown regulatory effects. Efforts to connect these variants with splicing quantitative trait loci (sQTLs) have provided functional insights, yet sQTLs reported by existing methods cannot explain many GWAS signals. We show current sQTL modeling approaches can be improved by considering alternative splicing representation, model calibration, and covariate integration. We then introduce MAJIQTL, a new pipeline for sQTL discovery. MAJIQTL includes two new statistical methods: a weighted multiple testing approach for sGene discovery and a model for sQTL effect size inference to improve variant prioritization. By applying MAJIQTL to GTEx, we find significantly more sGenes harboring sQTLs with functional significance. Notably, our analysis implicates the novel variant rs582283 in Alzheimer's disease. Using antisense oligonucleotides, we validate this variant's effect by blocking the implicated YBX3 binding site, leading to exon skipping in the gene MS4A3.

Dynamic convergence of autism disorder risk genes across neurodevelopment

Over a hundred risk genes underlie risk for autism spectrum disorder (ASD) but the extent to which they converge on shared downstream targets to increase ASD risk is unknown. To test the hypothesis that cellular context impacts the nature of convergence, here we apply a pooled CRISPR approach to target 29 ASD loss-of-function genes in human induced pluripotent stem cell (hiPSC)-derived neural progenitor cells, glutamatergic neurons, and GABAergic neurons. Two distinct approaches (gene-level and network-level analyses) demonstrate that convergence is greatest in mature glutamatergic neurons. Convergent effects are dynamic, varying in strength, composition, and biological role between cell types, increasing with functional similarity of the ASD genes examined, and driven by cell-type-specific gene co-expression patterns. Stratification of ASD genes yield targeted drug predictions capable of reversing gene-specific convergent signatures in human cells and ASD-related behaviors in zebrafish. Altogether, convergent networks downstream of ASD risk genes represent novel points of individualized therapeutic intervention.

Neuron type-specific proteomics reveals distinct Shank3 proteoforms in iSPNs and dSPNs lead to striatal synaptopathy in Shank3B-/- mice

Combinatorial expression of postsynaptic proteins underlies synapse diversity within and between neuron types. Thus, characterization of neuron-type-specific postsynaptic proteomes is key to obtaining a deeper understanding of discrete synaptic properties and how selective dysfunction manifests in synaptopathies. To overcome the limitations associated with bulk measures of synaptic protein abundance, we developed a biotin proximity protein tagging probe to characterize neuron-type-specific postsynaptic proteomes in vivo. We found Shank3 protein isoforms are differentially expressed by direct and indirect pathway spiny projection neurons (dSPNs and iSPNs). Investigation of Shank3B-/- mice lacking exons 13-16 within the Shank3 gene, reveal distinct Shank3 protein isoform expression in iSPNs and dSPNs. In Shank3B-/- striatum, Shank3E and Shank3NT are expressed by dSPNs but are undetectable in iSPNs. Proteomic analysis indicates significant and selective alterations in the postsynaptic proteome of Shank3B-/- iSPNs. Correspondingly, the deletion of exons 13-16 diminishes dendritic spine density, reduces spine head diameter, and hampers corticostriatal synaptic transmission in iSPNs. Remarkably, reintroducing Shank3E in adult Shank3B-/- iSPNs significantly rectifies the observed dendritic spine morphological and corticostriatal synaptic transmission deficits. We report unexpected cell-type specific synaptic protein isoform expression which could play a key causal role in specifying synapse diversity and selective synapse dysfunction in synaptopathies.

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.

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.

Identifying drug targets for schizophrenia through gene prioritization

Background

Schizophrenia genome-wide association studies (GWASes) have identified >250 significant loci and prioritized >100 disease-related genes. However, gene prioritization efforts have mostly been restricted to locus-based methods that ignore information from the rest of the genome.

Methods

To more accurately characterize genes involved in schizophrenia etiology, we applied a combination of highly-predictive tools to a published GWAS of 67,390 schizophrenia cases and 94,015 controls. We combined both locus-based methods (fine-mapped coding variants, distance to GWAS signals) and genome-wide methods (PoPS, MAGMA, ultra-rare coding variant burden tests). To validate our findings, we compared them with previous prioritization efforts, known neurodevelopmental genes, and results from the PsyOPS tool.

Results

We prioritized 62 schizophrenia genes, 41 of which were also highlighted by our validation methods. In addition to DRD2, the principal target of antipsychotics, we prioritized 9 genes that are targeted by approved or investigational drugs. These included drugs targeting glutamatergic receptors (GRIN2A and GRM3), calcium channels (CACNA1C and CACNB2), and GABAB receptor (GABBR2). These also included genes in loci that are shared with an addiction GWAS (e.g. PDE4B and VRK2).

Conclusions

We curated a high-quality list of 62 genes that likely play a role in the development of schizophrenia. Developing or repurposing drugs that target these genes may lead to a new generation of schizophrenia therapies. Rodent models of addiction more closely resemble the human disorder than rodent models of schizophrenia. As such, genes prioritized for both disorders could be explored in rodent addiction models, potentially facilitating drug development.

Multimodal analysis of RNA sequencing data powers discovery of complex trait genetics

Transcriptome data is commonly used to understand genome function via quantitative trait loci (QTL) mapping and to identify the molecular mechanisms driving genome wide association study (GWAS) signals through colocalization analysis and transcriptome-wide association studies (TWAS). While RNA sequencing (RNA-seq) has the potential to reveal many modalities of transcriptional regulation, such as various splicing phenotypes, such studies are often limited to gene expression due to the complexity of extracting and analyzing multiple RNA phenotypes. Here, we present Pantry (Pan-transcriptomic phenotyping), a framework to efficiently generate diverse RNA phenotypes from RNA-seq data and perform downstream integrative analyses with genetic data. Pantry currently generates phenotypes from six modalities of transcriptional regulation (gene expression, isoform ratios, splice junction usage, alternative TSS/polyA usage, and RNA stability) and integrates them with genetic data via QTL mapping, TWAS, and colocalization testing. We applied Pantry to Geuvadis and GTEx data, and found that 4,768 of the genes with no identified expression QTL in Geuvadis had QTLs in at least one other transcriptional modality, resulting in a 66% increase in genes over expression QTL mapping. We further found that QTLs exhibit modality-specific functional properties that are further reinforced by joint analysis of different RNA modalities. We also show that generalizing TWAS to multiple RNA modalities (xTWAS) approximately doubles the discovery of unique gene-trait associations, and enhances identification of regulatory mechanisms underlying GWAS signal in 42% of previously associated gene-trait pairs. We provide the Pantry code, RNA phenotypes from all Geuvadis and GTEx samples, and xQTL and xTWAS results on the web.

Alternative splicing in prostate cancer progression and therapeutic resistance

Prostate cancer (CaP) remains the second leading cause of cancer deaths in western men. CaP mortality results from diverse molecular mechanisms that mediate resistance to the standard of care treatments for metastatic disease. Recently, alternative splicing has been recognized as a hallmark of CaP aggressiveness. Alternative splicing events cause treatment resistance and aggressive CaP behavior and are determinants of the emergence of the two major types of late-stage treatment-resistant CaP, namely castration-resistant CaP (CRPC) and neuroendocrine CaP (NEPC). Here, we review recent multi-omics data that are uncovering the complicated landscape of alternative splicing events during CaP progression and the impact that different gene transcript isoforms can have on CaP cell biology and behavior. We discuss renewed insights in the molecular machinery by which alternative splicing occurs and contributes to the failure of systemic CaP therapies. The potential for alternative splicing events to serve as diagnostic markers and/or therapeutic targets is explored. We conclude by considering current challenges and promises associated with splicing-modulating therapies, and their potential for clinical translation into CaP patient care.

Dissecting Schizophrenia Biology Using Pleiotropy with Cognitive Genomics

Given the increasingly large number of loci discovered by psychiatric GWAS, specification of the key biological pathways underlying these loci has become a priority for the field. We have previously leveraged the pleiotropic genetic relationships between schizophrenia and two cognitive phenotypes (educational attainment and cognitive task performance) to differentiate two subsets of illness-relevant SNPs: (1) those with "concordant" alleles, which are associated with reduced cognitive ability/education and increased schizophrenia risk; and (2) those with "discordant" alleles linked to reduced educational and/or cognitive levels but lower schizophrenia susceptibility. In the present study, we extend our prior work, utilizing larger input GWAS datasets and a more powerful statistical approach to pleiotropic meta-analysis, the Pleiotropic Locus Exploration and Interpretation using Optimal test (PLEIO). Our pleiotropic meta-analysis of schizophrenia and the two cognitive phenotypes revealed 768 significant loci (159 novel). Among these, 347 loci harbored concordant SNPs, 270 encompassed discordant SNPs, and 151 "dual" loci contained concordant and discordant SNPs. Competitive gene-set analysis using MAGMA related concordant SNP loci with neurodevelopmental pathways (e.g., neurogenesis), whereas discordant loci were associated with mature neuronal synaptic functions. These distinctions were also observed in BrainSpan analysis of temporal enrichment patterns across developmental periods, with concordant loci containing more prenatally expressed genes than discordant loci. Dual loci were enriched for genes related to mRNA translation initiation, representing a novel finding in the schizophrenia literature.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Transcriptome-Wide Association Studies (TWAS): Methodologies, Applications, and Challenges

Transcriptome-wide association study (TWAS) methodologies aim to identify genetic effects on phenotypes through the mediation of gene transcription. In TWAS, in silico models of gene expression are trained as functions of genetic variants and then applied to genome-wide association study (GWAS) data. This post-GWAS analysis identifies gene-trait associations with high interpretability, enabling follow-up functional genomics studies and the development of genetics-anchored resources. We provide an overview of commonly used TWAS approaches, their advantages and limitations, and some widely used applications. © 2024 Wiley Periodicals LLC.

Accounting for isoform expression increases power to identify genetic regulation of gene expression

A core problem in genetics is molecular quantitative trait locus (QTL) mapping, in which genetic variants associated with changes in the molecular phenotypes are identified. One of the most-studied molecular QTL mapping problems is expression QTL (eQTL) mapping, in which the molecular phenotype is gene expression. It is common in eQTL mapping to compute gene expression by aggregating the expression levels of individual isoforms from the same gene and then performing linear regression between SNPs and this aggregated gene expression level. However, SNPs may regulate isoforms from the same gene in different directions due to alternative splicing, or only regulate the expression level of one isoform, causing this approach to lose power. Here, we examine a broader question: which genes have at least one isoform whose expression level is regulated by genetic variants? In this study, we propose and evaluate several approaches to answering this question, demonstrating that "isoform-aware" methods-those that account for the expression levels of individual isoforms-have substantially greater power to answer this question than standard "gene-level" eQTL mapping methods. We identify settings in which different approaches yield an inflated number of false discoveries or lose power. In particular, we show that calling an eGene if there is a significant association between a SNP and any isoform fails to control False Discovery Rate, even when applying standard False Discovery Rate correction. We show that similar trends are observed in real data from the GEUVADIS and GTEx studies, suggesting the possibility that similar effects are present in these consortia.