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.

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.

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.

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.

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.

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.

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.

Tumor mutational burden is a determinant of immune-mediated survival in breast cancer.

Mounting evidence supports a role for the immune system in breast cancer outcomes. The ability to distinguish highly immunogenic tumors susceptible to anti-tumor immunity from weakly immunogenic or inherently immune-resistant tumors would guide development of therapeutic strategies in breast cancer. Genomic, transcriptomic and clinical data from The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast cancer cohorts were used to examine statistical associations between tumor mutational burden (TMB) and the survival of patients whose tumors were assigned to previously-described prognostic immune subclasses reflecting favorable, weak or poor immune-infiltrate dispositions (FID, WID or PID, respectively). Tumor immune subclasses were associated with survival in patients with high TMB (TMB-Hi, P < 0.001) but not in those with low TMB (TMB-Lo, P = 0.44). This statistical relationship was confirmed in the METABRIC cohort (TMB-Hi, P = 0.047; TMB-Lo, P = 0.39), and also found to hold true in the more-indolent Luminal A tumor subtype (TMB-Hi, P = 0.011; TMB-Lo, P = 0.91). In TMB-Hi tumors, the FID subclass was associated with prolonged survival independent of tumor stage, molecular subtype, age and treatment. Copy number analysis revealed the reproducible, preferential amplification of chromosome 1q immune-regulatory genes in the PID immune subclass. These findings demonstrate a previously unappreciated role for TMB as a determinant of immune-mediated survival of breast cancer patients and identify candidate immune-regulatory mechanisms associated with immunologically cold tumors. Immune subtyping of breast cancers may offer opportunities for therapeutic stratification.

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.

A Phase I Study of CPI-613 in Combination with High-Dose Cytarabine and Mitoxantrone for Relapsed or Refractory Acute Myeloid Leukemia.

Purpose: CPI-613, a lipoate analogue that inhibits pyruvate dehydrogenase (PDH) and α-ketogluterate dehydrogenase (KGDH), has activity in patients with myeloid malignancies. This study explored the role of mitochondrial metabolism in chemotherapy response and determined the MTD, efficacy, and safety of CPI-613 combined with high-dose cytarabine and mitoxantrone in patients with relapsed or refractory acute myeloid leukemia.Experimental Design: The role of mitochondrial response to chemotherapy was assessed in cell lines and animal models. A phase I study of CPI-613 plus cytarabine and mitoxantrone was conducted in patients with relapsed or refractory AML.Results: Exposure to chemotherapy induced mitochondrial oxygen consumption that depended on PDH. CPI-613 sensitized AML cells to chemotherapy indicating that mitochondrial metabolism is a source of resistance. Loss of p53 did not alter response to CPI-613. The phase I study enrolled 67 patients and 62 were evaluable for response. The overall response rate was 50% (26CR+5CRi/62). Median survival was 6.7 months. In patients over 60 years old, the CR/CRi rate was 47% (15/32) with a median survival of 6.9 months. The response rate for patients with poor-risk cytogenetics also was encouraging with 46% (11/24 patients) achieving a CR or CRi. RNA sequencing analysis of a subset of baseline bone marrow samples revealed a gene expression signature consistent with the presence of B cells in the pretreatment marrow of responders.Conclusions: The addition of CPI-613 to chemotherapy is a promising approach in older patients and those with poor-risk cytogenetics. Clin Cancer Res; 24(9); 2060-73. ©2018 AACR.

Beta-catenin represses protein kinase D1 gene expression by non-canonical pathway through MYC/MAX transcription complex in prostate cancer.

Down regulation of Protein Kinase D1 (PrKD1), a novel serine threonine kinase, in prostate, gastric, breast and colon cancers in humans leads to disease progression. While the down regulation of PrKD1 by DNA methylation in gastric cancer and by nuclear beta-catenin in colon cancer has been shown, the regulatory mechanisms in other cancers are unknown. Because we had demonstrated that PrKD1 is the only known kinase to phosphorylate threonine 120 (T120) of beta-catenin in prostate cancer resulting in increased nuclear beta-catenin, we explored the role of beta-catenin in gene regulation of PrKD1. An initial CHIP assay identified potential binding sites for beta-catenin in and downstream of PrKD1 promoter and sequencing confirmed recruitment of beta-catenin to a 166 base pairs sequence upstream of exon 2. Co-transfection studies with PrKD1-promoter-reporter suggested that beta-catenin represses PrKD1 promoter. Efforts to identify transcription factors that mediate the co-repressor effects of beta-catenin identified recruitment of both MYC and its obligate heterodimer MAX to the same binding site as beta-catenin on the PrKD1 promoter site. Moreover, treatment with MYC inhibitor rescued the co-repressor effect of beta-catenin on PrKD1 gene expression. Prostate specific knock out of PrKD1 in transgenic mice demonstrated increased nuclear expression of beta-catenin validating the in vitro studies. Functional studies showed that nuclear translocation of beta-catenin as a consequence of PrKD1 down regulation, increases AR transcriptional activity with attendant downstream effects on androgen responsive genes. In silico human gene expression analysis confirmed the down regulation of PrKD1 in metastatic prostate cancer correlated inversely with the expression of MAX, but not MYC, and positively with MXD1, a competing heterodimer of MAX, suggesting that the dimerization of MAX with either MYC or MXD1 regulates PrKD1 gene expression. The study has identified a novel auto-repressive loop that perpetuates PrKD1 down regulation through beta-catenin/MYC/MAX protein complex.

Genetic regulation of adipose tissue transcript expression is involved in modulating serum triglyceride and HDL-cholesterol.

Dyslipidemia is a major contributor to the increased cardiovascular disease and mortality associated with obesity and type 2 diabetes. We hypothesized that variation in expression of adipose tissue transcripts is associated with serum lipid concentrations in African Americans (AAs), and common genetic variants regulate expression levels of these transcripts. Fasting serum lipid levels, genome-wide transcript expression profiles of subcutaneous adipose tissue, and genome-wide SNP genotypes were analyzed in a cohort of non-diabetic AAs (N=250). Serum triglyceride (TRIG) and high density lipoprotein-cholesterol (HDL-C) levels were associated (FDR<0.01) with expression level of 1021 and 1875 adipose tissue transcripts, respectively, but none associated with total cholesterol or LDL-C levels. Serum HDL-C-associated transcripts were enriched for salient biological pathways, including branched-chain amino acid degradation, and oxidative phosphorylation. Genes in immuno-inflammatory pathways were activated among individuals with higher serum TRIG levels. We identified significant cis-regulatory SNPs (cis-eSNPs) for 449 serum lipid-associated transcripts in adipose tissue. The cis-eSNPs of 12 genes were nominally associated (p<0.001) with serum lipid level in genome wide association studies in Global Lipids Genetics Consortium (GLGC) cohorts. Allelic effect direction of cis-eSNPs on expression of MARCH2, BEST1 and TMEM258 matched with effect direction of these SNP alleles on serum TRIG or HDL-C levels in GLGC cohorts. These data suggest that expressions of serum lipid-associated transcripts in adipose tissue are dependent on common cis-eSNPs in African Americans. Thus, genetically-mediated transcriptional regulation in adipose tissue may play a role in reducing HDL-C and increasing TRIG in serum.

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.

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.

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.

Effects of Pubertal Exposure to Dietary Soy on Estrogen Receptor Activity in the Breast of Cynomolgus Macaques.

Endogenous estrogens influence mammary gland development during puberty and breast cancer risk during adulthood. Early-life exposure to dietary or environmental estrogens may alter estrogen-mediated processes. Soy foods contain phytoestrogenic isoflavones (IF), which have mixed estrogen agonist/antagonist properties. Here, we evaluated mammary gland responses over time in pubertal female cynomolgus macaques fed diets containing either casein/lactalbumin (n = 12) or soy protein containing a human-equivalent dose of 120 mg IF/day (n = 17) for approximately 4.5 years spanning menarche. We assessed estrogen receptor (ER) expression and activity, promoter methylation of ERs and their downstream targets, and markers of estrogen metabolism. Expression of ERα and classical ERα response genes (TFF1, PGR, and GREB1) decreased with maturity, independent of diet. A significant inverse correlation was observed between TFF1 mRNA and methylation of CpG sites within the TFF1 promoter. Soy effects included lower ERß expression before menarche and lower mRNA for ERα and GREB1 after menarche. Expression of GATA-3, an epithelial differentiation marker that regulates ERα-mediated transcription, was elevated before menarche and decreased after menarche in soy-fed animals. Soy did not significantly alter expression of other ER activity markers, estrogen-metabolizing enzymes, or promoter methylation for ERs or ER-regulated genes. Our results demonstrate greater ER expression and activity during the pubertal transition, supporting the idea that this life stage is a critical window for phenotypic modulation by estrogenic compounds. Pubertal soy exposure decreases mammary ERα expression after menarche and exerts subtle effects on receptor activity and mammary gland differentiation. Cancer Prev Res; 9(5); 385-95. ©2016 AACR.