A PAK family kinase and the Hippo/Yorkie pathway modulate WNT signaling to functionally integrate body axes during regeneration.

Successful regeneration of missing tissues requires seamless integration of positional information along the body axes. Planarians, which regenerate from almost any injury, use conserved, developmentally important signaling pathways to pattern the body axes. However, the molecular mechanisms which facilitate cross talk between these signaling pathways to integrate positional information remain poorly understood. Here, we report a p21-activated kinase (smed-pak1) which functionally integrates the anterior-posterior (AP) and the medio-lateral (ML) axes. pak1 inhibits WNT/ß-catenin signaling along the AP axis and, functions synergistically with the ß-catenin-independent WNT signaling of the ML axis. Furthermore, this functional integration is dependent on warts and merlin-the components of the Hippo/Yorkie (YKI) pathway. Hippo/YKI pathway is a critical regulator of body size in flies and mice, but our data suggest the pathway regulates body axes patterning in planarians. Our study provides a signaling network integrating positional information which can mediate coordinated growth and patterning during planarian regeneration.

Conserved Genes in Highly Regenerative Metazoans Are Associated with Planarian Regeneration.

Metazoan species depict a wide spectrum of regeneration ability which calls into question the evolutionary origins of the underlying processes. Since species with high regeneration ability are widely distributed throughout metazoans, there is a possibility that the metazoan ancestor had an underlying common molecular mechanism. Early metazoans like sponges possess high regenerative ability, but, due to the large differences they have with Cnidaria and Bilateria regarding symmetry and neuronal systems, it can be inferred that this regenerative ability is different. We hypothesized that the last common ancestor of Cnidaria and Bilateria possessed remarkable regenerative ability which was lost during evolution. We separated Cnidaria and Bilateria into three classes possessing whole-body regenerating, high regenerative ability, and low regenerative ability. Using a multiway BLAST and gene phylogeny approach, we identified genes conserved in whole-body regenerating species and lost in low regenerative ability species and labeled them Cnidaria and Bilaterian regeneration genes. Through transcription factor analysis, we identified that Cnidaria and Bilaterian regeneration genes were associated with an overabundance of homeodomain regulatory elements. RNA interference of Cnidaria and Bilaterian regeneration genes resulted in loss of regeneration phenotype for HRJDa, HRJDb, DUF21, DISP3, and ARMR genes. We observed that DUF21 knockdown was highly lethal in the early stages of regeneration indicating a potential role in wound response. Also, HRJDa, HRJDb, DISP3, and ARMR knockdown showed loss of regeneration phenotype after second amputation. The results strongly correlate with their respective RNA-seq profiles. We propose that Cnidaria and Bilaterian regeneration genes play a major role in regeneration across highly regenerative Cnidaria and Bilateria.

Protocol for culturing and functionally manipulating planarian neoblasts using SiR-DNA-based flow cytometry.

Neoblasts are the only cells capable of proliferation in planarians. The traditional flow cytometry protocol using Hoechst inhibits the cell cycle. Here, we present a protocol for culturing and functionally manipulating planarian neoblasts using SiR-DNA-based flow cytometry. We describe steps for cell dissociation and staining, flow cytometry, and cell collection and culture. We then detail procedures for Nanoluciferase mRNA transfection. This protocol facilitates further investigations into the pluripotency and regeneration mechanisms within neoblasts. For complete details on the use and execution of this protocol, please refer to Lei et al.1.

Flatworm Transcriptomes Reveal Widespread Parasitism by Histophagous Ciliates.

Unicellular ciliates like Tetrahymena are best known as free-living bacteriovores, but many species are facultative or obligate parasites. These "histophages" feed on the tissues of hosts ranging from planarian flatworms to commercially important fish and the larvae of imperiled freshwater mussels. Here, we developed a novel bioinformatics pipeline incorporating the nonstandard ciliate genetic code and used it to search for Ciliophora sequences in 34 publicly available Platyhelminthes EST libraries. From 2,615,036 screened ESTs, we identified nearly 6,000 high-confidence ciliate transcripts, supporting parasitism of seven additional flatworm species. We also cultured and identified Tetrahymena from nine terrestrial and freshwater planarians, including invasive earthworm predators from the genus Bipalium and the widely studied regeneration models Dugesia japonica and Schmidtea mediterranea. A co-phylogenetic reconstruction provides strong evidence for the coevolution of histophagous Ciliophora with their Platyhelminthes hosts. We further report the antiprotozoal aminoglycoside paromomycin expels Tetrahymena from S. mediterranea, providing new opportunities to investigate the effects of this relationship on planarian biology. Together, our findings raise the possibility that invasive flatworms constitute a novel dispersal mechanism for Tetrahymena parasites and position the Platyhelminthes as an ideal model phylum for studying the ecology and evolution of histophagous ciliates.

Mutational profile of the regenerative process and de novo genome assembly of the planarian Schmidtea polychroa.

Planarians are organisms with a unique capacity to regenerate any part of their body. New tissues are generated in a process that requires many swift cell divisions. How costly is this process to an animal in terms of mutational load remains unknown. Using whole genome sequencing, we defined the mutational profile of the process of regeneration in the planarian species Schmidtea polychroa. We assembled de novo the genome of S. polychroa and analyzed mutations in animals that have undergone regeneration. We observed a threefold increase in the number of mutations and an altered mutational spectrum. High allele frequencies of subclonal mutations in regenerated animals suggested that most of the cells in the regenerated animal were descendants of a small number of stem cells with high expansion potential. We provide, for the first time, the draft genome assembly of S. polychroa, an estimation of the germline mutation rate for a planarian species and the mutational spectrum of the regeneration process of a living organism.

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.

Starvation resistance in planarians: multiple strategies to get a thrifty phenotype.

Starvation resistance is a life-saving mechanism for many organisms facing food availability fluctuation in the natural environment. Different strategies have been episodically identified for some model organisms, the first of which was the ability to suppress metabolic rate. Among the identified strategies, the ability of planarians to shrink their body under fasting conditions and revert the process after feeding (the growth-degrowth process) represents a fascinating mechanism to face long periods of fasting. The growth-degrowth process is strictly related to the capability of planarians to continuously maintain tissue homeostasis and body proportions even in challenging conditions, thanks to the presence of a population of pluripotent stem cells. Here, we take advantage of several previous studies describing the growth-degrowth process and of recent progress in the understanding of planarian homeostasis mechanisms, to identify tissue-selective transcriptional downregulation as a driving strategy for the development of a thrifty phenotype, and the p53 transcription factor as a player in adjusting tissue homeostasis in accordance with food availability.

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.

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.

A simplified and rapid in situ hybridization protocol for planarians.

Whole-mount in situ hybridization is a critical technique for analyzing gene expression in planarians. While robust in situ protocols have been developed, these protocols are laborious, making them challenging to incorporate in an academic setting, reducing throughput and increasing time to results. Here, the authors systematically tested modifications to all phases of the protocol with the goal of eliminating steps and reducing time without impacting quality. This modified protocol allows for whole-mount colorimetric in situ hybridization and multicolor fluorescence in situ hybridization to be completed in two days with a significant reduction in steps and hands-on processing time.

Evolutionary dynamics of whole-body regeneration across planarian flatworms.

Regenerative abilities vary dramatically across animals. Even amongst planarian flatworms, well-known for complete regeneration from tiny body fragments, some species have restricted regeneration abilities while others are almost entirely regeneration incompetent. Here, we assemble a diverse live collection of 40 planarian species to probe the evolution of head regeneration in the group. Combining quantification of species-specific head-regeneration abilities with a comprehensive transcriptome-based phylogeny reconstruction, we show multiple independent transitions between robust whole-body regeneration and restricted regeneration in freshwater species. RNA-mediated genetic interference inhibition of canonical Wnt signalling in RNA-mediated genetic interference-sensitive species bypassed all head-regeneration defects, suggesting that the Wnt pathway is linked to the emergence of planarian regeneration defects. Our finding that Wnt signalling has multiple roles in the reproductive system of the model species Schmidtea mediterranea raises the possibility that a trade-off between egg-laying, asexual reproduction by fission/regeneration and Wnt signalling drives regenerative trait evolution. Although quantitative comparisons of Wnt signalling levels, yolk content and reproductive strategy across our species collection remained inconclusive, they revealed divergent Wnt signalling roles in the reproductive system of planarians. Altogether, our study establishes planarians as a model taxon for comparative regeneration research and presents a framework for the mechanistic evolution of regenerative abilities.

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.

Neurotoxicity of Glyphosate to Planarian Dugesia japonica.

As one of the most widely used herbicides in agricultural industry, the residues of glyphosate (GLY) are frequent environmental pollutants. Freshwater planarian Dugesia japonica has been developed as a model for neurotoxicology. In this study, the effects of GLY on locomotion and feeding behavior, as well as neuroenzyme activities and mRNA expressions of D. japonica were determined. Additionally, histochemical localization was executed to explore the damage to the central nervous system (CNS) of planarians stressed by GLY. The results showed that the locomotor velocity, ingestion rate and the neuroenzyme activity were inhibited and the gene expressions were altered. Also, histo-architecture injury to CNS of planarians upon GLY exposure in a time-dependent manner was observed. Collectively, our results indicate that GLY can cause neurotoxicity to freshwater planarians representing as reduction in locomotor velocity and feeding rate by disturbing the neurotransmission systems and damaging the structure of CNS.

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.

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.

Isolation of planarian viable cells using fluorescence-activated cell sorting for advancing single-cell transcriptome analysis.

Preparing viable single cells is critical for conducting single-cell RNA sequencing (scRNA-seq) because the presence of ambient RNA from dead or damaged cells can interfere with data analysis. Here, we developed a method for isolating viable single cells from adult planarian bodies using fluorescence-activated cell sorting (FACS). This method was then applied to both adult pluripotent stem cells (aPSCs) and differentiating/differentiated cells. Initially, we employed a violet instead of ultraviolet (UV) laser to excite Hoechst 33342 to reduce cellular damage. After optimization of cell staining conditions and FACS compensation, we generated FACS profiles similar to those created using a previous method that employed a UV laser. Despite successfully obtaining high-quality RNA sequencing data for aPSCs, non-aPSCs produced low-quality RNA reads (i.e., <60% of cells possessing barcoding mRNAs). Subsequently, we identified an effective FACS gating condition that excluded low-quality cells and tissue debris without staining. This non-staining isolation strategy not only reduced post-dissociation time but also enabled high-quality scRNA-seq results for all cell types (i.e., >80%). Taken together, these findings imply that the non-staining FACS strategy may be beneficial for isolating viable cells not only from planarians but also from other organisms and tissues for scRNA-seq studies.

aristaless-like homeobox-3 is wound induced and promotes a low-Wnt environment required for planarian head regeneration.

The planarian Schmidtea mediterranea is a well-established model of adult regeneration, which is dependent on a large population of adult stem cells called neoblasts. Upon amputation, planarians undergo transcriptional wounding programs and coordinated stem cell proliferation to give rise to missing tissues. Interestingly, the Wnt signaling pathway is key to guiding what tissues are regenerated, yet less known are the transcriptional regulators that ensure proper activation and timing of signaling pathway components. Here, we have identified an aristaless-like homeobox transcription factor, alx-3, that is enriched in a population of putative neural-fated progenitor cells at homeostasis, and is also upregulated in stem cells and muscle cells at anterior-facing wounds upon amputation. Knockdown of alx-3 results in failure of head regeneration and patterning defects in amputated tail fragments. alx-3 is required for the expression of several early wound-induced genes, including the Wnt inhibitor notum, which is required to establish anterior polarity during regeneration. Together, these findings reveal a role for alx-3 as an early wound-response transcriptional regulator in both muscle cells and stem cells that is required for anterior regeneration by promoting a low-Wnt environment.

Genome-Wide Analysis of Planarian piRNAs.

In planarian flatworms, the piRNA pathway is operated by three PIWI proteins, termed SMEDWI-1, SMEDWI-2, and SMEDWI-3 (SMEDWI = Schmidtea mediterranea PIWI). The interplay between these three PIWI proteins and their associated small noncoding RNAs, termed piRNAs, fuels the outstanding regenerative abilities of planarians, enables tissue homeostasis, and, ultimately, ensures animal survival. As the molecular targets of PIWI proteins are determined by the sequences of their co-bound piRNAs, it is imperative to identify these sequences by next-generation sequencing applications. Following sequencing, the genomic targets and the regulatory potential of the isolated piRNA populations need to be uncovered. To that end, here we present a bioinformatics analysis pipeline for processing and systematic characterization of planarian piRNAs. The pipeline includes steps for the removal of PCR duplicates based on unique molecular identifier (UMI) sequences, and it accounts for piRNA multimapping to different loci in the genome. Importantly, our protocol also includes a fully automated pipeline that is freely available at GitHub. Together with the piRNA isolation and library preparation protocol (see accompanying chapter), the presented computational pipeline enables researchers to explore the functional role of the piRNA pathway in flatworm biology.