A regulatory network of T-box genes and the even-skipped homologue vab-7 controls patterning and morphogenesis in C. elegans

Spatiotemporal analysis of mRNA-protein relationships enhances transcriptome-based developmental inference

Elucidating the complex relationships between mRNA and protein expression at high spatiotemporal resolution is critical for unraveling multilevel gene regulation and enhancing mRNA-based developmental analyses. In this study, we conduct a single-cell analysis of mRNA and protein expression of transcription factors throughout C. elegans embryogenesis. Initially, cellular co-presence of mRNA and protein is low, increasing to a medium-high level (73%) upon factoring in delayed protein synthesis and long-term protein persistence. These factors substantially affect mRNA-protein concordance, leading to potential inaccuracies in mRNA-reliant gene detection and specificity characterization. Building on the learned relationship, we infer protein presence from mRNA expression and demonstrate its utility in identifying tissue-specific genes and elucidating relationships between genes and cells. This approach facilitates identifying the role of sptf-1/SP7 in neuronal lineage development. Collectively, this study provides insights into gene expression dynamics during rapid embryogenesis and approaches for improving the efficacy of transcriptome-based developmental analyses.

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.

A mosaic of conserved and novel modes of gene expression and morphogenesis in mesoderm and muscle formation of a larval bivalve

The mesoderm gives rise to several key morphological features of bilaterian animals including endoskeletal elements and the musculature. A number of regulatory genes involved in mesoderm and/or muscle formation (e.g., Brachyury (Bra), even-skipped (eve), Mox, myosin II heavy chain (mhc)) have been identified chiefly from chordates and the ecdysozoans Drosophila and Caenorhabditis elegans, but data for non-model protostomes, especially those belonging to the ecdysozoan sister clade, Lophotrochozoa (e.g., flatworms, annelids, mollusks), are only beginning to emerge. Within the lophotrochozoans, Mollusca constitutes the most speciose and diverse phylum. Interestingly, however, information on the morphological and molecular underpinnings of key ontogenetic processes such as mesoderm formation and myogenesis remains scarce even for prominent molluscan sublineages such as the bivalves. Here, we investigated myogenesis and developmental expression of Bra, eve, Mox, and mhc in the quagga mussel Dreissena rostriformis, an invasive freshwater bivalve and an emerging model in invertebrate evodevo. We found that all four genes are expressed during mesoderm formation, but some show additional, individual sites of expression during ontogeny. While Mox and mhc are involved in early myogenesis, eve is also expressed in the embryonic shell field and Bra is additionally present in the foregut. Comparative analysis suggests that Mox has an ancestral role in mesoderm and possibly muscle formation in bilaterians, while Bra and eve are conserved regulators of mesoderm development of nephrozoans (protostomes and deuterostomes). The fully developed Dreissena veliger larva shows a highly complex muscular architecture, supporting a muscular ground pattern of autobranch bivalve larvae that includes at least a velum muscle ring, three or four pairs of velum retractors, one or two pairs of larval retractors, two pairs of foot retractors, a pedal plexus, possibly two pairs of mantle retractors, and the muscles of the pallial line, as well as an anterior and a posterior adductor. As is typical for their molluscan kin, remodelling and loss of prominent larval features such as the velum musculature and various retractor systems appear to be also common in bivalves.

Supplementary information

The online version contains supplementary material available at 10.1007/s13127-022-00569-5.

A 4D single-cell protein atlas of transcription factors delineates spatiotemporal patterning during embryogenesis

Complex biological processes such as embryogenesis require precise coordination of cell differentiation programs across both space and time. Using protein-fusion fluorescent reporters and four-dimensional live imaging, we present a protein expression atlas of transcription factors (TFs) mapped onto developmental cell lineages during Caenorhabditis elegans embryogenesis, at single-cell resolution. This atlas reveals a spatiotemporal combinatorial code of TF expression, and a cascade of lineage-specific, tissue-specific and time-specific TFs that specify developmental states. The atlas uncovers regulators of embryogenesis, including an unexpected role of a skin specifier in neurogenesis and the critical function of an uncharacterized TF in convergent muscle differentiation. At the systems level, the atlas provides an opportunity to model cell state-fate relationships, revealing a lineage-dependent state diversity within functionally related cells and a winding trajectory of developmental state progression. Collectively, this single-cell protein atlas represents a valuable resource for elucidating metazoan embryogenesis at the molecular and systems levels.

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.

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.

Embryonic expression of priapulid Wnt genes

Posterior elongation of the developing embryo is a common feature of animal development. One group of genes that is involved in posterior elongation is represented by the Wnt genes, secreted glycoprotein ligands that signal to specific receptors on neighbouring cells and thereby establish cell-to-cell communication. In segmented animals such as annelids and arthropods, Wnt signalling is also likely involved in segment border formation and regionalisation of the segments. Priapulids represent unsegmented worms that are distantly related to arthropods. Despite their interesting phylogenetic position and their importance for the understanding of ecdysozoan evolution, priapulids still represent a highly underinvestigated group of animals. Here, we study the embryonic expression patterns of the complete sets of Wnt genes in the priapulids Priapulus caudatus and Halicryptus spinulosus. We find that both priapulids possess a complete set of 12 Wnt genes. At least in Priapulus, most of these genes are expressed in and around the posterior-located blastopore and thus likely play a role in posterior elongation. Together with previous work on the expression of other genetic factors such as caudal and even-skipped, this suggests that posterior elongation in priapulids is under control of the same (or very similar) conserved gene regulatory network as in arthropods.

Non-neuronal cell outgrowth in C. elegans

Cell outgrowth is a hallmark of some non-migratory developing cells during morphogenesis. Understanding the mechanisms that control cell outgrowth not only increases our knowledge of tissue and organ development, but can also shed light on disease pathologies that exhibit outgrowth-like behavior. C. elegans is a highly useful model for the analysis of genes and the function of their respective proteins. In addition, C. elegans also has several cells and tissues that undergo outgrowth during development. Here we discuss the outgrowth mechanisms of nine different C. elegans cells and tissues. We specifically focus on how these cells and tissues grow outward and the interactions they make with their environment. Through our own identification, and a meta-analysis, we also identify gene families involved in multiple cell outgrowth processes, which defined potential C. elegans core components of cell outgrowth, as well as identify a potential stepwise cell behavioral cascade used by cells undergoing outgrowth.

The Remarkably Diverse Family of T-Box Factors in Caenorhabditis elegans

The nematode Caenorhabditis elegans is a simple metazoan animal that is widely used as a model to understand the genetic control of development. The completely sequenced C. elegans genome contains 22 T-box genes, and they encode factors that show remarkable diversity in sequence, DNA-binding specificity, and function. Only three of the C. elegans T-box factors can be grouped into the conserved subfamilies found in other organisms, while the remaining factors are significantly diverged and unlike those in most other animals. While some of the C. elegans factors can bind canonical T-box binding elements, others bind and regulate target gene expression through distinct sequences. The nine genetically characterized T-box factors have varied functions in development and morphogenesis of muscle, hypodermal tissues, and neurons, as well as in early blastomere fate specification, cell migration, apoptosis, and sex determination, but the functions of most of the C. elegans T-box factors have not yet been extensively characterized. Like T-box factors in other animals, interaction with a Groucho-family corepressor and posttranslational SUMOylation have been shown to affect C. elegans T-box factor activity, and it is likely that additional mechanisms affecting T-box factor activity will be discovered using the effective genetic approaches in this organism.

Function of the C. elegans T-box factor TBX-2 depends on interaction with the UNC-37/Groucho corepressor

T-box transcription factors are important regulators of development in all animals, and altered expression of T-box factors has been identified in an increasing number of diseases and cancers. Despite these important roles, the mechanism of T-box factor activity is not well understood. We have previously shown that the Caenorhabditis elegans Tbx2 subfamily member TBX-2 functions as a transcriptional repressor to specify ABa-derived pharyngeal muscle, and that this function depends on SUMOylation. Here we show that TBX-2 function also depends on interaction with the Groucho-family corepressor UNC-37. TBX-2 interacts with UNC-37 in yeast two-hybrid assays via a highly conserved engrailed homology 1 (eh1) motif located near the TBX-2 C-terminus. Reducing unc-37 phenocopies tbx-2 mutants, resulting in a specific loss of anterior ABa-derived pharyngeal muscles and derepression of the tbx-2 promoter. Moreover, double mutants containing hypomorphic alleles of unc-37 and tbx-2 exhibit enhanced phenotypes, providing strong genetic evidence that unc-37 and tbx-2 share common functions in vivo. To test whether interaction with UNC-37 is necessary for TBX-2 activity, we developed a transgene rescue assay using a tbx-2 containing fosmid and found that mutating the tbx-2 eh1 motif reduced rescue of a tbx-2 null mutant. These results indicate that TBX-2 function in vivo depends on interaction with UNC-37. As many T-box factors contain eh1 motifs, we suggest that interaction with Groucho-family corepressors is a common mechanism contributing to their activity.

The Paired-box protein PAX-3 regulates the choice between lateral and ventral epidermal cell fates in C. elegans

The development of the single cell layer skin or hypodermis of Caenorhabditis elegans is an excellent model for understanding cell fate specification and differentiation. Early in C. elegans embryogenesis, six rows of hypodermal cells adopt dorsal, lateral or ventral fates that go on to display distinct behaviors during larval life. Several transcription factors are known that function in specifying these major hypodermal cell fates, but our knowledge of the specification of these cell types is sparse, particularly in the case of the ventral hypodermal cells, which become Vulval Precursor Cells and form the vulval opening in response to extracellular signals. Previously, the gene pvl-4 was identified in a screen for mutants with defects in vulval development. We found by whole genome sequencing that pvl-4 is the Paired-box gene pax-3, which encodes the sole PAX-3 transcription factor homolog in C. elegans. pax-3 mutants show embryonic and larval lethality, and body morphology abnormalities indicative of hypodermal cell defects. We report that pax-3 is expressed in ventral P cells and their descendants during embryogenesis and early larval stages, and that in pax-3 reduction-of-function animals the ventral P cells undergo a cell fate transformation and express several markers of the lateral seam cell fate. Furthermore, forced expression of pax-3 in the lateral hypodermal cells causes them to lose expression of seam cell markers. We propose that pax-3 functions in the ventral hypodermal cells to prevent these cells from adopting the lateral seam cell fate. pax-3 represents the first gene required for specification solely of the ventral hypodermal fate in C. elegans providing insights into cell type diversification.

The study of Priapulus caudatus reveals conserved molecular patterning underlying different gut morphogenesis in the Ecdysozoa

Background

The digestive systems of animals can become highly specialized in response to their exploration and occupation of new ecological niches. Although studies on different animals have revealed commonalities in gut formation, the model systems Caenorhabditis elegans and Drosophila melanogaster, which belong to the invertebrate group Ecdysozoa, exhibit remarkable deviations in how their intestines develop. Their morphological and developmental idiosyncrasies have hindered reconstructions of ancestral gut characters for the Ecdysozoa, and limit comparisons with vertebrate models. In this respect, the phylogenetic position, and slow evolving morphological and molecular characters of marine priapulid worms advance them as a key group to decipher evolutionary events that occurred in the lineages leading to C. elegans and D. melanogaster.

Results

In the priapulid Priapulus caudatus, the gut consists of an ectodermal foregut and anus, and a mid region of at least partial endodermal origin. The inner gut develops into a 16-cell primordium devoid of visceral musculature, arranged in three mid tetrads and two posterior duplets. The mouth invaginates ventrally and shifts to a terminal anterior position as the ventral anterior ectoderm differentially proliferates. Contraction of the musculature occurs as the head region retracts into the trunk and resolves the definitive larval body plan. Despite obvious developmental differences with C. elegans and D. melanogaster, the expression in P. caudatus of the gut-related candidate genes NK2.1, foxQ2, FGF8/17/18, GATA456, HNF4, wnt1, and evx demonstrate three distinct evolutionarily conserved molecular profiles that correlate with morphologically identified sub-regions of the gut.

Conclusions

The comparative analysis of priapulid development suggests that a midgut formed by a single endodermal population of vegetal cells, a ventral mouth, and the blastoporal origin of the anus are ancestral features in the Ecdysozoa. Our molecular data on P. caudatus reveal a conserved ecdysozoan gut-patterning program and demonstrates that extreme morphological divergence has not been accompanied by major molecular innovations in transcriptional regulators during digestive system evolution in the Ecdysozoa. Our data help us understand the origins of the ecdysozoan body plan, including those of C. elegans and D. melanogaster, and this is critical for comparisons between these two prominent model systems and their vertebrate counterparts.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0139-z) contains supplementary material, which is available to authorized users.

Caenorhabditis elegans TBX-2 Directly Regulates Its Own Expression in a Negative Autoregulatory Loop

T-box genes often exhibit dynamic expression patterns, and their expression levels can be crucial for normal function. Despite the importance of these genes, there is little known about T-box gene regulation. We have focused on the Caenorhabditis elegans gene tbx-2 to understand how T-box gene expression is regulated, and here we demonstrate TBX-2 itself directly represses its own expression in a negative autoregulatory loop. tbx-2 is essential for normal pharyngeal muscle development, and a tbx-2 promoter gfp fusion (Ptbx-2::gfp) is transiently expressed in the pharynx during embryogenesis and in a small number of head neurons in larvae and adults. Reduced tbx-2 function resulted in ectopic Ptbx-2::gfp expression in the seam cells and gut in larvae and adults. Mutation of potential T-box binding sites within the tbx-2 promoter resulted in a similar pattern of ectopic Ptbx-2::gfp expression, and chromatin immunoprecipitation analyses show TBX-2 binds these sites in vivo. This pattern of ectopic Ptbx-2::gfp expression in tbx-2 mutants was very similar to that observed in mutants affecting the NF-Y complex, and our results comparing tbx-2 and nfyb-1 single- and double mutants suggest TBX-2 and NF-Y function in a single pathway to repress the tbx-2 promoter. The tbx-2 promoter is the first direct target identified for TBX-2, and we used it to ask whether SUMOylation is essential for TBX-2 repression. RNAi knockdown of SUMOylation pathway components led to ectopic Ptbx-2::gfp expression in the seam cells and gut. Ectopic Ptbx-2::gfp also was observed in the syncytial hypodermis, suggesting either the tbx-2 promoter is repressed by other SUMOylation dependent mechanisms, or that decreased SUMOylation leads to stable changes in seam cell nuclei as they fuse with the syncytial hypodermis. We suggest negative autoregulation is an important mechanism that allows precise control of tbx-2 expression levels and may allow rapid changes in gene expression during development.

The NF-Y complex negatively regulates Caenorhabditis elegans tbx-2 expression

T-box genes are frequently expressed in dynamic patterns during animal development, but the mechanisms controlling expression of these genes are not well understood. The Caenorhabditis elegans T-box gene tbx-2 is essential for development of the ABa-derived pharyngeal muscles, specification of neural cell fate in the HSN/PHB lineage, and adaptation in olfactory neurons. The tbx-2 expression pattern is complex, and expression has been described in pharyngeal precursors and body wall muscles during embryogenesis, and amphid sensory neurons and pharyngeal neurons in adults. To examine mechanisms regulating tbx-2 gene expression, we performed an RNAi screen of transcription factor genes in strains containing a Ptbx-2::gfp reporter and identified the Nuclear Factor Y (NF-Y) complex as a negative regulator of tbx-2 expression. NF-Y is a heterotrimeric CCAAT-binding complex consisting of A-C subunits, and reduction of the NF-Y subunits nfya-1, nfyb-1, or nfyc-1 by RNAi or using mutants results in ectopic Ptbx-2::gfp expression in hypodermal seam cells and gut. Mutation of two CCAAT-boxes in the tbx-2 promoter results in a similar pattern of ectopic Ptbx-2::gfp expression, suggesting NF-Y directly represses the tbx-2 promoter. tbx-2 mRNA is moderately increased in nfya-1 null mutants, indicating NF-Y represses expression of endogenous tbx-2. Finally we identify and characterize a second-site mutation that enhances lethality of a temperature sensitive tbx-2 mutant and show that this mutation is a deletion in the nfyb-1 gene. Together, these results identify NF-Y as an important regulator of tbx-2 function in vivo.

Transcriptional regulation of gene expression in C. elegans

Protein coding gene sequences are converted to mRNA by the highly regulated process of transcription. The precise temporal and spatial control of transcription for many genes is an essential part of development in metazoans. Thus, understanding the molecular mechanisms underlying transcriptional control is essential to understanding cell fate determination during embryogenesis, post-embryonic development, many environmental interactions, and disease-related processes. Studies of transcriptional regulation in C. elegans exploit its genomic simplicity and physical characteristics to define regulatory events with single-cell and minute-time-scale resolution. When combined with the genetics of the system, C. elegans offers a unique and powerful vantage point from which to study how chromatin-associated proteins and their modifications interact with transcription factors and their binding sites to yield precise control of gene expression through transcriptional regulation.

Sex differences in body composition, fat storage, and gene expression profile in Caenorhabditis elegans in response to dietary restriction

The metabolic and health-promoting effects of dietary restriction (DR) have been extensively studied in several species. The response to DR with respect to sex is essentially unknown. To address this question, we used the model organism Caenorhabditis elegans to analyze body composition and gene expression in males and hermaphrodites in response to DR. Unexpectedly, DR increased the fat-to-fat-free mass ratio and enlarged lipid droplets in both sexes to a similar extent. These effects were linked to a downregulation of the lipase-like 5 (lipl-5) gene in both sexes at two developmental stages. By contrast, the reductions in body size, protein content, and total RNA content in response to DR were more pronounced in hermaphrodites than in males. Functional enrichment analysis of gene expression data showed a DR-induced downregulation of several embryogenesis-associated genes concomitant with an ongoing expression of sperm-associated genes in hermaphrodites. In conclusion, DR increases fat stores in both sexes of C. elegans in the form of large and possibly lipolysis-resistant lipid droplets and markedly alters the reproductive program in hermaphrodites but not in males.

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.

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.

Predicting mutation outcome from early stochastic variation in genetic interaction partners

Many mutations, including those that cause disease, only have a detrimental effect in a subset of individuals. The reasons for this are usually unknown, but may include additional genetic variation and environmental risk factors. However, phenotypic discordance remains even in the absence of genetic variation, for example between monozygotic twins, and incomplete penetrance of mutations is frequent in isogenic model organisms in homogeneous environments. Here we propose a model for incomplete penetrance based on genetic interaction networks. Using Caenorhabditis elegans as a model system, we identify two compensation mechanisms that vary among individuals and influence mutation outcome. First, feedback induction of an ancestral gene duplicate differs across individuals, with high expression masking the effects of a mutation. This supports the hypothesis that redundancy is maintained in genomes to buffer stochastic developmental failure. Second, during normal embryonic development we find that there is substantial variation in the induction of molecular chaperones such as Hsp90 (DAF-21). Chaperones act as promiscuous buffers of genetic variation, and embryos with stronger induction of Hsp90 are less likely to be affected by an inherited mutation. Simultaneously quantifying the variation in these two independent responses allows the phenotypic outcome of a mutation to be more accurately predicted in individuals. Our model and methodology provide a framework for dissecting the causes of incomplete penetrance. Further, the results establish that inter-individual variation in both specific and more general buffering systems combine to determine the outcome inherited mutations in each individual.

A Genome-Wide RNAi Screen for Factors Involved in Neuronal Specification in Caenorhabditis elegans

One of the central goals of developmental neurobiology is to describe and understand the multi-tiered molecular events that control the progression of a fertilized egg to a terminally differentiated neuron. In the nematode Caenorhabditis elegans, the progression from egg to terminally differentiated neuron has been visually traced by lineage analysis. For example, the two gustatory neurons ASEL and ASER, a bilaterally symmetric neuron pair that is functionally lateralized, are generated from a fertilized egg through an invariant sequence of 11 cellular cleavages that occur stereotypically along specific cleavage planes. Molecular events that occur along this developmental pathway are only superficially understood. We take here an unbiased, genome-wide approach to identify genes that may act at any stage to ensure the correct differentiation of ASEL. Screening a genome-wide RNAi library that knocks-down 18,179 genes (94% of the genome), we identified 245 genes that affect the development of the ASEL neuron, such that the neuron is either not generated, its fate is converted to that of another cell, or cells from other lineage branches now adopt ASEL fate. We analyze in detail two factors that we identify from this screen: (1) the proneural gene hlh-14, which we find to be bilaterally expressed in the ASEL/R lineages despite their asymmetric lineage origins and which we find is required to generate neurons from several lineage branches including the ASE neurons, and (2) the COMPASS histone methyltransferase complex, which we find to be a critical embryonic inducer of ASEL/R asymmetry, acting upstream of the previously identified miRNA lsy-6. Our study represents the first comprehensive, genome-wide analysis of a single neuronal cell fate decision. The results of this analysis provide a starting point for future studies that will eventually lead to a more complete understanding of how individual neuronal cell types are generated from a single-cell embryo.

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.

SUMOylation-dependent function of a T-box transcriptional repressor in Caenorhabditis elegans

T-box transcription factors are crucial developmental regulators, and they have not previously been associated with SUMOylation. In Caenorhabditis elegans, the Tbx2 subfamily member TBX-2 (T-box protein 2) is required for anterior pharyngeal muscle development. TBX-2 interacts with SUMOylation pathway enzymes, and loss of these enzymes phenocopies tbx-2 mutants. We hypothesize that TBX-2 functions as a SUMOylation-dependent transcriptional repressor. TBX-2 contains two consensus SUMOylation sites conserved in many T-box transcriptional repressors, and we suggest that the function of these T-box factors may similarly involve SUMOylation.

Role of T-box gene tbx-2 for anterior foregut muscle development in C. elegans

During organogenesis, pluripotent precursor cells acquire a defined identity such as muscle or nerve. The transition from naïve precursor towards the differentiated state is characterized by sequential waves of gene expression that are determined by regulatory transcription factors. A key question is how transcriptional circuitry dictates the succession of events that accompanies developmental competence, cell fate specification and differentiation. To address this question, we have examined how anterior muscles are established within the Caenorhabditis elegans foregut (pharynx). We find that the T-box transcription factor tbx-2 is essential to form anterior pharyngeal muscles from the ABa blastomere. In the absence of tbx-2 function, ABa-derived cells initiate development normally: they receive glp-1/Notch signaling cues, activate the T-box gene TBX-38 and express the organ selector gene PHA-4/FoxA. However, these cells subsequently arrest development, extinguish PHA-4 and fail to activate PHA-4 target genes. tbx-2 mutant cells do not undergo apoptosis and there is no evidence for adoption of an alternative fate. TBX-2 is expressed in ABa descendants and depends on activation by pha-4 and repression by components of glp-1/Notch signaling. Our analysis suggests that a positive feedback loop between tbx-2 and pha-4 is required for ABa-derived precursors to commit to pharyngeal muscle fate.

Combinatorial RNA interference in Caenorhabditis elegans reveals that redundancy between gene duplicates can be maintained for more than 80 million years of evolution

High-throughput combinatorial RNAi demonstrates that many duplicated genes in C. elegans can retain redundant functions for more than 80 million years

Background

Systematic analyses of loss-of-function phenotypes have been carried out for most genes in Saccharomyces cerevisiae, Caenorhabditis elegans, and Drosophila melanogaster. Although such studies vastly expand our knowledge of single gene function, they do not address redundancy in genetic networks. Developing tools for the systematic mapping of genetic interactions is thus a key step in exploring the relationship between genotype and phenotype.

Results

We established conditions for RNA interference (RNAi) in C. elegans to target multiple genes simultaneously in a high-throughput setting. Using this approach, we can detect the great majority of previously known synthetic genetic interactions. We used this assay to examine the redundancy of duplicated genes in the genome of C. elegans that correspond to single orthologs in S. cerevisiae or D. melanogaster and identified 16 pairs of duplicated genes that have redundant functions. Remarkably, 14 of these redundant gene pairs were duplicated before the divergence of C. elegans and C. briggsae 80-110 million years ago, suggesting that there has been selective pressure to maintain the overlap in function between some gene duplicates.

Conclusion

We established a high throughput method for examining genetic interactions using combinatorial RNAi in C. elegans. Using this technique, we demonstrated that many duplicated genes can retain redundant functions for more than 80 million years of evolution. This provides strong support for evolutionary models that predict that genetic redundancy between duplicated genes can be actively maintained by natural selection and is not just a transient side effect of recent gene duplication events.

Tbx5-dependent rheostatic control of cardiac gene expression and morphogenesis

Dominant mutations in the T-box transcription factor gene TBX5 cause Holt-Oram syndrome (HOS), an inherited human disease characterized by upper limb malformations and congenital heart defects (CHDs) of variable severity. We hypothesize that minor alterations in the dosage of Tbx5 directly influences severity of CHDs. Using a mouse allelic series, we show a sensitive inverse correlation between Tbx5 dosage and abnormal cardiac morphogenesis and gene expression. The CHDs found in mice harbouring a hypomorphic allele of Tbx5 (Tbx5(lox/+) mice) are less pronounced than those found in Tbx5 haploinsufficient mice (Tbx5(del/+)), and homozygous hypomorphic (Tbx5(lox/lox)) embryos have noticeably more advanced cardiac development than Tbx5 null (Tbx5(del/del)) embryos. Examination of target gene expression across the allelic series uncovers very fine sensitivity across the range of Tbx5 dosages, in which some genes respond dramatically differently to only 15% differences in Tbx5 mRNA levels. This analysis was expanded to a genome-wide level, which uncovered a Tbx5 dosage-sensitive genetic program involving a network of cardiac transcription factors, developmentally important cell-cell signaling molecules, and ion channel proteins. These results indicate an exquisite sensitivity of the developing heart to Tbx5 dosage and provide significant insight into the transcriptional and cellular mechanisms that are disrupted in CHDs.

The T-box factor TBX-2 and the SUMO conjugating enzyme UBC-9 are required for ABa-derived pharyngeal muscle in C. elegans

The C. elegans pharynx is produced from the embryonic blastomeres ABa and MS. Pharyngeal fate in the ABa lineage is specified by the combined activities of GLP-1/Notch-mediated signals and the TBX-37 and TBX-38 T-box transcription factors. Here, we show another T-box factor TBX-2 also functions in ABa-derived pharyngeal development. tbx-2 mutants arrest as L1 larvae lacking most or all ABa-derived pharyngeal muscles. In comparison, tbx-2 mutants retain ABa-derived marginal cells and pharyngeal muscles derived from MS. A tbx-2Colon, two colonsgfp translational fusion is expressed in a dynamic pattern in C. elegans embryos beginning near the 100-cell stage. Early expression is limited to a small number of cells, which likely include the ABa-derived pharyngeal precursors, while later expression is observed in body wall muscles and a subset of pharyngeal neurons. TBX-2 contains 2 consensus sumoylation sites, and it interacts in a yeast two-hybrid assay with the UBC-9 and GEI-17 components of the C. elegans SUMO-conjugating pathway. ubc-9(RNAi) has been previously shown to cause variable embryonic and larval arrest, and we find that, like tbx-2 mutants, ubc-9(RNAi) animals lack ABa-derived pharyngeal muscles. ubc-9(RNAi) also alters the subnuclear distribution of TBX-2::GFP fusion protein, suggesting that UBC-9 and TBX-2 interact in C. elegans. Together, these results indicate that TBX-2 and SUMO-conjugating enzymes are necessary for ABa-derived pharyngeal muscle, and we hypothesize that TBX-2 function requires sumoylation. Sumoylation is increasingly recognized as an important mechanism controlling activity of many nuclear factors, and these results provide the first evidence that T-box factor activity may require sumoylation.

A compendium of Caenorhabditis elegans regulatory transcription factors: a resource for mapping transcription regulatory networks

BACKGROUND

Transcription regulatory networks are composed of interactions between transcription factors and their target genes. Whereas unicellular networks have been studied extensively, metazoan transcription regulatory networks remain largely unexplored. Caenorhabditis elegans provides a powerful model to study such metazoan networks because its genome is completely sequenced and many functional genomic tools are available. While C. elegans gene predictions have undergone continuous refinement, this is not true for the annotation of functional transcription factors. The comprehensive identification of transcription factors is essential for the systematic mapping of transcription regulatory networks because it enables the creation of physical transcription factor resources that can be used in assays to map interactions between transcription factors and their target genes.

RESULTS

By computational searches and extensive manual curation, we have identified a compendium of 934 transcription factor genes (referred to as wTF2.0). We find that manual curation drastically reduces the number of both false positive and false negative transcription factor predictions. We discuss how transcription factor splice variants and dimer formation may affect the total number of functional transcription factors. In contrast to mouse transcription factor genes, we find that C. elegans transcription factor genes do not undergo significantly more splicing than other genes. This difference may contribute to differences in organism complexity. We identify candidate redundant worm transcription factor genes and orthologous worm and human transcription factor pairs. Finally, we discuss how wTF2.0 can be used together with physical transcription factor clone resources to facilitate the systematic mapping of C. elegans transcription regulatory networks.

CONCLUSION

wTF2.0 provides a starting point to decipher the transcription regulatory networks that control metazoan development and function.

Synthetic lethal analysis of Caenorhabditis elegans posterior embryonic patterning genes identifies conserved genetic interactions

Phenotypic robustness is evidenced when single-gene mutations do not result in an obvious phenotype. It has been suggested that such phenotypic stability results from 'buffering' activities of homologous genes as well as non-homologous genes acting in parallel pathways. One approach to characterizing mechanisms of phenotypic robustness is to identify genetic interactions, specifically, double mutants where buffering is compromised. To identify interactions among genes implicated in posterior patterning of the Caenorhabditis elegans embryo, we measured synthetic lethality following RNA interference of 22 genes in 15 mutant strains. A pair of homologous T-box transcription factors (tbx-8 and tbx-9) is found to interact in both C. elegans and C. briggsae, indicating that their compensatory function is conserved. Furthermore, a muscle module is defined by transitive interactions between the MyoD homolog hlh-1, another basic helix-loop-helix transcription factor, hnd-1, and the MADS-box transcription factor unc-120. Genetic interactions within a homologous set of genes involved in vertebrate myogenesis indicate broad conservation of the muscle module and suggest that other genetic modules identified in C. elegans will be conserved.

Identification of lineage-specific zygotic transcripts in early Caenorhabditis elegans embryos

During Caenorhabditis elegans embryogenesis, a maternally supplied transcription factor, SKN-1, is required for the specification of the mesendodermal precursor, EMS, in the 4-cell stage embryo. When EMS divides, it gives rise to a mesoderm-restricted precursor, MS, and an endoderm-restricted precursor, E. To systematically identify genes that function as key regulators of MS and/or E-derived tissues, we identified, by microarray analyses, genes that are newly transcribed within a short developmental window (approximately 30 min) encompassing the generation and fate specification of the MS and E blastomeres. By comparing total cDNAs generated from individual, carefully staged embryos, we identified 275 genes up-regulated in 12-cell embryos compared to 4-cell embryos. Fifty of these 275 genes are down-regulated in 12-cell skn-1 mutant embryos and are designated skn-1-dependent zygotic (sdz) genes. The spatial and temporal expression patterns in C. elegans embryos of 10 randomly selected sdz genes were analyzed by a nuclear GFP reporter driven by the endogenous 5' regulatory sequence of each gene. GFP expression, although absent at the 4-cell stage, was detected at the 12- to 16-cell stage for all 10 genes and was restricted to EMS-derived lineages for 7 of the 10. Among the seven lineage-specific genes, three genes are expressed equally in both MS and E lineages, two are expressed exclusively or predominantly in the MS lineage, and two are expressed exclusively in the E lineage. Depletion of skn-1 by RNAi abolishes the expression of all seven reporter transgenes in vivo, confirming that these genes are indeed skn-1 dependent. These results demonstrate the successful combination of single-staged embryo cDNAs, genetic mutants, and whole transcriptome microarray analysis to identify stage- and lineage-specific transcripts in early C. elegans embryos.

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

Maternal and zygotic activities of the homeodomain protein PAL-1 specify the identity and maintain the development of the multipotent C blastomere lineage in the C. elegans embryo. To identify PAL-1 regulatory target genes, we used microarrays to compare transcript abundance in wild-type embryos with mutant embryos lacking a C blastomere and to mutant embryos with extra C blastomeres. pal-1-dependent C-lineage expression was verified for select candidate target genes by reporter gene analysis, though many of the target genes are expressed in additional lineages as well. The set of validated target genes includes 12 transcription factors, an uncharacterized wingless ligand and five uncharacterized genes. Phenotypic analysis demonstrates that the identified PAL-1 target genes affect specification, differentiation and morphogenesis of C-lineage cells. In particular, we show that cell fate-specific genes (or tissue identity genes) and a posterior HOX gene are activated in lineage-specific fashion. Transcription of targets is initiated in four temporal phases, which together with their spatial expression patterns leads to a model of the regulatory network specified by PAL-1.