TelomereHunter - in silico estimation of telomere content and composition from cancer genomes

Multi-omic analysis of SDHB-deficient pheochromocytomas and paragangliomas identifies metastasis and treatment-related molecular profiles

Hereditary SDHB-mutant pheochromocytomas (PC) and paragangliomas (PG) are rare tumours with a high propensity to metastasize although their clinical behaviour is unpredictable. To characterize the genomic landscape of these tumours and identify metastasis biomarkers, we performed multi-omic analysis on 94 tumours from 79 patients using seven molecular methods. Sympathetic (chromaffin cell) and parasympathetic (non-chromaffin cell) PCPG had distinct molecular profiles reflecting their cell-of-origin and biochemical profile. TERT and ATRX-alterations were associated with metastatic PCPG and these tumours had an increased mutation load, and distinct transcriptional and telomeric features. Most PCPG had quiet genomes with some rare co-operative driver events observed, including EPAS1/HIF-2α mutations. Two mechanisms of acquired resistance to DNA alkylating chemotherapies were also detected - MGMT overexpression and mismatch repair-deficiency causing hypermutation. Our comprehensive multi-omic analysis of SDHB-mutant PCPG therefore identified features of metastatic disease and treatment response, expanding our understanding of these rare neuroendocrine tumours.

Cell-type-specific effects of age and sex on human cortical neurons

Altered transcriptional and epigenetic regulation of brain cell types may contribute to cognitive changes with advanced age. Using single-nucleus multi-omic DNA methylation and transcriptome sequencing (snmCT-seq) in frontal cortex from young adult and aged donors, we found widespread age- and sex-related variation in specific neuron types. The proportion of inhibitory SST- and VIP-expressing neurons was reduced in aged donors. Excitatory neurons had more profound age-related changes in their gene expression and DNA methylation than inhibitory cells. Hundreds of genes involved in synaptic activity, including EGR1, were less expressed in aged adults. Genes located in subtelomeric regions increased their expression with age and correlated with reduced telomere length. We further mapped cell-type-specific sex differences in gene expression and X-inactivation escape genes. Multi-omic single-nucleus epigenomes and transcriptomes provide new insight into the effects of age and sex on human neurons.

Characterization of telomere variant repeats using long reads enables allele-specific telomere length estimation

Telomeres are regions of repetitive DNA at the ends of linear chromosomes which protect chromosome ends from degradation. Telomere lengths have been extensively studied in the context of aging and disease, though most studies use average telomere lengths which are of limited utility. We present a method for identifying all 92 telomere alleles from long read sequencing data. Individual telomeres are identified using variant repeats proximal to telomere regions, which are unique across alleles. This high-throughput and high-resolution characterization of telomeres could be foundational to future studies investigating the roles of specific telomeres in aging and disease.

The evolution of lung adenocarcinoma precursors is associated with chromosomal instability and transition from innate to adaptive immune response/evasion

Studying lung adenocarcinoma (LUAD) early carcinogenesis is challenging, primarily due to the lack of LUAD precursors specimens. We amassed multi-omics data from 213 LUAD and LUAD precursors to identify molecular features underlying LUAD precancer evolution. We observed progressively increasing mutations, chromosomal aberrations, whole genome doubling and genomic instability from precancer to invasive LUAD, indicating aggravating chromosomal instability (CIN). Telomere shortening, a crucial genomic alteration linked to CIN, emerged at precancer stage. Moreover, later-stage lesions demonstrated increasing cancer stemness and decreasing alveolar identity, suggesting epithelial de-differentiation during early LUAD carcinogenesis. The innate immune cells progressively diminished from precancer to invasive LUAD, concomitant with a gradual recruitment of adaptive immune cells (except CD8+ and gamma-delta T cells that decreased in later stages) and upregulation of numerous immune checkpoints, suggesting LUAD precancer evolution is associated with a shift from innate to adaptive immune response and immune evasion mediated by various mechanisms.

Reprogramming Chromosome Ends by Functional Histone Acetylation

Cancers harness embryonic programs to evade aging and promote survival. Normally, sequences at chromosome ends called telomeres shorten with cell division, serving as a countdown clock to limit cell replication. Therefore, a crucial aspect of cancerous transformation is avoiding replicative aging by activation of telomere repair programs. Mouse embryonic stem cells (mESCs) activate a transient expression of the gene Zscan4, which correlates with chromatin de-condensation and telomere extension. Head and neck squamous cell carcinoma (HNSCC) cancers reactivate ZSCAN4, which in turn regulates the phenotype of cancer stem cells (CSCs). Our study reveals a new role for human ZSCAN4 in facilitating functional histone H3 acetylation at telomere chromatin. Next-generation sequencing indicates ZSCAN4 enrichment at telomere chromatin. These changes correlate with ZSCAN4-induced histone H3 acetylation and telomere elongation, while CRISPR/Cas9 knockout of ZSCAN4 leads to reduced H3 acetylation and telomere shortening. Our study elucidates the intricate involvement of ZSCAN4 and its significant contribution to telomere chromatin remodeling. These findings suggest that ZSCAN4 induction serves as a novel link between 'stemness' and telomere maintenance. Targeting ZSCAN4 may offer new therapeutic approaches to effectively limit or enhance the replicative lifespan of stem cells and cancer cells.

METTL3 drives telomere targeting of TERRA lncRNA through m6A-dependent R-loop formation: a therapeutic target for ALT-positive neuroblastoma

Telomerase-negative tumors maintain telomere length by alternative lengthening of telomeres (ALT), but the underlying mechanism behind ALT remains poorly understood. A proportion of aggressive neuroblastoma (NB), particularly relapsed tumors, are positive for ALT (ALT+), suggesting that a better dissection of the ALT mechanism could lead to novel therapeutic opportunities. TERRA, a long non-coding RNA (lncRNA) derived from telomere ends, localizes to telomeres in a R-loop-dependent manner and plays a crucial role in telomere maintenance. Here we present evidence that RNA modification at the N6 position of internal adenosine (m6A) in TERRA by the methyltransferase METTL3 is essential for telomere maintenance in ALT+ cells, and the loss of TERRA m6A/METTL3 results in telomere damage. We observed that m6A modification is abundant in R-loop enriched TERRA, and the m6A-mediated recruitment of hnRNPA2B1 to TERRA is critical for R-loop formation. Our findings suggest that m6A drives telomere targeting of TERRA via R-loops, and this m6A-mediated R-loop formation could be a widespread mechanism employed by other chromatin-interacting lncRNAs. Furthermore, treatment of ALT+ NB cells with a METTL3 inhibitor resulted in compromised telomere targeting of TERRA and accumulation of DNA damage at telomeres, indicating that METTL3 inhibition may represent a therapeutic approach for ALT+ NB.

Blastocyst telomere length predicts successful implantation after frozen-thawed embryo transfer

STUDY QUESTION

Do embryos with longer telomere length (TL) at the blastocyst stage have a higher capacity to survive after frozen-thawed embryo transfer (FET)?

SUMMARY ANSWER

Digitally estimated TL using low-pass whole genome sequencing (WGS) data from the preimplantation genetic testing for aneuploidy (PGT-A) process demonstrates that blastocyst TL is the most essential factor associated with likelihood of implantation.

WHAT IS KNOWN ALREADY

The lifetime TL is established in the early cleavage cycles following fertilization through a recombination-based lengthening mechanism and starts erosion beyond the blastocyst stage. In addition, a telomerase-mediated slow erosion of TL in human fetuses has been observed from a gestational age of 6-11 weeks. Finally, an abnormal shortening of telomeres is likely involved in embryo loss during early development.

STUDY DESIGN SIZE DURATION

Blastocyst samples were obtained from patients who underwent PGT-A and FET in an IVF center from March 2015 to May 2018. Digitally estimated mitochondrial copy number (mtCN) and TL were used to study associations with the implantation potential of each embryo.

PARTICIPANTS/MATERIALS SETTING AND METHODS

In total, 965 blastocysts from 232 cycles (164 patients) were available to investigate the biological and clinical relevance of TL. A WGS-based workflow was applied to determine the ploidy of each embryo. Data from low-pass WGS-PGT-A were used to estimate the mtCN and TL for each embryo. Single-variant and multi-variant logistic regression, decision tree, and random forest models were applied to study various factors in association with the implantation potential of each embryo.

MAIN RESULTS AND THE ROLE OF CHANCE

Of the 965 blastocysts originally available, only 216 underwent FET. While mtCN from the transferred embryos is significantly associated with the ploidy call of each embryo, mtCN has no role in impacting IVF outcomes after an embryo transfer in these women. The results indicate that mtCN is a marker of embryo aneuploidy. On the other hand, digitally estimated TL is the most prominent univariant factor and showed a significant positive association with pregnancy outcomes (P < 0.01, odds ratio 79.1). We combined several maternal and embryo parameters to study the joint effects on successful implantation. The machine learning models, namely decision tree and random forest, were trained and yielded classification accuracy of 0.82 and 0.91, respectively. Taken together, these results support the vital role of TL in governing implantation potential, perhaps through the ability to control embryo survival after transfer.

LIMITATIONS REASONS FOR CAUTION

The small sample size limits our study as only 216 blastocysts were transferred. The number was further reduced to 153 blastocysts, where pregnancy outcomes could be accurately traced. The other limitation of this study is that all data were collected from a single IVF center. The uniform and controlled operation of IVF cycles in a single center may cause selection bias.

WIDER IMPLICATIONS OF THE FINDINGS

We present novel findings to show that digitally estimated TL at the blastocyst stage is a predictor of pregnancy capacity after a FET cycle. As elective single-embryo transfer has become the mainstream direction in reproductive medicine, prioritizing embryos based on their implantation potential is crucial for clinical infertility treatment in order to reduce twin pregnancy rate and the time to pregnancy in an IVF center. The AI-powered, random forest prediction model established in this study thus provides a way to improve clinical practice and optimize the chances for people with fertility problems to achieve parenthood.

STUDY FUNDING/COMPETING INTERESTS

This study was supported by a grant from the National Science and Technology Council, Taiwan (MOST 108-2321-B-006-013 -). There were no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

Advances in the discovery and analyses of human tandem repeats

Long-read sequencing platforms provide unparalleled access to the structure and composition of all classes of tandemly repeated DNA from STRs to satellite arrays. This review summarizes our current understanding of their organization within the human genome, their importance with respect to disease, as well as the advances and challenges in understanding their genetic diversity and functional effects. Novel computational methods are being developed to visualize and associate these complex patterns of human variation with disease, expression, and epigenetic differences. We predict accurate characterization of this repeat-rich form of human variation will become increasingly relevant to both basic and clinical human genetics.

Experimental and Computational Approaches to Measure Telomere Length: Recent Advances and Future Directions

PURPOSE OF REVIEW

The length of telomeres, protective structures at the chromosome ends, is a well-established biomarker for pathological conditions including multisystemic syndromes called telomere biology disorders. Approaches to measure telomere length (TL) differ on whether they estimate average, distribution, or chromosome-specific TL, and each presents their own advantages and limitations.

RECENT FINDINGS

The development of long-read sequencing and publication of the telomere-to-telomere human genome reference has allowed for scalable and high-resolution TL estimation in pre-existing sequencing datasets but is still impractical as a dedicated TL test. As sequencing costs continue to fall and strategies for selectively enriching telomere regions prior to sequencing improve, these approaches may become a promising alternative to classic methods. Measurement methods rely on probe hybridization, qPCR or more recently, computational methods using sequencing data. Refinements of existing techniques and new approaches have been recently developed but a test that is accurate, simple, and scalable is still lacking.

ChArmTelo Enables Large-Scale Chromosome Arm-Level Telomere Analysis across Human Populations and in Cancer Patients

Telomeres are structures protecting chromosome ends. However, a scalable and cost-effective method to investigate chromosome arm-level (ChArm) telomeres (Telos) in large-scale projects is still lacking, hindering intensive investigation of high-resolution telomeres across cancers and other diseases. Here, ChArmTelo, the first computational toolbox to analyze telomeres at chromosome arm level in human and other animal species, using 10X linked-read and similar technologies, is presented. ChArmTelo currently consists of two algorithms, TeloEM and TeloKnow, for arm-level telomere length (TL) analysis. The algorithms are demonstrated by comprehensive analysis of chromosome arm-level telomere lengths (chArmTLs) in nearly 400 whole genome sequencing samples (WGS) from human populations and animals, including healthy and cancer samples. Notably, considerable performance improvement contributed by using the latest complete telomere-to-telomere reference genome (CHM13v2), compared to hg38, is shown. ChArmTelo reveals population-specific chArmTL differences and liver cancer signatures of chArmTLs and that DNA replication origin disruption may contribute to cancer by affecting TLs. Importantly, ChArmTelo can be readily applied to tens of thousands of cancer and healthy samples with published WGS data.

NCT/DKFZ MASTER handbook of interpreting whole-genome, transcriptome, and methylome data for precision oncology

Analysis of selected cancer genes has become an important tool in precision oncology but cannot fully capture the molecular features and, most importantly, vulnerabilities of individual tumors. Observational and interventional studies have shown that decision-making based on comprehensive molecular characterization adds significant clinical value. However, the complexity and heterogeneity of the resulting data are major challenges for disciplines involved in interpretation and recommendations for individualized care, and limited information exists on how to approach multilayered tumor profiles in clinical routine. We report our experience with the practical use of data from whole-genome or exome and RNA sequencing and DNA methylation profiling within the MASTER (Molecularly Aided Stratification for Tumor Eradication Research) program of the National Center for Tumor Diseases (NCT) Heidelberg and Dresden and the German Cancer Research Center (DKFZ). We cover all relevant steps of an end-to-end precision oncology workflow, from sample collection, molecular analysis, and variant prioritization to assigning treatment recommendations and discussion in the molecular tumor board. To provide insight into our approach to multidimensional tumor profiles and guidance on interpreting their biological impact and diagnostic and therapeutic implications, we present case studies from the NCT/DKFZ molecular tumor board that illustrate our daily practice. This manual is intended to be useful for physicians, biologists, and bioinformaticians involved in the clinical interpretation of genome-wide molecular information.

Distinct characteristics of two types of alternative lengthening of telomeres in mouse embryonic stem cells

Telomere length must be maintained in actively dividing cells to avoid cellular arrest or death. In the absence of telomerase activity, activation of alternative lengthening of telomeres (ALT) allows the maintenance of telomeric length and prolongs the cellular lifespan. Our previous studies have established two types of ALT survivors from mouse embryonic stem cells. The key differences between these ALT survivors are telomere-constituting sequences: non-telomeric sequences and canonical telomeric repeats, with each type of ALT survivors being referred to as type I and type II, respectively. We explored how the characteristics of the two types of ALT lines reflect their fates using multi-omics approaches. The most notable gene expression signatures of type I and type II ALT cell lines were chromatin remodelling and DNA repair, respectively. Compared with type II cells, type I ALT cells accumulated more mutations and demonstrated persistent telomere instability. These findings indicate that cells of the same origin have separate routes for survival, thus providing insights into the plasticity of crisis-suffering cells and cancers.

A pan-sarcoma landscape of telomeric content shows that alterations in RAD51B and GID4 are associated with higher telomeric content

Tumor cells need to activate a telomere maintenance mechanism, enabling limitless replication. The bulk of evidence supports that sarcomas predominantly use alternative lengthening of telomeres (ALT) mechanism, commonly associated with alterations in ATRX and DAXX. In our dataset, only 12.3% of sarcomas harbored alterations in these genes. Thus, we checked for the presence of other genomic determinants of high telomeric content in sarcomas. Our dataset consisted of 13555 sarcoma samples, sequenced as a part of routine clinical care on the FoundationOne®Heme platform. We observed a median telomeric content of 622.3 telomeric reads per GC-matched million reads (TRPM) across all samples. In agreement with previous studies, telomeric content was significantly higher in ATRX altered and POT1 altered sarcomas. We further observed that sarcomas with alterations in RAD51B or GID4 were enriched in samples with high telomeric content, specifically within uterus leiomyosarcoma for RAD51B and soft tissue sarcoma (not otherwise specified, nos) for GID4, Furthermore, RAD51B and POT1 alterations were mutually exclusive with ATRX and DAXX alterations, suggestive of functional redundancy. Our results propose a role played by RAD51B and GID4 in telomere elongation in sarcomas and open research opportunities for agents aimed at targeting this critical pathway in tumorigenesis.

Genomic landscape of TP53-mutated myeloid malignancies

TP53-mutated myeloid malignancies are associated with complex cytogenetics and extensive structural variants, which complicates detailed genomic analysis by conventional clinical techniques. We performed whole-genome sequencing (WGS) of 42 acute myeloid leukemia (AML)/myelodysplastic syndromes (MDS) cases with paired normal tissue to better characterize the genomic landscape of TP53-mutated AML/MDS. WGS accurately determines TP53 allele status, a key prognostic factor, resulting in the reclassification of 12% of cases from monoallelic to multihit. Although aneuploidy and chromothripsis are shared with most TP53-mutated cancers, the specific chromosome abnormalities are distinct to each cancer type, suggesting a dependence on the tissue of origin. ETV6 expression is reduced in nearly all cases of TP53-mutated AML/MDS, either through gene deletion or presumed epigenetic silencing. Within the AML cohort, mutations of NF1 are highly enriched, with deletions of 1 copy of NF1 present in 45% of cases and biallelic mutations in 17%. Telomere content is increased in TP53-mutated AMLs compared with other AML subtypes, and abnormal telomeric sequences were detected in the interstitial regions of chromosomes. These data highlight the unique features of TP53-mutated myeloid malignancies, including the high frequency of chromothripsis and structural variation, the frequent involvement of unique genes (including NF1 and ETV6) as cooperating events, and evidence for altered telomere maintenance.

Alternative lengthening of telomeres (ALT) in pediatric high-grade gliomas can occur without ATRX mutation and is enriched in patients with pathogenic germline mismatch repair (MMR) variants

BACKGROUND

To achieve replicative immortality, most cancers develop a telomere maintenance mechanism, such as reactivation of telomerase or alternative lengthening of telomeres (ALT). There are limited data on the prevalence and clinical significance of ALT in pediatric brain tumors, and ALT-directed therapy is not available.

METHODS

We performed C-circle analysis (CCA) on 579 pediatric brain tumors that had corresponding tumor/normal whole genome sequencing through the Open Pediatric Brain Tumor Atlas (OpenPBTA). We detected ALT in 6.9% (n = 40/579) of these tumors and completed additional validation by ultrabright telomeric foci in situ on a subset of these tumors. We used CCA to validate TelomereHunter for computational prediction of ALT status and focus subsequent analyses on pediatric high-grade gliomas (pHGGs) Finally, we examined whether ALT is associated with recurrent somatic or germline alterations.

RESULTS

ALT is common in pHGGs (n = 24/63, 38.1%), but occurs infrequently in other pediatric brain tumors (<3%). Somatic ATRX mutations occur in 50% of ALT+ pHGGs and in 30% of ALT- pHGGs. Rare pathogenic germline variants in mismatch repair (MMR) genes are significantly associated with an increased occurrence of ALT.

CONCLUSIONS

We demonstrate that ATRX is mutated in only a subset of ALT+ pHGGs, suggesting other mechanisms of ATRX loss of function or alterations in other genes may be associated with the development of ALT in these patients. We show that germline variants in MMR are associated with the development of ALT in patients with pHGG.

Integrated metabolic and genetic analysis reveals distinct features of human differentiated thyroid cancer

BACKGROUND

Differentiated thyroid cancer (DTC) affects thousands of lives worldwide each year. Typically, DTC is a treatable disease with a good prognosis. Yet, some patients are subjected to partial or total thyroidectomy and radioiodine therapy to prevent local disease recurrence and metastasis. Unfortunately, thyroidectomy and/or radioiodine therapy often worsen(s) quality of life and might be unnecessary in indolent DTC cases. On the other hand, the lack of biomarkers indicating a potential metastatic thyroid cancer imposes an additional challenge to managing and treating patients with this disease.

AIM

The presented clinical setting highlights the unmet need for a precise molecular diagnosis of DTC and potential metastatic disease, which should dictate appropriate therapy.

MATERIALS AND METHODS

In this article, we present a differential multi-omics model approach, including metabolomics, genomics, and bioinformatic models, to distinguish normal glands from thyroid tumours. Additionally, we are proposing biomarkers that could indicate potential metastatic diseases in papillary thyroid cancer (PTC), a sub-class of DTC.

RESULTS

Normal and tumour thyroid tissue from DTC patients had a distinct yet well-defined metabolic profile with high levels of anabolic metabolites and/or other metabolites associated with the energy maintenance of tumour cells. The consistency of the DTC metabolic profile allowed us to build a bioinformatic classification model capable of clearly distinguishing normal from tumor thyroid tissues, which might help diagnose thyroid cancer. Moreover, based on PTC patient samples, our data suggest that elevated nuclear and mitochondrial DNA mutational burden, intra-tumour heterogeneity, shortened telomere length, and altered metabolic profile reflect the potential for metastatic disease.

DISCUSSION

Altogether, this work indicates that a differential and integrated multi-omics approach might improve DTC management, perhaps preventing unnecessary thyroid gland removal and/or radioiodine therapy.

CONCLUSIONS

Well-designed, prospective translational clinical trials will ultimately show the value of this integrated multi-omics approach and early diagnosis of DTC and potential metastatic PTC.

Telomere Length Changes in Cancer: Insights on Carcinogenesis and Potential for Non-Invasive Diagnostic Strategies

Telomere dynamics play a crucial role in the maintenance of chromosome integrity; changes in telomere length may thus contribute to the development of various diseases including cancer. Understanding the role of telomeric DNA in carcinogenesis and detecting the presence of cell-free telomeric DNA (cf-telDNA) in body fluids offer a potential biomarker for novel cancer screening and diagnostic strategies. Liquid biopsy is becoming increasingly popular due to its undeniable benefits over conventional invasive methods. However, the organization and function of cf-telDNA in the extracellular milieu are understudied. This paper provides a review based on 3,398,017 cancer patients, patients with other conditions, and control individuals with the aim to shed more light on the inconsistent nature of telomere lengthening/shortening in oncological contexts. To gain a better understanding of biological factors (e.g., telomerase activation, alternative lengthening of telomeres) affecting telomere homeostasis across different types of cancer, we summarize mechanisms responsible for telomere length maintenance. In conclusion, we compare tissue- and liquid biopsy-based approaches in cancer assessment and provide a brief outlook on the methodology used for telomere length evaluation, highlighting the advances of state-of-the-art approaches in the field.

Large-Scale Detection of Telomeric Motif Sequences in Genomic Data Using TelFinder

Telomeres are regions of tandem repeated sequences at the ends of linear chromosomes that protect against DNA damage and chromosome fusion. Telomeres are associated with senescence and cancers and have attracted the attention of an increasing number of researchers. However, few telomeric motif sequences are known. Given the mounting interest in telomeres, an efficient computational tool for the de novo detection of the telomeric motif sequence of new species is needed since experimental-based methods are costly in terms of time and effort. Here, we report the development of TelFinder, an easy-to-use and freely available tool for the de novo detection of telomeric motif sequences from genomic data. The vast quantity of readily available genomic data makes it possible to apply this tool to any species of interest, which will undoubtedly inspire studies requiring telomeric repeat information and improve the utilization of these genomic data sets. We have tested TelFinder on telomeric sequences available in the Telomerase Database, and the detection accuracy reaches 90%. In addition, variation analyses in telomere sequences can be performed by TelFinder for the first time. The telomere variation preference of different chromosomes and even at the ends of the chromosome can provide clues regarding the underlying mechanisms of telomeres. Overall, these results shed new light on the divergent evolution of telomeres. IMPORTANCE Telomeres are reported to be highly correlated with the cell cycle and aging. As a result, research on telomere composition and evolution has become more and more urgent. However, using experimental methods to detect telomeric motif sequences is slow and costly. To combat this challenge, we developed TelFinder, a computational tool for the de novo detection of the telomere composition only using genomic data. In this study, we showed that a lot of complicated telomeric motifs could be identified by TelFinder only using genomic data. In addition, TelFinder can be used to check variation analyses in telomere sequences, which could lead to a deeper understanding of telomere sequences.

Colorectal cancer patient-derived organoids and cell lines harboring ATRX and/or DAXX mutations lack Alternative Lengthening of Telomeres (ALT)

Telomere maintenance is necessary to maintain cancer cell unlimited viability. However, the mechanisms maintaining telomere length in colorectal cancer (CRC) have not been extensively investigated. Telomere maintenance mechanisms (TMM) include the re-expression of telomerase or alternative lengthening of telomeres (ALT). ALT is genetically associated with somatic alterations in alpha-thalassemia/mental retardation X-linked (ATRX) and death domain-associated protein (DAXX) genes. Cells displaying ALT present distinctive features including C-circles made of telomeric DNA, long and heterogenous telomeric tracts, and telomeric DNA co-localized with promyelocytic leukemia (PML) bodies forming so-called ALT-associated PML bodies (APBs). Here, we identified mutations in ATRX and/or DAXX genes in an extensive collection of CRC samples including 119 patient-derived organoids (PDOs) and 232 established CRC cell lines. C-circles measured in CRC PDOs and cell lines showed low levels overall. We also observed that CRC PDOs and cell lines did not display a significant accumulation of APBs or long telomeres with no appreciable differences between wild-type and mutated ATRX/DAXX samples. Overall, our extensive analyses indicate that CRC is not prone to engage ALT, even when carrying genetic lesions in ATRX and/or DAXX, and support the notion that ATRX/DAXX genomic footprints are not reliable predictors of ALT.

Reprint of: The Diagnostic and Prognostic Utility of Incorporating DAXX, ATRX, and Alternative Lengthening of Telomeres (ALT) to the Evaluation of Pancreatic Neuroendocrine Tumors (PanNETs)

Pancreatic neuroendocrine tumors (PanNETs) are a heterogeneous group of neoplasms with increasing incidence and an ill-defined pathobiology. Although many PanNETs are indolent and remain stable for years, a subset may behave aggressively and metastasize widely. Thus, the increasing and frequent detection of PanNETs presents a treatment dilemma. Current prognostic systems are susceptible to interpretation errors, sampling issues, and do not accurately reflect the clinical behavior of these neoplasms. Hence, additional biomarkers are needed to improve the prognostic stratification of patients diagnosed with a PanNET. Recent studies have identified alterations in death domain-associated protein 6 (DAXX) and alpha-thalassemia/mental retardation X-linked (ATRX), as well as alternative lengthening of telomeres (ALT), as promising prognostic biomarkers. This review summarizes the identification, clinical utility, and specific nuances in testing for DAXX/ATRX by immunohistochemistry and ALT by telomere-specific fluorescence in situ hybridization in PanNETs. Furthermore, a discussion on diagnostic indications for DAXX, ATRX, and ALT status is provided to include the distinction between PanNETs and pancreatic neuroendocrine carcinomas (PanNECs), and determining pancreatic origin for metastatic neuroendocrine tumors in the setting of an unknown primary.

Genomic landscape of TP53 -mutated myeloid malignancies

TP53 -mutated myeloid malignancies are most frequently associated with complex cytogenetics. The presence of complex and extensive structural variants complicates detailed genomic analysis by conventional clinical techniques. We performed whole genome sequencing of 42 AML/MDS cases with paired normal tissue to characterize the genomic landscape of TP53 -mutated myeloid malignancies. The vast majority of cases had multi-hit involvement at the TP53 genetic locus (94%), as well as aneuploidy and chromothripsis. Chromosomal patterns of aneuploidy differed significantly from TP53 -mutated cancers arising in other tissues. Recurrent structural variants affected regions that include ETV6 on chr12p, RUNX1 on chr21, and NF1 on chr17q. Most notably for ETV6 , transcript expression was low in cases of TP53 -mutated myeloid malignancies both with and without structural rearrangements involving chromosome 12p. Telomeric content is increased in TP53 -mutated AML/MDS compared other AML subtypes, and telomeric content was detected adjacent to interstitial regions of chromosomes. The genomic landscape of TP53 -mutated myeloid malignancies reveals recurrent structural variants affecting key hematopoietic transcription factors and telomeric repeats that are generally not detected by panel sequencing or conventional cytogenetic analyses.

Key Points

WGS comprehensively determines TP53 mutation status, resulting in the reclassification of 12% of cases from mono-allelic to multi-hit Chromothripsis is more frequent than previously appreciated, with a preference for specific chromosomes ETV6 is deleted in 45% of cases, with evidence for epigenetic suppression in non-deleted cases NF1 is mutated in 48% of cases, with multi-hit mutations in 17% of these cases TP53 -mutated AML/MDS is associated with altered telomere content compared with other AMLs.

Whole-genome sequencing of chronic lymphocytic leukemia identifies subgroups with distinct biological and clinical features

The value of genome-wide over targeted driver analyses for predicting clinical outcomes of cancer patients is debated. Here, we report the whole-genome sequencing of 485 chronic lymphocytic leukemia patients enrolled in clinical trials as part of the United Kingdom's 100,000 Genomes Project. We identify an extended catalog of recurrent coding and noncoding genetic mutations that represents a source for future studies and provide the most complete high-resolution map of structural variants, copy number changes and global genome features including telomere length, mutational signatures and genomic complexity. We demonstrate the relationship of these features with clinical outcome and show that integration of 186 distinct recurrent genomic alterations defines five genomic subgroups that associate with response to therapy, refining conventional outcome prediction. While requiring independent validation, our findings highlight the potential of whole-genome sequencing to inform future risk stratification in chronic lymphocytic leukemia.

The diagnostic and prognostic utility of incorporating DAXX, ATRX, and alternative lengthening of telomeres to the evaluation of pancreatic neuroendocrine tumors

Pancreatic neuroendocrine tumors (PanNETs) are a heterogeneous group of neoplasms with increasing incidence and an ill-defined pathobiology. Although many PanNETs are indolent and remain stable for years, a subset may behave aggressively and metastasize widely. Thus, the increasing and frequent detection of PanNETs presents a treatment dilemma. Current prognostic systems are susceptible to interpretation errors, sampling issues, and do not accurately reflect the clinical behavior of these neoplasms. Hence, additional biomarkers are needed to improve the prognostic stratification of patients diagnosed with a PanNET. Recent studies have identified alterations in death domain-associated protein 6 (DAXX) and alpha-thalassemia/mental retardation X-linked (ATRX), as well as alternative lengthening of telomeres (ALT), as promising prognostic biomarkers. This review summarizes the identification, clinical utility, and specific nuances in testing for DAXX/ATRX by immunohistochemistry and ALT by telomere-specific fluorescence in situ hybridization in PanNETs. Furthermore, a discussion on diagnostic indications for DAXX, ATRX, and ALT status is provided to include the distinction between PanNETs and pancreatic neuroendocrine carcinomas (PanNECs), and determining pancreatic origin for metastatic neuroendocrine tumors in the setting of an unknown primary.

Reliable assessment of telomere maintenance mechanisms in neuroblastoma

BACKGROUND

Telomere maintenance mechanisms (TMM) are a hallmark of high-risk neuroblastoma, and are conferred by activation of telomerase or alternative lengthening of telomeres (ALT). However, detection of TMM is not yet part of the clinical routine, and consensus on TMM detection, especially on ALT assessment, remains to be achieved.

METHODS

Whole genome sequencing (WGS) data of 68 primary neuroblastoma samples were analyzed. Telomere length was calculated from WGS data or by telomere restriction fragment analysis (n = 39). ALT was assessed by C-circle assay (CCA, n = 67) and detection of ALT-associated PML nuclear bodies (APB) by combined fluorescence in situ hybridization and immunofluorescence staining (n = 68). RNA sequencing was performed (n = 64) to determine expression of TERT and telomeric long non-coding RNA (TERRA). Telomerase activity was examined by telomerase repeat amplification protocol (TRAP, n = 15).

RESULTS

Tumors were considered as telomerase-positive if they harbored a TERT rearrangement, MYCN amplification or high TERT expression (45.6%, 31/68), and ALT-positive if they were positive for APB and CCA (19.1%, 13/68). If all these markers were absent, tumors were considered TMM-negative (25.0%, 17/68). According to these criteria, the majority of samples were classified unambiguously (89.7%, 61/68). Assessment of additional ALT-associated parameters clarified the TMM status of the remaining seven cases with high likelihood: ALT-positive tumors had higher TERRA expression, longer telomeres, more telomere insertions, a characteristic pattern of telomere variant repeats, and were associated with ATRX mutations.

CONCLUSIONS

We here propose a workflow to reliably detect TMM in neuroblastoma. We show that unambiguous classification is feasible following a stepwise approach that determines both, activation of telomerase and ALT. The workflow proposed in this study can be used in clinical routine and provides a framework to systematically and reliably determine telomere maintenance mechanisms for risk stratification and treatment allocation of neuroblastoma patients.

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.

Inherited MUTYH mutations cause elevated somatic mutation rates and distinctive mutational signatures in normal human cells

Cellular DNA damage caused by reactive oxygen species is repaired by the base excision repair (BER) pathway which includes the DNA glycosylase MUTYH. Inherited biallelic MUTYH mutations cause predisposition to colorectal adenomas and carcinoma. However, the mechanistic progression from germline MUTYH mutations to MUTYH-Associated Polyposis (MAP) is incompletely understood. Here, we sequence normal tissue DNAs from 10 individuals with MAP. Somatic base substitution mutation rates in intestinal epithelial cells were elevated 2 to 4-fold in all individuals, except for one showing a 31-fold increase, and were also increased in other tissues. The increased mutation burdens were of multiple mutational signatures characterised by C > A changes. Different mutation rates and signatures between individuals are likely due to different MUTYH mutations or additional inherited mutations in other BER pathway genes. The elevated base substitution rate in normal cells likely accounts for the predisposition to neoplasia in MAP. Despite ubiquitously elevated mutation rates, individuals with MAP do not display overt evidence of premature ageing. Thus, accumulation of somatic mutations may not be sufficient to cause the global organismal functional decline of ageing.

The genomic and transcriptional landscape of primary central nervous system lymphoma

Primary lymphomas of the central nervous system (PCNSL) are mainly diffuse large B-cell lymphomas (DLBCLs) confined to the central nervous system (CNS). Molecular drivers of PCNSL have not been fully elucidated. Here, we profile and compare the whole-genome and transcriptome landscape of 51 CNS lymphomas (CNSL) to 39 follicular lymphoma and 36 DLBCL cases outside the CNS. We find recurrent mutations in JAK-STAT, NFkB, and B-cell receptor signaling pathways, including hallmark mutations in MYD88 L265P (67%) and CD79B (63%), and CDKN2A deletions (83%). PCNSLs exhibit significantly more focal deletions of HLA-D (6p21) locus as a potential mechanism of immune evasion. Mutational signatures correlating with DNA replication and mitosis are significantly enriched in PCNSL. TERT gene expression is significantly higher in PCNSL compared to activated B-cell (ABC)-DLBCL. Transcriptome analysis clearly distinguishes PCNSL and systemic DLBCL into distinct molecular subtypes. Epstein-Barr virus (EBV)+ CNSL cases lack recurrent mutational hotspots apart from IG and HLA-DRB loci. We show that PCNSL can be clearly distinguished from DLBCL, having distinct expression profiles, IG expression and translocation patterns, as well as specific combinations of genetic alterations.

Molecular Insights and Prognosis Associated With RBM8A in Glioblastoma

Background: Glioblastoma (GBM) is the most invasive brain tumors, and it is associated with high rates of recurrence and mortality. The purpose of this study was to investigate the expression of RBM8A in GBM and the potential influence of its expression on the disease.

Methods: Levels of RBM8A mRNA in GBM patients and controls were examined in The Cancer Genome Atlas (TCGA), GSE16011 and GSE90604 databases. GBM samples in TCGA were divided into RBM8A^high and RBM8A^low groups. Differentially expressed genes (DEGs) between GBM patients and controls were identified, as were DEGs between RBM8A^high and RBM8A^low groups. DEGs common to both of these comparisons were analyzed for coexpression and regression analyses. In addition, we identified potential effects of RBM8A on competing endogenous RNAs, immune cell infiltration, methylation modifications, and somatic mutations.

Results: RBM8A is expressed at significantly higher levels in GBM than control samples, and its level correlates with tumor purity. We identified a total of 488 mRNAs that differed between GBM and controls as well as between RBM8A^high and RBM8A^low groups, which enrichment analysis revealed to be associated mainly with neuroblast proliferation, and T cell immune responses. We identified 174 mRNAs that gave areas under the receiver operating characteristic curve >0.7 among coexpression module genes, of which 13 were significantly associated with overall survival of GBM patients. We integrated 11 candidate mRNAs through LASSO algorithm, then nomogram, risk score, and decision curve analyses were analyzed. We found that RBM8A may compete with DLEU1 for binding to miR-128-1-5p, and aberrant RBM8A expression was associations with tumor infiltration by immune cells. Some mRNAs associated with GBM prognosis also appear to be methylated or mutated.

Conclusions: Our study strongly links RBM8A expression to GBM pathobiology and patient prognosis. The candidate mRNAs identified here may lead to therapeutic targets against the disease.

Profiles of telomeric repeats in Insecta reveal diverse forms of telomeric motifs in Hymenopterans

Telomeres consist of highly conserved simple tandem telomeric repeat motif (TRM): (TTAGG)n in arthropods, (TTAGGG)n in vertebrates, and (TTTAGGG)n in most plants. TRM can be detected from chromosome-level assembly, which typically requires long-read sequencing data. To take advantage of short-read data, we developed an ultra-fast Telomeric Repeats Identification Pipeline and evaluated its performance on 91 species. With proven accuracy, we applied Telomeric Repeats Identification Pipeline in 129 insect species, using 7 Tbp of short-read sequences. We confirmed (TTAGG)n as the TRM in 19 orders, suggesting it is the ancestral form in insects. Systematic profiling in Hymenopterans revealed a diverse range of TRMs, including the canonical 5-bp TTAGG (bees, ants, and basal sawflies), three independent losses of tandem repeat form TRM (Ichneumonoids, hunting wasps, and gall-forming wasps), and most interestingly, a common 8-bp (TTATTGGG)n in Chalcid wasps with two 9-bp variants in the miniature wasp (TTACTTGGG) and fig wasps (TTATTGGGG). Our results identified extraordinary evolutionary fluidity of Hymenopteran TRMs, and rapid evolution of TRM and repeat abundance at all evolutionary scales, providing novel insights into telomere evolution.

Retinoblastoma protein regulates carcinogen susceptibility at heterochromatic cancer driver loci

Carcinogenic insult, such as UV light exposure, creates DNA lesions that evolve into mutations if left unrepaired. These resulting mutations can contribute to carcinogenesis and drive malignant phenotypes. Susceptibility to carcinogens (i.e., the propensity to form a carcinogen-induced DNA lesion) is regulated by both genetic and epigenetic factors. Importantly, carcinogen susceptibility is a critical contributor to cancer mutagenesis. It is known that mutations can be prevented by tumor suppressor regulation of DNA damage response pathways; however, their roles carcinogen susceptibility have not yet been reported. In this study, we reveal that the retinoblastoma (RB1) tumor suppressor regulates UV susceptibility across broad regions of the genome. In particular, centromere and telomere-proximal regions exhibit significant increases in UV lesion susceptibility when RB1 is deleted. Several cancer-related genes are located within genomic regions of increased susceptibility, including telomerase reverse transcriptase, TERT, thereby accelerating mutagenic potential in cancers with RB1 pathway alterations. These findings reveal novel genome stability mechanisms of a tumor suppressor and uncover new pathways to accumulate mutations during cancer evolution.

A pan-cancer landscape of telomeric content shows that RAD21 and HGF alterations are associated with longer telomeres

Background

Cancer cells can proliferate indefinitely through telomere maintenance mechanisms. These mechanisms include telomerase-dependent elongation, mediated by TERT activation, and alternative lengthening of telomeres (ALT), linked to loss of ATRX or DAXX.

Methods

We analyzed the telomeric content of 89,959 tumor samples within the Foundation Medicine dataset and investigated the genomic determinants of high telomeric content, linking them to clinical outcomes, when available.

Results

Telomeric content varied widely by disease type with leiomyosarcoma having the highest and Merkel cell carcinoma having the lowest telomeric content. In agreement with previous studies, telomeric content was significantly higher in samples with alterations in TERC, ATRX, and DAXX. We further identified that amplifications in two genes, RAD21 and HGF, were enriched in samples with high telomeric content, which was confirmed using the PCAWG/ICGC dataset. We identified the minimal amplified region associated with high telomeric content for RAD21 (8q23.1–8q24.12), which excludes MYC, and for HGF (7q21.11). Our results demonstrated that RAD21 and HGF exerted an additive telomere lengthening effect on samples with existing alterations in canonical genes previously associated with telomere elongation. Furthermore, patients with breast cancer who harbor RAD21 alterations had poor median overall survival and trended towards higher levels of Ki-67 staining.

Conclusions

This study highlights the importance of the role played by RAD21 (8q23.1–8q24.12) and HGF (7q21.11) in the lengthening of telomeres, supporting unlimited replication in tumors. These findings open avenues for work aimed at targeting this crucial pathway in tumorigenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01029-7.

qmotif: determination of telomere content from whole-genome sequence data

Motivation

Changes in telomere length have been observed in cancer and can be indicative of mechanisms involved in carcinogenesis. Most methods used to estimate telomere length require laboratory analysis of DNA samples. Here, we present qmotif, a fast and easy tool that determines telomeric repeat sequences content as an estimate of telomere length directly from whole-genome sequencing.

Results

qmotif shows similar results to quantitative PCR, the standard method for high-throughput clinical telomere length quantification. qmotif output correlates strongly with the output of other tools for determining telomere sequence content, TelSeq and TelomereHunter, but can run in a fraction of the time-usually under a minute.

Availability and implementation

qmotif is implemented in Java and source code is available at https://github.com/AdamaJava/adamajava, with instructions on how to build and use the application available from https://adamajava.readthedocs.io/en/latest/.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Telogator: a method for reporting chromosome-specific telomere lengths from long reads

MOTIVATION

Telomeres are the repetitive sequences found at the ends of eukaryotic chromosomes and are often thought of as a 'biological clock,' with their average length shortening during division in most cells. In addition to their association with senescence, abnormal telomere lengths are well known to be associated with multiple cancers, short telomere syndromes and as risk factors for a broad range of diseases. While a majority of methods for measuring telomere length will report average lengths across all chromosomes, it is known that aberrations in specific chromosome arms are biomarkers for certain diseases. Due to their repetitive nature, characterizing telomeres at this resolution is prohibitive for short read sequencing approaches, and is challenging still even with longer reads.

RESULTS

We present Telogator: a method for reporting chromosome-specific telomere length from long read sequencing data. We demonstrate Telogator's sensitivity in detecting chromosome-specific telomere length in simulated data across a range of read lengths and error rates. Telogator is then applied to 10 germline samples, yielding a high correlation with short read methods in reporting average telomere length. In addition, we investigate common subtelomere rearrangements and identify the minimum read length required to anchor telomere/subtelomere boundaries in samples with these haplotypes.

AVAILABILITY AND IMPLEMENTATION

Telogator is written in Python3 and is available at github.com/zstephens/telogator.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Novel Tumor-Specific Antigens for Immunotherapy Identified From Multi-omics Profiling in Thymic Carcinomas

Thymic carcinoma (TC) is the most aggressive thymic epithelial neoplasm. TC patients with microsatellite instability, whole-genome doubling, or alternative tumor-specific antigens from gene fusion are most likely to benefit from immunotherapies. However, due to the rarity of this disease, how to prioritize the putative biomarkers and what constitutes an optimal treatment regimen remains largely unknown. Therefore, we integrated genomic and transcriptomic analyses from TC patients and revealed that frameshift indels in KMT2C and CYLD frequently produce neoantigens. Moreover, a median of 3 fusion-derived neoantigens was predicted across affected patients, especially the CATSPERB-TC2N neoantigens that were recurrently predicted in TC patients. Lastly, potentially actionable alterations with early levels of evidence were uncovered and could be used for designing clinical trials. In summary, this study shed light on our understanding of tumorigenesis and presented new avenues for molecular characterization and immunotherapy in TC.

The telomere length landscape of prostate cancer

Replicative immortality is a hallmark of cancer, and can be achieved through telomere lengthening and maintenance. Although the role of telomere length in cancer has been well studied, its association to genomic features is less well known. Here, we report the telomere lengths of 392 localized prostate cancer tumours and characterize their relationship to genomic, transcriptomic and proteomic features. Shorter tumour telomere lengths are associated with elevated genomic instability, including single-nucleotide variants, indels and structural variants. Genes involved in cell proliferation and signaling are correlated with tumour telomere length at all levels of the central dogma. Telomere length is also associated with multiple clinical features of a tumour. Longer telomere lengths in non-tumour samples are associated with a lower rate of biochemical relapse. In summary, we describe the multi-level integration of telomere length, genomics, transcriptomics and proteomics in localized prostate cancer.

Neuroblastoma-Telomere maintenance, deregulated signaling transduction and beyond

The childhood malignancy neuroblastoma belongs to the group of embryonal tumors and originates from progenitor cells of the sympathoadrenal lineage. Treatment options for children with high-risk and relapsed disease are still very limited. In recent years, an ever-growing molecular diversity was identified using (epi)-genetic profiling of neuroblastoma tumors, indicating that molecularly targeted therapies could be a promising therapeutic option. In this review article, we summarize the various molecular subtypes and genetic events associated with neuroblastoma and describe recent advances in targeted therapies. We lay a strong emphasis on the importance of telomere maintenance mechanisms for understanding tumor progression and risk classification of neuroblastoma.

Increased somatic mutation burdens in normal human cells due to defective DNA polymerases

Mutation accumulation in somatic cells contributes to cancer development and is proposed as a cause of aging. DNA polymerases Pol ε and Pol δ replicate DNA during cell division. However, in some cancers, defective proofreading due to acquired POLE/POLD1 exonuclease domain mutations causes markedly elevated somatic mutation burdens with distinctive mutational signatures. Germline POLE/POLD1 mutations cause familial cancer predisposition. Here, we sequenced normal tissue and tumor DNA from individuals with germline POLE/POLD1 mutations. Increased mutation burdens with characteristic mutational signatures were found in normal adult somatic cell types, during early embryogenesis and in sperm. Thus human physiology can tolerate ubiquitously elevated mutation burdens. Except for increased cancer risk, individuals with germline POLE/POLD1 mutations do not exhibit overt features of premature aging. These results do not support a model in which all features of aging are attributable to widespread cell malfunction directly resulting from somatic mutation burdens accrued during life.

Integrated evaluation of telomerase activation and telomere maintenance across cancer cell lines

In cancer, telomere maintenance is critical for the development of replicative immortality. Using genome sequences from the Cancer Cell Line Encyclopedia and Genomics of Drug Sensitivity in Cancer Project, we calculated telomere content across 1299 cancer cell lines. We find that telomerase reverse transcriptase (TERT) expression correlates with telomere content in lung, central nervous system, and leukemia cell lines. Using CRISPR/Cas9 screening data, we show that lower telomeric content is associated with dependency of CST telomere maintenance genes. Increased dependencies of shelterin members are associated with wild-type TP53 status. Investigating the epigenetic regulation of TERT, we find widespread allele-specific expression in promoter-wildtype contexts. TERT promoter-mutant cell lines exhibit hypomethylation at PRC2-repressed regions, suggesting a cooperative global epigenetic state in the reactivation of telomerase. By incorporating telomere content with genomic features across comprehensively characterized cell lines, we provide further insights into the role of telomere regulation in cancer immortality.

Genomic and evolutionary classification of lung cancer in never smokers

Lung cancer in never smokers (LCINS) is a common cause of cancer mortality but its genomic landscape is poorly characterized. Here high-coverage whole-genome sequencing of 232 LCINS showed 3 subtypes defined by copy number aberrations. The dominant subtype (piano), which is rare in lung cancer in smokers, features somatic UBA1 mutations, germline AR variants and stem cell-like properties, including low mutational burden, high intratumor heterogeneity, long telomeres, frequent KRAS mutations and slow growth, as suggested by the occurrence of cancer drivers' progenitor cells many years before tumor diagnosis. The other subtypes are characterized by specific amplifications and EGFR mutations (mezzo-forte) and whole-genome doubling (forte). No strong tobacco smoking signatures were detected, even in cases with exposure to secondhand tobacco smoke. Genes within the receptor tyrosine kinase-Ras pathway had distinct impacts on survival; five genomic alterations independently doubled mortality. These findings create avenues for personalized treatment in LCINS.

aniFOUND: analysing the associated proteome and genomic landscape of the repaired nascent non-replicative chromatin

Specific capture of chromatin fractions with distinct and well-defined features has emerged as both challenging and a key strategy towards a comprehensive understanding of genome biology. In this context, we developed aniFOUND (accelerated native isolation of factors on unscheduled nascent DNA), an antibody-free method, which can label, capture, map and characterise nascent chromatin fragments that are synthesized in response to specific cues outside S-phase. We used the 'unscheduled' DNA synthesis (UDS) that takes place during the repair of UV-induced DNA lesions and coupled the captured chromatin to high-throughput analytical technologies. By mass-spectrometry we identified several factors with no previously known role in UVC-DNA damage response (DDR) as well as known DDR proteins. We experimentally validated the repair-dependent recruitment of the chromatin remodeller RSF1 and the cohesin-loader NIPBL at sites of UVC-induced photolesions. Developing aniFOUND-seq, a protocol for mapping UDS activity with high resolution, allowed us to monitor the landscape of UVC repair-synthesis events genome wide. We further resolved repair efficacy of the rather unexplored repeated genome, in particular rDNA and telomeres. In summary, aniFOUND delineates the proteome composition and genomic landscape of chromatin loci with specific features by integrating state-of-the-art 'omics' technologies to promote a comprehensive view of their function.

Higher-Order Chromatin Structures of Chromosomally Integrated HHV-6A Predict Integration Sites

Human herpesvirus -6A and 6B (HHV-6A/B) can integrate their genomes into the telomeres of human chromosomes. Viral integration can occur in several cell types, including germinal cells, resulting in individuals that harbor the viral genome in every cell of their body. The integrated genome is efficiently silenced but can sporadically reactivate resulting in various clinical symptoms. To date, the integration mechanism and the subsequent silencing of HHV-6A/B genes remains poorly understood. Here we investigate the genome-wide chromatin contacts of the integrated HHV-6A in latently-infected cells. We show that HHV-6A becomes transcriptionally silent upon infection of these cells over the course of seven days. In addition, we established an HHV-6-specific 4C-seq approach, revealing that the HHV-6A 3D interactome is associated with quiescent chromatin states in cells harboring integrated virus. Furthermore, we observed that the majority of virus chromatin interactions occur toward the distal ends of specific human chromosomes. Exploiting this finding, we established a 4C-seq method that accurately detects the chromosomal integration sites. We further implement long-read minION sequencing in the 4C-seq assay and developed a method to identify HHV-6A/B integration sites in clinical samples.

Alternative lengthening of telomeres in childhood neuroblastoma from genome to proteome

Telomere maintenance by telomerase activation or alternative lengthening of telomeres (ALT) is a major determinant of poor outcome in neuroblastoma. Here, we screen for ALT in primary and relapsed neuroblastomas (n = 760) and characterize its features using multi-omics profiling. ALT-positive tumors are molecularly distinct from other neuroblastoma subtypes and enriched in a population-based clinical sequencing study cohort for relapsed cases. They display reduced ATRX/DAXX complex abundance, due to either ATRX mutations (55%) or low protein expression. The heterochromatic histone mark H3K9me3 recognized by ATRX is enriched at the telomeres of ALT-positive tumors. Notably, we find a high frequency of telomeric repeat loci with a neuroblastoma ALT-specific hotspot on chr1q42.2 and loss of the adjacent chromosomal segment forming a neo-telomere. ALT-positive neuroblastomas proliferate slowly, which is reflected by a protracted clinical course of disease. Nevertheless, children with an ALT-positive neuroblastoma have dismal outcome.

Telomeres reforged with non-telomeric sequences in mouse embryonic stem cells

Telomeres are part of a highly refined system for maintaining the stability of linear chromosomes. Most telomeres rely on simple repetitive sequences and telomerase enzymes to protect chromosomal ends; however, in some species or telomerase-defective situations, an alternative lengthening of telomeres (ALT) mechanism is used. ALT mainly utilises recombination-based replication mechanisms and the constituents of ALT-based telomeres vary depending on models. Here we show that mouse telomeres can exploit non-telomeric, unique sequences in addition to telomeric repeats. We establish that a specific subtelomeric element, the mouse template for ALT (mTALT), is used for repairing telomeric DNA damage as well as for composing portions of telomeres in ALT-dependent mouse embryonic stem cells. Epigenomic and proteomic analyses before and after ALT activation reveal a high level of non-coding mTALT transcripts despite the heterochromatic nature of mTALT-based telomeres. After ALT activation, the increased HMGN1, a non-histone chromosomal protein, contributes to the maintenance of telomere stability by regulating telomeric transcription. These findings provide a molecular basis to study the evolution of new structures in telomeres.

Telomeres can be maintained by a telomerase-independent mechanism called an alternative lengthening of telomeres (ALT). Here the authors use mouse Terc (telomerase RNA) knockout embryonic cells and provide longitudinal analysis of ALT telomeres maintained with non-telomeric sequences.

Method comparison studies of telomere length measurement using qPCR approaches: A critical appraisal of the literature

Use of telomere length (TL) as a biomarker for various environmental exposures and diseases has increased in recent years. Various methods have been developed to measure telomere length. Polymerase chain reaction (PCR)-based methods remain wide-spread for population-based studies due to the high-throughput capability. While several studies have evaluated the repeatability and reproducibility of different TL measurement methods, the results have been variable. We conducted a literature review of TL measurement cross-method comparison studies that included a PCR-based method published between January 1, 2002 and May 25, 2020. A total of 25 articles were found that matched the inclusion criteria. Papers were reviewed for quality of methodologic reporting of sample and DNA quality, PCR assay characteristics, sample blinding, and analytic approaches to determine final TL. Overall, methodologic reporting was low as assessed by two different reporting guidelines for qPCR-based TL measurement. There was a wide range in the reported correlation between methods (as assessed by Pearson’s r) and few studies utilized the recommended intra-class correlation coefficient (ICC) for assessment of assay repeatability and methodologic comparisons. The sample size for nearly all studies was less than 100, raising concerns about statistical power. Overall, this review found that the current literature on the relation between TL measurement methods is lacking in validity and scientific rigor. In light of these findings, we present reporting guidelines for PCR-based TL measurement methods and results of analyses of the effect of assay repeatability (ICC) on statistical power of cross-sectional and longitudinal studies. Additional cross-laboratory studies with rigorous methodologic and statistical reporting, adequate sample size, and blinding are essential to accurately determine assay repeatability and replicability as well as the relation between TL measurement methods.

Bioinformatics workflows for clinical applications in precision oncology

High-throughput molecular profiling of tumors is a fundamental aspect of precision oncology, enabling the identification of genomic alterations that can be targeted therapeutically. In this context, a patient is matched to a specific drug or therapy based on the tumor's underlying genetic driver events rather than the histologic classification. This approach requires extensive bioinformatics methodology and workflows, including raw sequencing data processing and quality control, variant calling and annotation, integration of different molecular data types, visualization and finally reporting the data to physicians, cancer researchers and pharmacologists in a format that is readily interpretable for clinical decision making. This review comprises a broad overview of these bioinformatics aspects and discusses the multiple analytical, technical and interpretational challenges that remain to efficiently translate molecular findings into personalized treatment recommendations.

The landscape of chromothripsis across adult cancer types

Chromothripsis is a recently identified mutational phenomenon, by which a presumably single catastrophic event generates extensive genomic rearrangements of one or a few chromosome(s). Considered as an early event in tumour development, this form of genome instability plays a prominent role in tumour onset. Chromothripsis prevalence might have been underestimated when using low-resolution methods, and pan-cancer studies based on sequencing are rare. Here we analyse chromothripsis in 28 tumour types covering all major adult cancers (634 tumours, 316 whole-genome and 318 whole-exome sequences). We show that chromothripsis affects a substantial proportion of human cancers, with a prevalence of 49% across all cases. Chromothripsis generates entity-specific genomic alterations driving tumour development, including clinically relevant druggable fusions. Chromothripsis is linked with specific telomere patterns and univocal mutational signatures in distinct tumour entities. Longitudinal analysis of chromothriptic patterns in 24 matched tumour pairs reveals insights in the clonal evolution of tumours with chromothripsis.

The shattering of chromosomes is a dramatic early event in tumourigenesis and is termed chromothripsis. Here, the authors examine chromothripsis across 28 tumour types and show that 49% of cancers exhibit features of chromothripsis.

Genomic footprints of activated telomere maintenance mechanisms in cancer

Cancers require telomere maintenance mechanisms for unlimited replicative potential. They achieve this through TERT activation or alternative telomere lengthening associated with ATRX or DAXX loss. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, we dissect whole-genome sequencing data of over 2500 matched tumor-control samples from 36 different tumor types aggregated within the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium to characterize the genomic footprints of these mechanisms. While the telomere content of tumors with ATRX or DAXX mutations (ATRX/DAXXtrunc) is increased, tumors with TERT modifications show a moderate decrease of telomere content. One quarter of all tumor samples contain somatic integrations of telomeric sequences into non-telomeric DNA. This fraction is increased to 80% prevalence in ATRX/DAXXtrunc tumors, which carry an aberrant telomere variant repeat (TVR) distribution as another genomic marker. The latter feature includes enrichment or depletion of the previously undescribed singleton TVRs TTCGGG and TTTGGG, respectively. Our systematic analysis provides new insight into the recurrent genomic alterations associated with telomere maintenance mechanisms in cancer.