Front. Immunol.Frontiers in ImmunologyFront. Immunol.1664-3224Frontiers Media S.A.10.3389/fimmu.2023.1120742ImmunologyOriginal ResearchAnkylosing spondylitis and glaucoma in European population: A Mendelian randomization studyLiShengjie1ChenMinting1ZhangQing2FangMeijin1XiongWei1BaiLang1*1Department of Ophthalmology, Nanfang Hospital, Southern Medical University, Guangzhou, China2Department of Thoracic Surgery and Oncology, The First Affiliated Hospital of Guangzhou Medical University, China State Key Laboratory of Respiratory Disease and National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Disease, Guangzhou, China
Edited by: James Cheng-Chung Wei, Chung Shan Medical University Hospital, Taiwan
Reviewed by: Wuu-Tsun Perng, National Pingtung University of Science and Technology, Taiwan; Yi-Chiao Bai, Chung Shan Medical University, Taiwan
*Correspondence: Lang Bai, bailangsfy@126.com
This article was submitted to Autoimmune and Autoinflammatory Disorders : Autoimmune Disorders, a section of the journal Frontiers in Immunology
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Background
The relationship between ankylosing spondylitis (AS) and glaucoma in the European population remains unclear. In the present study, we applied a two-sample Mendelian randomization (MR) method to investigate their causal relationship.
Methods
MR analysis was conducted to validate the causal associations between AS with glaucoma using summary statistics from the genome-wide association studies of AS (9,069 cases and 13,578 control subjects) and glaucoma (8,591 cases and 210,201 control subjects). The inverse variance weighting method was performed to evaluate the causal relationship. The MR–Egger regression approach was applied to assess pleiotropy, while Cochran’s Q test was used to analyze heterogeneity. Subgroup analysis was performed according to primary open-angle glaucoma (POAG) and primary angle-closure glaucoma (PACG).
Results
The results of the MR study reveal a risk-increasing causal relationship between AS and glaucoma among European populations (OR = 1.35, 95%CI = 1.16–1.57, P = 8.81 × 10-5). Pleiotropy and heterogeneity were not found in our study. In the subgroup analysis, AS was also causal with POAG (OR = 1.48, 95%CI = 1.17–1.86, P = 8.80 × 10-4) and PACG (OR = 1.91, 95%CI = 1.03–3.51, P = 3.88 × 10-2).
Conclusion
The results of the MR analysis suggested a causal relationship between AS and glaucoma in the European population. Further studies are needed to identify the specific mechanism between these two diseases.
Glaucoma is a neurodegenerative disease characterized by the atrophy and depression of the optic disc and the defects of visual field and vision (1). Tham et al. reported a global prevalence of glaucoma of 3.54% in adults aged 40–80 years, of which the prevalence of primary open-angle glaucoma (POAG) was 3.05%, and it is estimated that, by 2040, the number of glaucoma patients worldwide will increase to 111.8 million (2). Glaucoma is the leading cause of irreversible blindness, with disease burden varying based on glaucoma subtypes and ethnicity. Without prompt treatment, it can result in irreversible visual impairment (3, 4). Pathological elevated intraocular pressure is the main risk factor for glaucoma (5). However, it cannot explain all glaucoma-related symptoms. Even after intraocular pressure (IOP) recovers to normal, some patients continue to acquire progressive visual field defects and optic atrophy. Therefore, early detection and diagnosis are undoubted of great significance to glaucoma.
Autoimmunity, inflammation, oxidative stress, mitochondrial dysfunction, and retinal ischemia have all been reported to be the risk factors for glaucoma (6–10). Several pieces of research have implicated antibodies and CD4+ T cells and the pathogenesis of glaucoma (11, 12). An increasing number of studies have unmasked the link between autoimmune diseases and glaucoma.
As an autoimmune disease, ankylosing spondylitis (AS) is a chronic arthritis characterized by the inflammation of the spine and sacroiliac joints (13, 14). AS affects the eyes through immune, cytokine, and infection mechanisms. There are 15.8% of AS patients who have the ocular disease (15). Uveitis is the most prevalent extra-articular manifestation, and there exists a correlation between certain ocular diseases and AS. Patients with optic neuritis are a high-risk population for autoimmune diseases, including AS (16, 17). However, due to limited data on AS cases, there are few studies on the association between AS and glaucoma (even primary glaucoma). Most studies on the relationship between AS and glaucoma focus only on the fact that AS causes uveitis, which further leads to secondary glaucoma. Clinicians often consider the possibility of AS in the face of uveitis patients in clinical practice. Lea et al. found that Europeans with the interleukin (IL)-38 polymorphism rs3811058 (IL-1F10.3) were more likely to develop AS (18). Recent studies have reported that IL-38 is involved in inflammatory autoimmune diseases, including AS and glaucoma (19). Furthermore, IL-36, IL-37, and IL-38 may be involved in the immune-mediated pathogenesis of chronic primary angle-closure glaucoma (PACG) and are associated with a mean deviation of visual fields (19, 20). Existing evidence suggests an association between them. Nonetheless, no studies have reported a direct causal relationship between AS and glaucoma.
Mendelian randomization (MR) analysis, which incorporates genetic variation as an instrumental variable (IV), has been frequently used to evaluate the potential causal association between risk factors and disease results (21, 22). Accordingly, MR analysis views genetic variation as a natural experiment entailing the random assignment of alleles associated with exposure during the time of gamete formation, in accordance with Mendel’s second law, which is typically independent of environmental risk factors and precedes risk factors and disease progression (23). More specifically, MR analysis avoids bias due to confounding or reverse causality as the genotype cannot be altered by the disease (24). Accordingly, MR analysis is less susceptible than traditional observational studies to confounding, reverse causality, and measurement error (25). By taking single-nucleotide polymorphisms (SNPs) as genetic variants, this approach bears a tremendous advantage in that the SNP-exposure effects and SNP-outcome effects are obtained from separate studies. With these summary data alone, even if they are not measured in the same set of samples, it is possible to estimate the causal influence of the exposure on the outcome (26). In the absence of observational data on both, we were still able to apply the two-sample MR analysis for the purpose of exploring the causal relationship between AS and glaucoma. Consequently, the aim of our study was to apply a two-sample MR analysis in order to identify whether AS is causally associated with the risk of glaucoma.
MethodsStudy design
The primary research of the genome-wide association studies (GWASs) from each nation had previously gained ethical consent, and the present analysis was based on publicly available summary-level data, which did not need additional approval. The publicly available summary data from the MR-Base platform were used in this study to assess the relationship between AS and glaucoma (26). To reduce the influence of population stratification, we restricted our study population to the individuals of European ancestry.
Selection of genetic instruments for ankylosing spondylitis
For the present study, we employed SNPs associated with AS at the genome-wide significance threshold from Cortes et al.’s GWASs (27). This GWAS comprised 9,069 cases and 13,578 controls from 22,647 subjects of European ancestry. The association between each AS-contained SNP [p< 5 × 10−8; linkage disequilibrium (LD) r2< 0.001] was evaluated, and no LD was found. If a requested SNP was absent, the proxy SNP (LD at r2>0.8) from the 1000 Genomes Project was employed as the surrogate effect allele. In this case, LD proxy lookups are automatically provided by MR-Base (26). In addition, we removed the SNPs with the presence of palindromes. Eventually, 24 SNPs as genetic variables were conducted into the final IV set (Supplementary Table S1, Supplementary Figures S1, S4, S7).
Data of genetic instruments for glaucoma
In order to avoid the inclusion of data from the same participant, the Data of Genetic Instruments for glaucoma was selected from the FinnGen consortium, with a total of 218,792 individuals of European ancestry (including 8,591 glaucoma cases and 210,201 controls) (28). Furthermore, POAG (3,412 cases and 210,201 controls) and PACG (588 cases and 210,201 controls) were included as histological subgroups for further comprehension of the causality between genetically predisposed AS and glaucoma subtypes. This study defines glaucoma by H40-H42 of the International Classification of Disease-10, and other studies have supported the use and accuracy of the International Classification of Disease coding for this diagnosis (28, 29) (Figure 1).
Overview of Mendelian randomization analysis and major assumptions. AS, ankylosing spondylitis; GWASs, genome-wide association studies; POAG, primary open-angle glaucoma; PACG, primary angle-closure glaucoma; *Confounder: smoking, myopia, diabetes.
Assessment of the Mendelian randomization assumptions
The following three assumptions must be fulfilled for the two-sample MR analysis to produce a correct interpretation: (I) SNPs are strongly associated with AS, usually, to satisfy this condition, with a p-value less than 5 × 10−8; (II) the IVs do not influence glaucoma through other pathways than the AS; and (III) the instrument variables do not associate with confounders. In addition, there is some potential pleiotropy in instrumental effects, which may indicate potential violations of the assumptions underlying two-sample MR. Therefore, other sensitivity analyses were carried out using multiple imputation methods to test the robustness of MR results. Accordingly, the first assumption can be quantified by the strength of an association between a genetic variant and the risk factor. On top of that, we calculated the F-statistic to measure the strength of IVs (30). This is related to the proportion of variance in the phenotype explained by the genetic variants (R2). The phenotype variance explained (R2) was obtained from GWAS databases or original literature. If the R2 were not provided in the GWAS databases or literature, we would assess the strength of instruments based on βexposure2 or SEexposure2. Subsequently, to verify the second assumption, we performed directional pleiotropy based on the intercept obtained from the MR–Egger regression analysis (31, 32). Concerning the third hypothesis, once these selected R2 were associated with confounding factors for AS and glaucoma, our MR analysis would prove unable to estimate the causal effect of AS on glaucoma accurately. Consequently, we manually searched in all the published GWASs whether AS-related SNPs have secondary phenotypes other than AS. The cut-off date was on 30th August 2022. In addition, smoking represents a known risk factor for greater disease activity and worse functioning in AS patients (33). Simultaneously, a number of studies have reported that smoking intensity is associated with the faster thinning of the retinal nerve fiber layer (34). In consequence, smoking may be a potential confounding factor in the relationship between AS and glaucoma (35). In addition, epidemiologic studies show that both myopia and diabetes may be risk factors for glaucoma (36–38). It cannot be ruled out that confounding factors such as myopia and diabetes may further affect glaucoma by influencing AS. Hence, we selected these three assumptions as possible confounding factors to study AS (Figure 1).
Statistical analysis
Inverse variance weighting (IVW) is generally regarded as the main method in MR analysis to test the causality between exposure and the outcome: in the absence of pleiotropy (and assuming that the instruments are valid), IVW returns the unbiased estimates of a causal effect so long as horizontal pleiotropy is balanced (26, 31). Accordingly, if there was just one SNP, we estimated the causal effect of AS on glaucoma using the Wald ratio method (22). In addition, we estimated the effects using the weighted median (WM) and MR–Egger methods (26). In this regard, the MR–Egger method has the advantage of allowing horizontal pleiotropic effects. Accordingly, if the results are consistent with the IVW method, then, confidence is gained in the conclusions. In this study, power calculations were performed based on a method suggested by Burgess et al. to confirm whether the IVs provided relatively accurate estimates of causal effects (39).
Finally, we performed Cochran’s Q test for the assessment of heterogeneity in order to assess whether the estimate of the causal effect of AS on glaucoma was consistent across each SNP (40). On top of that, a leave-one-out analysis was conducted to estimate whether the result was biased or driven by a single SNP, which sequentially excluded one SNP at a time, in order to assess the sensitivity of the results for individual variants (41).
MR analysis was performed in R (version 3.6.3) with the package “TwoSampleMR” (version 0.5.6).
ResultsMendelian randomization analysis of ankylosing spondylitis with risk of glaucoma and subtypes
Our results showed that AS was causally correlated with a markedly increased glaucoma risk among European populations. IVW was used to indicate that AS (OR = 1.35, 95%CI = 1.16–1.57, P = 8.81 × 10-5) was associated with the risk of glaucoma. We further validated the results of IVW by WM (OR = 1.42, 95%CI = 1.15–1.74, P = 1.08 × 10-3) and MR–Egger (OR = 1.33, 95%CI = 1.03–1.73, P = 3.97 × 10-2), and the results are similar. In addition, in the subgroup analysis, AS was found to be causally associated to POAG (OR = 1.48, 95%CI = 1.17–1.86, P = 8.80 × 10-4) and PACG (OR = 1.91, 95%CI = 1.03–3.51, P = 3.88 × 10-2), respectively, and these associations were consistent in the methods of WM and MR–Egger, with similar causal estimates in direction and amplitude (Table 1, Supplementary Figures S2, S5, S8).
Mendelian randomization (MR) results between ankylosing spondylitis (AS) and glaucoma.
Verification of the Mendelian randomization assumptionsF-statistic
We assessed the strength of instruments by an F-statistic. The F-statistic of AS was much larger than 10, indicating that the instrument accurately predicted glaucoma.
MR–Egger pleiotropy test
Furthermore, the impact of pleiotropy of the exposures given the intercept value could be negligible as no evidence for the presence of directional pleiotropy in the MR–Egger regression analysis was found (interrupt β<0.001, P = 0.92) (Table 2).
Causal effect from ankylosing spondylitis on potential glaucoma risk factors
For the third hypothesis, we found no evidence that the included AS-associated SNPs have secondary phenotypes other than AS during our manual search. Moreover, the overall correlations between AS and confounders (smoking, myopia, and diabetes) were not statistically significant (Table 3).
Under the supposition of the 24 selected SNPs explaining a total of 24.4% variance of AS, the MR analysis had 100% power at a significance level of 0.05 to detect the previously estimated causal effect size of AS (OR = 1.35), suggesting a more reliable estimation of the causal effect.
Heterogeneity and leave-one-out sensitivity analysis
According to the results of Cochran’s Q test, there was no heterogeneity among these SNPs (Supplementary Table S2). Leave-one-out sensitivity analysis could be found in the appendix. It seemed that no single SNP was strongly driving the overall effect of each exposure on glaucoma (Supplementary Figures S3, S6, S9).
Discussion
In our study, MR analysis revealed a significant causal relationship between AS and glaucoma, which was the first large-scale MR study to investigate the correlation between AS and glaucoma. This study provided a better understanding of potential risk factors for glaucoma in AS patients.
Ocular involvement is one of the most common extra-articular manifestations in AS patients. The pathogenesis of ocular disease in AS patients has yet to be elucidated. The fact that 15.8% of AS patients will develop anterior uveitis makes ophthalmologists pay more attention to AS in clinical practice, which is related to human leukocyte antigen (HLA)-B27 and the disease duration of AS (42). AS has an increased risk of developing chronic diseases such as diabetes mellitus and cardiovascular and cerebrovascular diseases (43). However, AS patients with acute anterior uveitis (AAU) did not increase the risk of subsequent major adverse cardiovascular events compared to non-AAU controls (44). The mechanism of how AS modifies glaucoma risk remains unclear. According to certain studies on the association between AS and glaucoma, corneal hysteresis (CH) reduced as the duration of AS increased. CH is considered a risk factor for all types of glaucoma and progressive glaucoma optic neuropathy because it represents the ability of the corneal stromal tissue to absorb and disseminate energy (45, 46). Furthermore, central corneal thickness (CCT) values reduced with increasing Bath Ankylosing Spondylitis Metrology Index score, which may result in lower IOP power during measurement, while thinner CCT with ocular hypertension is more prone to glaucoma (5, 47).
It is worth noting that the MR analysis conducted in our study basically satisfied three key assumptions (validated by sensitivity analysis), which can remove the influence of most confounding factors and avoid reverse causality. MR analysis has given us a better understanding of potential risk factors for glaucoma.
In the current clinical work, the conventional treatment of glaucoma has been directed toward controlling IOP, primarily because that is the proven modifiable risk factor for glaucoma. Topical drug therapy, which has been commonly used in the past, is being challenged by selective laser trabeculectomy, microinvasive glaucoma surgery, and sustained delivery methods. Surgery and drugs to reduce intraocular pressure and protect optic ganglion cells are effective ways to save the vision of glaucoma patients (48). However, early detection and intervention of glaucoma are still essential. The data analysis of medical therapy for the normotensive fellow eye of dogs previously diagnosed with primary glaucoma suggests that utilizing a carbonic anhydrase inhibitor and a prostaglandin analog would be reasonable (49). Individualized, refined preventive measures are beneficial for patients with high-risk factors for glaucoma. In other words, it is critical for patients with AS to assess their physical status and manage modifiable risk factors in glaucoma from a public health perspective. With the increasing availability of a wealth of genetic data, the expansion of GWASs could lead to the realization of early detection of glaucoma and aid in the diagnosis of AS, potentially enabling genetic-based treatments (4, 50).
The association between AS and glaucoma is influenced by a number of intermediary phenotypes, which must be taken into account. It is widely known that anterior uveitis is an important clinical feature of AS, and its complications may lead to secondary glaucoma. The MR analysis revealed a clear causal association between AS and uveitis, as well as an inseparable correlation between uveitis and secondary glaucoma, which is an important intermediary factor when analyzing the relationship between uveitis and secondary glaucoma. However, the AS–uveitis–secondary glaucoma axis may be a significant confounding effect in MR analysis on the relationship between AS and primary glaucoma. Based on the nature of the data we used, AS could not be stratified by uveitis. In this study, we were also unable to separate AS from confounding factors. Therefore, we performed an MR analysis between AS and primary glaucoma to strengthen the validation of the causal relationship between AS and glaucoma. The results suggested that uveitis could not fully explain the relationships between AS and glaucoma, while internal mechanisms have not been clearly unmasked. In the analysis of the relationship between AS and glaucoma, uveitis complicated by AS has provided a mechanistic explanation that does impact the numerical results, even though this is within the expected range.
Furthermore, we should keep in mind that there are still potential confounding factors in the study of AS and glaucoma, such as age, gender, medication use, and the duration of the type of glaucoma exposure. However, further analysis was difficult to conduct as the corresponding GWAS data were not available. In theory, we cannot rule out all possible confounds that may occur in MR analysis. The MR analysis (closest to a randomized controlled experiment) remains the most effective method for determining causality in the presence of confounders in the face of fewer research data and a lower likelihood of conducting a randomized controlled trial as it can minimize the effects of confounders and provides sufficient statistical power for causal estimation.
There are also limitations to MR analysis. Firstly, our MR analysis included only 24 SNPs in published GWAS. Secondly, the samples involved in our study were limited to European populations, which may not apply to Asian or African populations who have a higher incidence of glaucoma. Therefore, larger sample studies with more AS-related SNPs and more ethnic groups may help provide more valid conclusions about the causal relationship between AS and glaucoma. Furthermore, genetic variation in the analysis can explain only a small amount of individual diversity. In other words, in social epidemiological investigations, the mean of the relevant group must be evaluated rather than the mean of the individual. Finally, we were unable to stratify the population according to medication usage, disease duration, the history of uveitis, and other factors to avoid potential biases caused by the aforementioned confounding factors, and we cannot completely rule out the possible influence of pleiotropic effects on the results.
In conclusion, we conducted the first MR study and subgroup analysis to investigate the causal relationship between AS and glaucoma in the European population. Our study revealed a significant causal relationship between AS and glaucoma, especially POAG or PACG. We used MR analysis to mimic randomized controlled trials, which can avoid reversing causality and potential confounders common in conventional observational studies. However, the specific mechanism by which AS causes primary glaucoma requires further investigation.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
SL: Conceptualization, Data Curation, Methodology, Formal analysis, Writing- Original draft, Writing-Reviewing and Editing. MC: Data Curation, Methodology, Formal analysis, Writing- Original draft. QZ: Formal analysis, Writing-Reviewing and Editing. MF: Methodology, Formal analysis. WX: Methodology, Formal analysis. LB: Conceptualization, Supervision, Validation, Writing-Reviewing and Editing. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1120742/full#supplementary-material
ReferencesCohenLPPasqualeLR. Clinical characteristics and current treatment of glaucoma. Cold Spring Harbor Perspect Med (2014) 4(6). doi: 10.1101/cshperspect.a017236ThamYCLiXWongTYQuigleyHAAungTChengCY. Global prevalence of glaucoma and projections of glaucoma burden through 2040: A systematic review and meta-analysis. Ophthalmology (2014) 121(11):2081–90. doi: 10.1016/j.ophtha.2014.05.013ResnikoffSPascoliniDEtya'aleDKocurIPararajasegaramRPokharelGP. Global data on visual impairment in the year 2002. Bull World Health Organ (2004) 82(11):844–51.ZukermanRHarrisAVercellinAVSieskyBPasqualeLRCiullaTA. Molecular genetics of glaucoma: Subtype and ethnicity considerations. Genes (2020) 12(1). doi: 10.3390/genes12010055AugustoABLucaBAlessandroBJoaoBBChiaraBAndreiB. European Glaucoma society terminology and guidelines for glaucoma, 5th edition. Br J Ophthalmol (2021) 105(Suppl 1):1–169. doi: 10.1136/bjophthalmol-2021-egsguidelinesTezelGYangXLuoCKainADPowellDWKuehnMH. Oxidative stress and the regulation of complement activation in human glaucoma. Invest Ophthalmol Visual Sci (2010) 51(10):5071–82. doi: 10.1167/iovs.10-5289TezelG. The role of glia, mitochondria, and the immune system in glaucoma. Invest Ophthalmol Visual Sci (2009) 50(3):1001–12. doi: 10.1167/iovs.08-2717WeinrebRNAungTMedeirosFA. The pathophysiology and treatment of glaucoma: A review. Jama (2014) 311(18):1901–11. doi: 10.1001/jama.2014.3192Von Thun Und Hohenstein-BlaulNBellKPfeifferNGrusFH. Autoimmune aspects in glaucoma. Eur J Pharmacol (2016) 787:105–18. doi: 10.1016/j.ejphar.2016.04.031OsborneNNCassonRJWoodJPChidlowGGrahamMMelenaJ. Retinal ischemia: mechanisms of damage and potential therapeutic strategies. Prog Retinal Eye Res (2004) 23(1):91–147. doi: 10.1016/j.preteyeres.2003.12.001BellKFunkeSPfeifferNGrusFH. Serum and antibodies of glaucoma patients lead to changes in the proteome, especially cell regulatory proteins, in retinal cells. PloS One (2012) 7(10):e46910. doi: 10.1371/journal.pone.0046910ChenHChoKSVuTHKShenCHKaurMChenG. Commensal microflora-induced T cell responses mediate progressive neurodegeneration in glaucoma. Nat Commun (2018) 9(1):3209. doi: 10.1038/s41467-018-05681-9BraunJSieperJ. Ankylosing spondylitis. Lancet (Lond Engl) (2007) 369(9570):1379–90. doi: 10.1016/S0140-6736(07)60635-7ChenCWWeiJCGuJYuD. Editorial: Advances in pathogenesis, etiology, and therapies for ankylosing spondylitis. Front Immunol (2021) 12:822582. doi: 10.3389/fimmu.2021.822582ChenCHLinKCChenHALiaoHTLiangTHWangHP. Association of acute anterior uveitis with disease activity, functional ability and physical mobility in patients with ankylosing spondylitis: a cross-sectional study of Chinese patients in Taiwan. Clin Rheumatol (2007) 26(6):953–7. doi: 10.1007/s10067-006-0403-2MaKSLeeYHLinCJShihPCWeiJC. Management of extra-articular manifestations in spondyloarthritis. Int J Rheum Dis (2023) 26(2):183–6. doi: 10.1111/1756-185X.14485MaKSLeeCMChenPHYangYDongYWWangYH. Risk of autoimmune diseases following optic neuritis: A nationwide population-based cohort study. Front Med (2022) 9:903608. doi: 10.3389/fmed.2022.903608LeaWILeeYH. The associations between interleukin-1 polymorphisms and susceptibility to ankylosing spondylitis: A meta-analysis. Joint Bone Spine (2012) 79(4):370–4. doi: 10.1016/j.jbspin.2011.06.010XieLHuangZLiHLiuXZhengSSuW. IL-38: A new player in inflammatory autoimmune disorders. Biomolecules (2019) 9(8). doi: 10.3390/biom9080345ZhangJLSongXYChenYYNguyenTHAZhangJYBaoSS. Novel inflammatory cytokines (IL-36, 37, 38) in the aqueous humor from patients with chronic primary angle closure glaucoma. Int Immunopharmacol (2019) 71:164–8. doi: 10.1016/j.intimp.2019.03.016LawlorDAHarbordRMSterneJATimpsonNDavey SmithG. Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology. Stat Med (2008) 27(8):1133–63. doi: 10.1002/sim.3034BurgessSSmallDSThompsonSG. A review of instrumental variable estimators for mendelian randomization. Stat Methods Med Res (2017) 26(5):2333–55. doi: 10.1177/0962280215597579SmithGDEbrahimS. 'Mendelian randomization': Can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol (2003) 32(1):1–22. doi: 10.1093/ije/dyg070HolmesMVAla-KorpelaMSmithGD. Mendelian randomization in cardiometabolic disease: Challenges in evaluating causality. Nat Rev Cardiol (2017) 14(10):577–90. doi: 10.1038/nrcardio.2017.78Davey SmithGHemaniG. Mendelian randomization: Genetic anchors for causal inference in epidemiological studies. Hum Mol Genet (2014) 23(R1):R89–98. doi: 10.1093/hmg/ddu328HemaniGZhengJElsworthBWadeKHHaberlandVBairdD. The MR-base platform supports systematic causal inference across the human phenome. eLife (2018) 7. doi: 10.7554/eLife.34408CortesAHadlerJPointonJPRobinsonPCKaraderiTLeoP. Identification of multiple risk variants for ankylosing spondylitis through high-density genotyping of immune-related loci. Nat Genet (2013) 45(7):730–8. doi: 10.1038/ng.2667LiuKWuPChenBCaiYYuanRZouJ. Implicating causal brain magnetic resonance imaging in glaucoma using mendelian randomization. Front Med (2022) 9:956339. doi: 10.3389/fmed.2022.956339LiuKFanHHuHChengYLiuJYouZ. Genetic variation reveals the influence of steroid hormones on the risk of retinal neurodegenerative diseases. Front Endocrinol (2022) 13:1088557. doi: 10.3389/fendo.2022.1088557BurgessSThompsonSG. Avoiding bias from weak instruments in mendelian randomization studies. Int J Epidemiol (2011) 40(3):755–64. doi: 10.1093/ije/dyr036BowdenJDavey SmithGBurgessS. Mendelian randomization with invalid instruments: Effect estimation and bias detection through egger regression. Int J Epidemiol (2015) 44(2):512–25. doi: 10.1093/ije/dyv080VerbanckMChenCYNealeBDoR. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Genet (2018) 50(5):693–8. doi: 10.1038/s41588-018-0099-7ZhangHWanWLiuJDaiSZouYQianQ. Smoking quantity determines disease activity and function in Chinese patients with ankylosing spondylitis. Clin Rheumatol (2018) 37(6):1605–16. doi: 10.1007/s10067-018-4016-3NishidaTMahmoudinezhadGWeinrebRNBaxterSLEslaniMLiebmannJM. Smoking and progressive retinal nerve fibre layer thinning in glaucoma. Br J Ophthalmol (2022). doi: 10.1136/bjo-2022-321237Pérez-de-ArcelusMToledoEMartínez-GonzálezMMartín-CalvoNFernández-MonteroAMoreno-MontañésJ. Smoking and incidence of glaucoma: The SUN cohort. Medicine (2017) 96(1):e5761. doi: 10.1097/MD.0000000000005761HaAKimCYShimSRChangIBKimYK. Degree of myopia and glaucoma risk: A dose-response meta-analysis. Am J Ophthalmol (2022) 236:107–19. doi: 10.1016/j.ajo.2021.10.007HuZZhouFKamingaACXuH. Type 2 diabetes, fasting glucose, hemoglobin A1c levels and risk of primary open-angle glaucoma: A mendelian randomization study. Invest Ophthalmol Visual Sci (2022) 63(5):37. doi: 10.1167/iovs.63.5.37ShenLWalterSMellesRBGlymourMMJorgensonE. Diabetes pathology and risk of primary open-angle glaucoma: Evaluating causal mechanisms by using genetic information. Am J Epidemiol (2016) 183(2):147–55. doi: 10.1093/aje/kwv204BurgessS. Sample size and power calculations in mendelian randomization with a single instrumental variable and a binary outcome. Int J Epidemiol (2014) 43(3):922–9. doi: 10.1093/ije/dyu005BowdenJDel GrecoMFMinelliCDavey SmithGSheehanNThompsonJ. A framework for the investigation of pleiotropy in two-sample summary data mendelian randomization. Stat Med (2017) 36(11):1783–802. doi: 10.1002/sim.7221EmdinCAKheraAVKathiresanS. Mendelian randomization. Jama (2017) 318(19):1925–6. doi: 10.1001/jama.2017.17219WendlingD. Uveitis in seronegative arthritis. Curr Rheumatol Rep (2012) 14(5):402–8. doi: 10.1007/s11926-012-0267-1NurmohamedMTvan der Horst-BruinsmaIMaksymowychWP. Cardiovascular and cerebrovascular diseases in ankylosing spondylitis: Current insights. Curr Rheumatol Rep (2012) 14(5):415–21. doi: 10.1007/s11926-012-0270-6BaiYCLiuCHLeongPYLaiKLChenHHWeiJC. The risk of major adverse cardiovascular events in ankylosing spondylitis patients with a history of acute anterior uveitis: A nationwide, population based cohort study. Front Med (2022) 9:884800. doi: 10.3389/fmed.2022.884800DeolMTaylorDARadcliffeNM. Corneal hysteresis and its relevance to glaucoma. Curr Opin Ophthalmol (2015) 26(2):96–102. doi: 10.1097/ICU.0000000000000130CheeRISilvaFQEhrlichJRRadcliffeNM. Agreement of flicker chronoscopy for structural glaucomatous progression detection and factors associated with progression. Am J Ophthalmol (2013) 155(6):983–90.e1. doi: 10.1016/j.ajo.2013.01.005CabukKSÜstünEIAtalayKKirgizAAydinR. Corneal biomechanical features in patients with ankylosing spondylitis. Arquivos brasileiros oftalmol (2017) 80(5):309–12. doi: 10.5935/0004-2749.20170075CassonRJ. Medical therapy for glaucoma: A review. Clin Exp Ophthalmol (2022) 50(2):198–212. doi: 10.1111/ceo.13989PlummerCEBrasDGrozdanicSKomáromyAMMcLellanGMillerP. Prophylactic anti-glaucoma therapy in dogs with primary glaucoma: A practitioner survey of current medical protocols. Vet Ophthalmol (2021) 24(Suppl 1):96–108. doi: 10.1111/vop.12820LiZWuXLeoPJDe GuzmanEAkkocNBrebanM. Polygenic risk scores have high diagnostic capacity in ankylosing spondylitis. Ann rheum Dis (2021) 80(11):e187. doi: 10.1136/annrheumdis-2020-219446