Sustained ERK signaling promotes G2 cell cycle exit and primes cells for whole-genome duplication.

Whole-genome duplication (WGD) is a frequent event in cancer evolution that fuels chromosomal instability. WGD can result from mitotic errors or endoreduplication, yet the molecular mechanisms that drive WGD remain unclear. Here, we use live single-cell analysis to characterize cell-cycle dynamics upon aberrant Ras-ERK signaling. We find that sustained ERK signaling in human cells leads to reactivation of the APC/C in G2, resulting in tetraploid G0-like cells that are primed for WGD. This process is independent of DNA damage or p53 but dependent on p21. Transcriptomics analysis and live-cell imaging showed that constitutive ERK activity promotes p21 expression, which is necessary and sufficient to inhibit CDK activity and which prematurely activates the anaphase-promoting complex (APC/C). Finally, either loss of p53 or reduced ERK signaling allowed for endoreduplication, completing a WGD event. Thus, sustained ERK signaling-induced G2 cell cycle exit represents an alternative path to WGD.

Uncalled4 improves nanopore DNA and RNA modification detection via fast and accurate signal alignment.

Nanopore signal analysis enables detection of nucleotide modifications from native DNA and RNA sequencing, providing both accurate genetic/transcriptomic and epigenetic information without additional library preparation. Presently, only a limited set of modifications can be directly basecalled (e.g. 5-methylcytosine), while most others require exploratory methods that often begin with alignment of nanopore signal to a nucleotide reference. We present Uncalled4, a toolkit for nanopore signal alignment, analysis, and visualization. Uncalled4 features an efficient banded signal alignment algorithm, BAM signal alignment file format, statistics for comparing signal alignment methods, and a reproducible de novo training method for k-mer-based pore models, revealing potential errors in ONT's state-of-the-art DNA model. We apply Uncalled4 to RNA 6-methyladenine (m6A) detection in seven human cell lines, identifying 26% more modifications than Nanopolish using m6Anet, including in several genes where m6A has known implications in cancer. Uncalled4 is available open-source at github.com/skovaka/uncalled4.

IVT generation of guideRNAs for Cas9-enrichment Nanopore Sequencing.

Generating high-coverage sequencing coverage at select genomic loci has extensive applications in both research science and genetic medicine. Long-read sequencing technologies (e.g. nanopore sequencing) have expanded our ability to generate sequencing data in regions (e.g. repetitive elements) that are difficult to interrogate with short-read sequencing methods. In work presented here, we expand on our previous work using CRISPR/Cas9 for targeted nanopore sequencing by using in vitro transcribed guideRNAs, with 1100 guideRNAs in a single experiment. This approach decreases the cost per guideRNA, increases the number of guideRNAs that can be multiplexed in a single experiment, and provides a way to rapidly screen numerous guideRNAs for cutting efficiency. We apply this strategy in multiple patient-derived pancreatic cancer cell lines, demonstrating its ability to unveil structural variation in "deletion hotspots" around the tumor suppressor genes p16 (CDKN2A), and SMAD4.

Temozolomide-induced guanine mutations create exploitable vulnerabilities of guanine-rich DNA and RNA regions in drug-resistant gliomas.

Temozolomide (TMZ) is a chemotherapeutic agent that has been the first-line standard of care for the aggressive brain cancer glioblastoma (GBM) since 2005. Although initially beneficial, TMZ resistance is universal and second-line interventions are an unmet clinical need. Here, we took advantage of the known mechanism of action of TMZ to target guanines (G) and investigated G-rich G-quadruplex (G4) and splice site changes that occur upon TMZ resistance. We report that TMZ-resistant GBM has guanine mutations that disrupt the G-rich DNA G4s and splice sites that lead to deregulated alternative splicing. These alterations create vulnerabilities, which are selectively targeted by either the G4-stabilizing drug TMPyP4 or a novel splicing kinase inhibitor of cdc2-like kinase. Last, we show that the G4 and RNA binding protein EWSR1 aggregates in the cytoplasm in TMZ-resistant GBM cells and patient samples. Together, our findings provide insight into targetable vulnerabilities of TMZ-resistant GBM and present cytoplasmic EWSR1 as a putative biomarker.

Segmental duplications and their variation in a complete human genome.

Despite their importance in disease and evolution, highly identical segmental duplications (SDs) are among the last regions of the human reference genome (GRCh38) to be fully sequenced. Using a complete telomere-to-telomere human genome (T2T-CHM13), we present a comprehensive view of human SD organization. SDs account for nearly one-third of the additional sequence, increasing the genome-wide estimate from 5.4 to 7.0% [218 million base pairs (Mbp)]. An analysis of 268 human genomes shows that 91% of the previously unresolved T2T-CHM13 SD sequence (68.3 Mbp) better represents human copy number variation. Comparing long-read assemblies from human (n = 12) and nonhuman primate (n = 5) genomes, we systematically reconstruct the evolution and structural haplotype diversity of biomedically relevant and duplicated genes. This analysis reveals patterns of structural heterozygosity and evolutionary differences in SD organization between humans and other primates.

Epigenetically regulated digital signaling defines epithelial innate immunity at the tissue level.

To prevent damage to the host or its commensal microbiota, epithelial tissues must match the intensity of the immune response to the severity of a biological threat. Toll-like receptors allow epithelial cells to identify microbe associated molecular patterns. However, the mechanisms that mitigate biological noise in single cells to ensure quantitatively appropriate responses remain unclear. Here we address this question using single cell and single molecule approaches in mammary epithelial cells and primary organoids. We find that epithelial tissues respond to bacterial microbe associated molecular patterns by activating a subset of cells in an all-or-nothing (i.e. digital) manner. The maximum fraction of responsive cells is regulated by a bimodal epigenetic switch that licenses the TLR2 promoter for transcription across multiple generations. This mechanism confers a flexible memory of inflammatory events as well as unique spatio-temporal control of epithelial tissue-level immune responses. We propose that epigenetic licensing in individual cells allows for long-term, quantitative fine-tuning of population-level responses.

De novo genome assembly of the tobacco hornworm moth (Manduca sexta).

The tobacco hornworm, Manduca sexta, is a lepidopteran insect that is used extensively as a model system for studying insect biology, development, neuroscience, and immunity. However, current studies rely on the highly fragmented reference genome Msex_1.0, which was created using now-outdated technologies and is hindered by a variety of deficiencies and inaccuracies. We present a new reference genome for M. sexta, JHU_Msex_v1.0, applying a combination of modern technologies in a de novo assembly to increase continuity, accuracy, and completeness. The assembly is 470 Mb and is ∼20× more continuous than the original assembly, with scaffold N50 > 14 Mb. We annotated the assembly by lifting over existing annotations and supplementing with additional supporting RNA-based data for a total of 25,256 genes. The new reference assembly is accessible in annotated form for public use. We demonstrate that improved continuity of the M. sexta genome improves resequencing studies and benefits future research on M. sexta as a model organism.

Targeted nanopore sequencing by real-time mapping of raw electrical signal with UNCALLED.

Conventional targeted sequencing methods eliminate many of the benefits of nanopore sequencing, such as the ability to accurately detect structural variants or epigenetic modifications. The ReadUntil method allows nanopore devices to selectively eject reads from pores in real time, which could enable purely computational targeted sequencing. However, this requires rapid identification of on-target reads while most mapping methods require computationally intensive basecalling. We present UNCALLED ( https://github.com/skovaka/UNCALLED ), an open source mapper that rapidly matches streaming of nanopore current signals to a reference sequence. UNCALLED probabilistically considers k-mers that could be represented by the signal and then prunes the candidates based on the reference encoded within a Ferragina-Manzini index. We used UNCALLED to deplete sequencing of known bacterial genomes within a metagenomics community, enriching the remaining species 4.46-fold. UNCALLED also enriched 148 human genes associated with hereditary cancers to 29.6× coverage using one MinION flowcell, enabling accurate detection of single-nucleotide polymorphisms, insertions and deletions, structural variants and methylation in these genes.

Simultaneous profiling of chromatin accessibility and methylation on human cell lines with nanopore sequencing.

Probing epigenetic features on DNA has tremendous potential to advance our understanding of the phased epigenome. In this study, we use nanopore sequencing to evaluate CpG methylation and chromatin accessibility simultaneously on long strands of DNA by applying GpC methyltransferase to exogenously label open chromatin. We performed nanopore sequencing of nucleosome occupancy and methylome (nanoNOMe) on four human cell lines (GM12878, MCF-10A, MCF-7 and MDA-MB-231). The single-molecule resolution allows footprinting of protein and nucleosome binding, and determination of the combinatorial promoter epigenetic signature on individual molecules. Long-read sequencing makes it possible to robustly assign reads to haplotypes, allowing us to generate a fully phased human epigenome, consisting of chromosome-level allele-specific profiles of CpG methylation and chromatin accessibility. We further apply this to a breast cancer model to evaluate differential methylation and accessibility between cancerous and noncancerous cells.

Comprehensive analysis of structural variants in breast cancer genomes using single-molecule sequencing.

Improved identification of structural variants (SVs) in cancer can lead to more targeted and effective treatment options as well as advance our basic understanding of the disease and its progression. We performed whole-genome sequencing of the SKBR3 breast cancer cell line and patient-derived tumor and normal organoids from two breast cancer patients using Illumina/10x Genomics, Pacific Biosciences (PacBio), and Oxford Nanopore Technologies (ONT) sequencing. We then inferred SVs and large-scale allele-specific copy number variants (CNVs) using an ensemble of methods. Our findings show that long-read sequencing allows for substantially more accurate and sensitive SV detection, with between 90% and 95% of variants supported by each long-read technology also supported by the other. We also report high accuracy for long reads even at relatively low coverage (25×-30×). Furthermore, we integrated SV and CNV data into a unifying karyotype-graph structure to present a more accurate representation of the mutated cancer genomes. We find hundreds of variants within known cancer-related genes detectable only through long-read sequencing. These findings highlight the need for long-read sequencing of cancer genomes for the precise analysis of their genetic instability.

Telomere-to-telomere assembly of a complete human X chromosome.

After two decades of improvements, the current human reference genome (GRCh38) is the most accurate and complete vertebrate genome ever produced. However, no single chromosome has been finished end to end, and hundreds of unresolved gaps persist1,2. Here we present a human genome assembly that surpasses the continuity of GRCh382, along with a gapless, telomere-to-telomere assembly of a human chromosome. This was enabled by high-coverage, ultra-long-read nanopore sequencing of the complete hydatidiform mole CHM13 genome, combined with complementary technologies for quality improvement and validation. Focusing our efforts on the human X chromosome3, we reconstructed the centromeric satellite DNA array (approximately 3.1 Mb) and closed the 29 remaining gaps in the current reference, including new sequences from the human pseudoautosomal regions and from cancer-testis ampliconic gene families (CT-X and GAGE). These sequences will be integrated into future human reference genome releases. In addition, the complete chromosome X, combined with the ultra-long nanopore data, allowed us to map methylation patterns across complex tandem repeats and satellite arrays. Our results demonstrate that finishing the entire human genome is now within reach, and the data presented here will facilitate ongoing efforts to complete the other human chromosomes.

Characterization of Allele-Specific Regulation of Telomerase Reverse Transcriptase in Promoter Mutant Thyroid Cancer Cell Lines.

Background: Telomerase reverse transcriptase (TERT) promoter mutations play a role in carcinogenesis and are found in both tumors and cancer cell lines. TERT promoter methylation, transcription factor binding, chromatin remodeling, and alternative splicing are also known to play an integral role in TERT regulation. Methods: Using nanopore Cas9 targeted sequencing, we characterized allele-specific methylation in thyroid cancer cell lines heterozygous for the TERT promoter mutation. Furthermore, using chromatin immunoprecipitation followed by Sanger sequencing, we probed allele-specific binding of the transcription factors GABPA (GA binding protein transcription factor subunit alpha) and MYC, as well as the chromatin marks H3K4me3 and H3K27me3. Finally, using coding single nucleotide polymorphisms and the long-read sequencing, we examined complementary DNA for monoallelic expression (MAE). Results: We found the mutant TERT promoter allele to be significantly less methylated than wild type, while more methylated in the gene body in heterozygous TERT mutant cell lines. We demonstrated that the transcriptional activators GABPA and MYC bind only to the mutant TERT allele. In addition, the activating and repressive chromatin marks H3K4me3 and H3K27me3, respectively, bind mutant and wild-type alleles exclusively. Finally, in heterozygous mutant cell lines, TERT exhibits MAE from the mutant allele only. Conclusions: In summary, by employing new long-read sequencing methods, we were able to definitively demonstrate allele-specific DNA methylation, histone modifications, transcription factor binding, and the resulting monoallelic transcription in cell lines with heterozygous TERT mutations.

Targeted nanopore sequencing with Cas9-guided adapter ligation.

Despite recent improvements in sequencing methods, there remains a need for assays that provide high sequencing depth and comprehensive variant detection. Current methods1-4 are limited by the loss of native modifications, short read length, high input requirements, low yield or long protocols. In the present study, we describe nanopore Cas9-targeted sequencing (nCATS), an enrichment strategy that uses targeted cleavage of chromosomal DNA with Cas9 to ligate adapters for nanopore sequencing. We show that nCATS can simultaneously assess haplotype-resolved single-nucleotide variants, structural variations and CpG methylation. We apply nCATS to four cell lines, to a cell-line-derived xenograft, and to normal and paired tumor/normal primary human breast tissue. Median sequencing coverage was 675× using a MinION flow cell and 34× using the smaller Flongle flow cell. The nCATS sequencing requires only ~3 µg of genomic DNA and can target a large number of loci in a single reaction. The method will facilitate the use of long-read sequencing in research and in the clinic.

Transmission and clearance of potential procarcinogenic bacteria during fecal microbiota transplantation for recurrent Clostridioides difficile.

BACKGROUNDFecal microbiota transplantation (FMT) is an effective treatment for recurrent Clostridioides difficile infection (rCDI) in adults and children, but donor stool samples are currently screened for only a limited number of potential pathogens. We sought to determine whether putative procarcinogenic bacteria (enterotoxigenic Bacteroides fragilis, Fusobacterium nucleatum, and Escherichia coli harboring the colibactin toxin) could be durably transmitted from donors to patients during FMT.METHODSStool samples were collected from 11 pediatric rCDI patients and their respective FMT donors prior to FMT as well as from the patients at 2-10 weeks, 10-20 weeks, and 6 months after FMT. Bacterial virulence factors in stool DNA extracts and stool cultures were measured by quantitative PCR: Bacteroides fragilis toxin (bft), Fusobacterium adhesin A (fadA), and Escherichia coli colibactin (clbB).RESULTSFour of 11 patients demonstrated sustained acquisition of a procarcinogenic bacteria. Whole genome sequencing was performed on colony isolates from one of these donor/recipient pairs and demonstrated that clbB+ E. coli strains present in the recipient after FMT were identical to a strain present in the donor, confirming strain transmission. Conversely, 2 patients exhibited clearance of procarcinogenic bacteria following FMT from a negative donor.CONCLUSIONBoth durable transmission and clearance of procarcinogenic bacteria occurred following FMT, suggesting that additional studies on appropriate screening measures for FMT donors and the long-term consequences and/or benefits of FMT are warranted.FUNDINGCrohn's & Colitis Foundation, the Bloomberg~Kimmel Institute for Cancer Immunotherapy at Johns Hopkins University School of Medicine, the National Cancer Institute, and the Canadian Institutes of Health Research.

Epigenetic Changes Induced by Bacteroides fragilis Toxin.

Enterotoxigenic Bacteroides fragilis (ETBF) is a Gram-negative, obligate anaerobe member of the gut microbial community in up to 40% of healthy individuals. This bacterium is found more frequently in people with colorectal cancer (CRC) and causes tumor formation in the distal colon of multiple intestinal neoplasia (Apcmin/+ ) mice; tumor formation is dependent on ETBF-secreted Bacteroides fragilis toxin (BFT). Because of the extensive data connecting alterations in the epigenome with tumor formation, initial experiments attempting to connect BFT-induced tumor formation with methylation in colon epithelial cells (CECs) have been performed, but the effect of BFT on other epigenetic processes, such as chromatin structure, remains unexplored. Here, the changes in gene expression (transcriptome sequencing [RNA-seq]) and chromatin accessibility (assay for transposase-accessible chromatin using sequencing) induced by treatment of HT29/C1 cells with BFT for 24 and 48 h were examined. Our data show that several genes are differentially expressed after BFT treatment and that these changes relate to the interaction between bacteria and CECs. Further, sites of increased chromatin accessibility are associated with the location of enhancers in CECs and the binding sites of transcription factors in the AP-1/ATF family; they are also enriched for common differentially methylated regions (DMRs) in CRC. These data provide insight into the mechanisms by which BFT induces tumor formation and lay the groundwork for future in vivo studies to explore the impact of BFT on nuclear structure and function.

Characterization of human telomerase reverse transcriptase promoter methylation and transcription factor binding in differentiated thyroid cancer cell lines.

Telomerase reverse transcriptase (TERT) activation plays an important role in cancer development by enabling the immortalization of cells. TERT regulation is multifaceted, and its promoter methylation has been implicated in controlling expression through alteration in transcription factor binding. We have characterized TERT promoter methylation, transcription factor binding, and TERT expression levels in five differentiated thyroid cancer (DTC) cell lines and six normal thyroid tissue samples by targeted bisulfite sequencing, ChIP-qPCR, and qRT-PCR. DTC cell lines express varying levels of TERT and exhibit TERT promoter methylation patterns similar to patterns seen in other telomerase positive cancer cell lines. The minimal promoter immediately surrounding the transcription start site is hypomethylated, while further upstream portions show dense methylation. In contrast, the TERT promoter in normal thyroid tissue is largely unmethylated throughout and expresses TERT minimally. Transcription factor binding is also affected by TERT mutation status. The E-twenty-six (ETS) factor GABPA exhibits TERT binding in the TERT mutant DTC cells only, and allele-specific methylation patterns at the minimal promoter were observed as well, which may indicate allele-specific factor recruitment at the minimal promoter. Furthermore, we identified binding sites for activators MYC and GSC in the hypermethylated upstream region, pointing to its possible importance in TERT regulation. Overall, TERT expression and telomerase activity depend on the interplay of multiple regulatory mechanisms including TERT promoter methylation, mutation status, and recruitment of transcription factors. This work explores of the interplay between these regulatory mechanisms and offers insight into cellular control of active telomerase in human cancer.

MiR-203a is differentially expressed during branching morphogenesis and EMT in breast progenitor cells and is a repressor of peroxidasin.

MicroRNAs regulate developmental events such as branching morphogenesis, epithelial to mesenchymal transition (EMT) and its reverse process mesenchymal to epithelial transition (MET). In this study, we performed small RNA sequencing of a breast epithelial progenitor cell line (D492), and its mesenchymal derivative (D492M) cultured in three-dimensional microenvironment. Among the most downregulated miRNAs in D492M was miR-203a, a miRNA that plays an important role in epithelial differentiation. Increased expression of miR-203a was seen in D492, concomitant with increased complexity of branching. When miR-203a was overexpressed in D492M, a partial reversion towards epithelial phenotype was seen. Gene expression analysis of D492M and D492MmiR-203a revealed peroxidasin, a collagen IV cross-linker, as the most significantly downregulated gene in D492MmiR-203a. Collectively, we demonstrate that miR-203a expression temporally correlates with branching morphogenesis and is suppressed in D492M. Overexpression of miR-203a in D492M induces a partial MET and reduces the expression of peroxidasin. Furthermore, we demonstrate that miR-203a is a novel repressor of peroxidasin. MiR-203-peroxidasin axis may be an important regulator in branching morphogenesis, EMT/MET and basement membrane remodeling.

A Loss of Epigenetic Control Can Promote Cell Death through Reversing the Balance of Pathways in a Signaling Network.

Epigenetic control of regulatory networks is only partially understood. Expression of Insulin-like growth factor-II (IGF2) is controlled by genomic imprinting, mediated by silencing of the maternal allele. Loss of imprinting of IGF2 (LOI) is linked to intestinal and colorectal cancers, causally in murine models and epidemiologically in humans. However, the molecular underpinnings of the LOI phenotype are not clear. Surprisingly, in LOI cells, we find a reversal of the relative activities of two canonical signaling pathways triggered by IGF2, causing further rebalancing between pro- and anti-apoptotic signaling. A predictive mathematical model shows that this network rebalancing quantitatively accounts for the effect of receptor tyrosine kinase inhibition in both WT and LOI cells. This mechanism also quantitatively explains both the stable LOI phenotype and the therapeutic window for selective killing of LOI cells, and thus prevention of epigenetically controlled cancers. These findings suggest a framework for understanding epigenetically modified cell signaling.

Detecting DNA cytosine methylation using nanopore sequencing.

In nanopore sequencing devices, electrolytic current signals are sensitive to base modifications, such as 5-methylcytosine (5-mC). Here we quantified the strength of this effect for the Oxford Nanopore Technologies MinION sequencer. By using synthetically methylated DNA, we were able to train a hidden Markov model to distinguish 5-mC from unmethylated cytosine. We applied our method to sequence the methylome of human DNA, without requiring special steps for library preparation.

Nanopore sequencing detects structural variants in cancer.

Despite advances in sequencing, structural variants (SVs) remain difficult to reliably detect due to the short read length (<300 bp) of 2nd generation sequencing. Not only do the reads (or paired-end reads) need to straddle a breakpoint, but repetitive elements often lead to ambiguities in the alignment of short reads. We propose to use the long-reads (up to 20 kb) possible with 3rd generation sequencing, specifically nanopore sequencing on the MinION. Nanopore sequencing relies on a similar concept to a Coulter counter, reading the DNA sequence from the change in electrical current resulting from a DNA strand being forced through a nanometer-sized pore embedded in a membrane. Though nanopore sequencing currently has a relatively high mismatch rate that precludes base substitution and small frameshift mutation detection, its accuracy is sufficient for SV detection because of its long reads. In fact, long reads in some cases may improve SV detection efficiency. We have tested nanopore sequencing to detect a series of well-characterized SVs, including large deletions, inversions, and translocations that inactivate the CDKN2A/p16 and SMAD4/DPC4 tumor suppressor genes in pancreatic cancer. Using PCR amplicon mixes, we have demonstrated that nanopore sequencing can detect large deletions, translocations and inversions at dilutions as low as 1:100, with as few as 500 reads per sample. Given the speed, small footprint, and low capital cost, nanopore sequencing could become the ideal tool for the low-level detection of cancer-associated SVs needed for molecular relapse, early detection, or therapeutic monitoring.