Abstract
Objective
We aimed to evaluate the relationship between telomere length and systemic lupus erythematosus (SLE).
Methods
PUBMED and EMBASE databases were searched; meta-analyses were performed comparing telomere length in SLE patients and healthy controls, and on SLE patients in subgroups based on ethnicity, sample type, assay method and data type.
Results
Eight studies including 472 SLE patients and 365 controls were ultimately selected which showed that telomere length was significantly shorter in the SLE group than in the control group (standardized mean difference (SMD) = −0.835, 95% confidence interval (CI) = −1.291 to −0.380, p = 3.3 × 10−4). Stratification by ethnicity showed significantly shortened telomere length in the SLE group in Caucasian, Asian and mixed populations (SMD = −0.455, 95% CI = −0.763 to −0.147, p = 0.004; SMD = −0.887, 95% CI = −1.261 to −0.513, p = 3.4 × 10−4; SMD = −0.535, 95% CI = −0.923 to −0.147, p = 0.007; respectively). Furthermore, telomere length was significantly shorter in the SLE group than in the control group in whole blood and peripheral blood mononuclear cell groups (SMD = −0.361, 95% CI = −0.553 to −0.169, p = 2.3 × 10−4; SMD = −1.546, 95% CI = −2.583 to −0.510, p = 0.003; respectively); a similar trend was observed in leukocyte groups (SMD = −0.699, 95% CI = −1.511 to −0.114, p = 0.092). Meta-analyses based on assay method or data type revealed similar associations.
Conclusions
Our meta-analysis demonstrated that telomere length was significantly shorter in patients with SLE, regardless of ethnicity, sample type or assay method evaluated.
Introduction
Systemic lupus erythematosus (SLE) is a prototypic autoimmune disease, characterized by B cell hyperactivity, immune-complex deposition, high level of autoantibody production and multiple organ damage. Although the etiology of SLE has been known to be multifactorial, caused by interactions between genetic and environmental factors, findings including significant familial aggregation and multiple genetic linkages in SLE demonstrate that the etiology of SLE is genetic in nature. 1
Telomeres are DNA–protein complexes composed of TTAGGG-rich repeats capped at the ends of chromosomes in eukaryotic cells, and play an essential role in maintaining the stability and integrity of chromosomal architecture during duplication. 2 Telomeres shorten during cell divisions by loss of the hexameric repeats, resulting in shortening of telomeres. 3 Thus, telomere length reflects the replicative history of cells, and can be considered as a mitotic or biological clock. 3 Environmental factors such as inflammation and oxidative stress are associated with accelerated telomere shortening. 4 If cells reach a critically short size of telomere length, cellular senescence or apoptosis is triggered. 5 Telomere shortening has been involved in genetic instability of cells that are suspected to be implicated in the pathogenesis of various diseases, including autoimmune diseases, cancer, type 2 diabetes, neuropsychiatric disorders, as a possible risk factor.6–9
Studies have shown that shortened telomere length is associated with SLE, but, on the other hand, other reports found no such associations.10–17 These disparities are probably caused by small sample sizes, low statistical power, and/or clinical heterogeneity; therefore, to overcome the limitations of individual studies, resolve inconsistencies and reduce the likelihood that random errors are responsible for false-positive or false-negative associations, we turned to meta-analysis.18-20 The aim of the present study was to determine the relationship between telomere length and SLE by using a meta-analytic approach.
Materials and methods
Identification of eligible studies and data extraction
We performed a literature search for studies that examined telomere length in SLE patients and controls. PUBMED and EMBASE databases were searched to identify all available past articles (up to February 2016). The following keywords and subject terms were used in the search: ‘telomere’, ‘systemic lupus erythematosus’ and ‘SLE’. All references cited were also reviewed to identify additional studies not covered by the above-mentioned electronic databases. Studies were considered eligible if (1) they were case-control, cross-sectional or cohort studies; (2) they provided data on telomere length in case and control groups; and (3) they provided sufficient data on the association between telomere length and SLE. Studies were excluded if (1) they contained overlapping or insufficient data or (2) they were reviews or case reports. The following information was extracted from each study: primary author, year of publication, country, ethnicity, number of participants, sample type, assay method, and mean and standard deviation (SD) of telomere length. If the standard error of mean was reported, we calculated the SD using the appropriate statistical formula. When the data were in terms of median, range or p-value, we computed the mean and SD using previously described formulae.21,22 We conducted a sensitivity test on imputed values.
Evaluation of statistical associations
We performed a meta-analysis examining the relationship between telomere length and SLE. For continuity of data, results were presented as standardized mean differences (SMDs) and 95% confidence intervals (CIs). SMDs were calculated by dividing the mean difference between two groups by the pooled SD, and were used when different scales were integrated to measure the same concept. This measure compares case and control arms in terms of standardized scores. The magnitude of SMD was considered as follows: 0.2–0.5, small effect; 0.5–0.8, medium effect; ≥ 0.8, large effect. 23 We assessed within-study and between-study variations and heterogeneities using Cochran’s Q-statistics. 24 The heterogeneity test was used to assess the null hypothesis that all studies were evaluating the same effect. When the significant Q-statistic (p < 0.10) indicated heterogeneity across studies, the random effects model was used for the meta-analysis. 25 If not, the fixed effects model was used, which assumed that all studies estimated the same underlying effect, and considered within-study variation only. 24 We quantified the effect of heterogeneity using I2 = 100% × (Q–df)/Q, 26 where I2 measured the degree of inconsistency between studies and determined whether the percentage total variation across studies was due to heterogeneity rather than due to chance. I2 ranged between 0% and 100%; I2 values of 25%, 50% and 75% were referred to as low, moderate and high estimates, respectively. 26 Statistical manipulations were undertaken using the Comprehensive Meta-Analysis computer program (Biostat Inc., Englewood, NJ, USA).
Evaluation of heterogeneity, sensitivity test and publication bias
To examine potential sources of heterogeneity observed in the meta-analysis, meta-regression analysis was performed using the following variables: ethnicity, sample type, assay method, publication year, sample size and data type. A sensitivity test was performed to assess the influence of each individual study on the pooled odds ratio (OR) by omitting each study individually. Although funnel plots are often used to detect publication bias, they require diverse study types of varying sample sizes, and their interpretation involves subjective judgment. Therefore, we evaluated publication bias using Egger’s linear regression test, 27 which measured funnel plot asymmetry using a natural logarithm scale of ORs. When asymmetry was indicated, we used the ‘trim and fill’ method to adjust summary estimates for observed bias. 28 This method removes small studies until funnel plot symmetry is achieved by recalculating the center of the funnel before removed studies are replaced with their missing mirror-image counterparts. A revised summary estimate was then calculated using all original studies and hypothetical ‘filled’ studies.
Results
Studies included in the meta-analysis
Characteristics of individual studies included in the meta-analysis
SLE: systemic lupus erythematosus; SMD: standard mean difference; WB: whole blood; PBMC: peripheral blood mononuclear cell; qPCR: quantitative polymerase chain reaction; TRF: terminal restriction fragment; FISH: fluorescence in situ hybridization.
Magnitude of Cohen’s d effect size: 0.2–0.5: small effect; 0.5–0.8: medium effect; ≥0.8: large effect.
Not age-matched, but no statistical difference in age between SLE and control groups.
Meta-analysis of telomere length in SLE patients compared with controls
Meta-analysis was performed on all SLE patients, and on SLE patients in each ethnic group. Telomere length was significantly shorter in the SLE group than in the control group (SMD = −0.835, 95% CI = −1.291 to −0.380, p = 3.3 × 10−4) (Table 2, Figure 1). In addition, stratification by ethnicity showed a significantly shortened telomere length in SLE group in Caucasian, Asian and mixed populations (SMD = −0.455, 95% CI = −0.763 to −0.147, p = 0.004; SMD = −0.887, 95% CI = −1.261 to −0.513, p = 3.4 × 10−4; SMD = −0.535, 95% CI =−0.923 to −0.147, p = 0.007; respectively) (Table 2, Figure 2). A single African-American study showed a significantly shortened telomere length in SLE group (SMD = −2.798, 95% CI = −3.306 to −2.290, p < 1.0 × 10–8).
Meta-analysis of the relationship between telomere length and systemic lupus erythematosus. Meta-analysis of the association between telomere length and systemic lupus erythematosus SMD: standard mean difference; CI: confidence interval; WB: whole blood; PBMC: peripheral blood mononuclear cell; qPCR: quantitative polymerase chain reaction; TRF: terminal restriction fragment; FISH: fluorescence in situ hybridization; F: fixed effects model; R: random effects model. Magnitude of Cohen’s d effect size (SMD): 0.2–0.5: small effect; 0.5–0.8: medium effect; ≥0.8: large effect. Meta-analysis of the relationship between telomere length and systemic lupus erythematosus in specific ethnic groups.

Meta-analysis of telomere length in SLE patients compared with controls in subgroup
Meta-analysis was performed on SLE patients in each subgroup based on sample type, assay method and data type. Telomere length was significantly shorter in the SLE group than in the control group in whole blood and peripheral blood mononuclear cell groups (SMD = −0.361, 95% CI = −0.553 to −0.169, p = 2.3 × 10−4; SMD = −1.546, 95% CI = −2.583 to −0.510, p = 0.003; respectively), and a similar trend was observed in the leukocyte group (SMD = −0.699, 95% CI = −1.511 to −0.114, p = 0.092). Subgroup analysis by the assay method showed a significantly shortened telomere length in the SLE group by the quantitative polymerase chain reaction (qPCR) and terminal restriction fragment methods (SMD = −0.743, 95% CI = −1.478 to −0.009, p = 0.047; SMD = −1.163, 95% CI = −1.721 to −0.604, p = 4.6 × 10–5; respectively); a similar association was observed by the fluorescence in situ hybridization (FISH) method (SMD = −0.429, 95% CI = −0.859 to −0.000, p = 0.050). Stratification by the input data type revealed a significantly shortened telomere length in the original data and imputed data groups (Table 2).
Heterogeneity, sensitivity test and publication bias
Between-study heterogeneity was identified during the meta-analyses of telomere length in SLE patients (Table 2). Meta-regression analysis showed that ethnicity (p < 0.001), sample type (p < 0.001) and sample size (p = 0.012), but not publication year, assay method and data type (p > 0.05), had significant impacts on the heterogeneity in the meta-analysis of telomere length in SLE patients. Sensitivity analysis showed that no individual study significantly affected the pooled OR, indicating that the results of this meta-analysis are robust. It was difficult to correlate the funnel plot, which is usually used to detect publication bias, because the number of studies included in the analysis was relatively small. Egger’s regression test showed no evidence of publication bias (Egger’s regression test p-values > 0.1), but the funnel plot showed asymmetry; therefore, the ‘trim and fill’ method was used to adjust for publication bias (Figure 3). However, ORs that had been significant before this adjustment remained significant after the adjustment (SMD = −0.914, 95% CI = −1.363 to −0.465).
Funnel plot of studies that examined the association between the telomere length and systemic lupus erythematosus (Egger’s regression p-value = 0.327). The filled circle represents studies that showed publication bias. The diamonds at the bottom of the figure show summary effect estimates before (open) and after (filled) publication bias adjustment.
Discussion
Immune aging is associated with loss of immune functions, leading to autoimmunity, infection and malignancy.6–9 Critically shortened telomeres lead to cell senescence or apoptotic death and, thus, susceptibility to autoimmune diseases increases in the condition of accelerated immunosenescence. 5 Telomere shortening causes chromosomal instability that may lead to development of autoimmune diseases, including SLE. 7
Telomere length has been studied as a possible marker for SLE risk, showing mixed results. In this meta-analysis, we combined the evidence on the relationship between the telomere length and SLE. This meta-analysis of 12 comparative studies involving 472 SLE patients and 365 controls showed that telomere length was significantly shorter in the SLE group than in the control group. The significantly shortened telomere length in the SLE group was found in Caucasian, Asian and mixed populations. In addition, the telomere length is significantly shorter in patients with SLE, regardless of sample type or assay method evaluated. Accelerated immune aging, represented by telomere length reduction, has been involved in a variety of rheumatic diseases, including SLE. Immunosenescence by telomere shortening is related to increased reactivity to self-antigens, and loss of immune tolerance. Our meta-analysis showed that there was a significant association between shortened telomere length and SLE, independent of potential confounders such as age or ethnicity. It is well known that African-American patients with SLE have generally severe disease that is difficult to manage. A single African-American study showed a significantly shortened telomere length in SLE group (SMD = −2.798, 95% CI = −3.306 to −2.290, p < 1.0 × 10–8. Telomere length was lower in African-American patients with SLE than in Caucasian patients (SMD –2.798 vs. –0.445). Thus, the results of our study were consistent with the known finding.
It is not known whether this association is the cause or the consequence of shortened telomere length. There are some possible explanations. First, telomere shortening results from repeated mitotic divisions, exposure to oxidative stress, or DNA repair mechanisms. 4 SLE is associated with increased oxidative stress and can lead directly to accelerated telomere attrition. Immune activation observed in SLE may accelerate aging of the immune system. Telomere shortening is also influenced by environmental factors, lifestyle and socioeconomic factors. 29 Second, SLE patients may have congenitally shorter telomeres with decreased functionality, thus resulting in disease. 30 Telomere shortening is prevented via up-regulation of telomerase expression. 31 Genetic variation in telomerase, a regulator of telomere length, may be associated with shortened telomere length. 32 Third, initial autoantibodies developed in SLE may result in high leukocyte activation, leading to shortened telomeres; 33 although a mechanism by which either autoantibodies are produced or telomeres are congenitally shorter needs to be elucidated. A previous literature review investigated whether accelerated telomere shortening is a cause or a consequence of SLE. 30 The review indicated that accelerated telomere shortening may be a cause and not a consequence of disease development, because healthy relatives of patients have similarly short telomere length compared with patients and there seems to be a prelude of accelerated telomere shortening, which increases risk for disease development. 30
Multiple methods have been used for measurement of telomeres. 34 Terminal restriction fragment (TRF) is a convenient tool used for determining telomere length and, thus, is often called as the ‘gold standard’ method. Genomic DNA is exhaustively digested using frequent cutting restriction enzymes that lack recognition sites in the telomeric and subtelomeric regions. However, the restriction enzymes used can lead to the inclusion of subtelomeric DNA that is contiguous with the telomere, thereby possibly causing overestimation of the true telomere length. In addition, TRF requires a large amount of DNA and is labor-intensive. qPCR is a popular method for estimating telomere length because of the requirement for smaller amounts of DNA, its low cost and amenability for high-throughput testing. However, qPCR does not provide absolute length estimate. FISH uses fluorescent probes to quantify not only mean telomere lengths but also chromosome-specific telomere lengths, However, it requires mitotically active cells for metaphase chromosomes and cannot be used in cells that are not mitotically active.
The present study has some shortcomings that should be considered. First, most of the studies included in this meta-analysis had small sample sizes; thus, many of the individual studies that constitute this meta-analysis may have been underpowered. Second, the studies included in the meta-analysis were heterogeneous in demographic characteristics and clinical features. The heterogeneity and confounding factors such as disease activity and drugs used (e.g. immunosuppressive agents, hydroxychloroquine, corticosteroids) may have affected our results, which may be compounded by the limited information provided on clinical status and disease activity. This limited data did not allow further analysis, although we performed a sensitivity test and a meta-regression analysis.
Nevertheless, this meta-analysis has its strengths. The number of patients from individual studies ranged from five to 63, but our pooled analysis included 472 patients. Compared with individual studies, our study was able to provide more accurate data on the relation between the telomere length and SLE by increasing the statistical power and resolution through pooling of the results of independent analyses.
Our meta-analysis demonstrates that the telomere length is significantly shorter in patients with SLE, regardless of ethnicity, sample type and assay method evaluated, thereby suggesting that the shortened telomere length may be a risk factor for SLE or a consequence of SLE. Our meta-analysis indicated that telomere length likely plays an important role in the pathogenesis of SLE. Further studies are necessary to elucidate that a shortened telomere length directly contributes to the pathogenesis of SLE.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by a Korea University grant.
