Search

Transcriptomics data have been integrated with genome-wide association studies (GWASs) to help understand disease/trait molecular mechanisms. The utility of metabolomics, integrated with transcriptomics and disease GWASs, to understand molecular mechanisms for metabolite levels or diseases has not been thoroughly evaluated. We performed probabilistic transcriptome-wide association and locus-level colocalization analyses to integrate transcriptomics results for 49 tissues in 706 individuals from the GTEx project, metabolomics results for 1,391 plasma metabolites in 6,136 Finnish men from the METSIM study, and GWAS results for 2,861 disease traits in 260,405 Finnish individuals from the FinnGen study. We found that genetic variants that regulate metabolite levels were more likely to influence gene expression and disease risk compared to the ones that do not. Integrating transcriptomics with metabolomics results prioritized 397 genes for 521 metabolites, including 496 previously identified gene-metabolite pairs with strong functional connections and suggested 33.3% of such gene-metabolite pairs shared the same causal variants with genetic associations of gene expression. Integrating transcriptomics and metabolomics individually with FinnGen GWAS results identified 1,597 genes for 790 disease traits. Integrating transcriptomics and metabolomics jointly with FinnGen GWAS results helped pinpoint metabolic pathways from genes to diseases. We identified putative causal effects of UGT1A1/UGT1A4 expression on gallbladder disorders through regulating plasma (E,E)-bilirubin levels, of SLC22A5 expression on nasal polyps and plasma carnitine levels through distinct pathways, and of LIPC expression on age-related macular degeneration through glycerophospholipid metabolic pathways. Our study highlights the power of integrating multiple sets of molecular traits and GWAS results to deepen understanding of disease pathophysiology.

Summary In this study, we integrated large-scale GWAS summary data and used the predicted transcriptome-wide association study method to discover novel genes associated with osteoporosis. We identified 204 candidate genes, which provide novel clues for understanding the genetic mechanism of osteoporosis and indicate potential therapeutic targets. Introduction Osteoporosis is a highly polygenetic disease characterized by low bone mass and deterioration of the bone microarchitecture. Our objective was to discover novel candidate genes associated with osteoporosis. Methods To identify potential causal genes of the associated loci, we investigated trait-gene expression associations using the transcriptome-wide association study (TWAS) method. This method directly imputes gene expression effects from genome-wide association study (GWAS) data using a statistical prediction model trained on GTEx reference transcriptome data. We then performed a colocalization analysis to evaluate the posterior probability of biological patterns: associations characterized by a single causal variant or multiple distinct causal variants. Finally, a functional enrichment analysis of gene sets was performed using the VarElect and CluePedia tools, which assess the causal relationships between genes and a disease and search for potential gene’s functional pathways. The osteoporosis-associated genes were further confirmed based on the differentially expressed genes profiled from mRNA expression data of bone tissue. Results Our analysis identified 204 candidate genes, including 154 genes that have been previously associated with osteoporosis, 50 genes that have not been previously discovered. A biological function analysis found that 20 of the candidate genes were directly associated with osteoporosis. Further analysis of multiple gene expression profiles showed that 15 genes were differentially expressed in patients with osteoporosis. Among these, SLC11A2, MAP2K5, NFATC4, and HSP90B1 were enriched in four pathways, namely, mineral absorption pathway, MAPK signaling pathway, Wnt signaling pathway, and PI3K-Akt signaling pathway, which indicates a causal relationship with the occurrence of osteoporosis. Conclusions We demonstrated that transcriptome fine-mapping identifies more osteoporosis-related genes and provides key insight into the development of novel targeted therapeutics for the treatment of osteoporosis. Supplementary Information The online version contains supplementary material available at 10.1007/s00198-021-06024-z.

In transcriptome‐wide association studies (TWAS), gene expression values are predicted using genotype data and tested for association with a phenotype. The power of this approach to detect associations relies, at least in part, on the accuracy of the prediction. Here we compare the prediction accuracy of six different methods—LASSO, Ridge regression, Elastic net, Best Linear Unbiased Predictor, Bayesian Sparse Linear Mixed Model, and Random Forests—by performing cross‐validation using data from the Geuvadis Project. We also examine prediction accuracy (a) at different sample sizes, (b) when ancestry of the prediction model training and testing populations is different, and (c) when the tissue used to train the model is different from the tissue to be predicted. We find that, for most genes, the expression cannot be accurately predicted, but in general sparse statistical models tend to outperform polygenic models at prediction. Average prediction accuracy is reduced when the model training set size is reduced or when predicting across ancestries and is marginally reduced when predicting across tissues. We conclude that using sparse statistical models and the development of large reference panels across multiple ethnicities and tissues will lead to better prediction of gene expression, and thus may improve TWAS power.

Genome-wide association studies (GWAS) have identified over 60 genetic loci associated with immunoglobulin G (IgG) N-glycosylation; however, the causal genes and their abundance in relevant tissues are uncertain. Leveraging data from GWAS summary statistics for 8,090 Europeans, and large-scale expression quantitative trait loci (eQTL) data from the genotype-tissue expression of 53 types of tissues (GTEx v7), we derived a linkage disequilibrium score for the specific expression of genes (LDSC-SEG) and conducted a transcriptome-wide association study (TWAS). We identified 55 gene associations whose predicted levels of expression were significantly associated with IgG N-glycosylation in 14 tissues. Three working scenarios, i.e., tissue-specific, pleiotropic, and coassociated, were observed for candidate genetic predisposition affecting IgG N-glycosylation traits. Furthermore, pathway enrichment showed several IgG N-glycosylation-related pathways, such as asparagine N-linked glycosylation, N-glycan biosynthesis and transport to the Golgi and subsequent modification. Through phenome-wide association studies (PheWAS), most genetic variants underlying TWAS hits were found to be correlated with health measures (height, waist-hip ratio, systolic blood pressure) and diseases, such as systemic lupus erythematosus, inflammatory bowel disease, and Parkinson’s disease, which are related to IgG N-glycosylation. Our study provides an atlas of genetic regulatory loci and their target genes within functionally relevant tissues, for further studies on the mechanisms of IgG N-glycosylation and its related diseases.

Funder: University of Hawai'i at Mānoa; doi: http://dx.doi.org/10.13039/100008783, Background: Genome-wide association studies (GWAS) have identified over 56 susceptibility loci associated with Alzheimer’s disease (AD), but the genes responsible for these associations remain largely unknown. Methods: We performed a large transcriptome-wide association study (TWAS) leveraging modified UTMOST (Unified Test for MOlecular SignaTures) prediction models of ten brain tissues that are potentially related to AD to discover novel AD genetic loci and putative target genes in 71,880 (proxy) cases and 383,378 (proxy) controls of European ancestry. Results: We identified 53 genes with predicted expression associations with AD risk at Bonferroni correction threshold (P value < 3.38 × 10−6). Based on fine-mapping analyses, 21 genes at nine loci showed strong support for being causal. Conclusions: Our study provides new insights into the etiology and underlying genetic architecture of AD.

Genome-wide association study (GWAS) of Juvenile idiopathic arthritis (JIA) suffers from low power due to limited sample size and the interpretation challenge due to most signals located in non-coding regions. Gene-level analysis could alleviate these issues. Using GWAS summary statistics, we performed two typical gene-level analysis of JIA, transcriptome-wide association studies (TWAS) using FUnctional Summary-based ImputatiON (FUSION) and gene-based analysis using eQTL Multi-marker Analysis of GenoMic Annotation (eMAGMA), followed by comprehensive enrichment analysis. Among 33 overlapped significant genes from these two methods, 11 were previously reported, including TYK2 (PFUSION = 5.12 × 10−6, PeMAGMA = 1.94 × 10−7 for whole blood), IL-6R (PFUSION = 8.63 × 10−7, PeMAGMA = 2.74 × 10−6 for cells EBV-transformed lymphocytes), and Fas (PFUSION = 5.21 × 10−5, PeMAGMA = 1.08 × 10−6 for muscle skeletal). Some newly plausible JIA-associated genes are also reported, including IL-27 (PFUSION = 2.10 × 10−7, PeMAGMA = 3.93 × 10−8 for Liver), LAT (PFUSION = 1.53 × 10−4, PeMAGMA = 4.62 × 10−7 for Artery Aorta), and MAGI3 (PFUSION = 1.30 × 10−5, PeMAGMA = 1.73 × 10−7 for Muscle Skeletal). Enrichment analysis further highlighted 4 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and 10 Gene Ontology (GO) terms. Our findings can benefit the understanding of genetic determinants and potential therapeutic targets for JIA.

Contains fulltext : 237689.pdf (Publisher’s version ) (Open Access) Familial, sequencing, and genome-wide association studies (GWASs) and genetic correlation analyses have progressively unraveled the shared or pleiotropic germline genetics of breast and ovarian cancer. In this study, we aimed to leverage this shared germline genetics to improve the power of transcriptome-wide association studies (TWASs) to identify candidate breast cancer and ovarian cancer susceptibility genes. We built gene expression prediction models using the PrediXcan method in 681 breast and 295 ovarian tumors from The Cancer Genome Atlas and 211 breast and 99 ovarian normal tissue samples from the Genotype-Tissue Expression project and integrated these with GWAS meta-analysis data from the Breast Cancer Association Consortium (122,977 cases/105,974 controls) and the Ovarian Cancer Association Consortium (22,406 cases/40,941 controls). The integration was achieved through application of a pleiotropy-guided conditional/conjunction false discovery rate (FDR) approach in the setting of a TWASs. This identified 14 candidate breast cancer susceptibility genes spanning 11 genomic regions and 8 candidate ovarian cancer susceptibility genes spanning 5 genomic regions at conjunction FDR < 0.05 that were >1 Mb away from known breast and/or ovarian cancer susceptibility loci. We also identified 38 candidate breast cancer susceptibility genes and 17 candidate ovarian cancer susceptibility genes at conjunction FDR < 0.05 at known breast and/or ovarian susceptibility loci. The 22 genes identified by our cross-cancer analysis represent promising candidates that further elucidate the role of the transcriptome in mediating germline breast and ovarian cancer risk.

Summary Reliably generating rice varieties with low glycaemic index (GI) is an important nutritional intervention given the high rates of Type II diabetes incidences in Asia where rice is staple diet. We integrated a genome‐wide association study (GWAS) with a transcriptome‐wide association study (TWAS) to determine the genetic basis of the GI in rice. GWAS utilized 305 re‐sequenced diverse indica panel comprising ~2.4 million single nucleotide polymorphisms (SNPs) enriched in genic regions. A novel association signal was detected at a synonymous SNP in exon 2 of LOC_Os05g03600 for intermediate‐to‐high GI phenotypic variation. Another major hotspot region was predicted for contributing intermediate‐to‐high GI variation, involves 26 genes on chromosome 6 (GI6.1). These set of genes included GBSSI, two hydrolase genes, genes involved in signalling and chromatin modification. The TWAS and methylome sequencing data revealed cis‐acting functionally relevant genetic variants with differential methylation patterns in the hot spot GI6.1 region, narrowing the target to 13 genes. Conversely, the promoter region of GBSSI and its alternative splicing allele (G allele of Wx a) explained the intermediate‐to‐high GI variation. A SNP (C˃T) at exon‐10 was also highlighted in the preceding analyses to influence final viscosity (FV), which is independent of amylose content/GI . The low GI line with GC haplotype confirmed soft texture, while other two low GI lines with GT haplotype were characterized as hard and cohesive. The low GI lines were further confirmed through clinical in vivo studies. Gene regulatory network analysis highlighted the role of the non‐starch polysaccharide pathway in lowering GI.

BackgroundThe etiology of frontotemporal dementia (FTD) is poorly understood. To identify genes with predicted expression levels associated with FTD, we integrated summary statistics with external reference gene expression data using a transcriptome-wide association study approach.MethodsFUSION software was used to leverage FTD summary statistics (all FTD: n= 2154 cases, n= 4308 controls; behavioral variant FTD: n= 1337 cases, n= 2754 controls; semantic dementia: n= 308 cases, n= 616 controls; progressive nonfluent aphasia: n= 269 cases, n= 538 controls; FTD with motor neuron disease: n= 200 cases, n= 400 controls) from the International FTD-Genomics Consortium with 53 expression quantitative loci tissue type panels (n= 12,205; 5 consortia). Significance was assessed using a 5% false discovery rate threshold.ResultsWe identified 73 significant gene-tissue associations for FTD, representing 44 unique genes in 34 tissue types. Most significant findings were derived from dorsolateral prefrontal cortex splicing data (n= 19 genes, 26%). The 17q21.31 inversion locus contained 23 significant associations, representing 6 unique genes. Other top hits included SEC22B (a gene involved in vesicle trafficking), TRGV5, and ZNF302. A single gene finding (RAB38) was observed for behavioral variant FTD. For other clinical subtypes, no significant associations were observed.ConclusionsWe identified novel candidate genes (e.g., SEC22B) and previously reported risk regions (e.g., 17q21.31) for FTD. Most significant associations were observed in dorsolateral prefrontal cortex splicing data despite the modest sample size of this reference panel. This suggests that our findings are specific to FTD and are likely to be biologically relevant highlights of genes at different FTD risk loci that are contributing to the disease pathology.

BACKGROUND: The etiology of frontotemporal dementia (FTD) is poorly understood. To identify genes with predicted expression levels associated with FTD, we integrated summary statistics with external reference gene expression data using a transcriptome-wide association study approach. METHODS: FUSION software was used to leverage FTD summary statistics (all FTD: n = 2154 cases, n = 4308 controls; behavioral variant FTD: n = 1337 cases, n = 2754 controls; semantic dementia: n = 308 cases, n = 616 controls; progressive nonfluent aphasia: n = 269 cases, n = 538 controls; FTD with motor neuron disease: n = 200 cases, n = 400 controls) from the International FTD-Genomics Consortium with 53 expression quantitative loci tissue type panels (n = 12,205; 5 consortia). Significance was assessed using a 5% false discovery rate threshold. RESULTS: We identified 73 significant gene-tissue associations for FTD, representing 44 unique genes in 34 tissue types. Most significant findings were derived from dorsolateral prefrontal cortex splicing data (n = 19 genes, 26%). The 17q21.31 inversion locus contained 23 significant associations, representing 6 unique genes. Other top hits included SEC22B (a gene involved in vesicle trafficking), TRGV5, and ZNF302. A single gene finding (RAB38) was observed for behavioral variant FTD. For other clinical subtypes, no significant associations were observed. CONCLUSIONS: We identified novel candidate genes (e.g., SEC22B) and previously reported risk regions (e.g., 17q21.31) for FTD. Most significant associations were observed in dorsolateral prefrontal cortex splicing data despite the modest sample size of this reference panel. This suggests that our findings are specific to FTD and are likely to be biologically relevant highlights of genes at different FTD risk loci that are contributing to the disease pathology.

Summary Mismatch negativity (MMN) is a differential electrophysiological response measuring cortical adaptability to unpredictable stimuli. MMN is consistently attenuated in patients with psychosis. However, the genetics of MMN are uncharted, limiting the validation of MMN as a psychosis endophenotype. Here, we perform a transcriptome-wide association study of 728 individuals, which reveals 2 genes (FAM89A and ENGASE) whose expression in cortical tissues is associated with MMN. Enrichment analyses of neurodevelopmental expression signatures show that genes associated with MMN tend to be overexpressed in the frontal cortex during prenatal development but are significantly downregulated in adulthood. Endophenotype ranking value calculations comparing MMN and three other candidate psychosis endophenotypes (lateral ventricular volume and two auditory-verbal learning measures) find MMN to be considerably superior. These results yield promising insights into sensory processing in the cortex and endorse the notion of MMN as a psychosis endophenotype., Graphical abstract, Highlights • Expression of FAM89A and ENGASE is associated with reduced mismatch negativity (MMN) • Genes regulating neurotransmitter levels increase in MMN transcriptome-wide association • Genes influencing MMN are underexpressed in adult brain cortex • MMN ranks above verbal memory and ventricular volume as psychosis endophenotype, Bhat et al. identify two genes, FAM89A and ENGASE, whose expression in cortical tissue is negatively associated with mismatch negativity (MMN), an electrophysiological measure of cortical responses to unexpected sensory stimuli. They find enrichment of neurotransmission-regulating genes in these associations and endorse MMN as an endophenotype for psychosis.

Introduction:Humans and viruses have co-evolved for millennia resulting in a complex host genetic architecture. Understanding the genetic mechanisms of immune response to viral infection provides insight into disease etiology and informs public health interventions. Methods:We conducted a comprehensive study linking germline genetic variation and gene expression with antibody response to 28 antigens for 16 viruses using serological data from 7924 participants in the UK Biobank cohort. Using test results from 2010 UK Biobank subjects, we also investigated genetic determinants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Results:Signals in human leukocyte antigen (HLA) class II region dominated the landscape of viral antibody response, with 40 independent loci and 14 independent classical alleles, 7 of which exhibited pleiotropic effects across viral families. Genome-wide association analyses discovered 7 novel genetic loci associated with viral antibody response (P

Aim To identify novel candidate genes for pancreatic cancer. Methods We performed a transcriptome-wide association study (TWAS) analysis of pancreatic cancer (PC). GWAS summary data were driven from the published studies of PC, totally involving 558 542 SNPs in 1896 individuals with pancreatic cancer and 1939 healthy controls. FUSION software was applied to the PC GWAS summary data for tissue-related TWAS analysis, including whole blood, peripheral blood, adipose, and pancreas. The functional relevance of identified genes with PC was further validated by Oncomine, STRING, and CluePedia tool. Results Transcriptome-wide association study analysis identified 19 genes significantly associated with PC, such as LRP5L (P value = 5.21 × 10-5 ), SOX4 (P value = 3.2 × 10-4 ), and EGLN3 (P value = 6.2 × 10-3 ). KEGG pathway enrichment analysis detected several PC-associated pathways, such as One carbon pool by folate (P value = 1.60 × 10-16 ), Cell cycle (P value = 1.27 × 10-7 ), TGF-beta signaling pathway (P value = 4.64 × 10-6 ). Further comparing the 19 genes with previously identified overexpressed genes in PC patients found one overlapped gene SOX4. Conclusion We identified some novel candidate genes and pathways associated with PC. Our results provide novel clues for the genetic mechanism studies of pancreatic cancer.

Dupuytren’s disease is the common inherited tissue-specific fibrotic disorder. It’s characterized by progressive and irreversible fibroblastic proliferation affecting the palmar fascia of the hand, with an onset typically in the sixth decade of life. Although genome-wide association studies (GWAS) have identified 24 genomic regions associated with Dupuytren’s risk, the biological mechanisms driving signal at these regions remain elusive. We identify potential biological mechanisms for Dupuytren’s disease by integrating the most recent, largest GWAS (ncases = 3,871, ncontrols = 4,686) with eQTLs (47 tissue panels from five consortia, total n = 3,975) to perform a transcriptome-wide association study (TWAS). We identify 43 tissue-specific gene associations with Dupuytren’s risk, one of which resides at least 0.5 Mb away from the 24 risk regions previously identified. We also estimate the genome-wide genetic correlation between Dupuytren’s disease and 45 complex traits and find significant genetic correlations between Dupuytren’s disease and body mass index (, P = 1.6 × 10−6), type II diabetes (, P = 1.7 × 10−4), triglycerides (, P = 3.5 × 10−4), and high-density lipoprotein (, P = 4.1 × 10−4), which suggests a shared genetic etiology. We further refine the genome-wide genetic correlation signal to identify 8 regions significantly negatively correlated with BMI and 3 regions significantly correlated (1 positively and 2 negatively correlated) with HDL; none of these regions contained the novel gene association identified by TWAS. Our results are consistent with previous epidemiological findings which show that lower BMI increases risk for Dupuytren’s disease. These 12 novel risk regions provide new insight into the biological mechanisms of Dupuytren’s disease and serve as a starting point for functional validation.

Although genome-wide association studies (GWASs) have identified thousands of risk loci for many complex traits and diseases, the causal variants and genes at these loci remain largely unknown. Here, we introduce a method for estimating the local genetic correlation between gene expression and a complex trait and utilize it to estimate the genetic correlation due to predicted expression between pairs of traits. We integrated gene expression measurements from 45 expression panels with summary GWAS data to perform 30 multi-tissue transcriptome-wide association studies (TWASs). We identified 1,196 genes whose expression is associated with these traits; of these, 168reside more than 0.5 Mb away from any previously reported GWAS significant variant. We then used our approach to find 43 pairs of traits with significant genetic correlation at the level of predicted expression; of these, eight were not found through genetic correlation at the SNP level. Finally, we used bi-directional regression to find evidence that BMI causally influences triglyceride levels and that triglyceride levels causally influence low-density lipoprotein. Together, our results provide insight into the role of gene expression in the susceptibility of complex traits and diseases.