Genetic polymorphisms associated with telomere length and risk of developing myeloproliferative neoplasms

Genetics and epidemiology of mutational barcode-defined clonal hematopoiesis

Clonal hematopoiesis (CH) arises when a substantial proportion of mature blood cells is derived from a single hematopoietic stem cell lineage. Using whole-genome sequencing of 45,510 Icelandic and 130,709 UK Biobank participants combined with a mutational barcode method, we identified 16,306 people with CH. Prevalence approaches 50% in elderly participants. Smoking demonstrates a dosage-dependent impact on risk of CH. CH associates with several smoking-related diseases. Contrary to published claims, we find no evidence that CH is associated with cardiovascular disease. We provide evidence that CH is driven by genes that are commonly mutated in myeloid neoplasia and implicate several new driver genes. The presence and nature of a driver mutation alters the risk profile for hematological disorders. Nevertheless, most CH cases have no known driver mutations. A CH genome-wide association study identified 25 loci, including 19 not implicated previously in CH. Splicing, protein and expression quantitative trait loci were identified for CD164 and TCL1A.

Parental telomeres implications on immune senescence of newborns

Telomere, the biological chronometer, and its effect on the immune system considerably varies among individuals. During pregnancy, multiple risk factors affect telomere reprogramming during fetal life which can lead to health disparities in newborns. These changes may cause a long-term impact on the telomere genetics of the newborn and become a reason for lifelong health implications and immune senescence. Therefore, telomere shortening in parents due to genetic variation may act as a hallmark of immune senescence and aging in their newborns.

Predicted leukocyte telomere length and risk of myeloid neoplasms

Maintenance of telomere length has long been established to play a role in the biology of cancer and several studies suggest that it may be especially important in myeloid malignancies. To overcome potential bias in confounding and reverse causation of observational studies, we use both a polygenic risk score (PRS) and inverse-variance weighted (IVW) Mendelian randomization (MR) analyses to estimate the relationship between genetically predicted leukocyte telomere length (LTL) and acute myeloid leukemia (AML) risk in 498 cases and 2099 controls and myelodysplastic syndrome (MDS) risk in 610 cases and 1759 controls. Genetic instruments derived from four recent studies explaining 1.23-4.57% of telomere variability were considered. We used multivariable logistic regression to estimate odds ratios (OR, 95% confidence intervals [CI]) as the measure of association between individual single-nucleotide polymorphisms and myeloid malignancies. We observed a significant association between a PRS of longer predicted LTL and AML using three genetic instruments (OR = 4.03 per ~1200 base pair [bp] increase in LTL, 95% CI: 1.65, 9.85 using Codd et al. [Codd, V., Nelson, C.P., Albrecht, E., Mangino, M., Deelen, J., Buxton, J.L., Hottenga, J.J., Fischer, K., Esko, T., Surakka, I. et al. (2013) Identification of seven loci affecting mean telomere length and their association with disease. Nat. Genet., 45, 422-427 427e421-422.], OR = 3.48 per one-standard deviation increase in LTL, 95% CI: 1.74, 6.97 using Li et al. [Li, C., Stoma, S., Lotta, L.A., Warner, S., Albrecht, E., Allione, A., Arp, P.P., Broer, L., Buxton, J.L., Alves, A.D.S.C. et al. (2020) Genome-wide association analysis in humans links nucleotide metabolism to leukocyte telomere length. Am. J. Hum. Genet., 106, 389-404.] and OR = 2.59 per 1000 bp increase in LTL, 95% CI: 1.03, 6.52 using Taub et al. [Taub, M.A., Conomos, M.P., Keener, R., Iyer, K.R., Weinstock, J.S., Yanek, L.R., Lane, J., Miller-Fleming, T.W., Brody, J.A., Raffield, L.M. et al. (2022) Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed. Cell Genom., 2.] genetic instruments). MR analyses further indicated an association between LTL and AML risk (PIVW ≤ 0.049) but not MDS (all PIVW ≥ 0.076). Findings suggest variation in genes relevant to telomere function and maintenance may be important in the etiology of AML but not MDS.

A pleiotropic variant in DNAJB4 is associated with multiple myeloma risk

Pleiotropy, which consists of a single gene or allelic variant affecting multiple unrelated traits, is common across cancers, with evidence for genome-wide significant loci shared across cancer and noncancer traits. This feature is particularly relevant in multiple myeloma (MM) because several susceptibility loci that have been identified to date are pleiotropic. Therefore, the aim of this study was to identify novel pleiotropic variants involved in MM risk using 28 684 independent single nucleotide polymorphisms (SNPs) from GWAS Catalog that reached a significant association (P < 5 × 10^-8 ) with their respective trait. The selected SNPs were analyzed in 2434 MM cases and 3446 controls from the International Lymphoma Epidemiology Consortium (InterLymph). The 10 SNPs showing the strongest associations with MM risk in InterLymph were selected for replication in an independent set of 1955 MM cases and 1549 controls from the International Multiple Myeloma rESEarch (IMMEnSE) consortium and 418 MM cases and 147 282 controls from the FinnGen project. The combined analysis of the three studies identified an association between DNAJB4-rs34517439-A and an increased risk of developing MM (OR = 1.22, 95%CI 1.13-1.32, P = 4.81 × 10^-7 ). rs34517439-A is associated with a modified expression of the FUBP1 gene, which encodes a multifunctional DNA and RNA-binding protein that it was observed to influence the regulation of various genes involved in cell cycle regulation, among which various oncogenes and oncosuppressors. In conclusion, with a pleiotropic scan approach we identified DNAJB4-rs34517439 as a potentially novel MM risk locus.

Risk factors for clonal hematopoiesis of indeterminate potential and mosaic chromosomal alterations

Clonal hematopoiesis of indeterminate potential (CHIP) and mosaic chromosomal alterations (mCAs) of the autosomes, X, and Y chromosomes are aging-related somatic mutations detectable in peripheral blood. The presence of these acquired mutations predisposes otherwise healthy adults to increased risk of several chronic aging-related conditions including hematologic cancers, atherosclerotic cardiovascular diseases, other inflammatory conditions, and mortality. While the public health impact and disease associations of these blood-derived somatic mutations continue to expand, the inherited, behavioral/lifestyle, environmental risk factors and comorbid conditions that influence their occurrence and progression have been less well characterized. Age is the strongest risk factor for all types of CHIP and mCAs. CHIP and mCAs are generally more common in individuals of European than non-European ancestry. Evidence for a genetic predisposition has been strongest for mosaic loss of Y chromosome in men. Genome-wide association studies have recently begun to identify common and rare germline genetic variants associated with CHIP and mCAs. These loci include genes involving cell cycle regulation, cell proliferation/survival, hematopoietic progenitor cell regulation, DNA damage repair, and telomere maintenance. Some loci, such as TERT, ATM, TP53, CHEK2, and TCL1A, have overlapping associations with different types of CHIP, mCAs, and cancer predisposition. Various environmental or co-morbid contexts associated with presence or expansion of specific CHIP or mCA mutations are beginning to be elucidated, such as cigarette smoking, diet, cancer chemotherapy, particulate matter, and premature menopause. Further characterization of the germline genetic and environmental correlates of CHIP/mCAs may inform our ability to modify their progression and ultimately reduce the risk and burden of chronic diseases associated with these clonal somatic phenomena.

Genetic Variants of the TERT Gene and Telomere Length in Obstructive Sleep Apnea

Introduction: Obstructive sleep apnea (OSA) is a worldwide breathing disorder that has been diagnosed globally in almost 1 billion individuals aged 30−69 years. It is characterized by repeated upper airway collapses during sleep. Telomerase reverse transcriptase (TERT) is involved in the prevention of telomere shortening. This prospective, observational study aimed to investigate the relationship between single nucleotide polymorphisms (SNPs) of TERT and the severity of OSA, taking into account hypertension and diabetes prevalence. Methods: A total of 149 patients with OSA were diagnosed using one-night video-polysomnography based on the American Academy of Sleep Medicine guidelines. The TERT SNPs and telomere length (TL) were detected using real-time quantitative polymerase chain reaction. Results: Statistical analysis showed that there is no relationship between the rs2853669 and rs2736100 polymorphisms of TERT, and the severity of OSA (p > 0.05). Moreover, no relationship between TL and the severity of OSA was observed. The G allele in the locus of rs2736100 TERT was associated with hypertension prevalence and was more prevalent in hypertensives patients (46.00% vs. 24.49%, p = 0.011). The prevalence of hypertension was higher in patients with the C allele in the locus of rs2853669 than in patients without this allele (50.79% vs. 30.23%, p = 0.010). Moreover, a lower prevalence of diabetes was observed in homozygotes of rs2736100 TERT than in heterozygotes (5.63% vs. 15.38%, p = 0.039). Conclusion: This study showed no relationship between OSA and TERT SNPs. However, SNPs of the TERT gene (rs2736100 and rs2853669) were found to affect arterial hypertension and diabetes prevalence.

Transcriptomic Analysis of Conserved Telomere Maintenance Component 1 (CTC1) and Its Association with Leukemia

Telomere length (TEL) regulation is important for genome stability and is governed by the coordinated role of shelterin proteins, telomerase (TERT), and CST (CTC1/OBFC1/TEN1) complex. Previous studies have shown the association of telomerase expression with the risk of acute lymphoblastic leukemia (ALL). However, no data are available for CST association with the ALL. The current pilot study was designed to evaluate the CST expression levels in ALL. In total, 350 subjects were recruited, including 250 ALL cases and 100 controls. The subjects were stratified by age and categorized into pediatrics (1–18 years) and adults (19–54 years). TEL and expression patterns of CTC1, OBFC1, and TERT genes were determined by qPCR. The univariable logistic regression analysis was performed to determine the association of gene expression with ALL, and the results were adjusted for age and sex in multivariable analyses. Pediatric and adult cases did not reflect any change in telomere lengths relative to controls. However, expression of CTC1, OBFC1, and TERT genes were induced among ALL cases. Multivariable logistic regression analyses showed association of CTC1 with ALL in pediatric [β estimate (standard error (SE)= −0.013 (0.007), p = 0.049, and adults [0.053 (0.023), p = 0.025]. The association of CTC1 remained significant when taken together with OBFC1 and TERT in a multivariable model. Furthermore, CTC1 showed significant association with B-cell ALL [−0.057(0.017), p = 0.002) and T-cell ALL [−0.050 (0.018), p = 0.008] in pediatric group while no such association was noted in adults. Together, our findings demonstrated that telomere modulating genes, particularly CTC1, are strongly associated with ALL. Therefore, CTC1 can potentially be used as a risk biomarker for the identification of ALL in both pediatrics and adults.

A polymorphic variant in telomere maintenance is associated with worrisome features and high-risk stigmata development in IPMNs

Intraductal papillary mucinous neoplasms (IPMNs) are nonobligatory precursor lesions of pancreatic ductal adenocarcinoma (PDAC). The identification of molecular biomarkers able to predict the risk of progression of IPMNs toward malignancy is largely lacking and sorely needed. Telomere length (TL) is associated with the susceptibility of developing cancers, including PDAC. Moreover, several PDAC risk factors have been shown to be associated with IPMN transition to malignancy. TL is genetically determined, and the aim of this study was to use 11 SNPs, alone or combined in a score (teloscore), to estimate the causal relation between genetically determined TL and IPMNs progression. For this purpose, 173 IPMN patients under surveillance were investigated. The teloscore did not show any correlation, however, we observed an association between PXK-rs6772228-A and an increased risk of IPMN transition to malignancy (HR = 3.17; 95%CI 1.47-6.84; P = 3.24 × 10-3). This effect was also observed in a validation cohort of 142 IPMNs even though the association was not statistically significant. The combined analysis was consistent showing an association between PXK-rs6772228-A and increased risk of progression. The A allele of this SNP is strongly associated with shorter LTL that in turn have been reported to be associated with increased risk of developing PDAC. These results clearly highlight the importance of looking for genetic variants as potential biomarkers in this setting in order to further our understanding the etiopathogenesis of PDAC and suggest that genetically determined TL might be an additional marker of IPMN prognosis.

Germline-somatic JAK2 interactions are associated with clonal expansion in myelofibrosis

Myelofibrosis is a rare myeloproliferative neoplasm (MPN) with high risk for progression to acute myeloid leukemia. Our integrated genomic analysis of up to 933 myelofibrosis cases identifies 6 germline susceptibility loci, 4 of which overlap with previously identified MPN loci. Virtual karyotyping identifies high frequencies of mosaic chromosomal alterations (mCAs), with enrichment at myelofibrosis GWAS susceptibility loci and recurrently somatically mutated MPN genes (e.g., JAK2). We replicate prior MPN associations showing germline variation at the 9p24.1 risk haplotype confers elevated risk of acquiring JAK2^V617F mutations, demonstrating with long-read sequencing that this relationship occurs in cis. We also describe recurrent 9p24.1 large mCAs that selectively retained JAK2^V617F mutations. Germline variation associated with longer telomeres is associated with increased myelofibrosis risk. Myelofibrosis cases with high-frequency JAK2 mCAs have marked reductions in measured telomere length – suggesting a relationship between telomere biology and myelofibrosis clonal expansion. Our results advance understanding of the germline-somatic interaction at JAK2 and implicate mCAs involving JAK2 as strong promoters of clonal expansion of those mutated clones.

Myelofibrosis is a risk factor for the development of Acute Myeloid Leukaemia. Here, the authors carry out an integrated genomic investigation of 933 myelofibrosis patients, and identified interactions between germline and somatic variation in patients who required haematopoietic cell transplantation.

Mendelian randomization supports bidirectional causality between telomere length and clonal hematopoiesis of indeterminate potential

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program (n = 63,302) and UK Biobank (n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

Systematic review of Mendelian randomization studies on risk of cancer

BACKGROUND

We aimed to map and describe the current state of Mendelian randomization (MR) literature on cancer risk and to identify associations supported by robust evidence.

METHODS

We searched PubMed and Scopus up to 06/10/2020 for MR studies investigating the association of any genetically predicted risk factor with cancer risk. We categorized the reported associations based on a priori designed levels of evidence supporting a causal association into four categories, namely robust, probable, suggestive, and insufficient, based on the significance and concordance of the main MR analysis results and at least one of the MR-Egger, weighed median, MRPRESSO, and multivariable MR analyses. Associations not presenting any of the aforementioned sensitivity analyses were not graded.

RESULTS

We included 190 publications reporting on 4667 MR analyses. Most analyses (3200; 68.6%) were not accompanied by any of the assessed sensitivity analyses. Of the 1467 evaluable analyses, 87 (5.9%) were supported by robust, 275 (18.7%) by probable, and 89 (6.1%) by suggestive evidence. The most prominent robust associations were observed for anthropometric indices with risk of breast, kidney, and endometrial cancers; circulating telomere length with risk of kidney, lung, osteosarcoma, skin, thyroid, and hematological cancers; sex steroid hormones and risk of breast and endometrial cancer; and lipids with risk of breast, endometrial, and ovarian cancer.

CONCLUSIONS

Despite the large amount of research on genetically predicted risk factors for cancer risk, limited associations are supported by robust evidence for causality. Most associations did not present a MR sensitivity analysis and were thus non-evaluable. Future research should focus on more thorough assessment of sensitivity MR analyses and on more transparent reporting.

Association between telomere length and mitochondrial copy number and cancer risk in humans: A meta-analysis on more than 300,000 individuals

In the last decades the association of leukocyte telomere length (LTL) and mitochondrial copy number (mtDNAcn) with cancer risk has been the focus of many reports, however the relation is not yet completely understood. A meta-analysis of 112 studies including 64,184 cancer cases and 278,641 controls that analysed LTL and mtDNAcn in relation to cancer risk has been conducted to further our understanding of the topic. Stratified analyses for tumor type were also performed. Overall, no association was observed for all cancer combined neither for LTL nor mtDNAcn. Significant associations were detected for these biomarkers and specific cancer type; however, a large degree of heterogeneity was present, even within the same tumor type. Alternatives approaches based on polymorphic variants, such as polygenic risk scores and mendelian randomization, could be adopted to unravel the causal correlation of telomere length and mitochondrial copy number with cancer risk.

The Genetic Makeup of Myeloproliferative Neoplasms: Role of Germline Variants in Defining Disease Risk, Phenotypic Diversity and Outcome

Myeloproliferative neoplasms are hematologic malignancies typified by a substantial heritable component. Germline variants may affect the risk of developing a MPN, as documented by GWAS studies on large patient cohorts. In addition, once the MPN occurred, inherited host genetic factors can be responsible for tuning the disease phenotypic presentation, outcome, and response to therapy. This review covered the polymorphisms that have been variably associated to MPNs, discussing them in the functional perspective of the biological pathways involved. Finally, we reviewed host genetic determinants of clonal hematopoiesis, a pre-malignant state that may anticipate overt hematologic neoplasms including MPNs.

Clinical utility of targeted next-generation sequencing for the diagnosis of myeloid neoplasms with germline predisposition

Myeloid neoplasms (MN) with germline predisposition (MNGP) are likely to be more common than currently appreciated. Many of the genes involved in MNGP are also recurrently mutated in sporadic MN. Therefore, routine analysis of gene panels by next-generation sequencing provides an effective approach to detect germline variants with clinical significance in patients with hematological malignancies. Gene panel sequencing was performed in 88 consecutive and five nonconsecutive patients with MN diagnosis. Disease-causing germline mutations in CEBPα, ASXL1, TP53, MPL, GATA2, DDX41, and ETV6 genes were identified in nine patients. Six out of the nine patients with germline variants had a strong family history. These patients presented great heterogeneity in the age of diagnosis and phenotypic characteristics. In our study, there were families in which all the affected members presented the same subtype of disease, whereas members of other families presented various disease phenotypes. This intrafamiliar heterogeneity suggests that the acquisition of particular somatic variants may drive the evolution of the disease. This approach enabled high-throughput detection of MNGP in patients with MN diagnosis, which is of great relevance for both the patients themselves and the asymptomatic mutation carriers within the family. It is crucial to make a proper diagnosis of these patients to provide them with the most suitable treatment, follow-up, and genetic counseling.

Recent Advances in the Use of Molecular Analyses to Inform the Diagnosis and Prognosis of Patients with Polycythaemia Vera

Genetic studies in the past decade have improved our understanding of the molecular basis of the BCR-ABL1-negative myeloproliferative neoplasm (MPN) polycythaemia vera (PV). Such breakthroughs include the discovery of the JAK2V617F driver mutation in approximately 95% of patients with PV, as well as some very rare cases of familial hereditary MPN caused by inherited germline mutations. Patients with PV often progress to fibrosis or acute myeloid leukaemia, both associated with very poor clinical outcome. Moreover, thrombosis and major bleeding are the principal causes of morbidity and mortality. As a result of increasingly available and economical next-generation sequencing technologies, mutational studies have revealed the prognostic relevance of a few somatic mutations in terms of thrombotic risk and risk of transformation, helping to improve the risk stratification of patients with PV. Finally, knowledge of the molecular basis of PV has helped identify targets for directed therapy. The constitutive activation of the tyrosine kinase JAK2 is targeted by ruxolitinib, a JAK1/JAK2 tyrosine kinase inhibitor for PV patients who are resistant or intolerant to cytoreductive treatment with hydroxyurea. Other molecular mechanisms have also been revealed, and numerous agents are in various stages of development. Here, we will provide an update of the recent published literature on how molecular testing can improve the diagnosis and prognosis of patients with PV and present recent advances that may have prognostic value in the near future.

Genetically determined telomere length and multiple myeloma risk and outcome

Telomeres are involved in processes like cellular growth, chromosomal stability, and proper segregation to daughter cells. Telomere length measured in leukocytes (LTL) has been investigated in different cancer types, including multiple myeloma (MM). However, LTL measurement is prone to heterogeneity due to sample handling and study design (retrospective vs. prospective). LTL is genetically determined; genome-wide association studies identified 11 SNPs that, combined in a score, can be used as a genetic instrument to measure LTL and evaluate its association with MM risk. This approach has been already successfully attempted in various cancer types but never in MM. We tested the "teloscore" in 2407 MM patients and 1741 controls from the International Multiple Myeloma rESEarch (IMMeNSE) consortium. We observed an increased risk for longer genetically determined telomere length (gdTL) (OR = 1.69; 95% CI 1.36-2.11; P = 2.97 × 10^-6 for highest vs. lowest quintile of the score). Furthermore, in a subset of 1376 MM patients we tested the relationship between the teloscore and MM patients survival, observing a better prognosis for longer gdTL compared with shorter gdTL (HR = 0.93; 95% CI 0.86-0.99; P = 0.049). In conclusion, we report convincing evidence that longer gdTL is a risk marker for MM risk, and that it is potentially involved in increasing MM survival.

Telomere Length and Male Fertility

Over the past decade, telomeres have attracted increasing attention due to the role they play in human fertility. However, conflicting results have been reported on the possible association between sperm telomere length (STL) and leukocyte telomere length (LTL) and the quality of the sperm parameters. The aim of this study was to run a comprehensive study to investigate the role of STL and LTL in male spermatogenesis and infertility. Moreover, the association between the sperm parameters and 11 candidate single nucleotide polymorphisms (SNPs), identified in the literature for their association with telomere length (TL), was investigated. We observed no associations between sperm parameters and STL nor LTL. For the individual SNPs, we observed five statistically significant associations with sperm parameters: considering a p < 0.05. Namely, ACYP2˗rs11125529 and decreased sperm motility (p = 0.03); PXK˗rs6772228 with a lower sperm count (p = 0.02); NAF1˗rs7675998 with increased probability of having abnormal acrosomes (p = 0.03) and abnormal flagellum (p = 0.04); ZNF208˗rs8105767 and reduction of sperms with normal heads (p = 0.009). This study suggests a moderate involvement of telomere length in male fertility; however, in our analyses four SNPs were weakly associated with sperm variables, suggesting the SNPs to be pleiotropic and involved in other regulatory mechanisms independent of telomere homeostasis, but involved in the spermatogenic process.