A survey of cancer cell lines reveals highly structured and hierarchical relationships within and between DNA and mRNA that may be the result of selection

Transcriptional Analysis Reveals Evidence of Chronically Impeded ECM Turnover and Epithelium-to-Mesenchyme Transition in Scar Tissue Giving Rise to Marjolin's Ulcer

Marjolin's ulcer (MU) is an aggressive malignancy arising within chronic wounds. A major cause is unhealed burn injuries. This results in well-differentiated squamous cell carcinoma (SCC). This study aimed to elucidate transcriptional changes leading to malignancy by investigating differentially expressed genes in squamous cells present in a SCC compared with MU. MU tumor cells were isolated from histologically confirmed biopsy of SCC within an unhealed burn scar. Epithelial cells (ECs) adjacent to the tumor were co-isolated and a SCC cell line was commercially purchased. mRNA from all three samples was isolated and its expression was quantified using RNASeq. A threshold of log2fold change >2-fold in either direction was considered "differentially expressed." Gene expression analysis revealed distinct differences in gene expression in MU cells compared with EC (665 genes), EC and SCC (1673 genes). Enrichment analysis confirmed that pathways most affected included 1) elevation of genes associated with extracellular matrix organization/degradation, 2) activation of DNA damage, and 3) activation of cytokine signaling. Our analysis revealed two key insights about chronic wound microenvironment conducive to ulceration. First, in EC vs. MU comparison, downregulation of Collagen and Matrix metalloproteinase families suggests chronically impaired extracellular matrix turnover giving rise to a fibrotic microenvironment. Second, in SCC vs. MU comparison, dysregulation of cadherin-mediated cell-cell adhesions is suggestive of epithelial-to-mesenchymal transitions, similar to those during development. Acquisition of epithelial-to-mesenchymal transition may underlie the high metastatic rate in MU tumors. Taken together, this sheds light on mechanisms that underlie the divergent clinical features of these cutaneous cancers.

Potential of Aqueous Humor as a Surrogate Tumor Biopsy for Retinoblastoma

Importance

Retinoblastoma (Rb) is one of the first tumors to have a known genetic etiology. However, because biopsy of this tumor is contraindicated, it has not been possible to define the effects of secondary genetic changes on the disease course.

Objective

To investigate whether the aqueous humor (AH) of Rb eyes has sufficient tumor-derived DNA to perform genetic analysis of the tumor, including DNA copy number alterations.

Design, Setting, and Participants

This investigation was a case series study at a tertiary care hospital (Children's Hospital Los Angeles) with a large Rb treatment center. Cell-free DNA (cfDNA) was isolated from 6 AH samples from 3 children with Rb, including 2 after primary enucleation and 1 undergoing multiple intravitreous injections of melphalan for vitreous seeding. Samples were taken between December 2014 and September 2015.

Main Outcomes and Measures

Measurable levels of nucleic acids in the AH and identification of tumor-derived DNA copy number variation in the AH. The AH was analyzed for DNA, RNA, and micro-RNA using Qubit high-sensitivity kits. Cell-free DNA was isolated from the AH, and sequencing library protocols were optimized. Shallow whole-genome sequencing was performed on an Illumina platform, followed by genome-wide chromosomal copy number variation profiling to assess the presence of tumor DNA fractions in the AH cfDNA of the 3 patients. One child's cfDNA from the AH and tumor DNA were subjected to Sanger sequencing to isolate the RB1 mutation.

Results

Six AH samples were obtained from 3 Rb eyes in 3 children (2 male and 1 female; diagnosed at ages 7, 20, and 28 months). A corroborative pattern between the chromosomal copy number variation profiles of the AH cfDNA and tumor-derived DNA from the enucleated samples was identified. In addition, a nonsense RB1 mutation (Lys→STOP) from 1 child was also identified from the AH samples obtained during intravitreous injection of melphalan, which matched the tumor sample postsecondary enucleation. Sanger sequencing of the AH cfDNA and tumor DNA with polymerase chain reaction primers targeting RB1 gene c.1075A demonstrated this same RB1 mutation.

Conclusions and Relevance

In this study evaluating nucleic acids in the AH from Rb eyes undergoing salvage therapy with intravitreous injection of melphalan, the results suggest that the AH can serve as a surrogate tumor biopsy when Rb tumor tissue is not available. This novel method will allow for analyses of tumor-derived DNA in Rb eyes undergoing salvage therapy that have not been enucleated.

A comprehensive characterization of genome-wide copy number aberrations in colorectal cancer reveals novel oncogenes and patterns of alterations

To develop a comprehensive overview of copy number aberrations (CNAs) in stage-II/III colorectal cancer (CRC), we characterized 302 tumors from the PETACC-3 clinical trial. Microsatellite-stable (MSS) samples (n = 269) had 66 minimal common CNA regions, with frequent gains on 20 q (72.5%), 7 (41.8%), 8 q (33.1%) and 13 q (51.0%) and losses on 18 (58.6%), 4 q (26%) and 21 q (21.6%). MSS tumors have significantly more CNAs than microsatellite-instable (MSI) tumors: within the MSI tumors a novel deletion of the tumor suppressor WWOX at 16 q23.1 was identified (p<0.01). Focal aberrations identified by the GISTIC method confirmed amplifications of oncogenes including EGFR, ERBB2, CCND1, MET, and MYC, and deletions of tumor suppressors including TP53, APC, and SMAD4, and gene expression was highly concordant with copy number aberration for these genes. Novel amplicons included putative oncogenes such as WNK1 and HNF4A, which also showed high concordance between copy number and expression. Survival analysis associated a specific patient segment featured by chromosome 20 q gains to an improved overall survival, which might be due to higher expression of genes such as EEF1B2 and PTK6. The CNA clustering also grouped tumors characterized by a poor prognosis BRAF-mutant-like signature derived from mRNA data from this cohort. We further revealed non-random correlation between CNAs among unlinked loci, including positive correlation between 20 q gain and 8 q gain, and 20 q gain and chromosome 18 loss, consistent with co-selection of these CNAs. These results reinforce the non-random nature of somatic CNAs in stage-II/III CRC and highlight loci and genes that may play an important role in driving the development and outcome of this disease.

Predictive genes in adjacent normal tissue are preferentially altered by sCNV during tumorigenesis in liver cancer and may rate limiting

Background

In hepatocellular carcinoma (HCC) genes predictive of survival have been found in both adjacent normal (AN) and tumor (TU) tissues. The relationships between these two sets of predictive genes and the general process of tumorigenesis and disease progression remains unclear.

Methodology/Principal Findings

Here we have investigated HCC tumorigenesis by comparing gene expression, DNA copy number variation and survival using ∼250 AN and TU samples representing, respectively, the pre-cancer state, and the result of tumorigenesis. Genes that participate in tumorigenesis were defined using a gene-gene correlation meta-analysis procedure that compared AN versus TU tissues. Genes predictive of survival in AN (AN-survival genes) were found to be enriched in the differential gene-gene correlation gene set indicating that they directly participate in the process of tumorigenesis. Additionally the AN-survival genes were mostly not predictive after tumorigenesis in TU tissue and this transition was associated with and could largely be explained by the effect of somatic DNA copy number variation (sCNV) in cis and in trans. The data was consistent with the variance of AN-survival genes being rate-limiting steps in tumorigenesis and this was confirmed using a treatment that promotes HCC tumorigenesis that selectively altered AN-survival genes and genes differentially correlated between AN and TU.

Conclusions/Significance

This suggests that the process of tumor evolution involves rate-limiting steps related to the background from which the tumor evolved where these were frequently predictive of clinical outcome. Additionally treatments that alter the likelihood of tumorigenesis occurring may act by altering AN-survival genes, suggesting that the process can be manipulated. Further sCNV explains a substantial fraction of tumor specific expression and may therefore be a causal driver of tumor evolution in HCC and perhaps many solid tumor types.