Multiple regulatory elements with spatially and temporally distinct activities control the expression of the epithelial differentiation gene lin-26 in C. elegans

Mechanisms of lineage specification in Caenorhabditis elegans

The studies of cell fate and lineage specification are fundamental to our understanding of the development of multicellular organisms. Caenorhabditis elegans has been one of the premiere systems for studying cell fate specification mechanisms at single cell resolution, due to its transparent nature, the invariant cell lineage, and fixed number of somatic cells. We discuss the general themes and regulatory mechanisms that have emerged from these studies, with a focus on somatic lineages and cell fates. We next review the key factors and pathways that regulate the specification of discrete cells and lineages during embryogenesis and postembryonic development; we focus on transcription factors and include numerous lineage diagrams that depict the expression of key factors that specify embryonic founder cells and postembryonic blast cells, and the diverse somatic cell fates they generate. We end by discussing some future perspectives in cell and lineage specification.

Fluctuations in chromatin state at regulatory loci occur spontaneously under relaxed selection and are associated with epigenetically inherited variation in C. elegans gene expression

Some epigenetic information can be transmitted between generations without changes in the underlying DNA sequence. Changes in epigenetic regulators, termed epimutations, can occur spontaneously and be propagated in populations in a manner reminiscent of DNA mutations. Small RNA-based epimutations occur in C. elegans and persist for around 3-5 generations on average. Here, we explored whether chromatin states also undergo spontaneous change and whether this could be a potential alternative mechanism for transgenerational inheritance of gene expression changes. We compared the chromatin and gene expression profiles at matched time points from three independent lineages of C. elegans propagated at minimal population size. Spontaneous changes in chromatin occurred in around 1% of regulatory regions each generation. Some were heritable epimutations and were significantly enriched for heritable changes in expression of nearby protein-coding genes. Most chromatin-based epimutations were short-lived but a subset had longer duration. Genes subject to long-lived epimutations were enriched for multiple components of xenobiotic response pathways. This points to a possible role for epimutations in adaptation to environmental stressors.

A natural transdifferentiation event involving mitosis is empowered by integrating signaling inputs with conserved plasticity factors

Transdifferentiation, or direct cell reprogramming, is the conversion of one fully differentiated cell type into another. Whether core mechanisms are shared between natural transdifferentiation events when occurring with or without cell division is unclear. We have previously characterized the Y-to-PDA natural transdifferentiation in Caenorhabditis elegans, which occurs without cell division and requires orthologs of vertebrate reprogramming factors. Here, we identify a rectal-to-GABAergic transdifferentiation and show that cell division is required but not sufficient for conversion. We find shared mechanisms, including erasure of the initial identity, which requires the conserved reprogramming factors SEM-4/SALL, SOX-2, CEH-6/OCT, and EGL-5/HOX. We also find three additional and parallel roles of the Wnt signaling pathway: selection of a specific daughter, removal of the initial identity, and imposition of the precise final subtype identity. Our results support a model in which levels and antagonistic activities of SOX-2 and Wnt signaling provide a timer for the acquisition of final identity.

pop-1/TCF, ref-2/ZIC and T-box factors regulate the development of anterior cells in the C. elegans embryo

Patterning of the anterior-posterior axis is fundamental to animal development. The Wnt pathway plays a major role in this process by activating the expression of posterior genes in animals from worms to humans. This observation raises the question of whether the Wnt pathway or other regulators control the expression of the many anterior-expressed genes. We found that the expression of five anterior-specific genes in Caenorhabditis elegans embryos depends on the Wnt pathway effectors pop-1/TCF and sys-1/β-catenin. We focused further on one of these anterior genes, ref-2/ZIC, a conserved transcription factor expressed in multiple anterior lineages. Live imaging of ref-2 mutant embryos identified defects in cell division timing and position in anterior lineages. Cis-regulatory dissection identified three ref-2 transcriptional enhancers, one of which is necessary and sufficient for anterior-specific expression. This enhancer is activated by the T-box transcription factors TBX-37 and TBX-38, and surprisingly, concatemerized TBX-37/38 binding sites are sufficient to drive anterior-biased expression alone, despite the broad expression of TBX-37 and TBX-38. Taken together, our results highlight the diverse mechanisms used to regulate anterior expression patterns in the embryo.

A Maternal-Effect Toxin Affects Epithelial Differentiation and Tissue Mechanics in Caenorhabditis elegans

Selfish genetic elements that act as post-segregation distorters cause lethality in non-carrier individuals after fertilization. Two post-segregation distorters have been previously identified in Caenorhabditis elegans, the peel-1/zeel-1 and the sup-35/pha-1 elements. These elements seem to act as modification-rescue systems, also called toxin/antidote pairs. Here we show that the maternal-effect toxin/zygotic antidote pair sup-35/pha-1 is required for proper expression of apical junction (AJ) components in epithelia and that sup-35 toxicity increases when pathways that establish and maintain basal epithelial characteristics, die-1, elt-1, lin-26, and vab-10, are compromised. We demonstrate that pha-1(e2123) embryos, which lack the antidote, are defective in epidermal morphogenesis and frequently fail to elongate. Moreover, seam cells are frequently misshaped and mispositioned and cell bond tension is reduced in pha-1(e2123) embryos, suggesting altered tissue material properties in the epidermis. Several aspects of this phenotype can also be induced in wild-type embryos by exerting mechanical stress through uniaxial loading. Seam cell shape, tissue mechanics, and elongation can be restored in pha-1(e2123) embryos if expression of the AJ molecule DLG-1/Discs large is reduced. Thus, our experiments suggest that maternal-effect toxicity disrupts proper development of the epidermis which involves distinct transcriptional regulators and AJ components.

Engineering threshold-based selection systems

Using model organisms to identify novel therapeutic targets is frequently constrained by pre-existing genetic toolkits. To expedite positive selection for identification of novel downstream effectors, we engineered conditional expression of activated CED-10/Rac to disrupt Caenorhabditis elegans embryonic morphogenesis, titrated to 100% lethality. The strategy of engineering thresholds for positive selection using experimental animals was validated with pharmacological and genetic suppression and is generalizable to diverse molecular processes and experimental systems.

Multi-tissue patterning drives anterior morphogenesis of the C. elegans embryo

Complex structures derived from multiple tissue types are challenging to study in vivo, and our knowledge of how cells from different tissues are coordinated is limited. Model organisms have proven invaluable for improving our understanding of how chemical and mechanical cues between cells from two different tissues can govern specific morphogenetic events. Here we used Caenorhabditis elegans as a model system to show how cells from three different tissues are coordinated to give rise to the anterior lumen. While some aspects of pharyngeal morphogenesis have been well-described, it is less clear how cells from the pharynx, epidermis and neuroblasts coordinate to define the location of the anterior lumen and supporting structures. Using various microscopy and software approaches, we define the movements and patterns of these cells during anterior morphogenesis. Projections from the anterior-most pharyngeal cells (arcade cells) provide the first visible markers for the location of the future lumen, and facilitate patterning of the surrounding neuroblasts. These neuroblast patterns control the rate of migration of the anterior epidermal cells, whereas the epidermal cells ultimately reinforce and control the position of the future lumen, as they must join with the pharyngeal cells for their epithelialization. Our studies are the first to characterize anterior morphogenesis in C. elegans in detail and should lay the framework for identifying how these different patterns are controlled at the molecular level.

Game of Tissues: How the Epidermis Thrones C. elegans Shape

The versatility of epithelial cell structure is universally exploited by organisms in multiple contexts. Epithelial cells can establish diverse polarized axes within their tridimensional structure which enables them to flexibly communicate with their neighbors in a 360° range. Hence, these cells are central to multicellularity, and participate in diverse biological processes such as organismal development, growth or immune response and their misfunction ultimately impacts disease. During the development of an organism, the first task epidermal cells must complete is the formation of a continuous sheet, which initiates its own morphogenic process. In this review, we will focus on the C. elegans embryonic epithelial morphogenesis. We will describe how its formation, maturation, and spatial arrangements set the final shape of the nematode C. elegans. Special importance will be given to the tissue-tissue interactions, regulatory tissue-tissue feedback mechanisms and the players orchestrating the process.

DOT1L complex suppresses transcription from enhancer elements and ectopic RNAi in Caenorhabditis elegans

Methylation of histone H3 on lysine 79 (H3K79) by DOT1L is associated with actively transcribed genes. Earlier, we described that DOT-1.1, the Caenorhabditis elegans homolog of mammalian DOT1L, cooperates with the chromatin-binding protein ZFP-1 (AF10 homolog) to negatively modulate transcription of highly and widely expressed target genes. Also, the reduction of ZFP-1 levels has consistently been associated with lower efficiency of RNA interference (RNAi) triggered by exogenous double-stranded RNA (dsRNA), but the reason for this is not clear. Here, we demonstrate that the DOT1L complex suppresses transcription originating from enhancer elements and antisense transcription, thus potentiating the expression of enhancer-regulated genes. We also show that worms lacking H3K79 methylation do not survive, and this lethality is suppressed by a loss of caspase-3 or Dicer complex components that initiate gene silencing response to exogenous dsRNA. Our results suggest that ectopic elevation of endogenous dsRNA directly or indirectly resulting from global misregulation of transcription in DOT1L complex mutants may engage the Dicer complex and, therefore, limit the efficiency of exogenous RNAi.

A multicellular rosette-mediated collective dendrite extension

Coordination of neurite morphogenesis with surrounding tissues is crucial to the establishment of neural circuits, but the underlying cellular and molecular mechanisms remain poorly understood. We show that neurons in a C. elegans sensory organ, called the amphid, undergo a collective dendrite extension to form the sensory nerve. The amphid neurons first assemble into a multicellular rosette. The vertex of the rosette, which becomes the dendrite tips, is attached to the anteriorly migrating epidermis and carried to the sensory depression, extruding the dendrites away from the neuronal cell bodies. Multiple adhesion molecules including DYF-7, SAX-7, HMR-1 and DLG-1 function redundantly in rosette-to-epidermis attachment. PAR-6 is localized to the rosette vertex and dendrite tips, and promotes DYF-7 localization and dendrite extension. Our results suggest a collective mechanism of neurite extension that is distinct from the classical pioneer-follower model and highlight the role of mechanical cues from surrounding tissues in shaping neurites.

Cdc42 regulates junctional actin but not cell polarization in the Caenorhabditis elegans epidermis

During morphogenesis, adherens junctions (AJs) remodel to allow changes in cell shape and position while preserving adhesion. Here, we examine the function of Rho guanosine triphosphatase CDC-42 in AJ formation and regulation during Caenorhabditis elegans embryo elongation, a process driven by asymmetric epidermal cell shape changes. cdc-42 mutant embryos arrest during elongation with epidermal ruptures. Unexpectedly, we find using time-lapse fluorescence imaging that cdc-42 is not required for epidermal cell polarization or junction assembly, but rather is needed for proper junctional actin regulation during elongation. We show that the RhoGAP PAC-1/ARHGAP21 inhibits CDC-42 activity at AJs, and loss of PAC-1 or the interacting linker protein PICC-1/CCDC85A-C blocks elongation in embryos with compromised AJ function. pac-1 embryos exhibit dynamic accumulations of junctional F-actin and an increase in AJ protein levels. Our findings identify a previously unrecognized molecular mechanism for inhibiting junctional CDC-42 to control actin organization and AJ protein levels during epithelial morphogenesis.

The C. elegans mRNA decapping enzyme shapes morphology of cilia

Cilia and flagella are evolutionarily conserved organelles that protrude from cell surfaces. Most cilia and flagella are single rod-shaped but some cilia show a variety of shapes. For example, human airway epithelial cells are multiciliated, flagella of crayfish spermatozoon are star-like shaped, and fruit fly spermatozoon extends long flagella. In Caenorhabditis elegans, cilia display morphological diversity of shapes (single, dual rod-type and wing-like and highly-branched shapes). Here we show that DCAP-1 and DCAP-2, which are the homologues of mammalian DCP1 and DCP2 mRNA decapping enzymes, respectively, are involved in formation of dual rod-type and wing-like shaped cilia in C. elegans. mRNA decapping enzyme catalyzes hydrolysis of 5' cap structure of mRNA, which leads to degradation of mRNA. Rescue experiments showed that DCAP-2 acts not in glial cells surrounding cilia but in neurons. This is the first evidence to demonstrate that mRNA decapping is involved in ciliary shape formation.

Genome-wide discovery of active regulatory elements and transcription factor footprints in Caenorhabditis elegans using DNase-seq

Deep sequencing of size-selected DNase I–treated chromatin (DNase-seq) allows high-resolution measurement of chromatin accessibility to DNase I cleavage, permitting identification of de novo active cis-regulatory modules (CRMs) and individual transcription factor (TF) binding sites. We adapted DNase-seq to nuclei isolated from C. elegans embryos and L1 arrest larvae to generate high-resolution maps of TF binding. Over half of embryonic DNase I hypersensitive sites (DHSs) were annotated as noncoding, with 24% in intergenic, 12% in promoters, and 28% in introns, with similar statistics observed in L1 arrest larvae. Noncoding DHSs are highly conserved and enriched in marks of enhancer activity and transcription. We validated noncoding DHSs against known enhancers from myo-2, myo-3, hlh-1, elt-2, and lin-26/lir-1 and recapitulated 15 of 17 known enhancers. We then mined DNase-seq data to identify putative active CRMs and TF footprints. Using DNase-seq data improved predictions of tissue-specific expression compared with motifs alone. In a pilot functional test, 10 of 15 DHSs from pha-4, icl-1, and ceh-13 drove reporter gene expression in transgenic C. elegans. Overall, we provide experimental annotation of 26,644 putative CRMs in the embryo containing 55,890 TF footprints, as well as 15,841 putative CRMs in the L1 arrest larvae containing 32,685 TF footprints.

A CD47-associated super-enhancer links pro-inflammatory signalling to CD47 upregulation in breast cancer

CD47 is a cell surface molecule that inhibits phagocytosis of cells that express it by binding to its receptor, SIRPα, on macrophages and other immune cells. CD47 is expressed at different levels by neoplastic and normal cells. Here, to reveal mechanisms by which different neoplastic cells generate this dominant 'don't eat me' signal, we analyse the CD47 regulatory genomic landscape. We identify two distinct super-enhancers (SEs) associated with CD47 in certain cancer cell types. We show that a set of active constituent enhancers, located within the two CD47 SEs, regulate CD47 expression in different cancer cell types and that disruption of CD47 SEs reduces CD47 gene expression. Finally we report that the TNF-NFKB1 signalling pathway directly regulates CD47 by interacting with a constituent enhancer located within a CD47-associated SE specific to breast cancer. These results suggest that cancers can evolve SE to drive CD47 overexpression to escape immune surveillance.

Non-centrosomal epidermal microtubules act in parallel to LET-502/ROCK to promote C. elegans elongation

C. elegans embryonic elongation is a morphogenetic event driven by actomyosin contractility and muscle-induced tension transmitted through hemidesmosomes. A role for the microtubule cytoskeleton has also been proposed, but its contribution remains poorly characterized. Here, we investigate the organization of the non-centrosomal microtubule arrays present in the epidermis and assess their function in elongation. We show that the microtubule regulators γ-tubulin and NOCA-1 are recruited to hemidesmosomes and adherens junctions early in elongation. Several parallel approaches suggest that microtubule nucleation occurs from these sites. Disrupting the epidermal microtubule array by overexpressing the microtubule-severing protein Spastin or by inhibiting the C. elegans ninein homolog NOCA-1 in the epidermis mildly affected elongation. However, microtubules were essential for elongation when hemidesmosomes or the activity of the Rho kinase LET-502/ROCK were partially compromised. Imaging of junctional components and genetic analyses suggest that epidermal microtubules function together with Rho kinase to promote the transport of E-cadherin to adherens junctions and myotactin to hemidesmosomes. Our results indicate that the role of LET-502 in junctional remodeling is likely to be independent of its established function as a myosin II activator, but requires a microtubule-dependent pathway involving the syntaxin SYX-5. Hence, we propose that non-centrosomal microtubules organized by epidermal junctions contribute to elongation by transporting junction remodeling factors, rather than having a mechanical role.

Summary: During C. elegans embryonic elongation, microtubules nucleate at adjerens junctions and hemidesmosomes, and are important for the transport of junctional proteins.

Polarized Rac-dependent protrusions drive epithelial intercalation in the embryonic epidermis of C. elegans

Cell intercalation is a fundamental, coordinated cell rearrangement process that shapes tissues throughout animal development. Studies of intercalation within epithelia have focused almost exclusively on the localized constriction of specific apical junctions. Another widely deployed yet poorly understood alternative mechanism of epithelial intercalation relies on basolateral protrusive activity. Using the dorsal embryonic epidermis of Caenorhabditis elegans, we have investigated this alternative mechanism using high-resolution live cell microscopy and genetic analysis. We find that as dorsal epidermal cells migrate past one another they produce F-actin-rich protrusions polarized at their extending (medial) edges. These protrusions are controlled by the C. elegans Rac and RhoG orthologs CED-10 and MIG-2, which function redundantly to polarize actin polymerization upstream of the WAVE complex and WASP, respectively. We also identify UNC-73, the C. elegans ortholog of Trio, as a guanine nucleotide exchange factor (GEF) upstream of both CED-10 and MIG-2. Further, we identify a novel polarizing cue, CRML-1, which is the ortholog of human capping Arp2/3 myosin I linker (CARMIL), that localizes to the nonprotrusive lateral edges of dorsal cells. CRML-1 genetically suppresses UNC-73 function and, indirectly, actin polymerization. This network identifies a novel, molecularly conserved cassette that regulates epithelial intercalation via basolateral protrusive activity.

Glial development and function in the nervous system of Caenorhabditis elegans

The nematode, Caenorhabditis elegans, has served as a fruitful setting for understanding conserved biological processes. The past decade has seen the rise of this model organism as an important tool for uncovering the mysteries of the glial cell, which partners with neurons to generate a functioning nervous system in all animals. C. elegans affords unparalleled single-cell resolution in vivo in examining glia-neuron interactions, and similarities between C. elegans and vertebrate glia suggest that lessons learned from this nematode are likely to have general implications. Here, I summarize what has been gleaned over the past decade since C. elegans glia research became a concerted area of focus. Studies have revealed that glia are essential elements of a functioning C. elegans nervous system and play key roles in its development. Importantly, glial influence on neuronal function appears to be dynamic. Key questions for the field to address in the near- and long-term have emerged, and these are discussed within.

Multiple transcription factors directly regulate Hox gene lin-39 expression in ventral hypodermal cells of the C. elegans embryo and larva, including the hypodermal fate regulators LIN-26 and ELT-6

Background

Hox genes encode master regulators of regional fate specification during early metazoan development. Much is known about the initiation and regulation of Hox gene expression in Drosophila and vertebrates, but less is known in the non-arthropod invertebrate model system, C. elegans. The C. elegans Hox gene lin-39 is required for correct fate specification in the midbody region, including the Vulval Precursor Cells (VPCs). To better understand lin-39 regulation and function, we aimed to identify transcription factors necessary for lin-39 expression in the VPCs, and in particular sought factors that initiate lin-39 expression in the embryo.

Results

We used the yeast one-hybrid (Y1H) method to screen for factors that bound to 13 fragments from the lin-39 region: twelve fragments contained sequences conserved between C. elegans and two other nematode species, while one fragment was known to drive reporter gene expression in the early embryo in cells that generate the VPCs. Sixteen transcription factors that bind to eight lin-39 genomic fragments were identified in yeast, and we characterized several factors by verifying their physical interactions in vitro, and showing that reduction of their function leads to alterations in lin-39 levels and lin-39::GFP reporter expression in vivo. Three factors, the orphan nuclear hormone receptor NHR-43, the hypodermal fate regulator LIN-26, and the GATA factor ELT-6 positively regulate lin-39 expression in the embryonic precursors to the VPCs. In particular, ELT-6 interacts with an enhancer that drives GFP expression in the early embryo, and the ELT-6 site we identified is necessary for proper embryonic expression. These three factors, along with the factors ZTF-17, BED-3 and TBX-9, also positively regulate lin-39 expression in the larval VPCs.

Conclusions

These results significantly expand the number of factors known to directly bind and regulate lin-39 expression, identify the first factors required for lin-39 expression in the embryo, and hint at a positive feedback mechanism involving GATA factors that maintains lin-39 expression in the vulval lineage. This work indicates that, as in other organisms, the regulation of Hox gene expression in C. elegans is complicated, redundant and robust.

The PAF1 complex is involved in embryonic epidermal morphogenesis in Caenorhabditis elegans

The PAF1 complex (PAF1C) is an evolutionarily conserved protein complex involved in transcriptional regulation and chromatin remodeling. How the PAF1C is involved in animal development is still not well understood. Here, we report that, in the nematode Caenorhabditis elegans, the PAF1C is involved in epidermal morphogenesis in late embryogenesis. From an RNAi screen we identified the C. elegans ortholog of a component of the PAF1C, CTR-9, as a gene whose depletion caused various defects during embryonic epidermal morphogenesis, including epidermal cell positioning, ventral enclosure and epidermal elongation. RNAi of orthologs of other four components of the PAF1C (PAFO-1, LEO-1, CDC-73 and RTFO-1) caused similar epidermal defects. In these embryos, whereas the number and cell fate determination of epidermal cells were apparently unaffected, their position and shape were severely disorganized. PAFO-1::mCherry, mCherry::LEO-1 and GFP::RTFO-1 driven by the authentic promoters were detected in the nuclei of a wide range of cells. Nuclear localization of GFP::RTFO-1 was independent of other PAF1C components, while PAFO-1::mCherry and mCherry::LEO-1 dependent on other components except RTFO-1. Epidermis-specific expression of mCherry::LEO-1 rescued embryonic lethality of the leo-1 deletion mutant. Thus, although the PAF1C is universally expressed in C. elegans embryos, its epidermal function is crucial for the viability of this animal.

The landscape of RNA polymerase II transcription initiation in C. elegans reveals promoter and enhancer architectures

RNA polymerase transcription initiation sites are largely unknown in Caenorhabditis elegans. The initial 5′ end of most protein-coding transcripts is removed by trans-splicing, and noncoding initiation sites have not been investigated. We characterized the landscape of RNA Pol II transcription initiation, identifying 73,500 distinct clusters of initiation. Bidirectional transcription is frequent, with a peak of transcriptional pairing at 120 bp. We assign transcription initiation sites to 7691 protein-coding genes and find that they display features typical of eukaryotic promoters. Strikingly, the majority of initiation events occur in regions with enhancer-like chromatin signatures. Based on the overlap of transcription initiation clusters with mapped transcription factor binding sites, we define 2361 transcribed intergenic enhancers. Remarkably, productive transcription elongation across these enhancers is predominantly in the same orientation as that of the nearest downstream gene. Directed elongation from an upstream enhancer toward a downstream gene could potentially deliver RNA polymerase II to a proximal promoter, or alternatively might function directly as a distal promoter. Our results provide a new resource to investigate transcription regulation in metazoans.

C. elegans dystroglycan coordinates responsiveness of follower axons to dorsal/ventral and anterior/posterior guidance cues

Neural development in metazoans is characterized by the establishment of initial process tracts by pioneer axons and the subsequent extension of follower axons along these pioneer processes. Mechanisms governing the fidelity of follower extension along pioneered routes are largely unknown. In C. elegans, formation of the right angle-shaped lumbar commissure connecting the lumbar and preanal ganglia is an example of pioneer/follower dynamics. We find that the dystroglycan ortholog DGN-1 mediates the fidelity of follower lumbar commissure axon extension along the pioneer axon route. In dgn-1 mutants, the axon of the pioneer PVQ neuron faithfully establishes the lumbar commissure, but axons of follower lumbar neurons, such as PVC, frequently bypass the lumbar commissure and extend along an oblique trajectory directly toward the preanal ganglion. In contrast, disruption of the UNC-6/netrin guidance pathway principally perturbs PVQ ventral guidance to pioneer the lumbar commissure. Loss of DGN-1 in unc-6 mutants has a quantitatively similar effect on follower axon guidance regardless of PVQ axon route, indicating that DGN-1 does not mediate follower/pioneer adhesion. Instead, DGN-1 appears to block premature responsiveness of follower axons to a preanal ganglion-directed guidance cue, which mediates ventral-to-anterior reorientation of lumbar commissure axons. Deletion analysis shows that only the most N-terminal DGN-1 domain is required for these activities. These studies suggest that dystroglycan modulation of growth cone responsiveness to conflicting guidance cues is important for restricting follower axon extension to the tracts laid down by pioneers.

A regulatory network modeled from wild-type gene expression data guides functional predictions in Caenorhabditis elegans development

BACKGROUND

Complex gene regulatory networks underlie many cellular and developmental processes. While a variety of experimental approaches can be used to discover how genes interact, few biological systems have been systematically evaluated to the extent required for an experimental definition of the underlying network. Therefore, the development of computational methods that can use limited experimental data to define and model a gene regulatory network would provide a useful tool to evaluate many important but incompletely understood biological processes. Such methods can assist in extracting all relevant information from data that are available, identify unexpected regulatory relationships and prioritize future experiments.

RESULTS

To facilitate the analysis of gene regulatory networks, we have developed a computational modeling pipeline method that complements traditional evaluation of experimental data. For a proof-of-concept example, we have focused on the gene regulatory network in the nematode C. elegans that mediates the developmental choice between mesodermal (muscle) and ectodermal (skin) cell fates in the embryonic C lineage. We have used gene expression data to build two models: a knowledge-driven model based on gene expression changes following gene perturbation experiments, and a data-driven mathematical model derived from time-course gene expression data recovered from wild-type animals. We show that both models can identify a rich set of network gene interactions. Importantly, the mathematical model built only from wild-type data can predict interactions demonstrated by the perturbation experiments better than chance, and better than an existing knowledge-driven model built from the same data set. The mathematical model also provides new biological insight, including a dissection of zygotic from maternal functions of a key transcriptional regulator, PAL-1, and identification of non-redundant activities of the T-box genes tbx-8 and tbx-9.

CONCLUSIONS

This work provides a strong example for a mathematical modeling approach that solely uses wild-type data to predict an underlying gene regulatory network. The modeling approach complements traditional methods of data analysis, suggesting non-intuitive network relationships and guiding future experiments.

Some, but not all, retromer components promote morphogenesis of C. elegans sensory compartments

The endings of sensory receptor cells often lie within specialized compartments formed by glial cells. The main sensory organ of Caenorhabditis elegans, the amphid, provides a powerful setting for studying glial compartment morphogenesis. Our previous studies showed that amphid compartment size is controlled by opposing activities of the Nemo-like kinase LIT-1, which promotes compartment expansion, and the Patched-related protein DAF-6, which restricts compartment growth. From a genetic screen for mutations able to suppress the bloated sensory compartments of daf-6 mutants, we identified an allele of the sorting nexin gene snx-1. SNX-1 protein is a component of the retromer, a protein complex that facilitates recycling of transmembrane proteins from the endosome to the Golgi network. We find that snx-1 functions cell autonomously within glia to promote sensory compartment growth, and that SNX-1 protein is enriched near the surface of the sensory compartment. snx-1 interacts genetically with lit-1 and another regulator of compartment size, the Dispatched-related gene che-14. Mutations in snx-3 and vps-29, also retromer genes, can suppress daf-6 defects. Surprisingly, however, remaining retromer components seem not to be involved. Our results suggest that a novel assembly of retromer components is important for determining sensory compartment dimensions.

Sensory organ remodeling in Caenorhabditis elegans requires the zinc-finger protein ZTF-16

Neurons and glia display remarkable morphological plasticity, and remodeling of glia may facilitate neuronal shape changes. The molecular basis and control of glial shape changes is not well understood. In response to environmental stress, the nematode Caenorhabditis elegans enters an alternative developmental state, called dauer, in which glia and neurons of the amphid sensory organ remodel. Here, we describe a genetic screen aimed at identifying genes required for amphid glia remodeling. We previously demonstrated that remodeling requires the Otx-type transcription factor TTX-1 and its direct target, the receptor tyrosine kinase gene ver-1. We now find that the hunchback/Ikaros-like C2H2 zinc-finger factor ztf-16 is also required. We show that ztf-16 mutants exhibit pronounced remodeling defects, which are explained, at least in part, by defects in the expression of ver-1. Expression and cell-specific rescue studies suggest that ztf-16, like ttx-1, functions within glia; however, promoter deletion studies show that ztf-16 acts through a site on the ver-1 promoter that is independent of ttx-1. Our studies identify an important component of glia remodeling and suggest that transcriptional changes may underlie glial morphological plasticity in the sensory organs of C. elegans.

Neural maintenance roles for the matrix receptor dystroglycan and the nuclear anchorage complex in Caenorhabditis elegans

Recent studies in Caenorhabditis elegans have revealed specific neural maintenance mechanisms that protect soma and neurites against mispositioning due to displacement stresses, such as muscle contraction. We report that C. elegans dystroglycan (DG) DGN-1 functions to maintain the position of lumbar neurons during late embryonic and larval development. In the absence of DGN-1 the cell bodies of multiple lumbar neuron classes are frequently displaced anterior of their normal positions. Early but not later embryonic panneural expression of DGN-1 rescues positional maintenance, suggesting that dystroglycan is required for establishment of a critical maintenance pathway that persists throughout later developmental stages. Lumbar neural maintenance requires only a membrane-tethered N-terminal domain of DGN-1 and may involve a novel extracellular partner for dystroglycan. A genetic screen for similar lumbar maintenance mutants revealed a role for the nesprin/SYNE family protein ANC-1 as well as for the extracellular protein DIG-1, previously implicated in lumbar neuron maintenance. The involvement of ANC-1 reveals a previously unknown role for nucleus-cytoskeleton interactions in neural maintenance. Genetic analysis indicates that lumbar neuron position is maintained in late embryos by parallel DGN-1/DIG-1 and ANC-1-dependent pathways, and in larvae by separate DGN-1 and ANC-1 pathways. The effect of muscle paralysis on late embryonic- or larval-stage maintenance defects in mutants indicates that lumbar neurons are subject to both muscle contraction-dependent and contraction-independent displacement stresses, and that different maintenance pathways may protect against specific types of displacement stress.

Cell plasticity in Caenorhabditis elegans: from induced to natural cell reprogramming

Achieving controlled reprogramming of differentiated cells into a desired cell type would open new opportunities in stem-cell biology and regenerative medicine. Experimentation on cell reprogramming requires a model in which cell conversion can be induced and tracked individually. The tiny nematode, Caenorhabditis elegans, owing to its known cellular lineage, allows the study of direct cell type conversion with a single-cell resolution. Indeed, recent advances have shown that despite its invariant cell lineage, cellular identities can be reprogrammed, leading to cell conversion in vivo. In addition, natural transdifferentiation events occur in the worm, providing a powerful model for the study of cellular plasticity in a physiological cellular microenvironment. Here, we review pioneer studies on induced and naturally occurring reprogramming events in C. elegans and the new notions that have emerged.

Opposing activities of LIT-1/NLK and DAF-6/patched-related direct sensory compartment morphogenesis in C. elegans

Glial cells surround neuronal endings to create enclosed compartments required for neuronal function. This architecture is seen at excitatory synapses and at sensory neuron receptive endings. Despite the prevalence and importance of these compartments, how they form is not known. We used the main sensory organ of C. elegans, the amphid, to investigate this issue. daf-6/Patched-related is a glia-expressed gene previously implicated in amphid sensory compartment morphogenesis. By comparing time series of electron-microscopy (EM) reconstructions of wild-type and daf-6 mutant embryos, we show that daf-6 acts to restrict compartment size. From a genetic screen, we found that mutations in the gene lit-1/Nemo-like kinase (NLK) suppress daf-6. EM and genetic studies demonstrate that lit-1 acts within glia, in counterbalance to daf-6, to promote sensory compartment expansion. Although LIT-1 has been shown to regulate Wnt signaling, our genetic studies demonstrate a novel, Wnt-independent role for LIT-1 in sensory compartment size control. The LIT-1 activator MOM-4/TAK1 is also important for compartment morphogenesis and both proteins line the glial sensory compartment. LIT-1 compartment localization is important for its function and requires neuronal signals. Furthermore, the conserved LIT-1 C-terminus is necessary and sufficient for this localization. Two-hybrid and co-immunoprecipitation studies demonstrate that the LIT-1 C-terminus binds both actin and the Wiskott-Aldrich syndrome protein (WASP), an actin regulator. We use fluorescence light microscopy and fluorescence EM methodology to show that actin is highly enriched around the amphid sensory compartment. Finally, our genetic studies demonstrate that WASP is important for compartment expansion and functions in the same pathway as LIT-1. The studies presented here uncover a novel, Wnt-independent role for the conserved Nemo-like kinase LIT-1 in controlling cell morphogenesis in conjunction with the actin cytoskeleton. Our results suggest that the opposing daf-6 and lit-1 glial pathways act together to control sensory compartment size.

The nervous system of most animals consists of two related cell types, neurons and glia. A striking property of glia is their ability to ensheath neuronal cells, which can help increase the efficiency of synaptic communication between neurons. Sensory neuron receptive endings in the periphery, as well as excitatory synapses in the central nervous system, often lie within specialized compartments formed by glial processes. Despite the prevalence of these compartments, and their importance for neuronal function and signal transmission, little is known about how they form. We have used the amphid, the main sensory organ of the worm Caenorhabditis elegans, to investigate glial sensory compartment morphogenesis. We demonstrate that the glia-expressed gene daf-6/Patched-related acts to restrict the size of the sensory compartment, while the Nemo-like kinase lit-1 acts within glia in the opposite direction, to promote sensory compartment expansion. We show that LIT-1 localizes to the sensory compartment through a highly conserved domain. This domain can interact both with actin, which outlines the compartment, and with the regulator of actin polymerization WASP, which acts in the same pathway as lit-1. We postulate that Nemo-like kinases could have broader roles as regulators of cellular morphogenesis, in addition to their traditional role in regulating the Wnt signaling pathway.

A novel sperm-delivered toxin causes late-stage embryo lethality and transmission ratio distortion in C. elegans

The evolutionary fate of an allele ordinarily depends on its contribution to host fitness. Occasionally, however, genetic elements arise that are able to gain a transmission advantage while simultaneously imposing a fitness cost on their hosts. We previously discovered one such element in C. elegans that gains a transmission advantage through a combination of paternal-effect killing and zygotic self-rescue. Here we demonstrate that this element is composed of a sperm-delivered toxin, peel-1, and an embryo-expressed antidote, zeel-1. peel-1 and zeel-1 are located adjacent to one another in the genome and co-occur in an insertion/deletion polymorphism. peel-1 encodes a novel four-pass transmembrane protein that is expressed in sperm and delivered to the embryo via specialized, sperm-specific vesicles. In the absence of zeel-1, sperm-delivered PEEL-1 causes lethal defects in muscle and epidermal tissue at the 2-fold stage of embryogenesis. zeel-1 is expressed transiently in the embryo and encodes a novel six-pass transmembrane domain fused to a domain with sequence similarity to zyg-11, a substrate-recognition subunit of an E3 ubiquitin ligase. zeel-1 appears to have arisen recently, during an expansion of the zyg-11 family, and the transmembrane domain of zeel-1 is required and partially sufficient for antidote activity. Although PEEL-1 and ZEEL-1 normally function in embryos, these proteins can act at other stages as well. When expressed ectopically in adults, PEEL-1 kills a variety of cell types, and ectopic expression of ZEEL-1 rescues these effects. Our results demonstrate that the tight physical linkage between two novel transmembrane proteins has facilitated their co-evolution into an element capable of promoting its own transmission to the detriment of organisms carrying it.

Regulation of sperm gene expression by the GATA factor ELT-1

Cell fate specification is mediated primarily through the expression of cell-type-specific genes. The regulatory pathway that governs the sperm/egg decision in the hermaphrodite germ line of Caenorhabditis elegans has been well characterized, but the transcription factors that drive these developmental programs remain unknown. We report the identification of ELT-1, a GATA transcription factor that specifies hypodermal fate in the embryo, as a regulator of sperm-specific transcription in the germ line. Computational analysis identified a conserved bipartite sequence element that is found almost exclusively in the promoters of a number of sperm genes. ELT-1 was recovered in a yeast one-hybrid screen for factors that bind to that sperm consensus site. In vitro assays defined the sperm consensus sequence as an optimal binding site for ELT-1. We determined that expression of elt-1 is elevated in the sperm-producing germ line, and that ELT-1 is required for sperm function. Deletion of the ELT-1 binding site from a sperm promoter abrogates sperm-specific expression of a reporter transgene. This work demonstrates a role for the ELT-1 transcription factor in sperm, and provides a critical link between the germ line sex determination program and gamete-specific gene expression.

Myosin II regulation during C. elegans embryonic elongation: LET-502/ROCK, MRCK-1 and PAK-1, three kinases with different roles

Myosin II plays a central role in epithelial morphogenesis; however, its role has mainly been examined in processes involving a single cell type. Here we analyze the structure, spatial requirement and regulation of myosin II during C. elegans embryonic elongation, a process that involves distinct epidermal cells and muscles. We developed novel GFP probes to visualize the dynamics of actomyosin remodeling, and found that the assembly of myosin II filaments, but not actin microfilaments, depends on the myosin regulatory light chain (MLC-4) and essential light chain (MLC-5, which we identified herein). To determine how myosin II regulates embryonic elongation, we rescued mlc-4 mutants with various constructs and found that MLC-4 is essential in a subset of epidermal cells. We show that phosphorylation of two evolutionary conserved MLC-4 serine and threonine residues is important for myosin II activity and organization. Finally, in an RNAi screen for potential myosin regulatory light chain kinases, we found that the ROCK, PAK and MRCK homologs act redundantly. The combined loss of ROCK and PAK, or ROCK and MRCK, completely prevented embryonic elongation, but a constitutively active form of MLC-4 could only rescue a lack of MRCK. This result, together with systematic genetic epistasis tests with a myosin phosphatase mutation, suggests that ROCK and MRCK regulate MLC-4 and the myosin phosphatase. Moreover, we suggest that ROCK and PAK regulate at least one other target essential for elongation, in addition to MLC-4.

Determination of the fate and contribution of ex vivo expanded human bone marrow stem and progenitor cells for bone formation by 2.3ColGFP

Bone marrow transplantation can provide an effective cell-based strategy to enhance bone repair. However, the fate of implanted cells and the extent of their contribution to bone osteoinduction remain uncertain. To define the fate of bone marrow-derived cells and their contribution in vivo, we used a bone-specific collagen I promoter (2.3Col) driving green fluorescent protein (GFP) (2.3ColGFP) within a lentiviral vector. Prior to in vivo cell fate determination, we verified a high efficiency of lentiviral transduction in human bone marrow stromal cells (hBMSCs), without altering the proliferation or differentiation potential of these cells. We showed that the 2.3ColGFP marker responded to endogenous transcriptional regulation signals. In a mouse ossicle model, we demonstrated that the 2.3ColGFP marker is able to specifically define human bone marrow-derived stem cells that enter the osteoblast lineage in vivo. In addition, cells labeled with 2.3ColGFP with the donor origin, directly make a major contribution to bone formation. Furthermore, we also demonstrated in a calvarial defect model that a mixture of human bone marrow-derived populations, have stronger bone regenerative potential than that of hBMSCs, and an optimal dose is required for bone regeneration by the mixed populations.

DEX-1 and DYF-7 establish sensory dendrite length by anchoring dendritic tips during cell migration

Cells are devices whose structures delimit function. For example, in the nervous system, neuronal and glial shapes dictate paths of information flow. To understand how cells acquire their shapes, we examined the formation of a sense organ in C. elegans. Using time-lapse imaging, we found that sensory dendrites form by stationary anchoring of dendritic tips during cell-body migration. A genetic screen identified DEX-1 and DYF-7, extracellular proteins required for dendritic tip anchoring, which act cooperatively at the time and place of anchoring. DEX-1 and DYF-7 contain, respectively, zonadhesin and zona pellucida domains, and DYF-7 self-associates into multimers important for anchoring. Thus, unlike other dendrites, amphid dendritic tips are positioned by DEX-1 and DYF-7 without the need for long-range guidance cues. In sequence and function, DEX-1 and DYF-7 resemble tectorins, which anchor stereocilia in the inner ear, suggesting that a sensory dendrite anchor may have evolved into part of a mechanosensor.

Functional analysis of synaptotagmin gene regulatory regions in two distantly related ascidian species

We have studied the structure and function of a promoter region of the Halocynthia synaptotagmin (Hr-Syt) gene, which is abundantly expressed in neuronal cells. Our previous analysis suggested that the expression of Hr-Syt is regulated by at least one epidermal and two neuronal regulatory regions. In this study, the regulatory regions of Hr-Syt promoter were further characterized by using two species of ascidians, Halocynthia roretzi and Ciona intestinalis embryos. A putative GATA transcription factor binding site in the epidermal regulatory region has ectodermal enhancer activity in the Halocynthia embryo. Neuronal expression of Hr-Syt was regulated by multiple redundant enhancer regions. Among these enhancer regions, a 200-bp (-2900/-2700) region drove the reporter expression in neurons in both species of ascidian. Although the synaptotagmin promoter sequences did not show overall similarity between Hr-Syt and Ciona synaptotagmin (Ci-Syt), 5'-upsteream two short sequences of Ci-Syt have similarity to the -2766/-2732 region of the Hr-Syt promoter. The homeodomain binding sites in this region are required for the neuronal enhancer activity. These results suggest that GATA and homeodomain transcription factors regulate the expression of synaptotagmin.

A C. elegans genome-scale microRNA network contains composite feedback motifs with high flux capacity

MicroRNAs (miRNAs) and transcription factors (TFs) are primary metazoan gene regulators. Whereas much attention has focused on finding the targets of both miRNAs and TFs, the transcriptional networks that regulate miRNA expression remain largely unexplored. Here, we present the first genome-scale Caenorhabditis elegans miRNA regulatory network that contains experimentally mapped transcriptional TF --> miRNA interactions, as well as computationally predicted post-transcriptional miRNA --> TF interactions. We find that this integrated miRNA network contains 23 miRNA <--> TF composite feedback loops in which a TF that controls a miRNA is itself regulated by that same miRNA. By rigorous network randomizations, we show that such loops occur more frequently than expected by chance and, hence, constitute a genuine network motif. Interestingly, miRNAs and TFs in such loops are heavily regulated and regulate many targets. This "high flux capacity" suggests that loops provide a mechanism of high information flow for the coordinate and adaptable control of miRNA and TF target regulons.

Pairing of competitive and topologically distinct regulatory modules enhances patterned gene expression

Biological networks are inherently modular, yet little is known about how modules are assembled to enable coordinated and complex functions. We used RNAi and time series, whole-genome microarray analyses to systematically perturb and characterize components of a Caenorhabditis elegans lineage-specific transcriptional regulatory network. These data are supported by selected reporter gene analyses and comprehensive yeast one-hybrid and promoter sequence analyses. Based on these results, we define and characterize two modules composed of muscle- and epidermal-specifying transcription factors that function together within a single cell lineage to robustly specify multiple cell types. The expression of these two modules, although positively regulated by a common factor, is reliably segregated among daughter cells. Our analyses indicate that these modules repress each other, and we propose that this cross-inhibition coupled with their relative time of induction function to enhance the initial asymmetry in their expression patterns, thus leading to the observed invariant gene expression patterns and cell lineage. The coupling of asynchronous and topologically distinct modules may be a general principle of module assembly that functions to potentiate genetic switches.

Identification of muscle-specific regulatory modules in Caenorhabditis elegans

Transcriptional regulation is the major regulatory mechanism that controls the spatial and temporal expression of genes during development. This is carried out by transcription factors (TFs), which recognize and bind to their cognate binding sites. Recent studies suggest a modular organization of TF-binding sites, in which clusters of transcription-factor binding sites cooperate in the regulation of downstream gene expression. In this study, we report our computational identification and experimental verification of muscle-specific cis-regulatory modules in Caenorhabditis elegans. We first identified a set of motifs that are correlated with muscle-specific gene expression. We then predicted muscle-specific regulatory modules based on clusters of those motifs with characteristics similar to a collection of well-studied modules in other species. The method correctly identifies 88% of the experimentally characterized modules with a positive predictive value of at least 65%. The prediction accuracy of muscle-specific expression on an independent test set is highly significant (P<0.0001). We performed in vivo experimental tests of 12 predicted modules, and 10 of those drive muscle-specific gene expression. These results suggest that our method is highly accurate in identifying functional sequences important for muscle-specific gene expression and is a valuable tool for guiding experimental designs.

PHA-4/FoxA cooperates with TAM-1/TRIM to regulate cell fate restriction in the C. elegans foregut

A key question in development is how pluripotent progenitors are progressively restricted to acquire specific cell fates. Here we investigate how embryonic blastomeres in C. elegans develop into foregut (pharynx) cells in response to the selector gene PHA-4/FoxA. When pha-4 is removed from pharyngeal precursors, they exhibit two alternative responses. Before late-gastrulation (8E stage), these cells lose their pharyngeal identity and acquire an alternative fate such as ectoderm (Specification stage). After the Specification stage, mutant cells develop into aberrant pharyngeal cells (Morphogenesis/Differentiation stage). Two lines of evidence suggest that the Specification stage depends on transcriptional repression of ectodermal genes by pha-4. First, pha-4 exhibits strong synthetic phenotypes with the B class synMuv gene tam-1 (Tandam Array expression Modifier 1) and with a mediator of transcriptional repression, the NuRD complex (NUcleosome Remodeling and histone Deacetylase). Second, pha-4 associates with the promoter of the ectodermal regulator lin-26 and is required to repress lin-26 expression. We propose that restriction of early blastomeres to the pharyngeal fate depends on both repression of ectodermal genes and activation of pharyngeal genes by PHA-4.

MOM-5 frizzled regulates the distribution of DSH-2 to control C. elegans asymmetric neuroblast divisions

Asymmetric cell divisions produce all 302 neurons of the C. elegans hermaphrodite. Here, we describe a role for a C. elegans Dishevelled homolog, DSH-2, in an asymmetric neuroblast division. In dsh-2 mutants, neurons normally descended from the anterior neuroblast daughter of the ABpl/rpppa blast cell were frequently duplicated, while non-neuronal cells produced by the posterior daughter cell were often missing. These observations indicate that in the absence of dsh-2 function, the posterior daughter cell was transformed into a second anterior-like cell. Loss of mom-5, a C. elegans frizzled homolog, produced a similar phenotype. We also show that the DSH-2 protein localized to the cell cortex in most cells of the embryo. In the absence of MOM-5/Fz, DSH-2 was localized to the cytoplasm, suggesting that MOM-5 regulates asymmetric cell division by controlling the localization of DSH-2. Although all neurons in C. elegans are produced by an invariant pattern of cell divisions, our results indicate that cell signaling may contribute to asymmetric neuroblast division during embryogenesis.

Identification of elements essential for transcription in Brugia malayi promoters

Little is known concerning promoter structure in the filarial parasites. Recently, transient transfection methods have been developed for the human filarial parasite Brugia malayi. These methods have been employed to localize the promoter for the 70kDa heat shock protein (BmHSP70) to a region extending 394nt upstream from the initiating codon of the BmHSP70 open reading frame. Replacement mutagenesis was used to define the elements necessary for BmHSP70 promoter activity in detail. Four domains, ranging in size from six to 22 nucleotides, were found to be necessary for full promoter activity. The two most distal domains encoded a binding site for the heat shock transcription factor and a putative binding site for the GAGA transcription factor, motifs that are found in many other HSP70 promoters. However, none of the essential domains contained sequences typical of cis elements that are usually found in the core domain of a eukaryotic promoter. The largest essential domain was located at positions -53 to -32, and included the splice leader addition site. These data suggest that the regulatory domains of the BmHSP70 promoter were similar to those found in other eukaryotes, but that the core promoter domain exhibited features that appeared to be distinct from those found in most other well-characterized eukaryotic promoters. An analysis of two additional promoters of B.malayi highly transcribed genes suggests that they also lack features commonly found in most eukaryotic core promoters, suggesting that the unique features of the BmHSP70 core promoter are not confined to this gene.