Transcriptome age of individual cell types in Caenorhabditis elegans

Divergence in neuronal signaling pathways despite conserved neuronal identity among Caenorhabditis species

One avenue to better understand brain evolution is to map molecular patterns of evolutionary changes in neuronal cell types across entire nervous systems of distantly related species. Generating whole-animal single-cell transcriptomes of three nematode species from the Caenorhabditis genus, we observed a remarkable stability of neuronal-cell-type identities over more than 45 million years of evolution. Conserved patterns of combinatorial expression of homeodomain transcription factors are among the best classifiers of homologous neuron classes. Unexpectedly, we discover an extensive divergence in neuronal signaling pathways. Although identities of neurotransmitter-producing neurons (glutamate, acetylcholine, γ-aminobutyric acid [GABA], and several monoamines) remain stable, expression of ionotropic and metabotropic receptors for all these neurotransmitter systems shows substantial divergence, resulting in more than half of all neuron classes changing their capacity to be receptive to specific neurotransmitters. Neuropeptidergic signaling is also remarkably divergent, both at the level of neuropeptide expression and receptor expression, yet the overall dense network topology of the wireless neuropeptidergic connectome remains stable. Novel neuronal signaling pathways are suggested by our discovery of small secreted proteins that show no obvious hallmarks of conventional neuropeptides but show similar patterns of highly neuron-type-specific and highly evolvable expression profiles. In conclusion, by investigating the evolution of entire nervous systems at the resolution of single-neuron classes, we uncover patterns that may reflect basic principles governing evolutionary novelty in neuronal circuits.

Network properties constrain natural selection on gene expression in Caenorhabditis elegans

Gene regulatory networks (GRNs) integrate genetic and environmental signals to coordinate complex phenotypes and evolve through a balance of selection and drift. Using publicly available datasets from Caenorhabditis elegans, we investigated the extent of natural selection on transcript abundance by linking population-scale variation in gene expression to fecundity, a key fitness component. While the expression of most genes covaried only weakly with fitness, which is typical for polygenic traits, we identified seven transcripts under significant directional selection. These included nhr-114 and feh-1, implicating variation in nutrient-sensing and metabolic pathways as impacting fitness. Stronger directional selection on tissue-specific and older genes highlighted the germline and nervous system as focal points of adaptive change. Network position further constrained selection on gene expression; high-connectivity genes faced stronger stabilizing and directional selection, highlighting GRN architecture as a key factor in microevolutionary dynamics. The activity of transcription factors such as zip-3, which regulates mitochondrial stress responses, emerged as targets of selection, revealing potential links between energy homeostasis and fitness. Our findings demonstrate how GRNs mediate the interplay between selection and drift, shaping microevolutionary trajectories of gene expression and phenotypic diversity.

Molecular patterns of evolutionary changes throughout the whole nervous system of multiple nematode species

One avenue to better understand brain evolution is to map molecular patterns of evolutionary changes in neuronal cell types across entire nervous systems of distantly related species. Generating whole-animal single-cell transcriptomes of three nematode species from the Caenorhabditis genus, we observed a remarkable stability of neuronal cell type identities over more than 45 million years of evolution. Conserved patterns of combinatorial expression of homeodomain transcription factors are among the best classifiers of homologous neuron classes. Unexpectedly, we discover an extensive divergence in neuronal signaling pathways. While identities of neurotransmitter-producing neurons (glutamate, acetylcholine, GABA and several monoamines) remain stable, ionotropic and metabotropic receptors for all these neurotransmitter systems show substantial divergence, resulting in more than half of all neuron classes changing their capacity to be receptive to specific neurotransmitters. Neuropeptidergic signaling is also remarkably divergent, both at the level of neuropeptide expression and receptor expression, yet the overall dense network topology of the wireless neuropeptidergic connectome remains stable. Novel neuronal signaling pathways are suggested by our discovery of small secreted proteins that show no obvious hallmarks of conventional neuropeptides, but show similar patterns of highly neuron-type-specific and highly evolvable expression profiles. In conclusion, by investigating the evolution of entire nervous systems at the resolution of single neuron classes, we uncover patterns that may reflect basic principles governing evolutionary novelty in neuronal circuits.

A transcriptomic hourglass in brown algae

Complex multicellularity has emerged independently across a few eukaryotic lineages and is often associated with the rise of elaborate, tightly coordinated developmental processes1,2. How multicellularity and development are interconnected in evolution is a major question in biology. The hourglass model of embryonic evolution depicts how developmental processes are conserved during evolution, and predicts morphological and molecular divergence in early and late embryogenesis, bridged by a conserved mid-embryonic (phylotypic) period linked to the formation of the basic body plan3,4. Initially found in animal embryos5-8, molecular hourglass patterns have recently been proposed for land plants and fungi9,10. However, whether the hourglass pattern is an intrinsic feature of all complex multicellular eukaryotes remains unknown. Here we tested the presence of a molecular hourglass in the brown algae, a eukaryotic lineage that has evolved multicellularity independently from animals, fungi and plants1,11,12. By exploring transcriptome evolution patterns of brown algae with distinct morphological complexities, we uncovered an hourglass pattern during embryogenesis in morphologically complex species. Filamentous algae without canonical embryogenesis display transcriptome conservation in multicellular stages of the life cycle, whereas unicellular stages are more rapidly evolving. Our findings suggest that transcriptome conservation in brown algae is associated with cell differentiation stages, but is not necessarily linked to embryogenesis. Together with previous work in animals, plants and fungi, we provide further evidence for the generality of a developmental hourglass pattern across complex multicellular eukaryotes.

Single-cell phylotranscriptomics of developmental and cell type evolution

Single-cell phylotranscriptomics is an emerging tool to reveal the molecular and cellular mechanisms of evolution. We summarize its utility in studying the hourglass pattern of ontogenetic evolution and for understanding the evolutionary history of cell types. The developmental hourglass model suggests that the mid-embryonic stage is the most conserved period of development across species, which is supported by morphological and molecular studies. Single-cell phylotranscriptomic analysis has revealed previously underappreciated heterogeneity in transcriptome ages among lineages and cell types throughout development, and has identified the lineages and tissues that drive the whole-organism hourglass pattern. Single-cell transcriptome age analyses also provide important insights into the origin of germ layers, the different selective forces on tissues during adaptation, and the evolutionary relationships between cell types.

Conserved and specific gene expression patterns in the embryonic development of tardigrades

Tardigrades, commonly known as water bears, are enigmatic organisms characterized by their remarkable resilience to extreme environments despite their simple and compact body structure. To date, there is still much to understand about their evolutionary and developmental features contributing to their special body plan and abilities. This research provides preliminary insights on the conserved and specific gene expression patterns during embryonic development of water bears, focusing on the species Hypsibius exemplaris. The developmental dynamic expression analysis of the genes with various evolutionary age grades indicated that the mid-conserved stage of H. exemplaris corresponds to the period of ganglia and midgut development, with the late embryonic stage showing a transition from non-conserved to conserved state. Additionally, a comparison with Drosophila melanogaster highlighted the absence of certain pathway nodes in development-related pathways, such as Maml and Hairless, which are respectively the transcriptional co-activator and co-repressor of NOTCH regulated genes. We also employed Weighted Gene Co-expression Network Analysis (WGCNA) to investigate the expression patterns of tardigrade-specific genes during embryo development. Our findings indicated that the module containing the highest proportion of tardigrade-specific genes (TSGs) exhibits high expression levels before the mid-conserved stage, potentially playing a role in glutathione and lipid metabolism. These functions may be associated to the ecdysone synthesis and storage cell formation, which is unique to tardigrades.

Robustness and complexity

When a system robustly corrects component-level errors, the direct pressure on component performance declines. Components become less reliable, maintain more genetic variability, or drift neutrally, creating new forms of complexity. Examples include the hourglass pattern of biological development and the hourglass architecture for robustly complex systems in engineering.

oggmap: a Python package to extract gene ages per orthogroup and link them with single-cell RNA data

SUMMARY

For model species, single-cell RNA-based cell atlases are available. A good cell atlas includes all major stages in a species' ontogeny, and soon, they will be standard even for nonmodel species. Here, we propose a Python package called oggmap, which allows for the easy extraction of an orthomap (gene ages per orthogroup) for any given query species from OrthoFinder and other gene family data resources, like homologous groups from eggNOG or PLAZA. oggmap provides extracted gene ages for more than thousand eukaryotic species which can be further used to calculate gene age-weighted expression data from scRNA sequencing objects using the Python Scanpy toolkit. Not limited to one transcriptome evolutionary index, oggmap can visualize the individual gene category (e.g. age class, nucleotide diversity bin) and their corresponding expression profiles to investigate scRNA-based cell type assignments in an evolutionary context.

AVAILABILITY AND IMPLEMENTATION

oggmap source code is available at https://github.com/kullrich/oggmap, documentation is available at https://oggmap.readthedocs.io/en/latest/. oggmap can be installed via PyPi or directly used via a docker container.

daf-42 is an evolutionarily young gene essential for dauer development in Caenorhabditis elegans

Under adverse environmental conditions, nematodes arrest into dauer, an alternative developmental stage for diapause. Dauer endures unfavorable environments and interacts with host animals to access favorable environments, thus playing a critical role in survival. Here, we report that in Caenorhabditis elegans, daf-42 is essential for development into the dauer stage, as the null mutant of daf-42 exhibited a "no viable dauer" phenotype in which no viable dauers were obtained in any dauer-inducing conditions. Long-term time lapse microscopy of synchronized larvae revealed that daf-42 is involved in developmental changes from the pre-dauer L2d stage to the dauer stage. daf-42 encodes large, disordered proteins of various sizes that are expressed in and secreted from the seam cells within a narrow time window shortly before the molt into dauer stage. Transcriptome analysis showed that the transcription of genes involved in larval physiology and dauer metabolism is highly affected by the daf-42 mutation. Contrary to the notion that essential genes that control the life and death of an organism may be well conserved across diverse species, daf-42 is an evolutionarily young gene conserved only in the Caenorhabditis genus. Our study shows that dauer formation is a vital process that is controlled not only by conserved genes but also by newly emerged genes, providing important insights into evolutionary mechanisms.