Replication Study: Wnt activity defines colon cancer stem cells and is regulated by the microenvironment

  1. Anthony Essex
  2. Javier Pineda
  3. Grishma Acharya
  4. Hong Xin
  5. James Evans
  6. Reproducibility Project: Cancer Biology  Is a corresponding author
  1. PhenoVista Biosciences, United States
  2. Explora BioLabs Inc, United States

Abstract

As part of the Reproducibility Project: Cancer Biology we published a Registered Report (Evans et al., 2015), that described how we intended to replicate selected experiments from the paper ‘Wnt activity defines colon cancer stem cells and is regulated by the microenvironment’ (Vermeulen et al., 2010). Here, we report the results. Using three independent primary spheroidal colon cancer cultures that expressed a Wnt reporter construct we observed high Wnt activity was associated with the cell surface markers CD133, CD166, and CD29, but not CD24 and CD44, while the original study found all five markers were correlated with high Wnt activity (Figure 2F; Vermeulen et al., 2010). Clonogenicity was highest in cells with high Wnt activity and clonogenic potential of cells with low Wnt activity were increased by myofibroblast-secreted factors, including HGF. While the effects were in the same direction as the original study (Figure 6D; Vermeulen et al., 2010) whether statistical significance was reached among the different conditions varied. When tested in vivo, we did not find a difference in tumorigenicity between high and low Wnt activity, while the original study found cells with high Wnt activity were more effective in inducing tumors (Figure 7E; Vermeulen et al., 2010). Tumorigenicity, however, was increased with myofibroblast-secreted factors, which was in the same direction as the original study (Figure 7E; Vermeulen et al., 2010), but not statistically significant. Finally, we report meta-analyses for each results where possible.

https://doi.org/10.7554/eLife.45426.001

Introduction

The Reproducibility Project: Cancer Biology (RP:CB) is a collaboration between the Center for Open Science and Science Exchange that seeks to address concerns about reproducibility in scientific research by conducting replications of selected experiments from a number of high-profile papers in the field of cancer biology (Errington et al., 2014). For each of these papers, a Registered Report detailing the proposed experimental designs and protocols for the replications was peer reviewed and published prior to data collection. The present paper is a Replication Study that reports the results of the replication experiments detailed in the Registered Report (Evans et al., 2015) for a paper by Vermeulen et al. (2010) and uses a number of approaches to compare the outcomes of the original experiments and the replications.

Vermeulen et al. (2010) reported that colon cancer cell subpopulations with high Wnt activity correlated with markers of cancer stems cells (CSC) and displayed enhanced tumor initiating potential. Moreover, factors secreted from cancer-associated fibroblasts (CAFs), which are an important component of the stroma, such as HGF, were reported to play a role in the formation of the CSC niche and tumorigenicity by activating the Wnt signaling pathway (Vermeulen et al., 2010). This was demonstrated using in vitro clonogenicity and in vivo tumorigenicity assays, suggesting Wnt activity defines CSCs and is regulated by the microenvironment.

The Registered Report for the paper by Vermeulen et al. (2010) described the experiments to be replicated (Figures 2F, 6D and 7E), and summarized the current evidence for these findings (Evans et al., 2015). Since that publication additional studies have reported a relationship between Wnt activity, using a Wnt reporter like Vermeulen et al. (2010), and CSC properties in various malignancies, including colorectal, lung, gastric, and breast cancer (Jun et al., 2016; Su et al., 2015). Moreover, recent studies have also reported CSC properties from cells with high expression of Wnt target genes, such as LGR5 (Dame et al., 2018; Junttila et al., 2015; Shimokawa et al., 2017). Furthermore, recent studies have continued to examine the role of the microenvironment and cancer stemness. Niche factor requirements in colorectal tumors were found to decrease during tumorigenesis (Fujii et al., 2016; Kashfi et al., 2018). While a new modeling approach suggested stem cell functionality during colorectal tumor expansion was defined by secreted factors from CAFs rather than cell-intrinsic properties (Flanagan et al., 2018; Lenos et al., 2018).

Vermeulen et al. (2010) also reported CD133, the combination of CD29/CD24, and the combination of CD44/CD166 were correlated with high Wnt activity. CD133 has been suggested to mark CSCs in various tumor types, although the accuracy as a CSC biomarker has been highly controversial (Glumac and LeBeau, 2018). In colorectal cancer, variation in clonogenic potential with specific cell populations have been reported (LaBarge and Bissell, 2008), with CD133+ cells reported to be associated with the CSC population in two separate studies (O'Brien et al., 2007; Ricci-Vitiani et al., 2007), while Shmelkov et al. (2008) reported both CD133+ and CD133- populations were capable of forming colonospheres in vitro and were serially tumorigenic in mice. Variation has also been reported in independent studies that examined CD133 expression to define a clonogenic subfraction when examining the same cell line, HCT116 (Chen et al., 2011; Dittfeld et al., 2009). As such, other studies have reported other markers, such as CD44 and CD166, to be more robust in identifying colorectal CSCs (Dalerba et al., 2007; Ozawa et al., 2014). There is also variation of the significance of CD24 expression in colorectal cancer, while the significance of CD29 needs further investigation, although the presence of these molecules have been associated with CSC characteristics (Hatano et al., 2017; Izumi et al., 2015; Muraro et al., 2012).

The outcome measures reported in this Replication Study will be aggregated with those from the other Replication Studies to create a dataset that will be examined to provide evidence about reproducibility of cancer biology research, and to identify factors that influence reproducibility more generally.

Results and discussion

Generation and characterization of primary spheroidal cultures of colon cancer cells

To assess Wnt signaling activity in colon cancer stem cells (CSC), we transduced primary spheroidal cultures of colon cancer cells with the same Wnt reporter construct as the original study, which used a TCF/LEF-1 responsive promoter to drive expression of green fluorescent protein (GFP) (TOP-GFP; Reya et al., 2003). The experimental approach to generate TOP-GFP expressing CSC cultures was described in Protocol 1 of the Registered Report (Evans et al., 2015). We used three independent spheroidal cultures, one used in the original study (Co100) and two derived from primary human colorectal cancer tissues (CSC1 and E450). The three cultures were transduced with TOP-GFP and single-cell TOP-GFP cultures were isolated. This approach, similar to the original study, was done to exclude variation in lentiviral integration and copy number between cells (Vermeulen et al., 2008). The single-cell-derived TOP-GFP cultures displayed heterogeneity in Wnt signaling, similar to what was reported in the original study (Figure 1A,B).

Figure 1 with 1 supplement see all
Analysis of CSC marker expression in TOP-GFP cultures.

(A) Representative images of the three independent single-cell-cloned CSC cultures, lentivirally transduced with TOP-GFP. Phase contrast (top) and fluorescence microscopy (bottom) for each of the cultures indicated. Bar = 90 µm. (B) Single parameter histograms for GFP intensity for each of the TOP-GFP single-cell-cloned CSC cultures with the TOP-GFPlow (10% lowest) and TOP-GFPhigh (10% highest) populations indicated. (C) Single parameter histograms for the indicated cell surface markers for each of the indicated cultures. Gray denotes TOP-GFPlow (10% lowest) and green denotes TOP-GFPhigh (10% highest) populations. (D) Density plots for CD29/CD24 and CD44/CD166 from TOP-GFPlow (gray) and TOP-GFPhigh (green) populations of each culture. Additional details for this experiment can be found at https://osf.io/tfy28/.

https://doi.org/10.7554/eLife.45426.002

The three TOP-GFP cultures were then sorted into the highest and lowest 10% of TOP-GFP-expressing cells by fluorescence-activated cell sorting and analyzed for expression of the same cell surface markers reported in the original study. This experiment is similar to what was reported in Figure 2F of Vermeulen et al. (2010) and described in Protocol 2 in the Registered Report (Evans et al., 2015). We found the TOP-GFPhigh populations were more enriched for CD133+ or CD166+ cells compared to the TOP-GFPlow populations for each of the three cultures (Figure 1C,D). There were also more CD29+ cells in the TOP-GFPhigh populations for each of the three cultures, while the two TOP-GFP populations were mostly similar for CD24 expression. We also found the TOP-GFPlow populations from E450 and CSC1 cultures were more enriched for CD44+ cells, while both the populations displayed similar expression for Co100. The original study stated that CD133, the combination of CD29/CD24, and the combination of CD44/CD166 were correlated with the TOP-GFPhigh population (Vermeulen et al., 2010). However, since the degree that each of the markers correlated with the TOP-GFPhigh and TOP-GFPlow populations were not completely reported in the original study, it is difficult to directly compare to the results reported in this replication attempt. To summarize, for this experiment we found results that varied in direction relative to the original study.

Clonogenicity of TOP-GFP CSC cultures

The three TOP-GFP cultures were then used to assess the clonogenic potential of the cells using a limiting-dilution assay. Different TOP-GFP expressing fractions were examined to test if variation in TOP-GFP levels resulted in differential clonogenicity. Additionally, TOP-GFPlow fractions were treated with conditioned medium derived from myofibroblasts (MFCM) or hepatocyte growth factor (HGF), with and without a specific c-Met inhibitor (PHA665752), to test if myofibroblast-secreted factors increased the clonogenic potential. This experiment is similar to what was reported in Figure 6D of Vermeulen et al. (2010) and described in Protocol 3 in the Registered Report (Evans et al., 2015). We first performed a pilot assay of the dilution curve in untreated conditions and observed the clonogenic potential in the TOP-GFPhigh fractions were greater than the TOP-GFPlow fractions for each of the cultures tested (Figure 2—figure supplement 1). We then proceeded with the experiment to test all of the conditions specified in the Registered Report and reported in the original study. Similar to the pilot assay, we found the clonogenic potential of TOP-GFPhigh cells were greater than TOP-GFPlow cells for each of the three cultures (Figure 2). We also observed that the clonogenicity of TOP-GFPlow cells were increased in the presence of MCFM or HGF, which was reduced when PHA665752 was included, although to varying degrees across the different cultures. Interestingly, PHA665752 treatment on the whole population of TOP-GFP cells (TOP-GFPwhole) had varying effects on the clonogenicity among the different cultures tested. Both the Co100 and E450 cultures had decreased clonogenicity in the presence of PHA665752, while CSC1 cultures were increased. The original study reported the clonogenic potential of the TOP-GFPhigh fraction was greater compared to TOP-GFPlow cells, with the clonogenic potential of TOP-GFPlow cells enhanced with MFCM, or HGF, treatment, almost to the level of TOP-GFPhigh cells, which was blocked with PHA665752 (Vermeulen et al., 2010). PHA665752 was also reported to have no effect on the clonogenicity of TOP-GFPwhole cells (Vermeulen et al., 2010). The HGF and PHA665752 concentrations were the same between the original study and this replication attempt (25 ng/ml and 500 nM, respectively) as was the MFCM treatment that used a 1:2 dilution of MFCM diluted in CSC medium. As suggested during peer review of the Registered Report (Gilbertson, 2015), we also determined the concentration of HGF in MFCM, which was determined by enzyme-linked immunosorbent assay (ELISA) to be 0.61 ng/ml. The original study reported HGF production in MFCM was ~120 ng/ml (Figure 5E; Vermeulen et al., 2010) or approximately 200 times higher than what we observed. Other studies that measured the amount of HGF in MFCM using the same cell line (18Co) and timeline (24 hr) reported concentrations of ~0.4 ng/ml (Shao et al., 2006) and ~6 ng/ml (Woo et al., 2015). The variation of HGF production might be explained by differences in assay reagents, such as the generation of the standard curve (Jones et al., 1995) and variability in microplate surface properties (Lilyanna et al., 2018). Further, HGF production has been shown to be influenced by other soluble factors, such as transforming growth factor-β (TGF-β) and basic fibroblast growth factor (FGF2) (Neuss et al., 2004). The variation in HGF production in MFCM between the original study and this replication attempt might account for any observed differences in outcomes and should be taken into account when interpreting these results. Importantly, though, observing and reporting all outcomes are informative to establish the range of conditions under which a given phenotype can be observed (Bailoo et al., 2014).

Figure 2 with 1 supplement see all
Clonogenicity assay of TOP-GFP cultures.

A limiting-dilution assay was performed on the TOP-GFPlow, TOP-GFPhigh, or TOP-GFPwhole populations of the three indicated TOP-GFP cultures. Cells were left untreated, or treated with 25 ng/ml HGF, 1:2 dilution of MFCM, or 500 nM PHA-665752 (PHA), as indicated. The bar graphs present the clonogenic potential of each culture with error bars representing 95% confidence intervals (y-axis is log2 scale). This experiment was performed once for each culture. See Materials and methods and Registered Report (Evans et al., 2015) for details on limiting-dilution statistics and scheme. Planned contrast between TOP-GFPlow vs TOP-GFPhigh: E450 (χ2 = 39.8, uncorrected p=2.82×10−10, corrected p=1.69×10−9); CSC1 (χ2 = 4.82, uncorrected p=0.028, corrected p=0.169); Co100 (χ2 = 7.59, uncorrected p=0.0059, corrected p=0.035). Planned contrast between TOP-GFPlow vs TOP-GFPlow + HGF: E450 (χ2 = 1.49, uncorrected p=0.223, corrected p>0.99); CSC1 (χ2 = 0.337, uncorrected p=0.562, corrected p=0.99); Co100 (χ2 = 12.7, uncorrected p=3.70×10−4, corrected p=0.0022). Planned contrast between TOP-GFPlow vs TOP-GFPlow + MFCM: E450 (χ2 = 1.96, uncorrected p=0.162, corrected p=0.969); CSC1 (χ2 = 4.18, uncorrected p=0.041, corrected p=0.245); Co100 (χ2 = 28.7, uncorrected p=8.26×10−8, corrected p=4.96×10−7). Planned contrast between TOP-GFPlow + HGF vs TOP-GFPlow + HGF + PHA: E450 (χ2 = 0.376, uncorrected p=0.540, corrected p>0.99); CSC1 (χ2 = 34.0, uncorrected p=5.64×10−9, corrected p=3.39×10−8); Co100 (χ2 = 5.13, uncorrected p=0.024, corrected p=0.141). Planned contrast between TOP-GFPlow + MFCM vs TOP-GFPlow + MFCM + PHA: E450 (χ2 = 61.0, uncorrected p=5.71×10−15, corrected p=3.43×10−14); CSC1 (χ2 = 43.5, uncorrected p=4.14×10−11, corrected p=2.48×10−10); Co100 (χ2 = 17.6, uncorrected p=2.67×10−5, corrected p=1.60×10−4). Planned contrast between TOP-GFPwhole vs TOP-GFPwhole + PHA: E450 (χ2 = 68.3, uncorrected p=1.43×10−16, corrected p=8.56×10−16); CSC1 (χ2 = 72.2, uncorrected p=1.96×10−17, corrected p=1.17×10−16); Co100 (χ2 = 20.2, uncorrected p=6.91×10−6, corrected p=4.14×10−5). Additional details for this experiment can be found at https://osf.io/k9vce/.

https://doi.org/10.7554/eLife.45426.004

As outlined in the Registered Report (Evans et al., 2015), we planned to conduct six comparisons using the Bonferroni correction to adjust for multiple comparisons, making the a priori significance threshold 0.0083. We performed Extreme Limiting Dilution Analysis (ELDA) (Hu and Smyth, 2009) and tested for pairwise differences in frequency between groups (see Figure 2 figure legend). The sample sizes were determined a priori to detect the effects based on the originally reported data. We found that the test between TOP-GFPhigh and TOP-GFPlow cells was statistically significant for Co100 and E450 cultures, but not CSC1. Treatment with HGF, or MFCM, resulted in a statistically significant increase in clonogenicity in TOP-GFPlow cells from the Co100 culture, but not E450 or CSC1 cultures. The comparison of HGF treatment, with or without PHA665752, in TOP-GFPlow cells was statistically significant for the CSC1 culture, but not the Co100 or E450 cultures, while the comparison of MFCM treatment, with or without PHA665752, in TOP-GFPlow cells was statistically significant for all three cultures. Furthermore, the differences observed in TOP-GFPwhole cells with or without PHA665752 were statistically significant for all three cultures. To summarize, for this experiment we found results that were in the same direction as the original study, except for treatment of TOP-GFPwhole cells with PHA665752, and statistical significance that varied among the three cultures tested as well as the original study.

Tumorigenicity of TOP-GFP CSC culture

We also examined the frequency TOP-GFP cells form tumors when injected into nude mice. This experiment is similar to what was reported in Figure 7E of Vermeulen et al. (2010) and described in Protocol 4 in the Registered Report (Evans et al., 2015). While the original study included TOP-GFPlow cells co-injected with myofibroblasts and TOP-GFPwhole cells, this replication attempt was restricted to TOP-GFPhigh, TOP-GFPlow, and TOP-GFPlow cells co-injected with MFCM. The original study also reported results from two clones, while this replication attempt utilized a single clone. As stated in the Registered Report, we identified the clone to use as the one with the largest observed difference in clonogenicity between untreated TOP-GFPhigh and TOP-GFPlow cells, which, as described above, was the E450 culture. Different cell numbers were injected into female nude mice and blindly analyzed for tumor formation after nine weeks. We found the frequency of tumorigenicity was similar when TOP-GFPhigh cells (1 in every 3332, 95% CI [9174, 1210]) or TOP-GFPlow cells (1 in every 2744, 95% CI [7377, 1020]) were injected (Table 1), which was not a statistically significant difference (χ2 = 0.084, uncorrected p=0.772, corrected p>0.99). The addition of MFCM to TOP-GFPlow cells resulted in an increased frequency of tumorigenicity (1 in every 774, 95% CI [2268, 264]), which was not statistically significant when compared to untreated TOP-GFPlow cells (χ2 = 3.32, uncorrected p=0.069, corrected p=0.137). The original study reported for each of the two clones tested (C100.B5 and C100.G7) the TOP-GFPhigh fraction was more effective in inducing tumors (C100.B5:~1 in every 37, 95% CI [92, 15]; C100.G7:~1 in every 961, 95% CI [2498, 369]) than the TOP-GFPlow fraction (C100.B5:~1 in every 6939, 95% CI [18841, 2555]; C100.G7: frequency estimate unable to be determined) and that tumorigenicity was increased when TOP-GFPlow cells were co-injected with MFCM (C100.B5:~1 in every 310, 95% CI [789, 122]; C100.G7:~1 in every 2352, 95% CI [5236, 1056]) (Vermeulen et al., 2010). To summarize, for this experiment we found results that were in the same direction as the original study for the comparison between TOP-GFPlow with or without MFCM, but not for the comparison of TOP-GFPlow and TOP-GFPhigh, and not statistically significant where predicted.

Table 1
Tumorigenicity assay of TOP-GFP culture.

Cell numbers from the indicated populations were injected into female athymic nude mice. Cells were left untreated or treated with 1:2 dilution of MFCM for 2 hr before injection. The number of successful tumor initiations after nine weeks out of four injected mice for each condition is reported. Planned contrast between TOP-GFPlow vs TOP-GFPhigh (χ2 = 0.084, uncorrected p=0.772, corrected p>0.99). Planned contrast between TOP-GFPlow vs TOP-GFPlow + MFCM (χ2 = 3.32, uncorrected p=0.069, corrected p=0.137). Additional details for this experiment can be found at https://osf.io/j73xu/.

https://doi.org/10.7554/eLife.45426.006
LineCondition1010010005000
E450TOP-GFP Low0/40/42/43/4
TOP-GFP High0/40/43/42/4
TOP-GFP Low + MFCM0/42/42/44/4

Meta-analysis of original and replication effects

We performed a meta-analysis using a random-effects model, where possible, to combine each of the effects described above as pre-specified in the confirmatory analysis plan (Evans et al., 2015). To provide a standardized measure of the effect, a common effect size was calculated for each effect from the original and replication studies. Cohen’s ω is a standardized measure of the association between two variables, in this case the cells tested and clonogenic, or tumorigenic, frequency. The estimate of the effect size of one study, as well as the associated uncertainty (i.e. confidence interval), compared to the effect size of the other study provides another approach to compare the original and replication results (Errington et al., 2014; Valentine et al., 2011). Importantly, the width of the confidence interval for each study is a reflection of not only the confidence level (e.g. 95%), but also variability of the sample (e.g. SD) and sample size.

There were six comparisons of the in vitro clonogenicity assay, which were reported in Figure 2 of this study and Figure 6D of Vermeulen et al. (2010). The effect size point estimates of the original study for each of the effects was not within the 95% CI of the replication results, and vice versa (Figure 3A). Furthermore, the effect sizes were larger in the original study, compared to the three TOP-GFP cultures tested in this replication attempt, with the exception of treatment of TOP-GFPwhole cells with or without PHA665752. The meta-analyses were statistically significant for the comparison of untreated TOP-GFPlow and MFCM-treated TOP-GFPlow (p=0.047), but not for the other five comparisons (see Figure 3 figure legend). Additionally, for the comparison of untreated TOP-GFPlow and untreated TOP-GFPhigh as well as untreated TOP-GFPlow and HGF treated TOP-GFPhigh, the large CI of the meta-analyses along with statistically significant Cochran’s Q tests (p=0.0085 and p=0.0028, respectively) suggest heterogeneity between the original and replication studies.

Meta-analyses of each effect.

Effect size and 95% confidence interval are presented for Vermeulen et al. (2010), the results from this replication study (RP:CB), and a random effects meta-analysis of the effects. Cohen’s ω is a standardized measure of the association between the cells tested and clonogenic, or tumorigenic, frequency. The higher the value, the stronger the association, with an effect size of zero indicating there was no association. Sample sizes used in Vermeulen et al. (2010) and RP:CB are reported under the study name. (A) Comparison of clonogenic frequency between the indicated treated, or untreated, populations of TOP-GFP CSC cultures. TOP-GFPlow vs TOP-GFPhigh (meta-analysis p=0.094); TOP-GFPlow vs TOP-GFPlow + HGF (meta-analysis p=0.110); TOP-GFPlow vs TOP-GFPlow + MFCM (meta-analysis p=0.047); TOP-GFPlow + HGF vs TOP-GFPlow + HGF + PHA (meta-analysis p=0.218); TOP-GFPlow + MFCM vs TOP-GFPlow + MFCM + PHA (meta-analysis p=0.085); TOP-GFPwhole vs TOP-GFPwhole + PHA (meta-analysis p=0.498). (B) Comparison of frequency of tumorigenicity between the indicated treated, or untreated, populations of TOP-GFP CSC cultures injected into mice. TOP-GFPlow vs TOP-GFPhigh (meta-analysis p=0.330); TOP-GFPlow vs TOP-GFPlow + MFCM (meta-analysis p=0.033). Additional details for these meta-analyses can be found at https://osf.io/g4ewk/.

https://doi.org/10.7554/eLife.45426.007

There were two comparisons of the in vivo tumorigenicity assay, which were reported in Table 1 of this study and Figure 7E of Vermeulen et al. (2010). Similar to the clonogenicity assay, the effect sizes were larger in the original study compared to this replication attempt, and the point estimates of each study were not within the 95% CI of the other study (Figure 3B). The meta-analysis of the TOP-GFPlow and TOP-GFPhigh comparison was not statistically significant (p=0.330) with a large 95% CI and a statistically significant Cochran’s Q test (p=1.35×10−55) that suggests heterogeneity between the original and replication studies. The meta-analysis of untreated TOP-GFPlow and MFCM treated TOP-GFPlow was statistically significant (p=0.033), suggesting the null hypothesis that MFCM treatment does not impact tumorigenicity of TOP-GFPlow cells can be rejected; however, the large 95% CI and a statistically significant Cochran’s Q test (p=0.011) suggest heterogeneity between the original and replication studies.

This direct replication provides an opportunity to understand the present evidence of these effects. Any known differences, including reagents and protocol differences, were identified prior to conducting the experimental work and described in the Registered Report (Evans et al., 2015). However, this is limited to what was obtainable from the original paper and through communication with the original authors, which means there might be particular features of the original experimental protocol that could be critical, but unidentified. So while some aspects, such as the TOP-GFP reporter plasmid, cell surface markers, treatment conditions of cultures, Co100 culture, and mouse strain were maintained, others were unknown or not easily controlled for. These include variables such as cell line genetic drift (Hughes et al., 2007; Kleensang et al., 2016), genetic heterogeneity of mouse inbred strains (Casellas, 2011), the microbiome of recipient mice (Macpherson and McCoy, 2015), housing temperature in mouse facilities (Kokolus et al., 2013), lot variability of key reagents such as HGF and PHA-665752 (Leek et al., 2010), and similarities and differences in patient characteristics (Klevorn and Teague, 2016). Environmental differences such as husbandry staff, bedding type and source, light levels, and other intangibles, all of which, by necessity, differed between the studies also affect experimental outcomes with mice (Howard, 2002; Jensen and Ritskes-Hoitinga, 2007; Nevalainen, 2014; Sorge et al., 2014). The difference in HGF production in conditioned medium between the original study and this replication attempt, as described above, is another factor to consider. Also, differences in CSC features between the cultures used in this replication attempt, as well as the original study, is another important factor to consider. This includes the expression of the cell surface markers between the TOP-GFPlow and TOP-GFPhigh populations, particularly CD24 and CD44, which have been reported as Wnt target genes (Shulewitz et al., 2006; Wielenga et al., 1999). This could be due to differences in lentiviral integration and copy number of the TOP-GFP reporter, clonal artifacts, and genetic differences between the cancer cells the cultures were derived from as well as genetic drift during passaging of the cultures (Ben-David et al., 2018). Whether these or other factors influence the outcomes of this study is open to hypothesizing and further investigation, which is facilitated by direct replications and transparent reporting.

Materials and methods

Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional
information
Cell line (Homo sapiens)Co100doi:10.1038/ncb2048shared by
Medema lab, University of
Amsterdam
Cell line (H. sapiens, female)CSC1ProMab Biotechnologiescat# CC100103
Cell line (H. sapiens, female)E450this paper
Cell line (H. sapiens, female)18CoATCCcat# CRL-1459; RRID:CVCL_2379
Strain, strain background (Mus musculus, Athymic Nude, female)athymic nudeCharles RiverStrain code: 490; RRID:IMSR_CRL:490
Recombinant DNA reagentTOP-GFPdoi:10.1038/nature01593RRID:Addgene_14715shared by Medema lab, University of Amsterdam
Chemical compound, drugHGFSigma-Aldrichcat# H5791lot# MKBT3102V
Chemical compound, drugPHA-665752Sigma-Aldrichcat# PZ0147
OtherMatrigelCorningcat# 356230
AntibodyPE-conjugated anti-CD133Miltenyi Bioteccat# 130-098-826; clone: AC133; RRID:AB_26608821:100 dilution
AntibodyPE-conjugated anti-CD24BD Biosciencescat# 560991; clone ML5; RRID:AB_105630741:100 dilution
AntibodyAPC-conjugated anti-CD29BD Biosciencescat# 561794; clone: MAR4; RRID:AB_108981631:100 dilution
AntibodyPE-conjugated anti-CD166R and D Systemscat# FAB6561P; clone: 105902; RRID:AB_22238871:100 dilution
AntibodyAPC-conjugated anti-CD44BD Biosciencescat# 560890; clone: G44-26; RRID:AB_20339591:100 dilution
AntibodyPE-conjugated mouse IgG1 isotype controlMiltenyi Bioteccat# 130-098-106; clone: X-56; RRID:AB_26614631:100 dilution
AntibodyAPC-conjugated mouse IgG2b, κ isotype controlBD Biosciencescat# 555745; clone: 27–35; RRID:AB_3986121:100 dilution
AntibodyPE-conjugated mouse IgG2a, κ isotype controlBD Biosciencescat# 555574; clone: G155-178; RRID:AB_3959531:100 dilution
AntibodyAPC-conjugated mouse IgG1 isotype controlBD Biosciencescat# 555751; clone: MOPC-21; RRID:AB_3986131:100 dilution
Software, algorithmFACS Sortware sorterBD BiosciencesRRID:SCR_016722version 1.2.0.142
Software, algorithmHCS Studio Cell AnalysisThermo Fisher ScientificRRID:SCR_016787version 6.6.0
Software, algorithmFACSDivaBD BiosciencesRRID:SCR_016722version 6.1.3 or 8.0.1
Software, algorithmFlowJoTree Star, IncRRID:SCR_008520version 10
Software, algorithmR Project for statistical computinghttps://www.r-project.orgRRID:SCR_001905version 3.5.1

As described in the Registered Report (Evans et al., 2015), we attempted a replication of the experiments reported in Figures 2F, 6D, and 7E of Vermeulen et al. (2010). A detailed description of all protocols can be found in the Registered Report (Evans et al., 2015) and are described below with additional information not listed in the Registered Report, but needed during experimentation.

Cell culture

Request a detailed protocol

Three cultures of CSCs were isolated/obtained for this study. Co100 cells, which were used in the original study, were shared by Dr. Jan Paul Medema (University of Amsterdam). CSC1 cells were obtained commercially from primary human colorectal tumor tissue from a female Caucasian patient at the age of 65 (ProMab Biotechnologies, cat# CC100103; datasheet available at https://osf.io/det4j/). E450 cells were isolated as described in the Registered Report (Evans et al., 2015) from a freshly excised human colon adenocarcinoma tumor fragment from a female Caucasian patient at the age of 78. Of note, E450 was the only viable spheroidal culture that was successfully derived from twelve different colon tissue fragments that were attempted. This is slightly lower (8.3%) then the range of what was shared by the original authors during preparation of the Registered Report (10–20%: Evans et al., 2015) and previously reported efficiency rates (15%: Qureshi-Baig et al., 2016; 11%: Brattain et al., 1981; 33%: McBain et al., 1984), although methods to increase efficiency have since been reported (73%; Miyoshi et al., 2018). Patient phenotype (e.g. age, sex, ethnicity, diagnosis) for E450 are available at https://osf.io/ysf58/. Approval was obtained from Western Institutional Review Board (WIRB) (protocol MR #0701) and were in full compliance with good clinical practices as defined under the U.S. Food and Drug Administration (FDA) regulations, U.S. Department of Health and Human Services (HHS) regulations, and the International Conference on Harmonisation (ICH) guidelines. Shared samples and data were de-identified for this study. CSCs were maintained at 37°C in a humidified atmosphere at 5% CO2 in CSC medium (modified neurobasal A medium supplemented with 1X N2 supplement, lipid mixture-1 (1 ml/500 ml medium), 20 ng/ml fibroblast growth factor-basic, and 50 ng/ml epidermal growth factor) and passaged as described in the Registered Report (Evans et al., 2015) with additional details available at https://osf.io/dtbvp/. 18Co cells (ATCC, cat# CRL-1459, RRID:CVCL_2379) were maintained at 37°C in a humidified atmosphere at 5% CO2 in Eagle’s Minimum Essential Medium supplemented with 10% FBS, 500 U/ml penicillin, 500 U/ml streptomycin, and 1.25 µg/ml amphotericin B. Quality control data confirming the cells were free of mycoplasma contamination (MycoAlert Mycoplasma Detection kit; Lonza, cat# LT07-318) is available at https://osf.io/xzh9t/.

Lentiviral infection

Request a detailed protocol

TCF/LEF-1 responsive promoter to drive expression of green fluorescent protein (GFP) (TOP-GFP, RRID:Addgene_14715) was shared Dr. Jan Paul Medema (University of Amsterdam) with permission from Dr. Laurie Ailles (University Health Network; University of Toronto). Spheroidal cultures were transduced with lentiviral particles to express TOP-GFP which were produced by Cyagen Biosciences, Inc (Santa Clara, California) with a titer of 3.66 × 108 TU/ml as determined by quantitative PCR using a fragment in the WPRE region of the lentiviral vector amplified from genomic DNA of transduced HEK293 cells. Dissociated spheroidal cultures were each transduced for 24 hr with 20 µl concentrated lentivirus per 1 × 106 cells in 10 ml CSC medium supplemented with 8 µg/l polybrene before medium was replaced. Cells were cultured for 4 weeks before isolation of single-cell-derived cultures by fluorescence-activated cell sorting (FACS). FACS was performed on an Influx cell sorter (BD Biosciences) and analyzed with FACS Sortware sorter software (BD Biosciences, RRID:SCR_016722), version 1.2.0.142. Spheroids were dissociated as described in the Registered Report (Evans et al., 2015) and propidium iodide (PI) was added at 250 ng/ml immediately prior to sorting. Single, PI-negative, GFP-positive cells were sorted and deposited into individual wells of ultralow-adhesion 96-well plates containing 200 µl/well CSC medium. Four 96-well plates were tested for each culture with one viable single-cell clone generated from the CSC1 culture (0.26% efficiency), three clones from the E450 culture (0.78% efficiency, and eight clones from the Co100 culture (2.08% efficiency), which were near the range of what was shared by the original authors during preparation of the Registered Report (~1%: Evans et al., 2015). One clone was randomly selected from each culture for further analysis. Over a period of 13 weeks, spheroid cultures arising from single cells were gradually expanded into larger ultralow-adhesion flasks. Microscopy images of cultures were acquired with a CellInsight CX7 High-Content Screening (HCS) Platform (ThermoFisher Scientific) and HCS Studio Cell Analysis software (ThermoFisher Scientific, RRID:SCR_016787) version 6.6.0, build 8153.

Flow cytometry analysis of cell surface markers

Request a detailed protocol

Spheroid cultures were dissociated with trypsin and resuspended at a final concentration of 1 × 106 cells/ml in FACS buffer (PBS supplemented with 2% fetal bovine serum (FBS), 1X antifungal/antibiotic, and 2 mM EDTA). Cells were stained at 1:100 dilution with PE-conjugated monoclonal anti-CD133 (Miltenyi Biotec, cat# 130-098-826, clone AC133, RRID:AB_2660882), PE-conjugated monoclonal anti-CD24 (BD Biosciences, cat# 560991, clone ML5, RRID:AB_10563074) and APC-conjugated monoclonal anti-CD29 (BD Biosciences, cat# 561794, clone MAR4, RRID:AB_10898163), or PE-conjugated monoclonal anti-CD166 (R and D Systems, cat# FAB6561P, clone 105902, RRID:AB_2223887) and APC-conjugated monoclonal anti-CD44 (BD Biosciences, cat# 560890, clone G44-26, RRID:AB_2033959) and incubated at 4°C in the dark for 10 min. Cells were also stained with 1:100 dilution of control antibodies: PE-conjugated monoclonal mouse IgG1 isotype control (Miltenyi Biotec, cat# 130-098-106, clone X-56, RRID:AB_2661463), APC-conjugated monoclonal mouse IgG2b, κ isotype control (BD Biosciences, cat# 555745, clone 27–35, RRID:AB_398612), PE conjugated monoclonal mouse IgG2a, κ isotype control (BD Biosciences, cat# 555574, clone G155-178, RRID:AB_395953), or APC-conjugated monoclonal mouse IgG1 isotype control (BD Biosciences, cat# 555751, clone MOPC-21, RRID:AB_398613). Cells were washed by adding 20 times the reaction volume of FACS buffer and gently inverting tubes three times. Cells were centrifuged at 1000 RPM for 10 min, supernatant was carefully aspirated, and the cells were resuspended in 100 µl FACS buffer. Flow cytometry analysis was performed on a FACSAria II (BD Biosciences) and analyzed with FACSDiva software (BD Biosciences, RRID:SCR_016722), version 6.1.3. PI (250 ng/ml) was added to cells just before analysis. FACS data was imported into FlowJo software (Tree Star, Inc, RRID:SCR_008520), version 10, after which the scaled compensated values were exported as csv files. These values were then imported into Python 2.7 to perform rectangular gating. Cells were first gated using the forward scatter and propidium-iodide channels (i.e. cells negative for propidium-iodide were retained). Cells were subsequently gated for positive TOP-GFP expression. After this, cells below the 10th percentile and above the 90th percentile of TOP-GFP expression were compared. Gating strategy was described in the Registered Report with additional details available at https://osf.io/8c43g/ and a representative example depicted in Figure 1—figure supplement 1.

Conditioned medium

Request a detailed protocol

7.5 × 105 18 Co cells were seeded in a 75 cm2 flask and incubated overnight. The next day, cells were washed twice with PBS and incubated for 24 hr with 10 ml of CSC medium without EGF and FGF-basic. The next day the conditioned medium was collected and cleared by centrifugation for 5 min at 1400 RPM and used at 1:2 dilution in CSC medium for the assays described below. The level of HGF present in MFCM was determined by ELISA (Sigma-Aldrich, cat# RAB0212) according to manufacturer’s instructions with a standard curve. Data are available at https://osf.io/fpj4u/.

Limiting-dilution assay

Request a detailed protocol

An initial pilot experiment was performed to assess the potential for the three populations (TOP-GFPlow (10% lowest), TOP-GFPhigh (10% highest), and TOP-GFPwhole (total)), without treatment, on the three TOP-GFP CSC cultures (Co100, CSC1, E450). Cells from the different populations were deposited at 1, 2, 4, 6, 8, 12, 16, 20, and 24 cells per well in the number of wells indicated in the Registered Report with additional details available at https://osf.io/ydfrg/. Cells were deposited with an Influx cell sorter and analyzed with FACS Sortware sorter software, version 1.2.0.142. Cells were incubated at 37°C in a humidified atmosphere at 5% CO2, with culture medium replaced every 4 days. After 14 days of culture, the number of cultures with spheres, and the number of cells per sphere were quantified using automated high-content fluorescence imaging for GFP-positive and Hoechst-positive cells using a CellInsight CX7 High-Content Screening (HCS) Platform and HCS Studio Cell Analysis software. Spheres composed of two or more cells were used to determine clonal frequency which was evaluated by ELDA from the statmod R package (Hu and Smyth, 2009), version 1.4.30. Raw data are available at https://osf.io/ctqu2/ with data aggregated in csv format (https://osf.io/ydejb/). Pilot results reported in Figure 2—figure supplement 1. Based on these results, it was decided that the cell titration would remain the same for the confirmatory experiment, and that the E450 culture would be used for the in vivo tumorigenicity assay. Cells from the indicated TOP-GFP population were deposited into 96-well ultralow-adhesion plates with 100 µl of either CSC medium (untreated), CSC medium with 25 ng/ml HGF (Sigma-Aldrich, cat# H5791, lot# MKBT3102V), CSC medium with MFCM (1:2 dilution in CSC medium), CSC medium with 25 ng/ml HGF and 500 nM PHA-665752 (Sigma-Aldrich, cat# PZ0147), CSC medium with MFCM and 500 nM PHA-665752, or CSC medium with 500 nM PHA-665752. Cells were deposited with a FACSAria III (BD Biosciences) and analyzed with FACSDiva software, version 8.0.1. Cells were incubated at 37°C in a humidified atmosphere at 5% CO2, with the appropriate culture medium replaced every 4 days. After 14 days of culture, the number of cultures with spheres, and the number of cells per sphere were blindly quantified using automated high-content fluorescence imaging for GFP-positive and Hoechst-positive cells as described for the pilot assay. Raw data are available at https://osf.io/qwgx4/ with data aggregated in csv format (https://osf.io/26zp5/). Clonal frequency and statistical significance was determined by ELDA (Hu and Smyth, 2009).

In vivo tumorigenicity assay

Request a detailed protocol

All animal procedures were approved by the Explora BioLabs, Inc animal use committee (IACUC# SP17-009-005A) and were in accordance with Explora BioLabs, Inc policies on the care, welfare, and treatment of laboratory animals.

Nine-week old female Athymic Nude mice (Charles River, Strain code: 490, RRID:IMSR_CRL:490) were housed in sterile conditions under standard temperature, humidity, and timed lighting conditions with 12 hr light/dark cycles and acclimated to the housing environment for 3 days prior to the initiation of the study. They were housed on bedding material (Corn Cobb, cat# M-BTM-C8) that was changed bi-weekly. Animals were provided standard diet (Envigo, cat# 2920X (Irradiated Global 18% Soy Protein Extruded Rodent Diet)) and acidified water (pH 2.5–3.0) throughout the study period ad libitum. Body weights were measured on Day −1 for randomization and the 48 animals were stratified into 12 groups to obtain similar average body weight among groups. Following cell-injection, animal health, body weight, and tumor observation were recorded weekly and are available at https://osf.io/xs9up/. Tumor volumes were calculated from caliper measurements using the formula (volume = 1/2(length*width2). Experimental work was performed blinded to the identity of the sample the mice were injected with.

Mice were injected, on Day 0, with TOP-GFP transduced cultures (E450 culture) at 10, 100, 1000, or 5000 cells from the 10% lowest or 10% highest TOP-GFP intensities that were deposited, by FACS (FACSAria II with FACSDiva software, version 6.1.3), in a 96-well ultralow-adhesion plate and resuspended in 100 µl of CSC medium or MFCM (generated as described above and in the Registered Report [Evans et al., 2015]) and incubated at 37°C for 2 hr. After this incubation, plates were shipped to the facility that performed the mouse injection/monitoring (~30 min) where the cells and medium (100 µl) were mixed with growth factor reduced Matrigel (Corning, cat# 356230) at a 1:1 ratio and injected subcutaneously into the right flank of the female mice using a sterile 25 G needle and 1 ml syringe as described in the Registered Report (Evans et al., 2015). Mice were monitored for tumor formation for nine weeks after injection, the indicated study endpoint in the Registered Report. To explore if the frequency changed, the mice were monitored an additional 2 weeks; however we did not observe any new tumor initiations. Tumor-initiating cell frequency was determined by ELDA (Hu and Smyth, 2009).

Statistical analysis

Request a detailed protocol

Statistical analysis was performed with R software (RRID:SCR_001905), version 3.5.1 (R Development Core Team, 2018). All data, csv files, and analysis scripts are available on the OSF (https://osf.io/pgjhx/). Confirmatory statistical analysis was pre-registered (https://osf.io/rscpj/) before the experimental work began as outlined in the Registered Report (Evans et al., 2015). Data were checked to ensure assumptions of statistical tests were met. The fitted models to determine the stem cell frequency for different groups were compared using likelihood ratio tests using the asymptotic chi-square approximation to the log-ratio (Hu and Smyth, 2009). When described in the results, the Bonferroni correction, to account for multiple testings, was applied to the alpha error or the p-value. The Bonferroni corrected value was determined by divided the uncorrected value (0.05) by the number of tests performed. The confidence intervals for the Cohen’s ω estimates were determined using a Fisher’s z’ transformation (Rosenthal and DiMatteo, 2001). A meta-analysis of a common original and replication effect size was performed with a random effects model and the metafor R package (Viechtbauer, 2010), version 2.0–0 (https://osf.io/g4ewk/). The original study data of the stem cell frequency and 95% CI pertaining to Figure 6D was extracted a priori from the published figure and used to create simulated data sets with preserved sampling structure using ELDA (Hu and Smyth, 2009) during preparation of the experimental design, while the original study data pertaining to Figure 7E was published in the original paper (Vermeulen et al., 2010). The C100.B5 line from the original study was used in the meta-analysis for the tumorigenicity assay, but not the C100.G7 line because an estimate for TOP-GFPlow could not be calculated (i.e. estimate was infinity because of no observable responses). The summary data was published in the Registered Report (Evans et al., 2015) and used in the power calculations to determine the sample sizes for this study.

Data availability

Request a detailed protocol

Additional detailed experimental notes, data, and analysis are available on OSF (RRID:SCR_003238) (https://osf.io/pgjhx/; Essex et al., 2019). This includes the R Markdown file (https://osf.io/d6qp8/) that was used to compose this manuscript, which is a reproducible document linking the results in the article directly to the data and code that produced them (Hartgerink, 2017). Flow cytometry data for this study has also been deposited at Flow Repository (RRID:SCR_013779; Spidlen et al., 2012), where it is directly accessible at https://flowrepository.org/id/FR-FCM-ZYUG.

Deviations from registered report

Request a detailed protocol

We planned to isolate two CSC spheroidal cultures from patient samples, but due to only obtaining a single viable spheroidal culture from the 12 different colon tissue fragments that were attempted, we obtained another culture commercially that was also derived from primary human colorectal tumor tissue. The Registered Report indicated we would perform the flow cytometry analysis and clonogenicity assay on three different single-cell TOP-GFP clones from each of the three cultures, while the results reported in this study are from one random single-cell TOP-GFP clone from each culture. This was due to the CSC1 culture only producing one viable clone as stated in the ‘Lentiviral infection’ section above. We also did not perform the statistical analysis listed in Protocol 2 for the cell surface markers since the observed variation was from the same population of cells. Additional materials and instrumentation not listed in the Registered Report, but needed during experimentation are also listed.

Data availability

Additional detailed experimental notes, data, and analysis are available on OSF (RRID:SCR_003238) (https://osf.io/pgjhx/; Essex et al., 2019). This includes the R Markdown file (https://osf.io/d6qp8/) that was used to compose this manuscript, which is a reproducible document linking the results in the article directly to the data and code that produced them (Hartgerink, 2017). Flow cytometry data for this study has also been deposited at Flow Repository (RRID:SCR_013779; Spidlen et al., 2012), where it is directly accessible at https://flowrepository.org/id/FR-FCM-ZYUG.

The following data sets were generated
    1. Essex A
    2. Pineda J
    3. Acharya G
    4. Xin H
    5. Evans J
    6. Iorns E
    7. Tsui R
    8. Denis A
    9. Perfito N
    10. Errington TM
    (2019) FlowRepository
    Replication Study: Wnt activity defines colon cancer stem cells and is regulated by the microenvironment.

References

    1. Brattain MG
    2. Brattain DE
    3. Fine WD
    4. Khaled FM
    5. Marks ME
    6. Kimball PM
    7. Arcolano LA
    8. Danbury BH
    (1981)
    Initiation and characterization of cultures of human colonic carcinoma with different biological characteristics utilizing feeder layers of confluent fibroblasts
    Oncodevelopmental Biology and Medicine : The Journal of the International Society for Oncodevelopmental Biology and Medicine 2:355–366.
    1. McBain JA
    2. Weese JL
    3. Meisner LF
    4. Wolberg WH
    5. Willson JK
    (1984)
    Establishment and characterization of human colorectal cancer cell lines
    Cancer Research 44:5813–5821.
  1. Software
    1. R Development Core Team
    (2018)
    R: A language and environment for statistical computing
    R Foundation for Statistical Computing, Vienna, Austria.

Article and author information

Author details

  1. Anthony Essex

    PhenoVista Biosciences, San Diego, United States
    Contribution
    Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Conducted in vitro experimentation
    Competing interests
    PhenoVista Biosciences is a Science Exchange associated lab.
  2. Javier Pineda

    PhenoVista Biosciences, San Diego, United States
    Contribution
    Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Conducted in vitro experimentation
    Competing interests
    PhenoVista Biosciences is a Science Exchange associated lab.
  3. Grishma Acharya

    Explora BioLabs Inc, San Diego, United States
    Contribution
    Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Conducted in vivo experimentation
    Competing interests
    No competing interests declared
  4. Hong Xin

    Explora BioLabs Inc, San Diego, United States
    Contribution
    Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    No competing interests declared
  5. James Evans

    PhenoVista Biosciences, San Diego, United States
    Contribution
    Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Conducted in vitro experimentation
    Competing interests
    PhenoVista Biosciences is a Science Exchange associated lab.
  6. Reproducibility Project: Cancer Biology

    Contribution
    Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    1. tim@cos.io
    2. nicole@scienceexchange.com
    Competing interests
    EI, RT, NP: Employed by and hold shares in Science Exchange Inc.
    1. Elizabeth Iorns, Science Exchange, Palo Alto, United States
    2. Rachel Tsui, Science Exchange, Palo Alto, United States
    3. Alexandria Denis, Center for Open Science, Charlottesville, United States
    4. Nicole Perfito, Science Exchange, Palo Alto, United States
    5. Timothy M Errington, Center for Open Science, Charlottesville, United States
    6. Elizabeth Iorns, Science Exchange, Palo Alto, United States
    7. Rachel Tsui, Science Exchange, Palo Alto, United States
    8. Alexandria Denis, Center for Open Science, Charlottesville, United States
    9. Nicole Perfito, Science Exchange, Palo Alto, United States
    10. Timothy M Errington, Center for Open Science, Charlottesville, United States

Funding

The Reproducibility Project: Cancer Biology is funded by the Laura and John Arnold Foundation, provided to the Center for Open Science in collaboration with Science Exchange. The funder had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The Reproducibility Project: Cancer Biology would like to thank Dr Jan Paul Medema (University of Amsterdam) and Dr Giorgio Stassi (University of Palermo) for sharing critical information, data, and reagents, specifically the Co100 cells and TOP-GFP plasmid used in the original study. We want to thank Dr Laurie Ailles (University Health Network; University of Toronto) for providing permission for the TOP-GFP plasmid. We want to thank I Janette Delgadillo (Explora BioLabs, Inc) for providing study oversight for the animal experiments and Dr Lawrence R Blocher (Minerva Resource, Inc) for human colon adenocarcinoma tumor fragment collection. We would also like to thank the following companies for generously donating reagents to the Reproducibility Project: Cancer Biology; American Type and Tissue Collection (ATCC), Applied Biological Materials, BioLegend, Charles River Laboratories, Corning Incorporated, DDC Medical, EMD Millipore, Harlan Laboratories, LI-COR Biosciences, Mirus Bio, Novus Biologicals, Sigma-Aldrich, and System Biosciences (SBI).

Ethics

Human subjects: Approval was obtained from Western Institutional Review Board (WIRB) (protocol MR #0701) and were in full compliance with good clinical practices as defined under the US Food and Drug Administration (FDA) regulations, US Department of Health and Human Services (HHS) regulations, and the International Conference on Harmonisation (ICH) guidelines.

Animal experimentation: All animal procedures were approved by the Explora BioLabs, Inc animal use committee (IACUC# SP17-009-005A) and were in accordance with Explora BioLabs, Inc policies on the care, welfare, and treatment of laboratory animals.

Copyright

© 2019, Essex 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

  • 2,428
    views
  • 178
    downloads
  • 12
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

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)

  1. Anthony Essex
  2. Javier Pineda
  3. Grishma Acharya
  4. Hong Xin
  5. James Evans
  6. Reproducibility Project: Cancer Biology
(2019)
Replication Study: Wnt activity defines colon cancer stem cells and is regulated by the microenvironment
eLife 8:e45426.
https://doi.org/10.7554/eLife.45426

Share this article

https://doi.org/10.7554/eLife.45426

Further reading

    1. Cancer Biology
    Edited by Roger J Davis et al.
    Collection Updated

    Investigating reproducibility in preclinical cancer research.

    1. Cancer Biology
    2. Genetics and Genomics
    Yaroslav Kainov, Fursham Hamid, Eugene V Makeyev
    Research Article

    The expression of eukaryotic genes relies on the precise 3'-terminal cleavage and polyadenylation of newly synthesized pre-mRNA transcripts. Defects in these processes have been associated with various diseases, including cancer. While cancer-focused sequencing studies have identified numerous driver mutations in protein-coding sequences, noncoding drivers – particularly those affecting the cis-elements required for pre-mRNA cleavage and polyadenylation – have received less attention. Here, we systematically analysed somatic mutations affecting 3'UTR polyadenylation signals in human cancers using the Pan-Cancer Analysis of Whole Genomes (PCAWG) dataset. We found a striking enrichment of cancer-specific somatic mutations that disrupt strong and evolutionarily conserved cleavage and polyadenylation signals within tumour suppressor genes. Further bioinformatics and experimental analyses conducted as a part of our study suggest that these mutations have a profound capacity to downregulate the expression of tumour suppressor genes. Thus, this work uncovers a novel class of noncoding somatic mutations with significant potential to drive cancer progression.