Simon, Rachel A

Background: Telomeres are repetitive sequences at the end of linear chromosomes that preserve genetic information. Patients with telomere biology disorders (TBDs) are subject to multiorgan manifestations, which include the liver. However, due to the rare nature of the disorder, these manifestations are not fully characterized. Methods: Here, we describe the hepatic manifestations of patients with TBD at a large tertiary center within a specialized TBD clinic; 121 consecutive patients with TBD seen from 2017 to 2024 were included. The median age of the TBD diagnosis was 43 years old, with 60% male. The most common pathogenic/likely pathogenic variant was TERT (25%), followed by RTEL1 (10%) and TERC (7%). Results: Evidence of hepatic involvement in TBD was observed in 60% of patients (73 of 121), as determined by imaging studies or abnormal liver chemistries. Hepatic steatosis was the most prevalent imaging finding (56%), while a cholestatic pattern of liver injury was the predominant laboratory abnormality (73%). Biopsy-proven portosinusoidal vascular disorder was present in 4 patients. Findings of portal hypertension were observed in 25 patients (34%), with decompensation observed in 12 patients (16%). Traditional liver imaging modalities (computer tomography, ultrasound, and MRI) were compared with magnetic resonance elastography in 41 patients. Of these, 7 patients had abnormal magnetic resonance elastography with normal traditional imaging (17%), indicating higher sensitivity of magnetic resonance elastography to detect liver involvement. Conclusions: This study describes the liver manifestations within a cohort of patients with TBD with findings of portosinusoidal vascular disorder and hepatopulmonary syndrome similar to other cohorts and a noted higher prevalence of hepatic steatosis as compared with the general population. Liver stiffness measurements, integrated with imaging and laboratory testing, further improve characterization of hepatic manifestations in TBD.

Introduction: Telomere biology disorders (TBD) are multisystemic accelerated aging syndromes due to abnormal telomere length (TL) shortening and increased cancer risk. Pediatric-onset TBD is characterized by severely shortened TL and specific associations with germline pathogenic variants mainly in DKC1, TINF2, PARN genes and characteristic clinical features [bone marrow failure (BMF), abnormal nail/skin/oral mucosa]. Adult-onset TBD is genetically distinct (predominantly, TERT, TERC, and RTEL1 mutations) with varied organ involvement [interstitial lung disease (ILD), BMF, cirrhosis, or bone disease] and heterogenous outcomes. In our previous study of predominantly adult patients with short telomeres, multiple organ involvement (Mangaonkar A et al. Blood Cancer J 2021) predicted poor overall survival (OS). We sought to validate these findings in a larger, well annotated adult cohort of TBD patients. Methods: Patients with suspected TBD underwent comprehensive clinical and molecular evaluation through our TBD-dedicated clinic after regulatory approval. TL was assessed through clinical flow cytometry fluorescence in-situ hybridization (FlowFISH).Clinical and laboratory features, including pulmonary imaging, lung volumes/function testing, hepatic imaging (MR elastography/fibroscan) and bone density assessments were retrospectively abstracted. Germline and somatic genetic testing was performed on peripheral blood through targeted panels for known TBD and myeloid genes respectively. Cox proportional hazard model was used to identify significant variables that also retained independent prognostic impact. Results: Overall, 122 patients met criteria for TBD as defined by lymphocyte TL (LTL) <=10th percentile and/or a presence of a germline-TBD variant; 34 patients LTL <1st percentile, 84 1st and the 10th percentile (27 with TBD-germline variants) and 4 with only TBD-germline variants. Median age at diagnosis was 56 (18-77) years, with 72 (59%) males. Clinical features included presence of ILD (n=81, 66%), BMF (n=39, 32%), cirrhosis/hepatic fibrosis (n=19, 38%), bone disease osteopenia/osteoporosis/fractures, n=53, 70%), MDS/AML (n=11, 9%), non-hematological cancer (n=38, 31%), and head and neck cancer (n=2). Germline pathogenic variants were detected in 22 (18%), most frequently in TERT (n=10), RTEL1 (n=7) and TERC (n=2), while variants of uncertain significance (VUS) were seen in 35 (29%) patients, especially in TERT (n=20), TERC/RTEL1(n=6), NAF1 (n=4). Somatic variants in hematopoietic cells were detected in 37 (30%) patients with frequent alterations seen in DNA damage response and repair (PPM1D/TP53, n=11), and splicing component (U2AF1/SF3B1/ZRSR2, n=9) genes. At last follow-up [median of 3 (0-18) years], there were 26 (21%) deaths (most commonly, progressive ILD-18, cancer of lung/breast/pancreas-3, transplant-related complications-2). After censoring patients at the time of any transplantation (lung transplant-26, liver transplant-5, bone marrow transplant-8), the Kaplan-Meier estimate of median OS was not reached. Univariate survival analysis found adverse OS predicted by older age at diagnosis (P=0.006), presence of ILD (P=0.0002), hepatic fibrosis/cirrhosis (P=0.02), and somatic DNA repair variants (P=0.03). Multivariate model analysis found that only presence of ILD retained independent prognostic impact (HR=2, P=0.02). Among radiographic ILD subtypes, UIP was the most common (51%), followed by NSIP (10%), and others (6%), with median FEV1 of 2L (range: 0.5-4.6), FVC 2.4L (0.7-6.2), and TLC 4L (1.9-7.4). TBD-ILD patients were older (P<0.0001) with a male preponderance (P=0.03), 52 (64%) with multi-organ involvement, and had a less severe hematological phenotype [lower incidence of BMF (P<0.0001) and fetal hemoglobin elevation (P=0.01)], more frequent bone disease (P=0.0009), relatively longer telomeres in lymphocytes (P=0.04), and no particular association with any germline or somatic genetic alterations. Conclusions: Presence of ILD is an independent and significant predictor of adverse survival in adult TBD. Within adult TBD, patients with ILD represent a unique clinical subtype without any obvious genetic or pathologic correlates. Focusing on definitive therapeutic measures in this subgroup of adult TBD has the best chance to alter natural history and improve outcomes.

Background: Age-related clonal hematopoiesis (ARCH) is a precursor state to blood cancers and is associated with epigenetic age acceleration (EpiAgeAccel) (Aging Cell, 2021). While ARCH in combination with AgeAccel has the highest risk for all cause mortality and cardiovascular disease, not much is known about context-relevant clonal hematopoiesis (CRCH) - CH occurring in the context of prior chemo/radiation, bone marrow failure and germline predisposition syndromes. CRCH has a unique mutational signature enriched in genes in DNA damage response and repair (TP53 and PPM1D). To answer this question, we assessed EpiAge in a prospective cohort of patients with both ARCH and CRCH. Methods: Patients were prospectively recruited for CH and EpiAge assessment from the Mayo Clinic CH clinic between 2022 and 2024. Demographic and clinical data including chronologic age, BMI, sex, blood counts, and medical history were recorded at the time of enrollment. CH status was determined using error corrected sequencing on peripheral blood (PB) DNA. All somatic variants with a variant allele frequency (VAF) ≥0.5% were considered as CH. For EpiAge analysis PB DNA was analyzed using the Illumina MethylationEPIC array. The resulting beta values were used to obtain epigenetic age estimates using Horvath's clock (Genome Biol., 2013), Hannum's Clock (Mol Cell., 2013), PhenoAge (Aging, 2018), DunedinPACE (eLife, 2022), and predicted telomere length (Aging, 2019). Results: 106 (62%) of 170 patients assessed were found to have CH, with the most frequent mutations being TET2 (19%), PPM1D (17%), DNMT3A (16%), ASXL1 (9%), and TP53 (6%); median VAF 9.5%. Except for a higher median chronologic age (68 vs 60 years, p<0.001) in the CH group, baseline characteristics including BMI, comorbidities, baseline blood counts, sex and tobacco use history were similar between the two groups. The CH group demonstrated significant EpiAgeAccel by all three clocks: Horvath (p=0.001), Hannum (p<0.001), and PhenoAge (p<0.001); with an increase in biological age of 4.9 years (95% CI 2-7.7) by Horvath clock, 4.8 years (95% CI 2.2-7.5) by Hannum clock, and 7 years (95% CI 2.9-10.5) by PhenoAge, after adjusting for chronologic age, sex, smoking status, and cancer history. CH was also associated with increased DunedinPACE of aging by 0.13 biological year per chronologic year (bioyear/chronoyear) (95% CI 0.06-0.2, p<0.001), and shortened telomere lengths by 0.24kb (95% CI 0.11-0.35, p<0.001) after adjusting for chronologic age and excluding patients with telomere biology disorders. CH patients were then further classified into ARCH-63 (59%) and CRCH-43 (41%). While the median chronological age (69 vs 67 years) was similar between the 2 groups, CRCH patients had lower Hgb [10.9g/dL vs. 13g/dL (p=0.0001)] , lower platelet [146x109/L vs. 192x109/L (p=0.0132)] and lower WBC [4.4x109/L vs. 5.2x109/L (p=0.0345)]. The median BMI, comorbidities, and tobacco use history were similar between the 2 groups. The 3 most common ARCH mutations were TET2 (20%), DNMT3A (19%), and ASXL1 (14%), while, PPM1D (27%), TET2 (18%), and DNMT3A (14%) were frequent in CRCH. There was no significant difference found in median EpiAge via Horvath (68 vs 67.5years p=0.2826), Hannum (65 vs 66 years p=0.9333), or PhenoAge (55 vs 56 years p=0.618) between ARCH and CRCH. However, DunedinPACE of aging was higher in CRCH (1.27 vs 1.18 bioyear/chronoyear p=0.0148) and median telomere length was shorter in CRCH (6.35 vs 6.49kb p=0.0247), in comparison to ARCH. Multivariable regression analysis was completed for both ARCH and CRCH compared to non-CH, adjusting for chronologic age, smoking status, and male sex. While both ARCH and CRCH demonstrated EpiAgeAccel via Horvath(5.8years p=0.0001, 3.4 years p=0.051), Hannum(5years p=0.001, 4.5 years p=0.007), and PhenoAge (6.7years p=0.002, 7.6 years p=0.002) clocks, CRCH more than doubled DunedinPACE (0.2 bioyear/chronoyear, p<0.001) of aging compared to ARCH (0.08 bioyear/chronoyear, p=0.026). CRCH also demonstrated double the decrease in telomere length (0.34kb, p<0.001) compared to ARCH (0.17kb, p=0.015). Conclusion: Our prospective CH cohort study validates the association between CH and EpiAgeAccel using several different methylation clocks and for the first time, highlights increased pace of aging and shorter telomere lengths in CRCH versus ARCH. Correlations with cardiovascular disease and mortality outcomes are ongoing.

Background: It has been recently demonstrated that altered lipid metabolism involving increased free fatty acids (FFAs) uptake supports the growth and proliferation of MM cells (Panaroni et al. Blood 2022; Farell et al. Elife 2023). Whether the availability of FFAs in the bone marrow (BM) microenvironment is associated with the pathogenesis of MM from MGUS has not been previously evaluated. Thus, we evaluated if individual species of FAs or FA-containing lipids such as triacylglycerols (TAGs), diacylglycerols (DAGs), and monoacylglycerols (MAGs) are differentially present between the BM microenvironment of MGUS and MM patients. Methods: Lipid profiles of BM plasma samples (reflective of the BM microenvironment) from patients with a spectrum of plasma cell disorders (22 MGUS and 24 MM (15 newly diagnosed and 9 relapsed/refractory)) were acquired by liquid chromatography (LC)/tandem mass spectrometry (MS), and individual lipid species were quantified based on single-point quantification using 55 isotopically labeled internal standards. Concentration-based lipid species data was analyzed by the web-based tool MetaboAnalyst 5.0. Data was normalized by quantile normalization, log transformed and scaled by Pareto scaling. p values <0.05 were considered significant. Results: A total of 993 distinct lipid species were quantified by the MS-based targeted lipidomics assessment. Of these lipid species, the levels of individual mid-chain (≤50 carbon chain length) and long-chain (≥52 carbon chain length) polyunsaturated TAGs were the most differentially present in the BM plasma of MM compared to MGUS patients (Figure 1A). Even though a principal component analysis (PCA) involving all lipid species showed moderate separation of MM from MGUS patients, five MGUS patients (categorized as MGUS-other (MGO)) who clustered more towards MM patients compared to the remaining MGUS patients had a distinct pattern of individual mid and long-chain polyunsaturated TAG levels that were comparable to MM patients (Figure 1B). Furthermore, the proliferation rate of clonal plasma cells in the BM of these MM and MGUS as measured by the S-phase% was found to be directly correlated with the level of long-chain polyunsaturated TAGs. Conclusions: Mid and long-chain polyunsaturated TAGs are present in lower and higher levels, respectively, in the BM microenvironment of MM compared to MGUS patients. These findings continue to support the importance of FA-containing TAGs and other lipids in the pathogenesis of MM. Future confirmatory studies evaluating the biological significance of these differences are required to determine if they can be exploited for diagnostic and therapeutic purposes.

Introduction: Telomere biology disorders (TBD) are multisystem diseases often associated with bone marrow failure (BMF) and pulmonary fibrosis. TBD are clinically diagnosed when mean telomere lengths (TL) in blood cells measured by clinical flowFISH assays fall <1 st percentile, in comparison to age-matched reference controls. This information has limitations, e.g. TL in blood cells may not represent TL in other organs (e.g. lung, liver) and moreover, individuals with high clinical suspicion of TBD can show borderline average TL which complicates their diagnosis. The likely reason for this is that cell senescence can be triggered by shortening of a subset, or even a single chromosome arm telomere, a situation that may be missed when measuring average TL. Previous approaches to measure individual chromosome-arm telomere length (CA-TL) are only available for a small number of chromosome arms or can only estimate the number of shortened telomeres without matching them to specific chromosome arms. Taking advantage of recent developments in long-read next generation sequencing (NGS) and the availability of the recently published telomere-to-telomere (T2T) reference human genome, our group developed Telogator, (Stephens Z. et al, 2022) a bioinformatic tool that estimates CA-TL for most chromosome arms using long-read NGS data. Here we use Telogator to describe for the first time the CA-TL pattern in TBD patients in comparison with a control non-TBD cohort. Methods: TBD patients were identified through the Pre-Myeloid and Bone Marrow Failure Specialty Clinic at Mayo Clinic. All patients underwent flowFISH and genetic testing. Age-matched non-TBD patients identified at the Center for Individualized Medicine (Mayo Clinic) were used as reference controls. Informed consent was obtained from all individuals with approval from the Mayo Clinic Institutional Review Board. DNA from peripheral blood mononuclear cells was sequenced using PacBio Sequel II following the manufacturer indications for HiFi sequencing. BAM files were aligned to the T2T reference and CA-TL was calculated using Telogator as described in Stephens Z. et al. 2022. Average TL, CA-TL, and percent of short telomeres (both defined as <5 kb and <3.2 kb) were calculated. Comparisons of average TL and percent of short telomeres were performed using t-test, while Wilcoxon rank-sum test was used to compare CA-TL. Results: Four TBD patients [median age = 58.5 (52-60), 50% female] and 6 non-TBD controls [Median age = 56 (50-63), 33% female] were included in the analysis. All TBD patients presented average TL below the 1 st centile in both lymphocytes and granulocytes, and TBD-related phenotype (100% with pulmonary fibrosis, 50% with BMF, 25% with liver fibrosis). Three had a pathogenic variant in a TBD related genes [2 in TERT (NM_198253.3: c.2030G>A, p.(Gly677Asp) and c.2515_2533del, p.(Thr839Alafs*29)) and 1 in NHP2 (NM_017838.3: c.459_460delAT; p.(X154Argext29))]. The average TL calculated by Telogator showed high correlation with flowFISH values (r 2 = 0.9883, p = 0.0059) and was lower in TBD patients than controls (p = 0.0012). Next, when comparing individual CA-TL between groups, only half of CA-TL was statistically different in the TBD group compared to controls (1p, 3q, 4p, 5p, 5q, 6q, 7p, 7q, 8q, 9p, 9q, 10q, 11p, 12q, 12q, 16p, 16q, 18q, 19p, 19q, 20p, 20q, 21q - Figure 1). The percentage of telomeres <5kb and <3 kb was significantly increased in TBD patients compared to controls (p= 0.014 and 0.0333, respectively). Discussion: Using a novel informatic tool, we replicate previous results obtained using classic approaches and were able to expand on these studies by measuring specific CA-TL in TBD patients. We identified statistically significant differences at the CA-TL level in these patients when compared to age-matched controls, while other CA-TL were not statistically different. Future studies and the inclusion of additional TBD patients will help define the clinical relevance of the differential CA-TL shortening and its prevalence among TBD patients. While limited by the small number of samples in the study, we show the feasibility of using this approach to study TL with unprecedented resolution and show that Telogator offers important advantages compared to classic methods to study telomere biology that can be used in other relevant fields like aging and/or cancer. Our results pave the way to advance the mechanistic understanding of TBD.

INTRODUCTION Telomere biology disorders (TBDs) are multisystem diseases caused by accelerated telomere shortening that can present with bone marrow failure and pulmonary fibrosis, among other abnormalities. TBD diagnosis is considered when average blood telomere length (TL) measured by clinical flowFISH is <1 st percentile and/or a pathogenic variant in telomere maintenance genes (most commonly TERT, TERC, RTEL1) is identified. Pathogenic variants in genes coding for the telomerase protein ( TERT) and RNA component ( TERC) decrease telomerase activity. Similarly, variants in RTEL1 result in defective unwinding of the t-loop at telomere ends, increasing the levels of free t-circles in cells. Variants in these genes are often classified as variants of uncertain significance (VUSs by ACMG criteria) with additional efforts needed for variant curation. To address this, we describe results of functional testing performed in TBD patients with VUSs involving TERT, TERC, and RTEL1. METHODS TBD suspected patients were identified through the Pre-Myeloid Cancer and Bone Marrow Failure Clinic (Mayo Clinic) where they underwent clinical flowFISH and genetic testing for TBD related genes. Informed consent was obtained with approval from the Mayo Clinic Institutional Review Board. Viable peripheral blood and bone marrow samples were collected along with age matched healthy volunteers. Telomerase activity was measured in samples with TERT and TERC variants through telomeric repeat amplification protocol (TRAP) assays using the TRAPeze XL Telomerase Detection Kit (MilliporeSigma) following manufacturer instructions. HCT116 cells were used as positive controls while heat inactivated samples served as negative controls. Fluorescence was measured in a SpectraMax plate reader. Percentage of total product generated (TPG) relative to healthy age matched controls was calculated. T-circle detection assays were used on RTEL1 samples based on the protocol by Zhang et al. (2017). Final southern blotting was completed using the TeloTAGGG Telomere Length Assay Kit (Roche), and t-circle bands were visualized in a c600 GEL Imaging System (Azure Biosystems). RESULTS Ten patients (mean age = 53 years (27-71), 50% female) were included in our study. Patients 1 to 5 presented VUS in TERT (n=4) or TERC (n=1) and TL <1 st percentile in lymphocytes and/or granulocytes. Patient 6 carried a known pathogenic variant in RTEL1 and was included as a positive control. Patients 7 to 9 with RTEL1 variants presented TL close to the 10 th or between the 1 st-10 th percentile. No TL results were available for patient 10 but was tested due to the presence of a VUS in RTEL1 (Table 1). Similarly, while patient 7 carried a likely pathogenic variant, it was tested due to presenting lymphocyte TL between the 1 st-10 th percentile. TRAP results indicated decreased telomerase activity compared to age matched controls in all samples with TERT and TERC variants suggesting a pathogenic effect (52.%, 10.7%, 48.7%, 4.4% and 5.2% of TPG compared to age matched controls, respectively) (Table 1). On the other hand, t-circles were only seen in patient 6 (Figure 2). DISCUSSION We used TRAP and t-circle assays to explore the telomerase activity and presence of t-circles in patient samples with TERT/TERC and RTEL1 variants. In the case of TERT/TERC, we observed a reduction of telomerase activity in all samples. Together with blood TL below the 1 st percentile, this suggests that these variants are indeed deleterious (Table 1). However, since patient 4 presented two TERT VUSs, additional testing is needed to conclude how each variant impacts telomerase activity independently. In regard to RTEL1 variants, we observed t-circles in patient 6 as expected. No t-circles could be seen on all remaining samples, including patient 7 who carried a variant classified as likely pathogenic (Figure 1). This in-frame deletion has been described once in a TBD family (Cogan et al. 2015) but was not functionally tested. This result, together with the borderline TL seen in this patient, may suggest that some compensatory mechanism exists on this sample, or that the deleterious effect of the variant is not strong enough to be detected by our assay. Further studies are currently ongoing by our group to validate these results using a cell-based transfection disease model. In summary, we show that functional testing can be a helpful tool to determine pathogenicity of VUSs in TERT/TERC and RTEL1.