- Research article
- Open Access
- Open Peer Review
Genomic pathways modulated by Twist in breast cancer
BMC Cancervolume 17, Article number: 52 (2017)
The basic helix-loop-helix transcription factor TWIST1 (Twist) is involved in embryonic cell lineage determination and mesodermal differentiation. There is evidence to indicate that Twist expression plays a role in breast tumor formation and metastasis, but the role of Twist in dysregulating pathways that drive the metastatic cascade is unclear. Moreover, many of the genes and pathways dysregulated by Twist in cell lines and mouse models have not been validated against data obtained from larger, independant datasets of breast cancer patients.
We over-expressed the human Twist gene in non-metastatic MCF-7 breast cancer cells to generate the estrogen-independent metastatic breast cancer cell line MCF-7/Twist. These cells were inoculated in the mammary fat pad of female severe compromised immunodeficient mice, which subsequently formed xenograft tumors that metastasized to the lungs. Microarray data was collected from both in vitro (MCF-7 and MCF-7/Twist cell lines) and in vivo (primary tumors and lung metastases) models of Twist expression. Our data was compared to several gene datasets of various subtypes, classes, and grades of human breast cancers.
Our data establishes a Twist over-expressing mouse model of breast cancer, which metastasizes to the lung and replicates some of the ontogeny of human breast cancer progression. Gene profiling data, following Twist expression, exhibited novel metastasis driver genes as well as cellular maintenance genes that were synonymous with the metastatic process. We demonstrated that the genes and pathways altered in the transgenic cell line and metastatic animal models parallel many of the dysregulated gene pathways observed in human breast cancers.
Analogous gene expression patterns were observed in both in vitro and in vivo Twist preclinical models of breast cancer metastasis and breast cancer patient datasets supporting the functional role of Twist in promoting breast cancer metastasis. The data suggests that genetic dysregulation of Twist at the cellular level drives alterations in gene pathways in the Twist metastatic mouse model which are comparable to changes seen in human breast cancers. Lastly, we have identified novel genes and pathways that could be further investigated as targets for drugs to treat metastatic breast cancer.
Despite significant advancements in understanding the biology of breast cancer progression, there are still many aspects of this disease that are not well understood, such as the metastatic process, which is the largest cause of mortality in breast cancer patients. There has been a concerted effort to understand the mechanics of this process in order to provide new opportunities to design targeted interventions to inhibit or prevent metastatic progression. A number of in vitro and in vivo mouse preclinical models have been used to identify key components that are drivers of the metastatic process. However, given the heterogeneity of breast cancer and the relevance of the model used to the human disease, there are significant quantitative and qualitative differences between human breast tumors and preclinical mouse models. One possible alternative is to define gene-specific alterations and higher order pathways that are dysregulated both in human breast cancers and in preclinical mouse models and cell lines. If the dysregulated gene signatures identified from breast cancer patient samples can be functionally validated in cell lines and in xenograft models, this may help in understanding breast cancer progression. Simultaneously, if genetic changes that contribute to breast cancer progression can be catalogued using cell line and xenograft-based models, it may be possible to derive gene signatures that, at least partly, replicate the progression of breast cancer in humans. Importantly, this could be accomplished in cell lines that ordinarily do not express these genes. It may also be possible to altogether avoid mouse models in certain types of breast cancers if cell lines are demonstrated to be validated surrogates of the cancer.
In our quest to genetically define breast cancer progression, we identified the gene Twist (TWIST1), which promotes the metastatic process in humans. Twist is a transcription factor belonging to the basic helix-loop-helix (bHLH) family of proteins and is essential for normal vertebrate development [1–3]. Twist is overexpressed in breast [4, 5], gastric , and prostate cancers , as well as in melanomas , gliomas , osteosarcomas , and rhabdomyosarcomas . Functionally, Twist overexpression promotes breast cancer by increasing angiogenesis and chromosomal instability [4, 12], and by downregulating E-cadherin  and estrogen receptor . We have also demonstrated that Twist promotes the breast cancer stem cell phenotype by regulating CD24 . Even though the functional roles of Twist in cancer biogenesis are well documented, there is, however, little evidence presented in the literature characterizing the global dysregulation of genes and pathways brought about as a direct consequence of Twist overexpression [4, 5]. Here, we show by establishing a new mouse orthotopic xenograft model, that breast tumors overexpressing Twist are highly metastatic to the lungs. We also demonstrate that the genes and pathways altered in our Twist over expressing cell line and metastatic xenograft mouse model recapitulate some of the observations seen in human breast tumors. Our data establishes a new mouse model of metastatic breast cancer that can be used to test the efficacy of therapies for breast cancer treatment.
Cell lines and animal experiments
MCF-7 cell line was originally obtained from ATCC (Manassas, VA). MCF-7/Twist cell line was created earlier  and constitutively over-expresses Twist. One million MCF-7/Twist cells were injected into the second left mammary fat pad of female SCID mice which were 4–6 weeks old. Breast tumors were allowed to grow for a period of 2 months after which the tumors were resected in a sterile environment. MCF-7/Twist primary tumor explants were then grown in tissue culture till cell lines were obtained. Fibroblast contamination was removed by using low concentrations of G418 for short periods. After a further period of 2 months, mice were sacrificed and lungs were examined for metastases. Lungs were minced under sterile conditions and made into a single cell suspension by dispersing the tissue using collagenase and DNase I treatment. Single cells were then plated in cell culture dishes and expanded till enough MCF-7/Twist metastatic cells were obtained for RNA extraction. RNA was extracted from all cell lines in the early part of their growth for performing microarray analysis.
The RNA samples were converted into double stranded cDNA which was subsequently transcribed into biotinylated complementary RNA. After purification and fragmentation, the cRNA was hybridized to GeneChip Human Genome U133 Plus 2.0 array chips (Affymetrix, Santa Clara, CA), and scanned using an Affymetrix GeneChip Scanner 3000 with default parameters to obtain gene expression profiles. Affymetrix CEL file data was extracted and their data was normalized using RMA (Robust Multi-array Average), and averaged over biological replicates, creating the quantile-normalized log2 transcript signal values used in subsequent ANOVA analyses. Analysis was performed using Genomic Suite Software 6.6 (Partek, St. Louis, Missouri).
Genes whose expression differed by at least 1.5 fold from the median in at least 20% of the arrays were retained. We identified genes that were differentially expressed among the classes using a multivariate permutation test . We used the multivariate permutation test to provide 90% confidence that the false discovery rate was less than 10%. The false discovery rate is the proportion of the list of genes claimed to be differentially expressed that are false positives. The multivariate permutation test is non-parametric and does not require the assumption of Gaussian distributions.
Hierarchical clustering analysis
To confirm the relationship between the various classes, we performed unsupervised hierarchical clustering of classes using the UPGMA (Unweighted Pair Group Method with Arithmetic Mean) method.
To further study the gene expression differences among the different samples, we visualized the data by volcano plots. The log2 fold-change in expression is seen on the x-axis with red depicting mRNAs that are up-regulated and green depicting mRNAs that are down-regulated. The y-axis depicts the -log10 significance between the classes. The horizontal line depicts P = 0.05 and values above were considered statistically significant.
Experimental study design and bioinformatics analysis
To study genomic pathways that are disrupted by Twist expression in breast cancer, we designed an experimental plan that encompasses cell culture, a mouse xenograft model, and human breast cancer samples (Fig. 1). We stably overexpressed Twist in MCF-7 breast cancer cells to establish the MCF-7/Twist cell line. Functional analysis revealed that this cell line is aggressive, forms large primary tumors, and is metastatic when implanted orthotopically in the breast of SCID mice . Inoculation of MCF-7/Twist cells into the mammary fat pad of SCID mice gave rise to tumors in all mice (n = 8) within 1–2 months. These tumors were surgically resected two months after inoculation and expanded into cell lines. After an additional period of 1–2 months, all mice had developed lungs metastases, which were confirmed by histological analysis of the lungs. We thus established an orthotopic mouse xenograft model of breast cancer that mirrors the ontogeny of tumor formation and progression of lung metastasis in humans (Fig. 1). In order to analyze the genes and pathways that were dysregulated by Twist in our models, we studied RNA transcripts from the parental MCF-7 and the Twist over expressing MCF-7/Twist cell lines, primary orthotopic tumors and lung metastases in mice. RNA was analyzed by microarray analysis.
A total of 54,675 probe sets were obtained from microarray analysis after normalization. Of these, using the NetAffyx Annotation Release 35 dataset, probe sets of unknown origin (8,841), controls (62), and probes mapping to multiple gene assignments (2,796) were discarded. The remaining 42,976 probe sets were mapped to 21,655 genes.
To study the overall variation between the sample classes, we initially performed unsupervised PCA analysis (Fig. 2a). PCA analysis showed that the parental MCF-7 clustered with the MCF-7/Twist cell line, while the xenograft tumors clustered with the metastatic tumors, and the samples from breast cancers formed an independent tightly knit cluster. The cell lines MCF-7 and MCF-7/Twist clustered relatively further apart from one another indicating that over expression of Twist caused large genomic changes in the parental MCF-7 cells. From the close clustering of the orthotopic primary and metastatic tumors, we surmise that the gene alterations induced by Twist expression were attenuated to some degree in the in vivo mouse model.
Hierarchical clustering analysis was performed on the MCF-7 and MCF-7/Twist cell lines, on the mouse primary xenograft tumors, and lung metastases and on the human breast tumors (Fig. 2b). From the hierarchical clustering analysis, we determined the closest relationship observed was within human tumor samples and the mouse xenograft samples (Fig. 2b). We also observed the clustering of cell lines according to their propensity to metastasize with MCF-7 cell being farthest from MCF-7/Twist lung metastases.
Gene expression analysis
We compared mRNA expression between all three groups (MCF-7/Twist cells, primary tumors, and lung metastases) using parental MCF-7 cells as controls (Fig. 3a). There were a total of 1,451 probe sets that were common between all three datasets (MCF-7/Twist cells, MCF-7/Twist primary tumors and metastases) of which 516 were up-regulated and 935 were down-regulated (Fig. 3a). See Additional File 1 for the complete list of genes up and down-regulated in all three classes. We validated the top two genes by quantitative reverse transcription polymerase chain reaction (qRT-PCR) in these two datasets and the results are shown in Additional File 2.
Next, we analyzed the expression data for genes that were differentially expressed in the three groups (Fig. 3b). When comparing multidirectional gene expression in three classes, 12 probe sets were up-regulated in MCF-7/Twist while being down-regulated in the primary tumors and metastases. On the other hand, we found 7 probe sets down-regulated in MCF-7/Twist that were up-regulated in primary tumors and metastases. Subsequently, we compared multidirectional gene expression across two classes. 228 probe sets were up-regulated in the MCF-7/Twist cell line while being down-regulated in primary tumors. Similarly, 156 probe sets were up-regulated in Twist cells while down-regulated in metastatic tumors. Conversely, 98 probe sets were down-regulated in MCF-7/Twist cells while they were up regulated in primary tumors. Similarly, 87 probe sets were down-regulated in Twist cells while being up-regulated in metastatic tumors. Overall, we found that Twist caused larger changes in gene expression in the MCF-7/Twist cell line as compared to changes caused in the primary or metastatic tumor classes. Complete probe sets are listed in Additional File 3.
A total of 9,278 probe sets (6,478 genes) had their expression significantly altered (FC > 1.5, P < 0.05) in the MCF-7/Twist cell line (5,185 probe sets/3,497 genes up-regulated and 4,093 probe sets/2,981 genes down-regulated). The probe sets differentially expressed (up and down-regulated) in MCF-7/Twist cells are listed in Additional file 4. The top two genes in both directions were validated by qRT-PCR and the result is shown in Additional file 5.
In the primary tumors, 12,355 probe sets (10,182 genes) had altered expression compared to MCF-7 (5,887 probe sets/4,771 genes up-regulated and 6,648 probe sets/5,411 genes down-regulated) (See Additional File 6 for complete list). The qRT-PCR validation is displayed in Additional File 7.
In comparison to MCF-7 cells, the metastatic tumors had 15,085 probe sets (12,183 genes) that were significantly altered of which 7,225 probe sets/5,825 genes were up-regulated and 7,860 probe sets/6,358 genes were down-regulated (Additional File 8). The top genes from the up and down-regulated datasets were validated by qRT-PCR (Additional File 9).
As could be expected, bioinformatics analysis indicated that the largest overlap of 11,077 probe sets occurred between the primary and metastatic tumors (5,348 up- and 5,729 down-regulated). This is displayed in Fig. 3a, and the complete list is in Additional File 10. Additional File 11 displays the results of the validation of this data by qRT-PCR.
There were 1,849 probe sets in common between the MCF-7/Twist cell line and metastasis (635 up- and 1,214 down-regulated), the complete list of which is displayed in Additional File 12. Validation of the top 2 genes (up and down-regulated) was performed by qRT-PCR and is shown in Additional File 13.
Figure 3a also exhibits that the least overlapping gene signature was between the MCF-7/Twist cells and primary tumor classes, where 1,764 probe sets were regulated commonly (645 up- and 1,119 down-regulated). All the expressed probe sets are displayed in Additional File 14 and the validation performed by qRT-PCR is shown in Additional File 15.
Intriguingly, the MCF-7/Twist cell line had more down-regulated probe sets in common with the metastatic tumors and more up-regulated probe sets in common with the primary tumors.
Highly expressed transcripts in MCF-7/Twist cells include EDIL3, which has, so far, not been found to be involved in cancer development but is involved in immune system-dependent engulfment of apoptotic cells . Another high expressing transcript includes LDHB, which is involved in triple negative breast cancer . ESR1 (estrogen receptor alpha), GATA3 (GATA Binding Protein 3), and CDH1 (E-cadherin) are amongst the genes most down-regulated by Twist expression and shown to be involved in the progression of breast cancer earlier [13, 14, 19]. Genes that were up-regulated in primary and metastatic tumors included SETD5 (SET domain containing 5) and TFAP2A (transcription factor AP-2 alpha). Intriguingly, very little is known about SETD5 in breast cancer progression and metastasis , while TFAP2A is involved in both sporadic and hereditary breast cancer [21, 22]. SLC6A14 (solute carrier family 6, member 14) and TCN1 (transcobalamin 1) were down-regulated genes in primary and metastatic tumors. On comparing probe sets from all 3 classes, we found that LMO3 (LIM-only protein 3) was the highest significantly up-regulated gene, while SLC6A14 was down-regulated the most amongst the genes. Further details and complete gene lists can be found in additional files.
To determine the variation in gene expression compared to controls, we displayed our microarray datasets with the help of volcano plots (Fig. 4a). Genes that were significantly different between the two comparison classes are displayed above the horizontal line (P < 0.05). As seen from the volcano plots, the maximum changes induced by Twist overexpression was in the MCF-7/Twist cell line, followed by the primary tumors, and lastly the metastatic tumors.
Finally, we compared the genetic changes caused by Twist expression in the breast cell line and in the orthotopic mouse xenograft and metastasis model. As seen in the scatter plots (Fig. 4b), the highest correlation is seen between the primary and metastatic tumors in mouse while much lower correlations were seen when comparing the MCF-7/Twist cell line data with the primary and metastatic tumor data.
Genes dysregulated by Twist are observed in human breast cancer datasets
To study how our Twist overexpressing cell line and orthotopic mouse primary and metastatic mouse model compared with human breast cancers, we analyzed data from the van’t Veer , Sotiriou , and Perou datasets . From the van’t Veer “poor prognosis” dataset consisting of 50 genes, a total of 38, 37, and 40 genes were present in the MCF-7/Twist, primary, and metastatic datasets, respectively (Fig. 5a). Of these, 26 genes overlapped in all 3 classes. The van’t Veer “breast cancer up in metastasis gene set” was present in 36, 36, and 37 of 51 genes in MCF-7/Twist, primary, and metastatic datasets, respectively. In the other “breast cancer down in metastasis gene set” was present in 89, 87, and 93 of 114 genes in the respective datasets (Fig. 5b). Overall, the van’t Veer breast cancer metastasis dataset had 90 genes (22 up-regulated and 68 down-regulated) in common across all 3 of our datasets. From the van’t Veer “estrogen receptor positive breast cancer gene set” of 384 genes, we found an overlap of 281, 272, and 288 genes in our MCF-7/Twist, primary and metastatic datasets, respectively (Fig. 5c). This high degree of overlap corroborates our earlier work, which demonstrated that Twist expression correlates with and represses ER expression [4, 14]. Moreover, it is a validation of the role of Twist in breast cancer progression.
Additionally, we compared the Sotiriou dataset that has a gene signature that discriminates grade 3 from grade 1 breast tumors . A total of 99 genes from our datasets overlapped with the Sotiriou dataset (195 genes) (Fig. 5d). The MCF-7/Twist, primary tumor, and metastatic tumor datasets had 147, 151, and 142 genes, respectively, overlapping with the Sotiriou dataset.
Lastly, we compared our datasets with the Perou gene sets that classify breast tumors into intrinsic sub-types . As expected, we observed the least amount of overlapping genes when compared to the Luminal A and B datasets (Fig. 6a, b). Higher concordance (6 genes) between datasets was observed with the HER2 enriched dataset (Fig. 6c). As seen in Fig. 6d, the largest overlap was observed in the basal gene set with 8 genes being common between the 3 classes.
Rho GTPase and IL-6 signaling pathways are upregulated by Twist
Subsequently, when analyzed by IPA, the most significant pathways up-regulated in all three classes were signaling by Rho GTPases, IL-6, NGF, JAK/Stat, and CCK/Gastrin (Fig. 7a, b). On the other hand, PPAR, PTEN, PPARα/RXRα, dopamine receptor, and Wnt/β-catenin signaling were reduced. While comparing upstream regulators dysregulated in the three classes, we found RELA, CD40LG, STAT3, FOXM1, and CD24 the most upregulated. miR-122, CXCL12, COL18A1, miR-1, and TP73 were the top five upstream regulator molecules that were down-regulated in our dataset (Fig. 7c, d).
Metastases are the single largest cause of deaths from all cancers. The currently accepted model states that metastases arise when over a period of time a sub population of cells from the primary tumor acquire genetic alterations that confer upon them selective advantages that allow them to metastasize . This widely belief has been challenged earlier [27, 28] and has recently come under increasing scrutiny due to microarray analysis of tumors and their associated metastasis [23, 29]. These studies demonstrate that the so-called “gene signature” of metastatic tumors is not very different from primary tumors. This implies that the genetic alterations that lead to metastasis are acquired much earlier than is presently anticipated - perhaps during formation of the primary tumor itself. The lack of progress in treatment of metastases is partly caused by the scarcity of clinically relevant animal models of breast cancer and metastases [30, 31].
We have earlier demonstrated that the Twist gene is involved in the progression of breast cancers towards a highly aggressive phenotype [4, 12, 14]. In this work, we over expressed Twist in non-metastatic MCF-7 and orthotopically implanted these cells in mouse mammary fat pads. After resecting the primary tumors, we monitored the animals for the development of metastatic lesions, primarily in lungs. Unlike other studies that directly assayed gene changes from primary and metastatic sites, we grew the cells out in plastic for 1–2 weeks before microarray analysis. This allows for us to lessen the impact of mouse cells that could have infiltrated the orthotopic tumors. Moreover, this also eliminates any fibroblast cells that could be carried over from the metastatic sites in lungs. To lessen batch variation, we used the same gene arrays (Affymetrix Human Genome HGU133 Plus 2.0) and the same core facility for all the work in this study. We also analyzed human breast cancers in our study. However, due to the limited sample size, our study was not powered enough to draw conclusions from this part of the study. Instead, we compared our Twist model with currently accepted human gene signature databases ([23–25, 32].
As expected, our model showed a high species-dependent correlation. The primary and metastatic tumors clustered strongly both in PCA analysis as well as by unsupervised hierarchical clustering. Strong clustering was also seen in the breast cancer samples. The PCA distribution of MCF-7 and MCF-7/Twist showed a higher than expected degree of separation. This could be due to the very strong effect of Twist on the MCF-7 cell line, which becomes ER-negative after Twist overexpression . Over time, as tumors establish in mice, the effects of Twist may be ameliorated due to selection pressure from the various other dysregulated pathways of tumor growth.
In our model, the largest number of genes were dysregulated in the metastatic tumor class followed by the primary tumors, and least in the MCF-7/Twist class. This was expected since one could deduce that metastatic breast tumors are genetically the furthest away from the parental MCF-7 cell line, which formed the baseline of all the comparisons. We also observed that a high number of genes were commonly dysregulated in the primary and metastatic tumors. However, we also observed that MCF-7/Twist cells had a higher number of genes in common with the metastatic tumor class as compared to the primary tumors. This would indicate that the MCF-7/Twist cell line in vitro was a better representative of the metastatic process in mice than the primary tumors. When comparing the classes for genes that were expressed in opposite directions, as expected, we only found a significant number of genes when MCF-7/Twist was compared with the primary and metastatic tumor classes (either together or singly). This was further corroborated by quantile plots, which clearly indicated the close relationship between the primary breast tumors and the metastasis.
When comparing our datasets with the van’t Veer poor prognosis gene set, we observed that a large proportion of genes from all 3 classes overlapped with that dataset leading us to hypothesize that our model does recapitulate, in some measure, the progression of breast cancer in humans. Similarly, the Sotiriou gene set also exhibited a high degree of concordance with our dataset implying that the Twist model is a suitable model system for high-grade breast tumors. When comparing the Perou dataset of intrinsic sub-types of breast cancer, we observed that our datasets overlapped the most with the basal and HER-2 enriched datasets. A very low degree of overlap was seen with the Luminal A and Luminal B sub-types.
To model which biological functions and pathways were dysregulated by Twist overexpression, we employed IPA analysis. Overall, we saw an increase in transcription and translational activity, and involving genes related to cellular biogenesis and apoptosis. Also, enriched were genes in the protease family and extracellular matrix. This was expected since Twist is a master transcription factor for many pathways and is also involved in apoptosis and in modulating the extracellular matrix to cause invasion and metastasis.
In summary, we generated Twist-driven models of breast cancer progression, both in vitro and in vivo, that mimic the ontology of breast cancer formation and metastatic development in humans. These datasets will help in identifying Twist-associated dysregulated genes and pathways that can be targeted to prevent or decrease breast cancer metastases. Importantly, our work provides evidence that gene-specific functions can be ascertained by using preclinical models that can be associated with the ontogeny of human breast cancer formation and progression.
Robust multi-array Average, qRT-PCR (quantitative reverse transcription polymerase chain reaction)
Howard TD, Paznekas WA, Green ED, Chiang LC, Ma N. Ortiz de Luna RI, Garcia Delgado C, Gonzalez-Ramos M, Kline AD, jabs EW: mutations in TWIST, a basic helix-loop-helix transcription factor, in saethre-chotzen syndrome. Nat Genet. 1997;15(1):36–41.
Kunz J, Hudler M, Fritz B. Identification of a frameshift mutation in the gene TWIST in a family affected with robinow-sorauf syndrome. J Med Genet. 1999;36(8):650–2.
Gripp KW, Stolle CA, Celle L, McDonald-McGinn DM, Whitaker LA, Zackai EH. TWIST gene mutation in a patient with radial aplasia and craniosynostosis: further evidence for heterogeneity of Baller-Gerold syndrome. Am J Med Genet. 1999;82(2):170–6.
Mironchik Y, Winnard Jr PT, Vesuna F, Kato Y, Wildes F, Pathak AP, Kominsky S, Artemov D, Bhujwalla ZM, Van Diest P, et al. Twist overexpression induces in vivo angiogenesis and correlates with chromosomal instability in breast cancer. Cancer Res. 2005;65(23):10801–9.
Yang J, Mani SA, Donaher JL, Ramaswamy S, Itzykson RA, Come C, Savagner P, Gitelman I, Richardson A, Weinberg RA. Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis. Cell. 2004;117(7):927–39.
Rosivatz E, Becker I, Specht K, Fricke E, Luber B, Busch R, Hofler H, Becker KF. Differential expression of the epithelial-mesenchymal transition regulators snail, SIP1, and twist in gastric cancer. Am J Pathol. 2002;161(5):1881–91.
Kwok WK, Ling MT, Lee TW, Lau TC, Zhou C, Zhang X, Chua CW, Chan KW, Chan FL, Glackin C, et al. Up-regulation of TWIST in prostate cancer and its implication as a therapeutic target. Cancer Res. 2005;65(12):5153–62.
Hoek K, Rimm DL, Williams KR, Zhao H, Ariyan S, Lin A, Kluger HM, Berger AJ, Cheng E, Trombetta ES, et al. Expression profiling reveals novel pathways in the transformation of melanocytes to melanomas. Cancer Res. 2004;64(15):5270–82.
Elias MC, Tozer KR, Silber JR, Mikheeva S, Deng M, Morrison RS, Manning TC, Silbergeld DL, Glackin CA, Reh TA, et al. TWIST is expressed in human gliomas and promotes invasion. Neoplasia. 2005;7(9):824–37.
Entz-Werle N, Stoetzel C, Berard-Marec P, Kalifa C, Brugiere L, Pacquement H, Schmitt C, Tabone MD, Gentet JC, Quillet R, et al. Frequent genomic abnormalities at TWIST in human pediatric osteosarcomas. Int J Cancer. 2005;117(3):349–55.
Maestro R, Dei Tos AP, Hamamori Y, Krasnokutsky S, Sartorelli V, Kedes L, Doglioni C, Beach DH, Hannon GJ. Twist is a potential oncogene that inhibits apoptosis. Genes Dev. 1999;13(17):2207–17.
Vesuna F, Winnard Jr P, Glackin C, Raman V. Twist overexpression promotes chromosomal instability in the breast cancer cell line MCF-7. Cancer Genet Cytogenet. 2006;167(2):189–91.
Vesuna F, van Diest P, Chen JH, Raman V. Twist is a transcriptional repressor of E-cadherin gene expression in breast cancer. Biochem Biophys Res Commun. 2008;367(2):235–41.
Vesuna F, Lisok A, Kimble B, Domek J, Kato Y, van der Groep P, Artemov D, Kowalski J, Carraway H, van Diest P, et al. Twist contributes to hormone resistance in breast cancer by downregulating estrogen receptor-alpha. Oncogene. 2012;31(27):3223–34.
Vesuna F, Lisok A, Kimble B, Raman V. Twist modulates breast cancer stem cells by transcriptional regulation of CD24 expression. Neoplasia. 2009;11(12):1318–28.
Korn EL, Troendle JF, McShane LM, Simon R. Controlling the number of false discoveries: application to high-dimensional genomic data. J Stat Planning and Inference. 2004;124(2):379–98.
Hanayama R, Tanaka M, Miwa K, Nagata S. Expression of developmental endothelial locus-1 in a subset of macrophages for engulfment of apoptotic cells. J Immunol. 2004;172(6):3876–82.
McCleland ML, Adler AS, Shang Y, Hunsaker T, Truong T, Peterson D, Torres E, Li L, Haley B, Stephan JP, et al. An integrated genomic screen identifies LDHB as an essential gene for triple-negative breast cancer. Cancer Res. 2012;72(22):5812–23.
McCleskey BC, Penedo TL, Zhang K, Hameed O, Siegal GP, Wei S. GATA3 expression in advanced breast cancer: prognostic value and organ-specific relapse. Am J Clin Pathol. 2015;144(5):756–63.
Sowalsky AG, Xia Z, Wang L, Zhao H, Chen S, Bubley GJ, Balk SP, Li W. Whole transcriptome sequencing reveals extensive unspliced mRNA in metastatic castration-resistant prostate cancer. Mol Cancer Res: MCR. 2015;13(1):98–106.
Kang HJ, Lee MH, Kang HL, Kim SH, Ahn JR, Na H, Na TY, Kim YN, Seong JK, Lee MO. Differential regulation of estrogen receptor alpha expression in breast cancer cells by metastasis-associated protein 1. Cancer Res. 2014;74(5):1484–94.
Gaudet MM, Kuchenbaecker KB, Vijai J, Klein RJ, Kirchhoff T, McGuffog L, Barrowdale D, Dunning AM, Lee A, Dennis J, et al. Identification of a BRCA2-specific modifier locus at 6p24 related to breast cancer risk. PLoS Genet. 2013;9(3):e1003173.
Van’t Veer LJ, Dai H, Van De Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, Van Der Kooy K, Marton MJ, Witteveen AT, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002;415(6871):530–6.
Sotiriou C, Wirapati P, Loi S, Harris A, Fox S, Smeds J, Nordgren H, Farmer P, Praz V, Haibe-Kains B, et al. Gene expression profiling in breast cancer: understanding the molecular basis of Histologic grade to improve prognosis. J Natl Cancer Inst. 2006;98(4):262–72.
Parker JS, Mullins M, Cheang MC, Leung S, Voduc D, Vickery T, Davies S, Fauron C, He X, Hu Z, et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 2009;27(8):1160–7.
Fidler IJ, Kripke ML. Metastasis results from preexisting variant cells within a malignant tumor. Science. 1977;197(4306):893–5.
Vaage J. Metastasizing potentials of mouse mammary tumors and their metastases. Int J Cancer. 1988;41(6):855–8.
Milas L, Peters LJ, Ito H. Spontaneous metastasis: random or selective? Clin Exp Metastasis. 1983;1(4):309–15.
Ramaswamy S, Ross KN, Lander ES, Golub TR. A molecular signature of metastasis in primary solid tumors. Nat Genet. 2003;33(1):49–54.
Yi B, Williams PJ, Niewolna M, Wang Y, Yoneda T. Tumor-derived platelet-derived growth factor-BB plays a critical role in osteosclerotic bone metastasis in an animal model of human breast cancer. Cancer Res. 2002;62(3):917–23.
Kurebayashi J, McLeskey SW, Johnson MD, Lippman ME, Dickson RB, Kern FG. Quantitative demonstration of spontaneous metastasis by MCF-7 human breast cancer cells cotransfected with fibroblast growth factor 4 and LacZ. Cancer Res. 1993;53(9):2178–87.
Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, et al. Molecular portraits of human breast tumours. Nature. 2000;406(6797):747–52.
We would also like to thank Connie Talbot Jr. (JHMI Microarray Core) for his support in microarray analysis.
This study was supported by the National Institutes of Health by grants R01CA131250 and R01CA140226 to VR.
Availability of data and materials
The datasets supporting the conclusions of this article are available in the Gene Expression Omnibus (GEO) repository [GSE87705, http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE87705.
FV devised the study concept, acquired the data, interpreted and visualized the data, and drafted and revised the manuscript. YB performed the experiments. VR devised the study, interpreted the data, and drafted and revised the manuscript. All authors approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Animals were studied under protocols approved by the Johns Hopkins University Animal Care and Use Committee.
Ranked list of significant probe sets up and down-regulated in MCF-7/Twist cells, primary tumors, and metastases. Probe sets are ranked by average fold change across all 3 groups. (XLSX 148 kb)
Histogram displaying results of validation of genes up-regulated (LMO3 and BICC1) and down-regulated (DSCAM-AS1 and SLC6A14) in the MCF-7/Twist cell line, primary tumors, and metastases. (EPS 1432 kb)
Ranked list of probe sets expressed multi-directionally in all three groups. (XLSX 22 kb)
Ranked list of significant probe sets most differentially expressed in MCF-7/Twist cells. Probe sets are ranked by fold change. (XLSX 461 kb)
Histogram displaying results of qRT-PCR validation of genes up-regulated (EDIL3 and LDHB), and down-regulated (DSCAM-AS1 and NPY1R) in the MCF-7/Twist cell line. (EPS 1258 kb)
Ranked list of significant probe sets most differentially expressed in MCF-7/Twist primary tumors. Probe sets are ranked by fold change. (XLSX 601 kb)
Histogram of qRT-PCR validation of genes up-regulated (SETD5 and TFAP2A) and down-regulated (SLC6A14 and LINC01419) in MCF-7/Twist primary tumors. (EPS 1253 kb)
Ranked list of significant probe sets differentially expressed in MCF-7/Twist lung metastases. Probe sets are ranked by fold change. (XLSX 734 kb)
qRT-PCR validation results of genes up-regulated (SETD5 and TFAP2A) and down-regulated (SLC6A14 and HPGD) in MCF-7/Twist lung metastases. (EPS 1253 kb)
Ranked list of significant probe sets up and down-regulated in MCF-7/Twist primary tumors and MCF-7/Twist lung metastases. Probe sets are ranked by fold change. (XLSX 740 kb)
qRT-PCR analysis of genes up-regulated (SETD5 and TFAP2A), and down-regulated (NR2F1 and GJA1) in MCF-7/Twist primary xenograft tumors and in lung metastases. (EPS 1264 kb)
Ranked list of significant probe sets up and down-regulated in both MCF-7/Twist cells and MCF-7/Twist lung metastases. Probe sets are ranked by fold change. (XLSX 38 kb)
Histogram showing qRT-PCR validation of genes up-regulated (EIF1AY and PLAU) and down-regulated (WDR72 and AGR2) in MCF-7/Twist cell line and in lung metastases. (EPS 1266 kb)
Ranked list of significant probe sets up and down-regulated in both MCF-7/Twist cells and MCF-7/Twist primary tumors. Probe sets are ranked by fold change. (XLSX 32 kb)
Histogram displaying qRT-PCR validation results of up-regulated genes (POPDC3 and CXADR) and down-regulated genes (NPY1R and IGFBP5) in MCF-7/Twist cells and primary tumors. (EPS 1268 kb)