The time-resolved transcriptome of C. elegans

Transcriptional atlas of Daphnia magna

Transcriptomics studies are more likely to achieve predictive results when they rely on tissue- and cell-specific transcriptional data. Identification of cell types in novel model organisms by their transcriptional profiles is difficult without data on transcriptional differences among major tissues and anatomical features. Here we report the first dataset on tissue- and organ-specific transcriptomics in freshwater plankton crustacean Daphnia magna, reporting markers of embryos, hemocytes, gut, carapace, antennae-2, and head, as well as the remaining carcass. Embryos are the most transcriptionally different from adults' features, with antennae and carapace being the most differentiated among them. We demonstrate that transcriptional markers of embryos vs. adults and of various adult anatomical features can be used to provide validation and functional explanation to published differential expression in response to environmental factors like infection, hypoxia, toxicants, or kairomones; to annotate Daphnia single cell data; and to ask questions about transcriptional diversification within extended gene families.

FG repeats drive co-clustering of nuclear pores and P granules in the C. elegans germline

Condensates that accumulate small RNA biogenesis factors (nuage) are common in germ cells and often associate with nuclei. In the Caenorhabditis elegans germline, P granules overlay large clusters of nuclear pores and this organization has been proposed to facilitate surveillance of nascent transcripts by Argonaute proteins enriched in P granules. We report that co-clustering of nuclear pores and P granules depends on FG repeat-containing nucleoporins and FG repeats in the Vasa class helicase GLH-1. Worms with mutations that prevent this co-clustering are fertile under standard growth conditions and exhibit misregulation of only a minority of genes, including replication-dependent histones. Our observations suggest that association with nuclear pores, although non-essential for genome surveillance, may serve to tune mRNA flow through P granules and other nuage condensates.

Widespread destabilization of Caenorhabditis elegans microRNAs by the E3 ubiquitin ligase EBAX-1

MicroRNAs (miRNAs) associate with Argonaute (AGO) proteins to form complexes that direct mRNA repression. miRNAs are also the subject of regulation. For example, some miRNAs are destabilized through a pathway in which pairing to specialized transcripts recruits the ZSWIM8 E3 ubiquitin ligase, which polyubiquitinates AGO, leading to its degradation and exposure of the miRNA to cellular nucleases. Here, we found that 22 miRNAs in Caenorhabditis elegans are sensitive to loss of EBAX-1, the ZSWIM8 ortholog in nematodes, implying that these 22 miRNAs might be subject to this pathway of target-directed miRNA degradation (TDMD). The impact of EBAX-1 depended on the developmental stage, with the greatest effect on the miRNA pool (14.5%) observed in L1 larvae, and the greatest number of different miRNAs affected (17) observed in germline-depleted adults. The affected miRNAs included the miR-35-42 family, as well as other miRNAs among the least stable in the worm, suggesting that TDMD is a major miRNA-destabilization pathway in the worm. The excess miR-35-42 molecules that accumulated in ebax-1 mutants caused increased repression of their predicted target mRNAs and underwent 3' trimming over time. In general, however, miRNAs sensitive to EBAX-1 loss had no consistent pattern of either trimming or tailing. Replacement of the 3' region of miR-43 substantially reduced EBAX-1 sensitivity, a result that differed from that observed previously for miR-35. Together, these findings broaden the implied biological scope of TDMD-like regulation of miRNA stability in animals, and indicate that a role for miRNA 3' sequences is variable in the worm.

Multi-level transcriptional regulation of embryonic sex determination and dosage compensation by the X-signal element sex-1

The C. elegans nuclear hormone receptor sex-1 is known to be an embryonic X-signal element that represses xol-1, the sex-switch gene that is the master regulator of sex determination and dosage compensation. Several prior studies on sex-1 function have suggested that sex-1 may have additional downstream roles beyond the regulation of xol-1 expression. In this study we characterize some of these additional roles of sex-1 in regulating the dual processes of sex determination and dosage compensation during embryogenesis. Our study reveals that sex-1 acts on many of the downstream targets of xol-1 in a xol-1-independent manner. Further analysis of these shared but independently regulated downstream targets uncovered that sex-1 mediates the expression of hermaphrodite- and male-biased genes during embryogenesis. We validated sex-1 binding on one of these downstream targets, the male-developmental gene her-1. Our data suggests a model where sex-1 exhibits multi-level direct transcriptional regulation on several targets, including xol-1 and genes downstream of xol-1, to reinforce the appropriate expression of sex-biased transcripts in XX embryos. Furthermore, we found that xol-1 sex-1 double mutants show defects in dosage compensation. Our study provides evidence that misregulation of dpy-21, one of the components of the dosage compensation complex, and the subsequent misregulation of H4K20me1 enrichment on the X chromosomes, may contribute to this defect.

EDC-3 and EDC-4 regulate embryonic mRNA clearance and biomolecular condensate specialization

Animal development is dictated by the selective and timely decay of mRNAs in developmental transitions, but the impact of mRNA decapping scaffold proteins in development is unclear. This study unveils the roles and interactions of the DCAP-2 decapping scaffolds EDC-3 and EDC-4 in the embryonic development of C. elegans. EDC-3 facilitates the timely removal of specific embryonic mRNAs, including cgh-1, car-1, and ifet-1 by reducing their expression and preventing excessive accumulation of DCAP-2 condensates in somatic cells. We further uncover a role for EDC-3 in defining the boundaries between P bodies, germ granules, and stress granules. Finally, we show that EDC-4 counteracts EDC-3 and engenders the assembly of DCAP-2 with the GID (CTLH) complex, a ubiquitin ligase involved in maternal-to-zygotic transition (MZT). Our findings support a model where multiple RNA decay mechanisms temporally clear maternal and zygotic mRNAs throughout embryonic development.

Transcriptome analysis in C. elegans early embryos upon depletion of the Topoisomerase 2/condensin II axis

In C. elegans , the chromosome compaction factors topoisomerase 2 (TOP-2) and condensin II have been shown to globally repress transcription in multiple contexts. Our group has previously reported that TOP-2 and condensin II repress transcription in the C. elegans germline during larval starvation, oocyte maturation, and in germline progenitor cells of the early embryo. Here, we assess the transcriptome of early embryos treated with RNAi against TOP-2 and the condensin II subunit CAPG-2. We found 144 upregulated and 172 downregulated genes. Further analysis showed that the upregulated genes are mostly somatic, with a host of neuronal cells present in our tissue enrichment analysis.

Identification of shared gene expression programs activated in multiple modes of torpor across vertebrate clades

Torpor encompasses diverse adaptations to extreme environmental stressors such as hibernation, aestivation, brumation, and daily torpor. Here we introduce StrokeofGenus, an analytic pipeline that identifies distinct transcriptomic states and shared gene expression patterns across studies, tissues, and species. We use StrokeofGenus to study multiple and diverse forms of torpor from publicly-available RNA-seq datasets that span eight species and two classes. We identify three transcriptionally distinct states during the cycle of heterothermia: euthermia, torpor, and interbout arousal. We also identify torpor-specific gene expression patterns that are shared both across tissues and between species with over three hundred million years of evolutionary divergence. We further demonstrate the general sharing of gene expression patterns in multiple forms of torpor, implying a common evolutionary origin for this process. Although here we apply StrokeofGenus to analysis of torpor, it can be used to interrogate any other complex physiological processes defined by transient transcriptomic states.

Understanding Longevity: SIN-3 and DAF-16 Revealed as Independent Players in Lifespan Regulation

Aging is the process of gradual physio-biochemical deterioration. Although aging is inevitable, healthy aging is the key to individual and communal well-being. Therefore, it is essential to understand the regulation of aging. SIN-3/Sin-3 is a unique regulatory protein that regulates aging without DNA-binding activity. It functions by establishing multiple protein interactions. To understand the functional mechanism of this transcriptional regulator, the Caenorhabditis elegans protein interactome was assessed for SIN-3 interactions. DAF-16/FOXO emerged as one of the leading contenders for SIN-3-mediated regulation of aging. This study looks at the concerted role of SIN-3 and DAF-16 proteins in lifespan regulation. Phenotypic profiling for the mutants of these genes shows the functional accord between these 2 proteins with similar functions in stress response and vital biological processes. However, there were no significant physical interactions when checked for protein-protein interaction between SIN-3 and DAF-16 proteins. C. elegans genomics and transcriptomics data also indicated the possibilities of concerted gene regulation. This genetic regulation is more likely related to SIN-3 dominance on DAF-16 function. Overall, SIN-3 and DAF-16 proteins have strong functional interactions that ensure healthy aging. The influence of SIN-3 on DAF-16-mediated stress response is one of their convergence points in longevity regulation.

XOL-1 regulates developmental timing by modulating the H3K9 landscape in C. elegans early embryos

Sex determination in the nematode C. elegans is controlled by the master regulator XOL-1 during embryogenesis. Expression of xol-1 is dependent on the ratio of X chromosomes and autosomes, which differs between XX hermaphrodites and XO males. In males, xol-1 is highly expressed and in hermaphrodites, xol-1 is expressed at very low levels. XOL-1 activity is known to be critical for the proper development of C. elegans males, but its low expression was considered to be of minimal importance in the development of hermaphrodite embryos. Our study reveals that XOL-1 plays an important role as a regulator of developmental timing during hermaphrodite embryogenesis. Using a combination of imaging and bioinformatics techniques, we found that hermaphrodite embryos have an accelerated rate of cell division, as well as a more developmentally advanced transcriptional program when xol-1 is lost. Further analyses reveal that XOL-1 is responsible for regulating the timing of initiation of dosage compensation on the X chromosomes, and the appropriate expression of sex-biased transcriptional programs in hermaphrodites. We found that xol-1 mutant embryos overexpress the H3K9 methyltransferase MET-2 and have an altered H3K9me landscape. Some of these effects of the loss of xol-1 gene were reversed by the loss of met-2. These findings demonstrate that XOL-1 plays an important role as a developmental regulator in embryos of both sexes, and that MET-2 acts as a downstream effector of XOL-1 activity in hermaphrodites.

Epidermal cell fusion promotes the transition from an embryonic to a larval transcriptome in C. elegans

Cell fusion is a fundamental process in the development of multicellular organisms, yet its impact on gene regulation, particularly during crucial developmental stages, remains poorly understood. The Caenorhabditis elegans epidermis comprises 8-10 syncytial cells, with the largest integrating 139 individual nuclei through cell-cell fusion governed by the fusogenic protein EFF-1. To explore the effects of cell fusion on developmental progression and associated gene expression changes, we conducted transcriptomic analyses of eff-1 fusion-deficient mutants. Our RNAseq findings showed widespread transcriptomic changes that were enriched for epidermal genes and key molecular pathways involved in epidermal function during larval development. Subsequent single-molecule fluorescence in situ hybridization validated the altered expression of mRNA transcripts, confirming quantifiable changes in gene expression in the absence of embryonic epidermal fusion. These results underscore the significance of cell-cell fusion in shaping transcriptional programs during development and raise questions regarding the precise identities and specialized functions of different subclasses of nuclei within developing syncytial cells and tissues.

Sex-specific developmental gene expression atlas unveils dimorphic gene networks in C. elegans

Sex-specific traits and behaviors emerge during development by the acquisition of unique properties in the nervous system of each sex. However, the genetic events responsible for introducing these sex-specific features remain poorly understood. In this study, we create a comprehensive gene expression atlas of pure populations of hermaphrodites and males of the nematode Caenorhabditis elegans across development. We discover numerous differentially expressed genes, including neuronal gene families like transcription factors, neuropeptides, and G protein-coupled receptors. We identify INS-39, an insulin-like peptide, as a prominent male-biased gene expressed specifically in ciliated sensory neurons. We show that INS-39 serves as an early-stage male marker, facilitating the effective isolation of males in high-throughput experiments. Through complex and sex-specific regulation, ins-39 plays pleiotropic sexually dimorphic roles in various behaviors, while also playing a shared, dimorphic role in early life stress. This study offers a comparative sexual and developmental gene expression database for C. elegans. Furthermore, it highlights conserved genes that may underlie the sexually dimorphic manifestation of different human diseases.

The Caenorhabditis RNA-seq Browser: a web-based application for on-demand analysis of publicly available Caenorhabditis spp. bulk RNA-sequencing data

The Caenorhabditis RNA-seq Browser is an open-source Shiny web app that enables on-demand visualization and quantification of bulk RNA-sequencing data for five Caenorhabditis species: C. elegans , C. briggsae , C. brenneri , C. japonica , and C. remanei . The app is designed to allow researchers without previous coding experience to interactively explore publicly available Caenorhabditis RNA-sequencing data. Key app features include the ability to plot gene expression across life stages for user-specified gene sets, and modules for performing differential gene expression analyses. The Caenorhabditis RNA-seq Browser can be accessed online via shinyapps.io or can be installed locally in R from a GitHub repository.

RNA-seq analysis reveals changes in mRNA expression during development in Daphnia mitsukuri

Temporal transcriptional variation is a major contributor to functional evolution and the developmental process. Parthenogenetic water fleas of the genus Daphnia (Cladocera) provide an ideal model to characterize gene expression patterns across distinct developmental stages. Herein, we report RNA-seq data for female Daphnia mitsukuri at three developmental stages: the embryo, juvenile (three timepoints) and adult. Comparisons of gene expression patterns among these three developmental stages and weighted gene co-expression network analysis based on expression data across developmental stages identified sets of genes underpinning each of the developmental stages of D. mitsukuri. Specifically, highly expressed genes (HEGs) at the embryonic developmental stage were associated with cell proliferation, ensuring the necessary foundation for subsequent development; HEGs at the juvenile stages were associated with chemosensory perception, visual perception and neurotransmission, allowing individuals to enhance detection of potential environmental risks; HEGs at the adult stage were associated with antioxidative defensive systems, enabling adults to mount an efficient response to perceived environmental risks. Additionally, we found a significant overlap between expanded gene families of Daphnia species and HEGs at the juvenile stages, and these genes were associated with visual perception and neurotransmission. Our work provides a resource of developmental transcriptomes, and comparative analyses that characterize gene expression dynamics throughout development of Daphnia.

Characterization of the distribution and dynamics of chromatin states in the C. elegans germline reveals substantial H3K4me3 remodeling during oogenesis

Chromatin organization in the C. elegans germline is tightly regulated and critical for germ cell differentiation. Although certain germline epigenetic regulatory mechanisms have been identified, how they influence chromatin structure and ultimately gene expression remains unclear, in part because most genomic studies have focused on data collected from intact worms comprising both somatic and germline tissues. We therefore analyzed histone modification and chromatin accessibility data from isolated germ nuclei representing undifferentiated proliferating and meiosis I populations to define chromatin states. We correlated these states with overall transcript abundance, spatiotemporal expression patterns, and the function of small RNA pathways. Because the essential role of the germline is to transmit genetic information and establish gene expression in the early embryo, we compared epigenetic and transcriptomic profiles from undifferentiated germ cells to those of embryos to define the epigenetic changes during this developmental transition. The active histone modification H3K4me3 shows particularly dynamic remodeling as germ cells differentiate into oocytes, which suggests a mechanism for establishing early transcription of essential genes during zygotic genome activation. This analysis highlights the dynamism of the chromatin landscape across developmental transitions and provides a resource for future investigation into epigenetic regulatory mechanisms in germ cells.

Uncovering the essential roles of glutamate carboxypeptidase 2 orthologs in Caenorhabditis elegans

Human glutamate carboxypeptidase 2 (GCP2) from the M28B metalloprotease group is an important target for therapy in neurological disorders and an established tumor marker. However, its physiological functions remain unclear. To better understand general roles, we used the model organism Caenorhabditis elegans to genetically manipulate its three existing orthologous genes and evaluate the impact on worm physiology. The results of gene knockout studies showed that C. elegans GCP2 orthologs affect the pharyngeal physiology, reproduction, and structural integrity of the organism. Promoter-driven GFP expression revealed distinct localization for each of the three gene paralogs, with gcp-2.1 being most abundant in muscles, intestine, and pharyngeal interneurons, gcp-2.2 restricted to the phasmid neurons, and gcp-2.3 located in the excretory cell. The present study provides new insight into the unique phenotypic effects of GCP2 gene knockouts in C. elegans, and the specific tissue localizations. We believe that elucidation of particular roles in a non-mammalian organism can help to explain important questions linked to physiology of this protease group and in extension to human GCP2 involvement in pathophysiological processes.

Neuroscience Needs Network Science

The brain is a complex system comprising a myriad of interacting neurons, posing significant challenges in understanding its structure, function, and dynamics. Network science has emerged as a powerful tool for studying such interconnected systems, offering a framework for integrating multiscale data and complexity. To date, network methods have significantly advanced functional imaging studies of the human brain and have facilitated the development of control theory-based applications for directing brain activity. Here, we discuss emerging frontiers for network neuroscience in the brain atlas era, addressing the challenges and opportunities in integrating multiple data streams for understanding the neural transitions from development to healthy function to disease. We underscore the importance of fostering interdisciplinary opportunities through workshops, conferences, and funding initiatives, such as supporting students and postdoctoral fellows with interests in both disciplines. By bringing together the network science and neuroscience communities, we can develop novel network-based methods tailored to neural circuits, paving the way toward a deeper understanding of the brain and its functions, as well as offering new challenges for network science.

Neuroscience needs Network Science

The brain is a complex system comprising a myriad of interacting elements, posing significant challenges in understanding its structure, function, and dynamics. Network science has emerged as a powerful tool for studying such intricate systems, offering a framework for integrating multiscale data and complexity. Here, we discuss the application of network science in the study of the brain, addressing topics such as network models and metrics, the connectome, and the role of dynamics in neural networks. We explore the challenges and opportunities in integrating multiple data streams for understanding the neural transitions from development to healthy function to disease, and discuss the potential for collaboration between network science and neuroscience communities. We underscore the importance of fostering interdisciplinary opportunities through funding initiatives, workshops, and conferences, as well as supporting students and postdoctoral fellows with interests in both disciplines. By uniting the network science and neuroscience communities, we can develop novel network-based methods tailored to neural circuits, paving the way towards a deeper understanding of the brain and its functions.

Developmentally programmed histone H3 expression regulates cellular plasticity at the parental-to-early embryo transition

During metazoan development, the marked change in developmental potential from the parental germline to the embryo raises an important question regarding how the next life cycle is reset. As the basic unit of chromatin, histones are essential for regulating chromatin structure and function and, accordingly, transcription. However, the genome-wide dynamics of the canonical, replication-coupled (RC) histones during gametogenesis and embryogenesis remain unknown. In this study, we use CRISPR-Cas9-mediated gene editing in Caenorhabditis elegans to investigate the expression pattern and role of individual RC histone H3 genes and compare them to the histone variant, H3.3. We report a tightly regulated epigenome landscape change from the germline to embryos that are regulated through differential expression of distinct histone gene clusters. Together, this study reveals that a change from a H3.3- to H3-enriched epigenome during embryogenesis restricts developmental plasticity and uncovers distinct roles for individual H3 genes in regulating germline chromatin.

N-glycan antennal modifications are altered in Caenorhabditis elegans lacking the HEX-4 N-acetylgalactosamine-specific hexosaminidase

Simple organisms are often considered to have simple glycomes, but plentiful paucimannosidic and oligomannosidic glycans overshadow the less abundant N-glycans with highly variable core and antennal modifications; Caenorhabditis elegans is no exception. By use of optimized fractionation and assessing wildtype in comparison to mutant strains lacking either the HEX-4 or HEX-5 β-N-acetylgalactosaminidases, we conclude that the model nematode has a total N-glycomic potential of 300 verified isomers. Three pools of glycans were analyzed for each strain: either PNGase F released and eluted from a reversed-phase C18 resin with either water or 15% methanol or PNGase Ar released. While the water-eluted fractions were dominated by typical paucimannosidic and oligomannosidic glycans and the PNGase Ar-released pools by glycans with various core modifications, the methanol-eluted fractions contained a huge range of phosphorylcholine-modified structures with up to three antennae, sometimes with four N-acetylhexosamine residues in series. There were no major differences between the C. elegans wildtype and hex-5 mutant strains, but the hex-4 mutant strains displayed altered sets of methanol-eluted and PNGase Ar-released pools. In keeping with the specificity of HEX-4, there were more glycans capped with N-acetylgalactosamine in the hex-4 mutants, as compared with isomeric chito-oligomer motifs in the wildtype. Considering that fluorescence microscopy showed that a HEX-4::enhanced GFP fusion protein colocalizes with a Golgi tracker, we conclude that HEX-4 plays a significant role in late-stage Golgi processing of N-glycans in C. elegans. Furthermore, finding more "parasite-like" structures in the model worm may facilitate discovery of glycan-processing enzymes occurring in other nematodes.

wrk-1 and rig-5 control pioneer and follower axon navigation in the ventral nerve cord of Caenorhabditis elegans in a nid-1 mutant background

During nervous system development, neurons send out axons, which must navigate large distances to reach synaptic targets. Axons grow out sequentially. The early outgrowing axons, pioneers, must integrate information from various guidance cues in their environment to determine the correct direction of outgrowth. Later outgrowing follower axons can at least in part navigate by adhering to pioneer axons. In Caenorhabditis elegans, the right side of the largest longitudinal axon tract, the ventral nerve cord, is pioneered by the AVG axon. How the AVG axon navigates is only partially understood. In this study, we describe the role of two members of the IgCAM family, wrk-1 and rig-5, in AVG axon navigation. While wrk-1 and rig-5 single mutants do not show AVG navigation defects, both mutants have highly penetrant pioneer and follower navigation defects in a nid-1 mutant background. Both mutations increase the fraction of follower axons following the misguided pioneer axon. We found that wrk-1 and rig-5 act in different genetic pathways, suggesting that we identified two pioneer-independent guidance pathways used by follower axons. We assessed general locomotion, mechanosensory responsiveness, and habituation to determine whether axonal navigation defects impact nervous system function. In rig-5 nid-1 double mutants, we found no significant defects in free movement behavior; however, a subpopulation of animals shows minor changes in response duration habituation after mechanosensory stimulation. These results suggest that guidance defects of axons in the motor circuit do not necessarily lead to major movement or behavioral defects but impact more complex behavioral modulation.

Increasing Complexity of the N-Glycome During Caenorhabditis Development

Caenorhabditis elegans is a frequently employed genetic model organism and has been the object of a wide range of developmental, genetic, proteomic, and glycomic studies. Here, using an off-line MALDI-TOF-MS approach, we have analyzed the N-glycans of mixed embryos and liquid- or plate-grown L4 larvae. Of the over 200 different annotatable N-glycan structures, variations between the stages as well as the mode of cultivation were observed. While the embryonal N-glycome appears less complicated overall, the liquid- and plate-grown larvae differ especially in terms of methylation of bisecting fucose, α-galactosylation of mannose, and di-β-galactosylation of core α1,6-fucose. Furthermore, we analyzed the O-glycans by LC-electrospray ionization-MS following β-elimination; especially the embryonal O-glycomes included a set of phosphorylcholine-modified structures, previously not shown to exist in nematodes. However, the set of glycan structures cannot be clearly correlated with levels of glycosyltransferase transcripts in developmental RNA-Seq datasets, but there is an indication for coordinated expression of clusters of potential glycosylation-relevant genes. Thus, there are still questions to be answered in terms of how and why a simple nematode synthesizes such a diverse glycome.

Nematode gene annotation by machine-learning-assisted proteotranscriptomics enables proteome-wide evolutionary analysis

Nematodes encompass more than 24,000 described species, which were discovered in almost every ecological habitat, and make up >80% of metazoan taxonomic diversity in soils. The last common ancestor of nematodes is believed to date back to ∼650-750 million years, generating a large and phylogenetically diverse group to be explored. However, for most species high-quality gene annotations are incomprehensive or missing. Combining short-read RNA sequencing with mass spectrometry-based proteomics and machine-learning quality control in an approach called proteotranscriptomics, we improve gene annotations for nine genome-sequenced nematode species and provide new gene annotations for three additional species without genome assemblies. Emphasizing the sensitivity of our methodology, we provide evidence for two hitherto undescribed genes in the model organism Caenorhabditis elegans Extensive phylogenetic systems analysis using this comprehensive proteome annotation provides new insights into evolutionary processes of this metazoan group.

A lipid transfer protein ensures nematode cuticular impermeability

The cuticle of C. elegans is impermeable to chemicals, toxins, and pathogens. However, increased permeability is a desirable phenotype because it facilitates chemical uptake. Surface lipids contribute to the permeability barrier. Here, we identify the lipid transfer protein GMAP-1 as a critical element setting the permeability of the C. elegans cuticle. A gmap-1 deletion mutant increases cuticular permeability to sodium azide, levamisole, Hoechst, and DiI. Expressing GMAP-1 in the hypodermis or transiently in the adults is sufficient to rescue this gmap-1 permeability phenotype. GMAP-1 protein is secreted from the hypodermis to the aqueous fluid filling the space between collagen fibers of the cuticle. In vitro, GMAP-1 protein binds phosphatidylserine and phosphatidylcholine while in vivo, GMAP-1 sets the surface lipid composition and organization. Altogether, our results suggest GMAP-1 secreted by hypodermis shuttles lipids to the surface to form the permeability barrier of C. elegans.

•

GMAP-1 is secreted by the hypodermis toward the cuticle of Caenorhabditis elegans

•

GMAP-1 binds and shuttle phosphoglycerides

•

GMAP-1 sets the lipid composition of the cuticle

•

While healthy, gmap-1 mutant displays high cuticular permeability

AE-TPGG: a novel autoencoder-based approach for single-cell RNA-seq data imputation and dimensionality reduction

Single-cell RNA sequencing (scRNA-seq) technology has become an effective tool for high-throughout transcriptomic study, which circumvents the averaging artifacts corresponding to bulk RNA-seq technology, yielding new perspectives on the cellular diversity of potential superficially homogeneous populations. Although various sequencing techniques have decreased the amplification bias and improved capture efficiency caused by the low amount of starting material, the technical noise and biological variation are inevitably introduced into experimental process, resulting in high dropout events, which greatly hinder the downstream analysis. Considering the bimodal expression pattern and the right-skewed characteristic existed in normalized scRNA-seq data, we propose a customized autoencoder based on a two-part-generalized-gamma distribution (AE-TPGG) for scRNA-seq data analysis, which takes mixed discrete-continuous random variables of scRNA-seq data into account using a two-part model and utilizes the generalized gamma (GG) distribution, for fitting the positive and right-skewed continuous data. The adopted autoencoder enables AE-TPGG to captures the inherent relationship between genes. In addition to the ability of achieving low-dimensional representation, the AE-TPGG model also provides a denoised imputation according to statistical characteristic of gene expression. Results on real datasets demonstrate that our proposed model is competitive to current imputation methods and ameliorates a diverse set of typical scRNA-seq data analyses.

ELECTRONIC SUPPLEMENTARY MATERIAL

Supplementary material is available in the online version of this article at 10.1007/s11704-022-2011-y.

The unfolded protein response reverses the effects of glucose on lifespan in chemically-sterilized C. elegans

Metabolic diseases often share common traits, including accumulation of unfolded proteins in the endoplasmic reticulum (ER). Upon ER stress, the unfolded protein response (UPR) is activated to limit cellular damage which weakens with age. Here, we show that Caenorhabditis elegans fed a bacterial diet supplemented high glucose at day 5 of adulthood (HGD-5) extends their lifespan, whereas exposed at day 1 (HGD-1) experience shortened longevity. We observed a metabolic shift only in HGD-1, while glucose and infertility synergistically prolonged the lifespan of HGD-5, independently of DAF-16. Notably, we identified that UPR stress sensors ATF-6 and PEK-1 contributed to the longevity of HGD-5 worms, while ire-1 ablation drastically increased HGD-1 lifespan. Together, we postulate that HGD activates the otherwise quiescent UPR in aged worms to overcome ageing-related stress and restore ER homeostasis. In contrast, young animals subjected to HGD provokes unresolved ER stress, conversely leading to a detrimental stress response.

Dynamic Transcriptional Landscape of Grass Carp (Ctenopharyngodon idella) Reveals Key Transcriptional Features Involved in Fish Development

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

Defect-buffering cellular plasticity increases robustness of metazoan embryogenesis

Developmental processes are intrinsically robust so as to preserve a normal-like state in response to genetic and environmental fluctuations. However, the robustness and potential phenotypic plasticity of individual developing cells under genetic perturbations remain to be systematically evaluated. Using large-scale gene perturbation, live imaging, lineage tracing, and single-cell phenomics, we quantified the phenotypic landscape of C. elegans embryogenesis in >2,000 embryos following individual knockdown of over 750 conserved genes. We observed that cellular genetic systems are not sufficiently robust to single-gene perturbations across all cells; rather, gene knockdowns frequently induced cellular defects. Dynamic phenotypic analyses revealed many cellular defects to be transient, with cells exhibiting phenotypic plasticity that serves to alleviate, correct, and accommodate the defects. Moreover, potential developmentally related cell modules may buffer the phenotypic effects of individual cell position changes. Our findings reveal non-negligible contributions of cellular plasticity and multicellularity as compensatory strategies to increase developmental robustness.

Microbial byproducts determine reproductive fitness of free-living and parasitic nematodes

Trichuris nematodes reproduce within the microbiota-rich mammalian intestine and lay thousands of eggs daily, facilitating their sustained presence in the environment and hampering eradication efforts. Here, we show that bacterial byproducts facilitate the reproductive development of nematodes. First, we employed a pipeline using the well-characterized, free-living nematode C. elegans to identify microbial factors with conserved roles in nematode reproduction. A screen for E. coli mutants that impair C. elegans fertility identified genes in fatty acid biosynthesis and ethanolamine utilization pathways, including fabH and eutN. Additionally, Trichuris muris eggs displayed defective hatching in the presence of fabH- or eutN-deficient E. coli due to reduced arginine or elevated aldehydes, respectively. T. muris reared in gnotobiotic mice colonized with these E. coli mutants displayed morphological defects and failed to lay viable eggs. These findings indicate that microbial byproducts mediate evolutionarily conserved transkingdom interactions that impact the reproductive fitness of distantly related nematodes.

Nematode chromosomes

The nematode Caenorhabditis elegans has shed light on many aspects of eukaryotic biology, including genetics, development, cell biology, and genomics. A major factor in the success of C. elegans as a model organism has been the availability, since the late 1990s, of an essentially gap-free and well-annotated nuclear genome sequence, divided among 6 chromosomes. In this review, we discuss the structure, function, and biology of C. elegans chromosomes and then provide a general perspective on chromosome biology in other diverse nematode species. We highlight malleable chromosome features including centromeres, telomeres, and repetitive elements, as well as the remarkable process of programmed DNA elimination (historically described as chromatin diminution) that induces loss of portions of the genome in somatic cells of a handful of nematode species. An exciting future prospect is that nematode species may enable experimental approaches to study chromosome features and to test models of chromosome evolution. In the long term, fundamental insights regarding how speciation is integrated with chromosome biology may be revealed.

Intergenerational adaptations to stress are evolutionarily conserved, stress-specific, and have deleterious trade-offs

Despite reports of parental exposure to stress promoting physiological adaptations in progeny in diverse organisms, there remains considerable debate over the significance and evolutionary conservation of such multigenerational effects. Here, we investigate four independent models of intergenerational adaptations to stress in Caenorhabditis elegans - bacterial infection, eukaryotic infection, osmotic stress, and nutrient stress - across multiple species. We found that all four intergenerational physiological adaptations are conserved in at least one other species, that they are stress -specific, and that they have deleterious tradeoffs in mismatched environments. By profiling the effects of parental bacterial infection and osmotic stress exposure on progeny gene expression across species, we established a core set of 587 genes that exhibited a greater than twofold intergenerational change in expression in response to stress in C. elegans and at least one other species, as well as a set of 37 highly conserved genes that exhibited a greater than twofold intergenerational change in expression in all four species tested. Furthermore, we provide evidence suggesting that presumed adaptive and deleterious intergenerational effects are molecularly related at the gene expression level. Lastly, we found that none of the effects we detected of these stresses on C. elegans F1 progeny gene expression persisted transgenerationally three generations after stress exposure. We conclude that intergenerational responses to stress play a substantial and evolutionarily conserved role in regulating animal physiology and that the vast majority of the effects of parental stress on progeny gene expression are reversible and not maintained transgenerationally.

Identification of essential genes in Caenorhabditis elegans through whole genome sequencing of legacy mutant collections

It has been estimated that 15%–30% of the ∼20,000 genes in C. elegans are essential, yet many of these genes remain to be identified or characterized. With the goal of identifying unknown essential genes, we performed whole-genome sequencing on complementation pairs from legacy collections of maternal-effect lethal and sterile mutants. This approach uncovered maternal genes required for embryonic development and genes with apparent sperm-specific functions. In total, 58 putative essential genes were identified on chromosomes III–V, of which 52 genes are represented by novel alleles in this collection. Of these 52 genes, 19 (40 alleles) were selected for further functional characterization. The terminal phenotypes of embryos were examined, revealing defects in cell division, morphogenesis, and osmotic integrity of the eggshell. Mating assays with wild-type males revealed previously unknown male-expressed genes required for fertilization and embryonic development. The result of this study is a catalog of mutant alleles in essential genes that will serve as a resource to guide further study toward a more complete understanding of this important model organism. As many genes and developmental pathways in C. elegans are conserved and essential genes are often linked to human disease, uncovering the function of these genes may also provide insight to further our understanding of human biology.

Transcriptome analysis provides genome annotation and expression profiles in the central nervous system of Lymnaea stagnalis at different ages

BACKGROUND

The pond snail, Lymnaea stagnalis (L. stagnalis), has served as a valuable model organism for neurobiology studies due to its simple and easily accessible central nervous system (CNS). L. stagnalis has been widely used to study neuronal networks and recently gained popularity for study of aging and neurodegenerative diseases. However, previous transcriptome studies of L. stagnalis CNS have been exclusively carried out on adult L. stagnalis only. As part of our ongoing effort studying L. stagnalis neuronal growth and connectivity at various developmental stages, we provide the first age-specific transcriptome analysis and gene annotation of young (3 months), adult (6 months), and old (18 months) L. stagnalis CNS.

RESULTS

Using the above three age cohorts, our study generated 55-69 millions of 150 bp paired-end RNA sequencing reads using the Illumina NovaSeq 6000 platform. Of these reads, ~ 74% were successfully mapped to the reference genome of L. stagnalis. Our reference-based transcriptome assembly predicted 42,478 gene loci, of which 37,661 genes encode coding sequences (CDS) of at least 100 codons. In addition, we provide gene annotations using Blast2GO and functional annotations using Pfam for ~ 95% of these sequences, contributing to the largest number of annotated genes in L. stagnalis CNS so far. Moreover, among 242 previously cloned L. stagnalis genes, we were able to match ~ 87% of them in our transcriptome assembly, indicating a high percentage of gene coverage. The expressional differences for innexins, FMRFamide, and molluscan insulin peptide genes were validated by real-time qPCR. Lastly, our transcriptomic analyses revealed distinct, age-specific gene clusters, differentially expressed genes, and enriched pathways in young, adult, and old CNS. More specifically, our data show significant changes in expression of critical genes involved in transcription factors, metabolisms (e.g. cytochrome P450), extracellular matrix constituent, and signaling receptor and transduction (e.g. receptors for acetylcholine, N-Methyl-D-aspartic acid, and serotonin), as well as stress- and disease-related genes in young compared to either adult or old snails.

CONCLUSIONS

Together, these datasets are the largest and most updated L. stagnalis CNS transcriptomes, which will serve as a resource for future molecular studies and functional annotation of transcripts and genes in L. stagnalis.

A 4D single-cell protein atlas of transcription factors delineates spatiotemporal patterning during embryogenesis

Complex biological processes such as embryogenesis require precise coordination of cell differentiation programs across both space and time. Using protein-fusion fluorescent reporters and four-dimensional live imaging, we present a protein expression atlas of transcription factors (TFs) mapped onto developmental cell lineages during Caenorhabditis elegans embryogenesis, at single-cell resolution. This atlas reveals a spatiotemporal combinatorial code of TF expression, and a cascade of lineage-specific, tissue-specific and time-specific TFs that specify developmental states. The atlas uncovers regulators of embryogenesis, including an unexpected role of a skin specifier in neurogenesis and the critical function of an uncharacterized TF in convergent muscle differentiation. At the systems level, the atlas provides an opportunity to model cell state-fate relationships, revealing a lineage-dependent state diversity within functionally related cells and a winding trajectory of developmental state progression. Collectively, this single-cell protein atlas represents a valuable resource for elucidating metazoan embryogenesis at the molecular and systems levels.

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.

Comprehensive characterization of tissue-specific chromatin accessibility in L2 Caenorhabditis elegans nematodes

Recently developed single-cell technologies allow researchers to characterize cell states at ever greater resolution and scale. Caenorhabditis elegans is a particularly tractable system for studying development, and recent single-cell RNA-seq studies characterized the gene expression patterns for nearly every cell type in the embryo and at the second larval stage (L2). Gene expression patterns give insight about gene function and into the biochemical state of different cell types; recent advances in other single-cell genomics technologies can now also characterize the regulatory context of the genome that gives rise to these gene expression levels at a single-cell resolution. To explore the regulatory DNA of individual cell types in C. elegans, we collected single-cell chromatin accessibility data using the sci-ATAC-seq assay in L2 larvae to match the available single-cell RNA-seq data set. By using a novel implementation of the latent Dirichlet allocation algorithm, we identify 37 clusters of cells that correspond to different cell types in the worm, providing new maps of putative cell type-specific gene regulatory sites, with promise for better understanding of cellular differentiation and gene regulation.

High-Temporal-Resolution smFISH Method for Gene Expression Studies in Caenorhabditis elegans Embryos

Recent development in fluorescence-based molecular tools has contributed significantly to developmental studies, including embryogenesis. Many of these tools rely on multiple steps of sample manipulation, so obtaining large sample sizes presents a major challenge as it can be labor-intensive and time-consuming. However, large sample sizes are required to uncover critical aspects of embryogenesis, for example, subtle phenotypic differences or gene expression dynamics. This problem is particularly relevant for single-molecule fluorescence in situ hybridization (smFISH) studies in Caenorhabditis elegans embryogenesis. Microfluidics can help address this issue by allowing a large number of samples and parallelization of experiments. However, performing efficient reagent exchange on chip for large numbers of embryos remains a bottleneck. Here, we present a microfluidic pipeline for large-scale smFISH imaging of C. elegans embryos with minimized labor. We designed embryo traps and engineered a protocol allowing for efficient chemical exchange for hundreds of C. elegans embryos simultaneously. Furthermore, the device design and small footprint optimize imaging throughput by facilitating spatial registration and enabling minimal user input. We conducted the smFISH protocol on chip and demonstrated that image quality is preserved. With one device replacing the equivalent of 10 glass slides of embryos mounted manually, our microfluidic approach greatly increases throughput. Finally, to highlight the capability of our platform to perform longitudinal studies with high temporal resolution, we conducted a temporal analysis of par-1 gene expression in early C. elegans embryos. The method demonstrated here paves the way for systematic high-temporal-resolution studies that will benefit large-scale RNAi and drug screens and in systems beyond C. elegans embryos.

PhenoMIP: High-Throughput Phenotyping of Diverse Caenorhabditis elegans Populations via Molecular Inversion Probes

Whether generated within a lab setting or isolated from the wild, variant alleles continue to be an important resource for decoding gene function in model organisms such as Caenorhabditis elegans. With advances in massively parallel sequencing, multiple whole-genome sequenced (WGS) strain collections are now available to the research community. The Million Mutation Project (MMP) for instance, analyzed 2007 N2-derived, mutagenized strains. Individually, each strain averages ∼400 single nucleotide variants amounting to ∼80 protein-coding variants. The effects of these variants, however, remain largely uncharacterized and querying the breadth of these strains for phenotypic changes requires a method amenable to rapid and sensitive high-throughput analysis. Here we present a pooled competitive fitness approach to quantitatively phenotype subpopulations of sequenced collections via molecular inversion probes (PhenoMIP). We phenotyped the relative fitness of 217 mutant strains on multiple food sources and classified these into five categories. We also demonstrate on a subset of these strains, that their fitness defects can be genetically mapped. Overall, our results suggest that approximately 80% of MMP mutant strains may have a decreased fitness relative to the lab reference, N2. The costs of generating this form of analysis through WGS methods would be prohibitive while PhenoMIP analysis in this manner is accomplished at less than one-tenth of projected WGS costs. We propose methods for applying PhenoMIP to a broad range of population selection experiments in a cost-efficient manner that would be useful to the community at large.

Two MicroRNAs Are Sufficient for Embryonic Patterning in C. elegans

MicroRNAs (miRNAs) are a class of post-transcriptional repressors with diverse roles in animal development and physiology [1]. The Microprocessor complex, composed of Drosha and Pasha/DGCR8, is necessary for the biogenesis of all canonical miRNAs and essential for the early stages of animal embryogenesis [2, 3, 4, 5, 6, 7, 8]. However, the cause for this requirement is largely unknown. Animals often express hundreds of miRNAs, and it remains unclear whether the Microprocessor is required to produce one or few essential miRNAs or many individually non-essential miRNAs. Additionally, both Drosha and Pasha/DGCR8 bind and cleave a variety of non-miRNA substrates [9, 10, 11, 12, 13, 14, 15], and it is unknown whether these activities account for the Microprocessor’s essential requirement. To distinguish between these possibilities, we developed a system in C. elegans to stringently deplete embryos of Microprocessor activity. Using a combination of auxin-inducible degradation (AID) and RNA interference (RNAi), we achieved Drosha and Pasha/DGCR8 depletion starting in the maternal germline, resulting in Microprocessor and miRNA-depleted embryos, which fail to undergo morphogenesis or form organs. Using a Microprocessor-bypass strategy, we show that this early embryonic arrest is rescued by the addition of just two miRNAs, one miR-35 and one miR-51 family member, resulting in morphologically normal larvae. Thus, just two out of ∼150 canonical miRNAs are sufficient for morphogenesis and organogenesis, and the processing of these miRNAs accounts for the essential requirement for Drosha and Pasha/DGCR8 during the early stages of C. elegans embryonic development.

Video Abstract

•

Depletion of Drosha and Pasha results in embryos that fail to undergo morphogenesis

•

The mirtron pathway enables expression of miRNAs in the absence of Drosha and Pasha

•

Two miRNAs are sufficient to rescue embryogenesis in the absence of Drosha and Pasha

•

miR-35 and miR-51 play an unexplored, likely conserved role in animal development

Differential Expression Gene Explorer (DrEdGE): a tool for generating interactive online visualizations of gene expression datasets

SUMMARY

Differential Expression Gene Explorer (DrEdGE) is a web-based tool that guides genomicists through easily creating interactive online data visualizations, which colleagues can query according to their own conditions to discover genes, samples or patterns of interest. We demonstrate DrEdGE's features with three example websites generated from publicly available datasets-human neuronal tissue, mouse embryonic tissue and Caenorhabditis elegans whole embryos. DrEdGE increases the utility of large genomics datasets by removing technical obstacles to independent exploration.

AVAILABILITY AND IMPLEMENTATION

Freely available at http://dredge.bio.unc.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Ancestral function of Inhibitors-of-kappaB regulates Caenorhabditis elegans development

Mammalian IκB proteins (IκBs) exert their main function as negative regulators of NF-κB, a central signaling pathway controlling immunity and inflammation. An alternative chromatin role for IκBs has been shown to affect stemness and cell differentiation. However, the involvement of NF-κB in this function has not been excluded. NFKI-1 and IKB-1 are IκB homologs in Caenorhabditis elegans, which lacks NF-κB nuclear effectors. We found that nfki-1 and ikb-1 mutants display developmental defects that phenocopy mutations in Polycomb and UTX-1 histone demethylase, suggesting a role for C. elegans IκBs in chromatin regulation. Further supporting this possibility (1) we detected NFKI-1 in the nucleus of cells; (2) NFKI-1 and IKB-1 bind to histones and Polycomb proteins, (3) and associate with chromatin in vivo, and (4) mutations in nfki-1 and ikb-1 alter chromatin marks. Based on these results, we propose that ancestral IκB inhibitors modulate Polycomb activity at specific gene subsets with an impact on development.

Combinatorial Action of Temporally Segregated Transcription Factors

Combinatorial action of transcription factors (TFs) with partially overlapping expression is a widespread strategy to generate novel gene-expression patterns and, thus, cellular diversity. Known mechanisms underlying combinatorial activity require co-expression of TFs within the same cell. Here, we describe the mechanism by which two TFs that are never co-expressed generate a new, intersectional expression pattern in C. elegans embryos: lineage-specific priming of a gene by a transiently expressed TF generates a unique intersection with a second TF acting on the same gene four cell divisions later; the second TF is expressed in multiple cells but only activates transcription in those where priming occurred. Early induction of active transcription is necessary and sufficient to establish a competent state, maintained by broadly expressed regulators in the absence of the initial trigger. We uncover additional cells diversified through this mechanism. Our findings define a mechanism for combinatorial TF activity with important implications for generation of cell-type diversity.

•

Lineage-specific priming enables asymmetric gene expression in L/R neuron pairs

•

Transient, lineage-specific TFs prime a locus for later activation by a bilateral TF

•

An early active transcriptional state is necessary and sufficient for priming

•

Maintenance of asymmetric primed state occurs in a symmetric regulatory environment

Speciation and the developmental alarm clock

New species arise as the genomes of populations diverge. The developmental 'alarm clock' of speciation sounds off when sufficient divergence in genetic control of development leads hybrid individuals to infertility or inviability, the world awoken to the dawn of new species with intrinsic post-zygotic reproductive isolation. Some developmental stages will be more prone to hybrid dysfunction due to how molecular evolution interacts with the ontogenetic timing of gene expression. Considering the ontogeny of hybrid incompatibilities provides a profitable connection between 'evo-devo' and speciation genetics to better link macroevolutionary pattern, microevolutionary process, and molecular mechanisms. Here, we explore speciation alongside development, emphasizing their mutual dependence on genetic network features, fitness landscapes, and developmental system drift. We assess models for how ontogenetic timing of reproductive isolation can be predictable. Experiments and theory within this synthetic perspective can help identify new rules of speciation as well as rules in the molecular evolution of development.

Direct full-length RNA sequencing reveals unexpected transcriptome complexity during Caenorhabditis elegans development

Massively parallel sequencing of the polyadenylated RNAs has played a key role in delineating transcriptome complexity, including alternative use of an exon, promoter, 5′ or 3′ splice site or polyadenylation site, and RNA modification. However, reads derived from the current RNA-seq technologies are usually short and deprived of information on modification, compromising their potential in defining transcriptome complexity. Here, we applied a direct RNA sequencing method with ultralong reads using Oxford Nanopore Technologies to study the transcriptome complexity in Caenorhabditis elegans. We generated approximately six million reads using native poly(A)-tailed mRNAs from three developmental stages, with average read lengths ranging from 900 to 1100 nt. Around half of the reads represent full-length transcripts. To utilize the full-length transcripts in defining transcriptome complexity, we devised a method to classify the long reads as the same as existing transcripts or as a novel transcript using sequence mapping tracks rather than existing intron/exon structures, which allowed us to identify roughly 57,000 novel isoforms and recover at least 26,000 out of the 33,500 existing isoforms. The sets of genes with differential expression versus differential isoform usage over development are largely different, implying a fine-tuned regulation at isoform level. We also observed an unexpected increase in putative RNA modification in all bases in the coding region relative to the UTR, suggesting their possible roles in translation. The RNA reads and the method for read classification are expected to deliver new insights into RNA processing and modification and their underlying biology in the future.

A lineage-resolved molecular atlas of C. elegans embryogenesis at single-cell resolution

Caenorhabditis elegans is an animal with few cells but a wide diversity of cell types. In this study, we characterize the molecular basis for their specification by profiling the transcriptomes of 86,024 single embryonic cells. We identify 502 terminal and preterminal cell types, mapping most single-cell transcriptomes to their exact position in C. elegans' invariant lineage. Using these annotations, we find that (i) the correlation between a cell's lineage and its transcriptome increases from middle to late gastrulation, then falls substantially as cells in the nervous system and pharynx adopt their terminal fates; (ii) multilineage priming contributes to the differentiation of sister cells at dozens of lineage branches; and (iii) most distinct lineages that produce the same anatomical cell type converge to a homogenous transcriptomic state.

Evolutionary characteristics of intergenic transcribed regions indicate rare novel genes and widespread noisy transcription in the Poaceae

Extensive transcriptional activity occurring in intergenic regions of genomes has raised the question whether intergenic transcription represents the activity of novel genes or noisy expression. To address this, we evaluated cross-species and post-duplication sequence and expression conservation of intergenic transcribed regions (ITRs) in four Poaceae species. Among 43,301 ITRs across the four species, 34,460 (80%) are species-specific. ITRs found across species tend to be more divergent in expression and have more recent duplicates compared to annotated genes. To assess if ITRs are functional (under selection), machine learning models were established in Oryza sativa (rice) that could accurately distinguish between phenotype genes and pseudogenes (area under curve-receiver operating characteristic = 0.94). Based on the models, 584 (8%) and 4391 (61%) rice ITRs are classified as likely functional and nonfunctional with high confidence, respectively. ITRs with conserved expression and ancient retained duplicates, features that were not part of the model, are frequently classified as likely-functional, suggesting these characteristics could serve as pragmatic rules of thumb for identifying candidate sequences likely to be under selection. This study also provides a framework to identify novel genes using comparative transcriptomic data to improve genome annotation that is fundamental for connecting genotype to phenotype in crop and model systems.

H3.3K27M-induced chromatin changes drive ectopic replication through misregulation of the JNK pathway in C. elegans

Substitution of lysine 27 with methionine in histone H3.3 is a recently discovered driver mutation of pediatric high-grade gliomas. Mutant cells show decreased levels and altered distribution of H3K27 trimethylation (H3K27me3). How these chromatin changes are established genome-wide and lead to tumorigenesis remains unclear. Here we show that H3.3K27M-mediated alterations in H3K27me3 distribution result in ectopic DNA replication and cell cycle progression of germ cells in Caenorhabditis elegans. By genetically inducing changes in the H3.3 distribution, we demonstrate that both H3.3K27M and pre-existing H3K27me3 act locally and antagonistically on Polycomb Repressive Complex 2 (PRC2) in a concentration-dependent manner. The heterochromatin changes result in extensive gene misregulation, and genetic screening identified upregulation of JNK as an underlying cause of the germcell aberrations. Moreover, JNK inhibition suppresses the replicative fate in human tumor-derived H3.3K27M cells, thus establishing C. elegans as a powerful model for the identification of potential drug targets for treatment of H3.3K27M tumors.

Substitution of lysine 27 with methionine in histone H3.3 (H3.3K27M) is a driver mutation of pediatric high-grade gliomas. Here the authors show that H3.3K27M-mediated alterations in H3K27me3 distribution result in ectopic DNA replication and cell cycle progression of germ cells in Caenorhabditis elegans, through JNK pathway misregulation.

The C. elegans embryonic transcriptome with tissue, time, and alternative splicing resolution

We have used RNA-seq in Caenorhabditis elegans to produce transcription profiles for seven specific embryonic cell populations from gastrulation to the onset of terminal differentiation. The expression data for these seven cell populations, covering major cell lineages and tissues in the worm, reveal the complex and dynamic changes in gene expression, both spatially and temporally. Also, within genes, start sites and exon usage can be highly differential, producing transcripts that are specific to developmental periods or cell lineages. We have also found evidence of novel exons and introns, as well as differential usage of SL1 and SL2 splice leaders. By combining this data set with the modERN ChIP-seq resource, we are able to support and predict gene regulatory relationships. The detailed information on differences and similarities between gene expression in cell lineages and tissues should be of great value to the community and provides a framework for the investigation of expression in individual cells.

NITPicker: selecting time points for follow-up experiments

Background

The design of an experiment influences both what a researcher can measure, as well as how much confidence can be placed in the results. As such, it is vitally important that experimental design decisions do not systematically bias research outcomes. At the same time, making optimal design decisions can produce results leading to statistically stronger conclusions. Deciding where and when to sample are among the most critical aspects of many experimental designs; for example, we might have to choose the time points at which to measure some quantity in a time series experiment. Choosing times which are too far apart could result in missing short bursts of activity. On the other hand, there may be time points which provide very little information regarding the overall behaviour of the quantity in question.

Results

In this study, we develop a tool called NITPicker (Next Iteration Time-point Picker) for selecting optimal time points (or spatial points along a single axis), that eliminates some of the biases caused by human decision-making, while maximising information about the shape of the underlying curves. NITPicker uses ideas from the field of functional data analysis. NITPicker is available on the Comprehensive R Archive Network (CRAN) and code for drawing figures is available on Github (https://github.com/ezer/NITPicker).

Conclusions

NITPicker performs well on diverse real-world datasets that would be relevant for varied biological applications, including designing follow-up experiments for longitudinal gene expression data, weather pattern changes over time, and growth curves.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2717-5) contains supplementary material, which is available to authorized users.

A Model-Based Approach for Identifying Functional Intergenic Transcribed Regions and Noncoding RNAs

With advances in transcript profiling, the presence of transcriptional activities in intergenic regions has been well established. However, whether intergenic expression reflects transcriptional noise or activity of novel genes remains unclear. We identified intergenic transcribed regions (ITRs) in 15 diverse flowering plant species and found that the amount of intergenic expression correlates with genome size, a pattern that could be expected if intergenic expression is largely nonfunctional. To further assess the functionality of ITRs, we first built machine learning models using Arabidopsis thaliana as a model that accurately distinguish functional sequences (benchmark protein-coding and RNA genes) and likely nonfunctional ones (pseudogenes and unexpressed intergenic regions) by integrating 93 biochemical, evolutionary, and sequence-structure features. Next, by applying the models genome-wide, we found that 4,427 ITRs (38%) and 796 annotated ncRNAs (44%) had features significantly similar to benchmark protein-coding or RNA genes and thus were likely parts of functional genes. Approximately 60% of ITRs and ncRNAs were more similar to nonfunctional sequences and were likely transcriptional noise. The predictive framework established here provides not only a comprehensive look at how functional, genic sequences are distinct from likely nonfunctional ones, but also a new way to differentiate novel genes from genomic regions with noisy transcriptional activities.

Single-cell RNA-seq denoising using a deep count autoencoder

Single-cell RNA sequencing (scRNA-seq) has enabled researchers to study gene expression at a cellular resolution. However, noise due to amplification and dropout may obstruct analyses, so scalable denoising methods for increasingly large but sparse scRNA-seq data are needed. We propose a deep count autoencoder network (DCA) to denoise scRNA-seq datasets. DCA takes the count distribution, overdispersion and sparsity of the data into account using a negative binomial noise model with or without zero-inflation, and nonlinear gene-gene dependencies are captured. Our method scales linearly with the number of cells and can, therefore, be applied to datasets of millions of cells. We demonstrate that DCA denoising improves a diverse set of typical scRNA-seq data analyses using simulated and real datasets. DCA outperforms existing methods for data imputation in quality and speed, enhancing biological discovery.

Single-cell RNA sequencing is a powerful method to study gene expression, but noise in the data can obstruct analysis. Here the authors develop a denoising method based on a deep count autoencoder network that scales linearly with the number of cells, and therefore is compatible with large data sets.