Connectivity mapping using a combined gene signature from multiple colorectal cancer datasets identified candidate drugs including existing chemotherapies

BACKGROUND: While the discovery of new drugs is a complex, lengthy and costly process, identifying new uses for existing drugs is a cost-effective approach to therapeutic discovery. Connectivity mapping integrates gene expression profiling with advanced algorithms to connect genes, diseases and small molecule compounds and has been applied in a large number of studies to identify potential drugs, particularly to facilitate drug repurposing. Colorectal cancer (CRC) is a commonly diagnosed cancer with high mortality rates, presenting a worldwide health problem. With the advancement of high throughput omics technologies, a number of large scale gene expression profiling studies have been conducted on CRCs, providing multiple datasets in gene expression data repositories. In this work, we systematically apply gene expression connectivity mapping to multiple CRC datasets to identify candidate therapeutics to this disease.

RESULTS: We developed a robust method to compile a combined gene signature for colorectal cancer across multiple datasets. Connectivity mapping analysis with this signature of 148 genes identified 10 candidate compounds, including irinotecan and etoposide, which are chemotherapy drugs currently used to treat CRCs. These results indicate that we have discovered high quality connections between the CRC disease state and the candidate compounds, and that the gene signature we created may be used as a potential therapeutic target in treating the disease. The method we proposed is highly effective in generating quality gene signature through multiple datasets; the publication of the combined CRC gene signature and the list of candidate compounds from this work will benefit both cancer and systems biology research communities for further development and investigations.

Myc-dependent purine biosynthesis affects nucleolar stress and therapy response in prostate cancer

The androgen receptor is a key transcription factor contributing to the development of all stages of prostate cancer (PCa). In addition, other transcription factors have been associated with poor prognosis in PCa, amongst which c-Myc (MYC) is a well-established oncogene in many other cancers. We have previously reported that the AR promotes glycolysis and anabolic metabolism; many of these metabolic pathways are also MYC-regulated in other cancers. In this study, we report that in PCa cells de novo purine biosynthesis and the subsequent conversion to XMP is tightly regulated by MYC and independent of AR activity. We characterized two enzymes, PAICS and IMPDH2, within the pathway as PCa biomarkers in tissue samples and report increased efficacy of established anti-androgens in combination with a clinically approved IMPDH inhibitor, mycophenolic acid (MPA). Treatment with MPA led to a significant reduction in cellular guanosine triphosphate (GTP) levels accompanied by nucleolar stress and p53 stabilization. In conclusion, targeting purine biosynthesis provides an opportunity to perturb PCa metabolism and enhance tumour suppressive stress responses.

Molecular pathways: Estrogen pathway in colorectal cancer.

Worldwide, colorectal cancer has a higher incidence rate in men than in women, suggesting a protective role for sex hormones in the development of the disease. Preclinical data support a role for estrogen and its receptors in the initiation and progression of colorectal cancer and establishes that protective effects of estrogen are exerted through ERβ. Hormone replacement therapy (HRT) in postmenopausal women as well as consumption of soy reduces the incidence of colorectal cancer. In the Women's Health Initiative trial, use of HRT in postmenopausal women reduced the risk of colon cancer by 56% [95% confidence interval (CI), 0.38-0.81; P = 0.003]. A recent meta-analysis showed that in women, consumption of soy reduced the risk of colon cancer by 21% (95% CI, 0.03-0.35; P = 0.026). In this review, using the preclinical data, we translate the findings in the clinical trials and observational studies to define the role of estrogen in the prevention of colorectal cancer. We hypothesize that sometime during the tumorigenesis process ERβ expression in colonocytes is lost and the estrogen ligand, HRT, or soy products, exerts its effects through preventing this loss. Thus, in the adenoma-to-carcinoma continuum, timing of HRT is a significant determinant of the observed benefit from this intervention. We further argue that the protective effects of estrogen are limited to certain molecular subtypes. Successful development of estrogen modulators for prevention of colorectal cancer depends on identification of susceptible colorectal cancer population(s). Thus, research to better understand the estrogen pathway is fundamental for clinical delivery of these agents.

Histone deacetylase inhibitors suppress thymidylate synthase gene expression and synergize with the fluoropyrimidines in colon cancer cells

Despite recent therapeutic advances, the response rates to chemotherapy for patients with metastatic colon cancer remain at approximately 50% with the fluoropyrimidine, 5-fluorouracil (5-FU), continuing to serve as the foundation chemotherapeutic agent for the treatment of this disease. Previous studies have demonstrated that overexpression of thymidylate synthase (TS) is a key determinant of resistance to 5-FU-based chemotherapy. Therefore, there is a significant need to develop alternative therapeutic strategies to overcome TS-mediated resistance. In this study, we demonstrate that the histone deacetylase inhibitors (HDACi) vorinostat and LBH589 significantly downregulate TS gene expression in a panel of colon cancer cell lines. Downregulation of TS was independent of p53, p21 and HDAC2 expression and was achievable in vivo as demonstrated by mouse xenograft models. We provide evidence that HDACi treatment leads to a potent transcriptional repression of the TS gene. Combination of the fluoropyrimidines 5-FU or FUdR with both vorinostat and LBH589 enhanced cell cycle arrest and growth inhibition. Importantly, the downstream effects of TS inhibition were significantly enhanced by this combination including the inhibition of acute TS induction and the enhanced accumulation of the cytotoxic nucleotide intermediate dUTP. These data demonstrate that HDACi repress TS expression at the level of transcription and provides the first evidence suggesting a direct mechanistic link between TS downregulation and the synergistic interaction observed between HDACi and 5-FU. This study provides rationale for the continued clinical evaluation of HDACi in combination with 5-FU-based therapies as a strategy to overcome TS-mediated resistance.

A coactivator role of CARM1 in the dysregulation of β-catenin activity in colorectal cancer cell growth and gene expression

Aberrant activation of Wnt/β-catenin signaling, resulting in the expression of Wnt-regulated oncogenes, is recognized as a critical factor in the etiology of colorectal cancer. Occupancy of β-catenin at promoters of Wnt target genes drives transcription, but the mechanism of β-catenin action remains poorly understood. Here, we show that CARM1 (coactivator-associated arginine methyltransferase 1) interacts with β-catenin and positively modulates β-catenin-mediated gene expression. In colorectal cancer cells with constitutively high Wnt/β-catenin activity, depletion of CARM1 inhibits expression of endogenous Wnt/β-catenin target genes and suppresses clonal survival and anchorage-independent growth. We also identified a colorectal cancer cell line (RKO) with a low basal level of β-catenin, which is dramatically elevated by treatment with Wnt3a. Wnt3a also increased the expression of a subset of endogenous Wnt target genes, and CARM1 was required for the Wnt-induced expression of these target genes and the accompanying dimethylation of arginine 17 of histone H3. Depletion of β-catenin from RKO cells diminished the Wnt-induced occupancy of CARM1 on a Wnt target gene, indicating that CARM1 is recruited to Wnt target genes through its interaction with β-catenin and contributes to transcriptional activation by mediating events (including histone H3 methylation) that are downstream from the actions of β-catenin. Therefore, CARM1 is an important positive modulator of Wnt/β-catenin transcription and neoplastic transformation, and may thereby represent a novel target for therapeutic intervention in cancers involving aberrantly activated Wnt/β-catenin signaling.

A let-7 microRNA-binding site polymorphism in 3'-untranslated region of KRAS gene predicts response in wild-type KRAS patients with metastatic colorectal cancer treated with cetuximab monotherapy

PURPOSE: recent studies have found that KRAS mutations predict resistance to monoclonal antibodies targeting the epidermal growth factor receptor in metastatic colorectal cancer (mCRC). A polymorphism in a let-7 microRNA complementary site (lcs6) in the KRAS 3' untranslated region (UTR) is associated with an increased cancer risk in non-small-cell lung cancer and reduced overall survival (OS) in oral cancers. We tested the hypothesis whether this polymorphism may be associated with clinical outcome in KRAS wild-type (KRASwt) mCRC patients treated with cetuximab monotherapy.

PATIENTS AND METHODS: the presence of KRAS let-7 lcs6 polymorphism was evaluated in 130 mCRC patients who were enrolled in a phase II study of cetuximab monotherapy (IMCL-0144). Genomic DNA was extracted from dissected formalin-fixed paraffin-embedded tumor tissue, KRAS mutation status and polymorphism were assessed using direct sequencing and PCR restriction fragment length polymorphism technique.

RESULTS: KRAS let-7 lcs6 polymorphism was found to be related to object response rate (ORR) in mCRC patients whose tumors had KRASwt. The 12 KRASwt patients harboring at least a variant G allele (TG or GG) had a 42% ORR compared with a 9% ORR in 55 KRASwt patients with let-7 lcs6 TT genotype (P = 0.02, Fisher's exact test). KRASwt patients with TG/GG genotypes had trend of longer median progression-free survival (3.9 versus 1.3 months) and OS (10.7 versus 6.4 months) compared to those with TT genotypes.

CONCLUSIONS: these results are the first to indicate that the KRAS 3'UTR polymorphism may predict for cetuximab responsiveness in KRASwt mCRC patients, which warrants validation in other clinical trials.

Common cancer stem cell gene variants predict colon cancer recurrence

PURPOSE: Recent evidence suggests that cancer stem cells (CSC) are responsible for key elements of colon cancer progression and recurrence. Germline variants in CSC genes may result in altered gene function and/or activity, thereby causing interindividual differences in a patient's tumor recurrence capacity and chemoresistance. We investigated germline polymorphisms in a comprehensive panel of CSC genes to predict time to tumor recurrence (TTR) in patients with stage III and high-risk stage II colon cancer.

EXPERIMENTAL DESIGN: A total of 234 patients treated with 5-fluorouracil-based chemotherapy at the University of Southern California were included in this study. Whole blood samples were analyzed for germline polymorphisms in genes that have been previously associated with colon CSC (CD44, Prominin-1, DPP4, EpCAM, ALCAM, Msi-1, ITGB1, CD24, LGR5, and ALDH1A1) by PCR-RFLP or direct DNA-sequencing.

RESULTS: The minor alleles of CD44 rs8193 C>T, ALCAM rs1157 G>A, and LGR5 rs17109924 T>C were significantly associated with increased TTR (9.4 vs. 5.4 years; HR, 0.51; 95% CI: 0.35-0.93; P = 0.022; 11.3 vs. 5.7 years; HR, 0.56; 95% CI: 0.33-0.94; P = 0.024, and 10.7 vs. 5.7 years; HR, 0.33; 95% CI: 0.12-0.90; P = 0.023, respectively) and remained significant in the multivariate analysis stratified by ethnicity. In recursive partitioning, a specific gene variant profile including LGR5 rs17109924, CD44 rs8193, and ALDH1A1 rs1342024 represented a high-risk subgroup with a median TTR of 1.7 years (HR, 6.71, 95% CI: 2.71-16.63, P < 0.001).

CONCLUSION: This is the first study identifying common germline variants in colon CSC genes as independent prognostic markers for stage III and high-risk stage II colon cancer patients.

Single nucleotide polymorphisms in the APCS gene influence geographic atrophy secondary to VEGF inhibition therapy in neovascular macular degeneration.

Introduction: Age-related macular degeneration (AMD) is a leading cause of vision loss in the elderly mostly due to the development of neovascular AMD (nAMD) or geographic atrophy (GA). Intravitreal injections of anti-vascular endothelial growth factor (VEGF) agents are an effective therapeutic option for nAMD. Following anti-VEGF treatments, increased atrophy of the retinal pigment epithelium (RPE) and choriocapillaries that resembles GA has been reported. We sought to evaluate the underlying genetic influences that may contribute to this process. Methods: We selected 68 single nucleotide polymorphisms (SNPs) from genes previously identified as susceptibility factors in AMD, along with 43 SNPs from genes encoding the VEGF protein and its cognate receptors as this pathway is targeted by treatment. We enrolled 467 consecutive patients (Feb 2009 to October 2011) with nAMD who received anti-VEGF therapy. The acutely presenting eye was designated as the study eye and retinal tomograms graded for macular atrophy at study exit. Statistical analysis was performed using PLINK to identify SNPs with a P value < 0.01. Logistic regression models with macular atrophy as dependent variable were fitted with age, gender, smoking status, common genetic risk factors and the identified SNPs as explanatory variables. Results: Grading for macular atrophy was available in 304 study eyes and 70% (214) were classified as showing macular atrophy. In the unadjusted analysis we observed significant associations between macular atrophy and two independent SNPs in the APCS gene: rs6695377: odds ratio (OR) = 1.98; 95% confidence intervals (CI): 1.23, 3.19; P = 0.004; rs1446965: OR = 2.49, CI: 1.29, 4.82; P = 0.006 and these associations remained significant after adjustment for covariates. Conclusions: VEGF is a mitogen and growth factor for choroidal blood vessels and the RPE and its inhibition could lead to atrophy of these key tissues. Anti-VEGF treatment can interfere with ocular vascular maintenance and may be associated with RPE and choroidal atrophy. As such, these medications, which block the effects of VEGF, may influence the development of GA. The top associated SNPs are found in the APCS gene, a highly conserved glycoprotein that encodes Serum amyloid P (SAP) which opsonizes apoptotic cells. SAP can bind to and activate complement components via binding to C1q, a mechanism by which SAP may remove cellular debris, affecting regulation of the three complement pathways.

Natural Killer Cell-like signature observed in locally advanced rectal cancer after neoadjuvant chemoradiotherapy in a tumour regression grade dependent manner

Background: Around 10-15% of patients with locally advanced rectal cancer (LARC) undergo a pathologically complete response (TRG4) to neoadjuvant chemoradiotherapy; the rest of patients exhibit a spectrum of tumour regression (TRG1-3). Understanding therapy-related genomic alterations may help us to identify underlying biology or novel targets associated with response that could increase the efficacy of therapy in patients that do not benefit from the current standard of care.
Methods: 48 FFPE rectal cancer biopsies and matched resections were analysed using the WG-DASL HumanHT-12_v4 Beadchip array on the illumina iScan. Bioinformatic analysis was conducted in Partek genomics suite and R studio. Limma and glmnet packages were used to identify genes differentially expressed between tumour regression grades. Validation of microarray results will be carried out using IHC, RNAscope and RT-PCR.
Results: Immune response genes were observed from supervised analysis of the biopsies which may have predictive value. Differential gene expression from the resections as well as pre and post therapy analysis revealed induction of genes in a tumour regression dependent manner. Pathway mapping and Gene Ontology analysis of these genes suggested antigen processing and natural killer mediated cytotoxicity respectively. The natural killer-like gene signature was switched off in non-responders and on in the responders. IHC has confirmed the presence of Natural killer cells through CD56+ staining.
Conclusion: Identification of NK cell genes and CD56+ cells in patients responding to neoadjuvant chemoradiotherapy warrants further investigation into their association with tumour regression grade in LARC. NK cells are known to lyse malignant cells and determining whether their presence is a cause or consequence of response is crucial. Interrogation of the cytokines upregulated in our NK-like signature will help guide future in vitro models.