Cancer progression as a sequence of atavistic reversions.

It has long been recognized that cancer onset and progression represent a type of reversion to an ancestral quasi-unicellular phenotype. This general concept has been refined into the atavistic model of cancer that attempts to provide a quantitative analysis and testable predictions based on genomic data. Over the past decade, support for the multicellular-to-unicellular reversion predicted by the atavism model has come from phylostratigraphy. Here, we propose that cancer onset and progression involve more than a one-off multicellular-to-unicellular reversion, and are better described as a series of reversionary transitions. We make new predictions based on the chronology of the unicellular-eukaryote-to-multicellular-eukaryote transition. We also make new predictions based on three other evolutionary transitions that occurred in our lineage: eukaryogenesis, oxidative phosphorylation and the transition to adaptive immunity. We propose several modifications to current phylostratigraphy to improve age resolution to test these predictions. Also see the video abstract here: https://youtu.be/3unEu5JYJrQ.

Whole-genome sequencing of an aggressive BRAF wild-type papillary thyroid cancer identified EML4-ALK translocation as a therapeutic target.

BACKGROUND: Recent advances in the treatment of cancer have focused on targeting genomic aberrations with selective therapeutic agents. In radioiodine resistant aggressive papillary thyroid cancers, there remain few effective therapeutic options. A 62-year-old man who underwent multiple operations for papillary thyroid cancer and whose metastases progressed despite standard treatments provided tumor tissue. METHODS: We analyzed tumor and whole blood DNA by whole genome sequencing, achieving 80× or greater coverage over 94 % of the exome and 90 % of the genome. We determined somatic mutations and structural alterations. RESULTS: We found a total of 57 somatic mutations in 55 genes of the cancer genome. There was notably a lack of mutations in NRAS and BRAF, and no RET/PTC rearrangement. There was a mutation in the TRAPP oncogene and a loss of heterozygosity of the p16, p18, and RB1 tumor suppressor genes. The oncogenic driver for this tumor is a translocation involving the genes for anaplastic lymphoma receptor tyrosine kinase (ALK) and echinoderm microtubule associated protein like 4 (EML4). The EML4-ALK translocation has been reported in approximately 5 % of lung cancers, as well as in pediatric neuroblastoma, and is a therapeutic target for crizotinib. CONCLUSIONS: This is the first report of the whole genomic sequencing of a papillary thyroid cancer in which we identified an EML4-ALK translocation of a TRAPP oncogene mutation. These findings suggest that this tumor has a more distinct oncogenesis than BRAF mutant papillary thyroid cancer. Whole genome sequencing can elucidate an oncogenic context and expose potential therapeutic vulnerabilities in rare cancers.

New Validated Tool to Diagnose Breastfeeding Dysfunction.

BACKGROUND: Breastfeeding behaviors and experiences exist on a continuum. What differentiates normal from dysfunctional is defined by frequency and severity. No current validated tool addresses the subjective experience of dyads with a predictive score that can be followed over time. RESEARCH AIM: To create and validate a self-report tool to assess breastfeeding and evaluate its ability to predict risk of breastfeeding dysfunction. METHODS: This study used a cross-sectional design to determine the validity of a novel instrument to assess breastfeeding dysfunction. We gave the initial questionnaire to 2085 breastfeeding dyads. We assessed content validity by comparison with other tools. We used exploratory factor analysis with varimax rotation for concept identification and Cronbach's alpha for internal consistency. We employed logistic regression to assess the tool's ability to differentiate between normal breastfeeding and breastfeeding dysfunction. RESULTS: Factor analysis mapped 17 questions to four concepts to create a score (FLIP; flow, latch, injury [to the nipple], and post-feed behavior). Internal consistency and reliability of the scores in these concepts were acceptable (Cronbach's alpha ≥ 0.087 for all measures). A logistic regression model that controlled for infant age, with a breastfeeding dysfunction risk classification threshold of 60%, yielded a correct classification of 88.7%, with 93.1% sensitivity, 64.6% specificity, and a 6.5% false positive rate. CONCLUSIONS: The FLIP score was determined to be a valid and reliable instrument for quantifying the severity of breastfeeding dysfunction in children under 1 year old. Further studies will assess its usefulness in the management of breastfeeding dysfunction.

Preclinical investigation of nanoparticle albumin-bound paclitaxel as a potential treatment for adrenocortical cancer.

BACKGROUND: Traditional drug discovery methods have a limited role in rare cancers. We hypothesized that molecular technology including gene expression profiling could expose novel targets for therapy in adrenocortical carcinoma (ACC), a rare and lethal cancer. SPARC (secreted protein acidic rich in cysteine) is an albumin-binding matrix-associated protein that is proposed to act as a mechanism for the increased efficacy of a nanoparticle albumin-bound preparation of the antimicrotubular drug Paclitaxel (nab-paclitaxel). METHODS: The transcriptomes of 19 ACC tumors and 4 normal adrenal glands were profiled on Affymetrix U133 Plus2 expression microarrays to identify genes representing potential therapeutic targets. Immunohistochemical analysis for target proteins was performed on 10 ACC, 6 benign adenomas, and 1 normal adrenal gland. Agents known to inhibit selected targets were tested in comparison with mitotane in the 2 ACC cell lines (H295R and SW-13) in vitro and in mouse xenografts. RESULTS: SPARC expression is increased in ACC samples by 1.56 ± 0.44 (µ ± SD) fold. Paclitaxel and nab-paclitaxel show in vitro inhibition of H295R and SW-13 cells at IC50 concentrations of 0.33 µM and 0.0078 µM for paclitaxel and 0.35 µM and 0.0087 µM for nab-paclitaxel compared with mitotane concentrations of 15.9 µM and 46.4 µM. In vivo nab-paclitaxel treatment shows a greater decrease in tumor weight in both xenograft models than mitotane. CONCLUSIONS: Biological insights garnered through expression profiling of ACC tumors suggest further investigation into the use of nab-paclitaxel for the treatment of ACC.

Predicting drug sensitivity and resistance: profiling ABC transporter genes in cancer cells.

For analysis of multidrug resistance, a major barrier to effective cancer chemotherapy, we profiled mRNA expression of the 48 known human ABC transporters in 60 diverse cancer cell lines (the NCI-60) used by the National Cancer Institute to screen for anticancer activity. The use of real-time RT-PCR avoided artifacts commonly encountered with microarray technologies. By correlating the results with the growth inhibitory profiles of 1,429 candidate anticancer drugs tested against the cells, we identified which transporters are more likely than others to confer resistance to which agents. Unexpectedly, we also found and validated compounds whose activity is potentiated, rather than antagonized, by the MDR1 multidrug transporter. Such compounds may serve as leads for development.

Identification of a signature of evolutionarily conserved stress-induced mutagenesis in cancer.

The clustering of mutations observed in cancer cells is reminiscent of the stress-induced mutagenesis (SIM) response in bacteria. Bacteria deploy SIM when faced with DNA double-strand breaks in the presence of conditions that elicit an SOS response. SIM employs DinB, the evolutionary precursor to human trans-lesion synthesis (TLS) error-prone polymerases, and results in mutations concentrated around DNA double-strand breaks with an abundance that decays with distance. We performed a quantitative study on single nucleotide variant calls for whole-genome sequencing data from 1950 tumors, non-inherited mutations from 129 normal samples, and acquired mutations in 3 cell line models of stress-induced adaptive mutation. We introduce statistical methods to identify mutational clusters, quantify their shapes and tease out the potential mechanism that produced them. Our results show that mutations in both normal and cancer samples are indeed clustered and have shapes indicative of SIM. Clusters in normal samples occur more often in the same genomic location across samples than in cancer suggesting loss of regulation over the mutational process during carcinogenesis. Additionally, the signatures of TLS contribute the most to mutational cluster formation in both patient samples as well as experimental models of SIM. Furthermore, a measure of cluster shape heterogeneity was associated with cancer patient survival with a hazard ratio of 5.744 (Cox Proportional Hazard Regression, 95% CI: 1.824-18.09). Our results support the conclusion that the ancient and evolutionary-conserved adaptive mutation response found in bacteria is a source of genomic instability in cancer. Biological adaptation through SIM might explain the ability of tumors to evolve in the face of strong selective pressures such as treatment and suggests that the conventional 'hit it hard' approaches to therapy could prove themselves counterproductive.

Toward a pathway-centered approach for the treatment of adrenocortical carcinoma.

PURPOSE OF REVIEW: Adrenocortical carcinoma is an aggressive, lethal malignancy of the adrenal cortex. The rarity of the disease has stymied therapeutic development. Recent work toward understanding the molecular pathogenesis of the disease has identified several potential new diagnostic and therapeutic targets. RECENT FINDINGS: The molecular characterization of adrenocortical carcinoma has identified dysregulation of the Gap 2/mitosis transition and the insulin-like growth factor 1 receptor signaling cascade as two major pathways for therapeutic development. These studies have also highlighted an unappreciated heterogeneity of the disease at the gene level that nevertheless seems to converge onto common cellular pathways. Additionally, the characterization of Wnt signaling through ß-catenin in adrenal development, the demonstration of the involvement of BMP signaling in adrenocortical carcinoma growth regulation, and the discovery that ERCC1 expression levels can predict therapeutic response to platinum are just a few of the recent advances that promise to shed light on adrenocortical carcinoma biology. SUMMARY: Short-term, therapeutic development should target the Gap 2/mitosis transition and the downstream signaling of the insulin-like growth factor 1 receptor receptor. Long-term, additional characterization of patient samples, particularly at the sequence level, is required to fully understand adrenocortical carcinoma biology and apply that knowledge to clinical practice.

Genomic signatures of cancer: basis for individualized risk assessment, selective staging and therapy.

The development of DNA microarray and quantitative real-time polymerase chain reaction technologies has allowed for precise quantitation of RNA expression of hundreds to thousands of genes. These technologies have significantly impacted the study and understanding of cancer in terms of its molecular and genetic characteristics. In this review article, breast cancer, colon cancer, and adrenal carcinoma have been chosen to illustrate the principle and techniques of genomic profiling and to illustrate how such methods may be used to develop genomic signatures for personalized risk assessment and to individualize patient treatment.

Reverting to single-cell biology: The predictions of the atavism theory of cancer.

Cancer or cancer-like phenomena pervade multicellular life, implying deep evolutionary roots. Many of the hallmarks of cancer recapitulate unicellular modalities, suggesting that cancer initiation and progression represent a systematic reversion to simpler ancestral phenotypes in response to a stress or insult. This so-called atavism theory may be tested using phylostratigraphy, which can be used to assign ages to genes. Several research groups have confirmed that cancer cells tend to over-express evolutionary older genes, and rewire the architecture linking unicellular and multicellular gene networks. In addition, some of the elevated mutation rate - a well-known hallmark of cancer - is actually self-inflicted, driven by genes found to be homologs of the ancient SOS genes activated in stressed bacteria, and employed to evolve biological workarounds. These findings have obvious implications for therapy.

MYBPC3 Haplotype Linked to Hypertrophic Cardiomyopathy in Rhesus Macaques (Macaca mulatta).

In humans, abnormal thickening of the left ventricle of the heart clinically defines hypertrophic cardiomyopathy (HCM), a common inherited cardiovascular disorder that can precede a sudden cardiac death event. The wide range of clinical presentations in HCM obscures genetic variants that may influence an individual's susceptibility to sudden cardiac death. Although exon sequencing of major sarcomere genes can be used to detect high-impact causal mutations, this strategy is successful in only half of patient cases. The incidence of left ventricular hypertrophy (LVH) in a managed research colony of rhesus macaques provides an excellent comparative model in which to explore the genomic etiology of severe HCM and sudden cardiac death. Because no rhesus HCM-associated mutations have been reported, we used a next-generation genotyping assay that targets 7 sarcomeric rhesus genes within 63 genomic sites that are orthologous to human genomic regions known to harbor HCM disease variants. Amplicon sequencing was performed on 52 macaques with confirmed LVH and 42 unrelated, unaffected animals representing both the Indian and Chinese rhesus macaque subspecies. Bias-reduced logistic regression uncovered a risk haplotype in the rhesus MYBPC3 gene, which is frequently disrupted in both human and feline HCM; this haplotype implicates an intronic variant strongly associated with disease in either homozygous or carrier form. Our results highlight that leveraging evolutionary genomic data provides a unique, practical strategy for minimizing population bias in complex disease studies.

Genomic and expression profiling of adrenocortical carcinoma: application to diagnosis, prognosis and treatment.

Adrenocortical carcinoma (ACC) is an aggressive endocrine tumor with a poor 5-year survival rate of 10-20%. Current therapy is often ineffective and may be associated with intolerable side effects. Although ACC is extremely rare, recent advances in genomic and expression profiling, coupled with knowledge gained from the study of the inherited syndromes that increase ACC risk, are beginning to bring together a picture of a tumor type dependent on p53, the G2/M cell cycle transition and IGF2 stimulation. Nevertheless, ACC remains a heterogeneous disease. Only recently have sufficient tumors been characterized and results published to permit an exploration of this diversity. Advances in treatment will depend on exploiting those pathways already implicated in ACC, along with those yet to be identified, and testing those treatments in better models of the disease than the three cell lines that currently exist and are widely available to the community.

Asparagine synthetase as a causal, predictive biomarker for L-asparaginase activity in ovarian cancer cells.

L-Asparaginase (l-ASP), a bacterial enzyme used since the 1970s to treat acute lymphoblastic leukemia, selectively starves cells that cannot synthesize sufficient asparagine for their own needs. Molecular profiling of the NCI-60 cancer cell lines using five different microarray platforms showed strong negative correlations of asparagine synthetase (ASNS) expression and DNA copy number with sensitivity to l-ASP in the leukemia and ovarian cancer cell subsets. To assess whether the ovarian relationship is causal, we used RNA interference to silence ASNS in three ovarian lines and observed 4- to 5-fold potentiation of sensitivity to l-ASP with two of the lines. For OVCAR-8, the line that expresses the least ASNS, the potentiation was >500-fold. Significantly, that potentiation was >700-fold in the multidrug-resistant derivative OVCAR-8/ADR, showing that the causal relationship between ASNS expression and l-ASP activity survives development of classical multidrug resistance. Tissue microarrays confirmed low ASNS expression in a subset of clinical ovarian cancers as well as other tumor types. Overall, this pharmacogenomic/pharmacoproteomic study suggests the use of l-ASP for treatment of a subset of ovarian cancers (and perhaps other tumor types), with ASNS as a biomarker for patient selection.

Integrating data on DNA copy number with gene expression levels and drug sensitivities in the NCI-60 cell line panel.

Chromosome rearrangement, a hallmark of cancer, has profound effects on carcinogenesis and tumor phenotype. We used a panel of 60 human cancer cell lines (the NCI-60) as a model system to identify relationships among DNA copy number, mRNA expression level, and drug sensitivity. For each of 64 cancer-relevant genes, we calculated all 4,096 possible Pearson's correlation coefficients relating DNA copy number (assessed by comparative genomic hybridization using bacterial artificial chromosome microarrays) and mRNA expression level (determined using both cDNA and Affymetrix oligonucleotide microarrays). The analysis identified an association of ERBB2 overexpression with 3p copy number, a finding supported by data from human tumors and a mouse model of ERBB2-induced carcinogenesis. When we examined the correlation between DNA copy number for all 353 unique loci on the bacterial artificial chromosome microarray and drug sensitivity for 118 drugs with putatively known mechanisms of action, we found a striking negative correlation (-0.983; 95% bootstrap confidence interval, -0.999 to -0.899) between activity of the enzyme drug L-asparaginase and DNA copy number of genes near asparagine synthetase in the ovarian cancer cells. Previous analysis of drug sensitivity and mRNA expression had suggested an inverse relationship between mRNA levels of asparagine synthetase and L-asparaginase sensitivity in the NCI-60. The concordance of pharmacogenomic findings at the DNA and mRNA levels strongly suggests further study of L-asparaginase for possible treatment of a low-synthetase subset of clinical ovarian cancers. The DNA copy number database presented here will enable other investigators to explore DNA transcript-drug relationships in their own domains of research focus.

Ancient genes establish stress-induced mutation as a hallmark of cancer.

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.

Karyotypic complexity of the NCI-60 drug-screening panel.

We used spectral karyotyping to provide a detailed analysis of karyotypic aberrations in the diverse group of cancer cell lines established by the National Cancer Institute for the purpose of anticancer drug discovery. Along with the karyotypic description of these cell lines we defined and studied karyotypic complexity and heterogeneity (metaphase-to-metaphase variations) based on three separate components of genomic anatomy: (a) ploidy; (b) numerical changes; and (c) structural rearrangements. A wide variation in these parameters was evident in these cell lines, and different association patterns between them were revealed. Analysis of the breakpoints and other specific features of chromosomal changes across the entire set of cell lines or within particular lineages pointed to a striking lability of centromeric regions that distinguishes the epithelial tumor cell lines. We have also found that balanced translocations are as frequent in absolute number within the cell lines derived from solid as from hematopoietic tumors. Important similarities were noticed between karyotypic changes in cancer cell lines and that seen in primary tumors. This dataset offers insights into the causes and consequences of the destabilizing events and chromosomal instability that may occur during tumor development and progression. It also provides a foundation for investigating associations between structural genome anatomy and cancer molecular markers and targets, gene expression, gene dosage, and resistance or sensitivity to tens of thousands of molecular compounds.

Membrane transporters and channels: role of the transportome in cancer chemosensitivity and chemoresistance.

Membrane transporters and channels (collectively the transportome) govern cellular influx and efflux of ions, nutrients, and drugs. We used oligonucleotide arrays to analyze gene expression of the transportome in 60 human cancer cell lines used by the National Cancer Institute for drug screening. Correlating gene expression with the potencies of 119 standard anticancer drugs identified known drug-transporter interactions and suggested novel ones. Folate, nucleoside, and amino acid transporters positively correlated with chemosensitivity to their respective drug substrates. We validated the positive correlation between SLC29A1 (nucleoside transporter ENT1) expression and potency of nucleoside analogues, azacytidine and inosine-glycodialdehyde. Application of an inhibitor of SLC29A1, nitrobenzylmercaptopurine ribonucleoside, significantly reduced the potency of these two drugs, indicating that SLC29A1 plays a role in cellular uptake. Three ABC efflux transporters (ABCB1, ABCC3, and ABCB5) showed significant negative correlations with multiple drugs, suggesting a mechanism of drug resistance. ABCB1 expression correlated negatively with potencies of 19 known ABCB1 substrates and with Baker's antifol and geldanamycin. Use of RNA interference reduced ABCB1 mRNA levels and concomitantly increased sensitivity to these two drugs, as expected for ABCB1 substrates. Similarly, specific silencing of ABCB5 by small interfering RNA increased sensitivity to several drugs in melanoma cells, implicating ABCB5 as a novel chemoresistance factor. Ion exchangers, ion channels, and subunits of proton and sodium pumps variably correlated with drug potency. This study identifies numerous potential drug-transporter relationships and supports a prominent role for membrane transport in determining chemosensitivity. Measurement of transporter gene expression may prove useful in predicting anticancer drug response.

Diagnostic markers that distinguish colon and ovarian adenocarcinomas: identification by genomic, proteomic, and tissue array profiling.

Colon and ovarian cancers can be difficult to distinguish in the abdomen, and the distinction is important because it determines which drugs will be used for therapy. To identify molecular markers for that differential diagnosis, we developed a multistep protocol starting with the 60 human cancer cell lines used by the National Cancer Institute to screen for new anticancer agents. The steps included: (a) identification of candidate markers using cDNA microarrays; (b) verification of clone identities by resequencing; (c) corroboration of transcript levels using Affymetrix oligonucleotide chips; (d) quantitation of protein expression by "reverse-phase" protein microarray; and (e) prospective validation of candidate markers on clinical tumor sections in tissue microarrays. The two best candidates identified were villin for colon cancer cells and moesin for ovarian cancer cells. Because moesin stained stromal elements in both types of cancer, it would probably not have been identified as a marker if we had started with mRNA or protein profiling of bulk tumors. Villin appears at least as useful as the currently used colon cancer marker cytokeratin 20, and moesin also appears to have utility. The multistep process introduced here has the potential to produce additional markers for cancer diagnosis, prognosis, and therapy.

Pathway Implications of Aberrant Global Methylation in Adrenocortical Cancer.

CONTEXT: Adrenocortical carcinomas (ACC) are a rare tumor type with a poor five-year survival rate and limited treatment options. OBJECTIVE: Understanding of the molecular pathogenesis of this disease has been aided by genomic analyses highlighting alterations in TP53, WNT, and IGF signaling pathways. Further elucidation is needed to reveal therapeutically actionable targets in ACC. DESIGN: In this study, global DNA methylation levels were assessed by the Infinium HumanMethylation450 BeadChip Array on 18 ACC tumors and 6 normal adrenal tissues. A new, non-linear correlation approach, the discretization method, assessed the relationship between DNA methylation/gene expression across ACC tumors. RESULTS: This correlation analysis revealed epigenetic regulation of genes known to modulate TP53, WNT, and IGF signaling, as well as silencing of the tumor suppressor MARCKS, previously unreported in ACC. CONCLUSIONS: DNA methylation may regulate genes known to play a role in ACC pathogenesis as well as known tumor suppressors.

Vorinostat differentially alters 3D nuclear structure of cancer and non-cancerous esophageal cells.

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.