Cancer is sometimes depicted as a reversion to single cell behavior in cells adapted to live in a multicellular assembly. If this is the case, one would expect that mutation in cancer disrupts functional mechanisms that suppress cell-level traits detrimental to multicellularity. Such mechanisms should have evolved with or after the emergence of multicellularity. This leads to two related, but distinct hypotheses: 1) Somatic mutations in cancer will occur in genes that are younger than the emergence of multicellularity (1000 million years [MY]); and 2) genes that are frequently mutated in cancer and whose mutations are functionally important for the emergence of the cancer phenotype evolved within the past 1000 million years, and thus would exhibit an age distribution that is skewed to younger genes. In order to investigate these hypotheses we estimated the evolutionary ages of all human genes and then studied the probability of mutation and their biological function in relation to their age and genomic location for both normal germline and cancer contexts.

We observed that under a model of uniform random mutation across the genome, controlled for gene size, genes less than 500 MY were more frequently mutated in both cases. Paradoxically, causal genes, defined in the COSMIC Cancer Gene Census, were depleted in this age group. When we used functional enrichment analysis to explain this unexpected result we discovered that COSMIC genes with recessive disease phenotypes were enriched for DNA repair and cell cycle control. The non-mutated genes in these pathways are orthologous to those underlying stress-induced mutation in bacteria, which results in the clustering of single nucleotide variations. COSMIC genes were less common in regions where the probability of observing mutational clusters is high, although they are approximately 2-fold more likely to harbor mutational clusters compared to other human genes. Our results suggest this ancient mutational response to stress that evolved among prokaryotes was co-opted to maintain diversity in the germline and immune system, while the original phenotype is restored in cancer. Reversion to a stress-induced mutational response is a hallmark of cancer that allows for effectively searching “protected” genome space where genes causally implicated in cancer are located and underlies the high adaptive potential and concomitant therapeutic resistance that is characteristic of cancer.

The histone deacetylase (HDAC) inhibitor vorinostat has received significant attention in recent years as an ‘epigenetic’ drug used to treat solid tumors. However, its mechanisms of action are not entirely understood, particularly with regard to its interaction with the aberrations in 3D nuclear structure that accompany neoplastic progression. We investigated the impact of vorinostat on human esophageal epithelial cell lines derived from normal, metaplastic (pre-cancerous), and malignant tissue. Using a combination of novel optical computed tomography (CT)-based quantitative 3D absorption microscopy and conventional confocal fluorescence microscopy, we show that subjecting malignant cells to vorinostat preferentially alters their 3D nuclear architecture relative to non-cancerous cells. Optical CT (cell CT) imaging of fixed single cells showed that drug-treated cancer cells exhibit significant alterations in nuclear morphometry. Confocal microscopy revealed that vorinostat caused changes in the distribution of H3K9ac-marked euchromatin and H3K9me3-marked constitutive heterochromatin. Additionally, 3D immuno-FISH showed that drug-induced expression of the DNA repair gene MGMT was accompanied by spatial relocation toward the center of the nucleus in the nuclei of metaplastic but not in non-neoplastic cells. Our data suggest that vorinostat’s differential modulation of 3D nuclear architecture in normal and abnormal cells could play a functional role in its anti-cancer action.