Tracking cells in epithelial acini by light sheet microscopy reveals proximity effects in breast cancer initiation
Abstract
Cancer clone evolution takes place within tissue ecosystem habitats. But, how exactly tumors arise from a few malignant cells within an intact epithelium is a central, yet unanswered question. This is mainly due to the inaccessibility of this process to longitudinal imaging together with a lack of systems that model the progression of a fraction of transformed cells within a tissue. Here, we developed a new methodology based on primary mouse mammary epithelial acini, where oncogenes can be switched on in single cells within an otherwise normal epithelial cell layer. We combine this stochastic breast tumor induction model with inverted light-sheet imaging to study single-cell behavior for up to four days and analyze cell fates utilizing a newly developed image-data analysis workflow. The power of this integrated approach is illustrated by us finding that small local clusters of transformed cells form tumors while isolated transformed cells do not.
eLife digest
There are now drugs to treat many types of cancer, but questions still remain around how these diseases start in the first place. Researchers think that tumor growth begins when a single cell suffers damage to certain sites in its DNA that eventually cause it to divide uncontrollably. That damaged cell, and its descendants, go on to form a lump, or tumor. The trouble with proving this theory is that it is hard to watch it happening in real time. Doctors usually only meet people with cancer when their tumors start to cause health problems. By this point, the tumors contain millions of cells.
A way to watch the very beginnings of a cancer could reveal risk factors within a tissue that foster the growth of a tumor. But first, researchers need to test their theory about how the disease begins in the first place. One way to do this is to surround a single cancer cell with healthy cells and watch what happens next. To do this, Alladin, Chaible et al. took healthy cells from the breast tissue of mice and grew them in the laboratory into mini-organs called organoids. These organoids share a lot of features with actual mouse breast tissue; they can even make milk if given the right hormones.
Once the organoids were ready, Alladin, Chaible et al then started modifying a small number of single cells inside them by switching on genes called oncogenes, which are known to drive cancer formation in humans. Using fluorescent proteins and a sheet of laser light it was possible to watch what happened to the cells over time. This revealed that, even though all the oncogene-driven single cells received the same signals, not all of them started to divide uncontrollably. In fact, a single modified cell had a low chance of forming a tumor on its own. The more oncogene-driven cells there were near to each other, the more likely they were to form tumors. Alladin, Chaible et al. think that this is because the healthy tissue interacts with the modified, oncogene-driven cells to suppress tumor formation. It is only when a larger number of modified cells group together and start to communicate with each other that they can override the inhibitory messages of the healthy tissue.
How healthy tissue stops single modified cells from forming tumors is not yet clear. But, with this new mini-organ system, researchers now have the tools to investigate. In the future, this could lead to new strategies to stop cancer before it has a chance to get started.
Introduction
Organoid cultures grown from cell-lines or primary cells have been successfully employed to study molecular mechanisms during different stages of tumorigenesis (Clevers, 2016; Havas et al., 2017; Simian and Bissell, 2017). However, they are derived from primary material that usually allows oncogenic activation in all cells of the tissue and therefore cannot reproduce the localized expansion at a defined part of the tissue that is seen in the patient situation. To study tumor initiation, models need to be established wherein only a few tumorigenic cells expand in the context of its immediate non-tumor microenvironment (Tabassum and Polyak, 2015) and the interplay of these populations can be visualized.
Long-term imaging of primary organoids has been achieved via light-sheet microscopy (Drost et al., 2015; Serra et al., 2019; Verissimo et al., 2016) and these efforts have benefited from the lower phototoxicity provided by SPIM imaging. However, regarding fast dividing cells in tumor organoids, tracking single-cell dynamics necessitates high resolution imaging which in turn limits the time frame in which organoids can be imaged without phototoxic effects (Held et al., 2018). Conversely, imaging primary organoids for longer time periods requires an offset of temporal and cellular resolution that eventually cannot allow single-cell fate tracking (Dekkers et al., 2016).
Two models have been used to understand tumor progression and heterogeneity (Plaks et al., 2015); the (i) hierarchical model refers to tumor propagating cells as cancer stem cells (CSCs) (Kreso and Dick, 2014) and the (ii) stochastic model states that every cell within a tumor is equally likely to be the cell of origin and facilitate tumor initiation and progression (Vogelstein et al., 1988; Williams et al., 2007). These models for tumor propagation can be reconciled when considering that the population of oncogene driven cells have the capacity to interconvert between differentiated- and stem-like states (Chaffer and Weinberg, 2015; Gupta et al., 2011; Plaks et al., 2015), a flexibility that is supported by the microenvironment (McGovern et al., 2009; Quail et al., 2012). These concepts, however, have not been probed in the context of tumor initiation upon stochastic oncogene activation in an established cell-layer.
Here we present a novel model of breast tumorigenesis where only single cells express oncogenes within healthy mammary acini and are able to establish localized outgrowth. We thereby overcome the above-mentioned limitations of studying tumorigenesis events in the face of tissue-wide transformation. Furthermore, we perform long-term imaging of these stochastic breast tumor acini for the first time at a temporal resolution that allows us to follow single-cell fates. We also integrate this approach with an image analysis pipeline capable of segmenting cells in their dynamic progression towards tumorigenesis, so they can be tracked individually over time.
Results
For modeling tumorigenesis in breast tissue, we use an inducible mouse model of breast cancer (TetO-MYC/TetO-Neu/MMTV-rtTA transgenic mice) (Fry et al., 2016; Moody et al., 2002; Podsypanina et al., 2008) that has been shown to recapitulate hallmarks of human breast disease (Havas et al., 2017; Jechlinger et al., 2009; Figure 1a). In this tractable transgenic mouse model, the activity of two potent oncogenes – MYC and Neu (the rodent homolog for the human ERBB2 gene used in this mouse model)– are under the control of a Tet-O promoter that can only be activated in the presence of both the rtTA (reverse tetracycline-controlled transactivator) protein and the doxycycline compound. This three-part inducible system allows for temporal control of oncogenic expression by regulating doxycycline in the medium or animal diet and spatial control by the MMTV promoter that confines expression of the rtTA protein to the cells of the mammary lineage (Bockamp et al., 2002). We modify this tissue wide tumorigenesis model (tri-transgenic (T) model; Figure 1a, upper panels) to generate a stochastic system by retaining only the oncogenic constructs, making them TetO-MYC/TetO-Neu mice (bi-transgenic (B) model; Figure 1a, lower panels). The rtTA inducer gene is then lentivirally delivered to single cells within luminal epithelial acini derived from the bi-transgenic (B) mouse glands, granting only this subset of cells the competence to express transforming oncogenes and thereby preventing tissue wide transformation. Hence, this novel (B) model should allow acini to present with two cell populations –tumorigenic and non-tumorigenic– and could be used to study the dynamics of early tumorigenesis in an otherwise normal epithelium (Figure 1a, lower right panels).

Establishment of stochastic tumorigenesis in mammary acini.
(a) Schematic representation of the mouse models and the in vitro culture methods used. Acini are grown from single cells harvested from the mammary glands of either bi-trangenic (B) or tri-trangenic (T) mice, transduced with lentiviral particles in solution and re-seeded into 3D cultures. Doxycycline is added to the medium to induce the expression of oncogenes in cells expressing rtTA. B mice have the MYC and Neu oncogene constructs in their genome. These oncogenes are activated in single cells infected with the Inducer-reporter (pLenti-rtTA-GFP) lentiviral particles in the presence of doxycycline, modeling stochastic breast tumorigenesis (right panel). T mice have the rtTA transducer construct along with the oncogenes and all cells in T acini can be induced to express oncogenes in 3D culture in the presence of doxycycline. T mice infected with Reporter (pLenti-NULL-GFP) lentiviral particles are used as infection controls (left panel). Both viral particles mark single cells in the acini with H2B-GFP. (b) Representative immunofluorescence staining images of fixed 3D gels with B acini transduced with Inducer-reporter virus or T acini transduced with Reporter virus before induction (top), 24 hr post induction and (middle) and 72 hr post induction (bottom) with doxycycline. GFP expressing transduced cells (green), MYC oncoprotein (magenta), DAPI nuclear stain (blue). Scale bar, 10 µm.
Primary mammary epithelial cells derived from transgenic mice were seeded in 3D matrigel as single cells to be grown for 48 hr to form small acini consisting of around 15 cells that arrange around a small lumen. A small number of single cells in these acini were then transduced with lentiviral particles (Figure 1a, middle panel). In the tissue-wide tumorigenesis model(T), acini were transduced with the reporter virus (pLv-pGK-H2B-GFP) that marks a subset of cells with H2B-GFP, while tissue wide rtTA expression is driven in all cells (Figure 1a, left panels). To achieve stochastic tumorigenesis, bi-transgenic (B) acini were transduced with the inducer-reporter virus (pLv-pGK-rtTA-p2A-H2B-GFP) that expresses rtTA and reporter H2B-GFP in only these single cells within the normal epithelium (Figure 1a, right panels). Then, doxycycline was supplemented in the medium to induce tumorigenic growth in rtTA expressing cells. Immunofluorescent staining of 3D matrigel cultures, for both sets of doxycycline-induced transduced acini, was used to validate transgene specific protein expression of the MYC oncogene (human) in only the transduced cells of (B) acini as opposed to all cells of (T) acini (Figure 1b, anti-human MYC specific antibody). We FACS sorted H2B-GFP marked cells from both systems respectively for qPCR analysis of transgene specific MYC and Neu mRNA expression. This way we normalized doxycycline dosage in both systems to obtain the same mRNA expression levels of oncogenes using 800 ng/ml doxycycline in the tissue wide (T) system and 600 ng/ml doxycycline in the stochastic (B) system (Figure 1—figure supplement 1).
To verify oncogene induction patterns and to understand the impact on epithelial morphology in the acini, 3D matrigel cultures of both doxycycline-induced (T) and (B) transduced acini were analyzed by immunofluorescent staining for polarity marker distribution. Lentivirus-transduced structures from the tissue wide induction system (T) presented with filled lumens that contained both infected (GFP positive) as well as non-infected oncogene-expressing cells. These structures exhibited overall disrupted epithelial polarity and some random microlumen, as reported previously (Jechlinger et al., 2009; Figure 2a). For the stochastic tumor induction system (B), we observed 2 phenotypes after 96 hr of oncogene expression; some of the (B) acini showed exclusively H2B-GFP marked cells in hyperplastic areas that also displayed double cell rim morphology and disturbed epithelial polarity (Figures 2a and 96 hours ON DOX, left panels), while the other (B) acini did not show expansion of the few H2B-GFP positive cells (Figures 2a and 96 hours ON DOX, right panels). Intrigued by these 2 distinct phenotypes, we analyzed transduced (H2B-GFP expression) (B) acini for their oncogenic MYC expression (Figure 2b). We verified that multi layered regions in acini contained mainly H2B-GFP positive cells that expectedly stained for MYC (Figure 2b, middle panels). For normal appearing regions, this analysis showed that the inability of single H2B-GFP positive, transduced cells to expand over 96 hr was not due to the lack of MYC expression (Figure 2b, lower panels).

Characterization of stochastic tumorigenesis in mammary acini.
(a) Representative immunofluorescence staining images of fixed 3D gels with T acini transduced with Reporter virus (left panels) and B acini transduced with Inducer-reporter virus (right panels), before induction and 96 hr post induction with doxycycline. Polarity markers include alpha-6-Integrin (magenta) and ZO-1 (yellow). Transduced cells are marked with GFP (green) and nucleus is counterstained with DAPI (blue). Scale bar, 20 µm.(b) Representative immunofluorescence staining images of fixed 3D gels with B acini transduced with Inducer-reporter virus before induction (top panels) and 96 hr post induction (middle and bottom panels) with doxycycline. GFP expressing transduced cells (green), MYC oncoprotein (magenta), DAPI nuclear stain (blue). Scale bar, 50 µm.
To follow up these observations in more detail over time, we bred the nuclear reporter H2B-mCherry into the T and B mice to mark all the cells in the acinus for inverted light-sheet microscopy (Luxendo InVi SPIM, Figure 3—figure supplement 1). The InVI SPIM was adjusted for non-phototoxic, long-term imaging (up to 4 days, every 10 min with 1 µm z-spacing). The T/H2B-mCherry acini transduced with reporter virus proliferated swiftly upon doxycycline addition showing expansion of both the marked and unmarked cells (Figure 3a); a sturdy tumor phenotype developed, manifested by multi-cell-layered rims and pronounced proliferation-associated-apoptosis in all acini (Figure 3—Video 1). In contrast, B/H2B-mCherry acini transduced with inducer-reporter virus, displayed phenotypic variation upon induction of oncogenes in the transduced cells. Some acini showed fast clonal expansion of oncogene-expressing cells that form multilayer clusters in the acinus rim. This proliferative phenotype seems to stem from several transduced cells in vicinity to each other at the start of time-lapse imaging (Figure 3b, upper panel; Figure 3—Video 2). Again, more sparsely infected acini, did not sustain proliferation of the oncogene-expressing cells (Figure 3b, lower panel; Figure 3—Video 3). Taken together, the phenotypes for (B) and (T) acini are consistent with the respective behavior observed with the 96 hr end point immunofluorescent analysis (Figure 2a and b), verifying these prior observations and excluding phototoxic effects due to longitudinal SPIM imaging.

Light sheet imaging of stochastic tumorigenesis in mammary acini 3D images of selected timepoints during live-cell time-lapse microscopy of induced T acini transduced with Reporter virus.
(a) or B acini transduced with Inducer-reporter virus (b). All cells in the acini express H2B-mCherry (magenta) and only cells transduced with lentiviral particles express H2B-GFP (green). Imaging was started 24 hr after oncogenic induction with doxycycline. In (b) the upper panel shows the proliferative phenotype seen with stochastic transformation, whereas the lower panel shows the non-proliferative phenotype observed in some stochastically transformed acini. (Imaging conditions: H2B-mCherry 594 nm Ex, 610 LP Em; H2B-GFP 488 nm Ex and 497–554 nm Em). Scale bar, 20 µm. (c) Schematic representation of the big-image data analysis pipeline developed to analyze the light sheet microscopy images. Images are acquired in two channels (H2B-mCherry in magenta and H2B-GFP in green) at 10 min intervals for 3–4 days. Big Data Processor Fiji plugin is used to pre-process the raw images and CATS Fiji plugin is used for generation of pixel probability maps (Figure 3—figure supplement 2). Image pixels of the H2B-GFP images are classified into background (black), nucleus center (green), nucleus boundary (blue) classes by manual training. Processed raw images along with the probability maps from the nucleus center channel (green) are exported to Imaris for 3D visualization, nuclear segmentation and single-cell tracking.
To analyze the dual-color light-sheet movies (H2B-mCherry—all cells in the acinus , H2B-GFP—transduced cells within the acinus) on a single-cell level, we developed a big image data compatible image analysis pipeline that allows efficient visualization of longitudinal big image data, nuclear segmentation and single-cell tracking in 3D (Figure 3c, Figure 3—figure supplement 3 and Figure 4—figure supplement 2). The Big Data Tools plugin (Big Data Processor (Tischer, Norlin, and Pepperkok)) for Fiji (Schindelin et al., 2012) was used to visualize the 2D image stacks recorded over time. It uses lazy loading to stream the image files, which can be up to a few terabytes in size per acinus over 3–4 days of imaging. It also allows for easy cropping, binning, chromatic shift correction between channels and file format conversion to Imaris compatible formats. Once the raw data have been pre-processed, the CATS plug-in from Fiji (Tischer and Pepperkok) was used for segmenting the cell nuclei in the green channel images. The CATS plug-in uses machine learning algorithms to predict the probabilities of all pixels in the image and classify them into pre-defined classes. Three classes — background, nucleus, boundary — were trained manually on the CATS tool and the segmented images were exported to the commercially available Imaris software (Bitplane AG, 2020, Software available at http://bitplane.com) for 3D rendering and tracking. The nucleus class channel from the CATS tool was used for surface rendering and the surfaces in all time points were tracked over the course of the movie. Cell tracking allowed us to follow the clonal evolution for each transduced cell in the acinus over 3 days. Single transduced cells within one acinus show a difference in proliferation and cell fate as indicated in representative tracks (Figure 4a).

Proximity of transformed cells in a normal epithelium enhances tumor proliferation and establishment.
(a) Single-cell tracking results for every cell in a representative B acinus transduced with the Inducer-reporter virus. Top panel shows the acinus at the beginning of the time-lapse (24 hr post induction) with each transduced cell surface rendered with Imaris. The middle panel shows the lineage trees of each individual cell over the time-lapse recording. Lineage trees of single cells are grouped into proliferative (highlighted in red, orange) and non-proliferative (highlighted in blue) cell clusters. The bottom panel shows the acinus at the end of the time-lapse (~76 hr post induction with doxycycline). Color coding of each cell maintained in all panels. Scale bar, 15 µm. (b) Schematic representation of the 9 features of stochastically transformed cells extracted at the beginning of time-lapse imaging. These features were assessed for their impact on tumor cell proliferation within B acini transduced with the Inducer-reporter virus using logistic regression. Lower right panel: Coefficients (represented as odds ratios) of the three features included in the best logistic regression model, colored horizontal bars represent the 95% confidence interval of the estimate. ** indicates p-value (of having no effect)<0.01, * indicates p-value<0.05. The vertical grey line indicates the position of no effect. (c) Representative B mammary acini stochastically transduced with the Inducer-reporter virus and induced with doxycycline. Left panels show acini 24 hr post induction. Color highlights indicate clusters of transduced cells identified from hierarchical clustering (shown in middle panels) with proliferative clusters highlighted in orange and non-proliferative clusters highlighted in blue. Right panels show the same acini ~ 72–76 hr post induction. Scale bar, 20 µm.
To better understand the parameters that positively affect a transduced cell in the stochastic tumorigenesis model to start proliferating and establishing a tumor within a normal epithelium, we extracted the center of mass coordinates of all cells in the acinus at the start of the imaging and represented them as dots in a 3D space: depicted as magenta for normal cells and green for transduced cells in Figure 4b (upper right panel). Equipped with the tumor formation success of each stochastically transduced acinus, we hypothesized that tumors originate from groups of independently transduced oncogene-expressing cells in closer vicinity within the acinus. To explore this in detail, we defined clusters of cells in the stochastic model as a group of oncogene-expressing cells that are closer to each other than to other oncogene-expressing cells of the same acinus. Notably, these clusters contain both oncogene-expressing cells as well as normal cells of the acinus. The cluster volume was described as the sphere located at the center of mass of a cluster with a diameter equal to the distance between the two farthest oncogene-expressing cells of the cluster. Next, we extracted 9 features that represent reasonable factors that could influence tumor establishment in the epithelium such as size of the acinus (represented by number of cells in the acinus), density of cells in the acinus, infection efficiency (represented by number of infected cells in the acinus), number of cells in a defined cluster volume, number of independently transduced cells in a cluster volume, distances between all cells in a cluster volume, distances between transduced cells in a cluster volume, fraction of transduced cells in the cluster volume and number of cell-cell contacts between transduced cells in a cluster volume (Figure 4b, left panel).
We analyzed 20 acini and all the transduced cells in these acini (n = 150) at the start of the imaging to ascertain the effect of the immediate non-transformed cell microenvironment on tumor cell proliferation. Since clusters of oncogene-expressing cells could be tracked over time, each cluster was associated with either a tumor outcome, or not. To identify which features were linked to this outcome, we chose a model selection approach, based on information theory. We fitted a logistic regression model for all possible linear combinations of features and selected the best model based on the Akaike information criterion (with correction for small sample sizes) (Calcagno and Mazancourt, 2010). In this model only one feature, the ‘number of transduced cells in the cluster’, can predict whether a cluster can form a tumor. Each additional transduced cell in a cluster increases the odds of this cluster forming a tumor by 9 (Figure 4b, lower right panel). Since other models within a few information criterion units of the best one can be considered equally good at predicting tumor formation, we also looked at features included in models within 3 information criterion units from the best one. Relative importance of features can be estimated by the sum of the information criterion values of the models in which each feature appears (Buckland et al., 1997). Computing relative importance across the selected models also revealed ‘number of transduced cells in the cluster’ as the most important feature, followed by ‘number of cells in a cluster’ and ‘number of contacts between transduced cells in a cluster’. However, the only feature with significant contribution to tumor formation is ‘number of transduced cells in the cluster’. It is important to note that this result does not rule out contribution from other factors which more targeted experiments may be able to reveal.
We also analyzed the effect of whole acinar features at the beginning and end of the imaging period to try and detect acinus-specific advantages, like acinus size, cell density and cell proliferative rate. Acini with transduced cell clusters that displayed a tumor outcome were compared to acini where the cell clusters did not result in tumor outcomes (Figure 4—figure supplement 1a). Acini with tumors only differed significantly from those without tumor on the number and proliferation rate of transduced cells in line with the observation that clusters of transduced cells are the primary indicators of tumor outcome. To account for this and possibly other acinus-specific effects, we also evaluated a logistic regression model including the three most important features identified above plus a random effect for acinus of origin and interaction terms for the starting number of normal cells in the acinus. In this model again, the number of oncogene-expressing cells in the cluster appeared as the main contributing factor to tumorigenesis with an odds ratio of 6.7 (Figure 4—figure supplement 1b).
That tumors are mainly initiated by clusters of transduced cells is further illustrated by representative acini shown in Figure 4c; proliferating cells that establish tumor foci at the end of the long-term imaging cluster together at the start of the time lapse imaging and the non-proliferative cells are more sparsely located within the acinus. In line, we found examples of similar sized acini bearing the same number of transduced cells that showed differential cluster expansion depending on the ‘numbers of transduced cells within a cluster’ (Figure 4—figure supplement 2; comparison a to b). Consequently, and again depending on the number of cells within the specific cluster, we observed clusters behaving differently within the same organoid; the organoids shown in Figure 4a, Figure 4c (top and middle panels) and Figure 4—figure supplement 2a, display clusters of single cells that do not expand over time, whereas clusters containing several cells show massive expansion (beyond normal expected cells division rates due to growth of the acinus).
Discussion
Our results indicate that a proximity-controlled interaction or signaling network between different transformed cells might guide tumor outgrowth in a normal epithelium (Greaves and Maley, 2012; Reid et al., 2010). This might be due to the repressive effect that an intact polarized tissue layer exerts on single early stage cancerous cells (Lee and Vasioukhin, 2008). According to our modeling attempts the simple parameter of ‘number of transduced cells within a cluster’ is enough to impinge on overall epithelial integrity and most likely already accounts for sub-features like ‘distances between transduced cells in a cluster volume’ or ‘number of cell-cell contacts between transduced cells in a cluster volume’, which did not pan out as predictive factors on their own. These observations, however, do not exclude the importance of such traits, which should be explored in more detail for the different situations we report (proliferative versus non proliferative clusters). In addition, the increased probability for tumor outgrowth upon proximity effects of closely located transformed cells, might also be rooted in paracrine effects (e.g. via microRNAs [Kosaka et al., 2012] in the cluster, or hydrogen ion dynamics/pH [Reshkin et al., 2014]). The questions outlined above can now be assessed employing the presented 3D system and longitudinal imaging pipeline together with appropriate fluorescent-tagged reporter proteins.
Studies on loss of important polarity proteins have highlighted their function as non-canonical tumor suppressors in breast tumorigenesis (McCaffrey et al., 2012; Partanen et al., 2012); however, other reports (Macara and McCaffrey, 2013; Xue et al., 2013) cannot confirm these observations which were all obtained from tissue wide transformed model systems, leaving conclusions open for alternative studies. Similarly, it will interesting to trace the effect of prominent tumor factors in breast cancer, such as SCRIBL/MYC (Zhan et al., 2008) or the Hippo pathway (Calses et al., 2019) on overall tissue integrity and promotion of tumor growth in a none tissue wide expression setting. Clearly, to better interrogate the effect of oncogenic drivers on epithelial polarity in more detail, to understand deficiencies in cell-cell interactions towards tumorigenesis and to settle conflicting reports as mentioned above, there is a need for a more detailed analysis, employing a model system that does not show modification of all cells in the tissue.
The interaction of tumor cells with the immediate microenvironment has been subject of extensive studies with regards to immune cells (Binnewies et al., 2018) and other tumor associated cell types (Tabassum and Polyak, 2015); however, the interaction with the normal neighboring cells has not been explored in real time using an organotypic mammalian model system. Rather, questions of cell competition in heterogenous tissues have mainly been addressed in either 2D culture systems or in the Drosophila wing (Tamori and Deng, 2011). Alternatively, organ-specific mammalian cell types that self-organize in a manner similar to in vivo can now be studied in specialized 3D culture conditions in vitro (Lancaster and Knoblich, 2014). Furthermore, these cultures can be used to model disease and serve as an alternative system for drug testing that better recapitulate effects as compared to conventional 2D cell culture (Clevers, 2016; Simian and Bissell, 2017).
Our adaptation of such culture systems, to mark and at the same time enable single cells to express regulatable oncogenes, gave rise to the first stochastic breast tumor model using mammary luminal epithelial acini. Together with a newly established imaging and data processing pipeline, we have developed an integrated approach that allows us to follow cell fates and to interrogate cell-cell interactions of a tumor cell with the normal epithelium. We report that cells within a polarized, differentiated epithelial cell layer are able to establish tumorigenic outgrowth upon stochastic activation of strong oncogenic drivers, supporting the concept of interconversion between differentiated- and stem-like states (Chaffer and Weinberg, 2015; Gupta et al., 2011; Plaks et al., 2015) for tumor initiating cells. However, not all cells expressing the very powerful combination of MYC and Neu oncogenes can expand, demonstrating an inhibitory effect that needs to stem from the normal epithelial microenvironment in the presented system. The ability to distinguish marked tumor cells from the normal epithelium will now allow us to perform single-cell RNA sequencing analysis on select sorted cells. This will help delineate the signaling networks within the immediate tumor microenvironment and further clarify the molecular basis of these distinct processes.
The amenability of 3D multicellular systems to interference with small molecule inhibitors, viral shRNA vectors and genomic editing has the potential to further our understanding of the mechanisms important during tumor initiation. In the future, this will also permit drug screens on the heterogeneous populations of normal and malignant cells, as presented in the described stochastic breast tumor model. The ability to follow tumor regression and establishment of minimal residual disease within a tissue like acinus structure will help distinguish between overall cytotoxic drugs and tumor cell specific inhibitory drugs.
Taken together, we strongly believe that our integration of a true stochastic tumor model with the ability to image single-cell fates will successfully bridge the gap between genetically modified model systems and the clinical situation, helping gain novel insights on breast cancer.
Materials and methods
Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
---|---|---|---|---|
Strain, strain background (M. musculus, FVB/NJ strain) | TetO-MYC/TetO-Neu/MMTV-rtTA | Moody et al., 2002 Havas et al., 2017 | Tri-transgenic (T) strain | |
Strain, strain background (M. musculus, FVB/NJ strain) | TetO-MYC/TetO-Neu | Moody et al., 2002 Havas et al., 2017 | Bi-transgenic (B) strain | |
Strain, strain background (M. musculus, FVB/NJ strain) | R26-H2B-mCherry | Abe and Fujimori, 2013 | RIKEN, CDB0239K | |
Recombinant DNA reagent | pWPXL plasmid | Didier Trono; Addgene | RRID:Addgene_12257 | 2nd generation lentiviral transfer plasmid |
Recombinant DNA reagent | pLVPT-GDNF-rtTR-KRAB-2SM2 | Patrick Aebischer and Didier Trono; Addgene | RRID:Addgene_11647 | Template for pGK promoter |
Recombinant DNA reagent | pCMV-VSV-G | Bob Weinberg; Addgene | RRID:Addgene_8454 | Envelope protein for producing lentiviral and MuLV retroviral particles |
Recombinant DNA reagent | psPAX2 | Didier Trono; Addgene | RRID:Addgene_12260 | 2nd generation lentiviral packaging plasmid |
Tranfectedconstruct (M. musculus) | pLv-pGK-H2B-GFP | This paper | Reporter Virus | Drives expression of Reporter-GFP |
Transfected construct (M. musculus) | pLv-pGK-rtTA-p2A-H2B-GFP | This paper | Inducer-reporter Virus | Drives expression of rtTA and Reporter-GFP |
Antibody | Anti-c-MYC (D84C12) (Rabbit monoclonal) | Cell signalling Technologies | Cat. # 5605; RRID:AB_1903938 | IF dilution (1:900) |
Antibody | Anti-alpha6-integrin (clone NKI-GoH3) (mouse monoclonal) | Millipore | Cat. # MAB1378; RRID:AB_2128317 | IF dilution (1:80) |
Antibody | Anti-ZO1 (61–7300) (rabbit polyclonal) | Life Technologies | Cat. # 61–7300; RRID:AB_2533938 | IF dilution (1:500) |
Sequence-based reagent | Neu Fw | This paper | PCR primers | CGTTTTGTGGTCATCCAGAACG |
Sequence-based reagent | Neu Rv | This paper | PCR primers | CTTCAGCGTCTACCAGGTCACC |
Sequence-based reagent | MYC Fw | This paper | PCR primers | GCGACTCTGAGGAGGAACAAGA |
Sequence-based reagent | MYC Rv | This paper | PCR primers | CCAGCAGAAGGTGATCCAGACT |
Sequence-based reagent | mCherry Fw | This paper | PCR primers | GAGGCTGAAGCTGAAGGAC |
Sequence-based reagent | mCherry Rv | This paper | PCR primers | GATGGTGTAGTCCTCGTTGTG |
Sequence-based reagent | PUM11 Fw | This paper | PCR primers | AATGTGTGGCCGGATCTTGT |
Sequence-based reagent | PUM11 Rv | This paper | PCR primers | CCCACAGTGCCTTATACACCA |
Commercial assay or kit | RNA PureLink Mini Kit | Thermo Fisher | Cat. # 12183018A | RNA Extraction Kit |
Commercial assay or kit | SuperScript VILO cDNA Synthesis Kit | Thermo Fisher | Cat. # 11754050 | cDNA Synthesis Kit |
Software, algorithm | Big Data Processor (plugin) | Fiji; Tischer et al., 2019a | RRID:SCR_018484 | |
Software, algorithm | CATS (plugin) | Fiji; Tischer and Pepperkok, 2019b | RRID:SCR_018486 | |
Software, algorithm | Imaris | Imaris x64 | Bitplane; RRID:SCR_007370 | Software available at: http://bitplane.com |
Other | DAPI stain | Thermo Fisher | Cat# 62248; RRID:AB_2307445 | (1 mg/mL), (1:1000 dilution) |
Animals
The mouse strains TetO-MYC/MMTV rtTA (D'Cruz et al., 2001) and TetO-Neu/MMTV rtTA (Moody et al., 2002), that have been previously described, were bred in order to establish the tri-transgenic strain TetO-MYC/TetO-Neu/MMTV rtTA (T) or bi-transgenic strain TetO-MYC/TetO-Neu (B). Reporter H2B-mCherry was crossed into the B and T lines using a R26-H2B-mCherry line (Abe and Fujimori, 2013) (RIKEN, CDB0239K). All ten mammary glands were harvested (from virgin female mice between 8–10 weeks old), digested and singularized for establishing acinar cultures. All mice used in this study were housed according to the guidelines of the Federation of European Laboratory Animal Science Associations (FELASA).
Rational for the use of these oncogenes: ERBB2 is overexpressed in ~20% of breast cancers (Arteaga et al., 2012), MYC in 15–50% of human breast cancer (Blancato et al., 2004). The combination of MYC and ERBB2 is found in highly aggressive human breast cancer (Al-Kuraya et al., 2004):and - in fact- Neu (the rodent homolog for the human ERBB2 gene used in this mouse model) and MYC strongly accelerate tumor onset In the combined transgenic animals (average 45 days) as compared to single transgenic animals (MYC 155 days, Neu 99 days), In all cases tumors regress rapidly to non-palpable state following oncogene silencing.
Lentivirus cloning and production
Request a detailed protocolThe lentivirus design is based on pWPXL backbone, which was a gift from Didier Trono (Addgene #12257). The coding region from the original plasmid was excised using ClaI and NdeI in order to insert a new multiple cloning site (MCS). The pGK promoter was PCR amplified from pLVPT-GDNF-rtTR-KRAB-2SM2, which was a gift from Patrick Aebischer and Didier Trono (Addgene #11647) and cloned using XhoI and EcoRI restriction sites. For the plasmid pLenti-rtTA-GFP the synthetic region rtTA-p2A-H2B-GFP was cloned downstream of the pGK promoter using EcoRI and NheI sites. The plasmid pLenti-Null-GFP is derived from the pLenti-rtTA-GFP by removing the rtTA sequence, using the restriction sites EcoRI and BamHI, and retaining H2B-GFP in the coding region. For production of lentivirus particles, we seeded 1.6 × 107 HEK-293T cells (Lenti-X - Clontech Cat. # 632180) in 500 cm2 square dishes (Corning Cat. # 431110). After 24 hr, the cells were supplemented with medium containing 25 µM of chloroquine diphosphate (Sigma-Aldrich Cat. # C6628). After a 5 hr incubation, using 360 µg of polyethyleneimine (4 µg for each µg of plasmid), we transfect the cells with a mixture of endotoxin free plasmids: 20 µg pCMV-VSV-G (Addgene #8454); 30 µg psPAX2 (Addgene #12260); 40 µg transfer plasmids pLenti-rtTA-GFP or pLenti-Null-GFP. We harvested the medium after 48 hr, 72 hr and 96 hr after transfection. Concentration of the lentivirus from the collected medium was performed using an ultracentrifuge (Beckman Sw32 rotor) at 25,000 rpm for 2 hr at 4°C. The lentivirus pellet was resuspended in 1000 µl of HBBS buffer, aliquoted and stored at −80°C. The lentivirus titer was measured using FACS analysis as described by Kutner and colleagues (Kutner et al., 2009).
3D acinus cultures
Request a detailed protocolMammary glands harvested from mice (see above), were digested in order to prepare a single-cell solution. For this, the tissue was divided in four loosely capped 50 ml falcon, each supplemented with 5 ml serum-free media (DMEM/F12 supplemented with 25 mM HEPES and 1% Pen Strep(100 U/ml Penicillin; 100 µg/ml Streptomycin; ThermoFisher Cat. # 15140122)) and 750 U of Collagenase Type 3 (Worthington Biochemical Corp Cat. # LS004183), 20 µg of Liberase (Roche Cat. # 5401020001) and incubated overnight at 37°C and 5%CO2. The glands were then mechanically disrupted using a 5 ml pipette, and washed in PBS before being pelleted at 1000 rpm for 5 min. The cell pellet was resuspended in 5 ml of 0.25% Trypsin-EDTA and incubated for 45 min at 37°C and 5%CO2. The enzymatic reaction was then neutralized using 40 ml of serum supplemented medium (DMEM/F12 with 25 mM HEPES, 1% Pen Strep and 10% FBS Tetracycline Free certified (Biowest Cat. # S181T). The cells were pelleted again, resuspended in Mammary Epithelial Cell Basal Medium (PromoCell Cat. # C-21210) and seeded in collagen coated plates (Corning Cat. # 354400) overnight at 37°C and 5%CO2. This allows for epithelial cells to adhere to the surface of the plates while the other cell types float on top in the medium and can be easily removed by vacuum suction. The epithelial cells were detached from the collagen coated plates by incubating them with 0.25% Trypsin-EDTA for 5–7 min at 37°C and 5%CO2, following inactivation with serum supplemented media. The single-cell solution was pelleted, resuspended in MEBM and counted. We mixed 50,000 cells with 90 µl of Matrigel Matrix basement Membrane growth factor reduced phenol red free (Corning Cat. # 356231), and seeded this mixture into a 12 well plate (Corning Cat. # 3336) and incubated it for 30–40 min until the matrigel solidified. The gels were supplemented with 1.5 ml MEBM and allowed to grow at 37°C and 5%CO2.
For transduction, after 3 days of growth, the gels were mechanically disrupted and placed in a 15 ml falcon. Two disrupted gels were placed in one 15 ml falcon with 2 ml of MEBM supplemented with 25U of Collagenase type I and 5 µg of Liberase. Following incubation in this solution for 2 hr at 37°C and 5%CO2, when the matrigel was totally digested, the acini were washed 3 times with 15 ml of serum supplemented media and once with 15 ml of serum free media, and pelleted at 1000 rpm for 5 min. We then supplemented the acinus pellet (from two original gels) in 10 µl of MEBM and added 6 × 105 lentivirus particles to the solution. We then mixed this solution with 90 µl matrigel and plated it in 35 mm dishes (Greiner Bio-One Cat. # 627160) and placed in incubator for 30–40 min until the matrigel solidified. The gels were supplemented with 3 ml MEBM and incubated for 2 days at 37°C and 5%CO2 in order to allow for acinus recovery and lentiviral gene expression.
For induction of oncogenes in the cells of the acini, doxycycline (Sigma Cat. # D9891) was supplemented in the media. 800 ng/ml of doxycycline was used to induce T acini and 600 ng/ml was used for B acini. qPCR analysis was used to standardize the doxycycline dosage for B acini (see below).
qPCR analysis
Request a detailed protocolThe qPCR technique was performed following the MIQE guidelines, where the total RNA was isolated from the mammary gland acini using RNA PureLink Mini Kit (ThermoFisher Cat. # 12183018A) and 2.5 ug was reverse transcribed to cDNA using SuperScript VILO cDNA Synthesis Kit (ThermoFisher Cat. # 11754050). Using Primer3 software we designed specific primers for DNA intercalating fluorescent dye approach for the transgenes Neu (Forward: CGTTTTGTGGTCATCCAGAACG and Reverse: CTTCAGCGTCTACCAGGTCACC) and c- MYC (Forward: GCGACTCTGAGGAGGAACAAGA and Reverse: CCAGCAGAAGGTGATCCAGACT). As endogenous controls, mCherry (Forward: GAGGCTGAAGCTGAAGGAC and Reverse: GATGGTGTAGTCCTCGTTGTG) and Pum1 (Forward: AATGTGTGGCCGGATCTTGT and Reverse: CCCACAGTGCCTTATACACCA) were used. Primer efficiency was verified and established between 95% and 105% Each sample was analyzed in duplicate and non-template controls were used in each qPCR run. Analyses were carried out using a StepOne device (Applied Biosystems, USA). Analysis of relative gene expression data was performed according to the 2−ΔΔCq method and the results were expressed as fold change of ΔΔCq values obtained from the reference T800 acini (Figure 1—figure supplement 1).
Immunofluorescence staining
Request a detailed protocolMatrigel cultures were grown as described above and plated on Nunc Lab-Tek II (Thermo Cat. # 155382) chambers. At pre-defined timepoints, the gels were fixed using 4% PFA for 2–3 min, following 3 washes with PBS. The gels were blocked with 10% goat serum for 2 hr at room temperature, followed by incubation with primary antibodies was done overnight at 4°C. The remaining immunofluorescence staining was performed as per standard protocol for c-MYC (Cell Signaling Technologies, Cat. # D84C12, 1:900), alpha6-integrin (Millipore Cat. # MAB1378, dilution 1:80) and ZO1 (Life Technologies Cat. # 61–7300, dilution 1:500). The nuclei were counter stained with 1:1000 DAPI (ThermoFisher Cat. # 62248, 1 mg/ml, dilution 1:1000) and mounted in anti-fading mounting medium (VECTASHIELD Mounting Medium with DAPI (Vecto Cat. # H1500-10)). Please note that the c-MYC antibody (Cell Signaling Technologies, Cat. # D84C12) recognizes specifically the human protein, which is transgenically expressed and does not recognize endogenous mouse MYC protein.
Stained gels were imaged on Leica SP5 confocal microscope using 63x water lens and the LAS AF imaging software.
Light sheet microscopy
Sample holder preparation and mounting
Request a detailed protocolImaging was performed on the InVi SPIM inverted light-sheet microscope (Luxendo Light-Sheet, Bruker Corporation). Sample mounting for the InVi SPIM is suitable for 3D matrigel cultures that are used to grow and transduce mammary acini (see above). The sample holder is made of medical grade plastic (PEEK). A 25 μm thin membrane (FEP; Luxendo) with a refractive index matching that of water is glued to the upper surface of a groove in the sample holder with a biocompatible silicone glue (Silpuran 4200; Wacker), forming a trough with transparent bottom (Figure 3—figure supplement 1). Matrigel cultures were carefully cut with a scalpel into rectangular slivers and transferred onto the FEP membrane’s trough. Once the gel sliver was aligned in place, 20–30 µl of fresh matrigel drops were poured onto the gel sliver in the sample holder until there was a thin layer of liquid matrigel on top of the gel sliver. The setup was incubated for 20 min at 37°C in a 5% CO2 incubator to allow the matrigel layer on top to solidify. Once the gel was solidified, 600–800 µl of MEBM supplemented with/without doxycycline was added to the sample holder’s FEP sheet trough. Preferably, freshly mounted sample gels were allowed to settle overnight in the incubator to prevent any gel drift during imaging, when the holder is placed into the imaging chamber of the microscope. The imaging chamber acts as an incubator with environmental control and it has a reservoir for immersion medium, which is filled with water so that both objective lenses and the bottom of the sample holder are below the water surface (Figure 3—figure supplement 1).
Imaging configuration and conditions
Request a detailed protocolThe InVi SPIM is equipped with a Nikon CFI 10x/0.3NA water immersion lens for illumination and a Nikon CFI-75 25x/1.1NA water immersion lens for detection. For excitation of GFP and mCherry, 488 nm and 594 nm laser lines were used, respectively, while emission was selected using a 497–554 nm band pass filter and a 610 nm long pass filter, respectively. 3D image stacks were acquired with a light-sheet thickness of 4 µm, a final magnification of 62.5x, resulting in 104 nm pixel size. The In-Vi SPIM environmental control was set to 37°C, 5% CO2 and 95% humidity. A series of optimization experiments, involving different laser powers, exposure times and z-step sizes yielded laser powers of 13 µW for 488 nm and 36 µW for 594 nm, 100 millisecond exposure time per frame and 1 µm z-spacing between frames to be optimal for long-term imaging (96–120 hr) without photo-bleaching or photo-toxic effects on growth.
Images were recorded as 2D planes ranging from 100 to 500 in number, depending on the acinus size. Each 3D stack of planes was recorded in 2 channels - mCherry (all cells) and GFP (transduced cells). Depending on the duration of the time lapse imaging, 450–600 image stacks (equivalent to ~72–96 hr) were recorded per acinus at 10 min intervals.
Image analysis
Request a detailed protocolBig Data Processor (Tischer et al., 2019a), a Fiji plugin for lazy loading of big image data, was used to visualize the images in 2D slicing mode, crop stacks in x, y, z, and t, bin images (3 × 3×1 in x, y, z), perform chromatic shift correction between channels and convert. h5 files from the InVi SPIM into an Imaris compatible multi-resolution file format (.ims) for further analysis (Figure 3—figure supplement 3).
The oncogenic cells (H2B-GFP channel) displayed heterogeneous morphologies as well as varying intensity textures, making it difficult to segment them using conventional thresholding approaches. We thus used a trainable segmentation approach to convert the raw intensity values into pixel probability maps, using the Fiji plugin CATS (Tischer and Pepperkok, 2019b) (Context Aware Trainable Segmentation). Using the H2B-GFP channel images as input, we trained three pixel classes: background, nucleus center and nucleus boundary. For training we drew about 20(background), 120(nucleus center), 100(nucleus boundary) labels distributed across the different time-frames of the movie. After feature computation and training of a Random Forest classifier the whole dataset was processed on EMBL’s high performance computer cluster. The segmentation of one dataset -typically 100 timepoints- is distributed across few hundred jobs, each job using 32 GB RAM, 16 cores, and running for about 30 min. The nucleus center probability maps were then exported from CATS and added as an additional channel to the converted intensity data (Figure 3—figure supplement 3).
The data were then loaded into Imaris (Bitplane AG, 2020, Software available at http://bitplane.com) for 3D visualization and further processing. Using the Imaris’ Surfaces function, we segmented the nucleus center probability maps into objects. To do so, probability maps were manually thresholded, using a surface smoothening parameter of 0.3 µm; the minimum quality parameter for seed points was set to 0.1, and object splitting was applied for objects larger than 5.5 µm. Objects with volumes less than 20 µm3 were excluded. Next, all objects were tracked over time using Imaris’ Lineage tracking algorithm with a maximum distance between objects in subsequent time points limited to 10 µm and a maximum gap size between identification of the object in a particular track limited to 10 time points. Analysis of segmentation results is shown in Figure 4—figure supplement 1. Most of the errors in the object segmentation were false merges, where two cells were segmented as one. This kind of error is frequently not sustained in the previous or following time points and the maximum gap size parameter of the tracking algorithm thus frequently provides correct tracks nonetheless. The resulting lineage trees of proliferating tumor cells within the acinus were corrected manually within Imaris, for example, excluding apoptotic cells and auto-fluorescent debris.
Center of mass coordinates of each cell were measured and exported from Imaris for subsequent feature analysis (Figure 4b).
Feature analysis
Request a detailed protocolObservations suggest that tumors in acini originate from clusters of oncogene-expressing cells produced by independent transduction events. To identify these clusters, we computed the pairwise Euclidean distances between all oncogene-expressing cells in an acinus at the start of the experiment and applied hierarchical clustering with complete linkage. Clusters were identified automatically by cutting the branches of the trees using the dynamic tree cut algorithm (Langfelder et al., 2008). This defined a cluster as a group of oncogene-expressing cells that are closer to each other than to other oncogene-expressing cells of the same acinus. Note that a cluster can be composed of a single cell if this cell is comparatively isolated from other transduced cells. For each cluster we identified the following features as possibly linked to tumor formation: (1) number of cells in the acinus (2) cell density expressed as the ratio of number of cells to acinus surface area computed by assuming the acinus is a sphere with diameter equal to the distance between the two most distant cells (3) number of oncogene-expressing cells in the acinus (4) number of cells (including both oncogene-expressing and normal cells) in the cluster volume defined as the sphere centered at the center of mass of the cluster with diameter equal to the distance between the two farthest oncogene-expressing cells of the cluster (5) number of oncogene-expressing cells in the cluster (6) average pairwise distance between all cells in the cluster volume (7) average pairwise distance between oncogene-expressing cells in the cluster (8) fraction of oncogene-expressing cells in the cluster volume (9) number of contacts between oncogene-expressing cells in the cluster. Two cells are presumed in contact if they are less than the average cell diameter + 2 standard deviation apart.
Oncogene-expressing cells were tracked over time and a cluster was associated with a tumor outcome if any of its cells lead to tumor formation. Here, ‘tumor formation’ is defined as the phenotypic observation of increased proliferation rate and disrupted polarity observed in transduced cell clusters upon induction with doxycycline. Regions displaying ‘tumor formation’ within the acinus often form multi-layer cell clusters with fast-dividing cells and increased apoptosis. To identify which features were linked to this outcome, we took an information-theoretic approach to model selection. We fitted a logistic regression model for all possible linear combinations of features and selected the best model based on the Akaike information criterion (with correction for small sample sizes) (Calcagno and Mazancourt, 2010). This model included the following features: number of oncogene-expressing cells in the cluster, number of oncogene-expressing cells in the acinus, number of cells in the cluster and number of contacts between oncogene-expressing cells in the cluster. Of these features, only the number of oncogene-expressing cells in the cluster contributed significantly to tumor formation with an odds ratio of 8.96 (Figure 4b). Of note, all the models within 3 information criterion units of the best model also included the number of oncogene-expressing cells in the cluster as a significant contributor to tumorigenesis. Relative importance of different features can be estimated by the sum of the information criterion values of the models in which each feature appears (Buckland et al., 1997). Computing relative variable importance across all models also indicated this feature as the most important. To control for acinus of origin, we also evaluated a mixed effect model including the three most important fixed effect features (number of transduced cells in a cluster, number of cells in a cluster and number of contacts between transduced cells in a cluster), a random effect for acinus and interaction terms for the starting number of normal cells. In this case again, the number of oncogene-expressing cells in the cluster appeared as the main contributing factor to tumorigenesis with an odds ratio of 6.73 (Figure 4—figure supplement 2b).
Data availability
1) Entire image recordings (movies) of time-lapse panels in Figure 3a and 3b (3 Video files in total) have been provided as supplementary movie files. 2) We have uploaded the code for the Feature analysis of the nine acinar features described in Figure 4, as source code file "Feature_Analysis.Rmd". Refer to Supplement file 1 and Online Materials and Methods section for analysis summary. 3) We have uploaded the html file describing the source code as Supplementary file 1. 4) Three. xlsx files with 20 sheets each, one sheet for each acinus analyzed are provided as Supplementary files 2, 3, and 4. These contain the x,y,z coordinates for each cell in the respective acinus at the beginning of the SPIM recording (Supplementary file 2) and at the end (Supplementary file 3). Supplementary File 4 contains the "label" for each transduced cell (corresponding to the labels in Supplement File 2) for the acini at the beginning of the SPIM recording. These. xlsx files were input into the source code to carry out the acinus feature analysis described in Figure 4b and Figure 4 - figure supplement 1. 4) We have deposited the original imaging data for all acini recorded and analyzed (20 mammary acini) at the BioStudies archive at EMBL-EBI (https://www.ebi.ac.uk/biostudies/studies/S-BIAD13). A total of 390-450. h5 image files recorder from 2 channels on the microscope are uploaded for each acini (10 minute time intervals). Raw image data from the microscope was cropped to remove empty pixels, binned in x,y (3,3) and converted to 8-bit images using Big Data Processor Fiji Plug in (http://doi.org/10.5281/zenodo.2574702). This data repository also contains video files generated via Imaris for each acinus, showing fluorescence SPIM miscropscopy data (pre-processed raw files available in respective folders) in 2-color 3D projections (mcherry- magenta; GFP- green) for observing visual phenotypes.
-
ZenodoBigDataProcessor: Fiji plugin for big image data inspection and processing.https://doi.org/10.5281/zenodo.2574702
-
BioStudiesID S-BIAD13. Tracking the cells of tumor origin in breast organoids by light sheet microscopy - SPIM movie data.
References
-
Reporter mouse lines for fluorescence imagingDevelopment, Growth & Differentiation 55:390–405.https://doi.org/10.1111/dgd.12062
-
Treatment of HER2-positive breast Cancer: current status and future perspectivesNature Reviews Clinical Oncology 9:16–32.https://doi.org/10.1038/nrclinonc.2011.177
-
Of mice and models: improved animal models for biomedical researchPhysiological Genomics 11:115–132.https://doi.org/10.1152/physiolgenomics.00067.2002
-
Model selection: an integral part of inferenceBiometrics 53:603–618.https://doi.org/10.2307/2533961
-
Glmulti : an R package for easy automated model selection with (Generalized) Linear modelsJournal of Statistical Software 34:i12.https://doi.org/10.18637/jss.v034.i12
-
How does multistep tumorigenesis really proceed?Cancer Discovery 5:22–24.https://doi.org/10.1158/2159-8290.CD-14-0788
-
Characterizing responses to CFTR-modulating drugs using rectal organoids derived from subjects with cystic fibrosisScience Translational Medicine 8:ra344–ra384.https://doi.org/10.1126/scitranslmed.aad8278
-
Clinical applications of mouse models for breast Cancer engaging HER2/neuIntegrative Cancer Science and Therapeutics 3:593–603.https://doi.org/10.15761/ICST.1000210
-
Metabolic shifts in residual breast Cancer drive tumor recurrenceJournal of Clinical Investigation 127:2091–2105.https://doi.org/10.1172/JCI89914
-
Competitive interactions of Cancer cells and normal cells via secretory microRNAsJournal of Biological Chemistry 287:1397–1405.https://doi.org/10.1074/jbc.M111.288662
-
Evolution of the Cancer stem cell modelCell Stem Cell 14:275–291.https://doi.org/10.1016/j.stem.2014.02.006
-
Cell polarity and cancer--cell and tissue polarity as a non-canonical tumor suppressorJournal of Cell Science 121:1141–1150.https://doi.org/10.1242/jcs.016634
-
Cell polarity in morphogenesis and metastasisPhilosophical Transactions of the Royal Society B: Biological Sciences 368:20130012.https://doi.org/10.1098/rstb.2013.0012
-
Microenvironmental regulation of Cancer stem cell phenotypesCurrent Stem Cell Research & Therapy 7:197–216.https://doi.org/10.2174/157488812799859838
-
Barrett's oesophagus and oesophageal adenocarcinoma: time for a new synthesisNature Reviews Cancer 10:87–101.https://doi.org/10.1038/nrc2773
-
Role of pHi, and proton transporters in oncogene-driven neoplastic transformationPhilosophical Transactions of the Royal Society B: Biological Sciences 369:20130100.https://doi.org/10.1098/rstb.2013.0100
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Organoids: a historical perspective of thinking in three dimensionsJournal of Cell Biology 216:31–40.https://doi.org/10.1083/jcb.201610056
-
Cell competition and its implications for development and CancerJournal of Genetics and Genomics 38:483–495.https://doi.org/10.1016/j.jgg.2011.09.006
-
Genetic alterations during colorectal-tumor developmentNew England Journal of Medicine 319:525–532.https://doi.org/10.1056/NEJM198809013190901
-
Loss of Par3 promotes breast Cancer metastasis by compromising cell-cell cohesionNature Cell Biology 15:189–200.https://doi.org/10.1038/ncb2663
Decision letter
-
Lynne-Marie PostovitReviewing Editor; University of Alberta, Canada
-
Richard M WhiteSenior Editor; Memorial Sloan Kettering Cancer Center, United States
-
Anne RiosReviewer
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Acceptance summary:
This paper demonstrates a beautiful methodology for tracking tumour formation within 3D acinar structures, with single cell resolution. It also suggests that tumour formation may depend upon local interactions between transformed cells. The method described may be used to better understand the underpinnings of stochastic tumor evolution.
Decision letter after peer review:
Thank you for submitting your article "Tracking cells in a novel organoid model by light sheet microscopy reveals proximity effects in breast cancer initiation" for consideration by eLife. Your article has been reviewed by three peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Richard White as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Anne Rios (Reviewer #3).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
Summary:
This manuscript by Jechlinger and colleagues combines a novel transgenic organoid model of breast hyperplasia with light-sheet microscopy in order to demonstrate that hyperplastic outgrowths occur most readily when transformed cells are adjacent to each other. Specifically, the authors compared early progression (hyperplastic growth and depolarization) in organoids wherein all versus only a subset of epithelial cells harbor the Neu and MYC driver mutations. Aided by time lapsed light-sheet microscopy and computational analyses, they surmised that the number of transduced cells in a cluster dictates progression, suggesting that special cues are likely required for cancer progression to occur. This work was reviewed by 2 experts and myself. One of the reviewers was notably an expert in the imaging and characterization of organoids, and thus able to comment on the technical aspects of this study, which given its focus, are the most important.
Overall, the reviewers felt that the paper describes an interesting tool that is accompanied by beautiful images and a fairly compelling analysis. However, in order to better justify the biological conclusions, additional experiments would be needed. These would include an analysis of cells which outgrow (versus those that do not). Such studies may include single cell RNA sequencing to determine how cells in clusters differ from those which are more isolated, as well as a more thorough assessment of polarity and consideration of microenvironmental cues such the concentration of myo-epithelial signals. Hence, the main issues raised by the reviewers related to the claims that a distinct number of initiator cells can explain tumorigenic clone expansion and that a proximity-controlled interaction or signaling network between different transformed cells is required for tumor outgrowth in a normal epithelium. In order to substantiate these claims, a number of experiments would have to be done. I would suggest, rather, that the claims are tempered and that the suggested technical revisions are completed.
Essential revisions:
1) Heterogeneities associated with the growth of organoids in general should be considered with additional controls: The authors cannot claim that you need interaction/cooperation between different transformed cells when the “organoid” with few activated cells clearly show an overall growth slower or even non-existent compare to the one with expansion of tumorigenic clones (organoid size, cell compactness). Differences of growth and growth advantage for some “organoids” while culturing is quite common. The authors should think of a comparison of growth dynamics of untransduced tissue and provide the initial volume in both conditions as well as assess the impact of proliferation differences between different organoids.
2) The idea that the initial cluster size dictates tumor formation should be better supported: If indeed tumor formation is dependent on the initial cluster size, this effect should also be detectable at single organoid level: A single organoid should have clusters formed by several cells that will give rise to a tumor and clusters formed by only 1 cell that will not give rise to a tumor. The authors should define such examples, and fit a model with random effect, where each organoid will be considered as random effect. This could as well exclude the possibility that the results are due to the fact that organoids have a differential growth capacity.
3) The AI approach should be strengthened and/or better articulated: In the logistic regression model used to assess which feature is responsible for the tumor formation, authors should include interaction effects between different features, which is not considered by the linear model used. As indeed the size of the organoids could influence the number of cells per cluster, etc.
4) The conclusions must be either softened or expanded upon (with experiments such as those outlined in the summary): Considering that the conclusion of the authors would be correct and that the tumor formation is solely dependent on the cluster size, please provide an explanation why is it independent of the interactions between cells of the same cluster. It would be logical to assume that cells have to be in close contact to exert an effect on each other.
5) The semantics regarding "organoids" should be addressed: The authors call “organoids”, the in vitro 3D culture they use. Not all 3D spheroids culture can be called “organoids”, it is essential first to examine if the 3D structures grown in culture resemble the original breast tissue (e.g. cell composition, geometry). Why not use the well-defined protocol developed for mouse organoid culturing where both main mammary epithelial cell types can be modelled (luminal and myoepithelial) (Jamieson, Dekkers et al., 2016)? The authors should either prove that their cultures can be called organoids or use the right protocol. It is of great importance as they are claiming they use a stochastic model. Looking at shape and polarity, I presume their organoids are luminal-enriched.
6) Studies regarding polarity should be strengthened: 3D z-stack and rendering should have been performed to make conclusive interpretation in Figure 2A about disrupted polarity and random lumen. 2D often lead to a wrong interpretation of architectural phenotypes. Moreover, the authors should expand the polarity analyses by measuring a number of polarity markers (ZO-1, Par3).
7) The authors should refine terminology related to localized transformation and/or perform the experiments that would allow them to suggest that transformation is in fact localized: The authors also claimed a model of localized transformation while what they did is random and low rate transformation using a doxycycline inducible model with lentiviral delivery of rtTA. It is a confusing claim as localized transformation implies you control the location. This can be performed with optogenetic technology. Alternatively, the terminology used should be refined.
8) Please ensure the supplementary figures are uploaded.
9) The authors should not suggest that a big-data analysis was conducted and/or should apply more robust analytical techniques: Authors claim that they do big-data analysis, however only 20 organoids are segmented and tracked. The authors should either provide a more robust analysis or change the semantics that are used in the paper.
Reviewer #1:
This manuscript by Jechlinger and colleagues combines a novel transgenic organoid model of breast hyperplasia with light-sheet microscopy in order to demonstrate that hyperplastic outgrowths occur most readily when transformed cells are adjacent to each other. Specifically, the authors compared early progression (hyperplastic growth and depolarization) in organoids wherein all versus only a subset of epithelial cells harbor the Neu and MYC driver mutations. Aided by time lapsed light-sheet microscopy and computational analyses, they determined that the number of transduced cells in a cluster dictates progression, suggesting that special cues are likely required for cancer progression to occur. Overall, this is an interesting tool, accompanied by beautiful images and a fairly compelling analysis. However, in order to solidify what is still mostly a descriptive finding, an analysis of cells which outgrow (versus those that do not) would be needed.
As indicated by the authors, single cell RNA sequencing should be done to determine how cells in clusters differ from those which are more isolated. Alternatively, key factors may be analyzed perhaps with IHC, to determine if there are tangible differences.
Reviewer #2:
This is an interesting paper using real time imaging and a retroviral based strategy to induce both the MYC/activated Neu in single cell. The core observation is that clusters of oncogenes expressing mammary epithelial cells go on to form tumors whereas single cell isolates do not. The authors conclude that adjacent normal cells (untransduced cells) negatively inhibit tumor outgrowth. In general, the data presented support the author's conclusions. However, the authors should address the following issues before publication.
1) Although the data indicate that single cell expressing both oncogenes have normal polarity, using a number of polarity markers (ZO-1, Par3) can the authors exclude subtle alterations in polarity program.
2) What impact would co-cultivation of myoepithelial cells have on the observed tumor dynamics.
Reviewer #3:
Alladin et al. proposed a stochastic breast tumor doxycycline inducible system using mouse breast 3D cultures to establish lower cell transformation rate. They performed lightsheet time-lapse imaging together with tracking analysis to examine clone expansion overtime after transformation. The authors claimed that distinct number of initiator cells can explain tumorigenic clone expansion. In addition, they suggest that a proximity-controlled interaction or signaling network between different transformed cells is required for tumor outgrowth in a normal epithelium. I have serious concerns about this claim and even if the claim is correct, the paper does not provide any explanation for such a phenomenon. Please see my comments below, I hope it may help improve your manuscript.
– The authors cannot claim that you need interaction/cooperation between different transformed cells when the “organoid” with few activated cells clearly show an overall growth slower or even inexistent compare to the one with expansion of tumorigenic clones (organoid size, cell compactness). Differences of growth and growth advantage for some “organoids” while culturing is quite common. The authors should think of a comparison of growth dynamics of untransduced tissue and provide the initial volume in both conditions as well as assess the impact of proliferation differences between different organoids.
– If indeed tumor formation is dependent on the initial cluster size, this effect should also be detectable at single organoid level: A single organoid should have clusters formed by several cells that will give rise to a tumor and clusters formed by only 1 cell that will not give rise to a tumor. The authors should define such examples, and fit a model with random effect, where each organoid will be considered as random effect. This could as well exclude the possibility that the results are due to the fact that organoids have a differential growth capacity.
– In the logistic regression model used to assess which feature is responsible for the tumor formation, authors should include interaction effects between different features, which is not considered by the linear model used. As indeed the size of the organoids could influence the number of cells per cluster, etc.
– Considering that the conclusion of the authors would be correct and that the tumor formation is solely dependent on the cluster size, please provide an explanation why is it independent of the interactions between cells of the same cluster. It would be logical to assume that cells have to be in close contact to exert an effect on each other.
– The author call “organoids”, the in vitro 3D culture they use. Not all 3D spheroids culture can be called “organoids”, it is essential first to examine if the 3D structures grown in culture resemble the original breast tissue (e.g. cell composition, geometry). Why not use the well-defined protocol developed for mouse organoid culturing where both main mammary epithelial cell types can be modelled (luminal and myoepithelial) (Jamieson, Dekkers et al., 2016)? The authors should either prove that their cultures can be called organoids or use the right protocol. It is of great importance as they are claiming they use a stochastic model. Looking at shape and polarity, I presume their organoids are luminal-enriched.
– 3D z-stack and rendering should have been performed to make conclusive interpretation in Figure 2A about disrupted polarity and random lumen. 2D often lead to a wrong interpretation of architectural phenotypes.
– The authors also claimed a model of localized transformation while what they did is random and low rate transformation using a doxycycline inducible model with lentiviral delivery of rtTA. It is rather confusing claim as localized transformation implies you control the location. This can be performed with optogenetic technology.
– Important to mention: No supplementary figures were uploaded.
– Authors claim that they do big-data analysis, however only 20 organoids are segmented and tracked, this is misleading. The authors should either provide a more robust analysis or change the semantics that are used in the paper. Also, few terabytes per organoid imaged…I may advise to consider other imaging technologies such as confocal that will give better images for less data generated.
https://doi.org/10.7554/eLife.54066.sa1Author response
Essential revisions:
1) Heterogeneities associated with the growth of organoids in general should be considered with additional controls: The authors cannot claim that you need interaction/cooperation between different transformed cells when the “organoid” with few activated cells clearly show an overall growth slower or even non-existent compare to the one with expansion of tumorigenic clones (organoid size, cell compactness). Differences of growth and growth advantage for some “organoids” while culturing is quite common. The authors should think of a comparison of growth dynamics of untransduced tissue and provide the initial volume in both conditions as well as assess the impact of proliferation differences between different organoids.
We would like to thank the reviewers for bringing up this crucial point. We have addressed this question via an in-depth analysis of single organoid features, inclusion of 2 additional supplementary figures and new text in the manuscript:
a) Analysis of numbers and spatial coordinates of all (transduced and normal) cells of each acinus at the start and the end of imaging allowed us to “analyze the effect of whole acinar features at the beginning and end of the imaging period to try and detect acinus-specific advantages, like acinus size, cell density and cell proliferation rate. Acini with transduced cell clusters that displayed a tumor outcome were compared to acini where the cell clusters did not result in tumor outcomes (Figure 4—figure supplement 1A). Acini with tumors only differed significantly from those without tumor on the number and proliferation rate of transduced cells in line with the observation that clusters of transduced cells are the primary indicators of tumor outcome.”
b) We illustrate this observation further by inclusion of an example of “similar sized acini bearing the same number of transduced cells that showed differential cluster expansion depending on the numbers of transduced cells within a cluster” (Figure 4—figure supplement 2; comparison A to B).
c) To account for acinus specific effects, we “also evaluated a logistic regression model including the three most important features identified (above) plus a random effect for acinus of origin and interaction terms for the starting number of normal cells in the acinus. In this model again, the number of oncogene-expressing cells in the cluster appeared as the main contributing factor to tumorigenesis with an odds ratio of 6.7 (Figure 4—figure supplement 1B).”
d) We include 3D projection videos over time (uploaded to https://www.ebi.ac.uk/biostudies/studies/S-BIAD13?key=d65c53a7-253d-4e4e-82df-d27be3a836f) for all 20 analyzed acini (150 transduced cells in total) to give the reader a transparent representation of the data and the possibility to examine valid concerns like this one, with ease.
2) The idea that the initial cluster size dictates tumor formation should be better supported: If indeed tumor formation is dependent on the initial cluster size, this effect should also be detectable at single organoid level: A single organoid should have clusters formed by several cells that will give rise to a tumor and clusters formed by only 1 cell that will not give rise to a tumor. The authors should define such examples, and fit a model with random effect, where each organoid will be considered as random effect. This could as well exclude the possibility that the results are due to the fact that organoids have a differential growth capacity.
We entirely agree with the reviewer’s notion that a single acinus should give rise to different phenotypes (proliferative versus non-proliferative), depending on the respective cluster/microenvironment.
a) In fact, we present one example acinus (Figure 4A) that bears expanding clusters originating from several transduced cells in vicinity (red and orange circled) as well as non-expanding clusters (blue circles; single cells or duplets within the cluster). Cell linages are shown over the imaging period.
b) We also present 2 further cases in Figure 4C (upper panels and middle panels).
c) In addition, as suggested by the reviewers, we included a more in-depth analysis of two acini in Figure 4—figure supplement 2A, showing cell lineages over time.
3) The AI approach should be strengthened and/or better articulated: In the logistic regression model used to assess which feature is responsible for the tumor formation, authors should include interaction effects between different features, which is not considered by the linear model used. As indeed the size of the organoids could influence the number of cells per cluster, etc.
We strengthened the computational modelling approach and added onto the initially presented “logistic regression model for all possible linear combinations of features”, while also tuning down the interpretation of the results:
a) “Since other models within a few information criterion units of the best one can be considered equally good at predicting tumor formation, we also looked at features included in models within 3 information criterion units from the best one. […]. It is important to note that this result does not rule out contribution from other factors which more targeted experiments may be able to reveal.”
b) We included a mixed effect model for “a random effect for acinus of origin and interaction terms for the starting number of normal cells in the acinus.” (see rebuttal point 1c). Again, here we obtained "number of transduced cells in the cluster" as the most prominent determinant.
4) The conclusions must be either softened or expanded upon (with experiments such as those outlined in the summary): Considering that the conclusion of the authors would be correct and that the tumor formation is solely dependent on the cluster size, please provide an explanation why is it independent of the interactions between cells of the same cluster. It would be logical to assume that cells have to be in close contact to exert an effect on each other.
We do not intend to claim that tumor formation is solely dependent on the number of transduced cells in a cluster, rather that this factor is a potent indicator of tumorigenic outcome when compared to the features we analyzed in Figure 4B. We have endeavored to make this clearer by inserting the following statements:
a) “However, the only feature with significant contribution to tumor formation is "number of transduced cells in the cluster". It is important to note that this result does not rule out contribution from other factors which more targeted experiments may be able to reveal.”; see also rebuttal point 3b).
b) “According to our modeling attempts the simple parameter of “number of transduced cells within a cluster” is enough to impinge on overall epithelial integrity and most likely already accounts for sub-features like “distances between transduced cells in a cluster volume” or “number of cell-cell contacts between transduced cells in a cluster volume”, which did not pan out as predictive factors on their own. […] The questions outlined above can now be assessed employing the presented 3D system and longitudinal imaging pipeline together with appropriate fluorescent-tagged reporter proteins.”.
5) The semantics regarding "organoids" should be addressed: The authors call “organoids”, the in vitro 3D culture they use. Not all 3D spheroids culture can be called “organoids”, it is essential first to examine if the 3D structures grown in culture resemble the original breast tissue (e.g. cell composition, geometry). Why not use the well-defined protocol developed for mouse organoid culturing where both main mammary epithelial cell types can be modelled (luminal and myoepithelial) (Jamieson, Dekkers et al., 2016)? The authors should either prove that their cultures can be called organoids or use the right protocol. It is of great importance as they are claiming they use a stochastic model. Looking at shape and polarity, I presume their organoids are luminal-enriched.
We would like to take the opportunity to briefly describe the original primary mammary 3D system and culture conditions used in our lab (Jechlinger et al., 2009, Genes and Development; Havas et al., 2017, JCI). We have anecdotal evidence for K14/SMA positive cells in this 3D system and, hence, referred to these cultures as “organoids” in the aforementioned publications. These cultures permitted us to isolate hallmarks of tumor progression and residual disease that were successfully mapped to tissues and expression data of breast cancer patients following neoadjuvant therapy (harboring residual disease) (please refer to Havas et al., 2017 and Radic Shechter et al., 2020, bioRxiv). Building on these experiences and data we decided to stay with these conditions for correlation with a stochastic model (please see also Figure 1—figure supplement 1).
However, we are aware that new culture conditions allow for better compartmentalization of cell types within 3D organoid cultures. Since our mammary epithelial cell cultures are – in fact as the reviewer points out – predominantly luminal epithelial, we accede to rename our spheroids “luminal epithelial acini”/ simply “acini”. The terminology has been changed everywhere in the manuscript.
6) Studies regarding polarity should be strengthened: 3D z-stack and rendering should have been performed to make conclusive interpretation in Figure 2A about disrupted polarity and random lumen. 2D often lead to a wrong interpretation of architectural phenotypes. Moreover, the authors should expand the polarity analyses by measuring a number of polarity markers (ZO-1, Par3).
Despite the unfortunate work situation, we feel that we could address all essential revision points with the exception of this one. However, since we are already showing a 2D analysis on alpha6-integrin and ZO1 (Figure 2A) and considering the currently insecure process of obtaining reagents as well as restrictions on accessing the microscopy facility, we were hoping for your permission to forego this revision point. Having said this, EMBL has only recently begun to gradually re-open for pressing experiments and COVID19 related research and we will try our best to perform the immunofluorescence analysis as fast as possible, in case you deem this still necessary.
7) The authors should refine terminology related to localized transformation and/or perform the experiments that would allow them to suggest that transformation is in fact localized: The authors also claimed a model of localized transformation while what they did is random and low rate transformation using a doxycycline inducible model with lentiviral delivery of rtTA. It is a confusing claim as localized transformation implies you control the location. This can be performed with optogenetic technology. Alternatively, the terminology used should be refined.
We have removed this misleading terminology and understand that it could be interpreted wrongly.
8) Please ensure the supplementary figures are uploaded.
We have made sure that all figures (now 4 main figures and 6 supplementary figures) are present in the revised manuscript file (PDF version) and also sent as single PDFs as suggested in the guidelines.
9) The authors should not suggest that a big-data analysis was conducted and/or should apply more robust analytical techniques: Authors claim that they do big-data analysis, however only 20 organoids are segmented and tracked. The authors should either provide a more robust analysis or change the semantics that are used in the paper.
We agree that “big-data” analysis can be misleading. What we were trying to convey, is that the imaging data obtained during longitudinal SPIM imaging is fairly large in size and difficult to handle (more than a terabyte per acinus, Figure 3—figure supplement 2). We have substituted the misleading terminology with “big-image-data” to avoid confusion.
Reviewer #1:
[…]
As indicated by the authors, single cell RNA sequencing should be done to determine how cells in clusters differ from those which are more isolated. Alternatively, key factors may be analyzed perhaps with IHC, to determine if there are tangible differences.
Regarding single cell RNA sequencing, we have mentioned this possibility in the Discussion as an outlook and write above in response “(We) agree with the suggested revision points and have aimed to temper the claims as well as focused on technical revisions.”
As for IF analyses, we have to refer to our answer outlined to the Essential revisions point 6.
Reviewer #2:
[…]
2) What impact would co-cultivation of myoepithelial cells have on the observed tumor dynamics.
As we outlined above in the Essential revisions point 5, we have anecdotal evidence of myoepithelial cells in these culture conditions. Due to their sparseness, it is highly unlikely that cells are infected in the stochastic system and contribute to clonal expansion. In fact, we rely on this system for an easily accessible luminal epithelial layer for the viral transduction. The presence of myoepithelial cell enriched media or culture conditions that foster their outgrowth after the transduction step, represents a future aim to further explore tumor initiation within an otherwise healthy mammary organoid.
Reviewer #3:
[…]
– Authors claim that they do big-data analysis, however only 20 organoids are segmented and tracked, this is misleading. The authors should either provide a more robust analysis or change the semantics that are used in the paper. Also, few terabytes per organoid imaged…I may advise to consider other imaging technologies such as confocal that will give better images for less data generated.
Please refer to our answer outlined to the Essential revisions point 9.
Regarding the suggestion of employing confocal microscopy: This has been tried and was clearly resulting in phototoxic effects when imaging was carried out for even a few hours at 10-15 minute intervals (to not miss out on the directionality of cell division events which necessary to perform single cell tracing). Spinning disk confocal was less phototoxic but did not match the image clarity results obtained from light sheet microscopy (SPIM). Furthermore, resolution of images 2-3 cell-layers from the outer rim of the acinus towards the center was largely reduced with confocal technology as compared to light sheet imaging using low (non-phototoxic) laser doses.
https://doi.org/10.7554/eLife.54066.sa2Article and author information
Author details
Funding
No external funding was received for this work.
Acknowledgements
The authors want to thank Sylwia Gawrzak, Ksenija Radic, Rocio Sotillo, Robert Prevedel and Jan Ellenberg for critically reading the manuscript and Marta Garcia Montero for mouse husbandry. This study was technically supported by EMBL Advanced Light Microscopy Facility (ALMF) and the EMBL Laboratory for Animal Resources (LAR).
Ethics
Animal experimentation: Animals are treated at the European Molecular Laboratory in agreement with National and International laws and policies. All effort are made to use the minimal amount of animals as possible in accordance with Russell and Burch's (1959) principle of (3Rs) reduction and highest ethical standards. The IACUC (Institutional Animal Care and Use Committee) approved the work with these mice (approval # MJ160070).
Senior Editor
- Richard M White, Memorial Sloan Kettering Cancer Center, United States
Reviewing Editor
- Lynne-Marie Postovit, University of Alberta, Canada
Reviewer
- Anne Rios
Version history
- Received: November 29, 2019
- Accepted: June 16, 2020
- Version of Record published: July 21, 2020 (version 1)
Copyright
© 2020, Alladin et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 5,152
- Page views
-
- 611
- Downloads
-
- 22
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
Cancer stem cells (CSCs) undergo epithelial-mesenchymal transition (EMT) to drive metastatic dissemination in experimental cancer models. However, tumour cells undergoing EMT have not been observed disseminating into the tissue surrounding human tumour specimens, leaving the relevance to human cancer uncertain. We have previously identified both EpCAM and CD24 as CSC markers that, alongside the mesenchymal marker Vimentin, identify EMT CSCs in human oral cancer cell lines. This afforded the opportunity to investigate whether the combination of these three markers can identify disseminating EMT CSCs in actual human tumours. Examining disseminating tumour cells in over 12,000 imaging fields from 74 human oral tumours, we see a significant enrichment of EpCAM, CD24 and Vimentin co-stained cells disseminating beyond the tumour body in metastatic specimens. Through training an artificial neural network, these predict metastasis with high accuracy (cross-validated accuracy of 87-89%). In this study, we have observed single disseminating EMT CSCs in human oral cancer specimens, and these are highly predictive of metastatic disease.
-
- Cancer Biology
- Medicine
Esophageal cancer (EC) is a fatal digestive disease with a poor prognosis and frequent lymphatic metastases. Nevertheless, reliable biomarkers for EC diagnosis are currently unavailable. Accordingly, we have performed a comparative proteomics analysis on cancer and paracancer tissue-derived exosomes from eight pairs of EC patients using label-free quantification proteomics profiling and have analyzed the differentially expressed proteins through bioinformatics. Furthermore, nano-flow cytometry (NanoFCM) was used to validate the candidate proteins from plasma-derived exosomes in 122 EC patients. Of the 803 differentially expressed proteins discovered in cancer and paracancer tissue-derived exosomes, 686 were up-regulated and 117 were down-regulated. Intercellular adhesion molecule-1 (CD54) was identified as an up-regulated candidate for further investigation, and its high expression in cancer tissues of EC patients was validated using immunohistochemistry, real-time quantitative PCR (RT-qPCR), and western blot analyses. In addition, plasma-derived exosome NanoFCM data from 122 EC patients concurred with our proteomic analysis. The receiver operating characteristic (ROC) analysis demonstrated that the AUC, sensitivity, and specificity values for CD54 were 0.702, 66.13%, and 71.31%, respectively, for EC diagnosis. Small interference (si)RNA was employed to silence the CD54 gene in EC cells. A series of assays, including cell counting kit-8, adhesion, wound healing, and Matrigel invasion, were performed to investigate EC viability, adhesive, migratory, and invasive abilities, respectively. The results showed that CD54 promoted EC proliferation, migration, and invasion. Collectively, tissue-derived exosomal proteomics strongly demonstrates that CD54 is a promising biomarker for EC diagnosis and a key molecule for EC development.