Monitoring of telomere dynamics in peripheral blood leukocytes in relation to colorectal cancer patients' outcomes

doi:10.21203/rs.3.rs-1619264/v1

Download PDF

Article

Monitoring of telomere dynamics in peripheral blood leukocytes in relation to colorectal cancer patients' outcomes

https://doi.org/10.21203/rs.3.rs-1619264/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

We investigated the possible associations between leukocyte telomere length, therapy outcomes, and clinicopathological features in patients with colorectal cancer. Aditionally, telomerase reverse transcriptase (TERT) expression was evaluated.

Telomere length was measured using singleplex qPCR in 478 consecutive leukocyte DNA samples from 198 patients. Blood was drawn at diagnosis prior to any therapy, and then at 6-month intervals for 18 months.

Following diagnosis, the telomeres gradually shortened during the course of the treatment regardless of the patient's age. The most pronounced decrease was observed 12 months after the diagnosis (p < 0.0001). Based on tumor localization, the decrease in telomere length one year after the diagnosis followed different trajectories (p = 0.03). In patients treated with adjuvant therapy, telomere length was correlated with the time elapsed after completion of therapy (p = 0.03). TERT expression did not correlate with the telomere length; however, it was higher in women than men (1.35-fold, 95% CI 1.11–1.65, p = 0.003) and in smokers than non-smokers (1.27-fold, 95% CI 1.01–1.61, p = 0.04).

Leukocyte telomere length declines naturally during aging, but the accelerated shortening observed in our patients was age-independent. Telomere length manifestly reflected chemotherapy impact and could be linked to therapy toxicity.

leukocyte telomere length

longitudinal monitoring

colorectal cancer

TERT expression

adjuvant therapy

Telomeres, which are protective caps at the chromosomal ends, comprise of hexanucleotide tandem repeats in combination with the associated shelterin complex subunits. Telomere complexes protect chromosomal ends from DNA double-strand break repair machinery which can otherwise initiate inappropriate recombination events leading to extensive chromosomal rearrangements ¹. Telomeres, however, due to the intrinsic limitations of DNA replication undergo gradual age-dependent attritions ^2,3. Telomere function is also affected by oxidative DNA damage caused by endogenous and exogenous agents ^1,2. Excessive telomere shortening can trigger irreversible cell proliferation arrest, leading to the accumulation of senescencent cells and aging. Cancers develop mechanisms to bypass the replicative senescence and become immortal through upregulation of telomerase reverse transcriptase (TERT), which is a subunit of the holoenzyme that elongates telomeres ⁴.

Anticancer drugs perturb telomere length, either by their direct DNA-damaging activity or secondary-induced oxidative stress. Their effect depends on the treatment regimens and cancer types ^5–7. Currently, only a limited number of studies have focused on the effect on telomere length dynamics during the course of the therapy in solid tumors ⁵. Benitez-Buelga et al., observed a negative correlation between mean leukocyte telomere length (LTL) and the duration of cancer therapy in 489 patients with sporadic breast cancer ⁷. The shortening effect was significant in patients treated for > 90 days, and their LTL remained shortened until one year after the end of the treatment. However, other longitudinal studies on patients with solid tumors have not reported any association between therapy received and LTL ^8,9. To the best of our knowledge, only one study has investigated the therapeuic effect on LTL in colorectal cancer (CRC) patients, featuring, however, a very small (n = 32) and diverse cohort of solid cancer patients including only three CRC patients ¹⁰.

Upregulation of telomerase has been reported in hematopoietic cells following consecutive chemotherapy courses, and is interpreted as a mechanism to compensate for telomere attrition ¹¹. A study of childhood cancers by Franco et al. ¹² reported that telomerase activity does not change in any phase of treatment in solid tumors, but rather decreases sharply during the induction chemotherapy in acute lymphoblastic leukemia due to the elimination of the leukemic clone ¹².

Epidemiological studies have documented a complex association between disrupted LTL homeostasis (either extremely short or long LTL phenotypes) and different cancers ^2,13. This has resulted in an intensive investigation of LTL as a risk and prognostic marker for patients with cancer. Specifically, in CRC, there is a need for new biomarkers to predict therapeutic response, toxicity, and resistance. However, the question remains whether a variation in TERT expression and LTL, possibly affected by DNA-damaging anticancer drugs, can be the biological mechanism underlying the inadequate response to cancer therapy. Moreover, in CRC patients, the effect of treatment on the constitutive telomere length changes and TERT expression has not yet been investigated. Our study hypothesizes that LTL and TERT expression are surrogate markers of cancer therapy sensitivity, toxicity, and overall patient outcomes.

2.1 Population characteristics

One hundred and ninety-eight newly diagnosed, and histologically confirmed patients diagnosed with CRC (median age 66 years; range 32–88 years; 127 men; 71 women) were enrolled in the study. Peripheral blood samples were collected from all patients at Thomayer University Hospital in Prague, Czech Republic from June 2008 to May 2018. In total, 478 samples were collected from 198 patients. The characteristics of the study group are shown in Table 1. Blood specimens were drawn from patients at the baseline before any treatment (n = 137), followed by collection at the sixth (n = 156), twelfth (n = 135), and eighteenth (n = 50) month after diagnosis. The 6-month sampling encompassed, for the most part, patients (56%) for whom the anticancer therapy had been completed. An age-matched control group was not included, based on our previous results confirming that the lymphocyte telomere length of cancer patients is age-independent, in contrast to healthy controls ¹⁴.

Table 1. Patient´s characteristics. The table summarizes the characteristics of all patients who participated in the study. Number of patients do not always add up to 100% (n = 198) due to missing data for some attributes.

Studied cohort of patients	median age (years) [range]	66 [32-88]
		n = 198	%
Gender
	males	127	64.1
	female	71	35.9
Smoking status
	smokers	46	24.5
	non-smokers	142	75.5
Tumor site location^a
	proximal colon	35	18.0
	distal colon	86	44.3
	rectum	73	37.6
UICC TNM stage^b
	I + II	118	63.8
	III + IV	67	36.2
Microsatellite status
	stable	138	85.2
	instable	24	14.8
Therapy response
	good	64	70.3
	poor	27	29.7
Neoadjuvant therapy
	yes	66	33.3
	no	132	66.6
Adjuvant therapy
	yes	89	46.6
	no	102	53.4

Note: ^aPatients with CRC were categorized as having proximal colon (C18.0-18.4), distal colon (C18.5.-19.0), or rectal cancer (C20) according to the International Statistical Classification of Diseases and Related Health Problems 10^thRevision (ICD-10).

^bTumor-Node-Metastasis (TNM) staging system was developed by the Union for International Cancer Control (UICC) to classify the anatomical extent of tumor cancers.

In our study, we analyzed the oeffect of both neoadjuvant and adjuvant therapies on LTL. In general, for CRC, the type of treatment varies according to the tumor localization. The treatment strategy for colon cancer patients typically involves surgical removal, followed by adjuvant chemotherapy based on 5-fluorouracil (often in combination with oxaliplatin, irinotecan, or either of the above drugs with added folic acid; ¹⁵). In contrast, rectal cancer treatment is based on neoadjuvant (5-fluorouracil-based chemo)radiotherapy to downstage the tumor before surgical removal. Patients with locally advanced tumors may sometimes receive adjuvant chemotherapy (based on drugs described above for colon cancer treatment).

The cohort analyzed for the impact of 5-fluorouracil-based chemotherapy as standard care for CRC in the adjuvant setting comprised 42 individuals with a complete record of treatment given. To minimize the influence of other treatment strategies, only patients without previous neoadjuvant treatment were enrolled. The effect of neoadjuvant therapy was investigated in 23 patients with rectal cancer with complete treatment records. To eliminate the possible influence of other administered treatments, the cohort included only patients who did not received postoperative therapy. Neoadjuvant therapy for the patients involved (5-fluorouracil-based chemo)radiotherapy or chemotherapy with 5-fluorouracil used alone. The therapies prescribed for each patient are presented in Supplementary Table 1.

All participants provided written informed consent to participate in this study. Peripheral blood sampling was performed in accordance with the ethical standards outlined in the Declaration of Helsinki. Clinicopathological characteristics of all patients were collected from their medical records. Personal data were acquired using lifestyle questionnaires. Only collaborating clinicians could disclose the identities of the participants since they entered the anonymized study.

2.2 Peripheral Blood Collection, Processing, and DNA and RNA extraction

Blood was collected into Vacuette K3EDTA blood collection tubes by venipuncture. After sampling, fresh specimens were processed to produce mononuclear cells. These peripheral mononuclear cells were isolated from 2 ml of whole blood by Ficoll‑Paque PLUS (GE Healthcare Life Sciences) flotation, as recommended in the manufacturer's handbook, and frozen in TRI Reagent (Sigma-Aldrich). The remaining whole blood was stored at -80 °C.

Genomic DNA from whole blood (i.e. leukocytes) and total RNA from mononuclear cells were used as starting materials for qPCR-based analyses. DNA was extracted from 200 µL of blood, according to the DNeasy Blood and Tissue Kit protocol (Qiagen, Valencia, CA). RNA was extracted from Ficoll-isolated mononuclear cells using the standard TRI Reagent protocol.

DNA and RNA concentrations were determined using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific). RNA quality (RIN ≥ 8) was assessed using a Bioanalyzer 2100 (Agilent, Santa Clara, CA, USA) and the Agilent RNA 6000 Nano kit.

2.3 Measurement and Calculation of Telomere Length

Relative LTL was estimated using singleplex qPCR based on Cawthon ¹⁶ and later optimized by Gil and Coetzer ¹⁷. The method was adopted in collaboration with the Department of Biology and Medical Genetics, First Faculty of Medicine, Charles University, Prague, which has published articles on this topic ^18,19.

Briefly, LTL measurements were performed using the ViiA 7 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). qPCR reactions were ran on a MicroAmp Fast Optical 96-Well Reaction Plate (Life Technologies, Carlsbad, CA, USA) in triplicate. The 24 µL reactions consisted of PowerUp SYBR Green Master Mix (Life Technologies, Carlsbad, CA, USA), primers (300 nM for 36B4, 900 nM for telomeres), and 5 ng of DNA, and were prepared according to the manufacturer's instructions. To detect the reference 36B4 gene, which encodes the acidic ribosomal phosphoprotein P0, and telomere sequences, primers (Sigma-Aldrich, St. Louis, MO, USA) with the same annealing temperature were used. Primer sequences are listed in Supplementary Table 2. The reaction conditions were identical for both the targets. The holding stage was performed at 50 °C for 2 min and 95 °C for 2 min. This was followed by 40 cycles at 95 °C for 30 s and 54 °C for 1 min (ramp rate: 1.6 °C/s). The specificity of the primers was assessed by melting curve analysis ( Supplementary Figure 1) at the end of the amplification procedure. A calibration curve and interplate calibration were established using a 17 ng/µL reference DNA sample pooled from 30 healthy middle-aged hospital volunteers: 15 females and 15 males ranging from 41 to 55 years of age with no previous experience with cancer ²⁰. Participants were recruited at the Faculty Hospital Kralovské Vinohrady in Prague, Czech Republic. Amplification efficiency was verified using a 5-point calibration curve created through a 2‑fold serial dilution of the reference sample with three replicates at each point. The DNA amount ranged between 17 ng to 1.0625 ng per well. The reaction efficiency was calculated from the slope of the standard curves (see Supplementary Figure 2) using the equation: E = -1 + 10^(-1/slope) and evaluated using Viia7 software. The efficiency was 97% for telomeres and 101% for 36B4. The interplate and intraplate coefficients of variation were 1.03% and 0.35% for 36B4, and 1.92% and 0.24% for telomeres, respectively. Samples assayed in triplicate with the standard deviation (SD) ˃ 0.3 were omitted from the analysis, and their measurement was performed again. When two replicates were nearly identical and one was distant with a suspected dilution error, the outlier was omitted from the calculations.

LTL was estimated as the ratio of telomere repeat (T) copy number to single gene (S) copy number (T/S ratio) and calculated by the ΔΔCt method ²¹. The mean cycle threshold (Ct) of 36B4 was subtracted from the mean Ct of telomeres. The calculated ΔCt of each sample and ΔCt of the standard (derived from the ΔCts of standard replicates within the calibration curve) were substituted into the formula 2^{-(∆Ct(sample)-∆Ct(standard)}. The obtained T/S ratio is proportional to the average TL in the cell. A ratio above 1 indicated that the patient had longer telomeres than the reference sample.

2.4 Measurement and Calculation of TERT mRNA Expression

TERT mRNA expression levels were measured by two-step singleplex quantitative reverse‑transcription PCR (RT-qPCR) using the TaqMan chemistry.

Complementary DNA (cDNA) was synthesized using 300 ng οf total RNA with a High Capacity cDNA Reverse Transcription Kit (Life Technologies, Carlsbad, CA, USA) in 15 μL reactions. All components were combined and incubated in an MJ Research PTC-200 Thermal Cycler (Bio-Rad), following the manufacturer’s instructions. Subsequent qPCR reactions were carried out in a 384-well format using a MicroAmp Optical 384-Well Reaction plate (Life Technologies, Carlsbad, CA, USA) on the Viia7 Real-Time PCR System (Applied Biosystems). All reactions were performed in triplicates. To normalize TERT expression levels, endogenous control genes glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and beta-actin (ACTB), which are known to be stably expressed in mononuclear cells ²², were used. All 20 μl qPCR reactions included TaqMan Universal Master Mix II, with UNG (Applied Biosystems, Carlsbad, CA, USA), particular FAM/MGB TaqMan Gene Expression Assay (Hs00972650_m1 (TERT), Hs01060665_g1 (ACTB), or Hs02786624_g1 (GAPDH)), RNase-free water, and 2 μl of 3:1 diluted cDNA. Thermocycling conditions were 50 °C for 2 min and 95 °C for 10 min, followed by 40 cycles at 95 °C for 15 s and 60 °C for 1 min. The fluorescence signal was captured at the end of the 60 °C segments. Using the peripheral blood mononuclear cell RNA from a healthy donor, the interplate and intraplate coefficients of variation were 2.21% and 0.45% for TERT, 0.70% and 0.42% for GAPDH and 0.40% and 0.49% for ACTB. Technical replicates with the SD of the average Ct ˃ 0.3 were omitted from the study, and samples were measured in a new run.

The relative expression of TERT was calculated analogically to the previous description (see Section 2.3) using the mean Ct of the two housekeeping genes as a reference.

2.5 Statistical Methods

Statistical analysis was performed using STATISTICA (version 11Cz; TIBCO Software Inc., Palo Alto, CA, USA) and MATLAB (version 2019b; The MathWorks, Inc., Natick, MA, USA). All reported p-values were two-tailed and the level of statistical significance was set at α = 0.05.

Owing to their mostly normal distribution, LTL and log₂-transformed TERT expression values (i.e., -∆∆Ct values) were compared between relevant group pairs (according to sex, smoking habit, microsatellite instability, Union for International Cancer Control (UICC) tumor-node-metastasis (TNM) classification, and treatment response) using a two-sample t-test. The Mann-Whitney U test was used to compare the LTL of rectal cancer patients measured before and after neoadjuvant therapy (in an independent arrangement because of a small overlap between the patient groups). In the samples collected at the time of the diagnosis, an ordinal association between LTL and TERT expression was analyzed using the Kendall rank correlation coefficient. Kendall's τ rank test was also performed to investigate the relationship between age and LTL or TERT expression and the correlation between LTL and time from adjuvant therapy termination. One-way analysis of variance (ANOVA) was used to examine the differences in LTL and TERT expression in patients diagnosed with proximal colon, distal colon, and rectal cancer. To examine LTL changes over time, these diagnosis-related groups were further investigated using a repeated-measures ANOVA. Determination of whether LTL values of all patients changed over time was also performed using repeated-measures ANOVA as well. Lastly, repeated-measures ANOVA was used to compare LTL in patients with different TNM stages and treatment responses over time. Kruskal-Wallis ANOVA was used to compare the mean LTL values in time intervals from adjuvant therapy.

Overall survival (defined from diagnosis to death or censored at the time of the last follow-up) and its association with LTL were analyzed using the Kaplan-Meier method, Cox proportional-hazards model, and after stratification by LTL median, the Gehan-Wilcoxon test. The median follow-up time was 26.1 months (determined using the inverse Kaplan-Meier method).

2.6 Availability of data and materials

The datasets generated and/or analysed during the current study are not publicly available due patients' privacy but are available, in coded form, from the corresponding author on reasonable request.

3.1 Leukocyte telomere length dynamics in 12-month interval

The first three samplings (collected at diagnosis and 6 and 12 months after diagnosis) were available for 49 patients. Specifically, we monitored the overall LTL dynamics and dynamics across the patient subgroups. Monitoring of these 49 patients identified LTL shortening during the course of the treatment and follow-up examinations. In the whole group, there was a significant downward trend across the first three 6-month of observation periods independent of age, cancer stage, tumor localization, or treatment (repeated measures ANOVA, F(2, 96) = 16.93, n = 49, p < 0.0001). The Bonferroni post-hoc test revealed that the most pronounced decrease in LTL was between 6 and 12 months after the diagnosis (Bonferroni corrected p < 0.0001, Figure 1). Furthermore, patients analyzed in the subgroups categorized according to their tumor location showed different trajectories of LTL decreases. (repeated measures ANOVA F(4, 90) = 2.83, p = 0.03, proximal colon n = 11, distal colon n = 24, rectum n = 14, Figure 2). Concerning therapy response and cancer staging, we did not observe a difference in the rate of LTL decrease between good (n = 19) and poor responders (n = 4) (repeated-measures ANOVA, F(2, 42) = 1.30, p = 0.28), or between TNM I + II (n = 33) and III + IV (n = 13) (repeated-measures ANOVA, F(2, 88) = 0.52, p = 0.60). We also performed sampling 18 months after the diagnosis, but the four complete samples were drawn from 6 patients only.

3.2 Leukocyte telomere length in patients prior to therapy vs. clinicopathological data

We observed no correlation between LTL and age at diagnosis (Kendall's τ = 0.03, p = 0.63, n = 137), gender (females: 0.86 ± 0.35, n = 51; males: 0.86 ± 0.35, n = 86; two-sample t-test t(135) = -0.01, p = 0.99), or smoking status (smokers: 0.88 ± 0.34, n = 32; non-smokers: 0.87±0.35, n = 98; two-sample t-test t(128) = -0.13, p = 0.90). There was no difference in LTL between microsatelitte-stable (0.86 ± 0.34, n = 94) and unstable patients (0.82 ± 0.29, n = 22; two-sample t-test t(114) = 0.48, p = 0.63). Patients diagnosed with TNM stage I + II (0.84 ± 0.35, n = 82) and TNM stage III + IV (0.89 ± 0.34, n = 47) did not show statistically significant differences in LTL (two-sample t-test t(127) = -0.693, p = 0.49). Likewise, there was no difference in LTL in patients with tumors in the proximal colon, distal colon, and rectal locations (one-way ANOVA (F(2, 132) = 0.28, p = 0.76). The complete LTL results are listed in Table 2.

Table 2. Leukocyte telomere length in colorectal cancer patient subgroups at the time of the diagnosis. The table summarises all the results of LTL and stratifies the studied group of CRC patients based on the clinicopathological characteristics.

Clinicopathological characteristics of patients
		n	leukocyte telomere length mean ± SD	p value
Sampling collected at the time of the diagnosis		137	0.86 ± 0.35	-
Gender
	males	86	0.86 ± 0.35	0.99
	females	51	0.86 ± 0.35	0.99
Smoking status
	smokers	32	0.88 ± 0.34	0.90
	non-smokers	98	0.87 ± 0.35	0.90
Tumor site
	proximal colon	31	0.83 ± 0.36	0.76
	distal colon	69	0.87 ± 0.37
	rectum	36	0.83 ± 0.31
UICC TNM stage
	I + II	82	0.84 ± 0.35	0.49
	III + IV	47	0.89 ± 0:34	0.49
Microsatellite status
	stable	94	0.86 ± 0.34	0.63
	instable	22	0.82 ± 0.29	0.63
Therapy response
	good	41	0.92 ± 0.34	0.74
	poor	20	0.90 ± 0.30	0.74

3.3 Linking leukocyte telomere length to treatment response and therapy effect

In incident CRC cases, comparison of LTL between good (0.94 ± 0.34, n = 41) and poor (0.91 ± 0.30, n = 20) therapy responders revealed no difference (two-sample t-test t(59) = 0.34, p = 0.74, see Table 3). With regard to adjuvant therapy, LTL was inversely correlated with the time elapsed from the therapy completion (Kendall's τ = -0.22, p = 0.03), but the correlation was computed only by using n = 49 data points. Further, LTL was significantly non-uniform across the six time intervals defined with respect to the treatment (Kruskal-Wallis ANOVA H(5, 96) = 12.74, p = 0.03, Figure 3), although the post-hoc test did not confirm the differences between any of the time point pairs. Additionally, within the period < 90 days after the end of adjuvant therapy, we found that increase in LTL showed a borderline association with an increased overall survival (Cox proportional-hazards model, hazard ratio (HR) = 0.01, 95% confidence interval (CI) 0.00–0.01, p = 0.05, n = 16). However, stratification based on median LTL (cut-off value 0.94), in contrast, did not show that LTL positively affected overall survival of the patients (Gehan-Wilcoxon test, p = 0.06). When focusing on neoadjuvant therapy in rectal cancer cases only, LTL of the patients measured before (0.73±0.28, n = 7) and after the therapy administration (0.82±0.36, n = 32) did not differ (Mann-Whitney U = 94, p = 0.53). We also did not find any correlation with time after neoadjuvant therapy (Kendall's τ = -0.06, p = 0.65).

3.4 Leukocyte telomere length as predictor of overall survival rate

The 3-year overall survival rate of all CRC patients was 83.3% (95% CI 76.2–90.4 %). Investigatingly, overall survival according to telomere length measured in newly diagnosed patients suggested that LTL had no relation to the survival rate (Cox proportional-hazards model, HR = 1.36, 95% CI 0.33–5.59, p = 0.67, n = 137).

3.5 TERT expression in incident cases vs. clinicopathological data

We also measured TERT expression in mononuclear cells at the time of diagnosis. TERT expression did not correlate with LTL (Kendall's τ = 0.00, p = 1, n = 104). TERT expression was higher in females (n = 43) than males (n = 66) (two-sample t-test t(107) = -3.01, 1.35-fold, 95% CI 1.11–1.65, p = 0.003, Supplementary Figure 3) and smokers (n = 29) than non-smokers (n = 74) (two-sample t-test t(101) = -2.11, 1.27-fold, 95% CI 1.01–1.61, p = 0.04, Supplementary Figure 4). No correlation was observed between TERT expression and patient age (Kendall's τ = -0.08, p = 0.23, n = 109). Microsatelitte-unstable patients (n = 15) showed only negligibly lower TERT expression than microsatellite-stable patients (n = 75) (0.91-fold, 95% CI 0.67–1.23, two-sample t-test t(88) = 0.66, p = 0.51), as did TNM III + IV (n = 37) in comparison with TNM I + II (n = 66) (0.96-fold, 95% CI 0.79–1.17, two-sample t-test t(101) = 0.37, p = 0.72). When compared to patients with the tumors in the distal colon (n = 53), decrease in TERT expression was observed for the localization in proximal colon (n = 24, 1.12-fold, 95% CI 0.85–1.47) as well as in rectum (n = 31, 1.11-fold, 95% CI 0.86–1.44), however without statistically significant (one-way ANOVA F(2, 105) = 0.55, p = 0.58). Similarly, TERT was insignificantly decreased in poor therapy responders (n = 18) compared to good responders (0.87-fold, 95% CI 0.65–1.15, two-sample t-test t(45) = 0.96, p = 0.34).

In vitro studies have provided evidence that conventional chemotherapies may alter telomere homeostasis through their specific mechanisms of action ⁵. Thus, there is a biological plausibility that anticancer treatment may be a possible reason for permanent or transient telomere length changes, and these changes may reflect cancer patients' outcomes. Furthermore, telomere length presumably varies depending on the therapy regimen ⁷.

To our knowledge, only one longitudinal study has evaluated LTL before and after exposure to chemotherapy in CRC patients ¹⁰. In this study, LTL was shortened after standard-dose combination chemotherapy, and the degree of shortening corresponded to good therapy response and neutropenia severity. However, the study also included patients with other solid cancers (n = 29), and the number of patients with CRC was undersized (n = 3). Additionally, Garg et al. found that in CRC patients, a short LTL might predict mucositis as a side effect of of 5-FU in the adjuvant setting ²³.

In this study, we observed that CRC therapy exerted a moderate effect on LTL, particularly in association with adjuvant therapy. We recorded the LTL correlation with time elapsed after adjuvant therapy, but we did not observe any LTL differences between patients untreated and treated with neoadjuvant therapy or between good and poor therapy responders. The dynamics of LTL should be further considered in the context of white blood cell counts ²⁴, which was not available for the particular intervals of interest. One of the side effects of chemotherapy is that it can lower white blood cell levels, however, this effect is only short-lived, and the cells recover over time ²⁵. We hypothesized that the decrease and subsequent normalization of leukocyte counts after chemotherapy could affect the measured LTL as the newly produced leukocytes may have intact telomeres. The assessment of the telomere homeostasis with the treatment response among CRC patients must be conducted after a stringent control for phenotypic tumor heterogeneity and a large variety of treatment regimens. For instance, our previous study showed significant differences in telomere length measured in CRC tumors and adjacent mucosa, related to localization, microsatellite instability status, and staging ²⁰.

Of the limited number of studies investigating the effect of chemotherapy on telomere homeostasis in solid tumors, only two have suggested post-treatment LTL recovery ⁷. As Benitez-Buelga et al. reported that the recovery phase began approximately two years after diagnosis ⁷. There is a paucity of epidemiological evidence regarding this phenomenon in solid tumors. However, the results of studies on leukemia patients documented the post-treatment telomere recovery phase ⁵. Our study did not find any normalization of LTL after treatment in CRC patients monitored for approximately 12 months after diagnosis. Longer observations could not be obtained owing to the lack of samples collected 18 months after diagnosis.

In concordance with our results, the influence of tobacco smoking and sex on telomerase activity in peripheral blood mononuclear cells from healthy subjects was reported in a previous study, with higher telomerase activity in women and smokers ²⁶. This study also showed that the CpG sites in the TERT promoter are differentially methylated in smokers, suggesting that smoke-related oxidative stress has a possible epigenetic effect on TERT. In their study, smokers displayed longer telomeres than non-smokers, in addition to increased telomerase activity. The upregulated telomerase activity/TERT expression in females can be explained by estrogens stimulating telomerase expression ²⁷.

Overall, the main novelty of our work is that, as the first one, we investigated the effect of chemotherapy on LTL solely in patients with CRC. The study is unique because of its follow-up design, which enabled monitoring of LTL from patient diagnosis through treatment and follow-up intervals. However, the number of patients sampled every six months during the follow-up visits declined steadily. After the surgical procedure, patients most likely did not always attend the hospital regularly as planned; blood samples in some cases remained uncollected, and some patients moved to receive follow-up care at their place of residence or even died. We are aware that the most repeated samples, therefore, originate from patients with good therapeutic responses.

This study focused on mapping LTL in patients with CRC since their diagnosis, during receiving treatment, and follow-up care. Our research has demonstrated that LTL decreased 6 months after the approximate end of the therapy (i.e., 12 months from the diagnosis). Moreover, these LTL changes probably have different dynamics, depending on the anatomical localization of the tumors. Leukocytes also likely reflect the effect of adjuvant treatment, as we observed a correlation between LTL shortening and termination of adjuvant therapy. A large variance in treatment regimens suggests that it would be worthwhile to delve deeper into particular treatment strategies and study their possible effects separately. The higher TERT expression in women and smokers likely reflects the impact of estrogen hormones and tobacco smoking, respectively. Thus, estrogen- and smoke-induced oxidative stress may induce higher telomerase activity in peripheral blood mononuclear cells.

ACTB beta-actin

ANOVA one-way analysis of variance

cDNA complementary DNA

CI confidence interval

CRC colorectal cancer

Ct cycle treshold

GAPDH glyceraldehyde-3-phosphate dehydrogenase

HR hazard ratio

LTL leukocyte telomere length

SD standard deviation

T/S ratio telomere repeat (T) copy number to single gene (S) copy number

TERT telomerase reverse transcriptase

TNM Tumor-Node-Metastasis

UICC Union for International Cancer Control

Acknowledgments

We wish to thank the patients for their participation in the study.

Conflict of Interests

The authors declare no potential conflict of interests.

Ethics statement

All patients included in the study (recruited at the Thomayer University Hospital in Prague, Czech Republic) signed an informed consent form. The study was approved by ethics committees of the Institute of Experimental Medicine, The Czech Academy od Sciences, and the Institute for Clinical and Experimental Medicine and Thomayer Hospital. This study was performed in accordance with the principles of the Declaration of Helsinki.

Funding

We are grateful for the financial support from the Grant Agency of the Czech Republic (GACR 19-10543S, GACR 21-27902S), the Ministry of Health of the Czech Republic (AZV NV 18-03-00199 and AZV NU 21-03-00145), and the Charles University (project GA UK No. 120, UNCE/MED/006). This publication is part of a project that has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement No. 856620. This study was supported by the project No. CZ.02.1.01/0.0/0.0/16_019/0000787 „Fighting INfectious Diseases“, awarded by the MEYS CR.“, financed from EFRR.

Authors' contribution statement

KT: data curation (equal), funding acquisition (supporting), investigation (lead), vizualisation (equal), writing original draft (lead), writing review, and editing (equal). MK: data curation (equal), investigation (supporting), methodology (equal), vizualisation (equal), writing review, and editing (equal). AZ: methodology (equal), writing review, and editing (equal). MK: supervision (equal). VV: conceptualization (equal), writing review, and editing (equal). PS: resources (equal). LS: resources (equal). ML: resources (equal). KH: supervision (equal). VL: supervision (equal), PH: formal analysis (lead), writing review, and editing (equal). RK: supervision (equal), writing review, and editing (equal). LV: project administration (equal). PV: conceptualization (lead), funding acquisition (lead), project administration (equal), supervision (lead), writing review, and editing (equal).

Tomasova, K., Kroupa, M., Forsti, A., Vodicka, P. & Vodickova, L. Telomere maintenance in interplay with DNA repair in pathogenesis and treatment of colorectal cancer. Mutagenesis 35, 261–271 (2020).
Srinivas, N., Rachakonda, S. & Kumar, R. Telomeres and telomere length: a general overview. Cancers 12, 558 (2020).
Huang, Z. et al. Dynamics of leukocyte telomere length in adults aged 50 and older: a longitudinal population-based cohort study. GeroScience 43, 645–654 (2021).
Heidenreich, B. & Kumar, R. TERT promoter mutations in telomere biology. Mutat. Res. 771, 15–31 (2017).
Gallicchio, L., Gadalla, S. M., Murphy, J. D. & Simonds, N. I. The effect of cancer treatments on telomere length: a systematic review of the literature. J. Natl. Cancer Inst. 110, 1048–1058 (2018).
Lu, Y., Leong, W., Guérin, O., Gilson, E. & Ye, J. Telomeric impact of conventional chemotherapy. Front. Med. 7, 411–417 (2013).
Benitez-Buelga, C. et al. Impact of chemotherapy on telomere length in sporadic and familial breast cancer patients. Breast Cancer Res. Treat. 149, 385–394 (2015).
Duggan, C. et al. Change in peripheral blood leukocyte telomere length and mortality in breast cancer survivors. J. Natl. Cancer Inst. 106, dju035 (2014).
Brouwers, B. et al. The impact of adjuvant chemotherapy in older breast cancer patients on clinical and biological aging parameters. Oncotarget 7, 29977–29988 (2016).
Yoon, S. Y. et al. Telomere length shortening of peripheral blood mononuclear cells in solid-cancer patients undergoing standard-dose chemotherapy might be correlated with good treatment response and neutropenia severity. Acta Haematol. 118, 30–37 (2007).
Szyper-Kravitz, M. et al. Granulocyte colony-stimulating factor administration upregulates telomerase activity in CD34 ⁺ haematopoietic cells and may prevent telomere attrition after chemotherapy. Br. J. Haematol. 120, 329–336 (2003).
Franco, S. et al. Telomere dynamics in childhood leukemia and solid tumors: a follow-up study. Leukemia 17, 401–410 (2003).
Calado, R. T. &Young, N. S. Telomere Diseases. N. Engl. J. Med. 361, 2353–2365 (2009).
Vodenkova, S. et al. Chromosomal damage and telomere length in peripheral blood lymphocytes of cancer patients. Oncol. Rep. 44, 2219–2230 (2020).
Vodenkova, S. et al. 5-fluorouracil and other fluoropyrimidines in colorectal cancer: past, present and future. Pharmacol. Ther. 206, 107447 (2020).
Cawthon, R. M. Telomere length measurement by a novel monochrome multiplex quantitative PCR method. Nucleic Acids Res. 37, e21 (2009).
Gil. M. E. & Coetzer, T. L. Real-time quantitative PCR of telomere length. Mol. Biotechnol. 27, 169–172 (2004).
Zinkova, A., Brynychova, I., Svacina, A., Jirkovska, M. & Korabecna, M. Cell-free DNA from human plasma and serum differs in content of telomeric sequences and its ability to promote immune response. Sci. Rep. 7, 2591 (2017).
Zinkova, A. et al. Relative amount of telomeric sequences in terminal villi does not differ between normal term placentas and placentas from patients with well-controlled type 1 diabetes mellitus. Placenta 55, 1–4 (2017).
Kroupa, M. et al. Relationship of telomere length in colorectal cancer patients with cancer phenotype and patient prognosis. Br. J. Cancer 121, 344–350 (2019).
Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods San Diego Calif. 25, 402–408 (2001).
Chege, D. et al. Evaluation of a quantitative real-time PCR assay to measure HIV-specific mucosal CD8 + T cell responses in the cervix. PLoS ONE 5, e13077 (2010).
Garg, M. B. et al. Predicting 5-fluorouracil toxicity in colorectal cancer patients from peripheral blood cell telomere length: a multivariate analysis. Br. J. Cancer 107, 1525–1533 (2012).
Lin, Y. et al. Age-associated telomere attrition of lymphocytes in vivo is co-ordinated with changes in telomerase activity, composition of lymphocyte subsets and health conditions. Clin. Sci. 128, 367–377 (2015).
Wu, J. et al. Values of applying white blood cell counts in the prognostic evaluation of resectable colorectal cancer. Mol. Med. Rep. 19, 2330–2340 (2019).
Marcon, F. et al. Telomerase activity, telomere length and hTERT DNA methylation in peripheral blood mononuclear cells from monozygotic twins with discordant smoking habits. Environ. Mol. Mutagen. 58, 551–559 (2017).
Bayne, S. et al. Estrogen deficiency leads to telomerase inhibition, telomere shortening and reduced cell proliferation in the adrenal gland of mice. Cell Res. 18, 1141–1150 (2008).

Table 3 is not available with this version

No competing interests reported.

Download PDF

Version 1

posted

You are reading this latest preprint version

Monitoring of telomere dynamics in peripheral blood leukocytes in relation to colorectal cancer patients' outcomes

Status:

Version 1

Abstract

Figures

1 Introduction

2 Material And Methods

2.1 Population characteristics

2.2 Peripheral Blood Collection, Processing, and DNA and RNA extraction

2.3 Measurement and Calculation of Telomere Length

2.4 Measurement and Calculation of TERT mRNA Expression

2.5 Statistical Methods

3 Results

3.1 Leukocyte telomere length dynamics in 12-month interval

3.2 Leukocyte telomere length in patients prior to therapy vs. clinicopathological data

4 Discussion

5 Conclusion

Abbreviations

Declarations

References

Table 3

Additional Declarations

Supplementary Files

Status:

Version 1