Cell type transcriptome atlas for the planarian Schmidtea mediterranea

Planarians require ced-12/elmo-1 to clear dead cells by excretion through the gut

Cell corpse removal is a critical component of both development and homeostasis throughout the animal kingdom. Extensive research has revealed many of the mechanisms involved in corpse removal, typically involving engulfment and digestion by another cell; however, the dynamics of cell corpse clearance in adult tissues remain unclear. Here, we track cell death in the adult planarian Schmidtea mediterranea and find that, following light-induced cell death, pigment cell corpses transit to the gut and are excreted from the animal. Gut phagocytes, previously only known to phagocytose food, are required for pigment cells to enter the gut lumen. Finally, we show that the planarian ortholog of ced-12/engulfment and cell motility (ELMO) is required for corpse phagocytosis and removal through the gut. In total, we present a mechanism of cell clearance in an adult organism involving transit of dead cells to the gut, transport into the gut by phagocytes, and physical excretion of debris.

Cell cycle specific, differentially tagged ribosomal proteins to measure phase specific transcriptomes from asynchronously cycling cells

Asynchronously cycling cells pose a challenge to the accurate characterization of phase-specific gene expression. Current strategies, including RNAseq, survey the steady state gene expression across the cell cycle and are inherently limited by their inability to resolve dynamic gene regulatory networks. Single cell RNAseq (scRNAseq) can identify different cell cycle transcriptomes if enough cycling cells are present, however some cells are not amenable to scRNAseq. Therefore, we merged two powerful strategies, the CDT1 and GMNN degrons used in Fluorescent Ubiquitination-based Cell Cycle Indicator (FUCCI) cell cycle sensors and the ribosomal protein epitope tagging used in RiboTrap/Tag technologies to isolate cell cycle phase-specific mRNA for sequencing. The resulting cell cycle dependent, tagged ribosomal proteins (ccTaggedRP) were differentially expressed during the cell cycle, had similar subcellular locations as endogenous ribosomal proteins, incorporated into ribosomes and polysomes, and facilitated the recovery of cell cycle phase-specific RNA for sequencing. ccTaggedRP has broad applications to investigate phase-specific gene expression in complex cell populations.

A resource of single-cell gene expression profiles in a planarian Dugesia japonica

The freshwater planarian Dugesia japonica maintains an abundant heterogeneous cell population called neoblasts, which include adult pluripotent stem cells. Thus, it is an excellent model organism for stem cell and regeneration research. Recently, many single-cell RNA sequencing (scRNA-seq) databases of several model organisms, including other planarian species, have become publicly available; these are powerful and useful resources to search for gene expression in various tissues and cells. However, the only scRNA-seq dataset for D. japonica has been limited by the number of genes detected. Herein, we collected D. japonica cells, and conducted an scRNA-seq analysis. A novel, automatic, iterative cell clustering strategy produced a dataset of 3,404 cells, which could be classified into 63 cell types based on gene expression profiles. We introduced two examples for utilizing the scRNA-seq dataset in this study using D. japonica. First, the dataset provided results consistent with previous studies as well as novel functionally relevant insights, that is, the expression of DjMTA and DjP2X-A genes in neoblasts that give rise to differentiated cells. Second, we conducted an integrative analysis of the scRNA-seq dataset and time-course bulk RNA-seq of irradiated animals, demonstrating that the dataset can help interpret differentially expressed genes captured via bulk RNA-seq. Using the R package "Seurat" and GSE223927, researchers can easily access and utilize this dataset.

MRLC controls apoptotic cell death and functions to regulate epidermal development during planarian regeneration and homeostasis

Adult stem cells (ASCs) are pluripotent cells with the capacity to self-renew and constantly replace lost cells due to physiological turnover or injury. Understanding the molecular mechanisms of the precise coordination of stem cell proliferation and proper cell fate decision is important to regeneration and organismal homeostasis. The planarian epidermis provides a highly tractable model to study ASC complex dynamic due to the distinct spatiotemporal differentiation stages during lineage development. Here, we identified the myosin regulatory light chain (MRLC) homologue in the Dugesia japonica transcriptome. We found high expression levels of MRLC in wound region during regeneration and also expressed in late epidermal progenitors as an essential regulator of the lineage from neoblasts to mature epidermal cells. We investigated the function of MRLC using in situ hybridization, real-time polymerase chain reaction and double fluorescent and uncovered the potential mechanism. Knockdown of MRLC leads to a remarkable increase in cell death, causes severe abnormalities during regeneration and homeostasis and eventually leads to animal death. The global decrease in epidermal cell in MRLC RNAi animals induces accelerated epidermal proliferation and differentiation. Additionally, we find that MRLC is co-expressed with cdc42 and acts cooperatively to control the epidermal lineage development by affecting cell death. Our results uncover an important role of MRLC, as an inhibitor of apoptosis, involves in epidermal development.

Purification of Planarian Stem Cells Using a Draq5-Based FACS Approach

Planarians are flatworms that have the remarkable ability to regenerate entirely new animals. This regenerative ability requires abundant adult stem cells called neoblasts, which are relatively small in size, sensitive to irradiation and the only proliferative cells in the animal. Despite the lack of cell surface markers, fluorescence-activated cell sorting (FACS) protocols have been developed to discriminate and isolate neoblasts, based on DNA content. Here, we describe a protocol that combines staining of far-red DNA dye Draq5, Calcein-AM and DAPI, along with a shortened processing time. This profiling strategy can be used to functionally characterize the neoblast population in pharmacologically-treated or gene knockdown animals. Highly purified neoblasts can be analyzed with downstream assays, such as in situ hybridization and RNA sequencing.

SINGLE-CELL TRANSCRIPTOME ANALYSIS IN HEALTH AND DISEASE

ABSTRACT

The analysis of the single-cell transcriptome has emerged as a powerful tool to gain insights on the basic mechanisms of health and disease. It is widely used to reveal the cellular diversity and complexity of tissues at cellular resolution by RNA sequencing of the whole transcriptome from a single cell. Equally, it is applied to discover an unknown, rare population of cells in the tissue. The prime advantage of single-cell transcriptome analysis is the detection of stochastic nature of gene expression of the cell in tissue. Moreover, the availability of multiple platforms for the single-cell transcriptome has broadened its approaches to using cells of different sizes and shapes, including the capture of short or full-length transcripts, which is helpful in the analysis of challenging biological samples. And with the development of numerous packages in R and Python, new directions in the computational analysis of single-cell transcriptomes can be taken to characterize healthy versus diseased tissues to obtain novel pathological insights. Downstream analysis such as differential gene expression analysis, gene ontology term analysis, Kyoto Encyclopedia of Genes and Genomes pathway analysis, cell-cell interaction analysis, and trajectory analysis has become standard practice in the workflow of single-cell transcriptome analysis to further examine the biology of different cell types. Here, we provide a broad overview of single-cell transcriptome analysis in health and disease conditions currently applied in various studies.

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.

Marine Invertebrates One Cell at A Time: Insights from Single-Cell Analysis

Over the past decade, single-cell RNA-sequencing (scRNA-seq) has made it possible to study the cellular diversity of a broad range of organisms. Technological advances in single-cell isolation and sequencing have expanded rapidly, allowing the transcriptomic profile of individual cells to be captured. As a result, there has been an explosion of cell type atlases created for many different marine invertebrate species from across the tree of life. Our focus in this review is to synthesize current literature on marine invertebrate scRNA-seq. Specifically, we provide perspectives on key insights from scRNA-seq studies, including descriptive studies of cell type composition, how cells respond in dynamic processes such as development and regeneration, and the evolution of new cell types. Despite these tremendous advances, there also lie several challenges ahead. We discuss the important considerations that are essential when making comparisons between experiments, or between datasets from different species. Finally, we address the future of single-cell analyses in marine invertebrates, including combining scRNA-seq data with other 'omics methods to get a fuller understanding of cellular complexities. The full diversity of cell types across marine invertebrates remains unknown and understanding this diversity and evolution will provide rich areas for future study.

A niche-derived non-ribosomal peptide triggers planarian sexual development

Germ cells are regulated by local microenvironments (niches), which secrete instructive cues. Conserved developmental signaling molecules act as niche-derived regulatory factors, yet other types of niche signals remain to be identified. Single-cell RNA-sequencing of sexual planarians revealed niche cells expressing a non-ribosomal peptide synthetase (nrps). Inhibiting nrps led to loss of female reproductive organs and testis hyperplasia. Mass spectrometry detected the dipeptide β-alanyl-tryptamine (BATT), which is associated with reproductive system development and requires nrps and a monoamine-transmitter-synthetic enzyme (AADC) for its production. Exogenous BATT rescued the reproductive defects after nrps or aadc inhibition, restoring fertility. Thus, a non-ribosomal, monoamine-derived peptide provided by niche cells acts as a critical signal to trigger planarian reproductive development. These findings reveal an unexpected function for monoamines in niche-germ cell signaling. Furthermore, given the recently reported role for BATT as a male-derived factor required for reproductive maturation of female schistosomes, these results have important implications for the evolution of parasitic flatworms and suggest a potential role for non-ribosomal peptides as signaling molecules in other organisms.

LIM-HD transcription factors control axial patterning and specify distinct neuronal and intestinal cell identities in planarians

Adult planarians can regenerate the gut, eyes and even a functional brain. Proper identity and patterning of the newly formed structures require signals that guide and commit their adult stem cells. During embryogenesis, LIM-homeodomain (LIM-HD) transcription factors act in a combinatorial 'LIM code' to control cell fate determination and differentiation. However, our understanding about the role these genes play during regeneration and homeostasis is limited. Here, we report the full repertoire of LIM-HD genes in Schmidtea mediterranea. We found that lim homeobox (lhx) genes appear expressed in complementary patterns along the cephalic ganglia and digestive system of the planarian, with some of them being co-expressed in the same cell types. We have identified that Smed-islet1, -lhx1/5-1, -lhx2/9-3, -lhx6/8, -lmx1a/b-2 and -lmx1a/b-3 are essential to pattern and size the planarian brain as well as for correct regeneration of specific subpopulations of dopaminergic, serotonergic, GABAergic and cholinergic neurons, while Smed-lhx1/5.2 and -lhx2/9.2 are required for the proper expression of intestinal cell type markers, specifically the goblet subtype. LIM-HD are also involved in controlling axonal pathfinding (lhx6/8), axial patterning (islet1, lhx1/5-1, lmx1a/b-3), head/body proportions (islet2) and stem cell proliferation (lhx3/4, lhx2/9-3, lmx1a/b-2, lmx1a/b-3). Altogether, our results suggest that planarians might present a combinatorial LIM code that controls axial patterning and axonal growing and specifies distinct neuronal and intestinal cell identities.

Functional genomics in Spiralia

Our understanding of the mechanisms that modulate gene expression in animals is strongly biased by studying a handful of model species that mainly belong to three groups: Insecta, Nematoda and Vertebrata. However, over half of the animal phyla belong to Spiralia, a morphologically and ecologically diverse animal clade with many species of economic and biomedical importance. Therefore, investigating genome regulation in this group is central to uncovering ancestral and derived features in genome functioning in animals, which can also be of significant societal impact. Here, we focus on five aspects of gene expression regulation to review our current knowledge of functional genomics in Spiralia. Although some fields, such as single-cell transcriptomics, are becoming more common, the study of chromatin accessibility, DNA methylation, histone post-translational modifications and genome architecture are still in their infancy. Recent efforts to generate chromosome-scale reference genome assemblies for greater species diversity and optimise state-of-the-art approaches for emerging spiralian research systems will address the existing knowledge gaps in functional genomics in this animal group.

Spiralian genomics and the evolution of animal genome architecture

Recent developments in sequencing technologies have greatly improved our knowledge of phylogenetic relationships and genomic architectures throughout the tree of life. Spiralia, a diverse clade within Protostomia, is essential for understanding the evolutionary history of parasitism, gene conversion, nervous systems and animal body plans. In this review, we focus on the current hypotheses of spiralian phylogeny and investigate the impact of long-read sequencing on the quality of genome assemblies. We examine chromosome-level assemblies to highlight key genomic features that have driven spiralian evolution, including karyotype, synteny and the Hox gene organization. In addition, we show how chromosome rearrangement has influenced spiralian genomic structures. Although spiralian genomes have undergone substantial changes, they exhibit both conserved and lineage-specific features. We recommend increasing sequencing efforts and expanding functional genomics research to deepen insights into spiralian biology.

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.

Fate specification is spatially intermingled across planarian stem cells

Regeneration requires mechanisms for producing a wide array of cell types. Neoblasts are stem cells in the planarian Schmidtea mediterranea that undergo fate specification to produce over 125 adult cell types. Fate specification in neoblasts can be regulated through expression of fate-specific transcription factors. We utilize multiplexed error-robust fluorescence in situ hybridization (MERFISH) and whole-mount FISH to characterize fate choice distribution of stem cells within planarians. Fate choices are often made distant from target tissues and in a highly intermingled manner, with neighboring neoblasts frequently making divergent fate choices for tissues of different location and function. We propose that pattern formation is driven primarily by the migratory assortment of progenitors from mixed and spatially distributed fate-specified stem cells and that fate choice involves stem-cell intrinsic processes.

CRISPR/Cas9-based depletion of 16S ribosomal RNA improves library complexity of single-cell RNA-sequencing in planarians

BACKGROUND

Single-cell RNA-sequencing (scRNA-seq) relies on PCR amplification to retrieve information from vanishingly small amounts of starting material. To selectively enrich mRNA from abundant non-polyadenylated transcripts, poly(A) selection is a key step during library preparation. However, some transcripts, such as mitochondrial genes, can escape this elimination and overwhelm libraries. Often, these transcripts are removed in silico, but whether physical depletion improves detection of rare transcripts in single cells is unclear.

RESULTS

We find that a single 16S ribosomal RNA is widely enriched in planarian scRNA-seq datasets, independent of the library preparation method. To deplete this transcript from scRNA-seq libraries, we design 30 single-guide RNAs spanning its length. To evaluate the effects of depletion, we perform a side-by-side comparison of the effects of eliminating the 16S transcript and find a substantial increase in the number of genes detected per cell, coupled with virtually complete loss of the 16S RNA. Moreover, we systematically determine that library complexity increases with a limited number of PCR cycles following CRISPR treatment. When compared to in silico depletion of 16S, physically removing it reduces dropout rates, retrieves more clusters, and reveals more differentially expressed genes.

CONCLUSIONS

Our results show that abundant transcripts reduce the retrieval of informative transcripts in scRNA-seq and distort the analysis. Physical removal of these contaminants enables the detection of rare transcripts at lower sequencing depth, and also outperforms in silico depletion. Importantly, this method can be easily customized to deplete any abundant transcript from scRNA-seq libraries.

The BAF chromatin remodeling complex licenses planarian stem cells access to ectodermal and mesodermal cell fates

BACKGROUND

The flatworm planarian, Schmidtea mediterranea, has a large population of adult stem cells (ASCs) that replace any cell type during tissue turnover or regeneration. How planarian ASCs (called neoblasts) manage self-renewal with the ability to produce daughter cells of different cell lineages (multipotency) is not well understood. Chromatin remodeling complexes ultimately control access to DNA regions of chromosomes and together with specific transcription factors determine whether a gene is transcribed in a given cell type. Previous work in planarians determined that RNAi of core components of the BAF chromatin remodeling complex, brg1 and smarcc2, caused increased ASCs and failed regeneration, but how these cellular defects arise at the level of gene regulation in neoblasts is unknown.

RESULTS

Here, we perform ATAC and RNA sequencing on purified neoblasts, deficient for the BAF complex subunits brg-1 and smarcc2. The data demonstrate that the BAF complex promotes chromatin accessibility and facilitates transcription at target loci, as in other systems. Interestingly, we find that the BAF complex enables access to genes known to be required for the generation of mesoderm- and ectoderm-derived lineages, including muscle, parenchymal cathepsin, neural, and epithelial lineages. BAF complex knockdowns result in disrupted differentiation into these cell lineages and functional consequences on planarian regeneration and tissue turnover. Notably, we did not detect a role for the BAF complex in neoblasts making endodermal lineages.

CONCLUSIONS

Our study provides functional insights into how the BAF complex contributes to cell fate decisions in planarian ASCs in vivo.

Reappraisal of the Concept of Accelerated Aging in Neurodegeneration and Beyond

BACKGROUND

Genetic and epigenetic changes, oxidative stress and inflammation influence the rate of aging, which diseases, lifestyle and environmental factors can further accelerate. In accelerated aging (AA), the biological age exceeds the chronological age.

OBJECTIVE

The objective of this study is to reappraise the AA concept critically, considering its weaknesses and limitations.

METHODS

We reviewed more than 300 recent articles dealing with the physiology of brain aging and neurodegeneration pathophysiology.

RESULTS

(1) Application of the AA concept to individual organs outside the brain is challenging as organs of different systems age at different rates. (2) There is a need to consider the deceleration of aging due to the potential use of the individual structure-functional reserves. The latter can be restored by pharmacological and/or cognitive therapy, environment, etc. (3) The AA concept lacks both standardised terminology and methodology. (4) Changes in specific molecular biomarkers (MBM) reflect aging-related processes; however, numerous MBM candidates should be validated to consolidate the AA theory. (5) The exact nature of many potential causal factors, biological outcomes and interactions between the former and the latter remain largely unclear.

CONCLUSIONS

Although AA is commonly recognised as a perspective theory, it still suffers from a number of gaps and limitations that assume the necessity for an updated AA concept.

Chromatin analysis of adult pluripotent stem cells reveals a unique stemness maintenance strategy

Many highly regenerative organisms maintain adult pluripotent stem cells throughout their life, but how the long-term maintenance of pluripotency is accomplished is unclear. To decipher the regulatory logic of adult pluripotent stem cells, we analyzed the chromatin organization of stem cell genes in the planarian Schmidtea mediterranea. We identify a special chromatin state of stem cell genes, which is distinct from that of tissue-specific genes and resembles constitutive genes. Where tissue-specific promoters have detectable transcription factor binding sites, the promoters of stem cell-specific genes instead have sequence features that broadly decrease nucleosome binding affinity. This genic organization makes pluripotency-related gene expression the default state in these cells, which is maintained by the activity of chromatin remodelers ISWI and SNF2 in the stem cells.

Deciphering m6A dynamics at a single-base level during planarian anterior-posterior axis specification

BACKGROUND

The establishment of the anterior-posterior (A-P) axis is a crucial step during tissue repair and regeneration. Despite the association reported recently of N6-methyladenosine (m6A) with regeneration, the mechanism underlying the regulation of m6A in A-P axis specification during regeneration remains unknown. Herein, we deciphered the m6A landscape at a single-base resolution at multiple time points during A-P axis regeneration and constructed the de novo transcriptome assembly of the Dugesia japonica planarian.

RESULTS

Immunofluorescence staining and comparative analysis revealed that m6A is widespread across the planarian and dynamically regulated during regeneration along the A-P axis, exhibiting a strong spatiotemporal feature. The resulting datasets of m6A-modified genes identified 80 anterior-specific genes and 13 posterior-specific genes, respectively. In addition, we showed that YTHDC1 serves as the primary m6A reader to be involved in the m6A-mediated specification of A-P axis during regeneration in Dugesia japonica planarian.

CONCLUSIONS

Our study provides an RNA epigenetic explanation for the specification of the A-P axis during tissue regeneration in planarian.

Ets-1 transcription factor regulates glial cell regeneration and function in planarians

Glia play multifaceted roles in nervous systems in response to injury. Depending on the species, extent of injury and glial cell type in question, glia can help or hinder the regeneration of neurons. Studying glia in the context of successful regeneration could reveal features of pro-regenerative glia that could be exploited for new human therapies. Planarian flatworms completely regenerate their nervous systems after injury - including glia - and thus provide a strong model system for exploring glia in the context of regeneration. Here, we report that planarian glia regenerate after neurons, and that neurons are required for correct glial numbers and localization during regeneration. We also identify the planarian transcription factor-encoding gene ets-1 as a key regulator of glial cell maintenance and regeneration. Using ets-1 (RNAi) to perturb glia, we show that glial loss is associated with altered neuronal gene expression, impeded animal movement and impaired nervous system architecture - particularly within the neuropil. Importantly, our work reveals the inter-relationships of glia and neurons in the context of robust neural regeneration.

BMP suppresses WNT to integrate patterning of orthogonal body axes in adult planarians

Adult regeneration restores patterning of orthogonal body axes after damage in a post-embryonic context. Planarians regenerate using distinct body-wide signals primarily regulating each axis dimension: anteroposterior Wnts, dorsoventral BMP, and mediolateral Wnt5 and Slit determinants. How regeneration can coordinate perpendicular tissue axes without symmetry-breaking embryonic events is not fully understood. Here, we report that the planarian dorsoventral regulator bmp4 suppresses the posterior determinant wnt1 to provide patterning input to the anteroposterior axis. Double-FISH identified distinct anteroposterior domains within dorsal midline muscle that express either bmp4 or wnt1. Homeostatic inhibition bmp4 and smad1 expanded the wnt1 expression anteriorly, while elevation of BMP signaling through nog1;nog2 RNAi reduced the wnt1 expression domain and elevated bmp4 expression. Homeostatic BMP signal perturbation broadly affected anteroposterior identity as measured by expression of posterior Wnt pathway factors, and caused mislocalization of AP-regionalized pharynx progenitors, without strongly affecting expression domains of anterior regulators. Additionally, wnt1 inhibition elevated bmp4 expression in the tip of the tail. Therefore, dorsal BMP signals and posterior wnt1 mutually antagonize for patterning the tail. Furthermore, homeostatic bmp4 RNAi caused medial expansion of the lateral determinant wnt5 and reduced expression of the medial regulator slit. By contrast, nog1;nog2 RNAi restricted wnt5 expression. Double RNAi of bmp4 and wnt5 resulted in lateral ectopic eye phenotypes, suggesting bmp4 acts upstream of wnt5 to pattern the mediolateral axis. These results indicate bmp4 controls dorsoventral information and also, through suppression of Wnt signals, influences anteroposterior and mediolateral identity. Based on related functions across vertebrates and Cnidarians, Wnt and BMP cross-regulation could form an ancient mechanism for coordinating orthogonal axis patterning.

Ultrafast distant wound response is essential for whole-body regeneration

Injury induces systemic responses, but their functions remain elusive. Mechanisms that can rapidly synchronize wound responses through long distances are also mostly unknown. Using planarian flatworms capable of whole-body regeneration, we report that injury induces extracellular signal-regulated kinase (Erk) activity waves to travel at a speed 10-100 times faster than those in other multicellular tissues. This ultrafast propagation requires longitudinal body-wall muscles, elongated cells forming dense parallel tracks running the length of the organism. The morphological properties of muscles allow them to act as superhighways for propagating and disseminating wound signals. Inhibiting Erk propagation prevents tissues distant to the wound from responding and blocks regeneration, which can be rescued by a second injury to distal tissues shortly after the first injury. Our findings provide a mechanism for long-range signal propagation in large, complex tissues to coordinate responses across cell types and highlight the function of feedback between spatially separated tissues during whole-body regeneration.

New technologies to study helminth development and host-parasite interactions

How parasites develop and survive, and how they stimulate or modulate host immune responses are important in understanding disease pathology and for the design of new control strategies. Microarray analysis and bulk RNA sequencing have provided a wealth of data on gene expression as parasites develop through different life-cycle stages and on host cell responses to infection. These techniques have enabled gene expression in the whole organism or host tissue to be detailed, but do not take account of the heterogeneity between cells of different types or developmental stages, nor the spatial organisation of these cells. Single-cell RNA-seq (scRNA-seq) adds a new dimension to studying parasite biology and host immunity by enabling gene profiling at the individual cell level. Here we review the application of scRNA-seq to establish gene expression cell atlases for multicellular helminths and to explore the expansion and molecular profile of individual host cell types involved in parasite immunity and tissue repair. Studying host-parasite interactions in vivo is challenging and we conclude this review by briefly discussing the applications of organoids (stem-cell derived mini-tissues) to examine host-parasite interactions at the local level, and as a potential system to study parasite development in vitro. Organoid technology and its applications have developed rapidly, and the elegant studies performed to date support the use of organoids as an alternative in vitro system for research on helminth parasites.

CRISPR/Cas9-based depletion of 16S ribosomal RNA improves library complexity of single-cell RNA-sequencing

Background

Single-cell RNA-sequencing (scRNA-seq) relies on PCR amplification to retrieve information from vanishingly small amounts of starting material. To selectively enrich mRNA from abundant non-polyadenylated transcripts, poly(A) selection is a key step during library preparation. However, some transcripts, such as mitochondrial genes, can escape this elimination and overwhelm libraries. Often, these transcripts are removed in silico, but whether physical depletion improves detection of rare transcripts in single cells is unclear.

Results

We find that a single 16S ribosomal RNA is widely enriched in planarian scRNA-seq datasets, independent of the library preparation method. To deplete this transcript from scRNA-seq libraries, we design 30 single-guide RNAs spanning its length. To evaluate the effects of depletion, we perform a side-by-side comparison of the effects of eliminating the 16S transcript and find a substantial increase in the number of genes detected per cell, coupled with virtually complete loss of the 16S RNA. Moreover, we systematically determine that library complexity increases with a limited number of PCR cycles following CRISPR treatment. When compared to in silico depletion of 16S, physically removing it reduces dropout rates, retrieves more clusters, and reveals more differentially-expressed genes.

Conclusions

Our results show that abundant transcripts reduce the retrieval of informative transcripts in scRNA-seq and distort the analysis. Physical removal of these contaminants enables the detection of rare transcripts at lower sequencing depth, and also outperforms in silico depletion. Importantly, this method can be easily customized to deplete any abundant transcript from scRNA-seq libraries.

Stem cells: The cell that does it all

How do animals replace all their worn-out cells to maintain their tissues? A new study shows that, in the cnidarian Hydractinia symbiolongicarpus, a single adult stem cell is sufficient to generate the entire repertoire of somatic and germ line cells.

Spatiotemporal transcriptomic atlas reveals the dynamic characteristics and key regulators of planarian regeneration

Whole-body regeneration of planarians is a natural wonder but how it occurs remains elusive. It requires coordinated responses from each cell in the remaining tissue with spatial awareness to regenerate new cells and missing body parts. While previous studies identified new genes essential to regeneration, a more efficient screening approach that can identify regeneration-associated genes in the spatial context is needed. Here, we present a comprehensive three-dimensional spatiotemporal transcriptomic landscape of planarian regeneration. We describe a pluripotent neoblast subtype, and show that depletion of its marker gene makes planarians more susceptible to sub-lethal radiation. Furthermore, we identified spatial gene expression modules essential for tissue development. Functional analysis of hub genes in spatial modules, such as plk1, shows their important roles in regeneration. Our three-dimensional transcriptomic atlas provides a powerful tool for deciphering regeneration and identifying homeostasis-related genes, and provides a publicly available online spatiotemporal analysis resource for planarian regeneration research.

Acoel single-cell atlas reveals expression dynamics and heterogeneity of adult pluripotent stem cells

Adult pluripotent stem cell (aPSC) populations underlie whole-body regeneration in many distantly-related animal lineages, but how the underlying cellular and molecular mechanisms compare across species is unknown. Here, we apply single-cell RNA sequencing to profile transcriptional cell states of the acoel worm Hofstenia miamia during postembryonic development and regeneration. We identify cell types shared across stages and their associated gene expression dynamics during regeneration. Functional studies confirm that the aPSCs, also known as neoblasts, are the source of differentiated cells and reveal transcription factors needed for differentiation. Subclustering of neoblasts recovers transcriptionally distinct subpopulations, the majority of which are likely specialized to differentiated lineages. One neoblast subset, showing enriched expression of the histone variant H3.3, appears to lack specialization. Altogether, the cell states identified in this study facilitate comparisons to other species and enable future studies of stem cell fate potentials.

m6 A promotes planarian regeneration

Regeneration is the regrowth of damaged tissues or organs, a vital process in response to damages from primitive organisms to higher mammals. Planarian possesses active whole-body regenerative capability owing to its vast reservoir of adult stem cells, neoblasts, providing an ideal model to delineate the underlying mechanisms for regeneration. RNA N⁶ -methyladenosine (m⁶ A) modification participates in many biological processes, including stem cell self-renewal and differentiation, in particular the regeneration of haematopoietic stem cells and axons. However, how m⁶ A controls regeneration at the whole-organism level remains largely unknown. Here, we demonstrate that the depletion of m⁶ A methyltransferase regulatory subunit wtap abolishes planarian regeneration, potentially through regulating genes related to cell-cell communication and cell cycle. Single-cell RNA-seq (scRNA-seq) analysis unveils that the wtap knockdown induces a unique type of neural progenitor-like cells (NP-like cells), characterized by specific expression of the cell-cell communication ligand grn. Intriguingly, the depletion of m⁶ A-modified transcripts grn, cdk9 or cdk7 partially rescues the defective regeneration of planarian caused by wtap knockdown. Overall, our study reveals an indispensable role of m⁶ A modification in regulating whole-organism regeneration.

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 hours 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 suggest new long-term cycling populations. Focused analyses of non-skeletal muscle and the endoderm identified transcriptional profiles of understudied cell types and subpopulations, including the pneumatic duct, individual intestinal smooth muscle layers, spatially distinct pericyte subpopulations, and homologs of recently discovered human best4+ enterocytes. The transcriptional regulators of these populations remain unknown, so we reconstructed gene expression trajectories to suggest candidates. To enable additional discoveries, we make this comprehensive transcriptional atlas of early zebrafish development available through our website, Daniocell.

Differentiation trajectories of the Hydra nervous system reveal transcriptional regulators of neuronal fate

The small freshwater cnidarian polyp Hydra vulgaris uses adult stem cells (interstitial stem cells) to continually replace neurons throughout its life. This feature, combined with the ability to image the entire nervous system (Badhiwala et al., 2021; Dupre & Yuste, 2017) and availability of gene knockdown techniques (Juliano, Reich, et al., 2014; Lohmann et al., 1999; Vogg et al., 2022), makes Hydra a tractable model for studying nervous system development and regeneration at the whole-organism level. In this study, we use single-cell RNA sequencing and trajectory inference to provide a comprehensive molecular description of the adult nervous system. This includes the most detailed transcriptional characterization of the adult Hydra nervous system to date. We identified eleven unique neuron subtypes together with the transcriptional changes that occur as the interstitial stem cells differentiate into each subtype. Towards the goal of building gene regulatory networks to describe Hydra neuron differentiation, we identified 48 transcription factors expressed specifically in the Hydra nervous system, including many that are conserved regulators of neurogenesis in bilaterians. We also performed ATAC-seq on sorted neurons to uncover previously unidentified putative regulatory regions near neuron-specific genes. Finally, we provide evidence to support the existence of transdifferentiation between mature neuron subtypes and we identify previously unknown transition states in these pathways. All together, we provide a comprehensive transcriptional description of an entire adult nervous system, including differentiation and transdifferentiation pathways, which provides a significant advance towards understanding mechanisms that underlie nervous system regeneration.

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.

scTransSort: Transformers for Intelligent Annotation of Cell Types by Gene Embeddings

Single-cell transcriptomics is rapidly advancing our understanding of the composition of complex tissues and biological cells, and single-cell RNA sequencing (scRNA-seq) holds great potential for identifying and characterizing the cell composition of complex tissues. Cell type identification by analyzing scRNA-seq data is mostly limited by time-consuming and irreproducible manual annotation. As scRNA-seq technology scales to thousands of cells per experiment, the exponential increase in the number of cell samples makes manual annotation more difficult. On the other hand, the sparsity of gene transcriptome data remains a major challenge. This paper applied the idea of the transformer to single-cell classification tasks based on scRNA-seq data. We propose scTransSort, a cell-type annotation method pretrained with single-cell transcriptomics data. The scTransSort incorporates a method of representing genes as gene expression embedding blocks to reduce the sparsity of data used for cell type identification and reduce the computational complexity. The feature of scTransSort is that its implementation of intelligent information extraction for unordered data, automatically extracting valid features of cell types without the need for manually labeled features and additional references. In experiments on cells from 35 human and 26 mouse tissues, scTransSort successfully elucidated its high accuracy and high performance for cell type identification, and demonstrated its own high robustness and generalization ability.

Inhibition of ATM kinase rescues planarian regeneration after lethal radiation

As stem cells divide, they acquire mutations that can be passed on to daughter cells. To mitigate potentially deleterious outcomes, cells activate the DNA damage response (DDR) network, which governs several cellular outcomes following DNA damage, including repairing DNA or undergoing apoptosis. At the helm of the DDR are three PI3-like kinases including Ataxia-Telangiectasia Mutated (ATM). We report here that knockdown of ATM in planarian flatworms enables stem cells to withstand lethal doses of radiation which would otherwise induce cell death. In this context, stem cells circumvent apoptosis, replicate their DNA, and recover function using homologous recombination-mediated DNA repair. Despite radiation exposure, atm knockdown animals survive long-term and regenerate new tissues. These effects occur independently of ATM's canonical downstream effector p53. Together, our results demonstrate that in planarians, ATM promotes radiation-induced apoptosis. This acute, ATM-dependent apoptosis is a key determinant of long-term animal survival. Our results suggest that inhibition of ATM in these organisms could, therefore, potentially favor cell survival after radiation without obvious effects on stem cell behavior.

Ultrafast and long-range coordination of wound responses is essential for whole-body regeneration

Injury induces systemic, global responses whose functions remain elusive. In addition, mechanisms that rapidly synchronize wound responses through long distances across the organismal scale are mostly unknown. Using planarians, which have extreme regenerative ability, we report that injury induces Erk activity to travel in a wave-like manner at an unexpected speed (∼1 mm/h), 10-100 times faster than those measured in other multicellular tissues. This ultrafast signal propagation requires longitudinal body-wall muscles, elongated cells forming dense parallel tracks running the length of the organism. Combining experiments and computational models, we show that the morphological properties of muscles allow them to minimize the number of slow intercellular signaling steps and act as bidirectional superhighways for propagating wound signals and instructing responses in other cell types. Inhibiting Erk propagation prevents cells distant to the wound from responding and blocks regeneration, which can be rescued by a second injury to distal tissues within a narrow time window after the first injury. These results suggest that rapid responses in uninjured tissues far from wounds are essential for regeneration. Our findings provide a mechanism for long-range signal propagation in large and complex tissues to coordinate cellular responses across diverse cell types, and highlights the function of feedback between spatially separated tissues during whole-body regeneration.

Beyond the Genetic Code: A Tissue Code?

The genetic code determines how the precise amino acid sequence of proteins is specified by genomic information in cells. But what specifies the precise histologic organization of cells in plant and animal tissues is unclear. We now hypothesize that another code, the tissue code , exists at an even higher level of complexity which determines how tissue organization is dynamically maintained. Accordingly, we modeled spatial and temporal asymmetries of cell division and established that five simple mathematical laws ("the tissue code") convey a set of biological rules that maintain the specific organization and continuous self-renewal dynamics of cells in tissues. These laws might even help us understand wound healing, and how tissue disorganization leads to birth defects and tissue pathology like cancer.

Schmidtea mediterranea as a Model Organism to Study the Molecular Background of Human Motile Ciliopathies

Cilia and flagella are evolutionarily conserved organelles that form protrusions on the surface of many growth-arrested or differentiated eukaryotic cells. Due to the structural and functional differences, cilia can be roughly classified as motile and non-motile (primary). Genetically determined dysfunction of motile cilia is the basis of primary ciliary dyskinesia (PCD), a heterogeneous ciliopathy affecting respiratory airways, fertility, and laterality. In the face of the still incomplete knowledge of PCD genetics and phenotype-genotype relations in PCD and the spectrum of PCD-like diseases, a continuous search for new causative genes is required. The use of model organisms has been a great part of the advances in understanding molecular mechanisms and the genetic basis of human diseases; the PCD spectrum is not different in this respect. The planarian model (Schmidtea mediterranea) has been intensely used to study regeneration processes, and-in the context of cilia-their evolution, assembly, and role in cell signaling. However, relatively little attention has been paid to the use of this simple and accessible model for studying the genetics of PCD and related diseases. The recent rapid development of the available planarian databases with detailed genomic and functional annotations prompted us to review the potential of the S. mediterranea model for studying human motile ciliopathies.

Emergent dynamics of adult stem cell lineages from single nucleus and single cell RNA-Seq of Drosophila testes

Proper differentiation of sperm from germline stem cells, essential for production of the next generation, requires dramatic changes in gene expression that drive remodeling of almost all cellular components, from chromatin to organelles to cell shape itself. Here, we provide a single nucleus and single cell RNA-seq resource covering all of spermatogenesis in Drosophila starting from in-depth analysis of adult testis single nucleus RNA-seq (snRNA-seq) data from the Fly Cell Atlas (FCA) study. With over 44,000 nuclei and 6000 cells analyzed, the data provide identification of rare cell types, mapping of intermediate steps in differentiation, and the potential to identify new factors impacting fertility or controlling differentiation of germline and supporting somatic cells. We justify assignment of key germline and somatic cell types using combinations of known markers, in situ hybridization, and analysis of extant protein traps. Comparison of single cell and single nucleus datasets proved particularly revealing of dynamic developmental transitions in germline differentiation. To complement the web-based portals for data analysis hosted by the FCA, we provide datasets compatible with commonly used software such as Seurat and Monocle. The foundation provided here will enable communities studying spermatogenesis to interrogate the datasets to identify candidate genes to test for function in vivo.

Planarian LDB and SSDP proteins scaffold transcriptional complexes for regeneration and patterning

Sequence-specific transcription factors often function as components of large regulatory complexes. LIM-domain binding protein (LDB) and single-stranded DNA-binding protein (SSDP) function as core scaffolds of transcriptional complexes in animals and plants. Little is known about potential partners and functions for LDB/SSDP complexes in the context of tissue regeneration. In this work, we find that planarian LDB1 and SSDP2 promote tissue regeneration, with a particular function in mediolateral polarity reestablishment. We find that LDB1 and SSDP2 interact with one another and with characterized planarian LIM-HD proteins Arrowhead, Islet1, and Lhx1/5-1. SSDP2 and LDB1 also function with islet1 in polarity reestablishment and with lhx1/5-1 in serotonergic neuron maturation. Finally, we show new roles for LDB1 and SSDP2 in regulating gene expression in the planarian intestine and parenchyma; these functions may be LIM-HD-independent. Together, our work provides insight into LDB/SSDP complexes in a highly regenerative organism. Further, our work provides a strong starting point for identifying and characterizing potential binding partners of LDB1 and SSDP2 and for exploring roles for these proteins in diverse aspects of planarian physiology.

Djck1α Is Required for Proper Regeneration and Maintenance of the Medial Tissues in Planarians

CK1α (Casein kinase 1α) is a member of the casein kinase 1(CK1) family that is involved in diverse cellular processes, but its functions remain unclear in stem cell development. Freshwater planarians are capable of whole-body regeneration, making it a classic model for the study of regeneration, tissue homeostasis, and polarity in vivo. To investigate the roles of CK1α in regeneration and homeostasis progress, we characterize a homolog of CK1α from planarian Dugesia japonica. We find that Djck1α, which shows an enriched expression pattern in the nascent tissues, is widely expressed especially in the medial regions of planarians. Knockdown of CK1α by RNAi presents a thicker body due to dorsal hyperplasia, along with defects in the medial tissues including nerve proliferation, missing epidermis, intestine disturbance, and hyper-proliferation during the progression of regeneration and homeostasis. Moreover, we find that the ck1α RNAi animals exhibit expansion of the midline marker slit. The eye deficiency induced by slit RNAi can be rescued by ck1α and slit double RNAi. These results suggest that ck1α is required for the medial tissue regeneration and maintenance in planarian Dugesia japonica by regulating the expression of slit, which helps to further investigate the regulation of planarian mediolateral axis.

Construction of a cross-species cell landscape at single-cell level

Individual cells are basic units of life. Despite extensive efforts to characterize the cellular heterogeneity of different organisms, cross-species comparisons of landscape dynamics have not been achieved. Here, we applied single-cell RNA sequencing (scRNA-seq) to map organism-level cell landscapes at multiple life stages for mice, zebrafish and Drosophila. By integrating the comprehensive dataset of > 2.6 million single cells, we constructed a cross-species cell landscape and identified signatures and common pathways that changed throughout the life span. We identified structural inflammation and mitochondrial dysfunction as the most common hallmarks of organism aging, and found that pharmacological activation of mitochondrial metabolism alleviated aging phenotypes in mice. The cross-species cell landscape with other published datasets were stored in an integrated online portal-Cell Landscape. Our work provides a valuable resource for studying lineage development, maturation and aging.

Planarians to schistosomes: an overview of flatworm cell-types and regulators

Schistosomiasis remains a major neglected tropical disease that afflicts over 200 million people globally. Schistosomes, the aetiological agent of schistosomiasis, are parasitic flatworms that propagate between molluscan and mammalian hosts. Inside the mammalian host, schistosomes rapidly grow over 100-fold in size and develop into a sexually mature male or female that thrives in the bloodstream for several decades. Recent work has identified schistosome stem cells as the source that drives parasite transmission, reproduction and longevity. Moreover, studies have begun to uncover molecular programmes deployed by stem cells that are essential for tissue development and maintenance, parasite survival and immune evasion. Such programmes are reminiscent of neoblast-driven development and regeneration of planarians, the free-living flatworm relative of schistosomes. Over the last few decades, research in planarians has employed modern functional genomic tools that significantly enhanced our understanding of stem cell-driven animal development and regeneration. In this review, we take a broad stroke overview of major flatworm organ systems at the cellular and molecular levels. We summarize recent advances on genetic regulators that play critical roles in differentiation and maintenance of flatworm cell types. Finally, we provide perspectives on how investigation of basic parasite biology is critical to discovering new approaches to battle schistosomiasis.

BMP suppresses WNT to integrate patterning of orthogonal body axes in adult planarians

Adult regeneration restores patterning of orthogonal body axes after damage in a post-embryonic context. Planarians regenerate using distinct body-wide signals primarily regulating each axis dimension: anteroposterior Wnts, dorsoventral BMP, and mediolateral Wnt5 and Slit determinants. How regeneration can consistently form perpendicular tissue axes without symmetry-breaking embryonic events is unknown, and could either occur using fully independent, or alternatively, integrated signals defining each dimension. Here, we report that the planarian dorsoventral regulator bmp4 suppresses the posterior determinant wnt1 to pattern the anteroposterior axis. Double-FISH identified distinct anteroposterior domains within dorsal midline muscle that express either bmp4 or wnt1 . Homeostatic inhibition bmp4 and smad1 expanded the wnt1 expression anteriorly, while elevation of BMP signaling through nog1;nog2 RNAi reduced the wnt1 expression domain. BMP signal perturbation broadly affected anteroposterior identity as measured by expression of posterior Wnt pathway factors, without affecting head regionalization. Therefore, dorsal BMP signals broadly limit posterior identity. Furthermore, bmp4 RNAi caused medial expansion of the lateral determinant wnt5 and reduced expression of the medial regulator slit . Double RNAi of bmp4 and wnt5 resulted in lateral ectopic eye phenotypes, suggesting bmp4 acts upstream of wnt5 to pattern the mediolateral axis. Therefore, bmp4 acts at the top of a patterning hierarchy both to control dorsoventral information and also, through suppression of Wnt signals, to regulate anteroposterior and mediolateral identity. These results reveal that adult pattern formation involves integration of signals controlling individual orthogonal axes.

Author Summary

Systems that coordinate long-range communication across axes are likely critical for enabling tissue restoration in regenerative animals. While individual axis pathways have been identified, there is not yet an understanding of how signal integration allows repatterning across 3-dimensions. Here, we report an unanticipated linkage between anteroposterior, dorsoventral, and mediolateral systems in planarians through BMP signaling. We find that dorsally expressed BMP restricts posterior and lateral identity by suppressing distinct Wnt signals in adult planarians. These results demonstrate that orthogonal axis information is not fully independent and suggest a potentially ancient role of integrated axis patterning in generating stable 3-dimensional adult forms.

MorphoFeatures for unsupervised exploration of cell types, tissues, and organs in volume electron microscopy

Electron microscopy (EM) provides a uniquely detailed view of cellular morphology, including organelles and fine subcellular ultrastructure. While the acquisition and (semi-)automatic segmentation of multicellular EM volumes are now becoming routine, large-scale analysis remains severely limited by the lack of generally applicable pipelines for automatic extraction of comprehensive morphological descriptors. Here, we present a novel unsupervised method for learning cellular morphology features directly from 3D EM data: a neural network delivers a representation of cells by shape and ultrastructure. Applied to the full volume of an entire three-segmented worm of the annelid Platynereis dumerilii, it yields a visually consistent grouping of cells supported by specific gene expression profiles. Integration of features across spatial neighbours can retrieve tissues and organs, revealing, for example, a detailed organisation of the animal foregut. We envision that the unbiased nature of the proposed morphological descriptors will enable rapid exploration of very different biological questions in large EM volumes, greatly increasing the impact of these invaluable, but costly resources.

Wnt signaling in whole-body regeneration

Regeneration abilities are widespread among animals and select species can restore any body parts removed by wounds that sever the major body axes. This capability of whole-body regeneration as exemplified in flatworm planarians, Acoels, and Cnidarians involves initial responses to injury, the assessment of wound site polarization, determination of missing tissue and programming of blastema fate, and patterned outgrowth to restore axis content and proportionality. Wnt signaling drives many shared and conserved aspects of the biology of whole-body regeneration in the planarian species Schmidtea mediterranea and Dugesia japonica, in the Acoel Hofstenia miamia, and in Cnidarians Hydra and Nematostella. These overlapping mechanisms suggest whole-body regeneration might be an ancestral property across diverse animal taxa.

An RNA/DNA-Based Flow Cytometry Approach for Isolating Slow-Cycling Stem Cells

Flow cytometry methods for sorting specific populations of cells based on fluorescence or physical properties have been a widely used technique for decades. Flow cytometry has been particularly vital to the study of planarians, which remain refractory to transgenic transformation, as it has provided a work-around solution for studying stem cell biology and lineage relationships in the context of regeneration. Many flow cytometry applications have been published in planarians, beginning with broad Hoechst-based strategies for isolating cycling stem cells and progressing to more function-based approaches involving vital dyes and surface antibodies. In this protocol, we look to build on the classic DNA-labeling Hoechst staining strategy by adding pyronin Y staining to label RNA. While Hoechst labeling alone allows for the isolation of stem cells in the S/G2/M phases of the cell cycle, heterogeneity within the population of stem cells with 2 C DNA content is not resolved. By considering RNA levels, this protocol can further divide this population of stem cells into two groups: G1 stem cells with relatively high RNA content and a slow-cycling population with low RNA content, which we call RNA^low stem cells. In addition, we provide instruction for combining this RNA/DNA flow cytometry protocol with EdU labeling experiments and describe an optional step for incorporating immunostaining prior to cell sorting (in this case with the pluripotency marker TSPAN-1). This protocol adds a new staining strategy and examples of combinatorial flow cytometry approaches to the repertoire of flow cytometry techniques for studying planarian stem cells.

A Wnt11 and Dishevelled signaling pathway acts prior to injury to control wound polarization for the onset of planarian regeneration

Regeneration is initiated by wounding, but it is unclear how injury-induced signals precisely convey the identity of the tissues requiring replacement. In the planarian Schmidtea mediterranea, the first event in head regeneration is the asymmetric activation of the Wnt inhibitor notum in longitudinal body-wall muscle cells, preferentially at anterior-facing versus posterior-facing wound sites. However, the mechanism driving this early symmetry-breaking event is unknown. We identify a noncanonical Wnt11 and Dishevelled pathway regulating notum polarization, which opposes injury-induced notum-activating Wnt/β-catenin signals and regulates muscle orientation. Using expression analysis and experiments to define a critical time of action, we demonstrate that Wnt11 and Dishevelled signals act prior to injury and in a growth-dependent manner to orient the polarization of notum induced by wounding. In turn, injury-induced notum dictates polarization used in the next round of regeneration. These results identify a self-reinforcing feedback system driving the polarization of blastema outgrowth and indicate that regeneration uses pre-existing tissue information to determine the outcome of wound-induced signals.

Microtubule-associated protein 1B is implicated in stem cell commitment and nervous system regeneration in planarians

Microtubule-associated 1B (MAP1B) proteins are expressed at the nervous system level where they control cytoskeleton activity and regulate neurotransmitter release. Here, we report about the identification of a planarian MAP1B factor (DjMap1B) that is enriched in cephalic ganglia and longitudinal nerve cords but not in neoblasts, the plentiful population of adult stem cells present in planarians, thanks to which these animals can continuously cell turnover and regenerate any lost body parts. DjMap1B knockdown induces morphological anomalies in the nervous system and affects neoblast commitment. Our data put forward a correlation between a MAP1B factor and stem cells and suggest a function of the nervous system in non-cell autonomous control of planarian stem cells.

Embryonic origins of adult pluripotent stem cells

Although adult pluripotent stem cells (aPSCs) are found in many animal lineages, mechanisms for their formation during embryogenesis are unknown. Here, we leveraged Hofstenia miamia, a regenerative worm that possesses collectively pluripotent aPSCs called neoblasts and produces manipulable embryos. Lineage tracing and functional experiments revealed that one pair of blastomeres gives rise to cells that resemble neoblasts in distribution, behavior, and gene expression. In Hofstenia, aPSCs include transcriptionally distinct subpopulations that express markers associated with differentiated tissues; our data suggest that despite their heterogeneity, aPSCs are derived from one lineage, not from multiple tissue-specific lineages during development. Next, we combined single-cell transcriptome profiling across development with neoblast cell-lineage tracing and identified a molecular trajectory for neoblast formation that includes transcription factors Hes, FoxO, and Tbx. This identification of a cellular mechanism and molecular trajectory for aPSC formation opens the door for in vivo studies of aPSC regulation and evolution.