Spliceosome Profiling Visualizes Operations of a Dynamic RNP at Nucleotide Resolution.

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.

Enhancer activation requires trans-recruitment of a mega transcription factor complex.

Enhancers 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.

Pluripotency factors functionally premark cell-type-restricted enhancers in ES cells.

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.

Condensin I and II Complexes License Full Estrogen Receptor α-Dependent Enhancer Activation.

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.

Histone demethylase LSD1 regulates hematopoietic stem cells homeostasis and protects from death by endotoxic shock.

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.

Specifications of the ACMG/AMP standards and guidelines for mitochondrial DNA variant interpretation.

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 Consortium and obtained Expert Panel status from ClinGen. This group reviewed the 2015 American College of Medical Genetics and Association of Molecular Pathology 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 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 as specifications for utilization of mtDNA genomic databases and computational algorithms.

Tyrosine phosphorylation of histone H2A by CK2 regulates transcriptional elongation.

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.

Spemann organizer transcriptome induction by early beta-catenin, Wnt, Nodal, and Siamois signals in Xenopus laevis.

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.

lncRNA-dependent mechanisms of androgen-receptor-regulated gene activation programs.

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.

Functional roles of enhancer RNAs for oestrogen-dependent transcriptional activation.

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.

Epithelial cell integrin ß1 is required for developmental angiogenesis in the pituitary gland.

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.

Enhancer-bound LDB1 regulates a corticotrope promoter-pausing repression program.

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.

Chem-seq permits identification of genomic targets of drugs against androgen receptor regulation selected by functional phenotypic screens.

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.

Analysis of Clonal Composition in Human iPSC and ESC and Derived 2D and 3D Differentiated Cultures.

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.

Mitochondrial DNA haplogroup, genetic ancestry, and susceptibility to Ewing sarcoma.

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.

The spectrum of mitochondrial DNA (mtDNA) mutations in pediatric CNS tumors.

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.

Analysis of Brugada syndrome loci reveals that fine-mapping clustered GWAS hits enhances the annotation of disease-relevant variants.

Genome-wide association studies (GWASs) are instrumental in identifying loci harboring common single-nucleotide variants (SNVs) that affect human traits and diseases. GWAS hits emerge in clusters, but the focus is often on the most significant hit in each trait- or disease-associated locus. The remaining hits represent SNVs in linkage disequilibrium (LD) and are considered redundant and thus frequently marginally reported or exploited. Here, we interrogate the value of integrating the full set of GWAS hits in a locus repeatedly associated with cardiac conduction traits and arrhythmia, SCN5A-SCN10A. Our analysis reveals 5 common 7-SNV haplotypes (Hap1-5) with 2 combinations associated with life-threatening arrhythmia-Brugada syndrome (the risk Hap1/1 and protective Hap2/3 genotypes). Hap1 and Hap2 share 3 SNVs; thus, this analysis suggests that assuming redundancy among clustered GWAS hits can lead to confounding disease-risk associations and supports the need to deconstruct GWAS data in the context of haplotype composition.

RNA Sequence Analyses throughout the Course of Mouse Cardiac Laminopathy Identify Differentially Expressed Genes for Cell Cycle Control and Mitochondrial Function.

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.

Detection of mitochondrial DNA variants at low level heteroplasmy in pediatric CNS and extra-CNS solid tumors with three different enrichment methods.

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.

Super-Enhancer Redistribution as a Mechanism of Broad Gene Dysregulation in Repeatedly Drug-Treated Cancer Cells.

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.