Identification and characterization of intermediate states in mammalian neural crest cell epithelial to mesenchymal transition and delamination.

Epithelial to mesenchymal transition (EMT) is a cellular process that converts epithelial cells to mesenchymal cells with migratory potential in developmental and pathological processes. Although originally considered a binary event, EMT in cancer progression involves intermediate states between a fully epithelial and a fully mesenchymal phenotype, which are characterized by distinct combinations of epithelial and mesenchymal markers. This phenomenon has been termed epithelial to mesenchymal plasticity (EMP), however, the intermediate states remain poorly described and it's unclear whether they exist during developmental EMT. Neural crest cells (NCC) are an embryonic progenitor cell population that gives rise to numerous cell types and tissues in vertebrates, and their formation and delamination is a classic example of developmental EMT. However, whether intermediate states also exist during NCC EMT and delamination remains unknown. Through single-cell RNA sequencing of mouse embryos, we identified intermediate NCC states based on their transcriptional signature and then spatially defined their locations in situ in the dorsolateral neuroepithelium. Our results illustrate the importance of cell cycle regulation and functional role for the intermediate stage marker Dlc1 in facilitating mammalian cranial NCC delamination and may provide new insights into mechanisms regulating pathological EMP.

Identification and characterization of intermediate states in mammalian neural crest cell epithelial to mesenchymal transition and delamination.

Epithelial to mesenchymal transition (EMT) is a cellular process that converts epithelial cells to mesenchymal cells with migratory potential in both developmental and pathological processes. Although originally considered a binary event, EMT in cancer progression involves intermediate states between a fully epithelial and a fully mesenchymal phenotype, which are characterized by distinct combinations of epithelial and mesenchymal markers. This phenomenon has been termed epithelial to mesenchymal plasticity (EMP), however, the intermediate states remain poorly described and it's unclear whether they exist during developmental EMT. Neural crest cells (NCC) are an embryonic progenitor cell population that gives rise to numerous cell types and tissues in vertebrates, and their formation is a classic example of developmental EMT. An important feature of NCC development is their delamination from the neuroepithelium via EMT, following which NCC migrate throughout the embryo and undergo differentiation. NCC delamination shares similar changes in cellular state and structure with cancer cell invasion. However, whether intermediate states also exist during NCC EMT and delamination remains unknown. Through single cell RNA sequencing, we identified intermediate NCC states based on their transcriptional signature and then spatially defined their locations in situ in the dorsolateral neuroepithelium. Our results illustrate the progressive transcriptional and spatial transitions from premigratory to migratory cranial NCC during EMT and delamination. Of note gene expression and trajectory analysis indicate that distinct intermediate populations of NCC delaminate in either S phase or G2/M phase of the cell cycle, and the importance of cell cycle regulation in facilitating mammalian cranial NCC delamination was confirmed through cell cycle inhibition studies. Additionally, transcriptional knockdown revealed a functional role for the intermediate stage marker Dlc1 in regulating NCC delamination and migration. Overall, our work identifying and characterizing the intermediate cellular states, processes, and molecular signals that regulate mammalian NCC EMT and delamination furthers our understanding of developmental EMP and may provide new insights into mechanisms regulating pathological EMP.

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.

Paraspeckles interact with SWI/SNF subunit ARID1B to regulate transcription and splicing.

Paraspeckles are subnuclear RNA-protein structures that are implicated in important processes including cellular stress response, differentiation, and cancer progression. However, it is unclear how paraspeckles impart their physiological effect at the molecular level. Through biochemical analyses, we show that paraspeckles interact with the SWI/SNF chromatin-remodeling complex. This is specifically mediated by the direct interaction of the long-non-coding RNA NEAT1 of the paraspeckles with ARID1B of the cBAF-type SWI/SNF complex. Strikingly, ARID1B depletion, in addition to resulting in loss of interaction with the SWI/SNF complex, decreases the binding of paraspeckle proteins to chromatin modifiers, transcription factors, and histones. Functionally, the loss of ARID1B and NEAT1 influences the transcription and the alternative splicing of a common set of genes. Our findings reveal that dynamic granules such as the paraspeckles may leverage the specificity of epigenetic modifiers to impart their regulatory effect, thus providing a molecular basis for their function.

H3K36 methylation and DNA-binding both promote Ioc4 recruitment and Isw1b remodeler function.

The Isw1b chromatin-remodeling complex is specifically recruited to gene bodies to help retain pre-existing histones during transcription by RNA polymerase II. Recruitment is dependent on H3K36 methylation and the Isw1b subunit Ioc4, which contains an N-terminal PWWP domain. Here, we present the crystal structure of the Ioc4-PWWP domain, including a detailed functional characterization of the domain on its own as well as in the context of full-length Ioc4 and the Isw1b remodeler. The Ioc4-PWWP domain preferentially binds H3K36me3-containing nucleosomes. Its ability to bind DNA is required for nucleosome binding. It is also furthered by the unique insertion motif present in Ioc4-PWWP. The ability to bind H3K36me3 and DNA promotes the interaction of full-length Ioc4 with nucleosomes in vitro and they are necessary for its recruitment to gene bodies in vivo. Furthermore, a fully functional Ioc4-PWWP domain promotes efficient remodeling by Isw1b and the maintenance of ordered chromatin in vivo, thereby preventing the production of non-coding RNAs.

Metabolic regulation of telomere silencing by SESAME complex-catalyzed H3T11 phosphorylation.

Telomeres are organized into a heterochromatin structure and maintenance of silent heterochromatin is required for chromosome stability. How telomere heterochromatin is dynamically regulated in response to stimuli remains unknown. Pyruvate kinase Pyk1 forms a complex named SESAME (Serine-responsive SAM-containing Metabolic Enzyme complex) to regulate gene expression by phosphorylating histone H3T11 (H3pT11). Here, we identify a function of SESAME in regulating telomere heterochromatin structure. SESAME phosphorylates H3T11 at telomeres, which maintains SIR (silent information regulator) complex occupancy at telomeres and protects Sir2 from degradation by autophagy. Moreover, SESAME-catalyzed H3pT11 directly represses autophagy-related gene expression to further prevent autophagy-mediated Sir2 degradation. By promoting H3pT11, serine increases Sir2 protein levels and enhances telomere silencing. Loss of H3pT11 leads to reduced Sir2 and compromised telomere silencing during chronological aging. Together, our study provides insights into dynamic regulation of silent heterochromatin by histone modifications and autophagy in response to cell metabolism and aging.

Single-cell reconstruction with spatial context of migrating neural crest cells and their microenvironments during vertebrate head and neck formation.

The dynamics of multipotent neural crest cell differentiation and invasion as cells travel throughout the vertebrate embryo remain unclear. Here, we preserve spatial information to derive the transcriptional states of migrating neural crest cells and the cellular landscape of the first four chick cranial to cardiac branchial arches (BA1-4) using label-free, unsorted single-cell RNA sequencing. The faithful capture of branchial arch-specific genes led to identification of novel markers of migrating neural crest cells and 266 invasion genes common to all BA1-4 streams. Perturbation analysis of a small subset of invasion genes and time-lapse imaging identified their functional role to regulate neural crest cell behaviors. Comparison of the neural crest invasion signature to other cell invasion phenomena revealed a shared set of 45 genes, a subset of which showed direct relevance to human neuroblastoma cell lines analyzed after exposure to the in vivo chick embryonic neural crest microenvironment. Our data define an important spatio-temporal reference resource to address patterning of the vertebrate head and neck, and previously unidentified cell invasion genes with the potential for broad impact.

Translation of small downstream ORFs enhances translation of canonical main open reading frames.

In addition to canonical open reading frames (ORFs), thousands of translated small ORFs (containing less than 100 codons) have been identified in untranslated mRNA regions (UTRs) across eukaryotes. Small ORFs in 5' UTRs (upstream (u)ORFs) often repress translation of the canonical ORF within the same mRNA. However, the function of translated small ORFs in the 3' UTRs (downstream (d)ORFs) is unknown. Contrary to uORFs, we find that translation of dORFs enhances translation of their corresponding canonical ORFs. This translation stimulatory effect of dORFs depends on the number of dORFs, but not the length or peptide they encode. We propose that dORFs represent a new, strong, and universal translation regulatory mechanism in vertebrates.

An Atlas of Transcription Factors Expressed in Male Pupal Terminalia of Drosophila melanogaster.

During development, transcription factors and signaling molecules govern gene regulatory networks to direct the formation of unique morphologies. As changes in gene regulatory networks are often implicated in morphological evolution, mapping transcription factor landscapes is important, especially in tissues that undergo rapid evolutionary change. The terminalia (genital and anal structures) of Drosophila melanogaster and its close relatives exhibit dramatic changes in morphology between species. While previous studies have identified network components important for patterning the larval genital disc, the networks governing adult structures during pupal development have remained uncharted. Here, we performed RNA-seq in whole Drosophila melanogaster male terminalia followed by in situ hybridization for 100 highly expressed transcription factors during pupal development. We find that the male terminalia are highly patterned during pupal stages and that specific transcription factors mark separate structures and substructures. Our results are housed online in a searchable database (https://flyterminalia.pitt.edu/) as a resource for the community. This work lays a foundation for future investigations into the gene regulatory networks governing the development and evolution of Drosophila terminalia.

Set1-catalyzed H3K4 trimethylation antagonizes the HIR/Asf1/Rtt106 repressor complex to promote histone gene expression and chronological life span.

Aging is the main risk factor for many prevalent diseases. However, the molecular mechanisms regulating aging at the cellular level are largely unknown. Using single cell yeast as a model organism, we found that reducing yeast histone proteins accelerates chronological aging and increasing histone supply extends chronological life span. We sought to identify pathways that regulate chronological life span by controlling intracellular histone levels. Thus, we screened the histone H3/H4 mutant library to uncover histone residues and posttranslational modifications that regulate histone gene expression. We discovered 15 substitution mutations with reduced histone proteins and 5 mutations with increased histone proteins. Among these mutations, we found Set1 complex-catalyzed H3K4me3 promotes histone gene transcription and maintains normal chronological life span. Unlike the canonical functions of H3K4me3 in gene expression, H3K4me3 facilitates histone gene transcription by acting as a boundary to restrict the spread of the repressive HIR/Asf1/Rtt106 complex from histone gene promoters. Collectively, our study identified a novel mechanism by which H3K4me3 antagonizes the HIR/Asf1/Rtt106 repressor complex to promote histone gene expression and extend chronological life span.

ClusterMap: compare multiple single cell RNA-Seq datasets across different experimental conditions.

MOTIVATION: Single cell RNA-Seq (scRNA-Seq) facilitates the characterization of cell type heterogeneity and developmental processes. Further study of single cell profiles across different conditions enables the understanding of biological processes and underlying mechanisms at the sub-population level. However, developing proper methodology to compare multiple scRNA-Seq datasets remains challenging. RESULTS: We have developed ClusterMap, a systematic method and workflow to facilitate the comparison of scRNA-seq profiles across distinct biological contexts. Using hierarchical clustering of the marker genes of each sub-group, ClusterMap matches the sub-types of cells across different samples and provides 'similarity' as a metric to quantify the quality of the match. We introduce a purity tree cut method designed specifically for this matching problem. We use Circos plot and regrouping method to visualize the results concisely. Furthermore, we propose a new metric 'separability' to summarize sub-population changes among all sample pairs. In the case studies, we demonstrate that ClusterMap has the ability to provide us further insight into the different molecular mechanisms of cellular sub-populations across different conditions. AVAILABILITY AND IMPLEMENTATION: ClusterMap is implemented in R and available at https://github.com/xgaoo/ClusterMap. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MPTAC Determines APP Fragmentation via Sensing Sulfur Amino Acid Catabolism.

Metabolic disorder has been suggested to underlie Alzheimer's disease (AD). However, the decisive molecular linkages remain unclear. We discovered that human Molybdopterin Synthase Associating Complex, MPTAC, promotes sulfur amino acid catabolism to prevent oxidative damage from excess sulfur amino acids, which, in turn, advances fatty acid oxidation and acetyl coenzyme A (acetyl-CoA) synthesis. The association of MPTAC with Protein arginine (R) Methyltransferase 5 (PRMT5) complex and small nuclear ribonucleoprotein (SNRP) splicing factors enables SNRPs to sense metabolic states through their methylation. This promotes the splicing fidelity of amyloid precursor protein (APP) pre-mRNA and proper APP fragmentation, abnormalities of which have been observed in the platelets of AD patients. The functions of MPTAC are crucial to maintain expression of drebrin 1, which is required for synaptic plasticity, through prevention from oxidative damage. Thus, adjustment of sulfur amino acid catabolism by MPTAC prevents events that occur early in the onset of AD.

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.

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.

Histone H3 threonine 11 phosphorylation by Sch9 and CK2 regulates chronological lifespan by controlling the nutritional stress response.

Upon nutritional stress, the metabolic status of cells is changed by nutrient signaling pathways to ensure survival. Altered metabolism by nutrient signaling pathways has been suggested to influence cellular lifespan. However, it remains unclear how chromatin regulation is involved in this process. Here, we found that histone H3 threonine 11 phosphorylation (H3pT11) functions as a marker for nutritional stress and aging. Sch9 and CK2 kinases cooperatively regulate H3pT11 under stress conditions. Importantly, H3pT11 defective mutants prolonged chronological lifespan (CLS) by altering nutritional stress responses. Thus, the phosphorylation of H3T11 by Sch9 and CK2 links a nutritional stress response to chromatin in the regulation of CLS.

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.

Aubergine Controls Germline Stem Cell Self-Renewal and Progeny Differentiation via Distinct Mechanisms.

Piwi family protein Aubergine (Aub) maintains genome integrity in late germ cells of the Drosophila ovary through Piwi-associated RNA-mediated repression of transposon activities. Although it is highly expressed in germline stem cells (GSCs) and early progeny, it remains unclear whether it plays any roles in early GSC lineage development. Here we report that Aub promotes GSC self-renewal and GSC progeny differentiation. RNA-iCLIP results show that Aub binds the mRNAs encoding self-renewal and differentiation factors in cultured GSCs. Aub controls GSC self-renewal by preventing DNA-damage-induced Chk2 activation and by translationally controlling the expression of self-renewal factors. It promotes GSC progeny differentiation by translationally controlling the expression of differentiation factors, including Bam. Therefore, this study reveals a function of Aub in GSCs and their progeny, which promotes translation of self-renewal and differentiation factors by directly binding to its target mRNAs and interacting with translational initiation factors.

Composition and Function of Mutant Swi/Snf Complexes.

The 12-subunit Swi/Snf chromatin remodeling complex is conserved from yeast to humans. It functions to alter nucleosome positions by either sliding nucleosomes on DNA or evicting histones. Interestingly, 20% of all human cancers carry mutations in subunits of the Swi/Snf complex. Many of these mutations cause protein instability and loss, resulting in partial Swi/Snf complexes. Although several studies have shown that histone acetylation and activator-dependent recruitment of Swi/Snf regulate its function, it is less well understood how subunits regulate stability and function of the complex. Using functional proteomic and genomic approaches, we have assembled the network architecture of yeast Swi/Snf. In addition, we find that subunits of the Swi/Snf complex regulate occupancy of the catalytic subunit Snf2, thereby modulating gene transcription. Our findings have direct bearing on how cancer-causing mutations in orthologous subunits of human Swi/Snf may lead to aberrant regulation of gene expression by this complex.