- Open Access
HOXC8 regulates self-renewal, differentiation and transformation of breast cancer stem cells
© The Author(s). 2017
- Received: 27 October 2016
- Accepted: 23 January 2017
- Published: 16 February 2017
Homeobox genes are master regulators of cell fate during embryonic development and their expression is altered in cancer. By regulating the balance between cell proliferation and differentiation, they maintain homeostasis of normal tissues. Here, we screened the expression of homeobox genes in mammary stem cells to establish their role in stem cells transformation in breast cancer.
Using a Homeobox Genes PCR array, we screened 83 homeobox genes in normal cancer breast stem/progenitor cells isolated by flow cytometry. The candidate gene HOXC8 epigenetic regulation was studied by DNA methylation and miRNA expression analyses. Self-renewal and differentiation of HOXC8-overexpressing or knockdown cells were assessed by flow cytometry and mammosphere, 3D matrigel and soft agar assays. Clinical relevance of in vitro findings were validated by bioinformatics analysis of patient datasets from TCGA and METABRIC studies.
In this study we demonstrate altered expression of homeobox genes in breast cancer stem/progenitor cells. HOXC8 was consistently downregulated in stem/progenitor cells of all breast molecular subtypes, thus representing an interesting tumour suppressor candidate. We show that downregulated expression of HOXC8 is associated with DNA methylation at the gene promoter and expression of miR196 family members. Functional studies demonstrated that HOXC8 gain of function induces a decrease in the CD44+/CD24-/low cancer stem cell population and proportion of chemoresistant cells, with a concomitant increase in CD24+ differentiated cells. Increased HOXC8 levels also decrease the ability of cancer cells to form mammospheres and to grow in anchorage-independent conditions. Furthermore, loss of HOXC8 in non-tumorigenic mammary epithelial cells expands the cancer stem/progenitor cells pool, increases stem cell self-renewal, prevents differentiation induced by retinoic acid and induces a transformed phenotype.
Taken together, our study points to an important role of homeobox genes in breast cancer stem/progenitor cell function and establishes HOXC8 as a suppressor of stemness and transformation in the mammary gland lineage.
- Breast cancer
- Cancer stem cells
- Homeobox genes
Cancer stem/progenitor cell cells (CSC) are considered to play an important role in breast cancer complexity as they retain the fundamental features of normal stem cells, being able to self-renew and differentiate into the different cell types that comprise the heterogeneous tumour mass [1, 2]. Breast CSC have a degree of cellular plasticity and it is unclear whether they originate from transformation of normal mammary stem cells (MaSC) or by reprogramming of differentiated cells to a malignant phenotype .
Genes that are required for tissue development and normal stem cell plasticity may contribute to the molecular blueprint of breast CSC and tumorigenesis. Indeed, it is well established that many molecular pathways involved in normal development are altered during cancer progression . Amongst the developmental pathways, homeobox genes function as “master” gene regulators. They are characterised by a homeodomain that binds specific DNA elements and encode transcription factors involved in the regulation of cell growth and differentiation, thus establishing and maintaining cell identity and fate during development . Given their role, it has been postulated that a balanced expression of these genes is critical for stem cell function, especially for self-renewal, differentiation and quiescence . Within homeobox genes, the HOX gene family is important in early tissue differentiation. For instance, the HOXA gene cluster is silenced in pluripotent human embryonic stem cells and epigenetically reprogrammed and expressed when cells are induced to differentiate . They also regulate neural and hematopoietic stem cells function [6, 7].
Homeobox genes are expressed during mammary gland development and they appear to be regulated both by hormones and extracellular matrix remodelling . Examples include Hoxc6 and Msx2 being down-regulated and up-regulated by estrogen and progesterone, respectively. Msx2 is also regulated by the interaction between epithelial and stromal cells in the gland . Several homeobox genes are involved in proliferation and differentiation of mammary stem/progenitor cells. For instance Prrx1, Six1, Lbx1, Sox9, Msx, Zeb1, Dlx4 have been shown to regulate epithelial-to-mesenchymal transition (EMT) in mammary cells, thus promoting stem cell-like characteristics [9–15].
Aberrant expression of homeobox genes has been reported in several malignancies . Mis-expression of homeobox genes can lead to abnormal differentiation and proliferation, leading to a change in cell identity or homeotic transformation, therefore playing an important role in carcinogenesis . In cancer, homeobox genes function as “tumour modulators” as their deregulation normally involve either up-regulation of genes expressed in undifferentiated cells or down-regulation of genes expressed in differentiated tissue, thus acting either as oncogenes or tumour suppressor genes . Abnormal expression of homeobox genes leading to loss of differentiation is observed in breast cancer where expression of Hoxb6, Hoxb7, Hoxc6, Hoxc8, Hoxd4, Hoxd8, Hoxd9, Hoxd10, Hoxa5 is lost in mouse mammary tumours and Hoxa1, Hoxd3, Hoxd12, Msx1, Six1 and Oct3 are instead expressed in cancer, but absent or lowly expressed in the normal differentiated gland . This deregulated gene expression has been shown to be involved in neoplastic transformation by regulating cell cycle, apoptosis, angiogenesis, metastasis and cell adhesion . Epigenetic mechanisms, such as DNA methylation, histone modification and silencing by non-coding RNA are involved in the regulation of homeobox gene expression . Moreover, epigenetic silencing of these genes is an early event in breast carcinogenesis [18–20] and DNA methylation of homeobox genes is associated with specific breast cancer cell types and cancer molecular subtypes [21–23].
In this study, we screened the expression of homeobox genes in breast CSC and demonstrate that a large number of genes are differentially expressed in immortalised CSC compared to normal MaSC, in agreement with evidence suggesting an involvement of developmentally-regulated genes in early stages of cancer transformation. We also reported expression of homebox genes uniquely altered in different breast cancer subtypes, which could represent novel CSC biomarkers for patient stratification. Of the genes commonly altered across different breast cancer molecular subtypes, HOXC8 was chosen as a putative novel tumour suppressor gene. We show that HOXC8 downregulation in breast CSC is associated with epigenetic silencing. Downregulation of HOXC8 in breast cancer was also found by meta-analysis of breast cancer data from large cohort studies. Gain of function of HOXC8 reduced CSC self-renewal and the ability of cancer cells to grow in anchorage-independent conditions. Conversely, loss of function of HOXC8 in mammary normal cells induced CSC proliferation and colony formation. In addition, reduced HOXC8 expression impaired cell differentiation and response to retinoic acid. Taken together, our study shows that homeobox genes represent novel biomarkers of breast CSC and that HOXC8 functions as a novel tumour suppressor gene by regulating breast CSC proliferation and differentiation.
Cells and materials
Human mammary epithelial cells (HMEC) were obtained from Invitrogen and cultured with the proprietary HuMEC Ready Medium. The immortalised human mammary epithelial cell line MCF10A (ATCC) was cultured in HuMEC Ready Medium supplemented with 100 ng/ml cholera toxin. Breast cancer cell lines MCF-7, HCC1954, HCC1428 (all from ATCC), MDA-MB-468, MDA-MB-231, BT549 and Hs578T (all from the NCI-60 cell collection, CRN cell bank, University of Nottingham) were grown in RPMI medium supplemented with 10% foetal calf serum (FCS), 1% Penicillin/Streptomycin (Pen/Strep), 1% L-Glutamine, 1% sodium pyruvate, 1% non-essential amino acids (NEAA). HEK 293 T (ATCC) cells were grown in DMEM medium containing 10% FCS supplemented, 1% Pen/Strep, 1% L-Glutamine, 1% sodium pyruvate, 1% NEAA. Cell lines obtained from ATCC were used within few passages from the original stocks, whereas other cell lines were authenticated by STR profiling (Eurofins Genomics, Germany) using the Promega PowerPlex 21 PCR kit and matched against the ATCC STR database. All cell lines were tested form mycoplasma contamination using the EZ-PCR Mycoplasma Test Kit (Geneflow). All cell culture materials were from Invitrogen and chemicals from Sigma-Aldrich, unless otherwise stated.
For chemoresistance assay, cells were incubated with Paclitaxel (10nM), doxorubicin (500 nM) and 5-Fluorouracil (5-FU) (500 μM) for 5 days. Viable cells were trypsinised, collected and the CD44+/CD24-/low profile was analysed by flow cytometry.
For differentiation experiments, cells were treated with 1 μM all-trans-retinoic acid (ATRA) for 7 days.
Stem/progenitor cell populations were identified and isolated as CD44+/CD24-/low by FACS (Beckman Coulter MoFlo XDP sorter). The following antibodies and isotype controls were used: FITC-anti-CD44 (eBioscience 11-0441, 1:100 dilution), PE-anti-CD24 (eBioscience 12-0247, 1:20 dilution), APC-anti-CD44 (eBioscience 17-0441, 1:167 dilution), FITC-anti-CD24 (eBioscience 11-0247, 1:20 dilution), FITC-rat isotype control (eBioscience 11-4031, 1:100 dilution), PE-mouse isotype control (eBioscience 9012-4714, 1:20 dilution), APC-rat isotype control (eBioscience 17-4031, 1:167 dilution), FITC-mouse isotype control (eBioscience 11-4714, 1:20 dilution). Briefly, 5x105 cells were incubated with antibodies in the dark at 4 °C for 1 h. Cells were washed and re-suspended in complete medium for analysis. For cell sorting, 106 cells in 100 μL of complete medium were stained as described above. Cells were then treated with DNase I (Qiagen) in RDD buffer (10 mM Tris-HCL, 2.5 mM MgCl2, 0.5 mM CaCl2, pH7.6) for at least 15 min at room temperature, and passed through a 40 μm mesh (Becton Dickinson LTD) to avoid cell clumping. Data were analysed with the Weasel software.
Homeobox gene expression array
Total RNA was extracted with the RNeasy Mini Kit (Qiagen) and cDNA transcribed with the RT2 First Strand Kit (Qiagen). Expression of homeobox genes was analysed by real-time PCR (qRT-PCR) using the Homeobox (HOX) Genes RT2 Profiler PCR Array (Qiagen) on a Roche LightCycler® 480 System. Real-time PCR data analysis was performed by using the RT2 Profiler™ PCR Array Data Analysis (Qiagen) and based on the ΔΔCT method with normalization of the raw data to the housekeeping gene RPLP0 after analysis with the BestKeeper software .
Gene expression was analysed by qRT-PCR using TaqMan® chemistry (Applied Biosystems). Data were analysed by the ΔΔCT method with normalization to the housekeeping genes RPLP0. HOXC8 copy number was measured by qRT-PCR using genomic DNA (extracted with DNeasy Blood &Tissue kit, Qiagen) and LuminoCt® SYBR® Green qPCR ReadyMix. miRNAs were purified from cell pellets using the miRNeasy Mini Kit (Qiagen) and reverse transcribed with miScript II RT kit (Qiagen). qRT-PCR was performed using the LuminoCt® SYBR® Green qPCR ReadyMix. Data were analysed by the ΔΔCT method with normalization to the endogenous control RNU6B. For primers and assay used see Additional file 1: Table S1.
Genomic DNA was isolated using DNeasy Tissue kit (Qiagen). Bisulfite conversion (200 ng genomic DNA) was achieved using the EZ DNA methylation kit (Zymo Research) and 3 μl of bisulfite converted DNA was used for PCR reaction using the HotStarTaq masterm mix kit (Qiagen). Primers spanning the CpG island sequence were designed using the Methprimer software (http://www.urogene.org/cgi-bin/methprimer/methprimer.cgi) (Additional file 1: Table S1). Bisulfire converted and purified PCR products were directly sequenced.
Lentiviral plasmids for HOXC8 overexpression and knockdown
The HOXC8 plasmid was prepared by amplifying HOXC8 cDNA from HMEC cells using Phusion® high-fidelity DNA polymerase PCR master mix (New England Biolabs) and cloning between EcoRI and SpeI sites of a modified pSIN-EF2-Puro plasmid (pSin-EF2-Nanog-Pur was a gift from James Thomson, Addgene plasmid # 16578). HOXC8 and scrambled shRNA primers were constructed using previously reported sequences . shRNA harpins were ligated into AgeI/EcoRI-digested pLKO-Tet-On plasmid (Tet-pLKO-puro was a gift from Dmitri Wiederschain, Addgene plasmid # 21915). Plasmids were transformed in Stabl3 cells (Invitrogen) and purified using Qiagen plasmid purification kits (Qiagen). Viral particles were produced by transfecting HEK 293 T cells with shRNA or pSIN plasmids, together with psPAX2 and pMD2.G plamids (gift from Didier Trono, Addgene plasmids # 12260 and 12259) using FuGENE®6 (Promega). Virus-containing supernatant was collected at 72 h post-transfection, filtered through 0.45 μm filter and added of 6 μg/mL hexadimethrine bromide (polybrene) (Sigma-Aldrich). Stable transgenic cells were selected after 72 h using Puromycin (0.75–1 μg/ml). Tet-inducible shRNA was induced by treatment with 2ug/ml doxycycline.
Nuclear proteins were extracted with NucBuster™ Protein Extraction Kit (Calbiochem). Extracted proteins were loaded into a 12% or 15% Acrylamide gel (20 μg/lane), separated by SDS-PAGE electrophoresis and blotted onto a PVDF membrane. Membranes were blocked with 5% skimmed milk and then probed overnight at 4 °C with a rabbit anti-HOXC8 antibody (1:1,000, Sigma H1791) in the presence of 0.1% Tween 20 and 5% milk. Membranes were stripped and the incubated with a mouse anti-Lamin A/C antibody (1: 2,000, Cell Signalling 4777) overnight at 4 °C. Peroxidase conjugated donkey anti-rabbit and sheep anti-mouse (1:10,000; GE Healthcare NA934 and NA931, respectively) antibodies were incubated for 1 h at RT. ECL prime kit (GE Healthcare) was used to detect chemiluminescence.
Single cells were plated in ultra-low adherent flasks coated with Poly 2-hydroxyethyl methacrylate (PolyHEMA) (Sigma-Aldrich) at a density of 2x104 cells/mL. Cells were grown in mammosphere medium consisting of Dulbecco’s Modified Eagle Medium: Nutrient Mixture F-12 (DMEM/F12) supplemented with 20 ng/mL epidermal growth factor (EGF, R&D Systems) 20 ng/mL basic fibroblast growth factor (bFGF) (R&D Systems), 5 μg/mL insulin, and 1% pen/strep, B27, 4 μg/mL heparin, and 5 μg/mL hydrocortisone. Cells were incubated at 37 °C and 5% CO2 for 7 days to attain first generation mammospheres. Formed mammospheres were collected by gentle centrifugation at 200 × g for 4 min, dissociated enzymatically with trypsin and mechanically by pipetting. Serial passaging to attain consequent generations of mammospheres (up to four) was done by re-plating dissociated mammosphere-derived cells at a density of 5x103 cells/ml and culture for 7–10 days in mammosphere medium. For PKH26 labelled mammospheres, 2x106 cells were stained with 1 μM PKH26 for 5 min, washed and then and plated in mammosphere medium for the first generation. After 7 days, mammospheres were collected, dissociated as described above and sorted to isolate the PKH26 positivehigh population. Sorted cells were cultured again in mammosphere medium to obtain the second generation.
Soft agar assay
Cancer cells (1-3x104 cells/6 well) were seeded in 1 ml of 0.5% noble agar in complete RPMI medium overlaying 2 ml 1% agar in the same medium. MCF10A cells were instead seeded in HuMEC medium, as described above. After 2 weeks culture, cell colonies were stained with 0.05% crystal violet in 10% ethanol/PBS for 1 h and colonies ≥ of 100 μm counted under a MZ125 Leica stereomicroscope.
3D Matrigel assay
The 3D on-top Matrigel assay was performed according to Lee et al. . Cells (1.5x104) were plated into a 4 well plate covered with 120 μl of growth factor reduced Matrigel (BD Biosciences) and cells allowed to attached for 30 min. Cells were then cultured in 500 μl HuMEC Ready Medium supplemented with 100 ng/ml cholera toxin and 10% Matrigel for up to 7 days, with medium changed every other day.
Bioinformatics and statistics
For the analysis of the Homeobox gene array, hierarchical cluster analysis was performed by using the software ArrayMining (http://www.arraymining.net/). Venn diagrams were generated with the online tool Venny (http://bioinfogp.cnb.csic.es/tools/venny/index.html).
The analysis of copy number variation in clinical datasets was performed with the software OASIS (OASIS: A Web-based Platform for Exploratory Analysis of Cancer Genome and Transcriptome data, www.oasis-genomics.org). Expression analysis of homeobox genes in TCGA (RNASeq and PAM50) and METABRIC patient datasets was performed using the UCSC Cancer Browser (https://genome-cancer.ucsc.edu/) and the OASIS portal, respectively [27, 28]. For the analysis of HOXC8 in the TCGA dataset, clinical parameters and next-generation RNA sequencing data (RNAseqv2) were obtained from TCGA ccRCC and normal non-tumour breast tissue sample set [29, 30]. Patients were sub-divided depending on clinic-pathological parameters provided by TCGA. Normalized gene counts (gene counts rescaled according to library size) from each patient were compiled into a tab-delimited file for downstream analysis by RobinA implementation of the Bioconductor-edgeR software [31, 32]. Expression analysis of HOXC8 in the METABRIC dataset [33, 34] was also obtained from the OASIS portal. Methylation analysis of the TCGA dataset was performed by using the UCSC Cancer Browser and the MethHC browser (http://methhc.mbc.nctu.edu.tw/php/index.php) .
All data are expressed as mean ± standard deviation. Comparisons between two samples were done using Student’s t-test. One-way or Two-way analysis of variance (ANOVA) analyses followed by Bonferroni’s multiple comparison tests were used for multiple group comparisons. Statistical analyses were performed with GraphPad Prism with significance levels set at *P < 0.05, **P < 0.01, ***P < 0.001.
Expression of homeobox genes in breast CSC
In this study we investigated the expression of 83 different homeobox genes involved in morphogenesis, body pattern formation, embryonic development, and cellular differentiation in normal MaSC and breast CSC. We used cell lines representing normal mammary epithelium, immortalised mammary epithelium (IMM) and the breast cancer molecular subtypes luminal (LUM: ER+, PR+, HER2−), HER2-enriched (HER2-E: ER−, PR−, HER2+) and basal/triple negative (TN: ER−, PR−, HER2−). CD44+/CD24-/low stem cell/progenitor cells were isolated from normal and immortalised mammary cells (HMEC and MCF10A, respectively), LUM (MCF-7, HCC1428), HER2-E (MDA-MB-468, HCC1954) and TN (BT549, Hs578T and MDA-MB-231) cancer cells. Expression of the majority (77.11%) of the homeobox genes analysed was altered in CSC compared to MaSC (CD44+/CD24-/low in HMEC) across different cancer subtypes (≥2-fold change, with genes considered differentially expressed only if either up-regulated or down-regulated in at least 2 cell lines of the same molecular subtype).
A large number of genes were differentially expressed in immortalised stem cells, suggesting altered expression of homeobox genes is an early event during CSC transformation. The largest number of differentially expressed genes was observed in TN CSC. Differential expression of a large number of genes was consistent across cell lines of the same molecular subtypes, representing 54.8% of the genes in LUM, 41.8% in HER2-E, and 72.2% in TN. A minority of genes changed expression in either only one cell line within molecular subtypes or in different directions (either down- or upregulated). The majority of genes in the immortalised stem cells were downregulated, whereas upregulated genes were predominant in CSC across the different molecular subtypes (Additional file 2: Figure S1A). Of the downregulated genes in the immortal stem cells, the majority were shared by TN CSC, whereas upregulated genes were shared similarly by all CSC subtypes. Interestingly, the TN CSC showed the largest number of uniquely downregulated and upregulated homeobox genes (Additional file 2: Figure S1B, Figure S1C). When considering the cancer subtype-specific homeobox genes, three genes were unique to LUM CSC, 6 to HER2-E CSC, and 21 to TN CSC. Nine genes were instead common to all cancer subtypes (considering at least 2 cell lines) (Fig. 1b). Of these, six genes were either consistently upregulated (BARX1, DLX4, HOXB8, PITX2, SIX2) or downregulated (HOXC8).
For the luminal-specific genes, HOXA1 was found downregulated in all datasets, EN2 upregulated in the METABRIC, and HOXB1 expression showed no overlap. The genes in the HER2-E subtype corresponded well with patient data, with HOXC10, HOXC12, PITX1, SIX4 matching all datasets and HESX1 and PDHX the METABRIC and TCGA sets, respectively. Finally, In the TN-specific gene set, DLX5, HHEX, HOXA9 expression corresponded to all datasets, whereas HOXB9, HOXC6, ISL1, LBX2, VAX1 only to two datasets (TCGA/PAM50 or TCGA/METABRIC or METABRIC/PAM50). Interestingly, of the genes that overlapped with only one dataset (CDX2, DLX1, DLX6, HOXB2, OTP, PROP1), three matched with their PAM50 subtype (CDX2, DLX6 and HOXB2). However, the expression of a number of gene did not correspond to any dataset (DLX2, HOXC9, HOXD3, MEIS1, MEOX1, NKX3-1, PITX3), suggesting that the expression of these genes might possibly be specific to the particular cell lines analysed in this study.
HOXC8 is downregulated in breast CSC and tumour tissues
HOXC8 downregulation in breast CSC is associated with epigenetic silencing
Therefore we focused our attention on the epigenetic regulation of HOXC8. Expression of miR196 has been previously linked to regulation of HOXC8  and therefore we measured the expression of two miR196 family members in HER2-E and TN CSC, as these showed the most downregulated expression of HOXC8. The expression of miR196-a and miR196-b was significantly increased in CSC isolated from MDA-MB-231, but not in CSC isolated from other cell lines (Fig. 4b). Direct bisulphite sequencing analysis of the CpG island in the HOXC8 promoter upstream of the transcription start site revealed high levels of DNA methylation in all CSC compared to normal MaSC (Fig. 4c and Additional file 5: Figure S4). A significant increase in DNA methylation at the promoter region of HOXC8 was also observed in the breast TCGA tumour samples compared to normal tissue (Fig. 4d), thus demonstrating an important role of epigenetic modifications in the regulation of HOXC8 expression and the clinical relevance of the in vitro findings.
Gain of function of HOXC8 reduces CSC self-renewal and chemoresistance
Loss of function of HOXC8 increases stem cell self-renewal and transformation
Expression of HOXC8 is required for breast stem cell differentiation
The importance of developmental genes in the regulation of cancer stem cell self-renewal and plasticity is well established [1, 40]. This study investigated for the first time the expression of a panel of developmental regulatorsof the homeobox family in breast cancer stem/progenitor cells. We found that homeobox genes are deregulated at an early stage of stem cell transformation as well in different cancer molecular subtypes. We found that the majority of homeobox genes analysed were differentially expressed in transformed/non tumorigenic breast stem cells, with the majority being downregulated compared to normal MaSC. Epigenetic deregulation of homeobox genes is a frequent and early event in breast cancer. DNA methylation of the HOX cluster has been reported to be associated with mammary epithelial cell immortalisation, suggesting that modulation of early developmental genes plays a role in escaping cell senescence . In addition, alterations in DNA methylation of several homeobox genes is found in ductal carcinoma in situ (DCIS), and low-grade breast tumours  , with those found in DCIS being associated with disease progression . These findings suggest a model whereby stem cells could be the cell of origin of breast tumours, as developmental genes that are normally reversibly silenced by polycomb complexes in these cells  become permanently silenced via DNA methylation during carcinogenesis [40, 44, 45]. Our data support this hypothesis as the majority of downregulated homeobox genes in immortal stem cells were also found not to be expressed or expressed a low level in TN breast cancer, a molecular subtype particularly enriched in CSC. A number of homeobox genes were differentially expressed only in specific breast molecular subtypes. Particularly important is the finding that the highest number of unique differentially expressed homeobox genes were found in TN CSC as this molecular subtype of breast cancer is still poorly understood and does not respond to current conventional or targeted therapies [46, 47]. DNA methylation profile of different molecular subtypes identified high enrichment of homebox genes being hypermethylated specifically in ER− which were associated with high rate of tumour recurrence and therefore likely contributing to poor patient outcome . Expression of homeobox genes have been also found to be specific of TN breast cancer subtypes, with the MSL (mesenchymal-stem-like) subtype being particularly enriched . Therefore these newly identified CSC homeobox genes could represent novel biomarkers that could be used for patient stratification in precision medicine and warrant further investigation.
Among the differentially expressed genes identified in this study, HOXC8 was found consistently downregulated in CSC and therefore represented a novel putative tumour suppressor candidate. HOXC8 was especially downregulated in TN and HER2-E CSC and it was also found to be downregulated in clinical samples of TCGA and METABRIC datasets, strengthening its significance as a potential biomarker. In line with previous studies reporting the importance of epigenetic regulation of homeobox genes during development and carcinogenesis , we also found that HOXC8 silencing is associated with the expression of miR196 and DNA methylation at the gene promoter region in CSC. Regulation of HOXC8 via miRNA 196 has been reported in breast cancer whereby the ration of miR196 and HOXC8 correlates with cell migration and metastasis  and methylation of the HOXC8 promoter can induce gene silencing by inducing polymerase II stalling during transcription .
Our results demonstrate that HOXC8 is an important regulator of stem cell self-renewal and that HOXC8 can act as modulator of cell differentiation in the breast. Using a combination of strategies and cell culture models, we show that silencing of HOXC8 confers mammary cells with stemness potential and this can be reversed by its re-expression with induction of differentiation. The enhanced CSC population induced by downregulation of HOXC8 in mammary cells is also associated with an increase in transformed phenotype and this observation supports the notion that breast cell tumorigenicity is enhanced through acquisition of an undifferentiated state . HOX proteins are transcription factors that regulate cell proliferation and differentiation during embryonic development. In normal adult tissues, they maintain tissue homeostasis by defining cell fate and lineage commitment [4, 16, 17]. Several cancers demonstrate altered expression of HOX genes but their role in transformation of stem cells and establishment of CSC is not well-known. Recently, a role in CSC regulation has been strongly established for HOXA5 in breast and colorectal cancer, whereby HOXA5 can maintain homeostasis by suppressing stemness in these tissues [52, 53].
During development, HOXC8 is expressed in the neural tube and somatic mesoderm and its expression is essential for skeletal development [54, 55]. HOXC8 is also expressed in the normal mammary gland [8, 56, 57] and it participates in the initiation of mammary morphogenesis . Previous studies have reported an increased expression of HOXC8 during breast cancer progression that can induce metastasis through direct regulation of CDH11 and EMB genes [25, 59]. This different role of HOXC8 suggests that it is possible that this gene could act as a “modulator” of carcinogenesis and therefore play either a role as suppressor or driver of transformation depending on the physiological and pathological context. Consistent with this, HOXC8 has been implicated in different cancer types, including prostate, cervical, breast, oesophageal, and pancreatic cancer [36, 60–63]. In these cancers HOXC8 has been reported to act either as tumour suppressor gene or oncogene, suggesting that the function of HOXC8 may depend on the activation or repression of different HOXC8 targets whose expression can be tissue-specific. Several targets of HOXC8 have been characterised [57, 64] and their functional role in the context of HOXC8 regulation may provide new insights into the role of HOXC8 in different cancer types. Although an effective approach for inhibition of the interaction between HOX proteins and co-factors PBX is available , targeting specific HOX proteins still remains a challenge. Therefore, a better definition of HOX transcriptional regulators and HOX targets could offer new strategies for therapeutic intervention.
We have also shown that HOXC8 downregulation in mammary cells can lead to impaired differentiation in response to retinoic acid, and that the expression of HOXC8 can be induced by retinoids in normal breast cells. This role has also been recently demonstrated for HOXA5, suggesting a prominent role of HOX genes in the differentiation of the mammary gland . Retinoic acid has been shown to regulate breast CSC differentiation, reduce mammospheres formation and CSC frequency in combination with chemotherapy [66, 67]. Therefore our study adds to the evidence that HOX genes are important regulators of stem cell fate and that differentiation therapy should be considered as a viable approach to target breast CSC.
This study reports the screening of a large number of homeobox genes in breast CSC. We show that the expression of homeobox genes is altered in early-transformed breast cells and in different breast cancer molecular subtypes. Altered expression of homeobox genes was also supported by data obtained from large cohorts of patient through bioinformatics analysis. Therefore, this screening provides novel biomarkers that can be used in precision medicine and as possible targets for CSC-directed therapies. Of the altered genes, HOXC8 demonstrated to function as a possible novel tumour suppressor in breast CSC by regulating stem cell self-renewal, differentiation and transformation. Therefore, this study demonstrates the critical role of homeobox genes in stem cells and paves the way to larger functional studies to explore the networks involved in homeotic gene regulation in breast cancer.
The authors would like to thanks Adrian Robins, Nina Lane, David Onion, Nicola Croxall, and Scott Hulme for technical and FACS analysis support.
This study was supported by grants from the Biotechnology and Biological Sciences Research Council (BBSRC; grant number RS4949) Royal Society of London (grant number RA3810) and the University of Nottingham (grant number A2B090).
Availability of data and materials
The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.
MS performed experiments and contributed to the design of the study. RC performed miRNA analyses and assisted with bioinformatics analysis. BW performed gene knockdown experiments. JP critically reviewed the manuscript and assisted in study design. NM performed bioinformatics analyses and critically reviewed the manuscript. AG critically reviewed the manuscript and contributed to study design. CA performed experiments, designed the study and wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The study was approved by the University of Nottingham local ethical committee. Anonymised patient data from TGCA and METABRIC studies were available from the Genomic Data Commons Data Portal and OASIS Genomics Portal.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Shah M, Allegrucci C. Keeping an open mind: highlights and controversies of the breast cancer stem cell theory. Breast Cancer. 2012;4:155–66.PubMedPubMed CentralGoogle Scholar
- Valent P, Bonnet D, De Maria R, Lapidot T, Copland M, Melo JV, Chomienne C, Ishikawa F, Schuringa JJ, Stassi G, et al. Cancer stem cell definitions and terminology: the devil is in the details. Nat Rev Cancer. 2012;12:767–75.View ArticlePubMedGoogle Scholar
- Mark M, Rijli FM, Chambon P. Homeobox genes in embryogenesis and pathogenesis. Pediatr Res. 1997;42:421–9.View ArticlePubMedGoogle Scholar
- Shah N, Sukumar S. The Hox genes and their roles in oncogenesis. Nat Rev Cancer. 2010;10:361–71.View ArticlePubMedGoogle Scholar
- Atkinson SP, Koch CM, Clelland GK, Willcox S, Fowler JC, Stewart R, Lako M, Dunham I, Armstrong L. Epigenetic marking prepares the human HOXA cluster for activation during differentiation of pluripotent cells. Stem Cells. 2008;26:1174–85.View ArticlePubMedGoogle Scholar
- Barber BA, Rastegar M. Epigenetic control of Hox genes during neurogenesis, development, and disease. Ann Anat. 2010;192:261–74.View ArticlePubMedGoogle Scholar
- He H, Hua X, Yan J. Epigenetic regulations in hematopoietic Hox code. Oncogene. 2011;30:379–88.View ArticlePubMedGoogle Scholar
- Chen H, Sukumar S. Role of homeobox genes in normal mammary gland development and breast tumorigenesis. J Mammary Gland Biol Neoplasia. 2003;8:159–75.View ArticlePubMedGoogle Scholar
- Ocana OH, Corcoles R, Fabra A, Moreno-Bueno G, Acloque H, Vega S, Barrallo-Gimeno A, Cano A, Nieto MA. Metastatic colonization requires the repression of the epithelial-mesenchymal transition inducer Prrx1. Cancer Cell. 2012;22:709–24.View ArticlePubMedGoogle Scholar
- McCoy EL, Iwanaga R, Jedlicka P, Abbey NS, Chodosh LA, Heichman KA, Welm AL, Ford HL. Six1 expands the mouse mammary epithelial stem/progenitor cell pool and induces mammary tumors that undergo epithelial-mesenchymal transition. J Clin Invest. 2009;119:2663–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu M, Smolen GA, Zhang J, Wittner B, Schott BJ, Brachtel E, Ramaswamy S, Maheswaran S, Haber DA. A developmentally regulated inducer of EMT, LBX1, contributes to breast cancer progression. Genes Dev. 2009;23:1737–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Guo W, Keckesova Z, Donaher JL, Shibue T, Tischler V, Reinhardt F, Itzkovitz S, Noske A, Zurrer-Hardi U, Bell G, et al. Slug and Sox9 cooperatively determine the mammary stem cell state. Cell. 2012;148:1015–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Hu G, Lee H, Price SM, Shen MM, Abate-Shen C. Msx homeobox genes inhibit differentiation through upregulation of cyclin D1. Development. 2001;128:2373–84.PubMedGoogle Scholar
- Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, Farshid G, Vadas MA, Khew-Goodall Y, Goodall GJ. The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat Cell Biol. 2008;10:593–601.View ArticlePubMedGoogle Scholar
- Zhang L, Yang M, Gan L, He T, Xiao X, Stewart MD, Liu X, Yang L, Zhang T, Zhao Y, Fu J. DLX4 upregulates TWIST and enhances tumor migration, invasion and metastasis. Int J Biol Sci. 2012;8:1178–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Samuel S, Naora H. Homeobox gene expression in cancer: insights from developmental regulation and deregulation. Eur J Cancer. 2005;41:2428–37.View ArticlePubMedGoogle Scholar
- Abate-Shen C. Deregulated homeobox gene expression in cancer: cause or consequence? Nat Rev Cancer. 2002;2:777–85.View ArticlePubMedGoogle Scholar
- Tommasi S, Karm DL, Wu X, Yen Y, Pfeifer GP. Methylation of homeobox genes is a frequent and early epigenetic event in breast cancer. Breast Cancer Res. 2009;11:R14.View ArticlePubMedPubMed CentralGoogle Scholar
- Novak P, Jensen TJ, Garbe JC, Stampfer MR, Futscher BW. Stepwise DNA methylation changes are linked to escape from defined proliferation barriers and mammary epithelial cell immortalization. Cancer Res. 2009;69:5251–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Novak P, Jensen T, Oshiro MM, Watts GS, Kim CJ, Futscher BW. Agglomerative epigenetic aberrations are a common event in human breast cancer. Cancer Res. 2008;68:8616–25.View ArticlePubMedGoogle Scholar
- Bloushtain-Qimron N, Yao J, Snyder EL, Shipitsin M, Campbell LL, Mani SA, Hu M, Chen H, Ustyansky V, Antosiewicz JE, et al. Cell type-specific DNA methylation patterns in the human breast. Proc Natl Acad Sci U S A. 2008;105:14076–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Kamalakaran S, Varadan V, Giercksky Russnes HE, Levy D, Kendall J, Janevski A, Riggs M, Banerjee N, Synnestvedt M, Schlichting E, et al. DNA methylation patterns in luminal breast cancers differ from non-luminal subtypes and can identify relapse risk independent of other clinical variables. Mol Oncol. 2011;5:77–92.View ArticlePubMedGoogle Scholar
- Ronneberg JA, Fleischer T, Solvang HK, Nordgard SH, Edvardsen H, Potapenko I, Nebdal D, Daviaud C, Gut I, Bukholm I, et al. Methylation profiling with a panel of cancer related genes: association with estrogen receptor, TP53 mutation status and expression subtypes in sporadic breast cancer. Mol Oncol. 2011;5:61–76.View ArticlePubMedGoogle Scholar
- Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper--Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004;26:509–15.View ArticlePubMedGoogle Scholar
- Chao F, Zhang J, Zhang Y, Liu H, Yang C, Wang J, Guo Y, Wen X, Zhang K, Huang B, et al. Embigin, regulated by HOXC8, plays a suppressive role in breast tumorigenesis. Oncotarget. 2015;6:23496–509.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee GY, Kenny PA, Lee EH, Bissell MJ. Three-dimensional culture models of normal and malignant breast epithelial cells. Nat Methods. 2007;4:359–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Cline MS, Craft B, Swatloski T, Goldman M, Ma S, Haussler D, Zhu J. Exploring TCGA Pan-cancer data at the UCSC cancer genomics browser. Sci Rep. 2013;3:2652.View ArticlePubMedPubMed CentralGoogle Scholar
- Fernandez-Banet J, Esposito A, Coffin S, Horvath IB, Estrella H, Schefzick S, Deng S, Wang K, Aching K, Ding Y, et al. OASIS: web-based platform for exploring cancer multi-omics data. Nat Meth. 2016;13:9–10.View ArticleGoogle Scholar
- Ciriello G, Gatza ML, Beck AH, Wilkerson MD, Rhie SK, Pastore A, Zhang H, McLellan M, Yau C, Kandoth C, et al. Comprehensive molecular portraits of invasive lobular breast cancer. Cell. 2015;163:506–19.View ArticlePubMedPubMed CentralGoogle Scholar
- Cancer Genome Atlas Network. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490:61–70.View ArticleGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.View ArticlePubMedGoogle Scholar
- Lohse M, Bolger AM, Nagel A, Fernie AR, Lunn JE, Stitt M, Usadel B. RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic Acids Res. 2012;40:W622–627.View ArticlePubMedPubMed CentralGoogle Scholar
- Curtis C, Shah SP, Chin SF, Turashvili G, Rueda OM, Dunning MJ, Speed D, Lynch AG, Samarajiwa S, Yuan Y, et al. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012;486:346–52.PubMedPubMed CentralGoogle Scholar
- Pereira B, Chin SF, Rueda OM, Vollan HK, Provenzano E, Bardwell HA, Pugh M, Jones L, Russell R, Sammut SJ, et al. The somatic mutation profiles of 2,433 breast cancers refines their genomic and transcriptomic landscapes. Nat Commun. 2016;7:11479.View ArticlePubMedPubMed CentralGoogle Scholar
- Huang WY, Hsu SD, Huang HY, Sun YM, Chou CH, Weng SL, Huang HD. MethHC: a database of DNA methylation and gene expression in human cancer. Nucleic Acids Res. 2015;43:D856–861.View ArticlePubMedGoogle Scholar
- Li Y, Zhang M, Chen H, Dong Z, Ganapathy V, Thangaraju M, Huang S. Ratio of miR-196s to HOXC8 messenger RNA correlates with breast cancer cell migration and metastasis. Cancer Res. 2010;70:7894–904.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang YA, Shen K, Wang Y, Brooks SC. Retinoic acid signaling is required for proper morphogenesis of mammary gland. Dev Dyn. 2005;234:892–9.View ArticlePubMedGoogle Scholar
- Vermot J, Schuhbaur B, Le Mouellic H, McCaffery P, Garnier JM, Hentsch D, Brulet P, Niederreither K, Chambon P, Dolle P, Le Roux I. Retinaldehyde dehydrogenase 2 and Hoxc8 are required in the murine brachial spinal cord for the specification of Lim1+ motoneurons and the correct distribution of Islet1+ motoneurons. Development. 2005;132:1611–21.View ArticlePubMedGoogle Scholar
- Nadauld LD, Sandoval IT, Chidester S, Yost HJ, Jones DA. Adenomatous polyposis coli control of retinoic acid biosynthesis is critical for zebrafish intestinal development and differentiation. J Biol Chem. 2004;279:51581–9.View ArticlePubMedGoogle Scholar
- Shah M, Allegrucci C. Stem cell plasticity in development and cancer: epigenetic origin of cancer stem cells. Subcell Biochem. 2013;61:545–65.View ArticlePubMedGoogle Scholar
- Faryna M, Konermann C, Aulmann S, Bermejo JL, Brugger M, Diederichs S, Rom J, Weichenhan D, Claus R, Rehli M, et al. Genome-wide methylation screen in low-grade breast cancer identifies novel epigenetically altered genes as potential biomarkers for tumor diagnosis. FASEB J. 2012;26:4937–50.View ArticlePubMedGoogle Scholar
- Johnson KC, Koestler DC, Fleischer T, Chen P, Jenson EG, Marotti JD, Onega T, Kristensen VN, Christensen BC. DNA methylation in ductal carcinoma in situ related with future development of invasive breast cancer. Clin Epigenetics. 2015;7:75.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee TI, Jenner RG, Boyer LA, Guenther MG, Levine SS, Kumar RM, Chevalier B, Johnstone SE, Cole MF, Isono K, et al. Control of developmental regulators by Polycomb in human embryonic stem cells. Cell. 2006;125:301–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Ohm JE, McGarvey KM, Yu X, Cheng L, Schuebel KE, Cope L, Mohammad HP, Chen W, Daniel VC, Yu W, et al. A stem cell-like chromatin pattern may predispose tumor suppressor genes to DNA hypermethylation and heritable silencing. Nat Genet. 2007;39:237–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Schlesinger Y, Straussman R, Keshet I, Farkash S, Hecht M, Zimmerman J, Eden E, Yakhini Z, Ben-Shushan E, Reubinoff BE, et al. Polycomb-mediated methylation on Lys27 of histone H3 pre-marks genes for de novo methylation in cancer. Nat Genet. 2007;39:232–6.View ArticlePubMedGoogle Scholar
- Yadav BS, Chanana P, Jhamb S. Biomarkers in triple negative breast cancer: a review. World J Clin Oncol. 2015;6:252–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Eccles SA, Aboagye EO, Ali S, Anderson AS, Armes J, Berditchevski F, Blaydes JP, Brennan K, Brown NJ, Bryant HE, et al. Critical research gaps and translational priorities for the successful prevention and treatment of breast cancer. Breast Cancer Res. 2013;15:R92.View ArticlePubMedPubMed CentralGoogle Scholar
- Fackler MJ, Umbricht CB, Williams D, Argani P, Cruz LA, Merino VF, Teo WW, Zhang Z, Huang P, Visvananthan K, et al. Genome-wide methylation analysis identifies genes specific to breast cancer hormone receptor status and risk of recurrence. Cancer Res. 2011;71:6195–207.View ArticlePubMedPubMed CentralGoogle Scholar
- Lehmann BD, Bauer JA, Chen X, Sanders ME, Chakravarthy AB, Shyr Y, Pietenpol JA. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J Clin Invest. 2011;121:2750–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Tao Y, Xi S, Briones V, Muegge K. Lsh mediated RNA polymerase II stalling at HoxC6 and HoxC8 involves DNA methylation. PLoS One. 2010;5:e9163.View ArticlePubMedPubMed CentralGoogle Scholar
- Keller PJ, Lin AF, Arendt LM, Klebba I, Jones AD, Rudnick JA, DiMeo TA, Gilmore H, Jefferson DM, Graham RA, et al. Mapping the cellular and molecular heterogeneity of normal and malignant breast tissues and cultured cell lines. Breast Cancer Res. 2010;12:R87.View ArticlePubMedPubMed CentralGoogle Scholar
- Teo WW, Merino VF, Cho S, Korangath P, Liang X, Wu RC, Neumann NM, Ewald AJ, Sukumar S. HOXA5 determines cell fate transition and impedes tumor initiation and progression in breast cancer through regulation of E-cadherin and CD24. Oncogene. 2016;35(42):5539–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Ordonez-Moran P, Dafflon C, Imajo M, Nishida E, Huelsken J. HOXA5 counteracts stem cell traits by inhibiting Wnt signaling in colorectal cancer. Cancer Cell. 2015;28:815–29.View ArticlePubMedGoogle Scholar
- Shashikant CS, Ruddle FH. Combinations of closely situated cis-acting elements determine tissue-specific patterns and anterior extent of early Hoxc8 expression. Proc Natl Acad Sci U S A. 1996;93:12364–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Le Mouellic H, Lallemand Y, Brulet P. Homeosis in the mouse induced by a null mutation in the Hox-3.1 gene. Cell. 1992;69:251–64.View ArticlePubMedGoogle Scholar
- Lewis MT. Homeobox genes in mammary gland development and neoplasia. Breast Cancer Res. 2000;2:158–69.View ArticlePubMedPubMed CentralGoogle Scholar
- Kalyani R, Lee JY, Min H, Yoon H, Kim MH. Genes frequently coexpressed with Hoxc8 provide insight into the discovery of target genes. Mol Cells. 2016;39:395–402.View ArticlePubMedPubMed CentralGoogle Scholar
- Carroll LS, Capecchi MR. Hoxc8 initiates an ectopic mammary program by regulating Fgf10 and Tbx3 expression and Wnt/beta-catenin signaling. Development. 2015;142:4056–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Li Y, Chao F, Huang B, Liu D, Kim J, Huang S. HOXC8 promotes breast tumorigenesis by transcriptionally facilitating cadherin-11 expression. Oncotarget. 2014;5:2596–607.View ArticlePubMedPubMed CentralGoogle Scholar
- Waltregny D, Alami Y, Clausse N, de Leval J, Castronovo V. Overexpression of the homeobox gene HOXC8 in human prostate cancer correlates with loss of tumor differentiation. Prostate. 2002;50:162–9.View ArticlePubMedGoogle Scholar
- Alami Y, Castronovo V, Belotti D, Flagiello D, Clausse N. HOXC5 and HOXC8 expression are selectively turned on in human cervical cancer cells compared to normal keratinocytes. Biochem Biophys Res Commun. 1999;257:738–45.View ArticlePubMedGoogle Scholar
- Du YB, Dong B, Shen LY, Yan WP, Dai L, Xiong HC, Liang Z, Kang XZ, Qin B, Chen KN. The survival predictive significance of HOXC6 and HOXC8 in esophageal squamous cell carcinoma. J Surg Res. 2014;188:442–50.View ArticlePubMedGoogle Scholar
- Adwan H, Zhivkova-Galunska M, Georges R, Eyol E, Kleeff J, Giese NA, Friess H, Bergmann F, Berger MR. Expression of HOXC8 is inversely related to the progression and metastasis of pancreatic ductal adenocarcinoma. Br J Cancer. 2011;105:288–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Lei H, Juan AH, Kim MS, Ruddle FH. Identification of a Hoxc8-regulated transcriptional network in mouse embryo fibroblast cells. Proc Natl Acad Sci U S A. 2006;103:10305–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Morgan R, Boxall A, Harrington KJ, Simpson GR, Gillett C, Michael A, Pandha HS. Targeting the HOX/PBX dimer in breast cancer. Breast Cancer Res Treat. 2012;136:389–98.View ArticlePubMedGoogle Scholar
- Ginestier C, Wicinski J, Cervera N, Monville F, Finetti P, Bertucci F, Wicha MS, Birnbaum D, Charafe-Jauffret E. Retinoid signaling regulates breast cancer stem cell differentiation. Cell Cycle. 2009;8:3297–302.View ArticlePubMedPubMed CentralGoogle Scholar
- Sun R, Liu Y, Li SY, Shen S, Du XJ, Xu CF, Cao ZT, Bao Y, Zhu YH, Li YP, et al. Co-delivery of all-trans-retinoic acid and doxorubicin for cancer therapy with synergistic inhibition of cancer stem cells. Biomaterials. 2015;37:405–14.View ArticlePubMedGoogle Scholar