Mendelian Randomization Study Implies Causal Linkage Between Telomere Length and Juvenile Idiopathic Arthritis in a European Population.

BACKGROUND: Telomere maintenance is increasingly being considered as fundamental to the progression of immune-mediated inflammatory diseases. However, the causality underlying the purported relationship has not been fully elucidated. In the present work, we applied Mendelian randomization (MR) analysis to obtain estimates of the causal effect of telomere length (TL) on the risk of juvenile idiopathic arthritis (JIA) and JIA-associated iridocyclitis.
METHODS: Two-sample MR analysis was conducted using summary-level data from the largest genome-wide association studies concerning TL (78,592 individuals), JIA (6056 cases and 25,086 controls), and JIA-associated iridocyclitis (1430 cases and 9,2767 controls). All the participants were of European ancestry. The inverse variance weighted (IVW) method was applied to estimate the causal effects. Sensitivity analyses incorporating multiple complementary MR approaches were implemented to test the robustness of the association and examine potential bias from pleiotropy.
RESULTS: In our MR analysis, genetically predicted shorter TL was associated with an increased risk of JIA (IVW: odds ratio=1.68, 95% CI: 1.13-2.48, P=0.009), but not with the risk of JIA-associated iridocyclitis (IVW: odds ratio=1.75, 95% CI: 0.81-3.79, P=0.155). The other MR methods produced consistent results. Besides, a leave-one-out sensitivity analysis yielded similar findings and validated the robustness of the causal relationship. MR-Egger regression revealed no notable horizontal pleiotropy (intercept=0.046, P=0.175).
CONCLUSION: This work provides evidence of a negative association between TL and JIA risk, but not for the association between TL and the risk of JIA-associated iridocyclitis, in a European population. Future studies with larger sample sizes are warranted to elucidate the underlying role of TL in these diseases.

Introduction

Human telomeres are repetitive stretches of the hexameric nucleotide sequence TTAGGG located at the natural ends of linear chromosomes.1 Together with the associated protein complex capping the telomere, they are mainly dedicated to preserving DNA integrity and chromosomal stability in the human genome. In most somatic tissues, telomere attrition accompanied by the DNA damage response will ultimately contribute to cellular senescence or apoptosis,2,3 which is recognized as a primary hallmark of aging. Research suggests that telomeres may be inversely associated with the risk of immune-mediated inflammatory diseases, such as inflammatory bowel diseases, psoriasis, and rheumatoid arthritis.4–6 Although the etiology of these conditions remains poorly defined, advances in molecular research on telomere length (TL) indicate that their pathogenesis could be attributable to the expression of common inflammatory cytokines associated with cellular aging.7,8

Juvenile idiopathic arthritis (JIA) is a leading cause of rheumatic disease affecting children under the age of 16 years,9 with iridocyclitis as the most common extra-articular comorbidity, which involves the ocular iris and ciliary body.10 The clinical traits are recognized as complex immune-mediated inflammatory disorders, and are subject to a complicated association with the immune system dysfunction, and various environmental and genetic factors.11,12 Several population-based observational studies have found that JIA patients display significantly shorter TL compared to healthy controls.13–16 The relationship between JIA-associated iridocyclitis and TL has been seldom discussed. However, the case–control study design used in all retrieved research is prone to reverse causality and residual confounding. The difficult assays used for TL measurement, small sample sizes in observational studies, and lower statistical power are further limitations. Therefore, whether TL has a causal effect on the risk of JIA or JIA-associated iridocyclitis has not yet been reliably inferred.

To validate these findings and avoid the potential impact from confounding factors, Mendelian randomization (MR) analysis is widely adopted to estimate the causal effects of genetically predicted exposures on one or more complex outcomes.17 Given that germline genetic information is fixed at conception in a random process according to Mendel’s laws of inheritance, the MR method utilizes genetic variants such as single-nucleotide polymorphisms (SNPs) as instrumental variables (IVs) for risk factors of interest,18 which can produce a more reliable association than before. To our knowledge, there has been no previous MR investigation on the causal link between TL and JIA. Therefore, in the current study, we performed a two-sample MR using summary-level statistics from large genome-wide association study (GWAS) databases to evaluate whether TL is causally associated with JIA risk. Furthermore, we also examined the potential causal sequence of TL on JIA-associated iridocyclitis.

Materials and Methods

Data Sources

We adopted a standard two-sample MR framework in which IV-exposure and IV-outcome interactions were derived from independent and non-overlapping individual datasets.

Summary statistics on leukocyte TL were obtained from the most recent, largest GWAS meta-analysis, which included EPIC-InterAct (n=19,799), EPIC-CVD (n=11,915), and ENGAGE consortium (n=46,898).19 All participants were of European ancestry. Mean leukocyte TL was measured using an established quantitative polymerase chain reaction technique and expressed as a ratio of the telomere repeat number (T) to a single-copy gene (S). The detailed procedure is described in a previous study.20 The TL measurements were standardized by quantifying the standard curve or using calibration samples. The final reported result for each sample was the average of three independently measured T/S values.21

JIA GWAS summary statistics were acquired from the recently updated meta-analysis of 6056 patients with all JIA subtypes and 25,086 controls of European descent.22 Detailed and clear descriptions, including the study procedure, were presented in the two involved studies.22,23 All eligible JIA patients must meet the criteria created and approved by the board-certified rheumatologists of the International League of Associations for Rheumatology (ILAR). The JIA GWAS data were systematically investigated regarding the specificity and sharing of JIA susceptibility loci across ILAR subtypes; that is, all subtypes of JIA were included in our cases without particular exploration of each subtype. Besides, to explore the causal inference between TL and JIA comorbidity (iridocyclitis), summary statistics from the MRC IEU Open GWAS Project (), with 1430 JIA-associated iridocyclitis cases and 92,767 controls of European ancestry from the FinnGen Study, were included in this work.

The study was approved by the ethical committee of each institutional review board, and written informed consent was provided by all participants.

Genetic Instruments for Telomere Length

Three essential model assumptions of MR analysis should be fulfilled to guarantee valid IVs (Figure 1), namely: 1) the selected IVs are robustly associated with the exposure (TL); 2) the eligible IVs for exposure are not associated with any confounders; and 3) the selected IVs only affect the outcome exclusively through the exposure (TL).

Figure 1

Mendelian randomization model of leukocyte TL and risk of JIA and JIA-associated iridocyclitis. The study design is under the assumption that genetic variants are directly associated with TL, but not indirectly with TL and not directly with confounders, that is, the instrumental variables influence JIA and JIA-associated iridocyclitis only directly through TL.

Based on the three MR core assumptions, the top 20 SNPs (shown in Table 1) regarded as IVs for leukocyte TL were chosen at the established threshold of GWAS significance (P<5×10−8) and assessed through additive models adjusted for age, gender, and cohort-specific covariates.19 These SNPs were mainly located close to the genes coding for proteins involved in telomere homeostasis. Proxy SNPs identified through the online platform LDlink () were applied when the SNPs were unavailable for the outcome in the GWAS. Besides, palindromic SNPs were not included in our research. The strength of the eligible IVs was calculated by R2 and the F-statistic.

Table 1

Characteristics of SNPs Selected as the IV and Association with Leukocyte TL

No.	SNP	Chr.	Position	Proximal	EA	OA	EAF	Beta	SE	P-Value	R2	F-Statistic	Remove
1	rs2695242	1	226,594,038	PARP1	C	A	0.830	0.042	0.006	9.6E-11	0.0021	41.9610
proxy of rs3219104 (r2=1)
2	rs10936600	3	169,514,585	TERC	T	A	0.243	−0.086	0.006	7.18E-51	0.0019	225.2677	Palindromic
3	rs4683917	3	101,186,025	SENP7	A	C	0.579	−0.037	0.007	2.45E-08	0.0030	31.1529
proxy of rs55749605 (r2=1)
4	rs13137667	4	71,774,347	MOB1B	C	T	0.959	0.077	0.014	2.43E-08	0.0011	31.1630
5	rs4691895	4	164,048,199	NAF1	C	G	0.783	0.058	0.006	1.58E-21	0.0024	90.9546	Palindromic
6	rs7705526	5	1,285,974	TERT	A	C	0.328	0.082	0.006	5.34E-45	0.0026	198.3350
7	rs2853677	5	1,287,194	TERT	A	G	0.592	−0.064	0.005	3.35E-31	0.0021	135.1105
8	rs2736176	6	31,587,561	PRRC2A	C	G	0.313	0.034	0.005	3.53E-10	0.0026	39.4232	Palindromic
9	rs34991172	6	25,480,328	CARMIL1	G	T	0.068	−0.061	0.010	6.19E-09	0.0011	33.8236
10	rs59294613	7	124,554,267	POT1	A	C	0.293	−0.041	0.005	1.17E-13	0.0021	55.1500
11	rs9419958	10	105,675,946	STN1(OBFC1)	C	T	0.862	−0.064	0.007	5.05E-19	0.0016	79.5213
12	rs228595	11	108,105,593	ATM	A	G	0.417	−0.028	0.005	1.43E-08	0.0024	32.1985
13	rs2302588	14	73,404,752	DCAF4	C	G	0.100	0.048	0.008	1.68E-08	0.0017	31.8763	Palindromic
14	rs3785074	16	69,406,986	TERF2	G	A	0.263	0.035	0.006	4.64E-10	0.0024	38.8811
15	rs62053580	16	74,680,074	RFWD3	G	A	0.169	−0.039	0.007	4.06E-08	0.0017	30.1692
16	rs7194734	16	82,199,980	MPHOSPH6	T	C	0.782	−0.037	0.006	6.94E-10	0.0020	38.0989
17	rs8105767	19	22,215,441	ZNF208	G	A	0.289	0.039	0.005	5.42E-13	0.0025	52.1238
18	rs75691080	20	62,269,750	RTEL1/STMN3	T	C	0.091	−0.067	0.009	5.99E-14	0.0013	56.4560
19	rs73624724	20	62,436,398	RTEL1/ZBTB46	C	T	0.129	0.051	0.007	6.33E-12	0.0019	47.3029
20	rs34978822	20	62,291,599	RTEL1	G	C	0.015	−0.140	0.023	7.26E-10	0.0004	38.0094	Palindromic

Abbreviations: SNP, single-nucleotide polymorphism; IV, instrumental variable; TL, telomere length; EA, effect allele; OA, other allele; EAF, effect allele frequency; SE, standard error.

Statistical Analysis

Mendelian Randomization Analysis

Our two-sample MR analysis was conducted using the “TwoSampleMR” R package.24 The inverse variance weighted (IVW) method18 was applied to obtain the combined estimate of the causal effect. The conventional IVW method relies on the assumption that all of the SNPs obtained in the analysis are valid instruments. Besides, the maximum likelihood-based method25 was utilized to provide appropriately estimated confidence intervals (CIs) when weak IVs exist. For each SNP, causal effect estimates were generated for relative outcomes as odds ratios (ORs) per one SD unit increase in the putative risk factor using the Wald ratio.

Sensitivity Analysis

After MR analysis, Cochran’s Q test was applied to assess heterogeneity among the effects of individual TL-associated SNPs on JIA or JIA-associated iridocyclitis. Besides, several important sensitivity analyses were performed to verify MR model assumptions and evaluate the MR results. The MR-Egger regression was used to identify potential directional pleiotropy26 and an intercept P-value >0.05 indicated a lack of horizontal pleiotropy. The MR-RAPS (robust adjusted profile score) approach was used to validate the robustness of systemic and idiosyncratic pleiotropy.27 For each outlying SNP, the MR-PRESSO (MR pleiotropy residual sum and outlier) test was applied to assess the causal estimates.28 Besides, we applied radial variants of the IVW model for a better visualization of the causal estimate. Leave-one-SNP-out analysis was also used to assess to identify potentially influential SNPs and verify the reliability of the results. With the PhenoScanner database (), those SNPs that correlated with confounding factors were manually screened and omitted.

All statistical analyses in this study were conducted using the “TwoSampleMR” package in R software (version 4.1.0). All presented P-values were two sided and statistical significance was set at the 5% level.

Results

In our MR analysis, 20 leukocyte TL-related SNPs were described as IV models. Two SNPs that failed the identification in the JIA GWAS were replaced by their corresponding proxy SNPs. We searched the selected SNPs in the PhenoScanner dataset and did not identify any potential SNPs in association with other confounding factors. After harmonization and removal of palindromic SNPs with intermediate allele frequencies, the remaining 13 sentinel SNPs were taken as IVs for TL. The F-statistics for these IVs were greater than 10 (range: 30.17–198.34), indicating that these SNPs were sufficiently strong for MR analysis. These variants explained around 2.98% of the variance in leukocyte TL.

Genetically Determined TL Was Associated with JIA Risk

No obvious evidence was observed for the heterogeneity of these TL-associated SNPs in Cochran’s Q test (P=0.577). The IVW method under a fixed-effects model was implemented for primary estimation. As presented in Table 2, primary MR analysis indicated that genetically predicted shorter TL was associated with increased risk of JIA (OR=1.68, 95% CI: 1.13–2.48, P=0.009) (Figure 2). Similar effect estimates were observed through the maximum likelihood-based method and MR-PRESSO method. Notably, we used the MR-RAPS estimator and found that the results remained largely consistent with our primary findings (OR=1.69, 95% CI: 1.13–2.51, P=0.01). This indicates the robustness of the causal association between TL and JIA. Furthermore, we tried to find the causal effect between genetically determined TL and JIA-associated iridocyclitis. No strong evidence was yielded based on IVW, MR-PRESSO, MR-RAPS, or maximum likelihood-based methods (IVW, OR=1.75, 95% CI: 0.81–3.79, P=0.155). The detailed data from these methods are presented in the Supplementary Material, .

Table 2

MR Results of the Association Between Leukocyte TL and the Risk of JIA

	N	OR (95% CI)	P-Value	P-Value for Heterogeneity	P-Value for Intercept	P-Value for Distortion	N	OR (95% CI)	P-Value	P-Value for Heterogeneity	P-Value for Intercept	P-Value for Distortion
	Full Set						Proxy SNPs Not Identified in Outcome
IVW	13	1.68 (1.13–2.48)	0.009	0.577	–	–	15	1.59 (1.09–2.30)	0.015	0.664	–
Maximum likelihood	13	1.69 (1.14–2.51)	0.008	0.583	–	–	15	1.60 (1.10–2.32)	0.014	0.669	–	–
MR-PRESSO	13	1.68 (1.16–2.41)	0.017	–	–	0.549	15	1.59 (1.14–2.21)	0.017	–	–	0.648
MR-RAPS	13	1.69 (1.13–2.51)	0.010	–	–	–	15	1.59 (1.09–2.33)	0.016	–	–	–
MR-Egger intercept	13	–	–	–	0.175	–	15	–	–	–	0.314	–

Abbreviations: MR, Mendelian randomization; TL, telomere length; JIA, juvenile idiopathic arthritis; N, number of SNPs used in MR analysis; OR, odds ratio; CI, confidence interval; IVW, inverse variance weighted; PRESSO, pleiotropy residual sum and outlier; RAPS, robust adjusted profile score.

Figure 2

Scatter plot of the effect size and 95% CIs of each SNP on leukocyte TL and JIA risk. The horizontal axis represents genetic associations of each variant on leukocyte TL. The vertical axis indicates the genetic association of each SNP with risk of JIA.

Abbreviations: SNP, single-nucleotide polymorphism; TL, telomere length; JIA, juvenile idiopathic arthritis.

In an alternative approach, we utilized the proxy package to explore the potential causal effects of SNPs that were not identified in the outcome factors. Our MR analysis also yielded very comparable point estimates through multiple MR methods (IVW, OR=1.59, 95% CI: 1.09–2.30, P=0.015) (Table 2). Thus, genetic predisposition towards shorter TL was found to be a causative factor for a higher JIA risk. However, no significant estimates were discovered between TL and JIA-associated iridocyclitis (IVW, OR=1.65, 95% CI: 0.84–3.27, P=0.147). The detailed data are presented in the Supplementary Material, .

Pleiotropy and Sensitivity Analysis

The MR-Egger regression did not suggest any evidence of directional pleiotropy (intercept=0.046, SE=0.032, P=0.175). The result of MR-PRESSO analysis did not highlight any outliers for any IVs (P-distortion=0.549) (Table 2). In addition, the IVW radial MR approach did not delineate any evidence of outlying genetic variants (Figure 3 and the Supplementary Material, ). In the sensitivity analyses, the leave-one-out method demonstrated that the causal links between genetically determined TL and JIA or JIA-associated iridocyclitis were not driven by a single SNP, suggesting the stability of our MR analysis (shown in the Supplementary Material, and ).

Figure 3

Radial plots to visualize individual outlier SNPs in the MR estimates for the effect of leukocyte TL and risk of JIA. The radial curve displays the ratio estimate for each SNP, as well as the overall IVW estimate (in blue). Black dots show valid SNPs.

As suggested by Brion et al (), we also calculated the sample power to validate our analysis. Under the current situation, our study had 99% power to detect a causal effect of a 68% increase in JIA risk at a significance level of 0.05.

Discussion

To the best of our knowledge, this is the first two-sample MR study to evaluate the relationships between leukocyte TL and JIA risk, as well as between leukocyte TL and the risk of JIA-associated iridocyclitis. Having made an effort to fulfill stringent MR assumptions and utilize multiple complementary MR methodologies, our results demonstrated that shorter genetically predicted TL could increase the risk of JIA in populations of European ancestry. However, there was no evidence supporting the causal linkage between TL and JIA-associated iridocyclitis.

TL, a recognized complex heritable trait, is influenced by oxidative damage and replicative stress caused by genetic and epigenetic factors, as well as extrinsic factors such as socioeconomic status, working conditions, and stress.29,30 Longer TL in offspring could be attributed to higher paternal age at conception.31 Besides, an increase in telomerase activity was inversely correlated with improvement in depressive symptoms and disease severity.32 It has been reported that shortened TL is associated with a high level of psychosocial stress, which could be attributed to increased oxidative stress and reduced telomerase activity.33 Moreover, in an observational study on rheumatoid arthritis, the levels of oxidative stress biomarkers rose as disease activity increased and thereafter led to significant telomere shortening.34 In our MR study, we disregarded genetic instruments associated with the above traits and solely explored the relationship between TL and diseases (JIA and JIA-associated iridocyclitis) from a genetic perspective.

Our findings were consistent with the evidence from limited observational studies; that is, genetically predicted TL was associated with the risk of JIA. JIA patients showed significantly shorter TL than healthy controls.13,16 A previous study demonstrated that shorter TL was measured in naive CD45RA+CD4+T cells from JIA patients with persistent immune activation.15 However, the exact biological mechanisms of telomere shortening in JIA remain to be elucidated. Several possible explanations may be advanced. First, some studies have indicated that the loss of the telomere would lead to functional consequences after infection due to autoimmune deficiency.13,35 Second, a change in TL may drive the premature immunosenescence of CD4+ CD28− T cells, which have been implicated in immune system diseases and chronic immune-mediated diseases owing to their pro-inflammatory and cytotoxic characteristics.36,37 In addition, some studies have shown that shortened or dysfunctional telomeres are associated with hyperactivity of the transcription factor NF-kB and overexpression of inflammatory cytokines such as TNF-alpha and IL-6 in circulating macrophages.38,39 Shortened TL has also been linked to the upregulation of NLRP3 inflammasome and the ATM-YAP1-pro-IL-18 pathway in inflammatory bowel diseases.4,40 These findings imply that TL may be involved in the pathogenesis of immune-mediated inflammatory diseases.

Contrary to the inevitable bias in observational research, our two-sample MR study, with a large sample size and sufficient statistical power, can address the concern of unmeasured or uncontrolled confounding factors by reanalyzing GWAS summary data. Furthermore, our results are in keeping with a causal role of shortened TL in increasing the risk of multiple immune-mediated disorders such as rheumatoid arthritis and hip osteoarthritis, as discovered in a growing number of MR studies.6,41 Although the exact role of TL in JIA has not yet been clearly defined, these findings suggest that telomere attrition seems to be a common trait of these disorders. Thus, the modulation of TL, except for a better understanding of various TL analyses, may represent a potential therapeutic approach for taking control of a disease risk.

In the MR study, we discovered that TL was related to JIA risk, rather than to the risk of JIA-associated iridocyclitis, suggesting the specificity of its comorbidity. Unfortunately, no studies have been published to evaluate the role of TL in JIA-associated iridocyclitis until now. One study revealed longer TL in ocular birdshot uveitis patients compared to unrelated healthy controls.42 These contradictory phenomena may be attributed to deleterious lifestyle, distinct cellular turnover due to chronic inflammation, and other determinants in the disease process.3 Furthermore, the relatively small populations and low incidence of JIA-associated iridocyclitis would unquestionably result in the inconsistency of the results. Further studies are warranted to explore the biological mechanisms of TL in JIA-associated iridocyclitis.

Major strengths of this MR study include the robust genetic instruments identified from the largest GWAS, explaining 2.98% of the variance in TL. In order to reduce bias in the causal estimate, we conducted the two-sample MR, where IVs associated with JIA or JIA-associated iridocyclitis and the estimation of the IVs were derived from two independent individuals. Moreover, we used stringent criteria for the inclusion of SNPs as IVs in our analyses and undertook comprehensive assessment for confounding. The intercept of the MR-Egger regression test was close to zero, suggesting no strong evidence for directional pleiotropy. Besides, no outliers were highlighted by the MR-PRESSO test, which is often used to detect and correct the outliers in IVW linear regression. Therefore, the consistent results observed across multiple complementary approaches lend further support to our findings.

Despite multiple advantages of the stringent MR assumptions, several limitations of the current study are worthy of mention. First, since the included participants in this study were of European ancestry, our findings may not be extrapolated to other populations, and it is urgent that our results are verified in various races. Second, the full individual-level statistics of the most recent GWAS were not available, which made it difficult to perform further functional validation. Third, the effects of age, sex, and other confounding environmental factors on TL cannot be ignored. Only summary-level data was applied in our MR analysis, which limited further stratified analyses on these specific factors. Fourth, the uncertain incompatibility in sample size between JIA and JIA-associated iridocyclitis may lead to some bias in our MR analysis. Fifth, though the genetic effects on the exposure of interest are independent and fixed throughout the life course, the age incompatibility between TL in the adult population and JIA in pediatric individuals should be noticed in our MR study. Whereas various JIA subtypes were taken together into the evaluation of their causal inference with TL, the exploration on the association of genetically predicted TL with a specific JIA subtype is still a challenge in current MR analysis. Therefore, subgroup analysis based on a set of JIA subtypes is recommended in the further GWAS study. Finally, JIA and its comorbidities are considered complex clinical traits involving various components, including genetic, immunological, environmental, and socioeconomic factors, suggesting that genetic variants could temptingly and partially explain the results.

Conclusions

To our knowledge, our results provide novel evidence to support a potential effect of leukocyte TL on JIA risk. The relationship between genetically predicted leukocyte TL and JIA-associated iridocyclitis could not be confirmed in the present work. Further large-scale studies among different populations are warranted to validate our findings and elucidate the causal role of TL in these disorders.