Conservation of gene organization and trans-splicing in the glyceraldehyde-3-phosphate dehydrogenase-encoding genes of Caenorhabditis briggsae.

The genes encoding body-wall-specific glyceraldehyde-3-phosphate dehydrogenase from Caenorhabditis briggsae were sequenced and compared to the homologous genes from Caenorhabditis elegans. The direct tandem organization of these genes, gpd-2 and gpd-3, and the size and location of the two introns in each gene are the same in C. elegans and C. briggsae. Primer-extension studies demonstrated that the two genes in C. briggsae are trans-splice differentially with the same splice leader (SL) RNAs as are observed in C. elegans. The gdp-2 gene is trans-spliced with SL1 while gdp-3 is trans-spliced with SL2. Significant sequence conservation was observed within the promoter regions of each species and may indicate those regions responsible for body-wall-muscle-specific gene expression and/or differential trans-splicing. Comparisons of the sequences suggest that the tandem repeat of the genes has been subjected to concerted evolution and that C. briggsae and C. elegans diverged much earlier than would be anticipated based on morphological similarities alone. Finally, an open reading frame found several hundred nucleotides upstream from gpd-2, in both species, appears to be homologous to the ATP synthase subunit, ATPase inhibitor protein, from bovine mitochondria.     

Pattern of selective constraint in C. elegans and C. briggsae genomes.

Similarity between related genomes may carry information on selective constraint in each of them. We analysed patterns of similarity between several homologous regions of Caenorhabditis elegans and C. briggsae genomes. All homologous exons are quite similar. Alignments of introns and of intergenic sequences contain long gaps, segments where similarity is low and close to that between random sequences aligned using the same parameters, and segments of high similarity. Conservative estimates of the fractions of selectively constrained nucleotides are 72%, 17% and 18% for exons, introns and intergenic sequences, respectively. This implies that the total number of constrained nucleotides within non-coding sequences is comparable to that within coding sequences, so that at least one-third of nucleotides in C. elegans and C. briggsae genomes are under strong stabilizing selection.     

Cloning by synteny: identifying C. briggsae homologues of C. elegans genes.

Phylogenetic comparisons of gene and protein sequences between related species are often used to identify evolutionarily conserved elements that are important for gene expression, function, or regulation. However, homologoues may sometimes be difficult to identify by conventional low stringency hybridisation techniques, if they have undergone substantial sequence divergence. A new approach, cloning by synteny, is described that was used to identify the C. briggsae homologue of the C. elegans sex-determining gene tra-2. We show that four genes tra-2, ppp-1, art-1, and sod-1 are organised in a syntenic cluster and suggest that extensive conservation of gene linkage may exist between C. briggsae and C. elegans. We have also constructed a C. briggsae cDNA library to facilitate characterisation of these genes. Given the rapid progress in the physical mapping and sequencing of the C. elegans genome, cloning by synteny may provide the fastest method for identifying C. briggsae gene homologues, especially for genes encoding novel proteins.     

Functional elements and domains inferred from sequence comparisons of a heat shock gene in two nematodes

Caenorhabditis elegans and Caenorhabditis briggsae are two closely related nematode species that are nearly identical morphologically. Interspecific cross-hybridizing DNA appears to be restricted primarily to coding regions. We compared portions of the hsp-3 homologs, two grp 78-like genes, from C. elegans and C. briggsae and detected regions of DNA identity in the coding region, the 5' flanking DNAs, and the introns. The hsp-3 homologs share approximately 98% and 93% identity at the amino acid and nucleotide levels, respectively. Using the nucleotide substitution rate at the silent third position of the codons, we have estimated a lower limit for the date of divergence between C. elegans and C. briggsae to be approximately 23-32 million years ago. The 5' flanking DNAs and one of the introns contain elements that are highly conserved between C. elegans and C. briggsae. Some of the regions of nucleotide identity in the 5' flanking DNAs correspond to previously detected identities including viral enhancer sequences, a heat shock element, and an element present in the regulatory regions of mammalian grp78 and grp94 genes. We propose that a comparison of C. elegans and C. briggsae sequences will be useful in the detection of potential regulatory and structural elements.     

Conservation of the C.elegans tra-2 3'UTR translational control.

The Caenorhabditis elegans sex-determination gene, tra-2, is translationally regulated by two 28 nt elements (DREs) located in the 3'UTR that bind a factor called DRF. This regulation requires the laf-1 gene activity. We demonstrate that the nematode Caenorhabditis briggsae tra-2 gene and the human oncogene GLI are translationally regulated by elements that are functionally equivalent to DREs. Here, we rename the DREs to TGEs (tra-2 and GLI elements). Similarly to the C.elegans tra-2 TGEs, the C.briggsae tra-2 and GLI TGEs repress translation of a reporter transgene in a laf-1 dependent manner. Furthermore, they regulate poly(A) tail length and bind DRF. We also find that the C.elegans TGEs control translation and poly(A) tail length in C.briggsae and rodent cells. Moreover, these same organisms contain a factor that specifically associates with the C.elegans TGEs. These findings are consistent with the TGE control being present in C.briggsae and rodent cells. Three lines of evidence indicate that C.briggsae tra-2 and GLI are translationally controlled in vivo by TGEs. First, like C.elegans tra-2 TGEs, the C.briggsae tra-2 and GLI TGEs control translation and poly(A) tail lengths in C.briggsae and rodent cells, respectively. Second, the same factor in C.briggsae and mammalian cells that binds to the C.elegans tra-2 TGEs binds the C.briggsae tra-2 and GLI TGEs. Third, deletion of the GLI TGE increases GLI's ability to transform cells. These findings suggest that TGE control is conserved and regulates the expression of other mRNAs.     

Smoke without fire: most reported cases of intron gain in nematodes instead reflect intron losses

Identification of recently gained spliceosomal introns would provide crucial evidence in the continuing debate concerning the age and evolutionary significance of introns. A previously published genomic analysis reported to have identified 122 introns that had been gained since the divergence of the nematodes Caenorhabidits elegans and Caenorhabditis briggsae approximately 100 MYA. However, using newly available genomic sequence from additional Caenorhabditis species, we show that 74% (60/81) of the reported gains in C. elegans are present in a C. briggsae relative. This pattern indicates that these introns represent losses in C. briggsae, not gains in C. elegans. In addition, 61% (25/41) of the reported gains in C. briggsae are present in the more distant C. briggsae relative, in a pattern suggesting that additional reported gains in C. elegans and/or C. briggsae may in fact represent unrecognized losses. These results underscore the dominance of intron loss over intron gain in recent eukaryotic evolution, the pitfalls associated with parsimony in inferring intron gains, and the importance of genomic sequencing of clusters of closely related species for drawing accurate inferences about genome evolution.     

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.     

Evolution of regulatory elements producing a conserved gene expression pattern in Caenorhabditis

Natural selection acts at the level of function, not at the logistical level of how organisms achieve a particular function. Consequently, significant DNA sequence and regulatory differences can achieve the same function, such as a particular gene expression pattern. To investigate how regulatory features underlying a conserved function can evolve, we compared the regulation of a conserved gene expression pattern in the related species Caenorhabditis elegans and C. briggsae. We find that both C. elegans and C. briggsae express the ovo-related zinc finger gene lin-48 in the same pattern in hindgut cells. However, the regulation of this gene by the Pax-2/5/8 protein EGL-38 differs in two important ways. First, specific differences in the regulatory sequences of lin-48 result in the presence of two redundant EGL-38 response elements in C. elegans, whereas the redundancy is absent in C. briggsae. Second, there is a single egl-38 gene in C. briggsae. In contrast, the gene is duplicated in C. elegans, with only one copy retaining the ability to regulate lin-48 in vivo. These results illustrate molecular changes that can occur despite maintenance of conserved gene function in different species.     

A C. elegans gene encodes a protein homologous to mammalian calreticulin.

The gene encoding a C. elegans homologue of the mammalian reticuloplasmin, calreticulin, was cloned and sequenced and the amino-acid sequence of its product deduced. The coding region of the gene comprises three exons separated by introns of 95 and 55 nucleotides, followed by either 158 or 279 bases of 3' non-coding sequence before putative polyadenylation signals. The precursor protein of 395 residues includes an N-terminal signal sequence of 13 residues. The C-terminus has the ER retention signal HDEL preceded by a polyacidic zone similar to known mammalian calreticulins. The sequence shows a 61% identity with mouse calreticulin, increasing to 82% in the proline-rich region of the molecule. Comparison of the C. elegans sequence with the calreticulin-related antigen RAL-1 of Oncocerca volvulus shows 73% identity, excluding the calreticulin C-terminal region. The sequence of this region differs markedly from RAL-1 where the parasite protein has a polybasic stretch and no ER retention signal. The C. elegans gene described here and designated crt-1 was mapped to a region towards the left-hand end of Chromosome V on the physical map of the genome. Southern blotting of genomic DNA indicates that in C. elegans the calreticulin homologue exists in only one form as the product of a single gene.     

The spliceosomal snRNAs of Caenorhabditis elegans.

Nematodes are the only group of organisms in which both cis- and trans-splicing of nuclear mRNAs are known to occur. Most Caenorhabditis elegans introns are exceptionally short, often only 50 bases long. The consensus donor and acceptor splice site sequences found in other animals are used for both cis- and trans-splicing. In order to identify the machinery required for these splicing events, we have characterized the C. elegans snRNAs. They are similar in sequence and structure to those characterized in other organisms, and several sequence variations discovered in the nematode snRNAs provide support for previously proposed structure models. The C. elegans snRNAs are encoded by gene families. We report here the sequences of many of these genes. We find a highly conserved sequence, the proximal sequence element (PSE), about 65 bp upstream of all 21 snRNA genes thus far sequenced, including the SL RNA genes, which specify the snRNAs that provide the 5' exons in trans-splicing. The sequence of the C. elegans PSE is distinct from PSE's from other organisms.     

The glyceraldehyde-3-phosphate dehydrogenase gene family in the nematode, Caenorhabditis elegans: isolation and characterization of one of the genes

The isolation and genomic sequence of one of possibly four glyceraldehyde-3-phosphate dehydrogenase genes in the nematode, Caenorhabditis elegans is presented. The complete nucleotide sequence of the coding as well as the noncoding flanking regions of this gene has been determined. The deduced amino-acid sequence agrees with the sequence of typical glyceraldehyde-3-phosphate dehydrogenase enzymes and its molecular weight of 36,235 agrees with its size determined previously (Yarbrough, P. and Hecht, R. (1984) J. Biol. Chem. 259, 14711-14720). That this isolated gene encodes a nematode glyceraldehyde-3-phosphate dehydrogenase is additionally confirmed by demonstrating its immunoreactivity to an anti-nematode glyceraldehyde-3-phosphate dehydrogenase antibody after its expression as a fusion protein with dihydrofolate reductase. Codon utilization follows a pattern typical of other expressed nematode genes. The gene is split by two introns that are highly conserved in comparison to other introns observed in C. elegans. The placement of one of these introns is conserved with respect to the chicken glyceraldehyde-3-phosphate dehydrogenase gene. Within the 5' flanking sequence homology to actin and the homology 2 block of the major myosin gene (unc-54) is noted. It is of interest that the 3' flanking region contains a CAAAT box, followed by a TATAAT box, before an open reading frame of a closely linked gene that also contains a small AT-rich intron with the nematode consensus splice junction.     

Regulatory elements required for development of caenorhabditis elegans hermaphrodites are conserved in the tra-2 homologue of C. remanei, a male/female sister species

The Caenorhabditis elegans hermaphrodite is essentially a female that produces sperm. In C. elegans, tra-2 promotes female fates and must be repressed to achieve hermaphrodite spermatogenesis. In an effort to learn how mating systems evolve, we have cloned tra-2 from C. remanei, the closest gonochoristic relative of C. elegans. We found its structure to be similar to that of Ce-tra-2 but its sequence to be divergent. RNA interference demonstrates that Cr-tra-2 promotes female fates. Two sites of tra-2 regulation are required for the onset of hermaphrodite spermatogenesis in C. elegans. One, the MX region of TRA-2, is as well conserved in C. remanei as it is in C. briggsae (another male/hermaphrodite species), suggesting that this control is not unique to hermaphrodites. Another, the DRE/TGE element of the tra-2 3' UTR, was not detected by sequence analysis. However, gel-shift assays demonstrate that a factor in C. remanei can bind specifically to the Cr-tra-2 3' UTR, suggesting that this translational control is also conserved. We propose that both controls are general and do not constitute a novel "switch" that enables sexual mosaicism in hermaphrodites. However, subtle quantitative or qualitative differences in their employment may underlie differences in mating system seen in Caenorhabditis.