Developmental function and state transitions of a gene expression oscillator in Caenorhabditis elegans

C. elegans epicuticlins define specific compartments in the apical extracellular matrix and function in wound repair

The apical extracellular matrix (aECM) of external epithelia often contains lipid-rich outer layers that contribute to permeability barrier function. The external aECM of nematodes is known as the cuticle and contains an external lipid-rich layer - the epicuticle. Epicuticlins are a family of tandem repeat cuticle proteins of unknown function. Here, we analyze the localization and function of the three C. elegans epicuticlins (EPIC proteins). EPIC-1 and EPIC-2 localize to the surface of the cuticle near the outer lipid layer, as well as to interfacial cuticles and adult-specific struts. EPIC-3 is expressed in dauer larvae and localizes to interfacial aECM in the buccal cavity. Skin wounding in the adult induces epic-3 expression, and EPIC proteins localize to wound sites. Null mutants lacking EPIC proteins are viable with reduced permeability barrier function and normal epicuticle lipid mobility. Loss of function in EPIC genes modifies the skin blistering phenotypes of Bli mutants and reduces survival after skin wounding. Our results suggest EPIC proteins define specific cortical compartments of the aECM and promote wound repair.

The NHR-23-regulated putative protease inhibitor mlt-11 gene is necessary for C. elegans cuticle structure and function

C. elegans molting offers a powerful entry point to understanding developmentally programmed apical extracellular matrix remodeling. However, the gene regulatory network controlling this process remains poorly understood. Focusing on targets of NHR-23, a key transcription factor that drives molting, we confirmed the Kunitz family protease inhibitor gene mlt-11 as an NHR-23 target. Through reporter assays, we identified NHR-23-binding sites that are necessary and sufficient for epithelial expression. We generated a translational fusion and demonstrated that MLT-11 is localized to the cuticle and lined openings to the exterior (vulva, rectum, mouth). We created a set of strains expressing varied levels of MLT-11 by deleting endogenous cis-regulatory element sequences. Combined deletion of two cis-regulatory elements caused developmental delay, motility defects, and failure of the cuticle barrier. Inactivation of mlt-11 by RNAi produced even more pronounced defects. mlt-11 is necessary to pattern every layer of the adult cuticle, suggesting a broad patterning role prior to the formation of the mature cuticle. Together these studies provide an entry point into understanding how individual cis-regulatory elements function to coordinate expression of oscillating genes involved in molting and how MLT-11 ensures proper cuticle assembly.

Caenorhabditis elegans Hedgehog-related proteins are tissue- and substructure-specific components of the cuticle and precuticle

In Caenorhabditis elegans, expanded families of divergent Hedgehog-related and patched-related proteins promote numerous processes ranging from epithelial and sense organ development to pathogen responses to cuticle shedding during the molt cycle. The molecular functions of these proteins have been mysterious since nematodes lack a canonical Hedgehog signaling pathway. Here we show that Hedgehog-related proteins are components of the cuticle and precuticle apical extracellular matrices that coat, shape, and protect external epithelia. Of four Hedgehog-related proteins imaged, two (GRL-2 and GRL-18) stably associated with the cuticles of specific tubes and two (GRL-7 and WRT-10) labeled precuticle substructures such as furrows or alae. We found that wrt-10 mutations disrupt cuticle alae ridges, consistent with a structural role in matrix organization. We hypothesize that most nematode Hedgehog-related proteins are apical extracellular matrix components, a model that could explain many of the reported functions for this family. These results highlight ancient connections between Hedgehog proteins and the extracellular matrix and suggest that any signaling roles of C. elegans Hedgehog-related proteins will be intimately related to their matrix association.

The Caenorhabditis elegans cuticle and precuticle: a model for studying dynamic apical extracellular matrices in vivo

Apical extracellular matrices (aECMs) coat the exposed surfaces of animal bodies to shape tissues, influence social interactions, and protect against pathogens and other environmental challenges. In the nematode Caenorhabditis elegans, collagenous cuticle and zona pellucida protein-rich precuticle aECMs alternately coat external epithelia across the molt cycle and play many important roles in the worm's development, behavior, and physiology. Both these types of aECMs contain many matrix proteins related to those in vertebrates, as well as some that are nematode-specific. Extensive differences observed among tissues and life stages demonstrate that aECMs are a major feature of epithelial cell identity. In addition to forming discrete layers, some cuticle components assemble into complex substructures such as ridges, furrows, and nanoscale pillars. The epidermis and cuticle are mechanically linked, allowing the epidermis to sense cuticle damage and induce protective innate immune and stress responses. The C. elegans model, with its optical transparency, facilitates the study of aECM cell biology and structure/function relationships and all the myriad ways by which aECM can influence an organism.

Ghost Channels and Ghost Cycles Guiding Long Transients in Dynamical Systems

Dynamical descriptions and modeling of natural systems have generally focused on fixed points, with saddles and saddle-based phase-space objects such as heteroclinic channels or cycles being central concepts behind the emergence of quasistable long transients. Reliable and robust transient dynamics observed for real, inherently noisy systems is, however, not met by saddle-based dynamics, as demonstrated here. Generalizing the notion of ghost states, we provide a complementary framework that does not rely on the precise knowledge or existence of (un)stable fixed points, but rather on slow directed flows organized by ghost sets in ghost channels and ghost cycles. Moreover, we show that the appearance of these novel objects is an emergent property of a broad class of models typically used for description of natural systems.

Essential function of transmembrane transcription factor MYRF in promoting transcription of miRNA lin-4 during C. elegans development

Precise developmental timing control is essential for organism formation and function, but its mechanisms are unclear. In C. elegans, the microRNA lin-4 critically regulates developmental timing by post-transcriptionally downregulating the larval-stage-fate controller LIN-14. However, the mechanisms triggering the activation of lin-4 expression toward the end of the first larval stage remain unknown. We demonstrate that the transmembrane transcription factor MYRF-1 is necessary for lin-4 activation. MYRF-1 is initially localized on the cell membrane, and its increased cleavage and nuclear accumulation coincide with lin-4 expression timing. MYRF-1 regulates lin-4 expression cell-autonomously and hyperactive MYRF-1 can prematurely drive lin-4 expression in embryos and young first-stage larvae. The tandem lin-4 promoter DNA recruits MYRF-1GFP to form visible loci in the nucleus, suggesting that MYRF-1 directly binds to the lin-4 promoter. Our findings identify a crucial link in understanding developmental timing regulation and establish MYRF-1 as a key regulator of lin-4 expression.

Accelerated hermaphrodite maturation on male pheromones suggests a general principle of coordination between larval behavior and development

Environment in general and social signals in particular could alter development. In Caenorhabditis elegans, male pheromones hasten development of hermaphrodite larvae. We show that this involves acceleration of growth and both somatic and germline development during the last larval stage (L4). Larvae exposed to male pheromones spend more time in L3 and less in the quiescent period between L3 and L4. This behavioral alteration improves provision in early L4, likely allowing for faster development. Larvae must be exposed to male pheromones in late L3 for behavioral and developmental effects to occur. Latter portions of other larval stages also contain periods of heightened sensitivity to environmental signals. Behavior during the early part of the larval stages is biased toward exploration, whereas later the emphasis shifts to food consumption. We argue that this organization allows assessment of the environment to identify the most suitable patch of resources, followed by acquisition of sufficient nutrition and salient information for the developmental events in the next larval stage. Evidence from other species indicates that such coordination of behavior and development may be a general feature of larval development.

Epidermal cell fusion promotes the transition from an embryonic to a larval transcriptome in C. elegans

Cell fusion is a fundamental process in the development of multicellular organisms, yet its impact on gene regulation, particularly during crucial developmental stages, remains poorly understood. The Caenorhabditis elegans epidermis comprises 8-10 syncytial cells, with the largest integrating 139 individual nuclei through cell-cell fusion governed by the fusogenic protein EFF-1. To explore the effects of cell fusion on developmental progression and associated gene expression changes, we conducted transcriptomic analyses of eff-1 fusion-deficient mutants. Our RNAseq findings showed widespread transcriptomic changes that were enriched for epidermal genes and key molecular pathways involved in epidermal function during larval development. Subsequent single-molecule fluorescence in situ hybridization validated the altered expression of mRNA transcripts, confirming quantifiable changes in gene expression in the absence of embryonic epidermal fusion. These results underscore the significance of cell-cell fusion in shaping transcriptional programs during development and raise questions regarding the precise identities and specialized functions of different subclasses of nuclei within developing syncytial cells and tissues.

Recent advances in understanding microRNA function and regulation in C. elegans

MicroRNAs (miRNAs) were first discovered in C. elegans as essential post-transcriptional regulators of gene expression. Since their initial discovery, miRNAs have been implicated in numerous areas of physiology and disease in all animals examined. In recent years, the C. elegans model continues to contribute important advances to all areas of miRNA research. Technological advances in tissue-specific miRNA profiling and genome editing have driven breakthroughs in understanding biological functions of miRNAs, mechanism of miRNA action, and regulation of miRNAs. In this review, we highlight these new C. elegans findings from the past five to seven years.

Hypoxia-inducible factor induces cysteine dioxygenase and promotes cysteine homeostasis in Caenorhabditis elegans

Dedicated genetic pathways regulate cysteine homeostasis. For example, high levels of cysteine activate cysteine dioxygenase, a key enzyme in cysteine catabolism in most animal and many fungal species. The mechanism by which cysteine dioxygenase is regulated is largely unknown. In an unbiased genetic screen for mutations that activate cysteine dioxygenase (cdo-1) in the nematode Caenorhabditis elegans, we isolated loss-of-function mutations in rhy-1 and egl-9, which encode proteins that negatively regulate the stability or activity of the oxygen-sensing hypoxia inducible transcription factor (hif-1). EGL-9 and HIF-1 are core members of the conserved eukaryotic hypoxia response. However, we demonstrate that the mechanism of HIF-1-mediated induction of cdo-1 is largely independent of EGL-9 prolyl hydroxylase activity and the von Hippel-Lindau E3 ubiquitin ligase, the classical hypoxia signaling pathway components. We demonstrate that C. elegans cdo-1 is transcriptionally activated by high levels of cysteine and hif-1. hif-1-dependent activation of cdo-1 occurs downstream of an H2S-sensing pathway that includes rhy-1, cysl-1, and egl-9. cdo-1 transcription is primarily activated in the hypodermis where it is also sufficient to drive sulfur amino acid metabolism. Thus, the regulation of cdo-1 by hif-1 reveals a negative feedback loop that maintains cysteine homeostasis. High levels of cysteine stimulate the production of an H2S signal. H2S then acts through the rhy-1/cysl-1/egl-9 signaling pathway to increase HIF-1-mediated transcription of cdo-1, promoting degradation of cysteine via CDO-1.

The effect of developmental variation on expression QTLs in a multi parental Caenorhabditis elegans population

Regulation of gene expression plays a crucial role in developmental processes and adaptation to changing environments. expression quantitative trait locus (eQTL) mapping is a technique used to study the genetic regulation of gene expression using the transcriptomes of recombinant inbred lines (RILs). Typically, the age of the inbred lines at the time of RNA sampling is carefully controlled. This is necessary because the developmental process causes changes in gene expression, complicating the interpretation of eQTL mapping experiments. However, due to genetics and variation in ambient micro-environments, organisms can differ in their "developmental age," even if they are of the same chronological age. As a result, eQTL patterns are affected by developmental variation in gene expression. The model organism Caenorhabditis elegans is particularly suited for studying the effect of developmental variation on eQTL mapping patterns. In a span of days, C. elegans transitions from embryo through 4 larval stages to adult while undergoing massive changes to its transcriptome. Here, we use C. elegans to investigate the effect of developmental age variation on eQTL patterns and present a normalization procedure. We used dynamical eQTL mapping, which includes the developmental age as a cofactor, to separate the variation in development from genotypic variation and explain variation in gene expression levels. We compare classical single marker eQTL mapping and dynamical eQTL mapping using RNA-seq data of ∼200 multi-parental RILs of C. elegans. The results show that (1) many eQTLs are caused by developmental variation, (2) most trans-bands are developmental QTLs, and (3) dynamical eQTL mapping detects additional eQTLs not found with classical eQTL mapping. We recommend that correction for variation in developmental age should be strongly considered in eQTL mapping studies given the large impact of processes like development on the transcriptome.

C. elegans epicuticlins define specific compartments in the apical extracellular matrix and function in wound repair

The apical extracellular matrix (aECM) of external epithelia often contains lipid-rich outer layers that contribute to permeability barrier function. The external aECM of nematode is known as the cuticle and contains an external lipid-rich layer, the epicuticle. Epicuticlins are a family of tandem repeat proteins originally identified as components of the insoluble fraction of the cuticular aECM and thought to localize in or near epicuticle. However, there has been little in vivo analysis of epicuticlins. Here, we report the localization analysis of the three C. elegans epicuticlins (EPIC proteins) using fluorescent protein knock-ins to visualize endogenously expressed proteins, and further examine their in vivo function using genetic null mutants. By TIRF microscopy, we find that EPIC-1 and EPIC-2 localize to the surface of the cuticle in larval and adult stages in close proximity to the outer lipid layer. EPIC-1 and EPIC-2 also localize to interfacial cuticles and adult-specific cuticle struts. EPIC-3 expression is restricted to the stress-induced dauer stage, where it localizes to interfacial aECM in the buccal cavity. Strikingly, skin wounding in the adult induces epic-3 expression, and EPIC-3::mNG localizes to wound scars. Null mutants lacking one, two, or all three EPIC proteins display reduced survival after skin wounding yet are viable with low penetrance defects in epidermal morphogenesis. Our results suggest EPIC proteins define specific aECM compartments and have roles in wound repair.

C. elegans Hedgehog-related proteins are tissue- and substructure-specific components of the cuticle and pre-cuticle

In C. elegans, divergent Hedgehog-related (Hh-r) and Patched-related (PTR) proteins promote numerous processes ranging from epithelial and sense organ development to pathogen responses to cuticle shedding during the molt cycle. Here we show that Hh-r proteins are actual components of the cuticle and pre-cuticle apical extracellular matrices (aECMs) that coat, shape, and protect external epithelia. Different Hh-r proteins stably associate with the aECMs of specific tissues and with specific substructures such as furrows and alae. Hh-r mutations can disrupt matrix structure. These results provide a unifying model for the function of nematode Hh-r proteins and highlight ancient connections between Hh proteins and the extracellular matrix.

The Clockwork Embryo: Mechanisms Regulating Developmental Rate

Organismal development requires the reproducible unfolding of an ordered sequence of discrete steps (cell fate determination, migration, tissue folding, etc.) in both time and space. Here, we review the mechanisms that grant temporal specificity to developmental steps, including molecular clocks and timers. Individual timing mechanisms must be coordinated with each other to maintain the overall developmental sequence. However, phenotypic novelties can also arise through the modification of temporal patterns over the course of evolution. Two main types of variation in temporal patterning characterize interspecies differences in developmental time: allochrony, where the overall developmental sequence is either accelerated or slowed down while maintaining the relative duration of individual steps, and heterochrony, where the duration of specific developmental steps is altered relative to the rest. New advances in in vitro modeling of mammalian development using stem cells have recently enabled the revival of mechanistic studies of allochrony and heterochrony. In both cases, differences in the rate of basic cellular functions such as splicing, translation, protein degradation, and metabolism seem to underlie differences in developmental time. In the coming years, these studies should identify the genetic differences that drive divergence in developmental time between species.

Canalisation and plasticity on the developmental manifold of Caenorhabditis elegans

How do the same mechanisms that faithfully regenerate complex developmental programmes in spite of environmental and genetic perturbations also allow responsiveness to environmental signals, adaptation and genetic evolution? Using the nematode Caenorhabditis elegans as a model, we explore the phenotypic space of growth and development in various genetic and environmental contexts. Our data are growth curves and developmental parameters obtained by automated microscopy. Using these, we show that among the traits that make up the developmental space, correlations within a particular context are predictive of correlations among different contexts. Furthermore, we find that the developmental variability of this animal can be captured on a relatively low-dimensional phenotypic manifold and that on this manifold, genetic and environmental contributions to plasticity can be deconvolved independently. Our perspective offers a new way of understanding the relationship between robustness and flexibility in complex systems, suggesting that projection and concentration of dimension can naturally align these forces as complementary rather than competing.

An nhr-85::GFP::AID*::3xFLAG knock-in allele for investigation of molting and oscillatory gene regulation

C. elegans NHR-85 is a poorly characterized nuclear hormone receptor transcription factor with an emerging role in regulating microRNA expression to control developmental timing. We generated the first NHR-85 translational fusion by knocking a GFP::AID*::3xFLAG cassette into the endogenous locus to tag all known isoforms. nhr-85 ::GFP::AID*::3xFLAG animals have wild-type broodsizes and NHR-85 ::GFP peaks in expression at the start of the L4 stage in epithelial cells. NHR-85 is not expressed in the germline, suggesting that while it might cooperate with the NHR-23 transcription factor to control microRNA expression, NHR-23 promotes spermatogenesis independent of NHR-85 . This nhr-85 ::GFP::AID*::3xFLAG strain will be a valuable resource for studying when and where NHR-85 acts to promote developmental timing.

The proprotein convertase BLI-4 promotes collagen secretion prior to assembly of the Caenorhabditis elegans cuticle

Some types of collagens, including transmembrane MACIT collagens and C. elegans cuticle collagens, are N-terminally cleaved at a dibasic site that resembles the consensus for furin or other proprotein convertases of the subtilisin/kexin (PCSK) family. Such cleavage may release transmembrane collagens from the plasma membrane and affect extracellular matrix assembly or structure. However, the functional consequences of such cleavage are unclear and evidence for the role of specific PCSKs is lacking. Here, we used endogenous collagen fusions to fluorescent proteins to visualize the secretion and assembly of the first collagen-based cuticle in C. elegans and then tested the role of the PCSK BLI-4 in these processes. Unexpectedly, we found that cuticle collagens SQT-3 and DPY-17 are secreted into the extraembryonic space several hours before cuticle matrix assembly. Furthermore, this early secretion depends on BLI-4/PCSK; in bli-4 and cleavage-site mutants, SQT-3 and DPY-17 are not efficiently secreted and instead form large intracellular puncta. Their later assembly into cuticle matrix is reduced but not entirely blocked. These data reveal a role for collagen N-terminal processing in intracellular trafficking and the control of matrix assembly in vivo. Our observations also prompt a revision of the classic model for C. elegans cuticle matrix assembly and the pre-cuticle-to-cuticle transition, suggesting that cuticle layer assembly proceeds via a series of regulated steps and not simply by sequential secretion and deposition.

Periods of environmental sensitivity couple larval behavior and development

The typical life cycle in most animal phyla includes a larval period that bridges embryogenesis and adulthood1. Despite the great diversity of larval forms, all larvae grow, acquire adult morphology and function, while navigating their habitats to obtain resources necessary for development. How larval development is coordinated with behavior remains substantially unclear. Here, we describe features of the iterative organization of larval stages that serve to assess the environment and procure resources prior to costly developmental commitments. We found that male-excreted pheromones accelerate2-4 the onset of adulthood in C. elegans hermaphrodites by coordinately advancing multiple developmental events and growth during the last larval stage. The larvae are sensitive to the accelerating male pheromones only at the end of the penultimate larval stage, just before the acceleration begins. Other larval stages also contain windows of sensitivity to environmental inputs. Importantly, behaviors associated with search and consumption of food are distinct between early and late portions of larval stages. We infer that each larval stage in C. elegans is subdivided into two epochs: A) global assessment of the environment to identify the most suitable patch and B) consumption of sufficient food and acquisition of salient information for developmental events in the next stage. We predict that in larvae of other species behavior is also divided into distinct epochs optimized either for assessing the habitat or obtaining the resources. Thus, a major role of larval behavior is to coordinate the orderly progression of development in variable environments.

Effectors of anterior morphogenesis in C. elegans embryos

During embryogenesis the nascent Caenorhabditis elegans epidermis secretes an apical extracellular matrix (aECM) that serves as an external stabilizer, preventing deformation of the epidermis by mechanical forces exerted during morphogenesis. At present, the factors that contribute to aECM function are mostly unknown, including the aECM components themselves, their posttranslational regulators, and the pathways required for their secretion. Here we showed that two proteins previously linked to aECM function, SYM-3/FAM102A and SYM-4/WDR44, colocalize to intracellular and membrane-associated puncta and likely function in a complex. Proteomics experiments also suggested potential roles for SYM-3/FAM102A and SYM-4/WDR44 family proteins in intracellular trafficking. Nonetheless, we found no evidence to support a critical function for SYM-3 or SYM-4 in the apical deposition of two aECM components, NOAH-1 and FBN-1. Moreover, loss of a key splicing regulator of fbn-1, MEC-8/RBPMS2, had surprisingly little effect on the abundance or deposition of FBN-1. Using a focused screening approach, we identified 32 additional proteins that likely contribute to the structure and function of the embryonic aECM. We also characterized morphogenesis defects in embryos lacking mir-51 microRNA family members, which display a similar phenotype to mec-8; sym double mutants. Collectively, these findings add to our knowledge of factors controlling embryonic morphogenesis.

The proprotein convertase BLI-4 promotes collagen secretion during assembly of the Caenorhabditis elegans cuticle

Some types of collagens, including transmembrane MACIT collagens and C. elegans cuticle collagens, are N-terminally cleaved at a dibasic site that resembles the consensus for furin or other proprotein convertases of the subtilisin/kexin (PCSK) family. Such cleavage may release transmembrane collagens from the plasma membrane and affect extracellular matrix assembly or structure. However, the functional consequences of such cleavage are unclear and evidence for the role of specific PCSKs is lacking. Here, we used endogenous collagen fusions to fluorescent proteins to visualize the secretion and assembly of the first collagen-based cuticle in C. elegans and then tested the role of the PCSK BLI-4 in these processes. Unexpectedly, we found that cuticle collagens SQT-3 and DPY-17 are secreted into the extraembryonic space several hours before cuticle matrix assembly. Furthermore, this early secretion depends on BLI-4/PCSK; in bli-4 and cleavage-site mutants, SQT-3 and DPY-17 are not efficiently secreted and instead form large intracellular aggregates. Their later assembly into cuticle matrix is reduced but not entirely blocked. These data reveal a role for collagen N-terminal processing in intracellular trafficking and in the spatial and temporal restriction of matrix assembly in vivo . Our observations also prompt a revision of the classic model for C. elegans cuticle matrix assembly and the pre-cuticle-to-cuticle transition, suggesting that cuticle layer assembly proceeds via a series of regulated steps and not simply by sequential secretion and deposition.

NHR-23 activity is necessary for C. elegans developmental progression and apical extracellular matrix structure and function

Nematode molting is a remarkable process where animals must repeatedly build a new apical extracellular matrix (aECM) beneath a previously built aECM that is subsequently shed. The nuclear hormone receptor NHR-23 (also known as NR1F1) is an important regulator of C. elegans molting. NHR-23 expression oscillates in the epidermal epithelium, and soma-specific NHR-23 depletion causes severe developmental delay and death. Tissue-specific RNAi suggests that nhr-23 acts primarily in seam and hypodermal cells. NHR-23 coordinates the expression of factors involved in molting, lipid transport/metabolism and remodeling of the aECM. NHR-23 depletion causes dampened expression of a nas-37 promoter reporter and a loss of reporter oscillation. The cuticle collagen ROL-6 and zona pellucida protein NOAH-1 display aberrant annular localization and severe disorganization over the seam cells after NHR-23 depletion, while the expression of the adult-specific cuticle collagen BLI-1 is diminished and frequently found in patches. Consistent with these localization defects, the cuticle barrier is severely compromised when NHR-23 is depleted. Together, this work provides insight into how NHR-23 acts in the seam and hypodermal cells to coordinate aECM regeneration during development.

Pathways that affect anterior morphogenesis in C. elegans embryos

During embryogenesis the nascent Caenorhabditis elegans epidermis secretes an apical extracellular matrix (aECM) that serves as an external stabilizer, preventing deformation of the epidermis by mechanical forces exerted during morphogenesis. We showed that two conserved proteins linked to this process, SYM-3/FAM102A and SYM-4/WDR44, colocalize to intracellular and membrane-associated puncta and likely function together in a complex. Proteomics data also suggested potential roles for FAM102A and WDR44 family proteins in intracellular trafficking, consistent with their localization patterns. Nonetheless, we found no evidence to support a clear function for SYM-3 or SYM-4 in the apical deposition of two aECM components, FBN-1 and NOAH. Surprisingly, loss of MEC-8/RBPMS2, a conserved splicing factor and regulator of fbn-1 , had little effect on the abundance or deposition of FBN-1 to the aECM. Using a focused screening approach, we identified 32 additional proteins that likely contribute to the structure and function of the embryonic aECM. Lastly, we examined morphogenesis defects in embryos lacking mir-51 microRNA family members, which display a related embryonic phenotype to mec-8; sym double mutants. Collectively, our findings add to our knowledge of pathways controlling embryonic morphogenesis.

SUMMARY STATEMENT

We identify new proteins in apical ECM biology in C. elegans and provide evidence that SYM-3/FAM102A and SYM-4/WDR44 function together in trafficking but do not regulate apical ECM protein deposition.

Experimental considerations for study of C. elegans lysosomal proteins

Lysosomes are an important organelle required for the degradation of a range of cellular components. Lysosome function is critical for development and homeostasis as dysfunction can lead to inherited genetic disorders, cancer, and neurodegenerative and metabolic diseases. The acidic and protease-rich environment of lysosomes poses experimental challenges. Many fluorescent proteins are quenched or degraded, while specific red fluorescent proteins can be cleaved from translational fusion partners and accumulate. While studying MLT-11, a Caenorhabditis elegans molting factor that localizes to lysosomes and the cuticle, we sought to optimize several experimental parameters. We found that, in contrast to mNeonGreen fusions, mScarlet fusions to MLT-11 missed cuticular and rectal epithelial localization. Rapid sample lysis and denaturation were critical for preventing MLT-11 fragmentation while preparing lysates for western blots. Using a model lysosomal substrate (NUC-1), we found that rigid polyproline linkers and truncated mCherry constructs do not prevent cleavage of mCherry from NUC-1. We provide evidence that extended localization in lysosomal environments prevents the detection of FLAG epitopes in western blots. Finally, we optimize an acid-tolerant green fluorescent protein (Gamillus) for use in C. elegans. These experiments provide important experimental considerations and new reagents for the study of C. elegans lysosomal proteins.

C. elegans molting requires rhythmic accumulation of the Grainyhead/LSF transcription factor GRH-1

C. elegans develops through four larval stages that are rhythmically terminated by molts, that is, the synthesis and shedding of a cuticular exoskeleton. Each larval cycle involves rhythmic accumulation of thousands of transcripts, which we show here relies on rhythmic transcription. To uncover the responsible gene regulatory networks (GRNs), we screened for transcription factors that promote progression through the larval stages and identified GRH-1, BLMP-1, NHR-23, NHR-25, MYRF-1, and BED-3. We further characterize GRH-1, a Grainyhead/LSF transcription factor, whose orthologues in other animals are key epithelial cell-fate regulators. We find that GRH-1 depletion extends molt durations, impairs cuticle integrity and shedding, and causes larval death. GRH-1 is required for, and accumulates prior to, each molt, and preferentially binds to the promoters of genes expressed during this time window. Binding to the promoters of additional genes identified in our screen furthermore suggests that we have identified components of a core molting-clock GRN. Since the mammalian orthologues of GRH-1, BLMP-1 and NHR-23, have been implicated in rhythmic homeostatic skin regeneration in mouse, the mechanisms underlying rhythmic C. elegans molting may apply beyond nematodes.

Loss of the E3 ubiquitin ligases UBR-5 or HECD-1 restores Caenorhabditis elegans development in the absence of SWI/SNF function

SWItch/sucrose non-fermenting (SWI/SNF) complexes are a family of chromatin remodelers that are conserved across eukaryotes. Mutations in subunits of SWI/SNF cause a multitude of different developmental disorders in humans, most of which have no current treatment options. Here, we identify an alanine-to-valine-causing mutation in the SWI/SNF subunit snfc-5 (SMARCB1 in humans) that prevents embryonic lethality in Caenorhabditis elegans nematodes harboring a loss-of-function mutation in the SWI/SNF subunit swsn-1 (SMARCC1/2 in humans). Furthermore, we found that the combination of this specific mutation in snfc-5 and a loss-of-function mutation in either of the E3 ubiquitin ligases ubr-5 (UBR5 in humans) or hecd-1 (HECTD1 in humans) can restore development to adulthood in swsn-1 loss-of-function mutants that otherwise die as embryos. Using these mutant models, we established a set of 335 genes that are dysregulated in SWI/SNF mutants that arrest their development embryonically but exhibit near wild-type levels of expression in the presence of suppressor mutations that prevent embryonic lethality, suggesting that SWI/SNF promotes development by regulating some subset of these 335 genes. In addition, we show that SWI/SNF protein levels are reduced in swsn-1; snfc-5 double mutants and partly restored to wild-type levels in swsn-1; snfc-5; ubr-5 triple mutants, consistent with a model in which UBR-5 regulates SWI/SNF levels by tagging the complex for proteasomal degradation. Our findings establish a link between two E3 ubiquitin ligases and SWI/SNF function and suggest that UBR5 and HECTD1 could be potential therapeutic targets for the many developmental disorders caused by missense mutations in SWI/SNF subunits.

A specific type of Argonaute phosphorylation regulates binding to microRNAs during C. elegans development

Argonaute proteins are at the core of the microRNA-mediated gene silencing pathway essential for animals. In C. elegans, the microRNA-specific Argonautes ALG-1 and ALG-2 regulate multiple processes required for proper animal developmental timing and viability. Here we identified a phosphorylation site on ALG-1 that modulates microRNA association. Mutating ALG-1 serine 642 into a phospho-mimicking residue impairs microRNA binding and causes embryonic lethality and post-embryonic phenotypes that are consistent with alteration of microRNA functions. Monitoring microRNA levels in alg-1 phosphorylation mutant animals shows that microRNA passenger strands increase in abundance but are not preferentially loaded into ALG-1, indicating that the miRNA binding defects could lead to microRNA duplex accumulation. Our genetic and biochemical experiments support protein kinase A (PKA) KIN-1 as the putative kinase that phosphorylates ALG-1 serine 642. Our data indicate that PKA triggers ALG-1 phosphorylation to regulate its microRNA association during C. elegans development.

A spatiotemporal reconstruction of the C. elegans pharyngeal cuticle reveals a structure rich in phase-separating proteins

How the cuticles of the roughly 4.5 million species of ecdysozoan animals are constructed is not well understood. Here, we systematically mine gene expression datasets to uncover the spatiotemporal blueprint for how the chitin-based pharyngeal cuticle of the nematode Caenorhabditis elegans is built. We demonstrate that the blueprint correctly predicts expression patterns and functional relevance to cuticle development. We find that as larvae prepare to molt, catabolic enzymes are upregulated and the genes that encode chitin synthase, chitin cross-linkers, and homologs of amyloid regulators subsequently peak in expression. Forty-eight percent of the gene products secreted during the molt are predicted to be intrinsically disordered proteins (IDPs), many of which belong to four distinct families whose transcripts are expressed in overlapping waves. These include the IDPAs, IDPBs, and IDPCs, which are introduced for the first time here. All four families have sequence properties that drive phase separation and we demonstrate phase separation for one exemplar in vitro. This systematic analysis represents the first blueprint for cuticle construction and highlights the massive contribution that phase-separating materials make to the structure.

Arnold tongue entrainment reveals dynamical principles of the embryonic segmentation clock

Living systems exhibit an unmatched complexity, due to countless, entangled interactions across scales. Here, we aim to understand a complex system, that is, segmentation timing in mouse embryos, without a reference to these detailed interactions. To this end, we develop a coarse-grained approach, in which theory guides the experimental identification of the segmentation clock entrainment responses. We demonstrate period- and phase-locking of the segmentation clock across a wide range of entrainment parameters, including higher-order coupling. These quantifications allow to derive the phase response curve (PRC) and Arnold tongues of the segmentation clock, revealing its essential dynamical properties. Our results indicate that the somite segmentation clock has characteristics reminiscent of a highly non-linear oscillator close to an infinite period bifurcation and suggests the presence of long-term feedbacks. Combined, this coarse-grained theoretical-experimental approach reveals how we can derive simple, essential features of a highly complex dynamical system, providing precise experimental control over the pace and rhythm of the somite segmentation clock.

Feedback between a retinoid-related nuclear receptor and the let-7 microRNAs controls the pace and number of molting cycles in C. elegans

Animal development requires coordination among cyclic processes, sequential cell fate specifications, and once-a-lifetime morphogenic events, but the underlying timing mechanisms are not well understood. Caenorhabditis elegans undergoes four molts at regular 8 to 10 hour intervals. The pace of the cycle is governed by PERIOD/lin-42 and other as-yet unknown factors. Cessation of the cycle in young adults is controlled by the let-7 family of microRNAs and downstream transcription factors in the heterochronic pathway. Here, we characterize a negative feedback loop between NHR-23, the worm homolog of mammalian retinoid-related orphan receptors (RORs), and the let-7 family of microRNAs that regulates both the frequency and finite number of molts. The molting cycle is decelerated in nhr-23 knockdowns and accelerated in let-7(-) mutants, but timed similarly in let-7(-) nhr-23(-) double mutants and wild-type animals. NHR-23 binds response elements (ROREs) in the let-7 promoter and activates transcription. In turn, let-7 dampens nhr-23 expression across development via a complementary let-7-binding site (LCS) in the nhr-23 3' UTR. The molecular interactions between NHR-23 and let-7 hold true for other let-7 family microRNAs. Either derepression of nhr-23 transcripts by LCS deletion or high gene dosage of nhr-23 leads to protracted behavioral quiescence and extra molts in adults. NHR-23 and let-7 also coregulate scores of genes required for execution of the molts, including lin-42. In addition, ROREs and LCSs isolated from mammalian ROR and let-7 genes function in C. elegans, suggesting conservation of this feedback mechanism. We propose that this feedback loop unites the molting timer and the heterochronic gene regulatory network, possibly by functioning as a cycle counter.

Stereotyped behavioral maturation and rhythmic quiescence in C.elegans embryos

Systematic analysis of rich behavioral recordings is being used to uncover how circuits encode complex behaviors. Here, we apply this approach to embryos. What are the first embryonic behaviors and how do they evolve as early neurodevelopment ensues? To address these questions, we present a systematic description of behavioral maturation for Caenorhabditis elegans embryos. Posture libraries were built using a genetically encoded motion capture suit imaged with light-sheet microscopy and annotated using custom tracking software. Analysis of cell trajectories, postures, and behavioral motifs revealed a stereotyped developmental progression. Early movement is dominated by flipping between dorsal and ventral coiling, which gradually slows into a period of reduced motility. Late-stage embryos exhibit sinusoidal waves of dorsoventral bends, prolonged bouts of directed motion, and a rhythmic pattern of pausing, which we designate slow wave twitch (SWT). Synaptic transmission is required for late-stage motion but not for early flipping nor the intervening inactive phase. A high-throughput behavioral assay and calcium imaging revealed that SWT is elicited by the rhythmic activity of a quiescence-promoting neuron (RIS). Similar periodic quiescent states are seen prenatally in diverse animals and may play an important role in promoting normal developmental outcomes.

Real age prediction from the transcriptome with RAPToR

Transcriptomic data is often affected by uncontrolled variation among samples that can obscure and confound the effects of interest. This variation is frequently due to unintended differences in developmental stages between samples. The transcriptome itself can be used to estimate developmental progression, but existing methods require many samples and do not estimate a specimen's real age. Here we present real-age prediction from transcriptome staging on reference (RAPToR), a computational method that precisely estimates the real age of a sample from its transcriptome, exploiting existing time-series data as reference. RAPToR works with whole animal, dissected tissue and single-cell data for the most common animal models, humans and even for non-model organisms lacking reference data. We show that RAPToR can be used to remove age as a confounding factor and allow recovery of a signal of interest in differential expression analysis. RAPToR will be especially useful in large-scale single-organism profiling because it eliminates the need for accurate staging or synchronisation before profiling.

An unexpected role for the conserved ADAM-family metalloprotease ADM-2 in Caenorhabditis elegans molting

Molting is a widespread developmental process in which the external extracellular matrix (ECM), the cuticle, is remodeled to allow for organismal growth and environmental adaptation. Studies in the nematode Caenorhabditis elegans have identified a diverse set of molting-associated factors including signaling molecules, intracellular trafficking regulators, ECM components, and ECM-modifying enzymes such as matrix metalloproteases. C. elegans NEKL-2 and NEKL-3, two conserved members of the NEK family of protein kinases, are essential for molting and promote the endocytosis of environmental steroid-hormone precursors by the epidermis. Steroids in turn drive the cyclic induction of many genes required for molting. Here we report a role for the sole C. elegans ADAM–meltrin metalloprotease family member, ADM-2, as a mediator of molting. Loss of adm-2, including mutations that disrupt the metalloprotease domain, led to the strong suppression of molting defects in partial loss-of-function nekl mutants. ADM-2 is expressed in the epidermis, and its trafficking through the endo-lysosomal network was disrupted after NEKL depletion. We identified the epidermally expressed low-density lipoprotein receptor–related protein, LRP-1, as a candidate target of ADM-2 regulation. Whereas loss of ADM-2 activity led to the upregulation of apical epidermal LRP-1, ADM-2 overexpression caused a reduction in LRP-1 levels. Consistent with this, several mammalian ADAMs, including the meltrin ADAM12, have been shown to regulate mammalian LRP1 via proteolysis. In contrast to mammalian homologs, however, the regulation of LRP-1 by ADM-2 does not appear to involve the metalloprotease function of ADM-2, nor is proteolytic processing of LRP-1 strongly affected in adm-2 mutants. Our findings suggest a noncanonical role for an ADAM family member in the regulation of a lipoprotein-like receptor and lead us to propose that endocytic trafficking may be important for both the internalization of factors that promote molting as well as the removal of proteins that can inhibit the process.

The molecular and cellular features of molting in nematodes share many similarities with cellular and developmental processes that occur in mammals. This includes the degradation and reorganization of extracellular matrix materials that surround cells, as well as the intracellular machineries that allow cells to sample and modify their environments. In the current study, we found an unexpected function for a conserved protein that cleaves other proteins on the external surface of cells. Rather than promoting molting through extracellular matrix reorganization, however, the ADM-2 protease appears to function as a negative regulator of molting. This observation can be explained in part by data showing that ADM-2 negatively regulates a cell surface receptor required for molting. Surprisingly, it appears to do so through a mechanism that does not involve proteolysis. Our data provide insights into the mechanisms controlling molting and link several conserved proteins to show how they function together during development.

Deviations from temporal scaling support a stage-specific regulation for C. elegans postembryonic development

BACKGROUND

After embryonic development, Caenorhabditis elegans progress through for larval stages, each of them finishing with molting. The repetitive nature of C. elegans postembryonic development is considered an oscillatory process, a concept that has gained traction from regulation by a circadian clock gene homologue. Nevertheless, each larval stage has a defined duration and entails specific events. Since the overall duration of development is controlled by numerous factors, we have asked whether different rate-limiting interventions impact all stages equally.

RESULTS

We have measured the duration of each stage of development for over 2500 larvae, under varied environmental conditions known to alter overall developmental rate. We applied changes in temperature and in the quantity and quality of nutrition and analysed the effect of genetically reduced insulin signalling. Our results show that the distinct developmental stages respond differently to these perturbations. The changes in the duration of specific larval stages seem to depend on stage-specific events. Furthermore, our high-resolution measurement of the effect of temperature on the stage-specific duration of development has unveiled novel features of temperature dependence in C. elegans postembryonic development.

CONCLUSIONS

Altogether, our results show that multiple factors fine tune developmental timing, impacting larval stages independently. Further understanding of the regulation of this process will allow modelling the mechanisms that control developmental timing.

Temporal scaling in C. elegans larval development

SignificanceAn enduring mystery of development is how its timing is controlled, particularly for development after birth, where timing is highly flexible and depends on environmental conditions, such as food availability and diet. We followed timing of cell- and organism-level events in individual Caenorhabditis elegans larvae developing from hatching to adulthood, uncovering widespread variations in event timing, both between isogenic individuals in the same environment and when changing conditions and genotypes. However, in almost all cases, we found that events occurred at the same time, when time was rescaled by the duration of development measured in each individual. This observation of "temporal scaling" poses strong constraints on models to explain timing of larval development.

A life cycle alteration can correct molting defects in Caenorhabditis elegans

Molting is a widespread feature in the development of many invertebrates, including nematodes and arthropods. In Caenorhabditis elegans, the highly conserved protein kinases NEKL-2/NEK8/9 and NEKL-3/NEK6/7 (NEKLs) promote molting through their involvement in the uptake and intracellular trafficking of epidermal cargos. We found that the relative requirements for NEKL-2 and NEKL-3 differed at different life-cycle stages and under different environmental conditions. Most notably, the transition from the second to the third larval stage (L2→L3 molt) required a higher level of NEKL function than during several other life stages or when animals had experienced starvation at the L1 stage. Specifically, larvae that entered the pre-dauer L2d stage could escape molting defects when transiting to the (non-dauer) L3 stage. Consistent with this, mutations that promote entry into L2d suppressed nekl-associated molting defects, whereas mutations that inhibit L2d entry reduced starvation-mediated suppression. We further showed that loss or reduction of NEKL functions led to defects in the transcription of cyclically expressed molting genes, many of which are under the control of systemic steroid hormone regulation. Moreover, the timing and severity of these transcriptional defects correlated closely with the strength of nekl alleles and with their stage of arrest. Interestingly, transit through L2d rescued nekl-associated expression defects in suppressed worms, providing an example of how life-cycle decisions can impact subsequent developmental events. Given that NEKLs are implicated in the uptake of sterols by the epidermis, we propose that loss of NEKLs leads to a physiological reduction in steroid-hormone signaling and consequent defects in the transcription of genes required for molting.

The oscillating Mucin-type protein DPY-6 has a conserved role in nematode mouth and cuticle formation

Nematodes show an extraordinary diversity of mouth structures and strikingly different feeding strategies, which has enabled an invasion of all ecosystems. However, nearly nothing is known about the structural and molecular architecture of the nematode mouth (stoma). Pristionchus pacificus is an intensively studied nematode that exhibits unique life history traits, including predation, teeth-like denticle formation, and mouth-form plasticity. Here, we used a large-scale genetic screen to identify genes involved in mouth formation. We identified Ppa-dpy-6 to encode a Mucin-type hydrogel-forming protein that is macroscopically involved in the specification of the cheilostom, the anterior part of the mouth. We used a recently developed protocol for geometric morphometrics of miniature animals to characterize these defects further and found additional defects that affect mouth form, shape, and size resulting in an overall malformation of the mouth. Additionally, Ppa-dpy-6 is shorter than wild-type with a typical Dumpy phenotype, indicating a role in the formation of the external cuticle. This concomitant phenotype of the cheilostom and cuticle provides the first molecular support for the continuity of these structures and for the separation of the cheilostom from the rest of the stoma. In Caenorhabditis elegans, dpy-6 was an early mapping mutant but its molecular identity was only determined during genome-wide RNAi screens and not further investigated. Strikingly, geometric morphometric analysis revealed previously unrecognized cheilostom and gymnostom defects in Cel-dpy-6 mutants. Thus, the Mucin-type protein DPY-6 represents to the best of our knowledge, the first protein involved in nematode mouth formation with a conserved role in cuticle deposition. This study opens new research avenues to characterize the molecular composition of the nematode mouth, which is associated with extreme ecological diversification.

Towards a physical understanding of developmental patterning

The temporal coordination of events at cellular and tissue scales is essential for the proper development of organisms, and involves cell-intrinsic processes that can be coupled by local cellular signalling and instructed by global signalling, thereby creating spatial patterns of cellular states that change over time. The timing and structure of these patterns determine how an organism develops. Traditional developmental genetic methods have revealed the complex molecular circuits regulating these processes but are limited in their ability to predict and understand the emergent spatio-temporal dynamics. Increasingly, approaches from physics are now being used to help capture the dynamics of the system by providing simplified, generic descriptions. Combined with advances in imaging and computational power, such approaches aim to provide insight into timing and patterning in developing systems.

Single worm transcriptomics identifies a developmental core network of oscillating genes with deep conservation across nematodes

High-resolution spatial and temporal maps of gene expression have facilitated a comprehensive understanding of animal development and evolution. In nematodes, the small body size represented a major challenge for such studies, but recent advancements have helped overcome this limitation. Here, we have implemented single worm transcriptomics (SWT) in the nematode model organism Pristionchus pacificus to provide a high-resolution map of the developmental transcriptome. We selected 38 time points from hatching of the J2 larvae to young adults to perform transcriptome analysis over 60 h of postembryonic development. A mean sequencing depth of 4.5 million read pairs allowed the detection of more than 23,135 (80%) of all genes. Nearly 3000 (10%) genes showed oscillatory expression with discrete expression levels, phases, and amplitudes. Gene age analysis revealed an overrepresentation of ancient gene classes among oscillating genes, and around one-third of them have 1:1 orthologs in C. elegans. One important gene family overrepresented among oscillating genes is collagens. Several of these collagen genes are regulated by the developmental switch gene eud-1, indicating a potential function in the regulation of mouth-form plasticity, a key developmental process in this facultative predatory nematode. Together, our analysis provides (1) an updated protocol for SWT in nematodes that is applicable to many microscopic species, (2) a 1- to 2-h high-resolution catalog of P. pacificus gene expression throughout postembryonic development, and (3) a comparative analysis of oscillatory gene expression between the two model organisms P. pacificus and C. elegans and associated evolutionary dynamics.

An Epigenetic Priming Mechanism Mediated by Nutrient Sensing Regulates Transcriptional Output during C. elegans Development

Although precise tuning of gene expression levels is critical for most developmental pathways, the mechanisms by which the transcriptional output of dosage-sensitive molecules is established or modulated by the environment remain poorly understood. Here, we provide a mechanistic framework for how the conserved transcription factor BLMP-1/Blimp1 operates as a pioneer factor to decompact chromatin near its target loci during embryogenesis (hours prior to major transcriptional activation) and, by doing so, regulates both the duration and amplitude of subsequent target gene transcription during post-embryonic development. This priming mechanism is genetically separable from the mechanisms that establish the timing of transcriptional induction and functions to canalize aspects of cell-fate specification, animal size regulation, and molting. A key feature of the BLMP-1-dependent transcriptional priming mechanism is that chromatin decompaction is initially established during embryogenesis and maintained throughout larval development by nutrient sensing. This anticipatory mechanism integrates transcriptional output with environmental conditions and is essential for resuming normal temporal patterning after animals exit nutrient-mediated developmental arrests.

Gene expression oscillations in C. elegans underlie a new developmental clock

During C. elegans larval development, thousands of genes, accounting for >20% of the transcriptome, exhibit oscillatory expression with large amplitudes. The time of peaking varies for different genes, but expression generally peaks once per larval stage, with both the oscillation period and larval stage duration varying in concert with temperature. This and other evidence support the existence of a gene expression oscillator that functions as a developmental clock. In this article, we review what is known about the biology, architecture and possible mechanisms of this clock. We compare it to other oscillators, and highlight tools and approaches suited to its study. Finally, we point out implications of these wide-spread and dynamic changes of gene expression on any type of gene expression profiling experiment in C. elegans larvae and how such experiments need to be controlled.

Tissue-Specific Transcription Footprinting Using RNA PoI DamID (RAPID) in Caenorhabditis elegans

Differential gene expression across cell types underlies development and cell physiology in multicellular organisms. Caenorhabditis elegans is a powerful, extensively used model to address these biological questions. A remaining bottleneck relates to the difficulty to obtain comprehensive tissue-specific gene transcription data, since available methods are still challenging to execute and/or require large worm populations. Here, we introduce the RNA Polymerase DamID (RAPID) approach, in which the Dam methyltransferase is fused to a ubiquitous RNA polymerase subunit to create transcriptional footprints via methyl marks on the DNA of transcribed genes. To validate the method, we determined the polymerase footprints in whole animals, in sorted embryonic blastomeres and in different tissues from intact young adults by driving tissue-specific Dam fusion expression. We obtained meaningful transcriptional footprints in line with RNA-sequencing (RNA-seq) studies in whole animals or specific tissues. To challenge the sensitivity of RAPID and demonstrate its utility to determine novel tissue-specific transcriptional profiles, we determined the transcriptional footprints of the pair of XXX neuroendocrine cells, representing 0.2% of the somatic cell content of the animals. We identified 3901 candidate genes with putatively active transcription in XXX cells, including the few previously known markers for these cells. Using transcriptional reporters for a subset of new hits, we confirmed that the majority of them were expressed in XXX cells and identified novel XXX-specific markers. Taken together, our work establishes RAPID as a valid method for the determination of RNA polymerase footprints in specific tissues of C. elegans without the need for cell sorting or RNA tagging.