Gene family evolution underlies cell-type diversification in the hypothalamus of teleosts

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.

The decision of male medaka to mate or fight depends on two complementary androgen signaling pathways

Adult male animals typically court and attempt to mate with females, while attacking other males. Emerging evidence from mice indicates that neurons expressing the estrogen receptor ESR1 in behaviorally relevant brain regions play a central role in mediating these mutually exclusive behavioral responses to conspecifics. However, the findings in mice are unlikely to apply to vertebrates in general because, in many species other than rodents and some birds, androgens-rather than estrogens-have been implicated in male behaviors. Here, we report that male medaka (Oryzias latipes) lacking one of the two androgen receptor subtypes (Ara) are less aggressive toward other males and instead actively court them, while those lacking the other subtype (Arb) are less motivated to mate with females and conversely attack them. These findings indicate that, in male medaka, the Ara- and Arb-mediated androgen signaling pathways facilitate appropriate behavioral responses, while simultaneously suppressing inappropriate responses, to males and females, respectively. Notably, males lacking either receptor retain the ability to discriminate the sex of conspecifics, suggesting a defect in the subsequent decision-making process to mate or fight. We further show that Ara and Arb are expressed in intermingled but largely distinct populations of neurons, and stimulate the expression of different behaviorally relevant genes including galanin and vasotocin, respectively. Collectively, our results demonstrate that male teleosts make adaptive decisions to mate or fight as a result of the activation of one of two complementary androgen signaling pathways, depending on the sex of the conspecific that they encounter.

CACIMAR: cross-species analysis of cell identities, markers, regulations, and interactions using single-cell RNA sequencing data

Transcriptomic analysis across species is increasingly used to reveal conserved gene regulations which implicate crucial regulators. Cross-species analysis of single-cell RNA sequencing (scRNA-seq) data provides new opportunities to identify the cellular and molecular conservations, especially for cell types and cell type-specific gene regulations. However, few methods have been developed to analyze cross-species scRNA-seq data to uncover both molecular and cellular conservations. Here, we built a tool called CACIMAR, which can perform cross-species analysis of cell identities, markers, regulations, and interactions using scRNA-seq profiles. Based on the weighted sum models of the conserved features, we developed different conservation scores to measure the conservation of cell types, regulatory networks, and intercellular interactions. Using publicly available scRNA-seq data on retinal regeneration in mice, zebrafish, and chick, we demonstrated four main functions of CACIMAR. First, CACIMAR allows to identify conserved cell types even in evolutionarily distant species. Second, the tool facilitates the identification of evolutionarily conserved or species-specific marker genes. Third, CACIMAR enables the identification of conserved intracellular regulations, including cell type-specific regulatory subnetworks and regulators. Lastly, CACIMAR provides a unique feature for identifying conserved intercellular interactions. Overall, CACIMAR facilitates the identification of evolutionarily conserved cell types, marker genes, intracellular regulations, and intercellular interactions, providing insights into the cellular and molecular mechanisms of species evolution.

Barcoding Notch signaling in the developing brain

Developmental signaling inputs are fundamental for shaping cell fates and behavior. However, traditional fluorescent-based signaling reporters have limitations in scalability and molecular resolution of cell types. We present SABER-seq, a CRISPR-Cas molecular recorder that stores transient developmental signaling cues as permanent mutations in cellular genomes for deconstruction at later stages via single-cell transcriptomics. We applied SABER-seq to record Notch signaling in developing zebrafish brains. SABER-seq has two components: a signaling sensor and a barcode recorder. The sensor activates Cas9 in a Notch-dependent manner with inducible control while the recorder accumulates mutations that represent Notch activity in founder cells. We combine SABER-seq with an expanded juvenile brain atlas to define cell types whose fates are determined downstream of Notch signaling. We identified examples wherein Notch signaling may have differential impact on terminal cell fates. SABER-seq is a novel platform for rapid, scalable and high-resolution mapping of signaling activity during development.

Adult sex change leads to extensive forebrain reorganization in clownfish

Sexual differentiation of the brain occurs in all major vertebrate lineages but is not well understood at a molecular and cellular level. Unlike most vertebrates, sex-changing fishes have the remarkable ability to change reproductive sex during adulthood in response to social stimuli, offering a unique opportunity to understand mechanisms by which the nervous system can initiate and coordinate sexual differentiation. This study explores sexual differentiation of the forebrain using single nucleus RNA-sequencing in the anemonefish Amphiprion ocellaris, producing the first cellular atlas of a sex-changing brain. We uncover extensive sex differences in cell type-specific gene expression, relative proportions of cells, baseline neuronal excitation, and predicted inter-neuronal communication. Additionally, we identify the cholecystokinin, galanin, and estrogen systems as central molecular axes of sexual differentiation. Supported by these findings, we propose a model of neurosexual differentiation in the conserved vertebrate social decision-making network spanning multiple subtypes of neurons and glia, including neuronal subpopulations within the preoptic area that are positioned to regulate gonadal differentiation. This work deepens our understanding of sexual differentiation in the vertebrate brain and defines a rich suite of molecular and cellular pathways that differentiate during adult sex change in anemonefish.

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes

BACKGROUND

How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum have long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species.

RESULTS

The final assembled and filtered Berghia genome is comparable to other high-quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes) and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low.

CONCLUSIONS

Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

Single-cell analysis of shared signatures and transcriptional diversity during zebrafish development

During development, animals generate distinct cell populations with specific identities, functions, and morphologies. We mapped transcriptionally distinct populations across 489,686 cells from 62 stages during wild-type zebrafish embryogenesis and early larval development (3-120 h post-fertilization). Using these data, we identified the limited catalog of gene expression programs reused across multiple tissues and their cell-type-specific adaptations. We also determined the duration each transcriptional state is present during development and identify unexpected long-term cycling populations. Focused clustering and transcriptional trajectory analyses of non-skeletal muscle and endoderm identified transcriptional profiles and candidate transcriptional regulators of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and recently discovered best4+ cells. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

Single-cell transcriptomics reveals the brain evolution of web-building spiders

Spiders are renowned for their efficient capture of flying insects using intricate aerial webs. How the spider nervous systems evolved to cope with this specialized hunting strategy and various environmental clues in an aerial space remains unknown. Here we report a brain-cell atlas of >30,000 single-cell transcriptomes from a web-building spider (Hylyphantes graminicola). Our analysis revealed the preservation of ancestral neuron types in spiders, including the potential coexistence of noradrenergic and octopaminergic neurons, and many peptidergic neuronal types that are lost in insects. By comparing the genome of two newly sequenced plesiomorphic burrowing spiders with three aerial web-building spiders, we found that the positively selected genes in the ancestral branch of web-building spiders were preferentially expressed (42%) in the brain, especially in the three mushroom body-like neuronal types. By gene enrichment analysis and RNAi experiments, these genes were suggested to be involved in the learning and memory pathway and may influence the spiders' web-building and hunting behaviour. Our results provide key sources for understanding the evolution of behaviour in spiders and reveal how molecular evolution drives neuron innovation and the diversification of associated complex behaviours.

Single-cell transcriptomics refuels the exploration of spiralian biology

Spiralians represent the least studied superclade of bilaterian animals, despite exhibiting the widest diversity of organisms. Although spiralians include iconic organisms, such as octopus, earthworms and clams, a lot remains to be discovered regarding their phylogeny and biology. Here, we review recent attempts to apply single-cell transcriptomics, a new pioneering technology enabling the classification of cell types and the characterisation of their gene expression profiles, to several spiralian taxa. We discuss the methodological challenges and requirements for applying this approach to marine organisms and explore the insights that can be brought by such studies, both from a biomedical and evolutionary perspective. For instance, we show that single-cell sequencing might help solve the riddle of the homology of larval forms across spiralians, but also to better characterise and compare the processes of regeneration across taxa. We highlight the capacity of single-cell to investigate the origin of evolutionary novelties, as the mollusc shell or the cephalopod visual system, but also to interrogate the conservation of the molecular fingerprint of cell types at long evolutionary distances. We hope that single-cell sequencing will open a new window in understanding the biology of spiralians, and help renew the interest for these overlooked but captivating organisms.

A telencephalon cell type atlas for goldfish reveals diversity in the evolution of spatial structure and cell types

Teleost fish form the largest group of vertebrates and show a tremendous variety of adaptive behaviors, making them critically important for the study of brain evolution and cognition. The neural basis mediating these behaviors remains elusive. We performed a systematic comparative survey of the goldfish telencephalon. We mapped cell types using single-cell RNA sequencing and spatial transcriptomics, resulting in de novo molecular neuroanatomy parcellation. Glial cells were highly conserved across 450 million years of evolution separating mouse and goldfish, while neurons showed diversity and modularity in gene expression. Specifically, somatostatin interneurons, famously interspersed in the mammalian isocortex for local inhibitory input, were curiously aggregated in a single goldfish telencephalon nucleus but molecularly conserved. Cerebral nuclei including the striatum, a hub for motivated behavior in amniotes, had molecularly conserved goldfish homologs. We suggest elements of a hippocampal formation across the goldfish pallium. Last, aiding study of the teleostan everted telencephalon, we describe substantial molecular similarities between goldfish and zebrafish neuronal taxonomies.

Distinct and shared endothermic strategies in the heat producing tissues of tuna and other teleosts

Although most fishes are ectothermic, some, including tuna and billfish, achieve endothermy through specialized heat producing tissues that are modified muscles. How these heat producing tissues evolved, and whether they share convergent molecular mechanisms, remain unresolved. Here, we generated a high-quality genome from the mackerel tuna (Euthynnus affinis) and investigated the heat producing tissues of this fish by single-nucleus and bulk RNA sequencing. Compared with other teleosts, tuna-specific genetic variation is strongly associated with muscle differentiation. Single-nucleus RNA-seq revealed a high proportion of specific slow skeletal muscle cell subtypes in the heat producing tissues of tuna. Marker genes of this cell subtype are associated with the relative sliding of actin and myosin, suggesting that tuna endothermy is mainly based on shivering thermogenesis. In contrast, cross-species transcriptome analysis indicated that endothermy in billfish relies mainly on non-shivering thermogenesis. Nevertheless, the heat producing tissues of the different species do share some tissue-specific genes, including vascular-related and mitochondrial genes. Overall, although tunas and billfishes differ in their thermogenic strategies, they share similar expression patterns in some respects, highlighting the complexity of convergent evolution.

Benchmarking strategies for cross-species integration of single-cell RNA sequencing data

The growing number of available single-cell gene expression datasets from different species creates opportunities to explore evolutionary relationships between cell types across species. Cross-species integration of single-cell RNA-sequencing data has been particularly informative in this context. However, in order to do so robustly it is essential to have rigorous benchmarking and appropriate guidelines to ensure that integration results truly reflect biology. Here, we benchmark 28 combinations of gene homology mapping methods and data integration algorithms in a variety of biological settings. We examine the capability of each strategy to perform species-mixing of known homologous cell types and to preserve biological heterogeneity using 9 established metrics. We also develop a new biology conservation metric to address the maintenance of cell type distinguishability. Overall, scANVI, scVI and SeuratV4 methods achieve a balance between species-mixing and biology conservation. For evolutionarily distant species, including in-paralogs is beneficial. SAMap outperforms when integrating whole-body atlases between species with challenging gene homology annotation. We provide our freely available cross-species integration and assessment pipeline to help analyse new data and develop new algorithms.

A ligand-receptor interactome atlas of the zebrafish

Studies in zebrafish can unravel the functions of cellular communication and thus identify novel bench-to-bedside drugs targeting cellular communication signaling molecules. Due to the incomplete annotation of zebrafish proteome, the knowledge of zebrafish receptors, ligands, and tools to explore their interactome is limited. To address this gap, we de novo predicted the cellular localization of zebrafish reference proteome using deep learning algorithm. We combined the predicted and existing annotations on cellular localization of zebrafish proteins and created repositories of zebrafish ligands, membrane receptome, and interactome as well as associated diseases and targeting drugs. Unlike other tools, our interactome atlas is based on both the physical interaction data of zebrafish proteome and existing human ligand-receptor pair databases. The resources are available as R and Python scripts. DanioTalk provides a novel resource for researchers interested in targeting cellular communication in zebrafish, as we demonstrate in applications studying synapse and axo-glial interactome. DanioTalk methodology can be applied to build and explore the ligand-receptor atlas of other non-mammalian model organisms.

Evolutionary dynamics of mushroom body Kenyon cell types in hymenopteran brains from multifunctional type to functionally specialized types

Evolutionary dynamics of diversification of brain neuronal cell types that have underlain behavioral evolution remain largely unknown. Here, we compared transcriptomes and functions of Kenyon cell (KC) types that compose the mushroom bodies between the honey bee and sawfly, a primitive hymenopteran insect whose KCs likely have the ancestral properties. Transcriptome analyses show that the sawfly KC type shares some of the gene expression profile with each honey bee KC type, although unique gene expression profiles have also been acquired in each honey bee KC type. In addition, functional analysis of two sawfly genes suggested that the functions in learning and memory of the ancestral KC type were heterogeneously inherited among the KC types in the honey bee. Our findings strongly suggest that the functional evolution of KCs in Hymenoptera involved two previously hypothesized processes for evolution of cell function: functional segregation and divergence.

Genomic and transcriptomic analyses support a silk gland origin of spider venom glands

BACKGROUND

Spiders comprise a hyperdiverse lineage of predators with venom systems, yet the origin of functionally novel spider venom glands remains unclear. Previous studies have hypothesized that spider venom glands originated from salivary glands or evolved from silk-producing glands present in early chelicerates. However, there is insufficient molecular evidence to indicate similarity among them. Here, we provide comparative analyses of genome and transcriptome data from various lineages of spiders and other arthropods to advance our understanding of spider venom gland evolution.

RESULTS

We generated a chromosome-level genome assembly of a model spider species, the common house spider (Parasteatoda tepidariorum). Module preservation, GO semantic similarity, and differentially upregulated gene similarity analyses demonstrated a lower similarity in gene expressions between the venom glands and salivary glands compared to the silk glands, which questions the validity of the salivary gland origin hypothesis but unexpectedly prefers to support the ancestral silk gland origin hypothesis. The conserved core network in the venom and silk glands was mainly correlated with transcription regulation, protein modification, transport, and signal transduction pathways. At the genetic level, we found that many genes in the venom gland-specific transcription modules show positive selection and upregulated expressions, suggesting that genetic variation plays an important role in the evolution of venom glands.

CONCLUSIONS

This research implies the unique origin and evolutionary path of spider venom glands and provides a basis for understanding the diverse molecular characteristics of venom systems.

Evolution of homology: From archetype towards a holistic concept of cell type

The concept of homology lies in the heart of comparative biological science. The distinction between homology as structure and analogy as function has shaped the evolutionary paradigm for a century and formed the axis of comparative anatomy and embryology, which accept the identity of structure as a ground measure of relatedness. The advent of single-cell genomics overturned the classical view of cell homology by establishing a backbone regulatory identity of cell types, the basic biological units bridging the molecular and phenotypic dimensions, to reveal that the cell is the most flexible unit of living matter and that many approaches of classical biology need to be revised to understand evolution and diversity at the cellular level. The emerging theory of cell types explicitly decouples cell identity from phenotype, essentially allowing for the divergence of evolutionarily related morphotypes beyond recognition, as well as it decouples ontogenetic cell lineage from cell-type phylogeny, whereby explicating that cell types can share common descent regardless of their structure, function or developmental origin. The article succinctly summarizes current progress and opinion in this field and formulates a more generalistic view of biological cell types as avatars, transient or terminal cell states deployed in a continuum of states by the developmental programme of one and the same omnipotent cell, capable of changing or combining identities with distinct evolutionary histories or inventing ad hoc identities that never existed in evolution or development. It highlights how the new logic grounded in the regulatory nature of cell identity transforms the concepts of cell homology and phenotypic stability, suggesting that cellular evolution is inherently and massively network-like, with one-to-one homologies being rather uncommon and restricted to shallower levels of the animal tree of life.

Cell type diversity in a developing octopus brain

Octopuses are mollusks that have evolved intricate neural systems comparable with vertebrates in terms of cell number, complexity and size. The brain cell types that control their sophisticated behavioral repertoire are still unknown. Here, we profile the cell diversity of the paralarval Octopus vulgaris brain to build a cell type atlas that comprises mostly neural cells, but also multiple glial subtypes, endothelial cells and fibroblasts. We spatially map cell types to the vertical, subesophageal and optic lobes. Investigation of cell type conservation reveals a shared gene signature between glial cells of mouse, fly and octopus. Genes related to learning and memory are enriched in vertical lobe cells, which show molecular similarities with Kenyon cells in Drosophila. We construct a cell type taxonomy revealing transcriptionally related cell types, which tend to appear in the same brain region. Together, our data sheds light on cell type diversity and evolution in the octopus brain.

Molecular diversity and evolution of neuron types in the amniote brain

The existence of evolutionarily conserved regions in the vertebrate brain is well established. The rules and constraints underlying the evolution of neuron types, however, remain poorly understood. To compare neuron types across brain regions and species, we generated a cell type atlas of the brain of a bearded dragon and compared it with mouse datasets. Conserved classes of neurons could be identified from the expression of hundreds of genes, including homeodomain-type transcription factors and genes involved in connectivity. Within these classes, however, there are both conserved and divergent neuron types, precluding a simple categorization of the brain into ancestral and novel areas. In the thalamus, neuronal diversification correlates with the evolution of the cortex, suggesting that developmental origin and circuit allocation are drivers of neuronal identity and evolution.

A Primer for Single-Cell Sequencing in Non-Model Organisms

Single-cell sequencing technologies have led to a revolution in our knowledge of the diversity of cell types, connections between biological levels of organization, and relationships between genotype and phenotype. These advances have mainly come from using model organisms; however, using single-cell sequencing in non-model organisms could enable investigations of questions inaccessible with typical model organisms. This primer describes a general workflow for single-cell sequencing studies and considerations for using non-model organisms (limited to multicellular animals). Importantly, single-cell sequencing, when further applied in non-model organisms, will allow for a deeper understanding of the mechanisms between genotype and phenotype and the basis for biological variation.

Discordant Genome Assemblies Drastically Alter the Interpretation of Single-Cell RNA Sequencing Data Which Can Be Mitigated by a Novel Integration Method

Advances in sequencing and assembly technology have led to the creation of genome assemblies for a wide variety of non-model organisms. The rapid production and proliferation of updated, novel assembly versions can create vexing problems for researchers when multiple-genome assembly versions are available at once, requiring researchers to work with more than one reference genome. Multiple-genome assemblies are especially problematic for researchers studying the genetic makeup of individual cells, as single-cell RNA sequencing (scRNAseq) requires sequenced reads to be mapped and aligned to a single reference genome. Using the Astyanax mexicanus, this study highlights how the interpretation of a single-cell dataset from the same sample changes when aligned to its two different available genome assemblies. We found that the number of cells and expressed genes detected were drastically different when aligning to the different assemblies. When the genome assemblies were used in isolation with their respective annotations, cell-type identification was confounded, as some classic cell-type markers were assembly-specific, whilst other genes showed differential patterns of expression between the two assemblies. To overcome the problems posed by multiple-genome assemblies, we propose that researchers align to each available assembly and then integrate the resultant datasets to produce a final dataset in which all genome alignments can be used simultaneously. We found that this approach increased the accuracy of cell-type identification and maximised the amount of data that could be extracted from our single-cell sample by capturing all possible cells and transcripts. As scRNAseq becomes more widely available, it is imperative that the single-cell community is aware of how genome assembly alignment can alter single-cell data and their interpretation, especially when reviewing studies on non-model organisms.