Proteome Heterogeneity in Colorectal Cancer

Tumor heterogeneity is an important feature of colorectal cancer (CRC) manifested by dynamic changes in gene expression, protein expression, and availability of different tumor subtypes. Recent publications in the past 10 years have revealed proteome heterogeneity between different colorectal tumors and within the same tumor site. This paper reviews recent research works on the proteome heterogeneity in CRC, which includes the heterogeneity within a single tumor (intratumor heterogeneity), between different anatomical sites at the same organ, and between primary and metastatic sites (intertumor heterogeneity). The potential use of proteome heterogeneity in precision medicine and its implications in biomarker discovery and therapeutic outcomes will be discussed. Identification of the unique proteome landscape between and within individual tumors is imperative for understanding cancer biology and the management of CRC patients.     

Intratumor heterogeneity in evolutionary models of tumor progression

With rare exceptions, human tumors arise from single cells that have accumulated the necessary number and types of heritable alterations. Each such cell leads to dysregulated growth and eventually the formation of a tumor. Despite their monoclonal origin, at the time of diagnosis most tumors show a striking amount of intratumor heterogeneity in all measurable phenotypes; such heterogeneity has implications for diagnosis, treatment efficacy, and the identification of drug targets. An understanding of the extent and evolution of intratumor heterogeneity is therefore of direct clinical importance. In this article, we investigate the evolutionary dynamics of heterogeneity arising during exponential expansion of a tumor cell population, in which heritable alterations confer random fitness changes to cells. We obtain analytical estimates for the extent of heterogeneity and quantify the effects of system parameters on this tumor trait. Our work contributes to a mathematical understanding of intratumor heterogeneity and is also applicable to organisms like bacteria, agricultural pests, and other microbes.     

Intratumor heterogeneity: Nature and biological significance

Intratumor heterogeneity inherent in the majority of human cancers is a major obstacle for a highly efficient diagnosis and successful prognosis and treatment of these diseases. Being a result of clonal diversity within the same tumor, intratumor heterogeneity can be manifested in variability of genetic and epigenetic status, gene and protein expression, morphological structure, and other features of the tumor. It is most likely that the appearance of this diversity is a source for the adaptation of the tumor to changes in microenvironmental conditions and/or a tool for changing its malignant potential. In any case, both processes result in the appearance of cell clones with different undetermined sets of hallmarks. In this review, we describe the heterogeneity of molecular disorders in various human tumors and consider modern viewpoints of its development including genetic and non-genetic factors of heterogeneity origin and the role of cancer stem cells and clonal evolution. We also systematize data on the contribution of tumor diversity to progression of various tumors and the efficiency of their treatment. The main problems are indicated in the diagnosis and therapy of malignant tumors caused by intratumor heterogeneity and possible pathways for their solution. Moreover, we also suggest the key goals whose achievement promises to minimize the problem of intratumor heterogeneity and to identify new prognostic, predictive, and target markers for adequate and effective treatment of cancer.     

Intratumor heterogeneity of k-ras and p53 mutations among human colorectal adenomas containing early cancer.

The molecular pathways and the timing of genetic events during human colorectal carcinogenesis are still not fully understood. We have addressed the intratumor heterogeneity of the mutational status of the k-ras oncogene and of the p53 oncosuppressor gene during the adenoma-carcinoma sequence by investigating 26 human colorectal adenomas containing early cancer. An intratumor comparative analysis was obtained among the adenomatous and carcinomatous component pairs. Additionally, we have analyzed 17 adenomas having cancer in the near vicinity. The adenomatous components of the adenomas containing early cancer and the adenomas having cancer in the near vicinity had comparable frequencies for k-ras mutations (28 and 47%) but different for p53 mutations (52 and 7%, p-value = 0.01). Interestingly, the adenomatous and carcinomatous components of the adenomas containing early cancer were rarely heterogeneous for the k-ras mutational status (only in 13% of the cases) but were characterized by heterogeneity of the p53 status in 59% of the cases (p-value < 0.01). In addition, the mutations of p53 for the adenomatous components of the adenomas containing early cancer were statistically significantly associated with severe dysplasia (p-value = 0.01). Intratumor homogeneity of k-ras status during the human colorectal adenoma-carcinoma sequence suggests that the role of k-ras is more related to tumor initiation than to tumor progression. On the contrary, intratumor heterogeneity of p53 mutations indicates that the type of the p53 mutations may also be relevant for selection and expansion of new subclones leading to tumor progression.     

Integrated Multiregional Analysis Proposing a New Model of Colorectal Cancer Evolution

Understanding intratumor heterogeneity is clinically important because it could cause therapeutic failure by fostering evolutionary adaptation. To this end, we profiled the genome and epigenome in multiple regions within each of nine colorectal tumors. Extensive intertumor heterogeneity is observed, from which we inferred the evolutionary history of the tumors. First, clonally shared alterations appeared, in which C>T transitions at CpG site and CpG island hypermethylation were relatively enriched. Correlation between mutation counts and patients’ ages suggests that the early-acquired alterations resulted from aging. In the late phase, a parental clone was branched into numerous subclones. Known driver alterations were observed frequently in the early-acquired alterations, but rarely in the late-acquired alterations. Consistently, our computational simulation of the branching evolution suggests that extensive intratumor heterogeneity could be generated by neutral evolution. Collectively, we propose a new model of colorectal cancer evolution, which is useful for understanding and confronting this heterogeneous disease.     

DDX3X Induces Primary EGFR-TKI Resistance Based on Intratumor Heterogeneity in Lung Cancer Cells Harboring EGFR-Activating Mutations

The specific mechanisms how lung cancer cells harboring epidermal growth factor receptor (EGFR) activating mutations can survive treatment with EGFR-tyrosine kinase inhibitors (TKIs) until they eventually acquire treatment-resistance genetic mutations are unclear. The phenotypic diversity of cancer cells caused by genetic or epigenetic alterations (intratumor heterogeneity) confers treatment failure and may foster tumor evolution through Darwinian selection. Recently, we found DDX3X as the protein that was preferentially expressed in murine melanoma with cancer stem cell (CSC)-like phenotypes by proteome analysis. In this study, we transfected PC9, human lung cancer cells harboring EGFR exon19 deletion, with cDNA encoding DDX3X and found that DDX3X, an ATP-dependent RNA helicase, induced CSC-like phenotypes and the epithelial-mesenchymal transition (EMT) accompanied with loss of sensitivity to EGFR-TKI. DDX3X expression was associated with upregulation of Sox2 and increase of cancer cells exhibiting CSC-like phenotypes, such as anchorage-independent proliferation, strong expression of CD44, and aldehyde dehydrogenase (ALDH). The EMT with switching from E-cadherin to N-cadherin was also facilitated by DDX3X. Either ligand-independent or ligand-induced EGFR phosphorylation was inhibited in lung cancer cells that strongly expressed DDX3X. Lack of EGFR signal addiction resulted in resistance to EGFR-TKI. Moreover, we found a small nonadherent subpopulation that strongly expressed DDX3X accompanied by the same stem cell-like properties and the EMT in parental PC9 cells. The unique subpopulation lacked EGFR signaling and was highly resistant to EGFR-TKI. In conclusion, our data indicate that DDX3X may play a critical role for inducing phenotypic diversity, and that treatment targeting DDX3X may overcome primary resistance to EGFR-TKI resulting from intratumor heterogeneity.     

Deciphering intratumor heterogeneity using cancer genome analysis

Intratumor heterogeneity within individual cancer tissues underlies the numerous phenotypes of cancer. Tumor subclones ultimately affect therapeutic outcomes due to their distinct molecular features. Drug-resistant subclones are present at a low frequency in tissues at the time of biopsy, but can also arise as a result of acquired somatic mutations. A number of different approaches have been utilized to understand the nature of intratumor heterogeneity. Clonal analysis using whole exome or genome sequencing data can help monitor subclones in the context of tumor progression. Multiregional biopsies permit the molecular characterization of subclones within tumors. Deep sequencing has also provided researchers with the ability to measure the low allele fraction variant within a small number of cells. Ultimately, single-cell sequencing will enable the identification of every minor population within a tumor microenvironment. In the clinical context, the ability to identify and monitor the subclonal architecture of a tumor is valuable for the development of precise cancer therapeutic methods.     

The developing landscape of combinatorial therapies of immune checkpoint blockade with DNA damage repair inhibitors for the treatment of breast and ovarian cancers

In this era of precision medicine, with the help of biomarkers, immunotherapy has significantly improved prognosis of many patients with malignant tumor. Deficient mismatch repair (dMMR)/microsatellite instability (MSI) status is used as a biomarker in clinical practice to predict favorable response to immunotherapy and prognosis. MSI is an important characteristic which facilitates mutation and improves the likelihood of a favorable response to immunotherapy. However, many patients with dMMR/MSI still respond poorly to immunotherapies, which partly results from intratumor heterogeneity propelled by dMMR/MSI. In this review, we discuss how dMMR/MSI facilitates mutations in tumor cells and generates intratumor heterogeneity, especially through type II interferon (IFN-γ) signaling and tumor-infiltrating lymphocytes (TILs). We discuss the mechanism of immunotherapy from the perspective of dMMR/MSI, molecular pathways and TILs, and we discuss how intratumor heterogeneity hinders the therapeutic effect of immunotherapy. Finally, we summarize present techniques and strategies to look at the tumor as a whole to design personalized regimes and achieve favorable prognosis.

     

Intratumor Heterogeneity of ALK-Rearrangements and Homogeneity of EGFR-Mutations in Mixed Lung Adenocarcinoma

Background

Non Small Cell Lung Cancer is a highly heterogeneous tumor. Histologic intratumor heterogeneity could be ‘major’, characterized by a single tumor showing two different histologic types, and ‘minor’, due to at least 2 different growth patterns in the same tumor. Therefore, a morphological heterogeneity could reflect an intratumor molecular heterogeneity. To date, few data are reported in literature about molecular features of the mixed adenocarcinoma. The aim of our study was to assess EGFR-mutations and ALK-rearrangements in different intratumor subtypes and/or growth patterns in a series of mixed adenocarcinomas and adenosquamous carcinomas.

Methods

590 Non Small Cell Lung Carcinomas tumor samples were revised in order to select mixed adenocarcinomas with available tumor components. Finally, only 105 mixed adenocarcinomas and 17 adenosquamous carcinomas were included in the study for further analyses. Two TMAs were built selecting the different intratumor histotypes. ALK-rearrangements were detected through FISH and IHC, and EGFR-mutations were detected through IHC and confirmed by RT-PCR.

Results

10/122 cases were ALK-rearranged and 7 from those 10 showing an intratumor heterogeneity of the rearrangements. 12/122 cases were EGFR-mutated, uniformly expressing the EGFR-mutated protein in all histologic components.

Conclusion

Our data suggests that EGFR-mutations is generally homogeneously expressed. On the contrary, ALK-rearrangement showed an intratumor heterogeneity in both mixed adenocarcinomas and adenosquamous carcinomas. The intratumor heterogeneity of ALK-rearrangements could lead to a possible impact on the therapeutic responses and the disease outcomes.     

Unsupervised Deconvolution of Dynamic Imaging Reveals Intratumor Vascular Heterogeneity and Repopulation Dynamics

With the existence of biologically distinctive malignant cells originated within the same tumor, intratumor functional heterogeneity is present in many cancers and is often manifested by the intermingled vascular compartments with distinct pharmacokinetics. However, intratumor vascular heterogeneity cannot be resolved directly by most in vivo dynamic imaging. We developed multi-tissue compartment modeling (MTCM), a completely unsupervised method of deconvoluting dynamic imaging series from heterogeneous tumors that can improve vascular characterization in many biological contexts. Applying MTCM to dynamic contrast-enhanced magnetic resonance imaging of breast cancers revealed characteristic intratumor vascular heterogeneity and therapeutic responses that were otherwise undetectable. MTCM is readily applicable to other dynamic imaging modalities for studying intratumor functional and phenotypic heterogeneity, together with a variety of foreseeable applications in the clinic.     

The thin red line between the immune system and cancer evolution

The cancer immunoediting theory describes the dual ability of endogenous antitumor immunity to inhibit or promote progressing cancers. Tumor-specific neoantigens arising from somatic mutations serve as targets for the endogenous T-cell-mediated antitumor immunity and therefore possess a crucial role for tumor development. Additionally, targeting these molecules is conceptually appealing because neoantigens are not expressed in healthy tissue and therefore confer less toxicity and greater specificity when used in therapeutic interventions. Moreover, intratumor neo-antigenic heterogeneity is believed to play a pivotal role in the activation of adaptive immunity and in the efficacy of immunotherapies that are based on immune checkpoint inhibition. In this respect, mutual interactions between tumor cells and immune lymphocytes regulate the levels of antitumor immunity, but also shape tumor heterogeneity through the selective outgrowth of tumor subclones. Therefore, the exploration of the mechanistic pathways and the identification of the genomic aberrations underlying the clonal evolution of tumors is considered mandatory for improving the clinical outcomes of therapies, as it will assist in the selection of the appropriate therapeutic decisions so as to delay, avoid, or overcome resistance through the identification of the most effective therapeutic strategies.     

Intratumoral variations in DNA ploidy and s-phase fraction in human breast cancer.

To study intratumoral DNA ploidy heterogeneity and S-phase fraction (SPF) variability, we prospectively collected five different samples from 48 breast carcinomas and each sample was analysed separately by flow cytometry. Aneuploidy rate was 89.6% after analysis of four or five samples. DNA ploidy heterogeneity, i.e., different samples classified as either DNA euploid or DNA aneuploid in the same tumor was seen in 17%, and DNA index heterogeneity, i.e., tumor populations with different DNA indices (DIs) seen in different samples was 44%. A statistical model defining SPF heterogeneity is proposed. SPF heterogeneity as defined by us was 71%, and as expected the SPF heterogeneity rate increased significantly with increasing number of analysed samples. Four or more samples are needed to detect the most deviant (highest) SPF values. An unrecognized intratumor heterogeneity of DNA ploidy and SPF may partly explain the conflicting results reported in the literature on the above prognostic indicators.     

Variations in pO2 and pH response to hyperthermia: dependence on transplant site and duration of treatment.

It has been clearly established that changes in intratumor pO2 and pH occur following hyperthermia, and it has been hypothesized that these changes may, in some way, be related to the ultimate response (i.e., cure) of the lesion. The purpose of this study was twofold: first, to examine the changes in intratumor pH during the course of a hyperthermia treatment at biologically related end point "doses"; second, to examine the response of pO2 after treatment in a different lesion transplant site. During hyperthermia treatment of the tumor transplanted in the leg, intratumor pH was found to drop from a control value of 6.74 +/- 0.17 to 6.47 +/- 0.13 within 15 min following the start of treatment. The values then remained relatively constant throughout the remainder of the treatment (either 1 or 2 h at 43.5 degrees C). Following the subcurative (10% tumor cures at 30 days; 60 min at 43.5 degrees C) treatment the pH began to rise immediately, while after the higher dose (60% tumor cures at 30 days; 120 min at 43.5 degrees C) a slight rise in pH was followed by a continuous drop in pH for up to 4 h, as we have reported previously. Oxygen response in the two transplant sites (leg and flank) was found to be remarkably different even though the tumor cure rate was identical for a given hyperthermia "dose" in terms of time and temperature. In the leg, only very low levels of oxygen can be measured in the tumor 24 h after treatment with either "dose" studied (all measured pO2 values less than or equal to 5 mm Hg). In the flank, the tumor response is dependent on hyperthermia "dose." Only 28% of measured oxygen values are less than or equal to 5 mm Hg 24 h following a subcurative "dose," while 4 h following the higher "dose" there is a nonsignificant trend toward hypoxia (approximately 65% of values less than or equal to 5 mm Hg) with a subsequent shift toward reoxygenation. These latter observations are contrary to results reported previously and tend to contradict some current theories regarding the physiological mechanisms associated with hyperthermia treatment.     

GeneCount: genome-wide calculation of absolute tumor DNA copy numbers from array comparative genomic hybridization data

Absolute tumor DNA copy numbers can currently be achieved only on a single gene basis by using fluorescence in situ hybridization (FISH). We present GeneCount, a method for genome-wide calculation of absolute copy numbers from clinical array comparative genomic hybridization data. The tumor cell fraction is reliably estimated in the model. Data consistent with FISH results are achieved. We demonstrate significant improvements over existing methods for exploring gene dosages and intratumor copy number heterogeneity in cancers.     

Algorithms to Model Single Gene,Single Chromosome,and Whole Genome Copy Number Changes Jointly in Tumor Phylogenetics

We present methods to construct phylogenetic models of tumor progression at the cellular level that include copy number changes at the scale of single genes, entire chromosomes, and the whole genome. The methods are designed for data collected by fluorescence in situ hybridization (FISH), an experimental technique especially well suited to characterizing intratumor heterogeneity using counts of probes to genetic regions frequently gained or lost in tumor development. Here, we develop new provably optimal methods for computing an edit distance between the copy number states of two cells given evolution by copy number changes of single probes, all probes on a chromosome, or all probes in the genome. We then apply this theory to develop a practical heuristic algorithm, implemented in publicly available software, for inferring tumor phylogenies on data from potentially hundreds of single cells by this evolutionary model. We demonstrate and validate the methods on simulated data and published FISH data from cervical cancers and breast cancers. Our computational experiments show that the new model and algorithm lead to more parsimonious trees than prior methods for single-tumor phylogenetics and to improved performance on various classification tasks, such as distinguishing primary tumors from metastases obtained from the same patient population.     

Intratumor heterogeneity of biomarker expression in breast carcinomas

Small biopsy samples are used increasingly to assess the biomarker expression for prognostic information and for monitoring therapeutic responses prior to and during neoadjuvant therapy. The issue of intratumor heterogeneity of expression of biomarkers, however, has raised questions about the validity of the assessment of biomarker expression based on limited tissue samples. We examined immunohistochemically the expression of HER-2neu (p185^erbB-2), epidermal growth factor receptor (EGFR), Bcl-2, p53, and proliferating cell nuclear antigen (PCNA) in 30 breast carcinomas using archived, paraffin embedded tissue and determined the extent of intratumor heterogeneity. Each section was divided into four randomly oriented discrete regions, each containing a portion of the infiltrating carcinoma. For each tumor, the entire lesion and four regions were analyzed for the expression of these markers. Scores of both membrane and cytoplasmic staining of HER-2neu and EGFR, scores of cytoplasmic staining of Bcl-2, and scores of nuclear staining of both p53 and PCNA were recorded. The intensity of staining and the proportion of immunostained cells were determined. A semiquantitative immunoscore was calculated by determining the sum of the products of the intensity and corresponding proportion of stained tumor cells. We analyzed both invasive (IDC) and in situ (DCIS) carcinomas. The Wilcoxon signed-rank test was used for paired comparisons between overall and regional immunoscores and between overall and regional percentages of stained cells. Spearman's correlation coefficients were used to assess the level of agreement of overall biomarker expression with each of the regions. Generalized linear models were used to assess overall and pair-wise differences in the absolute values of percent changes between overall and regional expression of biomarkers. For IDCs, there were no statistically significant differences in the expression of the biomarkers in terms of either the percentage of cells staining or the immunoscores when comparing the entire tumor with each region except for the lower EGFR expression of arbitrarily selected region 1 and lower p53 expression of region 1 compared to that of the entire tumor section. For DCIS, there were no statistically significant differences in the expression of the biomarkers between the entire tumor and each region except in PCNA of region 2 compared to that of entire tumor section. Positive correlation of immunoscores was observed between the entire tumor and each region as well as across all four regions for IDC. Similar observations were noted with DCIS except for HER-2neu and PCNA. No statistically significant differences were observed in the absolute values of percent changes of biomarker expression between overall and the four regions for both DCIS and IDC. Therefore, no significant intratumor heterogeneity in the expression of HER-2neu, Bcl-2, and PCNA was observed in IDC. Minor regional variations were observed for EGFR and p53 in IDC. Similarly, no significant regional variation in the expression of markers was observed in DCIS except for PCNA.     

Targeted therapies in pancreatic cancer: Promises and failures

Pancreatic ductal adenocarcinoma (PDAC) is an incidence rate nearly equal to its mortality rate. The poor prognosis of the disease can be explained by the absence of effective biomarkers for screening and early detection, together with the aggressive behavior and resistance to the currently available chemotherapy. The therapeutic failure can also be attributed to the inter-/intratumor genetic heterogeneity and the abundance of tumor stroma that occupies the majority of the tumor mass. Gemcitabine is used in the treatment of PDAC; however, the response rate is less than 12%. A recent phase III trial revealed that the combination of oxaliplatin, irinotecan, fluorouracil, and leucovorin could be an option for the treatment of metastatic PDAC patients with good performance status, although these approaches can result in high toxicity level. Further investigations are required to develop innovative anticancer agents that either improve gemcitabine activity, within novel combinatorial approaches or acts with a better efficacy than gemcitabine. The aim of the current review is to give an overview of preclinical and clinical studies targeting key dysregulated signaling pathways in PDAC.