The embryonic muscle transcriptome of Caenorhabditis elegans

Advancing insights into microgravity induced muscle changes using Caenorhabditis elegans as a model organism

Spaceflight presents significant challenges to the physiological state of living organisms. This can be due to the microgravity environment experienced during long-term space missions, resulting in alterations in muscle structure and function, such as atrophy. However, a comprehensive understanding of the adaptive mechanisms of biological systems is required to devise potential solutions and therapeutic approaches for adapting to spaceflight conditions. This review examines the current understanding of the challenges posed by spaceflight on physiological changes, alterations in metabolism, dysregulation of pathways and the suitability and advantages of using the model organism Caenorhabditis elegans nematodes to study the effects of spaceflight. Research has shown that changes in the gene and protein composition of nematodes significantly occur across various larval stages and rearing environments, including both microgravity and Earth gravity settings, often mirroring changes observed in astronauts. Additionally, the review explores significant insights into the fundamental metabolic changes associated with muscle atrophy and growth, which could lead to the development of diagnostic biomarkers and innovative techniques to prevent and counteract muscle atrophy. These insights not only advance our understanding of microgravity-induced muscle atrophy but also lay the groundwork for the development of targeted interventions to mitigate its effects in the future.

Epithelial UNC-23 limits mechanical stress to maintain glia-neuron architecture in C. elegans

For an organ to maintain correct architecture and function, its diverse cellular components must coordinate their size and shape. Although cell-intrinsic mechanisms driving homotypic cell-cell coordination are known, it is unclear how cell shape is regulated across heterotypic cells. We find that epithelial cells maintain the shape of neighboring sense-organ glia-neuron units in adult Caenorhabditis elegans (C. elegans). Hsp co-chaperone UNC-23/BAG2 prevents epithelial cell shape from deforming, and its loss causes head epithelia to stretch aberrantly during animal movement. In the sense-organ glia, amphid sheath (AMsh), this causes progressive fibroblast growth factor receptor (FGFR)-dependent disruption of the glial apical cytoskeleton. Resultant glial cell shape alteration causes concomitant shape change in glia-associated neuron endings. Epithelial UNC-23 maintenance of glia-neuron shape is specific both spatially, within a defined anatomical zone, and temporally, in a developmentally critical period. As all molecular components uncovered are broadly conserved across central and peripheral nervous systems, we posit that epithelia may similarly regulate glia-neuron architecture cross-species.

Transcriptional landscape of a hypoxia response identifies cell-specific pathways for adaptation

How the HIF-1 (Hypoxia-Inducible) transcription factor drives and coordinates distinct responses to low oxygen across diverse cell types is poorly understood. We present a multi-tissue single-cell gene-expression atlas of the hypoxia response of the nematode Caenorhabditis elegans . This atlas highlights how cell-type-specific HIF-1 responses overlap and diverge among and within neuronal, intestinal, and muscle tissues. Using the atlas to guide functional analyses of candidate muscle-specific HIF-1 effectors, we discovered that HIF-1 activation drives downregulation of the tspo-1 ( TSPO, Translocator Protein) gene in vulval muscle cells to modulate a hypoxia-driven change in locomotion caused by contraction of body-wall muscle cells. We further showed that in human cardiomyocytes HIF-1 activation decreases levels of TSPO and thereby alters intracellular cholesterol transport and the mitochondrial network. We suggest that TSPO-1 is an evolutionarily conserved mediator of HIF-1-dependent modulation of muscle and conclude that our gene-expression atlas can help reveal how HIF-1 drives cell-specific adaptations to hypoxia.

PBRM-1/PBAF-regulated genes in a multipotent progenitor in Caenorhabditis elegans

The Caenorhabditis elegans somatic gonadal precursors (SGPs) are multipotent progenitors that generate all somatic cells of the adult reproductive system. The 2 SGPs originate in the mesodermal layer and are born through a division that produces one SGP and one head mesodermal cell (hmc). One hmc terminally differentiates, and the other dies by programmed cell death. The polybromo-associated BAF (PBAF) chromatin remodeling complex promotes the multipotent SGP fate. The complete loss of PBAF causes lethality, so we used a combination of Cre/lox recombination and GFP nanobody-directed protein degradation to eliminate PBRM-1, the signature subunit of the PBAF complex, from 83 mesodermal cells, including SGPs, body muscles, and the hmc. We used RNA sequencing to identify genes acting downstream of PBAF in these cells and identified 1,955 transcripts that were significantly differentially expressed between pbrm-1(-) and pbrm-1(+) in the mesoderm of L1 larvae. We found that genes involved in muscle cell function were overrepresented; most of these genes had lower expression in the absence of PBRM-1, suggesting that PBAF promotes muscle differentiation. Among the differentially expressed genes were 125 that are normally expressed at higher levels in SGP vs hmc and positively regulated by pbrm-1 and 53 that are normally expressed at higher levels in hmc vs SGP and are negatively regulated by pbrm-1; these are candidate regulators of the SGP/hmc fate decision. We validated one candidate gene using a fluorescent reporter; the hsp-12.3 reporter was derepressed in SGPs in pbrm-1 mutants, suggesting that hsp-12.3 expression is normally repressed by pbrm-1 in SGPs.

EAT-2 attenuates C. elegans development via metabolic remodeling in a chemically defined food environment

Dietary intake and nutrient composition regulate animal growth and development; however, the underlying mechanisms remain elusive. Our previous study has shown that either the mammalian deafness homolog gene tmc-1 or its downstream acetylcholine receptor gene eat-2 attenuates Caenorhabditis elegans development in a chemically defined food CeMM (C. elegans maintenance medium) environment, but the underpinning mechanisms are not well-understood. Here, we found that, in CeMM food environment, for both eat-2 and tmc-1 fast-growing mutants, several fatty acid synthesis and elongation genes were highly expressed, while many fatty acid β-oxidation genes were repressed. Accordingly, dietary supplementation of individual fatty acids, such as monomethyl branch chain fatty acid C17ISO, palmitic acid and stearic acid significantly promoted wild-type animal development on CeMM, and mutations in either C17ISO synthesis gene elo-5 or elo-6 slowed the rapid growth of eat-2 mutant. Tissue-specific rescue experiments showed that elo-6 promoted animal development mainly in the intestine. Furthermore, transcriptome and metabolome analyses revealed that elo-6/C17ISO regulation of C. elegans development may be correlated with up-regulating expression of cuticle synthetic and hedgehog signaling genes, as well as promoting biosynthesis of amino acids, amino acid derivatives and vitamins. Correspondingly, we found that amino acid derivative S-adenosylmethionine and its upstream metabolite methionine sulfoxide significantly promoted C. elegans development on CeMM. This study demonstrated that C17ISO, palmitic acid, stearic acid, S-adenosylmethionine and methionine sulfoxide inhibited or bypassed the TMC-1 and EAT-2-mediated attenuation of development via metabolic remodeling, and allowed the animals to adapt to the new nutritional niche.

An updated C. elegans nuclear body muscle transcriptome for studies in muscle formation and function

The body muscle is an important tissue used in organisms for proper viability and locomotion. Although this tissue is generally well studied and characterized, and many pathways have been elucidated throughout the years, we still lack a comprehensive understanding of its transcriptome and how it controls muscle development and function. Here, we have updated a nuclear FACS sorting-based methodology to isolate and sequence a high-quality muscle transcriptome from Caenorhabditis elegans mixed-stage animals. We have identified 2848 muscle-specific protein-coding genes, including 78 transcription factors and 206 protein-coding genes containing an RNA binding domain. We studied their interaction network, performed a detailed promoter analysis, and identified novel muscle-specific cis-acting elements. We have also identified 16 high-quality muscle-specific miRNAs, studied their function in vivo using fluorochrome-based analyses, and developed a high-quality C. elegans miRNA interactome incorporating other muscle-specific datasets produced by our lab and others.Our study expands our understanding of how muscle tissue functions in C. elegans andin turn provides results that can in the future be applied to humans to study muscular-related diseases.

Neuronal HLH-30/TFEB modulates peripheral mitochondrial fragmentation to improve thermoresistance in Caenorhabditis elegans

Transcription factor EB (TFEB) is a conserved master transcriptional activator of autophagy and lysosomal genes that modulates organismal lifespan regulation and stress resistance. As neurons can coordinate organism-wide processes, we investigated the role of neuronal TFEB in stress resistance and longevity. To this end, the Caenorhabditis elegans TFEB ortholog, hlh-30, was rescued panneuronally in hlh-30 loss of function mutants. While important in the long lifespan of daf-2 animals, neuronal HLH-30/TFEB was not sufficient to restore normal lifespan in short-lived hlh-30 mutants. However, neuronal HLH-30/TFEB rescue mediated robust improvements in the heat stress resistance of wildtype but not daf-2 animals. Notably, these mechanisms can be uncoupled, as neuronal HLH-30/TFEB requires DAF-16/FOXO to regulate longevity but not thermoresistance. Through further transcriptomics profiling and functional analysis, we discovered that neuronal HLH-30/TFEB modulates neurotransmission through the hitherto uncharacterized protein W06A11.1 by inducing peripheral mitochondrial fragmentation and organismal heat stress resistance in a non-cell autonomous manner. Taken together, this study uncovers a novel mechanism of heat stress protection mediated by neuronal HLH-30/TFEB.

Unconventional Myosins from Caenorhabditis elegans as a Probe to Study Human Orthologues

Unconventional myosins are a superfamily of actin-based motor proteins that perform a number of roles in fundamental cellular processes, including (but not limited to) intracellular trafficking, cell motility, endocytosis, exocytosis and cytokinesis. 40 myosins genes have been identified in humans, which belong to different 12 classes based on their domain structure and organisation. These genes are widely expressed in different tissues, and mutations leading to loss of function are associated with a wide variety of pathologies while over-expression often results in cancer. Caenorhabditis elegans (C. elegans) is a small, free-living, non-parasitic nematode. ~38% of the genome of C. elegans has predicted orthologues in the human genome, making it a valuable tool to study the function of human counterparts and human diseases. To date, 8 unconventional myosin genes have been identified in the nematode, from 6 different classes with high homology to human paralogues. The hum-1 and hum-5 (heavy chain of an unconventional myosin) genes encode myosin of class I, hum-2 of class V, hum-3 and hum-8 of class VI, hum-6 of class VII and hum-7 of class IX. The hum-4 gene encodes a high molecular mass myosin (307 kDa) that is one of the most highly divergent myosins and is a member of class XII. Mutations in many of the human orthologues are lethal, indicating their essential properties. However, a functional characterisation for many of these genes in C. elegans has not yet been performed. This article reviews the current knowledge of unconventional myosin genes in C. elegans and explores the potential use of the nematode to study the function and regulation of myosin motors to provide valuable insights into their role in diseases.

Molecular lesions in alleles of the Caenorhabditis elegans lin-11 gene

The LIM homeodomain transcription factor LIN-11 is a key regulator of vulva, uterine, and neuron development in C. elegans. Multiple alleles of lin-11 are available, but none had been sequenced. We found that the reference allele, n389, is a 15900 bp deletion that also affects two other protein-coding genes, ZC247.1 and ZC247.2. The frequently used n566 allele is a 288bp deletion located in an intron and affecting the splice acceptor site.

HLH-1 Modulates Muscle Proteostasis During Caenorhabditis elegans Larval Development

Muscle proteostasis is shaped by the myogenic transcription factor MyoD which regulates the expression of chaperones during muscle differentiation. Whether MyoD can also modulate chaperone expression in terminally differentiated muscle cells remains open. Here we utilized a temperature-sensitive (ts) conditional knockdown nonsense mutation in MyoD ortholog in C. elegans, HLH-1, to ask whether MyoD plays a role in maintaining muscle proteostasis post myogenesis. We showed that hlh-1 is expressed during larval development and that hlh-1 knockdown at the first, second, or third larval stages resulted in severe defects in motility and muscle organization. Motility defects and myofilament organization were rescued when the clearance of hlh-1(ts) mRNA was inhibited, and hlh-1 mRNA levels were restored. Moreover, hlh-1 knockdown modulated the expression of chaperones with putative HLH-1 binding sites in their promoters, supporting HLH-1 role in muscle maintenance during larval development. Finally, mild disruption of hlh-1 expression during development resulted in earlier dysregulation of muscle maintenance and function during adulthood. We propose that the differentiation transcription factor, HLH-1, contributes to muscle maintenance and regulates cell-specific chaperone expression post differentiation. HLH-1 may thus impact muscle proteostasis and potentially the onset and manifestation of sarcopenia.

Advances in our understanding of nematode ion channels as potential anthelmintic targets

Ion channels are specialized multimeric proteins that underlie cell excitability. These channels integrate with a variety of neuromuscular and biological functions. In nematodes, the physiological behaviors including locomotion, navigation, feeding and reproduction, are regulated by these protein entities. Majority of the antinematodal chemotherapeutics target the ion channels to disrupt essential biological functions. Here, we have summarized current advances in our understanding of nematode ion channel pharmacology. We review cys-loop ligand gated ion channels (LGICs), including nicotinic acetylcholine receptors (nAChRs), acetylcholine-chloride gated ion channels (ACCs), glutamate-gated chloride channels (GluCls), and GABA (γ-aminobutyric acid) receptors, and other ionotropic receptors (transient receptor potential (TRP) channels and potassium ion channels). We have provided an update on the pharmacological properties of these channels from various nematodes. This article catalogs the differences in ion channel composition and resulting pharmacology in the phylum Nematoda. This diversity in ion channel subunit repertoire and pharmacology emphasizes the importance of pursuing species-specific drug target research. In this review, we have provided an overview of recent advances in techniques and functional assays available for screening ion channel properties and their application.

Caenorhabditis elegans ETR-1/CELF has broad effects on the muscle cell transcriptome, including genes that regulate translation and neuroblast migration

Migration of neuroblasts and neurons from their birthplace is central to the formation of neural circuits and networks. ETR-1 is the Caenorhabditis elegans homolog of the CELF1 (CUGBP, ELAV-like family 1) RNA-processing factor involved in neuromuscular disorders. etr-1 regulates body wall muscle differentiation. Our previous work showed that etr-1 in muscle has a non-autonomous role in neuronal migration, suggesting that ETR-1 is involved in the production of a signal emanating from body wall muscle that controls neuroblast migration and that interacts with Wnt signaling. etr-1 is extensively alternatively-spliced, and we identified the viable etr-1(lq61) mutant, caused by a stop codon in alternatively-spliced exon 8 and only affecting etr-1 isoforms containing exon 8. We took advantage of viable etr-1(lq61) to identify potential RNA targets of ETR-1 in body wall muscle using a combination of fluorescence activated cell sorting (FACS) of body wall muscles from wild-type and etr-1(lq61) and subsequent RNA-seq. This analysis revealed genes whose splicing and transcript levels were controlled by ETR-1 exon 8 isoforms, and represented a broad spectrum of genes involved in muscle differentiation, myofilament lattice structure, and physiology. Genes with transcripts underrepresented in etr-1(lq61) included those involved in ribosome function and translation, similar to potential CELF1 targets identified in chick cardiomyocytes. This suggests that at least some targets of ETR-1 might be conserved in vertebrates, and that ETR-1 might generally stimulate translation in muscles. As proof-of-principle, a functional analysis of a subset of ETR-1 targets revealed genes involved in AQR and PQR neuronal migration. One such gene, lev-11/tropomyosin, requires ETR-1 for alternative splicing, and another, unc-52/perlecan, requires ETR-1 for the production of long isoforms containing 3' exons. In sum, these studies identified gene targets of ETR-1/CELF1 in muscles, which included genes involved in muscle development and physiology, and genes with novel roles in neuronal migration.

Temperature regulates synaptic subcellular specificity mediated by inhibitory glutamate signaling

Environmental factors such as temperature affect neuronal activity and development. However, it remains unknown whether and how they affect synaptic subcellular specificity. Here, using the nematode Caenorhabditis elegans AIY interneurons as a model, we found that high cultivation temperature robustly induces defects in synaptic subcellular specificity through glutamatergic neurotransmission. Furthermore, we determined that the functional glutamate is mainly released by the ASH sensory neurons and sensed by two conserved inhibitory glutamate-gated chloride channels GLC-3 and GLC-4 in AIY. Our work not only presents a novel neurotransmission-dependent mechanism underlying the synaptic subcellular specificity, but also provides a potential mechanistic insight into high-temperature-induced neurological defects.

In Vivo Simultaneous Analysis of Gene Expression by Dual-Color Luciferases in Caenorhabditis elegans

Both fluorescent and luminescent observation are widely used to examine real-time gene expression patterns in living organisms. Several fluuorescent and luminescent proteins with specific optical properties have been developed and applied for simultaneous, multi-color observation of more than two gene expression profiles. Compared to fluorescent proteins, however, the application of multi-color luminescent imaging in living organisms is still limited. In this study, we introduced two-color luciferases into the soil nematode C. elegans and performed simultaneous analysis of two gene expression profiles. Using a green-emitting luciferase Eluc (emerald luciferase) and red-emitting luciferase SLR (stable luciferase red), the expression patterns of two genes were simultaneously observed in single animals from embryonic to adult stages over its whole life span. In addition, dual gene activities were observed at the single embryo level, with the simultaneous observation of morphological changes. These are the first application of a two-color luciferase system into a whole animal and suggest that precise relationship of expression patterns of multiple genes of interest can be analyzed over the whole life of the animal, dependent on the changes in genetic and/or environmental conditions.

FLN-1/filamin is required to anchor the actomyosin cytoskeleton and for global organization of sub-cellular organelles in a contractile tissue

Actomyosin networks are organized in space, direction, size, and connectivity to produce coordinated contractions across cells. We use the C. elegans spermatheca, a tube composed of contractile myoepithelial cells, to study how actomyosin structures are organized. FLN-1/filamin is required for the formation and stabilization of a regular array of parallel, contractile, actomyosin fibers in this tissue. Loss of fln-1 results in the detachment of actin fibers from the basal surface, which then accumulate along the cell junctions and are stabilized by spectrin. In addition, actin and myosin are captured at the nucleus by the linker of nucleoskeleton and cytoskeleton complex (LINC) complex, where they form large foci. Nuclear positioning and morphology, distribution of the endoplasmic reticulum and the mitochondrial network are also disrupted. These results demonstrate that filamin is required to prevent large actin bundle formation and detachment, to prevent excess nuclear localization of actin and myosin, and to ensure correct positioning of organelles.

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.

mRNA profiling reveals significant transcriptional differences between a multipotent progenitor and its differentiated sister

BACKGROUND

The two Caenorhabditis elegans somatic gonadal precursors (SGPs) are multipotent progenitors that generate all somatic tissues of the adult reproductive system. The sister cells of the SGPs are two head mesodermal cells (hmcs); one hmc dies by programmed cell death and the other terminally differentiates. Thus, a single cell division gives rise to one multipotent progenitor and one differentiated cell with identical lineage histories. We compared the transcriptomes of SGPs and hmcs in order to learn the determinants of multipotency and differentiation in this lineage.

RESULTS

We generated a strain that expressed fluorescent markers specifically in SGPs (ehn-3A::tdTomato) and hmcs (bgal-1::GFP). We dissociated cells from animals after the SGP/hmc cell division, but before the SGPs had further divided, and subjected the dissociated cells to fluorescence-activated cell sorting to collect isolated SGPs and hmcs. We analyzed the transcriptomes of these cells and found that 5912 transcripts were significantly differentially expressed, with at least two-fold change in expression, between the two cell types. The hmc-biased genes were enriched with those that are characteristic of neurons. The SGP-biased genes were enriched with those indicative of cell proliferation and development. We assessed the validity of our differentially expressed genes by examining existing reporters for five of the 10 genes with the most significantly biased expression in SGPs and found that two showed expression in SGPs. For one reporter that did not show expression in SGPs, we generated a GFP knock-in using CRISPR/Cas9. This reporter, in the native genomic context, was expressed in SGPs.

CONCLUSIONS

We found that the transcriptional profiles of SGPs and hmcs are strikingly different. The hmc-biased genes are enriched with those that encode synaptic transmission machinery, which strongly suggests that it has neuron-like signaling properties. In contrast, the SGP-biased genes are enriched with genes that encode factors involved in transcription and translation, as would be expected from a cell preparing to undergo proliferative divisions. Mediators of multipotency are likely to be among the genes differentially expressed in SGPs.

Polycomb and Notch signaling regulate cell proliferation potential during Caenorhabditis elegans life cycle

Stable cell fate is an essential feature for multicellular organisms in which individual cells achieve specialized functions. Caenorhabditis elegans is a great model to analyze the determinants of cell fate stability because of its invariant lineage. We present a tractable cell fate challenge system that uses the induction of fate-specifying transcription factors. We show that wild-type differentiated animals are highly resistant to fate challenge. Removal of heterochromatin marks showed marked differences: the absence of histone 3 lysine 9 methylation (H3K9) has no effect on fate stability, whereas Polycomb homolog mes-2 mutants lacking H3K27 methylation terminally arrest larval development upon fate challenge. Unexpectedly, the arrest correlated with widespread cell proliferation rather than transdifferentiation. Using a candidate RNAi larval arrest-rescue screen, we show that the LIN-12Notch pathway is essential for hyperplasia induction. Moreover, Notch signaling appears downstream of food-sensing pathways, as dauers and first larval stage diapause animals are resistant to fate challenge. Our results demonstrate an equilibrium between proliferation and differentiation regulated by Polycomb and Notch signaling in the soma during the nematode life cycle.

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

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

STITCHER 2.0: primer design for overlapping PCR applications

Overlapping polymerase chain reaction (PCR) is a common technique used by researchers in very diverse fields that enables the user to ‘stitch’ individual pieces of DNA together. Previously, we have reported a web based tool called STITCHER that provides a platform for researchers to automate the design of primers for overlapping PCR applications. Here we present STITCHER 2.0, which represents a substantial update to STITCHER. STITCHER 2.0 is a newly designed web tool that automates the design of primers for overlapping PCR. Unlike STITCHER, STITCHER 2.0 considers diverse algorithmic parameters, and returns multiple result files that include a facility for the user to draw their own primers as well as comprehensive visual guides to the user’s input, output, and designed primers. These result files provide greater control and insight during experimental design and troubleshooting. STITCHER 2.0 is freely available to all users without signup or login requirements and can be accessed at the following webpage: www.ohalloranlab.net/STITCHER2.html.

A Differentiation Transcription Factor Establishes Muscle-Specific Proteostasis in Caenorhabditis elegans

Safeguarding the proteome is central to the health of the cell. In multi-cellular organisms, the composition of the proteome, and by extension, protein-folding requirements, varies between cells. In agreement, chaperone network composition differs between tissues. Here, we ask how chaperone expression is regulated in a cell type-specific manner and whether cellular differentiation affects chaperone expression. Our bioinformatics analyses show that the myogenic transcription factor HLH-1 (MyoD) can bind to the promoters of chaperone genes expressed or required for the folding of muscle proteins. To test this experimentally, we employed HLH-1 myogenic potential to genetically modulate cellular differentiation of Caenorhabditis elegans embryonic cells by ectopically expressing HLH-1 in all cells of the embryo and monitoring chaperone expression. We found that HLH-1-dependent myogenic conversion specifically induced the expression of putative HLH-1-regulated chaperones in differentiating muscle cells. Moreover, disrupting the putative HLH-1-binding sites on ubiquitously expressed daf-21(Hsp90) and muscle-enriched hsp-12.2(sHsp) promoters abolished their myogenic-dependent expression. Disrupting HLH-1 function in muscle cells reduced the expression of putative HLH-1-regulated chaperones and compromised muscle proteostasis during and after embryogenesis. In turn, we found that modulating the expression of muscle chaperones disrupted the folding and assembly of muscle proteins and thus, myogenesis. Moreover, muscle-specific over-expression of the DNAJB6 homolog DNJ-24, a limb-girdle muscular dystrophy-associated chaperone, disrupted the muscle chaperone network and exposed synthetic motility defects. We propose that cellular differentiation could establish a proteostasis network dedicated to the folding and maintenance of the muscle proteome. Such cell-specific proteostasis networks can explain the selective vulnerability that many diseases of protein misfolding exhibit even when the misfolded protein is ubiquitously expressed.

Molecular chaperones protect proteins from misfolding and aggregation. In multi-cellular organisms, the composition and expression levels of chaperones vary between tissues. However, little is known of how such differential expression is regulated. We hypothesized that the cellular differentiation that regulates the cell-type specific expression program may be involved in establishing a cell-type specific chaperone network. To test this possibility, we addressed the myogenic commitment transcription factor HLH-1 (CeMyoD) that converts embryonic cells to muscle cells in Caenorhabditis elegans. We demonstrated that HLH-1 regulates the expression of muscle chaperones during muscle differentiation. Moreover, we showed that HLH-1-dependent expression of chaperones is required for embryonic development and muscle function. We propose that cellular differentiation results in cell-specific differences in the chaperone network that may be detrimental in terms of the susceptibility of neurons and muscle cells to protein misfolding diseases.

Tissue enrichment analysis for C. elegans genomics

Background

Over the last ten years, there has been explosive development in methods for measuring gene expression. These methods can identify thousands of genes altered between conditions, but understanding these datasets and forming hypotheses based on them remains challenging. One way to analyze these datasets is to associate ontologies (hierarchical, descriptive vocabularies with controlled relations between terms) with genes and to look for enrichment of specific terms. Although Gene Ontology (GO) is available for Caenorhabditis elegans, it does not include anatomical information.

Results

We have developed a tool for identifying enrichment of C. elegans tissues among gene sets and generated a website GUI where users can access this tool. Since a common drawback to ontology enrichment analyses is its verbosity, we developed a very simple filtering algorithm to reduce the ontology size by an order of magnitude. We adjusted these filters and validated our tool using a set of 30 gold standards from Expression Cluster data in WormBase. We show our tool can even discriminate between embryonic and larval tissues and can even identify tissues down to the single-cell level. We used our tool to identify multiple neuronal tissues that are down-regulated due to pathogen infection in C. elegans.

Conclusions

Our Tissue Enrichment Analysis (TEA) can be found within WormBase, and can be downloaded using Python’s standard pip installer. It tests a slimmed-down C. elegans tissue ontology for enrichment of specific terms and provides users with a text and graphic representation of the results.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1229-9) contains supplementary material, which is available to authorized users.

Analysis of C. elegans muscle transcriptome using trans-splicing-based RNA tagging (SRT)

Current approaches to profiling tissue-specific gene expression in C. elegans require delicate manipulation and are difficult under certain conditions, e.g. from dauer or aging worms. We have developed an easy and robust method for tissue-specific RNA-seq by taking advantage of the endogenous trans-splicing process. In this method, transgenic worms are generated in which a spliced leader (SL) RNA gene is fused with a sequence tag and driven by a tissue-specific promoter. Only in the tissue of interest, the tagged SL RNA gene is transcribed and then trans-spliced onto mRNAs. The tag allows enrichment and sequencing of mRNAs from that tissue only. As a proof of principle, we profiled the muscle transcriptome, which showed high coverage and efficient enrichment of muscle specific genes, with low background noise. To demonstrate the robustness of our method, we profiled muscle gene expression in dauer larvae and aging worms, revealing gene expression changes consistent with the physiology of these stages. The resulting muscle transcriptome also revealed 461 novel RNA transcripts, likely muscle-expressed long non-coding RNAs. In summary, the splicing-based RNA tagging (SRT) method provides a convenient and robust tool to profile trans-spliced genes and identify novel transcripts in a tissue-specific manner, with a low false positive rate.

Robust Distal Tip Cell Pathfinding in the Face of Temperature Stress Is Ensured by Two Conserved microRNAS in Caenorhabditis elegans

Biological robustness, the ability of an organism to maintain a steady-state output as genetic or environmental inputs change, is critical for proper development. MicroRNAs have been implicated in biological robustness mechanisms through their post-transcriptional regulation of genes and gene networks. Previous research has illustrated examples of microRNAs promoting robustness as part of feedback loops and genetic switches and by buffering noisy gene expression resulting from environmental and/or internal changes. Here we show that the evolutionarily conserved microRNAs mir-34 and mir-83 (homolog of mammalian mir-29) contribute to the robust migration pattern of the distal tip cells in Caenorhabditis elegans by specifically protecting against stress from temperature changes. Furthermore, our results indicate that mir-34 and mir-83 may modulate the integrin signaling involved in distal tip cell migration by potentially targeting the GTPase cdc-42 and the beta-integrin pat-3. Our findings suggest a role for mir-34 and mir-83 in integrin-controlled cell migrations that may be conserved through higher organisms. They also provide yet another example of microRNA-based developmental robustness in response to a specific environmental stress, rapid temperature fluctuations.

Identification of Wnt Pathway Target Genes Regulating the Division and Differentiation of Larval Seam Cells and Vulval Precursor Cells in Caenorhabditis elegans

The evolutionarily conserved Wnt/β-catenin signaling pathway plays a fundamental role during metazoan development, regulating numerous processes including cell fate specification, cell migration, and stem cell renewal. Wnt ligand binding leads to stabilization of the transcriptional effector β-catenin and upregulation of target gene expression to mediate a cellular response. During larval development of the nematode Caenorhabditis elegans, Wnt/β-catenin pathways act in fate specification of two hypodermal cell types, the ventral vulval precursor cells (VPCs) and the lateral seam cells. Because little is known about targets of the Wnt signaling pathways acting during larval VPC and seam cell differentiation, we sought to identify genes regulated by Wnt signaling in these two hypodermal cell types. We conditionally activated Wnt signaling in larval animals and performed cell type-specific "mRNA tagging" to enrich for VPC and seam cell-specific mRNAs, and then used microarray analysis to examine gene expression compared to control animals. Two hundred thirty-nine genes activated in response to Wnt signaling were identified, and we characterized 50 genes further. The majority of these genes are expressed in seam and/or vulval lineages during normal development, and reduction of function for nine genes caused defects in the proper division, fate specification, fate execution, or differentiation of seam cells and vulval cells. Therefore, the combination of these techniques was successful at identifying potential cell type-specific Wnt pathway target genes from a small number of cells and at increasing our knowledge of the specification and behavior of these C. elegans larval hypodermal cells.

Regulation of structure and function of sarcomeric actin filaments in striated muscle of the nematode Caenorhabditis elegans

The nematode Caenorhabditis elegans has been used as a valuable system to study structure and function of striated muscle. The body wall muscle of C. elegans is obliquely striated muscle with highly organized sarcomeric assembly of actin, myosin, and other accessory proteins. Genetic and molecular biological studies in C. elegans have identified a number of genes encoding structural and regulatory components for the muscle contractile apparatuses, and many of them have counterparts in mammalian cardiac and skeletal muscles or striated muscles in other invertebrates. Applicability of genetics, cell biology, and biochemistry has made C. elegans an excellent system to study mechanisms of muscle contractility and assembly and maintenance of myofibrils. This review focuses on the regulatory mechanisms of structure and function of actin filaments in the C. elegans body wall muscle. Sarcomeric actin filaments in C. elegans muscle are associated with the troponin-tropomyosin system that regulates the actin-myosin interaction. Proteins that bind to the side and ends of actin filaments support ordered assembly of thin filaments. Furthermore, regulators of actin dynamics play important roles in initial assembly, growth, and maintenance of sarcomeres. The knowledge acquired in C. elegans can serve as bases to understand the basic mechanisms of muscle structure and function.

Comparative RNA-Seq analysis reveals pervasive tissue-specific alternative polyadenylation in Caenorhabditis elegans intestine and muscles

BACKGROUND

Tissue-specific RNA plasticity broadly impacts the development, tissue identity and adaptability of all organisms, but changes in composition, expression levels and its impact on gene regulation in different somatic tissues are largely unknown. Here we developed a new method, polyA-tagging and sequencing (PAT-Seq) to isolate high-quality tissue-specific mRNA from Caenorhabditis elegans intestine, pharynx and body muscle tissues and study changes in their tissue-specific transcriptomes and 3'UTRomes.

RESULTS

We have identified thousands of novel genes and isoforms differentially expressed between these three tissues. The intestine transcriptome is expansive, expressing over 30% of C. elegans mRNAs, while muscle transcriptomes are smaller but contain characteristic unique gene signatures. Active promoter regions in all three tissues reveal both known and novel enriched tissue-specific elements, along with putative transcription factors, suggesting novel tissue-specific modes of transcription initiation. We have precisely mapped approximately 20,000 tissue-specific polyadenylation sites and discovered that about 30% of transcripts in somatic cells use alternative polyadenylation in a tissue-specific manner, with their 3'UTR isoforms significantly enriched with microRNA targets.

CONCLUSIONS

For the first time, PAT-Seq allowed us to directly study tissue specific gene expression changes in an in vivo setting and compare these changes between three somatic tissues from the same organism at single-base resolution within the same experiment. We pinpoint precise tissue-specific transcriptome rearrangements and for the first time link tissue-specific alternative polyadenylation to miRNA regulation, suggesting novel and unexplored tissue-specific post-transcriptional regulatory networks in somatic cells.

Cell type-specific affinity purification of nuclei for chromatin profiling in whole animals

Analyzing cell differentiation during development in a complex organism requires the analysis of expression and chromatin profiles in individual cell types. Our laboratory has developed a simple and generally applicable strategy to purify specific cell types from whole organisms for simultaneous analysis of chromatin and expression. The method, termed INTACT for Isolation of Nuclei TAgged in specific Cell Types, depends on the expression of an affinity-tagged nuclear envelope protein in the cell type of interest. These nuclei can be affinity-purified from the total pool of nuclei and used as a source for RNA and chromatin. The method serves as a simple and scalable alternative to FACS sorting or laser capture microscopy to circumvent the need for expensive equipment and specialized skills. This chapter provides detailed protocols for the cell-type specific purification of nuclei from Caenorhabditis elegans.

Spatiotemporal transcriptomics reveals the evolutionary history of the endoderm germ layer

The concept of germ layers has been one of the foremost organizing principles in developmental biology, classification, systematics and evolution for 150 years (refs 1 - 3). Of the three germ layers, the mesoderm is found in bilaterian animals but is absent in species in the phyla Cnidaria and Ctenophora, which has been taken as evidence that the mesoderm was the final germ layer to evolve. The origin of the ectoderm and endoderm germ layers, however, remains unclear, with models supporting the antecedence of each as well as a simultaneous origin. Here we determine the temporal and spatial components of gene expression spanning embryonic development for all Caenorhabditis elegans genes and use it to determine the evolutionary ages of the germ layers. The gene expression program of the mesoderm is induced after those of the ectoderm and endoderm, thus making it the last germ layer both to evolve and to develop. Strikingly, the C. elegans endoderm and ectoderm expression programs do not co-induce; rather the endoderm activates earlier, and this is also observed in the expression of endoderm orthologues during the embryology of the frog Xenopus tropicalis, the sea anemone Nematostella vectensis and the sponge Amphimedon queenslandica. Querying the phylogenetic ages of specifically expressed genes reveals that the endoderm comprises older genes. Taken together, we propose that the endoderm program dates back to the origin of multicellularity, whereas the ectoderm originated as a secondary germ layer freed from ancestral feeding functions.

Integrative analysis of C. elegans modENCODE ChIP-seq data sets to infer gene regulatory interactions

The C. elegans modENCODE Consortium has defined in vivo binding sites for a large array of transcription factors by ChIP-seq. In this article, we present examples that illustrate how this compendium of ChIP-seq data can drive biological insights not possible with analysis of individual factors. First, we analyze the number of independent factors bound to the same locus, termed transcription factor complexity, and find that low-complexity sites are more likely to respond to altered expression of a single bound transcription factor. Next, we show that comparison of binding sites for the same factor across developmental stages can reveal insight into the regulatory network of that factor, as we find that the transcription factor UNC-62 has distinct binding profiles at different stages due to distinct cofactor co-association as well as tissue-specific alternative splicing. Finally, we describe an approach to infer potential regulators of gene expression changes found in profiling experiments (such as DNA microarrays) by screening these altered genes to identify significant enrichment for targets of a transcription factor identified in ChIP-seq data sets. After confirming that this approach can correctly identify the upstream regulator on expression data sets for which the regulator was previously known, we applied this approach to identify novel candidate regulators of transcriptional changes with age. The analysis revealed nine candidate aging regulators, of which three were previously known to have a role in longevity. We experimentally showed that two of the new candidate aging regulators can extend lifespan when overexpressed, indicating that this approach can identify novel functional regulators of complex processes.

Regulatory role of the second gelsolin-like domain of Caenorhabditis elegans gelsolin-like protein 1 (GSNL-1) in its calcium-dependent conformation and actin-regulatory activities

Caenorhabditis elegans gelsolin-like protein-1 (GSNL-1) is an unconventional member of the gelsolin family of actin-regulatory proteins. Unlike typical gelsolin-related proteins with three or six G domains, GSNL-1 has four gelsolin-like (G) domains (G1-G4) and exhibits calcium-dependent actin filament severing and capping activities. The first G domain (G1) of GSNL-1 is necessary for its actin-regulatory activities. However, how other domains in GSNL-1 participate in regulation of its functions is not understood. Here, we report biochemical evidence that the second G domain (G2) of GSNL-1 has a regulatory role in its calcium-dependent conformation and actin-regulatory activities. Comparison of the sequences of gelsolin-related proteins from various species indicates that sequences of G2 are highly conserved. Among the conserved residues in G2, we focused on D162 of GSNL-1, since equivalent residues in gelsolin and severin are part of the calcium-binding sites and is a pathogenic mutation site in human gelsolin causing familial amyloidosis, Finish-type. The D162N mutation does not alter the inactive and fully calcium-activated states of GSNL-1 for actin filament severing (at 20 nM GSNL-1) and capping activities (at 50 nM GSNL-1). However, under these conditions, the mutant shows reduced calcium sensitivity for activation. By contrast, the D162N mutation strongly enhances susceptibility of GSNL-1 to chymotrypsin digestion only at high calcium concentrations but not at low calcium concentrations. The mutation also reduces affinity of GSNL-1 with actin monomers. These results suggest that G2 of GSNL-1 functions as a regulatory domain for its calcium-dependent actin-regulatory activities by mediating conformational changes of the GSNL-1 molecule.

Identifying functional links between genes by evolutionary transcriptomics

The ability to determine gene expression profiles across distant species presents a unique opportunity to identify functional relationships between genes. In particular, transcriptome data may help to distinguish whether genes with similar expression profiles are functionally related or independent. Recent studies on the evolution of gene expression have revealed a striking amount of divergence across strains and species, a notion which has hitherto not been brought to bear on the problem of detecting functional relationships between genes. Here, we introduce evo-links, a method by which a pair of genes are linked if their expression profiles are consistently more similar within species, while their individual conservation across species is low. We show that genes connected through evo-links are more enriched in known functional interactions than genes linked by conventional correlation measures. The network of linked genes further allows the identification of gene communities which reflect distinct functional pathways. We classified communities into major cell-types and derived a temporal developmental map of tissue specification in the nematode C. elegans. This map shows the sequential activation of the endoderm, body wall muscle, and neuronal tissues, and later the pharynx. We propose that as comparative transcriptomics becomes increasingly feasible, evo-links offer a robust method to detect functional relationships and disentangle developmental pathways in data lacking spatial resolution.

The atypical calpains: evolutionary analyses and roles in Caenorhabditis elegans cellular degeneration

The calpains are physiologically important Ca²⁺-activated regulatory proteases, which are divided into typical or atypical sub-families based on constituent domains. Both sub-families are present in mammals, but our understanding of calpain function is based primarily on typical sub-family members. Here, we take advantage of the model organism Caenorhabditis elegans, which expresses only atypical calpains, to extend our knowledge of the phylogenetic evolution and function of calpains. We provide evidence that a typical human calpain protein with a penta EF hand, detected using custom profile hidden Markov models, is conserved in ancient metazoans and a divergent clade. These analyses also provide evidence for the lineage-specific loss of typical calpain genes in C. elegans and Ciona, and they reveal that many calpain-like genes lack an intact catalytic triad. Given the association between the dysregulation of typical calpains and human degenerative pathologies, we explored the phenotypes, expression profiles, and consequences of inappropriate reduction or activation of C. elegans atypical calpains. These studies show that the atypical calpain gene, clp-1, contributes to muscle degeneration and reveal that clp-1 activity is sensitive to genetic manipulation of [Ca²⁺]_i. We show that CLP-1 localizes to sarcomeric sub-structures, but is excluded from dense bodies (Z-disks). We find that the muscle degeneration observed in a C. elegans model of dystrophin-based muscular dystrophy can be suppressed by clp-1 inactivation and that nemadipine-A inhibition of the EGL-19 calcium channel reveals that Ca²⁺ dysfunction underlies the C. elegans MyoD model of myopathy. Taken together, our analyses highlight the roles of calcium dysregulation and CLP-1 in muscle myopathies and suggest that the atypical calpains could retain conserved roles in myofilament turnover.

Calpains are calcium activated non-lysosomal proteases that cleave proteins with exquisite selectivity. Proteins can be activated by calpain cleavage, because they are released from inhibitory constraints, or they can be targeted for further degradation to facilitate their normal physiological turnover or to promote cellular remodelling. Inappropriate calpain activity can lead to degenerative pathologies and cancers. Our understanding of calpain function is based primarily on typical calpains, which carry EF hand motifs that bind Ca²⁺ or mediate dimerization; however, typical and atypical calpains, which lack EF hand motifs, are both present in mammals. Hence, any therapeutic intervention designed to suppress degenerative conditions, particularly those caused by elevated Ca²⁺ levels, should also consider the potential involvement of atypical calpains. We have taken advantage of the model organism C. elegans, which only encodes atypical calpain proteins, to gain an understanding of the evolution and activities of these proteins. We show that the CLP-1 atypical calpain is normally expressed in muscle and localizes to sarcomeric sub-structures. We find that CLP-1 contributes to the muscle degeneration observed in a model of Duchenne muscular dystrophy. Our studies also highlight the importance of calcium dysregulation in promoting CLP-1 activity and muscle degeneration.

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.

Calpains mediate integrin attachment complex maintenance of adult muscle in Caenorhabditis elegans

Two components of integrin containing attachment complexes, UNC-97/PINCH and UNC-112/MIG-2/Kindlin-2, were recently identified as negative regulators of muscle protein degradation and as having decreased mRNA levels in response to spaceflight. Integrin complexes transmit force between the inside and outside of muscle cells and signal changes in muscle size in response to force and, perhaps, disuse. We therefore investigated the effects of acute decreases in expression of the genes encoding these multi-protein complexes. We find that in fully developed adult Caenorhabditis elegans muscle, RNAi against genes encoding core, and peripheral, members of these complexes induces protein degradation, myofibrillar and mitochondrial dystrophies, and a movement defect. Genetic disruption of Z-line– or M-line–specific complex members is sufficient to induce these defects. We confirmed that defects occur in temperature-sensitive mutants for two of the genes: unc-52, which encodes the extra-cellular ligand Perlecan, and unc-112, which encodes the intracellular component Kindlin-2. These results demonstrate that integrin containing attachment complexes, as a whole, are required for proper maintenance of adult muscle. These defects, and collapse of arrayed attachment complexes into ball like structures, are blocked when DIM-1 levels are reduced. Degradation is also blocked by RNAi or drugs targeting calpains, implying that disruption of integrin containing complexes results in calpain activation. In wild-type animals, either during development or in adults, RNAi against calpain genes results in integrin muscle attachment disruptions and consequent sub-cellular defects. These results demonstrate that calpains are required for proper assembly and maintenance of integrin attachment complexes. Taken together our data provide in vivo evidence that a calpain-based molecular repair mechanism exists for dealing with attachment complex disruption in adult muscle. Since C. elegans lacks satellite cells, this mechanism is intrinsic to the muscles and raises the question if such a mechanism also exists in higher metazoans.

Muscle is a dynamic tissue that grows in response to use and nutrition and shrinks in response to lack of use, poor nutrition, or disease. Loss of muscle mass is an important public health problem, but we understand little of the genes that regulate muscle shrinkage. We have found that, in adult worm muscle, attachment to the basement membrane is continuously required to prevent catastrophic sub-cellular defects that result in impaired ability of muscle to function. We have also identified a group of proteases that are activated when the attachment fails to be properly maintained. Conversely, when these proteases are lacking in adult muscle, the muscles fail to maintain attachment to the basement membrane. Thus, we have discovered a group of proteases that appear to act to maintain attachment to the basement membrane and therefore to maintain muscle itself. Because these worms lack satellite cells, this maintenance system is intrinsic to muscle, thus raising the question whether a similar or identical system also works in humans.

Cell-type-specific nuclei purification from whole animals for genome-wide expression and chromatin profiling

An understanding of developmental processes requires knowledge of transcriptional and epigenetic landscapes at the level of tissues and ultimately individual cells. However, obtaining tissue- or cell-type-specific expression and chromatin profiles for animals has been challenging. Here we describe a method for purifying nuclei from specific cell types of animal models that allows simultaneous determination of both expression and chromatin profiles. The method is based on in vivo biotin-labeling of the nuclear envelope and subsequent affinity purification of nuclei. We describe the use of the method to isolate nuclei from muscle of adult Caenorhabditis elegans and from mesoderm of Drosophila melanogaster embryos. As a case study, we determined expression and nucleosome occupancy profiles for affinity-purified nuclei from C. elegans muscle. We identified hundreds of genes that are specifically expressed in muscle tissues and found that these genes are depleted of nucleosomes at promoters and gene bodies in muscle relative to other tissues. This method should be universally applicable to all model systems that allow transgenesis and will make it possible to determine epigenetic and expression profiles of different tissues and cell types.

Caenorhabditis elegans cyclin D/CDK4 and cyclin E/CDK2 induce distinct cell cycle re-entry programs in differentiated muscle cells

Cell proliferation and differentiation are regulated in a highly coordinated and inverse manner during development and tissue homeostasis. Terminal differentiation usually coincides with cell cycle exit and is thought to engage stable transcriptional repression of cell cycle genes. Here, we examine the robustness of the post-mitotic state, using Caenorhabditis elegans muscle cells as a model. We found that expression of a G1 Cyclin and CDK initiates cell cycle re-entry in muscle cells without interfering with the differentiated state. Cyclin D/CDK4 (CYD-1/CDK-4) expression was sufficient to induce DNA synthesis in muscle cells, in contrast to Cyclin E/CDK2 (CYE-1/CDK-2), which triggered mitotic events. Tissue-specific gene-expression profiling and single molecule FISH experiments revealed that Cyclin D and E kinases activate an extensive and overlapping set of cell cycle genes in muscle, yet failed to induce some key activators of G1/S progression. Surprisingly, CYD-1/CDK-4 also induced an additional set of genes primarily associated with growth and metabolism, which were not activated by CYE-1/CDK-2. Moreover, CYD-1/CDK-4 expression also down-regulated a large number of genes enriched for catabolic functions. These results highlight distinct functions for the two G1 Cyclin/CDK complexes and reveal a previously unknown activity of Cyclin D/CDK-4 in regulating metabolic gene expression. Furthermore, our data demonstrate that many cell cycle genes can still be transcriptionally induced in post-mitotic muscle cells, while maintenance of the post-mitotic state might depend on stable repression of a limited number of critical cell cycle regulators.

During development, cells face the important decision whether to continue to proliferate, or to exit the cell-division cycle and fully differentiate. Improved insight into the molecular mechanisms that arrest the cell cycle during terminal differentiation is important for our understanding of normal development, as well as for cancer research and regenerative medicine. To investigate the arrested state of terminally differentiated cells, we examined muscle cells in the model organism C. elegans, which is known for its reproducible cell-division pattern. We found that expression of a single cell cycle kinase with its regulatory partner (Cyclin) induced many cell division genes in muscle. While Cyclin D and E kinases often act similarly, only Cyclin D with CDK-4 triggered DNA replication in muscle, and this combination induced a much broader transcriptional response than Cyclin E/CDK-2. Despite activation of a substantial cell cycle program, Cyclin/CDK expression did not induce complete muscle cell division and failed to induce some key cell cycle regulators. Our results highlight distinct activities of Cyclin D and Cyclin E kinases, and they indicate that cell-cycle gene expression remains remarkably flexible in differentiated cells. We propose that the post-mitotic state of differentiated cells is maintained by tight control of a few regulatory genes.

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.

Biochemical and cell biological analysis of actin in the nematode Caenorhabditis elegans

The nematode Caenorhabditis elegans has long been a useful model organism for muscle research. Its body wall muscle is obliquely striated muscle and exhibits structural similarities with vertebrate striated muscle. Actin is the core component of the muscle thin filaments, which are highly ordered in sarcomeric structures in striated muscle. Genetic studies have identified genes that regulate proper organization and function of actin filaments in C. elegans muscle, and sequence of the worm genome has revealed a number of conserved candidate genes that may regulate actin. To precisely understand the functions of actin-binding proteins, such genetic and genomic studies need to be complemented by biochemical characterization of these actin-binding proteins in vitro. This article describes methods for purification and biochemical characterization of actin from C. elegans. Although rabbit muscle actin is commonly used to characterize actin-binding proteins from many eukaryotic organisms, we detect several quantitative differences between C. elegans actin and rabbit muscle actin, highlighting that use of actin from an appropriate source is important in some cases. Additionally, we describe probes for cell biological analysis of actin in C. elegans.

Determining the sub-cellular localization of proteins within Caenorhabditis elegans body wall muscle

Determining the sub-cellular localization of a protein within a cell is often an essential step towards understanding its function. In Caenorhabditis elegans, the relatively large size of the body wall muscle cells and the exquisite organization of their sarcomeres offer an opportunity to identify the precise position of proteins within cell substructures. Our goal in this study is to generate a comprehensive “localizome” for C. elegans body wall muscle by GFP-tagging proteins expressed in muscle and determining their location within the cell. For this project, we focused on proteins that we know are expressed in muscle and are orthologs or at least homologs of human proteins. To date we have analyzed the expression of about 227 GFP-tagged proteins that show localized expression in the body wall muscle of this nematode (e.g. dense bodies, M-lines, myofilaments, mitochondria, cell membrane, nucleus or nucleolus). For most proteins analyzed in this study no prior data on sub-cellular localization was available. In addition to discrete sub-cellular localization we observe overlapping patterns of localization including the presence of a protein in the dense body and the nucleus, or the dense body and the M-lines. In total we discern more than 14 sub-cellular localization patterns within nematode body wall muscle. The localization of this large set of proteins within a muscle cell will serve as an invaluable resource in our investigation of muscle sarcomere assembly and function.

Isolation and culture of larval cells from C. elegans

Cell culture is an essential tool to study cell function. In C. elegans the ability to isolate and culture cells has been limited to embryonically derived cells. However, cells or blastomeres isolated from mixed stage embryos terminally differentiate within 24 hours of culture, thus precluding post-embryonic stage cell culture. We have developed an efficient and technically simple method for large-scale isolation and primary culture of larval-stage cells. We have optimized the treatment to maximize cell number and minimize cell death for each of the four larval stages. We obtained up to 7.8×10(4) cells per microliter of packed larvae, and up to 97% of adherent cells isolated by this method were viable for at least 16 hours. Cultured larval cells showed stage-specific increases in both cell size and multinuclearity and expressed lineage- and cell type-specific reporters. The majority (81%) of larval cells isolated by our method were muscle cells that exhibited stage-specific phenotypes. L1 muscle cells developed 1 to 2 wide cytoplasmic processes, while L4 muscle cells developed 4 to 14 processes of various thicknesses. L4 muscle cells developed bands of myosin heavy chain A thick filaments at the cell center and spontaneously contracted ex vivo. Neurons constituted less than 10% of the isolated cells and the majority of neurons developed one or more long, microtubule-rich protrusions that terminated in actin-rich growth cones. In addition to cells such as muscle and neuron that are high abundance in vivo, we were also able to isolate M-lineage cells that constitute less than 0.2% of cells in vivo. Our novel method of cell isolation extends C. elegans cell culture to larval developmental stages, and allows use of the wealth of cell culture tools, such as cell sorting, electrophysiology, co-culture, and high-resolution imaging of subcellular dynamics, in investigation of post-embryonic development and physiology.

A spatial and temporal map of C. elegans gene expression

The C. elegans genome has been completely sequenced, and the developmental anatomy of this model organism is described at single-cell resolution. Here we utilize strategies that exploit this precisely defined architecture to link gene expression to cell type. We obtained RNAs from specific cells and from each developmental stage using tissue-specific promoters to mark cells for isolation by FACS or for mRNA extraction by the mRNA-tagging method. We then generated gene expression profiles of more than 30 different cells and developmental stages using tiling arrays. Machine-learning-based analysis detected transcripts corresponding to established gene models and revealed novel transcriptionally active regions (TARs) in noncoding domains that comprise at least 10% of the total C. elegans genome. Our results show that about 75% of transcripts with detectable expression are differentially expressed among developmental stages and across cell types. Examination of known tissue- and cell-specific transcripts validates these data sets and suggests that newly identified TARs may exercise cell-specific functions. Additionally, we used self-organizing maps to define groups of coregulated transcripts and applied regulatory element analysis to identify known transcription factor- and miRNA-binding sites, as well as novel motifs that likely function to control subsets of these genes. By using cell-specific, whole-genome profiling strategies, we have detected a large number of novel transcripts and produced high-resolution gene expression maps that provide a basis for establishing the roles of individual genes in cellular differentiation.

A widespread distribution of genomic CeMyoD binding sites revealed and cross validated by ChIP-Chip and ChIP-Seq techniques

Identifying transcription factor binding sites genome-wide using chromatin immunoprecipitation (ChIP)-based technology is becoming an increasingly important tool in addressing developmental questions. However, technical problems associated with factor abundance and suitable ChIP reagents are common obstacles to these studies in many biological systems. We have used two completely different, widely applicable methods to determine by ChIP the genome-wide binding sites of the master myogenic regulatory transcription factor HLH-1 (CeMyoD) in C. elegans embryos. The two approaches, ChIP-seq and ChIP-chip, yield strongly overlapping results revealing that HLH-1 preferentially binds to promoter regions of genes enriched for E-box sequences (CANNTG), known binding sites for this well-studied class of transcription factors. HLH-1 binding sites were enriched upstream of genes known to be expressed in muscle, consistent with its role as a direct transcriptional regulator. HLH-1 binding was also detected at numerous sites unassociated with muscle gene expression, as has been previously described for its mouse homolog MyoD. These binding sites may reflect several additional functions for HLH-1, including its interactions with one or more co-factors to activate (or repress) gene expression or a role in chromatin organization distinct from direct transcriptional regulation of target genes. Our results also provide a comparison of ChIP methodologies that can overcome limitations commonly encountered in these types of studies while highlighting the complications of assigning in vivo functions to identified target sites.

Dynamic regulation of sarcomeric actin filaments in striated muscle

In striated muscle, the actin cytoskeleton is differentiated into myofibrils. Actin and myosin filaments are organized in sarcomeres and specialized for producing contractile forces. Regular arrangement of actin filaments with uniform length and polarity is critical for the contractile function. However, the mechanisms of assembly and maintenance of sarcomeric actin filaments in striated muscle are not completely understood. Live imaging of actin in striated muscle has revealed that actin subunits within sarcomeric actin filaments are dynamically exchanged without altering overall sarcomeric structures. A number of regulators for actin dynamics have been identified, and malfunction of these regulators often result in disorganization of myofibril structures or muscle diseases. Therefore, proper regulation of actin dynamics in striated muscle is critical for assembly and maintenance of functional myofibrils. Recent studies have suggested that both enhancers of actin dynamics and stabilizers of actin filaments are important for sarcomeric actin organization. Further investigation of the regulatory mechanism of actin dynamics in striated muscle should be a key to understanding how myofibrils develop and operate.

C. elegans twist gene expression in differentiated cell types is controlled by autoregulation through intron elements

The temporospatial regulation of genes encoding transcription factors is important during development. The hlh-8 gene encodes the C. elegans mesodermal transcription factor CeTwist. Elements in the hlh-8 promoter restrict gene expression to predominantly undifferentiated cells of the M lineage. We have discovered that hlh-8 expression in differentiated mesodermal cells is controlled by two well-conserved E box elements in the large first intron. Additionally, we found that these elements are bound in vitro by CeTwist and its transcription factor partner, CeE/DA. The E box driven expression is eliminated or diminished in an hlh-8 null allele or in hlh-2 (CeE/DA) RNAi, respectively. Expression of hlh-8 is also diminished in animals harboring an hlh-8 intron deletion allele. Altogether, our results support a model in which hlh-8 is initially expressed in the undifferentiated M lineage cells via promoter elements and then the CeTwist activates its own expression further (autoregulation) in differentiated cells derived from the M lineage via the intron elements. This model provides a mechanism for how a transcription factor may regulate distinct target genes in cells both before and after initiating the differentiation program. The findings could also be relevant to understanding human Twist gene regulation, which is currently not well understood.

Recent advances in understanding the molecular mechanisms regulating C. elegans transcription

We review recent studies that have advanced our understanding of the molecular mechanisms regulating transcription in the nematode C. elegans. Topics covered include: (i) general properties of C. elegans promoters; (ii) transcription factors and transcription factor combinations involved in cell fate specification and cell differentiation; (iii) new roles for general transcription factors; (iv) nucleosome positioning in C. elegans "chromatin"; and (v) some characteristics of histone variants and histone modifications and their possible roles in controlling C. elegans transcription.

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.

Time-lapse imaging and cell-specific expression profiling reveal dynamic branching and molecular determinants of a multi-dendritic nociceptor in C. elegans

Nociceptive neurons innervate the skin with complex dendritic arbors that respond to pain-evoking stimuli such as harsh mechanical force or extreme temperatures. Here we describe the structure and development of a model nociceptor, the PVD neuron of C. elegans, and identify transcription factors that control morphogenesis of the PVD dendritic arbor. The two PVD neuron cell bodies occupy positions on either the right (PVDR) or left (PVDL) sides of the animal in posterior-lateral locations. Imaging with a GFP reporter revealed a single axon projecting from the PVD soma to the ventral cord and an elaborate, highly branched arbor of dendritic processes that envelop the animal with a web-like array directly beneath the skin. Dendritic branches emerge in a step-wise fashion during larval development and may use an existing network of peripheral nerve cords as guideposts for key branching decisions. Time-lapse imaging revealed that branching is highly dynamic with active extension and withdrawal and that PVD branch overlap is prevented by a contact-dependent self-avoidance, a mechanism that is also employed by sensory neurons in other organisms. With the goal of identifying genes that regulate dendritic morphogenesis, we used the mRNA-tagging method to produce a gene expression profile of PVD during late larval development. This microarray experiment identified>2,000 genes that are 1.5X elevated relative to all larval cells. The enriched transcripts encode a wide range of proteins with potential roles in PVD function (e.g., DEG/ENaC and Trp channels) or development (e.g., UNC-5 and LIN-17/frizzled receptors). We used RNAi and genetic tests to screen 86 transcription factors from this list and identified eleven genes that specify PVD dendritic structure. These transcription factors appear to control discrete steps in PVD morphogenesis and may either promote or limit PVD branching at specific developmental stages. For example, time-lapse imaging revealed that MEC-3 (LIM homeodomain) is required for branch initiation in early larval development whereas EGL-44 (TEAD domain) prevents ectopic PVD branching in the adult. A comparison of PVD-enriched transcripts to a microarray profile of mammalian nociceptors revealed homologous genes with potentially shared nociceptive functions. We conclude that PVD neurons display striking structural, functional and molecular similarities to nociceptive neurons from more complex organisms and can thus provide a useful model system in which to identify evolutionarily conserved determinants of nociceptor fate.

Distinct roles of four gelsolin-like domains of Caenorhabditis elegans gelsolin-like protein-1 in actin filament severing, barbed end capping, and phosphoinositide binding

Caenorhabditis elegans gelsolin-like protein-1 (GSNL-1) is a new member of the gelsolin family of actin regulatory proteins [Klaavuniemi, T., Yamashiro, S., and Ono, S. (2008) J. Biol. Chem. 283, 26071-26080]. It is an unconventional gelsolin-related protein with four gelsolin-like (G) domains (G1-G4), unlike typical gelsolin-related proteins with three or six G domains. GSNL-1 severs actin filaments and caps the barbed end in a calcium-dependent manner similar to that of gelsolin. In contrast, GSNL-1 has properties different from those of gelsolin in that it remains bound to F-actin and does not nucleate actin polymerization. To understand the mechanism by which GSNL-1 regulates actin dynamics, we investigated the domain-function relationship of GSNL-1 by analyzing activities of truncated forms of GSNL-1. G1 and the linker between G1 and G2 were sufficient for actin filament severing, whereas G1 and G2 were required for barbed end capping. The actin severing activity of GSNL-1 was inhibited by phosphatidylinositol 4,5-bisphosphate (PIP2), and a PIP2-sensitive domain was mapped to G1 and G2. At least two actin-binding sites were detected: a calcium-dependent G-actin-binding site in G1 and a calcium-independent G- and F-actin-binding site in G3 and G4. These results reveal both conserved and different utilization of G domains between C. elegans GSNL-1 and mammalian gelsolin for actin regulatory functions.