Neotelomeres and telomere-spanning chromosomal arm fusions in cancer genomes revealed by long-read sequencing.

Alterations in the structure and location of telomeres are pivotal in cancer genome evolution. Here, we applied both long-read and short-read genome sequencing to assess telomere repeat-containing structures in cancers and cancer cell lines. Using long-read genome sequences that span telomeric repeats, we defined four types of telomere repeat variations in cancer cells: neotelomeres where telomere addition heals chromosome breaks, chromosomal arm fusions spanning telomere repeats, fusions of neotelomeres, and peri-centromeric fusions with adjoined telomere and centromere repeats. These results provide a framework for the systematic study of telomeric repeats in cancer genomes, which could serve as a model for understanding the somatic evolution of other repetitive genomic elements.

Neotelomeres and Telomere-Spanning Chromosomal Arm Fusions in Cancer Genomes Revealed by Long-Read Sequencing.

Alterations in the structure and location of telomeres are key events in cancer genome evolution. However, previous genomic approaches, unable to span long telomeric repeat arrays, could not characterize the nature of these alterations. Here, we applied both long-read and short-read genome sequencing to assess telomere repeat-containing structures in cancers and cancer cell lines. Using long-read genome sequences that span telomeric repeat arrays, we defined four types of telomere repeat variations in cancer cells: neotelomeres where telomere addition heals chromosome breaks, chromosomal arm fusions spanning telomere repeats, fusions of neotelomeres, and peri-centromeric fusions with adjoined telomere and centromere repeats. Analysis of lung adenocarcinoma genome sequences identified somatic neotelomere and telomere-spanning fusion alterations. These results provide a framework for systematic study of telomeric repeat arrays in cancer genomes, that could serve as a model for understanding the somatic evolution of other repetitive genomic elements.

Identifying and correcting repeat-calling errors in nanopore sequencing of telomeres.

Nanopore long-read sequencing is an emerging approach for studying genomes, including long repetitive elements like telomeres. Here, we report extensive basecalling induced errors at telomere repeats across nanopore datasets, sequencing platforms, basecallers, and basecalling models. We find that telomeres in many organisms are frequently miscalled. We demonstrate that tuning of nanopore basecalling models leads to improved recovery and analysis of telomeric regions, with minimal negative impact on other genomic regions. We highlight the importance of verifying nanopore basecalls in long, repetitive, and poorly defined regions, and showcase how artefacts can be resolved by improvements in nanopore basecalling models.

Clinical and molecular characterization of virus-positive and virus-negative Merkel cell carcinoma.

BACKGROUND: Merkel cell carcinoma (MCC) is a highly aggressive neuroendocrine carcinoma of the skin caused by either the integration of Merkel cell polyomavirus (MCPyV) and expression of viral T antigens or by ultraviolet-induced damage to the tumor genome from excessive sunlight exposure. An increasing number of deep sequencing studies of MCC have identified significant differences between the number and types of point mutations, copy number alterations, and structural variants between virus-positive and virus-negative tumors. However, it has been challenging to reliably distinguish between virus positive and UV damaged MCC. METHODS: In this study, we assembled a cohort of 71 MCC patients and performed deep sequencing with OncoPanel, a clinically implemented, next-generation sequencing assay targeting over 400 cancer-associated genes. To improve the accuracy and sensitivity for virus detection compared to traditional PCR and IHC methods, we developed a hybrid capture baitset against the entire MCPyV genome and software to detect integration sites and structure. RESULTS: Sequencing from this approach revealed distinct integration junctions in the tumor genome and generated assemblies that strongly support a model of microhomology-initiated hybrid, virus-host, circular DNA intermediate that promotes focal amplification of host and viral DNA. Using the clear delineation between virus-positive and virus-negative tumors from this method, we identified recurrent somatic alterations common across MCC and alterations specific to each class of tumor, associated with differences in overall survival. Finally, comparing the molecular and clinical data from these patients revealed a surprising association of immunosuppression with virus-negative MCC and significantly shortened overall survival. CONCLUSIONS: These results demonstrate the value of high-confidence virus detection for identifying molecular mechanisms of UV and viral oncogenesis in MCC. Furthermore, integrating these data with clinical data revealed features that could impact patient outcome and improve our understanding of MCC risk factors.

ViroPanel: Hybrid Capture and Massively Parallel Sequencing for Simultaneous Detection and Profiling of Oncogenic Virus Infection and Tumor Genome.

Precision cancer medicine aims to classify tumors by site, histology, and molecular testing to determine an individualized profile of cancer alterations. Viruses are a major contributor to oncogenesis, causing 12% to 20% of all human cancers. Several viruses are causal of specific types of cancer, promoting dysregulation of signaling pathways and resulting in carcinogenesis. In addition, integration of viral DNA into the host (human) genome is a hallmark of some viral species. Tests for the presence of viral infection used in the clinical setting most often use quantitative PCR or immunohistochemical staining. Both approaches have limitations and need to be interpreted/scored appropriately. In some cases, results are not binary (virus present/absent), and it is unclear what to do with a weakly or partially positive result. In addition, viral testing of cancers is performed separately from tests to detect human genomic alterations and can thus be time-consuming and use limited valuable specimen. We present a hybrid-capture and massively parallel sequencing approach to detect viral infection that is integrated with targeted genomic analysis to provide a more complete tumor profile from a single sample.

Unique, dual-indexed sequencing adapters with UMIs effectively eliminate index cross-talk and significantly improve sensitivity of massively parallel sequencing.

BACKGROUND: Sample index cross-talk can result in false positive calls when massively parallel sequencing (MPS) is used for sensitive applications such as low-frequency somatic variant discovery, ancient DNA investigations, microbial detection in human samples, or circulating cell-free tumor DNA (ctDNA) variant detection. Therefore, the limit-of-detection of an MPS assay is directly related to the degree of index cross-talk. RESULTS: Cross-talk rates up to 0.29% were observed when using standard, combinatorial adapters, resulting in 110,180 (0.1% cross-talk rate) or 1,121,074 (0.29% cross-talk rate) misassigned reads per lane in non-patterned and patterned Illumina flow cells, respectively. Here, we demonstrate that using unique, dual-matched indexed adapters dramatically reduces index cross-talk to ≤1 misassigned reads per flow cell lane. While the current study was performed using dual-matched indices, using unique, dual-unrelated indices would also be an effective alternative. CONCLUSIONS: For sensitive downstream analyses, the use of combinatorial indices for multiplexed hybrid capture and sequencing is inappropriate, as it results in an unacceptable number of misassigned reads. Cross-talk can be virtually eliminated using dual-matched indexed adapters. These results suggest that use of such adapters is critical to reduce false positive rates in assays that aim to identify low allele frequency events, and strongly indicate that dual-matched adapters be implemented for all sensitive MPS applications.

Synthetic mRNA with Superior Properties that Mimics the Intracellular Fates of Natural Histone mRNA.

Since DNA and histone levels must be closely balanced for cell survival, histone expressions are highly regulated. The regulation of replication-dependent histone expression is mainly achieved at the mRNA level, as the mRNAs are rapidly removed when DNA replication is inhibited during S-phase. Histone mRNA degradation initiates with addition of multiple uridines (oligouridylation) following the 3' stem-loop (SL) catalyzed by terminal uridyltransferase (TUTase). Previous studies showed that histone mRNA degradation occurs through both 5' → 3' and 3' → 5' processes, but the relative contributions are difficult to dissect due to lack of established protocols. The translational efficiency and stability of synthetic mRNA in both cultured cells and whole animals can be improved by structural modifications at the both 5' and 3' termini. In this chapter, we present methods of utilizing modified cap dinucleotide analogs to block 5' → 3' degradation of a reporter mRNA containing canonical histone mRNA 3' SL and monitoring how oligouridylation and 3' → 5' degradation occur. Protocols are presented for synthesis of reporter mRNA containing the histone 3' SL and modified cap analogs, monitoring mRNA stability and unidirectional degradation either from 5' or 3' termini, and detection of oligo(U) tracts from degradation products by either traditional or deep sequencing.

EnD-Seq and AppEnD: sequencing 3' ends to identify nontemplated tails and degradation intermediates.

Existing methods for detecting RNA intermediates resulting from exonuclease degradation are low-throughput and laborious. In addition, mapping the 3' ends of RNA molecules to the genome after high-throughput sequencing is challenging, particularly if the 3' ends contain post-transcriptional modifications. To address these problems, we developed EnD-Seq, a high-throughput sequencing protocol that preserves the 3' end of RNA molecules, and AppEnD, a computational method for analyzing high-throughput sequencing data. Together these allow determination of the 3' ends of RNA molecules, including nontemplated additions. Applying EnD-Seq and AppEnD to histone mRNAs revealed that a significant fraction of cytoplasmic histone mRNAs end in one or two uridines, which have replaced the 1-2 nt at the 3' end of mature histone mRNA maintaining the length of the histone transcripts. Histone mRNAs in fly embryos and ovaries show the same pattern, but with different tail nucleotide compositions. We increase the sensitivity of EnD-Seq by using cDNA priming to specifically enrich low-abundance tails of known sequence composition allowing identification of degradation intermediates. In addition, we show the broad applicability of our computational approach by using AppEnD to gain insight into 3' additions from diverse types of sequencing data, including data from small capped RNA sequencing and some alternative polyadenylation protocols.

Deep sequencing shows multiple oligouridylations are required for 3' to 5' degradation of histone mRNAs on polyribosomes.

Histone mRNAs are rapidly degraded when DNA replication is inhibited during S phase with degradation initiating with oligouridylation of the stem loop at the 3' end. We developed a customized RNA sequencing strategy to identify the 3' termini of degradation intermediates of histone mRNAs. Using this strategy, we identified two types of oligouridylated degradation intermediates: RNAs ending at different sites of the 3' side of the stem loop that resulted from initial degradation by 3'hExo and intermediates near the stop codon and within the coding region. Sequencing of polyribosomal histone mRNAs revealed that degradation initiates and proceeds 3' to 5' on translating mRNA and that many intermediates are capped. Knockdown of the exosome-associated exonuclease PM/Scl-100, but not the Dis3L2 exonuclease, slows histone mRNA degradation consistent with 3' to 5' degradation by the exosome containing PM/Scl-100. Knockdown of No-go decay factors also slowed histone mRNA degradation, suggesting a role in removing ribosomes from partially degraded mRNAs.

mRNAs containing the histone 3' stem-loop are degraded primarily by decapping mediated by oligouridylation of the 3' end.

Metazoan replication-dependent histone mRNAs are only present in S-phase, due partly to changes in their stability. These mRNAs end in a unique stem-loop (SL) that is required for both translation and cell-cycle regulation. Previous studies showed that histone mRNA degradation occurs through both 5'→3' and 3'→5' processes, but the relative contributions are not known. The 3' end of histone mRNA is oligouridylated during its degradation, although it is not known whether this is an essential step. We introduced firefly luciferase reporter mRNAs containing the histone 3' UTR SL (Luc-SL) and either a normal or hDcp2-resistant cap into S-phase HeLa cells. Both mRNAs were translated, and translation initially protected the mRNAs from degradation, but there was a lag of ∼40 min with the uncleavable cap compared to ∼8 min for the normal cap before rapid decay. Knockdown of hDcp2 resulted in a similar longer lag for Luc-SL containing a normal cap, indicating that 5'→3' decay is important in this system. Inhibition of DNA replication with hydroxyurea accelerated the degradation of Luc-SL. Knockdown of terminal uridyltransferase (TUTase) 4 but not TUTase 3 slowed the decay process, but TUTase 4 knockdown had no effect on destabilization of the mRNA by hydroxyurea. Both Luc-SL and its 5' decay intermediates were oligouridylated. Preventing oligouridylation by 3'-deoxyadenosine (cordycepin) addition to the mRNA slowed degradation, in the presence or absence of hydroxyurea, suggesting oligouridylation initiates degradation. The spectrum of oligouridylated fragments suggests the 3'→5' degradation machinery stalls during initial degradation, whereupon reuridylation occurs.

Circular RNAs are abundant, conserved, and associated with ALU repeats.

Circular RNAs composed of exonic sequence have been described in a small number of genes. Thought to result from splicing errors, circular RNA species possess no known function. To delineate the universe of endogenous circular RNAs, we performed high-throughput sequencing (RNA-seq) of libraries prepared from ribosome-depleted RNA with or without digestion with the RNA exonuclease, RNase R. We identified >25,000 distinct RNA species in human fibroblasts that contained non-colinear exons (a "backsplice") and were reproducibly enriched by exonuclease degradation of linear RNA. These RNAs were validated as circular RNA (ecircRNA), rather than linear RNA, and were more stable than associated linear mRNAs in vivo. In some cases, the abundance of circular molecules exceeded that of associated linear mRNA by >10-fold. By conservative estimate, we identified ecircRNAs from 14.4% of actively transcribed genes in human fibroblasts. Application of this method to murine testis RNA identified 69 ecircRNAs in precisely orthologous locations to human circular RNAs. Of note, paralogous kinases HIPK2 and HIPK3 produce abundant ecircRNA from their second exon in both humans and mice. Though HIPK3 circular RNAs contain an AUG translation start, it and other ecircRNAs were not bound to ribosomes. Circular RNAs could be degraded by siRNAs and, therefore, may act as competing endogenous RNAs. Bioinformatic analysis revealed shared features of circularized exons, including long bordering introns that contained complementary ALU repeats. These data show that ecircRNAs are abundant, stable, conserved and nonrandom products of RNA splicing that could be involved in control of gene expression.

The F-box protein Cdc4/Fbxw7 is a novel regulator of neural crest development in Xenopus laevis.

BACKGROUND: The neural crest is a unique population of cells that arise in the vertebrate ectoderm at the neural plate border after which they migrate extensively throughout the embryo, giving rise to a wide range of derivatives. A number of proteins involved in neural crest development have dynamic expression patterns, and it is becoming clear that ubiquitin-mediated protein degradation is partly responsible for this. RESULTS: Here we demonstrate a novel role for the F-box protein Cdc4/Fbxw7 in neural crest development. Two isoforms of Xenopus laevis Cdc4 were identified, and designated xCdc4alpha and xCdc4beta. These are highly conserved with vertebrate Cdc4 orthologs, and the Xenopus proteins are functionally equivalent in terms of their ability to degrade Cyclin E, an established vertebrate Cdc4 target. Blocking xCdc4 function specifically inhibited neural crest development at an early stage, prior to expression of c-Myc, Snail2 and Snail. CONCLUSIONS: We demonstrate that Cdc4, an ubiquitin E3 ligase subunit previously identified as targeting primarily cell cycle regulators for proteolysis, has additional roles in control of formation of the neural crest. Hence, we identify Cdc4 as a protein with separable but complementary functions in control of cell proliferation and differentiation.

ElrA binding to the 3'UTR of cyclin E1 mRNA requires polyadenylation elements.

The early cell divisions of Xenopus laevis and other metazoan embryos occur in the presence of constitutively high levels of the cell cycle regulator cyclin E1. Upon completion of the 12th cell division, a time at which many maternal proteins are downregulated by deadenylation and destabilization of their encoding mRNAs, maternal cyclin E1 protein is downregulated while its mRNA is polyadenylated and stable. We report here that stable polyadenylation of cyclin E1 mRNA requires three cis-acting elements in the 3' untranslated region; the nuclear polyadenylation sequence, a contiguous cytoplasmic polyadenylation element and an upstream AU-rich element. ElrA, the Xenopus homolog of HuR and a member of the ELAV gene family binds the cyclin E1 3'UTR with high affinity. Deletion of these elements dramatically reduces the affinity of ElrA for the cyclin E1 3'UTR, abolishes polyadenylation and destabilizes the mRNA. Together, these findings provide compelling evidence that ElrA functions in polyadenylation and stabilization of cyclin E1 mRNA via binding these elements.

Antisense knockdown of cyclin E does not affect the midblastula transition in Xenopus laevis embryos.

In Xenopus laevis embryos, cyclin E protein remains constitutively high throughout the first 12 cell cycles following fertilization until the onset of the midblastula transition (MBT) (after the 12(th) cell cycle) when it undergoes a dramatic reduction. The disappearance of cyclin E at the MBT occurs independently of active cell cycle progression, zygotic transcription, protein synthesis and the nuclear to cytoplasmic ratio. This has suggested that cyclin E is part of an autonomous maternal timer that regulates the onset of the MBT. To determine how constitutively high levels of cyclin E are maintained prior to the MBT and to investigate if the reduction in cyclin E protein affects the timing of the MBT, we have knocked down endogenous cyclin E mRNA using an N,N-diethyl-ethylene-diamine modified antisense oligonucleotide targeted to its open reading frame. We report that maintenance of high levels of cyclin E protein before the MBT is due to a balance between ongoing translation and proteolytic degradation. In support of our antisense experiments, polysome analysis demonstrates that cyclin E mRNA is associated with the translated fraction prior to the MBT. Furthermore, knockdown of cyclin E was not associated with defects in the timing of developmental events. Our data suggests that cyclin E is not required for the later cell cycles of embryonic development and that the pathway effecting downregulation of cyclin E rather then cyclin E degradation itself may be part of a maternal timer that affects the onset of the MBT.