Int J Med Sci 2026; 23(5):1738-1747. doi:10.7150/ijms.131188 This issue Cite

Research Paper

Potential influence of carbonic anhydrase 9 genetic variants and expression levels on the progression of diabetic retinopathy

I-Chia Liang1,2,3, Hsiang-Wen Chien3,4,5,6, Kai Wang3,4,5, Chia-Yi Lee7, Ying-Erh Chou8,9 Corresponding address, Shun-Fa Yang7,9 Corresponding address

1. Department of Ophthalmology, Tri-Service General Hospital, Taipei, Taiwan.
2. Department of Ophthalmology, College of Medicine, National Defense Medical Center, Taipei, Taiwan.
3. Department of Ophthalmology, Cathay General Hospital, Taipei, Taiwan.
4. Departments of Ophthalmology, Sijhih Cathay General Hospital, New Taipei City, Taiwan.
5. School of Medicine, College of Medicine, Fu Jen Catholic University, New Taipei, Taiwan.
6. School of Medicine, National Tsing Hua University, Hsinchu, Taiwan.
7. Institute of Medicine, Chung Shan Medical University, Taichung, Taiwan.
8. School of Medicine, Chung Shan Medical University, Taichung, Taiwan.
9. Department of Medical Research, Chung Shan Medical University Hospital, Taichung, Taiwan.

Received 2026-1-9; Accepted 2026-3-19; Published 2026-3-30

Citation:
Liang IC, Chien HW, Wang K, Lee CY, Chou YE, Yang SF. Potential influence of carbonic anhydrase 9 genetic variants and expression levels on the progression of diabetic retinopathy. Int J Med Sci 2026; 23(5):1738-1747. doi:10.7150/ijms.131188. https://www.medsci.org/v23p1738.htm
Other styles

File import instruction

Abstract

Graphic abstract

Diabetic retinopathy (DR) represents one of the most common microvascular complications of diabetes, and proliferative diabetic retinopathy (PDR) is the most vision-threatening form. Carbonic anhydrase IX (CA9), a hypoxia-inducible enzyme, has been implicated in several pathological processes, but its involvement in DR has not been clarified. In this study, three CA9 single nucleotide polymorphisms (rs3829078, rs2071676, and rs1048638) were genotyped in diabetic patients with and without DR. Clinical characteristics were compared between groups, and expression analyses were conducted using public databases and ARPE-19 cells under hyperglycemic and hypoxic conditions. No significant association was observed between CA9 variants and DR susceptibility in the overall cohort. However, among patients aged ≤60 years, carriers of the rs1048638 C/A and C/A + A/A genotypes exhibited a significantly increased risk of PDR. Expression quantitative trait locus (eQTL) data from the GTEx database revealed higher CA9 mRNA expression in tissues harboring the rs1048638 A allele. Moreover, both transcriptomic data and in vitro experiments demonstrated upregulation of CA9 in ARPE-19 cells exposed to high glucose and hypoxia. These findings suggest that the rs1048638 polymorphism may modulate CA9 expression and contribute to PDR development in younger diabetic patients through hypoxia- and glucose-related mechanisms. Collectively, our findings suggest that CA9 may serve as a potential risk biomarker associated with the progression of DR, particularly in younger patients.

Keywords: ARPE-19 cells, carbonic anhydrase IX, diabetic retinopathy, hypoxia, single nucleotide polymorphism, vascular dysfunction

1. Introduction

Diabetic retinopathy (DR) is one of the most common microvascular complications of diabetes mellitus. Non-proliferative diabetic retinopathy (NPDR) represents the early stage of the disease, whereas proliferative diabetic retinopathy (PDR) is an advanced form characterized by pathological neovascularization and is a leading cause of vision loss in adults [1]. The pathogenesis of DR is complex and multifactorial, involving a dynamic interplay among metabolic, hemodynamic, and genetic factors [2, 3]. The transition from NPDR to PDR is driven by progressive retinal ischemia, which induces the release of proangiogenic mediators, particularly vascular endothelial growth factor (VEGF), thereby promoting the formation of fragile, leaky neovessels [4]. Identifying patients at high risk for progression from NPDR to PDR is crucial, as those who develop PDR face a markedly increased likelihood of severe, vision-threatening complications and incur substantially higher medical costs [5]. Although clinical factors such as diabetes duration, poor glycemic control, and coexisting vascular complications (e.g., nephropathy or chronic non-healing ulcers) have been recognized as risk factors for DR progression, they do not fully account for the onset or severity of the disease [6, 7]. Genetic variability likely contributes to interindividual differences in DR susceptibility and progression, even among patients with similar clinical and metabolic profiles.

In PDR, hypoxia-inducible factor-1 alpha (HIF-1α), a transcription factor activated under hypoxic conditions, serves as a central regulator of genes mediating cellular adaptation to low oxygen tension [8, 9]. In response to hypoxia, HIF-1α directly binds to the promoter region of the VEGF gene and also upregulates the expression of carbonic anhydrase IX (CA9), a gene implicated in several cancer types [10-13]. Although hypoxia is the primary inducer of CA9 expression, proinflammatory cytokines such as interleukin-1 and -6 (IL-1, IL-6), tumor necrosis factor-alpha (TNF-α), and transforming growth factor-beta (TGF-β) have also been shown to modulate CA9 expression [14]. These cytokines likewise play crucial roles in the pathophysiology of PDR [15].

Carbonic anhydrases (CAs) constitute a large family of zinc metalloenzymes that catalyze the reversible hydration of carbon dioxide, generating bicarbonate and protons. The CA9 gene encodes carbonic anhydrase IX, a transmembrane protein that contributes to pH regulation and cell adhesion, processes critical for angiogenesis and the modulation of the tumor microenvironment [16]. Genetic variations in angiogenesis- and vascular-related genes have been shown to influence individual susceptibility to DR [17]. Single nucleotide polymorphisms (SNPs) within the CA9 gene have been associated with the development and progression of various cancers, potentially by altering gene expression, protein conformation, or molecular interactions [18-21]. While previous genetic studies have highlighted the roles of oxidative stress, apoptosis, inflammation, and autophagy in DR pathogenesis, the precise molecular mechanisms underlying disease initiation and progression remain unclear [22]. Given the established associations between CA9 polymorphisms and diverse clinicopathological features, it is plausible that these genetic variants may also modulate susceptibility to DR and influence its progression—an area that warrants further investigation.

2. Materials and Methods

2.1. Subjects

A total of 1,534 diabetic patients from the department of endocrinology at Chung Shan Medical University Hospital (Taichung, Taiwan) were recruited for this study. Prior to initiation, ethical approval was obtained from the Institutional Review Board of Chung Shan Medical University Hospital (IRB No. CS2-22190), and written informed consent was secured from all participants. Clinical data were retrieved from medical records, including the duration of diabetes, presence or absence of DR, glycated hemoglobin (HbA1c) levels, renal function, and lipid profiles. DR was defined as the presence of diabetes-related retinal microvascular abnormalities and was classified as NPDR or PDR. PDR was defined by the presence of retinal or optic disc neovascularization as previously described [23]. For the non-DR control group, the absence of DR was confirmed through dilated fundus examinations conducted by experienced ophthalmologists.

2.2. Specimen collection and genomic DNA extraction

Peripheral whole blood samples were collected in ethylenediaminetetraacetic acid (EDTA)-containing tubes, centrifuged immediately, and stored at -80°C [24]. Genomic DNA was extracted using the QIAamp DNA Blood Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer's protocol. Extracted DNA was dissolved in TE buffer (10 mM Tris, 1 mM EDTA; pH 7.8), and its concentration was determined by spectrophotometric measurement at 260 nm. The DNA samples were stored at -20 °C and used as templates for subsequent polymerase chain reaction (PCR) analyses [25].

2.3. Selection and Analysis of CA9 Genetic Polymorphisms

Among the more than 30 known single nucleotide polymorphisms (SNPs) in the CA9 gene, three SNPs, rs2071676 (+201, G/A) in exon 1, rs3829078 (+1081, A/G) in exon 7, and rs1048638 (+1584, C/A) in the 3′ untranslated region (UTR) of exon 11, were selected based on their potential functional relevance and previous reports indicating their involvement in hypoxia-related pathways and tumor vascular invasion [20, 26, 27]. Genotyping of these polymorphisms was performed using the ABI StepOne™ Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) and TaqMan SNP Genotyping Assays. Data were analyzed with SDS software version 3.0 (Applied Biosystems) [28].

2.4. Bioinformatics Analysis

To evaluate the association between the rs1048638 SNP and CA9 gene expression, data from the Genotype-Tissue Expression (GTEx) database were analyzed for correlations in aorta, coronary artery, and esophageal mucosal tissues [29]. Additionally, gene expression profiles were analyzed using datasets (GSE233164 and GSE151610) obtained from the Gene Expression Omnibus (GEO) database, as these datasets include gene expression data from samples relevant to diabetic retinopathy, enabling comparative analysis of gene expression patterns.

2.5. Cell culture

The ARPE-19 human retinal pigment epithelial cell line was obtained from the Bioresources Collection and Research Center (BCRC), Food Industry Research and Development Institute (Hsinchu, Taiwan). Cells were cultured in Dulbecco's Modified Eagle's Medium/Nutrient Mixture F-12 (DMEM/F12) supplemented with 10% fetal bovine serum (FBS) and 1% penicillin-streptomycin. Cells were maintained at 37 °C in a humidified atmosphere containing 5% CO₂. For experimental treatments, cells were exposed to high-glucose conditions (25 mM glucose) for 24 hours. Hypoxia was induced by incubation in a hypoxic chamber with 1% O₂ for 24 hours, depending on the experimental protocol.

2.6. Western blotting

ARPE-19 cells were exposed to either hypoxic or high-glucose conditions for 24 hours, washed with PBS, and lysed using ice-cold lysis buffer as previously described [30]. Total protein concentrations were determined, and 20 μg of protein from each sample was separated by SDS-PAGE and transferred onto a PVDF membrane. The membrane was blocked with 5% nonfat milk and incubated with anti-CA9 primary antibody at 4 °C overnight, followed by incubation with horseradish peroxidase (HRP)-conjugated secondary antibodies (Dako Corporation, Carpinteria, CA, USA) [31]. Protein bands were visualized using an enhanced chemiluminescence detection system, and band intensities were quantified using ImageJ software.

2.7. Real time PCR

Total RNA was extracted using the Total RNA Mini Kit (Geneaid, New Taipei City, Taiwan), and complementary DNA (cDNA) was synthesized using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA, USA). Real-time PCR amplification was conducted as described in previous studies, following standard protocols [32].

2.8. Statistical analysis

A p-value of < 0.05 was considered statistically significant. Differences in demographic and clinical characteristics between the NPDR and PDR groups were assessed using the Mann-Whitney U test and Fisher's exact test. Adjusted odds ratios (AORs) and 95% confidence intervals (CIs) for the associations between genotypes and DR risk were calculated using multiple logistic regression, controlling for potential confounders. Because subgroup analyses were performed according to age, the statistical power of these stratified analyses may be relatively limited due to reduced sample sizes and should therefore be interpreted with caution. All statistical analyses were performed using SAS software (version 9.1; SAS Institute, Cary, NC, USA) for Windows.

3. Results

3.1. Characteristics of the study participants

The clinical and laboratory characteristics of the two diabetic groups, patients without diabetic retinopathy (non-DR) and those with diabetic retinopathy (DR), are summarized in Table 1. No statistically significant differences were observed between the groups in terms of age, age at diabetes onset, or sex distribution. However, the duration of diabetes was significantly longer in the DR group compared to the non-DR group (9.05 vs. 11.30 years, p < 0.001). With respect to laboratory parameters, patients in the DR group exhibited significantly higher levels of glycated hemoglobin (HbA1c) and serum creatinine, as well as a lower estimated glomerular filtration rate (eGFR), compared to those in the non-DR group. No significant differences were found in lipid profile parameters, including total cholesterol, high-density lipoprotein (HDL), low-density lipoprotein (LDL), and triglycerides (TG).

3.2. Association of CA9 gene polymorphisms with diabetic retinopathy

The genotype distributions of three CA9 SNPs (rs3829078, rs2071676, and rs1048638) between the non-DR and DR groups are presented in Table 2. No significant differences were observed in the frequencies of these SNPs between the two groups. To further investigate the relationship between these polymorphisms and disease severity, the DR group was stratified into NPDR and PDR subgroups. As shown in Table 3, no statistically significant differences in the genotype distributions of the three CA9 SNPs were identified among the non-DR, NPDR, and PDR groups.

 Table 1 

Clinical and laboratory characteristics of diabetic patients with and without retinopathy.

VariableNon-Diabetic Retinopathy (N=875)Diabetic Retinopathy (N=659)p value
Age (years)60.88 ± 12.2961.99 ± 12.070.078
Age of onset (years)51.83 ± 11.5550.69 ± 11.500.057
Male gender [n (%)]473 (54.1%)350 (53.1%)0.713
Duration of diabetes (years)9.05 ± 7.2911.30 ± 8.13<0.001
HbA1c [% (mmol/mol)]7.07 ± 1.307.45 ± 1.42<0.001
Serum creatinine [mg/dL]0.97 ± 0.741.29 ± 1.38<0.001
Glomerular filtration rate [ml/min]82.04 ± 32.8874.11 ± 34.29<0.001
Total cholesterol [mmol/L]160.37 ± 45.32162.40 ± 43.450.381
HDL cholesterol [μmol/L]45.58 ± 12.1545.56 ± 14.210.979
LDL cholesterol [μmol/L]86.10 ± 30.5285.20 ± 31.680.577
Triglycerides, [μmol/L]144.30 ± 206.52144.07 ± 121.860.980
 Table 2 

Adjusted odds ratios (AORs) and 95% confidence intervals (CIs) for the association between CA9 genotypic frequencies and diabetic retinopathy.

VariableNo Diabetic Retinopathy (N=875)Diabetic Retinopathy (N=659)AOR (95% CI)p value
rs3829078
AA807 (92.2%)614 (93.2%)1.000
AG68 (7.8%)44 (6.7%)0.743 (0.491-1.123)0.158
GG0 (0%)1 (0.1%)------
AG+GG68 (7.8%)45 (6.8%)0.875 (0.712-1.074)0.200
rs2071676
GG216 (24.7%)160 (24.3%)1.000
GA464 (53.0%)339 (51.4%)0.963 (0.745-1.244)0.773
AA195 (22.3%)160 (24.3%)1.126 (0.833-1.522)0.441
GA+AA659 (75.3%)499 (75.7%)1.005 (0.890-1.135)0.930
rs1048638
CC764 (87.3%)577 (87.6%)1.000
CA107 (12.2%)78 (11.8%)0.978 (0.709-1.350)0.895
AA4 (0.5%)4 (0.6%)1.324 (0.303-5.792)0.709
CA+AA111 (12.7%)82 (12.4%)0.995 (0.850-1.166)0.952

The adjusted odds ratio (AOR) with their 95% confidence intervals were estimated by multiple logistic regression models after controlling for the duration of diabetes, HbA1c, serum creatinine levels and glomerular filtration rate.

 Table 3 

Adjusted odds ratios (AORs) and 95% confidence intervals (CIs) for the association between CA9 genotypic frequencies and non-proliferative or proliferative diabetic retinopathy.

No Diabetic Retinopathy (N=875)Non-Proliferative Diabetic Retinopathy (N=515)AOR (95% CI)Proliferative Diabetic Retinopathy (N=144)AOR (95% CI)
rs3829078
AA807 (92.2%)482 (93.6%)1.000132 (91.7%)1.000
AG68 (7.8%)32 (6.2%)0.728 (0.466-1.138)12 (8.3%)0.938 (0.464-1.898)
GG0 (0%)1 (0.2%)---0 (0.0%)---
AG+GG68 (7.8%)33 (6.4%)0.869 (0.697-1.083)12 (8.3%)0.938 (0.464-1.898)
rs2071676
GG216 (24.7%)126 (24.5%)1.00034 (23.6%)1.000
GA464 (53.0%)266 (51.7%)0.984 (0.750-1.291)73 (50.7%)0.901 (0.556-1.457)
AA195 (22.3%)123 (23.8%)1.109 (0.805-1.529)37 (25.7%)1.226 (0.708-2.120)
GA+AA659 (75.3%)389 (75.5%)1.010 (0.888-1.149)110 (76.4%)0.997 (0.795-1.251)
rs1048638
CC764 (87.3%)456 (88.5%)1.000121 (84.0%)1.000
CA107 (12.2%)55 (10.7%)0.885 (0.624-1.257)23 (16.0%)1.450 (0.847-2.484)
AA4 (0.5%)4 (0.8%)1.597 (0.371-6.866)0 (0.0%)---
CA+AA111 (12.7%)59 (11.5%)0.954 (0.804-1.132)23 (16.0%)1.187 (0.908-1.552)

The adjusted odds ratio (AOR) with their 95% confidence intervals were estimated by multiple logistic regression models after controlling for the duration of diabetes, HbA1c, serum creatinine levels and glomerular filtration rate.

 Table 4 

Adjusted odds ratios (AORs) and 95% confidence intervals (CIs) for the association between CA9 genotypic frequencies and diabetic retinopathy across different age groups.

Age ≤ 60 yearsAge > 60 years
No Diabetic Retinopathy
(N=384)
Diabetic Retinopathy (N=271)AOR (95% CI)No Diabetic Retinopathy
(N=491)
Diabetic Retinopathy (N=388)AOR (95% CI)
rs3829078
AA353 (91.9%)254 (93.7%)1.000454 (92.5%)360 (92.8%)1.000
AG31 (8.1%)17 (6.3%)0.659 (0.342-1.269)37 (7.5%)27 (7.0%)0.817 (0.478-1.395)
GG0 (0%)0 (0.0%)---0 (0%)1 (0.2%)---
AG+GG31 (8.1%)17 (6.3%)0.659 (0.342-1.269)37 (7.5%)28 (7.2%)0.925 (0.710-1.206)
rs2071676
GG94 (24.5%)68 (25.1%)1.000122 (24.8%)92 (23.7%)1.000
GA209 (54.4%)137 (50.6%)0.882 (0.594-1.309)255 (51.9%)202 (52.1%)1.036 (0.738-1.452)
AA81 (21.1%)66 (24.3%)1.172 (0.734-1.870)114 (23.3%)94 (24.2%)1.117 (0.751-1.661)
GA+AA290 (75.5%)203 (74.9%)0.981 (0.814-1.183)369 (75.2%)296 (76.3%)1.030 (0.877-1.209)
rs1048638
CC341 (88.8%)230 (84.9%)1.000423 (86.2%)347 (89.4%)1.000
CA42 (10.9%)40 (14.8%)1.292 (0.792-2.109)65 (13.2%)38 (9.8%)0.787 (0.509-1.216)
AA1 (0.3%)1 (0.3%)2.168 (0.132-35.572)3 (0.6%)3 (0.8%)1.014 (0.179-5.730)
CA+AA43 (11.2%)41 (15.1%)1.144 (0.898-1.458)68 (13.8%)41 (10.6%)0.893 (0.722-1.104)

The adjusted odds ratio (AOR) with their 95% confidence intervals were estimated by multiple logistic regression models after controlling for the duration of diabetes, HbA1c, serum creatinine levels and glomerular filtration rate.

Subsequently, the cohort was stratified by age into two groups: patients aged >60 years and those aged ≤60 years. As shown in Table 4, no significant associations were found between the CA9 polymorphisms and DR status in either age subgroup. However, when both age and disease severity were considered simultaneously, a significant association emerged. Among patients aged ≤60 years, carriers of the rs1048638 C/A genotype, as well as those with the combined C/A + A/A genotypes, showed a significantly higher risk of developing PDR compared with individuals carrying the C/C genotype. The adjusted odds ratio (AOR) for the C/A genotype was 1.509 (95% CI: 1.191-1.913, p = 0.035), and the 95% confidence interval does not include 1.0, supporting the statistical significance of this association. Similarly, individuals with the combined C/A + A/A genotypes had an increased risk of PDR, with an AOR of 1.475 (95% CI: 1.018-2.138, p = 0.040) (Table 5).

3.3. Functional and clinical relevance of CA9 rs1048638 in diabetic retinopathy

Given the observed association between rs1048638 and PDR, additional analyses using publicly available datasets were performed to explore the potential functional significance of this PDR-associated SNP. As shown in Figure 1A and 1B, expression quantitative trait locus (eQTL) analysis from the Genotype-Tissue Expression (GTEx) database demonstrated differential CA9 mRNA expression in aortic and esophageal tissues among individuals with distinct rs1048638 genotypes. To further elucidate the functional role of CA9 in ARPE-19 cells, we examined CA9 expression under high-glucose and hypoxic conditions using data obtained from the Gene Expression Omnibus (GEO) database. As shown in Figure 2A, data from GSE233164 revealed a significant upregulation of CA9 expression in ARPE-19 cells exposed to high glucose. Similarly, analysis of the GSE151610 dataset indicated that CA9 expression was markedly increased under hypoxic stress conditions (Figure 2B). The findings from the GEO datasets were further validated by our in vitro experiments. As shown in Figures 2C-2F, both CA9 protein and mRNA levels were significantly elevated in ARPE-19 cells following high-glucose treatment (Figure 2C-D) and hypoxic exposure (Figure 2E-F). To explore the downstream functional consequences of CA9 upregulation, we evaluated VEGF expression under hypoxic conditions and performed CA9 knockdown experiments. Analysis of the GSE151610 dataset indicated that VEGF expression was markedly increased under hypoxic stress (Figure 2G). Moreover, CA9-specific siRNA was used to silence CA9 in ARPE-19 cells under hypoxia, resulting in a significant reduction of VEGF mRNA levels (Figure 2H). These results suggest that CA9 contributes to the regulation of angiogenic signaling in retinal cells during hypoxic stress.

 Figure 1 

rs1048638 regulates CA9 expression. Expression quantitative trait locus (eQTL) analysis of rs1048638 in aortic and esophageal tissues was performed using data from the GTEx portal.

Int J Med Sci Image
 Table 5 

Adjusted odds ratios (AORs) and 95% confidence intervals (CIs) for the association between CA9 genotypic frequencies and non-proliferative or proliferative diabetic retinopathy in patients under 60 years of age.

No Diabetic Retinopathy (N=384)Non-Proliferative Diabetic Retinopathy (N=216)AOR (95% CI)Proliferative Diabetic Retinopathy (N=55)AOR (95% CI)
rs3829078
AA353 (91.9%)203 (94.0%)1.00051 (925.7%)1.000
AG31 (8.1%)13 (6.0%)0.669 (0.334-1.343)4 (7.3%)0.623 (0.181-2.142)
GG0 (0%)0 (0.0%)---0 (0.0%)---
AG+GG31 (8.1%)13 (6.0%)0.669 (0.334-1.343)4 (7.3%)0.623 (0.181-2.142)
rs2071676
GG94 (24.5%)57 (26.4%)1.00011 (20.0%)1.000
GA209 (54.4%)107 (49.5%)0.844 (0.558-1.277)30 (54.5%)1.112 (0.488-2.532)
AA81 (21.1%)52 (24.1%)1.119 (0.685-1.827)14 (25.5%)1.527 (0.598-3.901)
GA+AA290 (75.5%)159 (73.6%)0.959 (0.789-1.166)44 (80.0%)1.106 (0.747-1.638)
rs1048638
CC341 (88.8%)187 (86.6%)1.00043 (78.2%)1.000
CA42 (10.9%)28 (13.0%)1.216 (0.720-2.054)12 (21.8%)1.509 (1.191-1.913)a
AA1 (0.3%)1 (0.4%)2.579 (0.157-42.445)0 (0.0%)---
CA+AA43 (11.2%)29 (13.4%)1.115 (0.861-1.444)12 (21.8%)1.475 (1.018-2.138)b

The adjusted odds ratio (AOR) with their 95% confidence intervals were estimated by multiple logistic regression models after controlling for the duration of diabetes, HbA1c, serum creatinine levels and glomerular filtration rate. ap=0.035 bp=0.040

 Figure 2 

Validation of CA9 and VEGF gene expression in ARPE-19 cells using datasets from the Gene Expression Omnibus (GEO) and in vitro assay. (A) CA9 expression was significantly upregulated in ARPE-19 cells treated with high glucose (dataset: GSE233164). (B) Comparison of CA9 expression in ARPE-19 cells under normoxic and hypoxic conditions (dataset: GSE151610). (C-D) Effects of high-glucose treatment on CA9 protein and mRNA levels in ARPE-19 cells were analyzed by Western blotting and real-time PCR. (E-F) Effects of hypoxic stress (1% O₂) on CA9 protein and mRNA levels in ARPE-19 cells were assessed by Western blotting and real-time PCR. (G) VEGF expression under hypoxic stress in ARPE-19 cells (dataset: GSE151610). (H) ARPE-19 cells were transfected with CA9-specific siRNA and exposed to hypoxia (1% O₂). CA9 knockdown led to a significant reduction in VEGF mRNA levels compared with hypoxia-exposed cells. Data are presented as mean ± SD. *p < 0.001 versus control. # p < 0.001 versus hypoxia condition.

Int J Med Sci Image

4. Discussion

The development of DR is influenced by a complex interplay of genetic predisposition and acquired lifestyle and metabolic factors [33, 34]. In the present study, we observed that younger diabetic patients (aged ≤60 years) carrying the rs1048638 C/A or C/A + A/A genotypes of the CA9 gene had a significantly increased risk of developing PDR. Data from the GTEx database indicated elevated CA9 mRNA expression in tissues harboring the A allele of rs1048638, suggesting a regulatory role of this polymorphism in gene expression. Furthermore, transcriptomic analyses from the GEO database and in vitro experiments in ARPE-19 cells demonstrated upregulation of CA9 expression under high-glucose and hypoxic conditions. These findings collectively support a potential role for CA9 polymorphism in the pathogenesis of DR, particularly PDR, in younger diabetic individuals, potentially mediated through glucose- and hypoxia-responsive pathways.

The initiation and progression of DR are believed to be driven by retinal capillary dysfunction and subsequent tissue hypoxia, which stimulate angiogenesis and neovascularization [35, 36]. Carbonic anhydrases (CAs), including CA9, have been increasingly implicated in the pathophysiology of diabetic complications, and inhibitors of CA activity have shown therapeutic potential [37, 38]. The L-arginine/nitric oxide (NO) signaling axis, which is often disrupted in diabetes, plays a key role in maintaining vascular homeostasis and is tightly linked to hypoxic stress in the retina [39, 40]. Reduced NO bioavailability contributes to impaired vasodilation, decreased retinal blood flow, and aggravated hypoxia. Moreover, an imbalance between NO and reactive oxygen species (ROS) exacerbates oxidative stress, damaging retinal neurons and microvasculature. Dysregulation of this pathway also influences inflammatory signaling, further contributing to retinal injury [39, 41, 42].

Recent studies have identified the activation of the arginase pathway via nicotinamide adenine dinucleotide phosphate (NADPH) oxidase (NOX) as a major contributor to retinal vascular dysfunction in DR. NOX activation is known to modulate CA9 expression and, conversely, CA9 may influence NOX activity [39]. Notably, CA9 has been implicated in the regulation of the L-arginine/NO pathway and has been found to be altered in patients with type 1 diabetes [43]. Higher CA9 protein levels have been associated with more rapid DR progression in individuals with type 1 diabetes, compared to those with slower disease trajectories [40]. Experimental studies further demonstrate that CA inhibitors improve oxygenation in the optic nerve and retina [44, 45], and topical CA inhibitors have been shown to enhance blood flow velocity in major retinal arteries [36]. However, our experiments focused on ARPE-19 cells, DR pathogenesis also involves vascular endothelial cells, which play a central role in microvascular damage and neovascularization. Future studies should investigate the functional consequences of CA9 variants in retinal endothelial cells and other relevant cell types to clarify their contribution to DR progression.

Age is an important modifier in the development of diabetic complications. It reflects a multifactorial process shaped by genetic, metabolic, and environmental factors. Younger individuals with diabetes may have increased susceptibility to hypoxic injury due to higher metabolic demand, early-stage diabetic nephropathy, or impaired systemic oxygen delivery [46]. In this study, the age cutoff of 60 years, which has been commonly used in epidemiological studies of diabetic complications to distinguish younger and older patient populations, was applied to explore potential age-dependent genetic effects. In this context, genetic variants such as CA9 polymorphisms, previously implicated in tumorigenesis and hypoxia adaptation, may influence physiological responses to hypoxia and modulate disease risk across age groups [47, 48]. In our study, the rs1048638 C/A and A/A genotypes were associated with an increased risk of PDR among patients aged ≤60 years. However, given the moderate effect size (AOR ≈ 1.5) and the exploratory nature of the stratified analyses, these findings should be interpreted with caution. Multiple testing considerations were taken into account when evaluating the statistical significance of the stratified results, and further validation in larger and independent cohorts will be necessary to confirm the clinical relevance of this association.

In our study, the rs1048638 C/A and A/A genotypes were associated with a significantly increased risk of PDR in patients aged ≤60 years. This polymorphism may influence CA9 mRNA stability or translational efficiency, thereby affecting protein expression and cellular responses to hypoxic stress. In cervical cancer, the A allele of rs1048638 has been reported to disrupt the binding site for microRNA-34a (miR-34a), leading to elevated CA9 expression and potentially promoting tumor progression [49, 50]. Similar associations have been observed in urothelial carcinoma, where the CA/AA genotype correlates with more aggressive tumor behavior [49], as well as in colorectal and hepatocellular carcinomas [19, 27]. These effects are likely mediated through impaired miR-34a-mediated regulation. MiR-34a plays a multifaceted role in diabetes, where it regulates vascular integrity, pancreatic β-cell function, glucose metabolism, and other key physiological processes [51]. In retinal endothelial cells, miR-34a suppresses mitochondrial function and biogenesis, enhancing susceptibility to oxidative stress [52]. Under hyperglycemic conditions, miR-34a dysregulation contributes to vascular dysfunction and promotes stress-induced premature senescence, thereby accelerating retinal microvascular damage and cellular aging [53]. However, the present study did not investigate the interaction between rs1048638 and miR-34a in ARPE-19 cells or other retinal cell types. Therefore, the precise molecular mechanism by which this polymorphism may influence CA9 expression in the context of diabetic retinopathy remains unclear. Future studies are warranted to explore whether the rs1048638 variant affects miR-34a binding and CA9 regulation in retinal cells under hyperglycemic or hypoxic conditions.

Taken together, our findings suggest a potential modulatory effect of CA9 polymorphisms on the progression of DR in younger patients. However, several limitations should be acknowledged. First, the study was conducted in a single ethnic population, which may limit the generalizability of the findings. Second, gene expression analyses relied in part on the GTEx database, which comprises non-diabetic, non-retinal tissues; thus, the inferred expression patterns may not fully reflect the diabetic retina. Finally, the present study did not investigate the molecular mechanisms by which CA9 variants, such as rs1048638, may influence gene expression via microRNA interactions, including miR-34a. Future studies with larger, multiethnic cohorts, longitudinal designs, and mechanistic analyses in retinal cells are warranted to validate these findings and elucidate the functional role of CA9 polymorphisms in DR pathogenesis.

In conclusion, this study provides evidence that the CA9 rs1048638 polymorphism is associated with an increased risk of proliferative diabetic retinopathy in younger patients with diabetes. Individuals aged ≤60 years carrying the C/A or A/A genotypes demonstrated a significantly higher likelihood of developing PDR. Public transcriptomic datasets and in vitro experiments confirmed that CA9 expression is upregulated under hyperglycemic and hypoxic conditions, supporting its role in DR pathophysiology. The involvement of CA9 in hypoxia adaptation, the L-arginine/NO signaling pathway, and miR-34a regulation further underscores its potential contribution to retinal vascular dysfunction. These findings suggest a possible gene-environment interaction that warrants further exploration through well-powered, longitudinal studies across diverse populations. Ultimately, CA9 may serve as a candidate biomarker or therapeutic target for DR, particularly in younger at-risk individuals.

Competing Interests

The authors have declared that no competing interest exists.

References

1. Fong DS, Aiello L, Gardner TW, King GL, Blankenship G, Cavallerano JD. et al. Retinopathy in diabetes. Diabetes Care. 2004;27(Suppl 1):S84-7

2. Eshaq RS, Aldalati AMZ, Alexander JS, Harris NR. Diabetic retinopathy: Breaking the barrier. Pathophysiology. 2017;24:229-41

3. Mu L, Wang D, Dong Z, Wu J, Wu X, Su J. et al. Abnormal Levels of Serum Ferroptosis-Related Biomarkers in Diabetic Retinopathy. J Ophthalmol. 2022;2022:3353740

4. Aiello LP, Wong JS. Role of vascular endothelial growth factor in diabetic vascular complications. Kidney Int Suppl. 2000;77:S113-9

5. Schmier JK, Covert DW, Lau EC, Matthews GP. Medicare expenditures associated with diabetes and diabetic retinopathy. Retina. 2009;29:199-206

6. Harris Nwanyanwu K, Talwar N, Gardner TW, Wrobel JS, Herman WH, Stein JD. Predicting development of proliferative diabetic retinopathy. Diabetes Care. 2013;36:1562-8

7. Yang J, Fan XH, Guan YQ, Li Y, Sun W, Yang XZ. et al. MMP-2 gene polymorphisms in type 2 diabetes mellitus diabetic retinopathy. Int J Ophthalmol. 2010;3:137-40

8. Vadlapatla RK, Vadlapudi AD, Mitra AK. Hypoxia-inducible factor-1 (HIF-1): a potential target for intervention in ocular neovascular diseases. Curr Drug Targets. 2013;14:919-35

9. Lee D, Tomita Y, Miwa Y, Kunimi H, Nakai A, Shoda C. et al. Recent Insights into Roles of Hypoxia-Inducible Factors in Retinal Diseases. Int J Mol Sci. 2024 25

10. Olive PL, Aquino-Parsons C, MacPhail SH, Liao SY, Raleigh JA, Lerman MI. et al. Carbonic anhydrase 9 as an endogenous marker for hypoxic cells in cervical cancer. Cancer Res. 2001;61:8924-9

11. Kim SJ, Rabbani ZN, Dewhirst MW, Vujaskovic Z, Vollmer RT, Schreiber EG. et al. Expression of HIF-1alpha, CA IX, VEGF, and MMP-9 in surgically resected non-small cell lung cancer. Lung Cancer. 2005;49:325-35

12. Lin CW, Yang WE, Lee WJ, Hua KT, Hsieh FK, Hsiao M. et al. Lipocalin 2 prevents oral cancer metastasis through carbonic anhydrase IX inhibition and is associated with favourable prognosis. Carcinogenesis. 2016;37:712-22

13. Yang JS, Chou CH, Hsieh YH, Tang CH, Lu KH, Yang SF. Carbonic Anhydrase IX Induces Human Osteosarcoma Cell Metastasis by Activating HSPA6 Expression Through the AMPK Signalling Pathway. J Cell Mol Med. 2025;29:e70763

14. Kockar F, Yildrim H, Sagkan RI, Hagemann C, Soysal Y, Anacker J. et al. Hypoxia and cytokines regulate carbonic anhydrase 9 expression in hepatocellular carcinoma cells in vitro. World J Clin Oncol. 2012;3:82-91

15. El-Asrar AM. Role of inflammation in the pathogenesis of diabetic retinopathy. Middle East Afr J Ophthalmol. 2012;19:70-4

16. Lin CY, Wang SS, Yang CK, Li JR, Chen CS, Hung SC. et al. Genetic polymorphism and carbonic anhydrase 9 expression can predict nodal metastatic prostate cancer risk in patients with prostate-specific antigen levels </=10 ng/ml at initial biopsy. Urol Oncol. 2019;37:814 e9- e16

17. Rogers MS, D'Amato RJ. The effect of genetic diversity on angiogenesis. Exp Cell Res. 2006;312:561-74

18. Pastorekova S, Gillies RJ. The role of carbonic anhydrase IX in cancer development: links to hypoxia, acidosis, and beyond. Cancer Metastasis Rev. 2019;38:65-77

19. Huang HC, Shiu BH, Nassef Y, Huang CC, Chou YE, Ting WC. et al. Impact of carbonic anhydrase 9 gene polymorphism on the progression of colorectal cancer. J Cancer. 2022;13:2775-80

20. Wang SS, Liu YF, Ou YC, Chen CS, Li JR, Yang SF. Impacts of CA9 gene polymorphisms on urothelial cell carcinoma susceptibility and clinicopathologic characteristics in Taiwan. PLoS One. 2013;8:e82804

21. Lin CY, Wang SS, Yang CK, Li JR, Chen CS, Hung SC. et al. Genetic polymorphism and carbonic anhydrase 9 expression can predict nodal metastatic prostate cancer risk in patients with prostate-specific antigen levels ≤10 ng/ml at initial biopsy. Urol Oncol. 2019;37:814.e9-e16

22. Hu L, Liu Y, Wei C, Jin H, Mei L, Wu C. SERPINH1, Targeted by miR-29b, Modulated Proliferation and Migration of Human Retinal Endothelial Cells Under High Glucose Conditions. Diabetes Metab Syndr Obes. 2021;14:3471-83

23. Yao YP, Chien HW, Wang K, Yang YS, Su SC, Chang LC. et al. Genetic association of diabetic retinopathy with long noncoding RNA CDKN2B-AS1 gene polymorphism. Eur J Ophthalmol. 2025;35:637-44

24. Weng WC, Hsieh YH, Lin CY, Liu YF, Su SC, Wang SS. et al. Functional variants of the pentraxin 3 gene are associated with the metastasis and progression of prostate cancer. J Cell Mol Med. 2024;28:e70041

25. Yeh JC, Chen YT, Chou YE, Su SC, Chang LC, Chen YL. et al. Interactive effects of CDKN2B-AS1 gene polymorphism and habitual risk factors on oral cancer. J Cell Mol Med. 2023;27:3395-403

26. Fraga A, Ribeiro R, Coelho A, Vizcaino JR, Coutinho H, Lopes JM. et al. Genetic polymorphisms in key hypoxia-regulated downstream molecules and phenotypic correlation in prostate cancer. BMC Urol. 2017;17:12

27. Hua KT, Liu YF, Hsu CL, Cheng TY, Yang CY, Chang JS. et al. 3'UTR polymorphisms of carbonic anhydrase IX determine the miR-34a targeting efficiency and prognosis of hepatocellular carcinoma. Sci Rep. 2017;7:4466

28. Chen YT, Lin CW, Chou YE, Su SC, Chang LC, Lee CY. et al. Potential impact of ADAM-10 genetic variants with the clinical features of oral squamous cell carcinoma. J Cell Mol Med. 2023;27:1144-52

29. Wang Z, Jensen MA, Zenklusen JC. A Practical Guide to The Cancer Genome Atlas (TCGA). Methods Mol Biol. 2016;1418:111-41

30. Yang SF, Chen YS, Chien HW, Wang K, Lin CL, Chiou HL. et al. Melatonin attenuates epidermal growth factor-induced cathepsin S expression in ARPE-19 cells: Implications for proliferative vitreoretinopathy. J Pineal Res. 2020;68:e12615

31. Su SC, Lin CW, Chen MK, Lee YC, Su CW, Bai S. et al. Multimodal profiling of oral squamous cell carcinoma identifies genomic alterations and expression programs associated with betel quid chewing. Neoplasia. 2025;68:101218

32. Hu DN, Zhang R, Yao S, Iacob CE, Yang WE, Rosen R. et al. Cultured Human Uveal Melanocytes Express/secrete CXCL1 and CXCL2 Constitutively and Increased by Lipopolysaccharide via Activation of Toll-like Receptor 4. Curr Eye Res. 2021;46:1681-94

33. Baek SU, Park MS, Cho BJ, Park IW, Kwon S. Risk factors associated with progression of diabetic retinopathy in eyes treated with panretinal photocoagulation. Sci Rep. 2021;11:13850

34. Beaser RS, Turell WA, Howson A. Strategies to Improve Prevention and Management in Diabetic Retinopathy: Qualitative Insights from a Mixed-Methods Study. Diabetes Spectr. 2018;31:65-74

35. Arden GB, Wolf JE, Tsang Y. Does dark adaptation exacerbate diabetic retinopathy? Evidence and a linking hypothesis. Vision Res. 1998;38:1723-9

36. Siesky B, Harris A, Brizendine E, Marques C, Loh J, Mackey J. et al. Literature review and meta-analysis of topical carbonic anhydrase inhibitors and ocular blood flow. Surv Ophthalmol. 2009;54:33-46

37. Sulaiman AH, Ghassan ZI, Omar TN. Biochemical Evaluation of Carbonic Anhydrase and Some Antioxidant Markers in Patients with Diabetes Complications. Arch Razi Inst. 2022;77:169-78

38. Matough FA, Budin SB, Hamid ZA, Alwahaibi N, Mohamed J. The role of oxidative stress and antioxidants in diabetic complications. Sultan Qaboos Univ Med J. 2012;12:5-18

39. Caldwell RB, Zhang W, Romero MJ, Caldwell RW. Vascular dysfunction in retinopathy-an emerging role for arginase. Brain Res Bull. 2010;81:303-9

40. Lundkvist P, Gronberg A, Carlsson PO, Ludvigsson J, Espes D. Predictive biomarkers of rapidly developing insulin deficiency in children with type 1 diabetes. BMJ Open Diabetes Res Care. 2024 12

41. Patel C, Rojas M, Narayanan SP, Zhang W, Xu Z, Lemtalsi T. et al. Arginase as a mediator of diabetic retinopathy. Front Immunol. 2013;4:173

42. Goldstein IM, Ostwald P, Roth S. Nitric oxide: a review of its role in retinal function and disease. Vision Res. 1996;36:2979-94

43. Carmann C, Lilienthal E, Weigt-Usinger K, Schmidt-Choudhury A, Horster I, Kayacelebi AA. et al. The L-arginine/NO pathway, homoarginine, and nitrite-dependent renal carbonic anhydrase activity in young people with type 1 diabetes mellitus. Amino Acids. 2015;47:1865-74

44. Stefansson E, Jensen PK, Eysteinsson T, Bang K, Kiilgaard JF, Dollerup J. et al. Optic nerve oxygen tension in pigs and the effect of carbonic anhydrase inhibitors. Invest Ophthalmol Vis Sci. 1999;40:2756-61

45. Noergaard MH, Bach-Holm D, Scherfig E, Bang K, Jensen PK, Kiilgaard JF. et al. Dorzolamide increases retinal oxygen tension after branch retinal vein occlusion. Invest Ophthalmol Vis Sci. 2008;49:1136-41

46. Catrina SB, Zheng X. Hypoxia and hypoxia-inducible factors in diabetes and its complications. Diabetologia. 2021;64:709-16

47. Naval J, Alonso V, Herranz MA. Genetic polymorphisms and skin aging: the identification of population genotypic groups holds potential for personalized treatments. Clin Cosmet Investig Dermatol. 2014;7:207-14

48. Lamparello AJ, Namas RA, Schimunek L, Cohen M, El-Dehaibi F, Yin J. et al. An Aging-Related Single-Nucleotide Polymorphism is Associated With Altered Clinical Outcomes and Distinct Inflammatory Profiles in Aged Blunt Trauma Patients. Shock. 2020;53:146-55

49. Yang SF, Liu YF, Cheng CW, Yang WE, Lin WL, Ko JL. et al. Impact of microRNA-34a and polymorphism of its target gene CA9 on susceptibility to uterine cervical cancer. Oncotarget. 2017;8:77860-71

50. Shen HP, Hsiao YH, Yang SF, Liu YF, Ko JL, Wang PH. Single nucleotide polymorphisms and haplotypes of carbonic anhydrase 9 can predict invasive squamous cell carcinoma of uterine cervix. Int J Med Sci. 2018;15:587-94

51. Fomison-Nurse I, Saw EEL, Gandhi S, Munasinghe PE, Van Hout I, Williams MJA. et al. Diabetes induces the activation of pro-ageing miR-34a in the heart, but has differential effects on cardiomyocytes and cardiac progenitor cells. Cell Death Differ. 2018;25:1336-49

52. Thounaojam MC, Bartoli M. MicroRNA-34a and vascular senescence in diabetes. Aging (Albany NY). 2019;11:11799-800

53. Thounaojam MC, Jadeja RN, Warren M, Powell FL, Raju R, Gutsaeva D. et al. MicroRNA-34a (miR-34a) Mediates Retinal Endothelial Cell Premature Senescence through Mitochondrial Dysfunction and Loss of Antioxidant Activities. Antioxidants (Basel). 2019 8

Author contact

Corresponding address Corresponding authors: Shun-Fa Yang, PhD. Or Ying-Erh Chou, PhD. Institute of Medicine, Chung Shan Medical University, 110, Section 1, Chien-Kuo N. Road, Taichung, Taiwan, ROC. Fax: 886-4-24723229. E-mail: ysfedu.tw (Shun-Fa Yang); E-mail: intointo814com (Ying-Erh Chou).


Citation styles

APA
Liang, I.C., Chien, H.W., Wang, K., Lee, C.Y., Chou, Y.E., Yang, S.F. (2026). Potential influence of carbonic anhydrase 9 genetic variants and expression levels on the progression of diabetic retinopathy. International Journal of Medical Sciences, 23(5), 1738-1747. https://doi.org/10.7150/ijms.131188.

ACS
Liang, I.C.; Chien, H.W.; Wang, K.; Lee, C.Y.; Chou, Y.E.; Yang, S.F. Potential influence of carbonic anhydrase 9 genetic variants and expression levels on the progression of diabetic retinopathy. Int. J. Med. Sci. 2026, 23 (5), 1738-1747. DOI: 10.7150/ijms.131188.

NLM
Liang IC, Chien HW, Wang K, Lee CY, Chou YE, Yang SF. Potential influence of carbonic anhydrase 9 genetic variants and expression levels on the progression of diabetic retinopathy. Int J Med Sci 2026; 23(5):1738-1747. doi:10.7150/ijms.131188. https://www.medsci.org/v23p1738.htm

CSE
Liang IC, Chien HW, Wang K, Lee CY, Chou YE, Yang SF. 2026. Potential influence of carbonic anhydrase 9 genetic variants and expression levels on the progression of diabetic retinopathy. Int J Med Sci. 23(5):1738-1747.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See https://ivyspring.com/terms for full terms and conditions.
Popup Image