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.

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.

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.

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.

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.

Serine and SAM Responsive Complex SESAME Regulates Histone Modification Crosstalk by Sensing Cellular Metabolism.

Pyruvate kinase M2 (PKM2) is a key enzyme for glycolysis and catalyzes the conversion of phosphoenolpyruvate (PEP) to pyruvate, which supplies cellular energy. PKM2 also phosphorylates histone H3 threonine 11 (H3T11); however, it is largely unknown how PKM2 links cellular metabolism to chromatin regulation. Here, we show that the yeast PKM2 homolog, Pyk1, is a part of a novel protein complex named SESAME (Serine-responsive SAM-containing Metabolic Enzyme complex), which contains serine metabolic enzymes, SAM (S-adenosylmethionine) synthetases, and an acetyl-CoA synthetase. SESAME interacts with the Set1 H3K4 methyltransferase complex, which requires SAM synthesized from SESAME, and recruits SESAME to target genes, resulting in phosphorylation of H3T11. SESAME regulates the crosstalk between H3K4 methylation and H3T11 phosphorylation by sensing glycolysis and glucose-derived serine metabolism. This leads to auto-regulation of PYK1 expression. Thus, our study provides insights into the mechanism of regulating gene expression, responding to cellular metabolism via chromatin modifications.

Swi/Snf dynamics on stress-responsive genes is governed by competitive bromodomain interactions.

The Swi/Snf chromatin remodeling complex functions to alter nucleosome positions by either sliding nucleosomes on DNA or the eviction of histones. The presence of histone acetylation and activator-dependent recruitment and retention of Swi/Snf is important for its efficient function. It is not understood, however, why such mechanisms are required to enhance Swi/Snf activity on nucleosomes. Snf2, the catalytic subunit of the Swi/Snf remodeling complex, has been shown to be a target of the Gcn5 acetyltransferase. Our study found that acetylation of Snf2 regulates both recruitment and release of Swi/Snf from stress-responsive genes. Also, the intramolecular interaction of the Snf2 bromodomain with the acetylated lysine residues on Snf2 negatively regulates binding and remodeling of acetylated nucleosomes by Swi/Snf. Interestingly, the presence of transcription activators mitigates the effects of the reduced affinity of acetylated Snf2 for acetylated nucleosomes. Supporting our in vitro results, we found that activator-bound genes regulating metabolic processes showed greater retention of the Swi/Snf complex even when Snf2 was acetylated. Our studies demonstrate that competing effects of (1) Swi/Snf retention by activators or high levels of histone acetylation and (2) Snf2 acetylation-mediated release regulate dynamics of Swi/Snf occupancy at target genes.

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.

Dynamic transcriptional events in embryonic stem cells mediated by the super elongation complex (SEC).

Transcriptional regulation of developmentally controlled genes is at the heart of differentiation and organogenesis. In this study, we performed global genomic analyses in murine embryonic stem (ES) cells and in human cells in response to activation signals. We identified an essential role for the ELL (eleven-nineteen lysine-rich leukemia gene)/P-TEFb (positive transcription elongation factor)-containing super elongation complex (SEC) in the regulation of gene expression, including several genes bearing paused RNA polymerase II (Pol II). Paused Pol II has been proposed to be associated with loci that respond rapidly to environmental stimuli. However, our studies in ES cells also identified a requirement for SEC at genes without paused Pol II, which also respond dynamically to differentiation signals. Our findings suggest that SEC is a major class of active P-TEFb-containing complexes required for transcriptional activation in response to environmental cues such as differentiation signals.

Analysis of dynamic changes in retinoid-induced transcription and epigenetic profiles of murine Hox clusters in ES cells.

The clustered Hox genes, which are highly conserved across metazoans, encode homeodomain-containing transcription factors that provide a blueprint for segmental identity along the body axis. Recent studies have underscored that in addition to encoding Hox genes, the homeotic clusters contain key noncoding RNA genes that play a central role in development. In this study, we have taken advantage of genome-wide approaches to provide a detailed analysis of retinoic acid (RA)-induced transcriptional and epigenetic changes within the homeotic clusters of mouse embryonic stem cells. Although there is a general colinear response, our analyses suggest a lack of strict colinearity for several genes in the HoxA and HoxB clusters. We have identified transcribed novel noncoding RNAs (ncRNAs) and their cis-regulatory elements that function in response to RA and demonstrated that the expression of these ncRNAs from both strands represent some of the most rapidly induced transcripts in ES cells. Finally, we have provided dynamic analyses of chromatin modifications for the coding and noncoding genes expressed upon activation and suggest that active transcription can occur in the presence of chromatin modifications and machineries associated with repressed transcription state over the clusters. Overall, our data provide a resource for a better understanding of the dynamic nature of the coding and noncoding transcripts and their associated chromatin marks in the regulation of homeotic gene transcription during development.

Set2 methylation of histone H3 lysine 36 suppresses histone exchange on transcribed genes.

Set2-mediated methylation of histone H3 at Lys 36 (H3K36me) is a co-transcriptional event that is necessary for the activation of the Rpd3S histone deacetylase complex, thereby maintaining the coding region of genes in a hypoacetylated state. In the absence of Set2, H3K36 or Rpd3S acetylated histones accumulate on open reading frames (ORFs), leading to transcription initiation from cryptic promoters within ORFs. Although the co-transcriptional deacetylation pathway is well characterized, the factors responsible for acetylation are as yet unknown. Here we show that, in yeast, co-transcriptional acetylation is achieved in part by histone exchange over ORFs. In addition to its function of targeting and activating the Rpd3S complex, H3K36 methylation suppresses the interaction of H3 with histone chaperones, histone exchange over coding regions and the incorporation of new acetylated histones. Thus, Set2 functions both to suppress the incorporation of acetylated histones and to signal for the deacetylation of these histones in transcribed genes. By suppressing spurious cryptic transcripts from initiating within ORFs, this pathway is essential to maintain the accuracy of transcription by RNA polymerase II.

Improved transcription and translation with L-leucine stimulation of mTORC1 in Roberts syndrome.

BACKGROUND: Roberts syndrome (RBS) is a human developmental disorder caused by mutations in the cohesin acetyltransferase ESCO2. We previously reported that mTORC1 signaling was depressed and overall translation was reduced in RBS cells and zebrafish models for RBS. Treatment of RBS cells and zebrafish RBS models with L-leucine partially rescued mTOR function and protein synthesis, correlating with increased cell division and improved development. RESULTS: In this study, we use RBS cells to model mTORC1 repression and analyze transcription and translation with ribosome profiling to determine gene-level effects of L-leucine. L-leucine treatment partially rescued translational efficiency of ribosomal subunits, translation initiation factors, snoRNA production, and mitochondrial function in RBS cells, consistent with these processes being mTORC1 controlled. In contrast, other genes are differentially expressed independent of L-leucine treatment, including imprinted genes such as H19 and GTL2, miRNAs regulated by GTL2, HOX genes, and genes in nucleolar associated domains. CONCLUSIONS: Our study distinguishes between gene expression changes in RBS cells that are TOR dependent and those that are independent. Some of the TOR independent gene expression changes likely reflect the architectural role of cohesin in chromatin looping and gene expression. This study reveals the dramatic rescue effects of L-leucine stimulation of mTORC1 in RBS cells and supports that normal gene expression and translation requires ESCO2 function.

Promoter-driven splicing regulation in fission yeast.

The meiotic cell cycle is modified from the mitotic cell cycle by having a pre-meiotic S phase that leads to high levels of recombination, two rounds of nuclear division with no intervening DNA synthesis and a reductional pattern of chromosome segregation. Rem1 is a cyclin that is only expressed during meiosis in the fission yeast Schizosaccharomyces pombe. Cells in which rem1 has been deleted show decreased intragenic meiotic recombination and a delay at the onset of meiosis I (ref. 1). When ectopically expressed in mitotically growing cells, Rem1 induces a G1 arrest followed by severe mitotic catastrophes. Here we show that rem1 expression is regulated at the level of both transcription and splicing, encoding two proteins with different functions depending on the intron retention. We have determined that the regulation of rem1 splicing is not dependent on any transcribed region of the gene. Furthermore, when the rem1 promoter is fused to other intron-containing genes, the chimaeras show a meiotic-specific regulation of splicing, exactly the same as endogenous rem1. This regulation is dependent on two transcription factors of the forkhead family, Mei4 (ref. 2) and Fkh2 (ref. 3). Whereas Mei4 induces both transcription and splicing of rem1, Fkh2 is responsible for the intron retention of the transcript during vegetative growth and the pre-meiotic S phase.

C G composition in transposon-derived genes is increased in FXD with perturbed immune system.

Increasing incidence of Fragile X disorders (FXD) and of immune-mediated disorders in FXD suggests that additional factors besides FMR1 mutations contribute to the pathogenesis. Here, we discovered that the expression levels or splicing of specific transposon element (TE)-derived genes, regulating purine metabolism and immune responses against viral infections are altered in FXD. These genes include HLA genes clustered in chr6p21.3 and viral responsive genes in chr5q15. Remarkably, these TE-derived genes contain a low A T/C G suggesting base substitutions of A T to C G. The TE-derived genes with changed expression levels contained a higher content of 5'-CG-3' dinucleotides in FXD compared to healthy donors. This resembles the genomes of some RNA viruses, which maintain high contents of CG dinucleotides to sustain their latent infection exploiting antiviral responses. Thus, past viral infections may have persisted as TEs, provoking immune-mediated disorders in FXD.

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.

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.

Epigenetic control and inheritance of rDNA arrays.

Ribosomal RNA (rRNA) genes exist in multiple copies arranged in tandem arrays known as ribosomal DNA (rDNA). The total number of gene copies is variable, and the mechanisms buffering this copy number variation remain unresolved. We surveyed the number, distribution, and activity of rDNA arrays at the level of individual chromosomes across multiple human and primate genomes. Each individual possessed a unique fingerprint of copy number distribution and activity of rDNA arrays. In some cases, entire rDNA arrays were transcriptionally silent. Silent rDNA arrays showed reduced association with the nucleolus and decreased interchromosomal interactions, indicating that the nucleolar organizer function of rDNA depends on transcriptional activity. Methyl-sequencing of flow-sorted chromosomes, combined with long read sequencing, showed epigenetic modification of rDNA promoter and coding region by DNA methylation. Silent arrays were in a closed chromatin state, as indicated by the accessibility profiles derived from Fiber-seq. Removing DNA methylation restored the transcriptional activity of silent arrays. Array activity status remained stable through the iPS cell re-programming. Family trio analysis demonstrated that the inactive rDNA haplotype can be traced to one of the parental genomes, suggesting that the epigenetic state of rDNA arrays may be heritable. We propose that the dosage of rRNA genes is epigenetically regulated by DNA methylation, and these methylation patterns specify nucleolar organizer function and can propagate transgenerationally.

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.

Expression profiling of S. pombe acetyltransferase mutants identifies redundant pathways of gene regulation.

BACKGROUND: Histone acetyltransferase enzymes (HATs) are implicated in regulation of transcription. HATs from different families may overlap in target and substrate specificity. RESULTS: We isolated the elp3+ gene encoding the histone acetyltransferase subunit of the Elongator complex in fission yeast and characterized the phenotype of an Deltaelp3 mutant. We examined genetic interactions between Deltaelp3 and two other HAT mutants, Deltamst2 and Deltagcn5 and used whole genome microarray analysis to analyze their effects on gene expression. CONCLUSIONS: Comparison of phenotypes and expression profiles in single, double and triple mutants indicate that these HAT enzymes have overlapping functions. Consistent with this, overlapping specificity in histone H3 acetylation is observed. However, there is no evidence for overlap with another HAT enzyme, encoded by the essential mst1+ gene.