Integrative genetic analysis reveals new relationships between intraocular pressure, glaucoma and ischemic stroke risk: A study based on combined SNP-to-gene, mendelian randomization and pathway investigations
=================================================================================================================================================================================================================

* Qi Zhang
* Di Hu
* Zenan Lin
* Junhong Jiang

## Abstract

**Background** Ischemic stroke (IS) is a leading cause of death in elderly people. Previous studies on exploring the association between intraocular pressure (IOP), glaucoma and the IS risk provided inconsistent results and unclear elucidations.

**Methods** Here, multiple genetic approaches were employed to investigate the possible causality between these traits. First, we performed the traditional bidirectional mendelian randomization (MR) study to explore the causal relationship between IOP, glaucoma and IS. Second, the SNPs selected as instrumental variables for IOP and primary open-angle glaucoma (POAG) were mapped to relevant genes by the novel combined SNP-to-gene (cS2G) method. The genes with significant causal effects on IS were then introduced to the gene ontology (GO), pathway and colocalization analyses. Third, the partitioned heritability analysis was also performed to evaluate the genome complexity with the Linkage Disequilibrium Score (LDSC) tool. Fourth, we also performed single SNP mendelian randomization (SSMR) study to find the SNPs of IOP and glaucoma, which had significant causal influence on IS risk. Then, they were introduced to the cytogenetic investigation. The multiple variable MR (MVMR) was applied to assess the independence of the causal effect of the exposures.

**Results** The MR results supported the view that the elevated IOP and POAG may contribute to the IS risk, but not vice versa. By using the cS2G approach, we identified 31 and 3 genes which may play key roles in the IOP- and POAG-induced IS risk, respectively. The GO and pathway analyses indicated the olfactory pathway to be a crucial pathway in the IOP-associated mechanism. The colocalization study strengthened the causal implications of genes *CDKNA2A* and *CDKN2B-AS1* between POAG and IS. The partitioned heritability analysis showed that the most enriched categories for both IOP and POAG were regulatory-associated terms such as the Super-enhancer. The SSMR study demonstrated that the IOP-associated SNPs with causal effects on IS were located majorly in chromosomes 1 and 11, while the POAG-associated ones were mostly found in chromosomes 9 and 4. The MVMR confirmed that the causal effects of IOP and POAG were not independent from each other.

**Conclusions** This work provides novel evidences to support the causal implications between IOP, POAG and IS, and offered putative pathway and genes for managing IOP and POAG in IS.

Key words
*   intraocular pressure
*   ischemic stroke
*   causality
*   combined SNP-to-gene
*   mendelian randomization

## 1. Introduction

Ischemic stroke (IS) remains a significant global health concern, representing one of the major causes of disability and mortality 1. The pathogenesis of ischemic stroke involves the occlusion of cerebral blood vessels, resulting in the shortage of blood supply to affected brain areas and subsequent regional damage. While hypertension, diabetes mellitus, and smoking are acknowledged to contribute to stroke risk, emerging evidence suggest that factors beyond the conventional cardiovascular spectrum may also play a role in stroke pathophysiology 2, 3.

Recently, there has been growing interest in studying the potential association between intraocular pressure (IOP) or glaucoma, and IS 4, 5. IOP, a physical measure of the fluid pressure within the eye, is primarily regulated by the balance between aqueous humor production and outflow. Though the elevated IOP may not be detected in the normal tension glaucoma, it is a hallmark feature of other glaucoma types, and believed to be the leading cause for the progressive retinal ganglion cell loss and visual field defects 6. While the primary pathological consequence of elevated IOP is optic nerve damage and visual impairment, recent works have also indicated that elevated IOP may also be associated with systemic vascular dysfunction, including the alterations in cerebral blood flow dynamics 7.

In detail, the rationale for investigating the relationship between IOP and IS stems from the following aspects. First, both elevated IOP and IS share common risk factors, such as hypertension, diabetes mellitus, and metabolic syndrome 2, 8–10. Second, previous investigations have reported the anatomical and physiological connections between the ocular and cerebrovascular systems, indicating that alterations in IOP may contribute to the stroke risk through impacting the intracranial pressure 11. Third, previous population-based studies and clinical cohorts has provided preliminary insights into the plausible link between elevated IOP and increased stroke risk, although findings have been contradictory and warrant more in-depth investigations 12. Given the potential implications for stroke prevention and management, clarifying the relationship between IOP and ischemic stroke seems to be a promising research topic. A better understanding of the underlying mechanisms linking ocular and cerebrovascular health may provide novel insights into stroke pathophysiology and offer new therapeutic targets for stroke prevention and management.

Mendelian Randomization (MR) analysis is a powerful biomedical method that uses appropriate single nucleotide polymorphisms (SNPs) as instrumental variables (IVs) to identify the potential causal effects between the exposure and outcome traits. Since the conventional 2-sample MR was not able to provide the detailed portrait of the underlying mechanism, an up-to-date combined SNP-to-gene (cS2G) strategy13 was employed to demonstrate the key genes and pathways involved in the causality. By leveraging multiple genetic tools, herein, we investigated the causality between glaucoma, IOP and ischemic stroke, and further identified relevant genes and pathways for future biomarker research.

## 2. Methods

### 2.1 GWAS data source of intraocular pressure, glaucoma and ischemic stroke

Inn this work, numerous GWAS datasets for intraocular pressure, glaucoma and ischemic stroke were introduced for further analyses. The IOP (both of right and left eye) were assessed by the gold standard Goldmann method. Four glaucoma-associated datasets, i.e. primary open angle glaucoma (POAG), primary angle closure glaucoma (PACG), normal tension glaucoma (NTG) and undefined glaucoma) were also analysed in this study. Besides, a GWAS dataset on analysing IS was selected for further investigations (Malik, Chauhan et al. 2018). The GWAS summary statistics data of these datasets were obtained from the IEU OpenGWAS project ([https://gwas.mrcieu.ac.uk/](https://gwas.mrcieu.ac.uk/)). Their basic features (GWAS ID, sample size, ethnics and links) were summarized in Table S1.

### 2.2 Selection of genetic instrumental variables (IVs)

The summary statistics data of IOP contained 9851867 single nucleotide polymorphisms (SNPs) data. 85 SNPs with P values less than 5E-8 and valid Reference SNP identifications (rsIDs) were selected as significant SNPs (See supplementary Table S2.) By using R package TwoSampleMR (version 0.5.6) [31], we performed the linkage disequilibrium (LD) to acquire independent SNPs (r2 < 0.001 within 10Mb range). We set the parameter ‘pop’ to ‘EUR’ to use 1000 Genomes Project Phase 3 (EUR) as the reference panel. The absolute value of five F statistics (see Supplementary Table S2.) were all larger than 10 indicating sufficient statistical strength of the selected IVs [33]. The selection process of IVs for IS and various kinds of glaucoma subtypes were similar to the abovementioned procedure. Their IVs were summarized in Table S2-S7.

### 2.3 The univariable, bidirectional MR and multiple variable MR (MVMR)

The exposure and outcome datasets were harmonized. The proxy SNPs (R2 > 0.8) were obtained when no SNP was found in the outcome dataset. Then, five MR analytical methods (IVW (Inverse variance weighted, Weighted median, Weighted mode, MR Egger, Simple mode)) were employed to provide a comprehensive MR results. Each trait was introduced to the MR analysis as exposure and outcome, respectively. Besides, the MVMR analyses were performed on the traits which were found to have causal influence on the IS risk in the univariable MR investigations.

### 2.4 The Linkage disequilibrium score regression (LDSC) and partitioned heritability analysis

In the current work, the IOP of the right eye was selected for the MR study. However, the eyes are symmetric (namely right and left eye) organs, thus, there is also a GWAS dataset of the IOP of left eye. In order to prove that there is a strong correlation between the IOP of right and left eye, we performed a genetic correlation analysis between them with LDSC tool. The result indicates a strong genetic correlation (*P*= 3.9755e-86) between the IOP of right and left eye (See supplementary Table S8), which guarantees the sufficiency of the utilization of the GWAS data of right eye.

In the SNP partitioned heritability study, 96 previously developed baseline genomic molecular annotations (See Table S9; 1000G\_Phase3_baselineLD_v2.2_ldscores were downloaded from [https://doi.org/10.5281/zenodo.10515792](https://doi.org/10.5281/zenodo.10515792)) such as promoter, super-enhancers, conserved regions and etc 14. SNP-h2 partitioning implemented in the LDSC approach allowed us to evaluate simultaneously the functional annotations, and estimate their enrichments.

### 2.5 The cS2G, Gene ontology (GO) and pathway analyses

Under the guidance of the novel cS2G method 15, the genes correlated to the critical SNPs for IOP and glaucoma were identified. Then, they were introduced to the MR independently and those lead to significant MR causal effect on IS risk were selected for further colocalization, GO and pathway analyses.

### 2.6 Colocalization analysis

The genes of interest were introduced to the colocalization analysis. In detail, we performed the colocalization study in the specific region of each gene and tested whether identified relation between two phenotypical traits were driven by linkage disequilibrium within the given locus. Theoretically, the approximate Bayes Factor colocalization analyses was conducted under the following five hypotheses, 1. neither trait had an association in the given locus; 2. only trait 1 had an association in the locus; 3. only trait 2 had an association in the locus; 4. both traits had associations, but with different causal variants; 5. both traits had associations and shared one identical causal variant. In the end, five posterior probabilities (PH0, PH1, PH2, PH3, and PH4) were generated for each hypothesis 16.

### 2.7 Cytogenetic analysis

Using the SSMR (single SNP mendelian randomization), the single SNP was introduced to the MR separately, and those with a *P* value less than 0.05 were identified to be significant SNPs. Their cytogenetic information was investigated.

### 2.8 Statistical analysis

Using R packages TwoSampleMR (version 0.6.6),ieugwasr (1.0.1), gwasglue (0.0.0.9000), we performed the following analyses. The selected IVs’ summary statistics were harmonized with outcome GWAS dataset. The IVW method was regarded as the major analytical tool. As complementary sensitivity tests, the weighted median, weighted mode, simple mode, the MR-Egger regression, and leave-one-out methods were introduced to assess the heterogeneity and pleiotropy, and evaluate the robustness of identified causal associations. R packages coloc (5.2.3) and locuscomparer (1.0.0) were employed to conduct the colocalization analysis. The R package qqman (0.1.9) was used to demonstrate the mannhattan plot. The clusterProfiler (4.12.0) package was sued to perform the GO and pathway analysis. For the LDSC analysis, the LDSC (LD SCore v1.0.1) command tool was utilized under the Python (3.12.4) environment. A *P* value less than 0.05 was regarded as statistically significant.

## 3. Results

### 3.1 The IOP and POAG have causal effects on the IS risk

As shown in Figure 1., multiple genetic analytical approaches were employed to explore the associations between IOP, glaucoma and IS. In order to clarify the potential causality, each trait was investigated as exposure as well as outcome. However, no SNP in the PACG summary statistics had a significant *P* value less than 5E-8, it was analysed only as outcome in the current work. Standing in line with previously acknowledged view, the increase of IOP has an obvious causal effect on the risks of various kinds of glaucoma subtypes (i.e. primary open-angle glaucoma(POAG), primary angle-closure glaucoma(PACG), normal tension glaucoma (NTG)) (See Figure 2. and supplementary Table S10). Despite NTG, our result showed that the genetically proxied glaucoma, PACG and POAG also had causal influence on the IOP. (See Figure 2. and Table S11).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/26/2024.08.26.24312564/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2024/08/26/2024.08.26.24312564/F1)

Figure 1. 
The workflow for the study. Multiple genetic analytical tools, such as MR, cS2G, GO, colocalization, LDSC were employed to evaluate the associations between IOP, glaucoma and IS.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/26/2024.08.26.24312564/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2024/08/26/2024.08.26.24312564/F2)

Figure 2. 
The bidirectional MR conducted on IOP, various kinds of glaucoma types and IS. Scatter plots for MR analyses of the causal influence of IOP on IS (A), POAG on IS (B), IOP on glaucoma (C), IOP on POAG (D), POAG on IOP (F), IOP on PACG (G), IOP on NTG (H). The schema demonstrating the results of the bidirectional MR (E).

We detected no remarkable causal relationship between glaucoma, NTG, PACG and IS (See Figure 2, Table S12 and S13). Notably, our results indicated that the genetically proxied IOP and POAG had obvious causal effects on IS risk, but not vice versa (See Figure 2, Table S12-S15).

### 3.2 Using novel combined SNP-to-gene (cS2G) approach to identify significant genes

As the MR study suggested that the genetically predicted IOP and POAG had causal effects on IS risk, they were introduced into the further cS2G analyses. First, the SNPs with a *P* less than 5E-8 were selected (See Figure 3.A). Their corresponding genes were identified by searching the archive of the cS2G dataset. In order to obtain a relatively complete result, the cS2G datasets generated on two cohorts (i.e. UK Biobank and 1000 EUR genome) were simultaneously investigated. Finally, the corresponding genes from UK Biobank and 100 EUR genome were aggregated. By using this cS2G method, we identify 395 genes from 2406 SNPs of the IOP GWAS dataset. (see Table S16.) Among them, 378 were successfully introduced to the single gene MR on IS. The result indicated 31 genes to have significant causal effect on IS (*P* < 0.05). (See Table S17.) They were introduced into the Gene Ontology (cellular component, molecular function and biological process) and pathway analyses, and the result indicated the olfactory pathway to play a key role in the IOP-induced IS risk. Similarly, the same approach was introduced to study POAG. Three critical genes (i.e. *CDKN2A*, *CDKN2B-AS1* and *GAS7*) were found to potentially participate in the POAG-induced IS risk (See Table S18-S19).

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/26/2024.08.26.24312564/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2024/08/26/2024.08.26.24312564/F3)

Figure 3. 
A and D) The Manhattan plot demonstrating significant SNPs of IOP and POAG. B and E) By using the cS2G method, 378 and 36 genes were identified to be the IOP-and POAG-associated genes, respectively. The volcano plots highlighting the genes with significant causal effects on IS. C) The single gene MR analysis indicated *CDKN2A*, *CDKN2B-AS1* and *GAS7* to play key roles in the POAG-induced IS risk. The GO and pathway analyses with 31 IOP-associated genes indicated the olfactory pathway may contribute to the IOP-induced IS risk.

### 3.3 The Approximate Bayesian colocalization analyses

As shown in Figure 4., we gained strong evidences to support the 3rd Hypothesis for both *CDKN2A* (posterior probability = 0.92) and *CDKN2B-AS1* (posterior probability = 0.92), which means that POAG and IS had associations in the given locus of these two genes, however, with different causal variant. The study on *GAS7* demonstrated that only POAG had associations in the given region.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/26/2024.08.26.24312564/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2024/08/26/2024.08.26.24312564/F4)

Figure 4. 
The colocalization analyses on the POAG-associated genes *CDKN1A* (A), *CDKN2B-AS1* (B) and *GAS7* (C) and IS. The results demonstrated that both *CDKN2A* and *CDKN2B-AS1* had strong evidences (Posterior probability>0.8) to support the 3rd Hypothesis. However, *GAS7* lacked the evidence to support neither the 3rd nor 4th Hypothesis.

### 3.4 Partitioned heritability by multiple functional categories

In order to accurately articulate the complexity of the genome of POAG and IS, we partitioned their heritability into 96 functional categories, by using the baseline model provided by Finucane et. al.17. Intriguingly, both traits had statistically significant (less than Bonferoni-corrected *P*, i.e. 0.05/96) enrichments in three categories. In detail, SuperEnhancer\_HniszL2\_0 (Enrichment = 2.923, *P* = 0.0001), H3K27ac\_HniszL2\_0 (Enrichment =2.546, *P* = 1.51E-5) and H3K27ac\_PGC2L2_0 (Enrichment =4.288, *P* = 0.0001) were identified to be significant in POAG. Of notice, the SuperEnhancer_HniszL2_0 (Enrichment = 2.186, *P* = 1.26E-6) was also detected in ischemic stroke (See Figure 5. and Table S20-S21.).

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/26/2024.08.26.24312564/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2024/08/26/2024.08.26.24312564/F5)

Figure 5. 
The partitioned heritability analyses on POAG (A) and IS (B).

### 3.5. Single SNP mendelian randomization (SSMR)

SSMR revealed 278 and 507 IOP-associated variants with significant (*P* < 0.05) positive and negative causal estimates on IS, respectively (Figure 6A-B and Table S22). More than one third of the positive causal SNPs were on chromosomes 1, whereas over 50% of negative causal SNPs were on chromosomes 11.

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/26/2024.08.26.24312564/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2024/08/26/2024.08.26.24312564/F6)

Figure 6. 
The cytogenetic analyses of SNPs with significant causal effects on IS. (A and B) The plots demonstrating the top 10 IOP-associated SNPs with positive and negative causal influences on IS risk, respectively. (C and D) The pie charts showed the chromosome distribution of the significant IOP-associated SNPs with positive and negative causal estimates. (E and F) The plots demonstrating the top POAG-associated SNPs with positive and negative causal influences on IS risk, respectively. (G and H) The pie charts showed the chromosome distribution of the significant POAG-associated SNPs with positive and negative causal estimates.

In comparison to IOP, we found less POAG-associated variants with significant positive (n=81) and negative (n=8) causal estimates on IS. (Figure 6.E-F and Table S23). The positive causal SNPs were majorly located on chromosomes 9, and all of negative causal SNPs were on chromosomes 4.

### 3.6 MVMR

By performing the MVMR, we tested whether the POAG’ causal influence on IS was independent from IOP. Though both traits were proved to have significant causal effects on IS risk in the univariable MR, neither POAG nor IOP demonstrated remarkable influence in the MVMR (See Figure 7. and Table 24).

![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/26/2024.08.26.24312564/F7.medium.gif)

[Figure 7.](http://medrxiv.org/content/early/2024/08/26/2024.08.26.24312564/F7)

Figure 7. 
The forest plot of MVMR demonstrate that the causal effects on IS of IOP and POAG were not independent from each other.

## 4. Discussions

Ischemic stroke (IS) is a leading cause of disability and represents a major and growing public health problem worldwide. In this work, our integrative genetic analysis based on genome-wide association studies (GWAS) datasets revealed a significant positive genetic correlation of POAG with IS and confirmed for the first time there is the common genetic variation at both genome and transcriptome levels between the two disorders. Additionally, we further discovered that elevated IOP, a primary feature of glaucoma, may increase the risk of IS by the potential mechanistic link mediated through the olfactory pathway. The findings of this study provide compelling evidence to support a causal relationship between IOP, POAG and IS, which provide the novel insights into development of potential therapeutic targets for stroke prevention and treatment.

Glaucoma, as a potential risk factor for IS, has attracted increasing attention from researchers and clinicians. The correlation between glaucoma and stroke has been observed in previous observational studies. Ho and colleagues conducted a population-based, prospective cohort study with a large sample size of 4032 POAG patients and 20160 control participants from the Taiwan National Health Insurance Research Database, they found that 1.5-fold increased risk of stroke in POAG patients compared with population controls, when followed for 5 years18. A similar result in the US veteran population was observed by French et al.19, who found that the risk of stroke was higher in patients with POAG than in controls (6.4% vs 4.8%). Likewise, a recent meta-analysis involving 362267 participants showed an association of POAG with a higher risk of stroke5. Consistent with previous studies, our genetic analysis results show that POAG had obvious causal effects on IS risk.

While several studies found a correlation between POAG and IS, the exact mechanism and causal direction behind these observations have not yet been investigated clearly. In the present study, further analysis by co-localization showed that the variants in *CDKN1A* and *CDKN2B-AS1* two genes play key roles in the pathogenesis of POAG and IS. The cyclin-dependent kinase inhibitor 2A (*CDKN2A*) gene and the cyclin-dependent kinase inhibitor 2B antisense noncoding RNA (*CDKN2B-AS1*) gene are both located on chromosome 9p21, a region that has been found to be associated closely with a wide range of human diseases, including cancer20, cardiovascular disease21, 22 and neurodegenerative diseases23. *CDKN2A* gene encodes p16 and p14, two crucial cell cycle regulators, involved in cell cycle progression24, differentiation25, senescence26, and apoptosis27. While the *CDKN2B-AS1* gene, also known as *ANRIL*, transcribes into a long noncoding RNA in the antisense direction that is involved in modulating the nearby *CDKN2A/2B* genes by interacting with polycomb proteins, and subsequently participating in the alteration of cell cycle regulation28, 29. To date, the *CDKN2A* and *CDKN2B-AS1* genes has been reported to be associated with IS. Several human genetic research suggested that the *CDKN2A* gene polymorphism is associated with IS in different populations including Northern European30, Chinese31 and African Americans32 and West Africans33. Similarly, the genetic susceptibility of *CDKN2B-AS1* has been shown to be closely related to the occurrence and outcome of IS.34–36 On the other hand, *CDKN2A* and *CDKN2B-AS1* genes were also reported as the significantly mutated gene in the pathogenesis of POAG.29, 37, 38

Basic research also provided evidence for the potential roles of *CDKN2A* and *CDKN2B-AS1* in the pathogenesis of POAG and IS. In the animal models of glaucoma, multiple *CDKN2B-AS1* splice variants are observed in rat retina, furthermore, CDKN2A protein expression is significantly upregulated in rat retina, and this elevated expression corresponds to ongoing RGC death.29 *CDKN2A* encodes crucial cell cycle regulators that affect cell proliferation or senescence in RGCs by regulating cell cycle regulation, thereby contributing to POAG pathogenesis.39, 40 Mice with their *CDKN2B-AS*1 gene partially deleted show increased RGC vulnerability to apoptosis.40 Meanwhile, *CDKN2A* and *CDKN2B-AS1* were also shown to the regulate vascular smooth muscle cell proliferation28, and the accumulation of vascular smooth muscle cells can contribute to the vessel stenosis41, which in turn leads to an increased risk of IS.42 In the mice model of IS, the level of CDKN2B-AS1 exhibited a significant increase in the brain, and downregulation of CDKN2B-AS1 can reduce neuroinflammation by negatively regulating miR-671-5p to inhibit NF-κB.43 The latest GWAS results for POAG loci showed that *CDKN2A* and *CDKN2B-AS1* might be associated with vascular traits.37 These evidences may, partly, account for the genetic association between POAG and IS. A recent MR study found a positive causal effect value (OR=1.03) in the POAG-ischemic stroke association, however, the *P* value was above the 0.05 threshold.44 This might be owed to their relatively small sample size (n=63412, in comparison to the current dataset with a sample size of 214634 participants) of the POAG GWAS dataset, which may limit the statistical power of MR analysis.

Leveraging genome-wide association studies (GWAS) analysis, we identified significant genetic variants associated with both elevated IOP and increased risk of IS, supporting the notion of a shared genetic basis between these two phenotypes.45 Our results corroborate previous epidemiological studies suggesting an association between elevated IOP and ischemic stroke risk.46, 47 While the precise mechanisms underlying this association have remained elusive, our findings offer novel insights into the potential involvement of the olfactory pathway in mediating the observed relationship. Specifically, we identified genetic variants implicated in olfactory receptor activity and olfactory transduction pathways that were significantly associated with both IOP and ischemic stroke risk.48–52 This provides compelling evidence suggesting that alterations in olfactory function may represent a mechanistic link between elevated IOP and increased stroke risk.

The olfactory pathway serves as a critical interface between the external environment and the central nervous system, playing a fundamental role in odor detection and processing. Recent studies have highlighted the olfactory system’s role beyond sensory perception, implicating it in various physiological and pathological processes, including neurodegenerative diseases, cognitive decline, and cardiovascular health.53–58 Our findings extend this understanding by implicating olfactory dysfunction as a potential mediator linking elevated IOP to ischemic stroke risk. Several potential mechanisms may underlie the observed association between olfactory dysfunction and stroke risk. Firstly, olfactory dysfunction may serve as a marker of systemic vascular dysfunction, reflecting shared pathophysiological processes underlying both conditions.56, 59 Secondly, alterations in olfactory receptor activity and olfactory signaling pathways may directly impact cerebrovascular health through mechanisms such as inflammation, oxidative stress, and endothelial dysfunction.60–63 Finally, disruption of olfactory pathways may influence autonomic nervous system function, thereby modulating cardiovascular risk factors and contributing to stroke pathogenesis.56, 64

While our study provides evidence supporting the causal relationship between IOP, POAG and ischemic stroke risk, several limitations should be acknowledged. First, the observational nature of GWAS analysis precludes establishing causality definitively. Further mechanistic studies using animal models and experimental manipulations are warranted to validate our findings and elucidate underlying biological mechanisms. Second, previous studies have reported a potential association between the intracranial pressure (ICP) and IOP, thus, a mediated MR may help clarify whether the elevated IOP induces the IS risk through the ICP-associated pathway 65, 66. However, to our knowledge, a GWAS dataset on ICP has not been released yet. Third, the findings were concluded from the datasets with European participants, the generalizability of our results to diverse populations and ethnic groups warrants further investigation.

## 5. Conclusion

In conclusion, our study provides novel insights into the relationship between intraocular pressure, POAG, and ischemic stroke risk. By leveraging GWAS analysis, we have identified genetic variants implicating the olfactory pathway as a potential mediator of the elevated IOP-induced ischemic stroke risk, and the *CDKN2A* and *CDKN2B-AS1* genes of the observed association between POAG and ischemic stroke. Elucidating the mechanisms linking ocular and cerebrovascular health may offer new avenues for stroke prevention and management, highlighting the importance of interdisciplinary approaches in understanding complex disease pathogenesis.

## Supporting information

Supplementary Table S1-S24 [[supplements/312564_file10.zip]](pending:yes)

## Data Availability

The links of the GWAS data were described appropriately in the paper. The codes and detailed information required to replicate the results in this work are available from the corresponding authors upon reasonable request.

## Author Contributions

Conception, supervision and administration: Zenan Lin and Junhong Jiang; Data curation: Qi Zhang and Di Hu; Investigation: Qi Zhang and Di Hu; Methodology: Zenan Lin and Junhong Jiang; Writing – original draft: Qi Zhang and Di Hu; Writing – review & editing: Zenan Lin and Junhong Jiang.

## Funding

No fund was obtained for this work.

## Conflict of Interest

The authors declare no conflict of interests.

## Ethics Approval

The GWAS data were publicly available and approved by their original institutions. An ethics approval for the current work is not required.

## Availability of data and material

The links of the GWAS data were described appropriately in the paper. The codes and detailed information required to replicate the results in this work are available from the corresponding authors upon reasonable request.

## Acknowledgement

The authors appreciate the participants and investigators who made the GWAS datasets publicly available.

*   Received August 26, 2024.
*   Revision received August 26, 2024.
*   Accepted August 26, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/)

## References

1.  1.Feigin VL, Stark BA, Johnson CO, et al. Global, regional, and national burden of stroke and its risk factors, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019. The Lancet Neurology 2021; 20(10): 795–820.
    
    

2.  2.Boehme AK, Esenwa C, Elkind MSV. Stroke Risk Factors, Genetics, and Prevention. Circulation Research 2017; 120(3): 472–95.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6OToiMTIwLzMvNDcyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMjYvMjAyNC4wOC4yNi4yNDMxMjU2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

3.  3.Hameed S, Karim N, Wasay M, Venketasubramanian N. Emerging Stroke Risk Factors: A Focus on Infectious and Environmental Determinants. Journal of Cardiovascular Development and Disease, 2024. (accessed.
    
    

4.  4.He Q, Wang W, Xu D, et al. The Association of Glaucoma with Ischemic Stroke and Functional Outcome after Ischemic Stroke from the Perspective of Causality. Cerebrovascular Diseases 2024: 1–9.
    
    

5.  5.Wang M, Chen N, Sun BC, et al. Association between glaucoma and risk of stroke: A systematic review and meta-analysis. Frontiers in neurology 2022; 13: 1034976.
    
    

6.  6.Weinreb RN, Leung CKS, Crowston JG, et al. Primary open-angle glaucoma. Nature reviews Disease primers 2016; 2(1): 1–19.
    
    

7.  7.Tutaj M, Brown CM, Brys M, et al. Dynamic cerebral autoregulation is impaired in glaucoma. Journal of the Neurological Sciences 2004; 220(1): 49–54.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jns.2004.02.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15140605&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221633900007&link_type=ISI) 

8.  8.Yasukawa T, Hanyuda A, Yamagishi K, et al. Relationship between blood pressure and intraocular pressure in the JPHC-NEXT eye study. Scientific Reports 2022; 12(1): 17493.
    
    

9.  9.Hanyuda A, Sawada N, Yuki K, et al. Relationships of diabetes and hyperglycaemia with intraocular pressure in a Japanese population: the JPHC-NEXT Eye Study. Scientific Reports 2020; 10(1): 5355.
    
    

10. 10.Yi YH, Cho YH, Kim YJ, et al. Metabolic syndrome as a risk factor for high intraocular pressure: the Korea National Health and Nutrition Examination Survey 2008-2010. Diabetes Metab Syndr Obes 2019; 12: 131–7.
    
    

11. 11.Li Z, Yang Y, Lu Y, et al. Intraocular pressure vs intracranial pressure in disease conditions: A prospective cohort study (Beijing iCOP study). BMC Neurology 2012; 12(1): 66.
    
    

12. 12.Lin IC, Wang Y-H, Wang T-J, et al. Glaucoma, Alzheimer’s disease, and Parkinson’s disease: an 8-year population-based follow-up study. PloS one 2014; 9(10): e108938.
    
    

13. 13.Gazal S, Weissbrod O, Hormozdiari F, et al. Combining SNP-to-gene linking strategies to identify disease genes and assess disease omnigenicity. Nat Genet 2022; 54(6): 827–36.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-022-01087-y&link_type=DOI) 

14. 14.Bulik-Sullivan BK, Loh P-R, Finucane HK, et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nature Genetics 2015; 47(3): 291–5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3211&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

15. 15.Gazal S, Weissbrod O, Hormozdiari F, et al. Combining SNP-to-gene linking strategies to identify disease genes and assess disease omnigenicity. Nature Genetics 2022; 54(6): 827–36.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-022-01087-y&link_type=DOI) 

16. 16.Giambartolomei C, Vukcevic D, Schadt EE, et al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet 2014; 10(5): e1004383.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004383&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24830394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

17. 17.Finucane HK, Bulik-Sullivan B, Gusev A, et al. Partitioning heritability by functional annotation using genome-wide association summary statistics. Nature Genetics 2015; 47(11): 1228–35.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3404&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414678&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

18. 18.Ho JD, Hu CC, Lin HC. Open-angle glaucoma and the risk of stroke development: a 5-year population-based follow-up study. Stroke 2009; 40(8): 2685–90.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToic3Ryb2tlYWhhIjtzOjU6InJlc2lkIjtzOjk6IjQwLzgvMjY4NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzI2LzIwMjQuMDguMjYuMjQzMTI1NjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

19. 19.French DD, Margo CE. Open angle glaucoma and stroke. Ophthalmology 2010; 117(8): 1653.e3–4.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20682379&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

20. 20.Wang JZ, Patil V, Liu J, et al. Increased mRNA expression of CDKN2A is a transcriptomic marker of clinically aggressive meningiomas. Acta neuropathologica 2023; 146(1): 145–62.
    
    

21. 21.Hannou SA, Wouters K, Paumelle R, Staels B. Functional genomics of the CDKN2A/B locus in cardiovascular and metabolic disease: what have we learned from GWASs? Trends in endocrinology and metabolism: TEM 2015; 26(4): 176–84.
    
    

22. 22.Helgadottir A, Thorleifsson G, Manolescu A, et al. A common variant on chromosome 9p21 affects the risk of myocardial infarction. Science (New York, NY) 2007; 316(5830): 1491–3.
    
    

23. 23.Mori H, Funahashi Y, Yoshino Y, et al. Blood CDKN2A Gene Expression in Aging and Neurodegenerative Diseases. Journal of Alzheimer’s disease : JAD 2021; 82(4): 1737–44.
    
    

24. 24.Kamijo T, Zindy F, Roussel MF, et al. Tumor suppression at the mouse INK4a locus mediated by the alternative reading frame product p19ARF. Cell 1997; 91(5): 649–59.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0092-8674(00)80452-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9393858&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997YH96000012&link_type=ISI) 

25. 25.Monzon J, Liu L, Brill H, et al. CDKN2A mutations in multiple primary melanomas. The New England journal of medicine 1998; 338(13): 879–87.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM199803263381305&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9516223&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000072669800005&link_type=ISI) 

26. 26.Baker DJ, Childs BG, Durik M, et al. Naturally occurring p16(Ink4a)-positive cells shorten healthy lifespan. Nature 2016; 530(7589): 184–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature16932&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26840489&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

27. 27.Chakraborty S, Utter MB, Frias MA, Foster DA. Cancer cells with defective RB and CDKN2A are resistant to the apoptotic effects of rapamycin. Cancer letters 2021; 522: 164–70.
    
    

28. 28.Visel A, Zhu Y, May D, et al. Targeted deletion of the 9p21 non-coding coronary artery disease risk interval in mice. Nature 2010; 464(7287): 409–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature08801&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20173736&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275657100046&link_type=ISI) 

29. 29.Burdon KP, Macgregor S, Hewitt AW, et al. Genome-wide association study identifies susceptibility loci for open angle glaucoma at TMCO1 and CDKN2B-AS1. Nature genetics 2011; 43(6): 574–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.824&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21532571&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

30. 30.Wahlstrand B, Orho-Melander M, Delling L, et al. The myocardial infarction associated CDKN2A/CDKN2B locus on chromosome 9p21 is associated with stroke independently of coronary events in patients with hypertension. Journal of hypertension 2009; 27(4): 769–73.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/HJH.0b013e328326f7eb&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19293724&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265455800016&link_type=ISI) 

31. 31.Hu WL, Li SJ, Liu DT, et al. Genetic variants on chromosome 9p21 and ischemic stroke in Chinese. Brain research bulletin 2009; 79(6): 431–5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.brainresbull.2009.04.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19559344&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268614400016&link_type=ISI) 

32. 32.Carty CL, Keene KL, Cheng YC, et al. Meta-Analysis of Genome-Wide Association Studies Identifies Genetic Risk Factors for Stroke in African Americans. Stroke 2015; 46(8): 2063–8.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToic3Ryb2tlYWhhIjtzOjU6InJlc2lkIjtzOjk6IjQ2LzgvMjA2MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzI2LzIwMjQuMDguMjYuMjQzMTI1NjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

33. 33.Akinyemi R, Arnett DK, Tiwari HK, et al. Interleukin-6 (IL-6) rs1800796 and cyclin dependent kinase inhibitor (CDKN2A/CDKN2B) rs2383207 are associated with ischemic stroke in indigenous West African Men. Journal of the neurological sciences 2017; 379: 229–35.
    
    

34. 34.Ferreira LE, Secolin R, Lopes-Cendes I, Cabral NL, França PHC. Association and interaction of genetic variants with occurrence of ischemic stroke among Brazilian patients. Gene 2019; 695: 84–91.
    
    

35. 35.Fathy N, Kortam MA. Long Noncoding RNAs MALAT1 and ANRIL Gene Variants and the Risk of Cerebral Ischemic Stroke: An Association Study. ACS Chemical Neuroscience 2021; 12(8): 1351–62.
    
    

36. 36.Yang J, Gu L, Guo X, et al. LncRNA ANRIL Expression and ANRIL Gene Polymorphisms Contribute to the Risk of Ischemic Stroke in the Chinese Han Population. Cellular and Molecular Neurobiology 2018; 38(6): 1253–69.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10571-018-0593-6&link_type=DOI) 

37. 37. Lo Faro V, Bhattacharya A, Zhou W, et al. Novel ancestry-specific primary open-angle glaucoma loci and shared biology with vascular mechanisms and cell proliferation. Cell reports Medicine 2024; 5(2): 101430.
    
    

38. 38.Pasquale LR, Loomis SJ, Kang JH, et al. CDKN2B-AS1 genotype-glaucoma feature correlations in primary open-angle glaucoma patients from the United States. Am J Ophthalmol 2013; 155(2): 342–53.e5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajo.2012.07.023&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23111177&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

39. 39.Tawarayama H, Feng Q, Murayama N, Suzuki N, Nakazawa T. Cyclin-Dependent Kinase Inhibitor 2b Mediates Excitotoxicity-Induced Death of Retinal Ganglion Cells. Investigative ophthalmology & visual science 2019; 60(13): 4479–88.
    
    

40. 40.Teotia P, Van Hook MJ, Wichman CS, Allingham RR, Hauser MA, Ahmad I. Modeling Glaucoma: Retinal Ganglion Cells Generated from Induced Pluripotent Stem Cells of Patients with SIX6 Risk Allele Show Developmental Abnormalities. Stem Cells 2017; 35(11): 2239–52.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/stem.2675&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28792678&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

41. 41.Gibbons GH, Dzau VJ. The emerging concept of vascular remodeling. The New England journal of medicine 1994; 330(20): 1431–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJM199405193302008&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8159199&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1994NK97200008&link_type=ISI) 

42. 42.Bennett MR, Sinha S, Owens GK. Vascular Smooth Muscle Cells in Atherosclerosis. Circulation research 2016; 118(4): 692–702.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6OToiMTE4LzQvNjkyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMjYvMjAyNC4wOC4yNi4yNDMxMjU2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

43. 43.Deng L, Jiang J, Chen S, et al. Long Non-coding RNA ANRIL Downregulation Alleviates Neuroinflammation in an Ischemia Stroke Model via Modulation of the miR-671-5p/NF-κB Pathway. 2022; 47(7): 2002–15.
    
    

44. 44.Wang K, Lin X, Sheng S, Chen D, Liu X, Yao K. Association between glaucoma and stroke: A bidirectional mendelian randomization study. Advances in ophthalmology practice and research 2024; 4(3): 147–55.
    
    

45. 45.Choquet H, Thai KK, Yin J, et al. A large multi-ethnic genome-wide association study identifies novel genetic loci for intraocular pressure. Nature Communications 2017; 8(1): 2108.
    
    

46. 46.Su CW, Chang YC, Lin CL. Association of Neovascular Glaucoma with Risk of Stroke: A Population-Based Cohort Study. Journal of Ophthalmology 2017; 2017: 1851568.
    
    

47. 47.Liu Y, Yao B. Glaucoma in rural China (the Rural Epidemiology for Glaucoma in China (REG-China)): a national cross-sectional study. British Journal of Ophthalmology 2023; 107(10): 1458–66.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJqb3BodGhhbG1vbCI7czo1OiJyZXNpZCI7czoxMToiMTA3LzEwLzE0NTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOC8yNi8yMDI0LjA4LjI2LjI0MzEyNTY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

48. 48.Mainland JD, Keller A, Li YR, et al. The missense of smell: functional variability in the human odorant receptor repertoire. Nature neuroscience 2014; 17(1): 114–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.3598&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24316890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

49. 49.Wehling E, Naess H, Wollschlaeger D, et al. Olfactory dysfunction in chronic stroke patients. BMC neurology 2015; 15: 199.
    
    

50. 50.Zou J, Gao J, Sun L, et al. Perceived Taste and Olfactory Dysfunctions and Subsequent Stroke Risk. JACC Asia 2024; 4(6): 483–92.
    
    

51. 51.Gugleta K, Kochkorov A, Katamay R, et al. Olfactory function in primary open-angle glaucoma patients. Klinische Monatsblatter fur Augenheilkunde 2010; 227(4): 277–9.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20408073&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

52. 52.Dikmetas O, Aygün O, Kocabeyoglu S, et al. Smell Sensitivity in Primary Open-angle Glaucoma and Pseudoexfoliation Glaucoma. Journal of glaucoma 2022; 31(5): 300–4.
    
    

53. 53.Quail Z, Carter MM, Wei A, Li X. Management of cognitive decline in Alzheimer’s disease using a non-pharmacological intervention program: A case report. Medicine 2020; 99(21): e20128.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

54. 54.Roalf DR, Moberg MJ, Turetsky BI, et al. A quantitative meta-analysis of olfactory dysfunction in mild cognitive impairment. Journal of neurology, neurosurgery, and psychiatry 2017; 88(3): 226–32.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiam5ucCI7czo1OiJyZXNpZCI7czo4OiI4OC8zLzIyNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzI2LzIwMjQuMDguMjYuMjQzMTI1NjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

55. 55.Murphy C. Olfactory and other sensory impairments in Alzheimer disease. Nature reviews Neurology 2019; 15(1): 11–24.
    
    

56. 56.Orecchioni M, Kobiyama K, Winkels H. Olfactory receptor 2 in vascular macrophages drives atherosclerosis by NLRP3-dependent IL-1 production. Science 2022; 375(6577): 214–21.
    
    

57. 57.Mai Y, Zhang X, Li Z, et al. Olfaction is a Marker of Severity but Not Diagnosis in Anorexia Nervosa: a Systematic Review and Meta-Analysis. Neuropsychology review 2020; 30(2): 251–66.
    
    

58. 58.Siegel JK, Wroblewski KE, McClintock MK, Pinto JM. Olfactory dysfunction persists after smoking cessation and signals increased cardiovascular risk. International forum of allergy & rhinology 2019; 9(9): 977–85.
    
    

59. 59.Schubert CR, Cruickshanks KJ, Fischer ME, Klein BE, Klein R, Pinto AA. Inflammatory and vascular markers and olfactory impairment in older adults. Age and ageing 2015; 44(5): 878–82.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ageing/afv075&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26082178&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

60. 60.Wellford SA, Moseman EA. Olfactory immunology: the missing piece in airway and CNS defence. Nature Reviews Immunology 2024; 24(6): 381–98.
    
    

61. 61.Ashraf S, Frazier OH, Carranza S, McPherson DD, Taegtmeyer H. A Two-Step Transcriptome Analysis of the Human Heart Reveals Broad and Disease-Responsive Expression of Ectopic Olfactory Receptors. International Journal of Molecular Sciences 2023; 24(18).
    
    

62. 62.Wang B, Peng YJ, Su X, et al. Olfactory receptor 78 regulates erythropoietin and cardiorespiratory responses to hypobaric hypoxia. 2021; 130(4): 1122–32.
    
    

63. 63.Kim SH, Yoon YC, Lee AS, et al. Expression of human olfactory receptor 10J5 in heart aorta, coronary artery, and endothelial cells and its functional role in angiogenesis. Biochemical and biophysical research communications 2015; 460(2): 404–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbrc.2015.03.046&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25791473&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 

64. 64.Yang Z, Cheng J, Shang P, Sun JP, Yu X. Emerging roles of olfactory receptors in glucose metabolism. Trends in cell biology 2023; 33(6): 463–76.
    
    

65. 65.Lashutka MK, Chandra A, Murray HN, Phillips GS, Hiestand BC. The relationship of intraocular pressure to intracranial pressure. Annals of Emergency Medicine 2004; 43(5): 585–91.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.annemergmed.2003.12.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15111918&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221131600006&link_type=ISI) 

66. 66.Li Z, Yang Y, Lu Y, et al. Intraocular pressure vs intracranial pressure in disease conditions: a prospective cohort study (Beijing iCOP study). BMC Neurol 2012; 12: 66.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2377-12-66&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22862817&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F26%2F2024.08.26.24312564.atom)