Signal transduction during C. elegans vulval development: a NeverEnding story

Tissue-specific inhibition of protein sumoylation uncovers diverse SUMO functions during C. elegans vulval development

The sumoylation (SUMO) pathway is involved in a variety of processes during C. elegans development, such as gonadal and vulval fate specification, cell cycle progression and maintenance of chromosome structure. The ubiquitous expression and pleiotropic effects have made it difficult to dissect the tissue-specific functions of the SUMO pathway and identify its target proteins. To overcome these challenges, we have established tools to block protein sumoylation and degrade sumoylated target proteins in a tissue-specific and temporally controlled manner. We employed the auxin-inducible protein degradation system (AID) to down-regulate the SUMO E3 ligase GEI-17 or the SUMO ortholog SMO-1, either in the vulval precursor cells (VPCs) or in the gonadal anchor cell (AC). Our results indicate that the SUMO pathway acts in multiple tissues to control different aspects of vulval development, such as AC positioning, basement membrane (BM) breaching, VPC fate specification and morphogenesis. Inhibition of protein sumoylation in the VPCs resulted in abnormal toroid formation and ectopic cell fusions during vulval morphogenesis. In particular, sumoylation of the ETS transcription factor LIN-1 at K169 is necessary for the proper contraction of the ventral vulA toroids. Thus, the SUMO pathway plays several distinct roles throughout vulval development.

Many proteins are chemically modified after they have been synthesized. In particular, conjugation with the Small Ubiquitin-like Modifier (SUMO) regulates the functions and activities of a large number of proteins in animal and plant cells.

Here, we have used the Nematode Caenorhabditis elegans to study the various effects of SUMO protein modification on organ development. By applying a tissue-specific protein degradation system, we could selectively block the SUMO pathway in different tissues of the animals. We focused on the development of the egg-laying organ as a model, and found that the SUMO pathway acts in multiple tissues to regulate distinct cellular functions. Finally, we show that SUMO modification of one transcription factor, called LIN-1, is necessary for the proper morphogenesis of the organ. Our results indicate that the manifold effects of the SUMO pathway can be attributed to the combined action of a distinct number of SUMO modified proteins acting in different cell types.

The retromer complex regulates C. elegans development and mammalian ciliogenesis

The mammalian retromer consists of subunits VPS26 (either VPS26A or VPS26B), VPS29 and VPS35, and a loosely associated sorting nexin (SNX) heterodimer or a variety of other SNX proteins. Despite involvement in yeast and mammalian cell trafficking, the role of retromer in development is poorly understood, and its impact on primary ciliogenesis remains unknown. Using CRISPR/Cas9 editing, we demonstrate that vps-26-knockout worms have reduced brood sizes, impaired vulval development and decreased body length, all of which have been linked to ciliogenesis defects. Although preliminary studies did not identify worm ciliary defects, and impaired development limited additional ciliogenesis studies, we turned to mammalian cells to investigate the role of retromer in ciliogenesis. VPS35 localized to the primary cilium of mammalian cells, and depletion of VPS26, VPS35, VPS29, SNX1, SNX2, SNX5 or SNX27 led to decreased ciliogenesis. Retromer also coimmunoprecipitated with the centriolar protein, CP110 (also known as CCP110), and was required for its removal from the mother centriole. Herein, we characterize new roles for retromer in C. elegans development and in the regulation of ciliogenesis in mammalian cells, suggesting a novel role for retromer in CP110 removal from the mother centriole.

Nfya-1 functions as a substrate of ERK-MAP kinase during Caenorhabditis elegans vulval development

A common bridge between a linear cytoplasmic signal and broad nuclear regulation is the family of MAP kinases which can translocate to the nucleus upon activation by the cytoplasmic signal. One pathway which functions to activate the ERK family of MAP kinases is the Ras signaling pathway which functions at multiple times and locations during the development of Caenorhabditis elegans including the development of the excretory cell, germ cells, male tail, and vulva. It has been most extensively characterized during the development of the vulva which is formed from the vulval precursor cells (VPCs), a set of six equivalent, epithelial cells designated P3.p - P8.p. Although LIN-1 appears to be a primary target of ERK MAP kinase during vulval development, it is likely that other developmentally important molecules are also regulated by ERK-mediated phosphorylation. The identification of physiological substrates of MAP kinases has been aided by the identification of docking site domains in substrate proteins that contribute to high-affinity interactions with kinases. Our laboratory has identified the C. elegans protein, T08D10.1/NFYA-1, as a potential ERK MAP kinase substrate in this manner, and we have initiated a characterization of its role during Ras-mediated development. T08D10.1 possesses significant homology to the CCAAT-box DNA-binding domain of the vertebrate nuclear transcription factor-Y, alpha (NF-YA) family of proteins. NF-Y proteins act as part of a complex to regulate the transcription of a large number of genes, in particular, genes that function in the G1/S cell cycle transition. T08D10.1/NFYA-1 is predicted to code for a protein containing multiple potential phosphorylation sites for ERK MAP kinase and a D-domain docking site. We demonstrate through biochemical analysis of purified NFYA-1 protein that it can act in vitro as a high affinity substrate for activated ERK MAP kinase. Growth factor activation of the Ras pathway in a tissue culture system has negligible effect on the protein's transactivation potential, however, the DNA-binding activity of the protein is reduced after treatment with activated ERK-MAP kinase. We demonstrate through mutant analysis that nfya-1 acts to inhibit vulval development and functions downstream or in parallel to let-60/ras. Both the NF-Y complex and the Ras signaling pathway play a fundamental role in cell proliferation and oncogenesis and the connection between the two is an important insight into the mechanisms of cell fate specification and cellular response.

Nets-within-nets for modeling emergent patterns in ontogenetic processes

Ontogenesis is the development of an organism from its earliest stage to maturity, including homeostasis maintenance throughout adulthood despite environmental perturbations. Almost all cells of a multicellular organism share the same genomic information. Nevertheless, phenotypic diversity and complex supra-cellular architectures emerge at every level, starting from tissues and organs. This is possible thanks to a robust and dynamic interplay of regulative mechanisms. To study ontogenesis, it is necessary to consider different levels of regulation, both genetic and epigenetic. Each cell undergoes a specific path across a landscape of possible regulative states affecting both its structure and its functions during development. This paper proposes using the Nets-Within-Nets formalism, which combines Petri Nets' simplicity with the capability to represent and simulate the interplay between different layers of regulation connected by non-trivial and context-dependent hierarchical relations. In particular, this work introduces a modeling strategy based on Nets-Within-Nets that can model several critical processes involved in ontogenesis. Moreover, it presents a case study focusing on the first phase of Vulval Precursor Cells specification in C.Elegans. The case study shows that the proposed model can simulate the emergent morphogenetic pattern corresponding to the observed developmental outcome of that phase, in both the physiological case and different mutations. The model presented in the results section is available online at https://github.com/sysbio-polito/NWN_CElegans_VPC_model/.

Microfluidic-based imaging of complete Caenorhabditis elegans larval development

Several microfluidic-based methods for Caenorhabditis elegans imaging have recently been introduced. Existing methods either permit imaging across multiple larval stages without maintaining a stable worm orientation, or allow for very good immobilization but are only suitable for shorter experiments. Here, we present a novel microfluidic imaging method that allows parallel live-imaging across multiple larval stages, while maintaining worm orientation and identity over time. This is achieved through an array of microfluidic trap channels carefully tuned to maintain worms in a stable orientation, while allowing growth and molting to occur. Immobilization is supported by an active hydraulic valve, which presses worms onto the cover glass during image acquisition only. In this way, excellent quality images can be acquired with minimal impact on worm viability or developmental timing. The capabilities of the devices are demonstrated by observing the hypodermal seam and P-cell divisions and, for the first time, the entire process of vulval development from induction to the end of morphogenesis. Moreover, we demonstrate feasibility of on-chip RNAi by perturbing basement membrane breaching during anchor cell invasion.

Summary: Parallel microfluidic long-term imaging allows reliable long-term study of Caenorhabditis elegans development across multiple larval stages at high-resolution and with minimal effect on physiological development.

Quantifying cell transitions in C. elegans with data-fitted landscape models

Increasing interest has emerged in new mathematical approaches that simplify the study of complex differentiation processes by formalizing Waddington's landscape metaphor. However, a rational method to build these landscape models remains an open problem. Here we study vulval development in C. elegans by developing a framework based on Catastrophe Theory (CT) and approximate Bayesian computation (ABC) to build data-fitted landscape models. We first identify the candidate qualitative landscapes, and then use CT to build the simplest model consistent with the data, which we quantitatively fit using ABC. The resulting model suggests that the underlying mechanism is a quantifiable two-step decision controlled by EGF and Notch-Delta signals, where a non-vulval/vulval decision is followed by a bistable transition to the two vulval states. This new model fits a broad set of data and makes several novel predictions.

Programmed cell fusion in development and homeostasis

During multicellular organism development, complex structures are sculpted to form organs and tissues, which are maintained throughout adulthood. Many of these processes require cells to fuse with one another, or with themselves. These plasma membrane fusions merge endoplasmic cellular content across external, exoplasmic, space. In the nematode Caenorhabditis elegans, such cell fusions serve as a unique sculpting force, involved in the embryonic morphogenesis of the skin-like multinuclear hypodermal cells, but also in refining delicate structures, such as valve openings and the tip of the tail. During post-embryonic development, plasma membrane fusions continue to shape complex neuron structures and organs such as the vulva, while during adulthood fusion participates in cell and tissue repair. These processes rely on two fusion proteins (fusogens): EFF-1 and AFF-1, which are part of a broader family of structurally related membrane fusion proteins, encompassing sexual reproduction, viral infection, and tissue remodeling. The established capabilities of these exoplasmic fusogens are further expanded by new findings involving EFF-1 and AFF-1 in endocytic vesicle fission and phagosome sealing. Tight regulation by cell-autonomous and non-cell autonomous mechanisms orchestrates these diverse cell fusions at the correct place and time-these processes and their significance are discussed in this review.

A multi-layered and dynamic apical extracellular matrix shapes the vulva lumen in Caenorhabditis elegans

Biological tubes must develop and maintain their proper diameter to transport materials efficiently. These tubes are molded and protected in part by apical extracellular matrices (aECMs) that line their lumens. Despite their importance, aECMs are difficult to image in vivo and therefore poorly understood. The Caenorhabditis elegans vulva has been a paradigm for understanding many aspects of organogenesis. Here we describe the vulva luminal matrix, which contains chondroitin proteoglycans, Zona Pellucida (ZP) domain proteins, and other glycoproteins and lipid transporters related to those in mammals. Confocal and transmission electron microscopy revealed, with unprecedented detail, a complex and dynamic aECM. Different matrix factors assemble on the apical surfaces of each vulva cell type, with clear distinctions seen between Ras-dependent (1°) and Notch-dependent (2°) cell types. Genetic perturbations suggest that chondroitin and other aECM factors together generate a structured scaffold that both expands and constricts lumen shape.

Assessment and Maintenance of Unigametic Germline Inheritance for C. elegans

The recent work of Besseling and Bringmann (2016) identified a molecular intervention for C. elegans in which premature segregation of maternal and paternal chromosomes in the fertilized oocyte can produce viable animals exhibiting a non-Mendelian inheritance pattern. Overexpression in embryos of a single protein regulating chromosome segregation (GPR-1) provides a germline derived clonally from a single parental gamete. We present a collection of strains and cytological assays to consistently generate and track non-Mendelian inheritance. These tools allow reproducible and high-frequency (>80%) production of non-Mendelian inheritance, the facile and simultaneous homozygosis for all nuclear chromosomes in a single generation, the precise exchange of nuclear and mitochondrial genomes between strains, and the assessments of non-canonical mitosis events. We show the utility of these strains by demonstrating a rapid assessment of cell lineage requirements (AB versus P1) for a set of genes (lin-2, lin-3, lin-12, and lin-31) with roles in C. elegans vulval development.

Cell-Cell Contact and Receptor Tyrosine Kinase Signaling

The behavior of cells within tissues is governed by the activities of adhesion receptors that provide spatial cues and transmit forces through intercellular junctions, and by growth-factor receptors, particularly receptor tyrosine kinases (RTKs), that respond to biochemical signals from the environment. Coordination of these two activities is essential for the patterning and polarized migration of cells during morphogenesis and for homeostasis in mature tissues; loss of this coordination is a hallmark of developing cancer and driver of metastatic progression. Although much is known about the individual functions of adhesion and growth factor receptors, we have a surprisingly superficial understanding of the mechanisms by which their activities are coordinated.

Functional Dysregulation of CDC42 Causes Diverse Developmental Phenotypes

Exome sequencing has markedly enhanced the discovery of genes implicated in Mendelian disorders, particularly for individuals in whom a known clinical entity could not be assigned. This has led to the recognition that phenotypic heterogeneity resulting from allelic mutations occurs more commonly than previously appreciated. Here, we report that missense variants in CDC42, a gene encoding a small GTPase functioning as an intracellular signaling node, underlie a clinically heterogeneous group of phenotypes characterized by variable growth dysregulation, facial dysmorphism, and neurodevelopmental, immunological, and hematological anomalies, including a phenotype resembling Noonan syndrome, a developmental disorder caused by dysregulated RAS signaling. In silico, in vitro, and in vivo analyses demonstrate that mutations variably perturb CDC42 function by altering the switch between the active and inactive states of the GTPase and/or affecting CDC42 interaction with effectors, and differentially disturb cellular and developmental processes. These findings reveal the remarkably variable impact that dominantly acting CDC42 mutations have on cell function and development, creating challenges in syndrome definition, and exemplify the importance of functional profiling for syndrome recognition and delineation.

The hallmarks of cell-cell fusion

Cell-cell fusion is essential for fertilization and organ development. Dedicated proteins known as fusogens are responsible for mediating membrane fusion. However, until recently, these proteins either remained unidentified or were poorly understood at the mechanistic level. Here, we review how fusogens surmount multiple energy barriers to mediate cell-cell fusion. We describe how early preparatory steps bring membranes to a distance of ∼10 nm, while fusogens act in the final approach between membranes. The mechanical force exerted by cell fusogens and the accompanying lipidic rearrangements constitute the hallmarks of cell-cell fusion. Finally, we discuss the relationship between viral and eukaryotic fusogens, highlight a classification scheme regrouping a superfamily of fusogens called Fusexins, and propose new questions and avenues of enquiry.

Developmental ERK Signaling Goes Digital

Reporting in Developmental Cell, de la Cova et al. (2017) present a biosensor to measure ERK activity dynamics in C. elegans larvae. They find that fate decision signaling involves frequency-modulated, digital ERK activity pulses. These findings may explain how graded morphogen signals are converted into precise and robust cell fate patterns.

Studying Nonproliferative Roles for Egfr Signaling in Tissue Morphogenesis Using Dorsal Closure of the Drosophila Embryo

For several decades, genetic analysis in Drosophila has made important contributions to the understanding of signaling by Egfr. Egfr has been well characterized with regard to its oncogenic potential but is also being studied for its roles in organismal development. We have recently developed dorsal closure of the Drosophila embryo as a system for characterizing Egfr regulation of events that do not involve proliferation, as no cell divisions occur during this process. Dorsal closure is essentially a developmental wound healing event with parallels to vertebrate developmental epithelial fusions such as neural tube closure and palate fusion. We describe here a set of materials and protocols for studying Egfr signaling during dorsal closure, including assessing effects of altering Egfr signaling on other pathways, gene expression and, using live imaging, morphogenesis and programmed cell death. Although this "tool kit" is designed for looking at Egfr, it can be readily adapted to look at the participation of any signaling molecule in dorsal closure.

C. elegans SMA-10 regulates BMP receptor trafficking

Signal transduction of the conserved transforming growth factor-β (TGFβ) family signaling pathway functions through two distinct serine/threonine transmembrane receptors, the type I and type II receptors. Endocytosis orchestrates the assembly of signaling complexes by coordinating the entry of receptors with their downstream signaling mediators. Recently, we showed that the C. elegans type I bone morphogenetic protein (BMP) receptor SMA-6, part of the TGFβ family, is recycled through the retromer complex while the type II receptor, DAF-4 is recycled in a retromer-independent, ARF-6 dependent manner. From genetic screens in C. elegans aimed at identifying new modifiers of BMP signaling, we reported on SMA-10, a conserved LRIG (leucine-rich and immunoglobulin-like domains) transmembrane protein. It is a positive regulator of BMP signaling that binds to the SMA-6 receptor. Here we show that the loss of sma-10 leads to aberrant endocytic trafficking of SMA-6, resulting in its accumulation in distinct intracellular endosomes including the early endosome, multivesicular bodies (MVB), and the late endosome with a reduction in signaling strength. Our studies show that trafficking defects caused by the loss of sma-10 are not universal, but affect only a limited set of receptors. Likewise, in Drosophila, we find that the fly homolog of sma-10, lambik (lbk), reduces signaling strength of the BMP pathway, consistent with its function in C. elegans and suggesting evolutionary conservation of function. Loss of sma-10 results in reduced ubiquitination of the type I receptor SMA-6, suggesting a possible mechanism for its regulation of BMP signaling.

Wnt Signaling Polarizes C. elegans Asymmetric Cell Divisions During Development

Asymmetric cell division is a common mode of cell differentiation during the invariant lineage of the nematode, C. elegans. Beginning at the four-cell stage, and continuing throughout embryogenesis and larval development, mother cells are polarized by Wnt ligands, causing an asymmetric inheritance of key members of a Wnt/β-catenin signal transduction pathway termed the Wnt/β-catenin asymmetry pathway. The resulting daughter cells are distinct at birth with one daughter cell activating Wnt target gene expression via β-catenin activation of TCF, while the other daughter displays transcriptional repression of these target genes. Here, we seek to review the body of evidence underlying a unified model for Wnt-driven asymmetric cell division in C. elegans, identify global themes that occur during asymmetric cell division, as well as highlight tissue-specific variations. We also discuss outstanding questions that remain unanswered regarding this intriguing mode of asymmetric cell division.

Caenorhabditis elegans as a model for understanding ROS function in physiology and disease

ROS (reactive oxygen species) are potentially damaging by-products of aerobic metabolism which, unchecked, can have detrimental effects on cell function. However, it is now widely accepted that, at physiological levels, certain ROS play important roles in cell signaling, acting as second messengers to regulate cell choices that contribute to the development, adaptation and survival of plants and animals. Despite important recent advances in the biochemical tools available to study redox-signaling, the molecular mechanisms underlying most of these responses remain poorly understood, particularly in multicellular organisms. As we will review here, C. elegans has emerged as a powerful animal model to elucidate these and other aspects of redox biology.

Dynein-mediated trafficking negatively regulates LET-23 EGFR signaling

Epidermal Growth Factor Receptor (EGFR) signaling is essential for animal development and increased signaling underlies many human cancers. Identifying the genes and cellular processes that regulate EGFR signaling in vivo will help elucidate how this pathway can become inappropriately activated. Caenorhabditis elegans vulva development provides an in vivo model to genetically dissect EGFR signaling. Here we identified a mutation in dhc-1, the heavy chain of the cytoplasmic dynein minus-end directed microtubule motor, in a genetic screen for regulators of EGFR signaling. Despite the many cellular functions of dynein, DHC-1 is a strong negative regulator of EGFR signaling during vulva induction. DHC-1 is required in the signal-receiving cell, genetically functions upstream or in parallel to LET-23 EGFR. LET-23 EGFR accumulates in cytoplasmic foci in dhc-1 mutants consistent with mammalian cell studies whereby dynein has been shown to regulate late endosome trafficking of EGFR with the Rab7 GTPase. However, we found different distributions of LET-23 EGFR foci in rab-7 versus dhc-1 mutants, suggesting that dynein functions at an earlier step of LET-23 EGFR trafficking to the lysosome than RAB-7. Our results demonstrate an in vivo role for dynein in limiting LET-23 EGFR signaling via endosomal trafficking.

Natural Genetic Variation Influences Protein Abundances in C. elegans Developmental Signalling Pathways

Complex traits, including common disease-related traits, are affected by many different genes that function in multiple pathways and networks. The apoptosis, MAPK, Notch, and Wnt signalling pathways play important roles in development and disease progression. At the moment we have a poor understanding of how allelic variation affects gene expression in these pathways at the level of translation. Here we report the effect of natural genetic variation on transcript and protein abundance involved in developmental signalling pathways in Caenorhabditis elegans. We used selected reaction monitoring to analyse proteins from the abovementioned four pathways in a set of recombinant inbred lines (RILs) generated from the wild-type strains N2 (Bristol) and CB4856 (Hawaii) to enable quantitative trait locus (QTL) mapping. About half of the cases from the 44 genes tested showed a statistically significant change in protein abundance between various strains, most of these were however very weak (below 1.3-fold change). We detected a distant QTL on the left arm of chromosome II that affected protein abundance of the phosphatidylserine receptor protein PSR-1, and two separate QTLs that influenced embryonic and ionizing radiation-induced apoptosis on chromosome IV. Our results demonstrate that natural variation in C. elegans is sufficient to cause significant changes in signalling pathways both at the gene expression (transcript and protein abundance) and phenotypic levels.

The Mediator Kinase Module Restrains Epidermal Growth Factor Receptor Signaling and Represses Vulval Cell Fate Specification in Caenorhabditis elegans

Cell signaling pathways that control proliferation and determine cell fates are tightly regulated to prevent developmental anomalies and cancer. Transcription factors and coregulators are important effectors of signaling pathway output, as they regulate downstream gene programs. In Caenorhabditis elegans, several subunits of the Mediator transcriptional coregulator complex promote or inhibit vulva development, but pertinent mechanisms are poorly defined. Here, we show that Mediator's dissociable cyclin dependent kinase 8 (CDK8) module (CKM), consisting of cdk-8, cic-1/Cyclin C, mdt-12/dpy-22, and mdt-13/let-19, is required to inhibit ectopic vulval cell fates downstream of the epidermal growth factor receptor (EGFR)-Ras-extracellular signal-regulated kinase (ERK) pathway. cdk-8 inhibits ectopic vulva formation by acting downstream of mpk-1/ERK, cell autonomously in vulval cells, and in a kinase-dependent manner. We also provide evidence that the CKM acts as a corepressor for the Ets-family transcription factor LIN-1, as cdk-8 promotes transcriptional repression by LIN-1. In addition, we find that CKM mutation alters Mediator subunit requirements in vulva development: the mdt-23/sur-2 subunit, which is required for vulva development in wild-type worms, is dispensable for ectopic vulva formation in CKM mutants, which instead display hallmarks of unrestrained Mediator tail module activity. We propose a model whereby the CKM controls EGFR-Ras-ERK transcriptional output by corepressing LIN-1 and by fine tuning Mediator specificity, thus balancing transcriptional repression vs. activation in a critical developmental signaling pathway. Collectively, these data offer an explanation for CKM repression of EGFR signaling output and ectopic vulva formation and provide the first evidence of Mediator CKM-tail module subunit crosstalk in animals.

Genes conserved in bilaterians but jointly lost with Myc during nematode evolution are enriched in cell proliferation and cell migration functions

Animals use a stereotypical set of developmental genes to build body architectures of varying sizes and organizational complexity. Some genes are critical to developmental patterning, while other genes are important to physiological control of growth. However, growth regulator genes may not be as important in small-bodied “micro-metazoans” such as nematodes. Nematodes use a simplified developmental strategy of lineage-based cell fate specifications to produce an adult bilaterian body composed of a few hundreds of cells. Nematodes also lost the MYC proto-oncogenic regulator of cell proliferation. To identify additional regulators of cell proliferation that were lost with MYC, we computationally screened and determined 839 high-confidence genes that are conserved in bilaterians/lost in nematodes (CIBLIN genes). We find that 30 % of all CIBLIN genes encode transcriptional regulators of cell proliferation, epithelial-to-mesenchyme transitions, and other processes. Over 50 % of CIBLIN genes are unnamed genes in Drosophila, suggesting that there are many understudied genes. Interestingly, CIBLIN genes include many Myc synthetic lethal (MycSL) hits from recent screens. CIBLIN genes include key regulators of heparan sulfate proteoglycan (HSPG) sulfation patterns, and lysyl oxidases involved in cross-linking and modification of the extracellular matrix (ECM). These genes and others suggest the CIBLIN repertoire services critical functions in ECM remodeling and cell migration in large-bodied bilaterians. Correspondingly, CIBLIN genes are co-expressed with Myc in cancer transcriptomes, and include a preponderance of known determinants of cancer progression and tumor aggression. We propose that CIBLIN gene research can improve our understanding of regulatory control of cellular growth in metazoans.

An Atlas of Network Topologies Reveals Design Principles for Caenorhabditis elegans Vulval Precursor Cell Fate Patterning

The vulval precursor cell (VPC) fate patterning in Caenorhabditis elegans is a classic model experimental system for cell fate determination and patterning in development. Despite its apparent simplicity (six neighboring cells arranged in one dimension) and many experimental and computational efforts, the patterning strategy and mechanism remain controversial due to incomplete knowledge of the complex biology. Here, we carry out a comprehensive computational analysis and obtain a reservoir of all possible network topologies that are capable of VPC fate patterning under the simulation of various biological environments and regulatory rules. We identify three patterning strategies: sequential induction, morphogen gradient and lateral antagonism, depending on the features of the signal secreted from the anchor cell. The strategy of lateral antagonism, which has not been reported in previous studies of VPC patterning, employs a mutual inhibition of the 2° cell fate in neighboring cells. Robust topologies are built upon minimal topologies with basic patterning strategies and have more flexible and redundant implementations of modular functions. By simulated mutation, we find that all three strategies can reproduce experimental error patterns of mutants. We show that the topology derived by mapping currently known biochemical pathways to our model matches one of our identified functional topologies. Furthermore, our robustness analysis predicts a possible missing link related to the lateral antagonism strategy. Overall, we provide a theoretical atlas of all possible functional networks in varying environments, which may guide novel discoveries of the biological interactions in vulval development of Caenorhabditis elegans and related species.