The embryonic muscle transcriptome of Caenorhabditis elegans

Determination of the mobility of novel and established Caenorhabditis elegans sarcomeric proteins in vivo

A screen was instigated to identify novel protein components of the Caenorhabditis elegans sarcomere. The subcellular localisation of full-length GFP fusion proteins was examined, in transgenic animals, for 62 essentially uncharacterized genes thought to be expressed within bodywall muscle cells. Three genes, T03G6.3, C46G7.2 and K04A8.6, were identified for further study. K04A8.6::GFP only displayed a regular sarcomeric distribution sporadically. However, C46G7.2::GFP localised to the centre of A-bands and dense bodies and T03G6.3::GFP localised in the I-band, of the bodywall muscle sarcomeres, consistently. This success with such a small screen suggests that there are further minor components of the C. elegans sarcomere yet to be discovered. Fluorescence Recovery After Photobleaching (FRAP) was applied to live transgenic individuals to assess the mobility of T03G6.3 and C46G7.2 and other well-known constituents of the sarcomere in vivo. Proteins associated with the thin filaments showed dynamic exchange whilst those associated with thick filaments appeared more static. This is the first demonstration that there are sarcomeric proteins in C. elegans muscle cells in dynamic exchange and that the rates of exchange in vivo correspond in general terms with observations in other experimental systems.

GExplore: a web server for integrated queries of protein domains, gene expression and mutant phenotypes

Background

The majority of the genes even in well-studied multi-cellular model organisms have not been functionally characterized yet. Mining the numerous genome wide data sets related to protein function to retrieve potential candidate genes for a particular biological process remains a challenge.

Description

GExplore has been developed to provide a user-friendly database interface for data mining at the gene expression/protein function level to help in hypothesis development and experiment design. It supports combinatorial searches for proteins with certain domains, tissue- or developmental stage-specific expression patterns, and mutant phenotypes. GExplore operates on a stand-alone database and has fast response times, which is essential for exploratory searches. The interface is not only user-friendly, but also modular so that it accommodates additional data sets in the future.

Conclusion

GExplore is an online database for quick mining of data related to gene and protein function, providing a multi-gene display of data sets related to the domain composition of proteins as well as expression and phenotype data. GExplore is publicly available at: http://genome.sfu.ca/gexplore/

A universal nonmonotonic relationship between gene compactness and expression levels in multicellular eukaryotes

Analysis of gene architecture and expression levels of four organisms, Homo sapiens, Caenorhabditis elegans, Drosophila melanogaster, and Arabidopsis thaliana, reveals a surprising, nonmonotonic, universal relationship between expression level and gene compactness. With increasing expression level, the genes tend at first to become longer but, from a certain level of expression, they become more and more compact, resulting in an approximate bell-shaped dependence. There are two leading hypotheses to explain the compactness of highly expressed genes. The selection hypothesis predicts that gene compactness is predominantly driven by the level of expression, whereas the genomic design hypothesis predicts that expression breadth across tissues is the driving force. We observed the connection between gene expression breadth in humans and gene compactness to be significantly weaker than the connection between expression level and compactness, a result that is compatible with the selection hypothesis but not the genome design hypothesis. The initial gene elongation with increasing expression level could be explained, at least in part, by accumulation of regulatory elements enhancing expression, in particular, in introns. This explanation is compatible with the observed positive correlation between intron density and expression level of a gene. Conversely, the trend toward increasing compactness for highly expressed genes could be caused by selection for minimization of energy and time expenditure during transcription and splicing and for increased fidelity of transcription, splicing, and/or translation that is likely to be particularly critical for highly expressed genes. Regardless of the exact nature of the forces that shape the gene architecture, we present evidence that, at least, in animals, coding and noncoding parts of genes show similar architectonic trends.

C. elegans dysferlin homolog fer-1 is expressed in muscle, and fer-1 mutations initiate altered gene expression of muscle enriched genes

Mutations in the human dysferlin gene cause Limb Girdle Muscular Dystrophy 2B (LGMD2B). The Caenorhabditis elegans dysferlin homolog, fer-1, affects sperms development but is not known to be expressed in or have a functional roles outside of the male germline. Using several approaches, we show that fer-1 mRNA is present in C. elegans muscle cells but is absent from neurons. In mammals, loss of muscle-expressed dysferlin causes transcriptional deregulation of muscle expressed genes. To determine if similar alterations in gene expression are initiated in C. elegans due to loss of muscle-expressed fer-1, we performed whole genome Affymetrix microarray analysis of two loss-of-function fer-1 mutants. Both mutants gave rise to highly similar changes in gene expression and altered the expression of 337 genes. Using multiple analysis methods, we show that this gene set is enriched for genes known to regulate the structure and function of muscle. However, these transcriptional changes do not appear to be in response to gross sarcomeric damage, since genetically sensitized fer-1 mutants exhibit normal thin filament organization. Our data suggest that processes other than sarcomere stability may be affected by loss of fer-1 in C. elegans muscle. Therefore, C. elegans may be an attractive model system in which to explore new muscle-specific functions of the dysferlin protein and gain insights into the molecular pathogenesis of LGMD2B.

An ER-resident membrane protein complex regulates nicotinic acetylcholine receptor subunit composition at the synapse

Nicotinic acetylcholine receptors (nAChRs) are homo- or heteropentameric ligand-gated ion channels mediating excitatory neurotransmission and muscle activation. Regulation of nAChR subunit assembly and transfer of correctly assembled pentamers to the cell surface is only partially understood. Here, we characterize an ER transmembrane (TM) protein complex that influences nAChR cell-surface expression and functional properties in Caenorhabditis elegans muscle. Loss of either type I TM protein, NRA-2 or NRA-4 (nicotinic receptor associated), affects two different types of muscle nAChRs and causes in vivo resistance to cholinergic agonists. Sensitivity to subtype-specific agonists of these nAChRs is altered differently, as demonstrated by whole-cell voltage-clamp of dissected adult muscle, when applying exogenous agonists or after photo-evoked, channelrhodopsin-2 (ChR2) mediated acetylcholine (ACh) release, as well as in single-channel recordings in cultured embryonic muscle. These data suggest that nAChRs desensitize faster in nra-2 mutants. Cell-surface expression of different subunits of the 'levamisole-sensitive' nAChR (L-AChR) is differentially affected in the absence of NRA-2 or NRA-4, suggesting that they control nAChR subunit composition or allow only certain receptor assemblies to leave the ER.

The complexity of the secreted NPA and FAR lipid-binding protein families of Haemonchus contortus revealed by an iterative proteomics-bioinformatics approach

Two different classes of small nematode specific lipid-binding proteins, the nematode polyprotein allergens/antigens (NPAs) and the fatty acid- and retinol-binding (FAR) proteins, are secreted by helminth parasites. Until now, there was no evidence of the expression or secretion of these two families of proteins in Haemonchus contortus. In this study, we applied proteomic and bioinformatic tools in an iterative manner to reveal the expression and complexity of these proteins in the excretory/secretory products (ESP) of adult H. contortus at the protein and gene levels. Initial examination of the mass spectra of ESP fractions against standard databases returned nine peptides mapping to Ostertagia ostertagi NPA and FAR sequences. Searches of the H. contortus EST and genomic contig databases with the O. ostertagi and Caenorhabditis elegans homologues retrieved diverse sequences encoding H. contortus NPA and FAR proteins. H. contortus sequences were then integrated into a customized database and a new search of the mass spectra achieved a 10-fold improvement in coverage of the predicted H. contortus NPAs. The final analyses of the mass spectra achieved 49-60% coverage of H. contortus NPAs and 7-47% coverage of H. contortus FARs. Moreover, the diversity in structures of the encoding genes was revealed by assembling the genomic sequence data with predicted protein sequences confirmed by the peptide evidence. We predict there are at least one Hc-NPA gene and six Hc-FAR genes in H. contortus, and life stage gene expression of Hc-FAR-1 to -6 revealed unique transcription patterns for each of these genes.

The NK-2 class homeodomain factor CEH-51 and the T-box factor TBX-35 have overlapping function in C. elegans mesoderm development

The C. elegans MS blastomere, born at the 7-cell stage of embryogenesis, generates primarily mesodermal cell types, including pharynx cells, body muscles and coelomocytes. A presumptive null mutation in the T-box factor gene tbx-35, a target of the MED-1 and MED-2 divergent GATA factors, was previously found to result in a profound decrease in the production of MS-derived tissues, although the tbx-35(-) embryonic arrest phenotype was variable. We report here that the NK-2 class homeobox gene ceh-51 is a direct target of TBX-35 and at least one other factor, and that CEH-51 and TBX-35 share functions. Embryos homozygous for a ceh-51 null mutation arrest as larvae with pharynx and muscle defects, although these tissues appear to be specified correctly. Loss of tbx-35 and ceh-51 together results in a synergistic phenotype resembling loss of med-1 and med-2. Overexpression of ceh-51 causes embryonic arrest and generation of ectopic body muscle and coelomocytes. Our data show that TBX-35 and CEH-51 have overlapping function in MS lineage development. As T-box regulators and NK-2 homeodomain factors are both important for heart development in Drosophila and vertebrates, our results suggest that these regulators function in a similar manner in C. elegans to specify a major precursor of mesoderm.

Invertebrate models of age-related muscle degeneration

Functional and structural deterioration of muscles is an inevitable consequence of ageing in a wide variety of animal species. What underlies these changes is a complex network of interactions between the muscle-intrinsic and muscle-extrinsic factors, making it very difficult to distinguish between the cause and the consequence. Many of the genes, structures, and processes implicated in mammalian skeletal muscle ageing are preserved in invertebrate species Drosophila melanogaster and Caenorhabditis elegans. The absence in these organisms of mechanisms that promote muscle regeneration, and substantially different hormonal environment, warrant caution when extrapolating experimental data from studies conducted in invertebrates to mammalian species. The simplicity and accessibility of these models, however, offer ample opportunities for studying age-related myopathologies as well as investigating drugs and therapies to alleviate them.

An integrated strategy to study muscle development and myofilament structure in Caenorhabditis elegans

A crucial step in the development of muscle cells in all metazoan animals is the assembly and anchorage of the sarcomere, the essential repeat unit responsible for muscle contraction. In Caenorhabditis elegans, many of the critical proteins involved in this process have been uncovered through mutational screens focusing on uncoordinated movement and embryonic arrest phenotypes. We propose that additional sarcomeric proteins exist for which there is a less severe, or entirely different, mutant phenotype produced in their absence. We have used Serial Analysis of Gene Expression (SAGE) to generate a comprehensive profile of late embryonic muscle gene expression. We generated two replicate long SAGE libraries for sorted embryonic muscle cells, identifying 7,974 protein-coding genes. A refined list of 3,577 genes expressed in muscle cells was compiled from the overlap between our SAGE data and available microarray data. Using the genes in our refined list, we have performed two separate RNA interference (RNAi) screens to identify novel genes that play a role in sarcomere assembly and/or maintenance in either embryonic or adult muscle. To identify muscle defects in embryos, we screened specifically for the Pat embryonic arrest phenotype. To visualize muscle defects in adult animals, we fed dsRNA to worms producing a GFP-tagged myosin protein, thus allowing us to analyze their myofilament organization under gene knockdown conditions using fluorescence microscopy. By eliminating or severely reducing the expression of 3,300 genes using RNAi, we identified 122 genes necessary for proper myofilament organization, 108 of which are genes without a previously characterized role in muscle. Many of the genes affecting sarcomere integrity have human homologs for which little or nothing is known.

Muscular diseases affect many people worldwide. While we have learned much about the sarcomere, the basic building block of muscle cells, there are still numerous questions that remain to be answered. We must learn more about proteins expressed in muscle and how they interact so that better treatments for myopathies can be developed. The nematode Caenorhabditis elegans is a valuable model organism for the study of muscle due to similarities between worm body wall muscle and vertebrate muscle, along with its semi-transparent cuticle that allows for visualization of muscle structures in live animals. We have used transcriptional profiling methods to identify the majority of genes that are expressed in the embryonic body wall muscle cells of C. elegans. To gain insight into possible functions performed by these genes and their corresponding proteins, we examined animals and muscle cells for abnormalities after the targeted inactivation of about 3,300 genes. We identified 122 genes necessary for proper myofilament organization, 108 of which had no previously characterized role in muscle. This approach proved to be a rapid and sensitive means to identify genes that affect muscle differentiation and sarcomere assembly.

Global prediction of tissue-specific gene expression and context-dependent gene networks in Caenorhabditis elegans

Tissue-specific gene expression plays a fundamental role in metazoan biology and is an important aspect of many complex diseases. Nevertheless, an organism-wide map of tissue-specific expression remains elusive due to difficulty in obtaining these data experimentally. Here, we leveraged existing whole-animal Caenorhabditis elegans microarray data representing diverse conditions and developmental stages to generate accurate predictions of tissue-specific gene expression and experimentally validated these predictions. These patterns of tissue-specific expression are more accurate than existing high-throughput experimental studies for nearly all tissues; they also complement existing experiments by addressing tissue-specific expression present at particular developmental stages and in small tissues. We used these predictions to address several experimentally challenging questions, including the identification of tissue-specific transcriptional motifs and the discovery of potential miRNA regulation specific to particular tissues. We also investigate the role of tissue context in gene function through tissue-specific functional interaction networks. To our knowledge, this is the first study producing high-accuracy predictions of tissue-specific expression and interactions for a metazoan organism based on whole-animal data.

In animals, a crucial facet of any gene's function is the tissue or cell type in which that gene is expressed and the proteins that it interacts with in that cell. However, genome-wide identification of expression across the multitude of tissues of varying size and complexity is difficult to achieve experimentally. In this paper, we show that microararray data collected from whole animals can be analyzed to yield high-quality predictions of tissue-specific expression. These predictions are of better or comparable accuracy to tissue-specific expression determined from high-throughput experiments. Our results provide a global view of tissue-specific expression in Caenorhabditis elegans, allowing us to address the question of how expression patterns are regulated and to analyze how the functions of genes that are expressed in several tissues are influenced by the cellular context.

Structure and evolution of the C. elegans embryonic endomesoderm network

The specification of the Caenorhabditis elegans endomesoderm has been the subject of study for more than 15 years. Specification of the 4-cell stage endomesoderm precursor, EMS, occurs as a result of the activation of a transcription factor cascade that starts with SKN-1, coupled with input from the Wnt/beta-catenin asymmetry pathway through the nuclear effector POP-1. As development proceeds, transiently-expressed cell fate factors are succeeded by stable, tissue/organ-specific regulators. The pathway is complex and uses motifs found in all transcriptional networks. Here, the regulators that function in the C. elegans endomesoderm network are described. An examination of the motifs in the network suggests how they may have evolved from simpler gene interactions. Flexibility in the network is evident from the multitude of parallel functions that have been identified and from apparent changes in parts of the corresponding network in Caenorhabditis briggsae. Overall, the complexities of C. elegans endomesoderm specification build a picture of a network that is robust, complex, and still evolving.

Essential role of ADF/cofilin for assembly of contractile actin networks in the C. elegans somatic gonad

The somatic gonad of the nematode Caenorhabditis elegans contains a myoepithelial sheath, which surrounds oocytes and provides contractile forces during ovulation. Contractile apparatuses of the myoepithelial-sheath cells are non-striated and similar to those of smooth muscle. We report the identification of a specific isoform of actin depolymerizing factor (ADF)/cofilin as an essential factor for assembly of contractile actin networks in the gonadal myoepithelial sheath. Two ADF/cofilin isoforms, UNC-60A and UNC-60B, are expressed from the unc-60 gene by alternative splicing. RNA interference of UNC-60A caused disorganization of the actin networks in the myoepithelial sheath. UNC-60B, which is known to function in the body-wall muscle, was not necessary or sufficient for actin organization in the myoepithelial sheath. However, mutant forms of UNC-60B with reduced actin-filament-severing activity rescued the UNC-60A-depletion phenotype. UNC-60A has a much weaker filament-severing activity than UNC-60B, suggesting that an ADF/cofilin with weak severing activity is optimal for assembly of actin networks in the myoepithelial sheath. By contrast, strong actin-filament-severing activity of UNC-60B was required for assembly of striated myofibrils in the body-wall muscle. Our results suggest that an optimal level of actin-filament-severing activity of ADF/cofilin is required for assembly of actin networks in the somatic gonad.

Caenorhabditis elegans gelsolin-like protein 1 is a novel actin filament-severing protein with four gelsolin-like repeats

The gelsolin family of proteins is a major class of actin regulatory proteins that sever, cap, and nucleate actin filaments in a calcium-dependent manner and are involved in various cellular processes. Typically, gelsolin-related proteins have three or six repeats of gelsolin-like (G) domain, and each domain plays a distinct role in severing, capping, and nucleation. The Caenorhabditis elegans gelsolin-like protein-1 (gsnl-1) gene encodes an unconventional gelsolin-related protein with four G domains. Sequence alignment suggests that GSNL-1 lacks two G domains that are equivalent to fourth and fifth G domains of gelsolin. In vitro, GSNL-1 severed actin filaments and capped the barbed end in a calcium-dependent manner. However, unlike gelsolin, GSNL-1 remained bound to the side of F-actin with a submicromolar affinity and did not nucleate actin polymerization, although it bound to G-actin with high affinity. These results indicate that GSNL-1 is a novel member of the gelsolin family of actin regulatory proteins and provide new insight into functional diversity and evolution of gelsolin-related proteins.

Complementary RNA amplification methods enhance microarray identification of transcripts expressed in the C. elegans nervous system

Background

DNA microarrays provide a powerful method for global analysis of gene expression. The application of this technology to specific cell types and tissues, however, is typically limited by small amounts of available mRNA, thereby necessitating amplification. Here we compare microarray results obtained with two different methods of RNA amplification to profile gene expression in the C. elegans larval nervous system.

Results

We used the mRNA-tagging strategy to isolate transcripts specifically from C. elegans larval neurons. The WT-Ovation Pico System (WT-Pico) was used to amplify 2 ng of pan-neural RNA to produce labeled cDNA for microarray analysis. These WT-Pico-derived data were compared to microarray results obtained with a labeled aRNA target generated by two rounds of In Vitro Transcription (IVT) of 25 ng of pan-neural RNA. WT-Pico results in a higher fraction of present calls than IVT, a finding consistent with the proposal that DNA-DNA hybridization results in lower mismatch signals than the RNA-DNA heteroduplexes produced by IVT amplification. Microarray data sets from these samples were compared to a reference profile of all larval cells to identify transcripts with elevated expression in neurons. These results were validated by the high proportion of known neuron-expressed genes detected in these profiles and by promoter-GFP constructs for previously uncharacterized genes in these data sets. Together, the IVT and WT-Pico methods identified 2,173 unique neuron-enriched transcripts. Only about half of these transcripts (1,044), however, are detected as enriched by both IVT and WT-Pico amplification.

Conclusion

We show that two different methods of RNA amplification, IVT and WT-Pico, produce valid microarray profiles of gene expression in the C. elegans larval nervous system with a low rate of false positives. However, our results also show that each method of RNA amplification detects a unique subset of bona fide neural-enriched transcripts and thus a wider array of authentic neural genes are identified by the combination of these data sets than by the microarray profiles obtained with either method of RNA amplification alone. With its relative ease of implementation and greater sensitivity, WT-Pico is the preferred method of amplification for cases in which sample RNA is limiting.