m6A is required for resolving progenitor identity during planarian stem cell differentiation.

Regeneration and tissue homeostasis require accurate production of missing cell lineages. Cell production is driven by changes to gene expression, which is shaped by multiple layers of regulation. Here, we find that the ubiquitous mRNA base-modification, m6A, is required for proper cell fate choice and cellular maturation in planarian stem cells (neoblasts). We mapped m6A-enriched regions in 7,600 planarian genes and found that perturbation of the m6A pathway resulted in progressive deterioration of tissues and death. Using single-cell RNA sequencing of >20,000 cells following perturbation of the m6A pathway, we identified an increase in expression of noncanonical histone variants, and that inhibition of the pathway resulted in accumulation of undifferentiated cells throughout the animal in an abnormal transcriptional state. Analysis of >1,000 planarian gene expression datasets revealed that the inhibition of the chromatin modifying complex NuRD had almost indistinguishable consequences, unraveling an unappreciated link between m6A and chromatin modifications. Our findings reveal that m6A is critical for planarian stem cell homeostasis and gene regulation in tissue maintenance and regeneration.

Molecular characterization of the immediate wound response of the solitary ascidian Polycarpa mytiligera.

BACKGROUND: Injury response is key to successful regeneration. Yet, transcriptome analyses of injury response were performed only on a handful of regenerative organisms. Here, we studied the injury response of the solitary ascidian Polycarpa mytiligera, an emerging model system, capable of regenerating any body part. We used the siphon as a model for studying transcriptional changes following injury, and identified genes that were activated in the initial 24 hours post amputation (hpa). RESULTS: Highly conserved genes, such as bone morphogenetic protein-1 (BMP1), growth hormone secretagogue receptor (GHSR) and IL-17, were upregulated by 12 hpa, yet their expression was sustained only in non-regenerating tissue fragments. We optimized fluorescent in situ hybridization, and found that the majority of BMP1+ cells were localized to the rigid tunic that covers the animal. This highlights the importance of this tissue, particularly during injury response. BMP1 was overexpressed following injuries to other body regions, suggesting that it was a part of a common injury-induced program. CONCLUSION: Our study suggests that, initially, specific injury-induced genes were upregulated in P. mytiligera organs, yet, later, a unique transcriptional profile was observed only in regenerating tissues. These findings highlight the importance of studying diverse regenerating and non-regenerating organisms for complete understanding of regeneration.

Comparative transcriptomics across the prokaryotic tree of life.

Whole-transcriptome sequencing studies from recent years revealed an unexpected complexity in transcriptomes of bacteria and archaea, including abundant non-coding RNAs, cis-antisense transcription and regulatory untranslated regions (UTRs). Understanding the functional relevance of the plethora of non-coding RNAs in a given organism is challenging, especially since some of these RNAs were attributed to 'transcriptional noise'. To allow the search for conserved transcriptomic elements we produced comparative transcriptome maps for multiple species across the microbial tree of life. These transcriptome maps are detailed in annotations, comparable by gene families, and BLAST-searchable by user provided sequences. Our transcriptome collection includes 18 model organisms spanning 10 phyla/subphyla of bacteria and archaea that were sequenced using standardized RNA-seq methods. The utility of the comparative approach, as implemented in our web server, is demonstrated by highlighting genes with exceptionally long 5'UTRs across species, which correspond to many known riboswitches and further suggest novel putative regulatory elements. Our study provides a standardized reference transcriptome to major clinically and environmentally important microbial phyla. The viewer is available at http://exploration.weizmann.ac.il/TCOL, setting a framework for comparative studies of the microbial non-coding genome.

Genomic island variability facilitates Prochlorococcus-virus coexistence.

Prochlorococcus cyanobacteria are extremely abundant in the oceans, as are the viruses that infect them. How hosts and viruses coexist in nature remains unclear, although the presence of both susceptible and resistant cells may allow this coexistence. Combined whole-genome sequencing and PCR screening technology now enables us to investigate the effect of resistance on genome evolution and the genomic mechanisms behind the long-term coexistence of Prochlorococcus and their viruses. Here we present a genome analysis of 77 substrains selected for resistance to ten viruses, revealing mutations primarily in non-conserved, horizontally transferred genes that localize to a single hypervariable genomic island. Mutations affected viral attachment to the cell surface and imposed a fitness cost to the host, manifested by significantly lower growth rates or a previously unknown mechanism of more rapid infection by other viruses. The mutant genes are generally uncommon in nature yet some carry polymorphisms matching those found experimentally. These data are empirical evidence indicating that viral-attachment genes are preferentially located in genomic islands and that viruses are a selective pressure enhancing the diversity of both island genes and island gene content. This diversity emerges as a genomic mechanism that reduces the effective host population size for infection by a given virus, thus facilitating long-term coexistence between viruses and their hosts in nature.

Transcriptome-wide mapping of 5-methylcytidine RNA modifications in bacteria, archaea, and yeast reveals m5C within archaeal mRNAs.

The presence of 5-methylcytidine (m(5)C) in tRNA and rRNA molecules of a wide variety of organisms was first observed more than 40 years ago. However, detection of this modification was limited to specific, abundant, RNA species, due to the usage of low-throughput methods. To obtain a high resolution, systematic, and comprehensive transcriptome-wide overview of m(5)C across the three domains of life, we used bisulfite treatment on total RNA from both gram positive (B. subtilis) and gram negative (E. coli) bacteria, an archaeon (S. solfataricus) and a eukaryote (S. cerevisiae), followed by massively parallel sequencing. We were able to recover most previously documented m(5)C sites on rRNA in the four organisms, and identified several novel sites in yeast and archaeal rRNAs. Our analyses also allowed quantification of methylated m(5)C positions in 64 tRNAs in yeast and archaea, revealing stoichiometric differences between the methylation patterns of these organisms. Molecules of tRNAs in which m(5)C was absent were also discovered. Intriguingly, we detected m(5)C sites within archaeal mRNAs, and identified a consensus motif of AUCGANGU that directs methylation in S. solfataricus. Our results, which were validated using m(5)C-specific RNA immunoprecipitation, provide the first evidence for mRNA modifications in archaea, suggesting that this mode of post-transcriptional regulation extends beyond the eukaryotic domain.

Correlated occurrence and bypass of frame-shifting insertion-deletions (InDels) to give functional proteins.

Short insertions and deletions (InDels) comprise an important part of the natural mutational repertoire. InDels are, however, highly deleterious, primarily because two-thirds result in frame-shifts. Bypass through slippage over homonucleotide repeats by transcriptional and/or translational infidelity is known to occur sporadically. However, the overall frequency of bypass and its relation to sequence composition remain unclear. Intriguingly, the occurrence of InDels and the bypass of frame-shifts are mechanistically related - occurring through slippage over repeats by DNA or RNA polymerases, or by the ribosome, respectively. Here, we show that the frequency of frame-shifting InDels, and the frequency by which they are bypassed to give full-length, functional proteins, are indeed highly correlated. Using a laboratory genetic drift, we have exhaustively mapped all InDels that occurred within a single gene. We thus compared the naive InDel repertoire that results from DNA polymerase slippage to the frame-shifting InDels tolerated following selection to maintain protein function. We found that InDels repeatedly occurred, and were bypassed, within homonucleotide repeats of 3-8 bases. The longer the repeat, the higher was the frequency of InDels formation, and the more frequent was their bypass. Besides an expected 8A repeat, other types of repeats, including short ones, and G and C repeats, were bypassed. Although obtained in vitro, our results indicate a direct link between the genetic occurrence of InDels and their phenotypic rescue, thus suggesting a potential role for frame-shifting InDels as bridging evolutionary intermediates.

A vast collection of microbial genes that are toxic to bacteria.

In the process of clone-based genome sequencing, initial assemblies frequently contain cloning gaps that can be resolved using cloning-independent methods, but the reason for their occurrence is largely unknown. By analyzing 9,328,693 sequencing clones from 393 microbial genomes, we systematically mapped more than 15,000 genes residing in cloning gaps and experimentally showed that their expression products are toxic to the Escherichia coli host. A subset of these toxic sequences was further evaluated through a series of functional assays exploring the mechanisms of their toxicity. Among these genes, our assays revealed novel toxins and restriction enzymes, and new classes of small, non-coding toxic RNAs that reproducibly inhibit E. coli growth. Further analyses also revealed abundant, short, toxic DNA fragments that were predicted to suppress E. coli growth by interacting with the replication initiator DnaA. Our results show that cloning gaps, once considered the result of technical problems, actually serve as a rich source for the discovery of biotechnologically valuable functions, and suggest new modes of antimicrobial interventions.

The single-nucleotide resolution transcriptome of Pseudomonas aeruginosa grown in body temperature.

One of the hallmarks of opportunistic pathogens is their ability to adjust and respond to a wide range of environmental and host-associated conditions. The human pathogen Pseudomonas aeruginosa has an ability to thrive in a variety of hosts and cause a range of acute and chronic infections in individuals with impaired host defenses or cystic fibrosis. Here we report an in-depth transcriptional profiling of this organism when grown at host-related temperatures. Using RNA-seq of samples from P. aeruginosa grown at 28°C and 37°C we detected genes preferentially expressed at the body temperature of mammalian hosts, suggesting that they play a role during infection. These temperature-induced genes included the type III secretion system (T3SS) genes and effectors, as well as the genes responsible for phenazines biosynthesis. Using genome-wide transcription start site (TSS) mapping by RNA-seq we were able to accurately define the promoters and cis-acting RNA elements of many genes, and uncovered new genes and previously unrecognized non-coding RNAs directly controlled by the LasR quorum sensing regulator. Overall we identified 165 small RNAs and over 380 cis-antisense RNAs, some of which predicted to perform regulatory functions, and found that non-coding RNAs are preferentially localized in pathogenicity islands and horizontally transferred regions. Our work identifies regulatory features of P. aeruginosa genes whose products play a role in environmental adaption during infection and provides a reference transcriptional landscape for this pathogen.

RNA-seq analysis of small RNPs in Trypanosoma brucei reveals a rich repertoire of non-coding RNAs.

The discovery of a plethora of small non-coding RNAs (ncRNAs) has fundamentally changed our understanding of how genes are regulated. In this study, we employed the power of deep sequencing of RNA (RNA-seq) to examine the repertoire of ncRNAs present in small ribonucleoprotein particles (RNPs) of Trypanosoma brucei, an important protozoan parasite. We identified new C/D and H/ACA small nucleolar RNAs (snoRNAs), as well as tens of putative novel non-coding RNAs; several of these are processed from trans-spliced and polyadenylated transcripts. The RNA-seq analysis provided information on the relative abundance of the RNAs, and their 5'- and 3'-termini. The study demonstrated that three highly abundant snoRNAs are involved in rRNA processing and highlight the unique trypanosome-specific repertoire of these RNAs. Novel RNAs were studied using in situ hybridization, association in RNP complexes, and 'RNA walk' to detect interaction with their target RNAs. Finally, we showed that the abundance of certain ncRNAs varies between the two stages of the parasite, suggesting that ncRNAs may contribute to gene regulation during the complex parasite's life cycle. This is the first study to provide a whole-genome analysis of the large repertoire of small RNPs in trypanosomes.

Self-targeting by CRISPR: gene regulation or autoimmunity?

The recently discovered prokaryotic immune system known as CRISPR (clustered regularly interspaced short palindromic repeats) is based on small RNAs ('spacers') that restrict phage and plasmid infection. It has been hypothesized that CRISPRs can also regulate self gene expression by utilizing spacers that target self genes. By analyzing CRISPRs from 330 organisms we found that one in every 250 spacers is self-targeting, and that such self-targeting occurs in 18% of all CRISPR-bearing organisms. However, complete lack of conservation across species, combined with abundance of degraded repeats near self-targeting spacers, suggests that self-targeting is a form of autoimmunity rather than a regulatory mechanism. We propose that accidental incorporation of self nucleic acids by CRISPR can incur an autoimmune fitness cost, and this could explain the abundance of degraded CRISPR systems across prokaryotes.

A single-base resolution map of an archaeal transcriptome.

Organisms of the third domain of life, the Archaea, share molecular characteristics both with Bacteria and Eukarya. These organisms attract scientific attention as research models for regulation and evolution of processes such as transcription, translation, and RNA processing. We have reconstructed the primary transcriptome of Sulfolobus solfataricus P2, one of the most widely studied model archaeal organisms. Analysis of 625 million bases of sequenced cDNAs yielded a single-base-pair resolution map of transcription start sites and operon structures for more than 1000 transcriptional units. The analysis led to the discovery of 310 expressed noncoding RNAs, with an extensive expression of overlapping cis-antisense transcripts to a level unprecedented in any bacteria or archaea but resembling that of eukaryotes. As opposed to bacterial transcripts, most Sulfolobus transcripts completely lack 5'-UTR sequences, suggesting that mRNA/ncRNA interactions differ between Bacteria and Archaea. The data also reveal internal hotspots for transcript cleavage linked to RNA degradation and predict sequence motifs that promote RNA destabilization. This study highlights transcriptome sequencing as a key tool for understanding the mechanisms and extent of RNA-based regulation in Bacteria and Archaea.

Validation of two ribosomal RNA removal methods for microbial metatranscriptomics.

The predominance of rRNAs in the transcriptome is a major technical challenge in sequence-based analysis of cDNAs from microbial isolates and communities. Several approaches have been applied to deplete rRNAs from (meta)transcriptomes, but no systematic investigation of potential biases introduced by any of these approaches has been reported. Here we validated the effectiveness and fidelity of the two most commonly used approaches, subtractive hybridization and exonuclease digestion, as well as combinations of these treatments, on two synthetic five-microorganism metatranscriptomes using massively parallel sequencing. We found that the effectiveness of rRNA removal was a function of community composition and RNA integrity for these treatments. Subtractive hybridization alone introduced the least bias in relative transcript abundance, whereas exonuclease and in particular combined treatments greatly compromised mRNA abundance fidelity. Illumina sequencing itself also can compromise quantitative data analysis by introducing a G+C bias between runs.

Comparative transcriptomics of pathogenic and non-pathogenic Listeria species.

Listeria monocytogenes is a human, food-borne pathogen. Genomic comparisons between L. monocytogenes and Listeria innocua, a closely related non-pathogenic species, were pivotal in the identification of protein-coding genes essential for virulence. However, no comprehensive comparison has focused on the non-coding genome. We used strand-specific cDNA sequencing to produce genome-wide transcription start site maps for both organisms, and developed a publicly available integrative browser to visualize and analyze both transcriptomes in different growth conditions and genetic backgrounds. Our data revealed conservation across most transcripts, but significant divergence between the species in a subset of non-coding RNAs. In L. monocytogenes, we identified 113 small RNAs (33 novel) and 70 antisense RNAs (53 novel), significantly increasing the repertoire of ncRNAs in this species. Remarkably, we identified a class of long antisense transcripts (lasRNAs) that overlap one gene while also serving as the 5' UTR of the adjacent divergent gene. Experimental evidence suggests that lasRNAs transcription inhibits expression of one operon while activating the expression of another. Such a lasRNA/operon structure, that we named 'excludon', might represent a novel form of regulation in bacteria.

Methods for cell isolation and analysis of the highly regenerative tunicate Polycarpa mytiligera.

Background: Polycarpa mytiligera is the only molecularly characterized solitary ascidian capable of regenerating all organs and tissue types. The cellular basis for regeneration in P. mytiligera is largely unknown, and methods for isolating live cells from this species for functional analyses are unavailable. Results: Here, we developed a method for isolating live cells from P. mytiligera, overcoming major experimental challenges, including the dissociation of its thick body wall and native cellular autofluorescence. We demonstrated the applicability of our approach for tissue dissociation and cell analysis using three flow cytometry platforms, and by using broadly used non-species-specific cell labeling reagents. In addition to live cell isolation, proof-of-concept experiments showed that this approach was compatible with gene expression analysis of RNA extracted from the isolated cells, and with ex vivo analysis of phagocytosis. Conclusion: We presented efficient methods for cell purification from a highly regenerative ascidian, which could be transferable to diversity of non-model marine organisms. The ability to purify live cells will promote future studies of cell function in P. mytiligera regeneration.

Stemness Activity Underlying Whole Brain Regeneration in a Basal Chordate.

Understanding how neurons regenerate following injury remains a central challenge in regenerative medicine. Adult mammals have a very limited ability to regenerate new neurons in the central nervous system (CNS). In contrast, the basal chordate Polycarpa mytiligera can regenerate its entire CNS within seven days of complete removal. Transcriptome sequencing, cellular labeling, and proliferation in vivo essays revealed that CNS regeneration is mediated by a newly formed neural progeny and the activation of neurodevelopmental pathways that are associated with enhanced stem-cell activity. Analyzing the expression of 239 activated pathways enabled a quantitative understanding of gene-set enrichment patterns at key regeneration stages. The molecular and cellular mechanisms controlling the regenerative ability that this study reveals can be used to develop innovative approaches to enhancing neurogenesis in closely-related chordate species, including humans.

Planarian stem cells specify fate yet retain potency during the cell cycle.

Planarian whole-body regeneration is enabled by stem cells called neoblasts. At least some neoblasts are individually pluripotent. Neoblasts are also heterogeneous, with subpopulations of specialized neoblasts having different specified fates. Fate specification in neoblasts is regulated by fate-specific transcription factor (FSTF) expression. Here, we find that FSTF expression is common in neoblast S/G2/M cell-cycle phases but less common in G1. We find that specialized neoblasts can divide to produce progeny with asymmetric cell fates, suggesting that they could retain pluripotency. Furthermore, no known neoblast class was present in all neoblast colonies, suggesting that pluripotency is not the exclusive property of any known class. We tested this possibility with single-cell transplantations, which indicate that at least some specialized neoblasts are likely clonogenic. On the basis of these findings, we propose a model for neoblast pluripotency in which neoblasts can undergo specialization during the cell cycle without loss of potency.

foxF-1 Controls Specification of Non-body Wall Muscle and Phagocytic Cells in Planarians.

Planarians are flatworms capable of regenerating any missing body part in a process requiring stem cells and positional information. Muscle is a major source of planarian positional information and consists of several types of fibers with distinct regulatory roles in regeneration. The transcriptional regulatory programs used to specify different muscle fibers are poorly characterized. Using single-cell RNA sequencing, we define the transcriptomes of planarian dorsal-ventral muscle (DVM), intestinal muscle (IM), and pharynx muscle. This analysis identifies foxF-1, which encodes a broadly conserved Fox-family transcription factor, as a master transcriptional regulator of all non-body wall muscle. The transcription factors encoded by nk4 and gata4/5/6-2 specify two different subsets of DVM, lateral and medial, respectively, whereas gata4/5/6-3 specifies IM. These muscle types all express planarian patterning genes. Both lateral and medial DVM are required for medial-lateral patterning in regeneration, whereas medial DVM and IM have a role in maintaining and regenerating intestine morphology. In addition to the role in muscle, foxF-1 is required for the specification of multiple cell types with transcriptome similarities, including high expression levels of cathepsin genes. These cells include pigment cells, glia, and several other cells with unknown function. cathepsin+ cells phagocytose E. coli, suggesting these are phagocytic cells. In conclusion, we describe a regulatory program for planarian muscle cell subsets and phagocytic cells, both driven by foxF-1. FoxF proteins specify different mesoderm-derived tissues in other organisms, suggesting that FoxF regulates formation of an ancient and broadly conserved subset of mesoderm derivatives in the Bilateria.

Cell type transcriptome atlas for the planarian Schmidtea mediterranea.

The transcriptome of a cell dictates its unique cell type biology. We used single-cell RNA sequencing to determine the transcriptomes for essentially every cell type of a complete animal: the regenerative planarian Schmidtea mediterranea. Planarians contain a diverse array of cell types, possess lineage progenitors for differentiated cells (including pluripotent stem cells), and constitutively express positional information, making them ideal for this undertaking. We generated data for 66,783 cells, defining transcriptomes for known and many previously unknown planarian cell types and for putative transition states between stem and differentiated cells. We also uncovered regionally expressed genes in muscle, which harbors positional information. Identifying the transcriptomes for potentially all cell types for many organisms should be readily attainable and represents a powerful approach to metazoan biology.

Planarian Epidermal Stem Cells Respond to Positional Cues to Promote Cell-Type Diversity.

Successful regeneration requires that progenitors of different lineages form the appropriate missing cell types. However, simply generating lineages is not enough. Cells produced by a particular lineage often have distinct functions depending on their position within the organism. How this occurs in regeneration is largely unexplored. In planarian regeneration, new cells arise from a proliferative cell population (neoblasts). We used the planarian epidermal lineage to study how the location of adult progenitor cells results in their acquisition of distinct functional identities. Single-cell RNA sequencing of epidermal progenitors revealed the emergence of distinct spatial identities as early in the lineage as the epidermal neoblasts, with further pre-patterning occurring in their post-mitotic migratory progeny. Establishment of dorsal-ventral epidermal identities and functions, in response to BMP signaling, required neoblasts. Our work identified positional signals that activate regionalized transcriptional programs in the stem cell population and subsequently promote cell-type diversity in the epidermis.