Search

Transcriptional deregulation is a central event in the development of acute myeloid leukemia (AML). To identify potential disturbances in gene regulation, we conducted an unbiased screen of allele-specific expression (ASE) in 209 AML cases. The gene encoding GATA binding protein 2 (GATA2) displayed ASE more often than any other myeloid- or cancer-related gene. GATA2 ASE was strongly associated with CEBPA double mutations (DMs), with 95% of cases presenting GATA2 ASE. In CEBPA DM AML with GATA2 mutations, the mutated allele was preferentially expressed. We found that GATA2 ASE was a somatic event lost in complete remission, supporting the notion that it plays a role in CEBPA DM AML. Acquisition of GATA2 ASE involved silencing of 1 allele via promoter methylation and concurrent overactivation of the other allele, thereby preserving expression levels. Notably, promoter methylation was also lost in remission along with GATA2 ASE. In summary, we propose that GATA2 ASE is acquired by epigenetic mechanisms and is a prerequisite for the development of AML with CEBPA DMs. This finding constitutes a novel example of an epigenetic hit cooperating with a genetic hit in the pathogenesis of AML.

Transcriptional deregulation is a central event in the development of acute myeloid leukemia (AML). To identify potential disturbances in gene regulation, we conducted an unbiased screen of allele-specific expression (ASE) in 209 AML cases. The gene encoding GATA binding protein 2 (GATA2) displayed ASE more often than any other myeloid- or cancer-related gene. GATA2ASE was strongly associated with CEBPAdouble mutations (DMs), with 95% of cases presenting GATA2ASE. In CEBPADM AML with GATA2mutations, the mutated allele was preferentially expressed. We found that GATA2ASE was a somatic event lost in complete remission, supporting the notion that it plays a role in CEBPADM AML. Acquisition of GATA2ASE involved silencing of 1 allele via promoter methylation and concurrent overactivation of the other allele, thereby preserving expression levels. Notably, promoter methylation was also lost in remission along with GATA2ASE. In summary, we propose that GATA2ASE is acquired by epigenetic mechanisms and is a prerequisite for the development of AML with CEBPADMs. This finding constitutes a novel example of an epigenetic hit cooperating with a genetic hit in the pathogenesis of AML.

Ischemic stroke (IS), caused by obstruction of cerebral blood flow, is one of the leading causes of death. While neurologists agree on delineation of IS into three subtypes (cardioembolic stroke (CES), large artery stroke (LAS), and small vessel stroke (SVS)), several subtyping systems exist. The most commonly used systems are TOAST (Trial of Org 10172 in Acute Stroke Treatment) and CCS (Causative Classification System for Stroke), but agreement is only moderate. We have compared two approaches to combining the existing subtyping systems for a phenotype suited for a genome-wide association study (GWAS). We used the NINDS Stroke Genetics Network dataset (SiGN, 11,477 cases with CCS and TOAST subtypes and 28,026 controls). We defined two new phenotypes: the intersect, for which an individual must be assigned the same subtype by CCS and TOAST; and the union, for which an individual must be assigned a subtype by either CCS or TOAST. The union yields the largest sample size while the intersect yields a phenotype with less potential misclassification. We performed GWAS for all subtypes, using the original subtyping systems, the intersect, and the union as phenotypes. In each subtype, heritability was higher for the intersect compared with the other phenotypes. We observed stronger effects at known IS variants with the intersect compared with the other phenotypes. With the intersect, we identify rs10029218:G>A as an associated variant with SVS. We conclude that this approach increases the likelihood to detect genetic associations in ischemic stroke.

We present the experimental protocol and data analysis toolbox for multi-contact 4C (MC-4C), a new proximity ligation method tailored to study the higher-order chromatin contact patterns of selected genomic sites. Conventional chromatin conformation capture (3C) methods fragment proximity ligation products for efficient analysis of pairwise DNA contacts. By contrast, MC-4C is designed to preserve and collect large concatemers of proximity ligated fragments for long-molecule sequencing on an Oxford Nanopore or Pacific Biosciences platform. Each concatemer of proximity ligation products represents a snapshot topology of a different individual allele, revealing its multi-way chromatin interactions. By inverse PCR with primers specific for a fragment of interest (the viewpoint) and DNA size selection, sequencing is selectively targeted to thousands of different complex interactions containing this viewpoint. A tailored statistical analysis toolbox is able to generate background models and three-way interaction profiles from the same dataset. These profiles can be used to distinguish whether contacts between more than two regulatory sequences are mutually exclusive or, conversely, simultaneously occurring at chromatin hubs. The entire procedure can be completed in 2 w, and requires standard molecular biology and data analysis skills and equipment, plus access to a third-generation sequencing platform.

Background: Nanopore-based DNA sequencing relies on basecalling the electric current signal. Basecalling requires neural networks to achieve competitive accuracies. To improve sequencing accuracy further, new models are continuously proposed with new architectures. However, benchmarking is currently not standardized, and evaluation metrics and datasets used are defined on a per publication basis, impeding progress in the field. This makes it impossible to distinguish data from model driven improvements. Results: To standardize the process of benchmarking, we unified existing benchmarking datasets and defined a rigorous set of evaluation metrics. We benchmarked the latest seven basecaller models by recreating and analyzing their neural network architectures. Our results show that overall Bonito’s architecture is the best for basecalling. We find, however, that species bias in training can have a large impact on performance. Our comprehensive evaluation of 90 novel architectures demonstrates that different models excel at reducing different types of errors and using recurrent neural networks (long short-term memory) and a conditional random field decoder are the main drivers of high performing models. Conclusions: We believe that our work can facilitate the benchmarking of new basecaller tools and that the community can further expand on this work.

Introduction:An imbalance in gene expression of the different alleles of a gene (allele-specific expression or ASE) is an intriguing genome-wide and dynamic phenomenon that could have an important role in tumorigenesis. Little is known about the causes and consequences of the differential cis-and trans-regulation on resultant ASE, aside from the impact of copy number alterations (CNAs) or imprinting. Acute myeloid leukemia (AML) is an aggressive clonal hematopoietic stem cell malignancy, with a low mutational burden in comparison with other cancers, where the impact of ASE may constitute a novel mechanism of pathogenesis.

Nanopore sequencing is a rapidly maturing technology delivering long reads in real time on a portable instrument at low cost. Not surprisingly, the community has rapidly taken up this new way of sequencing and has used it successfully for a variety of research applications. A major limitation of nanopore sequencing is its high error rate, which despite recent improvements to the nanopore chemistry and computational tools still ranges between 5% and 15%. Here, we review computational approaches determining the nanopore sequencing error rate. Furthermore, we outline strategies for translation of raw sequencing data into base calls for detection of base modifications and for obtaining consensus sequences.

Chromatin folding contributes to the regulation of genomic processes such as gene activity. Existing conformation capture methods characterize genome topology through analysis of pairwise chromatin contacts in populations of cells but cannot discern whether individual interactions occur simultaneously or competitively. Here we present multi-contact 4C (MC-4C), which applies Nanopore sequencing to study multi-way DNA conformations of individual alleles. MC-4C distinguishes cooperative from random and competing interactions and identifies previously missed structures in subpopulations of cells. We show that individual elements of the β-globin superenhancer can aggregate into an enhancer hub that can simultaneously accommodate two genes. Neighboring chromatin domain loops can form rosette-like structures through collision of their CTCF-bound anchors, as seen most prominently in cells lacking the cohesin-unloading factor WAPL. Here, massive collision of CTCF-anchored chromatin loops is believed to reflect ‘cohesin traffic jams’. Single-allele topology studies thus help us understand the mechanisms underlying genome folding and functioning.

Alterations of genes encoding transcriptional regulators of lymphoid development are a hallmark of B-progenitor acute lymphoblastic leukemia (B-ALL) and most commonly involve PAX5, encoding the DNA-binding transcription factor paired-box 5. The majority of PAX5alterations in ALL are heterozygous, and key PAX5 target genes are expressed in leukemic cells, suggesting that PAX5 may be a haploinsufficient tumor suppressor. To examine the role of PAX5alterations in leukemogenesis, we performed mutagenesis screens of mice heterozygous for a loss-of-function Pax5allele. Both chemical and retroviral mutagenesis resulted in a significantly increased penetrance and reduced latency of leukemia, with a shift to B-lymphoid lineage. Genomic profiling identified a high frequency of secondary genomic mutations, deletions, and retroviral insertions targeting B-lymphoid development, including Pax5, and additional genes and pathways mutated in ALL, including tumor suppressors, Ras, and Janus kinase-signal transducer and activator of transcription signaling. These results show that in contrast to simple Pax5haploinsufficiency, multiple sequential alterations targeting lymphoid development are central to leukemogenesis and contribute to the arrest in lymphoid maturation characteristic of ALL. This cross-species analysis also validates the importance of concomitant alterations of multiple cellular growth, signaling, and tumor suppression pathways in the pathogenesis of B-ALL.

Alterations of genes encoding transcriptional regulators of lymphoid development are a hallmark of B-progenitor acute lymphoblastic leukemia (B-ALL) and most commonly involve PAX5, encoding the DNA-binding transcription factor paired-box 5. The majority of PAX5 alterations in ALL are heterozygous, and key PAX5 target genes are expressed in leukemic cells, suggesting that PAX5 may be a haploinsufficient tumor suppressor. To examine the role of PAX5 alterations in leukemogenesis, we performed mutagenesis screens of mice heterozygous for a loss-of-function Pax5 allele. Both chemical and retroviral mutagenesis resulted in a significantly increased penetrance and reduced latency of leukemia, with a shift to B-lymphoid lineage. Genomic profiling identified a high frequency of secondary genomic mutations, deletions, and retroviral insertions targeting B-lymphoid development, including Pax5, and additional genes and pathways mutated in ALL, including tumor suppressors, Ras, and Janus kinase-signal transducer and activator of transcription signaling. These results show that in contrast to simple Pax5 haploinsufficiency, multiple sequential alterations targeting lymphoid development are central to leukemogenesis and contribute to the arrest in lymphoid maturation characteristic of ALL. This cross-species analysis also validates the importance of concomitant alterations of multiple cellular growth, signaling, and tumor suppression pathways in the pathogenesis of B-ALL.

Melanoma is the most deadly form of skin cancer. Expression of oncogenic BRAFor NRAS, which are frequently mutated in human melanomas, promote the formation of nevi but are not sufficient for tumorigenesis. Even with germline mutated p53, these engineered melanomas present with variable onset and pathology, implicating additional somatic mutations in a multi-hit tumorigenic process. To decipher the genetics of these melanomas, we sequence the protein coding exons of 53 primary melanomas generated from several BRAFV600Eor NRASQ61Kdriven transgenic zebrafish lines. We find that engineered zebrafish melanomas show an overall low mutation burden, which has a strong, inverse association with the number of initiating germline drivers. Although tumors reveal distinct mutation spectrums, they show mostly C > T transitions without UV light exposure, and enrichment of mutations in melanogenesis, p53and MAPK signaling. Importantly, a recurrent amplification occurring with pre-configured drivers BRAFV600Eand p53-/-suggests a novel path of BRAFcooperativity through the protein kinase A pathway. This is the first analysis of a melanoma mutational landscape in the absence of UV light, where tumors manifest with remarkably low mutation burden and high heterogeneity. Genotype specific amplification of protein kinase A in cooperation with BRAFand p53mutation suggests the involvement of melanogenesis in these tumors. This work is important for defining the spectrum of events in BRAFor NRASdriven melanoma in the absence of UV light, and for informed exploitation of models such as transgenic zebrafish to better understand mechanisms leading to human melanoma formation.