Epigenetic regulation of mammalian Hedgehog signaling to the stroma determines the molecular subtype of bladder cancer
Abstract
In bladder, loss of mammalian Sonic Hedgehog (Shh) accompanies progression to invasive urothelial carcinoma, but the molecular mechanisms underlying this cancer-initiating event are poorly defined. Here, we show that loss of Shh results from hypermethylation of the CpG shore of the Shh gene, and that inhibition of DNA methylation increases Shh expression to halt the initiation of murine urothelial carcinoma at the early stage of progression. In full-fledged tumors, pharmacologic augmentation of Hedgehog (Hh) pathway activity impedes tumor growth, and this cancer-restraining effect of Hh signaling is mediated by the stromal response to Shh signals, which stimulates subtype conversion of basal to luminal-like urothelial carcinoma. Our findings thus provide a basis to develop subtype-specific strategies for the management of human bladder cancer.
https://doi.org/10.7554/eLife.43024.001eLife digest
In order to grow, cancer cells shut down or over-activate genes that normally maintain a cell’s health. The Sonic Hedgehog gene – named after a Japanese cartoon character – is associated with the cancer of several tissues, including the bladder. In 2014, researchers found that losing the Sonic Hedgehog gene, Shh for short, is necessary for bladder cancers to become aggressive: Shh signals prompt healthy cells near the tumor to inhibit the cancer cell growth, whilst aggressive bladder cancer cells turn off the Shh gene. Kim et al. – including many of the researchers involved in the 2014 work – now investigate how cancer cells switch off the Shh gene and what effect it has on bladder cancer cells and their surrounding tissue when turned back on.
DNA sequencing bladder cancer cells derived from human patients showed that there were no genetic deletions or mutations within the gene. However, the sequence and nearby regions of DNA did contain methylations – a chemical modification that generally switches genes off. When mice with early stages of bladder cancer were treated with a drug that inhibits methylation, the Shh gene turned back on, the bladder cancers stopped growing and the tumors stayed at an early stage of development. When the same drug was used on mice with aggressive bladder cancer, this caused non-cancer cells in the surrounding tissue to respond to Shh and send restraining signals back to the tumor. These signals eventually stopped cancer growth and converted the tumor into a less aggressive type of bladder cancer. Additionally, Kim et al. saw that blocking methylation had the same effect on human bladder cancer cells that had been transplanted into mice.
These results therefore indicate that Shh could be a new target for cancer treatments. For instance, drugs that decrease methylation and turn on the gene could be a way of managing cancer in patients with aggressive bladder cancers, which often show low activity of the gene. However, future studies are needed to understand what exactly happens within cancer cells during tumor conversion and to determine if this kind of intervention could have unintended consequences.
https://doi.org/10.7554/eLife.43024.002Introduction
Hedgehog (Hh) signaling has been recognized for its post-embryonic roles in the homeostatic maintenance of tissue integrity and the development of human malignancies (Ahn and Joyner, 2005; Goodrich et al., 1997; Shin et al., 2011; Taipale and Beachy, 2001). The initial identification of Hh pathway activity in human cancers, including basal cell carcinoma and medulloblastoma, has led to the development of the first FDA-approved drug targeting the Hh pathway for the treatment of human malignancy (Goodrich et al., 1997; Ruch and Kim, 2013; Sekulic et al., 2012; Tang et al., 2012), giving rise to a new field of pharmaceutical intervention (Teglund and Toftgård, 2010). Despite promising early preclinical studies (Olive et al., 2009; Yauch et al., 2008), recent studies investigating pancreatic, colon or ovarian cancers have shown that Hh pathway antagonism is not beneficial and clinical trials had to be halted in some cases because of accelerated cancer growth (Herter-Sprie et al., 2013; Kaye et al., 2012; Ruch and Kim, 2013). Consistent with the results of human trials, several recent studies have shown a protective role of Hh pathway activity in the progression of cancers that originate from endodermally derived tissues, including the bladder (Shin et al., 2014a; Shin et al., 2014b), pancreas (Lee et al., 2014; Rhim et al., 2014), colon (Gerling et al., 2016; Lee et al., 2016), and prostate (Yang et al., 2017). This tumor-restraining effect on a wide range of solid cancers is suggested to be exerted by the stromal response to Hh signals elicited from epithelial cancer cells.
In the bladder, an organ of endodermal origin, the Hh response is restricted to the stroma while Sonic hedgehog (Shh) protein is produced in basal epithelial cells. The regulatory circuit involving Hh signaling feedback between the bladder epithelium and supporting stromal cells is required for the proliferative response to injury during urothelial regeneration (Shin et al., 2011). Surprisingly, Shh expression is lost during the development of invasive urothelial carcinoma, even though the tumor is derived from Shh-expressing stem cells (Shin et al., 2014a), and genetic ablation of Hh signal response in stromal cells accelerates tumor progression at an early stage (Shin et al., 2014b). These studies have shown that loss of Shh expression invariably accompanies progression to invasive carcinoma and that suppression of the Hh response in the tumor stroma significantly accelerates the initiation of cancer, suggesting that Hh pathway activity protects against tumor progression at early stages of tumor development. Although the expression of Shh is invariably lost in both murine and human urothelial carcinoma, deep insights into the mechanisms underlying the tumor-cell-specific regulation of Shh at the early stage of carcinogenesis remain elusive. Interestingly, recent large-scale genomic studies in human bladder cancer (Cancer Genome Atlas Research Network, 2014a) have revealed that mutations of genes involved in epigenetic regulation are highly enriched in invasive urothelial carcinomas, while an extensive analysis by our group showed no mutational changes in the SHH gene; these results raise the possibility that epigenetic activities may be responsible for the loss of SHH expression during the initiation of urothelial carcinomas.
Technical advances in cancer genomics have permitted the subdivision of tumors into different molecular subtypes based on gene expression and mutational profiles (Cancer Genome Atlas Network and Cancer Genome Atlas, 2012; Cancer Genome Atlas Research Network, 2015). Recent large-scale genomic studies of gene expression in human urothelial carcinoma have revealed five distinct subtypes of bladder cancer (Cancer Genome Atlas Research Network, 2014a; Choi et al., 2014; Robertson et al., 2017), which can provide a strategic basis for developing personalized therapeutic interventions for individual patients with different molecular subtypes of urothelial carcinoma with genetic variability. An understanding of cellular and molecular dynamics, however, as cells evolve from the pre-cancerous state to distinct molecular subtypes of invasive carcinoma during tumor progression is required to develop more rationalized and precise treatment options for this malignancy and will require extensive experimental testing and validation through the integrative analysis of subtype-specific tumor initiation and progression beyond phenotypic analysis.
Interactions between epithelial cancer cells and the tumor stroma are important for the initiation and growth of human cancers (Calon et al., 2015; Isella et al., 2015; Mao et al., 2013). Our previous work has shown that Hh signaling to the stroma induces the expression of Bone morphogenetic proteins (BMPs), which impede bladder cancer progression. This anti-cancer effect of Hh-induced stromal expression of BMPs is mediated by urothelial differentiation of pre-cancerous cells at the early stage of tumor initiation (Shin et al., 2014b), indicating the importance of the stromal response elicited by tumor cells during tumor progression. Although the protective role of Hh response in the stroma at early stages of tumor initiation is interesting, it remains unknown whether the modulation of the Hh signaling to the stroma, especially increased activity of the Hh pathway, in full-fledged tumors would have a similar antitumor growth effect. This is of particular interest because most patients seen in a clinical setting are at a late stage of disease, with full-grown tumors. In this study, we elucidated the molecular basis for the loss of Shh during the development of bladder cancer and showed the role of the Hh signaling response in the tumor stroma in the determination of distinct molecular subtypes and the growth of full-fledged urothelial carcinoma.
Results
Loss of Shh expression in urothelial carcinoma results from hypermethylation of the CpG shore of the Shh gene
Having previously established the absence of SHH expression (Shin et al., 2014b) with a low incidence of genetic alterations (Figure 1—figure supplement 1A) and enrichment of mutations in genes involved in epigenetic regulation in human invasive urothelial carcinomas (Cancer Genome Atlas Research Network, 2014a), we compared the level of methylation in the regulatory region of Shh between wild-type bladders and N-butyl-N-4-hydroxybutyl nitrosamine (BBN)-induced urothelial carcinomas. By performing bisulfite sequencing analysis, we found significant increases in DNA methylation at the CpG shore upstream of the CpG island of the Shh promoter region (Figure 1—figure supplement 1B) in murine invasive urothelial carcinomas compared to that in wild-type bladders (Figure 1A,B). Pharmacological inhibition of DNA methyltransferase (DNMT) activity with 5’-azacitidine in BBN-induced urothelial carcinoma decreased the level of DNA methylation (Figure 1A,B), with significant increases in the expression of Shh (Figure 1C).
-
Figure 1—source data 1
- https://doi.org/10.7554/eLife.43024.005
In addition to primary murine tumors, we established 3D bladder tumor organoids derived from BBN-induced urothelial carcinoma (Figure 1—figure supplement 1C). Orthotopic transplantation of these tumor organoids revealed the histopathology of parental tumors (Figure 1—figure supplement 1D), suggesting that these organoids could recapitulate the pathology of the original BBN-induced urothelial carcinomas. To confirm that the loss of Shh expression was due to the increased methylation of the Shh gene, murine bladder organoids were cultured and treated with 5’-azacitidine, and the methylation status of the Shh regulatory region was analyzed. Consistent with the results from BBN-induced tumors, our analyses revealed that the CpG shore of the Shh promoter region in bladder tumor organoids were also hypermethylated, and following the treatment of 5’-azacitidine, the level of methylation was decreased (Figure 1D,E), with significant increases in the expression of Shh (Figure 1F). Our findings in BBN-induced urothelial carcinomas and bladder tumor organoids strongly suggested that the loss of Shh expression in invasive urothelial carcinomas results from the hypermethylation of the CpG shore in the Shh promoter region.
Pharmacological inhibition of DNA methylation halts the initiation of invasive urothelial carcinoma at the premalignant stage of progression through increased Hh signaling to the stroma
Our previous study showed that loss of the stromal Hh response triggers the initiation of invasive urothelial carcinoma and that elevated Hh signaling inhibits the development of bladder cancer at the early stage of progression (Shin et al., 2014a; Shin et al., 2014b). Given the experimental evidence that the expression of Shh is restored under 5’-azacitidine treatment (Figure 1C,F), we reasoned that inhibition of DNA methylation activity would impede the development of bladder cancer, especially at the early stage of tumor initiation. To test this possibility, we pharmacologically inhibited DNA methylation with 5’-azacitidine in our BBN-induced bladder cancer model. As previously established (Shin et al., 2014a), carcinoma in situ (CIS) became robust and widespread when mice were exposed to BBN for 4 months. This widespread CIS contains premalignant cells with Shh expression and represents a precursor lesion in muscle-invasive carcinoma. To investigate the effect of DNA methylation on bladder cancer initiation, mice were treated with BBN for 4 months, followed by the initiation of treatment with 5’-azacitidine at low dose for an additional 2 months (Tsai et al., 2012), while continuing exposure to BBN (Figure 1G). Without 5’-azacitidine treatment, BBN exposure for a total of 6 months resulted in the development of invasive carcinoma (Shin et al., 2014a; Shin et al., 2014b), whereas no invasive carcinoma was observed in mice treated with 5’-azacitidine during the final 2 months (Figure 1H, and Figure 1—figure supplement 1E,F), suggesting that inhibition of DNA methylation impeded tumor initiation if treatment was administered prior to the formation of invasive carcinoma.
To test whether the anticancer initiation effect of 5’-azacitidine was mediated by an increased stromal Hh response induced by increased expression of Shh in cancer cells, we attempted to rescue the tumor-restraining effect of 5’-azacitidine by genetically suppressing the Hh response in the stroma. As previously reported, Shh expression occurs in basal stem cells in the urothelium and responses to this signal are restricted to stromal cells (Shin et al., 2011). To genetically inactivate the stromal Hh response, we used Col1a2CreER;Smoflox/flox or Col1a2CreER;Gli2flox/flox mice expressing tamoxifen (TM)-inducible, stroma-specific CreER (Col1a2CreER) and carrying homozygous floxed alleles of essential Hh pathway components (Gli2 or Smoothened). These mice were exposed to BBN for 4 months and then injected with TM to genetically ablate Hh response in the stroma prior to the formation of invasive carcinoma. The mice were then continuously exposed to BBN for an additional two months in the presence of 5’-azacitidine (Figure 1I). In mice treated with TM to ablate the Hh response in the stroma, we found that the antitumor initiation effect of 5’-azacitidine was reversed and that invasive carcinomas appeared at 6 months, as in normal BBN-exposed mice, whereas no invasive carcinoma was observed in control animals (Figure 1J,K and Figure 1—figure supplement 1G,H). These findings suggested that DNA methylation of the Shh gene functions as a molecular basis for the loss of Shh expression in invasive urothelial carcinoma and confirmed the role of Hh signaling to the stroma in the initiation of bladder cancer at the early stage of disease, as previously reported (Shin et al., 2014a; Shin et al., 2014b).
Pharmacological inhibition of DNMT activity impedes the growth of urothelial carcinoma via increased activity of the BMP pathway, induced by the stromal Hh response
Although the increased stromal Hh response induced by inhibiting DNA methylation was shown to inhibit the transition of premalignant lesions to invasive carcinoma at the early stage of tumorigenesis, it remained unclear whether it exerted similar anticancer effects on the growth of mature urothelial carcinomas. To investigate this possibility, we used a recently established orthotopic transplantation model in which bladder cancer cells are intramurally injected into the wall of the bladder dome, allowing the transplanted tumor cells to propagate in physiologically relevant in vivo microenvironments (Shin et al., 2014a). Mice were orthotopically injected with BBN-induced tumor cells derived from isogenic mice and treated with 5’-azacitidine for 1.5 months after transplantation (Figure 2A). In the control group without inhibition of DNA methylation, tumor cells propagated and grew into full-fledged invasive carcinomas (Figure 2B and Figure 2—figure supplement 1A). In bladders from 5’-azacitidine-treated mice, however, no invasive carcinoma was observed (Figure 2B and Figure 2—figure supplement 1B), suggesting that inhibition of DNA methylation fully impeded the growth of bladder tumors in immunocompetent wild-type animals.
To investigate whether the anticancer propagation effect of 5’-azacitidine was mediated by the stromal Hh response, we combined our pharmacological approach using 5’-azacitidine with a genetic approach to genetically suppress the stromal Hh response while pharmacologically increasing the expression of Shh in tumor cells (Figure 2C). To genetically inactivate the stromal Hh response in recipient mice, we used the Col1a2CreER;Gli2flox/flox and Col1a2CreER;Smoflox/flox strains. After 5 days of recovery from TM injection, BBN-induced tumors, which are derived from mice with an isogenic background, were orthotopically transplanted, followed by treatment with 5’-azacitidine for 1.5 months (Figure 2C). In both strains with genetic ablation of the stromal Hh response, the anticancer growth effect of 5’-azacitidine disappeared (Figure 2D,E). These results strongly suggested that the effect of 5’-azacitidine on the propagation of tumor cells was mediated through the stromal Hh response elicited by Shh, whose expression is regulated epigenetically by cancer cells.
Next, we sought to determine whether the Hh signaling-mediated, antitumor propagation effect was regulated by Bmp, a secreted stromal factor whose expression is known to be regulated by the stromal Hh response in bladder (Shin et al., 2014b). Bmps are secreted stromal factors for urothelial differentiation (Mysorekar et al., 2009), and BMP pathway activity impedes bladder cancer progression prior to the formation of invasive carcinoma by stimulating urothelial differentiation (Shin et al., 2014b). However, the role of stromal Bmp in later stages of tumor development, especially in the tumor growth, is unknown. To determine whether stromal Hh response-regulated Bmp expression is involved in bladder cancer growth, we overexpressed Bmp4 in bladder tumor organoids derived from BBN-induced tumors (Figure 1—figure supplement 1C,D). The expression of Bmp4 in these organoids was increased by 10-fold compared with that in control organoids (Figure 2—figure supplement 1C). The resulting organoids with Bmp4 expression were orthotopically injected into Col1a2CreER;Smoflox/flox and Col1a2CreER;Gli2flox/flox mice, and the mice were then injected with TM to genetically ablate the stromal Hh response. These animals were subsequently treated with 5’-azacitidine for 1 month to increase Shh expression in tumor cells (Figure 3A). We found that, compared with wild-type bladder tumor organoids, Bmp4-expressing tumor organoids showed growth reductions after transplantation (Figure 3B,C and Figure 2—figure supplement 1D,E). In addition, after tumor organoids were cultured in the presence of the Bmp4 protein to increase the BMP response (Figure 2—figure supplement 1F), the Bmp4-treated organoids grew slowly, and the efficiency of organoid formation was significantly reduced (Figure 3D,E,F). Taken together, the tumor-suppressing effect of Bmp, particularly in vivo where tumor cells still express Shh (with 5’-azacitidine treatment) while the stromal Hh response is suppressed (using Col1a2CreER;Smoflox/flox and Col1a2CreER;Gli2flox/flox mice), supports a potential scenario of an increased reciprocal tumor-stromal signal feedback loop in which hypomethylation-induced Shh secretion by tumor cells activates the Hh response in bladder stroma, resulting in stromal expression of Bmps, which in turn signal back to tumor cells to impede their growth.
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.43024.009
Heightened activity of Hh signaling to the stroma induces a less aggressive luminal subtype of urothelial carcinoma
To investigate the cellular basis of the cancer-restraining effects of the stromal Hh response induced by Shh, which is regulated by DNA methylation in tumor cells, on the growth of bladder tumors, we performed the orthotopic transplantation of BBN induced tumors to nude mice. These mice were chosen to grow transplanted tumors under more permissive conditions because orthotopic transplantation of bladder tumors to wild-type mice in the presence of 5’-azacitidine led to complete blockade of tumor growth (Figure 2A,B); thus, the use of nude mice under more permissive conditions for tumor growth allowed us to overcome the difficulty in studying the basis of the anticancer effect of the hypomethylation-induced stromal Hh response on tumor growth. Indeed, transplanted bladder tumors in nude mice grew even under 5’-azacitidine treatment, but they gave rise to tumor lesions with smaller sizes than those without 5’-azacitidine treatment (Figure 4A–C and Figure 4—figure supplement 1A,B). These findings suggest that 5’-azacitidine treatment is still effective in suppressing tumor growth under immunocompromised conditions, which is consistent with our earlier results (Figure 2A,B).
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.43024.012
Next, we sought to investigate the molecular and cellular basis for the less aggressive growth of bladder tumors when Hh signaling activity is increased by inhibiting DNA methylation. As shown previously (Figures 2 and 3), the anticancer effects of Hh signaling appeared to be mediated by stromal BMP, whose signaling activity is known to be associated with the urothelial differentiation of basal cells into luminal cells (Mysorekar et al., 2009). As our previous study on the cellular origin of bladder cancer showed that urothelial carcinoma is derived from basal stem cells (Shin et al., 2014a), we hypothesized that increased Hh signaling activity might cause the tumor to differentiate into a less aggressive form of luminal subtype, leading to much slower growth upon 5’-azacitidine treatment. As previously reported (Fantini et al., 2018; Shin et al., 2014b), based on the expression level of basal markers and mutational profile, the invasive carcinomas produced in our BBN model are most similar to the basal subtype of human urothelial carcinoma (Figure 4—figure supplement 1C), which is the most aggressive form of bladder cancer (Choi et al., 2014; Vasconcelos-Nóbrega et al., 2012). Thus, we examined BBN-induced bladder tumors that were orthotopically grown in the presence of 5’-azacitidine to investigate the cellular differentiation of transplanted tumors. Our immunohistochemical analysis revealed that the expression of the luminal marker, Ck18, was increased in the tumors treated with 5’-azacitidine, while control bladders displayed basal phenotypes such as squamous differentiation and the expression of Ck5, a marker for the basal subtype (Figure 4D,E and Figure 4—figure supplement 1D,E). Moreover, our quantitative RT-PCR experiments showed that the expression level of luminal markers relative to basal markers was significantly increased in the transplanted bladder tumors under 5’-azacitidine treatment compared with that of the control group, which was accompanied by increased expression of Shh (Figure 4F and Figure 4—figure supplement 1F). Consistent with these results, our analysis of RNA-seq expression profiles revealed the strong luminal signature with relatively low expression of basal markers in the tumors orthotopically grown in the presence of 5’-azacitidine, whereas control allografts in the absence of 5’-azacitidine showed clear standard signature of basal subtype (Figure 4G and Figure 4—figure supplement 1G) (Damrauer et al., 2014). These results suggested that subtype conversion between basal and luminal-like subtypes resulting from the increased Hh response induced by epigenetically upregulated expression of Shh in tumor cells might account for the reduced tumor growth.
Stromal Hh response-induced conversion from the basal to luminal-like subtype requires BMP pathway activity
To investigate whether conversion from the basal to luminal-like subtype upon 5’-azacitidine treatment is mediated by elevated expression of Hh in tumor cells, we performed orthotopic transplantation of BBN-induced tumor organoids engineered to express shRNA targeting Shh (Figure 5A and Figure 5—figure supplement 1A). The resulting organoids maintained low levels of Shh expression after transplantation, even upon 5’-azacitidine treatment, whereas the control organoids showed increased expression of Shh (Figure 5—figure supplement 1B). The expression levels of luminal markers were significantly decreased in tumor organoids following the expression of the shRNA targeting Shh after transplantation, indicating the basal-like characteristics including squamous differentiation of the resulting tumors, even under 5’-azacitidine treatment (Figure 5B,C and Figure 5—figure supplement 1C,D). Consistent with these results, our analysis of RNA-seq expression profiles confirmed that the gene signature related with the luminal status is negatively enriched in the tumors expressing shRNA for Shh, whereas tumor expressing control shRNA showed standard signature of luminal-like subtype (Figure 5D and Figure 5—figure supplement 1F). Similar results were noted when tumor organoids expressing shRNA targeting Shh, marked by mCherry expression, were mixed and transplanted together with control organoids, marked by EGFP, into the same in vivo microenvironment (Figure 5F). The resulting allografts revealed that mCherry-labeled tumors developed into a more aggressive, rapidly growing basal-like subtype, whereas, in the same microenvironment, EGFP-labeled tumors developed into a less aggressive luminal-like subtype following treatment with 5’-azacitidine (Figure 5G and Figure 5—figure supplement 2A), demonstrating Hh-mediated conversion of the bladder tumor subtype.
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.43024.017
To further determine whether conversion between the basal and luminal-like subtypes requires Hh-mediated BMP signaling, whose activity is necessary for the suppression of tumor growth (Figure 3), we established transplant models in which BBN-induced tumor organoids transduced to express shRNA targeting Bmpr1a were orthotopically injected into the mouse bladder (Figure 5A and Figure 5—figure supplement 1A). The expression of Bmpr1a was significantly decreased in the resulting tumor organoids (Figure 5—figure supplement 1B), and transplantation of Bmpr1a knock-down tumor organoids in the presence of 5’-azacitidine gave rise to secondary tumors with decreased expression of luminal markers with squamous differentiation compared to tumor organoids with intact Bmpr1a (Figure 5B,C and Figure 5—figure supplement 1C,E). Consistent with these results, our analysis of RNA-seq expression profiles confirmed that the gene signature related with the luminal status is negatively enriched in the tumors expressing shRNA for Bmpr1a (Figure 5E and Figure 5—figure supplement 1G), whereas tumor expressing control shRNA showed standard signature of luminal-like subtype. Similar data were obtained when tumor organoids expressing shRNA targeting Bmpr1a, marked by mCherry expression, were mixed and transplanted together with control organoids, marked by EGFP, into the same in vivo microenvironment (Figure 5F). The mixed allografts showed that mCherry-labeled tumors developed into a more rapidly growing basal-like subtype than EGFP-labeled tumors with the luminal-like subtype, even under 5’-azacitidine treatment (Figure 5H and Figure 5—figure supplement 2B). Furthermore, we genetically ablated Bmpr1a from BBN-induced tumor organoids derived from Bmpr1aflox/flox mice (Mishina et al., 2002) by expressing Cre recombinase (Figure 5—figure supplement 3A,B). Consistent with the findings described above, the resulting organoids developed into basal invasive carcinoma, even with 5’-azacitidine treatment (Figure 5—figure supplement 3C,D).
Taken together, our results from various genetic experiments in combination with pharmacological approaches to manipulate Hh and BMP signaling feedback during the growth of bladder tumors strongly suggest that subtype conversion between basal and luminal-like subtypes depends on the reciprocal signaling feedback between tumor cells and the stroma involving epigenetically regulated epithelial Hh expression in tumor cells, stromal Hh response-induced Bmp expression, and the BMP response in tumor cells.
Increased methylation of the SHH gene induces the basal subtype of human urothelial carcinoma and promotes tumor growth through decreased activity of Hh/BMP signaling feedback between cancer cells and tumor stroma
To determine whether Hh/BMP signaling feedback between the tumor and stroma could control the growth of tumors and determines their subtype in human bladder cancer, we first examined multiple bladder cancer cell lines derived from human invasive urothelial carcinomas. Three invasive bladder cancer cell lines, J82, T24 and TCC-SUP, were chosen (Bubeník et al., 1973; Nayak et al., 1977; O'Toole et al., 1978) and analyzed for the level of methylation in the regulatory region of SHH by bisulfite sequencing. Consistent with our mouse experiments described above (Figure 1), all three cell lines showed a significant increase in DNA methylation at the CpG shore of the SHH promoter region (Figure 6A,B and Figure 6—figure supplement 1A) and demonstrated basal characteristics of urothelial carcinoma with low levels of SHH expression. Pharmacological inhibition of DNMT activity with 5’-azacitidine in these cell lines significantly decreased the level of DNA methylation (Figure 6A,B), which was associated with marked increases in the expression of SHH (Figure 6C). Interestingly, a cell line derived from the luminal-papillary subtype RT4 (Rigby and Franks, 1970), exhibiting high expression of SHH, showed no significant changes in the expression of SHH, even with the DNMT inhibition (data not shown).
-
Figure 6—source data 1
- https://doi.org/10.7554/eLife.43024.020
To investigate the functional roles of SHH expression in the growth of human bladder tumors and the effects of Hh/BMP signaling feedback between the tumor and stroma on the subtype determination of human invasive urothelial carcinoma, we established an orthotopic xenograft model in which one human invasive bladder cancer cell lines J82 was transplanted into immunocompromised mice (Nod/Scid/Rag2), followed by 5’-azacitidine treatment for 1 month (Figure 6D). In the control group without inhibition of DNA methylation, tumor cells developed into full-fledged invasive carcinomas (Figure 6E and Figure 6—figure supplement 1B). In bladders from 5’-azacitidine treated mice, however, much smaller lesions were observed (Figure 6E and Figure 6—figure supplement 1C), suggesting that inhibition of DNA methylation inhibited the growth of human bladder tumor. In addition, xenografts from 5’-azacitidine-treated animals showed characteristics of the luminal-like subtype, with increased expression of luminal markers (Figure 6F), consistent with our observations in the above-described mouse model (Figure 4).
To further confirm the requirement for Hh/BMP signaling feedback in the subtype determination of human bladder cancer, we performed xenograft transplantation of a human bladder cancer cell line that was engineered to express shRNA targeting SHH or BMPR1A (Figure 6G). After transplantation, the resulting tumors maintained a low level of SHH or BMPR1A expression (Figure 6H), and the expression levels of basal markers were significantly increased, with low expression levels of luminal markers being observed even upon 5’-azacitidine treatment (Figure 6I). These human data are consistent with our observations in the mouse model, showing that the Hh pathway response in the tumor stroma induced by elevated expression of Hh in tumor cells impedes the growth of bladder cancer through stromal BMP-induced subtype conversion.
Expression analysis of patient-derived urothelial carcinomas and large-scale transcriptional analyses suggest an association of the basal subtype of human urothelial carcinoma with poor patient survival, due to the loss of Hh signaling activity
To validate the role of the loss of SHH induced by hypermethylation in the determination of the molecular subtype of human bladder cancer, we examined 10 muscle-invasive urothelial carcinoma samples derived from patients (Table 1). We found that the invasive carcinoma samples from six patients showed characteristics of the basal subtype, with the expression of basal markers being significantly increased (Figure 7A, Figure 7—figure supplement 1). The invasive carcinoma samples from the remaining four patients exhibited a luminal phenotype, with elevated expression of luminal markers (Figure 7A, Figure 7—figure supplement 1). Consistent with the data from our BBN model and human bladder cancer cell lines (Figures 4 and 6), the luminal phenotype of the human primary bladder cancer samples was accompanied by a marked increase in the expression of SHH, whereas human primary tumors with the basal phenotype showed a significant decrease in SHH expression (Figure 7B).
-
Figure 7—source data 1
- https://doi.org/10.7554/eLife.43024.024
To further determine whether our findings in the BBN mouse model showing that the molecular subtype of bladder cancer is determined by the expression of Shh regulated by DNA methylation are valid in human urothelial carcinoma, we performed methylation analysis by bisulfite sequencing in six basal subtypes and three luminal subtypes of primary invasive carcinoma samples (Table 1), and then compared their methylation status to that of three benign urothelia. We found that the methylation of the SHH gene, especially in the CpG shore of the promoter region, was significantly increased in the basal subtype of human invasive carcinomas compared to that in benign urothelia (Figure 7C,D). Decreased methylation in luminal subtypes of invasive carcinoma relative to basal subtypes was also observed in these human samples (Figure 7C,D). Our analysis of primary human bladder cancer samples for the expression and associated methylation status of SHH was consistent with our findings from the murine BBN model of bladder cancer.
As the basal subtype of human urothelial carcinoma is the most aggressive form of bladder cancer, associated with lower median overall survival (Choi et al., 2014), we sought to determine whether the SHH expression level was associated with clinical outcome in human bladder cancer patients. We analyzed patient outcomes among a set of 41 muscle-invasive bladder tumors from Seoul National University Hospital (Table 2) that were clustered into two groups, showing low SHH expression (n = 31) or high SHH expression (n = 10). We found that the group with high SHH expression survived for a median of 50 months, compared to 39 months in the group with low SHH expression. The patient group with higher SHH expression therefore had better patient outcomes, presenting a 28% survival benefit compared to the group with lower SHH expression (log rank test p<0.05, Figure 7E).
In addition, we analyzed data from recently published large-scale studies of gene expression in human invasive carcinoma (Robertson et al., 2017) to assess the association of SHH expression with the molecular subtypes and patient survival. First, we analyzed the clinical outcomes of two subgroups stratified by Hh and BMP pathway genes (Figure 6—figure supplement 1D). We noted that the group with low SHH expression showed molecular characteristics of the basal subtype with poor patient outcomes compared to those of the high-SHH group with luminal characteristics (Figure 6—figure supplement 1E,F). Based on unsupervised hierarchical clustering, we further identified three distinct subgroups (Figure 7F): basal, luminal, and p53-like subtypes, as previously reported (Choi et al., 2014; Robertson et al., 2017). We noted that the basal cluster showed the lowest level of SHH expression (Figure 7F), consistent with our previous analysis and results from other studies (Robertson et al., 2017; Shin et al., 2014b), and the median overall survival of patients with the basal subtype was significantly lower than that of patients with the luminal subtype (log rank test p<0.0001, Figure 7G). It is important to note that only a small fraction of the patients carried genetic mutations in SHH and its downstream genes, suggesting that genetic alterations might not have been the major cause of the loss of SHH expression signatures in human invasive bladder cancer (Figure 1—figure supplement 1A). These results support our findings from the BBN mouse model and the analysis of primary human samples showing that the expression of Shh was regulated by epigenetic modifications, rather than mutational changes during the development and growth of invasive urothelial carcinoma (Figures 1, 6 and 7).
Taken together, our findings from extensive analyses of primary tumor samples and a TCGA dataset suggest a possible basis for the development of the basal subtype of human invasive carcinomas that lack active Hh signaling. Under this scenario, loss of SHH expression due to hypermethylation leads to the decreased expression of stromal BMP, which in turn stimulates the formation of the basal subtype of human bladder cancer, with poor clinical outcomes.
Discussion
The proposed events during the development of urothelial carcinoma, progressing from Shh-positive basal stem cells to distinct molecular subtypes, are summarized in Figure 8. Following the formation of widespread CIS, the clonal loss of Shh expression caused by DNA methylation leads to suppression of the stromal Hh response, which in turn increases the number of undifferentiated premalignant cells with proliferative advantages due to the decreased expression of stromal differentiation factors. This critical event facilitates the initiation and progression of tumors at early stages, ultimately resulting in the formation of the Shh-negative, basal subtype of invasive carcinoma. This model is strongly supported by our observations of halted initiation at the early stage and significantly decreased progression to invasive carcinoma following an increased stromal response to Hh signal upon the treatment with a potent DNMT inhibitor to epigenetically increase the expression of Shh. Once formed, the basal subtype of urothelial carcinoma remains Shh negative, consistent with the observed loss of Shh in murine and human invasive urothelial carcinoma with basal characteristics (Cancer Genome Atlas Research Network, 2014a; Cancer Genome Atlas Research Network, 2014b; Shin et al., 2014b). In contrast, invasive urothelial carcinoma of the luminal subtype develops through the distinct papillary lineage, as Shh expression is maintained, with a low CIS mutational profile. One key feature of our model is that a luminal-like subtype can be induced from the basal subtype of urothelial carcinoma by pharmacological modulation with 5’-azacitidine, to increase the expression of Shh in tumor cells, resulting in an increased stromal Hh response. This observation raises the possibility of additional pathways leading to the development of the luminal subtype of bladder cancer; a gain of Shh expression due to epigenetic plasticity in the basal subtype induces the luminal subtype of urothelial carcinoma. The Hh-dependent conversion of tumor subtypes is likely mediated by stromal Bmps, as supported by our observation that increased Bmp expression in tumor cells dramatically facilitates basal-to-luminal transition.
Epigenetic modifications, such as DNA methylation (Herman and Baylin, 2003) and histone modifications (Seligson et al., 2005) of genes associated with carcinogenesis, are commonly detected in human malignancies (Feinberg and Tycko, 2004; Jones and Baylin, 2002). Our previous findings of the invariable loss of Shh in murine and human urothelial carcinomas (Shin et al., 2014b), together with studies of other solid tumors, strongly support a role of Shh as a tumor suppressor during cancer progression (Lee et al., 2014; Rhim et al., 2014; Yang et al., 2017), raising the important questions concerning the molecular mechanisms by which the expression of Shh is lost in tumor cells at early stages of carcinoma development. A recent TCGA analysis of human bladder cancers revealed frequent mutations in genes associated with epigenetic regulation, but interestingly, showed no mutational changes in the SHH gene, suggesting the possibility that epigenetic changes are responsible for the loss of Shh expression during the development of urothelial carcinomas. Our study reveals that hypermethylation of the CpG shore, the upstream region with a reduced CG density of the CpG island in the Shh promoter region, is responsible for the loss of Shh, and that pharmacological inhibition of DNA methylation with agents such as 5’-azacitidine induces increased expression of Shh, thereby halting the progression of urothelial carcinoma. These findings are supported by other studies that show the importance of the CpG shore regions, rather than CpG islands, of genes involved in tumorigenesis (Ellis et al., 2009; Irizarry et al., 2009). It is, however, still not clear regarding the potential upstream regulators that control Shh methylation during the development of urothelial carcinoma. A previous TCGA analysis showed frequent mutations in methyltransferases and demethylases, such as MLL2, ARID1A, and KDM6A, raising the interesting possibility that the increased Shh methylation is due to the mutations in these genes associated with epigenetic regulation. Further studies will be necessary to investigate the specific roles of these genetic changes in the epigenetic regulation of the Shh gene. Furthermore, it would be interesting to determine whether other epigenetic events regulating chromatin organization, such as histone modifications, are also involved in the loss of Shh and whether the complex enhancer of the Shh gene, a cluster of three long-range, 840 kb regions (Lettice et al., 2003; Sagai et al., 2009), is associated with the development of urothelial carcinoma.
One of the most interesting findings of the present study is the remarkable ability of tumors to change their phenotypic characteristics without significant alterations in their genetic profiles. Conventionally, cancer has been considered as a disease of ‘sequential accumulation of mutations’, and tumor phenotypes are dictated by their genetic alterations. However, recent advances in cancer biology and genomics have modified this traditional view, demonstrating that cancer is not a single-factorial genetic disease and that tumors exhibit a great degree of plasticity. Our observation that the tumor subtypes are epigenetically switched through the dynamic interactions with the tumor stroma—mediated by Hh-BMP signaling feedback—provides a clear basis for this recent view regarding tumor plasticity. In the normal bladder, our previous study has indicated that the development and regeneration of the urothelium depend on the balanced proliferation of basal stem cells and their differentiation into luminal cells, which is controlled by the stromal production of multiple factors (Shin et al., 2011). It is interesting to speculate that in cases of bladder tumors, the hierarchical relationship between basal stem cells and luminal cells observed in the normal bladder is maintained as basal subtype is converted into luminal subtype of bladder cancer. Interestingly, all luminal-like tumors induced by 5’-azacitidine treatment in our study do not show histological features of the basal subtype (e.g. squamous differentiation), although some luminal-like tumors still maintain the low expression of certain basal markers. This suggests the interesting possibility that there may be the intermediate subtype between basal and luminal human bladder cancer, similar to the intermediate cells in the urothelium of the normal bladder. Moreover, in our rescue experiment, we noticed that suppressing Hh-BMP feedback signaling with individual shRNAs for Shh or Bmpr1a along with 5’-azacitidine treatment was insufficient to regain the gene expression of the strong basal signature based on the GSEA analyses (Figure 5D,E and Figure 5—figure supplement 1F,G), although the reverted tumors show increased expression of certain basal markers (e.g. CK5) and the strong histological characteristic of the basal subtype (e.g. squamous differentiation) (Figure 5B). These data imply the potential involvement of additional mechanisms regulated by DNA methylation—other than Hh-BMP feedback signaling—to support the complete conversion of bladder tumors into the luminal subtype, particularly the complete loss of the basal signature, because the tissue dynamics of normal urothelium are controlled by multiple stromal factors such as Wnt, Fgf, and Bmp (Shin et al., 2011; Shin et al., 2014b). Further studies are required to investigate this possibility. Taken together, our findings on subtype conversion provided further insights on the current understanding of tumor plasticity and tumor interaction with other tissue components for the development of more effective therapeutic interventions for urothelial carcinomas.
Early detection and prevention are important therapeutic options for the treatment of cancers (Tomasetti et al., 2017; Tomasetti and Vogelstein, 2015), and a better understanding of tumor initiation and progression at early stages is crucial for the development of the most effective prevention strategies for cancer therapy. Pharmacological manipulation with 5’-azacitidine to increase Shh expression, resulting in therapeutic activation of stromal Bmp to block bladder tumor initiation through the differentiating of premalignant cells, could have a potentially beneficial impact on the clinical management of human bladder cancer, especially at the early stage of onset or recurrence, which makes this malignancy difficult and more costly to treat than other cancers. In addition to clinical benefits in terms of cancer prevention, our findings regarding the tumor-restraining effects of the stromal Hh response induced with pharmacological agents in full-fledged tumors, carry significant therapeutic implications, especially for patients with late-stage of invasive urothelial carcinoma. A recent TCGA analysis revealed five different subtypes of muscle-invasive urothelial carcinoma with different patient outcomes (Cancer Genome Atlas Research Network, 2014a; Choi et al., 2014; Hedegaard et al., 2016; Robertson et al., 2017). In these works, human invasive urothelial carcinomas were subdivided into five distinct molecular subtypes: luminal-papillary, luminal-infiltrated, luminal, basal and neuronal. Our results show that the loss of Shh induced by DNA methylation of the CpG shore of the Shh promoter region leads to the initiation and, ultimately, growth of the basal subtype of urothelial carcinoma. With pharmacological induction of the expression of Shh in epithelial tumor cells of the basal subtype, which in turn increases the stromal Hh response, urothelial carcinoma develops into a less aggressive, more manageable luminal-like subtype of bladder cancer through stromally produced Bmp-mediated cellular differentiation. The luminal subtype of urothelial carcinoma appears to grow more slowly, to be associated with better patient outcomes and to be more susceptible to the immune system than the basal subtype (Figures 2 and 3). Interestingly, a recent study showed a phenotypic shift of luminal parental tumors to basal subtypes in patient-derived organoid culture, and showed that this change is reversed in organoid xenografts in which the tumor stroma is restored (Lee et al., 2018). In the present study, we observed epigenetically regulated Hh expression in cancer cells, which in turn induces a stromal response in the tumor microenvironment to determine the molecular subtype of bladder cancer. Our findings of Hh-dependent subtype switching through the tumor stroma may provide a strong rationale for the phenotypic plasticity observed in the human bladder tumor organoid model, which lacks a tumor stroma. The ability to convert molecular subtypes by modulation of Hh signaling with pharmacological agents such as 5’-azacitidine could therefore facilitate the development of personalized therapeutic options for urothelial carcinoma and permit potential combinatorial therapy. For example, treatment with 5’-azacitidine and/or FK506 (Shin et al., 2014b; Spiekerkoetter et al., 2013) to activate the Hh response-induced expression of cellular differentiation factors, such as BMPs, in the tumor stroma would induce the luminal-like subtype, which is more manageable with subsequent treatment, including immunotherapy (Massard et al., 2016; Powles et al., 2014).
In the present study, we speculate that several additional key events may be associated with Hh signaling during the initiation and growth of urothelial carcinomas not belonging to the basal and luminal subtypes, for which the mechanisms of progression are described above (Figure 8). Non-muscle-invasive carcinoma, for example, is derived from Shh-positive basal epithelial cells, which have been previously identified as cells of origin of bladder cancer (Shin et al., 2014a), via a papillary pathway by maintaining active Hh signaling, in which Shh is expressed in the epithelium and the Hh response occurs in the stroma. The Hh response in the stroma induces the expression of secreted factors such as Wnts and Fgfs, which in turn stimulate the proliferation of epithelial cells (Shin et al., 2011). With proliferative advantage of increased signaling feedback between the epithelium and stroma and the additional mutations of genes such as Fgfr3, papillary invasive tumors develop while maintaining Shh expression and low CIS signatures. The luminal–infiltrated subtype also follows the papillary lineage, but rather than mutations in Fgfr3, it accumulates mutations associated with epithelial-mesenchymal transition (EMT) while maintaining a moderate level of Shh expression and showing no mutations in p53. Finally, the neuronal subtype characterized by the expression of neuroendocrine and neuronal genes (Sjödahl et al., 2017) may originate from neuroendocrine cells, but not basal stem cells, and may develop in an Hh signaling-independent way. This speculation arises from the results of our bioinformatic analysis of the TCGA data for the levels of SHH and BMP in five distinct subtypes of invasive urothelial carcinoma (Figure 6—figure supplement 1G). Taken together, our study identifying several key events and their association with Hh signaling during the development of different subtypes of bladder cancer provides the basis for the design of precise therapeutic strategies for individual patients with genetic variability.
Materials and methods
Mice
For the genetic deletion experiments, Col1a2CreER (RRID:IMSR_JAX:029235) mice were mated with the Smoflox/flox (RRID:IMSR_JAX:004526) or Gli2flox/flox (RRID:IMSR_JAX:007926) strains to obtain Col1a2CreER;Smoflox/flox, Col1a2CreER;Gli2flox/flox mice, which were administered 8 mg of TM (Sigma) per 30 g of body weight on 3 consecutive days by oral gavage. Male mice at 8–10 weeks of age were used. For experiment involving 5’-azacitidine, mice were injected intraperitoneally with 1 mg of 5’-azacitidine (Sigma) per 1 kg of body weight daily. The duration of the dosing is described in each figure. In each experiment, mice in each cage were randomly selected for drug/TM or control treatments. Mouse procedures were performed under isoflurane anesthesia with a standard vaporizer. All procedures were performed under a protocol approved by the Institutional Animal Care and Use Committee at POSTECH (IACUC number: POSTECH-2017–0094).
BBN-induced bladder carcinogenesis
Request a detailed protocolA 0.1% concentration of BBN (TCI) was dissolved in drinking water, and BBN-containing water was provided to mice ad libitum for 4 to 6 months in a dark bottle. BBN-containing water was changed twice a week. Bladders were collected and analyzed after 4 to 6 months of BBN administration.
Methylation analysis using bisulfite sequencing of genomic DNA
Request a detailed protocolThe DNA methylation status of murine and human Shh was determined using bisulfite genomic DNA sequencing. For bisulfite conversion, 1 μg of genomic DNA was converted using the MethylEdge Bisulfite Conversion System (Promega), following the manufacturer’s instructions. The genomic sequence of the regulatory region of murine and human Shh was obtained from the NCBI nucleotide database (Mus musculus: NC_000071.6, Homo sapiens: NG_007504.2), and the CpG island and CpG shores in the regulatory region were identified by Methprimer 2.0 (Li and Dahiya, 2002) (RRID:SCR_010269). The 2 kb regions upstream and downstream of the CpG island were referred to as the ‘CpG upshore’ and ‘CpG downshore’ regions, respectively. For sequencing analysis, bisulfite-converted DNA was amplified by EpiTaq HS (TaKaRa), and PCR products containing the CpG island and CpG shore regions were subcloned into the pGEM-T easy vector (Promega). The region containing the CpG island and shores was divided into eight sub-regions, and each sub-region was amplified using specific primers designed for bisulfite-converted target sequences (summarized in Table 3). The sequencing data were assembled using SnapGene software (https://snapgene.com/, RRID:SCR_015053) and the MUSCLE: multiple sequence alignment tool (https://www.ebi.ac.uk/Tools/msa/muscle/, RRID:SCR_011812). The average degree of methylation was obtained from the analysis of 8–10 clones of each sub-region. The methylated CpG sites were counted and distinguished from the unmethylated CpG sites.
Bladder organoid culture
Request a detailed protocolBBN-induced bladder tumors were minced and then incubated in DMEM (Gibco) containing collagenase I, II (20 mg/ml each) and thermolysin (250 KU/ml) at 37°C for 2 hr, with 5 min trituration every 30 min. A single-cell suspension was obtained and filtered through 100 µm cell strainers (Falcon). After lysis of red blood cells in ACK lysis buffer (Gibco), the cells were washed with DMEM containing 10% fetal bovine serum (Millipore) and counted using a hemocytometer (Sigma). For bladder organoid culture, single tumor cells were overlaid in growth factor-reduced Matrigel (Corning) and incubated with advanced DMEM/F-12 (Gibco) supplemented with 10 mM HEPES (pH 7.4, Sigma), 10 mM Nicotinamide (Sigma), 1 mM N-acetyl-L-cysteine (Sigma), GlutaMAX (Gibco), 1% penicillin/streptomycin (Gibco), 50 ng/mL mouse EGF (Peprotech), 0.5X B-27 (Gibco), 1 μM A8301, and 10 μM Y-27632. For Bmp4 treatment, organoids were treated with recombinant Bmp4 protein (Peprotech) for 8 days, with medium transition every 2 days. For the knock-down experiments, bladder tumor organoids were infected with lentivirus containing shRNAs specific for mouse/human Shh or Bmpr1a (Table 4). For lentivirus production, 293 T cells were co-transfected with lentiviral plasmids and packaging vectors (pCMV.dR 8.74 and pMD2.G [RRID:Addgene_12259], Addgene) using TransIT-LT1 (Mirus Bio). The supernatant was collected 48 hr post-transfection and filtered through 0.45-μm-pore PES filter (Millipore). The viral titer was calculated in 3T3 cells by serial dilution of the virus-containing supernatant. For lentiviral infection, bladder organoids were incubated in lentivirus-containing medium with 8 μg/ml polybrene (Sigma) for 12 hr at 37°C. For selection of infected organoids, GFP- or mCherry-positive organoids were picked up from Matrigel under a fluorescence microscope.
Orthotopic transplantation of bladder tumors
Request a detailed protocolBladder tumors were dissociated into single cells as described above. Cells were resuspended in 80 μl DMEM containing 50% Matrigel (BD Bioscience) and then submucosally injected with 29-gauge insulin syringes into the anterior aspect of the bladder dome. The abdominal incisions and skin were then closed with a 4–0 nylon suture, and the surgical site was disinfected with alcohol. For orthotopic transplantation of tumor organoids, bladder tumor organoids were selected, resuspended in 50% organoid medium and 50% Matrigel and then transplanted into recipient mice.
Human bladder tumor samples and cancer cell lines
Request a detailed protocolFrozen human bladder tissue samples were obtained from the tissue bank of Seoul National University Hospital. For fresh bladder tumor samples, 0.5–1 cm3 specimens of fresh bladder tissue were obtained from patients undergoing cystectomy or TURB under a protocol approved by the SNUH Institutional Review Board (IRB Number: 1607-135-777). Informed consent and consent to publish was obtained from the patients. The cancer tissues were evaluated before being transported to POSTECH for further analysis. For experiments involving bladder cancer cell lines, J82 (RRID:CVCL_0359), T24 (RRID:CVCL_0554) and TCC (RRID:CVCL_1738) were used. All cell lines were authenticated using STR profiling method and were tested negative for mycoplasma contamination.
Quantitative RT-PCR
Request a detailed protocolHuman or mouse bladder samples were snap frozen in liquid nitrogen, then homogenized with a mortar and pestle, and RNA was extracted using the RNeasy Plus Mini Kit (Qiagen). The RNA samples were subsequently dissolved in RNase-free water, and their concentration and purity were determined with a spectrophotometer. The TAE/formamide electrophoresis method (Masek et al., 2005) was used for the analysis of RNA quality. For quantitative RT-PCR of mRNA transcripts, first-strand cDNA was synthesized using a High-Capacity cDNA reverse Transcriptase kit (Applied Biosystems) with oligo dT. Quantitative RT-PCR was performed using SYBR Green Supermix (Applied Biosystems) and a One-step cycler (Applied Biosystems), and gene expression was normalized to the housekeeping gene HPRT1.
Histological analysis
Request a detailed protocolTumor specimens were prefixed in 10% neutral-buffered formalin for 12 hr, then embedded in paraffin and sectioned into 4-µm-thick sections using a microtome. The slides were stained with hematoxylin and then counter-stained with eosin for histological analysis. For immunostaining, tumor specimens were embedded in OCT compound (Tissue-Tek) and sectioned into 10-µm-thick sections with a cryostat (Leica).
Immunofluorescence analysis of tissue sections
Request a detailed protocolBladder tumors were dissected from mice, fixed in 10% neutral-buffered for 3 hr, washed three times in PBS, incubated in 30% sucrose overnight, and embedded in OCT compound (Tissue-Tek). The slides were subsequently washed twice in PBS, blocked in 2% goat serum with 3% BSA in PBS containing 0.25% Triton X-100 for 1 hr, and incubated overnight at 4°C in a humidified chamber with primary antibodies diluted in blocking solution. Sections were washed three times with PBS containing 0.25% Triton X-100 and incubated for 1 hr at RT with the appropriate Alexa Fluor-conjugated secondary antibodies diluted 1:1000 in blocking solution. The slides were next washed three times with PBS and tissue sections were mounted with Prolong Gold mounting reagent (Invitrogen). All immunofluorescence images were analyzed by confocal microscopy (Leica SP5 or Olympus FV1000).
RNA-Seq library construction
Request a detailed protocolTotal RNA was extracted with TRIzol reagent (Thermo Fisher) according to the manufacturer’s instructions. RNA-seq libraries were constructed using the TruSeq sample Prep Kit V2 (Illumina). Quantity of RNA-seq library was determined by Nanodrop, and average quantity of RNA-seq libraries ranged from 30 to 50 ng/µl. RNA-seq libraries were sequenced using a NextSeq platform with single-end reads of 75 bases.
Differential gene expression and gene set enrichment analysis (GSEA) of RNA-seq data
Request a detailed protocolRaw reads from fastq files were aligned to the reference mouse genome assembly mm10 using Tophat allowing no mismatch (Trapnell et al., 2009). Differentially expressed genes were analyzed using Cufflinks tools (Trapnell et al., 2012). From all annotated genes, genes were removed if the average of rpkm across all sequenced samples is below 1.0, likely to have low depth to assign the genes. Gene set enrichment analyses were performed according to the instructions (RRID:SCR_003199). To generate a custom gene set for each luminal and basal marker, representative genes for each signature were obtained from previous study (Damrauer et al., 2014). The RNA-seq datasets used in the study have been deposited in NCBI GEO (Accession number: GSE129441).
Data analysis
Request a detailed protocolStatistical analysis was performed using GraphPad Prism software v.6. (RRID:SCR_015807). All data are presented as the mean ± SEM, and two group comparisons were conducted with a two-tailed Student’s test. A value of p<0.05 was considered statistically significant. For the analysis of TCGA data, the gene expression levels of muscle-invasive bladder cancer patients were downloaded from the TCGA data portal (https://portal.gdc.cancer.gov/). The FPKM expression values were log2(x + 1) transformed for convenient comparison of mRNA abundance estimates, where x denotes the FPKM value for each gene. The log-transformed expression values were normalized to z-scores for further analysis. Gene Cluster 3.0 was used for unsupervised hierarchical clustering (de Hoon et al., 2004), with default settings of an uncentered correlation and centroid linkage for the similarity metric and clustering method, respectively. Visualization of the mRNA cluster results was conducted using Java TreeView (Saldanha, 2004) (RRID:SCR_016916). To examine the clinical output of different mRNA clusters, survival analysis was conducted using the Oasis2 tool (Han et al., 2016). For the Kaplan-Meier survival plot, patients with overall survival of 5 years or less were considered for survival analysis. The Oncoprint format of mutation occurrences was plotted from cBioPortal (Cerami et al., 2012; Gao et al., 2013) (RRID:SCR_014555).
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 1, 3, 4, 5, 6 and 7.
-
NCBI Gene Expression OmnibusID GSE129441. Epigenetic regulation of mammalian Hedgehog signaling to the stroma determines the molecular subtype of human bladder cancer.
References
-
Established cell line of urinary bladder carcinoma (T24) containing tumour-specific antigenInternational Journal of Cancer 11:765–773.https://doi.org/10.1002/ijc.2910110327
-
Open source clustering softwareBioinformatics 20:1453–1454.https://doi.org/10.1093/bioinformatics/bth078
-
Epigenetics in cancer: targeting chromatin modificationsMolecular Cancer Therapeutics 8:1409–1420.https://doi.org/10.1158/1535-7163.MCT-08-0860
-
Gene silencing in cancer in association with promoter hypermethylationNew England Journal of Medicine 349:2042–2054.https://doi.org/10.1056/NEJMra023075
-
New cast for a new era: preclinical cancer drug development revisitedJournal of Clinical Investigation 123:3639–3645.https://doi.org/10.1172/JCI68340
-
Stromal contribution to the colorectal cancer transcriptomeNature Genetics 47:312–319.https://doi.org/10.1038/ng.3224
-
The fundamental role of epigenetic events in cancerNature Reviews Genetics 3:415–428.https://doi.org/10.1038/nrg816
-
Phase II, open-label, randomized, multicenter study comparing the efficacy and safety of olaparib, a poly (ADP-ribose) polymerase inhibitor, and pegylated liposomal doxorubicin in patients with BRCA1 or BRCA2 mutations and recurrent ovarian cancerJournal of Clinical Oncology 30:372–379.https://doi.org/10.1200/JCO.2011.36.9215
-
MethPrimer: designing primers for methylation PCRsBioinformatics 18:1427–1431.https://doi.org/10.1093/bioinformatics/18.11.1427
-
Stromal cells in tumor microenvironment and breast cancerCancer and Metastasis Reviews 32:303–315.https://doi.org/10.1007/s10555-012-9415-3
-
Denaturing RNA electrophoresis in TAE agarose gelsAnalytical Biochemistry 336:46–50.https://doi.org/10.1016/j.ab.2004.09.010
-
A cell line from an anaplastic transitional cell carcinoma of human urinary bladderBritish Journal of Cancer 35:142–151.https://doi.org/10.1038/bjc.1977.21
-
Java treeview--extensible visualization of microarray dataBioinformatics 20:3246–3248.https://doi.org/10.1093/bioinformatics/bth349
-
Efficacy and safety of vismodegib in advanced basal-cell carcinomaNew England Journal of Medicine 366:2171–2179.https://doi.org/10.1056/NEJMoa1113713
-
FK506 activates BMPR2, rescues endothelial dysfunction, and reverses pulmonary hypertensionThe Journal of Clinical Investigation 123:3600–3613.https://doi.org/10.1172/JCI65592
-
Inhibiting the hedgehog pathway in patients with the basal-cell nevus syndromeNew England Journal of Medicine 366:2180–2188.https://doi.org/10.1056/NEJMoa1113538
-
Hedgehog beyond medulloblastoma and basal cell carcinomaBiochimica Et Biophysica Acta (BBA) - Reviews on Cancer 1805:181–208.https://doi.org/10.1016/j.bbcan.2010.01.003
-
TopHat: discovering splice junctions with RNA-SeqBioinformatics 25:1105–1111.https://doi.org/10.1093/bioinformatics/btp120
-
Stromal hedgehog signaling maintains smooth muscle and hampers micro-invasive prostate cancerDisease Models & Mechanisms 10:39–52.https://doi.org/10.1242/dmm.027417
Article and author information
Author details
Funding
National Research Foundation of Korea (NRF-2017R1A2B4006043)
- Kunyoo Shin
National Research Foundation of Korea (NRF-2017M3C7A1047875)
- Kunyoo Shin
National Research Foundation of Korea (NRF-2017R1A5A1015366)
- Kunyoo Shin
National Research Foundation of Korea (Brain Korea 21 Program for Leading Universities & Students)
- Kunyoo Shin
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Phil Beachy for helpful discussion regarding this manuscript. This research was supported by grants from the National Research Foundation of Korea to KS: NRF-2017R1A2B4006043, NRF-2017M3C7A1047875, NRF-2017R1A5A1015366 and the BK21 Plus Program.
Ethics
Human subjects: Frozen human bladder tissue samples were obtained from the tissue bank of Seoul National University Hospital. For fresh bladder tumor samples, 0.5-1 cm3 specimens of fresh bladder tissue were obtained from patients undergoing cystectomy or TURB under a protocol approved by the SNUH Institutional Review Board (IRB Number: 1607-135-777). Informed consent and consent to publish was obtained from the patients.
Animal experimentation: Mouse procedures were performed under isoflurane anesthesia with a standard vaporizer. All procedures were performed under a protocol approved by the Institutional Animal Care and Use Committee (IACUC) at POSTECH (POSTECH-2017-0094).
Copyright
© 2019, Kim 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
-
- 3,522
- views
-
- 546
- downloads
-
- 17
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
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
- Evolutionary Biology
In growing cell populations such as tumours, mutations can serve as markers that allow tracking the past evolution from current samples. The genomic analyses of bulk samples and samples from multiple regions have shed light on the evolutionary forces acting on tumours. However, little is known empirically on the spatio-temporal dynamics of tumour evolution. Here, we leverage published data from resected hepatocellular carcinomas, each with several hundred samples taken in two and three dimensions. Using spatial metrics of evolution, we find that tumour cells grow predominantly uniformly within the tumour volume instead of at the surface. We determine how mutations and cells are dispersed throughout the tumour and how cell death contributes to the overall tumour growth. Our methods shed light on the early evolution of tumours in vivo and can be applied to high-resolution data in the emerging field of spatial biology.
-
- Cancer Biology
- Evolutionary Biology
In asexual populations that don’t undergo recombination, such as cancer, deleterious mutations are expected to accrue readily due to genome-wide linkage between mutations. Despite this mutational load of often thousands of deleterious mutations, many tumors thrive. How tumors survive the damaging consequences of this mutational load is not well understood. Here, we investigate the functional consequences of mutational load in 10,295 human tumors by quantifying their phenotypic response through changes in gene expression. Using a generalized linear mixed model (GLMM), we find that high mutational load tumors up-regulate proteostasis machinery related to the mitigation and prevention of protein misfolding. We replicate these expression responses in cancer cell lines and show that the viability in high mutational load cancer cells is strongly dependent on complexes that degrade and refold proteins. This indicates that the upregulation of proteostasis machinery is causally important for high mutational burden tumors and uncovers new therapeutic vulnerabilities.