The homeodomain protein PAL-1 specifies a lineage-specific regulatory network in the C. elegans embryo

Landscape and regulation of mRNA translation in the early C. elegans embryo

Animal embryos rely on regulated translation of maternally deposited mRNAs to drive early development. Using low-input ribosome profiling combined with RNA sequencing on precisely staged embryos, we measured mRNA translation during the first four cell cycles of C. elegans development. We uncovered stage-specific patterns of developmentally coordinated translational regulation. We confirmed that mRNA localization correlates with translational eLiciency, though initial translational repression in germline precursors occurs before P-granule association. Our analysis suggests that the RNA-binding protein OMA-1 represses the translation of its target mRNAs in a stage-specific manner, while indirectly promoting the translational eLiciency of other transcripts. These findings illuminate how post-transcriptional mechanisms shape the embryonic proteome to direct cell diLerentiation, with implications for understanding similar regulation across species where maternal factors guide early development.

A lineage-resolved cartography of microRNA promoter activity in C. elegans empowers multidimensional developmental analysis

Elucidating the expression of microRNAs in developing single cells is critical for functional discovery. Here, we construct scCAMERA (single-cell cartography of microRNA expression based on reporter assay), utilizing promoter-driven fluorescent reporters in conjunction with imaging and lineage tracing. The cartography delineates the transcriptional activity of 54 conserved microRNAs in lineage-resolved single cells throughout C. elegans embryogenesis. The combinatorial expression of microRNAs partitions cells into fine clusters reflecting their function and anatomy. Notably, the expression of individual microRNAs exhibits high cell specificity and divergence among family members. Guided by cellular expression patterns, we identify developmental functions of specific microRNAs, including miR-1 in pharynx development and physiology, miR-232 in excretory canal morphogenesis by repressing NHR-25/NR5A, and a functional synergy between miR-232 and miR-234 in canal development, demonstrating the broad utility of scCAMERA. Furthermore, integrative analysis reveals that tissue-specific fate determinants activate microRNAs to repress protein production from leaky transcripts associated with alternative, especially neuronal, fates, thereby enhancing the fidelity of developmental fate differentiation. Collectively, our study offers rich opportunities for multidimensional expression-informed analysis of microRNA biology in metazoans.

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.

Investigation of PAL-1 requirement in C. elegans physiology using the auxin-inducible degradation system

The C. elegans PAL-1 protein encodes a caudal-like transcription factor that is required for posterior development and was recently implicated in stress response. We generated a transgenic strain of C. elegans with AID*::3xFLAG::wrmScarlet cassette knocked in at the C-terminal end of the pal-1 locus to enable an auxin-inducible degradation of PAL-1 . We found that auxin-induced degradation of PAL-1 starting from the L1 larval stage does not affect body length development but renders the animal sterile and shortens lifespan. This pal-1 ::AID*::3xFLAG::wrmScarlet strain will be a valuable resource for studying the requirement of PAL-1 in a temporal and tissue-specific manner.

CCR4-NOT subunit CCF-1/CNOT7 promotes transcriptional activation to multiple stress responses in Caenorhabditis elegans

CCR4-NOT is a versatile eukaryotic protein complex that controls multiple steps in gene expression regulation from synthesis to decay. In yeast, CCR4-NOT has been implicated in stress response regulation, though this function in other organisms remains unclear. In a genome-wide RNAi screen, we identified a subunit of the CCR4-NOT complex, ccf-1, as a requirement for the C. elegans transcriptional response to cadmium and acrylamide stress. Using whole-transcriptome RNA sequencing, we show that the knockdown of ccf-1 attenuates the activation of a broad range of stress-protective genes in response to cadmium and acrylamide, including those encoding heat shock proteins and xenobiotic detoxification. Consistently, survival assays show that the knockdown of ccf-1 decreases C. elegans stress resistance and normal lifespan. A yeast 2-hybrid screen using a CCF-1 bait identified the homeobox transcription factor PAL-1 as a physical interactor. Knockdown of pal-1 inhibits the activation of ccf-1 dependent stress genes and reduces C. elegans stress resistance. Gene expression analysis reveals that knockdown of ccf-1 and pal-1 attenuates the activation of elt-2 and elt-3 under stress that encode master transcriptional co-regulators of stress response in the C. elegans, and that overexpression of ELT-2 can suppress ccf-1's requirement for gene transcription in a stress-dependent manner. Our findings reveal a new role for CCR4-NOT in the environmental stress response and define its role in stress resistance and longevity in C. elegans.

The role of RNA-binding proteins in orchestrating germline development in Caenorhabditis elegans

RNA passed from parents to progeny controls several aspects of early development. The germline of the free-living nematode Caenorhabditis elegans contains many families of evolutionarily conserved RNA-binding proteins (RBPs) that target the untranslated regions of mRNA transcripts to regulate their translation and stability. In this review, we summarize what is known about the binding specificity of C. elegans germline RNA-binding proteins and the mechanisms of mRNA regulation that contribute to their function. We examine the emerging role of miRNAs in translational regulation of germline and embryo development. We also provide an overview of current technology that can be used to address the gaps in our understanding of RBP regulation of mRNAs. Finally, we present a hypothetical model wherein multiple 3'UTR-mediated regulatory processes contribute to pattern formation in the germline to ensure the proper and timely localization of germline proteins and thus a functional reproductive system.

Intricate Genetic Programs Controlling Dormancy in Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) displays the remarkable ability to transition in and out of dormancy, a hallmark of the pathogen’s capacity to evade the immune system and exploit susceptible individuals. Uncovering the gene regulatory programs that underlie the phenotypic shifts in MTB during disease latency and reactivation has posed a challenge. We develop an experimental system to precisely control dissolved oxygen levels in MTB cultures in order to capture the transcriptional events that unfold as MTB transitions into and out of hypoxia-induced dormancy. Using a comprehensive genome-wide transcription factor binding map and insights from network topology analysis, we identify regulatory circuits that deterministically drive sequential transitions across six transcriptionally and functionally distinct states encompassing more than three-fifths of the MTB genome. The architecture of the genetic programs explains the transcriptional dynamics underlying synchronous entry of cells into a dormant state that is primed to infect the host upon encountering favorable conditions.

Mycobacterium tuberculosis (MTB) persists within the host by counteracting disparate stressors including hypoxia. Peterson et al. report a transcriptional program that coordinates sequential state transitions to drive MTB in and out of hypoxia-induced dormancy. Among varied properties, this program encodes advanced preparedness to infect the host in favorable conditions.

Accessible chromatin reveals regulatory mechanisms underlying cell fate decisions during early embryogenesis

This study was conducted to investigate epigenetic landscape across multiple species and identify transcription factors (TFs) and their roles in controlling cell fate decision events during early embryogenesis. We made a comprehensively joint-research of chromatin accessibility of five species during embryogenesis by integration of ATAC-seq and RNA-seq datasets. Regulatory roles of candidate early embryonic TFs were investigated. Widespread accessible chromatin in early embryos overlapped with putative cis-regulatory sequences. Sets of cell-fate-determining TFs were identified. YOX1, a key cell cycle regulator, were found to homologous to clusters of TFs that are involved in neuron and epidermal cell-fate determination. Our research provides an intriguing insight into evolution of cell-fate decision during early embryogenesis among organisms.

Genome-wide RNAi screen for regulators of UPRmt in Caenorhabditis elegans mutants with defects in mitochondrial fusion

Mitochondrial dynamics plays an important role in mitochondrial quality control and the adaptation of metabolic activity in response to environmental changes. The disruption of mitochondrial dynamics has detrimental consequences for mitochondrial and cellular homeostasis and leads to the activation of the mitochondrial unfolded protein response (UPRmt), a quality control mechanism that adjusts cellular metabolism and restores homeostasis. To identify genes involved in the induction of UPRmt in response to a block in mitochondrial fusion, we performed a genome-wide RNAi screen in Caenorhabditis elegans mutants lacking the gene fzo-1, which encodes the ortholog of mammalian Mitofusin, and identified 299 suppressors and 86 enhancers. Approximately 90% of these 385 genes are conserved in humans, and one third of the conserved genes have been implicated in human disease. Furthermore, many have roles in developmental processes, which suggests that mitochondrial function and the response to stress are defined during development and maintained throughout life. Our dataset primarily contains mitochondrial enhancers and non-mitochondrial suppressors of UPRmt, indicating that the maintenance of mitochondrial homeostasis has evolved as a critical cellular function, which, when disrupted, can be compensated for by many different cellular processes. Analysis of the subsets 'non-mitochondrial enhancers' and 'mitochondrial suppressors' suggests that organellar contact sites, especially between the ER and mitochondria, are of importance for mitochondrial homeostasis. In addition, we identified several genes involved in IP3 signaling that modulate UPRmt in fzo-1 mutants and found a potential link between pre-mRNA splicing and UPRmt activation.

DeltaNeTS+: elucidating the mechanism of drugs and diseases using gene expression and transcriptional regulatory networks

BACKGROUND

Knowledge on the molecular targets of diseases and drugs is crucial for elucidating disease pathogenesis and mechanism of action of drugs, and for driving drug discovery and treatment formulation. In this regard, high-throughput gene transcriptional profiling has become a leading technology, generating whole-genome data on the transcriptional alterations caused by diseases or drug compounds. However, identifying direct gene targets, especially in the background of indirect (downstream) effects, based on differential gene expressions is difficult due to the complexity of gene regulatory network governing the gene transcriptional processes.

RESULTS

In this work, we developed a network analysis method, called DeltaNeTS+, for inferring direct gene targets of drugs and diseases from gene transcriptional profiles. DeltaNeTS+ uses a gene regulatory network model to identify direct perturbations to the transcription of genes using gene expression data. Importantly, DeltaNeTS+ is able to combine both steady-state and time-course expression profiles, as well as leverage information on the gene network structure. We demonstrated the power of DeltaNeTS+ in predicting gene targets using gene expression data in complex organisms, including Caenorhabditis elegans and human cell lines (T-cell and Calu-3). More specifically, in an application to time-course gene expression profiles of influenza A H1N1 (swine flu) and H5N1 (avian flu) infection, DeltaNeTS+ shed light on the key differences of dynamic cellular perturbations caused by the two influenza strains.

CONCLUSION

DeltaNeTS+ is a powerful network analysis tool for inferring gene targets from gene expression profiles. As demonstrated in the case studies, by incorporating available information on gene network structure, DeltaNeTS+ produces accurate predictions of direct gene targets from a small sample size (~ 10 s). Integrating static and dynamic expression data with transcriptional network structure extracted from genomic information, as enabled by DeltaNeTS+, is crucial toward personalized medicine, where treatments can be tailored to individual patients. DeltaNeTS+ can be freely downloaded from http://www.github.com/cabsel/deltanetsplus .

Dual RNAseq analyses at soma and germline levels reveal evolutionary innovations in the elephantiasis-agent Brugia malayi, and adaptation of its Wolbachia endosymbionts

Brugia malayi is a human filarial nematode responsible for elephantiasis, a debilitating condition that is part of a broader spectrum of diseases called filariasis, including lymphatic filariasis and river blindness. Almost all filarial nematode species infecting humans live in mutualism with Wolbachia endosymbionts, present in somatic hypodermal tissues but also in the female germline which ensures their vertical transmission to the nematode progeny. These α-proteobacteria potentially provision their host with essential metabolites and protect the parasite against the vertebrate immune response. In the absence of Wolbachia wBm, B. malayi females become sterile, and the filarial nematode lifespan is greatly reduced. In order to better comprehend this symbiosis, we investigated the adaptation of wBm to the host nematode soma and germline, and we characterized these cellular environments to highlight their specificities. Dual RNAseq experiments were performed at the tissue-specific and ovarian developmental stage levels, reaching the resolution of the germline mitotic proliferation and meiotic differentiation stages. We found that most wBm genes, including putative effectors, are not differentially regulated between infected tissues. However, two wBm genes involved in stress responses are upregulated in the hypodermal chords compared to the germline, indicating that this somatic tissue represents a harsh environment to which wBm have adapted. A comparison of the B. malayi and C. elegans germline transcriptomes reveals a poor conservation of genes involved in the production of oocytes, with the filarial germline proliferative zone relying on a majority of genes absent from C. elegans. The first orthology map of the B. malayi genome presented here, together with tissue-specific expression enrichment analyses, indicate that the early steps of oogenesis are a developmental process involving genes specific to filarial nematodes, that likely result from evolutionary innovations supporting the filarial parasitic lifestyle.

Stressful development: integrating endoderm development, stress, and longevity

In addition to performing digestion and nutrient absorption, the intestine serves as one of the first barriers to the external environment, crucial for protecting the host from environmental toxins, pathogenic invaders, and other stress inducers. The gene regulatory network (GRN) governing embryonic development of the endoderm and subsequent differentiation and maintenance of the intestine has been well-documented in C. elegans. A key regulatory input that initiates activation of the embryonic GRN for endoderm and mesoderm in this animal is the maternally provided SKN-1 transcription factor, an ortholog of the vertebrate Nrf1 and 2, which, like C. elegans SKN-1, perform conserved regulatory roles in mediating a variety of stress responses across metazoan phylogeny. Other key regulatory factors in early gut development also participate in stress response as well as in innate immunity and aging and longevity. In this review, we discuss the intersection between genetic nodes that mediate endoderm/intestine differentiation and regulation of stress and homeostasis. We also consider how direct signaling from the intestine to the germline, in some cases involving SKN-1, facilitates heritable epigenetic changes, allowing transmission of adaptive stress responses across multiple generations. These connections between regulation of endoderm/intestine development and stress response mechanisms suggest that varying selective pressure exerted on the stress response pathways may influence the architecture of the endoderm GRN, thereby leading to genetic and epigenetic variation in early embryonic GRN regulatory events.

Detecting Ordinal Subcascades

Ordinal classifier cascades are constrained by a hypothesised order of the semantic class labels of a dataset. This order determines the overall structure of the decision regions in feature space. Assuming the correct order on these class labels will allow a high generalisation performance, while an incorrect one will lead to diminished results. In this way ordinal classifier systems can facilitate explorative data analysis allowing to screen for potential candidate orders of the class labels. Previously, we have shown that screening is possible for total orders of all class labels. However, as datasets might comprise samples of ordinal as well as non-ordinal classes, the assumption of a total ordering might be not appropriate. An analysis of subsets of classes is required to detect such hidden ordinal substructures. In this work, we devise a novel screening procedure for exhaustive evaluations of all order permutations of all subsets of classes by bounding the number of enumerations we have to examine. Experiments with multi-class data from diverse applications revealed ordinal substructures that generate new and support known relations.

Molecular Effects of Silver Nanoparticles on Monogenean Parasites: Lessons from Caenorhabditis elegans

The mechanisms of action of silver nanoparticles (AgNPs) in monogenean parasites of the genus Cichlidogyrus were investigated through a microarray hybridization approach using genomic information from the nematode Caenorhabditis elegans. The effects of two concentrations of AgNPs were explored, low (6 µg/L Ag) and high (36 µg/L Ag). Microarray analysis revealed that both concentrations of AgNPs activated similar biological processes, although by different mechanisms. Expression profiles included genes involved in detoxification, neurotoxicity, modulation of cell signaling, reproduction, embryonic development, and tegument organization as the main biological processes dysregulated by AgNPs. Two important processes (DNA damage and cell death) were mostly activated in parasites exposed to the lower concentration of AgNPs. To our knowledge, this is the first study providing information on the sub-cellular and molecular effects of exposure to AgNPs in metazoan parasites of fish.

Multilevel regulation of muscle-specific transcription factor hlh-1 during Caenorhabditis elegans embryogenesis

hlh-1 is a myogenic transcription factor required for body-wall muscle specification during embryogenesis in Caenorhabditis elegans. Despite its well-known role in muscle specification, comprehensive regulatory control upstream of hlh-1 remains poorly defined. Here, we first established a statistical reference for the spatiotemporal expression of hlh-1 at single-cell resolution up to the second last round of divisions for most of the cell lineages (from 4- to 350-cell stage) using 13 wild-type embryos. We next generated lineal expression of hlh-1 after RNA interference (RNAi) perturbation of 65 genes, which were selected based on their degree of conservation, mutant phenotypes, and known roles in development. We then compared the expression profiles between wild-type and RNAi embryos by clustering according to their lineal expression patterns using mean-shift and density-based clustering algorithms, which not only confirmed the roles of existing genes but also uncovered the potential functions of novel genes in muscle specification at multiple levels, including cellular, lineal, and embryonic levels. By combining the public data on protein-protein interactions, protein-DNA interactions, and genetic interactions with our RNAi data, we inferred regulatory pathways upstream of hlh-1 that function globally or locally. This work not only revealed diverse and multilevel regulatory mechanisms coordinating muscle differentiation during C. elegans embryogenesis but also laid a foundation for further characterizing the regulatory pathways controlling muscle specification at the cellular, lineal (local), or embryonic (global) level.

Caudal-dependent cell positioning directs morphogenesis of the C. elegans ventral epidermis

Strikingly, epithelial morphogenesis remains incomplete at the end of C. elegans embryonic development; newly hatched larvae undergo extensive remodelling of their ventral epidermis during the first larval stage (L1), when newly-born epidermal cells move ventrally to complete the epidermal syncytium. Prior to this remodelling, undivided lateral seam cells produce anterior adherens junction processes that are inherited by the anterior daughter cells following an asymmetric division during L1. These adherens junction processes provide the ventral migratory route for these anterior daughters. Here, we show that these processes are perturbed in pal-1/caudal mutant animals, resulting in their inheritance by posterior, seam-fated daughters. This causes aberrant migration of seam daughter cells, disrupting the ventral epidermis. Using 4D-lineaging, we demonstrate that this larval epidermal morphogenesis defect in pal-1 mutants can be traced directly back to an initial cell positioning defect in the embryo. pal-1 expression, driven by a single intronic enhancer, is required to correctly position the seam cells in embryos such that the appropriate cell junctions support the correct migratory paths of seam daughters later in development, irrespective of their fate. Thus, during ventral epithelial remodelling in C. elegans, we show that the position of migrating cells, specified by pal-1/caudal, appears to be more important than their fate in driving morphogenesis.

•

caudal/pal-1 is required to form the correct cell junctions during embryogenesis.

•

Correctly placed cell junctions direct larval ventral epithelial cell migration.

•

larval epithelial cell migration occurs independently of cell fate.

•

Embryonic epidermal expression of pal-1 is dependent on a single intronic enhancer.

Assessing phenotype order in molecular data

Biological entities are key elements of biomedical research. Their definition and their relationships are important in areas such as phylogenetic reconstruction, developmental processes or tumor evolution. Hypotheses about relationships like phenotype order are often postulated based on prior knowledge or belief. Evidence on a molecular level is typically unknown and whether total orders are reflected in the molecular measurements is unclear or not assessed. In this work we propose a method that allows a fast and exhaustive screening for total orders in large datasets. We utilise ordinal classifier cascades to identify discriminable molecular representations of the phenotypes. These classifiers are constrained by an order hypothesis and are highly sensitive to incorrect assumptions. Two new error bounds, which are introduced and theoretically proven, lead to a substantial speed-up and allow the application to large collections of many phenotypes. In our experiments we show that by exhaustively evaluating all possible candidate orders, we are able to identify phenotype orders that best coincide with the high-dimensional molecular profiles.

A universal transcriptomic signature of age reveals the temporal scaling of Caenorhabditis elegans aging trajectories

We collected 60 age-dependent transcriptomes for C. elegans strains including four exceptionally long-lived mutants (mean adult lifespan extended 2.2- to 9.4-fold) and three examples of lifespan-increasing RNAi treatments. Principal Component Analysis (PCA) reveals aging as a transcriptomic drift along a single direction, consistent across the vastly diverse biological conditions and coinciding with the first principal component, a hallmark of the criticality of the underlying gene regulatory network. We therefore expected that the organism's aging state could be characterized by a single number closely related to vitality deficit or biological age. The "aging trajectory", i.e. the dependence of the biological age on chronological age, is then a universal stochastic function modulated by the network stiffness; a macroscopic parameter reflecting the network topology and associated with the rate of aging. To corroborate this view, we used publicly available datasets to define a transcriptomic biomarker of age and observed that the rescaling of age by lifespan simultaneously brings together aging trajectories of transcription and survival curves. In accordance with the theoretical prediction, the limiting mortality value at the plateau agrees closely with the mortality rate doubling exponent estimated at the cross-over age near the average lifespan. Finally, we used the transcriptomic signature of age to identify possible life-extending drug compounds and successfully tested a handful of the top-ranking molecules in C. elegans survival assays and achieved up to a +30% extension of mean lifespan.

Systems biology of embryonic development: Prospects for a complete understanding of the Caenorhabditis elegans embryo

The convergence of developmental biology and modern genomics tools brings the potential for a comprehensive understanding of developmental systems. This is especially true for the Caenorhabditis elegans embryo because its small size, invariant developmental lineage, and powerful genetic and genomic tools provide the prospect of a cellular resolution understanding of messenger RNA (mRNA) expression and regulation across the organism. We describe here how a systems biology framework might allow large-scale determination of the embryonic regulatory relationships encoded in the C. elegans genome. This framework consists of two broad steps: (a) defining the "parts list"-all genes expressed in all cells at each time during development and (b) iterative steps of computational modeling and refinement of these models by experimental perturbation. Substantial progress has been made towards defining the parts list through imaging methods such as large-scale green fluorescent protein (GFP) reporter analysis. Imaging results are now being augmented by high-resolution transcriptome methods such as single-cell RNA sequencing, and it is likely the complete expression patterns of all genes across the embryo will be known within the next few years. In contrast, the modeling and perturbation experiments performed so far have focused largely on individual cell types or genes, and improved methods will be needed to expand them to the full genome and organism. This emerging comprehensive map of embryonic expression and regulatory function will provide a powerful resource for developmental biologists, and would also allow scientists to ask questions not accessible without a comprehensive picture. This article is categorized under: Invertebrate Organogenesis > Worms Technologies > Analysis of the Transcriptome Gene Expression and Transcriptional Hierarchies > Gene Networks and Genomics.

Cell Fate Maintenance and Reprogramming During the Oocyte-to-Embryo Transition

This chapter reviews our current understanding of the mechanisms that regulate reprogramming during the oocyte-to-embryo transition (OET). There are two major events reshaping the transcriptome during OET. One is the clearance of maternal transcripts in the early embryo, extensively reviewed by others. The other event, which is the focus of this chapter, is the embryonic (or zygotic) genome activation (EGA). The mechanisms controlling EGA can be broadly divided into transcriptional and posttranscriptional. The former includes the regulation of the basal transcription machinery, the regulation by specific transcription factors and chromatin modifications. The latter is performed mostly via specific RNA-binding proteins (RBPs). Different animal models have been used to decipher the regulation of EGA. These models are often biased for the specific type of regulation, which is why we discuss the models ranging from invertebrates to mammals. Whether these biases stem from incomplete understanding of EGA in these models, or reflect evolutionarily distinct solutions to EGA regulation, is a key unresolved problem in developmental biology. As the mechanisms controlling developmental reprogramming can, and in some cases have been shown to, function in differentiated cells subjected to induced reprogramming, our understanding of EGA regulation may have implications for the efficiency of induced reprogramming and, thus, for regenerative medicine.

A Transcriptional Lineage of the Early C. elegans Embryo

During embryonic development, cells must establish fates, morphologies, and behaviors in coordination with one another to form a functional body. A prevalent hypothesis for how this coordination is achieved is that each cell's fate and behavior is determined by a defined mixture of RNAs. Only recently has it become possible to measure the full suite of transcripts in a single cell. Here we quantify genome-wide mRNA abundance in each cell of the Caenorhabditis elegans embryo up to the 16-cell stage. We describe spatially dynamic expression, quantify cell-specific differential activation of the zygotic genome, and identify genes that were previously unappreciated as being critical for development. We present an interactive data visualization tool that allows broad access to our dataset. This genome-wide single-cell map of mRNA abundance, alongside the well-studied life history and fate of each cell, describes at a cellular resolution the mRNA landscape that guides development.

Probing and rearranging the transcription factor network controlling the C. elegans endoderm

The ELT-2 GATA factor is the predominant transcription factor regulating gene expression in the C. elegans intestine, following endoderm specification. We comment on our previous study (Wiesenfahrt et al., 2016) that investigated how the elt-2 gene is controlled by END-1, END-3 and ELT-7, the 3 endoderm specific GATA factors that lie upstream in the regulatory hierarchy. We also discuss the unexpected result that ELT-2, if expressed sufficiently early and at sufficiently high levels, can specify the C. elegans endoderm, replacing the normal functions of END-1 and END-3.

Combinatorial decoding of the invariant C. elegans embryonic lineage in space and time

Understanding how a single cell, the zygote, can divide and differentiate to produce the diverse animal cell types is a central goal of developmental biology research. The model organism Caenorhabditis elegans provides a system that enables a truly comprehensive understanding of this process across all cells. Its invariant cell lineage makes it possible to identify all of the cells in each individual and compare them across organisms. Recently developed methods automate the process of cell identification, allowing high-throughput gene expression characterization and phenotyping at single cell resolution. In this Review, we summarize the sequences of events that pattern the lineage including establishment of founder cell identity, the signaling pathways that diversify embryonic fate, and the regulators involved in patterning within these founder lineages before cells adopt their terminal fates. We focus on insights that have emerged from automated approaches to lineage tracking, including insights into mechanisms of robustness, context-specific regulation of gene expression, and temporal coordination of differentiation. We suggest a model by which lineage history produces a combinatorial code of transcription factors that act, often redundantly, to ensure terminal fate.

Quantitative Differences in Nuclear β-catenin and TCF Pattern Embryonic Cells in C. elegans

The Wnt signaling pathway plays a conserved role during animal development in transcriptional regulation of distinct targets in different developmental contexts but it remains unclear whether quantitative differences in the nuclear localization of effector proteins TCF and β-catenin contribute to context-specific regulation. We investigated this question in Caenorhabditis elegans embryos by quantifying nuclear localization of fluorescently tagged SYS-1/β-catenin and POP-1/TCF and expression of Wnt ligands at cellular resolution by time-lapse microscopy and automated lineage tracing. We identified reproducible, quantitative differences that generate a subset of Wnt-signaled cells with a significantly higher nuclear concentration of the TCF/β-catenin activating complex. Specifically, β-catenin and TCF are preferentially enriched in nuclei of daughter cells whose parents also had high nuclear levels of that protein, a pattern that could influence developmental gene expression. Consistent with this, we found that expression of synthetic reporters of POP-1-dependent activation is biased towards cells that had high nuclear SYS-1 in consecutive divisions. We identified new genes whose embryonic expression patterns depend on pop-1. Most of these require POP-1 for either transcriptional activation or repression, and targets requiring POP-1 for activation are more likely to be expressed in the cells with high nuclear SYS-1 in consecutive divisions than those requiring POP-1 for repression. Taken together, these results indicate that SYS-1 and POP-1 levels are influenced by the parent cell’s SYS-1/POP-1 levels and this may provide an additional mechanism by which POP-1 regulates distinct targets in different developmental contexts.

The Wnt signaling pathway is active during the development of all multi-cellular animals and also improperly re-activated in many cancers. Here, we use time-lapse microscopy to quantify the nuclear localization of several proteins in response to Wnt signaling throughout early embryonic development in the nematode worm, C. elegans. We find that cells that received a Wnt signal in the previous division respond more strongly to a Wnt signal in the next division, in part by localizing more of the regulator β-catenin to the nucleus. This causes the relative enrichment of Wnt pathway proteins in the nuclei of repeatedly signaled cells, which we show likely impacts the activation of Wnt target genes. This represents a novel mechanism for the regulation of Wnt pathway targets in development and disease.

Exploring the roles of basal transcription factor 3 in eukaryotic growth and development

Basal transcription factor 3 (BTF3) has been reported to play a significant part in the transcriptional regulation linking with eukaryotes growth and development. Alteration in the BTF3 gene expression patterns or variation in their activities adds to the explanation of different signaling pathways and regulatory networks. Moreover, BTF3s often respond to numerous stresses, and subsequently they are involved in regulation of various mechanisms. BTF3 proteins also function through protein-protein contact, which can assist us to identify the multifaceted processes of signaling and transcriptional regulation controlled by BTF3 proteins. In this review, we discuss current advances made in starting to explore the roles of BTF3 transcription factors in eukaryotes especially in plant growth and development.

Repression of somatic cell fate in the germline

Germ cells must transmit genetic information across generations, and produce gametes while also maintaining the potential to form all cell types after fertilization. Preventing the activation of somatic programs is, therefore, crucial to the maintenance of germ cell identity. Studies in Caenorhabditis elegans, Drosophila melanogaster, and mouse have revealed both similarities and differences in how somatic gene expression is repressed in germ cells, thereby preventing their conversion into somatic tissues. This review will focus on recent developments in our understanding of how global or gene-specific transcriptional repression, chromatin regulation, and translational repression operate in the germline to maintain germ cell identity and repress somatic differentiation programs.

E3 ubiquitin ligases promote progression of differentiation during C. elegans embryogenesis

Regulated choice between cell fate maintenance and differentiation provides decision points in development to progress toward more restricted cell fates or to maintain the current one. Caenorhabditis elegans embryogenesis follows an invariant cell lineage where cell fate is generally more restricted upon each cell division. EMS is a progenitor cell in the four-cell embryo that gives rise to the endomesoderm. We recently found that when ubiquitin-mediated protein degradation is compromised, the anterior daughter of EMS, namely MS, reiterates the EMS fate. This observation demonstrates an essential function of ubiquitin-mediated protein degradation in driving the progression of EMS-to-MS differentiation. Here we report a genome-wide screen of the ubiquitin pathway and extensive lineage analyses. The results suggest a broad role of E3 ligases in driving differentiation progression. First, we identified three substrate-binding proteins for two Cullin-RING ubiquitin ligase (CRL) E3 complexes that promote the progression from the EMS fate to MS, namely LIN-23/β-TrCP and FBXB-3 for the CRL1/SCF complex and ZYG-11/ZYG-11B for the CRL2 complex. Genetic analyses suggest these E3 ligases function through a multifunctional protein OMA-1 and the endomesoderm lineage specifier SKN-1 to drive differentiation. Second, we found that depletion of components of the CRL1/SCF complex induces fate reiteration in all major founder cell lineages. These data suggest that regulated choice between self-renewal and differentiation is widespread during C. elegans embryogenesis as in organisms with regulative development, and ubiquitin-mediated protein degradation drives the choice towards differentiation. Finally, bioinformatic analysis of time series gene expression data showed that expression of E3 genes is transiently enriched during time windows of developmental stage transitions. Transcription factors show similar enrichment, but not other classes of regulatory genes. Based on these findings we propose that ubiquitin-mediated protein degradation, like many transcription factors, function broadly as regulators driving developmental progression during embryogenesis in C. elegans.

Statistical stage transition detection method for small sample gene expression time series data

In terms of their internal (genetic) and external (phenotypic) states, living cells are always changing at varying rates. Periods of stable or low rate of change are often called States, Stages, or Phases, whereas high-rate periods are called Transitions or Transients. While states and transitions are observed phenotypically, such as cell differentiation, cancer progression, for example, are related with gene expression levels. On the other hand, stages of gene expression are definable based on changes of expression levels. Analyzing relations between state changes of phenotypes and stage transitions of gene expression levels is a general approach to elucidate mechanisms of life phenomena. Herein, we propose an algorithm to detect stage transitions in a time series of expression levels of a gene by defining statistically optimal division points. The algorithm shows detecting ability for simulated datasets. An annotation based analysis on detecting results for a dataset of initial development of Caenorhabditis elegans agrees with that are presented in the literature.

A dual role of the Wnt signaling pathway during aging in Caenorhabditis elegans

Wnt signaling is a major and highly conserved developmental pathway that guides many important events during embryonic and larval development. In adulthood, misregulation of Wnt signaling has been implicated in tumorigenesis and various age-related diseases. These effects occur through highly complicated cell-to-cell interactions mediated by multiple Wnt-secreted proteins. While they share a high degree of sequence similarity, their function is highly diversified. Although the role of Wnt ligands during development is well studied, very little is known about the possible actions of Wnt signaling in natural aging. In this study, Caenorhabditis elegans serves, for the first time, as a model system to determine the role of Wnt ligands in aging. Caenorhabditis elegans has five Wnt proteins, mom-2, egl-20, lin-44, cwn-1, and cwn-2. We show that all five Wnt ligands are expressed and active past the development stages. The ligand mom-2/Wnt plays a major detrimental role in longevity, whereas the function of lin-44/Wnt is beneficial for long life. Interestingly, no evidence was found for Wnt signaling being involved in cellular or oxidative stress responses during aging. Our results suggest that Wnt signaling regulates aging-intrinsic genetic pathways, opening a new research direction on the role of Wnt signaling in aging and age-related diseases.

Collaborative regulation of development but independent control of metabolism by two epidermis-specific transcription factors in Caenorhabditis elegans

Cell fate specification is typically initiated by a master regulator, which is relayed by tissue-specific regulatory proteins (usually transcription factors) for further enforcement of cell identities, but how the factors are coordinated among each other to "finish up" the specification remains poorly understood. Caenorhabditis elegans epidermis specification is initiated by a master regulator, ELT-1, that activates its targets, NHR-25 and ELT-3, two epidermis-specific transcription factors that are important for development but not for initial specification of epidermis, thus providing a unique paradigm for illustrating how the tissue-specific regulatory proteins work together to enforce cell fate specification. Here we addressed the question through contrasting genome-wide in vivo binding targets between NHR-25 and ELT-3. We demonstrate that the two factors bind discrete but conserved DNA motifs, most of which remain in proximity, suggesting formation of a complex between the two. In agreement with this, gene ontology analysis of putative target genes suggested differential regulation of metabolism but coordinated control of epidermal development between the two factors, which is supported by quantitative analysis of expression of their specific or common targets in the presence or absence of either protein. Functional validation of a subset of the target genes showed both activating and inhibitory roles of NHR-25 and ELT-3 in regulating their targets. We further demonstrated differential control of specification of AB and C lineage-derived epidermis. The results allow us to assemble a comprehensive gene network underlying C. elegans epidermis development that is likely to be widely used across species and provides insights into how tissue-specific transcription factors coordinate with one another to enforce cell fate specification initiated by its master regulator.

DREM 2.0: Improved reconstruction of dynamic regulatory networks from time-series expression data

Background

Modeling dynamic regulatory networks is a major challenge since much of the protein-DNA interaction data available is static. The Dynamic Regulatory Events Miner (DREM) uses a Hidden Markov Model-based approach to integrate this static interaction data with time series gene expression leading to models that can determine when transcription factors (TFs) activate genes and what genes they regulate. DREM has been used successfully in diverse areas of biological research. However, several issues were not addressed by the original version.

Results

DREM 2.0 is a comprehensive software for reconstructing dynamic regulatory networks that supports interactive graphical or batch mode. With version 2.0 a set of new features that are unique in comparison with other softwares are introduced. First, we provide static interaction data for additional species. Second, DREM 2.0 now accepts continuous binding values and we added a new method to utilize TF expression levels when searching for dynamic models. Third, we added support for discriminative motif discovery, which is particularly powerful for species with limited experimental interaction data. Finally, we improved the visualization to support the new features. Combined, these changes improve the ability of DREM 2.0 to accurately recover dynamic regulatory networks and make it much easier to use it for analyzing such networks in several species with varying degrees of interaction information.

Conclusions

DREM 2.0 provides a unique framework for constructing and visualizing dynamic regulatory networks. DREM 2.0 can be downloaded from: www.sb.cs.cmu.edu/drem.

Uptake of extracellular double-stranded RNA by SID-2

Ingested dsRNAs trigger RNA interference (RNAi) in many invertebrates, including the nematode Caenorhabditis elegans. Here we show that the C. elegans apical intestinal membrane protein SID-2 is required in C. elegans for the import of ingested dsRNA and that, when expressed in Drosophila S2 cells, SID-2 enables the uptake of dsRNAs. SID-2-dependent dsRNA transport requires an acidic extracellular environment and is selective for dsRNAs with at least 50 base pairs. Through structure-function analysis, we identify several SID-2 regions required for this activity, including three extracellular, positively charged histidines. Finally, we find that SID-2-dependent transport is inhibited by drugs that interfere with vesicle transport. Therefore, we propose that environmental dsRNAs are imported from the acidic intestinal lumen by SID-2 via endocytosis and are released from internalized vesicles in a secondary step mediated by the dsRNA channel SID-1. Similar multistep mechanisms may underlie the widespread observations of environmental RNAi.

Network inference of pal-1 lineage-specific regulation in the C. elegans embryo by structural equation modeling

The elucidation of spatial and temporal control during developmental stages is one of the central tasks for systems biology, and a variety of intracellular factors are known as regulators for specific gene expression. The activity information of those various factors is not directly reflected in their gene expression profiles. Hence, a method based on Structural Equation Modeling (SEM) is described. SEM can include the latent variables within the constructed model and infer the relationships among latent and observed variables, as a network model. An improved SEM approach for the construction of an optimal model is applied to infer the regulatory network for the determination of C lineage fate in C. elegans development. The inferred network model shows that the 13 analysed transcription factor genes were regulated by several other factors in addition to pal-1 expression. The other regulatory factors are those involved in protein accumulation and localization as important regulatory factors for normal development. Those regulatory factors were regulated sequentially in the network model. The regulation of the known pal-1 regulated genes was dependent on this sequential control of the regulatory factors. The interpretation of the network model shows insights to the complex regulation occurring during the C lineage determination by pal-1.

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.

Transcription factor redundancy and tissue-specific regulation: evidence from functional and physical network connectivity

Two major transcriptional regulators of Caenorhabditis elegans bodywall muscle (BWM) differentiation, hlh-1 and unc-120, are expressed in muscle where they are known to bind and regulate several well-studied muscle-specific genes. Simultaneously mutating both factors profoundly inhibits formation of contractile BWM. These observations were consistent with a simple network model in which the muscle regulatory factors drive tissue-specific transcription by binding selectively near muscle-specific targets to activate them. We tested this model by measuring the number, identity, and tissue-specificity of functional regulatory targets for each factor. Some joint regulatory targets (218) are BWM-specific and enriched for nearby HLH-1 binding. However, contrary to the simple model, the majority of genes regulated by one or both muscle factors are also expressed significantly in non-BWM tissues. We also mapped global factor occupancy by HLH-1, and created a genetic interaction map that identifies hlh-1 collaborating transcription factors. HLH-1 binding did not predict proximate regulatory action overall, despite enrichment for binding among BWM-specific positive regulatory targets of hlh-1. We conclude that these tissue-specific factors contribute much more broadly to the transcriptional output of muscle tissue than previously thought, offering a partial explanation for widespread HLH-1 occupancy. We also identify a novel regulatory connection between the BWM-specific hlh-1 network and the hlh-8/twist nonstriated muscle network. Finally, our results suggest a molecular basis for synthetic lethality in which hlh-1 and unc-120 mutant phenotypes are mutually buffered by joint additive regulation of essential target genes, with additional buffering suggested via newly identified hlh-1 interacting factors.

The Caenorhabditis elegans epidermis as a model skin. I: development, patterning, and growth

The skin of the nematode Caenorhabditis elegans is composed of a simple epidermal epithelium and overlying cuticle. The skin encloses the animal and plays central roles in body morphology and physiology; its simplicity and accessibility make it a tractable genetic model for several aspects of skin biology. Epidermal precursors are specified by a hierarchy of transcriptional regulators. Epidermal cells form on the dorsal surface of the embryo and differentiate to form the epidermal primordium, which then spreads out in a process of epiboly to enclose internal tissues. Subsequent elongation of the embryo into a vermiform larva is driven by cell shape changes and cell fusions in the epidermis. Most epidermal cells fuse in mid-embryogenesis to form a small number of multinucleate syncytia. During mid-embryogenesis the epidermis also becomes intimately associated with underlying muscles, performing a tendon-like role in transmitting muscle force. Post-embryonic development of the epidermis involves growth by addition of new cells to the syncytia from stem cell-like epidermal seam cells and by an increase in cell size driven by endoreplication of the chromosomes in epidermal nuclei.

Multidimensional regulation of gene expression in the C. elegans embryo

How cells adopt different expression patterns is a fundamental question of developmental biology. We quantitatively measured reporter expression of 127 genes, primarily transcription factors, in every cell and with high temporal resolution in C. elegans embryos. Embryonic cells are highly distinct in their gene expression; expression of the 127 genes studied here can distinguish nearly all pairs of cells, even between cells of the same tissue type. We observed recurrent lineage-regulated expression patterns for many genes in diverse contexts. These patterns are regulated in part by the TCF-LEF transcription factor POP-1. Other genes' reporters exhibited patterns correlated with tissue, position, and left–right asymmetry. Sequential patterns both within tissues and series of sublineages suggest regulatory pathways. Expression patterns often differ between embryonic and larval stages for the same genes, emphasizing the importance of profiling expression in different stages. This work greatly expands the number of genes in each of these categories and provides the first large-scale, digitally based, cellular resolution compendium of gene expression dynamics in live animals. The resulting data sets will be a useful resource for future research.

General trends in the utilization of structural factors contributing to biological complexity

During evolution, proteins containing newly emerged domains and the increasing proportion of multidomain proteins in the full genome-encoded proteome (GEP) have substantially contributed to increasing biological complexity. However, it is not known how these two potential structural factors are preferentially utilized at given physiological states. Here, we classified proteins according to domain number and domain age and explored the general trends across species for the utilization of proteins from GEP to various certain-state proteomes (CSPs, i.e., all the proteins expressed at certain physiological states). We found that multidomain proteins or only older domain-containing proteins are significantly overrepresented in CSPs compared with GEP, which is a trend that is stronger in multicellular organisms than in unicellular organisms. Interestingly, the strengths of overrepresentation decreased during evolution of multicellular eukaryotes. When comparing across CSPs, we found that multidomain proteins are more overrepresented in complex tissues than in simpler ones, whereas no difference among proteins with domains of different ages is evident between complex and simple tissues. Thus, biological complexity under certain conditions is more significantly realized by diverse domain organization than by the emergence of new types of domain. In addition, we found that multidomain or only older domain-containing proteins tend to evolve slowly and generally are under stronger purifying selection, which may partly result from their general overrepresentation trends in CSPs.

Networks in Caenorhabditis elegans

The network paradigm has become a pervasive theme in biology over the last decade, as increasingly large functional genomic datasets are being collected to interrogate regulatory influences, physical interactions, and genetic dependencies between genes, transcripts, and proteins. These 'molecular interaction' networks can be analyzed collectively and individually to define their global architecture and local patterns of connectivity. These structural features ultimately underlie functional properties such as robustness, modularity, component circuitry (e.g. feedback loops), dynamics, and responses to perturbations. This review focuses on recent progress in elucidating molecular interaction networks using different kinds of functional assays in the classical genetic model for animal development, the roundworm Caenorhabditis elegans, with representative examples to illustrate current directions in different areas of network biology.

Myogenic conversion and transcriptional profiling of embryonic blastomeres in Caenorhabditis elegans

Myogenesis has proven to be a powerful paradigm for understanding cell fate specification and differentiation in many model organisms. This includes the nematode Caenorhabditis elegans for which the genetic, cellular, and molecular tools have allowed an in-depth understanding of muscle development. One tool not yet available in C. elegans is a robust, pure and prolific cell culture system to study myogenesis. As an alternative, this chapter describes a method by which the cell fates of early, uncommitted blastomeres in the embryo are converted to a myogenic lineage. This technique permits the nearly synchronous induction of myogenesis in vivo with the potential to generate a nearly homogeneous population of cells. Coupled with the RNA isolation and cDNA amplification methods that are also described, one can now profile gene expression throughout myogenesis using any platform of choice (e.g. expression arrays, next generation sequencing). Although limited by the artificial nature of this developing mass of muscle inside the eggshell, blastomere conversion and transcriptional profiling is a very powerful tool to investigate changes in gene expression associated with myogenesis in C. elegans that is applicable to many different cell types. When coupled with next generation sequencing, the method has the potential to yield a very high-resolution map of changes in gene expression throughout myogenesis.

The branched actin nucleator Arp2/3 promotes nuclear migrations and cell polarity in the C. elegans zygote

Regulated movements of the nucleus are essential during zygote formation, cell migrations, and differentiation of neurons. The nucleus moves along microtubules (MTs) and is repositioned on F-actin at the cellular cortex. Two families of nuclear envelope proteins, SUN and KASH, link the nucleus to the actin and MT cytoskeletons during nuclear movements. However, the role of actin nucleators in nuclear migration and positioning is poorly understood. We show that the branched actin nucleator, Arp2/3, affects nuclear movements throughout embryonic and larval development in C. elegans, including nuclear migrations in epidermal cells and neuronal precursors. In one-cell embryos the migration of the male pronucleus to meet the female pronucleus after fertilization requires Arp2/3. Loss of Arp2/3 or its activators changes the dynamics of non-muscle myosin, NMY-2, and alters the cortical accumulation of posterior PAR proteins. Reduced establishment of the posterior microtubule cytoskeleton in Arp2/3 mutants correlates with reduced male pronuclear migration. The UNC-84/SUN nuclear envelope protein that links the nucleus to the MT and actin cytoskeleton is known to regulate later nuclear migrations. We show here it also positions the male pronucleus. These studies demonstrate a global role for Arp2/3 in nuclear migrations. In the C. elegans one-cell embryo Arp2/3 promotes the establishment of anterior/posterior polarity and promotes MT growth that propels the anterior migration of the male pronucleus. In contrast with previous studies emphasizing pulling forces on the male pronucleus, we propose that robust MT nucleation pushes the male pronucleus anteriorly to join the female pronucleus.

Overcoming redundancy: an RNAi enhancer screen for morphogenesis genes in Caenorhabditis elegans

Morphogenesis is an important component of animal development. Genetic redundancy has been proposed to be common among morphogenesis genes, posing a challenge to the genetic dissection of morphogenesis mechanisms. Genetic redundancy is more generally a challenge in biology, as large proportions of the genes in diverse organisms have no apparent loss of function phenotypes. Here, we present a screen designed to uncover redundant and partially redundant genes that function in an example of morphogenesis, gastrulation in Caenorhabditis elegans. We performed an RNA interference (RNAi) enhancer screen in a gastrulation-sensitized double-mutant background, targeting genes likely to be expressed in gastrulating cells or their neighbors. Secondary screening identified 16 new genes whose functions contribute to normal gastrulation in a nonsensitized background. We observed that for most new genes found, the closest known homologs were multiple other C. elegans genes, suggesting that some may have derived from rounds of recent gene duplication events. We predict that such genes are more likely than single copy genes to comprise redundant or partially redundant gene families. We explored this prediction for one gene that we identified and confirmed that this gene and five close relatives, which encode predicted substrate recognition subunits (SRSs) for a CUL-2 ubiquitin ligase, do indeed function partially redundantly with each other in gastrulation. Our results implicate new genes in C. elegans gastrulation, and they show that an RNAi-based enhancer screen in C. elegans can be used as an efficient means to identify important but redundant or partially redundant developmental genes.

Functional modularity of nuclear hormone receptors in a Caenorhabditis elegans metabolic gene regulatory network

We present the first gene regulatory network (GRN) that pertains to post-developmental gene expression. Specifically, we mapped a transcription regulatory network of Caenorhabditis elegans metabolic gene promoters using gene-centered yeast one-hybrid assays. We found that the metabolic GRN is enriched for nuclear hormone receptors (NHRs) compared with other gene-centered regulatory networks, and that these NHRs organize into functional network modules.

The NHR family has greatly expanded in nematodes; C. elegans has 284 NHRs, whereas humans have only 48. We show that the NHRs in the metabolic GRN have metabolic phenotypes, suggesting that they do not simply function redundantly.

The mediator subunit MDT-15 preferentially interacts with NHRs that occur in the metabolic GRN.

We describe an NHR circuit that responds to nutrient availability and propose a model for the evolution and organization of NHRs in C. elegans metabolic regulatory networks.

Physical and/or regulatory interactions between transcription factors (TFs) and their target genes are essential to establish body plans of multicellular organisms during development, and these interactions have been studied extensively in the context of GRNs. The precise control of differential gene expression is also of critical importance to maintain physiological homeostasis, and many metabolic disorders such as obesity and diabetes coincide with substantial changes in gene expression. Much work has focused on the GRNs that control metazoan development; however, the design principles and organization of the GRNs that control systems physiology remain largely unexplored.

In this study, we present the first gene-centered GRN that includes ∼70 genes involved in C. elegans metabolism and physiology, 100 TFs and more than 500 protein–DNA interactions between them. The resulting metabolic GRN is enriched for NHRs, compared with other gene-centered regulatory networks. NHRs are well-known regulators of lipid meta-qj;bolism in mammals. The transcriptional activity of NHRs can be modified by diffusible ligands, which allows these TFs to function as molecular sensors and rapidly alter the expression of their target genes. Interestingly, NHRs comprise the largest family of TFs in nematodes; the C. elegans genome encodes 284 NHRs, most of which are uncharacterized. Furthermore, their organization in GRNs has not yet been investigated. In our study, we show that the C. elegans NHRs that we retrieved in the metabolic GRN organize into network modules, and that most of these NHRs function to maintain lipid homeostasis in the nematode. Interestingly, network modularity has been proposed to facilitate rapid and robust changes in gene expression. Our results suggest that the C. elegans metabolic GRN may have evolved by combining NHR family expansion with the specific modular wiring of NHRs to enable the rapid adaptation of the animal to different environmental cues.

NHRs can interact with transcriptional cofactors such as chromatin remodeling complexes and Mediator components. For instance, the C. elegans Mediator subunit, MDT-15, can interact with NHR-49 to regulate the expression of its target genes. To find all the TFs that MDT-15 can interact with, we performed systematic yeast two-hybrid assays with MDT-15 versus 755 full-length TFs. We found that MDT-15 preferentially associates with NHRs, and specifically with those NHRs that confer a metabolic phenotype and that occur in the metabolic GRN. This illustrates the central role of MDT-15 in the regulation of metabolic gene expression.

Using a variety of genetic and biochemical approaches, we characterized NHR-86 in more detail. NHR-86 participates in one of the two NHR modules, and has a high-flux capacity; that is it has both a high incoming and a high outgoing degree. We obtained an nhr-86 mutant and generated an NHR-86 antibody, and showed that NHR-86 functions as an auto-repressor in vivo and that nhr-86 mutant animals store abnormally high levels of body fat.

Finally, we discovered a novel NHR circuit that responds to nutrient availability. In this circuit NHR-45 regulates the activity of nhr-178 promoter in two distinct physiologically important tissues: the intestine and the hypodermis. Both of these NHRs are required to maintain lipid homeostasis in C. elegans. The expression of nhr-178 is responsive to the nutritional status of the animal, which switches between ON and OFF states in the hypodermis. We found that NHR-45 activity is necessary to control this switch in the hypodermis. Interestingly, NHR-45 has opposite effects on the activity of the nhr-178 promoter in these tissues: NHR-45 activates this promoter in the intestine, but represses it in the hypodermis.

Altogether our study leads to a model in which the expansion of the NHR family, TFs that have the capacity to act as fast molecular sensors, is combined with a modular network organization to enable rapid and robust responses to various environmental cues.

Gene regulatory networks (GRNs) provide insights into the mechanisms of differential gene expression at a systems level. GRNs that relate to metazoan development have been studied extensively. However, little is still known about the design principles, organization and functionality of GRNs that control physiological processes such as metabolism, homeostasis and responses to environmental cues. In this study, we report the first experimentally mapped metazoan GRN of Caenorhabditis elegans metabolic genes. This network is enriched for nuclear hormone receptors (NHRs). The NHR family has greatly expanded in nematodes: humans have 48 NHRs, but C. elegans has 284, most of which are uncharacterized. We find that the C. elegans metabolic GRN is highly modular and that two GRN modules predominantly consist of NHRs. Network modularity has been proposed to facilitate a rapid response to different cues. As NHRs are metabolic sensors that are poised to respond to ligands, this suggests that C. elegans GRNs evolved to enable rapid and adaptive responses to different cues by a concurrence of NHR family expansion and modular GRN wiring.

Cell fate specification in the C. elegans embryo

Cell specification requires that particular subsets of cells adopt unique expression patterns that ultimately define the fates of their descendants. In C. elegans, cell fate specification involves the combinatorial action of multiple signals that produce activation of a small number of "blastomere specification" factors. These initiate expression of gene regulatory networks that drive development forward, leading to activation of "tissue specification" factors. In this review, the C. elegans embryo is considered as a model system for studies of cell specification. The techniques used to study cell fate in this species, and the themes that have emerged, are described.

Revealing developmental networks by comparative transcriptomics

Metazoan development relies upon the precise control of the genome's expression. This enables different cells in the animal to have different properties, despite having the same genetic material, and different animals to have different morphologies despite sharing developmental genes. However, near-identical organisms may have different overall genomic content, suggesting that the mechanisms by which evolution of the phenotype proceeds on a global level are not well understood. We review here recent works that have discovered a tremendous amount of variation between the developmental transcriptomes of both closely and distantly related organisms. It is evident that the evolution of regulatory programs occurs at a rapid rate comparable to that of other genomic processes. Distinguishing the selective pressures on each regulatory element will thus be crucial towards understanding its functional relevance. We propose that such a comparative approach is most suited to the identification of unifying principles in cell fate specification and differentiation in the animal embryo.

Asymmetric Wolbachia segregation during early Brugia malayi embryogenesis determines its distribution in adult host tissues

Wolbachia are required for filarial nematode survival and fertility and contribute to the immune responses associated with human filarial diseases. Here we developed whole-mount immunofluorescence techniques to characterize Wolbachia somatic and germline transmission patterns and tissue distribution in Brugia malayi, a nematode responsible for lymphatic filariasis. In the initial embryonic divisions, Wolbachia segregate asymmetrically such that they occupy only a small subset of cells in the developing embryo, facilitating their concentration in the adult hypodermal chords and female germline. Wolbachia are not found in male reproductive tissues and the absence of Wolbachia from embryonic germline precursors in half of the embryos indicates Wolbachia loss from the male germline may occur in early embryogenesis. Wolbachia rely on fusion of hypodermal cells to populate adult chords. Finally, we detect Wolbachia in the secretory canal lumen suggesting living worms may release bacteria and/or their products into their host.

Filarial diseases affect over 150 million people in tropical countries. They are caused by parasitic nematodes like Brugia malayi that rely on their endosymbiont Wolbachia for their survival and fertility. These bacteria are a recognized drug target in the search for treatments killing adult worms. To understand the transmission of Wolbachia from the embryonic to adult stages, we developed new techniques to track these bacteria at the cellular and tissue levels. These techniques include immunofluorescence in whole mount adult tissues and embryos. We found that Wolbachia segregate asymetrically in specific cells, in a lineage-specific manner during early Brugia embryogenesis, and rely on cell fusion to subsequently populate the adult hypodermal chords. From the chords, the Wolbachia can be secreted in the secretory-excretory canal, suggesting that in addition to dead worms releasing the bacteria in the human body, living worms may also secrete Wolbachia, whose role in stimulating the immune system in filarial pathologies is now well established.

A switch from a gradient to a threshold mode in the regulation of a transcriptional cascade promotes robust execution of meiosis in budding yeast

Tight regulation of developmental pathways is of critical importance to all organisms, and is achieved by a transcriptional cascade ensuring the coordinated expression of sets of genes. We aimed to explore whether a strong signal is required to enter and complete a developmental pathway, by using meiosis in budding yeast as a model. We demonstrate that meiosis in budding yeast is insensitive to drastic changes in the levels of its consecutive positive regulators (Ime1, Ime2, and Ndt80). Entry into DNA replication is not correlated with the time of transcription of the early genes that regulate this event. Entry into nuclear division is directly regulated by the time of transcription of the middle genes, as premature transcription of their activator NDT80, leads to a premature entry into the first meiotic division, and loss of coordination between DNA replication and nuclear division. We demonstrate that Cdk1/Cln3 functions as a negative regulator of Ime2, and that ectopic expression of Cln3 delays entry into nuclear division as well as NDT80 transcription. Because Ime2 functions as a positive regulator for premeiotic DNA replication and NDT80 transcription, as well as a negative regulator of Cdk/Cln, we suggest that a double negative feedback loop between Ime2 and Cdk1/Cln3 promotes a bistable switch from the cell cycle to meiosis. Moreover, our results suggest a regulatory mode switch that ensures robust meiosis as the transcription of the early meiosis-specific genes responds in a graded mode to Ime1 levels, whereas that of the middle and late genes as well as initiation of DNA replication, are regulated in a threshold mode.

The spatial dynamics of tissue-specific promoters during C. elegans development

To understand whether the spatial organization of the genome reflects the cell's differentiated state, we examined whether genes assume specific subnuclear positions during Caenorhabditis elegans development. Monitoring the radial position of developmentally controlled promoters in embryos and larval tissues, we found that small integrated arrays bearing three different tissue-specific promoters have no preferential position in nuclei of undifferentiated embryos. However, in differentiated cells, they shifted stably toward the nuclear lumen when activated, or to the nuclear envelope when silent. In contrast, large integrated arrays bearing the same promoters became heterochromatic and nuclear envelope-bound in embryos. Tissue-specific activation of promoters in these large arrays in larvae overrode the perinuclear anchorage. For transgenes that carry both active and inactive promoters, the inward shift of the active promoter was dominant. Finally, induction of master regulator HLH-1 prematurely induced internalization of a muscle-specific promoter array in embryos. Fluorescence in situ hybridization confirmed analogous results for the endogenous endoderm-determining gene pha-4. We propose that, in differentiated cells, subnuclear organization arises from the selective positioning of active and inactive developmentally regulated promoters. We characterize two forces that lead to tissue-specific subnuclear organization of the worm genome: large repeat-induced heterochromatin, which associates with the nuclear envelope like repressed genes in differentiated cells, and tissue-specific promoters that shift inward in a dominant fashion over silent promoters, when they are activated.