Cell Dispersal Influences Tumor Heterogeneity and Introduces a Bias in NGS Data Interpretation

Short and long distance cell dispersal can have a marked effect on tumor structure, high cellular motility could lead to faster cell mixing and lower observable intratumor heterogeneity. Here we evaluated a model for cell mixing that investigates how short-range dispersal and cell turnover will account for mutational proportions. We show that cancer cells can penetrate neighboring and distinct areas in a matter of days. In next generation sequencing runs, higher proportions of a given cell line generated frequencies with higher precision, while mixtures with lower amounts of each cell line had lower precision manifesting in higher standard deviations. When multiple cell lines were co-cultured, cellular movement altered observed mutation frequency by up to 18.5%. We propose that some of the shared mutations detected at low allele frequencies represent highly motile clones that appear in multiple regions of a tumor owing to dispersion throughout the tumor. In brief, cell movement will lead to a significant technical (sampling) bias when using next generation sequencing to determine clonal composition. A possible solution to this drawback would be to radically decrease detection thresholds and increase coverage in NGS analyses. © 2017 The Author(s)

Dynamics of cancer cell subpopulations in primary and metastatic colorectal tumors

Intratumor heterogeneity—heterogeneity of cancer cells within a single tumor—is considered one of the most problematic factors of treatment. Genetic heterogeneity, such as in somatic mutations and chromosome aberrations, is a common characteristic of human solid tumors and is probably the basis of biological heterogeneity. Using mutations in APC, TP53 and KRAS as markers to identify distinct colorectal cancer subpopulations, we analyzed a total of 42 primary colorectal cancer tissues and six paired liver metastases with multipoint microsampling, which enabled analysis of mutation patterns and allelic imbalances with a resolution of 0.01 mm2 (about 200 cells). There was usually more than one subpopulation in each primary tumor. Only two of 15 (13.3%) cases with three gene mutations and eight of 27 (29.6%) cases with two gene mutations had a single subpopulation. Cells with mutations in all of the examined genes usually constituted the major population. Multipoint microsampling of six primary and metastatic tumor pairs revealed that the majority of discrepancies in mutation patterns found with the bulk tissue analysis were due to loss of subpopulations in the metastatic tissues. In addition, multipoint microsampling uncovered substantial changes in subpopulations that were not detected with bulk tissue analysis. Specifically, the proportion of KRAS mutation-negative subpopulations increased in the metastatic tumors of four cases. Because KRAS mutation status is linked to cetuximab/panitumumab efficacy, subpopulation dynamics could lead to differences in response to cetuximab/panitumumab in primary versus metastatic tumors

Protein crystals in adenovirus type 5-infected cells: requirements for intranuclear crystallogenesis, structural and functional analysis

Intranuclear crystalline inclusions have been observed in the nucleus of epithelial cells infected with Adenovirus serotype 5 (Ad5) at late steps of the virus life cycle. Using immuno-electron microscopy and confocal microscopy of cells infected with various Ad5 recombinants modified in their penton base or fiber domains, we found that these inclusions represented crystals of penton capsomers, the heteromeric capsid protein formed of penton base and fiber subunits. The occurrence of protein crystals within the nucleus of infected cells required the integrity of the fiber knob and part of the shaft domain. In the knob domain, the region overlapping residues 489–492 in the FG loop was found to be essential for crystal formation. In the shaft, a large deletion of repeats 4 to 16 had no detrimental effect on crystal inclusions, whereas deletion of repeats 8 to 21 abolished crystal formation without altering the level of fiber protein expression. This suggested a crucial role of the five penultimate repeats in the crystallisation process. Chimeric pentons made of Ad5 penton base and fiber domains from different serotypes were analyzed with respect to crystal formation. No crystal was found when fiber consisted of shaft (S) from Ad5 and knob (K) from Ad3 (heterotypic S5-K3 fiber), but occurred with homotypic S3K3 fiber. However, less regular crystals were observed with homotypic S35-K35 fiber. TB5, a monoclonal antibody directed against the Ad5 fiber knob was found by immunofluorescence microscopy to react with high efficiency with the intranuclear protein crystals in situ. Data obtained with Ad fiber mutants indicated that the absence of crystalline inclusions correlated with a lower infectivity and/or lower yields of virus progeny, suggesting that the protein crystals might be involved in virion assembly. Thus, we propose that TB5 staining of Ad-infected 293 cells can be used as a prognostic assay for the viability and productivity of fiber-modified Ad5 vectors

Tumor innate immunity primed by specific interferon-stimulated endogenous retroviruses.

Mesenchymal tumor subpopulations secrete pro-tumorigenic cytokines and promote treatment resistance1-4. This phenomenon has been implicated in chemorefractory small cell lung cancer and resistance to targeted therapies5-8, but remains incompletely defined. Here, we identify a subclass of endogenous retroviruses (ERVs) that engages innate immune signaling in these cells. Stimulated 3 prime antisense retroviral coding sequences (SPARCS) are oriented inversely in 3' untranslated regions of specific genes enriched for regulation by STAT1 and EZH2. Derepression of these loci results in double-stranded RNA generation following IFN-γ exposure due to bi-directional transcription from the STAT1-activated gene promoter and the 5' long terminal repeat of the antisense ERV. Engagement of MAVS and STING activates downstream TBK1, IRF3, and STAT1 signaling, sustaining a positive feedback loop. SPARCS induction in human tumors is tightly associated with major histocompatibility complex class 1 expression, mesenchymal markers, and downregulation of chromatin modifying enzymes, including EZH2. Analysis of cell lines with high inducible SPARCS expression reveals strong association with an AXL/MET-positive mesenchymal cell state. While SPARCS-high tumors are immune infiltrated, they also exhibit multiple features of an immune-suppressed microenviroment. Together, these data unveil a subclass of ERVs whose derepression triggers pathologic innate immune signaling in cancer, with important implications for cancer immunotherapy

Aggregation and travelling wave dynamics in a two-population model of cancer cell growth and invasion

Funding: Engineering and Physical Sciences Research Council (UK) grant numbers EP/L504932/1 (VB), EP/K033689/1 (RE).Cells adhere to each other and to the extracellular matrix (ECM) through protein molecules on the surface of the cells. The breaking and forming of adhesive bonds, a process critical in cancer invasion and metas- tasis, can be influenced by the mutation of cancer cells. In this paper, we develop a nonlocal mathematical model describing cancer cell invasion and movement as a result of integrin-controlled cell-cell adhesion and cell-matrix adhesion, for two cancer cell populations with different levels of mutation. The partial differential equations for cell dynamics are coupled with ordinary differential equations describing the extracellular matrix (ECM) degradation and the production and decay of integrins. We use this model to investigate the role of cancer mutation on the possibility of cancer clonal competition with alternating dominance, or even competitive exclusion (phenomena observed experimentally). We discuss different possible cell aggregation patterns, as well as travelling wave patterns. In regard to the travelling waves, we investigate the effect of cancer mutation rate on the speed of cancer invasion.Publisher PDFPeer reviewe

Tumour Cell Heterogeneity.

The population of cells that make up a cancer are manifestly heterogeneous at the genetic, epigenetic, and phenotypic levels. In this mini-review, we summarise the extent of intra-tumour heterogeneity (ITH) across human malignancies, review the mechanisms that are responsible for generating and maintaining ITH, and discuss the ramifications and opportunities that ITH presents for cancer prognostication and treatment

Regional differences in prostaglandin E₂ metabolism in human colorectal cancer liver metastases

Background: Prostaglandin (PG) E₂ plays a critical role in colorectal cancer (CRC) progression, including epithelial-mesenchymal transition (EMT). Activity of the rate-limiting enzyme for PGE₂ catabolism (15-hydroxyprostaglandin dehydrogenase [15-PGDH]) is dependent on availability of NAD+. We tested the hypothesis that there is intra-tumoral variability in PGE₂ content, as well as in levels and activity of 15-PGDH, in human CRC liver metastases (CRCLM). To understand possible underlying mechanisms, we investigated the relationship between hypoxia, 15-PGDH and PGE₂ in human CRC cells in vitro. Methods: Tissue from the periphery and centre of 20 human CRCLM was analysed for PGE₂ levels, 15-PGDH and cyclooxygenase (COX)-2 expression, 15-PGDH activity, and NAD+/NADH levels. EMT of LIM1863 human CRC cells was induced by transforming growth factor (TGF) β. Results: PGE₂ levels were significantly higher in the centre of CRCLM compared with peripheral tissue (P = 0.04). There were increased levels of 15-PGDH protein in the centre of CRCLM associated with reduced 15-PGDH activity and low NAD+/NADH levels. There was no significant heterogeneity in COX-2 protein expression. NAD+ availability controlled 15-PGDH activity in human CRC cells in vitro. Hypoxia induced 15-PGDH expression in human CRC cells and promoted EMT, in a similar manner to PGE₂. Combined 15-PGDH expression and loss of membranous E-cadherin (EMT biomarker) were present in the centre of human CRCLM in vivo.Conclusions: There is significant intra-tumoral heterogeneity in PGE₂ content, 15-PGDH activity and NAD+ availability in human CRCLM. Tumour micro-environment (including hypoxia)-driven differences in PGE₂ metabolism should be targeted for novel treatment of advanced CRC

Parallel Evolution under Chemotherapy Pressure in 29 Breast Cancer Cell Lines Results in Dissimilar Mechanisms of Resistance

Background: Developing chemotherapy resistant cell lines can help to identify markers of resistance. Instead of using a panel of highly heterogeneous cell lines, we assumed that truly robust and convergent pattern of resistance can be identified in multiple parallel engineered derivatives of only a few parental cell lines. Methods: Parallel cell populations were initiated for two breast cancer cell lines (MDA-MB-231 and MCF-7) and these were treated independently for 18 months with doxorubicin or paclitaxel. IC50 values against 4 chemotherapy agents were determined to measure cross-resistance. Chromosomal instability and karyotypic changes were determined by cytogenetics. TaqMan RT-PCR measurements were performed for resistance-candidate genes. Pgp activity was measured by FACS. Results: All together 16 doxorubicin- and 13 paclitaxel-treated cell lines were developed showing 2-46 fold and 3-28 fold increase in resistance, respectively. The RT-PCR and FACS analyses confirmed changes in tubulin isofom composition, TOP2A and MVP expression and activity of transport pumps (ABCB1, ABCG2). Cytogenetics showed less chromosomes but more structural aberrations in the resistant cells. Conclusion: We surpassed previous studies by parallel developing a massive number of cell lines to investigate chemoresistance. While the heterogeneity caused evolution of multiple resistant clones with different resistance characteristics, the activation of only a few mechanisms were sufficient in one cell line to achieve resistance. © 2012 Tegze et al