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.

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.

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.

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.

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.

RARE-57. PEDIATRIC CHORDOMA: WHOLE EXOME SEQUENCING OF 11 PEDIATRIC CHORDOMA SAMPLES

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&gt;T MT-ND4 missense variants in three tumors resected at different time points from the same patient, and another heteroplasmic m.1023C&gt;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&gt;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.

TBIO-12. THE SPECTRUM OF MITOCHONDRIAl DNA (mtDNA) MUTATIONS IN PEDIATRIC CENTRAL NERVOUS SYSTEM (CNS) TUMORS

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.

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.

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.

Landscape of Germline and Somatic Mitochondrial DNA Mutations in Pediatric Malignancies

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.

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.

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.

Evolutionary trends and functional anatomy of the human expanded autophagy network

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.

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.

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

While recent studies indicated roles of long non-coding RNAs (lncRNAs) in physiologic aspects of cell-type determination and tissue homeostasis¹ yet their potential involvement in regulated gene transcription programs remain rather poorly understood. Androgen receptor (AR) regulates a large repertoire of genes central to the identity and behavior 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 and PCGEM1, bind successively to the AR and strongly enhance both ligand-dependent and ligand-independent AR-mediated gene activation programs and proliferation in prostate cancer cells. Binding of PRNCR1 to the C-terminally acetylated AR on enhancers and its association with DOT1L appear to be required for recruitment of the second lncRNA, PCGEM1, to the DOT1L-mediated methylated AR N-terminus. Unexpectedly, recognition of specific protein marks by PCGEM1-recruited Pygopus2 PHD domain proves to enhance selective looping of AR-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 AR, causing ligand-independent activation of the AR transcriptional program and cell proliferation. Conditionally-expressed short hairpin RNA (shRNA) targeting of these lncRNAs in castration-resistant prostate cancer (CRPC) cell lines strongly suppressed tumor xenograft growth in vivo. Together, these results suggest that these overexpressed lncRNAs can potentially serve as a required component of castration-resistance in prostatic tumors.

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.