Genome-wide searching with base-pairing kernel functions for noncoding RNAs: computational and expression analysis of snoRNA families in Caenorhabditis elegans

Despite the accumulating research on noncoding RNAs (ncRNAs), it is likely that we are seeing only the tip of the iceberg regarding our understanding of the functions and the regulatory roles served by ncRNAs in cellular metabolism, pathogenesis and host-pathogen interactions. Therefore, more powerful computational and experimental tools for analyzing ncRNAs need to be developed. To this end, we propose novel kernel functions, called base-pairing profile local alignment (BPLA) kernels, for analyzing functional ncRNA sequences using support vector machines (SVMs). We extend the local alignment kernels for amino acid sequences in order to handle RNA sequences by using STRAL's; scoring function, which takes into account sequence similarities as well as upstream and downstream base-pairing probabilities, thus enabling us to model secondary structures of RNA sequences. As a test of the performance of BPLA kernels, we applied our kernels to the problem of discriminating members of an RNA family from nonmembers using SVMs. The results indicated that the discrimination ability of our kernels is stronger than that of other existing methods. Furthermore, we demonstrated the applicability of our kernels to the problem of genome-wide search of snoRNA families in the Caenorhabditis elegans genome, and confirmed that the expression is valid in 14 out of 48 of our predicted candidates by using qRT-PCR. Finally, highly expressed six candidates were identified as the original target regions by DNA sequencing.     

Transformation of Caenorhabditis elegans with genes from parasitic nematodes

Our knowledge of many aspects of the molecular biology of animal parasitic nematodes has rapidly expanded in recent years but the classical genetic analysis of this group of organisms has yet to emerge as a viable discipline. For example, it is not possible to routinely perform crosses between single males and females to examine the genetic basis of even simple phenotypes such as anthelmintic resistance. This has meant that the function of many cloned parasite genes can only be inferred from sequence comparison with genes from other organisms where the function is known, or by correlation of DNA polymorphisms linked to the gene with phenotypic differences between strains or individuals. In the absence of classical genetic techniques, a molecular solution is to transform a suitable host with the gene of interest, but what defines a suitable host? Here, Warwick Grant describes recent work that aims to provide such a host.     

Trehalose metabolism genes in Caenorhabditis elegans and filarial nematodes

The sugar trehalose is claimed to be important in the physiology of nematodes where it may function in sugar transport, energy storage and protection against environmental stresses. In this study we investigated the role of trehalose metabolism in nematodes, using Caenorhabditis elegans as a model, and also identified complementary DNA clones putatively encoding genes involved in trehalose pathways in filarial nematodes. In C. elegans two putative trehalose-6-phosphate synthase (tps) genes encode the enzymes that catalyse trehalose synthesis and five putative trehalase (tre) genes encode enzymes catalysing hydrolysis of the sugar. We showed by RT-PCR or Northern analysis that each of these genes is expressed as mRNA at all stages of the C. elegans life cycle. Database searches and sequencing of expressed sequence tag clones revealed that at least one tps gene and two tre genes are expressed in the filarial nematode Brugia malayi, while one tps gene and at least one tre gene were identified for Onchocerca volvulus. We used the feeding method of RNA interference in C. elegans to knock down temporarily the expression of each of the tps and tre genes. Semiquantitative RT-PCR analysis confirmed that expression of each gene was silenced by RNA interference. We did not observe an obvious phenotype for any of the genes silenced individually but gas-chromatographic analysis showed >90% decline in trehalose levels when both tps genes were targeted simultaneously. This decline in trehalose content did not affect viability or development of the nematodes.     

Using Caenorhabditis elegans for functional analysis of genes of parasitic nematodes

Information on the functional genomics of Caenorhabditis elegans has increased significantly in the last few years with the development of RNA interference. In parasitic nematodes, RNA interference has shown some success in gene knockdown but optimisation of this technique will be required before it can be adopted as a reliable functional genomics tool. Comparative studies in C. elegans remain an appropriate alternative for studying the function and regulation of some parasite genes and will be extremely useful for fully exploiting the increasing parasite genome sequence data becoming available.     

Mutation and the experimental evolution of outcrossing in Caenorhabditis elegans

An understanding of the forces that contribute to the phylogenetically widespread phenomenon of sexual reproduction has posed a longstanding problem in evolutionary biology. Mutational theories contend that sex can be maintained when the deleterious mutation rate is sufficiently high, although empirical evidence is equivocal and experimental studies are rare. To test the influence of mutation on the evolution of obligate outcrossing, I introduced a genetic polymorphism for breeding system into populations of the nematode Caenorhabditis elegans with high- and low-mutation rate genetic backgrounds and tracked the change in frequency of females, hermaphrodites, and males over approximately 21 generations. Hermaphrodites invaded all populations, regardless of mutational background. However, experimental populations with elevated mutation rates experienced more outcrossing and greater retention of females. This provides experimental evidence consistent with deleterious mutational explanations for the evolution of sex in principle, but the action of other processes is required to explain the evolution of sex in entirety.     

Tissue-specific expression,heat inducibility,and biological roles of two hsp16 genes in Caenorhabditis elegans

In this report we have examined two new heat shock protein (HSP16) proteins in the nematode Caenorhabditis elegans encoded by the open reading frames F08H9.3 and F08H9.4. The F08H9.3 and F08H9.4 genes are oriented in the same direction next to each other on the chromosome, not sharing any promoter region, unlike other hsp16 genes that share common promoters in pairs. The F08H9.3 and F08H9.4 proteins were expressed in a tissue-specific manner, unlike the other four HSP16 proteins. F08H9.3 was expressed in the pharynx, and F08H9.4 in the excretory canal and a few neuronal cells. While F08H9.3 was weakly induced by heat shock only in the same tissue as under the normal condition, F08H9.4 was newly induced in the intestine. RNA interference experiments showed that these two proteins are required for survival under the heat shock condition.     

Comparisons of the complete sequences of two collagen genes from Caenorhabditis elegans

Several collagen genes have been isolated from the nematode Caenorhabditis elegans. The complete nucleotide sequences of two of these genes, col-1 and col-2, have been determined. These collagen genes differ from vertebrate collagen genes in that they contain only one or two introns, their triple-helical regions are interrupted by nonhelical amino acid sequences and they are smaller. A high degree of nucleotide and amino acid homology exists between col-1 and col-2. In particular, the regions around cysteines and lysines are most highly conserved. The C. elegans genome contains 50 or more collagen genes, the majority of which probably encode cuticle collagens; col-1 and col-2 apparently are members of this large family of cuticle collagen genes.     

Evolution of dnmt-2 and mbd-2-like genes in the free-living nematodes Pristionchus pacificus, Caenorhabditis elegans and Caenorhabditis briggsae

Whole genome sequencing of several metazoan model organisms provides a platform for studying genome evolution. How representative are the genomes of these model organisms for their respective phyla? Within nematodes, for example, the free-living soil nematode Caenorhabditis elegans is a highly derived species with unusual genomic characters, such as a reduced Hox cluster (Curr. Biol., 13, 37–40) and the absence of a Hedgehog signaling system. Here, we describe the recent loss of a DNA methyltransferase-2 gene (dnmt-2) in C.elegans. A dnmt-2-like gene is present in the satellite model organism Pristionchus pacificus, another free-living nematode that diverged from C.elegans 200–300 million years ago. In contrast, C.elegans, Caenorhabditis briggsae and P.pacificus all contain an mbd-2-like gene, which encodes another essential component of the methylation system of higher animals and fungi. Cel-mbd-2 is expressed throughout development and RNA interference (RNAi) experiments result in variable phenotypes. In contrast, Cbr-mbd-2 RNAi results in paralyzed larval or adult worms suggesting recent changes of gene function within the genus Caenorhabditis. We speculate that both genes were part of an ancestral DNA methylation system in nematodes and that gene loss and sequence divergence have abolished DNA methylation in C.elegans.     

Comparative analysis of expression of two p97 homologues in Caenorhabditis elegans

Caenorhabditis elegans possesses two p97/VCP/Cdc48p homologues, named CDC-48.1 (C06A1.1) and CDC-48.2 (C41C4.8), although their expression regulation and functional diversity have not yet been studied. We therefore investigated spatial and temporal expression patterns of two p97 homologues in this study. RT-PCR and Western blot analysis showed that the amount of cdc-48.1 was about twofold of that of cdc-48.2 in adults and that two p97 homologues were induced by ER stress. The amount of cdc-48.1 mRNA did not increase in the cdc-48.2 deletion mutant and vice versa. In situ hybridization showed that two p97 homologues are mainly expressed in germ cells. In vivo expression analysis by using GFP translational fusion constructs revealed that CDC-48.1::GFP was expressed from embryos through to adult worms, while CDC-48.2::GFP was expressed mainly in embryos. These results suggest that the expression of two p97 homologues of C. elegans is differently regulated and independent of each other.     

Cloning and analysis of three new homeobox genes from the nematode Caenorhabditis elegans

Three homeobox-containing genes from the nematode Caenorhabditis elegans are described. Two of them (ceh-11 and ceh-12) were isolated from a genomic library by hybridization at low stringency with the Ascaris lumbricoides homeobox AHB-1. The first clone contains a homeobox defining a new class of homeoboxes (ceh-11). This gene maps on the third chromosome of C. elegans, at the same locus as egl-5, a gene already known to be essential for the determination of specific neurons. In the second clone, sequence analysis revealed the existence of the third helix of a putative homeobox (ceh-12) which is interrupted by an intron located upstream of the codon for the amino acid 45 of the homeodomain. Using the ceh-11 homeobox as a probe, a third homeobox (ceh-13) was isolated from a cDNA library. As ceh-13 belongs to the labial class of homeoboxes, we conclude that, at the time when the nematode lineage diverged from the myriapod-insect and the vertebrate lineages, the duplication which led to the Antp and the labial families of homeoboxes had already taken place.     

Divergent structures of Caenorhabditis elegans cytochrome P450 genes suggest the frequent loss and gain of introns during the evolution of nematodes

The Caenorhabditis elegans genome contains more than 60 cytochrome P450 (CYP) genes. The exon-intron organizations of all of the available and potentially active C. elegans CYP genes were inferred by a newly developed program for predicting protein-coding exons based on the alignment of a genomic DNA sequence and a protein profile. From the predicted amino acid sequences, all of the C. elegans CYP genes except one were classified into three groups, which were closely related to the mammalian drug-metabolizing P450 gene families CYP2, CYP3, and CYP4. The gene structures were strikingly divergent within each group; 20, 10, and 5 unique gene organizations were identified among 40, 18, and 5 genes in the CYP2-, CYP3-, and CYP4-related groups, respectively. The degrees of divergence in gene organization were strongly correlated with those in the amino acid sequences of encoding proteins, and the minimum rate of change in an intron insertion site was estimated to be about 90 times less frequent than amino acid substitutions. Parsimonious analyses suggested that frequent loss and gain of introns has occurred during the evolution of CYP genes in each group after the divergence of nematodes, arthropods, and deuterostomia. Few, if any, incidents of intron sliding were evident, and a model that did not allow intron insertions was highly inconsistent with the observations. All of these findings are explained better by the intron-late view than by the intron-early view.     

Concerted evolution of two novel protein families in Caenorhabditis species

Among a large number of homologous gene clusters in C. elegans, two gene families that appear to undergo concerted evolution were studied in detail. Both gene families are nematode specific and encode small secreted proteins of unknown function. For both families in three Caenorhabditis species, concerted groups of genes are characterized by close genomic proximity and by genes in inverted orientation. The rate of protein evolution in one of the two families could be calibrated by comparison with a closely related nonconcerted singleton gene with one-to-one orthologs in all three species. This comparison suggests that protein evolution in concerted gene clusters is two- to sevenfold accelerated. A broader survey of clustered gene families, focused on adjacent inverted gene pairs, identified an additional seven families in which concerted evolution probably occurs. All nine identified families encode relatively small proteins, eight of them encode putative secreted proteins, and most of these have very unusual amino acid composition or sequence. I speculate that these genes encode rapidly evolving antimicrobial peptides.     

Prediction of structured non-coding RNAs in the genomes of the nematodes Caenorhabditis elegans and Caenorhabditis briggsae

We present a survey for non-coding RNAs and other structured RNA motifs in the genomes of Caenorhabditis elegans and Caenorhabditis briggsae using the RNAz program. This approach explicitly evaluates comparative sequence information to detect stabilizing selection acting on RNA secondary structure. We detect 3,672 structured RNA motifs, of which only 678 are known non-translated RNAs (ncRNAs) or clear homologs of known C. elegans ncRNAs. Most of these signals are located in introns or at a distance from known protein-coding genes. With an estimated false positive rate of about 50% and a sensitivity on the order of 50%, we estimate that the nematode genomes contain between 3,000 and 4,000 RNAs with evolutionary conserved secondary structures. Only a small fraction of these belongs to the known RNA classes, including tRNAs, snoRNAs, snRNAs, or microRNAs. A relatively small class of ncRNA candidates is associated with previously observed RNA-specific upstream elements.