Discovery, optimization and validation of an optimal DNA-binding sequence for the Six1 homeodomain transcription factor

Transcriptional Regulation of De Novo Lipogenesis by SIX1 in Liver Cancer Cells

De novo lipogenesis (DNL), a hallmark of cancer, facilitates tumor growth and metastasis. Therapeutic drugs targeting DNL are being developed. However, how DNL is directly regulated in cancer remains largely unknown. Here, transcription factor sine oculis homeobox 1 (SIX1) is shown to directly increase the expression of DNL-related genes, including ATP citrate lyase (ACLY), fatty acid synthase (FASN), and stearoyl-CoA desaturase 1 (SCD1), via histone acetyltransferases amplified in breast cancer 1 (AIB1) and lysine acetyltransferase 7 （HBO1/KAT7）, thus promoting lipogenesis. SIX1 expression is regulated by insulin/lncRNA DGUOK-AS1/microRNA-145-5p axis, which also modulates DNL-related gene expression as well as DNL. The DGUOK-AS1/microRNA-145-5p/SIX1 axis regulates liver cancer cell proliferation, invasion, and metastasis in vitro and in vivo. In patients with liver cancer, SIX1 expression is positively correlated with DGUOK-AS1 and SCD1 expression and is negatively correlated with microRNA-145-5p expression. DGUOK-AS1 is a good predictor of prognosis. Thus, the DGUOK-AS1/microRNA-145-5p/SIX1 axis strongly links DNL to tumor growth and metastasis and may become an avenue for liver cancer therapeutic intervention.

Altered binding affinity of SIX1-Q177R correlates with enhanced WNT5A and WNT pathway effector expression in Wilms tumor

Wilms tumors present as an amalgam of varying proportions of tissues located within the developing kidney, one being the nephrogenic blastema comprising multipotent nephron progenitor cells (NPCs). The recurring missense mutation Q177R in NPC transcription factors SIX1 and SIX2 is most correlated with tumors of blastemal histology and is significantly associated with relapse. Yet, the transcriptional regulatory consequences of SIX1/2-Q177R that might promote tumor progression and recurrence have not been investigated extensively. Utilizing multiple Wilms tumor transcriptomic datasets, we identified upregulation of the gene encoding non-canonical WNT ligand WNT5A in addition to other WNT pathway effectors in SIX1/2-Q177R mutant tumors. SIX1 ChIP-seq datasets from Wilms tumors revealed shared binding sites for SIX1/SIX1-Q177R within a promoter of WNT5A and at putative distal cis-regulatory elements (CREs). We demonstrate colocalization of SIX1 and WNT5A in Wilms tumor tissue and utilize in vitro assays that support SIX1 and SIX1-Q177R activation of expression from the WNT5A CREs, as well as enhanced binding affinity within the WNT5A promoter that may promote the differential expression of WNT5A and other WNT pathway effectors associated with SIX1-Q177R tumors.

Ginsenoside Rh4 Suppresses Metastasis of Gastric Cancer via SIX1-Dependent TGF-β/Smad2/3 Signaling Pathway

Gastric cancer (GC) is the leading causes of cancer-related death worldwide. Surgery remains the cornerstone of gastric cancer treatment, and new strategies with adjuvant chemotherapy are currently gaining more and more acceptance. Ginsenoside Rh4 has excellent antitumor activity. Conversely, the mechanisms involved in treatment of GC are not completely understood. In this study, we certified that Rh4 showed strong anti-GC efficiency in vitro and in vivo. MTT and colony formation assays were performed to exhibit that Rh4 significantly inhibited cellular proliferation and colony formation. Results from the wound healing assay, transwell assays, and Western blotting indicated that Rh4 restrained GC cell migration and invasion by reversing epithelial–mesenchymal transition (EMT). Further validation by proteomic screening, co-treatment with disitertide, and SIX1 signal silencing revealed that SIX1, a target of Rh4, induced EMT by activating the TGF-β/Smad2/3 signaling pathway. In summary, our discoveries demonstrated the essential basis of the anti-GC metastatic effects of Rh4 via suppressing the SIX1–TGF-β/Smad2/3 signaling axis, which delivers a new idea for the clinical treatment of GC.

Six1 promotes skeletal muscle thyroid hormone response through regulation of the MCT10 transporter

BACKGROUND

The Six1 transcription factor is implicated in controlling the development of several tissue types, notably skeletal muscle. Six1 also contributes to muscle metabolism and its activity is associated with the fast-twitch, glycolytic phenotype. Six1 regulates the expression of certain genes of the fast muscle program by directly stimulating their transcription or indirectly acting through a long non-coding RNA. We hypothesized that additional mechanisms of action of Six1 might be at play.

METHODS

A combined analysis of gene expression profiling and genome-wide location analysis data was performed. Results were validated using in vivo RNA interference loss-of-function assays followed by measurement of gene expression by RT-PCR and transcriptional reporter assays.

RESULTS

The Slc16a10 gene, encoding the thyroid hormone transmembrane transporter MCT10, was identified as a gene with a transcriptional enhancer directly bound by Six1 and requiring Six1 activity for full expression in adult mouse tibialis anterior, a predominantly fast-twitch muscle. Of the various thyroid hormone transporters, MCT10 mRNA was found to be the most abundant in skeletal muscle, and to have a stronger expression in fast-twitch compared to slow-twitch muscle groups. Loss-of-function of MCT10 in the tibialis anterior recapitulated the effect of Six1 on the expression of fast-twitch muscle genes and led to lower activity of a thyroid hormone receptor-dependent reporter gene.

CONCLUSIONS

These results shed light on the molecular mechanisms controlling the tissue expression profile of MCT10 and identify modulation of the thyroid hormone signaling pathway as an additional mechanism by which Six1 influences skeletal muscle metabolism.

Secondary structures, dynamics, and DNA binding of the homeodomain of human SIX1

Human sine oculis homeobox homolog (SIX) 1 contains a homeodomain (HD), which is important for binding to DNA. In this study, we carried out structural studies on the HD of human SIX1 using nuclear magnetic resonance (NMR) spectroscopy. Its secondary structures and dynamics in solution were explored. HD is well-structured in solution, and our study shows that it contains three α-helices. Dynamics study indicates that the N- and C-terminal residues of HD are flexible in solution. HD of human SIX1 exhibits molecular interactions with a short double-strand DNA sequence evidenced by the ¹ H-¹⁵ N-heteronuclear single quantum correlation (HSQC) and ¹⁹ F-NMR experiments. Our current study provides structural information for HD of human SIX1. Further studies indicate that this construct can be utilized to study SIX1 and DNA interactions.

Silencing SIX1 inhibits epithelial mesenchymal transition through regulating TGF-β/Smad2/3 signaling pathway in papillary thyroid carcinoma

OBJECTIVE

To investigate the sineoculis homeobox homolog 1 (SIX1) affect the epithelial mesenchymal transition (EMT) in papillary thyroid carcinoma (PTC) through regulating TGF-β/Smad2/3 signaling pathway.

METHODS

The SIX1 expression in cytological specimens, tissues or PTC cell lines was detected by qRT-PCR, western blotting or immunohistochemistry. A series of vitro experiments including flow cytometry, CCK-8, wound-healing and Transwell were used to evaluate the biological characteristics in a PTC cell line (NPA cells), which were divided into Blank, Negative control (NC), SIX1, SIX1-siRNA, LY-364947 (TGF-β/Smad2/3 pathway inhibitor) and SIX1 + LY-364947 groups. TGF-β/Smad2/3 pathway and EMT related protein expression were measured by qRT-PCR and western blotting.

RESULTS

SIX1 mRNA expression was increased in cytological specimens from PTC patients as compared with the non-toxic nodular goitre (NTG) patients. Moreover, compared with adjacent normal tissues, expressions of SIX1, N-cadherin and Vimentin were higher while E-cadherin was lower in PTC tissues; and SIX1 was positively correlated with N-cadherin and Vimentin but was negatively correlated with E-cadherin. Furthermore, the SIX1 expression was associated with histopathology, extrathyroidal extension (ETE), lymph node metastasis (LNM), pT stage, TNM stage, and distant metastasis. In addition, the expressions of TGFβ1, p-SMAD2/3, N-cadherin and Vimentin were downregulated in NPA cells after LY-364947 treatment with upregulated E-cadherin, decreased cell proliferation and metastasis, and enhanced cell apoptosis, which was reversed by SIX1 overexpression.

CONCLUSION

Silencing SIX1 can inhibit TGF-β/Smad2/3 pathway, thereby suppressing EMT in PTC, which may be a novel avenue for the treatment of PTC.

Myogenesis control by SIX transcriptional complexes

SIX homeoproteins were first described in Drosophila, where they participate in the Pax-Six-Eya-Dach (PSED) network with eyeless, eyes absent and dachsund to drive synergistically eye development through genetic and biochemical interactions. The role of the PSED network and SIX proteins in muscle formation in vertebrates was subsequently identified. Evolutionary conserved interactions with EYA and DACH proteins underlie the activity of SIX transcriptional complexes (STC) both during embryogenesis and in adult myofibers. Six genes are expressed throughout muscle development, in embryonic and adult proliferating myogenic stem cells and in fetal and adult post-mitotic myofibers, where SIX proteins regulate the expression of various categories of genes. In vivo, SIX proteins control many steps of muscle development, acting through feedforward mechanisms: in the embryo for myogenic fate acquisition through the direct control of Myogenic Regulatory Factors; in adult myofibers for their contraction/relaxation and fatigability properties through the control of genes involved in metabolism, sarcomeric organization and calcium homeostasis. Furthermore, during development and in the adult, SIX homeoproteins participate in the genesis and the maintenance of myofibers diversity.

SIX1 represses senescence and promotes SOX2-mediated cellular plasticity during tumorigenesis

Six1 is a developmental transcriptional regulator frequently overexpressed in human tumors. Recent results show that SIX1 also acts as a repressor of cell senescence, an antiproliferative response with a key role in tumor suppression, among other physiological and pathological settings. Here, we set to study the impact of SIX1 gain of function in transformation and tumorigenesis of fibroblasts, in connection with senescence. Using transcriptomic, histological, and functional analyses in murine tumors and cells of fibroblast origin, we show that SIX1 has a strong pro-tumorigenic action in this model, linked to the repression of a senescence-related gene signature and the induction of an undifferentiated phenotype mediated, at least in part, by the regulation of the stemness factor Sox2. Moreover, functional analyses with human glioma cell lines also show that SIX1 controls SOX2 expression, senescence and self-renewal in this model. Collectively, our results support a general link of SIX1 with senescence and SOX2-mediated cell plasticity in tumors.

PTE, a novel module to target Polycomb Repressive Complex 1 to the human cyclin D2 (CCND2) oncogene

Polycomb group proteins are essential epigenetic repressors. They form multiple protein complexes of which two kinds, PRC1 and PRC2, are indispensable for repression. Although much is known about their biochemical properties, how mammalian PRC1 and PRC2 are targeted to specific genes is poorly understood. Here, we establish the cyclin D2 (CCND2) oncogene as a simple model to address this question. We provide the evidence that the targeting of PRC1 to CCND2 involves a dedicated PRC1-targeting element (PTE). The PTE appears to act in concert with an adjacent cytosine-phosphate-guanine (CpG) island to arrange for the robust binding of PRC1 and PRC2 to repressed CCND2 Our findings pave the way to identify sequence-specific DNA-binding proteins implicated in the targeting of mammalian PRC1 complexes and provide novel link between polycomb repression and cancer.

Transcriptional Regulation of the Warburg Effect in Cancer by SIX1

Aerobic glycolysis (the Warburg effect) facilitates tumor growth, and drugs targeting aerobic glycolysis are being developed. However, how the Warburg effect is directly regulated is largely unknown. Here we show that transcription factor SIX1 directly increases the expression of many glycolytic genes, promoting the Warburg effect and tumor growth in vitro and in vivo. SIX1 regulates glycolysis through HBO1 and AIB1 histone acetyltransferases. Cancer-related SIX1 mutation increases its ability to promote aerobic glycolysis and tumor growth. SIX1 glycolytic function is directly repressed by microRNA-548a-3p, which is downregulated, inversely correlates with SIX1, and is a good predictor of prognosis in breast cancer patients. Thus, the microRNA-548a-3p/SIX1 axis strongly links aerobic glycolysis to carcinogenesis and may become a promising cancer therapeutic target.

Differential regulation of mouse and human nephron progenitors by the Six family of transcriptional regulators

Nephron endowment is determined by the self-renewal and induction of a nephron progenitor pool established at the onset of kidney development. In the mouse, the related transcriptional regulators Six1 and Six2 play non-overlapping roles in nephron progenitors. Transient Six1 activity prefigures, and is essential for, active nephrogenesis. By contrast, Six2 maintains later progenitor self-renewal from the onset of nephrogenesis. We compared the regulatory actions of Six2 in mouse and human nephron progenitors by chromatin immunoprecipitation followed by DNA sequencing (ChIP-seq). Surprisingly, SIX1 was identified as a SIX2 target unique to the human nephron progenitors. Furthermore, RNA-seq and immunostaining revealed overlapping SIX1 and SIX2 activity in 16 week human fetal nephron progenitors. Comparative bioinformatic analysis of human SIX1 and SIX2 ChIP-seq showed each factor targeted a similar set of cis-regulatory modules binding an identical target recognition motif. In contrast to the mouse where Six2 binds its own enhancers but does not interact with DNA around Six1, both human SIX1 and SIX2 bind homologous SIX2 enhancers and putative enhancers positioned around SIX1. Transgenic analysis of a putative human SIX1 enhancer in the mouse revealed a transient, mouse-like, pre-nephrogenic, Six1 regulatory pattern. Together, these data demonstrate a divergence in SIX-factor regulation between mouse and human nephron progenitors. In the human, an auto/cross-regulatory loop drives continued SIX1 and SIX2 expression during active nephrogenesis. By contrast, the mouse establishes only an auto-regulatory Six2 loop. These data suggest differential SIX-factor regulation might have contributed to species differences in nephron progenitor programs such as the duration of nephrogenesis and the final nephron count.

The Campylobacter jejuni Ferric Uptake Regulator Promotes Acid Survival and Cross-Protection against Oxidative Stress

Campylobacter jejuni is a prevalent cause of bacterial gastroenteritis in humans worldwide. The mechanisms by which C. jejuni survives stomach acidity remain undefined. In the present study, we demonstrated that the C. jejuni ferric uptake regulator (Fur) plays an important role in C. jejuni acid survival and acid-induced cross-protection against oxidative stress. A C. jejuni Δfur mutant was more sensitive to acid than the wild-type strain. Profiling of the acid stimulon of the C. jejuni Δfur mutant allowed us to uncover Fur-regulated genes under acidic conditions. In particular, Fur was found to upregulate genes involved in flagellar and cell envelope biogenesis upon acid stress, and mutants with deletions of these genes were found to be defective in surviving acid stress. Interestingly, prior acid exposure of C. jejuni cross-protected against oxidative stress in a catalase (KatA)- and Fur-dependent manner. Western blotting and reverse transcription-quantitative PCR revealed increased expression of KatA upon acid stress. Electrophoretic mobility shift assays (EMSAs) demonstrated that the binding affinity between Fur and the katA promoter is reduced in vitro under conditions of low pH, rationalizing the higher levels of expression of katA under acidic conditions. Strikingly, the Δfur mutant exhibited reduced virulence in both human epithelial cells and the Galleria mellonella infection model. Altogether, this is the first study showing that, in addition to its role in iron metabolism, Fur is an important regulator of C. jejuni acid responses and this function cross-protects against oxidative stress. Moreover, our results clearly demonstrate Fur's important role in C. jejuni pathogenesis.

Six family of homeobox genes and related mechanisms in tumorigenesis protocols

In recent years, the homeobox gene superfamily has been introduced as a master regulator in downstream target genes related to cell development and proliferation. An indispensable role of this family involved in organogenesis development has been widely demonstrated since expression of Six family led to a distinct increase in development of various organs. These functions of Six family genes are primarily based on structure as well as regulatory role in response to external or internal stimuli. In addition to these roles, mutation or aberrant expression of Six family plays a fundamental role in initiation of carcinogenesis, a multistep process including transformation, proliferation, angiogenesis, migration, and metastasis. This suggests that the Six superfamily members can be considered as novel target molecules to inhibit tumor growth and progression. This review focuses on the structure, function, and mechanisms of the Six family in cancer processes and possible strategies to apply these family members for diagnostic, prognostic, and therapeutic purposes.

Homeodomain proteins: an update

Here, we provide an update of our review on homeobox genes that we wrote together with Walter Gehring in 1994. Since then, comprehensive surveys of homeobox genes have become possible due to genome sequencing projects. Using the 103 Drosophila homeobox genes as example, we present an updated classification. In animals, there are 16 major classes, ANTP, PRD, PRD-LIKE, POU, HNF, CUT (with four subclasses: ONECUT, CUX, SATB, and CMP), LIM, ZF, CERS, PROS, SIX/SO, plus the TALE superclass with the classes IRO, MKX, TGIF, PBC, and MEIS. In plants, there are 11 major classes, i.e., HD-ZIP (with four subclasses: I to IV), WOX, NDX, PHD, PLINC, LD, DDT, SAWADEE, PINTOX, and the two TALE classes KNOX and BEL. Most of these classes encode additional domains apart from the homeodomain. Numerous insights have been obtained in the last two decades into how homeodomain proteins bind to DNA and increase their specificity by interacting with other proteins to regulate cell- and tissue-specific gene expression. Not only protein-DNA base pair contacts are important for proper target selection; recent experiments also reveal that the shape of the DNA plays a role in specificity. Using selected examples, we highlight different mechanisms of homeodomain protein-DNA interaction. The PRD class of homeobox genes was of special interest to Walter Gehring in the last two decades. The PRD class comprises six families in Bilateria, and tinkers with four different motifs, i.e., the PAIRED domain, the Groucho-interacting motif EH1 (aka Octapeptide or TN), the homeodomain, and the OAR motif. Homologs of the co-repressor protein Groucho are also present in plants (TOPLESS), where they have been shown to interact with small amphipathic motives (EAR), and in yeast (TUP1), where we find an EH1-like motif in MATα2.

Genome-wide association between Six4, MyoD, and the histone demethylase Utx during myogenesis

Adult skeletal muscles can regenerate after injury, due to the presence of satellite cells, a quiescent population of myogenic progenitor cells. Once activated, satellite cells repair the muscle damage by undergoing myogenic differentiation. The myogenic regulatory factors (MRFs) coordinate the process of progenitor differentiation in cooperation with other families of transcription factors (TFs). The Six1 and Six4 homeodomain TFs are expressed in developing and adult muscle and Six1 is critical for embryonic and adult myogenesis. However, the lack of a muscle developmental phenotype in Six4-null mice, which has been attributed to compensation by other Six family members, has discouraged further assessment of the role of Six4 during adult muscle regeneration. By employing genome-wide approaches to address the function of Six4 during adult skeletal myogenesis, we have identified a core set of muscle genes coordinately regulated in adult muscle precursors by Six4 and the MRF MyoD. Throughout the genome of differentiating adult myoblasts, the cooperation between Six4 and MyoD is associated with chromatin repressive mark removal by Utx, a demethylase of histone H3 trimethylated at lysine 27. Among the genes coordinately regulated by Six4 and MyoD are several genes critical for proper in vivo muscle regeneration, implicating a role of Six4 in this process. Using in vivo RNA interference of Six4, we expose an uncompensated function of this TF during muscle regeneration. Together, our results reveal a role for Six4 during adult muscle regeneration and suggest a widespread mechanism of cooperation between Six4 and MyoD.

The regulation of gene expression in hair cells

No genes have been discovered for which expression is limited only to inner ear hair cells. This is hardly surprising, since the number of mammalian genes is estimated to be 20-25,000, and each gene typically performs many tasks in various locations. Many genes are expressed in inner ear sensory cells and not in other cells of the labyrinth. However, these genes are also expressed in other locations, often in other sensory or neuronal cell types. How gene transcription is directed specifically to hair cells is unclear. Key transcription factors that act during development can specify cell phenotypes, and the hair cell is no exception. The transcription factor ATOH1 is well known for its ability to transform nonsensory cells of the developing inner ear into hair cells. And yet, ATOH1 also specifies different sensory cells at other locations, neuronal phenotypes in the brain, and epithelial cells in the gut. How it specifies hair cells in the inner ear, but alternate cell types in other locations, is not known. Studies of regulatory DNA and transcription factors are revealing mechanisms that direct gene expression to hair cells, and that determine the hair cell identity. The purpose of this review is to summarize what is known about such gene regulation in this key auditory and vestibular cell type.

Down-regulation of miR-29b is essential for pathogenesis of uterine leiomyoma

Uterine leiomyomata (LMs) are the most common tumor affecting the female reproductive organs. The most notable pathophysiologic feature of this tumor is the excessive accumulation of rigid extracellular matrix (ECM) composed mainly of collagen types I and III. It is believed that the rigidity of the collagen-rich ECM causes symptoms such as abnormal bleeding and pelvic pain/pressure. However, the molecular pathogenesis for this ECM-rich tumor has yet to be elucidated. We have established that miR-29b was consistently down-regulated in LM compared with myometrium (MM). Hence, the function of miR-29b in LM was examined in vivo using adult female ovariectomized NOD-scid IL2Rγ(null) mice for subrenal xenograft models. In LM xenografts, restoring miR-29b inhibited the accumulation of ECM and the development of solid tumors. Although the miR-29b knockdown in MM cells increased the expression of collagens, it did not transform MM cells into tumorigenic, indicating that the down-regulation of miR-29b is essential but not sufficient for LM tumorigenesis. In addition, 17β-estradiol and progesterone down-regulated miR-29b and up-regulated mRNAs for multiple collagens in LM xenografts. Thus, we conclude that ECM production in LMs is regulated by steroid hormones via down-regulation of miR-29b, which is one of the mechanisms underlying the excessive accumulation of ECM.

Six1 regulates MyoD expression in adult muscle progenitor cells

Quiescent satellite cells are myogenic progenitors that enable regeneration of skeletal muscle. One of the early events of satellite cell activation following myotrauma is the induction of the myogenic regulatory factor MyoD, which eventually induces terminal differentiation and muscle function gene expression. The purpose of this study was to elucidate the mechanism by which MyoD is induced during activation of satellite cells in mouse muscle undergoing regeneration. We show that Six1, a transcription factor essential for embryonic myogenesis, also regulates MyoD expression in muscle progenitor cells. Six1 knock-down by RNA interference leads to decreased expression of MyoD in myoblasts. Chromatin immunoprecipitation assays reveal that Six1 binds the Core Enhancer Region of MyoD. Further, transcriptional reporter assays demonstrate that Core Enhancer Region reporter gene activity in myoblasts and in regenerating muscle depends on the expression of Six1 and on Six1 binding sites. Finally, we provide evidence indicating that Six1 is required for the proper chromatin structure at the Core Enhancer Region, as well as for MyoD binding at its own enhancer. Together, our results reveal that MyoD expression in satellite cells depends on Six1, supporting the idea that Six1 plays an important role in adult myogenesis, in addition to its role in embryonic muscle formation.

Six1 promotes proliferation of pancreatic cancer cells via upregulation of cyclin D1 expression

Six1 is one of the transcription factors that act as master regulators of development and are frequently dysregulated in cancers. However, the role of Six1 in pancreatic cancer is not clear. Here we show that the relative expression of Six1 mRNA is increased in pancreatic cancer and correlated with advanced tumor stage. In vitro functional assays demonstrate that forced overexpression of Six1 significantly enhances the growth rate and proliferation ability of pancreatic cancer cells. Knockdown of endogenous Six1 decreases the proliferation of these cells dramatically. Furthermore, Six1 promotes the growth of pancreatic cancer cells in a xenograft assay. We also show that the gene encoding cyclin D1 is a direct transcriptional target of Six1 in pancreatic cancer cells. Overexpression of Six1 upregulates cyclin D1 mRNA and protein, and significantly enhances the activity of the cyclin D1 promoter in PANC-1 cells. We demonstrate that Six1 promotes cell cycle progression and proliferation by upregulation of cyclin D1. These data suggest that Six1 is overexpressed in pancreatic cancer and may contribute to the increased cell proliferation through upregulation of cyclin D1.

Structure-function analyses of the human SIX1-EYA2 complex reveal insights into metastasis and BOR syndrome

SIX1 interacts with EYA to form a bipartite transcription factor essential for development. Loss of function of this complex causes branchio-oto-renal syndrome (BOR), while re-expression of SIX1 or EYA promotes metastasis. Here we describe the 2.0 Å structure of SIX1 bound to EYA2, which suggests a novel DNA binding mechanism for SIX1 and provides a rationale for the effect of BOR syndrome mutations. The structure also reveals that SIX1 uses predominantly a single helix to interact with EYA. Substitution of a single amino acid in this helix is sufficient to disrupt the SIX1–EYA interaction, SIX1-mediated epithelial-mesenchymal transition and metastasis in mouse models. Given that SIX1 and EYA are co-overexpressed in many tumor types, our data indicate that targeting the SIX1–EYA complex may be a potent approach to inhibit tumor progression in multiple cancer types.