Comparative analysis of whole genome sequencing-based telomere length measurement techniques

High resolution long-read telomere sequencing reveals dynamic mechanisms in aging and cancer

Telomeres are the protective nucleoprotein structures at the end of linear eukaryotic chromosomes. Telomeres' repetitive nature and length have traditionally challenged the precise assessment of the composition and length of individual human telomeres. Here, we present Telo-seq to resolve bulk, chromosome arm-specific and allele-specific human telomere lengths using Oxford Nanopore Technologies' native long-read sequencing. Telo-seq resolves telomere shortening in five population doubling increments and reveals intrasample, chromosome arm-specific, allele-specific telomere length heterogeneity. Telo-seq can reliably discriminate between telomerase- and ALT-positive cancer cell lines. Thus, Telo-seq is a tool to study telomere biology during development, aging, and cancer at unprecedented resolution.

Telomere length and cancer risk: finding Goldilocks

Telomeres are the nucleoprotein complex at chromosome ends essential in genomic stability. Baseline telomere length (TL) is determined by rare and common germline genetic variants but shortens with age and is susceptible to certain environmental exposures. Cellular senescence or apoptosis are normally triggered when telomeres reach a critically short length, but cancer cells overcome these protective mechanisms and continue to divide despite chromosomal instability. Rare germline variants in telomere maintenance genes cause exceedingly short telomeres for age (< 1st percentile) and the telomere biology disorders, which are associated with elevated risks of bone marrow failure, myelodysplastic syndrome, acute myeloid leukemia, and squamous cell carcinoma of the head/neck and anogenital regions. Long telomeres due to rare germline variants in the same or different telomere maintenance genes are associated with elevated risks of other cancers, such as chronic lymphocytic leukemia or sarcoma. Early epidemiology studies of TL in the general population lacked reproducibility but new methods, including creation of a TL polygenic score using common variants, have found longer telomeres associated with excess risks of renal cell carcinoma, glioma, lung cancer, and others. It has become clear that when it comes to TL and cancer etiology, not too short, not too long, but "just right" telomeres are important in minimizing cancer risk.

Multiomics-Based Feature Extraction and Selection for the Prediction of Lung Cancer Survival

Lung cancer is a global health challenge, hindered by delayed diagnosis and the disease's complex molecular landscape. Accurate patient survival prediction is critical, motivating the exploration of various -omics datasets using machine learning methods. Leveraging multi-omics data, this study seeks to enhance the accuracy of survival prediction by proposing new feature extraction techniques combined with unbiased feature selection. Two lung adenocarcinoma multi-omics datasets, originating from the TCGA and CPTAC-3 projects, were employed for this purpose, emphasizing gene expression, methylation, and mutations as the most relevant data sources that provide features for the survival prediction models. Additionally, gene set aggregation was shown to be the most effective feature extraction method for mutation and copy number variation data. Using the TCGA dataset, we identified 32 molecular features that allowed the construction of a 2-year survival prediction model with an AUC of 0.839. The selected features were additionally tested on an independent CPTAC-3 dataset, achieving an AUC of 0.815 in nested cross-validation, which confirmed the robustness of the identified features.

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.

Leveraging whole-genome sequencing to estimate telomere length in plants

Changes in telomere length are increasingly used to indicate species' response to environmental stress across diverse taxa. Despite this broad use, few studies have explored telomere length in plants. Thus, evaluation of new approaches for measuring telomeres in plants is needed. Rapid advances in sequencing approaches and bioinformatic tools now allow estimation of telomere content from whole-genome sequencing (WGS) data, a proxy for telomere length. While telomere content has been quantified extensively using quantitative polymerase chain reaction (qPCR) and WGS in humans, no study to date has compared the effectiveness of WGS in estimating telomere length in plants relative to qPCR approaches. In this study, we use 100 Populus clones re-sequenced using short-read Illumina sequencing to quantify telomere length comparing three different bioinformatic approaches (Computel, K-seek and TRIP) in addition to qPCR. Overall, telomere length estimates varied across different bioinformatic approaches, but were highly correlated across methods for individual genotypes. A positive correlation was observed between WGS estimates and qPCR, however, Computel estimates exhibited the greatest correlation. Computel incorporates genome coverage into telomere length calculations, suggesting that genome coverage is likely important to telomere length quantification when using WGS data. Overall, telomere estimates from WGS provided greater precision and accuracy of telomere length estimates relative to qPCR. The findings suggest WGS is a promising approach for assessing telomere length and, as the field of telomere ecology evolves, may provide added value to assaying response to biotic and abiotic environments for plants needed to accelerate plant breeding and conservation management.

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.

A Multibreed Genome-Wide Association Study for Cattle Leukocyte Telomere Length

Telomeres are terminal DNA regions of chromosomes that prevent chromosomal fusion and degradation during cell division. In cattle, leukocyte telomere length (LTL) is associated with longevity, productive lifespan, and disease susceptibility. However, the genetic basis of LTL in this species is less studied than in humans. In this study, we utilized the whole-genome resequencing data of 239 animals from 17 cattle breeds for computational leukocyte telomere length estimation and subsequent genome-wide association study of LTL. As a result, we identified 42 significant SNPs, of which eight were found in seven genes (EXOC6B, PTPRD, RPS6KC1, NSL1, AGBL1, ENSBTAG00000052188, and GPC1) when using covariates for two major breed groups (Turano-Mongolian and European). Association analysis with covariates for breed effect detected 63 SNPs, including 13 in five genes (EXOC6B, PTPRD, RPS6KC1, ENSBTAG00000040318, and NELL1). The PTPRD gene, demonstrating the top signal in analysis with breed effect, was previously associated with leukocyte telomere length in cattle and likely is involved in the mechanism of alternative lengthening of telomeres. The single nucleotide variants found could be tested for marker-assisted selection to improve telomere-length-associated traits.

Telomeres are shorter in wild Saccharomyces cerevisiae isolates than in domesticated ones

Telomeres are ribonucleoproteins that cap chromosome-ends and their DNA length is controlled by counteracting elongation and shortening processes. The budding yeast Saccharomyces cerevisiae has been a leading model to study telomere DNA length control and dynamics. Its telomeric DNA is maintained at a length that slightly varies between laboratory strains, but little is known about its variation at the species level. The recent publication of the genomes of over 1,000 S. cerevisiae strains enabled us to explore telomere DNA length variation at an unprecedented scale. Here, we developed a bioinformatic pipeline (YeaISTY) to estimate telomere DNA length from whole-genome sequences and applied it to the sequenced S. cerevisiae collection. Our results revealed broad natural telomere DNA length variation among the isolates. Notably, telomere DNA length is shorter in those derived from wild rather than domesticated environments. Moreover, telomere DNA length variation is associated with mitochondrial metabolism, and this association is driven by wild strains. Overall, these findings reveal broad variation in budding yeast's telomere DNA length regulation, which might be shaped by its different ecological life-styles.

Telomere length and stroke recurrence after ischemic stroke and TIA

BACKGROUND AND OBJECTIVE

Shortening telomere length (TL), as an indicator of aging, has been associated with increased risk of cardiovascular disease and incident stroke. However, there are limited data relating to the association between TL and recurrent stroke.

METHODS

Patients from the Third China National Stroke Registry who had whole genome sequencing (WGS) were selected. TL was estimated using TelSeq based on binary sequence alignment/map files derived from WGS data. Cox proportional hazards regression models were performed to assess the association of TL with recurrent stroke.

RESULTS

A total of 8041 patients with ischemic stroke (IS) or transient ischemic attack (TIA) were included. Mean TL was 2.14 ± 0.82 kb. Patients in the lowest tertile of TL had higher incidence of stroke recurrence compared to those in the middle and highest tertile (6.4% vs 5.9% vs 5.2%), but the difference was not longer significant after adjusting for age, sex, cardiovascular risk factors and stroke severity. Similarly, when analyzing TL as a continuous variable, the HR per 1000 bp increase in TL was significant 0.88 (0.79-0.98), but after adjusting for co-variates, was no longer significant (0.91; 95% confidence interval (CI), 0.81-1.02). In patients aged > 65 years, but not in younger patients, after adjusting for co-variates, TL was significantly associated with stroke recurrence. Compared to the lowest tertile, HRs (95% CI) after adjustment for all co-variates for the middle and highest tertiles were 0.78 (0.55-1.10) and 0.67 (0.46-0.98), respectively, with p for trend of 0.03. In analyses using TL as a continuous variable, adjusted HR (95% CI) per 1000 bp increase in TL was 0.80 (0.66-0.96). However, there was no significant interaction between TL and age on risk of stroke recurrence (p for interaction = 0.09).

CONCLUSIONS

In Chinese IS or TIA patients, no independent association was found between TL and risk of stroke recurrence after adjusting for co-variates. We found a possible association in older patients but this needs replicating.

Germline TP53 mutations undergo copy number gain years prior to tumor diagnosis

Li-Fraumeni syndrome (LFS) is a hereditary cancer predisposition syndrome associated with germline TP53 pathogenic variants. Here, we perform whole-genome sequence (WGS) analysis of tumors from 22 patients with TP53 germline pathogenic variants. We observe somatic mutations affecting Wnt, PI3K/AKT signaling, epigenetic modifiers and homologous recombination genes as well as mutational signatures associated with prior chemotherapy. We identify near-ubiquitous early loss of heterozygosity of TP53, with gain of the mutant allele. This occurs earlier in these tumors compared to tumors with somatic TP53 mutations, suggesting the timing of this mark may distinguish germline from somatic TP53 mutations. Phylogenetic trees of tumor evolution, reconstructed from bulk and multi-region WGS, reveal that LFS tumors exhibit comparatively limited heterogeneity. Overall, our study delineates early copy number gains of mutant TP53 as a characteristic mutational process in LFS tumorigenesis, likely arising years prior to tumor diagnosis.

Identification of the Sequence and the Length of Telomere DNA

Telomeres are essential nucleoprotein structures at the very ends of linear eukaryote chromosomes. They shelter the terminal genome territories against degradation and prevent the natural chromosome ends from being recognized by repair mechanisms as double-strand DNA breaks.There are two basic characteristics of telomeric DNA, its sequence and its length. The telomere sequence is important as a "landing area" for specific telomere-binding proteins, which function as signals and moderate the interactions required for correct telomere function. While the sequence forms the proper "landing surface" of telomeric DNA, its length is similarly important. Too short or exceptionally long telomere DNA cannot perform its function properly. In this chapter, methods for the investigation of these two basic telomere DNA characteristics are described, namely, telomere motif identification and telomere length measurement.

Telomere Length: Implications for Atherogenesis

PURPOSE OF REVIEW

The purpose of the study is to explore the evidence linking telomere length with atherosclerotic ischemic disease.

RECENT FINDINGS

There has been a recent expansion in strategies for measuring telomere length, including analyzing genome sequence data and capitalizing on genomic loci that associate with telomere length. These, together with more established approaches, have been used to generate a more complete picture of telomere length relationships with ischemic disease. Whereas earlier meta-analyses suggested an association between short leukocyte telomeres and ischemic disease, several recent large population studies now provide particularly compelling data, including an association with cardiovascular mortality. In addition, whether short leukocyte telomeres might be causally related to ischemic disease has been interrogated using Mendelian randomization strategies, which point to shorter leukocyte telomeres as a determining risk factor. Importantly however, the wide, interindividual variability in telomere length still means that a single assessment of leukocyte telomere length in an individual does not reliably report on a biological aging process. In this regard, recent multi-tissue analyses of telomere length dynamics are providing both new mechanistic insights into how telomere length and shortening rates may participate in atherogenesis and risk prediction opportunities. The balance of evidence indicates that short leukocyte telomeres confer a risk for atherosclerotic cardiovascular disease. Moreover, an integrated analysis of telomere lengths in leukocytes and other tissues may provide a window into individualized telomere dynamics, raising new prospects for risk management.

Taking flight: An ecological, evolutionary and genomic perspective on bat telomeres

Over 20% of all living mammals are bats (order Chiroptera). Bats possess extraordinary adaptations including powered flight, laryngeal echolocation and a unique immune system that enables them to tolerate a diversity of viral infections without presenting clinical disease symptoms. They occupy multiple trophic niches and environments globally. Significant physiological and ecological diversity occurs across the order. Bats also exhibit extreme longevity given their body size with many species showing few signs of ageing. The molecular basis of this extended longevity has recently attracted attention. Telomere maintenance potentially underpins bats' extended healthspan, although functional studies are still required to validate the causative mechanisms. In this review, we detail the current knowledge on bat telomeres, telomerase expression, and how these relate to ecology, longevity and life-history strategies. Patterns of telomere shortening and telomerase expression vary across species, and comparative genomic analyses suggest that alternative telomere maintenance mechanisms evolved in the longest-lived bats. We discuss the unique challenges faced when working with populations of wild bats and highlight ways to advance the field including expanding long-term monitoring across species that display contrasting life-histories and occupy different environmental niches. We further review how new high quality, chromosome-level genome assemblies can enable us to uncover the molecular mechanisms governing telomere dynamics and how phylogenomic analyses can reveal the adaptive significance of telomere maintenance and variation in bats.

Leukocyte telomere length in children born following blastocyst-stage embryo transfer

Perinatal and childhood adverse outcomes associated with assisted reproductive technology (ART) has been reported, but it remains unknown whether the initial leukocyte telomere length (LTL), which is an indicator of age-related phenotypes in later life, is affected. Here, we estimated the LTLs of 1,137 individuals from 365 families, including 202 children conceived by ART and 205 children conceived spontaneously from two centers of the China National Birth Cohort, using whole-genome sequencing (WGS) data. One-year-old children conceived by ART had shorter LTLs than those conceived spontaneously (beta, -0.36; P = 1.29 × 10^-3) after adjusting for plurality, sex and other potential confounding factors. In particular, blastocyst-stage embryo transfer was associated with shorter LTL (beta, -0.54, P = 2.69 × 10^-3) in children conceived by ART. The association was validated in 586 children conceived by ART from five centers using different LTL quantification methods (that is, WGS or qPCR). Blastocyst-stage embryo transfer resulted in shorter telomere lengths in mice at postnatal day 1 (P = 2.10 × 10^-4) and mice at 6 months (P = 0.042). In vitro culturing of mice embryos did not result in shorter telomere lengths in the late cleavage stage, but it did suppress telomerase activity in the early blastocyst stage. Our findings demonstrate the need to evaluate the long-term consequences of ART, particularly for aging-related phenotypes, in children conceived by ART.

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.

Telomere Length as a marker of biological aging: A critical review of recent literature

INTRODUCTION

Aging is characterized as a syndrome of deleterious, progressive, universal, and irreversible function changes affecting every structural and functional aspect of the organism and accompanied by a generalized increase in mortality. Although a substantial number of candidates for biomarkers of aging have been proposed, none has been validated or universally accepted. Human telomeres constitute hexameric repetitive DNA sequence nucleoprotein complexes that cap chromosome ends, regulating gene expression and modulating stress-related pathways. Telomere length (TL) shortening is observed both in cellular senescence and advanced age, leading to the investigation of TL as a biomarker for aging and a risk factor indicator for the development and progression of the most common age-related diseases.

OBJECTIVE

The present review underlines the connection between TL and the pathophysiology of the diseases associated with telomere attrition.

METHODS

We performed a structured search of the PubMed database for peer-reviewed research of the literature regarding leukocyte TL and cardiovascular diseases (CVD), more specifically stroke and heart disease, and focused on the relevant articles published during the last 5 years. We also applied Hill's criteria of causation to strengthen this association.

RESULTS

We analyzed the recent literature regarding TL length, stroke, and CVD. Although approximately one-third of the available studies support the connection, the results of different studies seem to be rather conflicting as a result of different study designs, divergent methods of TL determination, small study samples, and patient population heterogeneity. After applying Hill's criteria, we can observe that the literature conforms to them weakly, with chronology being the only Hill criterion of causality that probably cannot be contested.

CONCLUSION

The present review attempted to examine the purported relation between leukocyte TL and age-related diseases such as CVD and more specific stroke and heart disease in view of the best established, comprehensive, medical and epidemiological criteria that have characterized the focused recent relevant research. Although several recommendations have been made that may contribute significantly to the field, a call for novel technical approaches and studies is mandatory to further elucidate the possible association.

Telomere length and brain aging: A systematic review and meta-analysis

The current evidence on the association of leukocyte telomere length (LTL) with age-related structural and cognitive changes in the brain is mixed. Herein conforming to PRISMA 2020 guidelines, we performed a systematic review and meta-analysis using data from 27 observational studies in non-demented individuals. We used effect size and p-value based meta-analysis methods considering marked heterogeneity among studies. We found that the longer LTL was associated with higher brain volume (β = 0.43, 95%CI: 0.36-0.50%, p = 0.008, N = 1102) and with higher global cognition (β = 0.01; 95%CI: 0.00-0.02, p = 0.03, N = 19609) by effect size based meta-analysis and with brain volume, hippocampal volume, global cognition, cognitive domains of attention/speed as well as executive functions by p-value based meta-analysis. No significant association of LTL with brain white matter hyperintensities was detected. Furthermore, the evidence strongly suggests a subgroup-specific canonical effect of telomeres, notably in older individuals and females. In conclusion, we provide meta-analytic evidence on the beneficial effect of telomeres on brain structure as well as cognition and advocate for a beneficial subgroup-specific effect that warrants further attention.

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.

Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed

Genetic studies on telomere length are important for understanding age-related diseases. Prior GWAS for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally-diverse individuals (European, African, Asian and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n=109,122 individuals. We identified 59 sentinel variants (p-value <5×10-9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated the independent signals colocalized with cell-type specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated our TL polygenic trait scores (PTS) were associated with increased risk of cancer-related phenotypes.

WALTER: an easy way to online evaluate telomere lengths from terminal restriction fragment analysis

BACKGROUND

Telomeres, nucleoprotein structures comprising short tandem repeats and delimiting the ends of linear eukaryotic chromosomes, play an important role in the maintenance of genome stability. Therefore, the determination of the length of telomeres is of high importance for many studies. Over the last years, new methods for the analysis of the length of telomeres have been developed, including those based on PCR or analysis of NGS data. Despite that, terminal restriction fragment (TRF) method remains the gold standard to this day. However, this method lacks universally accepted and precise tool capable to analyse and statistically evaluate TRF results.

RESULTS

To standardize the processing of TRF results, we have developed WALTER, an online toolset allowing rapid, reproducible, and user-friendly analysis including statistical evaluation of the data. Given its web-based nature, it provides an easily accessible way to analyse TRF data without any need to install additional software.

CONCLUSIONS

WALTER represents a major upgrade from currently available tools for the image processing of TRF scans. This toolset enables a rapid, highly reproducible, and user-friendly evaluation of almost any TRF scan including in-house statistical evaluation of the data. WALTER platform together with user manual describing the evaluation of TRF scans in detail and presenting tips and troubleshooting, as well as test data to demo the software are available at https://www.ceitec.eu/chromatin-molecular-complexes-jiri-fajkus/rg51/tab?tabId=125#WALTER and the source code at https://github.com/mlyc93/WALTER .

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.

Temporal Telomere and DNA Damage Responses in the Space Radiation Environment

Telomeres, repetitive terminal features of chromosomes essential for maintaining genome integrity, shorten with cell division, lifestyle factors and stresses, and environmental exposures, and so they provide a robust biomarker of health, aging, and age-related diseases. We assessed telomere length dynamics (changes over time) in three unrelated astronauts before, during, and after 1-year or 6-month missions aboard the International Space Station (ISS). Similar to our results for National Aeronautics and Space Administration's (NASA's) One-Year Mission twin astronaut (Garrett-Bakelman et al., 2019), significantly longer telomeres were observed during spaceflight for two 6-month mission astronauts. Furthermore, telomere length shortened rapidly after return to Earth for all three crewmembers and, overall, telomere length tended to be shorter after spaceflight than before spaceflight. Consistent with chronic exposure to the space radiation environment, signatures of persistent DNA damage responses were also detected, including mitochondrial and oxidative stress, inflammation, and telomeric and chromosomal aberrations, which together provide potential mechanistic insight into spaceflight-specific telomere elongation.

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.

Telomere Length Maintenance and Its Transcriptional Regulation in Lynch Syndrome and Sporadic Colorectal Carcinoma

Background: Activation of telomere maintenance mechanisms (TMMs) is a hallmark of most cancers, and is required to prevent genome instability and to establish cellular immortality through reconstitution of capping of chromosome ends. TMM depends on the cancer type. Comparative studies linking tumor biology and TMM have potential impact for evaluating cancer onset and development. Methods: We have studied alterations of telomere length, their sequence composition and transcriptional regulation in mismatch repair deficient colorectal cancers arising in Lynch syndrome (LS-CRC) and microsatellite instable (MSI) sporadic CRC (MSI s-CRC), and for comparison, in microsatellite stable (MSS) s-CRC and in benign colon mucosa. Our study applied bioinformatics analysis of whole genome DNA and RNA sequencing data and a pathway model to study telomere length alterations and the potential effect of the "classical" telomerase (TEL-) and alternative (ALT-) TMM using transcriptomic signatures. Results: We have found progressive decrease of mean telomere length in all cancer subtypes compared with reference systems. Our results support the view that telomere attrition is an early event in tumorigenesis. TMM gets activated in all tumors studied due to concerted overexpression of a large fraction of genes with direct relation to telomere function, where only a very small fraction of them showed recurrent mutations. TEL-related transcriptional state was dominating in all CRC subtypes, showing, however, subtype-specific activation patterns; while contribution of the ALT-TMM was slightly more prominent in the hypermutated MSI s-CRC and LS-CRC. TEL-TMM is mainly activated by over-expression of DKC1 and/or TERT genes and their interaction partners, where DKC1 is more prominent in MSS than in MSI s-CRC and can serve as a transcriptomic marker of TMM activity. Conclusions: Our results suggest that transcriptional patterns are indicative for TMM pathway activation with subtle differences between TEL and ALT mechanisms in a CRC subtype-specific fashion. Sequencing data potentially provide a suited measure to study alterations of telomere length and of underlying transcriptional regulation. Further studies are needed to improve this method.

A bioinformatics approach to identify telomere sequences

Conventional approaches to identify a telomere motif in a new genome are laborious and time-intensive. An efficient new methodology based on next-generation sequencing (NGS), de novo sequence repeat finder (SERF) and fluorescence in situ hybridization (FISH) is presented. Unlike existing heuristic approaches, SERF utilizes an exhaustive analysis of raw NGS reads or assembled contigs for rapid de novo detection of conserved tandem repeats representing telomere motifs. SERF was validated using the NGS data from Ipheion uniflorum and Allium cepa with known telomere motifs. The analysis program was then used on NGS data to investigate the telomere motifs in several additional plant species and together with FISH proved to be an efficient approach to identify as yet unknown telomere motifs.

Telomere sequence content can be used to determine ALT activity in tumours

The replicative immortality of human cancer cells is achieved by activation of a telomere maintenance mechanism (TMM). To achieve this, cancer cells utilise either the enzyme telomerase, or the Alternative Lengthening of Telomeres (ALT) pathway. These distinct molecular pathways are incompletely understood with respect to activation and propagation, as well as their associations with clinical outcomes. We have identified significant differences in the telomere repeat composition of tumours that use ALT compared to tumours that do not. We then employed a machine learning approach to stratify tumours according to telomere repeat content with an accuracy of 91.6%. Importantly, this classification approach is applicable across all tumour types. Analysis of pathway mutations that were under-represented in ALT tumours, across 1,075 tumour samples, revealed that the autophagy, cell cycle control of chromosomal replication, and transcriptional regulatory network in embryonic stem cells pathways are involved in the survival of ALT tumours. Overall, our approach demonstrates that telomere sequence content can be used to stratify ALT activity in cancers, and begin to define the molecular pathways involved in ALT activation.

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

BACKGROUND

Establishment of telomere maintenance mechanisms is a universal step in tumor development to achieve replicative immortality. These processes leave molecular footprints in cancer genomes in the form of altered telomere content and aberrations in telomere composition. To retrieve these telomere characteristics from high-throughput sequencing data the available computational approaches need to be extended and optimized to fully exploit the information provided by large scale cancer genome data sets.

RESULTS

We here present TelomereHunter, a software for the detailed characterization of telomere maintenance mechanism footprints in the genome. The tool is implemented for the analysis of large cancer genome cohorts and provides a variety of diagnostic diagrams as well as machine-readable output for subsequent analysis. A novel key feature is the extraction of singleton telomere variant repeats, which improves the identification and subclassification of the alternative lengthening of telomeres phenotype. We find that whole genome sequencing-derived telomere content estimates strongly correlate with telomere qPCR measurements (r = 0.94). For the first time, we determine the correlation of in silico telomere content quantification from whole genome sequencing and whole genome bisulfite sequencing data derived from the same tumor sample (r = 0.78). An analogous comparison of whole exome sequencing data and whole genome sequencing data measured slightly lower correlation (r = 0.79). However, this is considerably improved by normalization with matched controls (r = 0.91).

CONCLUSIONS

TelomereHunter provides new functionality for the analysis of the footprints of telomere maintenance mechanisms in cancer genomes. Besides whole genome sequencing, whole exome sequencing and whole genome bisulfite sequencing are suited for in silico telomere content quantification, especially if matched control samples are available. The software runs under a GPL license and is available at https://www.dkfz.de/en/applied-bioinformatics/telomerehunter/telomerehunter.html .

The telomere world and aging: Analytical challenges and future perspectives

Telomeres, the terminal nucleoprotein structures of eukaryotic chromosomes, play pleiotropic functions in cellular and organismal aging. Telomere length (TL) varies throughout life due to the influence of genetic factors and to a complex balancing between "shortening" and "elongation" signals. Telomerase, the only enzyme that can elongate a telomeric DNA chain, and telomeric repeat-containing RNA (TERRA), a long non-coding RNA involved in looping maintenance, play key roles in TL during life. Despite recent advances in the knowledge of TL, TERRA and telomerase activity (TA) biology and their measurement techniques, the experimental and theoretical issues involved raise a number of problems that should carefully be considered by researchers approaching the "telomere world". The increasing use of such parameters - hailed as promising clinically relevant biomarkers - has failed to be paralleled by the development of automated and standardized measurement technology. Consequently, associating given TL values to specific pathological conditions involves on the one hand technological issues and on the other clinical-biological issues related to the planning of clinically relevant association studies. Addressing these issues would help avoid major biases in association studies involving TL and a number of outcomes, especially those focusing on psychological and bio-behavioral variables. The main challenge in telomere research is the development of accurate and reliable measurement methods to achieve simple and sensitive TL, TERRA, and TA detection. The discovery of the localization of telomeres and TERRA in cellular and extracellular compartments had added an additional layer of complexity to the measurement of these age-related biomarkers. Since combined analysis of TL, TERRA and TA may well provide more exhaustive clinical information than a single parameter, we feel it is important for researchers in the various fields to become familiar with their most common measurement techniques and to be aware of the respective merits and drawbacks of these approaches.

Copy number signatures and mutational processes in ovarian carcinoma

The genomic complexity of profound copy number aberrations has prevented effective molecular stratification of ovarian cancers. Here, to decode this complexity, we derived copy number signatures from shallow whole-genome sequencing of 117 high-grade serous ovarian cancer (HGSOC) cases, which were validated on 527 independent cases. We show that HGSOC comprises a continuum of genomes shaped by multiple mutational processes that result in known patterns of genomic aberration. Copy number signature exposures at diagnosis predict both overall survival and the probability of platinum-resistant relapse. Measurement of signature exposures provides a rational framework to choose combination treatments that target multiple mutational processes.

Copy number signatures and mutational processes in ovarian carcinoma

The genomic complexity of profound copy number aberrations has prevented effective molecular stratification of ovarian cancers. Here, to decode this complexity, we derived copy number signatures from shallow whole-genome sequencing of 117 high-grade serous ovarian cancer (HGSOC) cases, which were validated on 527 independent cases. We show that HGSOC comprises a continuum of genomes shaped by multiple mutational processes that result in known patterns of genomic aberration. Copy number signature exposures at diagnosis predict both overall survival and the probability of platinum-resistant relapse. Measurement of signature exposures provides a rational framework to choose combination treatments that target multiple mutational processes.

Copy number signatures and mutational processes in ovarian carcinoma

The genomic complexity of profound copy number aberrations has prevented effective molecular stratification of ovarian cancers. Here, to decode this complexity, we derived copy number signatures from shallow whole-genome sequencing of 117 high-grade serous ovarian cancer (HGSOC) cases, which were validated on 527 independent cases. We show that HGSOC comprises a continuum of genomes shaped by multiple mutational processes that result in known patterns of genomic aberration. Copy number signature exposures at diagnosis predict both overall survival and the probability of platinum-resistant relapse. Measurement of signature exposures provides a rational framework to choose combination treatments that target multiple mutational processes.

Stress, Telomeres, and Psychopathology: Toward a Deeper Understanding of a Triad of Early Aging

Telomeres play an important part in aging and show relationships to lifetime adversity, particularly childhood adversity. Meta-analyses demonstrate reliable associations between psychopathology (primarily depression) and shorter telomere length, but the nature of this relationship has not been fully understood. Here, we review and evaluate the evidence for impaired telomere biology as a consequence of psychopathology or as a contributing factor, and the important mediating roles of chronic psychological stress and impaired allostasis. There is evidence for a triadic relationship among stress, telomere shortening, and psychiatric disorders that is positively reinforcing and unfolds across the life course and, possibly, across generations. We review the role of genetics and biobehavioral responses that may contribute to shorter telomere length, as well as the neurobiological impact of impaired levels of telomerase. These complex interrelationships are important to elucidate because they have implications for mental and physical comorbidity and, potentially, for the prevention and treatment of depression.

Direct Quantitative Monitoring of Homology-Directed DNA Repair of Damaged Telomeres

Homology-directed DNA repair (HDR) is an evolutionary conserved mechanism that is required for genome integrity and organismal fitness across species. While a myriad of different factors and mechanisms are able to execute HDR, all forms necessitate common steps of DNA damage recognition, homology search and capture, and assembly of a DNA polymerase complex to conduct templated DNA synthesis. The central question of what determines HDR mechanism utilization in mammalian cells has been limited by an inability to directly monitor the DNA damage response and products of repair as they arise from a defined genomic lesion. In this chapter, we describe several methodologies to delineate major steps of HDR during alternative lengthening of telomeres in human cells. This includes procedures to visualize interchromosomal telomere homology searches in real time and quantitatively detect HDR synthesis of nascent telomeres emanating from synchronous activation of telomere DNA double-strand breaks. We highlight the critical details of these methods and their applicability to monitoring HDR at telomeres in a broad variety of mammalian cell types.

Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data

Telomere length is a risk factor in disease and the dynamics of telomere length are crucial to our understanding of cell replication and vitality. The proliferation of whole genome sequencing represents an unprecedented opportunity to glean new insights into telomere biology on a previously unimaginable scale. To this end, a number of approaches for estimating telomere length from whole-genome sequencing data have been proposed. Here we present Telomerecat, a novel approach to the estimation of telomere length. Previous methods have been dependent on the number of telomeres present in a cell being known, which may be problematic when analysing aneuploid cancer data and non-human samples. Telomerecat is designed to be agnostic to the number of telomeres present, making it suited for the purpose of estimating telomere length in cancer studies. Telomerecat also accounts for interstitial telomeric reads and presents a novel approach to dealing with sequencing errors. We show that Telomerecat performs well at telomere length estimation when compared to leading experimental and computational methods. Furthermore, we show that it detects expected patterns in longitudinal data, repeated measurements, and cross-species comparisons. We also apply the method to a cancer cell data, uncovering an interesting relationship with the underlying telomerase genotype.

TERT structural rearrangements in metastatic pheochromocytomas

Pheochromocytomas (PC) and paragangliomas (PGL) are endocrine tumors for which the genetic and clinicopathological features of metastatic progression remain incompletely understood. As a result, the risk of metastasis from a primary tumor cannot be predicted. Early diagnosis of individuals at high risk of developing metastases is clinically important and the identification of new biomarkers that are predictive of metastatic potential is of high value. Activation of TERT has been associated with a number of malignant tumors, including PC/PGL. However, the mechanism of TERT activation in the majority of PC/PGL remains unclear. As TERT promoter mutations occur rarely in PC/PGL, we hypothesized that other mechanisms - such as structural variations - may underlie TERT activation in these tumors. From 35 PC and four PGL, we identified three primary PCs that developed metastases with elevated TERT expression, each of which lacked TERT promoter mutations and promoter DNA methylation. Using whole genome sequencing, we identified somatic structural alterations proximal to the TERT locus in two of these tumors. In both tumors, the genomic rearrangements led to the positioning of super-enhancers proximal to the TERT promoter, that are likely responsible for the activation of the normally tightly repressed TERT expression in chromaffin cells.

A method for measuring the distribution of the shortest telomeres in cells and tissues

Improved methods to measure the shortest (not just average) telomere lengths (TLs) are needed. We developed Telomere Shortest Length Assay (TeSLA), a technique that detects telomeres from all chromosome ends from <1 kb to 18 kb using small amounts of input DNA. TeSLA improves the specificity and efficiency of TL measurements that is facilitated by user friendly image-processing software to automatically detect and annotate band sizes, calculate average TL, as well as the percent of the shortest telomeres. Compared with other TL measurement methods, TeSLA provides more information about the shortest telomeres. The length of telomeres was measured longitudinally in peripheral blood mononuclear cells during human aging, in tissues during colon cancer progression, in telomere-related diseases such as idiopathic pulmonary fibrosis, as well as in mice and other organisms. The results indicate that TeSLA is a robust method that provides a better understanding of the shortest length of telomeres.

Short telomeres are a hallmark of senescence and can result in genomic instability as well as cancer progression. Here, the authors present TeSLA, a technique to accurately detect telomeres under 1 kb in length.