Dynamic compartmentalization of the pro-invasive transcription factor NHR-67 reveals a role for Groucho in regulating a proliferative-invasive cellular switch in C. elegans

A growing body of evidence suggests that cell division and basement membrane invasion are mutually exclusive cellular behaviors. How cells switch between proliferative and invasive states is not well understood. Here, we investigated this dichotomy in vivo by examining two cell types in the developing Caenorhabditis elegans somatic gonad that derive from equipotent progenitors, but exhibit distinct cell behaviors: the post-mitotic, invasive anchor cell and the neighboring proliferative, non-invasive ventral uterine (VU) cells. We show that the fates of these cells post-specification are more plastic than previously appreciated and that levels of NHR-67 are important for discriminating between invasive and proliferative behavior. Transcription of NHR-67 is downregulated following post-translational degradation of its direct upstream regulator, HLH-2 (E/Daughterless) in VU cells. In the nuclei of VU cells, residual NHR-67 protein is compartmentalized into discrete punctae that are dynamic over the cell cycle and exhibit liquid-like properties. By screening for proteins that colocalize with NHR-67 punctae, we identified new regulators of uterine cell fate maintenance: homologs of the transcriptional co-repressor Groucho (UNC-37 and LSY-22), as well as the TCF/LEF homolog POP-1. We propose a model in which the association of NHR-67 with the Groucho/TCF complex suppresses the default invasive state in non-invasive cells, which complements transcriptional regulation to add robustness to the proliferative-invasive cellular switch in vivo.


Introduction
Cellular proliferation and invasion are key aspects of development (reviewed in Medwig and Matus, 2017), and are also two of the defining hallmarks of cancer (reviewed in Hanahan and Weinberg, 2000).A growing body of evidence suggests that cell cycle progression and invasion through a basement membrane are mutually exclusive cellular behaviors in both development and disease states (reviewed in Kohrman and Matus, 2017).Switching between invasive and proliferative phenotypes has been observed in melanoma and recently in breast cancer (Hoek et al., 2008;Mondal et al., 2022), but how these cell states are regulated in the context of development is not well understood.To investigate how this dichotomy in cellular behavior is controlled in vivo, we used C. elegans, leveraging its highly stereotypical development (Sulston and Horvitz, 1977), as well as its genetic and optical tractability.During the development of the hermaphroditic reproductive system, the proximal granddaughters of the Z1 and Z4 somatic gonad progenitors, Z1. pp and Z4.aa, give rise to four cells that will adopt one of two cellular fates: a proliferative VU cell or the terminally differentiated, invasive anchor cell (AC) (Figure 1A; Kimble and Hirsh, 1979).The distal cells of this competency group, Z1. ppa and Z4.aap, quickly lose their bipotentiality and become VU cells (Seydoux et al., 1990).In contrast, the proximal cells, Z1. ppp and Z4.aaa, undergo a stochastic Notch-mediated cell fate decision, giving rise to another VU cell and the post-mitotic AC (Figure 1A and B;Greenwald et al., 1983;Seydoux and Greenwald, 1989).Following fate specification, the AC undergoes invasive differentiation and breaches the underlying basement membrane, connecting the uterus to the vulval epithelium to facilitate egg-laying (Figure 1B; Sherwood and Sternberg, 2003).The AC is the default cell fate of the somatic gonad, as disruption of Notch or Wnt signaling results in ectopic AC specification (Phillips et al., 2007;Seydoux and Greenwald, 1989).Therefore, AC fate must be actively repressed.
Our previous work has shown that AC invasion is dependent on G 0 cell cycle arrest, which is coordinated by the pro-invasive transcription factor NHR-67 (NR2E1/TLX) (Figure 1-figure supplement 1A; Matus et al., 2015).NHR-67 functions within a gene regulatory network comprised of four conserved transcription factors whose homologs have been implicated in several types of metastatic cancer (Liang and Wang, 2020;Milde-Langosch, 2005;Nelson et al., 2021;Wang and Baker, 2015).We previously reported that NHR-67 is regulated by a feed-forward loop formed by  and HLH-2 (E/Daughterless), which functions largely in parallel to a cell cycle-independent subcircuit controlled by FOS-1 (Fos) (Figure 1-figure supplement 1A; Medwig-Kinney et al., 2020).EGL-43, HLH-2, and NHR-67 are reiteratively used within the Z lineage of the somatic gonad, in that, they also function to independently regulate  signaling during the initial AC/VU cell fate decision (Medwig-Kinney et al., 2020).Despite its role in lateral inhibition between Z1. ppp and Z4.aaa, expression of LIN-12 is not absolutely required for VU cell fate (Sallee et al., 2015a).Cell cycle state also cannot explain the difference between AC and VU cell fates, as arresting VU cells in G 0 through ectopic expression of CKI-1 (p21/p27) does not make them invasive (Smith et al., 2022).Thus, the mechanisms responsible for maintaining AC and VU cellular identities following initial cell fate specification remain unclear.
Maintenance of differentiated cell identity is essential for ensuring tissue integrity during development and homeostasis, and the inability to restrict phenotypic plasticity is now being recognized as an integral part of cancer pathogenesis (Hanahan, 2022).In vitro studies have identified several factors that safeguard differentiated cell identity (reviewed in Brumbaugh et al., 2019).Despite its largely autonomous modality of development, C. elegans has emerged as an ideal model system to study cell fate maintenance in vivo.There have been several reports of cell fate transformations that occur naturally, including two epithelial-to-neural transdifferentiation events (Jarriault et al., 2008;Riva et al., 2022), or following fate challenges (reviewed in Rothman and Jarriault, 2019).In such contexts, several epigenetic factors, including chromatin remodelers and histone chaperones, have been identified for their roles in restricting cell fate reprogramming (Hajduskova et al., 2019;Kagias et al., 2012;Kolundzic et al., 2018;Patel et al., 2012;Rahe and Hobert, 2019;Zuryn et al., 2014).However, in some cases, ectopic expression of a specific transcription factor is sufficient to overcome these barriers, as was first shown through pioneering work in mouse embryonic fibroblasts (Davis et al., 1987).Indeed, there are several examples in C. elegans where ectopic expression of single lineage-specific transcription factors induces cell fate transformations (Fukushige and Krause, 2005;Gilleard and McGhee, 2001;Horner et al., 1998;Jin et al., 1994;Kiefer et al., 2007;Quintin et al., 2001;Richard et al., 2011;Riddle et al., 2013;Tursun et al., 2011;Zhu et al., 1998).Moreover, C. elegans uterine tissue may be particularly amenable to fate transformations, as ectopic expression of a single GATA transcription factor, ELT-7, is sufficient to induce transorganogenesis of the somatic gonad into the gut by reprogramming the mesodermally-derived tissue into endoderm (Riddle et al., 2016).Valuable insights have been made into how the function of fate-specifying transcription factors can be tuned through means such as autoregulation and dynamic heterodimerization (Leyva-Díaz and Hobert, 2019;Sallee et al., 2017).We are just beginning to understand how an additional layer of control over transcriptional regulators can be achieved through compartmentalization (Boija et al., 2018;Lim and Levine, 2021).
Here, in our endeavor to understand how AC and VU cellular fates are maintained, we identified two mechanisms that together modulate the invasive-proliferative switch in C. elegans.We found that high levels of NHR-67 expression are sufficient to drive invasive differentiation, and that NHR-67 is transcriptionally downregulated in the non-invasive VU cells following the post-translational degradation of its direct upstream regulator, HLH-2.Additionally, we observed that the remaining NHR-67 protein in the VU cells compartmentalizes into punctae that exhibit liquid-like properties including dynamic assembly, fusion, and dissolution over the cell cycle as well as rapid recovery kinetics after photobleaching.These NHR-67 punctae colocalize in vivo with UNC-37 and LSY-22, homologs of the transcriptional co-repressor Groucho, as well as with POP-1 (TCF/LEF), which is likely mediated through a direct interaction between UNC-37 and the intrinsically disordered C-terminal region of NHR-67.Through functional perturbations, we demonstrate that UNC-37, LSY-22, and POP-1 contribute to the repression of the default invasive state in VU cells.We propose a model in which NHR-67 compartmentalizes through its interaction with Groucho, which, combined with transcriptional downregulation of NHR-67, may suppress invasive differentiation.

NHR-67 expression levels are important for distinguishing AC and VU cell identity
Despite arising from initially equipotent cells, the differentiated AC and VU cells exhibit very distinct cellular behaviors.The AC terminally differentiates to invade the underlying basement membrane while the VU cells remain proliferative, undergoing several rounds of division before terminally differentiating.One potential explanation for this difference in cell behavior is asymmetric expression of pro-invasive transcription factors.To investigate this possibility, we examined endogenous expression levels of four transcription factors that function in the gene regulatory network coordinating    AC invasion (EGL-43, FOS-1, HLH-2, and NHR-67) using GFP-tagged alleles (Medwig-Kinney et al., 2020).While FOS-1 levels are enriched in the AC compared to the VU cells (Figure 1-figure supplement 1B and C), FOS-1 has no known role in cell cycle regulation, so we did not pursue this protein further (Medwig-Kinney et al., 2021).EGL-43 was also not a promising candidate, as it is expressed in both cell types at comparable levels, with VU cells exhibiting approximately 89% of AC expression (Figure 1-figure supplement 1B and C).In contrast, HLH-2 exhibits significant asymmetry in expression, as VU cells express merely 17% of HLH-2 levels observed in the AC (Figure 1C and D).Previous studies have shown that dimerization-driven degradation of HLH-2 is responsible for its downregulation in the VU cells (Benavidez et al., 2022;Karp and Greenwald, 2003;Sallee and Greenwald, 2015b).NHR-67 exhibits a similar pattern of expression with over threefold enrichment in the AC, consistent with prior observations of transgenic reporters (Figure 1C and D;Verghese et al., 2011).Given the known role of NHR-67 in regulating cell cycle and invasion, we hypothesized that its differential expression between the AC and VU cells could contribute to their distinct cellular behaviors.
To assess the potential role of NHR-67 in regulating uterine cell identities, we manipulated its expression levels.We found that strong depletion of NHR-67 through RNA interference (RNAi) treatment results in ACs adopting VU-like characteristics.During AC/VU cell fate specification, LIN-12 (Notch) normally becomes restricted to the VU cells while the Delta-like ligand LAG-2 (visualized by LAG-2::P2A::H2B::mTurquoise2 Medwig- Kinney et al., 2022) accumulates in the AC (Wilkinson et al., 1994).Here, we observe that NHR-67 deficient ACs not only proliferated and failed to invade, as reported previously (Matus et al., 2015), but also ectopically expressed membrane-localized Notch (visualized by LIN-12::mNeonGreen Pani et al., 2022;Figure 1E;Figure 1-figure supplement 2A and B).Notably, NHR-67-deficient ACs expressed both LIN-12 and LAG-2, potentially indicating an intermediate state between AC and VU cell fate (Figure 1E; Figure 1-figure supplement 2A and  B).Next, we ectopically expressed NHR-67 ubiquitously following AC/VU specification using a heat shock inducible transgene (hsp::NHR-67::2x-BFP) (Medwig-Kinney et al., 2020) and observed the presence of multiple invasive ACs at a low penetrance (approximately 5%, n>50), denoted by ectopic expression of an AC marker (cdh-3p::mCherry::moeABD) and expansion of the basement membrane gap (Figure 1F).As previous work has demonstrated that proliferative ACs cannot invade (Matus et al., 2015), we concluded that these invasive ectopic ACs most likely arose from the fate conversion of neighboring VU cells.Taken together, these pieces of evidence suggest that high and low levels of NHR-67 correlate to properties of AC and VU cell identities, respectively (Figure 1G).

NHR-67 is enriched in the AC through transcriptional regulation by HLH-2
Next, we investigated how NHR-67 expression levels become asymmetric between the AC and VU cells.We and others have previously shown that HLH-2 positively regulates NHR-67 expression in the context of the AC (Figure 1-figure supplement 1A; Bodofsky et al., 2018;Medwig-Kinney et al., 2020).If this regulatory interaction also exists in the context of the VU cells, it could explain why the relative expression pattern of NHR-67 in the AC and VU cells mirrors that of HLH-2.In support of this hypothesis, we found that the initial onset of HLH-2, which has been shown to be asymmetric in Z1. pp and Z4.aa (Attner et al., 2019), correlates to onset of an NHR-67 transgene (Figure 2-figure supplement 1A; Gerstein et al., 2010).To test whether HLH-2 degradation is responsible for NHR-67 downregulation in the VU, we drove ectopic expression of HLH-2 using a transgene under the control of a heat shock inducible promoter (hsp::HLH-2::2x-BFP) (Medwig-Kinney et al., 2020) and observed elevated NHR-67 expression in VU cells (43% increase; n>30) (Figure 2A and B).To control against potential dimerization-driven degradation of HLH-2 in the VU cells, which the heat shock inducible transgene would still be susceptible to, we disrupted UBA-1, an E1 ubiquitin-activating enzyme that has recently been shown to be necessary for HLH-2 degradation in VU cells (Benavidez et al., 2022).Following perturbation of UBA-1 through RNAi treatment, HLH-2 expression in the VU cells increased more than fourfold and NHR-67 expression increased by nearly 60% compared to the empty vector control (Figure 2-figure supplement 1B-D).Both experiments suggest that NHR-67 expression in the VU cells is at least partially regulated by levels of HLH-2.
It has previously been proposed that the interaction between HLH-2 and NHR-67 is direct.This is based on the identification of E binding motifs within a 276 bp region of the NHR-67 promoter that is required for NHR-67 expression in the uterine tissue and encompasses the location of several hypomorphic mutations (pf2, pf88, pf159) (Figure 2C; Bodofsky et al., 2018;Verghese et al., 2011).We performed a yeast one-hybrid assay by generating a bait strain containing this NHR-67 promoter region and pairing it with an HLH-2 Gal4-AD prey plasmid from an existing yeast one-hybrid library (Reece-Hoyes et al., 2005).Yeast growth on the selective SC-HIS-TRP plates containing the competitive inhibitor 3-aminotriazole (3-AT) suggests that HLH-2 is capable of binding directly to this 276 bp region of the NHR-67 promoter (Figure 2D).Together, these results support that direct transcriptional regulation of NHR-67 by HLH-2 contributes to the asymmetry in NHR-67 expression between the AC and VU cells.

NHR-67 dynamically compartmentalizes in VU cell nuclei
Upon closer examination of GFP-tagged NHR-67, it became evident that the AC and VU cells not only exhibit differences in overall NHR-67 levels, but also in localization of the protein.While NHR-67 localization is fairly uniform throughout the AC nucleus (excluding the nucleolus), we often observed discrete punctae throughout the nuclei of VU cells (Figure 3A and B).These punctae were observed with NHR-67 endogenously tagged with several different fluorescent proteins, including GFP, mNeo-nGreen, mScarlet-I, and TagRFP-T (Figure 3-figure supplement 1A and B) and in the absence of tissue fixation methods that can cause artificial puncta (Irgen-Gioro et al., 2022), suggesting that they are not an artifact of the fluorophore or sample preparation.To characterize the dynamics of these punctae during interphase states of the cell cycle, we paired GFP-tagged NHR-67 with a CDK activity sensor.The CDK activity sensor is comprised of a fragment of DNA Helicase B (DHB) fused to a fluorophore (2x-mKate2), expressed under a ubiquitous promoter (Figure 3C; Adikes et al., 2020).DHB contains a strong nuclear localization signal (NLS), flanked by four serine sites, as well as a weaker nuclear export signal (NES).As CDK activity increases over the cell cycle, the CDK sensor is translocated from the nucleus to the cytoplasm, allowing for correlation of its relative subcellular localization to the cell cycle state (Figure 3C; Adikes et al., 2020;Spencer et al., 2013).Time-lapse microscopy revealed that the number of NHR-67 punctae was dynamic over the course of the cell cycle, with punctae first appearing shortly after mitotic exit in the G1 phase, and then reducing in number to two large punctae prior to nuclear envelope breakdown before disappearing (Figure 3D and E).We collected additional recordings with finer time resolution and captured fusion, or condensation, of punctae prior to their dissolution (representative of 6 biological replicates) (Figure 3F).These punctae also exhibit relatively rapid diffusion kinetics, as observed by fluorescence recovery following photobleaching (t 1/2 =46 s; n=8) at a rate within the same order of magnitude as P granule proteins PGL-1 and PGL-3 (Figure 3G and H;Putnam et al., 2019).
Groucho homologs UNC-37 and LSY-22 associate with NHR-67 punctae and contribute to VU cell fate Next, we tested the extent to which NHR-67 punctae colocalize with homologs of other proteins known to compartmentalize in nuclei by pairing GFP-and mScarlet-I-tagged NHR-67 with other endogenously tagged alleles.As NHR-67 is a transcription factor, we speculated that its punctae may represent clustering around sites of active transcription, which would be consistent with data showing RNA Polymerase II and the Mediator complex can compartmentalize with transcription factors (Cho et al., 2018).To test this hypothesis, we co-visualized NHR-67 with a GFP-tagged allele of ama-1, the amanitin-binding subunit of RNA polymerase II (Hills-Muckey et al., 2022) and failed to observe significant colocalization between NHR-67 and AMA-1 punctae (Manders' overlap coefficient, M=0.066) compared to negative controls where a single channel was compared to its 90-degree rotation (M=0.108)(Figure 4A and B).Another possibility considered is that NHR-67 punctae could correlate to chromatin organization, as heterochromatin has been shown to be compartmentalized in the nucleus (Larson et al., 2017;Strom et al., 2017).However, we did not observe significant colocalization of NHR-67 with the endogenously tagged HP1 heterochromatin proteins (Patel and Hobert, 2017) HPL-1 (M=0.076) or HPL-2 (M=0.083)(Figure 4A and B).We next tested if NHR-67 colocalizes with the transcriptional co-repressor Groucho, as Groucho had recently been shown to compartmentalize in the nuclei of cells in Ciona (Treen et al., 2021).The C. elegans genome encodes one Groucho homolog, UNC-37, which we acquired an mNeonGreen-tagged allele of Ma et al., 2021, and a Groucho-like protein, LSY-22, which we tagged with TagRFP-T (Figure 4-figure supplement 1).We observed significant colocalization of NHR-67 punctae with both LSY-22 (M=0.686) and UNC-37 (M=0.741),comparable to colocalization measures in heterozygous NHR-67::mScarlet-I/NHR-67::GFP animals (M=0.651), which were used as positive controls (Figure 4A and B).This evidence suggests that NHR-67 punctae do not localize to sites of active transcription or chromatin compaction, but instead associate with transcriptional co-repressors.
Since the AC is the default state of the AC/VU cell fate decision (Seydoux and Greenwald, 1989), we hypothesized that the punctae including NHR-67, UNC-37, and LSY-22 may function in repressing invasive differentiation.To test this hypothesis, we depleted UNC-37 and LSY-22 utilizing the auxininducible degron (AID) protein degradation system, in which a protein of interest is tagged with an The online version of this article includes the following source data and figure supplement(s) for figure 3: Source data 1.Raw data of CDK sensor (DHB) ratios in ventral uterine (VU) cells over time, as reported in Figure 3D and E.
TCF/LEF homolog POP-1 associates with NHR-67 punctae and contributes to VU cell fate post-specification While UNC-37 and LSY-22 appear to be important for the maintenance of normal uterine cell fates, both genes are broadly expressed and exhibit comparable levels (<10% difference) between the AC and VU cells (Figure 5A and C; Figure 5-figure supplement 1A and B); therefore, we hypothesized that another factor must be involved.It had previously been reported that the sole TCF/LEF homolog in C. elegans, POP-1, forms a repressive complex with UNC-37 in the early embryo to restrict expression of the endoderm-determining gene, END-1 (Calvo et al., 2001).Additionally, POP-1 has a known role in the development of the somatic gonad, as perturbing its function results in ectopic ACs (Siegfried and Kimble, 2002).Examination of an eGFP-tagged pop-1 allele (van der Horst et al., 2019), showed significant enrichment in the VU cells (>25%) compared to the AC (Figure 5B and C; Figure 5-figure supplement 1A and B).We also observed that endogenous POP-1 forms punctae in the nuclei of VU cells, which had previously been observed during interphase in non-Wnt signaled embryonic cells (Maduro et al., 2002).These POP-1 punctae colocalize with NHR-67 (M=0.547),although to a lesser degree than UNC-37 and LSY-22, likely because the strong POP-1 fluorescence outside of punctae made them more difficult to segment (Figure 5D and E).Additionally, NHR-67(RNAi) treatment resulted in a significant increase in AC expression of eGFP::POP-1 compared to empty vector controls (225%, n>30), a pattern we observed following depletion of other transcription factors (Medwig-Kinney et al., 2020) and chromatin modifiers (Smith et al., 2022) required for AC arrest and invasion (Figure 5F and G;.This negative regulation of POP-1 by NHR-67 may explain why the proteins have opposite patterns of enrichment. It has previously been suggested that POP-1 may be functioning as an activator in the VU precursors Z1. ppa and Z4.aap based on the relative expression of a POP-1 transgene (Sallee et al., 2015a).This view is largely dependent on the notion that high levels of POP-1 correlate to repressive function and that low levels are conducive for activator roles (Shetty et al., 2005).In contrast, we did not find evidence of transcriptional activation by POP-1 in the AC/VU precursors nor their differentiated descendants using an established POPTOP (POP-1 and TCF optimal promoter) reporter, which contains seven copies of POP-1/TCF binding sites and the pes-10 minimal promoter (Figure 5H; of the auxin-inducible degron (AID) system, where At TIR1 mediates proteasomal degradation of AID-tagged proteins in the presence of auxin.(D) Representative images of phenotypes observed following individual AID-depletion of UNC-37 and LSY-22 compared to control animals without AID-tagged alleles.All animals compared here are expressing TIR1 ubiquitously (rpl-28p:: At TIR1::T2A::mCherry::HIS-11) and an anchor cell (AC) marker (cdh-3p::mCherry::moeABD). Insets depict different z-planes of the same image.(E) Quantification of AC marker (cdh-3p::mCherry::moeABD) expression in ectopic ACs resulting from AID-depletion of UNC-37 and LSY-22 compared to control AC and VU cells.
The online version of this article includes the following source data and figure supplement(s) for figure 4: Source data 1.Raw data of protein colocalization in ventral uterine (VU) cells, as reported in Figure 4A and B and Figure 5D and E.
Source data 2. Raw data of CDH-3 expression in the ectopic anchor cells (ACs) resulting from auxin-mediated depletion of AID-tagged LSY-22 or UNC-37, as reported in Figure 4D and E.      Figure 5-figure supplement 3A and B; Green et al., 2008).The growing consensus regarding the Wnt/β-catenin asymmetry pathway is that relative levels of POP-1 and β-catenin are more important than absolute protein levels of POP-1 (Phillips and Kimble, 2009).Our proposed model of POP-1 acting as a repressor in the proximal gonad is consistent with the finding that SYS-1 (β-catenin) expression is restricted to the distal gonad early in somatic gonad development and is not detectable in the AC or VU cells (Figure 5-figure supplement 3C; Phillips et al., 2007;Sallee et al., 2015a).It is also supported by recent evidence suggesting that UNC-37/LSY-22 mutant alleles phenocopy pop-1 knockdown, which produces ectopic distal tip cells (Bekas and Phillips, 2022).
One aspect that makes studying the repressive role of POP-1 in cell fate maintenance challenging is that its activator function is required for distal cell fate specification in the somatic gonad earlier in development.Loss of either POP-1 and SYS-1 results in a Sys (symmetrical sister cell) phenotype, where all somatic gonad cells adopt the default proximal fate and thereby give rise to ectopic ACs (Siegfried and Kimble, 2002;Siegfried et al., 2004).This likely occluded previous identification of the potential repressive role of POP-1 in maintaining VU cell fates.To achieve temporal control over POP-1 expression to tease apart its roles, we sought to use the AID system, but the insertion of the degron into the pop-1 locus disrupted the protein's function.Instead, we paired eGFP-tagged POP-1 with a uterine-specific anti-GFP nanobody (Smith et al., 2022;Wang et al., 2017).The anti-GFP nanobody is fused to ZIF-1 and serves as an adapter, recognizing GFP-tagged proteins and promoting their ubiquitination by the Cullin2-based E3 ubiquitin ligase, which ultimately targets them for degradation via the proteasome (Figure 6-figure supplement 1A; Wang et al., 2017).This anti-GFP nanobody, visualized by nuclear expression of mCherry, was not detectable prior to or even shortly after the AC/ VU cell fate decision, which allowed us to bypass disruption of initial cell specification (Figure 6figure supplement 1B).While this method only produced a mild knockdown of POP-1 in the VU cells, we still observed the ectopic AC phenotype at low penetrance (7%, n=60) (Figure 6-figure supplement 1C).To achieve stronger depletion, we used RNAi for further POP-1 perturbations.
To interrogate the phenotypic consequences of POP-1 perturbation, we utilized a strain expressing two markers of AC fate (cdh-3p::mCherry::moeABD and LAG-2::P2A::H2B::mTurquoise2).Following treatment with pop-1(RNAi), we observed several animals with two or more bright cdh-3/lag-2+ ACs, consistent with known phenotypes caused by cell fate misspecification in the somatic gonad (17%, n=30) (Figure 6A).We also observed animals with invasive cells that express AC markers at different levels (53%, n=30), suggesting that the cells did not adopt AC fate at the same time (Figure 6A).To test whether the subset of dim cdh-3/lag-2+ ACs are the result of VU-to-AC cell fate conversion, we visualized AC and VU fates simultaneously using the AC markers previously described along with an mNeonGreen-tagged allele of lag-1 (CSL), a protein downstream of Notch signaling whose expression becomes restricted to the VU cells following AC/VU cell fate specification.Following treatment with pop-1(RNAi), we found that a subset of ectopic ACs co-express AC markers and LAG-1, likely indicating an intermediate state between the two cell types (Figure 6-figure supplement 2).To visualize this process live, we used time-lapse microscopy and were able to capture ectopic ACs gradually upregulating LAG-2 (+51%, n=3) and downregulating LAG-1 (-16%, n=3) over time (Figure 6B and  C), consistent with VU-to-AC cell fate conversion.
To further characterize the protein-protein interaction between NHR-67 and UNC-37, we assessed the protein structure of NHR-67 using AlphaFold, an artificial intelligence-based protein structure prediction tool (Jumper et al., 2021;Varadi et al., 2022), and PONDR, a predictor of intrinsic disorder  (Peng and Zhang, 2006).Both identify an intrinsically disordered region (IDR) at the C-terminus of NHR-67 (Figure 7A and B).IDRs are low complexity domains that lack fixed three-dimensional structures and have been shown to support dynamic protein-protein interactions (Chong et al., 2018).To determine if the IDR of NHR-67 is important for facilitating its interaction with UNC-37, we repeated the yeast two-hybrid experiment using UNC-37 Gal4-AD prey, pairing it with different fragments of the NHR-67 protein: full-length, without its IDR (ΔIDR), and its IDR alone (Figure 7C and D).Yeast growth on the selective SC-HIS-TRP-LEU plates containing the competitive inhibitor 3-aminotriazole (3-AT) demonstrates that the 108 amino acid IDR sequence of NHR-67 is necessary and to bind with UNC-37 (Figure 7C and D).
Thus, we propose the following potential model of C. elegans uterine cell fate maintenance based on the data presented here and the known roles of Groucho and TCF proteins in regulating transcription.First, transcription of NHR-67 is directly regulated by HLH-2, resulting in its enrichment in the AC compared to the VU cells.In the AC, where NHR-67 levels are high and POP-1 is repressed, NHR-67 is free to activate genes promoting invasive differentiation.In the VU cells, where NHR-67 levels are low and POP-1 levels are high, POP-1 assembles with LSY-22, UNC-37, and NHR-67, either directly repressing NHR-67 targets or sequestering NHR-67 away from its targets (Figure 7E).It is possible that POP-1 negatively regulates NHR-67 at the transcriptional level as well, as the NHR-67 promoter contains seven putative TCF binding sites (Zacharias et al., 2015).

Discussion
In summary, here we provide evidence that the activity of the pro-invasive transcription factor, NHR-67, is simultaneously regulated by two distinct processes, which together modulate the proliferativeinvasive switch in C. elegans.We show that NHR-67 is a potent fate-specifying transcription factor, in that its expression is sufficient for the invasive differentiation of ACs in the somatic gonad.The compartmentalization of NHR-67 in the VU cells could serve as a potential mechanism to suppress its function in activating the pro-invasive program.We also discovered that NHR-67 forms nuclear foci in non-invasive cells, which exhibit liquid-like properties, indicated by observations of their condensation, dissolution, and relatively rapid recovery from photobleaching, similar to what has been described with P granules (Brangwynne et al., 2009).These NHR-67 punctae associate with Groucho homologs, UNC-37 and LSY-22, likely through a direct protein-protein interaction with UNC-37 mediated by the C-terminal IDR of NHR-67.We postulate that this association leads to protein condensation, as has recently been described in Ciona embryos (Treen et al., 2021).Furthermore, repression of the default invasive state appears to be dependent on the expression of the TCF/LEF homolog POP-1, which clarifies our understanding of the dual roles this protein plays during the development of the somatic gonad.It is also interesting to note that the dynamic punctae formed by POP-1 in non-Wnt signaled cells was first described 20 years ago (Maduro et al., 2002), but their function is only now being appreciated in light of recent advances in our understanding of the formation of higher-order associations in the nucleus.
With regard to protein compartmentalization in the nucleus, most research has been through the lens of transcriptional activation through RNA Polymerase II and the mediator complex (Boija et al., 2018;Cho et al., 2018;Sabari et al., 2018) or repression through HP1 heterochromatin proteins (Larson et al., 2017;Strom et al., 2017).Here, we report the second observed case of compartmentalization of Groucho proteins (Treen et al., 2021), which may suggest that Groucho proteins have evolutionarily conserved roles that require this type of subnuclear organization.The online version of this article includes the following source data and figure supplement(s) for figure 6: Source data 1.Raw data of LAG-2::P2A::H2B::mT2 and LAG-1::mNG expression during VU-to-AC transdifferentiation, as reported in Figure 6B and C.   Still, as this is one of the first studies into the compartmentalization of transcriptional repressors in vivo, there is much left to learn.For example, it is unknown whether DNA binding is necessary for nuclear puncta formation.The interaction between UNC-37 and NHR-67 does not appear to depend on DNA binding, as the C-terminal IDR region of NHR-67 (excluding its zinc finger domains) was sufficient for binding with UNC-37 in vitro, but it is possible that DNA binding is needed for oligomerization in vivo.Furthermore, it remains unclear if suppression of invasive differentiation is achieved by simply sequestering the pro-invasive transcription factor NHR-67 away from its transcriptional targets or through direct repression of transcription.If the latter, another question that arises is how the repressive complex gets recruited to specific genomic sites, since POP-1 and NHR-67 are both capable of binding to DNA, and whether repression is achieved through competition with transcriptional activators or recruitment of histone deacetylases.Direct targets of NHR-67 have not yet been discovered, which makes it difficult to investigate this specific aspect of the repressive mechanism at present.We see this as a promising avenue of future study as technologies advance, allowing for transcriptional profiling and target identification in specific tissues or cells (Gómez-Saldivar et al., 2020;Katsanos and Barkoulas, 2022).
In this work, we have also identified several perturbations (i.e.increasing levels of NHR-67, decreasing levels of UNC-37/LSY-22) that result in incompletely penetrant transdifferentiation phenotypes and/or intermediate cell fates.We foresee these being ideal cell fate challenge backgrounds in which to perform screens to identify regulators of cellular plasticity, as has been done in other contexts (Rahe and Hobert, 2019).Additionally, these induced fate transformations can be paired with tools to visualize and manipulate the cell cycle (Adikes et al., 2020) to determine if any cell cycle state is particularly permissive for cell fate plasticity.While G1 arrest has been shown to enhance the conversion of human fibroblasts to dopaminergic neurons (Jiang et al., 2015), mitosis is required for the natural K-to-DVB transdifferentiation event in C. elegans (Riva et al., 2022).As control of proliferation and invasion, as well as maintenance of differentiated cellular identities, are important for both homeostatic and disease states, it is our hope that this work will shed light on how cells switch between these states in the context of cancer growth and metastasis.

Materials and methods
C. elegans strains, culture, and nomenclature Methods for C. elegans culture and genetics were followed as previously described (Brenner, 1974).Developmental synchronization for experiments was achieved through alkaline hypochlorite treatment of gravid adults to isolate eggs (Porta-de-la- Riva et al., 2012).L1 stage animals were plated on nematode growth media plates and subsequently cultured at 20 °C or 25 °C.Heat shock-inducible transgenes were activated by incubating animals on plates sealed with Parafilm in a 33 °C water bath for 2-3 hr.In the text and in figures, promoter sequences are designated with a 'p' following the gene name and gene fusions are represented by a double-colon (::) symbol.

CRISPR/Cas9 injections
New alleles and single-copy transgenes were generated by homology-directed repair using CRISPRbased genome engineering.mScarlet::AID and mNeonGreen::AID were inserted into the C-terminus of the NHR-67 locus by injecting adult germlines with Cas9 guide-RNA ribonucleoprotein complexes and short single-stranded oligodeoxynucleotide donors, as previously described (Ghanta and Mello, 2020).Successful integration was identified through screening for fluorescence and PCR.The LSY-22 locus was edited by injecting a Cas9 guide RNA plasmid and repair template plasmid containing a selfexcising cassette with selectable markers to facilitate screening (Dickinson et al., 2015;Dickinson allows for yeast growth in the presence of competitive inhibitor 3-AT (20 mM).(E) Possible models of the roles of NHR-67, UNC-37, LSY-22, and POP-1 in the maintenance of anchor cell (AC) and ventral uterine (VU) cell fate.In the ventral uterine cells, the association of NHR-67 with the Groucho/TCF complex may result in the repression of NHR-67 targets (top) or the sequestration of NHR-67 away from its targets (bottom).
The online version of this article includes the following figure supplement(s) for figure 7: Figure supplement 1. NHR-67 exhibits protein-protein interaction with UNC-37.

Auxin inducible protein degradation
The auxin-inducible degron (AID) system was used for the strong depletion of proteins of interest (Zhang et al., 2015).AID-tagged alleles were paired with the Arabidopsis thaliana F-box protein, transport inhibitor response 1 ( At TIR1), and treated with the water-soluble auxin 1-Naphthaleneacetic acid (K-NAA) at 1 mM concentration (Martinez et al., 2020).Auxin was added to nematode growth media plates according to previously published protocols (Martinez and Matus, 2020), which were then seeded with OP50 E. coli.To achieve robust depletion, synchronized L1 stage animals were directly plated on auxin plates.

Live-cell imaging
With the exception of the FRAP experiments shown in Figure 3, all micrographs were collected on a Hamamatsu Orca EM-CCD camera mounted on an upright Zeiss AxioImager A2 with a Borealismodified CSU10 Yokagawa spinning disk scan head (Nobska Imaging) using 405 nm, 440 nm, 488 nm, 514 nm, and 561 nm Vortran lasers in a VersaLase merge and a Plan-Apochromat 100x/1.4(NA) Oil DIC objective.MetaMorph software (Molecular Devices) was used for microscopy automation.Several experiments were scored using epifluorescence visualized on a Zeiss Axiocam MRM camera, also mounted on an upright Zeiss AxioImager A2 and a Plan-Apochromat 100x/1.4(NA) Oil DIC objective.For static imaging, animals were mounted into a drop of M9 on a 5% Noble agar pad containing approximately 10 mM sodium azide anesthetic and topped with a coverslip.For long-term timelapse imaging, animals were first anesthetized in 5 mM levamisole diluted in M9 for approximately 20 min, then transferred to a 5% Noble agar pad and topped with a coverslip sealed with VALAP (Kelley et al., 2017).For short-term time-lapse imaging, the pre-anesthetization step was omitted, and animals were transferred directly into a drop of 5 mM levamisole solution on the slide.

Fluorescence recovery after photobleaching
FRAP experiments were performed using an Acal BFi UV Optimicroscan photostimulation device mounted on a spinning disk confocal system consisting of a Nikon Ti2 inverted microscope with Yokogawa CSU-W1 SoRa spinning disk.Data were acquired using a Hamamatsu ORCA Fusion camera, 60x 1.27 NA water immersion objection, SoRa disk, and 2.8 x SoRa magnifier.Single plane images were collected every 1 s.

Yeast one-hybrid
The 276 bp fragment of the NHR-67 promoter (Bodofsky et al., 2018) was cloned into the pMW2 vector, and linearized by BamHI digestion.Linearized plasmid was transformed into the Y1H yeast strain (as described in Reece-Hoyes and Walhout, 2018).Transformed yeast was plated on SC-HIS plates for three days before being transformed with the HLH-2 Gal4-AD plasmid.Three colonies from each transformation plate were streaked onto SC-HIS-TRP+3-aminotriazole (3-AT) plates.Protein-DNA interactions were determined by visible growth on 3-AT conditions with negative growth in empty vector controls after three days.Plates were imaged on a Fotodyne FOTO/Analyst Investigator/FX darkroom imaging station.

Yeast two-hybrid
Plasmids containing target proteins fused to GAL-4 DNA-binding-domain + LEU and GAL-4 Activation Domain + TRP were co-transformed into the pJ69-4a Y2H yeast strain as previously described (Reece-Hoyes and Walhout, 2018).Transformed yeast was plated on SC-TRP-LEU plates for three days.Three colonies from each transformation plate were streaked onto SC-HIS-TRP-LEU 3-AT plates.Protein interactions were determined by visible growth on 3-AT conditions with negative growth in empty vector controls after three days.Plates were imaged as described in the previous section.

Quantification of protein expression and cell cycle state
Image quantification was performed in Fiji/ImageJ (Schindelin et al., 2012).Protein expression was quantified by drawing a region of interest and measuring the mean gray value, then manually subtracting the mean gray value of a background region of similar area to account for camera noise.For nuclear-localized proteins, the region of interest was drawn around the nucleus.Membrane expression of LIN-12 was measured by drawing a line (1.5 pixel thickness) along a cell's basolateral surface, to distinguish between LIN-12 expression from the cell of interest and its neighbors.CDH-3 expression was measured by drawing a region of interest around the cell membrane excluding the nucleus (as the nuclear-localized TIR1 transgene was expressed in the same channel).The CDK sensor was quantified as previously described (Adikes et al., 2020).Following rolling ball subtraction (50 pixels), the mean gray value is measured in a region of interest drawn within the cytoplasm and one around the nucleus excluding the nucleolus.The cytoplasmic-to-nuclear ratio correlates to CDK activity and is used to assess the cell cycle state (Adikes et al., 2020;Spencer et al., 2013).Movies were collected by acquiring z-stacks at 5 min intervals.Samples were time-aligned relative to anaphase.Cells that did not undergo anaphase during the acquisition period were aligned based on their DHB ratios.Animals that were arrested in development (i.e. did not show evidence of progressing through the cell cycle) were excluded from the analysis.

Colocalization analyses
For colocalization analyses, single-plane images were collected to avoid z drift during acquisition and prevent photobleaching, which was often non-uniform between red and green fluorophores.Micrographs were subject to background subtraction (rolling ball radius = 50) followed by thresholding to segment punctae.Manders' overlap coefficients (M) were calculated by measuring the extent that segmented punctae of NHR-67 overlapped with that of other proteins using Just Another Colocalization Plugin (JACoP) in Fiji/ImageJ (Bolte and Cordelières, 2006;Schindelin et al., 2012).Heterozygous animals for nhr-67::mScarlet and nhr-67::GFP were used as positive controls.These images were then re-analyzed following a 90-degree rotation of one of the two channels being compared, resulting in random colocalization that served as a negative control.

Data visualization and statistical analyses
Representative images were processed using Fiji/ImageJ (Schindelin et al., 2012).Heat maps were generated using the Fire lookup table.A power analysis was performed prior to data collection to determine sample sizes (Cohen, 1992).Tests to determine the statistical significance of data were conducted in RStudio and plots were generated using the R package ggplot2 (Wickham, 2016).Error bars represent the mean ± standard deviation.Schematics of gene loci were generated using sequences from WormBase (Harris et al., 2020) Matus et al., 2015 qyIs225[cdh-3p::mCherry::moeABD] V; qyIs7[laminin::GFP] X.

Figure supplement 1 .
Figure supplement 1.Expression of pro-invasive transcription factors EGL-43 and FOS-1 in the somatic gonad.

Figure supplement 1 .
Figure supplement 1. Onset of expression and regulatory interaction between NHR-67 and HLH-2 in the somatic gonad.

Figure 3 .
Figure 3 over time, aligned to anaphase.(E) Representative time-lapse of NHR-67::GFP over the course of a cell cycle, with cell membranes indicated with dotted lines.(F) Time-lapse depicting NHR-67::GFP punctae fusion prior to cell division.Right panels are pseudo-colored.(G-H) Representative images (G) and quantification (H) depicting fluorescence recovery of NHR-67::GFP following photobleaching of individual punctae (arrow).

Figure supplement 2 .
Figure supplement 2. POP-1 function in ventral uterine (VU) cells is distinct from the activating role in distal somatic gonad.

Figure supplement 3 .
Figure supplement 3. POP-1 expression is regulated by the cell cycle-dependent pro-invasion pathway.