Defining characteristics and conservation of poorly annotated genes in Caenorhabditis elegans using WormCat 2.0.

Omics tools provide broad datasets for biological discovery. However, the computational tools for identifying important genes or pathways in RNA-seq, proteomics, or GWAS (Genome-Wide Association Study) data depend on Gene Ontogeny annotations and are biased toward well-described pathways. This limits their utility as poorly annotated genes, which could have novel functions, are often passed over. Recently, we developed an annotation and category enrichment tool for Caenorhabditis elegans genomic data, WormCat, which provides an intuitive visualization output. Unlike Gene Ontogeny-based enrichment tools, which exclude genes with no annotation information, WormCat 2.0 retains these genes as a special UNASSIGNED category. Here, we show that the UNASSIGNED gene category enrichment exhibits tissue-specific expression patterns and can include genes with biological functions identified in published datasets. Poorly annotated genes are often considered to be potentially species-specific and thus, of reduced interest to the biomedical community. Instead, we find that around 3% of the UNASSIGNED genes have human orthologs, including some linked to human diseases. These human orthologs themselves have little annotation information. A recently developed method that incorporates lineage relationships (abSENSE) indicates that the failure of BLAST to detect homology explains the apparent lineage specificity for many UNASSIGNED genes. This suggests that a larger subset could be related to human genes. WormCat provides an annotation strategy that allows the association of UNASSIGNED genes with specific phenotypes and known pathways. Building these associations in C. elegans, with its robust genetic tools, provides a path to further functional study and insight into these understudied genes.

WormCat: An Online Tool for Annotation and Visualization of Caenorhabditis elegans Genome-Scale Data.

The emergence of large gene expression datasets has revealed the need for improved tools to identify enriched gene categories and visualize enrichment patterns. While gene ontogeny (GO) provides a valuable tool for gene set enrichment analysis, it has several limitations. First, it is difficult to graph multiple GO analyses for comparison. Second, genes from some model systems are not well represented. For example, ∼30% of Caenorhabditis elegans genes are missing from the analysis in commonly used databases. To allow categorization and visualization of enriched C. elegans gene sets in different types of genome-scale data, we developed WormCat, a web-based tool that uses a near-complete annotation of the C. elegans genome to identify coexpressed gene sets and scaled heat map for enrichment visualization. We tested the performance of WormCat using a variety of published transcriptomic datasets, and show that it reproduces major categories identified by GO. Importantly, we also found previously unidentified categories that are informative for interpreting phenotypes or predicting biological function. For example, we analyzed published RNA-seq data from C. elegans treated with combinations of lifespan-extending drugs, where one combination paradoxically shortened lifespan. Using WormCat, we identified sterol metabolism as a category that was not enriched in the single or double combinations, but emerged in a triple combination along with the lifespan shortening. Thus, WormCat identified a gene set with potential. phenotypic relevance not found with previous GO analysis. In conclusion, WormCat provides a powerful tool for the analysis and visualization of gene set enrichment in different types of C. elegans datasets.

Stress-responsive and metabolic gene regulation are altered in low S-adenosylmethionine.

S-adenosylmethionine (SAM) is a donor which provides the methyl groups for histone or nucleic acid modification and phosphatidylcholine production. SAM is hypothesized to link metabolism and chromatin modification, however, its role in acute gene regulation is poorly understood. We recently found that Caenorhabditis elegans with reduced SAM had deficiencies in H3K4 trimethylation (H3K4me3) at pathogen-response genes, decreasing their expression and limiting pathogen resistance. We hypothesized that SAM may be generally required for stress-responsive transcription. Here, using genetic assays, we show that transcriptional responses to bacterial or xenotoxic stress fail in C. elegans with low SAM, but that expression of heat shock genes are unaffected. We also found that two H3K4 methyltransferases, set-2/SET1 and set-16/MLL, had differential responses to survival during stress. set-2/SET1 is specifically required in bacterial responses, whereas set-16/MLL is universally required. These results define a role for SAM in the acute stress-responsive gene expression. Finally, we find that modification of metabolic gene expression correlates with enhanced survival during stress.