Intraocular pressure, primary open-angle glaucoma and the risk of retinal vein occlusion: a Mendelian randomization mediation analysis

preprint OA: closed CC-BY-NC-ND-4.0
📄 Open PDF Full text JSON View at publisher

Abstract

Background The etiological connection between intraocular pressure (IOP) and the risk of retinal vein occlusion (RVO) remains elusive, particularly regarding whether this risk emanates from the direct influence of elevated intraocular pressure (IOP), irrespective of the presence of primary open-angle glaucoma (POAG), or if it arises as a consequence of the sequelae of POAG. Therefore, we conducted a Mendelian Randomization (MR) mediation analysis to elucidate the mediating role of POAG in the association between IOP and RVO. Methods We identified 47 single-nucleotide polymorphisms (SNPs) associated with IOP (P-value < 5×10 -8 ) leveraging data from a genome-wide association study (GWAS) (N = 97,653) obtained from the UK Biobank and 50 SNPs associated with POAG (P-value < 5×10 -8 ) from a GWAS meta-analysis (16,677 cases and 199,580 controls). We related these SNPs with RVO using a GWAS of 775 RVO cases and 376,502 controls from FinnGen. By utilizing univariable and multivariable MR analyses we calculated the total effect of IOP on RVO and estimated the degree to which POAG mediates this association. Results MR analyses showed that higher IOP is associated with higher RVO risk (odds ratio of RVO per 1 mmHg increase in IOP: 1.53; 95% confidence interval: 1.04 to 2.26; p-value = 0.03). Moreover, our MR mediation analysis suggested that 91.6% of the total effect of IOP on RVO risk was mediated through POAG. The primary results were consistent with estimates of pleiotropy-robust MR methods. Conclusion Our findings suggest that higher IOP increases the risk of RVO and that the majority of this effect is mediated through POAG.
Full text 38,386 characters · extracted from preprint-html · click to expand
Intraocular pressure, primary open-angle glaucoma and the risk of retinal vein occlusion: a Mendelian randomization mediation analysis | medRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-P4HH5NV'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search Intraocular pressure, primary open-angle glaucoma and the risk of retinal vein occlusion: a Mendelian randomization mediation analysis View ORCID Profile Andreas Katsimpris , View ORCID Profile Sebastian-Edgar Baumeister , Nafsika Voulgari , View ORCID Profile Hansjörg Baurecht , View ORCID Profile Stylianos Kandarakis , View ORCID Profile Michael Nolde doi: https://doi.org/10.1101/2024.01.31.24302094 Andreas Katsimpris 1 Princess Alexandra Eye Pavilion, University of Edinburgh , Edinburgh, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Andreas Katsimpris For correspondence: a.katsimpris{at}gna-gennimatas.gr nolde{at}uni-muenster.de Sebastian-Edgar Baumeister 2 Institute of Health Services Research in Dentistry, University of Münster , Münster, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Sebastian-Edgar Baumeister Nafsika Voulgari 1 Princess Alexandra Eye Pavilion, University of Edinburgh , Edinburgh, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site Hansjörg Baurecht 3 Department of Epidemiology and Preventive Medicine, University of Regensburg , Regensburg, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Hansjörg Baurecht Stylianos Kandarakis 4 First Department of Ophthalmology, National and Kapodistrian University of Athens , Athens, Greece Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Stylianos Kandarakis Michael Nolde 2 Institute of Health Services Research in Dentistry, University of Münster , Münster, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Michael Nolde For correspondence: a.katsimpris{at}gna-gennimatas.gr nolde{at}uni-muenster.de Abstract Full Text Info/History Metrics Supplementary material Data/Code Preview PDF Abstract Background The etiological connection between intraocular pressure (IOP) and the risk of retinal vein occlusion (RVO) remains elusive, particularly regarding whether this risk emanates from the direct influence of elevated intraocular pressure (IOP), irrespective of the presence of primary open-angle glaucoma (POAG), or if it arises as a consequence of the sequelae of POAG. Therefore, we conducted a Mendelian Randomization (MR) mediation analysis to elucidate the mediating role of POAG in the association between IOP and RVO. Methods We identified 47 single-nucleotide polymorphisms (SNPs) associated with IOP (P-value < 5×10 -8 ) leveraging data from a genome-wide association study (GWAS) (N = 97,653) obtained from the UK Biobank and 50 SNPs associated with POAG (P-value < 5×10 -8 ) from a GWAS meta-analysis (16,677 cases and 199,580 controls). We related these SNPs with RVO using a GWAS of 775 RVO cases and 376,502 controls from FinnGen. By utilizing univariable and multivariable MR analyses we calculated the total effect of IOP on RVO and estimated the degree to which POAG mediates this association. Results MR analyses showed that higher IOP is associated with higher RVO risk (odds ratio of RVO per 1 mmHg increase in IOP: 1.53; 95% confidence interval: 1.04 to 2.26; p-value = 0.03). Moreover, our MR mediation analysis suggested that 91.6% of the total effect of IOP on RVO risk was mediated through POAG. The primary results were consistent with estimates of pleiotropy-robust MR methods. Conclusion Our findings suggest that higher IOP increases the risk of RVO and that the majority of this effect is mediated through POAG. Introduction Retinal vein occlusion (RVO) is the second most common retinal vascular disease following diabetic retinopathy and an important cause of vision loss in individuals older than 80 years with a prevalence of 4.6% [ 1 ]. Risk factors for RVO are mainly related to atherosclerosis, but also conditions altering the rheologic properties in the retinal veins, such as hypercoagulability and vasculitis [ 2 ]. Glaucoma, an ocular disorder characterized by optic nerve damage and visual field loss, has been implicated in the pathogenesis of RVO [ 3 ]. Specifically, the association between RVO and primary open angle glaucoma (POAG) has been investigated in numerous studies, showing a significantly higher prevalence of RVO in patients with POAG compared to the general population [ 3 , 4 ]. Different pathogenic mechanisms have been proposed, implicating intraocular pressure which is a well-established and modifiable risk factor for glaucoma. On one hand, increased IOP has been hypothesized to collapse the retina vessel walls thus resulting in occlusion [ 5 , 6 ]. On the other hand, structural changes of the vessels at the glaucomatous optic discs might also make them more susceptible to thrombosis [ 7 – 10 ]. Whether RVO is caused by the effects of increased IOP independently from any glaucomatous optic disc abnormalities has not been determined. Mendelian randomization (MR) serves as a robust methodology utilizing genetic variants derived from genome-wide association studies (GWAS) as instrumental variables to elucidate the causal relationship between modifiable exposures and disease outcomes [ 11 ]. A recent MR study has identified a positive association between POAG and RVO [ 12 ]. However, it was not assessed whether this observed effect is intrinsic to POAG itself or if it is a consequence of the elevated intraocular pressure that is mediated through POAG. In addressing this gap, we conducted a MR mediation analysis to discern the mediating role of POAG in the relationship between IOP and RVO and sought to enhance our understanding of the nuanced interplay between these factors. Materials and Methods Study design MR employs genetic variants as instrumental variables to investigate the impact of risk factors on disease susceptibility [ 11 ]. This method mitigates susceptibility to confounding and reverse causation, as these genetic variants are randomly assigned at conception, creating a quasi-randomized exposure allocation analogous to a randomized trial [ 11 ]. In our investigation, a two-sample MR was executed using summary statistics of single-nucleotide polymorphisms (SNPs) obtained from GWAS focusing on IOP [ 13 ] and RVO [ 14 ]. The aim was to evaluate the influence of IOP on the risk of RVO. Moreover, we tried to assess the proportion of the effect of IOP on RVO that is potentially mediated through POAG, by employing a two-step MR for mediation analysis [ 15 ]. The study protocol was not pre-registered and our research adhered to the STROBE-MR guidelines [ 16 ] and “Guidelines for performing Mendelian randomization investigations” [ 17 ]. Data sources Summary data for corneal-compensated IOP were obtained from the UK Biobank (UKBB) cohort GWAS [ 13 ] (Supplementary Table 1). In particular, a subset comprising 97,653 individuals within the UKBB underwent ophthalmic evaluations, encompassing the quantification of corneal-compensated IOP measured in millimeters of mercury (mmHg) through the application of a non-contact tonometer. The preference for corneal-compensated IOP as the selected exposure phenotype over regular IOP stems from its designed capability to accommodate corneal biomechanical properties. Moreover, this metric has been previously employed in GWAS for IOP [ 18 ]. Summary data for RVO were obtained from the FinnGen consortium database’s R9 release, involving a cohort of 775 documented RVO cases and 376,502 controls [ 14 ]. The FinnGen GWAS participants were of European origin, and the identification of RVO cases adhered to diagnostic criteria established in either the International Classification of Diseases, Ninth Revision (ICD-9) or International Statistical Classification of Diseases, Tenth Revision (ICD-10) code. For the mediation MR analysis, GWAS summary statistics for POAG were also retrieved from a GWAS meta-analysis of 16,677 POAG cases and 199,580 controls of European ancestry [ 19 ] (Supplementary Table 1). POAG cases were selected based on the ICD-9/ICD-10 criteria [ 19 ]. The genotyping, quality control, and imputation methods applied to the GWAS data used in our study have been detailed elsewhere [ 14 , 19 , 20 ]. Selection of genetic variants as instrumental variables We selected SNPs identified in the IOP GWAS that achieved genome-wide significance (P-value < 5*10 -8 ), following clumping for linkage disequilibrium (LD) at r 2 < 0.001 over a 10mb window [ 17 ]. The MR-Steiger directionality test was employed to discern the directions of the causal relationship between IOP and RVO [ 21 ]. SNPs exhibiting a stronger correlation with the outcome than the exposure were systematically excluded, as were those demonstrating notable influence in the funnel plots and scatter plots. Ultimately, we identified 47 SNPs associated with IOP as instrumental variables. Furthermore, we quantified the proportion of variability in IOP that is explained by these 47 SNPs, through the summation of the coefficients of determination (R 2 ) derived from the associations between the selected SNPs and IOP. Employing a similar approach, we chose 50 SNPs from the GWAS on POAG for utilization in our MR mediation analysis. Statistical analysis After the process of data harmonization based on HapMap3 [ 22 ], the removal of strand-ambiguous variants, and the alignment of association estimates, we proceeded to calculate Wald ratios. We calculated these ratios by dividing the logarithm of odds ratio per allele for each SNP identified in the GWAS for RVO by its respective estimate obtained from the IOP GWAS. Subsequently, the total effect of IOP on RVO risk was assessed through a multiplicative random effects inverse variance weighted (IVW) meta-analysis of the Wald ratios [ 23 ]. We implemented a univariable two-sample MR methodology, utilizing summary data extracted from GWAS of IOP and RVO. Our two-sample MR analysis was based on three key assumptions. Firstly, the imperative that the genetic variants selected should be associated with the targeted risk factor, a principle known as the “relevance” assumption [ 24 ]. We adhered to this assumption by choosing SNPs as instrumental variables that reached the threshold of genome-wide significance (P-value < 5*10 -8 ). Furthermore, we evaluated the strength of our instrumental variables by assessing the F-statistic of the selected SNPs, concurrently examining the proportion of variance in the exposure that they accounted for [ 25 ]. Secondly, the selected genetic instruments should not be correlated with factors that could potentially confound the relationship between the exposure and outcome, also known as the “exchangeability” assumption [ 24 ]. Thirdly, the “exclusion restriction” assumption [ 24 ], where the genetic instruments should not influence the outcome except through the risk factor of interest. While the “exchangeability” and “exclusion restriction” assumptions are inherently unverifiable, we executed sensitivity analyses to identify potential violations of these MR assumptions. PhenoScanner [ 26 ] was employed to assess whether any of our selected genetic instruments were associated with phenotypes that could serve as potential confounders in our analysis. In instances where pleiotropic pathways were identified, we applied multivariable MR to adjust for these effects [ 27 ]. Moreover, we examined heterogeneity among our selected SNPs using the Cochran Q heterogeneity test and IGX 2 [ 28 ] to detect pleiotropy, utilized several MR methods proposed to enhance robustness in instances where genetic variants exhibit pleiotropy [ 29 ] (MR Egger regression, penalized weighted median, IVW radial regression, and MR-Pleiotropy Residual Sum and Outlier (MR-PRESSO)) and performed a leave-one-out analysis to ascertain if the IVW estimate was influenced by a singular SNP. To further explore the mediating effect of IOP on RVO through POAG, we conducted a two-step MR for mediation analysis, as outlined in previous literature [ 15 ]. This process involves the calculation of two distinct MR estimates: firstly, calculating the effect of IOP on POAG through a univariable MR model, and secondly, assessing the effect of POAG on RVO using a multivariable MR model that includes an adjustment for IOP. The IOP estimate on POAG multiplied with the adjusted POAG estimate on RVO provided the indirect effect of IOP on RVO, mediated through POAG. We also quantified the proportion of the total effect of IOP on RVO explained by the mediator (POAG) by dividing the calculated indirect effect of IOP on RVO by the total effect (adjusted direct and indirect effects of IOP on RVO). The delta method was employed to derive 95% confidence intervals (95%CI) for the indirect effect [ 30 ]. Furthermore, considering the inherent requirement for independence between the SNPs selected as instruments for the exposure (IOP) and mediator (POAG) in MR mediation analysis, we ensured that the selected SNPs from the IOP and POAG GWAS datasets were non-overlapping [ 15 ]. This precautionary step was essential to uphold the validity and integrity of the MR for mediation analysis, preventing potential bias arising from SNP overlap. We conducted all analyses with R version 4.2.1 [ 31 ] using the MVMR (0.3), TwoSampleMR (0.5.6), MendelianRandomization (0.5.1), MRPRESSO (1.0) and cause (1.2.0) packages. Results The 47 SNPs selected from the IOP GWAS (Supplementary Figure 1) accounted for 2.64% of the variability in IOP, with all SNPs displaying F-statistics of ≥29.96 (Supplementary Table 2). Employing the IVW method, genetically predicted IOP was positively associated with RVO risk (OR = 1.53 per 1mmHg increase in IOP; 95%CI = 1.04 to 2.26; P-value = 0.03) ( Figure 2 and Supplementary Figure 2). Results from pleiotropy-robust MR methods aligned with the IVW analysis estimate ( Figure 2 ). We found no associations of our selected SNPs with RVO risk factors apart from POAG (Supplementary Table 3) and, thus, we refrained from conducting multivariable MR to correct for potential correlated horizontal pleiotropy. Download figure Open in new tab Figure 1. Mendelian randomization estimates for the effect of intraocular pressure on retinal vein occlusion. Estimates are reported as changes in odds of retinal vein occlusion per 1 mmHg increase in intraocular pressure. SNP, single nucleotide polymorphism; CI, confidence interval; MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier. Download figure Open in new tab Figure 2. Directed acyclic graphs of the mediation analysis with Mendelian randomization. The indirect effect of intraocular pressure (IOP) on retinal vein occlusion (RVO) can be calculated by multiplying α times β, where α is the effect of IOP on primary open-anlge glaucoma (POAG), and β the effect of POAG on RVO. The proportion mediated can be estimated by dividing the indirect effect by the total effect of IOP on RVO. All estimates are shown as the difference in the logarithm of odds of the outcome, per 1 unit increase of the exposure (continuous variables: IOP) or between the two exposure groups (binary variables: POAG). The Wald ratios for IOP with RVO did not exhibit heterogeneity (Supplementary Table 4) and the Cochran’s Q heterogeneity test yielded a value of 37.89 (p-value = 0.797). Additionally, no evidence for directional pleiotropy was found, since the intercepts from the MR-Egger analyses did not deviate from zero (Supplementary Table 4). The leave-one-SNP-out analysis did not identify any SNPs exerting significant influence on the IVW estimate for the association between IOP and RVO risk (Supplementary Table 5). An illustrative depiction of the MR mediation analysis can be seen in Figure 2 . We found that 91.6% of the total effect of IOP on RVO was mediated through POAG (Supplementary Table 6). Discussion In this MR mediation analysis, we utilized genetic data to disentangle the causal pathway between IOP and RVO, and to also explore the mediating role of POAG in this relationship. Our results support a positive causal effect of IOP on RVO which is mostly mediated through POAG. Many studies have examined the relationship between glaucoma, specifically POAG, and risk for RVO [ 4 ]. Yin et al demonstrated in their meta-analysis that POAG is significantly associated with RVO risk (OR: 5.03; 95% CI: 3.97 to 6.37) [ 3 ]. Moreover, in their subgroup analysis, POAG and chronic open angle-glaucoma were also correlated with central RVO (OR: 13.30; 95% CI: 3.34 to 53.20) and branch RVO (OR: 2.14; 95% CI: 1.09 to 4.20), but conclusions are guarded due to the small number of studies included in the meta-analyses. Similarly, Na et al. determined that the RVO incidence rate for the open-angle glaucoma patients was 528.95 per 100,000 person-years (95% CI: 515.46 to 542.79), 3.27 times higher than in the general population [ 32 ]. Another Korean nationwide, population-based 11-year longitudinal study showed that patients with POAG were under higher risk for developing RVO (HR: 3.27, 95% CI: 2.55-4.19) [ 33 ], while in the Beaver Dam Eye Study, no significant difference was found in the age-adjusted 5-year incidence of RVO between participants with and without POAG [ 34 ]. Different potential mechanisms have been proposed to interpret the association between POAG and RVO, and IOP has been shown to be a major contributing factor affecting retinal hemodynamics. It has been suggested that increased IOP may compress blood vessels, leading to proliferation of the vein intima and subsequently to the collapse of retinal vessel walls [ 5 ]. The effects of this mechanical compression are more pronounced at the level of the lamina cribrosa where the venous pressure is lowest. In older individuals with impaired retinal circulation and compromised autoregulatory vascular mechanisms, even a slight increase in IOP could be significant enough to decrease ocular blood flow, leading to venous slowing or stasis [ 6 ]. Luntz et al. strongly support an etiological relationship between increased IOP and central RVO having shown that the incidence of RVO in ocular hypertension is similar to that of POAG [ 35 ]. However, the Ocular Hypertension Treatment study, which included only participants with ocular hypertension, did not find a significant difference in RVO incidence between the observation and medication group [ 36 ]. The effect of IOP on retinal vessels might be further accentuated by anatomical characteristics of the glaucomatous optic nerve head. The posterior bowing of the lamina cribrosa present in glaucomatous optic discs confers less glial support to the retinal vessels thus making them more vulnerable to changes in IOP [ 8 ]. Also, higher cup-to-disc ratio has been correlated with the development of RVO, postulating that optic disc cupping induces structural abnormalities to the vessel at the disc, thus altogether making them susceptible to occlusion, venous stasis, and, consequently, to thrombosis, according to Virchow’s triad [ 7 , 9 , 10 ]. On top of that, according to the vascular theory of glaucoma pathogenesis, patients with POAG have narrower retinal arteries and veins compared to normal individuals, predisposing them to greater vessel wall compression following a rise in IOP [ 37 ]. Moreover, it has been hypothesized that RVO and POAG could potentially arise from the same vascular abnormality given their shared vascular risk factors. Optic disc hemorrhages frequently encountered in POAG and RVO indicate small vein occlusions, raising the question of a common underlying pathway [ 38 ]. The key strength of this study is that it stands out as the first study to utilize the approach of MR mediation analysis to disentangle the complex causal pathway between IOP, POAG and RVO. Another strength of it lies in the consistency of association estimates obtained from pleiotropy-robust methods, aligning with the IVW estimate and indicating an absence of model violations. Nevertheless, certain limitations warrant consideration. Firstly, our investigation primarily focused on delineating the mediating role of POAG in the association between IOP and RVO, excluding exploration into other glaucoma subtypes such as primary angle-closure glaucoma. Secondly, our MR models assumed a linear relationship between the identified risk factors and the observed outcomes, whereas the true association among IOP, POAG, and RVO may manifest as non-linear. Thirdly, it is imperative to exercise caution when extrapolating the genetic associations derived from European populations to other ethnic groups due to potential population-specific variations. Fourthly, we did not assess the associations of IOP with the two main subtypes of RVO, namely branch RVO and central RVO, since no GWAS datasets are available for these phenotypes. In conclusion, our MR analyses infer an elevated risk of RVO associated with higher IOP. By utilizing MR mediation analysis, our findings additionally suggest that the majority of the effect of IOP on RVO risk is mediated through POAG. To advance our understanding of this intricate interplay, further investigations through population-based prospective studies, as well as experimental studies, are warranted to comprehensively explore the complex pathway involving IOP, POAG, and RVO. Author contributions KA, BSE and NM contributed to the study conception and design, drafted the manuscript and analyzed the data. All authors critically revised the manuscript for important intellectual content, provided administrative, technical, or material support and approved the final version. Data availability The summary statistics for the intraocular for the UKBB GWAS are available at https://pan.ukbb.broadinstitute.org (access date: 2023/10/12). The retinal vein occlusion summary statistics for the FinnGen GWAS are available at https://www.finngen.fi/en/access_results (access date: 2023/10/12). The primary open-angle glaucoma summary statistics are available at https://www.ebi.ac.uk/gwas/publications/33627673 (access date: 2023/10/12) Additional information Competing Interests None declared Data Availability All data produced in the present work are contained in the manuscript Acknowledgments We want to acknowledge the participants and investigators of the FinnGen and UK Biobank studies. References 1. ↵ Mitchell , P. , Smith , W. & Chang , A . Prevalence and Associations of Retinal Vein Occlusion in Australia: The Blue Mountains Eye Study . Arch. Ophthalmol . 114 , 1243 – 1247 , doi: 10.1001/archopht.1996.01100140443012 ( 1996 ). OpenUrl CrossRef PubMed Web of Science 2. ↵ Kolar , P . Risk factors for central and branch retinal vein occlusion: a meta-analysis of published clinical data . J Ophthalmol 2014 , 724780 , doi: 10.1155/2014/724780 ( 2014 ). OpenUrl CrossRef PubMed 3. ↵ Yin , X. , Li , J. , Zhang , B. & Lu , P . Association of glaucoma with risk of retinal vein occlusion: A meta-analysis . Acta Ophthalmol. (Copenh .) 97 , 652 – 659 , doi: 10.1111/aos.14141 ( 2019 ). OpenUrl CrossRef 4. ↵ Jabbehdari , S. , Yazdanpanah , G. , Cantor , L. B. & Hajrasouliha , A. R . A narrative review on the association of high intraocular pressure and glaucoma in patients with retinal vein occlusion . Ann Transl Med 10 , 1072 , doi: 10.21037/atm-22-2730 ( 2022 ). OpenUrl CrossRef 5. ↵ Verhoeff FH . The effect of chronic glaucoma on the central retinal vessels . Arch Ophthalmol . 1913 ; 42 : 145 – 52 . OpenUrl 6. ↵ Frucht , J. , Shapiro , A. & Merin , S . Intraocular pressure in retinal vein occlusion . Br. J. Ophthalmol . 68 , 26 – 28 , doi: 10.1136/bjo.68.1.26 ( 1984 ). OpenUrl Abstract / FREE Full Text 7. ↵ Beaumont , P. E. & Kang , H. K . Cup-to-disc ratio, intraocular pressure, and primary open-angle glaucoma in retinal venous occlusion 1 . Ophthalmology 109 , 282 – 286 , doi: 10.1016/S0161-6420(01)00922-8 ( 2002 ). OpenUrl CrossRef PubMed 8. ↵ Behrman , S. RETINAL VEIN OBSTRUCTION . Br. J. Ophthalmol . 46 , 336 – 342 , doi: 10.1136/bjo.46.6.336 ( 1962 ). OpenUrl FREE Full Text 9. ↵ Dobree , J. H . Venous obstruction and neovascularization at the disc in chronic glaucoma . Trans. Ophthalmol. Soc. U. K . 77 , 229 – 237 ; discussion 237-228 ( 1957 ). OpenUrl 10. ↵ Rubinstein , K. & Jones , E. B . Retinal vein occlusion: long-term prospects: 10 years’ follow-up of 143 patients . Br. J. Ophthalmol . 60 , 148 – 150 , doi: 10.1136/bjo.60.2.148 ( 1976 ). OpenUrl FREE Full Text 11. ↵ Sanderson , E. et al. Mendelian randomization . Nature Reviews Methods Primers 2 , 6 , doi: 10.1038/s43586-021-00092-5 ( 2022 ). OpenUrl CrossRef PubMed 12. ↵ Zixun Wang , Zhiqing Li . Causal Association between Glaucoma and Risk of Retinal Vascular Occlusion: A Mendelian Randomization Study , 10 August 2023, PREPRINT (Version 1) available at Research Square , doi: 10.21203/rs.3.rs-3223120/v1 . OpenUrl CrossRef 13. ↵ Pan-UKB team. https://pan.ukbb.broadinstitute.org . 2020 . 14. ↵ Kurki , M. I. et al. FinnGen provides genetic insights from a well-phenotyped isolated population . Nature 613 , 508 – 518 , doi: 10.1038/s41586-022-05473-8 ( 2023 ). OpenUrl CrossRef PubMed 15. ↵ Carter , A. R. et al. Mendelian randomisation for mediation analysis: current methods and challenges for implementation . Eur. J. Epidemiol . 36 , 465 – 478 , doi: 10.1007/s10654-021-00757-1 ( 2021 ). OpenUrl CrossRef PubMed 16. ↵ Skrivankova , V. W. et al. Strengthening the reporting of observational studies in epidemiology using mendelian randomisation (STROBE-MR): explanation and elaboration . BMJ 375 , n2233 , doi: 10.1136/bmj.n2233 ( 2021 ). OpenUrl FREE Full Text 17. ↵ Burgess , S. et al. Guidelines for performing Mendelian randomization investigations: update for summer 2023 [version 3; peer review: 2 approved] . Wellcome Open Research 4 , doi: 10.12688/wellcomeopenres.15555.3 ( 2023 ). OpenUrl CrossRef 18. ↵ Simcoe , M. J. , Khawaja , A. P. , Hysi , P. G. & Hammond , C. J . Genome-wide association study of corneal biomechanical properties identifies over 200 loci providing insight into the genetic etiology of ocular diseases . Hum. Mol. Genet . 29 , 3154 – 3164 , doi: 10.1093/hmg/ddaa155 ( 2020 ). OpenUrl CrossRef 19. ↵ Gharahkhani , P. et al. Genome-wide meta-analysis identifies 127 open-angle glaucoma loci with consistent effect across ancestries . Nature Communications 12 , 1258 , doi: 10.1038/s41467-020-20851-4 ( 2021 ). OpenUrl CrossRef 20. ↵ Bycroft , C. et al. The UK Biobank resource with deep phenotyping and genomic data . Nature 562 , 203 – 209 , doi: 10.1038/s41586-018-0579-z ( 2018 ). OpenUrl CrossRef PubMed 21. ↵ Hemani , G. , Tilling , K. & Davey Smith , G . Orienting the causal relationship between imprecisely measured traits using GWAS summary data . PLoS Genet . 13 , e1007081 , doi: 10.1371/journal.pgen.1007081 ( 2017 ). OpenUrl CrossRef PubMed 22. ↵ Frazer , K. A. et al. A second generation human haplotype map of over 3.1 million SNPs . Nature 449 , 851 – 861 , doi: 10.1038/nature06258 ( 2007 ). OpenUrl CrossRef PubMed Web of Science 23. ↵ Burgess , S. , Dudbridge , F. & Thompson , S. G . Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods . Stat. Med . 35 , 1880 - 1906 ( 2015 ). 24. ↵ Labrecque , J. & Swanson , S. A . Understanding the Assumptions Underlying Instrumental Variable Analyses: a Brief Review of Falsification Strategies and Related Tools . Curr Epidemiol Rep 5 , 214 – 220 , doi: 10.1007/s40471-018-0152-1 ( 2018 ). OpenUrl CrossRef PubMed 25. ↵ Lawlor , D. A. , Harbord , R. M. , Sterne , J. A. , Timpson , N. & Davey Smith , G . Mendelian randomization: using genes as instruments for making causal inferences in epidemiology . Stat. Med . 27 , 1133 – 1163 , doi: 10.1002/sim.3034 ( 2008 ). OpenUrl CrossRef PubMed 26. ↵ Kamat , M. A. et al. PhenoScanner V2: an expanded tool for searching human genotype- phenotype associations . Bioinformatics 35 , 4851 – 4853 , doi: 10.1093/bioinformatics/btz469 ( 2019 ). OpenUrl CrossRef PubMed 27. ↵ Sanderson , E. Multivariable Mendelian Randomization and Mediation . Cold Spring Harb. Perspect. Med . 11 , doi: 10.1101/cshperspect.a038984 ( 2021 ). OpenUrl Abstract / FREE Full Text 28. ↵ Hemani , G. , Bowden , J. & Davey Smith , G . Evaluating the potential role of pleiotropy in Mendelian randomization studies . Hum. Mol. Genet . 27 , R195 – R208 , doi: 10.1093/hmg/ddy163 ( 2018 ). OpenUrl CrossRef PubMed 29. ↵ Slob , E. A. W. & Burgess , S . A comparison of robust Mendelian randomization methods using summary data . Genet. Epidemiol . 44 , 313 – 329 , doi: 10.1002/gepi.22295 ( 2020 ). OpenUrl CrossRef 30. ↵ Ogasawara , H . Asymptotic standard errors of estimated standard errors in structural equation modelling . Br. J. Math. Stat. Psychol . 55 , 213 – 229 , doi: 10.1348/000711002760554552 ( 2002 ). OpenUrl CrossRef PubMed 31. ↵ R Core Team ( 2022 ). R : A language and environment for statistical computing . R Foundation for Statistical Computing, Vienna, Austria . URL https://www.R-project.org/ . 32. ↵ Na , K. I. , Jeoung , J. W. , Kim , Y. K. , Lee , W. J. & Park , K. H . Incidence of retinal vein occlusion in open-angle glaucoma: a nationwide, population-based study using the Korean Health Insurance Review and Assessment Database . Clin. Experiment. Ophthalmol . 46 , 637 – 644 , doi: 10.1111/ceo.13157 ( 2018 ). OpenUrl CrossRef 33. ↵ Park , H. L. , Jung , Y. , Han , K. , Lee , M. Y. & Park , C. K. Health care claims for primary open-angle glaucoma and retinal vein occlusion from an 11-year nationwide dataset . Sci. Rep . 7 , 8038 , doi: 10.1038/s41598-017-07890-6 ( 2017 ). OpenUrl CrossRef 34. ↵ Klein , R. , Klein , B. E. , Moss , S. E. & Meuer , S. M . The epidemiology of retinal vein occlusion: the Beaver Dam Eye Study . Trans. Am. Ophthalmol. Soc . 98 , 133 – 141 ; discussion 141-133 ( 2000 ). OpenUrl PubMed 35. ↵ Luntz , M. H. & Schenker , H. I . Retinal vascular accidents in glaucoma and ocular hypertension . Surv. Ophthalmol . 25 , 163 – 167 , doi: 10.1016/0039-6257(80)90093-4 ( 1980 ). OpenUrl CrossRef PubMed Web of Science 36. ↵ Barnett , E. M. , Fantin , A. , Wilson , B. S. , Kass , M. A. & Gordon , M. O . The incidence of retinal vein occlusion in the ocular hypertension treatment study . Ophthalmology 117 , 484 – 488 , doi: 10.1016/j.ophtha.2009.08.022 ( 2010 ). OpenUrl CrossRef PubMed 37. ↵ Sonnsjö , B. & Krakau , C. E . Arguments for a vascular glaucoma etiology . Acta Ophthalmol. (Copenh .) 71 , 433 – 444 , doi: 10.1111/j.1755-3768.1993.tb04615.x ( 1993 ). OpenUrl CrossRef PubMed 38. ↵ Hayreh , S. S . Progress in the understanding of the vascular etiology of glaucoma . Curr. Opin. Ophthalmol . 5 ( 1994 ). View the discussion thread. Back to top Previous Next Posted February 01, 2024. Download PDF Supplementary Material Data/Code Email Thank you for your interest in spreading the word about medRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Intraocular pressure, primary open-angle glaucoma and the risk of retinal vein occlusion: a Mendelian randomization mediation analysis Message Subject (Your Name) has forwarded a page to you from medRxiv Message Body (Your Name) thought you would like to see this page from the medRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Intraocular pressure, primary open-angle glaucoma and the risk of retinal vein occlusion: a Mendelian randomization mediation analysis Andreas Katsimpris , Sebastian-Edgar Baumeister , Nafsika Voulgari , Hansjörg Baurecht , Stylianos Kandarakis , Michael Nolde medRxiv 2024.01.31.24302094; doi: https://doi.org/10.1101/2024.01.31.24302094 Share This Article: Copy Citation Tools Intraocular pressure, primary open-angle glaucoma and the risk of retinal vein occlusion: a Mendelian randomization mediation analysis Andreas Katsimpris , Sebastian-Edgar Baumeister , Nafsika Voulgari , Hansjörg Baurecht , Stylianos Kandarakis , Michael Nolde medRxiv 2024.01.31.24302094; doi: https://doi.org/10.1101/2024.01.31.24302094 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Ophthalmology Subject Areas All Articles Addiction Medicine (573) Allergy and Immunology (865) Anesthesia (303) Cardiovascular Medicine (4457) Dentistry and Oral Medicine (445) Dermatology (383) Emergency Medicine (610) Endocrinology (including Diabetes Mellitus and Metabolic Disease) (1517) Epidemiology (15244) Forensic Medicine (30) Gastroenterology (1132) Genetic and Genomic Medicine (6620) Geriatric Medicine (669) Health Economics (1002) Health Informatics (4557) Health Policy (1372) Health Systems and Quality Improvement (1615) Hematology (543) HIV/AIDS (1272) Infectious Diseases (except HIV/AIDS) (15936) Intensive Care and Critical Care Medicine (1106) Medical Education (624) Medical Ethics (147) Nephrology (670) Neurology (6634) Nursing (346) Nutrition (999) Obstetrics and Gynecology (1148) Occupational and Environmental Health (957) Oncology (3348) Ophthalmology (980) Orthopedics (369) Otolaryngology (421) Pain Medicine (436) Palliative Medicine (130) Pathology (665) Pediatrics (1696) Pharmacology and Therapeutics (693) Primary Care Research (714) Psychiatry and Clinical Psychology (5463) Public and Global Health (9257) Radiology and Imaging (2210) Rehabilitation Medicine and Physical Therapy (1371) Respiratory Medicine (1198) Rheumatology (598) Sexual and Reproductive Health (716) Sports Medicine (532) Surgery (714) Toxicology (99) Transplantation (289) Urology (265) (function(){function c(){var b=a.contentDocument||a.contentWindow.document;if(b){var d=b.createElement('script');d.innerHTML="window.__CF$cv$params={r:'a031838a0a5f8e2e',t:'MTc4MDAxNTkxOQ=='};var a=document.createElement('script');a.src='/cdn-cgi/challenge-platform/scripts/jsd/main.js';document.getElementsByTagName('head')[0].appendChild(a);";b.getElementsByTagName('head')[0].appendChild(d)}}if(document.body){var a=document.createElement('iframe');a.height=1;a.width=1;a.style.position='absolute';a.style.top=0;a.style.left=0;a.style.border='none';a.style.visibility='hidden';document.body.appendChild(a);if('loading'!==document.readyState)c();else if(window.addEventListener)document.addEventListener('DOMContentLoaded',c);else{var e=document.onreadystatechange||function(){};document.onreadystatechange=function(b){e(b);'loading'!==document.readyState&&(document.onreadystatechange=e,c())}}}})();

Text is read by the "Ask this paper" AI Q&A widget below. Extraction quality varies by source — PMC NXML preserves structure cleanly, OA-HTML may include some navigation residue, and OA-PDF can have broken hyphenation. The publisher copy (via DOI) is the canonical version.

My notes (saved in your browser only)

Ask this paper AI returns verbatim quotes from the full text · source: preprint-html

Answers must be backed by verbatim quotes from this paper's full text. Hallucinated quotes are dropped automatically; if no verbatim passage answers the question, we say so. How this works

Citation neighborhood (no data yet)

We don't have any in-corpus citations linked to this paper yet. This is a recent paper (2024) — citers typically take a year or two to land, and the OpenAlex reference graph may still be filling in.

Source provenance

europepmc
last seen: 2026-05-20T01:45:00.602351+00:00
unpaywall
last seen: 2026-05-22T02:00:06.705733+00:00
License: CC-BY-NC-ND-4.0