WormBase: a comprehensive data resource for Caenorhabditis biology and genomics

ADBP-1 regulates ADR-2 nuclear localization to control editing substrate selection

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by ADAR enzymes, is a prevalent and conserved RNA modification. While A-to-I RNA editing is essential in mammals, in Caenorhabditis elegans, it is not, making them invaluable for RNA editing research. In C. elegans, ADR-2 is the sole catalytic A-to-I editing enzyme, and ADR-1 is an RNA editing regulator. ADAR localization is well-studied in humans but not well-established in C. elegans. In this study, we examine the cellular and tissue-specific localization of ADR-2. We show that while ADR-2 is present in most cells in the embryo, at later developmental stages, its expression is both tissue- and cell-type-specific. Additionally, both ADARs are mainly in the nucleus. ADR-2 is adjacent to the chromosomes during the cell cycle. We show that the nuclear localization of endogenous ADR-2 depends on ADBP-1, not ADR-1. In adbp-1 mutant worms, ADR-2 is mislocalized, while ADR-1 is not, leading to decreased editing levels and de-novo editing, mostly in exons, suggesting that ADR-2 is also functional in the cytoplasm. Besides, mutated ADBP-1 affects gene expression. Furthermore, we show that ADR-2 targets adenosines with different surrounding nucleotides in exons and introns. Our findings indicate that ADR-2 cellular localization is highly regulated and affects its function.

ADBP-1 regulates ADR-2 nuclear localization to control editing substrate selection

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by ADAR enzymes, is a prevalent and conserved RNA modification. While A-to-I RNA editing is essential in mammals, in Caenorhabditis elegans , it is not, making them invaluable for RNA editing research. In C. elegans , ADR-2 is the sole catalytic A-to-I editing enzyme, and ADR-1 is an RNA editing regulator. ADAR localization is well-studied in humans but not well-established in C. elegans . In this study, we examine the cellular and tissue-specific localization of ADR-2. We show that while ADR-2 is present in most cells in the embryo, at later developmental stages, its expression is both tissue- and cell-type-specific. Additionally, both ADARs are mainly in the nucleus. ADR-2 is adjacent to the chromosomes during the cell cycle. We show that the nuclear localization of endogenous ADR-2 depends on ADBP-1, not ADR-1. In adbp-1 mutant worms, ADR-2 is mislocalized, while ADR-1 is not, leading to decreased editing levels and de-novo editing, mostly in exons, suggesting that ADR-2 is also functional in the cytoplasm. Besides, mutated ADBP-1 affects gene expression. Furthermore, we show that ADR-2 targets adenosines with different surrounding nucleotides in exons and introns. Our findings indicate that ADR-2 cellular localization is highly regulated and affects its function.

WormBase 2024: status and transitioning to Alliance infrastructure

WormBase has been the major repository and knowledgebase of information about the genome and genetics of Caenorhabditis elegans and other nematodes of experimental interest for over 2 decades. We have 3 goals: to keep current with the fast-paced C. elegans research, to provide better integration with other resources, and to be sustainable. Here, we discuss the current state of WormBase as well as progress and plans for moving core WormBase infrastructure to the Alliance of Genome Resources (the Alliance). As an Alliance member, WormBase will continue to interact with the C. elegans community, develop new features as needed, and curate key information from the literature and large-scale projects.

A bioinformatics approach to elucidate conserved genes and pathways in C. elegans as an animal model for cardiovascular research

Cardiovascular disease (CVD) is a collective term for disorders of the heart and blood vessels. The molecular events and biochemical pathways associated with CVD are difficult to study in clinical settings on patients and in vitro conditions. Animal models play a pivotal and indispensable role in CVD research. Caenorhabditis elegans, a nematode species, has emerged as a prominent experimental organism widely utilized in various biomedical research fields. However, the specific number of CVD-related genes and pathways within the C. elegans genome remains undisclosed to date, limiting its in-depth utilization for investigations. In the present study, we conducted a comprehensive analysis of genes and pathways related to CVD within the genomes of humans and C. elegans through a systematic bioinformatic approach. A total of 1113 genes in C. elegans orthologous to the most significant CVD-related genes in humans were identified, and the GO terms and pathways were compared to study the pathways that are conserved between the two species. In order to infer the functions of CVD-related orthologous genes in C. elegans, a PPI network was constructed. Orthologous gene PPI network analysis results reveal the hubs and important KRs: pmk-1, daf-21, gpb-1, crh-1, enpl-1, eef-1G, acdh-8, hif-1, pmk-2, and aha-1 in C. elegans. Modules were identified for determining the role of the orthologous genes at various levels in the created network. We also identified 9 commonly enriched pathways between humans and C. elegans linked with CVDs that include autophagy (animal), the ErbB signaling pathway, the FoxO signaling pathway, the MAPK signaling pathway, ABC transporters, the biosynthesis of unsaturated fatty acids, fatty acid metabolism, glutathione metabolism, and metabolic pathways. This study provides the first systematic genomic approach to explore the CVD-associated genes and pathways that are present in C. elegans, supporting the use of C. elegans as a prominent animal model organism for cardiovascular diseases.

Conserved Cardiovascular Network: Bioinformatics Insights into Genes and Pathways for Establishing Caenorhabditis elegans as an Animal Model for Cardiovascular Diseases

Cardiovascular disease (CVD) is a collective term for disorders of the heart and blood vessels. The molecular events and biochemical pathways associated with CVD are difficult to study in clinical settings on patients and in vitro conditions. Animal models play a pivotal and indispensable role in cardiovascular disease (CVD) research. Caenorhabditis elegans , a nematode species, has emerged as a prominent experimental organism widely utilised in various biomedical research fields. However, the specific number of CVD-related genes and pathways within the C. elegans genome remains undisclosed to date, limiting its in-depth utilisation for investigations. In the present study, we conducted a comprehensive analysis of genes and pathways related to CVD within the genomes of humans and C. elegans through a systematic bioinformatic approach. A total of 1113 genes in C. elegans orthologous to the most significant CVD-related genes in humans were identified, and the GO terms and pathways were compared to study the pathways that are conserved between the two species. In order to infer the functions of CVD-related orthologous genes in C. elegans, a PPI network was constructed. Orthologous gene PPI network analysis results reveal the hubs and important KRs: pmk-1, daf-21, gpb-1, crh-1, enpl-1, eef-1G, acdh-8, hif-1, pmk-2, and aha-1 in C. elegans. Modules were identified for determining the role of the orthologous genes at various levels in the created network. We also identified 9 commonly enriched pathways between humans and C. elegans linked with CVDs that include autophagy (animal), the ErbB signalling pathway, the FoxO signalling pathway, the MAPK signalling pathway, ABC transporters, the biosynthesis of unsaturated fatty acids, fatty acid metabolism, glutathione metabolism, and metabolic pathways. This study provides the first systematic genomic approach to explore the CVD-associated genes and pathways that are present in C. elegans, supporting the use of C. elegans as a prominent animal model organism for cardiovascular diseases.

Ivermectin-induced gene expression changes in adult Parascaris univalens and Caenorhabditis elegans: a comparative approach to study anthelminthic metabolism and resistance in vitro

Background

The nematode Parascaris univalens is one of the most prevalent parasitic pathogens infecting horses but anthelmintic resistance undermines treatment approaches. The molecular mechanisms underlying drug activity and resistance remain poorly understood in this parasite since experimental in vitro models are lacking. The aim of this study was to evaluate the use of Caenorhabditis elegans as a model for P. univalens drug metabolism/resistance studies by a comparative gene expression approach after in vitro exposure to the anthelmintic drug ivermectin (IVM).

Methods

Twelve adult P. univalens worms in groups of three were exposed to ivermectin (IVM, 10^–13 M, 10^–11 M, 10^–9 M) or left unexposed for 24 h at 37 °C, and total RNA, extracted from the anterior end of the worms, was sequenced using Illumina NovaSeq. Differentially expressed genes (DEGs) involved in metabolism, transportation, or gene expression with annotated Caernorhabditis elegans orthologues were identified as candidate genes to be involved in IVM metabolism/resistance. Similarly, groups of 300 adult C. elegans worms were exposed to IVM (10^–9 M, 10^–8 M and 10^–7 M) or left unexposed for 4 h at 20 °C. Quantitative RT-PCR of RNA extracted from the C. elegans worm pools was used to compare against the expression of selected P. univalens candidate genes after drug treatment.

Results

After IVM exposure, 1085 DEGs were found in adult P. univalens worms but the relative gene expression changes were small and large variabilities were found between different worms. Fifteen of the DEGs were chosen for further characterization in C. elegans after comparative bioinformatics analyses. Candidate genes, including the putative drug target lgc-37, responded to IVM in P. univalens, but marginal to no responses were observed in C. elegans despite dose-dependent behavioral effects observed in C. elegans after IVM exposure. Thus, the overlap in IVM-induced gene expression in this small set of genes was minor in adult worms of the two nematode species.

Conclusion

This is the first time to our knowledge that a comparative gene expression approach has evaluated C. elegans as a model to understand IVM metabolism/resistance in P. univalens. Genes in P. univalens adults that responded to IVM treatment were identified. However, identifying conserved genes in P. univalens and C. elegans involved in IVM metabolism/resistance by comparing gene expression of candidate genes proved challenging. The approach appears promising but was limited by the number of genes studied (n = 15). Future studies comparing a larger number of genes between the two species may result in identification of additional candidate genes involved in drug metabolism and/or resistance.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05260-4.

Identification of Essential Genes in Caenorhabditis elegans with Lethal Mutations Maintained by Genetic Balancers

Genetic balancer systems, which allow effective capture and maintenance of lethal mutations stably, play an important role in identifying essential genes. Whole-genome sequencing (WGS) followed by bioinformatics analysis, combined with genetic mapping data analysis, allows for an efficient and economical means of identifying genomic mutations in essential genes. Using this approach, we successfully identified 104 essential genes on ChrI, ChrIII, and ChrV in C. elegans. In this report, we described a protocol that sequences the genome of prebalanced Caenorhabditis elegans (C. elegans) strains to carry lethal mutations and identifies candidate causal mutations and candidate essential genes using a robust bioinformatics procedure.

Towards a reporting guideline for developmental and reproductive toxicology testing in C. elegans and other nematodes

Implementation of reliable methodologies allowing Reduction, Refinement, and Replacement (3Rs) of animal testing is a process that takes several decades and is still not complete. Reliable methods are essential for regulatory hazard assessment of chemicals where differences in test protocol can influence the test outcomes and thus affect the confidence in the predictive value of the organisms used as an alternative for mammals. Although test guidelines are common for mammalian studies, they are scarce for non-vertebrate organisms that would allow for the 3Rs of animal testing. Here, we present a set of 30 reporting criteria as the basis for such a guideline for Developmental and Reproductive Toxicology (DART) testing in the nematode Caenorhabditis elegans. Small organisms like C. elegans are upcoming in new approach methodologies for hazard assessment; thus, reliable and robust test protocols are urgently needed. A literature assessment of the fulfilment of the reporting criteria demonstrates that although studies describe methodological details, essential information such as compound purity and lot/batch number or type of container is often not reported. The formulated set of reporting criteria for C. elegans testing can be used by (i) researchers to describe essential experimental details (ii) data scientists that aggregate information to assess data quality and include data in aggregated databases (iii) regulators to assess study data for inclusion in regulatory hazard assessment of chemicals.

The genome-wide rate and spectrum of EMS-induced heritable mutations in the microcrustacean Daphnia: on the prospect of forward genetics

Forward genetic screening using the alkylating mutagen ethyl methanesulfonate (EMS) is an effective method for identifying phenotypic mutants of interest, which can be further genetically dissected to pinpoint the causal genetic mutations. An accurate estimate of the rate of EMS-induced heritable mutations is fundamental for determining the mutant sample size of a screening experiment that aims to saturate all the genes in a genome with mutations. This study examines the genome-wide EMS-induced heritable base-substitutions in three species of the freshwater microcrustacean Daphnia to help guide screening experiments. Our results show that the 10 mM EMS treatment induces base substitutions at an average rate of 1.17 × 10^-6/site/generation across the three species, whereas a significantly higher average mutation rate of 1.75 × 10^-6 occurs at 25 mM. The mutation spectrum of EMS-induced base substitutions at both concentration is dominated by G:C to A:T transitions. Furthermore, we find that female Daphnia exposed to EMS (F₀ individuals) can asexually produce unique mutant offspring (F₁) for at least 3 consecutive broods, suggestive of multiple broods as F₁ mutants. Lastly, we estimate that about 750 F₁s are needed for all genes in the Daphnia genome to be mutated at least once with a 95% probability. We also recommend 4-5 F₂s should be collected from each F₁ mutant through sibling crossing so that all induced mutations could appear in the homozygous state in the F₂ population at 70-80% probability.

Significance of probiotics in remodeling the gut consortium to enhance the immunity of Caenorhabditis elegans

In the recent past, Caenorhabditis elegans has emerged as one of the leading nematode models for studying host-microbe interactions on molecular, cellular, or organismal levels. In general, morphological and functional similarities of the gut of C. elegans with respect to that of human has brought in speculations on the study of the intestinal microbiota. On the other hand, probiotics have proved their efficacy in metabolism, development, and pathogenesis thereby inducing an immune response in C. elegans. Nurturing C. elegans with probiotics has led to immunomodulatory effects in the intestinal microbiota, proposing C. elegans as one of the in vivo screening criteria to select potential probiotic bacteria for host health-promoting factors. The major prospect of these probiotics is to exert longevity toward the host in diverse environmental conditions. The extent of research on probiotic metabolism has shed light on mechanisms of the immunomodulatory effect exerted by the nematode model. This review discusses various aspects of the effects of probiotics in improving the health and mechanisms involved in conferring immunity in C. elegans.

Deep Collaborative Filtering for Prediction of Disease Genes

Accurate prioritization of potential disease genes is a fundamental challenge in biomedical research. Various algorithms have been developed to solve such problems. Inductive Matrix Completion (IMC) is one of the most reliable models for its well-established framework and its superior performance in predicting gene-disease associations. However, the IMC method does not hierarchically extract deep features, which might limit the quality of recovery. In this case, the architecture of deep learning, which obtains high-level representations and handles noises and outliers presented in large-scale biological datasets, is introduced into the side information of genes in our Deep Collaborative Filtering (DCF) model. Further, for lack of negative examples, we also exploit Positive-Unlabeled (PU) learning formulation to low-rank matrix completion. Our approach achieves substantially improved performance over other state-of-the-art methods on diseases from the Online Mendelian Inheritance in Man (OMIM) database. Our approach is 10 percent more efficient than standard IMC in detecting a true association, and significantly outperforms other alternatives in terms of the precision-recall metric at the top-k predictions. Moreover, we also validate the disease with no previously known gene associations and newly reported OMIM associations. The experimental results show that DCF is still satisfactory for ranking novel disease phenotypes as well as mining unexplored relationships. The source code and the data are available at https://github.com/xzenglab/DCF.

Caenorhabditis elegans as a model system for human diseases

The nematode Caenorhabditis elegans offers unique advantages that enable a comprehensive delineation of the cellular and molecular mechanisms underlying devastating human pathologies such as stroke, ischemia and age-associated neurodegenerative disorders. Genetic models of human diseases that closely simulate several disease-related phenotypes have been established in the worm. These models allow the implementation of multidisciplinary approaches, in addition to large-scale genetic and pharmacological screenings, designed to elucidate the molecular mechanisms mediating pathogenesis and to identify targets and drugs for emergent therapeutic interventions. Such strategies have already provided valuable insights, highly relevant to human health and quality of life. This article considers the potential of C. elegans as a versatile platform for systematic dissection of the molecular basis of human disease, focusing on neurodegenerative disorders.

Transcriptional variation and divergence of host-finding behaviour in Steinernema carpocapsae infective juveniles

Background

Steinernema carpocapsae is an entomopathogenic nematode that employs nictation and jumping behaviours to find potential insect hosts. Here we aimed to investigate the transcriptional basis of variant host-finding behaviours in the infective juvenile (IJ) stage of three S. carpocapsae strains (ALL, Breton and UK1), with a focus on neuronal genes known to influence behaviour in other nematode species. Identifying gene expression changes that correlate with variant host-finding behaviours will further our understanding of nematode biology.

Results

RNA-seq analysis revealed that whilst up to 28% of the S. carpocapsae transcriptome was differentially expressed (P < 0.0001) between strains, remarkably few of the most highly differentially expressed genes (> 2 log2 fold change, P < 0.0001) were from neuronal gene families. S. carpocapsae Breton displays increased chemotaxis toward the laboratory host Galleria mellonella, relative to the other strains. This correlates with the up-regulation of four srsx chemosensory GPCR genes, and a sodium transporter gene, asic-2, relative to both ALL and UK1 strains. The UK1 strain exhibits a decreased nictation phenotype relative to ALL and Breton strains, which correlates with co-ordinate up-regulation of neuropeptide like protein 36 (nlp-36), and down-regulation of an srt family GPCR gene, and a distinct asic-2-like sodium channel paralogue. To further investigate the link between transcriptional regulation and behavioural variation, we sequenced microRNAs across IJs of each strain. We have identified 283 high confidence microRNA genes, yielding 321 predicted mature microRNAs in S. carpocapsae, and find that up to 36% of microRNAs are differentially expressed (P < 0.0001) between strains. Many of the most highly differentially expressed microRNAs (> 2 log2 fold, P < 0.0001) are predicted to regulate a variety of neuronal genes that may contribute to variant host-finding behaviours. We have also found evidence for differential gene isoform usage between strains, which alters predicted microRNA interactions, and could contribute to the diversification of behaviour.

Conclusions

These data provide insight to the transcriptional basis of behavioural variation in S. carpocapsae, supporting efforts to understand the molecular basis of complex behaviours in nematodes.

WormBase: A Model Organism Database

WormBase is an open-access model-organism database that provides current and accurate genetic information of C. elegans and related nematodes. Users can search WormBase by several fields, including a gene or human disease. A special feature of the database is the inclusion of micropublications, peer-reviewed data that may go unpublished in traditional venues.

Prediction and validation of the structural features of Ov58GPCR, an immunogenic determinant of Onchocerca volvulus

Onchocerciasis is a severely debilitating yet neglected tropical disease (NTD) that creates social stigma, generates and perpetuates poverty, and leads ultimately in some cases to irreversible unilateral or bilateral blindness if untreated. Consequently, the disease is a major impediment to socioeconomic development. Many control programs have been launched for the disease with moderate successes achieved. This mitigated hit is partially due to the lingering need for reliable, non-invasive and easily applicable tools for mapping endemic regions and post-elimination surveillance. In this work, bioinformatics analyses combined with immunological assays were applied in a bid to develop potential tools for diagnosis and assessing the success of drug treatment programs. We report that (i) the O. volvulus antigen, Ov58GPCR is a G-protein coupled receptor (GPCR) conserved in related nematodes, (ii) synthetic peptides predicted to be in the extracellular domain (ECD) of Ov58GPCR are indeed immunogenic epitopes in actively-infected individuals, (iii) synthetic peptide cocktails discriminate between actively-infected individuals, treated individuals and healthy African controls, (iv) polyclonal antibodies against one of the peptides or against the bacterially-expressed ECD reacted specifically with the native antigen of O. volvulus total and surface extracts, (v) Ov58GPCR is transcribed in both larvae and adult parasite stages, (vi) IgG and IgE responses to the recombinant ECD decline with ivermectin treatment. All these findings suggest that the extracellular domain and synthetic peptides of Ov58GPCR, as well as the specific immune response generated could be harnessed in the context of disease diagnosis and surveillance.

QsRNA-seq: a method for high-throughput profiling and quantifying small RNAs

The ability to profile and quantify small non-coding RNAs (sRNAs), specifically microRNAs (miRNAs), using high-throughput sequencing is challenging because of their small size. We developed QsRNA-seq, a method for preparation of sRNA libraries for high-throughput sequencing that overcomes this difficulty by enabling a gel-free separation of fragments shorter than 100 nt that differ only by 20 nt in length. The method allows the use of unique molecular identifiers for quantification and is more amenable to automation than gel-based methods. We show that QsRNA-seq gives very accurate, comprehensive, and reproducible results by looking at miRNAs in Caenorhabditis elegans embryos and larvae.

ODG: Omics database generator - a tool for generating, querying, and analyzing multi-omics comparative databases to facilitate biological understanding

BACKGROUND

Rapid generation of omics data in recent years have resulted in vast amounts of disconnected datasets without systemic integration and knowledge building, while individual groups have made customized, annotated datasets available on the web with few ways to link them to in-lab datasets. With so many research groups generating their own data, the ability to relate it to the larger genomic and comparative genomic context is becoming increasingly crucial to make full use of the data.

RESULTS

The Omics Database Generator (ODG) allows users to create customized databases that utilize published genomics data integrated with experimental data which can be queried using a flexible graph database. When provided with omics and experimental data, ODG will create a comparative, multi-dimensional graph database. ODG can import definitions and annotations from other sources such as InterProScan, the Gene Ontology, ENZYME, UniPathway, and others. This annotation data can be especially useful for studying new or understudied species for which transcripts have only been predicted, and rapidly give additional layers of annotation to predicted genes. In better studied species, ODG can perform syntenic annotation translations or rapidly identify characteristics of a set of genes or nucleotide locations, such as hits from an association study. ODG provides a web-based user-interface for configuring the data import and for querying the database. Queries can also be run from the command-line and the database can be queried directly through programming language hooks available for most languages. ODG supports most common genomic formats as well as generic, easy to use tab-separated value format for user-provided annotations.

CONCLUSIONS

ODG is a user-friendly database generation and query tool that adapts to the supplied data to produce a comparative genomic database or multi-layered annotation database. ODG provides rapid comparative genomic annotation and is therefore particularly useful for non-model or understudied species. For species for which more data are available, ODG can be used to conduct complex multi-omics, pattern-matching queries.

The Skp1 Homologs SKR-1/2 Are Required for the Caenorhabditis elegans SKN-1 Antioxidant/Detoxification Response Independently of p38 MAPK

SKN-1/Nrf are the primary antioxidant/detoxification response transcription factors in animals and they promote health and longevity in many contexts. SKN-1/Nrf are activated by a remarkably broad-range of natural and synthetic compounds and physiological conditions. Defining the signaling mechanisms that regulate SKN-1/Nrf activation provides insights into how cells coordinate responses to stress. Nrf2 in mammals is regulated in part by the redox sensor repressor protein named Keap1. In C. elegans, the p38 MAPK cascade in the intestine activates SKN-1 during oxidative stress by promoting its nuclear accumulation. Interestingly, we find variation in the kinetics of p38 MAPK activation and tissues with SKN-1 nuclear accumulation among different pro-oxidants that all trigger strong induction of SKN-1 target genes. Using genome-wide RNAi screening, we identify new genes that are required for activation of the core SKN-1 target gene gst-4 during exposure to the natural pro-oxidant juglone. Among 10 putative activators identified in this screen was skr-1/2, highly conserved homologs of yeast and mammalian Skp1, which function to assemble protein complexes. Silencing of skr-1/2 inhibits induction of SKN-1 dependent detoxification genes and reduces resistance to pro-oxidants without decreasing p38 MAPK activation. Global transcriptomics revealed strong correlation between genes that are regulated by SKR-1/2 and SKN-1 indicating a high degree of specificity. We also show that SKR-1/2 functions upstream of the WD40 repeat protein WDR-23, which binds to and inhibits SKN-1. Together, these results identify a novel p38 MAPK independent signaling mechanism that activates SKN-1 via SKR-1/2 and involves WDR-23.

Glucose or Altered Ceramide Biosynthesis Mediate Oxygen Deprivation Sensitivity Through Novel Pathways Revealed by Transcriptome Analysis in Caenorhabditis elegans

Individuals with type 2 diabetes display metabolic abnormalities, such as hyperglycemia, increased free fatty acids, insulin resistance, and altered ceramide levels, that contribute to vascular dysfunctions and compromised oxygen delivery. Caenorhabditis elegans fed a glucose-supplemented diet or with altered ceramide metabolism, due to a hyl-2 mutation, are sensitive to oxygen deprivation (anoxia). Our experiments showed that the combination of these factors further decreased the anoxia survival. RNA-sequencing analysis was performed to assess how a glucose-supplemented diet and/or a hyl-2 mutation altered the transcriptome. Comparison analysis of transcripts associated with anoxia-sensitive animals [hyl-2(tm2031) mutation or a glucose diet] revealed 199 common transcripts encoded by genes with known or predicted functions involving innate immunity, cuticle function (collagens), or xenobiotic and endobiotic phase I and II detoxification system. Use of RNA interference (RNAi) to target gene products of the xenobiotic and endobiotic phase I and II detoxification system (UDP-glycosyltransferase and Cytochrome p450 genes; ugt-15, ugt-18, ugt-19, ugt-41, ugt-63, cyp-13A12, cyp-25A1, and cyp-33C8) increased anoxia survival in wild-type animals fed a standard diet. Anoxia sensitivity of the hyl-2(tm2031) animals was suppressed by RNAi of cyp-25A1 or cyp-33C8 genes. A glucose diet fed to the P0 hermaphrodite decreased the anoxia survival of its F1 embryos; however, the RNAi of ugt-63 and cyp-33C8 suppressed anoxia sensitivity. These studies provide evidence that the detoxification system impacts oxygen deprivation responses and that C. elegans can be used to model the conserved detoxification system.

Large-Scale Functional RNAi Screen in C. elegans Identifies TGF-β and Notch Signaling Pathways as Modifiers of CACNA1A

Variants in CACNA1A that encodes the pore-forming α₁-subunit of human voltage-gated Cav2.1 (P/Q-type) Ca²⁺ channels cause several autosomal-dominant neurologic disorders, including familial hemiplegic migraine type 1, episodic ataxia type 2, and spinocerebellar ataxia type 6. To identify modifiers of incoordination in movement disorders, we performed a large-scale functional RNAi screen, using the Caenorhabditis elegans strain CB55, which carries a truncating mutation in the unc-2 gene, the worm ortholog for the human CACNA1A. The screen was carried out by the feeding method in 96-well liquid culture format, using the ORFeome v1.1 feeding library, and time-lapse imaging of worms in liquid culture was used to assess changes in thrashing behavior. We looked for genes that, when silenced, either ameliorated the slow and uncoordinated phenotype of unc-2, or interacted to produce a more severe phenotype. Of the 350 putative hits from the primary screen, 37 genes consistently showed reproducible results. At least 75% of these are specifically expressed in the C. elegans neurons. Functional network analysis and gene ontology revealed overrepresentation of genes involved in development, growth, locomotion, signal transduction, and vesicle-mediated transport. We have expanded the functional network of genes involved in neurodegeneration leading to cerebellar ataxia related to unc-2/CACNA1A, further confirming the involvement of the transforming growth factor β pathway and adding a novel signaling cascade, the Notch pathway.

Glucose induces sensitivity to oxygen deprivation and modulates insulin/IGF-1 signaling and lipid biosynthesis in Caenorhabditis elegans

Diet is a central environmental factor that contributes to the phenotype and physiology of individuals. At the root of many human health issues is the excess of calorie intake relative to calorie expenditure. For example, the increasing amount of dietary sugars in the human diet is contributing to the rise of obesity and type 2 diabetes. Individuals with obesity and type 2 diabetes have compromised oxygen delivery, and thus it is of interest to investigate the impact a high-sugar diet has on oxygen deprivation responses. By utilizing the Caenorhabditis elegans genetic model system, which is anoxia tolerant, we determined that a glucose-supplemented diet negatively impacts responses to anoxia and that the insulin-like signaling pathway, through fatty acid and ceramide synthesis, modulates anoxia survival. Additionally, a glucose-supplemented diet alters lipid localization and initiates a positive chemotaxis response. Use of RNA-sequencing analysis to compare gene expression responses in animals fed either a standard or glucose-supplemented diet revealed that glucose impacts the expression of genes involved with multiple cellular processes including lipid and carbohydrate metabolism, stress responses, cell division, and extracellular functions. Several of the genes we identified show homology to human genes that are differentially regulated in response to obesity or type 2 diabetes, suggesting that there may be conserved gene expression responses between C. elegans fed a glucose-supplemented diet and a diabetic and/or obesity state observed in humans. These findings support the utility of the C. elegans model for understanding the molecular mechanisms regulating dietary-induced metabolic diseases.

DAF-16/FOXO and EGL-27/GATA promote developmental growth in response to persistent somatic DNA damage

Genome maintenance defects cause complex disease phenotypes characterized by developmental failure, cancer susceptibility, and premature aging. It remains poorly understood how DNA damage responses function during organismal development and maintain tissue functionality when DNA damage accumulates with aging. Here we show that the FoxO transcription factor DAF-16 is activated in response to DNA damage during development while the DNA damage responsiveness of DAF-16 declines with aging. We find that in contrast to its established role in mediating starvation arrest, DAF-16 alleviates DNA damage-induced developmental arrest and even in the absence of DNA repair promotes developmental growth and enhances somatic tissue functionality. We demonstrate that the GATA transcription factor EGL-27 co-regulates DAF-16 target genes in response to DNA damage and together with DAF-16 promotes developmental growth. We propose that EGL-27/GATA activity specifies DAF-16 mediated DNA damage responses to enable developmental progression and to prolong tissue functioning when DNA damage persists.

A C. elegans homolog of the Cockayne syndrome complementation group A gene

Cockayne syndrome (CS) is a debilitating and complex disorder that results from inherited mutations in the CS complementation genes A and B, CSA and CSB. The links between the molecular functions of the CS genes and the complex pathophysiology of CS are as of yet poorly understood and are the subject of intense debate. While mouse models reflect the complexity of CS, studies on simpler genetic models might shed new light on the consequences of CS mutations. Here we describe a functional homolog of the human CSA gene in Caenorhabditis elegans. Similar to its human counterpart, mutations in the nematode csa-1 gene lead to developmental growth defects as a consequence of DNA lesions.

circBase: a database for circular RNAs

Recently, several laboratories have reported thousands of circular RNAs (circRNAs) in animals. Numerous circRNAs are highly stable and have specific spatiotemporal expression patterns. Even though a function for circRNAs is unknown, these features make circRNAs an interesting class of RNAs as possible biomarkers and for further research. We developed a database and website, "circBase," where merged and unified data sets of circRNAs and the evidence supporting their expression can be accessed, downloaded, and browsed within the genomic context. circBase also provides scripts to identify known and novel circRNAs in sequencing data. The database is freely accessible through the web server at http://www.circbase.org/.

Inductive matrix completion for predicting gene-disease associations

MOTIVATION

Most existing methods for predicting causal disease genes rely on specific type of evidence, and are therefore limited in terms of applicability. More often than not, the type of evidence available for diseases varies-for example, we may know linked genes, keywords associated with the disease obtained by mining text, or co-occurrence of disease symptoms in patients. Similarly, the type of evidence available for genes varies-for example, specific microarray probes convey information only for certain sets of genes. In this article, we apply a novel matrix-completion method called Inductive Matrix Completion to the problem of predicting gene-disease associations; it combines multiple types of evidence (features) for diseases and genes to learn latent factors that explain the observed gene-disease associations. We construct features from different biological sources such as microarray expression data and disease-related textual data. A crucial advantage of the method is that it is inductive; it can be applied to diseases not seen at training time, unlike traditional matrix-completion approaches and network-based inference methods that are transductive.

RESULTS

Comparison with state-of-the-art methods on diseases from the Online Mendelian Inheritance in Man (OMIM) database shows that the proposed approach is substantially better-it has close to one-in-four chance of recovering a true association in the top 100 predictions, compared to the recently proposed Catapult method (second best) that has <15% chance. We demonstrate that the inductive method is particularly effective for a query disease with no previously known gene associations, and for predicting novel genes, i.e. genes that are previously not linked to diseases. Thus the method is capable of predicting novel genes even for well-characterized diseases. We also validate the novelty of predictions by evaluating the method on recently reported OMIM associations and on associations recently reported in the literature.

AVAILABILITY

Source code and datasets can be downloaded from http://bigdata.ices.utexas.edu/project/gene-disease.

FixingTIM: interactive exploration of sequence and structural data to identify functional mutations in protein families

BACKGROUND

Knowledge of the 3D structure and functionality of proteins can lead to insight into the associated cellular processes, speed up the creation of pharmaceutical products, and develop drugs that are more effective in combating disease.

METHODS

We present the design and implementation of a visual mining and analysis tool to help identify protein mutations across a family of structural models and to help discover the effect of these mutations on protein function. We integrate 3D structure and sequence information in a common visual interface; multiple linked views and a computational backbone allow comparison at the molecular and atomic levels, while a novel trend-image visual abstraction allows for the sorting and mining of large collections of sequences and of their residues.

RESULTS

We evaluate our approach on the triosephosphate isomerase (TIM) family structural models and sequence data and show that our tool provides an effective, scalable way to navigate a family of proteins, as well as a means to inspect the structure and sequence of individual proteins.

CONCLUSIONS

The TIM application shows that our tool can assist in the navigation of families of proteins, as well as in the exploration of individual protein structures. In conjunction with domain expert knowledge, this interactive tool can help provide biophysical insight into why specific mutations affect function and potentially suggest additional modifications to the protein that could be used to rescue functionality.

BLMP-1/Blimp-1 regulates the spatiotemporal cell migration pattern in C. elegans

Spatiotemporal regulation of cell migration is crucial for animal development and organogenesis. Compared to spatial signals, little is known about temporal signals and the mechanisms integrating the two. In the Caenorhabditis elegans hermaphrodite, the stereotyped migration pattern of two somatic distal tip cells (DTCs) is responsible for shaping the gonad. Guidance receptor UNC-5 is necessary for the dorsalward migration of DTCs. We found that BLMP-1, similar to the mammalian zinc finger transcription repressor Blimp-1/PRDI-BF1, prevents precocious dorsalward turning by inhibiting precocious unc-5 transcription and is only expressed in DTCs before they make the dorsalward turn. Constitutive expression of blmp-1 when BLMP-1 would normally disappear delays unc-5 transcription and causes turn retardation, demonstrating the functional significance of blmp-1 down-regulation. Correct timing of BLMP-1 down-regulation is redundantly regulated by heterochronic genes daf-12, lin-29, and dre-1, which regulate the temporal fates of various tissues. DAF-12, a steroid hormone receptor, and LIN-29, a zinc finger transcription factor, repress blmp-1 transcription, while DRE-1, the F-Box protein of an SCF ubiquitin ligase complex, binds to BLMP-1 and promotes its degradation. We have therefore identified a gene circuit that integrates the temporal and spatial signals and coordinates with overall development of the organism to direct cell migration during organogenesis. The tumor suppressor gene product FBXO11 (human DRE-1 ortholog) also binds to PRDI-BF1 in human cell cultures. Our data suggest evolutionary conservation of these interactions and underscore the importance of DRE-1/FBXO11-mediated BLMP-1/PRDI-BF1 degradation in cellular state transitions during metazoan development.

The migratory path of DTCs determines the shape of the C. elegans gonad. How the spatiotemporal migration pattern is regulated is not clear. We identified a conserved transcription factor BLMP-1 as a central component of a gene regulatory circuit required for the spatiotemporal control of DTC migration. BLMP-1 levels regulate the timing of the DTC dorsal turn, as high levels delay the turn and low levels result in an early turn. We identify and characterize upstream regulators that control BLMP-1 levels. These regulators function in two ways, i.e. by destabilization of BLMP-1 through ubiquitin-mediated proteolysis and by transcriptional repression of the blmp-1 gene to down-regulate BLMP-1. Interestingly, blmp-1 also negatively controls these regulators. Our data suggest that a dietary signal input acts together with a double-negative feedback loop to switch DTCs from the “blmp-1-on” to the “blmp-1-off” state, promoting their dorsal turn. Furthermore, we show that some protein interactions in the circuit are conserved in C. elegans and humans. Our work defines a novel function of the conserved blmp-1 gene in the temporal control of cell migration, and establishes a gene regulatory circuit that integrates the temporal and spatial inputs to direct cell migration during organogenesis.

High-throughput capturing and characterization of mutations in essential genes of Caenorhabditis elegans

BACKGROUND

Essential genes are critical for the development of all organisms and are associated with many human diseases. These genes have been a difficult category to study prior to the availability of balanced lethal strains. Despite the power of targeted mutagenesis, there are limitations in identifying mutations in essential genes. In this paper, we describe the identification of coding regions for essential genes mutated using forward genetic screens in Caenorhabditis elegans. The lethal mutations described here were isolated and maintained by a wild-type allele on a rescuing duplication.

RESULTS

We applied whole genome sequencing to identify the causative molecular lesion resulting in lethality in existing C. elegans mutant strains. These strains are balanced and can be easily maintained for subsequent characterization. Our method can be effectively used to analyze mutations in a large number of essential genes. We describe here the identification of 64 essential genes in a region of chromosome I covered by the duplication sDp2. Of these, 42 are nonsense mutations, six are splice signal mutations, one deletion, and 15 are non-synonymous mutations. Many of the essential genes in this region function in cell cycle, transcriptional regulation, and RNA processing.

CONCLUSIONS

The essential genes identified here are represented by mutant strains, many of which have more than one mutant allele. The genetic resource can be utilized to further our understanding of essential gene function and will be applicable to the study of C. elegans development, conserved cellular function, and ultimately lead to improved human health.

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.

POMO--Plotting Omics analysis results for Multiple Organisms

Background

Systems biology experiments studying different topics and organisms produce thousands of data values across different types of genomic data. Further, data mining analyses are yielding ranked and heterogeneous results and association networks distributed over the entire genome. The visualization of these results is often difficult and standalone web tools allowing for custom inputs and dynamic filtering are limited.

Results

We have developed POMO (http://pomo.cs.tut.fi), an interactive web-based application to visually explore omics data analysis results and associations in circular, network and grid views. The circular graph represents the chromosome lengths as perimeter segments, as a reference outer ring, such as cytoband for human. The inner arcs between nodes represent the uploaded network. Further, multiple annotation rings, for example depiction of gene copy number changes, can be uploaded as text files and represented as bar, histogram or heatmap rings. POMO has built-in references for human, mouse, nematode, fly, yeast, zebrafish, rice, tomato, Arabidopsis, and Escherichia coli. In addition, POMO provides custom options that allow integrated plotting of unsupported strains or closely related species associations, such as human and mouse orthologs or two yeast wild types, studied together within a single analysis. The web application also supports interactive label and weight filtering. Every iterative filtered result in POMO can be exported as image file and text file for sharing or direct future input.

Conclusions

The POMO web application is a unique tool for omics data analysis, which can be used to visualize and filter the genome-wide networks in the context of chromosomal locations as well as multiple network layouts. With the several illustration and filtering options the tool supports the analysis and visualization of any heterogeneous omics data analysis association results for many organisms. POMO is freely available and does not require any installation or registration.

Knockdown of the C. elegans kinome identifies kinases required for normal protein homeostasis, mitochondrial network structure, and sarcomere structure in muscle

BACKGROUND

Kinases are important signalling molecules for modulating cellular processes and major targets of drug discovery programs. However, functional information for roughly half the human kinome is lacking. We conducted three kinome wide, >90%, RNAi screens and epistasis testing of some identified kinases against known intramuscular signalling systems to increase the functional annotation of the C. elegans kinome and expand our understanding of kinome influence upon muscle protein degradation.

RESULTS

96 kinases were identified as required for normal protein homeostasis, 74 for normal mitochondrial networks and 50 for normal sarcomere structure. Knockdown of kinases required only for normal protein homeostasis and/or mitochondrial structure was significantly less likely to produce a developmental or behavioural phenotype than knockdown of kinases required for normal sarcomere structure and/or other sub-cellular processes. Lastly, assessment of kinases for which knockdown produced muscle protein degradation against the known regulatory pathways in C. elegans muscle revealed that close to half of kinase knockdowns activated autophagy in a MAPK dependent fashion.

CONCLUSIONS

Roughly 40% of kinases studied, 159 of 397, are important in establishing or maintaining muscle cell health, with most required for both. For kinases where decreased expression triggers protein degradation, autophagy is most commonly activated. These results increase the annotation of the C. elegans kinome to roughly 75% and enable future kinome research. As 33% of kinases identified have orthologues expressed in human muscle, our results also enable testing of whether identified kinases function similarly in maintaining human muscle homeostasis.

Prediction and validation of gene-disease associations using methods inspired by social network analyses

Correctly identifying associations of genes with diseases has long been a goal in biology. With the emergence of large-scale gene-phenotype association datasets in biology, we can leverage statistical and machine learning methods to help us achieve this goal. In this paper, we present two methods for predicting gene-disease associations based on functional gene associations and gene-phenotype associations in model organisms. The first method, the Katz measure, is motivated from its success in social network link prediction, and is very closely related to some of the recent methods proposed for gene-disease association inference. The second method, called Catapult (Combining dATa Across species using Positive-Unlabeled Learning Techniques), is a supervised machine learning method that uses a biased support vector machine where the features are derived from walks in a heterogeneous gene-trait network. We study the performance of the proposed methods and related state-of-the-art methods using two different evaluation strategies, on two distinct data sets, namely OMIM phenotypes and drug-target interactions. Finally, by measuring the performance of the methods using two different evaluation strategies, we show that even though both methods perform very well, the Katz measure is better at identifying associations between traits and poorly studied genes, whereas Catapult is better suited to correctly identifying gene-trait associations overall.

The authors want to thank Jon Laurent and Kris McGary for some of the data used, and Li and Patra for making their code available. Most of Ambuj Tewari's contribution to this work happened while he was a postdoctoral fellow at the University of Texas at Austin.

Using Multiple Phenotype Assays and Epistasis Testing to Enhance the Reliability of RNAi Screening and Identify Regulators of Muscle Protein Degradation

RNAi is a convenient, widely used tool for screening for genes of interest. We have recently used this technology to screen roughly 750 candidate genes, in C. elegans, for potential roles in regulating muscle protein degradation in vivo. To maximize confidence and assess reproducibility, we have only used previously validated RNAi constructs and have included time courses and replicates. To maximize mechanistic understanding, we have examined multiple sub-cellular phenotypes in multiple compartments in muscle. We have also tested knockdowns of putative regulators of degradation in the context of mutations or drugs that were previously shown to inhibit protein degradation by diverse mechanisms. Here we discuss how assaying multiple phenotypes, multiplexing RNAi screens with use of mutations and drugs, and use of bioinformatics can provide more data on rates of potential false positives and negatives as well as more mechanistic insight than simple RNAi screening.

DGEM — A Microarray Gene Expression Database for Primary Human Disease Tissues

Gene expression patterns can reflect gene regulations in human tissues under normal or pathologic conditions. Gene expression profiling data from studies of primary human disease samples are particularly valuable since these studies often span many years in order to collect patient clinical information and achieve a large sample size. Disease-to-Gene Expression Mapper (DGEM) provides a beneficial community resource to access and analyze these data; it currently includes Affymetrix oligonucleotide array datasets for more than 40 human diseases and 1400 samples. The data are normalized to the same scale and stored in a relational database. A statistical-analysis pipeline was implemented to identify genes abnormally expressed in disease tissues or genes whose expressions are associated with clinical parameters such as cancer patient survival. Data-mining results can be queried through a web-based interface at http://dgem.dhcp.iupui.edu/. The query tool enables dynamic generation of graphs and tables that are further linked to major gene and pathway resources that connect the data to relevant biology, including Entrez Gene and Kyoto Encyclopedia of Genes and Genomes (KEGG). In summary, DGEM provides scientists and physicians a valuable tool to study disease mechanisms, to discover potential disease biomarkers for diagnosis and prognosis, and to identify novel gene targets for drug discovery. The source code is freely available for non-profit use, on request to the authors.

A novel mutation in β integrin reveals an integrin-mediated interaction between the extracellular matrix and cki-1/p27KIP1

The cell-extracellular matrix (ECM) interaction plays an essential role in maintaining tissue shapes and regulates cell behaviors such as cell adhesion, differentiation and proliferation. The mechanism by which the ECM influences the cell cycle in vivo is poorly understood. Here we demonstrate that the β integrin PAT-3 regulates the localization and expression of CKI-1, a C. elegans homologue of the cyclin dependent kinase inhibitor p27(KIP1). In nematodes expressing wild type PAT-3, CKI-1::GFP localizes primarily to nucleoli in hypodermal cells, whereas in animals expressing mutant pat-3 with a defective splice junction, CKI-1::GFP appears clumped and disorganized in nucleoplasm. RNAi analysis links cell adhesion genes to the regulation of CKI-1. RNAi of unc-52/perlecan, ina-1/α integrin, pat-4/ILK, and unc-97/PINCH resulted in abnormal CKI-1::GFP localization. Additional RNAi experiments revealed that the SCF E3 ubiquitin-ligase complex genes, skpt-1/SKP2, cul-1/CUL1 and lin-23/F-box, are required for the proper localization and expression of CKI-1, suggesting that integrin signaling and SCF E3 ligase work together to regulate the cellular distribution of CKI-1. These data also suggest that integrin plays a major role in maintaining proper CKI-1/p27(KIP1) levels in the cell. Perturbed integrin signaling may lead to the inhibition of SCF ligase activity, mislocalization and elevation of CKI-1/p27(KIP1). These results suggest that adhesion signaling is crucial for cell cycle regulation in vivo.

Large-scale functional RNAi screen in C. elegans identifies genes that regulate the dysfunction of mutant polyglutamine neurons

BACKGROUND

A central goal in Huntington's disease (HD) research is to identify and prioritize candidate targets for neuroprotective intervention, which requires genome-scale information on the modifiers of early-stage neuron injury in HD.

RESULTS

Here, we performed a large-scale RNA interference screen in C. elegans strains that express N-terminal huntingtin (htt) in touch receptor neurons. These neurons control the response to light touch. Their function is strongly impaired by expanded polyglutamines (128Q) as shown by the nearly complete loss of touch response in adult animals, providing an in vivo model in which to manipulate the early phases of expanded-polyQ neurotoxicity. In total, 6034 genes were examined, revealing 662 gene inactivations that either reduce or aggravate defective touch response in 128Q animals. Several genes were previously implicated in HD or neurodegenerative disease, suggesting that this screen has effectively identified candidate targets for HD. Network-based analysis emphasized a subset of high-confidence modifier genes in pathways of interest in HD including metabolic, neurodevelopmental and pro-survival pathways. Finally, 49 modifiers of 128Q-neuron dysfunction that are dysregulated in the striatum of either R/2 or CHL2 HD mice, or both, were identified.

CONCLUSIONS

Collectively, these results highlight the relevance to HD pathogenesis, providing novel information on the potential therapeutic targets for neuroprotection in HD.

Silencing of FAD synthase gene in Caenorhabditis elegans upsets protein homeostasis and impacts on complex behavioral patterns

BACKGROUND

FAD synthase is a ubiquitous enzyme that catalyses the last step of FAD biosynthesis, allowing for the biogenesis of several flavoproteins. In humans different isoforms are generated by alternative splicing, isoform 1 being localized in mitochondria. Homology searching in Caenorabditis elegans leads to the identification of two human FAD synthase homologues, coded by the single copy gene R53.1.

METHODS

The C. elegans R53.1 gene was silenced by feeding. The expression level of transcripts was established by semi-quantitative RT-PCR. Overall protein composition was evaluated by two-dimensional electrophoresis. Enzymatic activities were measured by spectrophotometry and oxygen consumption by polarography on isolated mitochondria.

RESULTS

From R53.1 two transcripts are generated by trans-splicing. Reducing by 50% the transcription efficiency of R53.1 by RNAi results in a 50% reduction in total flavin with decrease in ATP content and increase in ROS level. Significant phenotypical changes are noticed in knock-down nematodes. Among them, a significant impairment in locomotion behaviour possibly due to altered cholinergic transmission. At biochemical level, impairment of flavoenzyme activities and of some KCN-insensitive oxygen-consuming enzymes is detected. At proteomic level, at least 15 abundant proteins are affected by R53.1 gene silencing, among which superoxide dismutases.

CONCLUSION AND GENERAL SIGNIFICANCE

For the first time we addressed the existence of different isoforms of FAD-metabolizing enzymes in nematodes. A correlation between FAD synthase silencing and flavoenzyme derangement, energy shortage and redox balance impairment is apparent. In this aspect R53.1-interfered nematodes could provide an animal model system for studying human pathologies with alteration in flavin homeostasis/flavoenzyme biogenesis.

Soybean Knowledge Base (SoyKB): a web resource for soybean translational genomics

BACKGROUND

Soybean Knowledge Base (SoyKB) is a comprehensive all-inclusive web resource for soybean translational genomics. SoyKB is designed to handle the management and integration of soybean genomics, transcriptomics, proteomics and metabolomics data along with annotation of gene function and biological pathway. It contains information on four entities, namely genes, microRNAs, metabolites and single nucleotide polymorphisms (SNPs).

METHODS

SoyKB has many useful tools such as Affymetrix probe ID search, gene family search, multiple gene/metabolite search supporting co-expression analysis, and protein 3D structure viewer as well as download and upload capacity for experimental data and annotations. It has four tiers of registration, which control different levels of access to public and private data. It allows users of certain levels to share their expertise by adding comments to the data. It has a user-friendly web interface together with genome browser and pathway viewer, which display data in an intuitive manner to the soybean researchers, producers and consumers.

CONCLUSIONS

SoyKB addresses the increasing need of the soybean research community to have a one-stop-shop functional and translational omics web resource for information retrieval and analysis in a user-friendly way. SoyKB can be publicly accessed at http://soykb.org/.

Calpains mediate integrin attachment complex maintenance of adult muscle in Caenorhabditis elegans

Two components of integrin containing attachment complexes, UNC-97/PINCH and UNC-112/MIG-2/Kindlin-2, were recently identified as negative regulators of muscle protein degradation and as having decreased mRNA levels in response to spaceflight. Integrin complexes transmit force between the inside and outside of muscle cells and signal changes in muscle size in response to force and, perhaps, disuse. We therefore investigated the effects of acute decreases in expression of the genes encoding these multi-protein complexes. We find that in fully developed adult Caenorhabditis elegans muscle, RNAi against genes encoding core, and peripheral, members of these complexes induces protein degradation, myofibrillar and mitochondrial dystrophies, and a movement defect. Genetic disruption of Z-line– or M-line–specific complex members is sufficient to induce these defects. We confirmed that defects occur in temperature-sensitive mutants for two of the genes: unc-52, which encodes the extra-cellular ligand Perlecan, and unc-112, which encodes the intracellular component Kindlin-2. These results demonstrate that integrin containing attachment complexes, as a whole, are required for proper maintenance of adult muscle. These defects, and collapse of arrayed attachment complexes into ball like structures, are blocked when DIM-1 levels are reduced. Degradation is also blocked by RNAi or drugs targeting calpains, implying that disruption of integrin containing complexes results in calpain activation. In wild-type animals, either during development or in adults, RNAi against calpain genes results in integrin muscle attachment disruptions and consequent sub-cellular defects. These results demonstrate that calpains are required for proper assembly and maintenance of integrin attachment complexes. Taken together our data provide in vivo evidence that a calpain-based molecular repair mechanism exists for dealing with attachment complex disruption in adult muscle. Since C. elegans lacks satellite cells, this mechanism is intrinsic to the muscles and raises the question if such a mechanism also exists in higher metazoans.

Muscle is a dynamic tissue that grows in response to use and nutrition and shrinks in response to lack of use, poor nutrition, or disease. Loss of muscle mass is an important public health problem, but we understand little of the genes that regulate muscle shrinkage. We have found that, in adult worm muscle, attachment to the basement membrane is continuously required to prevent catastrophic sub-cellular defects that result in impaired ability of muscle to function. We have also identified a group of proteases that are activated when the attachment fails to be properly maintained. Conversely, when these proteases are lacking in adult muscle, the muscles fail to maintain attachment to the basement membrane. Thus, we have discovered a group of proteases that appear to act to maintain attachment to the basement membrane and therefore to maintain muscle itself. Because these worms lack satellite cells, this maintenance system is intrinsic to muscle, thus raising the question whether a similar or identical system also works in humans.

Multiple phenotypes resulting from a mutagenesis screen for pharynx muscle mutations in Caenorhabditis elegans

We describe a novel screen to isolate pharyngeal cell morphology mutants in Caenorhabditis elegans using myo-2::GFP to rapidly identify abnormally shaped pharynxes in EMS (Ethyl Methanesulfonate) mutagenized worms. We observed over 83 C. elegans lines with distinctive pharyngeal phenotypes in worms surviving to the L1 larval stage, with phenotypes ranging from short pharynx, unattached pharynx, missing cells, asymmetric morphology, and non-adherent pharynx cells. Thirteen of these mutations have been chromosomally mapped using Single Nucleotide Polymorphisms (SNPs) and deficiency strain complementation. Our studies have focused on genetically mapping and functionally testing two phenotypes, the short pharynx and the loss of muscle cohesion phenotypes. We have also identified new alleles of sma-1, and our screen suggests many genes directing pharynx assembly and structure may be either pharynx specific or less critical in other tissues.

Identification and functional clustering of genes regulating muscle protein degradation from amongst the known C. elegans muscle mutants

Loss of muscle mass via protein degradation is an important clinical problem but we know little of how muscle protein degradation is regulated genetically. To gain insight our labs developed C. elegans into a model for understanding the regulation of muscle protein degradation. Past studies uncovered novel functional roles for genes affecting muscle and/or involved in signalling in other cells or tissues. Here we examine most of the genes previously identified as the sites of mutations affecting muscle for novel roles in regulating degradation. We evaluate genomic (RNAi knockdown) approaches and combine them with our established genetic (mutant) and pharmacologic (drugs) approaches to examine these 159 genes. We find that RNAi usually recapitulates both organismal and sub-cellular mutant phenotypes but RNAi, unlike mutants, can frequently be used acutely to study gene function solely in differentiated muscle. In the majority of cases where RNAi does not produce organismal level phenotypes, sub-cellular defects can be detected; disrupted proteostasis is most commonly observed. We identify 48 genes in which mutation or RNAi knockdown causes excessive protein degradation; myofibrillar and/or mitochondrial morphologies are also disrupted in 19 of these 48 cases. These 48 genes appear to act via at least three sub-networks to control bulk degradation of protein in muscle cytosol. Attachment to the extracellular matrix regulates degradation via unidentified proteases and affects myofibrillar and mitochondrial morphology. Growth factor imbalance and calcium overload promote lysosome based degradation whereas calcium deficit promotes proteasome based degradation, in both cases myofibrillar and mitochondrial morphologies are largely unaffected. Our results provide a framework for effectively using RNAi to identify and functionally cluster novel regulators of degradation. This clustering allows prioritization of candidate genes/pathways for future mechanistic studies.

Guanylyl cyclase structure, function and regulation

Nitric oxide, bicarbonate, natriuretic peptides (ANP, BNP and CNP), guanylins, uroguanylins and guanylyl cyclase activating proteins (GCAPs) activate a family of enzymes variously called guanyl, guanylyl or guanylate cyclases that catalyze the conversion of guanosine triphosphate to cyclic guanosine monophosphate (cGMP) and pyrophosphate. Intracellular cyclic GMP is a second messenger that modulates: platelet aggregation, neurotransmission, sexual arousal, gut peristalsis, blood pressure, long bone growth, intestinal fluid secretion, lipolysis, phototransduction, cardiac hypertrophy and oocyte maturation. This review briefly discusses the discovery of cGMP and guanylyl cyclases, then nitric oxide, nitric oxide synthase and soluble guanylyl cyclase are described in slightly greater detail. Finally, the structure, function, and regulation of the individual mammalian single membrane-spanning guanylyl cyclases GC-A, GC-B, GC-C, GC-D, GC-E, GC-F and GC-G are described in greatest detail as determined by biochemical, cell biological and gene-deletion studies.

Identifying insertion mutations by whole-genome sequencing

Insertion mutagenesis via mobile genetic element is a common technique for the analysis of gene function in model organisms. Next-generation sequencing offers an attractive approach for localizing the site of insertion, but alignment-based mapping of mobile genetic elements is challenging. A computational method for identifying insertion sites is reported herein. The technique was validated by mapping transposons in both bacterial and nematode species. The approach should be extensible to other systems that employ mobile genetic elements to generate mutations.

Identifying novel genes in C. elegans using SAGE tags

Background

Despite extensive efforts devoted to predicting protein-coding genes in genome sequences, many bona fide genes have not been found and many existing gene models are not accurate in all sequenced eukaryote genomes. This situation is partly explained by the fact that gene prediction programs have been developed based on our incomplete understanding of gene feature information such as splicing and promoter characteristics. Additionally, full-length cDNAs of many genes and their isoforms are hard to obtain due to their low level or rare expression. In order to obtain full-length sequences of all protein-coding genes, alternative approaches are required.

Results

In this project, we have developed a method of reconstructing full-length cDNA sequences based on short expressed sequence tags which is called sequence tag-based amplification of cDNA ends (STACE). Expressed tags are used as anchors for retrieving full-length transcripts in two rounds of PCR amplification. We have demonstrated the application of STACE in reconstructing full-length cDNA sequences using expressed tags mined in an array of serial analysis of gene expression (SAGE) of C. elegans cDNA libraries. We have successfully applied STACE to recover sequence information for 12 genes, for two of which we found isoforms. STACE was used to successfully recover full-length cDNA sequences for seven of these genes.

Conclusions

The STACE method can be used to effectively reconstruct full-length cDNA sequences of genes that are under-represented in cDNA sequencing projects and have been missed by existing gene prediction methods, but their existence has been suggested by short sequence tags such as SAGE tags.

Modeling human diseases in Caenorhabditis elegans

Genes linked to human diseases often function in evolutionarily conserved pathways, which can be readily dissected in simple model organisms. Because of its short lifespan and well-known biology, coupled with a completely sequenced genome that shares extensive homology with that of mammals, Caenorhabditis elegans is one of the most versatile and powerful model organisms. Research in C. elegans has been instrumental for the elucidation of molecular pathways implicated in many human diseases. In this review, we introduce C. elegans as a model organism for biomedical research and we survey recent relevant findings that shed light on the basic molecular determinants of human disease pathophysiology. The nematode holds promise of providing clear leads towards the identification of potential targets for the development of new therapeutic interventions against human diseases.

Predicted sub-populations in a marine shrimp proteome as revealed by combined EST and cDNA data from multiple Penaeus species

Background

Many species of marine shrimp in the Family Penaeidae, viz. Penaeus (Litopenaeus) vannamei, Penaeus monodon, Penaeus (Fenneropenaeus) chinensis, and Penaeus (Marsupenaeus) japonicus, are animals of economic importance in the aquaculture industry. Yet information about their DNA and protein sequences is lacking. In order to predict their collective proteome, we combined over 270,000 available EST and cDNA sequences from the 4 shrimp species with all protein sequences of Drosophila melanogaster and Caenorhabditis elegans. EST data from 4 other crustaceans, the crab Carcinus maenas, the lobster Homarus americanus (Decapoda), the water flea Daphnia pulex, and the brine shrimp Artemia franciscana were also used.

Findings

Similarity searches from EST collections of the 4 shrimp species matched 64% of the protein sequences of the fruit fly, but only 45% of nematode proteins, indicating that the shrimp proteome content is more similar to that of an insect than a nematode. Combined results with 4 additional non-shrimp crustaceans increased matching to 78% of fruit fly and 56% of nematode proteins, suggesting that present shrimp EST collections still lack sequences for many conserved crustacean proteins. Analysis of matching data revealed the presence of 4 EST groups from shrimp, namely sequences for proteins that are both fruit fly-like and nematode-like, fruit fly-like only, nematode-like only, and non-matching. Gene ontology profiles of proteins for the 3 matching EST groups were analyzed. For non-matching ESTs, a small fraction matched protein sequences from other species in the UniProt database, including other crustacean-specific proteins.

Conclusions

Shrimp ESTs indicated that the shrimp proteome is comprised of sub-populations of proteins similar to those common to both insect and nematode models, those present specifically in either model, or neither. Combining small EST collections from related species to compensate for their small size allowed prediction of conserved expressed protein components encoded by their uncharacterized genomes. The organized data should be useful for transferring annotation data from model species into shrimp data and for further studies on shrimp proteins with particular functions or groups.

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.

Elucidating gene regulatory mechanisms for sperm function through the integration of classical and systems approaches in C. elegans

From worms to mammals, successful spermatogenesis depends on a gene expression profile that balances activating and repressive mechanisms. Besides developmental control of specific spermatogenic genes, male fertility requires temporal shifts in global gene expression and dramatic changes in chromatin structure and condensation. Recent studies are beginning to elucidate the molecular processes that both drive these temporal changes in gene expression and underlie fertility. In this review, we provide an overview of relevant C. elegans studies that have laid the groundwork for modern approaches. Next, we highlight recent studies that investigate how gene expression in C. elegans is modulated during spermatogenesis. These studies use large-scale genomic profiling in combination with bioinformatics, genetics, biochemistry, and in vitro methods to target specific stages or processes during sperm formation. Such studies are beginning to elucidate the multiple layers of gene regulation required during spermatogenesis, i.e., transcriptional, post-transcriptional, and epigenetic. Moreover, knowledge of how C. elegans coordinately regulates gene expression during spermatogenesis promises to provide key insights into parallel processes in mammals that are vital for fertility.

Phylogenomics of ligand-gated ion channels predicts monepantel effect

The recently launched veterinary anthelmintic drench for sheep (Novartis Animal Health Inc., Switzerland) containing the nematocide monepantel represents a new class of anthelmintics: the amino-acetonitrile derivatives (AADs), much needed in view of widespread resistance to the classical drugs. Recently, it was shown that the ACR-23 protein in Caenorhabditis elegans and a homologous protein, MPTL-1 in Haemonchus contortus, are potential targets for AAD action. Both proteins belong to the DEG-3 subfamily of acetylcholine receptors, which are thought to be nematode-specific, and different from those targeted by the imidazothiazoles (e.g. levamisole). Here we provide further evidence that Cel-ACR-23 and Hco-MPTL-1-like subunits are involved in the monepantel-sensitive phenotype. We performed comparative genomics of ligand-gated ion channel genes from several nematodes and subsequently assessed their sensitivity to anthelmintics. The nematode species in the Caenorhabditis genus, equipped with ACR-23/MPTL-1-like receptor subunits, are sensitive to monepantel (EC₅₀<1.25 µM), whereas the related nematodes Pristionchus pacificus and Strongyloides ratti, which lack an ACR-23/MPTL-1 homolog, are insensitive (EC₅₀>43 µM). Genome sequence information has long been used to identify putative targets for therapeutic intervention. We show how comparative genomics can be applied to predict drug sensitivity when molecular targets of a compound are known or suspected.

Increased use of anthelmintics has contributed to the emergence of drug-resistant nematodes, causing serious problems for more than one billion sheep worldwide. The last class of compounds indicated for livestock was introduced 28 years ago. Recently, however, Novartis AH developed a new anthelmintic active against drug-resistant nematodes of sheep, the amino-acetonitrile derivative (AAD) monepantel. We have previously indirectly shown that the AADs have a novel mode of action involving acetylcholine receptor subunits: the ACR-23 protein in Caenorhabditis elegans and a homologous protein, MPTL-1 in Haemonchus contortus. To better understand the mode of action of the AADs, we performed comparative genomics of all ligand-gated ion channel genes from a range of organisms, including members from all nematode clades. We confirmed that MPTL-1 belongs to a unique, nematode-specific sub-family of receptor subunits. We also found that some nematode species lack ACR-23/MPTL-1 and predicted them to be monepantel insensitive. We challenged this hypothesis in a panel of drug tests: several species of Caenorhabditis nematodes equipped with ACR-23/MPTL-1-like receptor subunits were found susceptible to monepantel, whereas Pristionchus pacificus, closely related to these worms but lacking an ACR-23/MPTL-1 homolog, was tolerant. The parasitic nematode Strongyloides ratti, which has only a remote homolog of DES-2 and ACR-23/MPTL-1, was also tolerant to monepantel. This confirms our prediction and highlights how comparative genomic data can be used to predict a drug effect.