Cardiomyocyte-derived circulating extracellular vesicles allow a non-invasive liquid biopsy of myocardium in health and disease.

The ability to track disease without tissue biopsy in patients is a major goal in biology and medicine. Here, we identify and characterize cardiomyocyte-derived extracellular vesicles in circulation (EVs; "cardiovesicles") through comprehensive studies of induced pluripotent stem cell-derived cardiomyocytes, genetic mouse models, and state-of-the-art mass spectrometry and low-input transcriptomics. These studies identified two markers ( POPDC2 , CHRNE ) enriched on cardiovesicles for biotinylated antibody-based immunocapture. Captured cardiovesicles were enriched in canonical cardiomyocyte transcripts/pathways with distinct profiles based on human disease type (heart failure, myocardial infarction). In paired myocardial tissue-plasma from patients, highly expressed genes in cardiovesicles were largely cardiac-enriched (vs. "bulk" EVs, which were more organ non-specific) with high expression in myocardial tissue by single nuclear RNA-seq, largely in cardiomyocytes. These results demonstrate the first "liquid" biopsy discovery platform to interrogate cardiomyocyte states non-invasively in model systems and in human disease, allowing non-invasive characterization of cardiomyocyte biology for discovery and therapeutic applications.

Multiancestry transferability of a polygenic risk score for diverticulitis.

OBJECTIVE: Polygenic risk scores (PRS) for diverticular disease must be evaluated in diverse cohorts. We sought to explore shared genetic predisposition across the phenome and to assess risk stratification in individuals genetically similar to European, African and Admixed-American reference samples. METHODS: A 44-variant PRS was applied to the All of Us Research Program. Phenome-wide association studies (PheWAS) identified conditions linked with heightened genetic susceptibility to diverticular disease. To evaluate the PRS in risk stratification, logistic regression models for symptomatic and for severe diverticulitis were compared with base models with covariates of age, sex, body mass index, smoking and principal components. Performance was assessed using area under the receiver operating characteristic curves (AUROC) and Nagelkerke's R2. RESULTS: The cohort comprised 181 719 individuals for PheWAS and 50 037 for risk modelling. PheWAS identified associations with diverticular disease, connective tissue disease and hernias. Across ancestry groups, one SD PRS increase was consistently associated with greater odds of severe (range of ORs (95% CI) 1.60 (1.27 to 2.02) to 1.86 (1.42 to 2.42)) and of symptomatic diverticulitis ((95% CI) 1.27 (1.10 to 1.46) to 1.66 (1.55 to 1.79)) relative to controls. European models achieved the highest AUROC and Nagelkerke's R2 (AUROC (95% CI) 0.78 (0.75 to 0.81); R2 0.25). The PRS provided a maximum R2 increase of 0.034 and modest AUROC improvement. CONCLUSION: Associations between a diverticular disease PRS and severe presentations persisted in diverse cohorts when controlling for known risk factors. Relative improvements in model performance were observed, but absolute change magnitudes were modest.

Myeloid-Specific JAK2 Contributes to Inflammation and Salt Sensitivity of Blood Pressure.

BACKGROUND: Salt sensitivity of blood pressure (SSBP), characterized by acute changes in blood pressure with changes in dietary sodium intake, is an independent risk factor for cardiovascular disease and mortality in people with and without hypertension. We previously found that elevated sodium concentration activates antigen-presenting cells (APCs), resulting in high blood pressure, but the mechanisms are unknown. Here, we hypothesized that APC-specific JAK2 (Janus kinase 2) through STAT3 (signal transducer and activator of transcription 3) and SMAD3 (small mothers against decapentaplegic homolog 3) contributes to SSBP. METHODS: We performed bulk or single-cell transcriptomic analyses following in vitro monocytes exposed to high salt and in vivo high sodium treatment in humans using a rigorous salt-loading/depletion protocol to phenotype SSBP. We also used a myeloid cell-specific CD11c+ JAK2 knockout mouse model and measured blood pressure with radiotelemetry after N-omega-nitro-L-arginine-methyl ester and a high salt diet treatment. We used flow cytometry for immunophenotyping and measuring cytokine levels. Fluorescence in situ hybridization and immunohistochemistry were performed to spatially visualize the kidney's immune cells and cytokine levels. Echocardiography was performed to assess cardiac function. RESULTS: We found that high salt treatment upregulates gene expression of the JAK/STAT/SMAD pathway while downregulating inhibitors of this pathway, such as suppression of cytokine signaling and cytokine-inducible SH2, in human monocytes. Expression of the JAK2 pathway genes mirrored changes in blood pressure after salt loading and depletion in salt-sensitive but not salt-resistant humans. Ablation of JAK2, specifically in CD11c+ APCs, attenuated salt-induced hypertension in mice with SSBP. Mechanistically, we found that SMAD3 acted downstream of JAK2 and STAT3, leading to increased production of highly reactive isolevuglandins and proinflammatory cytokine IL (interleukin)-6 in renal APCs, which activate T cells and increase production of IL-17A, IL-6, and TNF-α (tumor necrosis factor-alpha). CONCLUSIONS: Our findings reveal the APC JAK2 signaling pathway as a potential target for the diagnosis and treatment of SSBP in humans.

Metabolic gene function discovery platform GeneMAP identifies SLC25A48 as necessary for mitochondrial choline import.

Organisms maintain metabolic homeostasis through the combined functions of small-molecule transporters and enzymes. While many metabolic components have been well established, a substantial number remains without identified physiological substrates. To bridge this gap, we have leveraged large-scale plasma metabolome genome-wide association studies (GWAS) to develop a multiomic Gene-Metabolite Association Prediction (GeneMAP) discovery platform. GeneMAP can generate accurate predictions and even pinpoint genes that are distant from the variants implicated by GWAS. In particular, our analysis identified solute carrier family 25 member 48 (SLC25A48) as a genetic determinant of plasma choline levels. Mechanistically, SLC25A48 loss strongly impairs mitochondrial choline import and synthesis of its downstream metabolite betaine. Integrative rare variant and polygenic score analyses in UK Biobank provide strong evidence that the SLC25A48 causal effects on human disease may in part be mediated by the effects of choline. Altogether, our study provides a discovery platform for metabolic gene function and proposes SLC25A48 as a mitochondrial choline transporter.

Transcriptome-wide association study and Mendelian randomization in pancreatic cancer identifies susceptibility genes and causal relationships with type 2 diabetes and venous thromboembolism.

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 <0.10), including six genes not yet reported for PDAC risk (PPIP5K2, TFR2, HNF4G, LRRC10B, PRC1 and FBXL20) and ten previously reported genes (INHBA, SMC2, ABO, PDX1, MTMR6, ACOT2, PGAP3, STARD3, GSDMB, ADAM33). MR provided support for a causal effect of PDAC on T2D using genetic instruments in the HNF4G and PDX1 loci, and unidirectional causality of VTE on PDAC involving the ABO locus (OR 2.12, P < 1e-7). No evidence of a causal effect of PDAC on VTE was found. INTERPRETATION: These analyses identified candidate susceptibility genes and disease relationships for PDAC that warrant further investigation. HNF4G and PDX1 may induce PDAC-associated diabetes, whereas ABO may induce the causative effect of VTE on PDAC. FUNDING: National Institutes of Health (USA).

Genetic correlation and causal associations between psychiatric disorders and lung cancer risk.

BACKGROUND: Patients with certain psychiatric disorders have increased lung cancer incidence. However, establishing a causal relationship through traditional epidemiological methods poses challenges. METHODS: Available summary statistics of genome-wide association studies of cigarette smoking, lung cancer, and eight psychiatric disorders, including attention deficit/hyperactivity disorder (ADHD), autism, depression, major depressive disorder, bipolar disorder, insomnia, neuroticism, and schizophrenia (range N: 46,350-1,331,010) were leveraged to estimate genetic correlations using Linkage Disequilibrium Score Regression and assess causal effect of each psychiatric disorder on lung cancer using two-sample Mendelian randomization (MR) models, comprising inverse-variance weighted (IVW), weighted median, MR-Egger, pleiotropy residual sum and outlier testing (MR-PRESSO), and a constrained maximum likelihood approach (cML-MR). RESULTS: Significant positive correlations were observed between each psychiatric disorder and both smoking and lung cancer (all FDR < 0.05), except for the correlation between autism and lung cancer. Both univariable and the cML-MA MR analyses demonstrated that liability to schizophrenia, depression, ADHD, or insomnia was associated with an increased risk of overall lung cancer. Genetic liability to insomnia was linked specifically to squamous cell carcinoma (SCC), while genetic liability to ADHD was associated with an elevated risk of both SCC and small cell lung cancer (all P < 0.05). The later was further supported by multivariable MR analyses, which accounted for smoking. LIMITATIONS: Participants were constrained to European ancestry populations. Causal estimates from binary psychiatric disorders may be biased. CONCLUSION: Our findings suggest appropriate management of several psychiatric disorders, particularly ADHD, may potentially reduce the risk of developing lung cancer.

Novel ancestry-specific primary open-angle glaucoma loci and shared biology with vascular mechanisms and cell proliferation.

Primary open-angle glaucoma (POAG), a leading cause of irreversible blindness globally, shows disparity in prevalence and manifestations across ancestries. We perform meta-analysis across 15 biobanks (of the Global Biobank Meta-analysis Initiative) (n = 1,487,441: cases = 26,848) and merge with previous multi-ancestry studies, with the combined dataset representing the largest and most diverse POAG study to date (n = 1,478,037: cases = 46,325) and identify 17 novel significant loci, 5 of which were ancestry specific. Gene-enrichment and transcriptome-wide association analyses implicate vascular and cancer genes, a fifth of which are primary ciliary related. We perform an extensive statistical analysis of SIX6 and CDKN2B-AS1 loci in human GTEx data and across large electronic health records showing interaction between SIX6 gene and causal variants in the chr9p21.3 locus, with expression effect on CDKN2A/B. Our results suggest that some POAG risk variants may be ancestry specific, sex specific, or both, and support the contribution of genes involved in programmed cell death in POAG pathogenesis.

Alternative polyadenylation quantitative trait methylation mapping in human cancers provides clues into the molecular mechanisms of APA.

Genetic variants are involved in the orchestration of alternative polyadenylation (APA) events, while the role of DNA methylation in regulating APA remains unclear. We generated a comprehensive atlas of APA quantitative trait methylation sites (apaQTMs) across 21 different types of cancer (1,612 to 60,219 acting in cis and 4,448 to 142,349 in trans). Potential causal apaQTMs in non-cancer samples were also identified. Mechanistically, we observed a strong enrichment of cis-apaQTMs near polyadenylation sites (PASs) and both cis- and trans-apaQTMs in proximity to transcription factor (TF) binding regions. Through the integration of ChIP-signals and RNA-seq data from cell lines, we have identified several regulators of APA events, acting either directly or indirectly, implicating novel functions of some important genes, such as TCF7L2, which is known for its involvement in type 2 diabetes and cancers. Furthermore, we have identified a vast number of QTMs that share the same putative causal CpG sites with five different cancer types, underscoring the roles of QTMs, including apaQTMs, in the process of tumorigenesis. DNA methylation is extensively involved in the regulation of APA events in human cancers. In an attempt to elucidate the potential underlying molecular mechanisms of APA by DNA methylation, our study paves the way for subsequent experimental validations into the intricate biological functions of DNA methylation in APA regulation and the pathogenesis of human cancers. To present a comprehensive catalog of apaQTM patterns, we introduce the Pancan-apaQTM database, available at https://pancan-apaqtm-zju.shinyapps.io/pancanaQTM/.

Integrative genetic analysis identifies FLVCR1 as a plasma-membrane choline transporter in mammals.

Genome-wide association studies (GWASs) of serum metabolites have the potential to uncover genes that influence human metabolism. Here, we combined an integrative genetic analysis that associates serum metabolites to membrane transporters with a coessentiality map of metabolic genes. This analysis revealed a connection between feline leukemia virus subgroup C cellular receptor 1 (FLVCR1) and phosphocholine, a downstream metabolite of choline metabolism. Loss of FLVCR1 in human cells strongly impairs choline metabolism due to the inhibition of choline import. Consistently, CRISPR-based genetic screens identified phospholipid synthesis and salvage machinery as synthetic lethal with FLVCR1 loss. Cells and mice lacking FLVCR1 exhibit structural defects in mitochondria and upregulate integrated stress response (ISR) through heme-regulated inhibitor (HRI) kinase. Finally, Flvcr1 knockout mice are embryonic lethal, which is partially rescued by choline supplementation. Altogether, our findings propose FLVCR1 as a major choline transporter in mammals and provide a platform to discover substrates for unknown metabolite transporters.

A Local Genetic Correlation Analysis Provides Biological Insights Into the Shared Genetic Architecture of Psychiatric and Substance Use Phenotypes.

BACKGROUND: Global genetic correlation analysis has provided valuable insight into the shared genetic basis between psychiatric and substance use disorders. However, little is known about which regions disproportionately contribute to the global correlation. METHODS: We used Local Analysis of [co]Variant Annotation to calculate bivariate local genetic correlations across 2495 approximately equal-sized, semi-independent genomic regions for 20 psychiatric and substance use phenotypes. We performed a transcriptome-wide association study using expression weights from the prefrontal cortex to identify risk genes for each phenotype, followed by probabilistic fine-mapping to prioritize credible causal genes within each bivariate locus. RESULTS: We detected 80 significant (p < 2.08 × 10-6) bivariate local genetic correlations across 61 loci. The expression effect directions for risk genes within each bivariate locus were largely consistent with the local correlation coefficients, suggesting that genetically regulated gene expression may be used in the functional interpretation of local genetic correlations. Probabilistic fine-mapping identified several genes that may drive pleiotropic mechanisms for genetically correlated phenotypes. For example, we confirmed a local genetic correlation between schizophrenia and smoking behavior at 15q25 and prioritized PSMA4 as the most credible gene candidate underlying both phenotypes. CONCLUSIONS: Our study reveals previously unreported local bivariate genetic correlations between psychiatric and substance use phenotypes, which we fine-mapped to identify shared credible causal genes underlying genetically correlated phenotypes.

A transcriptome-wide association study identifies novel candidate susceptibility genes for prostate cancer risk.

A large proportion of heritability for prostate cancer risk remains unknown. Transcriptome-wide association study combined with validation comparing overall levels will help to identify candidate genes potentially playing a role in prostate cancer development. Using data from the Genotype-Tissue Expression Project, we built genetic models to predict normal prostate tissue gene expression using the statistical framework PrediXcan, a modified version of the unified test for molecular signatures and Joint-Tissue Imputation. We applied these prediction models to the genetic data of 79 194 prostate cancer cases and 61 112 controls to investigate the associations of genetically determined gene expression with prostate cancer risk. Focusing on associated genes, we compared their expression in prostate tumor vs normal prostate tissue, compared methylation of CpG sites located at these loci in prostate tumor vs normal tissue, and assessed the correlations between the differentiated genes' expression and the methylation of corresponding CpG sites, by analyzing The Cancer Genome Atlas (TCGA) data. We identified 573 genes showing an association with prostate cancer risk at a false discovery rate (FDR) ≤ 0.05, including 451 novel genes and 122 previously reported genes. Of the 573 genes, 152 showed differential expression in prostate tumor vs normal tissue samples. At loci of 57 genes, 151 CpG sites showed differential methylation in prostate tumor vs normal tissue samples. Of these, 20 CpG sites were correlated with expression of 11 corresponding genes. In this TWAS, we identified novel candidate susceptibility genes for prostate cancer risk, providing new insights into prostate cancer genetics and biology.

CD36 maintains the gastric mucosa and associates with gastric disease.

The gastric epithelium is often exposed to injurious elements and failure of appropriate healing predisposes to ulcers, hemorrhage, and ultimately cancer. We examined the gastric function of CD36, a protein linked to disease and homeostasis. We used the tamoxifen model of gastric injury in mice null for Cd36 (Cd36-/-), with Cd36 deletion in parietal cells (PC-Cd36-/-) or in endothelial cells (EC-Cd36-/-). CD36 expresses on corpus ECs, on PC basolateral membranes, and in gastrin and ghrelin cells. Stomachs of Cd36-/- mice have altered gland organization and secretion, more fibronectin, and inflammation. Tissue respiration and mitochondrial efficiency are reduced. Phospholipids increased and triglycerides decreased. Mucosal repair after injury is impaired in Cd36-/- and EC-Cd36-/-, not in PC-Cd36-/- mice, and is due to defect of progenitor differentiation to PCs, not of progenitor proliferation or mature PC dysfunction. Relevance to humans is explored in the Vanderbilt BioVu using PrediXcan that links genetically-determined gene expression to clinical phenotypes, which associates low CD36 mRNA with gastritis, gastric ulcer, and gastro-intestinal hemorrhage. A CD36 variant predicted to disrupt an enhancer site associates (p < 10-17) to death from gastro-intestinal hemorrhage in the UK Biobank. The findings support role of CD36 in gastric tissue repair, and its deletion associated with chronic diseases that can predispose to malignancy.

SLC25A39 is necessary for mitochondrial glutathione import in mammalian cells.

Glutathione (GSH) is a small-molecule thiol that is abundant in all eukaryotes and has key roles in oxidative metabolism1. Mitochondria, as the major site of oxidative reactions, must maintain sufficient levels of GSH to perform protective and biosynthetic functions2. GSH is synthesized exclusively in the cytosol, yet the molecular machinery involved in mitochondrial GSH import remains unknown. Here, using organellar proteomics and metabolomics approaches, we identify SLC25A39, a mitochondrial membrane carrier of unknown function, as a regulator of GSH transport into mitochondria. Loss of SLC25A39 reduces mitochondrial GSH import and abundance without affecting cellular GSH levels. Cells lacking both SLC25A39 and its paralogue SLC25A40 exhibit defects in the activity and stability of proteins containing iron-sulfur clusters. We find that mitochondrial GSH import is necessary for cell proliferation in vitro and red blood cell development in mice. Heterologous expression of an engineered bifunctional bacterial GSH biosynthetic enzyme (GshF) in mitochondria enables mitochondrial GSH production and ameliorates the metabolic and proliferative defects caused by its depletion. Finally, GSH availability negatively regulates SLC25A39 protein abundance, coupling redox homeostasis to mitochondrial GSH import in mammalian cells. Our work identifies SLC25A39 as an essential and regulated component of the mitochondrial GSH-import machinery.

Multilayer modelling of the human transcriptome and biological mechanisms of complex diseases and traits.

Here, we performed a comprehensive intra-tissue and inter-tissue multilayer network analysis of the human transcriptome. We generated an atlas of communities in gene co-expression networks in 49 tissues (GTEx v8), evaluated their tissue specificity, and investigated their methodological implications. UMAP embeddings of gene expression from the communities (representing nearly 18% of all genes) robustly identified biologically-meaningful clusters. Notably, new gene expression data can be embedded into our algorithmically derived models to accelerate discoveries in high-dimensional molecular datasets and downstream diagnostic or prognostic applications. We demonstrate the generalisability of our approach through systematic testing in external genomic and transcriptomic datasets. Methodologically, prioritisation of the communities in a transcriptome-wide association study of the biomarker C-reactive protein (CRP) in 361,194 individuals in the UK Biobank identified genetically-determined expression changes associated with CRP and led to considerably improved performance. Furthermore, a deep learning framework applied to the communities in nearly 11,000 tumors profiled by The Cancer Genome Atlas across 33 different cancer types learned biologically-meaningful latent spaces, representing metastasis (p < 2.2 × 10-16) and stemness (p < 2.2 × 10-16). Our study provides a rich genomic resource to catalyse research into inter-tissue regulatory mechanisms, and their downstream consequences on human disease.

Deep Learning Enables Fast and Accurate Imputation of Gene Expression.

A question of fundamental biological significance is to what extent the expression of a subset of genes can be used to recover the full transcriptome, with important implications for biological discovery and clinical application. To address this challenge, we propose two novel deep learning methods, PMI and GAIN-GTEx, for gene expression imputation. In order to increase the applicability of our approach, we leverage data from GTEx v8, a reference resource that has generated a comprehensive collection of transcriptomes from a diverse set of human tissues. We show that our approaches compare favorably to several standard and state-of-the-art imputation methods in terms of predictive performance and runtime in two case studies and two imputation scenarios. In comparison conducted on the protein-coding genes, PMI attains the highest performance in inductive imputation whereas GAIN-GTEx outperforms the other methods in in-place imputation. Furthermore, our results indicate strong generalization on RNA-Seq data from 3 cancer types across varying levels of missingness. Our work can facilitate a cost-effective integration of large-scale RNA biorepositories into genomic studies of disease, with high applicability across diverse tissue types.

A Transcriptome-Wide Association Study Identifies Candidate Susceptibility Genes for Pancreatic Cancer Risk.

Pancreatic cancer is among the most well-characterized cancer types, yet a large proportion of the heritability of pancreatic cancer risk remains unclear. Here, we performed a large transcriptome-wide association study to systematically investigate associations between genetically predicted gene expression in normal pancreas tissue and pancreatic cancer risk. Using data from 305 subjects of mostly European descent in the Genotype-Tissue Expression Project, we built comprehensive genetic models to predict normal pancreas tissue gene expression, modifying the UTMOST (unified test for molecular signatures). These prediction models were applied to the genetic data of 8,275 pancreatic cancer cases and 6,723 controls of European ancestry. Thirteen genes showed an association of genetically predicted expression with pancreatic cancer risk at an FDR ≤ 0.05, including seven previously reported genes (INHBA, SMC2, ABO, PDX1, RCCD1, CFDP1, and PGAP3) and six novel genes not yet reported for pancreatic cancer risk [6q27: SFT2D1 OR (95% confidence interval (CI), 1.54 (1.25-1.89); 13q12.13: MTMR6 OR (95% CI), 0.78 (0.70-0.88); 14q24.3: ACOT2 OR (95% CI), 1.35 (1.17-1.56); 17q12: STARD3 OR (95% CI), 6.49 (2.96-14.27); 17q21.1: GSDMB OR (95% CI), 1.94 (1.45-2.58); and 20p13: ADAM33 OR (95% CI): 1.41 (1.20-1.66)]. The associations for 10 of these genes (SFT2D1, MTMR6, ACOT2, STARD3, GSDMB, ADAM33, SMC2, RCCD1, CFDP1, and PGAP3) remained statistically significant even after adjusting for risk SNPs identified in previous genome-wide association study. Collectively, this analysis identified novel candidate susceptibility genes for pancreatic cancer that warrant further investigation. SIGNIFICANCE: A transcriptome-wide association analysis identified seven previously reported and six novel candidate susceptibility genes for pancreatic cancer risk.

Metabolic coessentiality mapping identifies C12orf49 as a regulator of SREBP processing and cholesterol metabolism.

Coessentiality mapping has been useful to systematically cluster genes into biological pathways and identify gene functions1-3. Here, using the debiased sparse partial correlation (DSPC) method3, we construct a functional coessentiality map for cellular metabolic processes across human cancer cell lines. This analysis reveals 35 modules associated with known metabolic pathways and further assigns metabolic functions to unknown genes. In particular, we identify C12orf49 as an essential regulator of cholesterol and fatty acid metabolism in mammalian cells. Mechanistically, C12orf49 localizes to the Golgi, binds membrane-bound transcription factor peptidase, site 1 (MBTPS1, site 1 protease) and is necessary for the cleavage of its substrates, including sterol regulatory element binding protein (SREBP) transcription factors. This function depends on the evolutionarily conserved uncharacterized domain (DUF2054) and promotes cell proliferation under cholesterol depletion. Notably, c12orf49 depletion in zebrafish blocks dietary lipid clearance in vivo, mimicking the phenotype of mbtps1 mutants. Finally, in an electronic health record (EHR)-linked DNA biobank, C12orf49 is associated with hyperlipidaemia through phenome analysis. Altogether, our findings reveal a conserved role for C12orf49 in cholesterol and lipid homeostasis and provide a platform to identify unknown components of other metabolic pathways.

Electronic health record phenotypes associated with genetically regulated expression of CFTR and application to cystic fibrosis.

PURPOSE: The increasing use of electronic health records (EHRs) and biobanks offers unique opportunities to study Mendelian diseases. We described a novel approach to summarize clinical manifestations from patient EHRs into phenotypic evidence for cystic fibrosis (CF) with potential to alert unrecognized patients of the disease. METHODS: We estimated genetically predicted expression (GReX) of cystic fibrosis transmembrane conductance regulator (CFTR) and tested for association with clinical diagnoses in the Vanderbilt University biobank (N = 9142 persons of European descent with 71 cases of CF). The top associated EHR phenotypes were assessed in combination as a phenotype risk score (PheRS) for discriminating CF case status in an additional 2.8 million patients from Vanderbilt University Medical Center (VUMC) and 125,305 adult patients including 25,314 CF cases from MarketScan, an independent external cohort. RESULTS: GReX of CFTR was associated with EHR phenotypes consistent with CF. PheRS constructed using the EHR phenotypes and weights discovered by the genetic associations improved discriminative power for CF over the initially proposed PheRS in both VUMC and MarketScan. CONCLUSION: Our study demonstrates the power of EHRs for clinical description of CF and the benefits of using a genetics-informed weighing scheme in construction of a phenotype risk score. This research may find broad applications for phenomic studies of Mendelian disease genes.

Genomic Variants of Cytarabine Sensitivity Associated with Treatment-Related Mortality in Pediatric AML: A Report from the Children's Oncology Group.

PURPOSE: Cytarabine is an effective treatment for AML with associated toxicities including treatment related mortality (TRM). The purpose is to determine the clinical relevance of SNPs identified through the use of HapMap lymphoblastoid cell-based models, in predicting cytarabine response and toxicity in AML. EXPERIMENTAL DESIGN: We tested clinical significance of SNPs associated with cytarabine sensitivity in children with AML treated on Children's Oncology Group regimens (CCG 2941/2961). Endpoints included overall survival (OS), event-free survival (EFS), and TRM. Patients who received bone marrow transplant were excluded. We tested 124 SNPs associated with cytarabine sensitivity in HapMap cell lines in 348 children to determine whether any associated with treatment outcomes. In addition, we tested five SNPs previously associated with TRM in children with AML in our independent dataset of 385 children. RESULTS: Homozygous variant genotypes of rs2025501 and rs6661575 had increased in vitro cellular sensitivity to cytarabine and were associated with increased TRM. TRM was particularly increased in children with variant genotype randomized to high-dose cytarabine (rs2025501: P = 0.0024 and rs6661575 P = 0.0188). In analysis of previously reported SNPs, only the variant genotype rs17202778 C/C was significantly associated with TRM (P < 0.0001). CONCLUSIONS: We report clinical importance of two SNPs not previously associated with cytarabine toxicity. Moreover, we confirm that SNP rs17202778 significantly impacts TRM in pediatric AML. Cytarabine sensitivity genotypes may predict TRM and could be used to stratify to standard versus high-dose cytarabine regimens, warranting further study in prospective AML trials.

Post-GWAS analysis of six substance use traits improves the identification and functional interpretation of genetic risk loci.

BACKGROUND: Little is known about the functional mechanisms through which genetic loci associated with substance use traits ascertain their effect. This study aims to identify and functionally annotate loci associated with substance use traits based on their role in genetic regulation of gene expression. METHODS: We evaluated expression Quantitative Trait Loci (eQTLs) from 13 brain regions and whole blood of the Genotype-Tissue Expression (GTEx) database, and from whole blood of the Depression Genes and Networks (DGN) database. The role of single eQTLs was examined for six substance use traits: alcohol consumption (N = 537,349), cigarettes per day (CPD; N = 263,954), former vs. current smoker (N = 312,821), age of smoking initiation (N = 262,990), ever smoker (N = 632,802), and cocaine dependence (N = 4,769). Subsequently, we conducted a gene level analysis of gene expression on these substance use traits using S-PrediXcan. RESULTS: Using an FDR-adjusted p-value <0.05 we found 2,976 novel candidate genetic loci for substance use traits, and identified genes and tissues through which these loci potentially exert their effects. Using S-PrediXcan, we identified significantly associated genes for all substance traits. DISCUSSION: Annotating genes based on transcriptomic regulation improves the identification and functional characterization of candidate loci and genes for substance use traits.