On the interpretation of transcriptome-wide association studies.

Transcriptome-wide association studies (TWAS) aim to detect relationships between gene expression and a phenotype, and are commonly used for secondary analysis of genome-wide association study (GWAS) results. Results from TWAS analyses are often interpreted as indicating a genetic relationship between gene expression and a phenotype, but this interpretation is not consistent with the null hypothesis that is evaluated in the traditional TWAS framework. In this study we provide a mathematical outline of this TWAS framework, and elucidate what interpretations are warranted given the null hypothesis it actually tests. We then use both simulations and real data analysis to assess the implications of misinterpreting TWAS results as indicative of a genetic relationship between gene expression and the phenotype. Our simulation results show considerably inflated type 1 error rates for TWAS when interpreted this way, with 41% of significant TWAS associations detected in the real data analysis found to have insufficient statistical evidence to infer such a relationship. This demonstrates that in current implementations, TWAS cannot reliably be used to investigate genetic relationships between gene expression and a phenotype, but that local genetic correlation analysis can serve as a potential alternative.

Refining the scope of genetic influences on alcohol misuse through environmental stratification and gene-environment interaction.

BACKGROUND: Gene-environment interaction (G × E) is likely an important influence shaping individual differences in alcohol misuse (AM), yet it has not been extensively studied in molecular genetic research. In this study, we use a series of genome-wide gene-environment interaction (GWEIS) and in silico annotation methods with the aim of improving gene identification and biological understanding of AM. METHODS: We carried out GWEIS for four AM phenotypes in the large UK Biobank sample (N = 360,314), with trauma exposure and socioeconomic status (SES) as moderators of the genetic effects. Exploratory analyses compared stratified genome-wide association (GWAS) and GWEIS modeling approaches. We applied functional annotation, gene- and gene-set enrichment, and polygenic score analyses to interpret the GWEIS results. RESULTS: GWEIS models showed few genetic variants with significant interaction effects across gene-environment pairs. Enrichment analyses identified moderation by SES of the genes NOXA1, DLGAP1, and UBE2L3 on drinking quantity and the gene IFIT1B on drinking frequency. Except for DLGAP1, these genes have not previously been linked to AM. The most robust results (GWEIS interaction p = 4.59e-09) were seen for SES moderating the effects of variants linked to immune-related genes on a pattern of drinking with versus without meals. CONCLUSIONS: Our results highlight several genes and a potential mechanism of immune system functioning behind the moderating effect of SES on the genetic influences on AM. Although GWEIS seems to be a preferred approach over stratified GWAS, modeling G × E effects at the molecular level remains a challenge even in large samples. Understanding these effects will require substantial effort and more in-depth phenotypic measurement.

Local patterns of genetic sharing challenge the boundaries between neuropsychiatric and insulin resistance-related conditions.

The co-occurrence of insulin resistance (IR)-related metabolic conditions with neuropsychiatric disorders is a complex public health challenge. Evidence of the genetic links between these phenotypes is emerging, but little is currently known about the genomic regions and biological functions that are involved. To address this, we performed Local Analysis of [co]Variant Association (LAVA) using large-scale (N=9,725-933,970) genome-wide association studies (GWASs) results for three IR-related conditions (type 2 diabetes mellitus, obesity, and metabolic syndrome) and nine neuropsychiatric disorders. Subsequently, positional and expression quantitative trait locus (eQTL)-based gene mapping and downstream functional genomic analyses were performed on the significant loci. Patterns of negative and positive local genetic correlations (|rg|=0.21-1, pFDR<0.05) were identified at 109 unique genomic regions across all phenotype pairs. Local correlations emerged even in the absence of global genetic correlations between IR-related conditions and Alzheimer's disease, bipolar disorder, and Tourette's syndrome. Genes mapped to the correlated regions showed enrichment in biological pathways integral to immune-inflammatory function, vesicle trafficking, insulin signalling, oxygen transport, and lipid metabolism. Colocalisation analyses further prioritised 10 genetically correlated regions for likely harbouring shared causal variants, displaying high deleterious or regulatory potential. These variants were found within or in close proximity to genes, such as SLC39A8 and HLA-DRB1, that can be targeted by supplements and already known drugs, including omega-3/6 fatty acids, immunomodulatory, antihypertensive, and cholesterol-lowering drugs. Overall, our findings underscore the complex genetic landscape of IR-neuropsychiatric multimorbidity, advocating for an integrated disease model and offering novel insights for research and treatment strategies in this domain.

The Genetic Architectures of Functional and Structural Connectivity Properties within Cerebral Resting-State Networks.

Functional connectivity within resting-state networks (RSN-FC) is vital for cognitive functioning. RSN-FC is heritable and partially translates to the anatomic architecture of white matter, but the genetic component of structural connections of RSNs (RSN-SC) and their potential genetic overlap with RSN-FC remain unknown. Here, we perform genome-wide association studies (N discovery = 24,336; N replication = 3412) and annotation on RSN-SC and RSN-FC. We identify genes for visual network-SC that are involved in axon guidance and synaptic functioning. Genetic variation in RSN-FC impacts biological processes relevant to brain disorders that previously were only phenotypically associated with RSN-FC alterations. Correlations of the genetic components of RSNs are mostly observed within the functional domain, whereas less overlap is observed within the structural domain and between the functional and structural domains. This study advances the understanding of the complex functional organization of the brain and its structural underpinnings from a genetics viewpoint.

An integrated framework for local genetic correlation analysis.

Genetic correlation (rg) analysis is used to identify phenotypes that may have a shared genetic basis. Traditionally, rg is studied globally, considering only the average of the shared signal across the genome, although this approach may fail when the rg is confined to particular genomic regions or in opposing directions at different loci. Current tools for local rg analysis are restricted to analysis of two phenotypes. Here we introduce LAVA, an integrated framework for local rg analysis that, in addition to testing the standard bivariate local rgs between two phenotypes, can evaluate local heritabilities and analyze conditional genetic relations between several phenotypes using partial correlation and multiple regression. Applied to 25 behavioral and health phenotypes, we show considerable heterogeneity in the bivariate local rgs across the genome, which is often masked by the global rg patterns, and demonstrate how our conditional approaches can elucidate more complex, multivariate genetic relations.

Exploring the genetic overlap between twelve psychiatric disorders.

The widespread comorbidity among psychiatric disorders demonstrated in epidemiological studies1-5 is mirrored by non-zero, positive genetic correlations from large-scale genetic studies6-10. To identify shared biological processes underpinning this observed phenotypic and genetic covariance and enhance molecular characterization of general psychiatric disorder liability11-13, we used several strategies aimed at uncovering pleiotropic, that is, cross-trait-associated, single-nucleotide polymorphisms (SNPs), genes and biological pathways. We conducted cross-trait meta-analysis on 12 psychiatric disorders to identify pleiotropic SNPs. The meta-analytic signal was driven by schizophrenia, hampering interpretation and joint biological characterization of the cross-trait meta-analytic signal. Subsequent pairwise comparisons of psychiatric disorders identified substantial pleiotropic overlap, but mainly among pairs of psychiatric disorders, and mainly at less stringent P-value thresholds. Only annotations related to evolutionarily conserved genomic regions were significant for multiple (9 out of 12) psychiatric disorders. Overall, identification of shared biological mechanisms remains challenging due to variation in power and genetic architecture between psychiatric disorders.

Genome-wide gene-environment interactions in neuroticism: an exploratory study across 25 environments.

Gene-environment interactions (GxE) are often suggested to play an important role in the aetiology of psychiatric phenotypes, yet so far, only a handful of genome-wide environment interaction studies (GWEIS) of psychiatric phenotypes have been conducted. Representing the most comprehensive effort of its kind to date, we used data from the UK Biobank to perform a series of GWEIS for neuroticism across 25 broadly conceptualised environmental risk factors (trauma, social support, drug use, physical health). We investigated interactions on the level of SNPs, genes, and gene-sets, and computed interaction-based polygenic risk scores (PRS) to predict neuroticism in an independent sample subset (N = 10,000). We found that the predictive ability of the interaction-based PRSs did not significantly improve beyond that of a traditional PRS based on SNP main effects from GWAS, but detected one variant and two gene-sets showing significant interaction signal after correction for the number of analysed environments. This study illustrates the possibilities and limitations of a comprehensive GWEIS in currently available sample sizes.

Systematic assessment of variability in the proteome of iPSC derivatives.

The use of induced pluripotent stem cells (iPSC) to model human complex diseases is gaining popularity as it allows investigation of human cells that are otherwise sparsely available. However, due to its laborious and cost intensive nature, iPSC research is often plagued by limited sample size and putative large variability between clones, decreasing statistical power for detecting experimental effects. Here, we investigate the source and magnitude of variability in the proteome of parallel differentiated astrocytes using mass spectrometry. We compare three possible sources of variability: inter-donor variability, inter- and intra-clonal variability, at different stages of maturation. We show that the interclonal variability is significantly smaller than the inter-donor variability, and that including more donors has a much larger influence on statistical power than adding more clones per donor. Our results provide insight into the sources of variability at protein level between iPSC samples derived in parallel and will aid in optimizing iPSC studies.