Search

Summary: Background: Tissue-specific analysis of the transcriptome is critical to elucidating the molecular basis of complex traits, but central tissues are often not accessible. We propose a methodology, Multi-mOdal-based framework to bridge the Transcriptome between PEripheral and Central tissues (MOTPEC). Methods: Multi-modal regulatory elements in peripheral blood are incorporated as features for gene expression prediction in 48 central tissues. To demonstrate the utility, we apply it to the identification of BMI-associated genes and compare the tissue-specific results with those derived directly from surrogate blood. Findings: MOTPEC models demonstrate superior performance compared with both baseline models in blood and existing models across the 48 central tissues. We identify a set of BMI-associated genes using the central tissue MOTPEC-predicted transcriptome data. The MOTPEC-based differential gene expression (DGE) analysis of BMI in the central tissues (including brain caudate basal ganglia and visceral omentum adipose tissue) identifies 378 genes overlapping the results from a TWAS of BMI, while only 162 overlapping genes are identified using gene expression in blood. Cellular perturbation analysis further supports the utility of MOTPEC for identifying trait-associated gene sets and narrowing the effect size divergence between peripheral blood and central tissues. Interpretation: The MOTPEC framework improves the gene expression prediction accuracy for central tissues and enhances the identification of tissue-specific trait-associated genes. Funding: This research is supported by the National Natural Science Foundation of China 82204118 (D.Z.), the seed funding of the Key Laboratory of Intelligent Preventive Medicine of Zhejiang Province (2020E10004), the National Institutes of Health (NIH) Genomic Innovator Award R35HG010718 (E.R.G.), NIH/NHGRI R01HG011138 (E.R.G.), NIH/NIA R56AG068026 (E.R.G.), NIH Office of the Director U24OD035523 (E.R.G.), and NIH/NIGMS R01GM140287 (E.R.G.).

Summary: Background: Two important questions regarding the genetics of pancreatic adenocarcinoma (PDAC) are 1. Which germline genetic variants influence the incidence of this cancer; and 2. Whether PDAC causally predisposes to associated non-malignant phenotypes, such as type 2 diabetes (T2D) and venous thromboembolism (VTE). Methods: In this study of 8803 patients with PDAC and 67,523 controls, we first performed a large-scale transcriptome-wide association study to investigate the association between genetically determined gene expression in normal pancreas tissue and PDAC risk. Secondly, we used Mendelian Randomization (MR) to analyse the causal relationships among PDAC, T2D (74,124 cases and 824,006 controls) and VTE (30,234 cases and 172,122 controls). Findings: Sixteen genes showed an association with PDAC risk (FDR

Deep Mutational Scanning (DMS) has enabled multiplexed measurement of mutational effects on protein properties, including kinematics and self-organization, with unprecedented resolution. However, potential bottlenecks of DMS characterization include experimental design, data quality, and depth of mutational coverage. Here, we apply deep learning to comprehensively model the mutational effect of the Alzheimer’s Disease associated peptide Aβ42 on aggregation-related biochemical traits from DMS measurements. Among tested neural network architectures, Convolutional Neural Networks and Recurrent Neural Networks are found to be the most cost-effective models with high performance even under insufficiently-sampled DMS studies. While sequence features are essential for satisfactory prediction from neural networks, geometric-structural features further enhance the prediction performance. Notably, we demonstrate how mechanistic insights into phenotype may be extracted from the neural networks themselves suitably designed. This methodological benefit is particularly relevant for biochemical systems displaying a strong coupling between structure and phenotype such as the conformation of Aβ42 aggregate and nucleation, as shown here using a Graph Convolutional Neural Network (GCN) developed from the protein atomic structure input. In addition to accurate imputation of missing values (which here ranged up to 55% of all phenotype values at key residues), the mutationally-defined nucleation phenotype generated from a GCN shows improved resolution for identifying known disease-causing mutations relative to the original DMS phenotype. Our study suggests that neural network derived sequence-phenotype mapping can be exploited not only to provide direct support for protein engineering or genome editing but also to facilitate therapeutic design with the gained perspectives from biological modeling.

Summary: Polygenic risk scores (PRSs) have been widely explored in precision medicine. However, few studies have thoroughly investigated their best practices in global populations across different diseases. We here utilized data from Global Biobank Meta-analysis Initiative (GBMI) to explore methodological considerations and PRS performance in 9 different biobanks for 14 disease endpoints. Specifically, we constructed PRSs using pruning and thresholding (P + T) and PRS-continuous shrinkage (CS). For both methods, using a European-based linkage disequilibrium (LD) reference panel resulted in comparable or higher prediction accuracy compared with several other non-European-based panels. PRS-CS overall outperformed the classic P + T method, especially for endpoints with higher SNP-based heritability. Notably, prediction accuracy is heterogeneous across endpoints, biobanks, and ancestries, especially for asthma, which has known variation in disease prevalence across populations. Overall, we provide lessons for PRS construction, evaluation, and interpretation using GBMI resources and highlight the importance of best practices for PRS in the biobank-scale genomics era.

Summary: Deep mutational scanning (DMS) experiments have been performed on SARS-CoV-2’s spike receptor-binding domain (RBD) and human angiotensin-converting enzyme 2 (ACE2) zinc-binding peptidase domain—both central players in viral infection and evolution and antibody evasion—quantifying how mutations impact biochemical phenotypes. We modeled biochemical phenotypes from massively parallel assays, using neural networks trained on protein sequence mutations in the virus and human host. Neural networks were significantly predictive of binding affinity, protein expression, and antibody escape, learning complex interactions and higher-order features that are difficult to capture with conventional methods from structural biology. Integrating the physicochemical properties of amino acids, such as hydrophobicity and long-range non-bonded energy per atom, significantly improved prediction (empirical p

Summary: We conducted PrediXcan analysis of hydrocephalus risk in ten neurological tissues and whole blood. Decreased expression of MAEL in the brain was significantly associated (Bonferroni-adjusted p < 0.05) with hydrocephalus. PrediXcan analysis of brain imaging and genomics data in the independent UK Biobank (N = 8,428) revealed that MAEL expression in the frontal cortex is associated with white matter and total brain volumes. Among the top differentially expressed genes in brain, we observed a significant enrichment for gene-level associations with these structural phenotypes, suggesting an effect on disease risk through regulation of brain structure and integrity. We found additional support for these genes through analysis of the choroid plexus transcriptome of a murine model of hydrocephalus. Finally, differential protein expression analysis in patient cerebrospinal fluid recapitulated disease-associated expression changes in neurological tissues, but not in whole blood. Our findings provide convergent evidence highlighting the importance of tissue-specific pathways and mechanisms in the pathophysiology of hydrocephalus.