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

Meta-Research: Understudied genes are lost in a leaky pipeline between genome-wide assays and reporting of results

Present-day publications on human genes primarily feature genes that already appeared in many publications prior to completion of the Human Genome Project in 2003. These patterns persist despite the subsequent adoption of high-throughput technologies, which routinely identify novel genes associated with biological processes and disease. Although several hypotheses for bias in the selection of genes as research targets have been proposed, their explanatory powers have not yet been compared. Our analysis suggests that understudied genes are systematically abandoned in favor of better-studied genes between the completion of -omics experiments and the reporting of results. Understudied genes remain abandoned by studies that cite these -omics experiments. Conversely, we find that publications on understudied genes may even accrue a greater number of citations. Among 45 biological and experimental factors previously proposed to affect which genes are being studied, we find that 33 are significantly associated with the choice of hit genes presented in titles and abstracts of -omics studies. To promote the investigation of understudied genes, we condense our insights into a tool, find my understudied genes (FMUG), that allows scientists to engage with potential bias during the selection of hits. We demonstrate the utility of FMUG through the identification of genes that remain understudied in vertebrate aging. FMUG is developed in Flutter and is available for download at fmug.amaral.northwestern.edu as a MacOS/Windows app.

A molecular atlas of adult C. elegans motor neurons reveals ancient diversity delineated by conserved transcription factor codes

Motor neurons (MNs) constitute an ancient cell type targeted by multiple adult-onset diseases. It is therefore important to define the molecular makeup of adult MNs in animal models and extract organizing principles. Here, we generate a comprehensive molecular atlas of adult Caenorhabditis elegans MNs and a searchable database. Single-cell RNA sequencing of 13,200 cells reveals that ventral nerve cord MNs cluster into 29 molecularly distinct subclasses. Extending C. elegans Neuronal Gene Expression Map and Network (CeNGEN) findings, all MN subclasses are delineated by distinct expression codes of either neuropeptide or transcription factor gene families. Strikingly, combinatorial codes of homeodomain transcription factor genes succinctly delineate adult MN diversity in both C. elegans and mice. Further, molecularly defined MN subclasses in C. elegans display distinct patterns of connectivity. Hence, our study couples the connectivity map of the C. elegans motor circuit with a molecular atlas of its constituent MNs and uncovers organizing principles and conserved molecular codes of adult MN diversity.

Combinatorial transcriptomic and genetic dissection of insulin/IGF-1 signaling-regulated longevity in Caenorhabditis elegans

Classical genetic analysis is invaluable for understanding the genetic interactions underlying specific phenotypes, but requires laborious and subjective experiments to characterize polygenic and quantitative traits. Contrarily, transcriptomic analysis enables the simultaneous and objective identification of multiple genes whose expression changes are associated with specific phenotypes. Here, we conducted transcriptomic analysis of genes crucial for longevity using datasets with daf-2/insulin/IGF-1 receptor mutant Caenorhabditis elegans. Our analysis unraveled multiple epistatic relationships at the transcriptomic level, in addition to verifying genetically established interactions. Our combinatorial analysis also revealed transcriptomic changes associated with longevity conferred by daf-2 mutations. In particular, we demonstrated that the extent of lifespan changes caused by various mutant alleles of the longevity transcription factor daf-16/FOXO matched their effects on transcriptomic changes in daf-2 mutants. We identified specific aging-regulating signaling pathways and subsets of structural and functional RNA elements altered by different genes in daf-2 mutants. Lastly, we elucidated the functional cooperation between several longevity regulators, based on the combination of transcriptomic and molecular genetic analysis. These data suggest that different biological processes coordinately exert their effects on longevity in biological networks. Together our work demonstrates the utility of transcriptomic dissection analysis for identifying important genetic interactions for physiological processes, including aging and longevity.

Modeling neurodevelopmental disorder-associated human AGO1 mutations in Caenorhabditis elegans Argonaute alg-1

MicroRNAs (miRNA) associate with Argonaute (AGO) proteins and repress gene expression by base pairing to sequences in the 3' untranslated regions of target genes. De novo coding variants in the human AGO genes AGO1 and AGO2 cause neurodevelopmental disorders (NDD) with intellectual disability, referred to as Argonaute syndromes. Most of the altered amino acids are conserved between the miRNA-associated AGO in Homo sapiens and Caenorhabditis elegans, suggesting that the human mutations could disrupt conserved functions in miRNA biogenesis or activity. We genetically modeled four human AGO1 mutations in C. elegans by introducing identical mutations into the C. elegans AGO1 homologous gene, alg-1. These alg-1 NDD mutations cause phenotypes in C. elegans indicative of disrupted miRNA processing, miRISC (miRNA silencing complex) formation, and/or target repression. We show that the alg-1 NDD mutations are antimorphic, causing developmental and molecular phenotypes stronger than those of alg-1 null mutants, likely by sequestrating functional miRISC components into non-functional complexes. The alg-1 NDD mutations cause allele-specific disruptions in mature miRNA profiles, accompanied by perturbation of downstream gene expression, including altered translational efficiency and/or messenger RNA abundance. The perturbed genes include those with human orthologs whose dysfunction is associated with NDD. These cross-clade genetic studies illuminate fundamental AGO functions and provide insights into the conservation of miRNA-mediated post-transcriptional regulatory mechanisms.

Meta-Research: understudied genes are lost in a leaky pipeline between genome-wide assays and reporting of results

Present-day publications on human genes primarily feature genes that already appeared in many publications prior to completion of the Human Genome Project in 2003. These patterns persist despite the subsequent adoption of high-throughput technologies, which routinely identify novel genes associated with biological processes and disease. Although several hypotheses for bias in the selection of genes as research targets have been proposed, their explanatory powers have not yet been compared. Our analysis suggests that understudied genes are systematically abandoned in favor of better-studied genes between the completion of -omics experiments and the reporting of results. Understudied genes remain abandoned by studies that cite these -omics experiments. Conversely, we find that publications on understudied genes may even accrue a greater number of citations. Among 45 biological and experimental factors previously proposed to affect which genes are being studied, we find that 33 are significantly associated with the choice of hit genes presented in titles and abstracts of - omics studies. To promote the investigation of understudied genes we condense our insights into a tool, find my understudied genes (FMUG), that allows scientists to engage with potential bias during the selection of hits. We demonstrate the utility of FMUG through the identification of genes that remain understudied in vertebrate aging. FMUG is developed in Flutter and is available for download at fmug.amaral.northwestern.edu as a MacOS/Windows app.

In vivo neuroprotective capacity of a Dunaliella salina extract - comprehensive transcriptomics and metabolomics study

In this study, an exhaustive chemical characterization of a Dunaliella salina (DS) microalga extract obtained using supercritical fluids has been performed, and its neuroprotective capacity has been evaluated in vivo using an Alzheimer's disease (AD) transgenic model of Caenorhabditis elegans (strain CL4176). More than 350 compounds were annotated in the studied DS extract, with triacylglycerols, free fatty acids (FAs), carotenoids, apocarotenoids and glycerol being the most abundant. DS extract significantly protects C. elegans in a dose-dependent manner against Aβ-peptide paralysis toxicity, after 32 h, 53% of treated worms at 50 µg/mL were not paralyzed. This concentration was selected to further evaluate the transcriptomics and metabolomics changes after 26 h by using advanced analytical methodologies. The RNA-Seq data showed an alteration of 150 genes, mainly related to the stress and detoxification responses, and the retinol and lipid metabolism. The comprehensive metabolomics and lipidomics analyses allowed the identification of 793 intracellular metabolites, of which 69 were significantly altered compared to non-treated control animals. Among them, different unsaturated FAs, lysophosphatidylethanolamines, nucleosides, dipeptides and modified amino acids that have been previously reported as beneficial during AD progression, were assigned. These compounds could explain the neuroprotective capacity observed, thus, providing with new evidences of the protection mechanisms of this promising extract.

Immunity-linked genes are stimulated by a membrane stress pathway linked to Golgi function and the ARF-1 GTPase

Infection response and other immunity-linked genes (ILGs) were first named in Caenorhabditis elegans-based expression after pathogen challenge, but many are also up-regulated when lipid metabolism is perturbed. Why pathogen attack and metabolic changes both increase ILGs is unclear. We find that ILGs are activated when phosphatidylcholine (PC) levels change in membranes of secretory organelles in C. elegans. RNAi targeting of the ADP-ribosylation factor arf-1, which disrupts the Golgi and secretory function, also activates ILGs. Low PC limits ARF-1 function, suggesting a mechanism for ILG activation via lipid metabolism, as part of a membrane stress response acting outside the ER. RNAi of selected ILGs uncovered defects in the secretion of two GFP reporters and the accumulation of a pathogen-responsive complement C1r/C1s, Uegf, Bmp1 (CUB) domain fusion protein. Our data argue that up-regulation of some ILGs is a coordinated response to changes in trafficking and may act to counteract stress on secretory function.

Transcriptional suppression of sphingolipid catabolism controls pathogen resistance in C. elegans

Sphingolipids are required for diverse biological functions and are degraded by specific catabolic enzymes. However, the mechanisms that regulate sphingolipid catabolism are not known. Here we characterize a transcriptional axis that regulates sphingolipid breakdown to control resistance against bacterial infection. From an RNAi screen for transcriptional regulators of pathogen resistance in the nematode C. elegans, we identified the nuclear hormone receptor nhr-66, a ligand-gated transcription factor homologous to human hepatocyte nuclear factor 4. Tandem chromatin immunoprecipitation-sequencing and RNA sequencing experiments revealed that NHR-66 is a transcriptional repressor, which directly targets sphingolipid catabolism genes. Transcriptional de-repression of two sphingolipid catabolic enzymes in nhr-66 loss-of-function mutants drives the breakdown of sphingolipids, which enhances host susceptibility to infection with the bacterial pathogen Pseudomonas aeruginosa. These data define transcriptional control of sphingolipid catabolism in the regulation of cellular sphingolipids, a process that is necessary for pathogen resistance.

The Y-ome Conundrum: Insights into Uncharacterized Genes and Approaches for Functional Annotation

The ever-increasing availability of genome sequencing data has revealed a substantial number of uncharacterized genes without known functions across various organisms. The first comprehensive genome sequencing of E. coli K12 revealed that more than 50% of its open reading frames corresponded to transcripts with no known functions. The group of protein-coding genes without a functional description and/or a recognized pathway, beginning with the letter "Y", is classified as the "y-ome". Several efforts have been made to elucidate the functions of these genes and to recognize their role in biological processes. This review provides a brief update on various strategies employed when studying the y-ome, such as high-throughput experimental approaches, comparative omics, metabolic engineering, gene expression analysis, and data integration techniques. Additionally, we highlight recent advancements in functional annotation methods, including the use of machine learning, network analysis, and functional genomics approaches. Novel approaches are required to produce more precise functional annotations across the genome to reduce the number of genes with unknown functions.

A molecular atlas of adult C. elegans motor neurons reveals ancient diversity delineated by conserved transcription factor codes

Motor neurons (MNs) constitute an ancient cell type targeted by multiple adult-onset diseases. It is therefore important to define the molecular makeup of adult MNs in animal models and extract organizing principles. Here, we generated a comprehensive molecular atlas of adult Caenorhabditis elegans MNs and a searchable database (http://celegans.spinalcordatlas.org). Single-cell RNA-sequencing of 13,200 cells revealed that ventral nerve cord MNs cluster into 29 molecularly distinct subclasses. All subclasses are delineated by unique expression codes of either neuropeptide or transcription factor gene families. Strikingly, we found that combinatorial codes of homeodomain transcription factor genes define adult MN diversity both in C. elegans and mice. Further, molecularly defined MN subclasses in C. elegans display distinct patterns of connectivity. Hence, our study couples the connectivity map of the C. elegans motor circuit with a molecular atlas of its constituent MNs, and uncovers organizing principles and conserved molecular codes of adult MN diversity.

Aminopeptidase MNP-1 triggers intestine protease production by activating daf-16 nuclear location to degrade pore-forming toxins in Caenorhabditis elegans

Pore-forming toxins (PFTs) are effective tools for pathogens infection. By disrupting epithelial barriers and killing immune cells, PFTs promotes the colonization and reproduction of pathogenic microorganisms in their host. In turn, the host triggers defense responses, such as endocytosis, exocytosis, or autophagy. Bacillus thuringiensis (Bt) bacteria produce PFT, known as crystal proteins (Cry) which damage the intestinal cells of insects or nematodes, eventually killing them. In insects, aminopeptidase N (APN) has been shown to act as an important receptor for Cry toxins. Here, using the nematode Caenorhabditis elegans as model, an extensive screening of APN gene family was performed to analyze the potential role of these proteins in the mode of action of Cry5Ba against the nematode. We found that one APN, MNP-1, participate in the toxin defense response, since the mnp-1(ok2434) mutant showed a Cry5Ba hypersensitive phenotype. Gene expression analysis in mnp-1(ok2434) mutant revealed the involvement of two protease genes, F19C6.4 and R03G8.6, that participate in Cry5Ba degradation. Finally, analysis of the transduction pathway involved in F19C6.4 and R03G8.6 expression revealed that upon Cry5Ba exposure, the worms up regulated both protease genes through the activation of the FOXO transcription factor DAF-16, which was translocated into the nucleus. The nuclear location of DAF-16 was found to be dependent on mnp-1 under Cry5Ba treatment. Our work provides evidence of new host responses against PFTs produced by an enteric pathogenic bacterium, resulting in activation of host intestinal proteases that degrade the PFT in the intestine.

Status of genome function annotation in model organisms and crops

Since the entry into genome-enabled biology several decades ago, much progress has been made in determining, describing, and disseminating the functions of genes and their products. Yet, this information is still difficult to access for many scientists and for most genomes. To provide easy access and a graphical summary of the status of genome function annotation for model organisms and bioenergy and food crop species, we created a web application (https://genomeannotation.rheelab.org) to visualize, search, and download genome annotation data for 28 species. The summary graphics and data tables will be updated semi-annually, and snapshots will be archived to provide a historical record of the progress of genome function annotation efforts. Clear and simple visualization of up-to-date genome function annotation status, including the extent of what is unknown, will help address the grand challenge of elucidating the functions of all genes in organisms.

Non-canonical pattern recognition of a pathogen-derived metabolite by a nuclear hormone receptor identifies virulent bacteria in C. elegans

Distinguishing infectious pathogens from harmless microorganisms is essential for animal health. The mechanisms used to identify infectious microbes are not fully understood, particularly in metazoan hosts that eat bacteria as their food source. Here, we characterized a non-canonical pattern-recognition system in Caenorhabditis elegans (C. elegans) that assesses the relative threat of virulent Pseudomonas aeruginosa (P. aeruginosa) to activate innate immunity. We discovered that the innate immune response in C. elegans was triggered by phenazine-1-carboxamide (PCN), a toxic metabolite produced by pathogenic strains of P. aeruginosa. We identified the nuclear hormone receptor NHR-86/HNF4 as the PCN sensor in C. elegans and validated that PCN bound to the ligand-binding domain of NHR-86/HNF4. Activation of NHR-86/HNF4 by PCN directly engaged a transcriptional program in intestinal epithelial cells that protected against P. aeruginosa. Thus, a bacterial metabolite is a pattern of pathogenesis surveilled by nematodes to identify a pathogen in its bacterial diet.

Modeling neurodevelopmental disorder-associated hAGO1 mutations in C. elegans Argonaute ALG-1

MicroRNAs (miRNA) are endogenous non-coding RNAs important for post-transcriptional regulation of gene expression. miRNAs associate with Argonaute proteins to bind to the 3' UTR of target genes and confer target repression. Recently, multiple de novo coding variants in the human Argonaute gene AGO1 ( hAGO1 ) have been reported to cause a neurodevelopmental disorder (NDD) with intellectual disability (ID). Most of the altered amino acids are conserved between the miRNA-associated Argonautes in H. sapiens and C. elegans , suggesting the hAGO1 mutations could disrupt evolutionarily conserved functions in the miRNA pathway. To investigate how the hAGO1 mutations may affect miRNA biogenesis and/or functions, we genetically modeled four of the hAGO1 de novo variants (referred to as NDD mutations) by introducing the identical mutations to the C. elegans hAGO1 homolog, alg-1 . This array of mutations caused distinct effects on C. elegans miRNA functions, miRNA populations, and downstream gene expression, indicative of profound alterations in aspects of miRNA processing and miRISC formation and/or activity. Specifically, we found that the alg-1 NDD mutations cause allele-specific disruptions in mature miRNA profiles both in terms of overall abundances and association with mutant ALG-1. We also observed allele-specific profiles of gene expression with altered translational efficiency and/or mRNA abundance. The sets of perturbed genes include human homologs whose dysfunction is known to cause NDD. We anticipate that these cross-clade genetic studies may advance the understanding of fundamental Argonaute functions and provide insights into the conservation of miRNA-mediated post-transcriptional regulatory mechanisms.

Cell context-dependent CFI-1/ARID3 functions control neuronal terminal differentiation

AT-rich interaction domain 3 (ARID3) transcription factors are expressed in the nervous system, but their mechanisms of action are largely unknown. Here, we provide, in vivo, a genome-wide binding map for CFI-1, the sole C. elegans ARID3 ortholog. We identify 6,396 protein-coding genes as putative direct targets of CFI-1, most of which encode neuronal terminal differentiation markers. In head sensory neurons, CFI-1 directly activates multiple terminal differentiation genes, thereby acting as a terminal selector. In motor neurons, however, CFI-1 acts as a direct repressor, continuously antagonizing three transcriptional activators. By focusing on the glr-4/GRIK4 glutamate receptor locus, we identify proximal CFI-1 binding sites and histone methyltransferase activity as necessary for glr-4 repression. Rescue assays reveal functional redundancy between core and extended DNA-binding ARID domains and a strict requirement for REKLES, the ARID3 oligomerization domain. Altogether, this study uncovers cell-context-dependent mechanisms through which a single ARID3 protein controls the terminal differentiation of distinct neuron types.

S-adenosylmethionine synthases specify distinct H3K4me3 populations and gene expression patterns during heat stress

Methylation is a widely occurring modification that requires the methyl donor S-adenosylmethionine (SAM) and acts in regulation of gene expression and other processes. SAM is synthesized from methionine, which is imported or generated through the 1-carbon cycle (1 CC). Alterations in 1 CC function have clear effects on lifespan and stress responses, but the wide distribution of this modification has made identification of specific mechanistic links difficult. Exploiting a dynamic stress-induced transcription model, we find that two SAM synthases in Caenorhabditis elegans, SAMS-1 and SAMS-4, contribute differently to modification of H3K4me3, gene expression and survival. We find that sams-4 enhances H3K4me3 in heat shocked animals lacking sams-1, however, sams-1 cannot compensate for sams-4, which is required to survive heat stress. This suggests that the regulatory functions of SAM depend on its enzymatic source and that provisioning of SAM may be an important regulatory step linking 1 CC function to phenotypes in aging and stress.

Neuronal temperature perception induces specific defenses that enable C. elegans to cope with the enhanced reactivity of hydrogen peroxide at high temperature

Hydrogen peroxide is the most common reactive chemical that organisms face on the microbial battlefield. The rate with which hydrogen peroxide damages biomolecules required for life increases with temperature, yet little is known about how organisms cope with this temperature-dependent threat. Here, we show that Caenorhabditis elegans nematodes use temperature information perceived by sensory neurons to cope with the temperature-dependent threat of hydrogen peroxide produced by the pathogenic bacterium Enterococcus faecium. These nematodes preemptively induce the expression of specific hydrogen peroxide defenses in response to perception of high temperature by a pair of sensory neurons. These neurons communicate temperature information to target tissues expressing those defenses via an insulin/IGF1 hormone. This is the first example of a multicellular organism inducing their defenses to a chemical when they sense an inherent enhancer of the reactivity of that chemical.