Epigenetic MLH1 silencing concurs with mismatch repair deficiency in sporadic, naturally occurring colorectal cancer in rhesus macaques.

BACKGROUND: Naturally occurring colorectal cancers (CRC) in rhesus macaques share many features with their human counterparts and are useful models for cancer immunotherapy; but mechanistic data are lacking regarding the comparative molecular pathogenesis of these cancers. METHODS: We conducted state-of-the-art imaging including CT and PET, clinical assessments, and pathological review of 24 rhesus macaques with naturally occurring CRC. Additionally, we molecularly characterized these tumors utilizing immunohistochemistry (IHC), microsatellite instability assays, DNAseq, transcriptomics, and developed a DNA methylation-specific qPCR assay for MLH1, CACNA1G, CDKN2A, CRABP1, and NEUROG1, human markers for CpG island methylator phenotype (CIMP). We furthermore employed Monte-Carlo simulations to in-silico model alterations in DNA topology in transcription-factor binding site-rich promoter regions upon experimentally demonstrated DNA methylation. RESULTS: Similar cancer histology, progression patterns, and co-morbidities could be observed in rhesus as reported for human CRC patients. IHC identified loss of MLH1 and PMS2 in all cases, with functional microsatellite instability. DNA sequencing revealed the close genetic relatedness to human CRCs, including a similar mutational signature, chromosomal instability, and functionally-relevant mutations affecting KRAS (G12D), TP53 (R175H, R273*), APC, AMER1, ALK, and ARID1A. Interestingly, MLH1 mutations were rarely identified on a somatic or germline level. Transcriptomics not only corroborated the similarities of rhesus and human CRCs, but also demonstrated the significant downregulation of MLH1 but not MSH2, MSH6, or PMS2 in rhesus CRCs. Methylation-specific qPCR suggested CIMP-positivity in 9/16 rhesus CRCs, but all 16/16 exhibited significant MLH1 promoter hypermethylation. DNA hypermethylation was modelled to affect DNA topology, particularly propeller twist and roll profiles. Modelling the DNA topology of a transcription factor binding motif (TFAP2A) in the MLH1 promoter that overlapped with a methylation-specific probe, we observed significant differences in DNA topology upon experimentally shown DNA methylation. This suggests a role of transcription factor binding interference in epigenetic silencing of MLH1 in rhesus CRCs. CONCLUSIONS: These data indicate that epigenetic silencing suppresses MLH1 transcription, induces the loss of MLH1 protein, abrogates mismatch repair, and drives genomic instability in naturally occurring CRC in rhesus macaques. We consider this spontaneous, uninduced CRC in immunocompetent, treatment-naïve rhesus macaques to be a uniquely informative model for human CRC.

Prognostic attributes of immune signatures in soft tissue sarcomas show differential dependencies on tumor mutational burden.

BACKGROUND: Cellular and intrinsic markers of sarcoma immunogenicity are poorly understood. To gain insight into whether tumor-immune interactions correlate with clinical aggressiveness, the authors examined the prognostic significance of immune gene signatures in combination with tumor mutational burden (TMB) and cancer-testis antigen (CTA) expression. METHODS: RNA sequencing and clinical data of 259 soft tissue sarcomas from The Cancer Genome Atlas project were used to investigate associations between published immune gene signatures and patient overall survival (OS) in the contexts of TMB, as computed from whole-exome sequencing data, and CTA gene expression. Multivariate Cox proportional hazards regression models and log-rank tests were used to assess survival associations. RESULTS: Immune signature scores that reflected in part the intratumoral abundance of cytotoxic T cells showed significant positive associations with OS. However, the prognostic power of the T-cell signatures was highly dependent on TMB-high status, consistent with protective effects of tumor-infiltrating T cells in tumors with elevated antigenicity. In TMB-low tumors, a signature of infiltrating plasma B cells was significantly and positively associated with OS, independent of T-cell signature status. Although tumor subtypes based on differential expression patterns of CTA genes showed different survival associations within leiomyosarcoma and myxofibrosarcoma histologies, neither CTA nor histologic subtype interacted with the T-cell-survival association. CONCLUSIONS: Signatures of T-cell and plasma B-cell infiltrates were associated with a survival benefit in soft tissue sarcomas. TMB, but not CTA expression, influenced the prognostic power of T-cell-associated, but not plasma B-cell-associated, survival. LAY SUMMARY: Clinical data and RNA analysis of 259 soft tissue sarcomas from The Cancer Genome Atlas project were used to investigate associations between five published gene immune cell expression signatures and survival in the context of tumor mutations. Activated T cells had a significant positive association with patient survival. Although high tumor mutation burden was associated with good survival, the prognostic power of T-cell signatures was highly dependent on tumor mutational status, consistent with protective effects of tumor-infiltrating T cells in tumors with high levels of antigens. In low tumor mutation-bearing tumors, plasma B cells were positively associated with survival.

Skeletal muscle extracellular matrix remodeling with worsening glycemic control in nonhuman primates.

Type 2 diabetes (T2D) development may be mediated by skeletal muscle (SkM) function, which is responsible for >80% of circulating glucose uptake. The goals of this study were to assess changes in global- and location-level gene expression, remodeling proteins, fibrosis, and vascularity of SkM with worsening glycemic control, through RNA sequencing, immunoblotting, and immunostaining. We evaluated SkM samples from health-diverse African green monkeys (Cholorcebus aethiops sabaeus) to investigate these relationships. We assessed SkM remodeling at the molecular level by evaluating unbiased transcriptomics in age-, sex-, weight-, and waist circumference-matched metabolically healthy, prediabetic (PreT2D) and T2D monkeys (n = 13). Our analysis applied novel location-specific gene differences and shows that extracellular facing and cell membrane-associated genes and proteins are highly upregulated in metabolic disease. We verified transcript patterns using immunohistochemical staining and protein analyses of matrix metalloproteinase 16 (MMP16), tissue inhibitor of metalloproteinase 2 (TIMP2), and VEGF. Extracellular matrix (ECM) functions to support intercellular communications, including the coupling of capillaries to muscle cells, which was worsened with increasing blood glucose. Multiple regression modeling from age- and health-diverse monkeys (n = 33) revealed that capillary density was negatively predicted by only fasting blood glucose. The loss of vascularity in SkM co-occurred with reduced expression of hypoxia-sensing genes, which is indicative of a disconnect between altered ECM and reduced endothelial cells, and known perfusion deficiencies present in PreT2D and T2D. This report supports that rising blood glucose values incite ECM remodeling and reduce SkM capillarization, and that targeting ECM would be a rational approach to improve health with metabolic disease.

Prognostic Molecular Classification of Appendiceal Mucinous Neoplasms Treated with Cytoreductive Surgery and Hyperthermic Intraperitoneal Chemotherapy.

BACKGROUND: Appendiceal mucinous neoplasm (AMN) with peritoneal metastasis is a rare but deadly disease with few prognostic or therapy-predictive biomarkers to guide treatment decisions. Here, we investigated the prognostic and biological attributes of gene expression-based AMN molecular subtypes. METHODS: AMN specimens (n = 138) derived from a population-based subseries of patients treated at our institution with cytoreductive surgery and hyperthermic intraperitoneal chemotherapy (CRS/HIPEC) between 05/2000 and 05/2013 were analyzed for gene expression using a custom-designed NanoString 148-gene panel. Signed non-negative matrix factorization (sNMF) was used to define a gene signature capable of delineating robustly-classified AMN molecular subtypes. The sNMF class assignments were evaluated by topology learning, reverse-graph embedding and cross-cohort performance analysis. RESULTS: Three molecular subtypes of AMN were discerned by the expression patterns of 17 genes with roles in cancer progression or anti-tumor immunity. Tumor subtype assignments were confirmed by topology learning. AMN subtypes were termed immune-enriched (IE), oncogene-enriched (OE) and mixed (M) as evidenced by their gene expression patterns, and exhibited significantly different post-treatment survival outcomes. Genes with specialized immune functions, including markers of T-cells, natural killer cells, B-cells, and cytolytic activity showed increased expression in the low-risk IE subtype, while genes implicated in the promotion of cancer growth and progression were more highly expressed in the high-risk OE subtype. In multivariate analysis, the subtypes demonstrated independent prediction power for post-treatment survival. CONCLUSIONS: Our findings suggest a greater role for the immune system in AMN than previously recognized. AMN subtypes may have clinical utility for predicting CRS/HIPEC treatment outcomes.

APOL1 Renal-Risk Variants Induce Mitochondrial Dysfunction.

APOL1 G1 and G2 variants facilitate kidney disease in blacks. To elucidate the pathways whereby these variants contribute to disease pathogenesis, we established HEK293 cell lines stably expressing doxycycline-inducible (Tet-on) reference APOL1 G0 or the G1 and G2 renal-risk variants, and used Illumina human HT-12 v4 arrays and Affymetrix HTA 2.0 arrays to generate global gene expression data with doxycycline induction. Significantly altered pathways identified through bioinformatics analyses involved mitochondrial function; results from immunoblotting, immunofluorescence, and functional assays validated these findings. Overexpression of APOL1 by doxycycline induction in HEK293 Tet-on G1 and G2 cells led to impaired mitochondrial function, with markedly reduced maximum respiration rate, reserve respiration capacity, and mitochondrial membrane potential. Impaired mitochondrial function occurred before intracellular potassium depletion or reduced cell viability occurred. Analysis of global gene expression profiles in nondiseased primary proximal tubule cells from black patients revealed that the nicotinate phosphoribosyltransferase gene, responsible for NAD biosynthesis, was among the top downregulated transcripts in cells with two APOL1 renal-risk variants compared with those without renal-risk variants; nicotinate phosphoribosyltransferase also displayed gene expression patterns linked to mitochondrial dysfunction in HEK293 Tet-on APOL1 cell pathway analyses. These results suggest a pivotal role for mitochondrial dysfunction in APOL1-associated kidney disease.

SOX7 Target Genes and Their Contribution to Its Tumor Suppressive Function.

SOX7 is a transcription factor and acts as a tumor suppressor, but its target genes in cancers are poorly explored. We revealed SOX7-mediated gene expression profile in breast cancer cells using microarray chips and discovered multiple altered signaling pathways. When combinatorially analyzing the microarray data with a gene array dataset from 759 breast cancer patients, we identified four genes as potential targets of SOX7 and validated them by quantitative PCR and chromatin immunoprecipitation assays. Among these four genes, we determined that SOX7-activated SPRY1 and SLIT2, and SOX7-repressed TRIB3 and MTHFD2 could all differentially contribute to SOX7-mediated tumor suppression. Overall, we identified multiple cancer-related pathways mediated by SOX7 and for the first time revealed SOX7-regulated target genes in a cancer-relevant context.

Mapping adipose and muscle tissue expression quantitative trait loci in African Americans to identify genes for type 2 diabetes and obesity.

Relative to European Americans, type 2 diabetes (T2D) is more prevalent in African Americans (AAs). Genetic variation may modulate transcript abundance in insulin-responsive tissues and contribute to risk; yet, published studies identifying expression quantitative trait loci (eQTLs) in African ancestry populations are restricted to blood cells. This study aims to develop a map of genetically regulated transcripts expressed in tissues important for glucose homeostasis in AAs, critical for identifying the genetic etiology of T2D and related traits. Quantitative measures of adipose and muscle gene expression, and genotypic data were integrated in 260 non-diabetic AAs to identify expression regulatory variants. Their roles in genetic susceptibility to T2D, and related metabolic phenotypes, were evaluated by mining GWAS datasets. eQTL analysis identified 1971 and 2078 cis-eGenes in adipose and muscle, respectively. Cis-eQTLs for 885 transcripts including top cis-eGenes CHURC1, USMG5, and ERAP2 were identified in both tissues. 62.1 % of top cis-eSNPs were within ±50 kb of transcription start sites and cis-eGenes were enriched for mitochondrial transcripts. Mining GWAS databases revealed association of cis-eSNPs for more than 50 genes with T2D (e.g. PIK3C2A, RBMS1, UFSP1), gluco-metabolic phenotypes (e.g. INPP5E, SNX17, ERAP2, FN3KRP), and obesity (e.g. POMC, CPEB4). Integration of GWAS meta-analysis data from AA cohorts revealed the most significant association for cis-eSNPs of ATP5SL and MCCC1 genes, with T2D and BMI, respectively. This study developed the first comprehensive map of adipose and muscle tissue eQTLs in AAs (publically accessible at https://mdsetaa.phs.wakehealth.edu ) and identified genetically regulated transcripts for delineating genetic causes of T2D, and related metabolic phenotypes.

Conservation of immune gene signatures in solid tumors and prognostic implications.

BACKGROUND: Tumor-infiltrating leukocytes can either limit cancer growth or facilitate its spread. Diagnostic strategies that comprehensively assess the functional complexity of tumor immune infiltrates could have wide-reaching clinical value. In previous work we identified distinct immune gene signatures in breast tumors that reflect the relative abundance of infiltrating immune cells and exhibited significant associations with patient outcomes. Here we hypothesized that immune gene signatures agnostic to tumor type can be identified by de novo discovery of gene clusters enriched for immunological functions and possessing internal correlation structure conserved across solid tumors from different anatomic sites. METHODS: We assembled microarray expression datasets encompassing 5,295 tumors of the breast, colon, lung, ovarian and prostate. Unsupervised clustering methods were used to determine number and composition of gene clusters within each dataset. Immune-enriched gene clusters (signatures) identified by gene ontology enrichment were analyzed for internal correlation structure and conservation across tumors then compared against expression profiles of: 1) flow-sorted leukocytes from peripheral blood and 2) >300 cancer cell lines from solid and hematologic cancers. Cox regression analysis was used to identify signatures with significant associations with clinical outcome. RESULTS: We identified nine distinct immune-enriched gene signatures conserved across all five tumor types. The signatures differentiated specific leukocyte lineages with moderate discernment overall, and naturally organized into six discrete groups indicative of admixed lineages. Moreover, seven of the signatures exhibit minimal and uncorrelated expression in cancer cell lines, suggesting that these signatures derive predominantly from infiltrating immune cells. All nine immune signatures achieved statistically significant associations with patient prognosis (p<0.05) in one or more tumor types with greatest significance observed in breast and skin cancers. Several signatures indicative of myeloid lineages exhibited poor outcome associations that were most apparent in brain and colon cancers. CONCLUSIONS: These findings suggest that tumor infiltrating immune cells can be differentiated by immune-specific gene expression patterns that quantify the relative abundance of multiple immune infiltrates across a range of solid tumor types. That these markers of immune involvement are significantly associated with patient prognosis in diverse cancers suggests their clinical utility as pan-cancer markers of tumor behavior and immune responsiveness.

Re-Sequencing of the APOL1-APOL4 and MYH9 Gene Regions in African Americans Does Not Identify Additional Risks for CKD Progression.

BACKGROUND: In African Americans (AAs), APOL1 G1 and G2 nephropathy risk variants are associated with non-diabetic end-stage kidney disease (ESKD) in an autosomal recessive pattern. Additional risk and protective genetic variants may be present near the APOL1 loci, since earlier age ESKD is observed in some AAs with one APOL1 renal-risk variant, and because the adjacent gene MYH9 is associated with nephropathy in populations lacking G1 and G2 variants. METHODS: Re-sequencing was performed across a ∼275 kb region encompassing the APOL1-APOL4 and MYH9 genes in 154 AA cases with non-diabetic ESKD and 38 controls without nephropathy who were heterozygous for a single APOL1 G1 or G2 risk variant. RESULTS: Sequencing identified 3,246 non-coding single nucleotide polymorphisms (SNPs), 55 coding SNPs, and 246 insertion/deletions. No new coding variations were identified. Eleven variants, including a rare APOL3 Gln58Ter null variant (rs11089781), were genotyped in a replication panel of 1,571 AA ESKD cases and 1,334 controls. After adjusting for APOL1 G1 and G2 risk effects, these variations were not significantly associated with ESKD. In subjects with <2 APOL1 G1 and/or G2 alleles (849 cases; 1,139 controls), the APOL3 null variant was nominally associated with ESKD (recessive model, OR 1.81; p = 0.026); however, analysis in 807 AA cases and 634 controls from the Family Investigation of Nephropathy and Diabetes did not replicate this association. CONCLUSION: Additional common variants in the APOL1-APOL4-MYH9 region do not contribute significantly to ESKD risk beyond the APOL1 G1 and G2 alleles.

Complexes of vesicular stomatitis virus matrix protein with host Rae1 and Nup98 involved in inhibition of host transcription.

Vesicular stomatitis virus (VSV) suppresses antiviral responses in infected cells by inhibiting host gene expression at multiple levels, including transcription, nuclear cytoplasmic transport, and translation. The inhibition of host gene expression is due to the activity of the viral matrix (M) protein. Previous studies have shown that M protein interacts with host proteins Rae1 and Nup98 that have been implicated in regulating nuclear-cytoplasmic transport. However, Rae1 function is not essential for host mRNA transport, raising the question of how interaction of a viral protein with a host protein that is not essential for gene expression causes a global inhibition at multiple levels. We tested the hypothesis that there may be multiple M protein-Rae1 complexes involved in inhibiting host gene expression at multiple levels. Using size exclusion chromatography and sedimentation velocity analysis, it was determined that Rae1 exists in high, intermediate, and low molecular weight complexes. The intermediate molecular weight complexes containing Nup98 interacted most efficiently with M protein. The low molecular weight form also interacted with M protein in cells that overexpress Rae1 or cells in which Nup98 expression was silenced. Silencing Rae1 expression had little if any effect on nuclear accumulation of host mRNA in VSV-infected cells, nor did it affect VSV's ability to inhibit host translation. Instead, silencing Rae1 expression reduced the ability of VSV to inhibit host transcription. M protein interacted efficiently with Rae1-Nup98 complexes associated with the chromatin fraction of host nuclei, consistent with an effect on host transcription. These results support the idea that M protein-Rae1 complexes serve as platforms to promote the interaction of M protein with other factors involved in host transcription. They also support the idea that Rae1-Nup98 complexes play a previously under-appreciated role in regulation of transcription.

Comparative toxicogenomic responses of mercuric and methyl-mercury.

BACKGROUND: Mercury is a ubiquitous environmental toxicant that exists in multiple chemical forms. A paucity of information exists regarding the differences or similarities by which different mercurials act at the molecular level. RESULTS: Transcriptomes of mixed-stage C. elegans following equitoxic sub-, low- and high-toxicity exposures to inorganic mercuric chloride (HgCl2) and organic methylmercury chloride (MeHgCl) were analyzed. In C. elegans, the mercurials had highly different effects on transcription, with MeHgCl affecting the expression of significantly more genes than HgCl2. Bioinformatics analysis indicated that inorganic and organic mercurials affected different biological processes. RNAi identified 18 genes that were important in C. elegans response to mercurial exposure, although only two of these genes responded to both mercurials. To determine if the responses observed in C. elegans were evolutionarily conserved, the two mercurials were investigated in human neuroblastoma (SK-N-SH), hepatocellular carcinoma (HepG2) and embryonic kidney (HEK293) cells. The human homologs of the affected C. elegans genes were then used to test the effects on gene expression and cell viability after using siRNA during HgCl2 and MeHgCl exposure. As was observed with C. elegans, exposure to the HgCl2 and MeHgCl had different effects on gene expression, and different genes were important in the cellular response to the two mercurials. CONCLUSIONS: These results suggest that, contrary to previous reports, inorganic and organic mercurials have different mechanisms of toxicity. The two mercurials induced disparate effects on gene expression, and different genes were important in protecting the organism from mercurial toxicity.

Phase II trial of cytarabine and mitoxantrone with devimistat in acute myeloid leukemia.

Devimistat is a TCA cycle inhibitor. A previously completed phase I study of devimistat in combination with cytarabine and mitoxantrone in patients with relapsed or refractory AML showed promising response rates. Here we report the results of a single arm phase II study (NCT02484391). The primary outcome of feasibility of maintenance devimistat following induction and consolidation with devimistat in combination with high dose cytarabine and mitoxantrone was not met, as maintenance devimistat was only administered in 2 of 21 responders. The secondary outcomes of response (CR + CRi) and median survival were 44% (21/48) and 5.9 months respectively. There were no unexpected toxicities observed. An unplanned, post-hoc analysis of the phase I and II datasets suggests a trend of a dose response in older but not younger patients. RNA sequencing data from patient samples reveals an age-related decline in mitochondrial gene sets. Devimistat impairs ATP synthesis and we find a correlation between mitochondrial membrane potential and sensitivity to chemotherapy. Devimistat also induces mitochondrial reactive oxygen species and turnover consistent with mitophagy. We find that pharmacological or genetic inhibition of mitochondrial fission or autophagy sensitizes cells to devimistat. These findings suggest that an age related decline in mitochondrial quality and autophagy may be associated with response to devimistat however this needs to be confirmed in larger cohorts with proper trial design.

Bulk and Single-Cell Profiling of Breast Tumors Identifies TREM-1 as a Dominant Immune Suppressive Marker Associated With Poor Outcomes.

BACKGROUND: Triggering receptor expressed on myeloid cells (TREM)-1 is a key mediator of innate immunity previously associated with the severity of inflammatory disorders, and more recently, the inferior survival of lung and liver cancer patients. Here, we investigated the prognostic impact and immunological correlates of TREM1 expression in breast tumors. METHODS: Breast tumor microarray and RNAseq expression profiles (n=4,364 tumors) were analyzed for associations between gene expression, tumor immune subtypes, distant metastasis-free survival (DMFS) and clinical response to neoadjuvant chemotherapy (NAC). Single-cell (sc)RNAseq was performed using the 10X Genomics platform. Statistical associations were assessed by logistic regression, Cox regression, Kaplan-Meier analysis, Spearman correlation, Student's t-test and Chi-square test. RESULTS: In pre-treatment biopsies, TREM1 and known TREM-1 inducible cytokines (IL1B, IL8) were discovered by a statistical ranking procedure as top genes for which high expression was associated with reduced response to NAC, but only in the context of immunologically "hot" tumors otherwise associated with a high NAC response rate. In surgical specimens, TREM1 expression varied among tumor molecular subtypes, with highest expression in the more aggressive subtypes (Basal-like, HER2-E). High TREM1 significantly and reproducibly associated with inferior distant metastasis-free survival (DMFS), independent of conventional prognostic markers. Notably, the association between high TREM1 and inferior DMFS was most prominent in the subset of immunogenic tumors that exhibited the immunologically hot phenotype and otherwise associated with superior DMFS. Further observations from bulk and single-cell RNAseq analyses indicated that TREM1 expression was significantly enriched in polymorphonuclear myeloid-derived suppressor cells (PMN-MDSCs) and M2-like macrophages, and correlated with downstream transcriptional targets of TREM-1 (IL8, IL-1B, IL6, MCP-1, SPP1, IL1RN, INHBA) which have been previously associated with pro-tumorigenic and immunosuppressive functions. CONCLUSIONS: Together, these findings indicate that increased TREM1 expression is prognostic of inferior breast cancer outcomes and may contribute to myeloid-mediated breast cancer progression and immune suppression.

Dysregulated Pyrimidine Biosynthesis Contributes to 5-FU Resistance in SCLC Patient-Derived Organoids but Response to a Novel Polymeric Fluoropyrimidine, CF10.

Chemo-immunotherapy is central to the treatment of small cell lung cancer (SCLC). Despite modest progress made with the addition of immunotherapy, current cytotoxic regimens display minimal survival benefit and new treatments are needed. Thymidylate synthase (TS) is a well-validated anti-cancer drug target, but conventional TS inhibitors display limited clinical efficacy in refractory or recurrent SCLC. We performed RNA-Seq analysis to identify gene expression changes in SCLC biopsy samples to provide mechanistic insight into the potential utility of targeting pyrimidine biosynthesis to treat SCLC. We identified systematic dysregulation of pyrimidine biosynthesis, including elevated TYMS expression that likely contributes to the lack of efficacy for current TS inhibitors in SCLC. We also identified E2F1-3 upregulation in SCLC as a potential driver of TYMS expression that may contribute to tumor aggressiveness. To test if TS inhibition could be a viable strategy for SCLC treatment, we developed patient-derived organoids (PDOs) from human SCLC biopsy samples and used these to evaluate both conventional fluoropyrimidine drugs (e.g., 5-fluorouracil), platinum-based drugs, and CF10, a novel fluoropyrimidine polymer with enhanced TS inhibition activity. PDOs were relatively resistant to 5-FU and while moderately sensitive to the front-line agent cisplatin, were relatively more sensitive to CF10. Our studies demonstrate dysregulated pyrimidine biosynthesis contributes to drug resistance in SCLC and indicate that a novel approach to target these pathways may improve outcomes.

Transcriptomic Features of T Cell-Barren Tumors Are Conserved Across Diverse Tumor Types.

Background: Understanding how tumors subvert immune destruction is essential to the development of cancer immunotherapies. New evidence suggests that tumors limit anti-tumor immunity by exploiting transcriptional programs that regulate intratumoral trafficking and accumulation of effector cells. Here, we investigated the gene expression profiles that distinguish immunologically "cold" and "hot" tumors across diverse tumor types. Methods: RNAseq profiles of tumors (n = 8,920) representing 23 solid tumor types were analyzed using immune gene signatures that quantify CD8+ T cell abundance. Genes and pathways associated with a low CD8+ T cell infiltration profile (CD8-Low) were identified by correlation, differential expression, and statistical ranking methods. Gene subsets were evaluated in immunotherapy treatment cohorts and functionally characterized in cell lines and mouse tumor models. Results: Among different cancer types, we observed highly significant overlap of genes enriched in CD8-Low tumors, which included known immunomodulatory genes (e.g., BMP7, CMTM4, KDM5B, RCOR2) and exhibited significant associations with Wnt signaling, neurogenesis, cell-cell junctions, lipid biosynthesis, epidermal development, and cancer-testis antigens. Analysis of mutually exclusive gene clusters demonstrated that different transcriptional programs may converge on the T cell-cold phenotype as well as predict for response and survival of patients to Nivo treatment. Furthermore, we confirmed that a top-ranking candidate belonging to the TGF-ß superfamily, BMP7, negatively regulates CD8+ T cell abundance in immunocompetent murine tumor models, with and without anti-PD-L1 treatment. Conclusions: This study presents the first evidence that solid tumors of diverse anatomical origin acquire conserved transcriptional alterations that may be operative in the T cell-cold state. Our findings demonstrate the potential clinical utility of CD8-Low tumor-associated genes for predicting patient immunotherapy outcomes and point to novel mechanisms with potential for broad therapeutic exploitation.

MAP3K7 and CHD1 Are Novel Mediators of Resistance to Oncolytic Vesicular Stomatitis Virus in Prostate Cancer Cells.

A key principle of oncolytic viral therapy is that many cancers develop defects in their antiviral responses, making them more susceptible to virus infection. However, some cancers display resistance to viral infection. Many of these resistant cancers constitutively express interferon-stimulated genes (ISGs). The goal of these experiments was to determine the role of two tumor suppressor genes, MAP3K7 and CHD1, in viral resistance and ISG expression in PC3 prostate cancer cells resistant to oncolytic vesicular stomatitis virus (VSV). MAP3K7 and CHD1 are often co-deleted in aggressive prostate cancers. Silencing expression of MAP3K7 and CHD1 in PC3 cells increased susceptibility to the matrix (M) gene mutant M51R-VSV, as shown by increased expression of viral genes, increased yield of progeny virus, and reduction of tumor growth in nude mice. Silencing MAP3K7 alone had a greater effect on virus susceptibility than did silencing CHD1. Silencing MAP3K7 and CHD1 decreased constitutive expression of ISG mRNAs and proteins, whereas silencing MAP3K7 alone decreased expression of ISG proteins, but actually increased expression of ISG mRNAs. These results suggest a role for the protein product of MAP3K7, transforming growth factor ß-activated kinase 1 (TAK1), in regulating translation of ISG mRNAs and a role of CHD1 in maintaining the transcription of ISGs.

Solute Carrier Family 37 Member 2 (SLC37A2) Negatively Regulates Murine Macrophage Inflammation by Controlling Glycolysis.

Increased flux of glucose through glycolysis is a hallmark of inflammatory macrophages and is essential for optimal effector functions. Solute carrier (SLC) 37A2 is an endoplasmic reticulum-anchored phosphate-linked glucose-6-phosphate transporter that is highly expressed in macrophages and neutrophils. We demonstrate that SLC37A2 plays a pivotal role in murine macrophage inflammatory activation and cellular metabolic rewiring. Toll-like receptor (TLR) 4 stimulation by lipopolysaccharide (LPS) rapidly increases macrophage SLC37A2 protein expression. SLC37A2 deletion reprograms macrophages to a hyper-glycolytic process and accelerates LPS-induced inflammatory cytokine production, which partially depends on nicotinamide adenine dinucleotide (NAD+) biosynthesis. Blockade of glycolysis normalizes the differential expression of pro-inflammatory cytokines between control and SLC37A2 deficient macrophages. Conversely, overexpression of SLC37A2 lowers macrophage glycolysis and significantly reduces LPS-induced pro-inflammatory cytokine expression. In conclusion, our study suggests that SLC37A2 dampens murine macrophage inflammation by down-regulating glycolytic reprogramming as a part of macrophage negative feedback system to curtail acute innate activation.

APOL1 Kidney-Risk Variants Induce Mitochondrial Fission.

INTRODUCTION: APOL1 G1 and G2 nephropathy-risk variants cause mitochondrial dysfunction and contribute to kidney disease. Analyses were performed to determine the genetic regulation of APOL1 and elucidate potential mechanisms in APOL1-nephropathy. METHODS: A global gene expression analysis was performed in human primary renal tubule cell lines derived from 50 African American individuals. Follow-up gene knock out, cell-based rescue, and microscopy experiments were performed. RESULTS: APOL1 genotypes did not alter APOL1 expression levels in the global gene expression analysis. Expression quantitative trait locus (eQTL) analysis in polyinosinic-polycytidylic acid (poly IC)-stimulated renal tubule cells revealed that single nucleotide polymorphism (SNP) rs513349 adjacent to BAK1 was a trans eQTL for APOL1 and a cis eQTL for BAK1; APOL1 and BAK1 were co-expressed in cells. BAK1 knockout in a human podocyte cell line resulted in diminished APOL1 protein, supporting a pivotal effect for BAK1 on APOL1 expression. Because BAK1 is involved in mitochondrial dynamics, mitochondrial morphology was examined in primary renal tubule cells and HEK293 Tet-on cells of various APOL1 genotypes. Mitochondria in APOL1 wild-type (G0G0) tubule cells maintained elongated morphology when stimulated by low-dose poly IC, whereas those with G1G1, G2G2, and G1G2 genotypes appeared to fragment. HEK293 Tet-on cells overexpressing APOL1 G0, G1, and G2 were created; G0 cells appeared to promote mitochondrial fusion, whereas G1 and G2 induced mitochondrial fission. The mitochondrial dynamic regulator Mdivi-1 significantly preserved cell viability and mitochondrial cristae structure and reversed mitochondrial fission induced by overexpression of G1 and G2. CONCLUSION: Results suggest the mitochondrial fusion/fission pathway may be a therapeutic target in APOL1-nephropathy.

Discernment of possible mechanisms of hepatotoxicity via biological processes over-represented by co-expressed genes.

BACKGROUND: Hepatotoxicity is a form of liver injury caused by exposure to stressors. Genomic-based approaches have been used to detect changes in transcription in response to hepatotoxicants. However, there are no straightforward ways of using co-expressed genes anchored to a phenotype or constrained by the experimental design for discerning mechanisms of a biological response. RESULTS: Through the analysis of a gene expression dataset containing 318 liver samples from rats exposed to hepatotoxicants and leveraging alanine aminotransferase (ALT), a serum enzyme indicative of liver injury as the phenotypic marker, we identified biological processes and molecular pathways that may be associated with mechanisms of hepatotoxicity. Our analysis used an approach called Coherent Co-expression Biclustering (cc-Biclustering) for clustering a subset of genes through a coherent (consistency) measure within each group of samples representing a subset of experimental conditions. Supervised biclustering identified 87 genes co-expressed and correlated with ALT in all the samples exposed to the chemicals. None of the over-represented pathways related to liver injury. However, biclusters with subsets of samples exposed to one of the 7 hepatotoxicants, but not to a non-toxic isomer, contained co-expressed genes that represented pathways related to a stress response. Unsupervised biclustering of the data resulted in 1) four to five times more genes within the bicluster containing all the samples exposed to the chemicals, 2) biclusters with co-expression of genes that discerned 1,4 dichlorobenzene (a non-toxic isomer at low and mid doses) from the other chemicals, pathways and biological processes that underlie liver injury and 3) a bicluster with genes up-regulated in an early response to toxic exposure. CONCLUSION: We obtained clusters of co-expressed genes that over-represented biological processes and molecular pathways related to hepatotoxicity in the rat. The mechanisms involved in the response of the liver to the exposure to 1,4-dichlorobenzene suggest non-genotoxicity whereas the exposure to the hepatotoxicants could be DNA damaging leading to overall genomic instability and activation of cell cycle check point signaling. In addition, key pathways and biological processes representative of an inflammatory response, energy production and apoptosis were impacted by the hepatotoxicant exposures that manifested liver injury in the rat.

Sources of variation in baseline gene expression levels from toxicogenomics study control animals across multiple laboratories.

BACKGROUND: The use of gene expression profiling in both clinical and laboratory settings would be enhanced by better characterization of variance due to individual, environmental, and technical factors. Meta-analysis of microarray data from untreated or vehicle-treated animals within the control arm of toxicogenomics studies could yield useful information on baseline fluctuations in gene expression, although control animal data has not been available on a scale and in a form best served for data-mining. RESULTS: A dataset of control animal microarray expression data was assembled by a working group of the Health and Environmental Sciences Institute's Technical Committee on the Application of Genomics in Mechanism Based Risk Assessment in order to provide a public resource for assessments of variability in baseline gene expression. Data from over 500 Affymetrix microarrays from control rat liver and kidney were collected from 16 different institutions. Thirty-five biological and technical factors were obtained for each animal, describing a wide range of study characteristics, and a subset were evaluated in detail for their contribution to total variability using multivariate statistical and graphical techniques. CONCLUSION: The study factors that emerged as key sources of variability included gender, organ section, strain, and fasting state. These and other study factors were identified as key descriptors that should be included in the minimal information about a toxicogenomics study needed for interpretation of results by an independent source. Genes that are the most and least variable, gender-selective, or altered by fasting were also identified and functionally categorized. Better characterization of gene expression variability in control animals will aid in the design of toxicogenomics studies and in the interpretation of their results.