Two new functions in the WormBase Enrichment Suite

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.

C. elegans: a prominent platform for modeling and drug screening in neurological disorders

INTRODUCTION

Human neurodevelopmental and neurodegenerative diseases (NDevDs and NDegDs, respectively) encompass a broad spectrum of disorders affecting the nervous system with an increasing incidence. In this context, the nematode C. elegans, has emerged as a benchmark model for biological research, especially in the field of neuroscience.

AREAS COVERED

The authors highlight the numerous advantages of this tiny worm as a model for exploring nervous system pathologies and as a platform for drug discovery. There is a particular focus given to describing the existing models of C. elegans for the study of NDevDs and NDegDs. Specifically, the authors underscore their strong applicability in preclinical drug development. Furthermore, they place particular emphasis on detailing the common techniques employed to explore the nervous system in both healthy and diseased states.

EXPERT OPINION

Drug discovery constitutes a long and expensive process. The incorporation of invertebrate models, such as C. elegans, stands as an exemplary strategy for mitigating costs and expediting timelines. The utilization of C. elegans as a platform to replicate nervous system pathologies and conduct high-throughput automated assays in the initial phases of drug discovery is pivotal for rendering therapeutic options more attainable and cost-effective.

A pair of RNA binding proteins inhibit ion transporter expression to maintain lifespan

Regulation of lifespan by transcription factors has been well established. More recently, a role for RNA binding proteins (RBPs) in regulating lifespan has also emerged. In both cases, a major challenge is to determine which regulatory targets are functionally responsible for the observed lifespan phenotype. We recently identified a pair of neuronal RBPs, exc-7/ELAVL and mbl-1/Muscleblind, which in Caenorhabditis elegans display synthetic (nonadditive) lifespan defects: single mutants do not affect lifespan, but exc-7; mbl-1 double mutants have strongly reduced lifespan. Such a strong synthetic phenotype represented an opportunity to use transcriptomics to search for potential causative targets that are synthetically regulated. Focus on such genes would allow us to narrow our target search by ignoring the hundreds of genes altered only in single mutants, and provide a shortlist of synthetically regulated candidate targets that might be responsible for the double mutant phenotype. We identified a small handful of genes synthetically dysregulated in double mutants and systematically tested each candidate gene for functional contribution to the exc-7; mbl-1 lifespan phenotype. We identified 1 such gene, the ion transporter nhx-6, which is highly upregulated in double mutants. Overexpression of nhx-6 causes reduced lifespan, and deletion of nhx-6 in an exc-7; mbl-1 background partially restores both lifespan and healthspan. Together, these results reveal that a pair of RBPs mediate lifespan in part by inhibiting expression of an ion transporter, and provide a template for how synthetic phenotypes (including lifespan) can be dissected at the transcriptomic level to reveal potential causative genes.

A spatiotemporally resolved atlas of mRNA decay in the C. elegans embryo reveals differential regulation of mRNA stability across stages and cell types

During embryonic development, cells undergo dynamic changes in gene expression that are required for appropriate cell fate specification. Although both transcription and mRNA degradation contribute to gene expression dynamics, patterns of mRNA decay are less well-understood. Here we directly measured spatiotemporally resolved mRNA decay rates transcriptome-wide throughout C. elegans embryogenesis by transcription inhibition followed by bulk and single-cell RNA-sequencing. This allowed us to calculate mRNA half-lives within specific cell types and developmental stages and identify differentially regulated mRNA decay throughout embryonic development. We identified transcript features that are correlated with mRNA stability and found that mRNA decay rates are associated with distinct peaks in gene expression over time. Moreover, we provide evidence that, on average, mRNA is more stable in the germline compared to in the soma and in later embryonic stages compared to in earlier stages. This work suggests that differential mRNA decay across cell states and time helps to shape developmental gene expression, and it provides a valuable resource for studies of mRNA turnover regulatory mechanisms.

Starvation-induced changes in somatic insulin/IGF-1R signaling drive metabolic programming across generations

Exposure to adverse nutritional and metabolic environments during critical periods of development can exert long-lasting effects on health outcomes of an individual and its descendants. Although such metabolic programming has been observed in multiple species and in response to distinct nutritional stressors, conclusive insights into signaling pathways and mechanisms responsible for initiating, mediating, and manifesting changes to metabolism and behavior across generations remain scarce. By using a starvation paradigm in Caenorhabditis elegans, we show that starvation-induced changes in dauer formation-16/forkhead box transcription factor class O (DAF-16/FoxO) activity, the main downstream target of insulin/insulin-like growth factor 1 (IGF-1) receptor signaling, are responsible for metabolic programming phenotypes. Tissue-specific depletion of DAF-16/FoxO during distinct developmental time points demonstrates that DAF-16/FoxO acts in somatic tissues, but not directly in the germline, to both initiate and manifest metabolic programming. In conclusion, our study deciphers multifaceted and critical roles of highly conserved insulin/IGF-1 receptor signaling in determining health outcomes and behavior across generations.

A comprehensive survey of C. elegans argonaute proteins reveals organism-wide gene regulatory networks and functions

Argonaute (AGO) proteins associate with small RNAs to direct their effector function on complementary transcripts. The nematode Caenorhabditis elegans contains an expanded family of 19 functional AGO proteins, many of which have not been fully characterized. In this work, we systematically analyzed every C. elegans AGO using CRISPR-Cas9 genome editing to introduce GFP::3xFLAG tags. We have characterized the expression patterns of each AGO throughout development, identified small RNA binding complements, and determined the effects of ago loss on small RNA populations and developmental phenotypes. Our analysis indicates stratification of subsets of AGOs into distinct regulatory modules, and integration of our data led us to uncover novel stress-induced fertility and pathogen response phenotypes due to ago loss.

The nematode Oscheius tipulae as a genetic model for programmed DNA elimination

Programmed DNA elimination (PDE) is a notable exception to the paradigm of genome integrity. In metazoa, PDE often occurs coincident with germline to somatic cell differentiation. During PDE, portions of genomic DNA are lost, resulting in reduced somatic genomes. Prior studies have described the sequences lost, as well as chromosome behavior, during metazoan PDE. However, a system for studying the mechanisms and consequences of PDE in metazoa is lacking. Here, we present a functional and genetic model for PDE in the free-living Rhabditidae nematode Oscheius tipulae, a family that also includes Caenorhabditis elegans. O. tipulae was recently suggested to eliminate DNA. Using staged embryos and DNA FISH, we showed that O. tipulae PDE occurs during embryogenesis at the 8-16 cell stages. We identified a conserved motif, named Sequence For Elimination (SFE), for all 12 break sites on the six chromosomes at the junctions of retained and eliminated DNA. SFE mutants exhibited a "fail-to-eliminate" phenotype only at the modified sites. END-seq revealed that breaks can occur at multiple positions within the SFE, with extensive end resection followed by telomere addition to both retained and eliminated ends. We identified many functional SFEs at the chromosome ends through END-seq in the wild-type embryos, genome sequencing of SFE mutants, and comparative genomics of 23 wild isolates. We suggest that these alternative SFEs provide flexibility in the sequences eliminated and a fail-safe mechanism for PDE. These studies establish O. tipulae as a new, attractive model for studying the mechanisms and consequences of PDE in a metazoan.

Alternative somatic and germline gene-regulatory strategies during starvation-induced developmental arrest

Nutrient availability governs growth and quiescence, and many animals arrest development when starved. Using C. elegans L1 arrest as a model, we show that gene expression changes deep into starvation. Surprisingly, relative expression of germline-enriched genes increases for days. We conditionally degrade the large subunit of RNA polymerase II using the auxin-inducible degron system and analyze absolute expression levels. We find that somatic transcription is required for survival, but the germline maintains transcriptional quiescence. Thousands of genes are continuously transcribed in the soma, though their absolute abundance declines, such that relative expression of germline transcripts increases given extreme transcript stability. Aberrantly activating transcription in starved germ cells compromises reproduction, demonstrating important physiological function of transcriptional quiescence. This work reveals alternative somatic and germline gene-regulatory strategies during starvation, with the soma maintaining a robust transcriptional response to support survival and the germline maintaining transcriptional quiescence to support future reproductive success.

Webster et al. show that the transcriptional response to starvation is mounted early in larval somatic cells supporting survival but that it wanes over time. In contrast, they show that the germline remains transcriptionally quiescent deep into starvation, supporting reproductive potential, while maintaining its transcriptome via transcript stability.

C. elegans monitor energy status via the AMPK pathway to trigger innate immune responses against bacterial pathogens

Pathogen recognition and the triggering of host innate immune system are critical to understanding pathogen-host interaction. Cellular surveillance systems have been identified as an important strategy for the identification of microbial infection. In the present study, using Bacillus thuringiensis-Caenorhabditis elegans as a model, we found an approach for surveillance systems to sense pathogens. We report that Bacillus thuringiensis Cry5Ba, a typical pore-forming toxin, caused mitochondrial damage and energy imbalance by triggering potassium ion leakage, instead of directly targeting mitochondria. Interestingly, we find C. elegans can monitor intracellular energy status to trigger innate immune responses via AMP-activated protein kinase (AMPK), secreting multiple effectors to defend against pathogenic attacks. Our study indicates that the imbalance of energy status is a prevalent side effect of pathogen infection. Furthermore, the AMPK-dependent surveillance system may serve as a practicable strategy for the host to recognize and defense against pathogens.

Bacillus thuringiensis toxin Cry5Ba triggers potassium ion leakage, causing mitochondrial damage and energy imbalance. C. elegans can monitor this intracellular energy imbalance via AMP-activated protein kinase to trigger innate immune responses.

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.

Caenorhabditis elegans transposable elements harbor diverse transcription factor DNA-binding sites

Transposable elements are powerful agents of evolution that can diversify transcriptional programs by distributing transcription factor DNA-binding sites throughout genomes. To investigate the extent that transposable elements provide transcription factor-binding motifs in Caenorhabditis elegans, we determined the genomic positions of DNA-binding motifs for 201 different transcription factors. Surprisingly, we found that almost all examined transcription factors have binding motifs that reside within transposable elements, and all types of transposable elements have at least 1 instance of a transcription factor motif, demonstrating that transposable elements provide previously unappreciated numbers of transcription factor-binding motifs to the C. elegans genome. After determining the occurrence of transcription factor motifs in transposable elements relative to the rest of the genome, we identified DNA-binding motifs for 45 different transcription factors that are greater than 20-fold enriched within transposable elements compared to what would be expected by chance. Consistent with potential functional roles for these transposable element-enriched transcription factor-binding sequences, we determined that all transcription factor motif types found in transposable elements have instances of residing within accessible chromatin sites associated with transcription factor binding. The overwhelming majority of transcription factor-binding motifs located within transposable elements associate with their cognate transcription factors, suggesting extensive binding of transcription factors to sequences within transposable elements. In addition, transposable elements with accessible or transcription factor-bound motifs reside in the putative promoter regions of approximately 12% of all protein-coding genes, providing widespread possibilities for influencing gene expression. This work represents the first comprehensive analysis of transposable element–transcription factor interactions in C. elegans and demonstrates that transposable element-provided transcription factor-binding sites are prevalent in this important model organism.

Targeted DamID in C. elegans reveals a direct role for LIN-22 and NHR-25 in antagonizing the epidermal stem cell fate

Transcription factors are key players in gene networks controlling cell fate specification during development. In multicellular organisms, they display complex patterns of expression and binding to their targets, hence, tissue specificity is required in the characterization of transcription factor-target interactions. We introduce here targeted DamID (TaDa) as a method for tissue-specific transcription factor target identification in intact Caenorhabditis elegans animals. We use TaDa to recover targets in the epidermis for two factors, the HES1 homolog LIN-22, and the NR5A1/2 nuclear hormone receptor NHR-25. We demonstrate a direct link between LIN-22 and the Wnt signaling pathway through repression of the Frizzled receptor lin-17. We report a direct role for NHR-25 in promoting cell differentiation via repressing the expression of stem cell-promoting GATA factors. Our results expand our understanding of the epidermal gene network and highlight the potential of TaDa to dissect the architecture of tissue-specific gene regulatory networks.

WormBase in 2022-data, processes, and tools for analyzing Caenorhabditis elegans

WormBase (www.wormbase.org) is the central repository for the genetics and genomics of the nematode Caenorhabditis elegans. We provide the research community with data and tools to facilitate the use of C. elegans and related nematodes as model organisms for studying human health, development, and many aspects of fundamental biology. Throughout our 22-year history, we have continued to evolve to reflect progress and innovation in the science and technologies involved in the study of C. elegans. We strive to incorporate new data types and richer data sets, and to provide integrated displays and services that avail the knowledge generated by the published nematode genetics literature. Here, we provide a broad overview of the current state of WormBase in terms of data type, curation workflows, analysis, and tools, including exciting new advances for analysis of single-cell data, text mining and visualization, and the new community collaboration forum. Concurrently, we continue the integration and harmonization of infrastructure, processes, and tools with the Alliance of Genome Resources, of which WormBase is a founding member.

A genetic screen in C. elegans reveals roles for KIN17 and PRCC in maintaining 5' splice site identity

Pre-mRNA splicing is an essential step of eukaryotic gene expression carried out by a series of dynamic macromolecular protein/RNA complexes, known collectively and individually as the spliceosome. This series of spliceosomal complexes define, assemble on, and catalyze the removal of introns. Molecular model snapshots of intermediates in the process have been created from cryo-EM data, however, many aspects of the dynamic changes that occur in the spliceosome are not fully understood. Caenorhabditis elegans follow the GU-AG rule of splicing, with almost all introns beginning with 5’ GU and ending with 3’ AG. These splice sites are identified early in the splicing cycle, but as the cycle progresses and “custody” of the pre-mRNA splice sites is passed from factor to factor as the catalytic site is built, the mechanism by which splice site identity is maintained or re-established through these dynamic changes is unclear. We performed a genetic screen in C. elegans for factors that are capable of changing 5’ splice site choice. We report that KIN17 and PRCC are involved in splice site choice, the first functional splicing role proposed for either of these proteins. Previously identified suppressors of cryptic 5’ splicing promote distal cryptic GU splice sites, however, mutations in KIN17 and PRCC instead promote usage of an unusual proximal 5’ splice site which defines an intron beginning with UU, separated by 1nt from a GU donor. We performed high-throughput mRNA sequencing analysis and found that mutations in PRCC, and to a lesser extent KIN17, changed alternative 5’ splice site usage at native sites genome-wide, often promoting usage of nearby non-consensus sites. Our work has uncovered both fine and coarse mechanisms by which the spliceosome maintains splice site identity during the complex assembly process.

Pre-messenger RNA splicing is an important regulator of eukaryotic gene expression, changing the content, frame, and functionality of both coding and non-coding transcripts. Our understanding of how the spliceosome chooses where to cut has focused on the initial identification of splice sites. However, our results suggest that the spliceosome also relies on other components in later steps to maintain the identity of the splice donor sites. We are currently in the midst of a “resolution revolution”, with ever-clearer cryo-EM snapshots of stalled complexes, allowing researchers to visualize moments in time in the splicing cycle. These models are illuminating, but do not always elucidate mechanistic functioning of a highly dynamic ribonucleoprotein complex. Therefore, our lab takes a complementary approach, using the power of genetics in a multicellular animal to gain functional insights into the spliceosome. Using a C.elegans genetic screen, we have found novel functional splicing roles for two proteins, KIN17 and PRCC. Mutations in PRCC in particular promote nearby alternative 5’ splice sites at native loci. This work improves our understanding of how the spliceosome maintains the identity of where to cut the pre-mRNA, and thus how genes are expressed and used in multicellular animals.

MOET: a web-based gene set enrichment tool at the Rat Genome Database for multiontology and multispecies analyses

Biological interpretation of a large amount of gene or protein data is complex. Ontology analysis tools are imperative in finding functional similarities through overrepresentation or enrichment of terms associated with the input gene or protein lists. However, most tools are limited by their ability to do ontology-specific and species-limited analyses. Furthermore, some enrichment tools are not updated frequently with recent information from databases, thus giving users inaccurate, outdated or uninformative data. Here, we present MOET or the Multi-Ontology Enrichment Tool (v.1 released in April 2019 and v.2 released in May 2021), an ontology analysis tool leveraging data that the Rat Genome Database (RGD) integrated from in-house expert curation and external databases including the National Center for Biotechnology Information (NCBI), Mouse Genome Informatics (MGI), The Kyoto Encyclopedia of Genes and Genomes (KEGG), The Gene Ontology Resource, UniProt-GOA, and others. Given a gene or protein list, MOET analysis identifies significantly overrepresented ontology terms using a hypergeometric test and provides nominal and Bonferroni corrected P-values and odds ratios for the overrepresented terms. The results are shown as a downloadable list of terms with and without Bonferroni correction, and a graph of the P-values and number of annotated genes for each term in the list. MOET can be accessed freely from https://rgd.mcw.edu/rgdweb/enrichment/start.html.

An autoregulation loop in fust-1 for circular RNA regulation in Caenorhabditis elegans

Many circular RNAs (circRNAs) are differentially expressed in different tissues or cell types, suggestive of specific factors that regulate their biogenesis. Here, taking advantage of available mutation strains of RNA-binding proteins (RBPs) in Caenorhabditis elegans, I performed a screening of circRNA regulation in 13 conserved RBPs. Among them, loss of FUST-1, the homolog of Fused in Sarcoma (FUS), caused downregulation of multiple circRNAs. By rescue experiments, I confirmed FUST-1 as a circRNA regulator. Through RNA sequencing using circRNA-enriched samples, circRNAs targets regulated by FUST-1 were identified globally, with hundreds of them significantly altered. Furthermore, I showed that FUST-1 regulates circRNA formation with only small to little effect on the cognate linear mRNAs. When recognizing circRNA pre-mRNAs, FUST-1 can affect both exon-skipping and circRNA in the same genes. Moreover, I identified an autoregulation loop in fust-1, where FUST-1, isoform a (FUST-1A) promotes the skipping of exon 5 of its own pre-mRNA, which produces FUST-1, isoform b (FUST-1B) with different N-terminal sequences. FUST-1A is the functional isoform in circRNA regulation. Although FUST-1B has the same functional domains as FUST-1A, it cannot regulate either exon-skipping or circRNA formation. This study provided an in vivo investigation of circRNA regulation, which will be helpful to understand the mechanisms that govern circRNA formation.

DREAM represses distinct targets by cooperating with different THAP domain proteins

The DREAM (dimerization partner [DP], retinoblastoma [Rb]-like, E2F, and MuvB) complex controls cellular quiescence by repressing cell-cycle and other genes, but its mechanism of action is unclear. Here, we demonstrate that two C. elegans THAP domain proteins, LIN-15B and LIN-36, co-localize with DREAM and function by different mechanisms for repression of distinct sets of targets. LIN-36 represses classical cell-cycle targets by promoting DREAM binding and gene body enrichment of H2A.Z, and we find that DREAM subunit EFL-1/E2F is specific for LIN-36 targets. In contrast, LIN-15B represses germline-specific targets in the soma by facilitating H3K9me2 promoter marking. We further find that LIN-36 and LIN-15B differently regulate DREAM binding. In humans, THAP proteins have been implicated in cell-cycle regulation by poorly understood mechanisms. We propose that THAP domain proteins are key mediators of Rb/DREAM function.

Region-specific H3K9me3 gain in aged somatic tissues in Caenorhabditis elegans

Epigenetic alterations occur as organisms age, and lead to chromatin deterioration, loss of transcriptional silencing and genomic instability. Dysregulation of the epigenome has been associated with increased susceptibility to age-related disorders. In this study, we aimed to characterize the age-dependent changes of the epigenome and, in turn, to understand epigenetic processes that drive aging phenotypes. We focused on the aging-associated changes in the repressive histone marks H3K9me3 and H3K27me3 in C. elegans. We observed region-specific gain and loss of both histone marks, but the changes are more evident for H3K9me3. We further found alteration of heterochromatic boundaries in aged somatic tissues. Interestingly, we discovered that the most statistically significant changes reflected H3K9me3-marked regions that are formed during aging, and are absent in developing worms, which we termed "aging-specific repressive regions" (ASRRs). These ASRRs preferentially occur in genic regions that are marked by high levels of H3K9me2 and H3K36me2 in larval stages. Maintenance of high H3K9me2 levels in these regions have been shown to correlate with a longer lifespan. Next, we examined whether the changes in repressive histone marks lead to de-silencing of repetitive DNA elements, as reported for several other organisms. We observed increased expression of active repetitive DNA elements but not global re-activation of silent repeats in old worms, likely due to the distributed nature of repetitive elements in the C. elegans genome. Intriguingly, CELE45, a putative short interspersed nuclear element (SINE), was greatly overexpressed at old age and upon heat stress. SINEs have been suggested to regulate transcription in response to various cellular stresses in mammals. It is likely that CELE45 RNAs also play roles in stress response and aging in C. elegans. Taken together, our study revealed significant and specific age-dependent changes in repressive histone modifications and repetitive elements, providing important insights into aging biology.

Gene expression profiling of epidermal cell types in C. elegans using Targeted DamID

The epidermis of Caenorhabditis elegans is an essential tissue for survival because it contributes to the formation of the cuticle barrier as well as facilitating developmental progression and animal growth. Most of the epidermis consists of the hyp7 hypodermal syncytium, the nuclei of which are largely generated by the seam cells, which exhibit stem cell-like behaviour during development. How seam cell progenitors differ transcriptionally from the differentiated hypodermis is poorly understood. Here, we introduce Targeted DamID (TaDa) in C. elegans as a method for identifying genes expressed within a tissue of interest without cell isolation. We show that TaDa signal enrichment profiles can be used to identify genes transcribed in the epidermis and use this method to resolve differences in gene expression between the seam cells and the hypodermis. Finally, we predict and functionally validate new transcription and chromatin factors acting in seam cell development. These findings provide insights into cell type-specific gene expression profiles likely associated with epidermal cell fate patterning.

Reverse complementary matches simultaneously promote both back-splicing and exon-skipping

Background

Circular RNAs (circRNAs) play diverse roles in different biological and physiological environments and are always expressed in a tissue-specific manner. Especially, circRNAs are enriched in the brain tissues of almost all investigated species, including humans, mice, Drosophila, etc. Although circRNAs were found in C. elegans, the neuron-specific circRNA data is not available yet. Exon-skipping is found to be correlated to circRNA formation, but the mechanisms that link them together are not clear.

Results

Here, through large-scale neuron isolation from the first larval (L1) stage of C. elegans followed by RNA sequencing with ribosomal RNA depletion, the neuronal circRNA data in C. elegans were obtained. Hundreds of novel circRNAs were annotated with high accuracy. circRNAs were highly expressed in the neurons of C. elegans and were positively correlated to the levels of their cognate linear mRNAs. Disruption of reverse complementary match (RCM) sequences in circRNA flanking introns effectively abolished circRNA formation. In the zip-2 gene, deletion of either upstream or downstream RCMs almost eliminated the production of both the circular and the skipped transcript. Interestingly, the 13-nt RCM in zip-2 is highly conserved across five nematode ortholog genes, which show conserved exon-skipping patterns. Finally, through in vivo one-by-one mutagenesis of all the splicing sites and branch points required for exon-skipping and back-splicing in the zip-2 gene, I showed that back-splicing still happened without exon-skipping, and vice versa.

Conclusions

Through protocol optimization, total RNA obtained from sorted neurons is increased to hundreds of nanograms. circRNAs highly expressed in the neurons of C. elegans are more likely to be derived from genes also highly expressed in the neurons. RCMs are abundant in circRNA flanking introns, and RCM-deletion is an efficient way to knockout circRNAs. More importantly, these RCMs are not only required for back-splicing but also promote the skipping of exon(s) to be circularized. Finally, RCMs in circRNA flanking introns can directly promote both exon-skipping and back-splicing, providing a new explanation for the correlation between them.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07910-w.

Interspecies RNA Interactome of Pathogen and Host in a Heritable Defensive Strategy

Communication with bacteria deeply impacts the life history traits of their hosts. Through specific molecules and metabolites, bacteria can promote short- and long-term phenotypic and behavioral changes in the nematode Caenorhabditis elegans. The chronic exposure of C. elegans to pathogens promotes the adaptive behavior in the host’s progeny called pathogen-induced diapause formation (PIDF). PIDF is a pathogen avoidance strategy induced in the second generation of animals infected and can be recalled transgenerationally. This behavior requires the RNA interference machinery and specific nematode and bacteria small RNAs (sRNAs). In this work, we assume that RNAs from both species co-exist and can interact with each other. Under this principle, we explore the potential interspecies RNA interactions during PIDF-triggering conditions, using transcriptomic data from the holobiont. We study two transcriptomics datasets: first, the dual sRNA expression of Pseudomonas aeruginosa PAO1 and C. elegans in a transgenerational paradigm for six generations and second, the simultaneous expression of sRNAs and mRNA in intergenerational PIDF. We focus on those bacterial sRNAs that are systematically overexpressed in the intestines of animals compared with sRNAs expressed in host-naïve bacteria. We selected diverse in silico methods that represent putative mechanisms of RNA-mediated interspecies interaction. These interactions are as follows: heterologous perfect and incomplete pairing between bacterial RNA and host mRNA; sRNAs of similar sequence expressed in both species that could mimic each other; and known or predicted eukaryotic motifs present in bacterial transcripts. We conclude that a broad spectrum of tools can be applied for the identification of potential sRNA and mRNA targets of the interspecies RNA interaction that can be subsequently tested experimentally.

Powerful detection of polygenic selection and evidence of environmental adaptation in US beef cattle

Selection on complex traits can rapidly drive evolution, especially in stressful environments. This polygenic selection does not leave intense sweep signatures on the genome, rather many loci experience small allele frequency shifts, resulting in large cumulative phenotypic changes. Directional selection and local adaptation are changing populations; but, identifying loci underlying polygenic or environmental selection has been difficult. We use genomic data on tens of thousands of cattle from three populations, distributed over time and landscapes, in linear mixed models with novel dependent variables to map signatures of selection on complex traits and local adaptation. We identify 207 genomic loci associated with an animal's birth date, representing ongoing selection for monogenic and polygenic traits. Additionally, hundreds of additional loci are associated with continuous and discrete environments, providing evidence for historical local adaptation. These candidate loci highlight the nervous system's possible role in local adaptation. While advanced technologies have increased the rate of directional selection in cattle, it has likely been at the expense of historically generated local adaptation, which is especially problematic in changing climates. When applied to large, diverse cattle datasets, these selection mapping methods provide an insight into how selection on complex traits continually shapes the genome. Further, understanding the genomic loci implicated in adaptation may help us breed more adapted and efficient cattle, and begin to understand the basis for mammalian adaptation, especially in changing climates. These selection mapping approaches help clarify selective forces and loci in evolutionary, model, and agricultural contexts.

hsp-90 and unc-45 depletion induce characteristic transcriptional signatures in coexpression cliques of C. elegans

Nematode development is characterized by progression through several larval stages. Thousands of genes were found in large scale RNAi-experiments to block this development at certain steps, two of which target the molecular chaperone HSP-90 and its cofactor UNC-45. Aiming to define the cause of arrest, we here investigate the status of nematodes after treatment with RNAi against hsp-90 and unc-45 by employing an in-depth transcriptional analysis of the arrested larvae. To identify misregulated transcriptional units, we calculate and validate genome-wide coexpression cliques covering the entire nematode genome. We define 307 coexpression cliques and more than half of these can be related to organismal functions by GO-term enrichment, phenotype enrichment or tissue enrichment analysis. Importantly, hsp-90 and unc-45 RNAi induce or repress many of these cliques in a coordinated manner, and then several specifically regulated cliques are observed. To map the developmental state of the arrested nematodes we define the expression behaviour of each of the cliques during development from embryo to adult nematode. hsp-90 RNAi can be seen to arrest development close to the L4 larval stage with further deviations in daf-16 regulated genes. unc-45 RNAi instead leads to arrested development at young adult stage prior to the programmatic downregulation of sperm-cell specific genes. In both cases processes can be defined to be misregulated upon depletion of the respective chaperone. With most of the defined gene cliques showing concerted behaviour at some stage of development from embryo to late adult, the “clique map” together with the clique-specific GO-terms, tissue and phenotype assignments will be a valuable tool in understanding concerted responses on the genome-wide level in Caenorhabditis elegans.

SID-2 negatively regulates development likely independent of nutritional dsRNA uptake

RNA interference (RNAi) is a gene regulatory mechanism based on RNA-RNA interaction conserved through eukaryotes. Surprisingly, many animals can take-up human-made double stranded RNA (dsRNA) from the environment to initiate RNAi suggesting a mechanism for dsRNA-based information exchange between organisms and their environment. However, no naturally occurring example has been identified since the discovery of the phenomenon 22 years ago. Therefore it remains enigmatic why animals are able to take up dsRNA. Here, we explore other possible functions by performing phenotypic studies of dsRNA uptake deficient sid-2 mutants in Caenorhabditis elegans. We find that SID-2 does not have a nutritional role in feeding experiments using genetic sensitized mutants. Furthermore, we use robot assisted imaging to show that sid-2 mutants accelerate growth rate and, by maternal contribution, body length at hatching. Finally, we perform transcriptome and lipidome analysis showing that sid-2 has no effect on energy storage lipids, but affects signalling lipids and the embryo transcriptome. Overall, these results suggest that sid-2 has mild effects on development and is unlikely functioning in the nutritional uptake of dsRNA. These findings broaden our understanding of the biological role of SID-2 and motivate studies identifying the role of environmental dsRNA uptake.

Antagonistic fungal enterotoxins intersect at multiple levels with host innate immune defences

Animals and plants need to defend themselves from pathogen attack. Their defences drive innovation in virulence mechanisms, leading to never-ending cycles of co-evolution in both hosts and pathogens. A full understanding of host immunity therefore requires examination of pathogen virulence strategies. Here, we take advantage of the well-studied innate immune system of Caenorhabditis elegans to dissect the action of two virulence factors from its natural fungal pathogen Drechmeria coniospora. We show that these two enterotoxins have strikingly different effects when expressed individually in the nematode epidermis. One is able to interfere with diverse aspects of host cell biology, altering vesicle trafficking and preventing the key STAT-like transcription factor STA-2 from activating defensive antimicrobial peptide gene expression. The second increases STA-2 levels in the nucleus, modifies the nucleolus, and, potentially as a consequence of a host surveillance mechanism, causes increased defence gene expression. Our results highlight the remarkably complex and potentially antagonistic mechanisms that come into play in the interaction between co-evolved hosts and pathogens.

Chronic temperature stress inhibits reproduction and disrupts endocytosis via chaperone titration in Caenorhabditis elegans

Background

Temperature influences biology at all levels, from altering rates of biochemical reactions to determining sustainability of entire ecosystems. Although extended exposure to elevated temperatures influences organismal phenotypes important for human health, agriculture, and ecology, the molecular mechanisms that drive these responses remain largely unexplored. Prolonged, mild temperature stress (48 h at 28 °C) has been shown to inhibit reproduction in Caenorhabditis elegans without significantly impacting motility or viability.

Results

Analysis of molecular responses to chronic stress using RNA-seq uncovers dramatic effects on the transcriptome that are fundamentally distinct from the well-characterized, acute heat shock response (HSR). While a large portion of the genome is differentially expressed ≥ 4-fold after 48 h at 28 °C, the only major class of oogenesis-associated genes affected is the vitellogenin gene family that encodes for yolk proteins (YPs). Whereas YP mRNAs decrease, the proteins accumulate and mislocalize in the pseudocoelomic space as early as 6 h, well before reproduction declines. A trafficking defect in a second, unrelated fluorescent reporter and a decrease in pre-synaptic neuronal signaling indicate that the YP mislocalization is caused by a generalized defect in endocytosis. Molecular chaperones are involved in both endocytosis and refolding damaged proteins. Decreasing levels of the major HSP70 chaperone, HSP-1, causes similar YP trafficking defects in the absence of stress. Conversely, increasing chaperone levels through overexpression of the transcription factor HSF-1 rescues YP trafficking and restores neuronal signaling.

Conclusions

These data implicate chaperone titration during chronic stress as a molecular mechanism contributing to endocytic defects that influence multiple aspects of organismal physiology. Notably, HSF-1 overexpression improves recovery of viable offspring after exposure to stress. These findings provide important molecular insights into understanding organismal responses to temperature stress as well as phenotypes associated with chronic protein misfolding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01008-1.

Large genetic diversity and strong positive selection in F-box and GPCR genes among the wild isolates of Caenorhabditis elegans

The F-box and chemosensory GPCR (csGPCR) gene families are greatly expanded in nematodes, including the model organism Caenorhabditis elegans, compared with insects and vertebrates. However, the intraspecific evolution of these two gene families in nematodes remain unexamined. In this study, we analyzed the genomic sequences of 330 recently sequenced wild isolates of C. elegans using a range of population genetics approaches. We found that F-box and csGPCR genes, especially the Srw family csGPCRs, showed much more diversity than other gene families. Population structure analysis and phylogenetic analysis divided the wild strains into eight non-Hawaiian and three Hawaiian subpopulations. Some Hawaiian strains appeared to be more ancestral than all other strains. F-box and csGPCR genes maintained a great amount of the ancestral variants in the Hawaiian subpopulation and their divergence among the non-Hawaiian subpopulations contributed significantly to population structure. F-box genes are mostly located at the chromosomal arms and high recombination rate correlates with their large polymorphism. Moreover, using both neutrality tests and extended haplotype homozygosity analysis, we identified signatures of strong positive selection in the F-box and csGPCR genes among the wild isolates, especially in the non-Hawaiian population. Accumulation of high-frequency-derived alleles in these genes was found in non-Hawaiian population, leading to divergence from the ancestral genotype. In summary, we found that F-box and csGPCR genes harbor a large pool of natural variants, which may be subjected to positive selection. These variants are mostly mapped to the substrate-recognition domains of F-box proteins and the extracellular and intracellular regions of csGPCRs, possibly resulting in advantages during adaptation by affecting protein degradation and the sensing of environmental cues, respectively.

FHOD formin and SRF promote post-embryonic striated muscle growth through separate pathways in C. elegans

Previous work with cultured cells has shown transcription of muscle genes by serum response factor (SRF) can be stimulated by actin polymerization driven by proteins of the formin family. However, it is not clear if endogenous formins similarly promote SRF-dependent transcription during muscle development in vivo. We tested whether formin activity promotes SRF-dependent transcription in striated muscle in the simple animal model, Caenorhabditis elegans. Our lab has shown FHOD-1 is the only formin that directly promotes sarcomere formation in the worm's striated muscle. We show here FHOD-1 and SRF homolog UNC-120 both support muscle growth and also muscle myosin II heavy chain A expression. However, while a hypomorphic unc-120 allele blunts expression of a set of striated muscle genes, these genes are largely upregulated or unchanged by absence of FHOD-1. Instead, pharmacological inhibition of the proteasome restores myosin protein levels in worms lacking FHOD-1, suggesting elevated proteolysis accounts for their myosin deficit. Interestingly, proteasome inhibition does not restore normal muscle growth to fhod-1(Δ) mutants, suggesting formin contributes to muscle growth by some alternative mechanism. Overall, we find SRF does not depend on formin to promote muscle gene transcription in a simple in vivo system.

Intergenerational Pathogen-Induced Diapause in Caenorhabditis elegans Is Modulated by mir-243

Persistent infection of the bacterivore nematode C. elegans with bacteria such as P. aeruginosa and S. enterica makes the worm diapause or hibernate. By doing this, the worm closes its mouth, avoiding infection. This response takes two generations to be implemented. In this work, we looked for genes expressed upon infection that could mediate the worm diapause triggered by pathogens. We identify mir-243-3p as the only transcript commonly upregulated when animals feed on P. aeruginosa and S. enterica for two consecutive generations. Moreover, we demonstrate that mir-243-3p is required for pathogen-induced dauer formation, a new function that has not been previously described for this microRNA (miRNA). We also find that the transcriptional activators DAF-16, PQM-1, and CRH-2 are necessary for the expression of mir-243 under pathogenesis. Here we establish a relationship between a small RNA and a developmental change that ensures the survival of a percentage of the progeny.

The interaction and communication between bacteria and their hosts modulate many aspects of animal physiology and behavior. Dauer entry as a response to chronic exposure to pathogenic bacteria in Caenorhabditis elegans is an example of a dramatic survival response. This response is dependent on the RNA interference (RNAi) machinery, suggesting the involvement of small RNAs (sRNAs) as effectors. Interestingly, dauer formation occurs after two generations of interaction with two unrelated moderately pathogenic bacteria. Therefore, we sought to discover the identity of C. elegans RNAs involved in pathogen-induced diapause. Using transcriptomics and differential expression analysis of coding and long and small noncoding RNAs, we found that mir-243-3p (the mature form of mir-243) is the only transcript continuously upregulated in animals exposed to both Pseudomonas aeruginosa and Salmonella enterica for two generations. Phenotypic analysis of mutants showed that mir-243 is required for dauer formation under pathogenesis but not under starvation. Moreover, DAF-16, a master regulator of defensive responses in the animal and required for dauer formation was found to be necessary for mir-243 expression. This work highlights the role of a small noncoding RNA in the intergenerational defensive response against pathogenic bacteria and interkingdom communication.

The Ancestral Caenorhabditis elegans Cuticle Suppresses rol-1

Genetic background commonly modifies the effects of mutations. We discovered that worms mutant for the canonical rol-1 gene, identified by Brenner in 1974, do not roll in the genetic background of the wild strain CB4856. Using linkage mapping, association analysis and gene editing, we determined that N2 carries an insertion in the collagen gene col-182 that acts as a recessive enhancer of rol-1 rolling. From population and comparative genomics, we infer the insertion is derived in N2 and related laboratory lines, likely arising during the domestication of Caenorhabditis elegans, and breaking a conserved protein. The ancestral version of col-182 also modifies the phenotypes of four other classical cuticle mutant alleles, and the effects of natural genetic variation on worm shape and locomotion. These results underscore the importance of genetic background and the serendipity of Brenner's choice of strain.

ngn-1/neurogenin Activates Transcription of Multiple Terminal Selector Transcription Factors in the Caenorhabditis elegans Nervous System

Proper nervous system development is required for an organism's survival and function. Defects in neurogenesis have been linked to neurodevelopmental disorders such as schizophrenia and autism. Understanding the gene regulatory networks that orchestrate neural development, specifically cascades of proneural transcription factors, can better elucidate which genes are most important during early neurogenesis. Neurogenins are a family of deeply conserved factors shown to be both necessary and sufficient for the development of neural subtypes. However, the immediate downstream targets of neurogenin are not well characterized. The objective of this study was to further elucidate the role of ngn-1/neurogenin in nervous system development and to identify its downstream transcriptional targets, using the nematode Caenorhabditis elegans as a model for this work. We found that ngn-1 is required for axon outgrowth, nerve ring architecture, and neuronal cell fate specification. We also showed that ngn-1 may have roles in neuroblast migration and epithelial integrity during embryonic development. Using RNA sequencing and comparative transcriptome analysis, we identified eight transcription factors (hlh-34/NPAS1, unc-42/PROP1, ceh-17/PHOX2A, lim-4/LHX6, fax-1/NR2E3, lin-11/LHX1, tlp-1/ZNF503, and nhr-23/RORB) whose transcription is activated, either directly or indirectly, by ngn-1 Our results show that ngn-1 has a role in transcribing known terminal regulators that establish and maintain cell fate of differentiated neural subtypes and confirms that ngn-1 functions as a proneural transcription factor in C. elegans neurogenesis.

Transcriptomic Profiling of DAF-7/TGFβ Pathway Mutants in C. elegans

The transforming growth factor beta superfamily encompasses a large family of ligands that are well conserved across many organisms. They are regulators of a number of physiological and pathological processes. The model nematode Caenorhabditis elegans has been instrumental in identifying key components of the transforming growth factor beta (TGFβ) pathway. In C. elegans, the TGFβ homolog DAF-7 signals through the DAF-1 Type I and DAF-4 Type II receptors to phosphorylate downstream R-SMADs DAF-8 and DAF-14. These R-SMADs translocate into the nucleus to inhibit Co-SMAD DAF-3. Many of the roles of the canonical DAF-7 pathway, involving both DAF-1 and DAF-3, have been identified using targeted genetic studies. Few have assessed the global transcriptomic changes in response to these genes, especially in adult animals. In this study, we performed RNA sequencing on wild type, daf-1, and daf-1; daf-3 adult hermaphrodites. To assess the overall trends of the data, we identified differentially expressed genes (DEGs) and performed gene ontology analysis to identify the types of downstream genes that are differentially expressed. Hierarchical clustering showed that the daf-1; daf-3 double mutants are transcriptionally more similar to wild type than daf-1 mutants. Analysis of the DEGs showed a disproportionally high number of genes whose expression is increased in daf-1 mutants, suggesting that DAF-1 acts as a general repressor of gene expression in wild type animals. Gene ontology analysis of the DEGs produced many significantly enriched terms, including Molting Cycle, Response to Topologically Incorrect Protein, and Response to Biotic Stimulus. Understanding the direct and indirect targets of the DAF-7 TGFβ pathway through this RNA-seq dataset can provide insight into novel roles of the multifunctional signaling pathway, as well as identify novel genes that may participate in previously reported functions of TGFβ signaling.