Distinctiveness and continuity in transcriptome and connectivity in the anterior-posterior axis of the paraventricular nucleus of the thalamus.

The paraventricular nucleus of the thalamus (PVT) projects axons to multiple areas, mediates a wide range of behaviors, and exhibits regional heterogeneity in both functions and axonal projections. Still, questions regarding the cell types present in the PVT and the extent of their differences remain inadequately addressed. We applied single-cell RNA sequencing to depict the transcriptomic characteristics of mouse PVT neurons. We found that one of the most significant variances in the PVT transcriptome corresponded to the anterior-posterior axis. While the single-cell transcriptome classified PVT neurons into five types, our transcriptomic and histological analyses showed continuity among the cell types. We discovered that anterior and posterior subpopulations had nearly non-overlapping projection patterns, while another population showed intermediate patterns. In addition, these subpopulations responded differently to appetite-related neuropeptides, with their activation showing opposing effects on food consumption. Our studies unveiled the contrasts and the continuity of PVT neurons that underpin their function.

WUSCHEL-RELATED HOMEOBOX 13 suppresses de novo shoot regeneration via cell fate control of pluripotent callus.

Plants can regenerate their bodies via de novo establishment of shoot apical meristems (SAMs) from pluripotent callus. Only a small fraction of callus cells is eventually specified into SAMs but the molecular mechanisms underlying fate specification remain obscure. The expression of WUSCHEL (WUS) is an early hallmark of SAM fate acquisition. Here, we show that a WUS paralog, WUSCHEL-RELATED HOMEOBOX 13 (WOX13), negatively regulates SAM formation from callus in Arabidopsis thaliana. WOX13 promotes non-meristematic cell fate via transcriptional repression of WUS and other SAM regulators and activation of cell wall modifiers. Our Quartz-Seq2-based single cell transcriptome revealed that WOX13 plays key roles in determining cellular identity of callus cell population. We propose that reciprocal inhibition between WUS and WOX13 mediates critical cell fate determination in pluripotent cell population, which has a major impact on regeneration efficiency.

Single-cell transcriptomics uncovers EGFR signaling-mediated gastric progenitor cell differentiation in stomach homeostasis.

Defects in gastric progenitor cell differentiation are associated with various gastric disorders, including atrophic gastritis, intestinal metaplasia, and gastric cancer. However, the mechanisms underlying the multilineage differentiation of gastric progenitor cells during healthy homeostasis remain poorly understood. Here, using a single-cell RNA sequencing method, Quartz-Seq2, we analyzed the gene expression dynamics of progenitor cell differentiation toward pit cell, neck cell, and parietal cell lineages in healthy adult mouse corpus tissues. Enrichment analysis of pseudotime-dependent genes and a gastric organoid assay revealed that EGFR-ERK signaling promotes pit cell differentiation, whereas NF-κB signaling maintains gastric progenitor cells in an undifferentiated state. In addition, pharmacological inhibition of EGFR in vivo resulted in a decreased number of pit cells. Although activation of EGFR signaling in gastric progenitor cells has been suggested as one of the major inducers of gastric cancers, our findings unexpectedly identified that EGFR signaling exerts a differentiation-promoting function, not a mitogenic function, in normal gastric homeostasis.

A common epigenetic mechanism across different cellular origins underlies systemic immune dysregulation in an idiopathic autism mouse model.

Immune dysregulation plays a key role in the pathogenesis of autism. Changes occurring at the systemic level, from brain inflammation to disturbed innate/adaptive immune in the periphery, are frequently observed in patients with autism; however, the intrinsic mechanisms behind them remain elusive. We hypothesize a common etiology may lie in progenitors of different types underlying widespread immune dysregulation. By single-cell RNA sequencing (sc-RNA seq), we trace the developmental origins of immune dysregulation in a mouse model of idiopathic autism. It is found that both in aorta-gonad-mesonephros (AGM) and yolk sac (YS) progenitors, the dysregulation of HDAC1-mediated epigenetic machinery alters definitive hematopoiesis during embryogenesis and downregulates the expression of the AP-1 complex for microglia development. Subsequently, these changes result in the dysregulation of the immune system, leading to gut dysbiosis and hyperactive microglia in the brain. We further confirm that dysregulated immune profiles are associated with specific microbiota composition, which may serve as a biomarker to identify autism of immune-dysregulated subtypes. Our findings elucidate a shared mechanism for the origin of immune dysregulation from the brain to the gut in autism and provide new insight to dissecting the heterogeneity of autism, as well as the therapeutic potential of targeting immune-dysregulated autism subtypes.

Local states of chromatin compaction at transcription start sites control transcription levels.

The 'open' and 'compact' regions of chromatin are considered to be regions of active and silent transcription, respectively. However, individual genes produce transcripts at different levels, suggesting that transcription output does not depend on the simple open-compact conversion of chromatin, but on structural variations in chromatin itself, which so far have remained elusive. In this study, weakly crosslinked chromatin was subjected to sedimentation velocity centrifugation, which fractionated the chromatin according to its degree of compaction. Open chromatin remained in upper fractions, while compact chromatin sedimented to lower fractions depending on the level of nucleosome assembly. Although nucleosomes were evenly detected in all fractions, histone H1 was more highly enriched in the lower fractions. H1 was found to self-associate and crosslinked to histone H3, suggesting that H1 bound to H3 interacts with another H1 in an adjacent nucleosome to form compact chromatin. Genome-wide analyses revealed that nearly the entire genome consists of compact chromatin without differences in compaction between repeat and non-repeat sequences; however, active transcription start sites (TSSs) were rarely found in compact chromatin. Considering the inverse correlation between chromatin compaction and RNA polymerase binding at TSSs, it appears that local states of chromatin compaction determine transcription levels.

Developmental excitation-inhibition imbalance underlying psychoses revealed by single-cell analyses of discordant twins-derived cerebral organoids.

Despite extensive genetic and neuroimaging studies, detailed cellular mechanisms underlying schizophrenia and bipolar disorder remain poorly understood. Recent progress in single-cell RNA sequencing (scRNA-seq) technologies enables identification of cell-type-specific pathophysiology. However, its application to psychiatric disorders is challenging because of methodological difficulties in analyzing human brains and the confounds due to a lifetime of illness. Brain organoids derived from induced pluripotent stem cells (iPSCs) of the patients are a powerful avenue to investigate the pathophysiological processes. Here, we generated iPSC-derived cerebral organoids from monozygotic twins discordant for psychosis. scRNA-seq analysis of the organoids revealed enhanced GABAergic specification and reduced cell proliferation following diminished Wnt signaling in the patient, which was confirmed in iPSC-derived forebrain neuronal cells. Two additional monozygotic twin pairs discordant for schizophrenia also confirmed the excess GABAergic specification of the patients' neural progenitor cells. With a well-controlled genetic background, our data suggest that unbalanced specification of excitatory and inhibitory neurons during cortical development underlies psychoses.

Benchmarking single-cell RNA-sequencing protocols for cell atlas projects.

Single-cell RNA sequencing (scRNA-seq) is the leading technique for characterizing the transcriptomes of individual cells in a sample. The latest protocols are scalable to thousands of cells and are being used to compile cell atlases of tissues, organs and organisms. However, the protocols differ substantially with respect to their RNA capture efficiency, bias, scale and costs, and their relative advantages for different applications are unclear. In the present study, we generated benchmark datasets to systematically evaluate protocols in terms of their power to comprehensively describe cell types and states. We performed a multicenter study comparing 13 commonly used scRNA-seq and single-nucleus RNA-seq protocols applied to a heterogeneous reference sample resource. Comparative analysis revealed marked differences in protocol performance. The protocols differed in library complexity and their ability to detect cell-type markers, impacting their predictive value and suitability for integration into reference cell atlases. These results provide guidance both for individual researchers and for consortium projects such as the Human Cell Atlas.

Primed to Naive-Like Conversion of the Common Marmoset Embryonic Stem Cells.

Mammalian pluripotent stem cells are thought to exist in two states: naive and primed. Generally, unlike those in rodents, pluripotent stem cells in primates, including humans, are regarded as being in the primed pluripotent state. Recently, several groups reported the existence of naive pluripotent stem cells in humans. In this study, we report the conversion of primed state embryonic stem cells from common marmoset, a New World monkey, to the naive state using transgenes. The cells showed typical naive state features, including dome-like colony morphology, growth factor requirement, gene expression profile, X chromosome activation state, and energy metabolic status. Moreover, interspecies chimeric embryo formation ability with mouse embryos was increased in the naive state. This technique can be applied in basic medical research using nonhuman primates, such as preclinical use of naive pluripotent stem cells and generating genetically modified primates.

Cell type-specific transcriptome analysis unveils secreted signaling molecule genes expressed in apical epithelial cap during appendage regeneration.

Wound epidermis (WE) and the apical epithelial cap (AEC) are believed to trigger regeneration of amputated appendages such as limb and tail in amphibians by producing certain secreted signaling molecules. To date, however, only limited information about the molecular signatures of these epidermal structures is available. Here we used a transgenic Xenopus laevis line harboring the enhanced green fluorescent protein (egfp) gene under control of an es1 gene regulatory sequence to isolate WE/AEC cells by performing fluorescence-activated cell sorting during the time course of tail regeneration (day 1, day 2, day 3 and day 4 after amputation). Time-course transcriptome analysis of these isolated WE/AEC cells revealed that more than 8,000 genes, including genes involved in signaling pathways such as those of reactive oxygen species, fibroblast growth factor (FGF), canonical and non-canonical Wnt, transforming growth factor ß (TGF ß) and Notch, displayed dynamic changes of their expression during tail regeneration. Notably, this approach enabled us to newly identify seven secreted signaling molecule genes (mdk, fstl, slit1, tgfß1, bmp7.1, angptl2 and egfl6) that are highly expressed in tail AEC cells. Among these genes, five (mdk, fstl, slit1, tgfß1 and bmp7.1) were also highly expressed in limb AEC cells but the other two (angptl2 and egfl6) are specifically expressed in tail AEC cells. Interestingly, there was no expression of fgf8 in tail WE/AEC cells, whose expression and pivotal role in limb AEC cells have been reported previously. Thus, we identified common and different properties between tail and limb AEC cells.

Strategies for Converting RNA to Amplifiable cDNA for Single-Cell RNA Sequencing Methods.

This review describes the features of molecular biology techniques for single-cell RNA sequencing (scRNA-seq), including methods developed in our laboratory. Existing scRNA-seq methods require the conversion of first-strand cDNA to amplifiable cDNA followed by whole-transcript amplification. There are three primary strategies for this conversion: poly-A tagging, template switching, and RNase H-DNA polymerase I-mediated second-strand cDNA synthesis for in vitro transcription. We discuss the merits and limitations of these strategies and describe our Reverse Transcription with Random Displacement Amplification technology that allows for direct first-strand cDNA amplification from RNA without the need for conversion to an amplifiable cDNA. We believe that this review provides all users of single-cell transcriptome technologies with an understanding of the relationship between the quantitative performance of various methods and their molecular features.

Polysaccharide hydrolase of the hadal zone amphipods Hirondellea gigas.

Hirondellea species are common inhabitants in the hadal region deeper than 7,000 m. We found that Hirondellea gigas thrived in the Challenger Deep possessed polysaccharide hydrolases as digestive enzymes. To obtain various enzymes of other H. gigas, we captured amphipods from the Japan Trench, and Izu-Ogasawara (Bonin) Trench. A phylogenetic analysis based on the cytochrome oxidase I gene showed close relationships among amphipods, despite the geographic distance between the localities. However, several differences in enzymatic properties were observed in these H. gigas specimens. We also carried out RNA sequencing of H. gigas from the Izu-Ogasawara Trench. The cellulase gene of H. gigas was highly homologous to cellobiohydrolase of Glucosyl Hydrolase family 7 (GH7). On the other hand, enzymatic properties of H. gigas's cellulase were different from those of typical GH7 cellobiohydrolase. Thus, these results indicate that hadal-zone amphipod can be good candidates as the new enzyme resource.

Quartz-Seq2: a high-throughput single-cell RNA-sequencing method that effectively uses limited sequence reads.

High-throughput single-cell RNA-seq methods assign limited unique molecular identifier (UMI) counts as gene expression values to single cells from shallow sequence reads and detect limited gene counts. We thus developed a high-throughput single-cell RNA-seq method, Quartz-Seq2, to overcome these issues. Our improvements in the reaction steps make it possible to effectively convert initial reads to UMI counts, at a rate of 30-50%, and detect more genes. To demonstrate the power of Quartz-Seq2, we analyzed approximately 10,000 transcriptomes from in vitro embryonic stem cells and an in vivo stromal vascular fraction with a limited number of reads.

Single-cell full-length total RNA sequencing uncovers dynamics of recursive splicing and enhancer RNAs.

Total RNA sequencing has been used to reveal poly(A) and non-poly(A) RNA expression, RNA processing and enhancer activity. To date, no method for full-length total RNA sequencing of single cells has been developed despite the potential of this technology for single-cell biology. Here we describe random displacement amplification sequencing (RamDA-seq), the first full-length total RNA-sequencing method for single cells. Compared with other methods, RamDA-seq shows high sensitivity to non-poly(A) RNA and near-complete full-length transcript coverage. Using RamDA-seq with differentiation time course samples of mouse embryonic stem cells, we reveal hundreds of dynamically regulated non-poly(A) transcripts, including histone transcripts and long noncoding RNA Neat1. Moreover, RamDA-seq profiles recursive splicing in >300-kb introns. RamDA-seq also detects enhancer RNAs and their cell type-specific activity in single cells. Taken together, we demonstrate that RamDA-seq could help investigate the dynamics of gene expression, RNA-processing events and transcriptional regulation in single cells.

Chromatin remodeler CHD7 regulates the stem cell identity of human neural progenitors.

Multiple congenital disorders often present complex phenotypes, but how the mutation of individual genetic factors can lead to multiple defects remains poorly understood. In the present study, we used human neuroepithelial (NE) cells and CHARGE patient-derived cells as an in vitro model system to identify the function of chromodomain helicase DNA-binding 7 (CHD7) in NE-neural crest bifurcation, thus revealing an etiological link between the central nervous system (CNS) and craniofacial anomalies observed in CHARGE syndrome. We found that CHD7 is required for epigenetic activation of superenhancers and CNS-specific enhancers, which support the maintenance of the NE and CNS lineage identities. Furthermore, we found that BRN2 and SOX21 are downstream effectors of CHD7, which shapes cellular identities by enhancing a CNS-specific cellular program and indirectly repressing non-CNS-specific cellular programs. Based on our results, CHD7, through its interactions with superenhancer elements, acts as a regulatory hub in the orchestration of the spatiotemporal dynamics of transcription factors to regulate NE and CNS lineage identities.

Constrained vertebrate evolution by pleiotropic genes.

Despite morphological diversification of chordates over 550 million years of evolution, their shared basic anatomical pattern (or 'bodyplan') remains conserved by unknown mechanisms. The developmental hourglass model attributes this to phylum-wide conserved, constrained organogenesis stages that pattern the bodyplan (the phylotype hypothesis); however, there has been no quantitative testing of this idea with a phylum-wide comparison of species. Here, based on data from early-to-late embryonic transcriptomes collected from eight chordates, we suggest that the phylotype hypothesis would be better applied to vertebrates than chordates. Furthermore, we found that vertebrates' conserved mid-embryonic developmental programmes are intensively recruited to other developmental processes, and the degree of the recruitment positively correlates with their evolutionary conservation and essentiality for normal development. Thus, we propose that the intensively recruited genetic system during vertebrates' organogenesis period imposed constraints on its diversification through pleiotropic constraints, which ultimately led to the common anatomical pattern observed in vertebrates.

Quartz-Seq: a highly reproducible and sensitive single-cell RNA sequencing method, reveals non-genetic gene-expression heterogeneity.

Development of a highly reproducible and sensitive single-cell RNA sequencing (RNA-seq) method would facilitate the understanding of the biological roles and underlying mechanisms of non-genetic cellular heterogeneity. In this study, we report a novel single-cell RNA-seq method called Quartz-Seq that has a simpler protocol and higher reproducibility and sensitivity than existing methods. We show that single-cell Quartz-Seq can quantitatively detect various kinds of non-genetic cellular heterogeneity, and can detect different cell types and different cell-cycle phases of a single cell type. Moreover, this method can comprehensively reveal gene-expression heterogeneity between single cells of the same cell type in the same cell-cycle phase.

CDC-48/p97 is required for proper meiotic chromosome segregation via controlling AIR-2/Aurora B kinase localization in Caenorhabditis elegans.

CDC-48/p97 is a AAA (ATPases associated with diverse cellular activities) chaperone involved in protein conformational changes such as the disassembly of protein complexes. We previously reported that Caenorhabditis elegans CDC-48.1 and CDC-48.2 (CDC-48s) are essential for the progression of meiosis I metaphase. Here, we report that CDC-48s are required for proper chromosome segregation during meiosis in C. elegans. In wild-type worms, at the diakinesis phase, phosphorylation of histone H3, one of the known substrates of aurora B kinase (AIR-2), on meiosis I chromatids correlated with AIR-2 localization at the cohesion sites of homologous chromatids. Conversely, depletion of CDC-48s resulted in a significant expansion of signals for AIR-2 and phosphorylated histone H3 over the entire length of meiotic chromosomes, leading to defective chromosome segregation, while the total amount of AIR-2 in lysates was not changed by the depletion of CDC-48s. The defective segregation of meiotic chromosomes caused by the depletion of CDC-48s was suppressed by the simultaneous depletion of AIR-2 and is similar to that observed following the depletion of protein phosphatase 1 (PP1) phosphatases. However, the amount and localization of PP1 were not changed by the depletion of CDC-48s. These results suggest that CDC-48s control the restricted localization of AIR-2 to the cohesion sites of homologous chromatids in meiosis I.

Recent advances in p97/VCP/Cdc48 cellular functions.

p97/VCP/Cdc48 is one of the best-characterized type II AAA (ATPases associated with diverse cellular activities) ATPases. p97 is suggested to be a ubiquitin-selective chaperone and its key function is to disassemble protein complexes. p97 is involved in a wide variety of cellular activities. Recently, novel functions, namely autophagy and mitochondrial quality control, for p97 have been uncovered. p97 was identified as a causative factor for inclusion body myopathy associated with Paget disease of bone and frontotemporal dementia (IBMPFD) and more recently as a causative factor for amyotrophic lateral sclerosis (ALS). In this review, we will summarize and discuss recent progress and topics in p97 functions and the relationship to its associated diseases.