Sherlock: detecting gene-disease associations by matching patterns of expression QTL and GWAS

From GWASs toward Mechanistic Understanding with Case Studies in Dermatogenetics

Many human skin diseases result from the complex interplay of genetic and environmental mechanisms that are largely unknown. GWASs have yielded insight into the genetic aspect of complex disease by highlighting regions of the genome or specific genetic variants associated with disease. Leveraging this information to identify causal genes and cell types will provide insight into fundamental biology, inform diagnostics, and aid drug discovery. However, the etiological mechanisms from genetic variant to disease are still unestablished in most cases. There now exists an unprecedented wealth of data and computational methods for variant interpretation in a functional context. It can be challenging to decide where to start owing to a lack of consensus on the best way to identify causal genetic mechanisms. This article highlights 3 key aspects of genetic variant interpretation: prioritizing causal genes, cell types, and pathways. We provide a practical overview of the main methods and datasets, giving examples from recent atopic dermatitis studies to provide a blueprint for variant interpretation. A collection of resources, including brief description and links to the packages and web tools, is provided for researchers looking to start in silico follow-up genetic analysis of associated genetic variants.

An overview of detecting gene-trait associations by integrating GWAS summary statistics and eQTLs

Detecting genes that affect specific traits (such as human diseases and crop yields) is important for treating complex diseases and improving crop quality. A genome-wide association study (GWAS) provides new insights and directions for understanding complex traits by identifying important single nucleotide polymorphisms. Many GWAS summary statistics data related to various complex traits have been gathered recently. Studies have shown that GWAS risk loci and expression quantitative trait loci (eQTLs) often have a lot of overlaps, which makes gene expression gradually become an important intermediary to reveal the regulatory role of GWAS. In this review, we review three types of gene-trait association detection methods of integrating GWAS summary statistics and eQTLs data, namely colocalization methods, transcriptome-wide association study-oriented approaches, and Mendelian randomization-related methods. At the theoretical level, we discussed the differences, relationships, advantages, and disadvantages of various algorithms in the three kinds of gene-trait association detection methods. To further discuss the performance of various methods, we summarize the significant gene sets that influence high-density lipoprotein, low-density lipoprotein, total cholesterol, and triglyceride reported in 16 studies. We discuss the performance of various algorithms using the datasets of the four lipid traits. The advantages and limitations of various algorithms are analyzed based on experimental results, and we suggest directions for follow-up studies on detecting gene-trait associations.

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.

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.

Joint-tissue integrative analysis identifies high-risk genes for Parkinson's disease

The loss of dopaminergic neurons in the substantia nigra and the abnormal accumulation of synuclein proteins and neurotransmitters in Lewy bodies constitute the primary symptoms of Parkinson's disease (PD). Besides environmental factors, scholars are in the early stages of comprehending the genetic factors involved in the pathogenic mechanism of PD. Although genome-wide association studies (GWAS) have unveiled numerous genetic variants associated with PD, precisely pinpointing the causal variants remains challenging due to strong linkage disequilibrium (LD) among them. Addressing this issue, expression quantitative trait locus (eQTL) cohorts were employed in a transcriptome-wide association study (TWAS) to infer the genetic correlation between gene expression and a particular trait. Utilizing the TWAS theory alongside the enhanced Joint-Tissue Imputation (JTI) technique and Mendelian Randomization (MR) framework (MR-JTI), we identified a total of 159 PD-associated genes by amalgamating LD score, GTEx eQTL data, and GWAS summary statistic data from a substantial cohort. Subsequently, Fisher's exact test was conducted on these PD-associated genes using 5,152 differentially expressed genes sourced from 12 PD-related datasets. Ultimately, 29 highly credible PD-associated genes, including CTX1B, SCNA, and ARSA, were uncovered. Furthermore, GO and KEGG enrichment analyses indicated that these genes primarily function in tissue synthesis, regulation of neuron projection development, vesicle organization and transportation, and lysosomal impact. The potential PD-associated genes identified in this study not only offer fresh insights into the disease's pathophysiology but also suggest potential biomarkers for early disease detection.

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.

Trans -eQTL hotspots shape complex traits by modulating cellular states

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

Alternative splicing in mouse brains affected by psychological stress is enriched in the signaling, neural transmission and blood-brain barrier pathways

Psychological stress increases the risk of major psychiatric disorders. Psychological stress on mice was reported to induce differential gene expression (DEG) in mice brain regions. Alternative splicing is a fundamental aspect of gene expression and has been associated with psychiatric disorders but has not been investigated in the stressed brain yet. This study investigated changes in gene expression and splicing under psychological stress, the related pathways, and possible relationship with psychiatric disorders. RNA-seq raw data of 164 mouse brain samples from 3 independent datasets with stressors including chronic social defeat stress (CSDS), early life stress (ELS), and two-hit stress of combined CSDS and ELS were collected. There were more changes in splicing than in gene expression in the ventral hippocampus and medial prefrontal cortex, but stress-induced changes of individual genes by differential splicing and differential expression could not be replicated. In contrast, pathway analyses produced robust findings: stress-induced differentially spliced genes (DSGs) were reproducibly enriched in neural transmission and blood-brain barrier systems, and DEGs were reproducibly enriched in stress response-related functions. The hub genes of DSG-related PPI networks were enriched in synaptic functions. The corresponding human homologs of stress-induced DSGs were robustly enriched in AD-related DSGs as well as BD and SCZ in GWAS. These results suggested that stress-induced DSGs from different datasets belong to the same biological system throughout the stress response process, resulting in consistent stress response effects.

Harnessing Genomics and Transcriptomics Approaches to Improve Female Fertility in Beef Cattle-A Review

Female fertility is the foundation of the cow-calf industry, impacting both efficiency and profitability. Reproductive failure is the primary reason why beef cows are sold in the U.S. and the cause of an estimated annual gross loss of USD 2.8 billion. In this review, we discuss the status of the genomics, transcriptomics, and systems genomics approaches currently applied to female fertility and the tools available to cow-calf producers to maximize genetic progress. We highlight the opportunities and limitations associated with using genomic and transcriptomic approaches to discover genes and regulatory mechanisms related to beef fertility. Considering the complex nature of fertility, significant advances in precision breeding will rely on holistic, multidisciplinary approaches to further advance our ability to understand, predict, and improve reproductive performance. While these technologies have advanced our knowledge, the next step is to translate research findings from bench to on-farm applications.

Whole-genome sequencing analysis of suicide deaths integrating brain-regulatory eQTLs data to identify risk loci and genes

Recent large-scale genome-wide association studies (GWAS) have started to identify potential genetic risk loci associated with risk of suicide; however, a large portion of suicide-associated genetic factors affecting gene expression remain elusive. Dysregulated gene expression, not assessed by GWAS, may play a significant role in increasing the risk of suicide death. We performed the first comprehensive genomic association analysis prioritizing brain expression quantitative trait loci (eQTLs) within regulatory regions in suicide deaths from the Utah Suicide Genetic Risk Study (USGRS). 440,324 brain-regulatory eQTLs were obtained by integrating brain eQTLs, histone modification ChIP-seq, ATAC-seq, DNase-seq, and Hi-C results from publicly available data. Subsequent genomic analyses were conducted in whole-genome sequencing (WGS) data from 986 suicide deaths of non-Finnish European (NFE) ancestry and 415 ancestrally matched controls. Additional independent USGRS suicide deaths with genotyping array data (n = 4657) and controls from the Genome Aggregation Database were explored for WGS result replication. One significant eQTL locus, rs926308 (p = 3.24e-06), was identified. The rs926308-T is associated with lower expression of RFPL3S, a gene important for neocortex development and implicated in arousal. Gene-based analyses performed using Sherlock Bayesian statistical integrative analysis also detected 20 genes with expression changes that may contribute to suicide risk. From analyzing publicly available transcriptomic data, ten of these genes have previous evidence of differential expression in suicide death or in psychiatric disorders that may be associated with suicide, including schizophrenia and autism (ZNF501, ZNF502, CNN3, IGF1R, KLHL36, NBL1, PDCD6IP, SNX19, BCAP29, and ARSA). Electronic health records (EHR) data was further merged to evaluate if there were clinically relevant subsets of suicide deaths associated with genetic variants. In summary, our study identified one risk locus and ten genes associated with suicide risk via gene expression, providing new insight into possible genetic and molecular mechanisms leading to suicide.

Joint analysis of GWAS and multi-omics QTL summary statistics reveals a large fraction of GWAS signals shared with molecular phenotypes

Molecular quantitative trait loci (xQTLs) are often harnessed to prioritize genes or functional elements underpinning variant-trait associations identified from genome-wide association studies (GWASs). Here, we introduce OPERA, a method that jointly analyzes GWAS and multi-omics xQTL summary statistics to enhance the identification of molecular phenotypes associated with complex traits through shared causal variants. Applying OPERA to summary-level GWAS data for 50 complex traits (n = 20,833-766,345) and xQTL data from seven omics layers (n = 100-31,684) reveals that 50% of the GWAS signals are shared with at least one molecular phenotype. GWAS signals shared with multiple molecular phenotypes, such as those at the MSMB locus for prostate cancer, are particularly informative for understanding the genetic regulatory mechanisms underlying complex traits. Future studies with more molecular phenotypes, measured considering spatiotemporal effects in larger samples, are required to obtain a more saturated map linking molecular intermediates to GWAS signals.

Mendelian randomization and transcriptomic analysis reveal an inverse causal relationship between Alzheimer's disease and cancer

BACKGROUND

Alzheimer's disease (AD) and cancer are common age-related diseases, and epidemiological evidence suggests an inverse relationship between them. However, investigating the potential mechanism underlying their relationship remains insufficient.

METHODS

Based on genome-wide association summary statistics for 42,034 AD patients and 609,951 cancer patients from the GWAS Catalog using the two-sample Mendelian randomization (MR) method. Moreover, we utilized two-step MR to identify metabolites mediating between AD and cancer. Furthermore, we employed colocalization analysis to identify genes whose upregulation is a risk factor for AD and demonstrated the genes' upregulation to be a favorable prognostic factor for cancer by analyzing transcriptomic data for 33 TCGA cancer types.

RESULTS

Two-sample MR analysis revealed a significant causal influence for increased AD risk on reduced cancer risk. Two-step MR analysis identified very low-density lipoprotein (VLDL) as a key mediator of the negative cause-effect relationship between AD and cancer. Colocalization analysis uncovered PVRIG upregulation to be a risk factor for AD. Transcriptomic analysis showed that PVRIG expression had significant negative correlations with stemness scores, and positive correlations with antitumor immune responses and overall survival in pan-cancer and multiple cancer types.

CONCLUSION

AD may result in lower cancer risk. VLDL is a significant intermediate variable linking AD with cancer. PVRIG abundance is a risk factor for AD but a protective factor for cancer. This study demonstrates a causal influence for AD on cancer and provides potential molecular connections between both diseases.

Prioritization of genes associated with type 2 diabetes mellitus for functional studies

Existing therapies for type 2 diabetes mellitus (T2DM) show limited efficacy or have adverse effects. Numerous genetic variants associated with T2DM have been identified, but progress in translating these findings into potential drug targets has been limited. Here, we describe the tools and platforms available to identify effector genes from T2DM-associated coding and non-coding variants and prioritize them for functional studies. We discuss QSER1 and SLC12A8 as examples of genes that have been identified as possible T2DM candidate genes using these tools and platforms. We suggest further approaches, including the use of sequencing data with increased sample size and ethnic diversity, single-cell omics data for analyses, glycaemic trait associations to predict gene function and, potentially, human induced pluripotent stem cell 'village' cultures, to strengthen current gene functionalization workflows. Effective prioritization of T2DM-associated genes for experimental validation could expedite our understanding of the genetic mechanisms responsible for T2DM to facilitate the use of precision medicine in its treatment.

Brain transcriptome-wide association study implicates novel risk genes underlying schizophrenia risk

BACKGROUND

To identify risk genes whose expression are regulated by the reported risk variants and to explore the potential regulatory mechanism in schizophrenia (SCZ).

METHODS

We systematically integrated three independent brain expression quantitative traits (eQTLs) (CommonMind, GTEx, and BrainSeq Phase 2, a total of 1039 individuals) and GWAS data (56 418 cases and 78 818 controls), with the use of transcriptome-wide association study (TWAS). Diffusion magnetic resonance imaging was utilized to quantify the integrity of white matter bundles and determine whether polygenic risk of novel genes linked to brain structure was present in patients with first-episode antipsychotic SCZ.

RESULTS

TWAS showed that eight risk genes (CORO7, DDAH2, DDHD2, ELAC2, GLT8D1, PCDHA8, THOC7, and TYW5) reached transcriptome-wide significance (TWS) level. These findings were confirmed by an independent integrative approach (i.e. Sherlock). We further conducted conditional analyses and identified the potential risk genes that driven the TWAS association signal in each locus. Gene expression analysis showed that several TWS genes (including CORO7, DDAH2, DDHD2, ELAC2, GLT8D1, THOC7 and TYW5) were dysregulated in the dorsolateral prefrontal cortex of SCZ cases compared with controls. TWS genes were mainly expressed on the surface of glutamatergic neurons, GABAergic neurons, and microglia. Finally, SCZ cases had a substantially greater TWS genes-based polygenic risk (PRS) compared to controls, and we showed that fractional anisotropy of the cingulum-hippocampus mediates the influence of TWS genes PRS on SCZ.

CONCLUSIONS

Our findings identified novel SCZ risk genes and highlighted the importance of the TWS genes in frontal-limbic dysfunctions in SCZ, indicating possible therapeutic targets.

Identification of MKNK1 and TOP3A as ovarian endometriosis risk-associated genes using integrative genomic analyses and functional experiments

The risk of endometriosis (EM), which is a common complex gynaecological disease, is related to genetic predisposition. However, it is unclear how genetic variants confer the risk of EM. Here, via Sherlock integrative analysis, we combined large-scale genome-wide association studies (GWAS) summary statistics on EM (N = 245,494) with a blood-based eQTL dataset (N = 1490) to identify EM risk-related genes. For validation, we leveraged two independent eQTL datasets (N = 769) for integration with the GWAS data. Thus, we prioritised 14 genes, including GIMAP4, TOP3A, and NMNAT3, which showed significant association with susceptibility to EM. We also utilised two independent methods, Multi-marker Analysis of GenoMic Annotation and S-PrediXcan, to further validate the EM risk-associated genes. Moreover, protein-protein interaction network analysis showed the 14 genes were functionally connected. Functional enrichment analyses further demonstrated that these genes were significantly enriched in metabolic and immune-related pathways. Differential gene expression analysis showed that in peripheral blood samples from patients with ovarian EM, TOP3A, MKNK1, SIPA1L2, and NUCB1 were significantly upregulated, while HOXB2, GIMAP5, and MGMT were significantly downregulated compared with their expression levels in samples from the controls. Immunohistochemistry further confirmed the increased expression levels of MKNK1 and TOP3A in the ectopic and eutopic endometrium compared to normal endometrium, while HOBX2 was downregulated in the endometrium of women with ovarian EM. Finally, in ex vivo functional experiments, MKNK1 knockdown inhibited ectopic endometrial stromal cells (EESCs) migration and invasion. TOP3A knockdown inhibited EESCs proliferation, migration, and invasion, while promoting their apoptosis. Convergent lines of evidence suggested that MKNK1 and TOP3A are novel EM risk-related genes.

Mendelian Randomization Study Using Dopaminergic Neuron-Specific eQTL Identifies Novel Risk Genes for Schizophrenia

Multiple integrative studies have been performed to identify the potential target genes of the non-coding schizophrenia (SCZ) risk variants. However, all the integrative studies used expression quantitative trait loci (eQTL) data from bulk tissues. Considering the cell type-specific regulatory effect of many genetic variants, it is important to conduct integrative studies using cell type-specific eQTL data. Here, we conduct a Mendelian randomization (MR) study by integrating genome-wide associations of SCZ (74,776 cases and 101,023 controls) and eQTL data (N = 215) from dopaminergic neurons, which were differentiated from human-induced pluripotent stem cell (iPSC) lines. For eQTL from young post-mitotic dopaminergic neurons (differentiation of iPSC for 30 days, D30), we identified 34 genes whose genetically regulated expression in dopaminergic neurons may have a causal role in SCZ. Among which, ARL3 showed the most significant associations with SCZ. For eQTL from more mature dopaminergic neurons (D52), we identified 37 potential SCZ causal genes, and ARL3 and GNL3 showed the most significant associations. Only 12 genes showed significant associations with SCZ in both D30 and D52 eQTL datasets, indicating the time point-specific genetic regulatory effects in young post-mitotic dopaminergic neurons and more mature dopaminergic neurons. Comparing the results from dopaminergic neurons with bulk brain tissues prioritized 2 high-confidence risk genes, including DDHD2 and GALNT10. Our study identifies multiple risk genes whose genetically regulated expression in dopaminergic neurons may have a causal role in SCZ. Further mechanistic investigation will provide pivotal insights into SCZ pathophysiology.

Exome-wide screening identifies novel rare risk variants for bone mineral density

Bone mineral density (BMD) is an independent risk factor of osteoporosis-related fractures. We performed gene-based burden tests to assess the association between rare variants and BMD, and identified several BMD candidate genes.

PURPOSE

BMD is highly heritable and a major predictor of osteoporotic fractures, but its genetic basis remains unclear. We aimed to identify rare risk variants contributing to BMD.

METHODS

Utilizing the newly released UK Biobank 200,643 exome dataset, we conducted a gene-based exome-wide association study in males and females, respectively. First, 100,639 males and 117,338 females with BMD values were included in the polygenic risk scores (PRS) analysis. Among individuals with lower 30% PRS, cases were individuals with top 10% BMD, and individuals with bottom 10% BMD were the controls. Considering the effects of vitamin D (VD), individuals with the highest 30% VD concentration were selected for VD-BMD analysis. After quality control, 741 males and 697 females were included in the BMD analysis, and 717 males and 708 females were included in the VD-BMD analysis. The variants were annotated by ANNOVAR software, then BMD and VD-BMD qualified variants were imported into the SKAT R-package to perform gene-based burden tests, respectively.

RESULTS

The gene-based burden test of the exonic variants identified genome-wide candidate associations in ANKRD18A (P = 1.60 × 10^-5, P_{Bonferroni adjust} = 2.11 × 10^-3), C22orf31 (P = 3.49 × 10^-4, P_{Bonferroni adjust} = 3.17 × 10^-2), and SPATC1L (P = 1.09 × 10^-5, P_{Bonferroni adjust} = 8.80 × 10^-3). For VD-BMD analysis, three genes were associated with BMD, such as NIPAL1 (P = 1.06 × 10^-3, P_{Bonferroni adjust} = 3.91 × 10^-2).

CONCLUSIONS

Our study suggested that rare variants contribute to BMD, providing new sights for broadening the genetic structure of BMD.

Multi-omics data integration methods and their applications in psychiatric disorders

To study mental illness and health, in the past researchers have often broken down their complexity into individual subsystems (e.g., genomics, transcriptomics, proteomics, clinical data) and explored the components independently. Technological advancements and decreasing costs of high throughput sequencing has led to an unprecedented increase in data generation. Furthermore, over the years it has become increasingly clear that these subsystems do not act in isolation but instead interact with each other to drive mental illness and health. Consequently, individual subsystems are now analysed jointly to promote a holistic understanding of the underlying biological complexity of health and disease. Complementing the increasing data availability, current research is geared towards developing novel methods that can efficiently combine the information rich multi-omics data to discover biologically meaningful biomarkers for diagnosis, treatment, and prognosis. However, clinical translation of the research is still challenging. In this review, we summarise conventional and state-of-the-art statistical and machine learning approaches for discovery of biomarker, diagnosis, as well as outcome and treatment response prediction through integrating multi-omics and clinical data. In addition, we describe the role of biological model systems and in silico multi-omics model designs in clinical translation of psychiatric research from bench to bedside. Finally, we discuss the current challenges and explore the application of multi-omics integration in future psychiatric research. The review provides a structured overview and latest updates in the field of multi-omics in psychiatry.

Incorporating genome-wide and transcriptome-wide association studies to identify genetic elements of longissimus dorsi muscle in Huaxi cattle

Locating the genetic variation of important livestock and poultry economic traits is essential for genetic improvement in breeding programs. Identifying the candidate genes for the productive ability of Huaxi cattle was one crucial element for practical breeding. Based on the genotype and phenotype data of 1,478 individuals and the RNA-seq data of 120 individuals contained in 1,478 individuals, we implemented genome-wide association studies (GWAS), transcriptome-wide association studies (TWAS), and Fisher's combined test (FCT) to identify the candidate genes for the carcass trait, the weight of longissimus dorsi muscle (LDM). The results indicated that GWAS, TWAS, and FCT identified seven candidate genes for LDM altogether: PENK was located by GWAS and FCT, PPAT was located by TWAS and FCT, and XKR4, MTMR3, FGFRL1, DHRS4, and LAP3 were only located by one of the methods. After functional analysis of these candidate genes and referring to the reported studies, we found that they were mainly functional in the progress of the development of the body and the growth of muscle cells. Combining advanced breeding techniques such as gene editing with our study will significantly accelerate the genetic improvement for the future breeding of Huaxi cattle.

GediNET for discovering gene associations across diseases using knowledge based machine learning approach

The most common approaches to discovering genes associated with specific diseases are based on machine learning and use a variety of feature selection techniques to identify significant genes that can serve as biomarkers for a given disease. More recently, the integration in this process of prior knowledge-based approaches has shown significant promise in the discovery of new biomarkers with potential translational applications. In this study, we developed a novel approach, GediNET, that integrates prior biological knowledge to gene Groups that are shown to be associated with a specific disease such as a cancer. The novelty of GediNET is that it then also allows the discovery of significant associations between that specific disease and other diseases. The initial step in this process involves the identification of gene Groups. The Groups are then subjected to a Scoring component to identify the top performing classification Groups. The top-ranked gene Groups are then used to train a Machine Learning Model. The process of Grouping, Scoring and Modelling (G-S-M) is used by GediNET to identify other diseases that are similarly associated with this signature. GediNET identifies these relationships through Disease-Disease Association (DDA) based machine learning. DDA explores novel associations between diseases and identifies relationships which could be used to further improve approaches to diagnosis, prognosis, and treatment. The GediNET KNIME workflow can be downloaded from: https://github.com/malikyousef/GediNET.git or https://kni.me/w/3kH1SQV_mMUsMTS .

Best practices for multi-ancestry, meta-analytic transcriptome-wide association studies: Lessons from the Global Biobank Meta-analysis Initiative

The Global Biobank Meta-analysis Initiative (GBMI), through its diversity, provides a valuable opportunity to study population-wide and ancestry-specific genetic associations. However, with multiple ascertainment strategies and multi-ancestry study populations across biobanks, GBMI presents unique challenges in implementing statistical genetics methods. Transcriptome-wide association studies (TWASs) boost detection power for and provide biological context to genetic associations by integrating genetic variant-to-trait associations from genome-wide association studies (GWASs) with predictive models of gene expression. TWASs present unique challenges beyond GWASs, especially in a multi-biobank, meta-analytic setting. Here, we present the GBMI TWAS pipeline, outlining practical considerations for ancestry and tissue specificity, meta-analytic strategies, and open challenges at every step of the framework. We advise conducting ancestry-stratified TWASs using ancestry-specific expression models and meta-analyzing results using inverse-variance weighting, showing the least test statistic inflation. Our work provides a foundation for adding transcriptomic context to biobank-linked GWASs, allowing for ancestry-aware discovery to accelerate genomic medicine.

Interpretation of the role of germline and somatic non-coding mutations in cancer: expression and chromatin conformation informed analysis

BACKGROUND

There has been extensive scrutiny of cancer driving mutations within the exome (especially amino acid altering mutations) as these are more likely to have a clear impact on protein functions, and thus on cell biology. However, this has come at the neglect of systematic identification of regulatory (non-coding) variants, which have recently been identified as putative somatic drivers and key germline risk factors for cancer development. Comprehensive understanding of non-coding mutations requires understanding their role in the disruption of regulatory elements, which then disrupt key biological functions such as gene expression.

MAIN BODY

We describe how advancements in sequencing technologies have led to the identification of a large number of non-coding mutations with uncharacterized biological significance. We summarize the strategies that have been developed to interpret and prioritize the biological mechanisms impacted by non-coding mutations, focusing on recent annotation of cancer non-coding variants utilizing chromatin states, eQTLs, and chromatin conformation data.

CONCLUSION

We believe that a better understanding of how to apply different regulatory data types into the study of non-coding mutations will enhance the discovery of novel mechanisms driving cancer.

Functional variant rs2270363 on 16p13.3 confers schizophrenia risk by regulating NMRAL1

Recent genome-wide association studies have reported multiple schizophrenia risk loci, yet the functional variants and their roles in schizophrenia remain to be characterized. Here we identify a functional single nucleotide polymorphism (rs2270363: G>A) at the schizophrenia risk locus 16p13.3. rs2270363 lies in the E-box element of the promoter of NMRAL1 and disrupts binding of the basic helix-loop-helix leucine zipper family proteins, including USF1, MAX and MXI1. We validated the regulatory effects of rs2270363 using reporter gene assays and electrophoretic mobility shift assay. Besides, expression quantitative trait loci analysis showed that the risk allele (A) of rs2270363 was significantly associated with elevated NMRAL1 expression in the human brain. Transcription factors knockdown and CRISPR-Cas9-mediated editing further confirmed the regulatory effects of the genomic region containing rs2270363 on NMRAL1. Intriguingly, NMRAL1 was significantly downregulated in the brain of schizophrenia patients compared with healthy subjects, and knockdown of Nmral1 expression affected proliferation and differentiation of mouse neural stem cells, as well as genes and pathways associated with brain development and synaptic transmission. Of note, Nmral1 knockdown resulted in significant decrease of dendritic spine density, revealing the potential pathophysiological mechanisms of NMRAL1 in schizophrenia. Finally, we independently confirmed the association between rs2270363 and schizophrenia in the Chinese population and found that the risk allele of rs2270363 was the same in European and Chinese populations. These lines of evidence suggest that rs2270363 may confer schizophrenia risk by regulating NMRAL1, a gene whose expression dysregulation might be involved in the pathogenesis of schizophrenia by affecting neurodevelopment and synaptic plasticity.

Predicting causal genes from psychiatric genome-wide association studies using high-level etiological knowledge

Genome-wide association studies have discovered hundreds of genomic loci associated with psychiatric traits, but the causal genes underlying these associations are often unclear, a research gap that has hindered clinical translation. Here, we present a Psychiatric Omnilocus Prioritization Score (PsyOPS) derived from just three binary features encapsulating high-level assumptions about psychiatric disease etiology - namely, that causal psychiatric disease genes are likely to be mutationally constrained, be specifically expressed in the brain, and overlap with known neurodevelopmental disease genes. To our knowledge, PsyOPS is the first method specifically tailored to prioritizing causal genes at psychiatric GWAS loci. We show that, despite its extreme simplicity, PsyOPS achieves state-of-the-art performance at this task, comparable to a prior domain-agnostic approach relying on tens of thousands of features. Genes prioritized by PsyOPS are substantially more likely than other genes at the same loci to have convergent evidence of direct regulation by the GWAS variant according to both DNA looping assays and expression or splicing quantitative trait locus (QTL) maps. We provide examples of genes hundreds of kilobases away from the lead variant, like GABBR1 for schizophrenia, that are prioritized by all three of PsyOPS, DNA looping and QTLs. Our results underscore the power of incorporating high-level knowledge of trait etiology into causal gene prediction at GWAS loci, and comprise a resource for researchers interested in experimentally characterizing psychiatric gene candidates.

Gene-Based Association Tests Using New Polygenic Risk Scores and Incorporating Gene Expression Data

Recently, gene-based association studies have shown that integrating genome-wide association studies (GWAS) with expression quantitative trait locus (eQTL) data can boost statistical power and that the genetic liability of traits can be captured by polygenic risk scores (PRSs). In this paper, we propose a new gene-based statistical method that leverages gene-expression measurements and new PRSs to identify genes that are associated with phenotypes of interest. We used a generalized linear model to associate phenotypes with gene expression and PRSs and used a score-test statistic to test the association between phenotypes and genes. Our simulation studies show that the newly developed method has correct type I error rates and can boost statistical power compared with other methods that use either gene expression or PRS in association tests. A real data analysis figure based on UK Biobank data for asthma shows that the proposed method is applicable to GWAS.

The genetic basis for adult onset glaucoma: Recent advances and future directions

Glaucoma, a diverse group of eye disorders that results in the degeneration of retinal ganglion cells, is the world's leading cause of irreversible blindness. Apart from age and ancestry, the major risk factor for glaucoma is increased intraocular pressure (IOP). In primary open-angle glaucoma (POAG), the anterior chamber angle is open but there is resistance to aqueous outflow. In primary angle-closure glaucoma (PACG), crowding of the anterior chamber angle due to anatomical alterations impede aqueous drainage through the angle. In exfoliation syndrome and exfoliation glaucoma, deposition of white flaky material throughout the anterior chamber directly interfere with aqueous outflow. Observational studies have established that there is a strong hereditable component for glaucoma onset and progression. Indeed, a succession of genome wide association studies (GWAS) that were centered upon single nucleotide polymorphisms (SNP) have yielded more than a hundred genetic markers associated with glaucoma risk. However, a shortcoming of GWAS studies is the difficulty in identifying the actual effector genes responsible for disease pathogenesis. Building on the foundation laid by GWAS studies, research groups have recently begun to perform whole exome-sequencing to evaluate the contribution of protein-changing, coding sequence genetic variants to glaucoma risk. The adoption of this technology in both large population-based studies as well as family studies are revealing the presence of novel, protein-changing genetic variants that could enrich our understanding of the pathogenesis of glaucoma. This review will cover recent advances in the genetics of primary open-angle glaucoma, primary angle-closure glaucoma and exfoliation glaucoma, which collectively make up the vast majority of all glaucoma cases in the world today. We will discuss how recent advances in research methodology have uncovered new risk genes, and how follow up biological investigations could be undertaken in order to define how the risk encoded by a genetic sequence variant comes into play in patients. We will also hypothesise how data arising from characterising these genetic variants could be utilized to predict glaucoma risk and the manner in which new therapeutic strategies might be informed.

Comprehensive and integrative analyses identify TYW5 as a schizophrenia risk gene

BACKGROUND

Identifying the causal genes at the risk loci and elucidating their roles in schizophrenia (SCZ) pathogenesis remain significant challenges. To explore risk variants associated with gene expression in the human brain and to identify genes whose expression change may contribute to the susceptibility of SCZ, here we report a comprehensive integrative study on SCZ.

METHODS

We systematically integrated the genetic associations from a large-scale SCZ GWAS (N = 56,418) and brain expression quantitative trait loci (eQTL) data (N = 175) using a Bayesian statistical framework (Sherlock) and Summary data-based Mendelian Randomization (SMR). We also measured brain structure of 86 first-episode antipsychotic-naive schizophrenia patients and 152 healthy controls with the structural MRI.

RESULTS

Both Sherlock (P = 3. 38 × 10^-6) and SMR (P = 1. 90 × 10^-8) analyses showed that TYW5 mRNA expression was significantly associated with risk of SCZ. Brain-based studies also identified a significant association between TYW5 protein abundance and SCZ. The single-nucleotide polymorphism rs203772 showed significant association with SCZ and the risk allele is associated with higher transcriptional level of TYW5 in the prefrontal cortex. We further found that TYW5 was significantly upregulated in the brain tissues of SCZ cases compared with controls. In addition, TYW5 expression was also significantly higher in neurons induced from pluripotent stem cells of schizophrenia cases compared with controls. Finally, combining analysis of genotyping and MRI data showed that rs203772 was significantly associated with gray matter volume of the right middle frontal gyrus and left precuneus.

CONCLUSIONS

We confirmed that TYW5 is a risk gene for SCZ. Our results provide useful information toward a better understanding of the genetic mechanism of TYW5 in risk of SCZ.

Integrative analysis prioritizes the relevant genes and risk factors for chronic venous disease

OBJECTIVES

Chronic venous disease (CVD) refers to a range of symptoms resulting from long-term morphological and functional abnormalities of the venous system. However, the mechanism of CVD development remains largely unknown. Here we aim to provide more information on CVD pathogenesis, prevention strategies and therapy development through the integrative analysis of large-scale genetic data.

METHODS

Genetic data were obtained from publicly accessible databases. We utilized different approaches, including FUMA, DEPICT, Sherlock, SMR/HEIDIS, DEPICT and NetWAS to identify possible causal genes for CVD. Candidate genes were prioritized to further literature-based review. The differential expression of prioritized genes was validated by microarray from the Gene Expression Omnibus (GEO), a public genomics data repository" and Real-time quantitative PCR (qPCR) of varicose veins (VVs) specimens. The causal relationships between risk factors and CVD were assessed using the Two-sample Mendelian randomization (MR) approach.

RESULTS

We identified 46 lead single-nucleotide polymorphisms (SNPs) and 26 plausible causal genes for CVD. Microarray data indicated differential expression of possible causal genes in CVD when compared to controls. The expression levels of WDR92, RSPO3, LIMA, ABCB10, DNAJC7, C1S, CXCL1 were significantly down-regulated (P<0.05). PHLDA1 and SERPINE1 were significantly upregulated (P<0.05). Dysregulated expression of WDR92, RSPO3 and CASZ1 was also found in varicose vein specimens by qPCR. Two-sample MR suggested causative effects of BMI (OR, 1.008, 95%CI, 1.005-1.010), standing height (OR, 1.009, 95%CI, 1.007-1.011), college degree (OR, 0.983, 95%CI, 0.991-0.976), insulin (OR, 0.858, 95%CI, 0.794-0.928) and metformin (OR, 0.944, 95%CI, 0.904-0.985) on CVD.

CONCLUSIONS

Our study integrates genetic and gene expression data to make an effective risk gene prediction and etiological inferences for CVD. Prioritized candidate genes provide more insights into CVD pathogenesis, and the causative effects of risk factors on CVD that deserve further investigation.

A flexible summary statistics-based colocalization method with application to the mucin cystic fibrosis lung disease modifier locus

Mucus obstruction is a central feature in the cystic fibrosis (CF) airways. A genome-wide association study (GWAS) of lung disease by the CF Gene Modifier Consortium (CFGMC) identified a significant locus containing two mucin genes, MUC20 and MUC4. Expression quantitative trait locus (eQTL) analysis using human nasal epithelia (HNE) from 94 CF-affected Canadians in the CFGMC demonstrated MUC4 eQTLs that mirrored the lung association pattern in the region, suggesting that MUC4 expression may mediate CF lung disease. Complications arose, however, with colocalization testing using existing methods: the locus is complex and the associated SNPs span a 0.2 Mb region with high linkage disequilibrium (LD) and evidence of allelic heterogeneity. We previously developed the Simple Sum (SS), a powerful colocalization test in regions with allelic heterogeneity, but SS assumed eQTLs to be present to achieve type I error control. Here we propose a two-stage SS (SS2) colocalization test that avoids a priori eQTL assumptions, accounts for multiple hypothesis testing and the composite null hypothesis, and enables meta-analysis. We compare SS2 to published approaches through simulation and demonstrate type I error control for all settings with the greatest power in the presence of high LD and allelic heterogeneity. Applying SS2 to the MUC20/MUC4 CF lung disease locus with eQTLs from CF HNE revealed significant colocalization with MUC4 (p = 1.31 × 10⁻⁵) rather than with MUC20. The SS2 is a powerful method to inform the responsible gene(s) at a locus and guide future functional studies. SS2 has been implemented in the application LocusFocus.

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

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

Fine mapping with epigenetic information and 3D structure

Since 2005, thousands of genome-wide association studies (GWAS) have been published, identifying hundreds of thousands of genetic variants that increase risk of complex traits such as autoimmune diseases. This wealth of data has the potential to improve patient care, through personalized medicine and the identification of novel drug targets. However, the potential of GWAS for clinical translation has not been fully achieved yet, due to the fact that the functional interpretation of risk variants and the identification of causal variants and genes are challenging. The past decade has seen the development of great advances that are facilitating the overcoming of these limitations, by utilizing a plethora of genomics and epigenomics tools to map and characterize regulatory elements and chromatin interactions, which can be used to fine map GWAS loci, and advance our understanding of the biological mechanisms that cause disease.

Cis-regulation of antisense noncoding RNA at the JAZF1 locus in type 2 diabetes

BACKGROUND

Several genomic loci of type 2 diabetes (T2D) nominated in genome-wide association studies (GWASs) have been suggested to regulate metabolism in muscle. However, a large portion of the genetic risk and the underlying regulation remain unexplained. This study aimed to localize the potentially functional regions or genes at juxtaposed with another zinc finger protein 1 (JAZF1) locus and interpret their possible biological mechanisms in the muscle of T2D.

METHODS AND RESULTS

With a cross-population meta-analysis of 7 GWASs, we identified a linkage disequilibrium (LD) block within intron 1 of JAZF1 that was significantly associated with T2D (FDR < 0.05). The colocalization analysis showed a significant association between genetically determined expression of JAZF1 in skeletal muscle and T2D with a strong probability of colocalization (PP4=75.09%). This region also encodes the upstream regulatory region (URR) of the antisense noncoding RNA JAZF1-AS1. Expression-QTL (e-QTL) analysis detected a regulatory SNP within this LD block, rs864745, that is associated with the expression of JAZF1-AS1 and JAZF1. With in vitro cloning, we further reported the role of JAZF1-AS1 in cis-regulating JAZF1 by directly forming RNA double strands. Downregulation of JAZF1, caused by JAZF1-AS1 depletion, inhibited the glucose uptake and lipid oxidation in skeletal muscle.

CONCLUSIONS

This study proposes a strategy to identify a novel T2D gene at the reported locus and generated a model in which polymorphisms at JAZF1 influence T2D risk through antisense-mediated gene regulation.

Secondary analyses for genome-wide association studies using expression quantitative trait loci

Genome-wide association studies (GWAS) have successfully identified thousands of single nucleotide polymorphisms (SNPs) associated with complex traits; however, the identified SNPs account for a fraction of trait heritability, and identifying the functional elements through which genetic variants exert their effects remains a challenge. Recent evidence suggests that SNPs associated with complex traits are more likely to be expression quantitative trait loci (eQTL). Thus, incorporating eQTL information can potentially improve power to detect causal variants missed by traditional GWAS approaches. Using genomic, transcriptomic, and platelet phenotype data from the Genetic Study of Atherosclerosis Risk family-based study, we investigated the potential to detect novel genomic risk loci by incorporating information from eQTL in the relevant target tissues (i.e., platelets and megakaryocytes) using established statistical principles in a novel way. Permutation analyses were performed to obtain family-wise error rates for eQTL associations, substantially lowering the genome-wide significance threshold for SNP-phenotype associations. In addition to confirming the well known association between PEAR1 and platelet aggregation, our eQTL-focused approach identified a novel locus (rs1354034) and gene (ARHGEF3) not previously identified in a GWAS of platelet aggregation phenotypes. A colocalization analysis showed strong evidence for a functional role of this eQTL.

Proteome-wide Association Study Provides Insights Into the Genetic Component of Protein Abundance in Psychiatric Disorders

BACKGROUND

Genome-wide association studies have identified multiple risk variants for psychiatric disorders. Nevertheless, how the risk variants confer risk of psychiatric disorders remains largely unknown.

METHODS

We performed proteome-wide association studies to identify genes whose cis-regulated protein abundance change in the human brain were associated with psychiatric disorders.

RESULTS

By integrating genome-wide associations of four common psychiatric disorders and two independent brain proteomes (n = 376 and n = 152, respectively) from the dorsolateral prefrontal cortex, we identified 61 genes (including 48 genes for schizophrenia, 12 genes for bipolar disorder, 5 genes for depression, and 2 genes for attention-deficit/hyperactivity disorder) whose genetically regulated protein abundance levels were associated with risk of psychiatric disorders. Comparison with transcriptome-wide association studies identified 18 overlapping genes that showed significant associations with psychiatric disorders at both proteome-wide and transcriptome-wide levels, suggesting that genetic risk variants likely confer risk of psychiatric disorders by regulating messenger RNA expression and protein abundance of these genes.

CONCLUSIONS

Our study not only provides new insights into the genetic component of protein abundance in psychiatric disorders but also highlights several high-confidence risk proteins (including CNNM2 and CTNND1) for schizophrenia and depression. These high-confidence risk proteins represent promising therapeutic targets for future drug development.

Joint-Tissue Integrative Analysis Identified Hundreds of Schizophrenia Risk Genes

Genome-wide association studies (GWAS) have identified a large number of schizophrenia risk variants, and most of them are mapped to noncoding regions. By leveraging multiple joint-tissue gene expression data and GWAS data, we herein performed a transcriptome-wide association study (TWAS) and Mendelian randomization (MR) analysis and identified 144 genes whose mRNA levels were related to genetic risk of schizophrenia. Most of these genes exhibited diametrically opposite trends of expression in prenatal and postnatal brain tissues, despite that their expression levels in dorsolateral prefrontal cortex (DLPFC) tissues did not significantly differ between schizophrenics and healthy controls. We then found significant enrichment of these genes in dopamine-related pathways that were repeatedly implicated in schizophrenia pathogenesis and in the action of antipsychotic drugs. Gene expression analysis using single cell RNA-sequencing (scRNA-seq) data of mid-gestation fetal brains further revealed enrichment of these genes in glutamatergic excitatory neurons and cycling progenitors. These lines of evidence, in consistency with previous findings, confirmed the polygenic nature of schizophrenia and highlighted involvement of early neurodevelopment aberrations in this disorder. Further investigations using advanced algorithms in both bulk brain tissues and in single cells and at different developmental stages are necessary to characterize transcriptomic features of schizophrenia pathogenesis along brain development.

Identification of an IL-22-Dependent Gene Signature as a Pharmacodynamic Biomarker

Interleukin-22 (IL-22) plays a role in epithelial barrier function and repair, and may provide benefits in conditions like inflammatory bowel disease. However, limited human data are available to assess the clinical effect of IL-22 administration. This study used a human intestinal cell line to identify an IL-22-dependent gene signature that could serve as a pharmacodynamic biomarker for IL-22 therapy. The response to IL-22Fc (UTTR1147A, an Fc-stabilized version of IL-22) was assessed in HT-29 cells by microarray, and the selected responsive genes were confirmed by qPCR. HT-29 cells demonstrated dose-dependent increases in STAT3 phosphorylation and multiple gene expression changes in response to UTTR1147A. Genes were selected that were upregulated by UTTR1147A, but to a lesser extent by IL-6, which also signals via STAT3. IL-1R1 was highly upregulated by UTTR1147A, and differential gene expression patterns were observed in response to IL-22Fc in the presence of IL-1β. An IL-22-dependent gene signature was identified that could serve as a pharmacodynamic biomarker in intestinal biopsies to support the clinical development of an IL-22 therapeutic. The differential gene expression pattern in the presence of IL-1β suggests that an inflammatory cytokine milieu in the disease setting could influence the clinical responses to IL-22.

Demystifying emerging bulk RNA-Seq applications: the application and utility of bioinformatic methodology

Significant innovations in next-generation sequencing techniques and bioinformatics tools have impacted our appreciation and understanding of RNA. Practical RNA sequencing (RNA-Seq) applications have evolved in conjunction with sequence technology and bioinformatic tools advances. In most projects, bulk RNA-Seq data is used to measure gene expression patterns, isoform expression, alternative splicing and single-nucleotide polymorphisms. However, RNA-Seq holds far more hidden biological information including details of copy number alteration, microbial contamination, transposable elements, cell type (deconvolution) and the presence of neoantigens. Recent novel and advanced bioinformatic algorithms developed the capacity to retrieve this information from bulk RNA-Seq data, thus broadening its scope. The focus of this review is to comprehend the emerging bulk RNA-Seq-based analyses, emphasizing less familiar and underused applications. In doing so, we highlight the power of bulk RNA-Seq in providing biological insights.

eQTpLot: a user-friendly R package for the visualization of colocalization between eQTL and GWAS signals

BACKGROUND

Genomic studies increasingly integrate expression quantitative trait loci (eQTL) information into their analysis pipelines, but few tools exist for the visualization of colocalization between eQTL and GWAS results. Those tools that do exist are limited in their analysis options, and do not integrate eQTL and GWAS information into a single figure panel, making the visualization of colocalization difficult.

RESULTS

To address this issue, we developed the intuitive and user-friendly R package eQTpLot. eQTpLot takes as input standard GWAS and cis-eQTL summary statistics, and optional pairwise LD information, to generate a series of plots visualizing colocalization, correlation, and enrichment between eQTL and GWAS signals for a given gene-trait pair. With eQTpLot, investigators can easily generate a series of customizable plots clearly illustrating, for a given gene-trait pair: 1) colocalization between GWAS and eQTL signals, 2) correlation between GWAS and eQTL p-values, 3) enrichment of eQTLs among trait-significant variants, 4) the LD landscape of the locus in question, and 5) the relationship between the direction of effect of eQTL signals and the direction of effect of colocalizing GWAS peaks. These clear and comprehensive plots provide a unique view of eQTL-GWAS colocalization, allowing for a more complete understanding of the interaction between gene expression and trait associations.

CONCLUSIONS

eQTpLot provides a unique, user-friendly, and intuitive means of visualizing eQTL and GWAS signal colocalization, incorporating novel features not found in other eQTL visualization software. We believe eQTpLot will prove a useful tool for investigators seeking a convenient and customizable visualization of eQTL and GWAS data colocalization.

AVAILABILITY AND IMPLEMENTATION

the eQTpLot R package and tutorial are available at https://github.com/RitchieLab/eQTpLot.

Mendelian randomization analysis identified genes potentially pleiotropically associated with periodontitis

OBJECTIVE

To prioritize genes that were pleiotropically or potentially causally associated with periodontitis.

METHODS

We applied the summary data-based Mendelian randomization (SMR) method integrating genome-wide association study (GWAS) for periodontitis and expression quantitative trait loci (eQTL) data to identify genes that were pleiotropically associated with periodontitis. We performed separate SMR analysis using CAGE eQTL data and GTEx eQTL data. SMR analysis were done for participants of European and East Asian ancestries, separately.

RESULTS

We identified multiple genes showing pleiotropic association with periodontitis in participants of European ancestry and participants of East Asian ancestry. PDCD2 (corresponding probe: ILMN_1758915) was the top hit showing pleotropic association with periodontitis in the participants of European ancestry using CAGE eQTL data, and BX093763 (corresponding probe: ILMN_1899903) and AC104135.3 (corresponding probe: ENSG00000204792.2) were the top hits in the participants of East Asian ancestry using CAGE eQTL data and GTEx eQTL data, respectively.

CONCLUSION

We identified multiple genes that may be involved in the pathogenesis of periodontitis in participants of European ancestry and participants of East Asian ancestry. Our findings provided important leads to a better understanding of the mechanisms underlying periodontitis and revealed potential therapeutic targets for the effective treatment of periodontitis.

Integrative genomics analysis reveals a 21q22.11 locus contributing risk to COVID-19

The systematic identification of host genetic risk factors is essential for the understanding and treatment of coronavirus disease 2019 (COVID-19). By performing a meta-analysis of two independent genome-wide association summary datasets (N = 680 128), a novel locus at 21q22.11 was identified to be associated with COVID-19 infection (rs9976829 in IFNAR2-IL10RB, odds ratio = 1.16, 95% confidence interval = 1.09-1.23, P = 2.57 × 10-6). The rs9976829 represents a strong splicing quantitative trait locus for both IFNAR2 and IL10RB genes, especially in lung tissue (P = 1.8 × 10-24). Integrative genomics analysis of combining genome-wide association study with expression quantitative trait locus data showed the expression variations of IFNAR2 and IL10RB have prominent effects on COVID-19 in various types of tissues, especially in lung tissue. The majority of IFNAR2-expressing cells were dendritic cells (40%) and plasmacytoid dendritic cells (38.5%), and IL10RB-expressing cells were mainly nonclassical monocytes (29.6%). IFNAR2 and IL10RB are targeted by several interferons-related drugs. Together, our results uncover 21q22.11 as a novel susceptibility locus for COVID-19, in which individuals with G alleles of rs9976829 have a higher probability of COVID-19 susceptibility than those with non-G alleles.

Advances and challenges in quantitative delineation of the genetic architecture of complex traits

Background

Genome-wide association studies (GWAS) have been widely adopted in studies of human complex traits and diseases.

Results

This review surveys areas of active research: quantifying and partitioning trait heritability, fine mapping functional variants and integrative analysis, genetic risk prediction of phenotypes, and the analysis of sequencing studies that have identified millions of rare variants. Current challenges and opportunities are highlighted.

Conclusion

GWAS have fundamentally transformed the field of human complex trait genetics. Novel statistical and computational methods have expanded the scope of GWAS and have provided valuable insights on the genetic architecture underlying complex phenotypes.

Novel Variance-Component TWAS method for studying complex human diseases with applications to Alzheimer's dementia

Transcriptome-wide association studies (TWAS) have been widely used to integrate transcriptomic and genetic data to study complex human diseases. Within a test dataset lacking transcriptomic data, traditional two-stage TWAS methods first impute gene expression by creating a weighted sum that aggregates SNPs with their corresponding cis-eQTL effects on reference transcriptome. Traditional TWAS methods then employ a linear regression model to assess the association between imputed gene expression and test phenotype, thereby assuming the effect of a cis-eQTL SNP on test phenotype is a linear function of the eQTL's estimated effect on reference transcriptome. To increase TWAS robustness to this assumption, we propose a novel Variance-Component TWAS procedure (VC-TWAS) that assumes the effects of cis-eQTL SNPs on phenotype are random (with variance proportional to corresponding reference cis-eQTL effects) rather than fixed. VC-TWAS is applicable to both continuous and dichotomous phenotypes, as well as individual-level and summary-level GWAS data. Using simulated data, we show VC-TWAS is more powerful than traditional TWAS methods based on a two-stage Burden test, especially when eQTL genetic effects on test phenotype are no longer a linear function of their eQTL genetic effects on reference transcriptome. We further applied VC-TWAS to both individual-level (N = ~3.4K) and summary-level (N = ~54K) GWAS data to study Alzheimer's dementia (AD). With the individual-level data, we detected 13 significant risk genes including 6 known GWAS risk genes such as TOMM40 that were missed by traditional TWAS methods. With the summary-level data, we detected 57 significant risk genes considering only cis-SNPs and 71 significant genes considering both cis- and trans- SNPs, which also validated our findings with the individual-level GWAS data. Our VC-TWAS method is implemented in the TIGAR tool for public use.

Integrative genomics analysis identifies five promising genes implicated in insomnia risk based on multiple omics datasets

In recent decades, many genome-wide association studies on insomnia have reported numerous genes harboring multiple risk variants. Nevertheless, the molecular functions of these risk variants conveying risk to insomnia are still ill-studied. In the present study, we integrated GWAS summary statistics (N=386,533) with two independent brain expression quantitative trait loci (eQTL) datasets (N=329) to determine whether expression-associated SNPs convey risk to insomnia. Furthermore, we applied numerous bioinformatics analyses to highlight promising genes associated with insomnia risk. By using Sherlock integrative analysis, we detected 449 significant insomnia-associated genes in the discovery stage. These identified genes were significantly overrepresented in six biological pathways including Huntington’s disease (P=5.58 × 10⁻⁵), Alzheimer’s disease (P=5.58 × 10⁻⁵), Parkinson’s disease (P=6.34 × 10⁻⁵), spliceosome (P=1.17 × 10⁻⁴), oxidative phosphorylation (P=1.09 × 10⁻⁴), and wnt signaling pathways (P=2.07 × 10⁻⁴). Further, five of these identified genes were replicated in an independent brain eQTL dataset. Through a PPI network analysis, we found that there existed highly functional interactions among these five identified genes. Three genes of LDHA (P=0.044), DALRD3 (P=5.0 × 10⁻⁵), and HEBP2 (P=0.032) showed significantly lower expression level in brain tissues of insomnic patients than that in controls. In addition, the expression levels of these five genes showed prominently dynamic changes across different time points between behavioral states of sleep and sleep deprivation in mice brain cortex. Together, the evidence of the present study strongly suggested that these five identified genes may represent candidate genes and contributed risk to the etiology of insomnia.

Agonal Factors Distort Gene-Expression Patterns in Human Postmortem Brains

Agonal factors, the conditions that occur just prior to death, can impact the molecular quality of postmortem brains, influencing gene expression results. Our study used gene expression data of 262 samples from ROSMAP with the detailed terminal state recorded for each donor, such as fever, infection, and unconsciousness. Fever and infection were the primary contributors to brain gene expression changes, brain cell-type-specific gene expression, and cell proportion changes. Furthermore, we also found that previous studies of gene expression in postmortem brains were confounded by agonal factors. Therefore, correction for agonal factors is important in the step of data preprocessing. Our analyses revealed fever and infection contributing to gene expression changes in postmortem brains and emphasized the necessity of study designs that document and account for agonal factors.

Integration of summary data from GWAS and eQTL studies identified novel risk genes for coronary artery disease

Several genetic loci have been reported to be significantly associated with coronary artery disease (CAD) by multiple genome-wide association studies (GWAS). Nevertheless, the biological and functional effects of these genetic variants on CAD remain largely equivocal. In the current study, we performed an integrative genomics analysis by integrating large-scale GWAS data (N = 459,534) and 2 independent expression quantitative trait loci (eQTL) datasets (N = 1890) to determine whether CAD-associated risk single nucleotide polymorphisms (SNPs) exert regulatory effects on gene expression. By using Sherlock Bayesian, MAGMA gene-based, multidimensional scaling (MDS), functional enrichment, and in silico permutation analyses for independent technical and biological replications, we highlighted 4 susceptible genes (CHCHD1, TUBG1, LY6G6C, and MRPS17) associated with CAD risk. Based on the protein-protein interaction (PPI) network analysis, these 4 genes were found to interact with each other. We detected a remarkably altered co-expression pattern among these 4 genes between CAD patients and controls. In addition, 3 genes of CHCHD1 (P = .0013), TUBG1 (P = .004), and LY6G6C (P = .038) showed significantly different expressions between CAD patients and controls. Together, we provide evidence to support that these identified genes such as CHCHD1 and TUBG1 are indicative factors of CAD.

Shared associations identify causal relationships between gene expression and immune cell phenotypes

Genetic mapping studies have identified thousands of associations between common variants and hundreds of human traits. Translating these associations into mechanisms is complicated by two factors: they fall into gene regulatory regions; and they are rarely mapped to one causal variant. One way around these limitations is to find groups of traits that share associations, using this genetic link to infer a biological connection. Here, we assess how many trait associations in the same locus are due to the same genetic variant, and thus shared; and if these shared associations are due to causal relationships between traits. We find that only a subset of traits share associations, with many due to causal relationships rather than pleiotropy. We therefore suggest that simply observing overlapping associations at a genetic locus is insufficient to infer causality; direct evidence of shared associations is required to support mechanistic hypotheses in genetic studies of complex traits.

Colocalization of association signals at nicotinic acetylcholine receptor genes between schizophrenia and smoking traits

BACKGROUND

Epidemiological data suggest that smoking may be a risk factor for schizophrenia (SCZ), but more evidence is needed. Two regions coding nicotinic acetylcholine receptor (nAchR) subunits, atCHRNA2 and the CHRNA5/A3/B4 cluster, were associated with SCZ in genome-wide association studies (GWAS). Additionally, a signal at CHRNA4 is near significance. CHRNA2 was also associated with cannabis use disorder (CUD). These regions were also associated with smoking behaviors. If tobacco is a risk factor, the GWAS signals at smoking behaviors and SCZ have to be due to the same causal variants, i.e. they have to colocalize, although colocalization does not necessarily imply causality. Here, we present colocalization analysis at these loci between SCZ and smoking behaviors.

METHODS

The Bayesian approach implemented in coloc was used to check for posterior probability of colocalization versus independent signals at the three loci with some evidence of association with SCZ and smoking behaviors, using GWAS summary statistics. Colocalization was also assessed for positive control traits and CUD. Three different sensibility analyses were used to confirm the results. A visualization tool, LocusCompare, was used to facilitate interpretation of the coloc results.

RESULTS

Evidence for colocalization of GWAS signals between SCZ and smoking behaviors was found for CHRNA2. Evidence for independent causal variants was found for the other two loci. CUD GWAS signal at CHRNA2 colocalizes with SCZ and smoking behaviors.

CONCLUSIONS

Overall, the results indicate that the association between some nAchR subunit genes and SCZ cannot be solely explained by their effect on smoking behaviors.

MOSTWAS: Multi-Omic Strategies for Transcriptome-Wide Association Studies

Traditional predictive models for transcriptome-wide association studies (TWAS) consider only single nucleotide polymorphisms (SNPs) local to genes of interest and perform parameter shrinkage with a regularization process. These approaches ignore the effect of distal-SNPs or other molecular effects underlying the SNP-gene association. Here, we outline multi-omics strategies for transcriptome imputation from germline genetics to allow more powerful testing of gene-trait associations by prioritizing distal-SNPs to the gene of interest. In one extension, we identify mediating biomarkers (CpG sites, microRNAs, and transcription factors) highly associated with gene expression and train predictive models for these mediators using their local SNPs. Imputed values for mediators are then incorporated into the final predictive model of gene expression, along with local SNPs. In the second extension, we assess distal-eQTLs (SNPs associated with genes not in a local window around it) for their mediation effect through mediating biomarkers local to these distal-eSNPs. Distal-eSNPs with large indirect mediation effects are then included in the transcriptomic prediction model with the local SNPs around the gene of interest. Using simulations and real data from ROS/MAP brain tissue and TCGA breast tumors, we show considerable gains of percent variance explained (1-2% additive increase) of gene expression and TWAS power to detect gene-trait associations. This integrative approach to transcriptome-wide imputation and association studies aids in identifying the complex interactions underlying genetic regulation within a tissue and important risk genes for various traits and disorders.