Search

Based on current studies, the incidence of Ewing sarcoma (ES) varies significantly by race and ethnicity, with the disease being most common in patients of European ancestry. However, race/ethnicity has generally been self-reported rather than formally evaluated at a population level using DNA evidence. Additionally, mitochondrial dysfunction is a hallmark of ES, yet there have been no reported studies of mitochondrial genetics in ES. Thus, we evaluated both the mitochondrial and nuclear ancestries of 420 pediatric ES patients in the United States using whole-genome sequencing. We found that the mitochondrial DNA (mtDNA) genomes of only six (1.4 %) patients belonged to African L haplogroups, while those of 90 % of the patients belonged to macrohaplogroup R, which includes haplogroup H, the most common maternal lineage in Europe. Compared to the general US population, European haplogroups were significantly enriched in ES patients (p 2.2e-16) and the African haplogroups are significantly impoverished (p 4.6e-16). Using the ancestry informative markers defined in a National Genographic study, the vast majority of patients exhibited significant nuclear ancestry originating from the Mediterranean, Northern Europe, and Southwest Asia, including all six patients with African L mtDNAs. Very few had primarily African nuclear ancestry. This is the first genomic epidemiology study to simultaneously interrogate the mitochondrial and nuclear ancestries of ES patients. While supporting previous findings of enriched European ancestry in ES patients, these results also suggest alternative hypotheses for the significant contribution of mitochondrial ancestry in ES patients, as well as the protective role of African ancestry.

Supplementary Figure 6 show the distributions of the heteroplasmy levels of mtDNA mutations, depending on the gene and the location affected.

Little is known about the spectrum of mitochondrial DNA (mtDNA) mutations across pediatric malignancies. In this study, we analyzed matched tumor and normal whole genome sequencing data from 616 pediatric patients with hematopoietic malignancies, solid tumors, and brain tumors. We identified 391 mtDNA mutations in 284 tumors including 45 loss-of-function mutations, which clustered at four statistically significant hotspots in MT-COX3, MT-ND4, and MT-ND5, and at a mutation hotspot in MT-tRNA-MET. A skewed ratio (4.83) of nonsynonymous versus synonymous (dN/dS) mtDNA mutations with high statistical significance was identified on the basis of Monte Carlo simulations in the tumors. In comparison, opposite ratios of 0.44 and 0.93 were observed in 616 matched normal tissues and in 249 blood samples from children without cancer, respectively. mtDNA mutations varied by cancer type and mtDNA haplogroup. Collectively, these results suggest that deleterious mtDNA mutations play a role in the development and progression of pediatric cancers.Significance:This pan-cancer mtDNA study establishes the landscape of germline and tumor mtDNA mutations and identifies hotspots of tumor mtDNA mutations to pinpoint key mitochondrial functions in pediatric malignancies.

Supplementary Figure 3 illustrates how the observed dN/dS ratio varies by the heteroplasmy.

Supplementary Figure 2 illustrates that individual tumor cells showing the same mtDNA variants at different heteroplasmy levels as exhibited by the bulk tumor samples.

Figure Legends for Supplementary Figures 1-6

Supplementary Figure 1 illustrates single-cell ATAC-seq data contains sufficient mtDNA reads for interrogating mtDNA variants.

Supplementary Figure 5 illustrates that mtDNA mutations in the mtDNA hypermutated pediatric tumor samples tend to occur at low heteroplasmy levels and are more likely to be benign in nature.

Supplementary Tables 1-18 shows 1) mtDNA mutations found in 616 pediatric cancer samples, 2) calculation of expected dN/dS ratio, 3) mtDNA mutations found in 2202 tumor samples studied by two previous studies, 4) the number of nuclear and mtDNA somatic mutations found in 134 randomly selected tumor samples, 5) mtDNA genome coverage for all tumor samples studies, 6) mtDNA genome coverage for all matched normal samples studies, 7) samples with low mtDNA coverage, 8) mtDNA genome coverage for non-cancer children, 9) mtDNA variants found in matched normal samples, 10) mtDNA variants found in non-cancer samples, 11) samples with known pathogenic mtDNA variants, 12) mitochondrial haplogroups for all patients studies, 13) number of cancer patients by mtDNA haplogroup, 14) number of patients of each haplogroup by cancer subtype, 15) permutation results from meta-analysis, 16) meta-analysis of long homopolymer regions, 17) meta-analysis results of mtDNA tRNA mutations, 18) mtDNA mutations found in 5 mtDNA hypermutated cancer samples.

Mitochondrial DNA (mtDNA) variant pathogenicity interpretation has special considerations given unique features of the mtDNA genome, including maternal inheritance, variant heteroplasmy, threshold effect, absence of splicing, and contextual effects of haplogroups. Currently there are insufficient standardized criteria for mtDNA variant assessment, which leads to inconsistencies in clinical variant pathogenicity reporting. An international working group of mtDNA experts was assembled within the Mitochondrial Disease Sequence Data Resource (MSeqDR) Consortium and obtained Expert Panel status from ClinGen. This group reviewed the 2015 American College of Medical Genetics (ACMG) and Association of Molecular Pathology (AMP) standards and guidelines that are widely used for clinical interpretation of DNA sequence variants and provided further specifications for additional and specific guidance related to mtDNA variant classification. These Expert Panel based consensus specifications allow for consistent consideration of the unique aspects of the mtDNA genome that directly influence variant assessment, including addressing mtDNA genome composition and structure, haplogroups and phylogeny, maternal inheritance, heteroplasmy, and functional analyses unique to mtDNA, as well specifications for utilization of mtDNA genomic databases and computational algorithms.

The mitochondrial genome is small, 16.5 kb, and yet complex to study due to an abundance of mitochondria in any given cell or tissue. Mitochondrial DNA (mtDNA) mutations have been previously described in cancer, many of which were detected at low heteroplasmy. In this study we enriched the mitochondrial genome in primary pediatric tumors for detection of mtDNA variants. We completed mtDNA enrichment using REPLI-g, Agilent SureSelect, and long-range polymerase chain reaction (LRPCR) followed by next generation sequencing (NGS) on Illumina platforms. Primary tumor and germline genomic DNA from a variety of pediatric central nervous system (CNS) and extra-CNS solid tumors were analyzed by the three different methods. Although all three methods performed equally well for detecting variants at high heteroplasmy or homoplasmy, only LRPCR and SureSelect-based enrichment methods provided consistent results for variants that were present at less than five percent heteroplasmy. We then applied both LRPCR and SureSelect to three successive samples from a patient with multiply-recurrent gliofibroma and detected a low-level novel mutation as well as a change in heteroplasmy levels of a synonymous variant that was correlated with progression of disease. Implication This study demonstrates that LRPCR and SureSelect enrichment, but not REPLI-g, followed by NGS are accurate methods for studying the mtDNA variations at low heteroplasmy, which may be applied to studying mtDNA mutations in cancer.

Little is known about the spectrum of mitochondrial DNA (mtDNA) mutations across pediatric malignancies. In this study, we analyzed matched tumor and normal whole genome sequencing data from 616 pediatric patients with hematopoietic malignancies, solid tumors, and brain tumors. We identified 391 mtDNA mutations in 284 tumors including 45 loss-of-function mutations, which clustered at four statistically significant hotspots in MT-COX3, MT-ND4, and MT-ND5, and at a mutation hotspot in MT-tRNA-MET. A skewed ratio (4.83) of nonsynonymous versus synonymous (dN/dS) mtDNA mutations with high statistical significance was identified on the basis of Monte Carlo simulations in the tumors. In comparison, opposite ratios of 0.44 and 0.93 were observed in 616 matched normal tissues and in 249 blood samples from children without cancer, respectively. mtDNA mutations varied by cancer type and mtDNA haplogroup. Collectively, these results suggest that deleterious mtDNA mutations play a role in the development and progression of pediatric cancers. Significance: This pan-cancer mtDNA study establishes the landscape of germline and tumor mtDNA mutations and identifies hotspots of tumor mtDNA mutations to pinpoint key mitochondrial functions in pediatric malignancies.

Background We previously established the landscape of mitochondrial DNA (mtDNA) mutations in 23 subtypes of pediatric malignancies, characterized mtDNA mutation profiles among these subtypes, and provided statistically significant evidence for a contributory role of mtDNA mutations to pediatric malignancies. Methods To further delineate the spectrum of mtDNA mutations in pediatric central nervous system (CNS) tumors, we analyzed 545 tumor-normal paired whole-genome sequencing datasets from the Children’s Brain Tumor Tissue Consortium. Results Germline mtDNA variants were used to determine the haplogroup, and maternal ancestry, which was not significantly different among tumor types. Among 166 (30.5%) tumors we detected 220 somatic mtDNA mutations, primarily missense mutations (36.8%), as well as 22 loss-of-function mutations. Different pediatric CNS tumor subtypes had distinct mtDNA mutation profiles. The number of mtDNA mutations per tumor ranged from 0.20 (dysembryoplastic neuroepithelial tumor [DNET]) to 0.75 (meningiomas). The average heteroplasmy was 10.7%, ranging from 4.6% in atypical teratoid/rhabdoid tumor (AT/RT) to 26% in diffuse intrinsic pontine glioma. High-grade gliomas had a significant higher number of mtDNA mutations per sample than low-grade gliomas (0.6 vs 0.27) (P = .004), with almost twice as many missense mtDNA mutations per sample (0.24 vs 0.11), and higher average heteroplasmy levels (16% vs 10%). Recurrent mtDNA mutations may represent hotspots which may serve as biologic markers of disease. Conclusions Our findings demonstrate varying contributions of mtDNA mutations in different subtypes of CNS tumors. Sequencing the mtDNA genome may ultimately be used to characterize CNS tumors at diagnosis and monitor disease progression.

Background We previously established the landscape of mitochondrial DNA (mtDNA) mutations in 23 subtypes of pediatric malignancies, characterized mtDNA mutation profiles among these subtypes, and provided statistically significant evidence for a contributory role of mtDNA mutations to pediatric malignancies. Methods To further delineate the spectrum of mtDNA mutations in pediatric CNS tumors, we analyzed 545 tumor-normal paired whole genome sequencing data sets from the Children’s Brain Tumor Tissue Consortium. Results Germline mtDNA variants were used to determine the haplogroup, and maternal ancestry, which was not significantly different among tumor types. Among 166 (30.5%) tumors we detected 220 somatic mtDNA mutations, primarily missense mutations (36.8%), as well as 22 loss-of-function mutations. Different pediatric CNS tumor subtypes had distinct mtDNA mutation profiles. The number of mtDNA mutations per tumor ranged from 0.20 (dysembryoplastic neuroepithelial tumor) to 0.75 (meningiomas). The average heteroplasmy was 10.7%, ranging from 4.6% in atypical teratoid/rhabdoid tumor (AT/RT) to 26% in diffuse intrinsic pontine glioma. High-grade gliomas had a significant higher number of mtDNA mutations per sample than low-grade gliomas (0.6 vs. 0.27) (p = 0.004), with almost twice as many missense mtDNA mutations per sample (0.24 vs. 0.11), and higher average heteroplasmy levels (16% vs. 10%). Recurrent mtDNA mutations may represent hotspots which may serve as biologic markers of disease. Conclusions Our findings demonstrate varying contributions of mtDNA mutations in different subtypes of CNS tumors. Sequencing the mtDNA genome may ultimately be used to characterize CNS tumors at diagnosis and monitor disease progression.

Human induced pluripotent and embryonic stem cell cultures (hiPSC/hESC) are phenotypically heterogeneous and prone to clonal deviations during subculturing and differentiation. Clonal deviations often emerge unnoticed, but they can change the biology of the cell culture with a negative impact on experimental reproducibility. Here, we describe a computational workflow to profile the bulk clonal composition in a hiPSC/hESC culture that can also be used to infer clonal deviations. This workflow processes data obtained with two versions of the same method. The two versions-epigenetic and transcriptomic-rely on a mechanism of stochastic H3K4me3 deposition during hiPSC/hESC derivation. This mechanism generates a signature of ten or more H3K4me3-enriched clustered protocadherin (PCDH) promoters distinct in every single cell. The aggregate of single-cell signatures provides an identificatory feature in every hiPSC/hESC line. This feature is stably transmitted to the cell progeny of the culture even after differentiation unless there is a clonal deviation event that changes the internal balance of single-cell signatures. H3K4me3 signatures can be profiled by chromatin immunoprecipitation and next-generation sequencing (ChIP-seq). Alternatively, an equivalent PCDH-expression version can be profiled by RNA-seq in PCDH-expressing hiPSC/hESC-derived cells (such as neurons, astrocytes, and cardiomyocytes; and, in long-term cultures, such as cerebral organoids). Notably, our workflow can also distinguish genetically identical hiPSC/hESC lines derived from the same patient or generated in the same editing process. Together, we propose a method to improve data sharing and reproducibility in the hiPSC and hESC fields.

Summary: Cisplatin is an antineoplastic drug administered at suboptimal and intermittent doses to avoid life-threatening effects. Although this regimen shortly improves symptoms in the short term, it also leads to more malignant disease in the long term. We describe a multilayered analysis ranging from chromatin to translation—integrating chromatin immunoprecipitation sequencing (ChIP-seq), global run-on sequencing (GRO-seq), RNA sequencing (RNA-seq), and ribosome profiling—to understand how cisplatin confers (pre)malignant features by using a well-established ovarian cancer model of cisplatin exposure. This approach allows us to segregate the human transcriptome into gene modules representing distinct regulatory principles and to characterize that the most cisplatin-disrupted modules are associated with underlying events of super-enhancer plasticity. These events arise when cancer cells initiate without ultimately ending the program of drug-stimulated death. Using a PageRank-based algorithm, we predict super-enhancer regulator ISL1 as a driver of this plasticity and validate this prediction by using CRISPR/dCas9-KRAB inhibition (CRISPRi) and CRISPR/dCas9-VP64 activation (CRISPRa) tools. Together, we propose that cisplatin reprograms cancer cells when inducing them to undergo near-to-death experiences. : Ovarian tumors acquire new phenotypic features when exposed to anticancer therapy, which helps them endure the struggle for survival and increase malignancy. Using a six-layered multiomic approach, Ma et al. find that selective super-enhancer redistribution is a major driver of the acquisition of new phenotypic features. Keywords: multi-omic approach, gene dysregulation, super-enhancer, cisplatin, drug resistance, ovarian cancer, GRO-seq, ribosome profiling, Seurat, PageRank

A localized transcriptome at the synapse facilitates synapse-, stimulus- and transcript-specific local protein synthesis in response to neuronal activity. While enzyme-mediated mRNA modifications are known to regulate cellular mRNA turnover, the role of these modifications in regulating synaptic RNA has not been studied. We established low-input m6A-sequencing of synaptosomal RNA to determine the chemically modified local transcriptome in healthy adult mouse forebrains and identified 4,469 selectively enriched m6A sites in 2,921 genes as the synaptic m6A epitranscriptome (SME). The SME is functionally enriched in synthesis and modulation of tripartite synapses and in pathways implicated in neurodevelopmental and neuropsychiatric diseases. Interrupting m6A-mediated regulation via knockdown of readers in hippocampal neurons altered expression of SME member Apc, resulting in synaptic dysfunction including immature spine morphology and dampened excitatory synaptic transmission concomitant with decreased clusters of postsynaptic density-95 (PSD-95) and decreased surface expression of AMPA receptor subunit GluA1. Our findings indicate that chemical modifications of synaptic mRNAs critically contribute to synaptic function.

Leukemia cells rely on two nucleotide biosynthetic pathways, de novo and salvage, to produce dNTPs for DNA replication. Here, using metabolomic, proteomic, and phosphoproteomic approaches, we show that inhibition of the replication stress sensing kinase ataxia telangiectasia and Rad3-related protein (ATR) reduces the output of both de novo and salvage pathways by regulating the activity of their respective rate-limiting enzymes, ribonucleotide reductase (RNR) and deoxycytidine kinase (dCK), via distinct molecular mechanisms. Quantification of nucleotide biosynthesis in ATR-inhibited acute lymphoblastic leukemia (ALL) cells reveals substantial remaining de novo and salvage activities, and could not eliminate the disease in vivo. However, targeting these remaining activities with RNR and dCK inhibitors triggers lethal replication stress in vitro and long-term disease-free survival in mice with B-ALL, without detectable toxicity. Thus the functional interplay between alternative nucleotide biosynthetic routes and ATR provides therapeutic opportunities in leukemia and potentially other cancers., Leukemic cells depend on the nucleotide synthesis pathway to proliferate. Here the authors use metabolomics and proteomics to show that inhibition of ATR reduced the activity of these pathways thus providing a valuable therapeutic target in leukemia.

To explore the role of mitochondrial DNA mutations in pediatric CNS tumors, we analyzed 749 tumor-normal paired whole genome sequencing data sets from the Children’s Brain Tumor Tissue Consortium (CBTTC). We detected 307 somatic mtDNA mutations in 222 CNS tumors (29.6%). Most frequently observed were missense mutations (38.1%). We also detected 34 loss-of-function mutations. Different pediatric CNS tumor subtypes have distinct mtDNA mutation profiles. For categorical comparisons, we analyzed subtypes with at least 15 samples. The highest number of mtDNA mutations per tumor sample was in meningiomas (0.85), while atypical teratoid rhabdoid tumors (ATRTs) had the lowest number per sample (0.18). High-grade gliomas had a higher number of mtDNA mutations per sample than low-grade gliomas (0.56 vs. 0.31) (p = 0.0011), with almost twice as many missense mtDNA mutations per sample (0.22 vs. 0.13) (p < 0.001), and higher average heteroplasmy levels (11% vs. 9%). The average heteroplasmy was 10.1%, ranging from 15.6% in medulloblastoma to 6.36% in schwannoma suggesting that these are clonal alterations and not artifacts. Intriguingly, the two chordoma patients in the CBTTC database had an identical heteroplasmic m.10971G>A MT-ND4 nonsense mutation. Similarly, our patient with recurrent gliofibroma harbored the same somatic MT-ND4 synonymous variant (m.10700A>G) detected at 53% heteroplasmy in the initial tumor, 79% in the first recurrence, and 97% in the second recurrence. Although the functional consequences of these alterations are not yet understood, our findings suggest that sequencing the mtDNA genome may be used to characterize CNS tumors at diagnosis and monitor disease progression.

In the mammalian genome, the clustered protocadherin (cPcdh) locus is a paradigm of stochastic gene expression with the potential to generate a unique cPcdh combination in every neuron. Here, we report a chromatin-based mechanism emerging during the transition from the naive to the primed states of cell pluripotency that reduces by orders of magnitude the combinatorial potential in the human cPcdh locus. This mechanism selectively increases the frequency of stochastic selection of a small subset of cPcdh genes after neuronal differentiation in monolayers, months-old organoids, and engrafted cells in the rat spinal cord. Signs of these frequent selections can be observed in the brain throughout fetal development and disappear after birth, unless there is a condition of delayed maturation such as Down Syndrome. We therefore propose that a pattern of limited cPcdh diversity is maintained while human neurons still retain fetal-like levels of maturation., Editorial SUMMARY: Short and long-term cultures of human stem cell-derived neurons reveal that a pattern of restricted selection of clustered protocadherin isoforms, pre-established in pluripotent cells, distinguishes immature from mature neurons.

Lamin A/C (LMNA) gene mutations are a known cause of familial dilated cardiomyopathy, but the precise mechanisms triggering disease progression remain unknown. We hypothesize that analysis of differentially expressed genes (DEGs) throughout the course of Lmna knockout (Lmna−/−)-induced cardiomyopathy may reveal novel Lmna-mediated alterations of signaling pathways leading to dilated cardiomyopathy. Although Lmna was the only DEG down-regulated at 1 week of age, we identified 730 and 1004 DEGs in Lmna−/− mice at 2 weeks and 1 month of age, respectively. At 2 weeks, Lmna−/− mice demonstrated both down- and up-regulation of the key genes involving cell cycle control, mitochondrial dysfunction, and oxidative phosphorylation, as well as down-regulated genes governing DNA damage repair and up-regulated genes involved in oxidative stress response, cell survival, and cardiac hypertrophy. At 1 month, the down-regulated genes included those involved in oxidative phosphorylation, mitochondrial dysfunction, nutrient metabolism, cardiac β-adrenergic signaling, action potential generation, and cell survival. We also found 96 overlapping DEGs at both ages involved in oxidative phosphorylation, mitochondrial function, and calcium signaling. Impaired oxidative phosphorylation was observed at early disease stage, even before the appearance of disease phenotypes, and worsened with disease progression, suggesting its importance in the pathogenesis and progression of LMNA cardiomyopathy. Reduction of oxidative stress might therefore prevent or delay the development from Lmna mutation to LMNA cardiomyopathy.

As a key component of the vertebrate neuroendocrine system, the pituitary gland relies on the progressive and coordinated development of distinct hormone-producing cell types and an invading vascular network. The molecular mechanisms that drive formation of the pituitary vasculature, which is necessary for regulated synthesis and secretion of hormones that maintain homeostasis, metabolism, and endocrine function, remain poorly understood. Here, we report that expression of integrin β1 in embryonic pituitary epithelial cells is required for angiogenesis in the developing mouse pituitary gland. Deletion of pituitary epithelial integrin β1 before the onset of angiogenesis resulted in failure of invading endothelial cells to recruit pericytes efficiently, whereas deletion later in embryogenesis led to decreased vascular density and lumen formation. In both cases, lack of epithelial integrin β1 was associated with a complete absence of vasculature in the pituitary gland at birth. Within pituitary epithelial cells, integrin β1 directs a large transcriptional program that includes components of the extracellular matrix and associated signaling factors that are linked to the observed non-cell-autonomous effects on angiogenesis. We conclude that epithelial integrin β1 functions as a critical and canonical regulator of developmental angiogenesis in the pituitary gland, thus providing insight into the long-standing systems biology conundrum of how vascular invasion is coordinated with tissue development.

Chordoma is a rare tumor and while SMARCB1 alterations have been observed in poorly differentiated chordomas, conventional chordomas are not well understood. We interrogated nuclear and mitochondrial genomes of 11 chordoma samples from 7 children. Frozen tumor tissue DNA was extracted and whole exome libraries generated using Agilent SureSelect Human All Exon V6 kit plus mtDNA genome capture kit. Libraries were sequenced using Illumina Nextseq 500. MuTect2, VarDict and LUBA variant callers were used with allele frequency cutoff 2%. Potential germline variants were filtered bioinformatically. In total, 656±74 high-confidence somatic variants, including 368±43 nonsynonymous variants per sample were detected. Of 2,607 combined unique nonsynonymous variants, 95% were missense. Remaining high impact variants were frameshift (37%), stop gain (39%), splice acceptor/donor (22%), start and stop loss (2%). Of the unique nonsynonymous variants, 137 fall within Cosmic Cancer Census Genes, including high impact variants in SETD2, MLLT4. No previously reported TBXT, CDKN2A, PI3K, LYST mutations identified. Tumor Mutation Burden/Megabase was 10±1. The mitochondrial analysis revealed heteroplasmic m.11727C>T MT-ND4 missense variants in three tumors resected at different time points from the same patient, and another heteroplasmic m.1023C>T rRNA mutation from the primary and recurrent tumors of another patient. Intriguingly, two Children’s Brain Tumor Tissue Consortium patients with chordoma had identical heteroplasmic m.10971G>A MT-ND4 nonsense mutations. Pediatric chordomas appear to lack somatic nuclear mutations. Observing recurrent mitochondrial mutations across multiple tumors from the same and/or different patients is striking, suggesting they may be implicated in tumorigenesis and be potential diagnostic markers.

Anticancer drugs with narrow therapeutic window, such as cisplatin, are administrated at suboptimal and intermittent doses to limit their adverse effects in patients. This regimen may allow survival of cancer cells that problematically also gain malignancy as result of the treatment. Here, we describe a multilayered analysis ranging from chromatin to translation (based on ChIP-seq, GRO-seq, RNA-seq, and ribosome profiling) to comprehensively understand how cisplatin confers new cellular properties in a well-established ovarian cancer model of cisplatin exposure. This approach allowed us to segregate the human transcriptome into seventeen modules underlying distinct regulatory principles. The most robustly disrupted modules are primarily associated with super-enhancer reprogramming, induced after cells repeatedly initiate but fail to complete the process of drug-stimulated death (corroborated by CRISPRi/a technology after PageRank-based prediction of candidate mediators). We therefore propose that cisplatin provides survivor cells with new properties, paradoxically, when they undergo a drug-stimulated near to death experience.

© 2018 Elsevier Inc. Tools to understand how the spliceosome functions in vivo have lagged behind advances in the structural biology of the spliceosome. Here, methods are described to globally profile spliceosome-bound pre-mRNA, intermediates, and spliced mRNA at nucleotide resolution. These tools are applied to three yeast species that span 600 million years of evolution. The sensitivity of the approach enables the detection of canonical and non-canonical events, including interrupted, recursive, and nested splicing. This application of statistical modeling uncovers independent roles for the size and position of the intron and the number of introns per transcript in substrate progression through the two catalytic stages. These include species-specific inputs suggestive of spliceosome-transcriptome coevolution. Further investigations reveal the ATP-dependent discard of numerous endogenous substrates after spliceosome assembly in vivo and connect this discard to intron retention, a form of splicing regulation. Spliceosome profiling is a quantitative, generalizable global technology used to investigate an RNP central to eukaryotic gene expression. Quantitative in vivo spliceosome profiling in three divergent yeast species reveals unannotated splicing events and determinants of catalytic efficiency and assesses non-optimal substrates.

Hematopoietic stem cells (HSCs) maintain a quiescent state during homeostasis, but with acute infection, they exit the quiescent state to increase the output of immune cells, the so-called "emergency hematopoiesis." However, HSCs' response to severe infection during septic shock and the pathological impact remain poorly elucidated. Here, we report that the histone demethylase KDM1A/LSD1, serving as a critical regulator of mammalian hematopoiesis, is a negative regulator of the response to inflammation in HSCs during endotoxic shock typically observed during acute bacterial or viral infection. Inflammation-induced LSD1 deficiency results in an acute expansion of a pathological population of hyperproliferative and hyperinflammatory myeloid progenitors, resulting in a septic shock phenotype and acute death. Unexpectedly, in vivo administration of bacterial lipopolysaccharide (LPS) to wild-type mice results in acute suppression of LSD1 in HSCs with a septic shock phenotype that resembles that observed following induced deletion of LSD1 The suppression of LSD1 in HSCs is caused, at least in large part, by a cohort of inflammation-induced microRNAs. Significantly, reconstitution of mice with bone marrow progenitor cells expressing inhibitors of these inflammation-induced microRNAs blocked the suppression of LSD1 in vivo following acute LPS administration and prevented mortality from endotoxic shock. Our results indicate that LSD1 activators or miRNA antagonists could serve as a therapeutic approach for life-threatening septic shock characterized by dysfunction of HSCs.

SummaryTools to understand how the spliceosome functionsin vivohave lagged behind advances in its structural biology. We describe methods to globally profile spliceosome-bound precursor, intermediates and products at nucleotide resolution. We apply these tools to three divergent yeast species that span 600 million years of evolution. The sensitivity of the approach enables detection of novel cases of non-canonical catalysis including interrupted, recursive and nested splicing. Employing statistical modeling to understand the quantitative relationships between RNA features and the data, we uncover independent roles for intron size, position and number in substrate progression through the two catalytic stages. These include species-specific inputs suggestive of spliceosome-transcriptome coevolution. Further investigations reveal ATP-dependent discard of numerous endogenous substrates at both the precursor and lariat-intermediate stages and connect discard to intron retention, a form of splicing regulation. Spliceosome profiling is a quantitative, generalizable global technology to investigate an RNP central to eukaryotic gene expression.HighlightsMeasurement of spliceosome-bound precursor and intermediate in three speciesNon-canonical splicing events revealedStatistical modeling uncovers substrate features that predict catalytic efficiencyDiscard of suboptimal substrates occursin vivoand predicts intron-retained mRNAs

A localized transcriptome at the synapse facilitates synapse-, stimulus-, and transcript-specific synthesis of the local proteome in response to neuronal activity. While enzyme-mediated mRNA modifications have been shown to regulate cellular mRNA turnover and translation, the role of these modifications in regulating synaptic RNA has not been studied. We established low-input m6A-seq of synaptosomal RNA to determine the chemically modified local transcriptome in healthy adult mouse forebrain and identified 4,329 selectively enriched m6A RNA peaks in 2,987 genes, which we refer to as the synaptic m6A epitranscriptome (SME). SME is functionally enriched in synthesis and modulation of tripartite synapses, and in pathways implicated in neurodevelopmental and neuropsychiatric diseases. Interrupting m6A-mediated regulation via knockdown of reader YTHDF1 in hippocampal neurons alters expression of SME member Apc, and causes synaptic malfunctions manifesting immature spine morphology and dampened excitatory synaptic transmission concomitant with decreased PSD-95 clustering and GluA1 surface expression. Our findings indicate that chemical modifications of synaptic mRNAs critically contribute to synaptic function.

A localized transcriptome at the synapse facilitates synapse-, stimulus- and transcript-specific local protein synthesis in response to neuronal activity. While enzyme-mediated mRNA modifications are known to regulate cellular mRNA turnover, the role of these modifications in regulating synaptic RNA has not been studied. We established low-input m

The earliest event in Xenopus development is the dorsal accumulation of nuclear β-catenin under the influence of cytoplasmic determinants displaced by fertilization. In this study, a genome-wide approach was used to examine transcription of the 43,673 genes annotated in the Xenopus laevis genome under a variety of conditions that inhibit or promote formation of the Spemann organizer signaling center. Loss of function of β-catenin with antisense morpholinos reproducibly reduced the expression of 247 mRNAs at gastrula stage. Interestingly, only 123 β-catenin targets were enriched on the dorsal side and defined an early dorsal β-catenin gene signature. These genes included several previously unrecognized Spemann organizer components. Surprisingly, only 3 of these 123 genes overlapped with the late Wnt signature recently defined by two other groups using inhibition by Dkk1 mRNA or Wnt8 morpholinos, which indicates that the effects of β-catenin/Wnt signaling in early development are exquisitely regulated by stage-dependent mechanisms. We analyzed transcriptome responses to a number of treatments in a total of 46 RNA-seq libraries. These treatments included, in addition to β-catenin depletion, regenerating dorsal and ventral half-embryos, lithium chloride treatment, and the overexpression of Wnt8, Siamois, and Cerberus mRNAs. Only some of the early dorsal β-catenin signature genes were activated at blastula whereas others required the induction of endomesoderm, as indicated by their inhibition by Cerberus overexpression. These comprehensive data provide a rich resource for analyzing how the dorsal and ventral regions of the embryo communicate with each other in a self-organizing vertebrate model embryo.

© 2014 Macmillan Publishers Limited. Post-translational histone modifications have a critical role in regulating transcription, the cell cycle, DNA replication and DNA damage repair. The identification of new histone modifications critical for transcriptional regulation at initiation, elongation or termination is of particular interest. Here we report a new layer of regulation in transcriptional elongation that is conserved from yeast to mammals. This regulation is based on the phosphorylation of a highly conserved tyrosine residue, Tyr 57, in histone H2A and is mediated by the unsuspected tyrosine kinase activity of casein kinase 2 (CK2). Mutation of Tyr 57 in H2A in yeast or inhibition of CK2 activity impairs transcriptional elongation in yeast as well as in mammalian cells. Genome-wide binding analysis reveals that CK2α, the catalytic subunit of CK2, binds across RNA-polymerase-II-transcribed coding genes and active enhancers. Mutation of Tyr 57 causes a loss of H2B mono-ubiquitination as well as H3K4me3 and H3K79me3, histone marks associated with active transcription. Mechanistically, both CK2 inhibition and the H2A(Y57F) mutation enhance H2B deubiquitination activity of the Spt-Ada-Gcn5 acetyltransferase (SAGA) complex, suggesting a critical role of this phosphorylation in coordinating the activity of the SAGA complex during transcription. Together, these results identify a new component of regulation in transcriptional elongation based on CK2-dependent tyrosine phosphorylation of the globular domain of H2A.

Understanding the mechanisms by which compounds discovered using cell-based phenotypic screening strategies might exert their effects would be highly augmented by new approaches exploring their potential interactions with the genome. For example, altered androgen receptor (AR) transcriptional programs, including castration resistance and subsequent chromosomal translocations, play key roles in prostate cancer pathological progression, making the quest for identification of new therapeutic agents and an understanding of their actions a continued priority. Here we report an approach that has permitted us to uncover the sites and mechanisms of action of a drug, referred to as "SD70," initially identified by phenotypic screening for inhibitors of ligand and genotoxic stress-induced translocations in prostate cancer cells. Based on synthesis of a derivatized form of SD70 that permits its application for a ChIP-sequencing-like approach, referred to as "Chem-seq," we were next able to efficiently map the genome-wide binding locations of this small molecule, revealing that it largely colocalized with AR on regulatory enhancers. Based on these observations, we performed the appropriate global analyses to ascertain that SD70 inhibits the androgen-dependent AR program, and prostate cancer cell growth, acting, at least in part, by functionally inhibiting the Jumonji domain-containing demethylase, KDM4C. Global location of candidate drugs represents a powerful strategy for new drug development by mapping genome-wide location of small molecules, a powerful adjunct to contemporary drug development strategies.

Although recent studies have indicated roles of long non-coding RNAs (lncRNAs) in physiological aspects of cell-type determination and tissue homeostasis, their potential involvement in regulated gene transcription programs remains rather poorly understood. The androgen receptor regulates a large repertoire of genes central to the identity and behaviour of prostate cancer cells, and functions in a ligand-independent fashion in many prostate cancers when they become hormone refractory after initial androgen deprivation therapy. Here we report that two lncRNAs highly overexpressed in aggressive prostate cancer, PRNCR1 (also known as PCAT8) and PCGEM1, bind successively to the androgen receptor and strongly enhance both ligand-dependent and ligand-independent androgen-receptor-mediated gene activation programs and proliferation in prostate cancer cells. Binding of PRNCR1 to the carboxy-terminally acetylated androgen receptor on enhancers and its association with DOT1L appear to be required for recruitment of the second lncRNA, PCGEM1, to the androgen receptor amino terminus that is methylated by DOT1L. Unexpectedly, recognition of specific protein marks by PCGEM1-recruited pygopus 2 PHD domain enhances selective looping of androgen-receptor-bound enhancers to target gene promoters in these cells. In 'resistant' prostate cancer cells, these overexpressed lncRNAs can interact with, and are required for, the robust activation of both truncated and full-length androgen receptor, causing ligand-independent activation of the androgen receptor transcriptional program and cell proliferation. Conditionally expressed short hairpin RNA targeting these lncRNAs in castration-resistant prostate cancer cell lines strongly suppressed tumour xenograft growth in vivo. Together, these results indicate that these overexpressed lncRNAs can potentially serve as a required component of castration-resistance in prostatic tumours.

Enhancers for embryonic stem (ES) cell-expressed genes and lineage-determining factors are characterized by conventional marks of enhancer activation in ES cells1-3, but it remains unclear whether enhancers destined to regulate cell-type-restricted transcription units might also have distinct signatures in ES cells. Here we show that cell-type-restricted enhancers are 'premarked' and activated as transcription units by the binding of one or two ES cell transcription factors, although they do not exhibit traditional enhancer epigenetic marks in ES cells, thus uncovering the initial temporal origins of cell-type-restricted enhancers. This premarking is required for future cell-type-restricted enhancer activity in the differentiated cells, with the strength of the ES cell signature being functionally important for the subsequent robustness of cell-type-restricted enhancer activation. We have experimentally validated this model in macrophage-restricted enhancers and neural precursor cell (NPC)-restricted enhancers using ES cell-derived macrophages or NPCs, edited to contain specific ES cell transcription factor motif deletions. DNA hydroxyl-methylation of enhancers in ES cells, determined by ES cell transcription factors, may serve as a potential molecular memory for subsequent enhancer activation in mature macrophages. These findings suggest that the massive repertoire of cell-type-restricted enhancers are essentially hierarchically and obligatorily premarked by binding of a defining ES cell transcription factor in ES cells, dictating the robustness of enhancer activation in mature cells.

In the version of this article initially published, a Supplementary Fig. 6f was cited in the last paragraph of the Results. No such panel exists; the citation has been deleted. The error has been corrected in the HTML and PDF versions of the article.

Enhancers instruct spatio-temporally specific gene expression in a manner tightly linked to higher-order chromatin architecture. Critical chromatin architectural regulators condensin I and condensin II play non-redundant roles controlling mitotic chromosomes. But the chromosomal locations of condensins and their functional roles in interphase are poorly understood. Here we report that both condensin complexes exhibit an unexpected, dramatic estrogen-induced recruitment to estrogen receptor α (ER-α)-bound eRNA(+) active enhancers in interphase breast cancer cells, exhibiting non-canonical interaction with ER-α via its DNA-binding domain (DBD). Condensins positively regulate ligand-dependent enhancer activation at least in part by recruiting an E3 ubiquitin ligase, HECTD1, to modulate the binding of enhancer-associated coactivators/corepressors, including p300 and RIP140, permitting full eRNA transcription, formation of enhancer:promoter looping, and the resultant coding gene activation. Collectively, our results reveal an important, unanticipated transcriptional role of interphase condensins in modulating estrogen-regulated enhancer activation and coding gene transcriptional program.

© 2015 Taylor and Francis Group, LLC. All eukaryotic cells utilize autophagy for protein and organelle turnover, thus assuring subcellular quality control, homeostasis, and survival. In order to address recent advances in identification of human autophagy associated genes, and to describe autophagy on a system-wide level, we established an autophagy-centered gene interaction network by merging various primary data sets and by retrieving respective interaction data. The resulting network (‘AXAN’) was analyzed with respect to subnetworks, e.g. the prime gene subnetwork (including the core machinery, signaling pathways and autophagy receptors) and the transcription subnetwork. To describe aspects of evolution within this network, we assessed the presence of protein orthologs across 99 eukaryotic model organisms. We visualized evolutionary trends for prime gene categories and evolutionary tracks for selected AXAN genes. This analysis confirms the eukaryotic origin of autophagy core genes while it points to a diverse evolutionary history of autophagy receptors. Next, we used module identification to describe the functional anatomy of the network at the level of pathway modules. In addition to obvious pathways (e.g., lysosomal degradation, insulin signaling) our data unveil the existence of context-related modules such as Rho GTPase signaling. Last, we used a tripartite, image-based RNAi – screen to test candidate genes predicted to play a role in regulation of autophagy. We verified the Rho GTPase, CDC42, as a novel regulator of autophagy-related signaling. This study emphasizes the applicability of system-wide approaches to gain novel insights into a complex biological process and to describe the human autophagy pathway at a hitherto unprecedented level of detail.

Substantial evidence supports the hypothesis that enhancers are critical regulators of cell-type determination, orchestrating both positive and negative transcriptional programs; however, the basic mechanisms by which enhancers orchestrate interactions with cognate promoters during activation and repression events remain incompletely understood. Here we report the required actions of LIM domain-binding protein 1 (LDB1)/cofactor of LIM homeodomain protein 2/nuclear LIM interactor, interacting with the enhancer-binding protein achaete-scute complex homolog 1, to mediate looping to target gene promoters and target gene regulation in corticotrope cells. LDB1-mediated enhancer:promoter looping appears to be required for both activation and repression of these target genes. Although LDB1-dependent activated genes are regulated at the level of transcriptional initiation, the LDB1-dependent repressed transcription units appear to be regulated primarily at the level of promoter pausing, with LDB1 regulating recruitment of metastasis-associated 1 family, member 2, a component of the nucleosome remodeling deacetylase complex, on these negative enhancers, required for the repressive enhancer function. These results indicate that LDB1-dependent looping events can deliver repressive cargo to cognate promoters to mediate promoter pausing events in a pituitary cell type.

SummaryEnhancers provide critical information directing cell-type-specific transcriptional programs, regulated by binding of signal-dependent transcription factors and their associated cofactors. Here, we report that the most strongly activated estrogen (E2)-responsive enhancers are characterized by trans-recruitment and in situ assembly of a large 1–2 MDa complex of diverse DNA-binding transcription factors by ERα at ERE-containing enhancers. We refer to enhancers recruiting these factors as mega transcription factor-bound in trans (MegaTrans) enhancers. The MegaTrans complex is a signature of the most potent functional enhancers and is required for activation of enhancer RNA transcription and recruitment of coactivators, including p300 and Med1. The MegaTrans complex functions, in part, by recruiting specific enzymatic machinery, exemplified by DNA-dependent protein kinase. Thus, MegaTrans-containing enhancers represent a cohort of functional enhancers that mediate a broad and important transcriptional program and provide a molecular explanation for transcription factor clustering and hotspots noted in the genome.

The functional importance of gene enhancers in regulated gene expression is well established. In addition to widespread transcription of long non-coding RNAs (lncRNAs) in mammalian cells, bidirectional ncRNAs are transcribed on enhancers, and are thus referred to as enhancer RNAs (eRNAs). However, it has remained unclear whether these eRNAs are functional or merely a reflection of enhancer activation. Here we report that in human breast cancer cells 17β-oestradiol (E2)-bound oestrogen receptor α (ER-α) causes a global increase in eRNA transcription on enhancers adjacent to E2-upregulated coding genes. These induced eRNAs, as functional transcripts, seem to exert important roles for the observed ligand-dependent induction of target coding genes, increasing the strength of specific enhancer-promoter looping initiated by ER-α binding. Cohesin, present on many ER-α-regulated enhancers even before ligand treatment, apparently contributes to E2-dependent gene activation, at least in part by stabilizing E2/ER-α/eRNA-induced enhancer-promoter looping. Our data indicate that eRNAs are likely to have important functions in many regulated programs of gene transcription.