A new model isolates glioblastoma clonal interactions and reveals unexpected modes for regulating motility, proliferation, and drug resistance

Signaling dynamics in coexisting monoclonal cell subpopulations unveil mechanisms of resistance to anti-cancer compounds

BACKGROUND

Tumor heterogeneity is a main contributor of resistance to anti-cancer targeted agents though it has proven difficult to study. Unfortunately, model systems to functionally characterize and mechanistically study dynamic responses to treatment across coexisting subpopulations of cancer cells remain a missing need in oncology.

METHODS

Using single cell cloning and expansion techniques, we established monoclonal cell subpopulations (MCPs) from a commercially available epidermal growth factor receptor (EGFR)-mutant non-small cell lung cancer cell line. We then used this model sensitivity to the EGFR inhibitor osimertinib across coexisting cell populations within the same tumor. Pathway-centered signaling dynamics associated with response to treatment and morphological characteristics of the MCPs were assessed using Reverse Phase Protein Microarray. Signaling nodes differentially activated in MCPs less sensitive to treatment were then pharmacologically inhibited to identify target signaling proteins putatively implicated in promoting drug resistance.

RESULTS

MCPs demonstrated highly heterogeneous sensitivities to osimertinib. Cell viability after treatment increased > 20% compared to the parental line in selected MCPs, whereas viability decreased by 75% in other MCPs. Reduced treatment response was detected in MCPs with higher proliferation rates, EGFR L858R expression, activation of EGFR binding partners and downstream signaling molecules, and expression of epithelial-to-mesenchymal transition markers. Levels of activation of EGFR binding partners and MCPs' proliferation rates were also associated with response to c-MET and IGFR inhibitors.

CONCLUSIONS

MCPs represent a suitable model system to characterize heterogeneous biomolecular behaviors in preclinical studies and identify and functionally test biological mechanisms associated with resistance to targeted therapeutics.

Comparative single-cell RNA-sequencing profiling of BMP4-treated primary glioma cultures reveals therapeutic markers

BACKGROUND

Glioblastoma (GBM) is the most aggressive primary brain tumor. Its cellular composition is very heterogeneous, with cells exhibiting stem-cell characteristics (GSCs) that co-determine therapy resistance and tumor recurrence. Bone Morphogenetic Protein (BMP)-4 promotes astroglial and suppresses oligodendrocyte differentiation in GSCs, processes associated with superior patient prognosis. We characterized variability in cell viability of patient-derived GBM cultures in response to BMP4 and, based on single-cell transcriptome profiling, propose predictive positive and early-response markers for sensitivity to BMP4.

METHODS

Cell viability was assessed in 17 BMP4-treated patient-derived GBM cultures. In two cultures, one highly-sensitive to BMP4 (high therapeutic efficacy) and one with low-sensitivity, response to treatment with BMP4 was characterized. We applied single-cell RNA-sequencing, analyzed the relative abundance of cell clusters, searched for and identified the aforementioned two marker types, and validated these results in all 17 cultures.

RESULTS

High variation in cell viability was observed after treatment with BMP4. In three cultures with highest sensitivity for BMP4, a substantial new cell subpopulation formed. These cells displayed decreased cell proliferation and increased apoptosis. Neuronal differentiation was reduced most in cultures with little sensitivity for BMP4. OLIG1/2 levels were found predictive for high sensitivity to BMP4. Activation of ribosomal translation (RPL27A, RPS27) was up-regulated within one day in cultures that were very sensitive to BMP4.

CONCLUSION

The changes in composition of patient-derived GBM cultures obtained after treatment with BMP4 correlate with treatment efficacy. OLIG1/2 expression can predict this efficacy, and upregulation of RPL27A and RPS27 are useful early-response markers.

Morphological Phenotyping of Organotropic Brain- and Bone-Seeking Triple Negative Metastatic Breast Tumor Cells

Triple negative breast cancer (TNBC) follows a non-random pattern of metastasis to the bone and brain tissue. Prior work has found that brain-seeking breast tumor cells display altered proteomic profiles, leading to alterations in pathways related to cell signaling, cell cycle, metabolism, and extracellular matrix remodeling. Given the unique microenvironmental characteristics of brain and bone tissue, we hypothesized that brain- or bone-seeking TNBC cells may have altered morphologic or migratory phenotypes from each other, or from the parental TNBC cells, as a function of the biochemical or mechanical microenvironment. In this study, we utilized TNBC cells (MDA-MB-231) that were conditioned to metastasize solely to brain (MDA-BR) or bone (MDA-BO) tissue. We quantified characteristics such as cell morphology, migration, and stiffness in response to cues that partially mimic their final metastatic niche. We have shown that MDA-BO cells have a distinct protrusive morphology not found in MDA-P or MDA-BR. Further, MDA-BO cells migrate over a larger area when on a collagen I (abundant in bone tissue) substrate when compared to fibronectin (abundant in brain tissue). However, migration in highly confined environments was similar across the cell types. Modest differences were found in the stiffness of MDA-BR and MDA-BO cells plated on collagen I vs. fibronectin-coated surfaces. Lastly, MDA-BO cells were found to have larger focal adhesion area and density in comparison with the other two cell types. These results initiate a quantitative profile of mechanobiological phenotypes in TNBC, with future impacts aiming to help predict metastatic propensities to organ-specific sites in a clinical setting.

Genetic diversity through social heterosis can increase virulence in RNA viral infections and cancer progression

In viral infections and cancer tumours, negative health outcomes often correlate with increasing genetic diversity. Possible evolutionary processes for such relationships include mutant lineages escaping host control or diversity, per se, creating too many immune system targets. Another possibility is social heterosis where mutations and replicative errors create clonal lineages varying in intrinsic capability for successful dispersal; improved environmental buffering; resource extraction or effective defence against immune systems. Rather than these capabilities existing in one genome, social heterosis proposes complementary synergies occur across lineages in close proximity. Diverse groups overcome host defences as interacting 'social genomes' with group genetic tool kits exceeding limited individual plasticity. To assess the possibility of social heterosis in viral infections and cancer progression, we conducted extensive literature searches for examples consistent with general and specific predictions from the social heterosis hypothesis. Numerous studies found supportive patterns in cancers across multiple tissues and in several families of RNA viruses. In viruses, social heterosis mechanisms probably result from long coevolutionary histories of competition between pathogen and host. Conversely, in cancers, social heterosis is a by-product of recent mutations. Investigating how social genomes arise and function in viral quasi-species swarms and cancer tumours may lead to new therapeutic approaches.

4-Acetylantrocamol LT3 Inhibits Glioblastoma Cell Growth and Downregulates DNA Repair Enzyme O6-Methylguanine-DNA Methyltransferase

Glioblastoma multiforme (GBM) is a deadly malignant brain tumor that is resistant to most clinical treatments. Novel therapeutic agents that are effective against GBM are required. Antrodia cinnamomea has shown antiproliferative effects in GBM cells. However, the exact mechanisms and bioactive components remain unclear. Thus, the present study aimed to investigate the effect and mechanism of 4-acetylantrocamol LT3 (4AALT3), a new ubiquinone from Antrodia cinnamomeamycelium, in vitro. U87 and U251 cell lines were treated with the indicated concentration of 4AALT3. Cell viability, cell colony-forming ability, migration, and the expression of proteins in well-known signaling pathways involved in the malignant properties of glioblastoma were then analyzed by CCK-8, colony formation, wound healing, and western blotting assays, respectively. We found that 4AALT3 significantly decreased cell viability, colony formation, and cell migration in both in vitro models. The epidermal growth factor receptor (EGFR), phosphatidylinositol-3-kinase (PI3K)/AKT/mammalian target of rapamycin (mTOR), Hippo/yes-associated protein (YAP), and cAMP-response element binding protein (CREB) pathways were suppressed by 4AALT3. Moreover, 4AALT3 decreased the level of DNA repair enzyme O⁶-methylguanine-DNA methyltransferase and showed a synergistic effect with temozolomide. Our findings provide the basis for exploring the beneficial effect of 4AALT3 on GBM in vivo.

Group Behavior and Emergence of Cancer Drug Resistance

Drug resistance is a major impediment in cancer. Although it is generally thought that acquired drug resistance is due to genetic mutations, emerging evidence indicates that nongenetic mechanisms also play an important role. Resistance emerges through a complex interplay of clonal groups within a heterogeneous tumor and the surrounding microenvironment. Traits such as phenotypic plasticity, intercellular communication, and adaptive stress response, act in concert to ensure survival of intermediate reversible phenotypes, until permanent, resistant clones can emerge. Understanding the role of group behavior, and the underlying nongenetic mechanisms, can lead to more efficacious treatment designs and minimize or delay emergence of resistance.

Deciphering neural heterogeneity through cell lineage tracing

Understanding how an adult brain reaches an appropriate size and cell composition from a pool of progenitors that proliferates and differentiates is a key question in Developmental Neurobiology. Not only the control of final size but also, the proper arrangement of cells of different embryonic origins is fundamental in this process. Each neural progenitor has to produce a precise number of sibling cells that establish clones, and all these clones will come together to form the functional adult nervous system. Lineage cell tracing is a complex and challenging process that aims to reconstruct the offspring that arise from a single progenitor cell. This tracing can be achieved through strategies based on genetically modified organisms, using either genetic tracers, transfected viral vectors or DNA constructs, and even single-cell sequencing. Combining different reporter proteins and the use of transgenic mice revolutionized clonal analysis more than a decade ago and now, the availability of novel genome editing tools and single-cell sequencing techniques has vastly improved the capacity of lineage tracing to decipher progenitor potential. This review brings together the strategies used to study cell lineages in the brain and the role they have played in our understanding of the functional clonal relationships among neural cells. In addition, future perspectives regarding the study of cell heterogeneity and the ontogeny of different cell lineages will also be addressed.

If Artificial In Vitro Microenvironment Can Influence Tumor Drug Resistance Network via Modulation of lncRNA Expression?-Comparative Analysis of Glioblastoma-Derived Cell Culture Models and Initial Tumors In Vivo

The tumor resistance of glioblastoma cells in vivo is thought to be enhanced by their heterogeneity and plasticity, which are extremely difficult to curb in vitro. The external microenvironment shapes the molecular profile of tumor culture models, thus influencing potential therapy response. Our study examines the expression profile of selected lncRNAs involved in tumor resistance network in three different glioblastoma-derived models commonly utilized for testing drug response in vitro. Differential expression analysis revealed significant divergence in lncRNA profile between parental tumors and tumor-derived cell cultures in vitro, including the following particles: MALAT1, CASC2, H19, TUSC7, XIST, RP11-838N2.4, DLX6-AS1, GLIDR, MIR210HG, SOX2-OT. The examined lncRNAs influence the phenomenon of tumor resistance via their downstream target genes through a variety of processes: multi-drug resistance, epithelial-mesenchymal transition, autophagy, cell proliferation and viability, and DNA repair. A comparison of in vivo and in vitro expression identified differences in the levels of potential lncRNA targets, with the highest discrepancies detected for the MDR1, LRP1, BCRP and MRP1 genes. Co-expression analyses confirmed the following interrelations: MALAT1-TYMS, MALAT1-MRP5, H19-ZEB1, CASC2-VIM, CASC2-N-CAD; they additionally suggest the possibility of MALAT1-BCRP, MALAT1-mTOR and TUSC7-PTEN interconnections in glioblastoma. Although our results clearly demonstrate that the artificial ex vivo microenvironment changes the profile of lncRNAs related to tumor resistance, it is difficult to anticipate the final phenotypic effect, since this phenomenon is a complex one that involves a network of molecular interactions underlying a variety of cellular processes.

The Organoid Era Permits the Development of New Applications to Study Glioblastoma

Glioblastoma (GB) is the most frequent and aggressive type of glioma. The lack of reliable GB models, together with its considerable clinical heterogeneity, has impaired a comprehensive investigation of the mechanisms that lead to tumorigenesis, cancer progression, and response to treatments. Recently, 3D cultures have opened the possibility to overcome these challenges and cerebral organoids are emerging as a leading-edge tool in GB research. The opportunity to easily engineer brain organoids via gene editing and to perform co-cultures with patient-derived tumor spheroids has enabled the analysis of cancer development in a context that better mimics brain tissue architecture. Moreover, the establishment of biobanks from GB patient-derived organoids represents a crucial starting point to improve precision medicine therapies. This review exemplifies relevant aspects of 3D models of glioblastoma, with a specific focus on organoids and their involvement in basic and translational research.

Detection of chemotherapy-resistant patient-derived acute lymphoblastic leukemia clones in murine xenografts using cellular barcodes

Clonal heterogeneity fuels leukemia evolution, therapeutic resistance, and relapse. Upfront detection of therapy-resistant leukemia clones at diagnosis may allow adaptation of treatment and prevention of relapse, but this is hampered by a paucity of methods to identify and trace single leukemia-propagating cells and their clonal offspring. Here, we tested methods of cellular barcoding analysis, to trace the in vivo competitive dynamics of hundreds of patient-derived leukemia clones upon chemotherapy-mediated selective pressure. We transplanted Nod/Scid/Il2Rγ^-/- (NSG) mice with barcoded patient-derived or SupB15 acute lymphoblastic leukemia (ALL) cells and assessed clonal responses to dexamethasone, methotrexate, and vincristine, longitudinally and across nine anatomic locations. We illustrate that chemotherapy reduces clonal diversity in a drug-dependent manner. At end-stage disease, methotrexate-treated patient-derived xenografts had significantly fewer clones compared with placebo-treated mice (100 ± 10 vs. 160 ± 15 clones, p = 0.0005), while clonal complexity in vincristine- and dexamethasone-treated xenografts was unaffected (115 ± 33 and 150 ± 7 clones, p = NS). Using tools developed to assess differential gene expression, we determined whether these clonal patterns resulted from random clonal drift or selection. We identified 5 clones that were reproducibly enriched in methotrexate-treated patient-derived xenografts, suggestive of pre-existent resistance. Finally, we found that chemotherapy-mediated selection resulted in a more asymmetric distribution of leukemia clones across anatomic sites. We found that cellular barcoding is a powerful method to trace the clonal dynamics of human patient-derived leukemia cells in response to chemotherapy. In the future, integration of cellular barcoding with single-cell sequencing technology may allow in-depth characterization of therapy-resistant leukemia clones and identify novel targets to prevent relapse.

Unveiling functional heterogeneity in breast cancer multicellular tumor spheroids through single-cell RNA-seq

Heterogeneity is an intrinsic characteristic of cancer. Even in isogenic tumors, cell populations exhibit differential cellular programs that overall supply malignancy and decrease treatment efficiency. In this study, we investigated the functional relationship among cell subtypes and how this interdependency can promote tumor development in a cancer cell line. To do so, we performed single-cell RNA-seq of MCF7 Multicellular Tumor Spheroids as a tumor model. Analysis of single-cell transcriptomes at two-time points of the spheroid growth, allowed us to dissect their functional relationship. As a result, three major robust cellular clusters, with a non-redundant complementary composition, were found. Meanwhile, one cluster promotes proliferation, others mainly activate mechanisms to invade other tissues and serve as a reservoir population conserved over time. Our results provide evidence to see cancer as a systemic unit that has cell populations with task stratification with the ultimate goal of preserving the hallmarks in tumors.