Identification of a nematode chemosensory gene family

Anatomy and transcriptomics of the common jingle shell (Bivalvia, Anomiidae) support a sensory function for bivalve tentacles

Animals have evolved numerous mechanisms to perceive and interact with the environment that can be translated into different sensory modalities. However, the genomic and phenotypic features that support sensory functions remain enigmatic for many invertebrates, such as bivalves, an ecologically and economically important taxonomic group. No repertoire of sensory genes has been characterized in bivalves, representing a significant knowledge gap in molluscan sensory biology. Here, we gather multiple lines of evidence to explore the specialized sensory function of bivalve tentacles in the common jingle shell, Anomia simplex. In addition to applying microscopy techniques, we performed transcriptome sequencing of dissected tentacles using phylogenetically-informed annotation to identify candidate receptors. Our results demonstrate the expression of candidate GPCRs, including one opsin type, five small-molecule receptors, and 11 chemosensory-related receptors, supporting the involvement of sensory neurons in the organ, likely in association with the ciliated receptor cells observed along the tentacle surface. In addition, we identified seven ionotropic receptors as putative chemosensory receptors and one member of the Piezo mechanosensitive ion channel, which might be involved in touch sensation by ciliated sensory receptors. Our results provide the first evidence of putative sensory genes expressed in a bivalve sensory organ, representing an important starting point to investigate chemosensation in this class.

G protein-coupled receptor-based thermosensation determines temperature acclimatization of Caenorhabditis elegans

Animals must sense and acclimatize to environmental temperatures for survival, yet their thermosensing mechanisms other than transient receptor potential (TRP) channels remain poorly understood. We identify a trimeric G protein-coupled receptor (GPCR), SRH-40, which confers thermosensitivity in sensory neurons regulating temperature acclimatization in Caenorhabditis elegans. Systematic knockdown of 1000 GPCRs by RNAi reveals GPCRs involved in temperature acclimatization, among which srh-40 is highly expressed in the ADL sensory neuron, a temperature-responsive chemosensory neuron, where TRP channels act as accessorial thermoreceptors. In vivo Ca2+ imaging demonstrates that an srh-40 mutation reduced the temperature sensitivity of ADL, resulting in supranormal temperature acclimatization. Ectopically expressing SRH-40 in a non-warmth-sensing gustatory neuron confers temperature responses. Moreover, temperature-dependent SRH-40 activation is reconstituted in Drosophila S2R+ cells. Overall, SRH-40 may be involved in thermosensory signaling underlying temperature acclimatization. We propose a dual thermosensing machinery through a GPCR and TRP channels in a single sensory neuron.

Transcriptome reveals the toxicity difference of dimethyl disulfide by contact and fumigation on Meloidogyne incognita through calcium channel-mediated oxidative phosphorylation

The prevention and control of root-knot nematode disease has been posing a severe challenge worldwide. Fumigant dimethyl disulfide (DMDS) has excellent biological activity against nematodes. However, DMDS displays significant differences in contact and fumigation toxicity on nematodes. The specific regulatory mechanisms of DMDS on nematodes were investigated by characterizing the ultrastructure of nematodes, examining the physiological and biochemical indicators, and conducting transcriptome high-throughput sequencing. As indicated by the results, DMDS fumigation exhibited the biological activity of against M. incognita 121 times higher than DMDS contact. DMDS contact destroyed nematode body wall cells. Besides, DMDS fumigation destroyed the structure of pseudocoelom. DMDS treatment expedited the oxygen consumption of nematode while inhibiting acetylcholinesterase activity. As indicated by the analysis of vital signaling pathways based on transcriptome, DMDS based on the contact mode penetrated directly into the nematode through the body wall and subsequently affected calcium channels in the body wall and muscle, disrupting their structure; it serves as an uncoupling agent to interfere with ATP synthase. Moreover, DMDS based on the fumigation mode entered the body through the respiratory pathway of olfactory perception-oxygen exchange and subsequently affected calcium channels in the nerve; eventually, DMDS acted on complex IV or complex I.

Arrestin-mediated desensitization enables intraneuronal olfactory discrimination in Caenorhabditis elegans

In the mammalian olfactory system, cross-talk between olfactory signals is minimized through physical isolation: individual neurons express one or few olfactory receptors among those encoded in the genome. Physical isolation allows for segregation of stimuli during signal transduction; however, in the nematode worm Caenorhabditis elegans, ∼1,300 olfactory receptors are primarily expressed in only 32 neurons, precluding this strategy. Here, we report genetic and behavioral evidence that β-arrestin-mediated desensitization of olfactory receptors, working downstream of the kinase GRK-1, enables discrimination between intraneuronal olfactory stimuli. Our findings suggest that C. elegans exploits β-arrestin desensitization to maximize responsiveness to novel odors, allowing for behaviorally appropriate responses to olfactory stimuli despite the large number of olfactory receptors signaling in single cells. This represents a fundamentally different solution to the problem of olfactory discrimination than that which evolved in mammals, allowing for economical use of a limited number of sensory neurons.

Chemosensory signal transduction in Caenorhabditis elegans

Chemosensory neurons translate perception of external chemical cues, including odorants, tastants, and pheromones, into information that drives attraction or avoidance motor programs. In the laboratory, robust behavioral assays, coupled with powerful genetic, molecular and optical tools, have made Caenorhabditis elegans an ideal experimental system in which to dissect the contributions of individual genes and neurons to ethologically relevant chemosensory behaviors. Here, we review current knowledge of the neurons, signal transduction molecules and regulatory mechanisms that underlie the response of C. elegans to chemicals, including pheromones. The majority of identified molecules and pathways share remarkable homology with sensory mechanisms in other organisms. With the development of new tools and technologies, we anticipate that continued study of chemosensory signal transduction and processing in C. elegans will yield additional new insights into the mechanisms by which this animal is able to detect and discriminate among thousands of chemical cues with a limited sensory neuron repertoire.

Signatures of conserved and unique molecular features in Afrotheria

Afrotheria is a clade of African-origin species with striking dissimilarities in appearance and habitat. In this study, we compared whole proteome sequences of six Afrotherian species to obtain a broad viewpoint of their underlying molecular make-up, to recognize potentially unique proteomic signatures. We find that 62% of the proteomes studied here, predominantly involved in metabolism, are orthologous, while the number of homologous proteins between individual species is as high as 99.5%. Further, we find that among Afrotheria, L. africana has several orphan proteins with 112 proteins showing < 30% sequence identity with their homologues. Rigorous sequence searches and complementary approaches were employed to annotate 156 uncharacterized protein sequences and 28 species-specific proteins. For 122 proteins we predicted potential functional roles, 43 of which we associated with protein- and nucleic-acid binding roles. Further, we analysed domain content and variations in their combinations within Afrotheria and identified 141 unique functional domain architectures, highlighting proteins with potential for specialized functions. Finally, we discuss the potential relevance of highly represented protein families such as MAGE-B2, olfactory receptor and ribosomal proteins in L. africana and E. edwardii, respectively. Taken together, our study reports the first comparative study of the Afrotherian proteomes and highlights salient molecular features.

Transcriptomic analysis of hookworm Ancylostoma ceylanicum life cycle stages reveals changes in G-protein coupled receptor diversity associated with the onset of parasitism

Free-living nematodes respond to variable and unpredictable environmental stimuli whereas parasitic nematodes exist in a more stable host environment. A positive correlation between the presence of environmental stages in the nematode life cycle and an increasing number of G-protein coupled receptors (GPCRs) reflects this difference in free-living and parasitic lifestyles. As hookworm larvae move from the external environment into a host, they detect uncharacterized host components, initiating a signalling cascade that results in the resumption of development and eventual maturation. Previous studies suggest this process is mediated by GPCRs in amphidial neurons. Here we set out to uncover candidate GPCRs required by a hookworm to recognise its host. First, we identified all potential Ancylostoma ceylanicum GPCRs encoded in the genome. We then used life cycle stage-specific RNA-seq data to identify differentially expressed GPCRs between the free-living infective L3 (iL3) and subsequent parasitic stages to identify receptors involved in the transition to parasitism. We reasoned that GPCRs involved in host recognition and developmental activation would be expressed at higher levels in the environmental iL3 stage than in subsequent stages. Our results support the model that a decrease in GPCR diversity occurs as the larvae develop from the free-living iL3 stage to the parasitic L3 (pL3) in the host over 24-72 h. We find that overall GPCR expression and diversity is highest in the iL3 compared with subsequent parasitic stages. By 72 h, there was an approximately 50% decrease in GPCR richness associated with the moult from the pL3 to the L4. Taken together, our data uncover a negative correlation between GPCR diversity and parasitic development in hookworm. Finally, we demonstrate proof of principal that Caenorhabditis elegans can be used as a heterologous system to examine the expression pattern of candidate host signal chemoreceptors (CRs) from hookworm. We observe expression of candidate host signal CRs in C. elegans, demonstrating that C. elegans can be effectively used as a surrogate to identify expressed hookworm genes. We present several preliminary examples of this strategy and confirm a candidate CR as neuronally expressed.

Repertoires of G protein-coupled receptors for Ciona-specific neuropeptides

Neuropeptides play pivotal roles in various biological events in the nervous, neuroendocrine, and endocrine systems, and are correlated with both physiological functions and unique behavioral traits of animals. Elucidation of functional interaction between neuropeptides and receptors is a crucial step for the verification of their biological roles and evolutionary processes. However, most receptors for novel peptides remain to be identified. Here, we show the identification of multiple G protein-coupled receptors (GPCRs) for species-specific neuropeptides of the vertebrate sister group, Ciona intestinalis Type A, by combining machine learning and experimental validation. We developed an original peptide descriptor-incorporated support vector machine and used it to predict 22 neuropeptide-GPCR pairs. Of note, signaling assays of the predicted pairs identified 1 homologous and 11 Ciona-specific neuropeptide-GPCR pairs for a 41% hit rate: the respective GPCRs for Ci-GALP, Ci-NTLP-2, Ci-LF-1, Ci-LF-2, Ci-LF-5, Ci-LF-6, Ci-LF-7, Ci-LF-8, Ci-YFV-1, and Ci-YFV-3. Interestingly, molecular phylogenetic tree analysis revealed that these receptors, excluding the Ci-GALP receptor, were evolutionarily unrelated to any other known peptide GPCRs, confirming that these GPCRs constitute unprecedented neuropeptide receptor clusters. Altogether, these results verified the neuropeptide-GPCR pairs in the protochordate and evolutionary lineages of neuropeptide GPCRs, and pave the way for investigating the endogenous roles of novel neuropeptides in the closest relatives of vertebrates and the evolutionary processes of neuropeptidergic systems throughout chordates. In addition, the present study also indicates the versatility of the machine-learning-assisted strategy for the identification of novel peptide-receptor pairs in various organisms.

An atlas of Caenorhabditis elegans chemoreceptor expression

One goal of modern day neuroscience is the establishment of molecular maps that assign unique features to individual neuron types. Such maps provide important starting points for neuron classification, for functional analysis, and for developmental studies aimed at defining the molecular mechanisms of neuron identity acquisition and neuron identity diversification. In this resource paper, we describe a nervous system-wide map of the potential expression sites of 244 members of the largest gene family in the C. elegans genome, rhodopsin-like (class A) G-protein-coupled receptor (GPCR) chemoreceptors, using classic gfp reporter gene technology. We cover representatives of all sequence families of chemoreceptor GPCRs, some of which were previously entirely uncharacterized. Most reporters are expressed in a very restricted number of cells, often just in single cells. We assign GPCR reporter expression to all but two of the 37 sensory neuron classes of the sex-shared, core nervous system. Some sensory neurons express a very small number of receptors, while others, particularly nociceptive neurons, coexpress several dozen GPCR reporter genes. GPCR reporters are also expressed in a wide range of inter- and motorneurons, as well as non-neuronal cells, suggesting that GPCRs may constitute receptors not just for environmental signals, but also for internal cues. We observe only one notable, frequent association of coexpression patterns, namely in one nociceptive amphid (ASH) and two nociceptive phasmid sensory neurons (PHA, PHB). We identified GPCRs with sexually dimorphic expression and several GPCR reporters that are expressed in a left/right asymmetric manner. We identified a substantial degree of GPCR expression plasticity; particularly in the context of the environmentally-induced dauer diapause stage when one third of all tested GPCRs alter the cellular specificity of their expression within and outside the nervous system. Intriguingly, in a number of cases, the dauer-specific alterations of GPCR reporter expression in specific neuron classes are maintained during postdauer life and in some case new patterns are induced post-dauer, demonstrating that GPCR gene expression may serve as traits of life history. Taken together, our resource provides an entry point for functional studies and also offers a host of molecular markers for studying molecular patterning and plasticity of the nervous system.

C. elegans avoids toxin-producing Streptomyces using a seven transmembrane domain chemosensory receptor

Predators and prey co-evolve, each maximizing their own fitness, but the effects of predator–prey interactions on cellular and molecular machinery are poorly understood. Here, we study this process using the predator Caenorhabditis elegans and the bacterial prey Streptomyces, which have evolved a powerful defense: the production of nematicides. We demonstrate that upon exposure to Streptomyces at their head or tail, nematodes display an escape response that is mediated by bacterially produced cues. Avoidance requires a predicted G-protein-coupled receptor, SRB-6, which is expressed in five types of amphid and phasmid chemosensory neurons. We establish that species of Streptomyces secrete dodecanoic acid, which is sensed by SRB-6. This behavioral adaptation represents an important strategy for the nematode, which utilizes specialized sensory organs and a chemoreceptor that is tuned to recognize the bacteria. These findings provide a window into the molecules and organs used in the coevolutionary arms race between predator and potential prey.

Evolutionary History of Chemosensory-Related Gene Families across the Arthropoda

Chemosensory-related gene (CRG) families have been studied extensively in insects, but their evolutionary history across the Arthropoda had remained relatively unexplored. Here, we address current hypotheses and prior conclusions on CRG family evolution using a more comprehensive data set. In particular, odorant receptors were hypothesized to have proliferated during terrestrial colonization by insects (hexapods), but their association with other pancrustacean clades and with independent terrestrial colonizations in other arthropod subphyla have been unclear. We also examine hypotheses on which arthropod CRG family is most ancient. Thus, we reconstructed phylogenies of CRGs, including those from new arthropod genomes and transcriptomes, and mapped CRG gains and losses across arthropod lineages. Our analysis was strengthened by including crustaceans, especially copepods, which reside outside the hexapod/branchiopod clade within the subphylum Pancrustacea. We generated the first high-resolution genome sequence of the copepod Eurytemora affinis and annotated its CRGs. We found odorant receptors and odorant binding proteins present only in hexapods (insects) and absent from all other arthropod lineages, indicating that they are not universal adaptations to land. Gustatory receptors likely represent the oldest chemosensory receptors among CRGs, dating back to the Placozoa. We also clarified and confirmed the evolutionary history of antennal ionotropic receptors across the Arthropoda. All antennal ionotropic receptors in E. affinis were expressed more highly in males than in females, suggestive of an association with male mate-recognition behavior. This study is the most comprehensive comparative analysis to date of CRG family evolution across the largest and most speciose metazoan phylum Arthropoda.

Natural Marine and Synthetic Xenobiotics Get on Nematode's Nerves: Neuro-Stimulating and Neurotoxic Findings in Caenorhabditis elegans

Marine algae release a plethora of organic halogenated compounds, many of them with unknown ecological impact if environmentally realistic concentrations are applied. One major compound is dibromoacetic acid (DBAA) which was tested for neurotoxicity in the invertebrate model organism Caenorhabditis elegans (C. elegans). This natural compound was compared with the widespread synthetic xenobiotic tetrabromobisphenol-A (TBBP-A) found in marine sediments and mussels. We found a neuro-stimulating effect for DBAA; this is contradictory to existing toxicological reports of mammals that applied comparatively high dosages. For TBBP-A, we found a hormetic concentration-effect relationship. As chemicals rarely occur isolated in the environment, a combination of both organobromines was also examined. Surprisingly, the presence of DBAA increased the toxicity of TBBP-A. Our results demonstrated that organohalogens have the potential to affect single organisms especially by altering the neurological processes, even with promoting effects on exposed organisms.

Genome-wide variations in a natural isolate of the nematode Caenorhabditis elegans

Background

Increasing genetic and phenotypic differences found among natural isolates of C. elegans have encouraged researchers to explore the natural variation of this nematode species.

Results

Here we report on the identification of genomic differences between the reference strain N2 and the Hawaiian strain CB4856, one of the most genetically distant strains from N2. To identify both small- and large-scale genomic variations (GVs), we have sequenced the CB4856 genome using both Roche 454 (~400 bps single reads) and Illumina GA DNA sequencing methods (101 bps paired-end reads). Compared to previously described variants (available in WormBase), our effort uncovered twice as many single nucleotide variants (SNVs) and increased the number of small InDels almost 20-fold. Moreover, we identified and validated large insertions, most of which range from 150 bps to 1.2 kb in length in the CB4856 strain. Identified GVs had a widespread impact on protein-coding sequences, including 585 single-copy genes that have associated severe phenotypes of reduced viability in RNAi and genetics studies. Sixty of these genes are homologs of human genes associated with diseases. Furthermore, our work confirms previously identified GVs associated with differences in behavioural and biological traits between the N2 and CB4856 strains.

Conclusions

The identified GVs provide a rich resource for future studies that aim to explain the genetic basis for other trait differences between the N2 and CB4856 strains.

A survey of putative secreted and transmembrane proteins encoded in the C. elegans genome

BACKGROUND

Almost half of the Caenorhabditis elegans genome encodes proteins with either a signal peptide or a transmembrane domain. Therefore a substantial fraction of the proteins are localized to membranes, reside in the secretory pathway or are secreted. While these proteins are of interest to a variety of different researchers ranging from developmental biologists to immunologists, most of secreted proteins have not been functionally characterized so far.

RESULTS

We grouped proteins containing a signal peptide or a transmembrane domain using various criteria including evolutionary origin, common domain organization and functional categories. We found that putative secreted proteins are enriched for small proteins and nematode-specific proteins. Many secreted proteins are predominantly expressed in specific life stages or in one of the two sexes suggesting stage- or sex-specific functions. More than a third of the putative secreted proteins are upregulated upon exposure to pathogens, indicating that a substantial fraction may have a role in immune response. Slightly more than half of the transmembrane proteins can be grouped into broad functional categories based on sequence similarity to proteins with known function. By far the largest groups are channels and transporters, various classes of enzymes and putative receptors with signaling function.

CONCLUSION

Our analysis provides an overview of all putative secreted and transmembrane proteins in C. elegans. This can serve as a basis for selecting groups of proteins for large-scale functional analysis using reverse genetic approaches.

Cross-Genome Clustering of Human and C. elegans G-Protein Coupled Receptors

G-protein coupled receptors (GPCRs) are one of the largest groups of membrane proteins and are popular drug targets. The work reported here attempts to perform cross-genome phylogeny on GPCRs from two widely different taxa, human versus C. elegans genomes and to address the issues on evolutionary plasticity, to identify functionally related genes, orthologous relationship, and ligand binding properties through effective bioinformatic approaches. Through RPS blast around 1106 nematode GPCRs were given chance to associate with previously established 8 types of human GPCR profiles at varying E-value thresholds and resulted 32 clusters were illustrating co-clustering and class-specific retainsionship. In the significant thresholds, 81% of the C. elegans GPCRs were associated with 32 clusters and 27 C. elegans GPCRs (2%) inferred for orthology. 177 hypothetical proteins were observed in cluster association and could be reliably associated with one of 32 clusters. Several nematode-specific GPCR clades were observed suggesting lineage-specific functional recruitment in response to environment.

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.

Genome-wide comparative gene family classification

Correct classification of genes into gene families is important for understanding gene function and evolution. Although gene families of many species have been resolved both computationally and experimentally with high accuracy, gene family classification in most newly sequenced genomes has not been done with the same high standard. This project has been designed to develop a strategy to effectively and accurately classify gene families across genomes. We first examine and compare the performance of computer programs developed for automated gene family classification. We demonstrate that some programs, including the hierarchical average-linkage clustering algorithm MC-UPGMA and the popular Markov clustering algorithm TRIBE-MCL, can reconstruct manual curation of gene families accurately. However, their performance is highly sensitive to parameter setting, i.e. different gene families require different program parameters for correct resolution. To circumvent the problem of parameterization, we have developed a comparative strategy for gene family classification. This strategy takes advantage of existing curated gene families of reference species to find suitable parameters for classifying genes in related genomes. To demonstrate the effectiveness of this novel strategy, we use TRIBE-MCL to classify chemosensory and ABC transporter gene families in C. elegans and its four sister species. We conclude that fully automated programs can establish biologically accurate gene families if parameterized accordingly. Comparative gene family classification finds optimal parameters automatically, thus allowing rapid insights into gene families of newly sequenced species.

Large synteny blocks revealed between Caenorhabditis elegans and Caenorhabditis briggsae genomes using OrthoCluster

BACKGROUND

Accurate identification of synteny blocks is an important step in comparative genomics towards the understanding of genome architecture and expression. Most computer programs developed in the last decade for identifying synteny blocks have limitations. To address these limitations, we recently developed a robust program called OrthoCluster, and an online database OrthoClusterDB. In this work, we have demonstrated the application of OrthoCluster in identifying synteny blocks between the genomes of Caenorhabditis elegans and Caenorhabditis briggsae, two closely related hermaphrodite nematodes.

RESULTS

Initial identification and analysis of synteny blocks using OrthoCluster enabled us to systematically improve the genome annotation of C. elegans and C. briggsae, identifying 52 potential novel genes in C. elegans, 582 in C. briggsae, and 949 novel orthologous relationships between these two species. Using the improved annotation, we have detected 3,058 perfect synteny blocks that contain no mismatches between C. elegans and C. briggsae. Among these synteny blocks, the majority are mapped to homologous chromosomes, as previously reported. The largest perfect synteny block contains 42 genes, which spans 201.2 kb in Chromosome V of C. elegans. On average, perfect synteny blocks span 18.8 kb in length. When some mismatches (interruptions) are allowed, synteny blocks ("imperfect synteny blocks") that are much larger in size are identified. We have shown that the majority (80%) of the C. elegans and C. briggsae genomes are covered by imperfect synteny blocks. The largest imperfect synteny block spans 6.14 Mb in Chromosome X of C. elegans and there are 11 synteny blocks that are larger than 1 Mb in size. On average, imperfect synteny blocks span 63.6 kb in length, larger than previously reported.

CONCLUSIONS

We have demonstrated that OrthoCluster can be used to accurately identify synteny blocks and have found that synteny blocks between C. elegans and C. briggsae are almost three-folds larger than previously identified.

Copy number variation in the genomes of twelve natural isolates of Caenorhabditis elegans

Background

Copy number variation is an important component of genetic variation in higher eukaryotes. The extent of natural copy number variation in C. elegans is unknown outside of 2 highly divergent wild isolates and the canonical N2 Bristol strain.

Results

We have used array comparative genomic hybridization (aCGH) to detect copy number variation in the genomes of 12 natural isolates of Caenorhabditis elegans. Deletions relative to the canonical N2 strain are more common in these isolates than duplications, and indels are enriched in multigene families on the autosome arms. Among the strains in our study, the Hawaiian and Madeiran strains (CB4856 and JU258) carry the largest number of deletions, followed by the Vancouver strain (KR314). Overall we detected 510 different deletions affecting 1136 genes, or over 5% of the genes in the canonical N2 genome. The indels we identified had a median length of 2.7 kb. Since many deletions are found in multiple isolates, deletion loci were used as markers to derive an unrooted tree to estimate genetic relatedness among the strains.

Conclusion

Copy number variation is extensive in C. elegans, affecting over 5% of the genes in the genome. The deletions we have detected in natural isolates of C. elegans contribute significantly to the number of deletion alleles available to researchers. The relationships between strains are complex and different regions of the genome possess different genealogies due to recombination throughout the natural history of the species, which may not be apparent in studies utilizing smaller numbers of genetic markers.

Divergent mechanisms controlling hypoxic sensitivity and lifespan by the DAF-2/insulin/IGF-receptor pathway

Organisms and their cells vary greatly in their tolerance of low oxygen environments (hypoxia). A delineation of the determinants of hypoxia tolerance is incomplete, despite intense interest for its implications in diseases such as stroke and myocardial infarction. The insulin/IGF-1 receptor (IGFR) signaling pathway controls survival of Caenorhabditis elegans from a variety of stressors including aging, hyperthermia, and hypoxia. daf-2 encodes a C. elegans IGFR homolog whose primary signaling pathway modulates the activity of the FOXO transcription factor DAF-16. DAF-16 regulates the transcription of a large number of genes, some of which have been shown to control aging. To identify genes that selectively regulate hypoxic sensitivity, we compared the whole-organismal transcriptomes of three daf-2 reduction-of-function alleles, all of which are hypoxia resistant, thermotolerant, and long lived, but differ in their rank of severities for these phenotypes. The transcript levels of 172 genes were increased in the most hypoxia resistant daf-2 allele, e1370, relative to the other alleles whereas transcripts from only 10 genes were decreased in abundance. RNAi knockdown of 6 of the 10 genes produced a significant increase in organismal survival after hypoxic exposure as would be expected if down regulation of these genes by the e1370 mutation was responsible for hypoxia resistance. However, RNAi knockdown of these genes did not prolong lifespan. These genes definitively separate the mechanisms of hypoxic sensitivity and lifespan and identify biological strategies to survive hypoxic injury.

Polymorphic segmental duplication in the nematode Caenorhabditis elegans

BACKGROUND

The nematode Caenorhabditis elegans was the first multicellular organism to have its genome fully sequenced. Over the last 10 years since the original publication in 1998, the C. elegans genome has been scrutinized and the last gaps were filled in November 2002, which present a unique opportunity for examining genome-wide segmental duplications.

RESULTS

Here, we performed analysis of the C. elegans genome in search for segmental duplications using a new tool -- OrthoCluster -- we have recently developed. We detected 3,484 duplicated segments -- duplicons -- ranging in size from 234 bp to 108 Kb. The largest pair of duplicons, 108 kb in length located on the left arm of Chromosome V, was further characterized. They are nearly identical at the DNA level (99.7% identity) and each duplicon contains 26 putative protein coding genes. Genotyping of 76 wild-type strains obtained from different labs in the C. elegans community revealed that not all strains contain this duplication. In fact, only 29 strains carry this large segmental duplication, suggesting a very recent duplication event in the C. elegans genome.

CONCLUSION

This report represents the first demonstration that the C. elegans laboratory wild-type N2 strains has acquired large-scale differences.

Evolution of the Caenorhabditis elegans genome

A fundamental problem in genome biology is to elucidate the evolutionary forces responsible for generating nonrandom patterns of genome organization. As the first metazoan to benefit from full-genome sequencing, Caenorhabditis elegans has been at the forefront of research in this area. Studies of genomic patterns, and their evolutionary underpinnings, continue to be augmented by the recent push to obtain additional full-genome sequences of related Caenorhabditis taxa. In the near future, we expect to see major advances with the onset of whole-genome resequencing of multiple wild individuals of the same species. In this review, we synthesize many of the important insights to date in our understanding of genome organization and function that derive from the evolutionary principles made explicit by theoretical population genetics and molecular evolution and highlight fertile areas for future research on unanswered questions in C. elegans genome evolution. We call attention to the need for C. elegans researchers to generate and critically assess nonadaptive hypotheses for genomic and developmental patterns, in addition to adaptive scenarios. We also emphasize the potential importance of evolution in the gonochoristic (female and male) ancestors of the androdioecious (hermaphrodite and male) C. elegans as the source for many of its genomic and developmental patterns.

The Caenorhabditis chemoreceptor gene families

BACKGROUND

Chemoreceptor proteins mediate the first step in the transduction of environmental chemical stimuli, defining the breadth of detection and conferring stimulus specificity. Animal genomes contain families of genes encoding chemoreceptors that mediate taste, olfaction, and pheromone responses. The size and diversity of these families reflect the biology of chemoperception in specific species.

RESULTS

Based on manual curation and sequence comparisons among putative G-protein-coupled chemoreceptor genes in the nematode Caenorhabditis elegans, we identified approximately 1300 genes and 400 pseudogenes in the 19 largest gene families, most of which fall into larger superfamilies. In the related species C. briggsae and C. remanei, we identified most or all genes in each of the 19 families. For most families, C. elegans has the largest number of genes and C. briggsae the smallest number, suggesting changes in the importance of chemoperception among the species. Protein trees reveal family-specific and species-specific patterns of gene duplication and gene loss. The frequency of strict orthologs varies among the families, from just over 50% in two families to less than 5% in three families. Several families include large species-specific expansions, mostly in C. elegans and C. remanei.

CONCLUSION

Chemoreceptor gene families in Caenorhabditis species are large and evolutionarily dynamic as a result of gene duplication and gene loss. These dynamics shape the chemoreceptor gene complements in Caenorhabditis species and define the receptor space available for chemosensory responses. To explain these patterns, we propose the gray pawn hypothesis: individual genes are of little significance, but the aggregate of a large number of diverse genes is required to cover a large phenotype space.

Beta-secretase: structure, function, and evolution

The most popular current hypothesis is that Alzheimer's disease (AD) is caused by aggregates of the amyloid peptide (Abeta), which is generated by cleavage of the Abeta protein precursor (APP) by beta-secretase (BACE-1) followed by gamma-secretase. BACE-1 cleavage is limiting for the production of Abeta, making it a particularly good drug target for the generation of inhibitors that lower Abeta. A landmark discovery in AD was the identification of BACE-1 (a.k.a. Memapsin-2) as a novel class of type I transmembrane aspartic protease. Although BACE-2, a homologue of BACE-1, was quickly identified, follow up studies using knockout mice demonstrated that BACE-1 was necessary and sufficient for most neuronal Abeta generation. Despite the importance of BACE-1 as a drug target, development has been slow due to the incomplete understanding of its function and regulation and the difficulties in developing a brain penetrant drug that can specifically block its large catalytic pocket. This review summarizes the biological properties of BACE-1 and attempts to use phylogenetic perspectives to understand its function. The article also addresses the challenges in discovering a selective drug-like molecule targeting novel mechanisms of BACE-1 regulation.

GPCR genes are preferentially retained after whole genome duplication

One of the most interesting questions in biology is whether certain pathways have been favored during evolution, and if so, what properties could cause such a preference. Due to the lack of experimental evidence, whether select gene families have been preferentially retained over time after duplication in metazoan organisms remains unclear. Here, by syntenic mapping of nonchemosensory G protein-coupled receptor genes (nGPCRs which represent half the receptome for transmembrane signaling) in the vertebrate genomes, we found that, as opposed to the 8–15% retention rate for whole genome duplication (WGD)-derived gene duplicates in the entire genome of pufferfish, greater than 27.8% of WGD-derived nGPCRs which interact with a nonpeptide ligand were retained after WGD in pufferfish Tetraodon nigroviridis. In addition, we show that concurrent duplication of cognate ligand genes by WGD could impose selection of nGPCRs that interact with a polypeptide ligand. Against less than 2.25% probability for parallel retention of a pair of WGD-derived ligands and a pair of cognate receptor duplicates, we found a more than 8.9% retention of WGD-derived ligand-nGPCR pairs–threefold greater than one would surmise. These results demonstrate that gene retention is not uniform after WGD in vertebrates, and suggest a Darwinian selection of GPCR-mediated intercellular communication in metazoan organisms.

Identification of ciliary and ciliopathy genes in Caenorhabditis elegans through comparative genomics

Comparative genomic analysis of three nematode species identifies 93 genes that encode putative components of the ciliated neurons in C. elegans and are subject to the same regulatory control.

Background

The recent availability of genome sequences of multiple related Caenorhabditis species has made it possible to identify, using comparative genomics, similarly transcribed genes in Caenorhabditis elegans and its sister species. Taking this approach, we have identified numerous novel ciliary genes in C. elegans, some of which may be orthologs of unidentified human ciliopathy genes.

Results

By screening for genes possessing canonical X-box sequences in promoters of three Caenorhabditis species, namely C. elegans, C. briggsae and C. remanei, we identified 93 genes (including known X-box regulated genes) that encode putative components of ciliated neurons in C. elegans and are subject to the same regulatory control. For many of these genes, restricted anatomical expression in ciliated cells was confirmed, and control of transcription by the ciliogenic DAF-19 RFX transcription factor was demonstrated by comparative transcriptional profiling of different tissue types and of daf-19(+) and daf-19(-) animals. Finally, we demonstrate that the dye-filling defect of dyf-5(mn400) animals, which is indicative of compromised exposure of cilia to the environment, is caused by a nonsense mutation in the serine/threonine protein kinase gene M04C9.5.

Conclusion

Our comparative genomics-based predictions may be useful for identifying genes involved in human ciliopathies, including Bardet-Biedl Syndrome (BBS), since the C. elegans orthologs of known human BBS genes contain X-box motifs and are required for normal dye filling in C. elegans ciliated neurons.

Rapid birth-death evolution specific to xenobiotic cytochrome P450 genes in vertebrates

Genes vary greatly in their long-term phylogenetic stability and there exists no general explanation for these differences. The cytochrome P450 (CYP450) gene superfamily is well suited to investigating this problem because it is large and well studied, and it includes both stable and unstable genes. CYP450 genes encode oxidase enzymes that function in metabolism of endogenous small molecules and in detoxification of xenobiotic compounds. Both types of enzymes have been intensively studied. My analysis of ten nearly complete vertebrate genomes indicates that each genome contains 50-80 CYP450 genes, which are about evenly divided between phylogenetically stable and unstable genes. The stable genes are characterized by few or no gene duplications or losses in species ranging from bony fish to mammals, whereas unstable genes are characterized by frequent gene duplications and losses (birth-death evolution) even among closely related species. All of the CYP450 genes that encode enzymes with known endogenous substrates are phylogenetically stable. In contrast, most of the unstable genes encode enzymes that function as xenobiotic detoxifiers. Nearly all unstable CYP450 genes in the mouse and human genomes reside in a few dense gene clusters, forming unstable gene islands that arose by recurrent local gene duplication. Evidence for positive selection in amino acid sequence is restricted to these unstable CYP450 genes, and sites of selection are associated with substrate-binding regions in the protein structure. These results can be explained by a general model in which phylogenetically stable genes have core functions in development and physiology, whereas unstable genes have accessory functions associated with unstable environmental interactions such as toxin and pathogen exposure. Unstable gene islands in vertebrates share some functional properties with bacterial genomic islands, though they arise by local gene duplication rather than horizontal gene transfer.

Comparative genomics and adaptive selection of the ATP-binding-cassette gene family in caenorhabditis species

ABC transporters constitute one of the largest gene families in all species. They are mostly involved in transport of substrates across membranes. We have previously demonstrated that the Caenorhabditis elegans ABC family shows poor one-to-one gene orthology with other distant model organisms. To address the evolution dynamics of this gene family among closely related species, we carried out a comparative analysis of the ABC family among the three nematode species C. elegans, C. briggsae, and C. remanei. In contrast to the previous observations, the majority of ABC genes in the three species were found in orthologous trios, including many tandemly duplicated ABC genes, indicating that the gene duplication took place before speciation. Species-specific expansions of ABC members are rare and mostly observed in subfamilies A and B. C. briggsae and C. remanei orthologous ABC genes tend to cluster on trees, with those of C. elegans as an outgroup, consistent with their proposed species phylogeny. Comparison of intron/exon structures of the highly conserved ABCE subfamily members also indicates a closer relationship between C. briggsae and C. remanei than between either of these species and C. elegans. A comparison between insect and mammalian species indicates lineage-specific duplications or deletions of ABC genes, while the family size remains relatively constant. Sites undergoing positive selection within subfamily D, which are implicated in very-long-chain fatty acid transport, were identified. The evolution of these sites might be driven by the changes in food source with time.

Efficient high-resolution deletion discovery in Caenorhabditis elegans by array comparative genomic hybridization

We have developed array Comparative Genomic Hybridization for Caenorhabditis elegans as a means of screening for novel induced deletions in this organism. We designed three microarrays consisting of overlapping 50-mer probes to annotated exons and micro-RNAs, the first with probes to chromosomes X and II, the second with probes to chromosome II alone, and a third to the entire genome. These arrays were used to reliably detect both a large (50 kb) multigene deletion and a small (1 kb) single-gene deletion in homozygous and heterozygous samples. In one case, a deletion breakpoint was resolved to fewer than 50 bp. In an experiment designed to identify new mutations we used the X:II and II arrays to detect deletions associated with lethal mutants on chromosome II. One is an 8-kb deletion targeting the ast-1 gene on chromosome II and another is a 141-bp deletion in the gene C06A8.1. Others span large sections of the chromosome, up to >750 kb. As a further application of array Comparative Genomic Hybridization in C. elegans we used the whole-genome array to detect the extensive natural gene content variation (almost 2%) between the N2 Bristol strain and the strain CB4856, a strain isolated in Hawaii and JU258, a strain isolated in Madeira.

Comparative chemosensation from receptors to ecology

Odour perception is initiated by specific interactions between odorants and a large repertoire of receptors in olfactory neurons. During the past few years, considerable progress has been made in tracing olfactory perception from the odorant receptor protein to the activity of olfactory neurons to higher processing centres and, ultimately, to behaviour. The most complete picture is emerging for the simplest olfactory system studied--that of the fruitfly Drosophila melanogaster. Comparison of rodent, insect and nematode olfaction reveals surprising differences and unexpected similarities among chemosensory systems.

Unusual DNA structures associated with germline genetic activity in Caenorhabditis elegans

We describe a surprising long-range periodicity that underlies a substantial fraction of C. elegans genomic sequence. Extended segments (up to several hundred nucleotides) of the C. elegans genome show a strong bias toward occurrence of AA/TT dinucleotides along one face of the helix while little or no such constraint is evident on the opposite helical face. Segments with this characteristic periodicity are highly overrepresented in intron sequences and are associated with a large fraction of genes with known germline expression in C. elegans. In addition to altering the path and flexibility of DNA in vitro, sequences of this character have been shown by others to constrain DNA::nucleosome interactions, potentially producing a structure that could resist the assembly of highly ordered (phased) nucleosome arrays that have been proposed as a precursor to heterochromatin. We propose a number of ways that the periodic occurrence of An/Tn clusters could reflect evolution and function of genes that express in the germ cell lineage of C. elegans.

Genomics in C. elegans: so many genes, such a little worm

The Caenorhabditis elegans genome sequence is now complete, fully contiguous telomere to telomere and totaling 100,291,840 bp. The sequence has catalyzed the collection of systematic data sets and analyses, including a curated set of 19,735 protein-coding genes--with >90% directly supported by experimental evidence--and >1300 noncoding RNA genes. High-throughput efforts are under way to complete the gene sets, along with studies to characterize gene expression, function, and regulation on a genome-wide scale. The success of the worm project has had a profound effect on genome sequencing and on genomics more broadly. We now have a solid platform on which to build toward the lofty goal of a true molecular understanding of worm biology with all its implications including those for human health.

How sensitive is a nose?

Odor sensitivity may not be due to odor-receptor (OR) binding but rather may be due to emergent properties of transduction pathways and the anatomical convergence of olfactory neurons. A recent study suggests that odor-OR interactions are brief and infrequently activate heterotrimeric GTP-binding proteins (G proteins); in contrast, visual receptors have long-lasting activation states and activate many G proteins. These differences may reflect strategies that evolved to accommodate very different signals, and the mechanisms described might be applicable for receptors across phyla. However, whereas visual receptors (rhodopsin) appeared before protostome-deuterostome separation, ORs may be independently derived in different phyla. Alternatively, phylum-distinct ORs may share common ancestry but be influenced by diversifying selection. Phylum-distinct ORs may imply phylum-specific OR mechanisms, whereas common ancestry may imply common mechanisms. Nonetheless, most animals detect a similar repertoire of olfactory signals, and OR mechanisms may be convergent on those signals independent of receptor relatedness. Thus, recent insights into the molecular characteristics of odor perception in frogs may well be relevant to such processes as how mosquitoes detect host odors for a malaria-transmitting blood meal.

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.

Analysis of homologous gene clusters in Caenorhabditis elegans reveals striking regional cluster domains

An algorithm for detecting local clusters of homologous genes was applied to the genome of Caenorhabditis elegans. Clusters of two or more homologous genes are abundant, totaling 1391 clusters containing 4607 genes, over one-fifth of all genes in C. elegans. Cluster genes are distributed unevenly in the genome, with the large majority located on autosomal chromosome arms, regions characterized by higher genetic recombination and more repeat sequences than autosomal centers and the X chromosome. Cluster genes are transcribed at much lower levels than average and very few have gross phenotypes as assayed by RNAi-mediated reduction of function. The molecular identity of cluster genes is unusual, with a preponderance of nematode-specific gene families that encode putative secreted and transmembrane proteins, and enrichment for genes implicated in xenobiotic detoxification and innate immunity. Gene clustering in Drosophila melanogaster is also substantial and the molecular identity of clustered genes follows a similar pattern. I hypothesize that autosomal chromosome arms in C. elegans undergo frequent local gene duplication and that these duplications support gene diversification and rapid evolution in response to environmental challenges. Although specific gene clusters have been documented in C. elegans, their abundance, genomic distribution, and unusual molecular identities were previously unrecognized.

Automatic gene collection system for genome-scale overview of G-protein coupled receptors in eukaryotes

We have developed an automatic system for identifying GPCR (G-protein coupled receptor) genes from various kinds of genomes, which is finally deposited in the SEVENS database (http://sevens.cbrc.jp/), by integrating such software as a gene finder, a sequence alignment tool, a motif and domain assignment tool, and a transmembrane helix predictor. SEVENS enables us to perform a genome-scale overview of the "GPCR universe" using sequences that are identified with high accuracy (99.4% sensitivity and 96.6% specificity). Using this system, we surveyed the complete genomes of 7 eukaryotes and 224 prokaryotes, and found that there are 4 to 1016 GPCR genes in the 7 eukaryotes, and only a total of 16 GPCR genes in all the prokaryotes. Our preliminary results indicate that 11 subfamilies of the Class A family, the Class 2(B) family, the Class 3(C) family and the fz/smo family are commonly found among human, fly, and nematode genomes. We also analyzed the chromosomal locations of the GPCR genes with the Kolmogorov-Smirnov test, and found that species-specific families, such as olfactory, taste, and chemokine receptors in human and nematode chemoreceptor in worm, tend to form clusters extensively, whereas no significant clusters were detected in fly and plant genomes.