Exploring the Caenorhabditis elegans and Drosophila melanogaster genomes to understand neuropeptide and peptidase function

Comparison of peptidase gene families in the newly released Drosophila melanogaster and Caenorhabditis elegans genomes highlights important differences in peptidase distributions with relevance to the evolution of both form and function in these two organisms and can help to identify the most appropriate model when using comparative studies relevant to the human condition.     

MetaBlasts: tracing protein tyrosine phosphatase gene family roots from Man to Drosophila melanogaster and Caenorhabditis elegans genomes

At increasing speed, sequencing data are being made public from both complex and simple life forms. Although biomedical interests tend to focus on mammalian genes, only simple organisms allow rapid genetic manipulation and functional analysis. A prerequisite for the meaningful extrapolation of gene functional studies from invertebrates to man is that the orthologs under study are unambiguously linked. However, identifying orthologs is not trivial, especially where large gene families are involved. We present here an automated sequence analysis procedure that allows a rapid visualization of most likely ortholog pairs. We illustrate the utility of this approach for the human gene family of protein tyrosine phosphatases (PTPs) as compared with the full set of Caenorhabditis elegans and Drosophila melanogaster conceptual ORFs. The approach is based on a reciprocal series of BLAST searches, which are automatically stored and represented in an HTML-formatted table. We have used this 'MetaBlast' approach to compile lists of human PTPs and their worm and fly orthologs. Many of these PTP orthologs had not been previously identified as such.     

Germ plasm in Caenorhabditis elegans, Drosophila and Xenopus

Special cytoplasm, called germ plasm, that is essential for the differentiation of germ cells is localized in a particular region of Caenorhabditis elegans, Drosophila and Xenopus eggs. The mode of founder cell formation of germline, the origin and behavior of the germline granules, and the molecules localized in germline cells are compared in these organisms. The common characteristics of the organisms are mainly as follows. First, the founder cells of germline are established before the intiation of gastrulation. Second, the germline granules or their derivatives are always present in germline cells or germ cells throughout the life cycle in embryos, larvae, and adults. Lastly, among the proteins localized in the germ plasm, only Vasa protein or its homolog is detected in the germline cells or germ cells throughout the life cycle. As the protein of vasa homolog has been reported to be also localized in the germline-specific structure or nuage in some of the organisms without the germ plasm, the possibility that the mechanism for differentiation of primordial germ cells is basically common in all organisms with or without the germ plasm is discussed.     

Distilling artificial recombinants from large sets of complete mtDNA genomes

Background

Large-scale genome sequencing poses enormous problems to the logistics of laboratory work and data handling. When numerous fragments of different genomes are PCR amplified and sequenced in a laboratory, there is a high immanent risk of sample confusion. For genetic markers, such as mitochondrial DNA (mtDNA), which are free of natural recombination, single instances of sample mix-up involving different branches of the mtDNA phylogeny would give rise to reticulate patterns and should therefore be detectable.

Methodology/Principal Findings

We have developed a strategy for comparing new complete mtDNA genomes, one by one, to a current skeleton of the worldwide mtDNA phylogeny. The mutations distinguishing the reference sequence from a putative recombinant sequence can then be allocated to two or more different branches of this phylogenetic skeleton. Thus, one would search for two (or three) near-matches in the total mtDNA database that together best explain the variation seen in the recombinants. The evolutionary pathway from the mtDNA tree connecting this pair together with the recombinant then generate a grid-like median network, from which one can read off the exchanged segments.

Conclusions

We have applied this procedure to a large collection of complete human mtDNA sequences, where several recombinants could be distilled by our method. All these recombinant sequences were subsequently corrected by de novo experiments – fully concordant with the predictions from our data-analytical approach.     

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.     

Limited microsynteny between the genomes of Pristionchus pacificus and Caenorhabditis elegans

Nematodes are an attractive group of organisms for studying the evolution of developmental processes. Pristionchus pacificus was established as a satellite organism for comparing vulva development and other processes to Caenorhabditis elegans. The generation of a genetic linkage map of P.pacificus has provided a first insight into the structure and organization of the genome of this species. Pristionchus pacificus and C.elegans are separated from one another by >100 000 000 years such that the structure of the genomes of these two nematodes might differ substantially. To evaluate the amount of synteny between the two genomes, we have obtained 126 kb of continuous genomic sequence of P.pacificus, flanking the developmental patterning gene pal-1. Of the 20 predicted open reading frames in this interval, 11 have C.elegans orthologs. Ten of these 11 orthologs are located on C.elegans chromosome III, indicating the existence of synteny. However, most of these genes are distributed over a 12 Mb interval of the C.elegans genome and only three pairs of genes show microsynteny. Thus, intrachromosomal rearrange ments occur frequently in nematodes, limiting the likelihood of identifying orthologous genes of P.pacificus and C.elegans based on positional information within the two genomes.     

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.     

Counting mutagenized genomes and optimizing genetic screens in Caenorhabditis elegans

In genetic screens, the number of mutagenized gametes examined is an important parameter for evaluating screen progress, the number of genes of a given mutable phenotype, gene size, cost, and labor. Since genetic screens often entail examination of thousands or tens of thousands of animals, strategies for optimizing genetics screens are important for minimizing effort while maximizing the number of mutagenized gametes examined. To date, such strategies have not been described for genetic screens in the nematode Caenorhabditis elegans. Here we review general principles of genetic screens in C. elegans, and use a modified binomial strategy to obtain a general expression for the number of mutagenized gametes examined in a genetic screen. We use this expression to calculate optimal screening parameters for a large range of genetic screen types. In addition, we developed a simple online genetic-screen-optimization tool that can be used independently of this paper. Our results demonstrate that choosing the optimal F2-to-F1 screening ratio can significantly improve screen efficiency.     

Hill-Robertson interference is a minor determinant of variations in codon bias across Drosophila melanogaster and Caenorhabditis elegans genomes

According to population genetics models, genomic regions with lower crossing-over rates are expected to experience less effective selection because of Hill-Robertson interference (HRi). The effect of genetic linkage is thought to be particularly important for a selection of weak intensity such as selection affecting codon usage. Consistent with this model, codon bias correlates positively with recombination rate in Drosophila melanogaster and Caenorhabditis elegans. However, in these species, the G+C content of both noncoding DNA and synonymous sites correlates positively with recombination, which suggests that mutation patterns and recombination are associated. To remove this effect of mutation patterns on codon bias, we used the synonymous sites of lowly expressed genes that are expected to be effectively neutral sites. We measured the differences between codon biases of highly expressed genes and their lowly expressed neighbors. In D. melanogaster we find that HRi weakly reduces selection on codon usage of genes located in regions of very low recombination; but these genes only comprise 4% of the total. In C. elegans we do not find any evidence for the effect of recombination on selection for codon bias. Computer simulations indicate that HRi poorly enhances codon bias if the local recombination rate is greater than the mutation rate. This prediction of the model is consistent with our data and with the current estimate of the mutation rate in D. melanogaster. The case of C. elegans, which is highly self-fertilizing, is discussed. Our results suggest that HRi is a minor determinant of variations in codon bias across the genome.     

Regulation of Caenorhabditis elegans male mate searching behavior by the nuclear receptor DAF-12

Coordination of animal behavior with reproductive status is often achieved through elaboration of hormones by the gonad. In the nematode Caenorhabditis elegans, adult males explore their environment to locate mates. Mate searching is regulated by presence of mates, nutritional status, and a signal from the gonad. Here we show that the gonadal signal acts via the nuclear receptor DAF-12, a protein known to regulate several C. elegans life-history traits. DAF-12 has both activational and organizational functions to stimulate exploratory behavior and acts downstream of the gonadal signal, outside of the gonad. DAF-12 acts upstream of sensory input from mating partners and physiological signals indicating nutritional status. Mate searching was rescued in germ-line ablated animals, but not if both germ line and somatic gonad were ablated, by a precursor of the DAF-12 ligand, dafachronic acid (DA). The results are interpreted to suggest that the germ line produces a DA precursor that is converted to DA outside of the germ line, possibly in the somatic gonad. As it does in other pathways in which it functions, in regulation of male mate searching behavior DAF-12 acts at a choice point between alternatives favoring reproduction (mate searching) vs. survival (remaining on food).     

C-Type lectin-like domains in Caenorhabditis elegans: predictions from the complete genome sequence

Protein modules related to the C-type carbohydrate-recognition domains of animal lectins are found in at least 125 proteins encoded in the Caenorhabditis elegans genome. Within these proteins, 183 C-type lectin-like domains (CTLDs) have been identified. The proteins have been classified based on the overall arrangement of modules within the polypeptides and based on sequence similarity between the CTLDs. The C.elegans proteins generally have different domain organization from known mammalian proteins containing CTLDs. Most of the CTLDs are divergent in sequence from those in mammalian proteins. However, 19 show conservation of most of the amino acid residues that ligate Ca(2+)to form a carbohydrate-binding site in vertebrate C-type carbohydrate-recognition domains. Seven of these domains are particularly similar in sequence to mannose- and N-acetylglucosamine-binding domains in the vicinity of this Ca(2+)site.     

Neutral evolution of ten types of mariner transposons in the genomes of Caenorhabditis elegans and Caenorhabditis briggsae

Ten types of mariner transposable elements (232 individual sequences) are present in the completed genomic DNA sequence of Caenorhabditis elegans and the partial sequence of Caenorhabditis briggsae. We analyze these replicated instances of mariner evolution and find that elements of a type have evolved within their genomes under no selection on their transposase genes. Seven of the ten reconstructed ancestral mariners carry defective transposase genes. Selection has acted during the divergence of some ancestral elements. The neutrally-evolving mariners are used to analyze the pattern of molecular evolution in Caenorhabditis. There is a significant mutational bias against transversions and significant variation in rates of change across sites. Deletions accumulate at a rate of 0.034 events/bp per substitution/site, with an average size of 166 bp (173 gaps observed). Deletions appear to obliterate preexisting deletions over time, creating larger gaps. Insertions accumulate at a rate of 0.019 events/bp per substitution/site, with an average size of 151 bp (61 events). Although the rate of deletion is lower than most estimates in other species, the large size of deletions causes rapid elimination of neutral DNA: a mariners half-life (the time by which half an elements sequence should have been deleted) is ~0.1 subsitutions/site. This high rate of DNA deletion may explain the compact nature of the nematode genome. When this work was done, both authors were affiliated with the University of Illinois at Urbana-Champaign. Dr. Witherspoon is now working in the private sector, Dr. Robertson remains affiliated with the University of Illinois.     

New BEL-like LTR-retrotransposons in Fugu rubripes, Caenorhabditis elegans, and Drosophila melanogaster

The BEL group of retroelements is present in greater numbers, variety and taxonomic range than may have been thought previously. In addition to the insects, nematodes and schistosomes, BEL-like elements are present in echinoderms, urochordates, and at least two highly diverged species of fish. We describe one new full-length BEL-like element in Fugu that we call Suzu, another in Drosophila that we call Tinker, and seven new families in C. elegans. Many of the C. elegans elements have an unusual insertion at the 5' end. The previously known Roo, TRAM and Telemac are also BEL-like retrotransposons. Some BEL-like elements have captured an envelope gene, probably from other retroelements in some cases but from a phlebovirus in one case.     

Cadherin superfamily proteins in Caenorhabditis elegans and Drosophila melanogaster

The ability to form selective cell-cell adhesions is an essential property of metazoan cells. Members of the cadherin superfamily are important regulators of this process in both vertebrates and invertebrates. With the advent of genome sequencing projects, determination of the full repertoire of cadherins available to an organism is possible and here we present the identification and analysis of the cadherin repertoires in the genomes of Caenorhabditis elegans and Drosophila melanogaster. Hidden Markov models of cadherin domains were matched to the protein sequences obtained from the translation of the predicted gene sequences. Matches were made to 21 C. elegans and 18 D. melanogaster sequences. Experimental and theoretical work on C. elegans sequences, and data from ESTs, show that three pairs of genes, and two triplets, should be merged to form five single genes. It also produced sequence changes at one or both of the 5' and 3' termini of half the sequences. In D. melanogaster it is probable that two of the cadherin genes should also be merged together and that three cadherin genes should be merged with other neighbouring genes.Of the 15 cadherin proteins found in C. elegans, 13 have the features of cell surface proteins, signal sequences and transmembrane helices; the other two have only signal sequences. Of the 17 in D. melanogaster, 11 at present have both features and another five have transmembrane helices. The evidence currently available suggests about one-third of the cadherins in the two organisms can be grouped into subfamilies in which all, or parts of, the molecules are conserved. Each organism also has a approximately 980 residue protein (CDH-11 and CG11059) with two cadherin domains and whose sequences match well over their entire length two proteins from human brain. Two proteins in C. elegans, HMR-1A and HMR-1B, and three in D. melanogaster, CadN, Shg and CG7527, have cytoplasmic domains homologous to those of the classical cadherin genes of chordates but their extracellular regions have different domain structures. Other common subclasses include the seven-helix membrane cadherins, Fat-like protocadherins and the Ret-like cadherins. At present, the remaining cadherins have no obvious similarities in their extracellular domain architecture or homologies to their cytoplasmic domains and may, therefore, represent species-specific or phylum-specific molecules.     

The Caenorhabditis elegans NR4A nuclear receptor is required for spermatheca morphogenesis

The gene nhr-6 encodes the Caenorhabditis elegans ortholog of the NR4A nuclear receptor. We determined the biological functions of NHR-6 through the isolation and characterization of a deletion allele of nhr-6, lg6001. We demonstrate that nhr-6 has an essential role in the development of the C. elegans somatic gonad. Specifically, nhr-6 is required for the development of the hermaphrodite spermatheca, a somatic gonad organ that serves as the site of sperm storage and oocyte fertilization. Using a variety of spermatheca cell markers, we have determined that loss of nhr-6 function causes severe morphological defects in the spermatheca and associated spermathecal valves. This appears to be due to specific requirements for nhr-6 in regulating cell proliferation and cell differentiation during development of these structures. The improper development of these structures in nhr-6(lg6001) mutants leads to defects in ovulation and significantly reduced fecundity of C. elegans hermaphrodites. The phenotypes of nhr-6(lg6001) mutants are consistent with a role for nhr-6 in organogenesis, similar to the functions of its mammalian homologs.     

Mitochondrial pseudogenes in the nuclear genomes of Drosophila

Mitochondrial pseudogenes in nuclear chromosomes (numts) have been detected in the genomes of a diverse range of eukaryotic species. However, the numt content of different genomes and their properties is not uniform, and study of these differences provides insight into the mechanisms and dynamics of genome evolution in different organisms. In the genus Drosophila, numts have previously only been identified on a genome-wide scale in the melanogaster subgroup. The present study extends the identification to 11 species of the Drosophila genus. We identify a total of 302 numts and show that the numt complement is highly variable in Drosophilids, ranging from just 4 in D. melanogaster to 67 in D. willistoni, broadly correlating with genome size. Many numts have undergone large-scale rearrangements in the nucleus, including interruptions, inversions, deletions and duplications of sequence of variable size. Estimating the age of the numts in the nucleus by phylogenetic tree reconstruction reveals the vast majority of numts to be recent gains, 90% having arisen on terminal branches of the species tree. By identifying paralogs and counting duplications among the extant numts we estimate that 23% of extant numts arose through post-insertion duplications. We estimate genus average rates of insertion of 0.75 per million years, and a duplication rate of 0.010 duplications per numt per million years.