Legends for Supplemental Figures S1-S5 and Supplemental Table S1.

Supplementary Table 3 from Integrative Genomics Identifies Distinct Molecular Classes of Neuroblastoma and Shows That Multiple Genes Are Targeted by Regional Alterations in DNA Copy Number

Supplementary Table 2 from Integrative Genomics Identifies Distinct Molecular Classes of Neuroblastoma and Shows That Multiple Genes Are Targeted by Regional Alterations in DNA Copy Number

Olig 2 is significantly reduced in Snf5F/F p53L/L GFAP-Cre adult mutant mice.

Description of additional methods and procedures used in the study.

Supplementary Table 3 from Integrative Genomics Identifies Distinct Molecular Classes of Neuroblastoma and Shows That Multiple Genes Are Targeted by Regional Alterations in DNA Copy Number

Primer sequences of qRT-PCR.

NG2 expression in Snf5F/F p53L/L GFAP-Cre mutant brain tissue.

Supplementary Table 1 from Integrative Genomics Identifies Distinct Molecular Classes of Neuroblastoma and Shows That Multiple Genes Are Targeted by Regional Alterations in DNA Copy Number

Snf5 inactivation causes neurodegeneration.

Supplementary Figures 1-2, Tables 1-2 from Serial Transcriptome Analysis and Cross-Species Integration Identifies <i>Centromere-Associated Protein E</i> as a Novel Neuroblastoma Target

Gene Set Enrichment Analysis shows close resemblance between human and mouse brain tumors.

Twelvepatients without therapy-related leukemia were studied after completing TOP2 poison chemotherapy in a high-risk neuroblastoma regimen. One patient harbored an inv(11) that was a KMT2A rearrangement. The KMT2A-MAML2 transcript was expressed at low level. The patient was prospectively followed. The inv(11) was undetectable in ensuing samples. Leukemia never developed after a 12.8-year follow-up period. Enriched etoposide-induced TOP2A cleavage in the relevant MAML2 genomic region supports a TOP2A DNA damage mechanism. After completing TOP2 poison chemotherapies, covert KMT2A-R clones may occur in a small minority of patients; however, not all KMT2A rearrangements herald a therapy-related leukemia diagnosis.

Single-cell gene expression analyses of mammalian tissues have uncovered profound stage-specific molecular regulatory phenomena that have changed the understanding of unique cell types and signaling pathways critical for lineage determination, morphogenesis, and growth. We discuss here the case for a Pediatric Cell Atlas as part of the Human Cell Atlas consortium to provide single-cell profiles and spatial characterization of gene expression across human tissues and organs. Such data will complement adult and developmentally focused HCA projects to provide a rich cytogenomic framework for understanding not only pediatric health and disease but also environmental and genetic impacts across the human lifespan.

Activins regulate numerous processes including inflammation and are synthesized as precursors consisting of a long N-terminal pro-region and a mature protein. Genomic human databases currently list three activin A (Act A) variants termed X1, X2 and X3. The X3 variant is the shortest, lacks N-terminal segments present in X1 and X2, and has been the focus of most past literature. Here, we asked whether these variants are expressed by human cells and tissues and what structural features are contained within their pro-regions. Human monocytic-like cells THP1 and U937 expressed X1 and X2 variants after exposure to phorbol ester or granulocyte-macrophage colony-stimulating factor, while X2 transcripts were present in placenta. Expression vectors encoding full length X2 or X3 variants resulted in production and secretion of biologically active Act A from cultured cells. Previous studies reported a putative HS-binding domain (HBD) in the X3 pro-region. Here, we identified a novel HBD with consensus HS-binding motifs near the N-terminal end of X1 and X2 proregions. Peptides encompassing this new domain interacted with substrate-bound HS with nanomolar affinity, while peptides from putative X3 HBD did not. In good agreement, full length X2 pro-region interacted with heparin-agarose, while the X3 pro-region did not. In sum, our study reveals that Act A variants are expressed by inflammatory cells and placenta and yield biological activity. The high affinity HBD in X1 and X2 pro-region and its absence in X3 could greatly influence overall Act A distribution, availability and activity in physiological and pathological circumstances.

Type II topoisomerases orchestrate proper DNA topology, and they are the targets of anti-cancer drugs that cause treatment-related leukemias with balanced translocations. Here, we develop a high-throughput sequencing technology to define TOP2 cleavage sites at single-base precision, and use the technology to characterize TOP2A cleavage genome-wide in the human K562 leukemia cell line. We find that TOP2A cleavage has functionally conserved local sequence preferences, occurs in cleavage cluster regions (CCRs), and is enriched in introns and lincRNA loci. TOP2A CCRs are biased toward the distal regions of gene bodies, and TOP2 poisons cause a proximal shift in their distribution. We find high TOP2A cleavage levels in genes involved in translocations in TOP2 poison–related leukemia. In addition, we find that a large proportion of genes involved in oncogenic translocations overall contain TOP2A CCRs. The TOP2A cleavage of coding and lincRNA genes is independently associated with both length and transcript abundance. Comparisons to ENCODE data reveal distinct TOP2A CCR clusters that overlap with marks of transcription, open chromatin, and enhancers. Our findings implicate TOP2A cleavage as a broad DNA damage mechanism in oncogenic translocations as well as a functional role of TOP2A cleavage in regulating transcription elongation and gene activation.

Reduced mitochondrial DNA copy number, mitochondrial DNA mutations or disruption of electron transfer chain complexes induce mitochondria-to-nucleus retrograde signaling, which induces global change in nuclear gene expression ultimately contributing to various human pathologies including cancer. Recent studies suggest that these mitochondrial changes cause transcriptional reprogramming of nuclear genes although the mechanism of this cross talk remains unclear. Here, we provide evidence that mitochondria-to-nucleus retrograde signaling regulates chromatin acetylation and alters nuclear gene expression through the heterogeneous ribonucleoprotein A2 (hnRNAP2). These processes are reversed when mitochondrial DNA content is restored to near normal cell levels. We show that the mitochondrial stress-induced transcription coactivator hnRNAP2 acetylates Lys 8 of H4 through an intrinsic histone lysine acetyltransferase (KAT) activity with Arg 48 and Arg 50 of hnRNAP2 being essential for acetyl-CoA binding and acetyltransferase activity. H4K8 acetylation at the mitochondrial stress-responsive promoters by hnRNAP2 is essential for transcriptional activation. We found that the previously described mitochondria-to-nucleus retrograde signaling-mediated transformation of C2C12 cells caused an increased expression of genes involved in various oncogenic processes, which is retarded in hnRNAP2 silenced or hnRNAP2 KAT mutant cells. Taken together, these data show that altered gene expression by mitochondria-to-nucleus retrograde signaling involves a novel hnRNAP2-dependent epigenetic mechanism that may have a role in cancer and other pathologies.

Background The ability to capture and sequence large contiguous DNA fragments represents a significant advancement towards the comprehensive characterization of complex genomic regions. While emerging sequencing platforms are capable of producing several kilobases-long reads, the fragment sizes generated by current DNA target enrichment technologies remain a limiting factor, producing DNA fragments generally shorter than 1 kbp. The DNA enrichment methodology described herein, Region-Specific Extraction (RSE), produces DNA segments in excess of 20 kbp in length. Coupling this enrichment method to appropriate sequencing platforms will significantly enhance the ability to generate complete and accurate sequence characterization of any genomic region without the need for reference-based assembly. Results RSE is a long-range DNA target capture methodology that relies on the specific hybridization of short (20-25 base) oligonucleotide primers to selected sequence motifs within the DNA target region. These capture primers are then enzymatically extended on the 3’-end, incorporating biotinylated nucleotides into the DNA. Streptavidin-coated beads are subsequently used to pull-down the original, long DNA template molecules via the newly synthesized, biotinylated DNA that is bound to them. We demonstrate the accuracy, simplicity and utility of the RSE method by capturing and sequencing a 4 Mbp stretch of the major histocompatibility complex (MHC). Our results show an average depth of coverage of 164X for the entire MHC. This depth of coverage contributes significantly to a 99.94 % total coverage of the targeted region and to an accuracy that is over 99.99 %. Conclusions RSE represents a cost-effective target enrichment method capable of producing sequencing templates in excess of 20 kbp in length. The utility of our method has been proven to generate superior coverage across the MHC as compared to other commercially available methodologies, with the added advantage of producing longer sequencing templates amenable to DNA sequencing on recently developed platforms. Although our demonstration of the method does not utilize these DNA sequencing platforms directly, our results indicate that the capture of long DNA fragments produce superior coverage of the targeted region.

Aims/hypothesis One of the most strongly associated type 2 diabetes loci reported to date resides within the TCF7L2 gene. Previous studies point to the T allele of rs7903146 in intron 3 as the causal variant at this locus. We aimed to identify the actual gene(s) under the influence of this variant.Methods Using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein-9 nuclease, we generated a 1.4 kb deletion of the genomic region harbouring rs7903146 in the HCT116 cell line, followed by global gene expression analysis. We then carried out a combination of circularised chromosome conformation capture (4C) and Capture C in cell lines, HCT116 and NCM460 in order to ascertain which promoters of these perturbed genes made consistent physical contact with this genomic region.Results We observed 99 genes with significant differential expression (false discovery rate [FDR] cut-off: 10%) and an effect size of at least twofold. The subsequent promoter contact analyses revealed just one gene, ACSL5, which resides in the same topologically associating domain as TCF7L2. The generation of additional, smaller deletions (66 bp and 104 bp) comprising rs7903146 showed consistently reduced ACSL5 mRNA levels across all three deletions of up to 30-fold, with commensurate loss of acyl-CoA synthetase long-chain family member 5 (ACSL5) protein. Notably, the deletion of this single-nucleotide polymorphism region abolished significantly detectable chromatin contacts with the ACSL5 promoter. We went on to confirm that contacts between rs7903146 and the ACSL5 promoter regions were conserved in human colon tissue. ACSL5 encodes ACSL5, an enzyme with known roles in fatty acid metabolism.Conclusions/interpretation This 'variant to gene mapping' effort implicates the genomic location harbouring rs7903146 as a regulatory region for ACSL5.