miR-430 regulates zygotic mRNA during zebrafish embryogenesis.

BACKGROUND: Early embryonic developmental programs are guided by the coordinated interplay between maternally inherited and zygotically manufactured RNAs and proteins. Although these processes happen concomitantly and affecting gene function during this period is bound to affect both pools of mRNAs, it has been challenging to study their expression dynamics separately. RESULTS: By employing SLAM-seq, a nascent mRNA labeling transcriptomic approach, in a developmental time series we observe that over half of the early zebrafish embryo transcriptome consists of maternal-zygotic genes, emphasizing their pivotal role in early embryogenesis. We provide an hourly resolution of de novo transcriptional activation events and follow nascent mRNA trajectories, finding that most de novo transcriptional events are stable throughout this period. Additionally, by blocking microRNA-430 function, a key post transcriptional regulator during zebrafish embryogenesis, we directly show that it destabilizes hundreds of de novo transcribed mRNAs from pure zygotic as well as maternal-zygotic genes. This unveils a novel miR-430 function during embryogenesis, fine-tuning zygotic gene expression. CONCLUSION: These insights into zebrafish early embryo transcriptome dynamics emphasize the significance of post-transcriptional regulators in zygotic genome activation. The findings pave the way for future investigations into the coordinated interplay between transcriptional and post-transcriptional landscapes required for the establishment of animal cell identities and functions.

Identification and Localization of Cell Types in the Mouse Olfactory Bulb Using Slide-SeqV2.

Spatial transcriptomics maps RNA molecules to the location in a tissue where they are expressed. Here we document the use of Slide-SeqV2 to visualize gene expression in the mouse olfactory bulb (OB). This approach relies on spatially identified beads to locate and quantify individual transcripts. The expression profiles associated with the beads are used to identify and localize individual cell types in an unbiased manner. We demonstrate the various cell types and subtypes with distinct spatial locations in the olfactory bulb that are identified using Slide-SeqV2.

Cell-type profiling of the sympathetic nervous system using spatial transcriptomics and spatial mapping of mRNA.

BACKGROUND: The molecular identification of neural progenitor cell populations that connect to establish the sympathetic nervous system (SNS) remains unclear. This is due to technical limitations in the acquisition and spatial mapping of molecular information to tissue architecture. RESULTS: To address this, we applied Slide-seq spatial transcriptomics to intact fresh frozen chick trunk tissue transversely cryo-sectioned at the developmental stage prior to SNS formation. In parallel, we performed age- and location-matched single cell (sc) RNA-seq and 10× Genomics Visium to inform our analysis. Downstream bioinformatic analyses led to the unique molecular identification of neural progenitor cells within the peripheral sympathetic ganglia (SG) and spinal cord preganglionic neurons (PGNs). We then successfully applied the HiPlex RNAscope fluorescence in situ hybridization and multispectral confocal microscopy to visualize 12 gene targets in stage-, age- and location-matched chick trunk tissue sections. CONCLUSIONS: Together, these data demonstrate a robust strategy to acquire and integrate single cell and spatial transcriptomic information, resulting in improved resolution of molecular heterogeneities in complex neural tissue architectures. Successful application of this strategy to the developing SNS provides a roadmap for functional studies of neural connectivity and platform to address complex questions in neural development and regeneration.

Off-target piRNA gene silencing in Drosophila melanogaster rescued by a transposable element insertion.

Transposable elements (TE) are selfish genetic elements that can cause harmful mutations. In Drosophila, it has been estimated that half of all spontaneous visible marker phenotypes are mutations caused by TE insertions. Several factors likely limit the accumulation of exponentially amplifying TEs within genomes. First, synergistic interactions between TEs that amplify their harm with increasing copy number are proposed to limit TE copy number. However, the nature of this synergy is poorly understood. Second, because of the harm posed by TEs, eukaryotes have evolved systems of small RNA-based genome defense to limit transposition. However, as in all immune systems, there is a cost of autoimmunity and small RNA-based systems that silence TEs can inadvertently silence genes flanking TE insertions. In a screen for essential meiotic genes in Drosophila melanogaster, a truncated Doc retrotransposon within a neighboring gene was found to trigger the germline silencing of ald, the Drosophila Mps1 homolog, a gene essential for proper chromosome segregation in meiosis. A subsequent screen for suppressors of this silencing identified a new insertion of a Hobo DNA transposon in the same neighboring gene. Here we describe how the original Doc insertion triggers flanking piRNA biogenesis and local gene silencing. We show that this local gene silencing occurs in cis and is dependent on deadlock, a component of the Rhino-Deadlock-Cutoff (RDC) complex, to trigger dual-strand piRNA biogenesis at TE insertions. We further show how the additional Hobo insertion leads to de-silencing by reducing flanking piRNA biogenesis triggered by the original Doc insertion. These results support a model of TE-mediated gene silencing by piRNA biogenesis in cis that depends on local determinants of transcription. This may explain complex patterns of off-target gene silencing triggered by TEs within populations and in the laboratory. It also provides a mechanism of sign epistasis among TE insertions, illuminates the complex nature of their interactions and supports a model in which off-target gene silencing shapes the evolution of the RDC complex.

Tumor-initiating stem cell shapes its microenvironment into an immunosuppressive barrier and pro-tumorigenic niche.

Tumor-initiating stem cells (TSCs) are critical for drug resistance and immune escape. However, the mutual regulations between TSC and tumor microenvironment (TME) remain unclear. Using DNA-label retaining, single-cell RNA sequencing (scRNA-seq), and other approaches, we investigated intestinal adenoma in response to chemoradiotherapy (CRT), thus identifying therapy-resistant TSCs (TrTSCs). We find bidirectional crosstalk between TSCs and TME using CellPhoneDB analysis. An intriguing finding is that TSCs shape TME into a landscape that favors TSCs for immunosuppression and propagation. Using adenoma-organoid co-cultures, niche-cell depletion, and lineaging tracing, we characterize a functional role of cyclooxygenase-2 (Cox-2)-dependent signaling, predominantly occurring between tumor-associated monocytes and macrophages (TAMMs) and TrTSCs. We show that TAMMs promote TrTSC proliferation through prostaglandin E2 (PGE2)-PTGER4(EP4) signaling, which enhances ß-catenin activity via AKT phosphorylation. Thus, our study shows that the bidirectional crosstalk between TrTSC and TME results in a pro-tumorigenic and immunosuppressive contexture.

Multiple Niche Compartments Orchestrate Stepwise Germline Stem Cell Progeny Differentiation.

The niche controls stem cell self-renewal and progenitor differentiation for maintaining adult tissue homeostasis in various organisms. However, it remains unclear whether the niche is compartmentalized to control stem cell self-renewal and stepwise progeny differentiation. In the Drosophila ovary, inner germarial sheath (IGS) cells form a niche for controlling germline stem cell (GSC) progeny differentiation. In this study, we have identified four IGS subpopulations, which form linearly arranged niche compartments for controlling GSC maintenance and multi-step progeny differentiation. Single-cell analysis of the adult ovary has identified four IGS subpopulations (IGS1-IGS4), the identities and cellular locations of which have been further confirmed by fluorescent in situ hybridization. IGS1 and IGS2 physically interact with GSCs and mitotic cysts to control GSC maintenance and cyst formation, respectively, whereas IGS3 and IGS4 physically interact with 16-cell cysts to regulate meiosis, oocyte development, and cyst morphological change. Finally, one follicle cell progenitor population has also been transcriptionally defined for facilitating future studies on follicle stem cell regulation. Therefore, this study has structurally revealed that the niche is organized into multiple compartments for orchestrating stepwise adult stem cell development and has also provided useful resources and tools for further functional characterization of the niche in the future.

Drosophila YBX1 homolog YPS promotes ovarian germ line stem cell development by preferentially recognizing 5-methylcytosine RNAs.

5-Methylcytosine (m5C) is a RNA modification that exists in tRNAs and rRNAs and was recently found in mRNAs. Although it has been suggested to regulate diverse biological functions, whether m5C RNA modification influences adult stem cell development remains undetermined. In this study, we show that Ypsilon schachtel (YPS), a homolog of human Y box binding protein 1 (YBX1), promotes germ line stem cell (GSC) maintenance, proliferation, and differentiation in the Drosophila ovary by preferentially binding to m5C-containing RNAs. YPS is genetically demonstrated to function intrinsically for GSC maintenance, proliferation, and progeny differentiation in the Drosophila ovary, and human YBX1 can functionally replace YPS to support normal GSC development. Highly conserved cold-shock domains (CSDs) of YPS and YBX1 preferentially bind to m5C RNA in vitro. Moreover, YPS also preferentially binds to m5C-containing RNAs, including mRNAs, in germ cells. The crystal structure of the YBX1 CSD-RNA complex reveals that both hydrophobic stacking and hydrogen bonds are critical for m5C binding. Overexpression of RNA-binding-defective YPS and YBX1 proteins disrupts GSC development. Taken together, our findings show that m5C RNA modification plays an important role in adult stem cell development.

Multisite Evaluation of Next-Generation Methods for Small RNA Quantification.

Small RNAs (smRNAs) are important regulators of many biologic processes and are now most frequently characterized using Illumina sequencing. However, although standard RNA sequencing library preparation has become routine in most sequencing facilities, smRNA sequencing library preparation has historically been challenging because of high input requirements, laborious protocols involving gel purifications, inability to automate, and a lack of benchmarking standards. Additionally, studies have suggested that many of these methods are nonlinear and do not accurately reflect the amounts of smRNAs in vivo. Recently, a number of new kits have become available that permit lower input amounts and less laborious, gel-free protocol options. Several of these new kits claim to reduce RNA ligase-dependent sequence bias through novel adapter modifications and to lessen adapter-dimer contamination in the resulting libraries. With the increasing number of smRNA kits available, understanding the relative strengths of each method is crucial for appropriate experimental design. In this study, we systematically compared 9 commercially available smRNA library preparation kits as well as NanoString probe hybridization across multiple study sites. Although several of the new methodologies do reduce the amount of artificially over- and underrepresented microRNAs (miRNAs), we observed that none of the methods was able to remove all of the bias in the library preparation. Identical samples prepared with different methods show highly varied levels of different miRNAs. Even so, many methods excelled in ease of use, lower input requirement, fraction of usable reads, and reproducibility across sites. These differences may help users select the most appropriate methods for their specific question of interest.

scRNA-Seq reveals distinct stem cell populations that drive hair cell regeneration after loss of Fgf and Notch signaling.

Loss of sensory hair cells leads to deafness and balance deficiencies. In contrast to mammalian hair cells, zebrafish ear and lateral line hair cells regenerate from poorly characterized support cells. Equally ill-defined is the gene regulatory network underlying the progression of support cells to differentiated hair cells. scRNA-Seq of lateral line organs uncovered five different support cell types, including quiescent and activated stem cells. Ordering of support cells along a developmental trajectory identified self-renewing cells and genes required for hair cell differentiation. scRNA-Seq analyses of fgf3 mutants, in which hair cell regeneration is increased, demonstrates that Fgf and Notch signaling inhibit proliferation of support cells in parallel by inhibiting Wnt signaling. Our scRNA-Seq analyses set the foundation for mechanistic studies of sensory organ regeneration and is crucial for identifying factors to trigger hair cell production in mammals. The data is searchable and publicly accessible via a web-based interface.

A Population of Navigator Neurons Is Essential for Olfactory Map Formation during the Critical Period.

In the developing brain, heightened plasticity during the critical period enables the proper formation of neural circuits. Here, we identify the "navigator" neurons, a group of perinatally born olfactory sensory neurons, as playing an essential role in establishing the olfactory map during the critical period. The navigator axons project circuitously in the olfactory bulb and traverse multiple glomeruli before terminating in perspective glomeruli. These neurons undergo a phase of exuberant axon growth and exhibit a shortened lifespan. Single-cell transcriptome analyses reveal distinct molecular signatures for the navigators. Extending their lifespan prolongs the period of exuberant growth and perturbs axon convergence. Conversely, a genetic ablation experiment indicates that, despite postnatal neurogenesis, only the navigators are endowed with the ability to establish a convergent map. The presence and the proper removal of the navigator neurons are both required to establish tight axon convergence into the glomeruli.

Author Correction: Suppression of m6A reader Ythdf2 promotes hematopoietic stem cell expansion.

In the initial published version of this article, there was an inadvertent omission from the Acknowledgements that this work was supported by Stowers Institute for Medical Research (SIMR-1004) and NIH National Cancer Institute grant to University of Kansas Cancer Center (P30 CA168524). This omission does not affect the description of the results or the conclusions of this work.

Suppression of m6A reader Ythdf2 promotes hematopoietic stem cell expansion.

Transplantation of hematopoietic stem cells (HSCs) from human umbilical cord blood (hUCB) holds great promise for treating a broad spectrum of hematological disorders including cancer. However, the limited number of HSCs in a single hUCB unit restricts its widespread use. Although extensive efforts have led to multiple methods for ex vivo expansion of human HSCs by targeting single molecules or pathways, it remains unknown whether it is possible to simultaneously manipulate the large number of targets essential for stem cell self-renewal. Recent studies indicate that N6-methyladenosine (m6A) modulates the expression of a group of mRNAs critical for stem cell-fate determination by influencing their stability. Among several m6A readers, YTHDF2 is recognized as promoting targeted mRNA decay. However, the physiological functions of YTHDF2 in adult stem cells are unknown. Here we show that following the conditional knockout (KO) of mouse Ythdf2 the numbers of functional HSC were increased without skewing lineage differentiation or leading to hematopoietic malignancies. Furthermore, knockdown (KD) of human YTHDF2 led to more than a 10-fold increase in the ex vivo expansion of hUCB HSCs, a fivefold increase in colony-forming units (CFUs), and more than an eightfold increase in functional hUCB HSCs in the secondary serial of a limiting dilution transplantation assay. Mapping of m6A in RNAs from mouse hematopoietic stem and progenitor cells (HSPCs) as well as from hUCB HSCs revealed its enrichment in mRNAs encoding transcription factors critical for stem cell self-renewal. These m6A-marked mRNAs were recognized by Ythdf2 and underwent decay. In Ythdf2 KO HSPCs and YTHDF2 KD hUCB HSCs, these mRNAs were stabilized, facilitating HSC expansion. Knocking down one of YTHDF2's key targets, Tal1 mRNA, partially rescued the phenotype. Our study provides the first demonstration of the function of YTHDF2 in adult stem cell maintenance and identifies its important role in regulating HSC ex vivo expansion by regulating the stability of multiple mRNAs critical for HSC self-renewal, thus identifying potential for future clinical applications.

Rapid Low-Cost Assembly of the Drosophila melanogaster Reference Genome Using Low-Coverage, Long-Read Sequencing.

Accurate and comprehensive characterization of genetic variation is essential for deciphering the genetic basis of diseases and other phenotypes. A vast amount of genetic variation stems from large-scale sequence changes arising from the duplication, deletion, inversion, and translocation of sequences. In the past 10 years, high-throughput short reads have greatly expanded our ability to assay sequence variation due to single nucleotide polymorphisms. However, a recent de novo assembly of a second Drosophila melanogaster reference genome has revealed that short read genotyping methods miss hundreds of structural variants, including those affecting phenotypes. While genomes assembled using high-coverage long reads can achieve high levels of contiguity and completeness, concerns about cost, errors, and low yield have limited widespread adoption of such sequencing approaches. Here we resequenced the reference strain of D. melanogaster (ISO1) on a single Oxford Nanopore MinION flow cell run for 24 hr. Using only reads longer than 1 kb or with at least 30x coverage, we assembled a highly contiguous de novo genome. The addition of inexpensive paired reads and subsequent scaffolding using an optical map technology achieved an assembly with completeness and contiguity comparable to the D. melanogaster reference assembly. Comparison of our assembly to the reference assembly of ISO1 uncovered a number of structural variants (SVs), including novel LTR transposable element insertions and duplications affecting genes with developmental, behavioral, and metabolic functions. Collectively, these SVs provide a snapshot of the dynamics of genome evolution. Furthermore, our assembly and comparison to the D. melanogaster reference genome demonstrates that high-quality de novo assembly of reference genomes and comprehensive variant discovery using such assemblies are now possible by a single lab for under $1,000 (USD).

Prospectively Isolated Tetraspanin+ Neoblasts Are Adult Pluripotent Stem Cells Underlying Planaria Regeneration.

Proliferating cells known as neoblasts include pluripotent stem cells (PSCs) that sustain tissue homeostasis and regeneration of lost body parts in planarians. However, the lack of markers to prospectively identify and isolate these adult PSCs has significantly hampered their characterization. We used single-cell RNA sequencing (scRNA-seq) and single-cell transplantation to address this long-standing issue. Large-scale scRNA-seq of sorted neoblasts unveiled a novel subtype of neoblast (Nb2) characterized by high levels of PIWI-1 mRNA and protein and marked by a conserved cell-surface protein-coding gene, tetraspanin 1 (tspan-1). tspan-1-positive cells survived sub-lethal irradiation, underwent clonal expansion to repopulate whole animals, and when purified with an anti-TSPAN-1 antibody, rescued the viability of lethally irradiated animals after single-cell transplantation. The first prospective isolation of an adult PSC bridges a conceptual dichotomy between functionally and molecularly defined neoblasts, shedding light on mechanisms governing in vivo pluripotency and a source of regeneration in animals. VIDEO ABSTRACT.

Retinoid-Sensitive Epigenetic Regulation of the Hoxb Cluster Maintains Normal Hematopoiesis and Inhibits Leukemogenesis.

Hox genes modulate the properties of hematopoietic stem cells (HSCs) and reacquired Hox expression in progenitors contributes to leukemogenesis. Here, our transcriptome and DNA methylome analyses revealed that Hoxb cluster and retinoid signaling genes are predominantly enriched in LT-HSCs, and this coordinate regulation of Hoxb expression is mediated by a retinoid-dependent cis-regulatory element, distal element RARE (DERARE). Deletion of the DERARE reduced Hoxb expression, resulting in changes to many downstream signaling pathways (e.g., non-canonical Wnt signaling) and loss of HSC self-renewal and reconstitution capacity. DNA methyltransferases mediate DNA methylation on the DERARE, leading to reduced Hoxb cluster expression. Acute myeloid leukemia patients with DNMT3A mutations exhibit DERARE hypomethylation, elevated HOXB expression, and adverse outcomes. CRISPR-Cas9-mediated specific DNA methylation at DERARE attenuated HOXB expression and alleviated leukemogenesis. Collectively, these findings demonstrate pivotal roles for retinoid signaling and the DERARE in maintaining HSCs and preventing leukemogenesis by coordinate regulation of Hoxb genes.

Single-cell transcriptome analysis of avian neural crest migration reveals signatures of invasion and molecular transitions.

Neural crest cells migrate throughout the embryo, but how cells move in a directed and collective manner has remained unclear. Here, we perform the first single-cell transcriptome analysis of cranial neural crest cell migration at three progressive stages in chick and identify and establish hierarchical relationships between cell position and time-specific transcriptional signatures. We determine a novel transcriptional signature of the most invasive neural crest Trailblazer cells that is consistent during migration and enriched for approximately 900 genes. Knockdown of several Trailblazer genes shows significant but modest changes to total distance migrated. However, in vivo expression analysis by RNAscope and immunohistochemistry reveals some salt and pepper patterns that include strong individual Trailblazer gene expression in cells within other subregions of the migratory stream. These data provide new insights into the molecular diversity and dynamics within a neural crest cell migratory stream that underlie complex directed and collective cell behaviors.

Aneuploidy as a cause of impaired chromatin silencing and mating-type specification in budding yeast.

Aneuploidy and epigenetic alterations have long been associated with carcinogenesis, but it was unknown whether aneuploidy could disrupt the epigenetic states required for cellular differentiation. In this study, we found that ~3% of random aneuploid karyotypes in yeast disrupt the stable inheritance of silenced chromatin during cell proliferation. Karyotype analysis revealed that this phenotype was significantly correlated with gains of chromosomes III and X. Chromosome X disomy alone was sufficient to disrupt chromatin silencing and yeast mating-type identity as indicated by a lack of growth response to pheromone. The silencing defect was not limited to cryptic mating type loci and was associated with broad changes in histone modifications and chromatin localization of Sir2 histone deacetylase. The chromatin-silencing defect of disome X can be partially recapitulated by an extra copy of several genes on chromosome X. These results suggest that aneuploidy can directly cause epigenetic instability and disrupt cellular differentiation.

Sequencing of the TBX6 Gene in Families with Familial Idiopathic Scoliosis.

STUDY DESIGN: A hypothesis-driven study was conducted in a familial cohort to determine the potential association between variants within the TBX6 gene and Familial Idiopathic Scoliosis (FIS). OBJECTIVE: To determine if variants within exons of the TBX6 gene segregate with the FIS phenotype within a sample of families with FIS. SUMMARY OF BACKGROUND DATA: Idiopathic Scoliosis (IS) is a structural curvature of the spine whose underlying genetic etiology has not been established. IS has been reported to occur at a higher rate than expected in family members of individuals with congenital scoliosis (CS), suggesting that the two diseases might have a shared etiology. The TBX6 gene on chromosome 16p, essential to somite development, has been associated with CS in a Chinese population. Previous studies have identified linkage to this locus in families with FIS, and specifically with rs8060511, located in an intron of the TBX6 gene. METHODS: Parent-offspring trios from 11 families (13 trios, 42 individuals) with FIS were selected for Sanger sequencing of the TBX6 gene. Trios were selected from a large population of families with FIS in which a genome-wide scan had resulted in linkage to 16p. RESULTS: Sequencing analyses of the subset of families resulted in the identification of five coding variants. Three of the five variants were novel; the remaining two variants were previously characterized and account for 90% of the observed variants in these trios. In all cases, there was no correlation between transmission of the TBX6 variant allele and FIS phenotype. However, an analysis of regulatory markers in osteoblasts showed that rs8060511 is in a putative enhancer element. CONCLUSIONS: Although this study did not identify any TBX6 coding variants that segregate with FIS, we identified a variant that is located in a potential TBX6 enhancer element. Therefore, further investigation of the region is needed.

Evaluation of commercially available RNA amplification kits for RNA sequencing using very low input amounts of total RNA.

This article includes supplemental data. Please visit http://www.fasebj.org to obtain this information.Multiple recent publications on RNA sequencing (RNA-seq) have demonstrated the power of next-generation sequencing technologies in whole-transcriptome analysis. Vendor-specific protocols used for RNA library construction often require at least 100 ng total RNA. However, under certain conditions, much less RNA is available for library construction. In these cases, effective transcriptome profiling requires amplification of subnanogram amounts of RNA. Several commercial RNA amplification kits are available for amplification prior to library construction for next-generation sequencing, but these kits have not been comprehensively field evaluated for accuracy and performance of RNA-seq for picogram amounts of RNA. To address this, 4 types of amplification kits were tested with 3 different concentrations, from 5 ng to 50 pg, of a commercially available RNA. Kits were tested at multiple sites to assess reproducibility and ease of use. The human total reference RNA used was spiked with a control pool of RNA molecules in order to further evaluate quantitative recovery of input material. Additional control data sets were generated from libraries constructed following polyA selection or ribosomal depletion using established kits and protocols. cDNA was collected from the different sites, and libraries were synthesized at a single site using established protocols. Sequencing runs were carried out on the Illumina platform. Numerous metrics were compared among the kits and dilutions used. Overall, no single kit appeared to meet all the challenges of small input material. However, it is encouraging that excellent data can be recovered with even the 50 pg input total RNA.