Explore the vast range of titles available.

Transcriptional regulation and posttranscriptional processing underlie many cellular and organismal phenotypes. We used RNA sequence data generated by Genotype-Tissue Expression (GTEx) project to investigate the patterns of transcriptome variation across individuals and tissues. Tissues exhibit characteristic transcriptional signatures that show stability in postmortem samples. These signatures are dominated by a relatively small number of genes–which is most clearly seen in blood–though few are exclusive to a particular tissue and vary more across tissues than individuals. Genes exhibiting high interindividual expression variation include disease candidates associated with sex, ethnicity, and age. Primary transcription is the major driver of cellular specificity, with splicing playing mostly a complementary role; except for the brain, which exhibits a more divergent splicing program. Variation in splicing, despite its stochasticity, may play in contrast a comparatively greater role in defining individual phenotypes.

Somatic cells can accumulate mutations over the course of an individual's lifetime. This generates cells that differ genetically at specific loci within the genome. To explore how this genetic diversity in individuals contributes to disease, Yizhak et al. developed a method to detect mutations from RNA sequencing data (see the Perspective by Tomasetti). Applying this method to Cancer Genome Atlas samples and normal samples from the Genotype-Tissue Expression (GTEx) project generated a tissue-specific study of mutation accumulation. Somatic mutations were detected in nearly all individuals and across many normal human tissues in genomic regions called cancer hotspots and in genes that play a role in cancer. Interestingly, the skin, lung, and esophagus exhibited the most mutations, suggesting that the environment generates many human mutations. Science , this issue p. eaaw0726 ; see also p. 938 “Normal” skin and other human tissues include macroscopic clonal expansions that contain genes associated with cancer risk. How somatic mutations accumulate in normal cells is poorly understood. A comprehensive analysis of RNA sequencing data from ~6700 samples across 29 normal tissues revealed multiple somatic variants, demonstrating that macroscopic clones can be found in many normal tissues. We found that sun-exposed skin, esophagus, and lung have a higher mutation burden than other tested tissues, which suggests that environmental factors can promote somatic mosaicism. Mutation burden was associated with both age and tissue-specific cell proliferation rate, highlighting that mutations accumulate over both time and number of cell divisions. Finally, normal tissues were found to harbor mutations in known cancer genes and hotspots. This study provides a broad view of macroscopic clonal expansion in human tissues, thus serving as a foundation for associating clonal expansion with environmental factors, aging, and risk of disease.

We apply integrative approaches to expression quantitative loci (eQTLs) from 44 tissues from the Genotype-Tissue Expression project and genome-wide association study data. About 60% of known trait-associated loci are in linkage disequilibrium with a cis -eQTL, over half of which were not found in previous large-scale whole blood studies. Applying polygenic analyses to metabolic, cardiovascular, anthropometric, autoimmune, and neurodegenerative traits, we find that eQTLs are significantly enriched for trait associations in relevant pathogenic tissues and explain a substantial proportion of the heritability (40–80%). For most traits, tissue-shared eQTLs underlie a greater proportion of trait associations, although tissue-specific eQTLs have a greater contribution to some traits, such as blood pressure. By integrating information from biological pathways with eQTL target genes and applying a gene-based approach, we validate previously implicated causal genes and pathways, and propose new variant and gene associations for several complex traits, which we replicate in the UK BioBank and BioVU. Integration of expression quantitative trait locus (eQTL) data from the Genotype-Tissue Expression project with genome-wide association study data shows that eQTLs are enriched for trait associations in disease-relevant tissues.

The resources generated by the GTEx consortium offer unprecedented opportunities to advance our understanding of the biology of human diseases. Here, we present an in-depth examination of the phenotypic consequences of transcriptome regulation and a blueprint for the functional interpretation of genome-wide association study-discovered loci. Across a broad set of complex traits and diseases, we demonstrate widespread dose-dependent effects of RNA expression and splicing. We develop a data-driven framework to benchmark methods that prioritize causal genes and find no single approach outperforms the combination of multiple approaches. Using colocalization and association approaches that take into account the observed allelic heterogeneity of gene expression, we propose potential target genes for 47% (2519 out of 5385) of the GWAS loci examined.

Sekar Kathiresan and colleagues perform a genome-wide association test for coronary artery disease (CAD) using data from the UK Biobank. They identify 15 new loci and perform phenome-wide association scanning, implicating insulin resistance pathways and transendothelial migration of leukocytes in CAD. UK Biobank is among the world's largest repositories for phenotypic and genotypic information in individuals of European ancestry 1 . We performed a genome-wide association study in UK Biobank testing ∼9 million DNA sequence variants for association with coronary artery disease (4,831 cases and 115,455 controls) and carried out meta-analysis with previously published results. We identified 15 new loci, bringing the total number of loci associated with coronary artery disease to 95 at the time of analysis. Phenome-wide association scanning showed that CCDC92 likely affects coronary artery disease through insulin resistance pathways, whereas experimental analysis suggests that ARHGEF26 influences the transendothelial migration of leukocytes.

DNA methylation provides a crucial epigenetic mark linking genetic variations to environmental influence. We have analyzed array-based DNA methylation profiles of 160 human retinas with co-measured RNA-seq and >8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 methylation quantitative trait loci and 12,505 expression quantitative trait loci) and 13,747 DNA methylation loci affecting gene expression, with over one-third specific to the retina. Methylation and expression quantitative trait loci show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration. Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of macular degeneration pathology by genotype–environment interaction in retina. Here, the authors perform genome-wide mapping of DNA methylation and expression quantitative trait loci, revealing associations among genotype, epigenome and transcriptome, uncovering genes and gene-environment interactions contributing to age-related macular degeneration (AMD).

Background In primary open-angle glaucoma (POAG), lowering intraocular pressure (IOP) is the only proven way of slowing vision loss. Schlemm’s canal (SC) is a hybrid vascular and lymphatic vessel that mediates aqueous humour drainage from the anterior ocular chamber. Animal studies support the importance of SC endothelial angiopoietin-TEK signalling, and more recently TIE1 signalling, in maintaining normal IOP. However, human genetic support for a causal role of TIE1 and TEK signalling in lowering IOP is currently lacking. Methods GWAS summary statistics were obtained for plasma soluble TIE1 (sTIE1) protein levels ( N  = 35,559), soluble TEK (sTEK) protein levels ( N  = 35,559), IOP ( N  = 139,555) and POAG ( N cases  = 16,677, N controls  = 199,580). Mendelian randomization (MR) was performed to estimate the association of genetically proxied TIE1 and TEK protein levels with IOP and POAG liability. Where significant MR estimates were obtained, genetic colocalization was performed to assess the probability of a shared causal variant (PP shared ) versus distinct (PP distinct ) causal variants underlying TIE1/TEK signalling and the outcome. Publicly available single-nucleus RNA-sequencing data were leveraged to investigate differential expression of TIE1 and TEK in the human ocular anterior segment. Results Increased genetically proxied TIE1 signalling and TEK signalling associated with a reduction in IOP (− 0.21 mmHg per SD increase in sTIE1, 95% CI = − 0.09 to − 0.33 mmHg, P  = 6.57 × 10 –4 , and − 0.14 mmHg per SD decrease in sTEK, 95% CI = − 0.03 to − 0.25 mmHg, P  = 0.011), but not with POAG liability. Colocalization analysis found that the probability of a shared causal variant was greater for TIE1 and IOP than for TEK and IOP (PP shared /(PP distinct  + PP shared ) = 0.98 for TIE1 and 0.30 for TEK). In the anterior segment, TIE1 and TEK were preferentially expressed in SC, lymphatic, and vascular endothelium. Conclusions This study provides novel human genetic support for a causal role of both TIE1 and TEK signalling in regulating IOP. Here, combined evidence from cis- MR and colocalization analyses provide stronger support for TIE1 than TEK as a potential IOP-lowering therapeutic target.

Psychotropic medications have been reported as a risk factor for angle closure disease. However, the interaction between background genetic risk for primary angle closure glaucoma (PACG) and susceptibility to angle closure disease among psychotropic medication users has not been investigated. Here we demonstrate the utility of a genome-wide polygenic risk score (PRS) in identifying and risk-stratifying subjects with PACG and investigate the association between PACG genetic burden and exposure to psychotropic medications on prevalent angle closure. This analysis used the UK Biobank dataset, a prospective cohort study of 502,506 UK residents. We constructed a PACG PRS for participants using genome-wide association study summary statistics from a multiethnic meta-analysis using the Lassosum method. Among the 441,054 participants, 959 (0.22%) were identified as PACG cases. Individuals with PACG had higher PRS compared to those without PACG (0.24±1.03 SD vs. 0.00±1.00 SD, p<0.001) and PACG prevalence increased with each decile of higher PRS. Among individuals using psychotropic medication, those with PACG had higher average PRS (0.31±1.00 SD vs. 0.00±1.00 SD, p<0.001) and were more likely to have a PRS in upper deciles of polygenic risk (p = 0.04). At each decile of PRS, psychotropic medication use was associated with increased risk of PACG. These effects were more pronounced and significant in higher deciles. We demonstrate the utility of a PRS for identifying individuals at higher risk of PACG. Additionally, we demonstrate an important relationship where the association between psychotropic medications use and PACG diagnosis varies across the polygenic risk spectrum.

How do aberrations in widely expressed genes lead to tissue‐selective hereditary diseases? Previous attempts to answer this question were limited to testing a few candidate mechanisms. To answer this question at a larger scale, we developed “Tissue Risk Assessment of Causality by Expression” (TRACE), a machine learning approach to predict genes that underlie tissue‐selective diseases and selectivity‐related features. TRACE utilized 4,744 biologically interpretable tissue‐specific gene features that were inferred from heterogeneous omics datasets. Application of TRACE to 1,031 disease genes uncovered known and novel selectivity‐related features, the most common of which was previously overlooked. Next, we created a catalog of tissue‐associated risks for 18,927 protein‐coding genes ( https://netbio.bgu.ac.il/trace/ ). As proof‐of‐concept, we prioritized candidate disease genes identified in 48 rare‐disease patients. TRACE ranked the verified disease gene among the patient's candidate genes significantly better than gene prioritization methods that rank by gene constraint or tissue expression. Thus, tissue selectivity combined with machine learning enhances genetic and clinical understanding of hereditary diseases. Synopsis An interpretable machine‐learning framework predicts disease genes for tissue‐selective hereditary diseases. The framework highlights known and novel tissue‐selectivity features and enhances genetic diagnosis of rare‐disease patients. An interpretable machine‐learning (ML) framework uses thousands of gene features to predict tissue‐associated disease genes. ML models highlight known and novel tissue‐selectivity mechanisms. An online catalogue of tissue‐associated risks for 18,927 protein‐coding genes in eight tissues is presented. The framework and catalogue enhance genetic diagnosis of rare‐disease patients. Graphical Abstract An interpretable machine‐learning framework predicts disease genes for tissue‐selective hereditary diseases. The framework highlights known and novel tissue‐selectivity features and enhances genetic diagnosis of rare‐disease patients.

Congenital heart disease (CHD) occurs in ∼1% of newborns. CHD arises from many distinct etiologies, ranging from genetic or genomic variation to exposure to teratogens, which elicit diverse cell and molecular responses during cardiac development. To systematically explore the relationships between CHD risk factors and responses, we compiled and integrated comprehensive datasets from studies of CHD in humans and model organisms. We examined two alternative models of potential functional relationships between genes in these datasets: direct convergence, in which CHD risk factors significantly and directly impact the same genes and molecules and functional convergence, in which risk factors significantly impact different molecules that participate in a discrete heart development network. We observed no evidence for direct convergence. In contrast, we show that CHD risk factors functionally converge in protein networks driving the development of specific anatomical structures (e.g., outflow tract, ventricular septum, and atrial septum) that are malformed by CHD. This integrative analysis of CHD risk factors and responses suggests a complex pattern of functional interactions between genomic variation and environmental exposures that modulate critical biological systems during heart development.