Next Article in Journal
The Genetic and Molecular Analyses of RAD51C and RAD51D Identifies Rare Variants Implicated in Hereditary Ovarian Cancer from a Genetically Unique Population
Next Article in Special Issue
Targeted Long-Read Bisulfite Sequencing Identifies Differences in the TERT Promoter Methylation Profiles between TERT Wild-Type and TERT Mutant Cancer Cells
Previous Article in Journal
Impact of the Tumor Microenvironment for Esophageal Tumor Development—An Opportunity for Prevention?
Previous Article in Special Issue
A Comparative Assessment of Replication Stress Markers in the Context of Telomerase
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relative Telomere Length Change in Colorectal Carcinoma and Its Association with Tumor Characteristics, Gene Expression and Microsatellite Instability

1
Institute for Population and Precision Health, Department of Public Health Sciences, Biological Sciences Division, The University of Chicago, Chicago, IL 60637, USA
2
Department of Pathology, Jahurul Islam Medical College, Kishoregonj 2336, Bangladesh
3
Department of Pathology, The Laboratory, Dhaka 1205, Bangladesh
4
Department of Surgery, Bangabandhu Sheikh Mujib Medical University, Dhaka 1000, Bangladesh
5
Department of Pathology, Cox’s Bazar Medical College, Cox’s Bazar 4700, Bangladesh
*
Author to whom correspondence should be addressed.
Cancers 2022, 14(9), 2250; https://doi.org/10.3390/cancers14092250
Submission received: 3 April 2022 / Revised: 27 April 2022 / Accepted: 28 April 2022 / Published: 30 April 2022
(This article belongs to the Special Issue The Dual Roles of Telomeres and Telomerase in Aging and Cancer)

Abstract

:

Simple Summary

Telomeres are specialized repeat DNA sequences at the ends of each chromosome. The primary function of telomeres is to prevent genomic instability. Telomeres shorten with every cell division. The role of relative telomere length (RTL) change in colorectal cancer (CRC) is not fully understood. We studied RTL in CRC tissue and adjacent normal tissue from the same patients (n = 165). We found significant shortening of RTL in cancer tissue, irrespective of age group, gender, tumor pathology, location and microsatellite instability (MSI) status. However, shortening was more pronounced in low-grade tumors and in the presence of MSI. Gene expression data suggested an association of telomere shortening with rapid cell division (upregulation of DNA replication and cyclin-dependent kinase genes), as well as downregulation of apoptosis genes. CRC tissue showed upregulation of some telomere maintenance genes.

Abstract

We compared tumor and adjacent normal tissue samples from 165 colorectal carcinoma (CRC) patients to study change in relative telomere length (RTL) and its association with different histological and molecular features. To measure RTL, we used a Luminex-based assay. We observed shorter RTL in the CRC tissue compared to paired normal tissue (RTL 0.722 ± SD 0.277 vs. 0.809 ± SD 0.242, p = 0.00012). This magnitude of RTL shortening (by ~0.08) in tumor tissue is equivalent to RTL shortening seen in human leukocytes over 10 years of aging measured by the same assay. RTL was shorter in cancer tissue, irrespective of age group, gender, tumor pathology, location and microsatellite instability (MSI) status. RTL shortening was more prominent in low-grade CRC and in the presence of microsatellite instability (MSI). In a subset of patients, we also examined differential gene expression of (a) telomere-related genes, (b) genes in selected cancer-related pathways and (c) genes at the genome-wide level in CRC tissues to determine the association between gene expression and RTL changes. RTL shortening in CRC was associated with (a) upregulation of DNA replication genes, cyclin dependent-kinase genes (anti-tumor suppressor) and (b) downregulation of “caspase executor”, reducing apoptosis.

1. Introduction

Telomeres are the terminal DNA–protein complex of chromosomes composed of repeat sequences of hexamers (TTAGGG) [1]. Proteins associated with telomeres, such as the shelterin complex, protect chromosomal ends from double-strand DNA breaks [2]. The primary function of telomeres is to prevent genomic instability [3]. Telomeres shorten with every cell division, ultimately limiting replication and potentially triggering apoptotic cell death. Telomere shortening correlates with age. The role of telomere dysfunction in colorectal carcinoma (CRC) is not fully understood. In case–control studies, it was found that relative telomere length (RTL) was longer in the colon tissue of individuals with adenoma compared to that of healthy individuals [4,5], suggesting longer RTL as a risk factor for adenoma and/or CRC. However, most studies comparing tumor tissue to paired adjacent normal tissue reveal shorter telomeres in CRC [6,7,8,9], although in some studies, no shortening observed [10].
Telomerase is the enzyme required for telomere maintenance [11] and is controlled by the telomerase reverse transcriptase (TERT) gene. TERT is downregulated in the majority of human tissues, but the telomerase RNA component (TERC) template is constitutively expressed [12]. Telomerase activity reappears in immortalized cell lines and in about 85% of human tumors, which has led to studies on the usefulness of telomerase for cancer diagnostics [13]. The exact mechanism of telomere length maintenance in tumor tissue is still undetermined but may be attributed to the genes related to telomerase activity, shelterin–telosome protein complex, histone binding, alternative telomere lengthening mechanism (ATL) or non-canonical telomere maintenance [14].
In this study, we examined the change in RTL in CRC tissue compared to corresponding normal colon tissue and its association with different histological and molecular changes of CRC. We also looked at differential gene expression of (a) the telomere-related genes directly or indirectly related to telomere function, (b) genes involved in selected cancer-related pathways and (c) genes at genome-wide level in CRC tissues to determine the association between gene expression and changes in RTL in paired (tumor–normal) tissues in CRC.

2. Materials and Methods

For this study, we included 330 paired (tumor and adjacent normal) samples from 165 CRC patients (m = 96, f = 69). Among them, 33 had right-sided CRC (cecum, 8; ascending colon, 17; hepatic flexure, 4; transverse colon, 4), and 132 had left-sided CRC (descending colon, 8; sigmoid colon, 15; rectosigmoid junction, 7; rectum, 102). A total 67 of the patients (39 male and 28 female) were aged ≤ 40 years. Fresh frozen samples were collected from 165 CRC patients from the department of Pathology, Bangabandhu Sheikh Mujib Medical University (BSMMU), Dhaka, Bangladesh, at different times spanning between December 2009 and May 2016. The patients were at different stages of CRC (Stage-1: 36; Stage-2: 42; Stage-3: 87). From each patient, the specimens were collected from the surgically resected tumor and the surrounding unaffected part of the colon about 5–10 cm away from the tumor mass. A surgical pathology fellow collected all samples from the operating room immediately after the surgical resection. Histopathology examination was performed on H&E-stained slides in routinely processed paraffin-impregnated tissue blocks. The slides were examined independently by two pathologists, and there was concordance in all 165 cases. For staging and grading of the CRC, the World Health Organization Classification of tumors was followed [15]. From each individual, we obtained a pair of tumor and normal tissues, which were frozen immediately and shipped on dry ice to the molecular genomics lab at the University of Chicago for subsequent DNA extraction and molecular assay.
For each patient, we also abstracted key demographic and clinical data and tumor characteristics from hospital medical records. Written informed consent was obtained from all participants. The research protocol was approved by the “Ethical Review Committee, Bangabandhu Sheikh Mujib Medical University”, Dhaka, Bangladesh (BSMMU/2010/10096), and by the “Biological Sciences Division, University of Chicago Hospital Institutional Review Board”, Chicago, IL, USA (10-264-E).

2.1. DNA and RNA Extraction and Quality Control

DNA was extracted from fresh frozen tissue using a Puregene Core kit (Qiagen, Maryland, USA). An electropherogram from Agilent BioAnalyzer with Agilent DNA 12000 chips showed the fragment size to be >10,000 bp. RNA was extracted from RNAlater preserved colonic tissue using a Ribopure tissue kit (Ambion, Austin, TX, USA, Cat# AM1924).

2.2. Relative Telomere Length (RTL) Measurement

For RTL measurement, we used a Luminex-based assay using QuantiGene Plex chemistry (Invitrogen, Fisher Scientific). The details of the assay were described previously [16,17]. Briefly, the assay requires ~50 ng of DNA, which is hybridized to sequence-specific probes for the telomere repeat sequence (TEL) and reference gene sequence (ALK). The TEL and ALK gene signals are amplified using branched DNA technology and detected using Luminex technology. We used custom-designed probes to measure abundance of the telomere repeat sequence. The 24-mer probe was targeted against 4 repeats—“TTAGGGTTAGGGTTAGGGTTAGGG”. As a reference single gene, we used ALK, which showed a very stable copy number (CN = 2) in all the DNA samples detected by oligonucleotide-based microarray SNP chips from our previous study. The result is a ratio; hence, there is no unit. The assay precision was good to excellent, with an intraclass correlation coefficient (ICC) of 0.91 (95% CI 0.86–0.94) [17]. The RTL assay failed in 17 samples (9 CRC and 8 normal tissue) out of the total 330 samples tested. This failure rate (5.1%) was similar to that reported in a previous large-scale study using the same Luminex-based RTL measurement assay [18].
Details of the assay and statistical analysis are presented in the Supplemental Methods.

3. Results

The different patient characteristics (shown in Table 1) suggested that the CRC was similar in male and female patients in terms of age group (age ≤ 40 years, vs. >40 years), location of the tumor, pathological staging, grade, KRAS mutation status and BRAF mutation status. However, MSI was slightly more frequent among the male participants (m = 30.2% vs. f = 17.4%, p = 0.06, Pearson chi-square test). In this study from Bangladesh (a southeast Asian country), in both male and female patients, in ~40% of the CRC cases were diagnosed at or below the age of 40 years. In our study, CRC in the ≤40 yr age group patients was more frequently associated with advanced stage at the time of diagnosis (65.7% Stage-3 in age ≤ 40 years compared to 43.9% Stage-3 in age > 40 years group, p = 0.006, Pearson chi-square test).
Comparison of RTL in CRC tissue and corresponding adjacent normal colonic tissue: The change in RTL (delta) in a patient was calculated as RTL in CRC tissue minus the RTL in the corresponding normal tissue. The frequency distributions of delta RTL, as well as paired RTL values (CRC and normal), are shown in Figure S1A,B, respectively. The RTL change in CRC tissue may also be shown as the ratio of RTL in CRC/RTL in normal tissue. However, the result was similar in respect to the biological interpretation (see Figure 1C and Figure S1D, where we also adjusted for age and sex in the ANOVA model). We report both the RTL values for CRC and normal tissue to demonstrate the magnitude of difference. On average, we observed relative telomere shortening in the CRC tissue compared to paired adjacent normal colon tissue from the same patient (RTL 0.722 ± SD 0.277 vs. 0.809 ± SD 0.242, p= 0.00012, paired t-test). This magnitude of telomere shortening (by ~0.087) in tumor tissue measured by Luminex assay is equivalent to telomere shortening in human leukocytes older than ~10 year of age [17]. In our previous study, in an independent set of DNA samples from 505 individuals between 21 years and 70 years of age (m = 329, f = 176), leukocyte RTL was measured by the same assay. Linear regression analysis suggested that with an increase in age by 10 years, the leukocyte RTL decreased by 0.08 (95% CI 0.06–0.09) [17].
In stratified analysis, we looked at the differences of RTL in CRC and normal tissue (RTL in CRC minus RTL in normal tissue) according to different clinico-morphological and molecular factors (shown in Table 2). A negative mean difference indicates shortening, and a positive value indicates lengthening of telomeres in tumor tissue. “Person-to-person variation” is an important factor that explains the variation in RTL (see Figure 1). We observed statistically significant shortening of telomeres in cancer tissue, irrespective of age group, gender, location (right/left sided), pathological diagnosis, staging, grading, mutational status (KRAS rs112445441 and BRAFv600E) or MSI status (Table 2). RTL shortening with aging was also seen for both normal and CRC tissue (ranging from 41 to 70 years of age) in this study. Essentially parallel lines were obtained for normal and CRC tissue (see Figure 2), confirming age-related RTL shortening. The shortening was not statistically significant in patients with Stage-1 disease but was significant in Stage-2 and Stage-3 cases. Stage-2 and Stage-3 may represent more chronic cases than Stage-1 and therefore an increased number of cell divisions in the cancer tissue and more telomere shortening.
Using ANOVA models including an interaction term (Tissue × Factor), we examined if the telomere shortening in tumor tissue was significantly different in presence or absence of any particular factor. The RTL was shorter in CRC irrespective of MSI status. However in the tumor tissue, there was significantly more telomere shortening in the patients with MSI compared to the patients with MSS tumor (see Figure 3A, ANOVA interaction term p = 0.0059). Similarly, telomere shortening in CRC tissue was greater in patients with low-grade tumor compared to those with high-grade tumor (see Figure 3B, ANOVA interaction term p = 0.0043). In the same line, the telomere shortening was more pronounced in Stage-2 and Stage-3, than it was in Stage-1, but the interaction p-values just crossed the conventional level of significance (ANOVA interaction p = 0.07).

Gene Expression of Telomere-Related Genes

In the next step, we explored the differential expression of genes related to telomere maintenance. We utilized our previous gene expression data (microarray, Illumina HT12 v4 chips) on the first 71 CRC cases and corresponding adjacent normal colonic mucosa (n = 70, one sample failed on microarray) from the same series of 165 cases for which we present the telomere data. Patient characteristics of the 71 CRC patients are shown in Table S1.
Depending on the difference of RTL (RTL in CRC vs. RTL in normal tissue), we divided the patients into two groups: (a) patients with delta values < 0, termed “patients with shortening of telomere” (n = 52); and (b) patients with delta >= 0, termed “no shortening of Telomere” (n = 19). One of the most well-known genes for telomerase activity is telomerase reverse transcriptase (TERT). The microarray we used had two probes in the TERT gene. One of the two probes showed 1.05-fold (95% CI 1.017–1.090, p = 0.0037) overexpression of TERT in CRC tissue compared to paired normal colonic mucosa (see Figure S2A). This slight degree of overexpression was not different between patients with “telomere shortening” and patients with “no telomere shortening” (see Figure S2B).
The complete list of telomere-related genes and paired comparisons of their expression in CRC and normal tissue is shown in Table S2. Considering the telomere length biology, we looked for four categories of genes:
  • Telomere maintenance genes: Among the patients with telomere shortening, gene set ANOVA analysis suggested that the “telomere maintenance” group of genes was, on average, 1.26-fold (95% CI 1.23–1.29) overexpressed in tumor tissue compared to corresponding normal colon tissue (p = 1.06 × 10−68). (Figure 4A). Even among patients without telomere shortening, these telomere maintenance genes were also overexpressed in tumor tissue compared to corresponding normal colon tissue, but by 1.17-fold (95% CI 1.13–1.1.22, p = 2.24 × 10−17) (Figure 4B). This difference in magnitude of overexpression was significantly higher in patients with “telomere shortening” (“tissue × TEL-shortening” interaction, p = 0.002). In other words, on average, the telomere maintenance genes were overexpressed in tumor tissue compared to paired normal colonic tissue, irrespective of the presence or absence of telomere shortening; however, the magnitude of overexpression was significantly higher if the patient had telomere shortening compared to those without telomere shortening.
2.
Alternative Lengthening of Telomere genes: A list of genes is shown in Table S2. Although statistically significant, the overall overexpression of this group of genes was only 1.03-fold (95% CI 1.01–1.05, p = 0.0012) in CRC tissue compared to normal tissue in patients with telomere shortening (see Figure S3A). Similar analysis for the same genes in patients with no telomere shortening did not show any statistically significant differential expression (fold change =1.03 (95% CI −1.004 to 1.06, p = 0.089)) in CRC compared to normal tissue (Figure S3B).
3.
Non-canonical Telomere maintenance genes: A list of genes is shown in Table S2. These genes were, on average, 1.03-fold (95% CI 1.01–1.05, p = 0.0022) overexpressed in tumor tissue compared to corresponding normal colon tissue in patients with “telomere shortening” (see Figure S4A). Similar analysis showed that these genes were also overexpressed to a similar extent by 1.03-fold (95% CI 1.01–1.06, p = 0.0052) in tumor tissue compared to corresponding normal colon tissue (see Figure S4B) in patients without telomere shortening. In other words, there was no difference in differential expression of this group of genes in the presence or absence of telomere shortening (interaction, p = 0.76).
4.
Shelterin–telomere protein-related genes: These genes were not differentially expressed in CRC compared to normal colon tissue in patients with “telomere shortening” (p = 0.19, Figure S5A). In patients with no telomere shortening, the fold change was minimal −1.03 (95% CI −1.001 to −1.06, p = 0.041, see Figure S5B).
Telomere shortening and selected group of cancer related genes: In light of well-known biological processes involved in carcinogenesis, we examined whether telomere shortening was associated with differential expression of some selected group of genes involved in different molecular mechanisms known to be altered in cancer, such as DNA repair, apoptosis, caspase initiator, caspase executor, tumor suppressor, antitumor suppressor, etc. (see Table S3 for list). Using the “tissue x TEL-shortening” interaction term in gene set ANOVA, we first detected the functional group(s) of genes where the magnitude of differential expression was different in patients with or without telomere shortening. Then, we looked at the paired analysis stratified by telomere-shortening status. We found that there was a difference of magnitude in differential expression of some functional groups of genes, such as “caspase executor” (interaction, p = 1.34 × 10−6); “antitumor suppressor” genes, such as cyclin-dependent kinase (CDK) (interaction, p = 0.003) and “DNA repair” (interaction, p = 0.01) in the presence or absence of telomere shortening.
Stratified paired analysis suggested that in patients with telomere shortening, CRC tissue exhibited more overexpression of the CDK group of genes (1.52-fold (95% CI 1.46–1.59), Figure 5A) than the overexpression seen in patients with no RTL shortening (1.35-fold (95% CI 1.26−1.43, Figure 5B). Similarly, the DNA repair genes were also more overexpressed in patients with telomere shortening (1.09-fold (95%CI 1.08–1.11) than in patients without telomere shortening (1.05-fold (95% CI 1.03–1.08). Considering the fact that (a) overexpression of CDK genes (antitumor suppressor) promotes tumorigenesis and (b) DNA damage in the tumor is associated with the cellular response of upregulation in DNA repair genes, the above-mentioned results suggest an underlying association of tumorigenesis (hence, increased cell division) and telomere shortening.
On the other hand, compared to paired normal tissue, CRC tissue in patients with telomere shortening had more downregulation of caspase executor genes (−1.24-fold (95% CI −1.20 to −1.28)) than the patients without telomere shortening (−1.06-fold (95%CI −1.12 to −1.01), Figure 6A,B, respectively). This suggests an association between impaired apoptosis (tumor growth) and telomere shortening seen in CRC.
Genome-wide Differential Gene Expression at “individual gene” level: In addition to the genes known to be related to telomere length or related to known cancer pathways, we also looked at all other genes at the genome-wide level that are differentially expressed in tumor tissue compared to paired normal colonic mucosa. By including an interaction term, “tissue × telomere shortening”, in ANOVA models, we tried to identify the individual genes, the differential expression of which (CRC vs. normal tissue) differs significantly in the presence or absence of telomere shortening. The p-value of the interaction term in ANOVA models suggested that there was a total of 3080 probes (covering 2683 genes), the differential expression (CRC vs. normal tissue) of which differed significantly in magnitude in the presence or absence of telomere shortening.
Gene ontology (GO) enrichment analysis of this list of genes is shown in Table S4. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway terms that were significantly enriched (enrichment p-value ≤ 0.05) within the list are shown in the Figure 7. The most significant enrichment was seen in genes related to “DNA replication”. This may suggest an association between DNA replication and telomere shortening.
Genome-wide Differential Gene Expression at “group of genes” level: In the next step, using gene set ANOVA, we tried to identify the “group of genes” (involved in different KEGG pathways), the overexpression or under-expression of which may be significantly different in magnitude in presence or absence of telomere shortening. We tested KEGG pathway terms (groups of genes). The top groups that showed difference in overall differential expression in the presence of telomere shortening are presented in Table S5. Interestingly, the most significant group of genes was “DNA replication” (interaction p-value = 1.2 × 10−11). The other groups include “retinol metabolism”, “nucleotide excision repair”, ”base excision repair”, “mismatch repair”, etc. When we stratified the paired analysis by telomere shortening, we found that, on average, the DNA replication genes were 1.26-fold (95% CI 1.24–1.27) overexpressed in CRC compared to paired normal colonic mucosa in the presence of telomere shortening, whereas the same group of genes was 1.16-fold (95% CI 1.14–1.18) overexpressed in CRC compared to paired normal colonic mucosa in patients without telomere shortening (Figure 8A,B, respectively). This may further emphasize the association of the DNA replication process with telomere shortening. The majority of the genes in the DNA replication group (57 out of 60) were over-expressed to some degree in CRC compared to paired normal colonic mucosa. Examples of DNA replication genes are shown in Figure S6A–E, which clearly shows that the magnitude of overexpression was greater in presence of telomere shortening.
Association of Telomere shortening, Tumor Grade and DNA replication: Our data showed that telomere shortening occurred more in low-grade tumors; we also found an association between upregulation of DNA replication genes and telomere shortening. Therefore, in the next step, using gene set ANOVA with interaction terms in the model, explored whether the overexpression of DNA replication genes was different in magnitude in low-grade and high-grade tumors. The interaction p-value for “tissue × grade” (p = 2.04 × 10−26) suggested that the DNA replication genes were more overexpressed in low-grade CRC than in high-grade CRC compared to paired normal samples (see Table 3). Given the fact that many of anticancer drugs work on proliferating cells, this gene expression finding may be consistent with the fact that anticancer drugs work better in low-grade than in high-grade CRC. We also included the interaction terms “tissue x telomere shortening” and “tissue × grade” in the same model as explanatory variables for the gene expression of functional groups (such as DNA replication) to explore whether the association between overexpression of DNA replication genes in CRC and telomere shortening was independent of tumor grade. The significant interaction p-values for both interaction terms suggest that the overexpression of DNA replication genes (presumably responsible for repeated cell division) was independently associated with telomere shortening in CRC tissue compared to paired normal colon tissue. When we looked for differential expression of telomere-related genes, we also noticed more overexpression of in low-grade CRC (see Table 3).
Association of Telomere shortening, Tumor Stage and DNA replication: Stratified analysis by stage (see Table 3) suggested that on average, DNA replication genes were significantly overexpressed in CRC tissue compared to paired normal tissue in all stages. However, there was significantly more overexpression in Stage-1 (Table 3). When we looked for differential expression of telomere-related genes in CRC, we found that the telomere maintenance genes were significantly more overexpressed in Stage-1 CRC. In other words, our data suggested that in Stage-1, there is more overexpression of DNA replication genes (which may result in more RTL shortening), but there is also more overexpression of telomere maintenance (which may compensate for the RTL shortening). In addition, chronicity may also play a role in the slightly more telomere shortening that occurs in Stage-2 and Stage-3 compared to Stage-1, although DNA replication genes were slightly more overexpressed in Stage-1.
Association of Telomere shortening, MSI status and DNA replication: Our data showed that telomere shortening was more pronounced in MSI tumors; we also found an association between upregulation of DNA replication genes and telomere shortening. Therefore, in the next step, we explored whether the overexpression of DNA replication genes was different in magnitude in MSS and MSI tumors. We found overexpression of DNA replication genes both in MSS and MSI tumors compared to paired normal tissue without any difference in magnitude (see Table 3). In other words, the difference in RTL shortening (more in MSI tumors) could not be attributed to a difference in upregulation of DNA replication genes. When we looked for differential expression of telomere maintenance genes stratified by MSI status, we found significant overexpression of telomere maintenance genes both in MSS and MSI tumors compared to paired normal tissue without any difference in magnitude, depending on MSI status (Table 3).
Association of Telomere shortening and Expression of Immune Target genes: Previously, we showed that MSI status could be used to identify CRC patients who would potentially benefit from particular immune checkpoint inhibitor drugs [19]. We explored whether the differential expression of such target genes in CRC tissue (compared to paired normal tissue) differs by telomere shortening status. Figure 9 shows a similar pattern of differential expression in CRC tissue with or without telomere shortening. There was no significant difference (interaction p = 0.66). This suggests that telomere shortening may not be helpful to identify patients who may potentially benefit from immune checkpoint inhibitors.
Association of Telomere shortening and Methylation status of Telomere-related genes: In addition to looking at the gene expression pattern, we also looked at the DNA methylation pattern. We utilized our previous genome-wide methylation data (microarray, Illumina Infinium methylation 450 K chips) on the first 125 CRC cases and corresponding adjacent normal colonic mucosa (n = 125) from the same series of 165 cases for which we presented the telomere data. In the methylation array we used, there was a total of 236 methylation probes covering the telomere-related genes. None of the probes showed differential methylation in tumor tissue compared to corresponding normal tissue exceeding 0.14 in either direction (hypo- or hypermethylation). The magnitude of differential methylation (beta of tumor vs. beta of corresponding normal) was low regardless of p-value (see Figure 10). We also used stratification by telomere shortening. The differential methylation of probes in (a) telomere maintenance genes (Figure S7), (b) alternate lengthening of telomere related genes (Figure S8), (c) non-canonical telomere maintenance genes (Figure S9) and (d) shelterin complex genes (Figure S10) are shown in additional figures. In general, telomere maintenance genes were very slightly hypomethylated in tumor tissue compared to corresponding normal colonic tissue both in the presence or absence of telomere shortening. This correlated also with the gene expression finding.

4. Discussion

A large-scale study including 31 cancer types suggests telomere shortening in a vast majority of cancers [20]. Cancer cells require telomere maintenance mechanisms for unlimited replicative potential. They may achieve this through TERT activation or alternative telomere lengthening associated with ATRX or DAXX loss. A recent large study involving 2500 matched tumor–control samples from 36 different tumor types suggests that whereas the telomere content of tumors with ATRX or DAXX mutations is increased, tumors with TERT modifications show a moderate decrease in telomere content [21]. It has also been proposed that telomere attrition could repress the outgrowth of cancer [22,23,24]. The role of telomere length in the etiology of CRC is not fully understood. A large study was conducted in Singapore enrolling more than 2600 subjects, among which 776 had CRC. Researchers measured telomere length from leucocytes using the qPCR method. Their findings support the notion that longer telomeres in leukocytes may be associated with a higher risk of CRC, particularly rectal cancer [25].
Fernandez-Rozadilla et al. compared telomere length of colon tissue from healthy controls with that of individuals with CRC, polyps and past or family history of polyps or CRC [4]. They found that the telomere length of normal mucosa of cases was higher than that of the controls. They found no correlation with location of tissue. Blood telomere length showed no association with case control status [4]. Peacock et al. examined the relationship between telomere length in normal colon tissue and the prevalence of colorectal adenoma, a precursor to CRC [5]. Forty cases with adenomas detected by colonoscopy and 45 controls were studied. RTL was measured using quantitative real-time PCR. RTL was significantly longer in colon tissue of individuals with adenomas compared to that of healthy individuals (p = 0.008) [5].
In some studies, investigators tested fresh frozen CRC tissue and corresponding adjacent normal colon mucosa and found a shorter telomere length in cancer tissue compared to the healthy tissue [8,9,26,27,28,29]. Suraweera et al. carried out a large study wherein they analyzed the elative telomere length of 90 adenoma, 419 CRC and their corresponding normal mucosa. Both adenoma and CRC tissue had significantly shorter RTL than the corresponding healthy mucosa [30]. They also found that female gender, proximal location and the TERT rs2736100 G allele were independently associated with longer age-adjusted RTL in normal mucosa of these patients. The age-adjusted RTL in normal and cancer tissue were independently associated with tumor stage. The authors did not find any association of tumor status of MSI, CIMP, TP53, KRAS and BRAF with RTL in normal or tumor tissue. However, they reported that near-tetraploid CRCs exhibited significantly longer RTL [30].
Rampazzo et al. studied 118 (CRC) samples (53 right colon, 30 left colon and 35 rectum) and corresponding adjacent non-cancerous tissues for telomere length, p53 mutation and MSI [9]. Telomeres were significantly shorter in CRCs than in adjacent tissues, regardless of tumor stage and grade, site or genetic alterations (p = 0.0001). Telomere length in CRCs did not differ with tumor progression or p53 status. In CRCs carrying the wild-type p53, telomeres were significantly shorter in tumors with MSI than in those with stable microsatellites (p = 0.027). Furthermore, telomere length differed according to tumor location, being longer in rectal cancers (p = 0.03).
O’Sullivan et al. measured telomere length by Q FISH and qPCR in colon adenoma [3]. Telomere length in larger colon adenoma lesions (>2 cm) was significantly shorter than that in normal adjacent (p = 0.004) or normal distant (p = 0.05) tissue from the same individuals. However, this association was not found in cases with colon adenoma of less than 2 cm in size.
A study was done on 42 MSS Stage-2 CRC tissue samples, where whole-exome sequencing data from paired normal colon and tumor tissue showed shorter telomere in tumor tissue [7]. The number of reads from the tumor sample was normalized to build the tumor telomere length ratio (TTLR). The functional gene set enrichment analysis showed pathways related to immune response (TCTN3, RNF7, ATP6VOE1, SNX3 and UBAP1) significantly associated with TLLR. However, telomere length had no correlation with gene expression changes of most of the genes responsible for telomere maintenance, TERT, RAP1, DKC1, TERF1, TERF2, POT1, TERF2IP and TPP1 [7].
Kroupa et al. [6] performed a large study with many important findings that are reproduced in our current study. They collected CRC tissue and adjacent healthy tissue from 721 individuals. RTL was measured by qPCR. RTL in tumor tissues was shorter than that in the adjacent mucosa (p < 0.0001). Shorter RTL was also observed in early-stage tumors compared to the advanced stages (p = 0.001). RTL was shorter in tumors at the proximal colon than at the distal colon (p < 0.0001). The authors compared paired normal and tumor tissue in patients with MSI and MSS, where patients with MSI had significantly shorter telomeres (p = 0.001). A similar finding was reported by Takagi et al. [31]. Telomere shortening was also shown in mucinous tumor histology (p < 0.0001) compared to tubular carcinoma [6].
One of the limitations of the current study is the lack of data on telomerase enzyme activity. The core components of telomerase are the reverse transcriptase (TERT) and telomerase RNA component (TERC), which provide the template for the synthesis of telomeric DNA. We looked at the differential expression of these genes. The role of telomerase in cancer was reviewed recently [32].Telomerase is associated with a set of accessory proteins, including dyskerin and nucleolar protein 10 (NOP10). We found that DKC1 (dyskerin pseudouridine synthase 1) was overexpressed in CRC. This is consistent with results reported by Nersisyan et al. [29], who studied gene expression in paired samples (tumor vs. healthy) in Lynch syndrome, sporadic CRC with MSI and MSS [29]. The DKC1 gene provides instructions for making a protein called dyskerin. This protein is involved in maintaining telomeric structure, which may be responsible for prolonged cell life. Alder et al. showed that mutation of the DKC1 gene was associated with early-age morbidity due to failure of telomere maintenance [33].
Telomerase is also known to be expressed in certain non-malignant cells, although generally at low levels. Engelhardt et al. investigated telomerase activity and telomere length in premalignant, malignant, inflammatory and normal colon specimens [34]. Telomerase activity was evaluated in 130 frozen specimens from human colon cancer (n= 50), adjacent normal colon tissue (n= 50), colon polyps (n = 20) and colitis (n = 10) using a modified telomeric repeat amplification protocol assay, whereas telomere length was assessed by Southern blot. High to moderate levels of telomerase activity were detected in 90% of colorectal tumors. Weakly positive activity was detected in 10%. None of the normal tissues exhibited telomerase activity. In polyps and colitis, telomerase activity was found but was 25- to 54-fold lower than that detected in colon cancer (p < 0.001). Late-stage tumors demonstrated increased telomerase activity compared to early-stage tumors [34].
Katayama et al. examined telomerase activity and telomere length in gastric and colon cancer [10]. Although telomerase activity was detected in 1/12 (8%) cases of gastric polyp and in 2/9 (22%) cases of colorectal polyps, its positivity in gastric cancer and CRC was 7/10 (70%) and 21/26 (81%; p < 0.0003 and p < 0.0001, respectively) [10].
Gertler et al. found significant positive correlations between telomere length and hTERT expression in both non-cancer colorectal mucosa (r 0.54; p 0.001) and CRC (r 0.52; p 0.001) [27].
Garcia-Aranda et al. studied 91 primary CRCs and their corresponding normal tissue samples for telomere length and telomerase activity [26]. Most tumors (81.3%) displayed telomerase activity. Telomeres in CRC specimens were significantly shorter compared to the normal, adjacent specimens (p = 0.02). Moreover, tumors that demonstrated shortened telomeres displayed higher TRF1 levels than tumors without telomere shortening.
Bergstrand et al. studied Hoyeraal–Hreidarsson syndrome (HHS), a multisystem disorder characterized by bone marrow failure, developmental defects and very short telomeres, which is caused by germline mutations in genes related to telomere biology [35]. WRAP53 is an essential component of the telomerase holoenzyme complex, a ribonucleoprotein complex required for telomere synthesis. A mutation of this gene is responsible for the HHS. We found the WDR79 (WRAP53) gene to be overexpressed in CRC tissue compared to normal tissue, and the overexpression was more pronounced in the presence of MSI, as well as in presence of RTL shortening (Figure S11).
Fu et al. found that overexpression of the NAT10 gene is associated with telomere shortening [36]. The shorter telomeres in the CRC tissue and NAT10 overexpression found in our study are consistent with this finding.
Considering the fact that the RTL varies from person to person, examination of paired samples (tumor and adjacent normal tissue) from the same patient may be the source of the most robust result with respect to an association between cancer and RTL. Our current study strongly suggests that in general, there is shortening of the telomere in CRC tissue compared to adjacent normal colonic tissue. However, in a few cases, there is lengthening of telomere as well. Our data suggest that the shortening of the telomere is the consequence of rapid cell division in CRC, which is supported by its significant association with upregulation of “DNA replication” genes, upregulation of “antitumor suppressor genes (cyclin dependent kinase genes)”, upregulation of “DNA repair” genes and downregulation of “caspase executor” (reducing apoptosis). Our data also suggest that cancer cells exhibit upregulation of some “telomere maintenance” genes, irrespective of telomere shortening, although the upregulation was greater in the presence of telomere shortening. We admit that we have data on neither telomerase activity nor somatic mutation of telomere-related genes. We also do not have clinical follow-up data to test the prognostic significance of telomere shortening in CRC. Telomere shortening was more noticeable in low-grade tumors (which was associated with more upregulation of “DNA replication” genes) and in MSI tumors (which was not related to upregulation of “DNA replication” genes). The current study cannot determine the cause or effect relationship of MSI and telomere shortening. Both may represent genomic instability. A recent study by Takagi et al. also found association between MSI and telomere shortening in CRC [31].
To our knowledge, this is the first study utilizing clinical samples from a southeast Asian country. In addition, this is one of the few studies incorporating genome-wide differential gene expression, as well as differential methylation data, from a subset of the same patients to study the association with RTL change in CRC. About 40% of the patients were ≤40 years of age, who had more severe disease than the patients of >40 years. Some of our findings support the previous findings in different populations. Unlike MSI, which may be used as a surrogate marker for upregulation of some immune target genes (hence, a potential candidate for certain checkpoint inhibitors), this telomere shortening may not be a marker for the ICI group of drugs but may be a potential marker of increased cell division and hence better candidates for traditional antineoplastic drugs working on proliferating cells in general. Future studies are needed to support this assumption.

5. Conclusions

In CRC, we found significant shortening of RTL in cancer tissue compared to the paired normal tissue, irrespective of age group, gender, tumor pathology, location and MSI status. RTL shortening was more frequent in low-grade CRC and in the presence of MSI. In CRC, RTL shortening is associated with (a) upregulation of DNA replication genes, cyclin dependent kinase genes (antitumor suppressor), and (b) downregulation of “caspase executor genes”, reducing apoptosis.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14092250/s1, Methods, statistical analysis [17,18,19,37,38,39,40,41,42,43,44,45,46,47,48,49,50], Figure S1: Histogram of Delta: RTL in CRC-RTL in corresponding Normal, Figure S2: Expression of TERT in CRC tissue, Figure S3: Alternative Lengthening of Telomere genes, Figure S4: Non-canonical Telomere maintenance genes, Figure S5: Shelterin–telomere protein-related genes, Figure S6: DNA replication genes, Figure S7: Telomere maintenance genes, Figure S8: Alternate lengthening of telomere related genes, Figure S9: Non-canonical telomere maintenance genes, Figure S10: shelterin complex genes, Figure S11: Gene expression of WDR79. Table S1: Patient characteristics of the 71 CRC patients with Gene expression data available, Table S2: TEL maintenanceGenes_GX, Table S3: _ANOVA_Paired_GX_Cancer, Table S4: Enrichment_KEGG_Intera, Table S5: GO-ANOVA_GX_KEGG_InteracTELshor.

Author Contributions

M.G.K. conceived and designed the study, developed the telomere assay, performed data analysis and wrote the manuscript. M.R. and R.P. collected the samples and carried out the histopathology. M.K. organized and supervised the tissue collection and was responsible for histopathology. Z.H. performed the surgery and contributed the tissue samples. A.M. helped in literature search and manuscript preparation. B.L.P. helped in manuscript preparation. H.A. supported and coordinated the study and helped in manuscript preparation. F.J. designed and developed the telomere assay, carried out the microarray assays and drafted the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partially supported by NIH funds P30ES027792 and R24ES028532.

Institutional Review Board Statement

This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board, Bangabandhu Sheikh Mujib Medical University”, Dhaka, Bangladesh (BSMMU/2010/10096), and by the “Biological Sciences Division, University of Chicago Hospital Institutional Review Board”, Chicago, IL, USA (10-264-E).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All supporting data are presented in the tables presented in the main manuscript and as additional supplemental material.

Acknowledgments

We acknowledge the support and help from all the patients included this study. We thank the University of Chicago Research, Bangladesh (URB) staff for handling and shipment of all the study material to the University of Chicago molecular genomics laboratory.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Blackburn, E.H. Structure and function of telomeres. Nature 1991, 350, 569–573. [Google Scholar] [CrossRef] [PubMed]
  2. De Lange, T. Shelterin-Mediated Telomere Protection. Annu. Rev. Genet. 2018, 52, 223–247. [Google Scholar] [CrossRef] [PubMed]
  3. O’Sullivan, R.J.; Karlseder, J. Telomeres: Protecting chromosomes against genome instability. Nat. Rev. Mol. Cell Biol. 2010, 11, 171–181. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Fernandez-Rozadilla, C.; Kartsonaki, C.; Woolley, C.; McClellan, M.; Whittington, D.; Horgan, G.; Leedham, S.; Kriaucionis, S.; East, J.; Tomlinson, I. Telomere length and genetics are independent colorectal tumour risk factors in an evaluation of biomarkers in normal bowel. Br. J. Cancer 2018, 118, 727–732. [Google Scholar] [CrossRef] [Green Version]
  5. Peacock, S.D.; Massey, T.E.; Vanner, S.J.; King, W.D. Telomere length in the colon is related to colorectal adenoma prevalence. PLoS ONE 2018, 13, e0205697. [Google Scholar] [CrossRef]
  6. Kroupa, M.; Rachakonda, S.K.; Liska, V.; Srinivas, N.; Urbanova, M.; Jiraskova, K.; Schneiderova, M.; Vycital, O.; Vymetalkova, V.; Vodickova, L.; et al. Relationship of telomere length in colorectal cancer patients with cancer phenotype and patient prognosis. Br. J. Cancer 2019, 121, 344–350. [Google Scholar] [CrossRef]
  7. Lopez-Doriga, A.; Valle, L.; Alonso, M.H.; Aussó, S.; Closa, A.; Sanjuan, X.; Barquero, D.; Rodríguez-Moranta, F.; Sanz-Pamplona, R.; Moreno, V. Telomere length alterations in microsatellite stable colorectal cancer and association with the immune response. Biochim. Biophys. Acta Mol. Basis Dis. 2018, 1864, 2992–3000. [Google Scholar] [CrossRef]
  8. Nakamura, K.; Furugori, E.; Esaki, Y.; Arai, T.; Sawabe, M.; Okayasu, I.; Fujiwara, M.; Kammori, M.; Mafune, K.; Kato, M.; et al. Correlation of telomere lengths in normal and cancers tissue in the large bowel. Cancer Lett. 2000, 158, 179–184. [Google Scholar] [CrossRef]
  9. Rampazzo, E.; Bertorelle, R.; Serra, L.; Terrin, L.; Candiotto, C.; Pucciarelli, S.; Del Bianco, P.; Nitti, D.; De Rossi, A. Relationship between telomere shortening, genetic instability, and site of tumour origin in colorectal cancers. Br. J. Cancer 2010, 102, 1300–1305. [Google Scholar] [CrossRef] [Green Version]
  10. Katayama, S.; Shiota, G.; Oshimura, M.; Kawasaki, H. Clinical usefulness of telomerase activity and telomere length in the preoperative diagnosis of gastric and colorectal cancer. J. Cancer Res. Clin. Oncol. 1999, 125, 405–410. [Google Scholar] [CrossRef]
  11. Palm, W.; de Lange, T. How shelterin protects mammalian telomeres. Annu. Rev. Genet. 2008, 42, 301–334. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Shay, J.W.; Bacchetti, S. A survey of telomerase activity in human cancer. Eur. J. Cancer 1997, 33, 787–791. [Google Scholar] [CrossRef]
  13. Nakamura, T.M.; Morin, G.B.; Chapman, K.B.; Weinrich, S.L.; Andrews, W.H.; Lingner, J.; Harley, C.B.; Cech, T.R. Telomerase catalytic subunit homologs from fission yeast and human. Science 1997, 277, 955–959. [Google Scholar] [CrossRef]
  14. Monteagudo, M.; Martínez, P.; Leandro-García, L.J.; Martínez-Montes, Á.M.; Calsina, B.; Pulgarín-Alfaro, M.; Díaz-Talavera, A.; Mellid, S.; Letón, R.; Gil, E.; et al. Analysis of Telomere Maintenance Related Genes Reveals NOP10 as a New Metastatic-Risk Marker in Pheochromocytoma/Paraganglioma. Cancers 2021, 13, 4758. [Google Scholar] [CrossRef] [PubMed]
  15. Hamilton, S.R. Tumours of the colon and rectum. In World Health Organization Classification of Tumours-Pathology and Genetics of Tumours of the Digestive System; International Agency for Research on Cancer (IARC): Lyon, France, 2000. [Google Scholar]
  16. Kibriya, M.G.; Jasmine, F.; Roy, S.; Ahsan, H.; Pierce, B. Measurement of telomere length: A new assay using QuantiGene chemistry on a Luminex platform. Cancer Epidemiol. Biomark. Prev. 2014, 23, 2667–2672. [Google Scholar] [CrossRef] [Green Version]
  17. Jasmine, F.; Shinkle, J.; Sabarinathan, M.; Ahsan, H.; Pierce, B.L.; Kibriya, M.G. A novel pooled-sample multiplex luminex assay for high-throughput measurement of relative telomere length. Am. J. Hum. Biol. 2018, 30, e23118. [Google Scholar] [CrossRef]
  18. Demanelis, K.; Jasmine, F.; Chen, L.S.; Chernoff, M.; Tong, L.; Delgado, D.; Zhang, C.; Shinkle, J.; Sabarinathan, M.; Lin, H.; et al. Determinants of telomere length across human tissues. Science 2020, 369, eaaz6876. [Google Scholar] [CrossRef]
  19. Jasmine, F.; Haq, Z.; Kamal, M.; Raza, M.; da Silva, G.; Gorospe, K.; Paul, R.; Strzempek, P.; Ahsan, H.; Kibriya, M.G. Interaction between Microsatellite Instability (MSI) and Tumor DNA Methylation in the Pathogenesis of Colorectal Carcinoma. Cancers 2021, 13, 4956. [Google Scholar] [CrossRef]
  20. Barthel, F.P.; Wei, W.; Tang, M.; Martinez-Ledesma, E.; Hu, X.; Amin, S.B.; Akdemir, K.C.; Seth, S.; Song, X.; Wang, Q.; et al. Systematic analysis of telomere length and somatic alterations in 31 cancer types. Nat. Genet. 2017, 49, 349–357. [Google Scholar] [CrossRef]
  21. Sieverling, L.; Hong, C.; Koser, S.D.; Ginsbach, P.; Kleinheinz, K.; Hutter, B.; Braun, D.M.; Cortés-Ciriano, I.; Xi, R.; Kabbe, R.; et al. Genomic footprints of activated telomere maintenance mechanisms in cancer. Nat. Commun. 2020, 11, 733. [Google Scholar] [CrossRef] [Green Version]
  22. Maciejowski, J.; de Lange, T. Telomeres in cancer: Tumour suppression and genome instability. Nat. Rev. Mol. Cell Biol. 2017, 18, 175–186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Schmutz, I.; Mensenkamp, A.R.; Takai, K.K.; Haadsma, M.; Spruijt, L.; de Voer, R.M.; Choo, S.S.; Lorbeer, F.K.; van Grinsven, E.J.; Hockemeyer, D.; et al. TINF2 is a haploinsufficient tumor suppressor that limits telomere length. eLife 2020, 9, e61235. [Google Scholar] [CrossRef] [PubMed]
  24. Günes, C.; Rudolph, K.L. Telomere dysfunction puts the brakes on oncogene-induced cancers. EMBO J. 2012, 31, 2833–2834. [Google Scholar] [CrossRef] [PubMed]
  25. Luu, H.N.; Qi, M.; Wang, R.; Adams-Haduch, J.; Miljkovic, I.; Opresko, P.L.; Jin, A.; Koh, W.P.; Yuan, J.M. Association Between Leukocyte Telomere Length and Colorectal Cancer Risk in the Singapore Chinese Health Study. Clin. Transl. Gastroenterol. 2019, 10, e00043. [Google Scholar] [CrossRef]
  26. Garcia-Aranda, C.; de Juan, C.; Diaz-Lopez, A.; Sanchez-Pernaute, A.; Torres, A.J.; Diaz-Rubio, E.; Balibrea, J.L.; Benito, M.; Iniesta, P. Correlations of telomere length, telomerase activity, and telomeric-repeat binding factor 1 expression in colorectal carcinoma. Cancer 2006, 106, 541–551. [Google Scholar] [CrossRef]
  27. Gertler, R.; Rosenberg, R.; Stricker, D.; Friederichs, J.; Hoos, A.; Werner, M.; Ulm, K.; Holzmann, B.; Nekarda, H.; Siewert, J.R. Telomere length and human telomerase reverse transcriptase expression as markers for progression and prognosis of colorectal carcinoma. J. Clin. Oncol. 2004, 22, 1807–1814. [Google Scholar] [CrossRef]
  28. Hastie, N.D.; Dempster, M.; Dunlop, M.G.; Thompson, A.M.; Green, D.K.; Allshire, R.C. Telomere reduction in human colorectal carcinoma and with ageing. Nature 1990, 346, 866–868. [Google Scholar] [CrossRef]
  29. Nersisyan, L.; Hopp, L.; Loeffler-Wirth, H.; Galle, J.; Loeffler, M.; Arakelyan, A.; Binder, H. Telomere Length Maintenance and Its Transcriptional Regulation in Lynch Syndrome and Sporadic Colorectal Carcinoma. Front. Oncol. 2019, 9, 1172. [Google Scholar] [CrossRef] [Green Version]
  30. Suraweera, N.; Mouradov, D.; Li, S.; Jorissen, R.N.; Hampson, D.; Ghosh, A.; Sengupta, N.; Thaha, M.; Ahmed, S.; Kirwan, M.; et al. Relative telomere lengths in tumor and normal mucosa are related to disease progression and chromosome instability profiles in colorectal cancer. Oncotarget 2016, 7, 36474–36488. [Google Scholar] [CrossRef] [Green Version]
  31. Takagi, S.; Kinouchi, Y.; Hiwatashi, N.; Nagashima, F.; Chida, M.; Takahashi, S.; Negoro, K.; Shimosegawa, T.; Toyota, T. Relationship between microsatellite instability and telomere shortening in colorectal cancer. Dis. Colon Rectum 2000, 43, S12–S17. [Google Scholar] [CrossRef]
  32. Robinson, N.J.; Schiemann, W.P. Telomerase in Cancer: Function, Regulation, and Clinical Translation. Cancers 2022, 14, 808. [Google Scholar] [CrossRef] [PubMed]
  33. Alder, J.K.; Parry, E.M.; Yegnasubramanian, S.; Wagner, C.L.; Lieblich, L.M.; Auerbach, R.; Auerbach, A.D.; Wheelan, S.J.; Armanios, M. Telomere phenotypes in females with heterozygous mutations in the dyskeratosis congenita 1 (DKC1) gene. Hum. Mutat. 2013, 34, 1481–1485. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Engelhardt, M.; Drullinsky, P.; Guillem, J.; Moore, M.A. Telomerase and telomere length in the development and progression of premalignant lesions to colorectal cancer. Clin. Cancer Res. 1997, 3, 1931–1941. [Google Scholar] [PubMed]
  35. Bergstrand, S.; Böhm, S.; Malmgren, H.; Norberg, A.; Sundin, M.; Nordgren, A.; Farnebo, M. Biallelic mutations in WRAP53 result in dysfunctional telomeres, Cajal bodies and DNA repair, thereby causing Hoyeraal-Hreidarsson syndrome. Cell Death Dis. 2020, 11, 238. [Google Scholar] [CrossRef] [PubMed]
  36. Fu, D.; Collins, K. Purification of human telomerase complexes identifies factors involved in telomerase biogenesis and telomere length regulation. Mol. Cell 2007, 28, 773–785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Cawthon, R.M. Telomere measurement by quantitative PCR. Nucleic Acids Res. 2002, 30, e47. [Google Scholar] [CrossRef] [PubMed]
  38. Cawthon, R.M. Telomere length measurement by a novel monochrome multiplex quantitative PCR method. Nucleic Acids Res. 2009, 37, e21. [Google Scholar] [CrossRef] [Green Version]
  39. Kibriya, M.G.; Jasmine, F.; Roy, S.; Ahsan, H.; Pierce, B.L. Novel Luminex Assay for Telomere Repeat Mass Does Not Show Well Position Effects Like qPCR. PLoS ONE 2016, 11, e0155548. [Google Scholar] [CrossRef]
  40. Harley, C.B.; Futcher, A.B.; Greider, C.W. Telomeres shorten during ageing of human fibroblasts. Nature 1990, 345, 458–460. [Google Scholar] [CrossRef]
  41. Kimura, M.; Stone, R.C.; Hunt, S.C.; Skurnick, J.; Lu, X.; Cao, X.; Harley, C.B.; Aviv, A. Measurement of telomere length by the Southern blot analysis of terminal restriction fragment lengths. Nat. Protoc. 2010, 5, 1596–1607. [Google Scholar] [CrossRef]
  42. Pierce, B.L.; Jasmine, F.; Roy, S.; Zhang, C.; Aviv, A.; Hunt, S.C.; Ahsan, H.; Kibriya, M.G. Telomere length measurement by a novel Luminex-based assay: A blinded comparison to Southern blot. Int. J. Mol. Epidemiol. Genet 2016, 7, 18–23. [Google Scholar] [PubMed]
  43. Eisenberg, D.T. Telomere length measurement validity: The coefficient of variation is invalid and cannot be used to compare quantitative polymerase chain reaction and Southern blot telomere length measurement techniques. Int. J. Epidemiol. 2016, 45, 1295–1298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Eisenberg, D.T.; Kuzawa, C.W.; Hayes, M.G. Improving qPCR telomere length assays: Controlling for well position effects increases statistical power. Am. J. Hum. Biol. 2015, 27, 570–575. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Janavicius, R.; Matiukaite, D.; Jakubauskas, A.; Griskevicius, L. Microsatellite instability detection by high-resolution melting analysis. Clin. Chem. 2010, 56, 1750–1757. [Google Scholar] [CrossRef] [Green Version]
  46. Kibriya, M.G.; Raza, M.; Jasmine, F.; Roy, S.; Paul-Brutus, R.; Rahaman, R.; Dodsworth, C.; Rakibuz-Zaman, M.; Kamal, M.; Ahsan, H. A genome-wide DNA methylation study in colorectal carcinoma. BMC Med. Genom. 2011, 4, 50. [Google Scholar] [CrossRef] [Green Version]
  47. Deschoolmeester, V.; Baay, M.; Wuyts, W.; Van Marck, E.; Van Damme, N.; Vermeulen, P.; Lukaszuk, K.; Lardon, F.; Vermorken, J.B. Detection of microsatellite instability in colorectal cancer using an alternative multiplex assay of quasi-monomorphic mononucleotide markers. J. Mol. Diagn. 2008, 10, 154–159. [Google Scholar] [CrossRef] [Green Version]
  48. Findeisen, P.; Kloor, M.; Merx, S.; Sutter, C.; Woerner, S.M.; Dostmann, N.; Benner, A.; Dondog, B.; Pawlita, M.; Dippold, W.; et al. T25 repeat in the 3′ untranslated region of the CASP2 gene: A sensitive and specific marker for microsatellite instability in colorectal cancer. Cancer Res. 2005, 65, 8072–8078. [Google Scholar] [CrossRef] [Green Version]
  49. Eisenhart, C. The assumptions underlying the analysis of variance. Biometrics 1947, 3, 1–21. [Google Scholar] [CrossRef]
  50. Downey, T. Analysis of a multifactor microarray study using Partek genomics solution. Methods Enzym. 2006, 411, 256–270. [Google Scholar] [CrossRef]
Figure 1. Variation of RTL in the data. The source of variations in the RTL that can be explained in the ANOVA models are shown. The upper row (A,B) shows unpaired analysis, and the lower row (C,D) shows result from paired analysis (taking the “Person ID” in the ANOVA model). The mean F-ratio (F-statistics for the factor/F-statistics for the model error) representing the significance of the factor in the ANOVA model is shown in the bar graphs (A,C). The sums of squares in the ANOVA model representing the proportion of the variation explained by the factors are shown as pie charts (B,D). When the “Person ID” is not included in the ANOVA model (unpaired analysis shown in B), 96% of the variation in RTL cannot be explained (error, 96.28). However, when the “Person ID” is incorporated into the ANOVA model (paired analysis shown in D), the error reduces to 25%, and the person-to-person variation explains 71.93% of the variation in RTL.
Figure 1. Variation of RTL in the data. The source of variations in the RTL that can be explained in the ANOVA models are shown. The upper row (A,B) shows unpaired analysis, and the lower row (C,D) shows result from paired analysis (taking the “Person ID” in the ANOVA model). The mean F-ratio (F-statistics for the factor/F-statistics for the model error) representing the significance of the factor in the ANOVA model is shown in the bar graphs (A,C). The sums of squares in the ANOVA model representing the proportion of the variation explained by the factors are shown as pie charts (B,D). When the “Person ID” is not included in the ANOVA model (unpaired analysis shown in B), 96% of the variation in RTL cannot be explained (error, 96.28). However, when the “Person ID” is incorporated into the ANOVA model (paired analysis shown in D), the error reduces to 25%, and the person-to-person variation explains 71.93% of the variation in RTL.
Cancers 14 02250 g001
Figure 2. Comparison of RTL in CRC tissue (in blue) and corresponding normal tissue (in red) by age group (A), location of the tumor (B), tumor stage (C) and gender (D).
Figure 2. Comparison of RTL in CRC tissue (in blue) and corresponding normal tissue (in red) by age group (A), location of the tumor (B), tumor stage (C) and gender (D).
Cancers 14 02250 g002
Figure 3. Comparison of RTL in CRC tissue (in blue) and corresponding normal tissue (in red) by MSI status (A) and by tumor grade (B). RTL shortening in CRC tissue was greater in MSI tumors, as well as in low-grade tumors.
Figure 3. Comparison of RTL in CRC tissue (in blue) and corresponding normal tissue (in red) by MSI status (A) and by tumor grade (B). RTL shortening in CRC tissue was greater in MSI tumors, as well as in low-grade tumors.
Cancers 14 02250 g003
Figure 4. Differential gene expression of telomere maintenance genes in CRC tissue (in blue) compared to corresponding normal tissue (in red). The gene probes (multiple probes for some genes) are shown on the x-axis, and the Log2-transformed expression value is shown on the y-axis. A difference of 1 on the y-axis represents a two-fold overexpression. Data from patients with RTL shortening are shown on the left (A), and data from patients without RTL shortening are shown on the right (B). The average magnitude of overexpression was significantly higher (p = 0.002) if the patient had telomere shortening (1.26 fold (95% CI 1.23–1.29)) compared to those without telomere shortening (1.17 fold (95% CI 1.13–1.1.22)).
Figure 4. Differential gene expression of telomere maintenance genes in CRC tissue (in blue) compared to corresponding normal tissue (in red). The gene probes (multiple probes for some genes) are shown on the x-axis, and the Log2-transformed expression value is shown on the y-axis. A difference of 1 on the y-axis represents a two-fold overexpression. Data from patients with RTL shortening are shown on the left (A), and data from patients without RTL shortening are shown on the right (B). The average magnitude of overexpression was significantly higher (p = 0.002) if the patient had telomere shortening (1.26 fold (95% CI 1.23–1.29)) compared to those without telomere shortening (1.17 fold (95% CI 1.13–1.1.22)).
Cancers 14 02250 g004
Figure 5. Differential gene expression of antitumor suppressor genes (cyclin D kinase) in CRC tissue (in blue) compared to corresponding normal tissue (in red). The gene probes (multiple probes for some genes) are shown on the x-axis, and the Log2-transformed expression value is shown on the y-axis. CDC2 is also known as CDK-1. Data from patients with RTL shortening are shown on the left (A), and data from patients without RTL shortening are shown on the right (B). The average magnitude of overexpression was significantly higher (p = 0.003) if the patient had telomere shortening (1.52-fold (95% CI 1.46–1.59), shown in A) compared to those without telomere shortening (1.35-fold (95%CI 1.26−1.43), shown in B).
Figure 5. Differential gene expression of antitumor suppressor genes (cyclin D kinase) in CRC tissue (in blue) compared to corresponding normal tissue (in red). The gene probes (multiple probes for some genes) are shown on the x-axis, and the Log2-transformed expression value is shown on the y-axis. CDC2 is also known as CDK-1. Data from patients with RTL shortening are shown on the left (A), and data from patients without RTL shortening are shown on the right (B). The average magnitude of overexpression was significantly higher (p = 0.003) if the patient had telomere shortening (1.52-fold (95% CI 1.46–1.59), shown in A) compared to those without telomere shortening (1.35-fold (95%CI 1.26−1.43), shown in B).
Cancers 14 02250 g005
Figure 6. Differential gene expression of caspase executor genes in CRC tissue (in blue) compared to corresponding normal tissue (in red). The gene probes (multiple probes for some genes) are shown on the x-axis, and the Log2-transformed expression value is shown on the y-axis. Data from patients with RTL shortening are shown on the left (A), and data from patients without RTL shortening are shown on the right (B). The average magnitude of downregulation was significantly more (p = 1.34 × 10−6) if the patient had telomere shortening (−1.24-fold (95% CI −1.20 to −1.28), shown in A) compared to those without telomere shortening (−1.06-fold (95%CI −1.12 to −1.01), shown in B).
Figure 6. Differential gene expression of caspase executor genes in CRC tissue (in blue) compared to corresponding normal tissue (in red). The gene probes (multiple probes for some genes) are shown on the x-axis, and the Log2-transformed expression value is shown on the y-axis. Data from patients with RTL shortening are shown on the left (A), and data from patients without RTL shortening are shown on the right (B). The average magnitude of downregulation was significantly more (p = 1.34 × 10−6) if the patient had telomere shortening (−1.24-fold (95% CI −1.20 to −1.28), shown in A) compared to those without telomere shortening (−1.06-fold (95%CI −1.12 to −1.01), shown in B).
Cancers 14 02250 g006
Figure 7. Enrichment score of the list of genes that had a significantly different magnitude of differential expression in CRC tissue than normal tissue, depending on the presence or absence of MSI.
Figure 7. Enrichment score of the list of genes that had a significantly different magnitude of differential expression in CRC tissue than normal tissue, depending on the presence or absence of MSI.
Cancers 14 02250 g007
Figure 8. Differential expression of genes involved in “DNA replication” in CRC tissue (in blue) compared to corresponding normal tissue (in red). The gene probes (multiple probes for some genes) are shown on the x-axis, and the Log2-transformed expression value is shown on the y-axis. Data from patients with RTL shortening are shown on the left (A), and data from patients without RTL shortening are shown on the right (B).
Figure 8. Differential expression of genes involved in “DNA replication” in CRC tissue (in blue) compared to corresponding normal tissue (in red). The gene probes (multiple probes for some genes) are shown on the x-axis, and the Log2-transformed expression value is shown on the y-axis. Data from patients with RTL shortening are shown on the left (A), and data from patients without RTL shortening are shown on the right (B).
Cancers 14 02250 g008
Figure 9. Immune checkpoint inhibitors (ICIs) are targeted to certain immunotarget genes. Differential expression of those immunotarget genes in CRC tissue (in blue) compared to corresponding normal tissue (in red) is shown. Data from patients with RTL shortening are shown on the left (A), and data from patients without RTL shortening are shown on the right (B). The magnitude of differential expression was similar in patients with or without RTL shortening.
Figure 9. Immune checkpoint inhibitors (ICIs) are targeted to certain immunotarget genes. Differential expression of those immunotarget genes in CRC tissue (in blue) compared to corresponding normal tissue (in red) is shown. Data from patients with RTL shortening are shown on the left (A), and data from patients without RTL shortening are shown on the right (B). The magnitude of differential expression was similar in patients with or without RTL shortening.
Cancers 14 02250 g009
Figure 10. Volcano plot showing the probe-level differential methylation data of telomere-related genes in CRC tissue compared to corresponding normal tissue. Magnitude of difference (beta value of CRC minus beta value of normal tissue) is shown on the x-axis, and the p-value from paired t-test is shown on the y-axis. The location of the probes in relation to the CpG island are color-coded.
Figure 10. Volcano plot showing the probe-level differential methylation data of telomere-related genes in CRC tissue compared to corresponding normal tissue. Magnitude of difference (beta value of CRC minus beta value of normal tissue) is shown on the x-axis, and the p-value from paired t-test is shown on the y-axis. The location of the probes in relation to the CpG island are color-coded.
Cancers 14 02250 g010
Table 1. Characteristics of the 165 CRC patients.
Table 1. Characteristics of the 165 CRC patients.
CharacteristicCategoryMale (n = 96)Female (n = 69)p-Value
Age (years)mean45.644.860.773
(SD)13.8213.98
Age category≤40 years39280.582
>40 years5741
LocationLeft75570.556
Right2112
HistopathologyAdenocarcinoma78630.171
Mucinous adenocarcinoma176
Squamous cell carcinoma1
StageStage-121150.987
Stage-22418
Stage-35136
GradeLow43350.528
High5334
MicrosatelliteMSI29120.069
MSS6757
CEA (ng/mL)mean42.5945.490.844
(SD)95.1486.9
KRASWild65500.435
Mutant2919
BRAFV600EWild87640.327
Mutant65
Table 2. RTL in CRC and paired normal tissue.
Table 2. RTL in CRC and paired normal tissue.
CharacteristicnRTL in CRC TissueRTL in Normal TissueDelta = Tumor-Normalp Value
MeanSDMeanSDMean95% Lower BOUND95% Upper Bound
Gender
Male880.73±0.270.81±0.24−0.08−0.14−0.020.009
Female640.71±0.290.81±0.25−0.09−0.16−0.030.004
Age
≥40910.72±0.260.81±0.23−0.09−0.15−0.04<0.001
<40610.73±0.30.81±0.26−0.08−0.160.0020.055
KRAS
Mutant450.67±0.310.80±0.24−0.12−0.21−0.040.007
Wild1060.75±0.260.82±0.24−0.07−0.12−0.020.007
BRAF
Mutant110.69±0.210.820.18−0.13−0.330.070.190
Wild1390.73±0.280.81±0.24−0.08−0.13−0.04<0.001
MSI Status
MSI380.67±0.260.85±0.28−0.18−0.26−0.10<0.001
MSS1140.74±0.280.80±0.23−0.06−0.11−0.0060.029
Histopathological Diagnosis
Adenocarcinoma1300.72±0.280.81±0.24−0.09−0.14−0.05<0.001
Mucinous Adenocarcinoma210.75±0.260.81±0.28−0.06−0.200.090.444
Squamous Cell Carcinoma10.97 0.84
Stage
Stage-1340.74±0.340.76±0.21−0.02−0.140.090.692
Stage-2370.70±0.270.82±0.29−0.12−0.19−0.050.001
Stage-3810.73±0.250.82±0.23−0.10−0.16−0.040.002
Grade
High-grade780.75±0.260.79±0.24−0.03−0.080.030.313
Low-grade740.69±0.290.84±0.25−0.15−0.21−0.08<0.001
Location
Left1220.74±0.290.82±0.25−0.08−0.13−0.030.003
Right300.65±0.240.28±0.20−0.13−0.02−0.030.009
Table 3. Overexpression of DNA replication genes and telomere maintenance genes in CRC compared to paired normal tissue stratified by tumor grade, tumor stage and MSI status.
Table 3. Overexpression of DNA replication genes and telomere maintenance genes in CRC compared to paired normal tissue stratified by tumor grade, tumor stage and MSI status.
StratificationDNA Replication GenesTelomere Maintenance Genes
Fold Change (95%CI)Fold Change (95%CI)
Grade
Low-Grade1.27(1.25–1.28)1.28(1.25–1.31)
High-Grade1.11(1.09–1.14)1.11(1.06–1.16)
p2.04 × 10−26 2.97 × 10−8
Stage
Stage-11.31(1.28–1.34)1.42(1.35–1.49)
Stage-21.22(1.20–1.25)1.19(1.15–1.24)
Stage-31.22(1.20–1.23)1.22(1.18–1.25)
p8.58 × 10−7 5.54 × 10−8
MSI status
MSS1.23(1.22–1.25)1.25(1.22–1.28)
MSI1.23(1.21–1.26)1.22(1.18–1.26)
p0.801001 0.34
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kibriya, M.G.; Raza, M.; Kamal, M.; Haq, Z.; Paul, R.; Mareczko, A.; Pierce, B.L.; Ahsan, H.; Jasmine, F. Relative Telomere Length Change in Colorectal Carcinoma and Its Association with Tumor Characteristics, Gene Expression and Microsatellite Instability. Cancers 2022, 14, 2250. https://doi.org/10.3390/cancers14092250

AMA Style

Kibriya MG, Raza M, Kamal M, Haq Z, Paul R, Mareczko A, Pierce BL, Ahsan H, Jasmine F. Relative Telomere Length Change in Colorectal Carcinoma and Its Association with Tumor Characteristics, Gene Expression and Microsatellite Instability. Cancers. 2022; 14(9):2250. https://doi.org/10.3390/cancers14092250

Chicago/Turabian Style

Kibriya, Muhammad G., Maruf Raza, Mohammed Kamal, Zahidul Haq, Rupash Paul, Andrew Mareczko, Brandon L. Pierce, Habibul Ahsan, and Farzana Jasmine. 2022. "Relative Telomere Length Change in Colorectal Carcinoma and Its Association with Tumor Characteristics, Gene Expression and Microsatellite Instability" Cancers 14, no. 9: 2250. https://doi.org/10.3390/cancers14092250

APA Style

Kibriya, M. G., Raza, M., Kamal, M., Haq, Z., Paul, R., Mareczko, A., Pierce, B. L., Ahsan, H., & Jasmine, F. (2022). Relative Telomere Length Change in Colorectal Carcinoma and Its Association with Tumor Characteristics, Gene Expression and Microsatellite Instability. Cancers, 14(9), 2250. https://doi.org/10.3390/cancers14092250

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop