Search

Supplementary Table S5 from BET Bromodomain Inhibition Blocks an AR-Repressed, E2F1-Activated Treatment-Emergent Neuroendocrine Prostate Cancer Lineage Plasticity Program

Supplementary Figure S3 from BET Bromodomain Inhibition Blocks an AR-Repressed, E2F1-Activated Treatment-Emergent Neuroendocrine Prostate Cancer Lineage Plasticity Program

Supplementary Data from BET Bromodomain Inhibition Blocks an AR-Repressed, E2F1-Activated Treatment-Emergent Neuroendocrine Prostate Cancer Lineage Plasticity Program

Supplementary Table S2 S4 S6 S7 from BET Bromodomain Inhibition Blocks an AR-Repressed, E2F1-Activated Treatment-Emergent Neuroendocrine Prostate Cancer Lineage Plasticity Program

Purpose:Lineage plasticity in prostate cancer—most commonly exemplified by loss of androgen receptor (AR) signaling and a switch from a luminal to alternate differentiation program—is now recognized as a treatment resistance mechanism. Lineage plasticity is a spectrum, but neuroendocrine prostate cancer (NEPC) is the most virulent example. Currently, there are limited treatments for NEPC. Moreover, the incidence of treatment-emergent NEPC (t-NEPC) is increasing in the era of novel AR inhibitors. In contradistinction to de novo NEPC, t-NEPC tumors often express the AR, but AR's functional role in t-NEPC is unknown. Furthermore, targetable factors that promote t-NEPC lineage plasticity are also unclear.Experimental Design:Using an integrative systems biology approach, we investigated enzalutamide-resistant t-NEPC cell lines and their parental, enzalutamide-sensitive adenocarcinoma cell lines. The AR is still expressed in these t-NEPC cells, enabling us to determine the role of the AR and other key factors in regulating t-NEPC lineage plasticity.Results:AR inhibition accentuates lineage plasticity in t-NEPC cells—an effect not observed in parental, enzalutamide-sensitive adenocarcinoma cells. Induction of an AR-repressed, lineage plasticity program is dependent on activation of the transcription factor E2F1 in concert with the BET bromodomain chromatin reader BRD4. BET inhibition (BETi) blocks this E2F1/BRD4-regulated program and decreases growth of t-NEPC tumor models and a subset of t-NEPC patient tumors with high activity of this program in a BETi clinical trial.Conclusions:E2F1 and BRD4 are critical for activating an AR-repressed, t-NEPC lineage plasticity program. BETi is a promising approach to block this program.

Supplementary Table S8-S12 from BET Bromodomain Inhibition Blocks an AR-Repressed, E2F1-Activated Treatment-Emergent Neuroendocrine Prostate Cancer Lineage Plasticity Program

Human genome-wide association studies found single-nucleotide polymorphisms (SNPs) near LYPLAL1 (Lysophospholipase-like protein 1) that have sex-specific effects on fat distribution and metabolic traits. To determine whether altering LYPLAL1 affects obesity and metabolic disease, we created and characterized a mouse knockout (KO) of Lyplal1. We fed the experimental group of mice a high-fat, high-sucrose (HFHS) diet for 23 weeks, and the controls were fed regular chow diet. Here, we show that CRISPR-Cas9 whole-body Lyplal1 KO mice fed an HFHS diet showed sex-specific differences in weight gain and fat accumulation as compared to chow diet. Female, not male, KO mice weighed less than WT mice, had reduced body fat percentage, had white fat mass, and had adipocyte diameter not accounted for by changes in the metabolic rate. Female, but not male, KO mice had increased serum triglycerides, decreased aspartate, and decreased alanine aminotransferase. Lyplal1 KO mice of both sexes have reduced liver triglycerides and steatosis. These diet-specific effects resemble the effects of SNPs near LYPLAL1 in humans, suggesting that LYPLAL1 has an evolutionary conserved sex-specific effect on adiposity. This murine model can be used to study this novel gene-by-sex-by-diet interaction to elucidate the metabolic effects of LYPLAL1 on human obesity.

The SH-SY5Y neuroblastoma cells are a widely used in vitro model approximating neurons for testing the target engagement of therapeutics designed for neurodegenerative diseases and pain disorders. However, their potential as a model for receptor-mediated delivery and uptake of novel modalities, such as antibody-drug conjugates, remains understudied. Investigation of the SH-SY5Y cell surfaceome will aid in greater in vitro to in vivo correlation of delivery and uptake, thereby accelerating drug discovery. So far, the majority of studies have focused on total cell proteomics from undifferentiated and differentiated SH-SY5Y cells. While some studies have investigated the expression of specific proteins in neuroblastoma tissue, a global approach for comparison of neuroblastoma cell surfaceome to the brain and dorsal root ganglion (DRG) neurons remains uninvestigated. Furthermore, an isoform-specific evaluation of cell surface proteins expressed on neuroblastoma cells remains unexplored. In this study, we define a bioinformatic workflow for the identification of high-confidence surface proteins expressed on brain and DRG neurons using tissue proteomic and transcriptomic data. We then delineate the SH-SY5Y cell surfaceome by surface proteomics and show that it significantly overlaps with the human brain and DRG neuronal surface proteome. We find that, for 32% of common surface proteins, SH-SY5Y-specific major isoforms are alternatively spliced, maintaining their protein-coding ability, and are predicted to localize to the cell surface. Validation of these isoforms using surface proteomics confirms a SH-SY5Y-specific alternative NRCAM (neuron-glia related cell adhesion molecule) isoform, which is absent in typical brain neurons, but present in neuroblastomas, making it a receptor of interest for neuroblastoma-specific therapeutics.

Psoriasis and type 2 diabetes (T2D) are complex conditions with significant impact on health. Psoriasis patients have higher risk of type 2 diabetes (~1.5 Odds Ratio) and vice versa, controlling for body mass index (BMI), yet there has been limited study comparing their genetic architecture. We hypothesized there are shared genetic components between psoriasis and T2D. Trans-disease meta-analysis (TDMA) was applied to 8,016,731 well-imputed genetic markers from large-scale meta-analyses of psoriasis (11,024 cases and 16,336 controls) and T2D adjusted for BMI (74,124 cases and 824,006 controls). We confirmed our findings in a hospital-based study (42,112 patients) and tested for causal relationships with multi-variable Mendelian randomization. Mendelian randomization identified a causal relationship between psoriasis and T2D (p=1.6x10(−4), OR=1.01), and highlighted the impact of BMI. TDMA further revealed 4 genome-wide significant loci (p

Serum liver enzyme concentrations are the most frequently-used laboratory markers of liver disease, a major cause of mortality. We conduct a meta-analysis of genome-wide association studies of liver enzymes from UK BioBank and BioBank Japan. We identified 160 previously-unreported independent alanine aminotransferase, 190 aspartate aminotransferase, and 199 alkaline phosphatase genome-wide significant associations, with some affecting multiple different enzymes. Associated variants implicate genes that demonstrate diverse liver cell type expression and promote a range of metabolic and liver diseases. These findings provide insight into the pathophysiology of liver and other metabolic diseases that are associated with serum liver enzyme concentrations., Serum liver enzymes are used as markers of liver disease, their concentration influenced in part by genetic factors. Here the authors meta-analyse genome-wide association studies on the UK Biobank and BioBank Japan to evaluate the association of three liver enzymes with liver and other metabolic diseases.

Our systematic literature collection and annotation identified 106 chemical drugs and 31 antibodies effective against the infection of at least one human coronavirus (including SARS-CoV, SAR-CoV-2, and MERS-CoV) in vitro or in vivo in an experimental or clinical setting. A total of 163 drug protein targets were identified, and 125 biological processes involving the drug targets were significantly enriched based on a Gene Ontology (GO) enrichment analysis. The Coronavirus Infectious Disease Ontology (CIDO) was used as an ontological platform to represent the anti-coronaviral drugs, chemical compounds, drug targets, biological processes, viruses, and the relations among these entities. In addition to new term generation, CIDO also adopted various terms from existing ontologies and developed new relations and axioms to semantically represent our annotated knowledge. The CIDO knowledgebase was systematically analyzed for scientific insights. To support rational drug design, a “Host-coronavirus interaction (HCI) checkpoint cocktail” strategy was proposed to interrupt the important checkpoints in the dynamic HCI network, and ontologies would greatly support the design process with interoperable knowledge representation and reasoning.

Genetic studies of bone mineral density (BMD) largely have been conducted in European populations. We therefore conducted a meta-analysis of six independent African ancestry cohorts to determine whether previously reported BMD loci identified in European populations were transferable to African ancestry populations. We included nearly 5000 individuals with both genetic data and assessments of BMD. Genotype imputation was conducted using the 1000G reference panel. We assessed single-nucleotide polymorphism (SNP) associations with femoral neck and lumbar spine BMD in each cohort separately, then combined results in fixed effects (or random effects if study heterogeneity was high, I2 index >60) inverse variance weighted meta-analyses. In secondary analyses, we conducted locus-based analyses of rare variants using SKAT-O. Mean age ranged from 12 to 68 years. One cohort included only men and another cohort included only women; the proportion of women in the other four cohorts ranged from 52% to 63%. Of 56 BMD loci tested, one locus, 6q25 (C6orf97, p = 8.87 × 10-4 ), was associated with lumbar spine BMD and two loci, 7q21 (SLC25A13, p = 2.84 × 10-4 ) and 7q31 (WNT16, p = 2.96 × 10-5 ), were associated with femoral neck BMD. Effects were in the same direction as previously reported in European ancestry studies and met a Bonferroni-adjusted p value threshold, the criteria for transferability to African ancestry populations. We also found associations that met locus-specific Bonferroni-adjusted p value thresholds in 11q13 (LRP5, p

Carotid artery atherosclerotic disease (CAAD) is a risk factor for stroke. We used a genome‐wide association (GWAS) approach to discover genetic variants associated with CAAD in participants in the electronic Medical Records and Genomics (eMERGE) Network. We identified adult CAAD cases with unilateral or bilateral carotid artery stenosis and controls without evidence of stenosis from electronic health records at eight eMERGE sites. We performed GWAS with a model adjusting for age, sex, study site, and genetic principal components of ancestry. In eMERGE we found 1793 CAAD cases and 17,958 controls. Two loci reached genome‐wide significance, on chr6 in LPA (rs10455872, odds ratio [OR] (95% confidence interval [CI]) = 1.50 (1.30–1.73), p = 2.1 × 10−8) and on chr7, an intergenic single nucleotide variant (SNV; rs6952610, OR (95% CI) = 1.25 (1.16–1.36), p = 4.3 × 10−8). The chr7 association remained significant in the presence of the LPA SNV as a covariate. The LPA SNV was also associated with coronary heart disease (CHD; 4199 cases and 11,679 controls) in this study (OR (95% CI) = 1.27 (1.13–1.43), p = 5 × 10−5) but the chr7 SNV was not (OR (95% CI) = 1.03 (0.97–1.09), p = .37). Both variants replicated in UK Biobank. Elevated lipoprotein(a) concentrations ([Lp(a)]) and LPA variants associated with elevated [Lp(a)] have previously been associated with CAAD and CHD, including rs10455872. With electronic health record phenotypes in eMERGE and UKB, we replicated a previously known association and identified a novel locus associated with CAAD.

Human genome-wide association studies found SNPs near LYPLAL1 that have sex-specific effects on fat distribution and metabolic traits. To determine whether altering LYPLAL1 affects obesity and metabolic disease we created and characterized a mouse knockout of Lyplal1. Here we show that CRISPR-Cas9 whole-body Lyplal1 knockout (KO) mice fed a high fat, high sucrose (HFHS) diet showed sex-specific differences in weight gain and fat accumulation. Female, not male, KO mice weighed less than WT mice, had reduced body fat percentage, white fat mass, and adipocyte diameter not accounted for by changes in metabolic rate. Female, but not male, KO mice had increased serum triglycerides, decreased aspartate, and alanine aminotransferase. Lyplal1 KO mice of both sexes have reduced liver triglycerides and steatosis. These diet-specific effects resemble the effects of SNPs near LYPLAL1 in humans, suggesting that LYPLAL1 has an evolutionary conserved sex-specific effect on adiposity. This murine model can be used to study this novel gene-by-sex-by-diet interaction to elucidate the metabolic effects of LYPLAL1 on human obesity.

Up to 25% of patients with nonalcoholic fatty liver disease (NAFLD) are not obese but may have a fat or muscle composition that predisposes them to NAFLD. Our aim was to determine whether body composition parameters associate with NAFLD and to identify genetic contributors to this association. This study included two cohorts. The first included 2,249 participants from the Framingham Heart Study who underwent a computed tomography scan to evaluate hepatic steatosis, dual‐energy x‐ray absorptiometry testing to assess body composition, and clinical examination. Body composition parameters were normalized to total body weight. A subset of participants underwent genotyping with an Affymetrix 550K single‐nucleotide polymorphism array. The second cohort, Michigan Genomics Initiative, included 19,239 individuals with genotyping on the Illumina HumanCoreExome v.12.1 array and full electronic health record data. Using sex‐stratified multivariable linear regression, greater central body fat associated with increased hepatic steatosis while greater lower extremity body fat associated with decreased hepatic steatosis. Greater appendicular lean mass was associated with decreased hepatic steatosis in men but not in women. A polygenic risk score for lipodystrophy (regional or global loss of adipose tissue) was associated with increased hepatic steatosis, increased liver fibrosis, and decreased lower extremity fat mass. Conclusion: Greater central body fat associated with increased hepatic steatosis, while greater lower extremity body fat and, in men, greater appendicular lean mass were associated with decreased hepatic steatosis. A genetic risk score for lipodystrophy was associated with NAFLD and liver fibrosis. Our results suggest that buffering of excess energy by peripheral fat and muscle may protect against NAFLD and liver fibrosis in the general population.

Focal adhesion (FA)-stimulated reorganization of the F-actin cytoskeleton regulates cellular size, shape, and mechanical properties. However, FA cross-talk with the intermediate filament cytoskeleton is poorly understood. Genetic ablation of the FA-associated scaffold protein Hic-5 in mouse cancer-associated fibroblasts (CAFs) promoted a dramatic collapse of the vimentin network, which was rescued following EGFP-Hic-5 expression. Vimentin collapse correlated with a loss of detergent-soluble vimentin filament precursors and decreased vimentin S72/S82 phosphorylation. Additionally, fluorescence recovery after photobleaching analysis indicated impaired vimentin dynamics. Microtubule (MT)-associated EB1 tracking and Western blotting of MT posttranslational modifications indicated no change in MT dynamics that could explain the vimentin collapse. However, pharmacological inhibition of the RhoGTPase Cdc42 in Hic-5 knockout CAFs rescued the vimentin collapse, while pan-formin inhibition with SMIFH2 promoted vimentin collapse in Hic-5 heterozygous CAFs. Our results reveal novel regulation of vimentin organization/dynamics by the FA scaffold protein Hic-5 via modulation of RhoGTPases and downstream formin activity.

Simple Summary Studying rare outcomes in cancer is challenging because observation of the rare event may require a very high number of patients or experimental animals. Here, we propose a new, predictive approach to understanding the biological mechanisms underlying such rare events in cancer treatment outcomes. We take as a case-study, the treatment of metastatic, castration-resistant prostate cancer with the live cell vaccine, sipuleucel-t (Provenge). Only a fraction of patients benefited from Provenge; that is, clinical success is a rare event. It remains an open question why Provenge conferred such a modest survival benefit. Our modeling paradigm captures the inherent heterogeneity that characterizes individuals in a population, and provides an explanation for the observed clinical outcomes of treatment with Provenge. Our approach readily generalizes to a range of emerging cancer immunotherapies, and more generally, to predicting and understanding how a population responds to any intervention targeting a human disease. Abstract Sipuleucel-T (Provenge) is the first live cell vaccine approved for advanced, hormonally refractive prostate cancer. However, survival benefit is modest and the optimal combination or schedule of sipuleucel-T with androgen depletion remains unknown. We employ a nonlinear dynamical systems approach to modeling the response of hormonally refractive prostate cancer to sipuleucel-T. Our mechanistic model incorporates the immune response to the cancer elicited by vaccination, and the effect of androgen depletion therapy. Because only a fraction of patients benefit from sipuleucel-T treatment, inter-individual heterogeneity is clearly crucial. Therefore, we introduce our novel approach, Standing Variations Modeling, which exploits inestimability of model parameters to capture heterogeneity in a deterministic model. We use data from mouse xenograft experiments to infer distributions on parameters critical to tumor growth and to the resultant immune response. Sampling model parameters from these distributions allows us to represent heterogeneity, both at the level of the tumor cells and the individual (mouse) being treated. Our model simulations explain the limited success of sipuleucel-T observed in practice, and predict an optimal combination regime that maximizes predicted efficacy. This approach will generalize to a range of emerging cancer immunotherapies.

BACKGROUND: High macular pigment optical density (MPOD) has been associated with improved eye health and better cognitive functions. Genetic variations have been associated with MPOD in adults. However, these associations between genetic variations and MPOD have not been studied in children. OBJECTIVES: This was a secondary analysis of the FK2 (Fitness Improves Thinking in Kids 2) trial (n = 134, 41% male). The aim was to determine differences in MPOD among children (aged 7–9 y) based on genetic variants that either are biologically relevant to lutein (L) and zeaxanthin (Z) accumulation or have been associated with MPOD in adults. METHODS: MPOD was measured using customized heterochromatic flicker photometry via a macular densitometer. DXA was used to assess whole-body and visceral adiposity. DNA was extracted from saliva samples and was genotyped for 26 hypothesis-driven single nucleotide polymorphisms and 75 ancestry-informative markers (AIMs). Habitual diet history was obtained via 3-d food logs completed by parents (n = 88). General linear models were used to compare MPOD between different genotypes. Principal component analysis was performed for the AIMs to account for ethnic heterogeneity. RESULTS: Children carrying ≥1 minor allele on β-carotene-15,15′-monooxygenase (BCO1)-rs7501331 (T allele) (P = 0.045), cluster of differentiation 36(CD36)-rs1527483 (T allele) (P = 0.038), or CD36-rs3173798 (C allele) (P = 0.001) had significantly lower MPOD (range: 14.1%–26.4%) than those who were homozygotes for the major alleles. MPOD differences based on CD36-rs3173798 genotypes persisted after adjustment for dietary L and Z intake. CONCLUSIONS: The findings indicate that genetic variations of CD36 and BCO1 contribute to MPOD in children. The influence of genetic variation in CD36-rs3173798 persisted after adjusting for variation in dietary intake. This trial was registered at clinicaltrials.gov as NCT01619826.

Purpose: Lineage plasticity in prostate cancer—most commonly exemplified by loss of androgen receptor (AR) signaling and a switch from a luminal to alternate differentiation program—is now recognized as a treatment resistance mechanism. Lineage plasticity is a spectrum, but neuroendocrine prostate cancer (NEPC) is the most virulent example. Currently, there are limited treatments for NEPC. Moreover, the incidence of treatment-emergent NEPC (t-NEPC) is increasing in the era of novel AR inhibitors. In contradistinction to de novo NEPC, t-NEPC tumors often express the AR, but AR's functional role in t-NEPC is unknown. Furthermore, targetable factors that promote t-NEPC lineage plasticity are also unclear. Experimental Design: Using an integrative systems biology approach, we investigated enzalutamide-resistant t-NEPC cell lines and their parental, enzalutamide-sensitive adenocarcinoma cell lines. The AR is still expressed in these t-NEPC cells, enabling us to determine the role of the AR and other key factors in regulating t-NEPC lineage plasticity. Results: AR inhibition accentuates lineage plasticity in t-NEPC cells—an effect not observed in parental, enzalutamide-sensitive adenocarcinoma cells. Induction of an AR-repressed, lineage plasticity program is dependent on activation of the transcription factor E2F1 in concert with the BET bromodomain chromatin reader BRD4. BET inhibition (BETi) blocks this E2F1/BRD4-regulated program and decreases growth of t-NEPC tumor models and a subset of t-NEPC patient tumors with high activity of this program in a BETi clinical trial. Conclusions: E2F1 and BRD4 are critical for activating an AR-repressed, t-NEPC lineage plasticity program. BETi is a promising approach to block this program.

The pharmacogenomic network responsible for the rapid antidepressant action of ketamine and concomitant adverse events in patients has been poorly defined. Integrative, multi-scale biological data analytics helps explain ketamine’s action. Using a validated computational pipeline, candidate ketamine-response genes and regulatory RNAs from published literature, binding affinity studies, and single nucleotide polymorphisms (SNPs) from genomewide association studies (GWAS), we identified 108 SNPs associated with 110 genes and regulatory RNAs. All of these SNPs are classified as enhancers, and additional chromatin interaction mapping in human neural cell lines and tissue shows enhancer-promoter interactions involving other network members. Pathway analysis and gene set optimization identified three composite sub-networks within the broader ketamine pharmacogenomic network. Expression patterns of ketamine network genes within the postmortem human brain are concordant with ketamine neurocircuitry based on the results of 24 published functional neuroimaging studies. The ketamine pharmacogenomic network is enriched in forebrain regions known to be rapidly activated by ketamine, including cingulate cortex and frontal cortex, and is significantly regulated by ketamine (p=6.26E-33; Fisher’s exact test). The ketamine pharmacogenomic network can be partitioned into distinct enhancer sub-networks associated with: (1) glutamate neurotransmission, chromatin remodeling, smoking behavior, schizophrenia, pain, nausea, vomiting, and post-operative delirium; (2) neuroplasticity, depression, and alcohol consumption; and (3) pharmacokinetics. The component sub-networks explain the diverse action mechanisms of ketamine and its analogs. These results may be useful for optimizing pharmacotherapy in patients diagnosed with depression, pain or related stress disorders.One Sentence SummaryThe ketamine network in the human brain consists of sub-networks associated with glutamate neurotransmission, neuroplasticity, and pharmacokinetics.

BACKGROUND AND AIMS: Cirrhosis is characterized by extensive fibrosis of the liver and is a major cause of liver-related mortality. Cirrhosis is partially heritable but genetic contributions to cirrhosis have not been systemically explored. Here, we carry out association analyses with cirrhosis in two large biobanks and determine the effects of cirrhosis associated variants on multiple human disease/traits. METHODS: We carried out a genome-wide association analysis of cirrhosis as a diagnosis in UK BioBank (UKBB; 1088 cases vs. 407 873 controls) and then tested top-associating loci for replication with cirrhosis in a hospital-based cohort from the Michigan Genomics Initiative (MGI; 875 cases of cirrhosis vs. 30 346 controls). For replicating variants or variants previously associated with cirrhosis that also affected cirrhosis in UKBB or MGI, we determined single nucleotide polymorphism effects on all other diagnoses in UKBB (PheWAS), common metabolic traits/diseases and serum/plasma metabolites. RESULTS: Unbiased genome-wide association study identified variants in/near PNPLA3 and HFE, and candidate variant analysis identified variants in/near TM6SF2, MBOAT7, SERPINA1, HSD17B13, STAT4 and IFNL4 that reproducibly affected cirrhosis. Most affected liver enzyme concentrations and/or aspartate transaminase-to-platelet ratio index. PheWAS, metabolic trait and serum/plasma metabolite association analyses revealed effects of these variants on lipid, inflammatory and other processes including new effects on many human diseases and traits. CONCLUSIONS: We identified eight loci that reproducibly associate with population-based cirrhosis and define their diverse effects on human diseases and traits.

Diverticular disease is common and has a high morbidity. Treatments are limited owing to the poor understanding of its pathophysiology. Here, to elucidate its etiology, we performed a genome-wide association study of diverticular disease (27,444 cases; 382,284 controls) from the UK Biobank and tested for replication in the Michigan Genomics Initiative (2,572 cases; 28,649 controls). We identified 42 loci associated with diverticular disease; 39 of these loci are novel. Using data-driven expression-prioritized integration for complex traits (DEPICT), we show that genes in these associated regions are significantly enriched for expression in mesenchymal stem cells and multiple connective tissue cell types and are co-expressed with genes that have a role in vascular and mesenchymal biology. Genes in these associated loci have roles in immunity, extracellular matrix biology, cell adhesion, membrane transport and intestinal motility. Phenome-wide association analysis of the 42 variants shows a common etiology of diverticular disease with obesity and hernia. These analyses shed light on the genomic landscape of diverticular disease. Genome-wide analyses identify 42 risk loci for diverticular disease, 39 of which are new. Genes in associated regions are enriched for expression in connective tissue cell types and are coexpressed with genes involved in vascular and mesenchymal biology.

Numerous pharmacogenetic clinical guidelines and recommendations have been published, but barriers have hindered the clinical implementation of pharmacogenetics. The Translational Pharmacogenetics Program (TPP) of the National Institutes of Health (NIH) Pharmacogenomics Research Network was established in 2011 to catalog and contribute to the development of pharmacogenetic implementations at eight US healthcare systems, with the goal to disseminate real-world solutions for the barriers to clinical pharmacogenetic implementation. The TPP collected and normalized pharmacogenetic implementation metrics through June 2015, including gene-drug pairs implemented, interpretations of alleles and diplotypes, numbers of tests performed and actionable results, and workflow diagrams. TPP participant institutions developed diverse solutions to overcome many barriers, but the use of Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines provided some consistency among the institutions. The TPP also collected some pharmacogenetic implementation outcomes (scientific, educational, financial, and informatics), which may inform healthcare systems seeking to implement their own pharmacogenetic testing programs.

Background: Fracture is a complex trait, affected by both genetic and environmental factors. A meta-analysis of genome-wide association studies (GWASs) identified multiple bone mineral density (BMD) and fracture-associated loci.Objective: We conducted a study to evaluate whether fracture genetic risk score (Fx-GRS) and bone mineral density genetic risk score (BMD-GRS) modify the association between the intake of calcium with vitamin D (CaD) and fracture risk.Design: Data from 5823 white postmenopausal women from the Women's Health Initiative CaD randomized trial were included. Participants received 1000 mg elemental Ca with 400 IU vitamin D3/d or placebo (median follow-up: 6.5 y). Total fracture was defined as first fracture of any type. We computed the Fx-GRS with 16 fracture- and BMD-associated variants, and the BMD-GRS with 50 BMD-associated variants. We used Cox regression and a case-only approach to test for multiplicative interaction. Additive interaction was assessed with the relative excess risk due to interaction (RERI). We analyzed genetic risk score as a continuous variable and a categorical variable based on quartile (quartile 1, quartiles 2-3, and quartile 4).Results: We observed no interaction between the Fx-GRS and CaD on fracture risk; however, we observed a significant multiplicative interaction between the BMD-GRS and CaD assignment (P-interaction = 0.01). In addition, there was a significant negative additive interaction between placebo assignment and higher BMD-GRS: quartiles 2-3, PRERI = 0.03; quartile 4, PRERI = 0.03. In a stratified analysis, the protective effect of CaD on fracture risk was observed in women in the lowest BMD-GRS quartile (HR: 0.60, 95% CI: 0.44, 0.81) but not in women with a higher BMD-GRS.Conclusions: We observed significant effects of CaD intake on fracture risk only in women with the lowest genetic predisposition to low BMD. Future large-scale studies with functional characterization of GWAS findings are warranted to assess the utility of genetic risk score in analysis of risks and benefits of CaD for bone.

Context Glycogen storage diseases are rare. Increased glycogen in the liver results in increased attenuation. Objective Investigate the association and function of a noncoding region associated with liver attenuation but not histologic nonalcoholic fatty liver disease. Design Genetics of Obesity-associated Liver Disease Consortium. Setting Population-based. Main Outcome Computed tomography measured liver attenuation. Results Carriers of rs4841132-A (frequency 2%-19%) do not show increased hepatic steatosis; they have increased liver attenuation indicative of increased glycogen deposition. rs4841132 falls in a noncoding RNA LOC157273 ~190 kb upstream of PPP1R3B. We demonstrate that rs4841132-A increases PPP1R3B through a cis genetic effect. Using CRISPR/Cas9 we engineered a 105-bp deletion including rs4841132-A in human hepatocarcinoma cells that increases PPP1R3B, decreases LOC157273, and increases glycogen perfectly mirroring the human disease. Overexpression of PPP1R3B or knockdown of LOC157273 increased glycogen but did not result in decreased LOC157273 or increased PPP1R3B, respectively, suggesting that the effects may not all occur via affecting RNA levels. Based on electronic health record (EHR) data, rs4841132-A associates with all components of the metabolic syndrome (MetS). However, rs4841132-A associated with decreased low-density lipoprotein (LDL) cholesterol and risk for myocardial infarction (MI). A metabolic signature for rs4841132-A includes increased glycine, lactate, triglycerides, and decreased acetoacetate and beta-hydroxybutyrate. Conclusions These results show that rs4841132-A promotes a hepatic glycogen storage disease by increasing PPP1R3B and decreasing LOC157273. rs4841132-A promotes glycogen accumulation and development of MetS but lowers LDL cholesterol and risk for MI. These results suggest that elevated hepatic glycogen is one cause of MetS that does not invariably promote MI.

BackgroundIn comparison to the non-pregnant state, the first trimester of pregnancy is characterized by systemic adaptation of the mother. The extent to which these adaptive processes are reflected in the maternal blood metabolome is not well characterized.ObjectiveTo determine the differences between the plasma metabolome of non-pregnant and pregnant women before 16 weeks gestation.Study designThis study included plasma samples from 21 non-pregnant women and 50 women with a normal pregnancy (8-16 weeks of gestation). Combined measurements by ultrahigh performance liquid chromatography/tandem mass spectrometry and by gas chromatography/mass spectrometry generated molecular abundance measurements for each sample. Molecular species detected in at least 10 samples were included in the analysis. Differential abundance was inferred based on false discovery adjusted p-values (FDR) from Mann-Whitney-Wilcoxon U tests ResultsOverall, 637 small molecules met the inclusion criteria and were tested for association with pregnancy; 44% (281/637) of small molecules had significantly different abundance, of which 81% (229/281) were less abundant in pregnant than in non-pregnant women. Eight percent (14/169) of the metabolites that remained significant in the adjusted analysis also changed as a function of gestational age. A pathway analysis revealed enrichment in steroid metabolites related to sex hormones, caffeine metabolites, lysolipids, dipeptides, and polypeptide bradykinin derivatives (all, FDR < 0.1).ConclusionsThis high-throughput mass spectrometry study identified: 1) differences between pregnant vs. non-pregnant women in the abundance of 44% of the profiled plasma metabolites, including known and novel molecules and pathways; and 2) specific metabolites that changed with gestational age.

Dementia with Lewy Bodies (DLB) is the second most common neurodegenerative disorder in the elderly. The development and progression of DLB remain unclear. In this study we used next generation sequencing to assess RNA expression profiles and cellular processes associated with DLB in the anterior cingulate cortex, a brain region affected by DLB pathology. The expression measurements were made in autopsy brain tissues from 8 DLB subjects and 10 age-matched controls using AmpliSeq technology with ion torrent sequencing. The analysis of RNA expression profiles revealed 490 differentially expressed genes, among which 367 genes were down-regulated and 123 were up-regulated. Functional enrichment analysis of genes differentially expressed in DLB indicated downregulation of genes associated with myelination, neurogenesis, and regulation of nervous system development. miRNA binding sites enriched in these mRNAs yielded a list of candidate miRNAs participating in DLB pathophysiology. Our study provides a comprehensive picture of gene expression landscape in DLB, identifying key cellular processes associated with DLB pathology.

The accumulation of excess fat in the liver (hepatic steatosis) in the absence of heavy alcohol consumption causes nonalcoholic fatty liver disease (NAFLD), which has become a global epidemic. Identifying metabolic risk factors that interact with the genetic risk of NAFLD is important for reducing disease burden. We tested whether serum glucose, insulin, insulin resistance, triglyceride (TG), low-density lipoprotein cholesterol, high-density lipoprotein cholesterol, body mass index (BMI), and waist-to-hip ratio adjusted for BMI interact with genetic variants in or near the patatin-like phospholipase domain containing 3 (PNPLA3) gene, the glucokinase regulatory protein (GCKR) gene, the neurocan/transmembrane 6 superfamily member 2 (NCAN/TM6SF2) gene, and the lysophospholipase-like 1 (LYPLAL1) gene to exacerbate hepatic steatosis, estimated by liver attenuation. We performed association analyses in 10 population-based cohorts separately and then meta-analyzed results in up to 14,751 individuals (11,870 of European ancestry and 2,881 of African ancestry). We found that PNPLA3-rs738409 significantly interacted with insulin, insulin resistance, BMI, glucose, and TG to increase hepatic steatosis in nondiabetic individuals carrying the G allele. Additionally, GCKR-rs780094 significantly interacted with insulin, insulin resistance, and TG. Conditional analyses using the two largest European ancestry cohorts in the study showed that insulin levels accounted for most of the interaction of PNPLA3-rs738409 with BMI, glucose, and TG in nondiabetic individuals. Insulin, PNPLA3-rs738409, and their interaction accounted for at least 8% of the variance in hepatic steatosis in these two cohorts. Conclusion: Insulin resistance, either directly or through the resultant elevated insulin levels, more than other metabolic traits, appears to amplify the PNPLA3-rs738409-G genetic risk for hepatic steatosis. Improving insulin resistance in nondiabetic individuals carrying PNPLA3-rs738409-G may preferentially decrease hepatic steatosis.

Non-alcoholic fatty liver disease (NAFLD) is a common cause of chronic liver disease. A single nucleotide polymorphism (SNPs), rs6834314, was associated with serum liver enzymes in the general population, presumably reflecting liver fat or injury. We studied rs6834314 and its nearest gene, HSD17B13 (17-beta hydroxysteroid dehydrogenase 13) to identify associations with histological features of NAFLD, and to characterize the functional role of HSD17B13 in NAFLD pathogenesis. The minor allele of rs6834314 was significantly associated with increased steatosis, but decreased inflammation, ballooning, Mallory-Denk bodies, and liver enzyme levels in 768 adult Caucasians with biopsy-proven NAFLD, and with cirrhosis in the general population. We found two plausible causative variants in the HSD17B13 gene. rs72613567, a splice-site SNP in high linkage with rs6834314 (r(2)=0.94) generates novel splice variants and shows a similar pattern of association with NAFLD histology. Its minor allele generates simultaneous expression of exon 6-skipping and G-nucleotide insertion variants. Another SNP, rs62305723 (encoding a P260S mutation) is significantly associated with decreased ballooning and inflammation. Hepatic expression of HSD17B13 is 5.9-fold higher (p=0.003) in patients with NAFLD. HSD17B13 is targeted to lipid droplets, requiring the conserved AA22–28 sequence and AA71–106 region. The protein has a retinol dehydrogenase (RDH) activity, with enzymatic activity dependent on lipid droplet targeting and co-factor binding site. The exon 6-deletion, G-insertion, and newly-described naturally-occurring P260S mutation, all confer loss of enzymatic activity. In conclusion, we demonstrate the association of variants in HSD17B13 with specific features of NAFLD histology, and identify the enzyme as a lipid droplet-associated RDH. Our data suggest that HSD17B13 plays a role in NAFLD through its enzymatic activity.

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course.

Human alveolar macrophages (HAM) are primary bacterial niche and immune response cells during Mycobacterium tuberculosis (M.tb) infection, and human blood monocyte-derived macrophages (MDM) are a model for investigating M.tb-macrophage interactions. Here, we use a targeted RNA-Seq method to measure transcriptome-wide changes in RNA expression patterns of freshly obtained HAM (used within 6 h) and 6 day cultured MDM upon M.tb infection over time (2, 24 and 72 h), in both uninfected and infected cells from three donors each. The Ion AmpliSeq™ Transcriptome Human Gene Expression Kit (AmpliSeq) uses primers targeting 18,574 mRNAs and 2,228 non-coding RNAs (ncRNAs) for a total of 20,802 transcripts. AmpliSeqTM yields highly precise and reproducible gene expression profiles (R2 >0.99). Taking advantage of AmpliSeq's reproducibility, we establish well-defined quantitative RNA expression patterns of HAM versus MDM, including significant M.tb-inducible genes, in networks and pathways that differ in part between MDM and HAM. A similar number of expressed genes are detected at all time-points between uninfected MDM and HAM, in common pathways including inflammatory and immune functions, but canonical pathway differences also exist. In particular, at 2 h, multiple genes relevant to the immune response are preferentially expressed in either uninfected HAM or MDM, while the HAM RNA profiles approximate MDM profiles over time in culture, highlighting the unique RNA expression profile of freshly obtained HAM. MDM demonstrate a greater transcriptional response than HAM upon M.tb infection, with 2 to >10 times more genes up- or down-regulated. The results identify key genes involved in cellular responses to M.tb in two different human macrophage types. Follow-up bioinformatics analysis indicates that approximately 30% of response genes have expression quantitative trait loci (eQTLs in GTEx), common DNA variants that can influence host gene expression susceptibility or resistance to M.tb, illustrated with the TREM1 gene cluster and IL-10.

The 4D nucleome has the potential to render challenges in neuropsychiatric pharmacogenomics more tractable. The epigenome roadmap consortium has demonstrated the critical role that noncoding regions of the human genome play in determination of human phenotype. Chromosome conformation capture methods have revealed the 4D organization of the nucleus, bringing interactions between distant regulatory elements into close spatial proximity in a periodic manner. These functional interactions have the potential to elucidate mechanisms of CNS drug response and side effects that previously have been unrecognized. This perspective assesses recent advances likely to reveal novel pharmacodynamic regulatory pathways in human brain, charting a future new avenue of pharmacogenomics research, using the spatial and temporal architecture of the human epigenome as its foundation.

Genetic factors strongly influence risk of common human diseases and treatment outcomes but the causative variants remain largely unknown; this gap has been called the ‘missing heritability’. We propose several hypotheses that in combination have the potential to narrow the gap. First, given a multi-stage path from wellness to disease, we propose that common variants under positive evolutionary selection represent normal variation and gate the transition between wellness and an ‘off-well’ state, revealing adaptations to changing environmental conditions. In contrast, genome-wide association studies (GWAS) focus on deleterious variants conveying disease risk, accelerating the path from off-well to illness and finally specific diseases, while common ‘normal’ variants remain hidden in the noise. Second, epistasis (dynamic gene-gene interactions) likely assumes a central role in adaptations and evolution; yet, GWAS analyses currently are poorly designed to reveal epistasis. As gene regulation is germane to adaptation, we propose that epistasis among common normal regulatory variants, or between common variants and less frequent deleterious variants, can have strong protective or deleterious phenotypic effects. These gene-gene interactions can be highly sensitive to environmental stimuli and could account for large differences in drug response between individuals. Residing largely outside the protein-coding exome, common regulatory variants affect either transcription of coding and non-coding RNAs (regulatory SNPs, or rSNPs) or RNA functions and processing (structural RNA SNPs, or srSNPs). Third, with the vast majority of causative variants yet to be discovered, GWAS rely on surrogate markers, a confounding factor aggravated by the presence of more than one causative variant per gene and by epistasis. We propose that the confluence of these factors may be responsible to large extent for the observed heritability gap.

Approximately 20% of all HIV-1 mother-to-child transmission (MTCT) occurs in utero (IU). In a chronic HIV infection, HIV-1 exists as a complex swarm of genetic variants, and following IU MTCT, viral genomic diversity is restricted through a mechanism that remains to be described. The 5' U3R region of the HIV-1 long terminal repeat (LTR) contains multiple transcription factor (TF) binding sites and regulates viral transcription. In this study, we tested the hypothesis that sequence polymorphisms in the U3R region of LTR are associated with IU MTCT. To this end, we used single template amplification to isolate 517 U3R sequences from maternal, placental, and infant plasma derived from 17 HIV-infected Malawian women: eight whose infants remained HIV uninfected (NT) and nine whose infants became HIV infected IU. U3R sequences show pairwise diversities ranging from 0.2% to 2.3%. U3R sequences from one participant contained two, three, or four putative NF-κB binding sites. Phylogenetic reconstructions indicated that U3R sequences from eight of nine IU participants were consistent with placental compartmentalization of HIV-1 while only one of eight NT cases was consistent with such compartmentalization. Specific TF sequence polymorphisms were not significantly associated with IU MTCT. To determine if replication efficiency of the U3R sequences was associated with IU MTCT, we cloned 90 U3R sequences and assayed promoter activity in multiple cell lines. Although we observed significant, yet highly variable promoter activity and TAT induction of promoter activity in the cell lines tested, there was no association between measured promoter activity and MTCT status. Thus, we were unable to detect a promoter genotype or phenotype associated with IU MTCT.

Background: Whole genome analyses of peripheral blood has demonstrated that acquired somatic mutations in peripheral blood also known as clonal hematopoiesis of indeterminate potential (CHIP) is present in up to 10% of individuals older than 60 years and associated with increased risk of cardiovascular mortality and hematologic malignancies (Jaiswal et al, Genovese et al, NEJM 2014). We have demonstrated that CHIP is associated specifically with increased risk of leukemia (Desai et al, Nat. Medicine 2018). CHIP has also been detected in 25% of individuals with concomitant advanced solid malignancies using targeted deep sequencing (Coombs et al, 2017). However, the role of antecedent CHIP in risk of solid malignancies has not been established due to lack of non-cancer controls and absence of CHIP data before the diagnosis of cancer in the published literature. The role of antecedent CHIP in cancer specific and non-cancer specific mortality after diagnosis of cancer is also not known. Methods: We analyzed whole genome sequencing data (30X coverage) from 10,089 participants in the Women's Health Initiative (WHI) enrolled in the Trans-Omics for Precision Medicine (TOPMED) consortium to assess the relationship between antecedent CHIP measured at baseline entry into the study, and cancer risk and mortality after cancer. Cox proportional hazards regression models were used to evaluate the relationship between CHIP and cancer risk, all-cause mortality as well as mortality due to cancer, cardiovascular disease (CVD) and other causes. All statistical tests are two-sided with an alpha level of 0.05. Results: Of 10,089 eligible participants, the overall rate of CHIP mutations was 9.3 % and the most common CHIP mutations included DNMT3A (55%), TET2 (19%), and ASXL1 (7.1%). CHIP was associated with both current and past history of smoking and the best-fit model suggested a dose dependent relationship by pack years. Over 13±6 years of follow-up, there were 2,337 women diagnosed with at least one cancer including (at 7±5 years of follow-up) 774 with breast cancer, 319 with colon, 355 with lung and 282 with hematologic malignancies. Among patients with cancer, there were 813 cancer specific deaths and 3,946 non-cancer related deaths. Among participants with and without a previous history of cancer, CHIP was detected in 8.4% and 9.3% of participants respectively. PPM1D or TP53 (n=50, of whom 8 had a previous history of malignancy) mutations were detected in 3.6% and 1.7 % of total participants with CH. 29% of participants with TP53 or PPM1D CHIP developed a future cancer , compared to 26% in non TP53/PPM1D CHIP and 21% without CHIP. There was no significant difference in TP53 CHIP in patients with or without previous history of cancer. After excluding 697 participants with a previous history of malignancy, antecedent CHIP (prior to the diagnosis of cancer) was not significantly associated with risk of solid malignancy overall (Hazards Ratio (HR), 95% confidence interval (CI) 1.07, (0.94−1.23) p=0.31). However, there was a borderline increased risk of breast cancer (HR, 95% CI, 1.23, (1.003−1.25) P=0.04) and no association with lung (HR 1.23, P=0.15) or colon cancer (HR 0.96, P=0.86). Among participants with solid malignancies, CHIP was associated with an increased cancer specific mortality (HR, 95% CI, 1.24 (1.02−1.51), p=0.03) but no association with mortality due to CVD post diagnosis of solid malignancy (HR 0.63, P=0.51) was seen. Any antecedent CHIP mutation was associated with increased risk of hematological malignancies (HR, 95% CI, 1.76, (1.25−2.48), p Conclusion: This is the first report to relate antecedent CHIP with solid tumor risk and mortality in a multi-site study of post-menopausal women. Antecedent CHIP was associated with increased solid cancer specific mortality but not with risk of solid malignancy or CVD mortality post cancer diagnosis. It is possible that deeper, targeted sequencing may identify an association between antecedent CHIP and the incidence of solid tumors; further investigation is warranted. Disclosures Desai: Sanofi: Consultancy; Astellas: Honoraria; Astex: Research Funding; Cellerant: Consultancy; Celgene: Consultancy. Lee:Helsinn: Consultancy; Jazz Pharmaceuticals, Inc: Consultancy; Roche Molecular Systems: Consultancy; AstraZeneca Pharmaceuticals: Consultancy; Karyopharm Therapeutics: Consultancy; Ai Therapeutics: Research Funding. Ritchie:Celgene, Incyte, Novartis, Pfizer: Consultancy; Ariad, Celgene, Incyte, Novartis: Speakers Bureau; AStella, Bristol-Myers Squibb, Novartis, NS Pharma, Pfizer: Research Funding; Celgene, Novartis: Other: travel support; Jazz Pharmaceuticals: Research Funding; Celgene: Other: Advisory board; Pfizer: Other: Advisory board, travel support; agios: Other: Advisory board; Tolero: Other: Advisory board; Genentech: Other: Advisory board. Guzman:Samus Therapeutics: Patents & Royalties: intellectual rights to the PU-FITC assay; SeqRx: Consultancy; Cellectis: Research Funding. Roboz:AbbVie: Consultancy, Membership on an entity's Board of Directors or advisory committees; Actinium: Consultancy, Membership on an entity's Board of Directors or advisory committees; Agios: Consultancy, Membership on an entity's Board of Directors or advisory committees; Amphivena: Consultancy, Membership on an entity's Board of Directors or advisory committees; Argenx: Consultancy, Membership on an entity's Board of Directors or advisory committees; Astex: Consultancy, Membership on an entity's Board of Directors or advisory committees; Astellas: Consultancy, Membership on an entity's Board of Directors or advisory committees; Bayer: Consultancy, Membership on an entity's Board of Directors or advisory committees; Celgene: Consultancy, Membership on an entity's Board of Directors or advisory committees; Celltrion: Consultancy, Membership on an entity's Board of Directors or advisory committees; Daiichi Sankyo: Consultancy, Membership on an entity's Board of Directors or advisory committees; Eisai: Consultancy, Membership on an entity's Board of Directors or advisory committees; Janssen: Consultancy, Membership on an entity's Board of Directors or advisory committees; Jazz: Consultancy, Membership on an entity's Board of Directors or advisory committees; Novartis: Consultancy, Membership on an entity's Board of Directors or advisory committees; MEI Pharma: Consultancy, Membership on an entity's Board of Directors or advisory committees; Orsenix: Consultancy, Membership on an entity's Board of Directors or advisory committees; Otsuka: Consultancy, Membership on an entity's Board of Directors or advisory committees; Pfizer: Consultancy, Membership on an entity's Board of Directors or advisory committees; Roche/Genentech: Consultancy, Membership on an entity's Board of Directors or advisory committees; Sandoz: Consultancy, Membership on an entity's Board of Directors or advisory committees; Takeda: Consultancy, Membership on an entity's Board of Directors or advisory committees; Trovagene: Consultancy, Membership on an entity's Board of Directors or advisory committees.

It has been proposed that supertree approaches should be applied to large multilocus datasets to achieve computational tractability. Large datasets such as those derived from phylogenomics studies can be broken into many locus-specific tree searches and the resulting trees can be stitched together via a supertree method. Using simulated data, workers have reported that they can rapidly construct a supertree that is comparable to the results of heuristic tree search on the entire dataset. To test this assertion with organismal data, we compare tree length under the parsimony criterion and computational time for 20 multilocus datasets using supertree (SuperFine and SuperTriplets) and supermatrix (heuristic search in TNT) approaches. Tree length and computational times were compared among methods using the Wilcoxon matched-pairs signed rank test. Supermatrix searches produced significantly shorter trees than either supertree approach (SuperFine or SuperTriplets; P 0.0002 in both cases). Moreover, the processing time of supermatrix search was significantly lower than SuperFine+locus-specific search (P 0.01) but roughly equivalent to that of SuperTriplets+locus-specific search (P 0.4, not significant). In conclusion, we show by using real rather than simulated data that there is no basis, either in time tractability or in tree length, for use of supertrees over heuristic tree search using a supermatrix for phylogenomics.

Background Alterations in gene expression are key events in disease etiology and risk. Poor reproducibility in detecting differentially expressed genes across studies suggests individual genes may not be sufficiently informative for complex diseases, such as myocardial infarction (MI). Rather, dysregulation of the ‘molecular network’ may be critical for pathogenic processes. Such a dynamic network can be built from pairwise non-linear interactions. Results We investigate non-linear interactions represented in mRNA expression profiles that integrate genetic background and environmental factors. Using logistic regression, we test the association of individual GWAS-based candidate genes and non-linear interaction terms (between these mRNA expression levels) with MI. Based on microarray data in CATHGEN (CATHeterization in GENetics) and FHS (Framingham Heart Study), we find individual genes and pairs of mRNAs, encoded by 41 MI candidate genes, with significant interaction terms in the logistic regression model. Two pairs replicate between CATHGEN and FHS (CNNM2|GUCY1A3 and CNNM2|ZEB2). Analysis of RNAseq data from GTEx (Genotype-Tissue Expression) shows that 20 % of these disease-associated RNA pairs are co-expressed, further prioritizing significant interactions. Because edges in sparse co-expression networks formed solely by the 41 candidate genes are unlikely to represent direct physical interactions, we identify additional RNAs as links between network pairs of candidate genes. This approach reveals additional mRNAs and interaction terms significant in the context of MI, for example, the path CNNM2|ACSL5|SCARF1|GUCY1A3, characterized by the common themes of magnesium and lipid processing. Conclusions The results of this study support a role for non-linear interactions between genes in MI and provide a basis for further study of MI systems biology. mRNA expression profiles encoded by a limited number of candidate genes yield sparse networks of MI-relevant interactions that can be expanded to include additional candidates by co-expression analysis. The non-linear interactions observed here inform our understanding of the clinical relevance of gene-gene interactions in the pathophysiology of MI, while providing a new strategy in developing clinical biomarker panels. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3075-6) contains supplementary material, which is available to authorized users.