A Multicentric Study of ACE2 rs2285666 Polymorphism: The Association with COVID-19 in the Iranian Population

authors:

avatar Hadi Razavi Nikoo ORCID 1 , 2 , avatar Seyed Majid Fatahi 3 , avatar Mahdieh Taheri 2 , avatar Mohammad Jamshidi 4 , avatar Ali Safarzadeh 5 , avatar Ashkan Alamdary 6 , avatar Hossein Teimouri 2 , avatar Bahman Aghcheli 2 , avatar Alijan Tabarraei 2 , avatar Sadegh Ali Azimi 2 , avatar Ali Asghar Kiani 4 , avatar Mehdi Ajorloo 4 , 6 , *

Laboratory Sciences Research Center, Golestan University of Medical Sciences, Gorgan, Iran
Department of Microbiology, Faculty of Medicine, Golestan University of Medical Sciences, Gorgan, Iran
Student Research Committee, Lorestan University of Medical Sciences, Khorramabad, Iran
Molecular and Cellular Research Laboratory, Department of Clinical Laboratory Sciences, Lorestan University of Medical Sciences, Khorramabad, Iran
Department of Biology, University of Padova, Padova, Italy
Blood Transfusion Research Center, High Institute for Research and Education in Transfusion Medicine, Tehran, Iran

how to cite: Razavi Nikoo H, Fatahi S M, Taheri M, Jamshidi M, Safarzadeh A, et al. A Multicentric Study of ACE2 rs2285666 Polymorphism: The Association with COVID-19 in the Iranian Population. Jundishapur J Microbiol. 2024;17(10):e153765. https://doi.org/10.5812/jjm-153765.

Abstract

Background:

Clinical studies have shown that variations in the ACE2 gene can influence susceptibility to and the severity of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infections. Although the evidence is not conclusive, recent studies have suggested a potential link between the ACE2 rs2285666 polymorphism and COVID-19.

Objectives:

This multicentric clinical study, conducted in Lorestan and Golestan provinces, aimed to examine the association between the ACE2 rs2285666 polymorphism and COVID-19 in Iranians of various ethnicities, including Fars, Lur, Turkmen, and Balooch.

Methods:

A total of 372 participants were evaluated, with 201 testing positive for SARS-CoV-2 and 171 negative. The PCR-based restriction fragment length polymorphism (PCR-RFLP) technique, using the AluI enzyme, was employed to identify the polymorphism. Demographic and clinical data were analyzed using IBM SPSS. Hazards were assessed through odds ratio analysis, while Hardy-Weinberg equilibrium (HWE) and genotype variations were evaluated using SNPSTATS.

Results:

PCR-RFLP analysis indicated that the AA genotype may increase susceptibility to COVID-19. Among female COVID-19-positive patients, 56%, 29%, and 15% exhibited GG, GA, and AA genotypes, respectively, compared to 61%, 35%, and 4% in the control group. In male individuals, the frequencies of G and A genotypes were 89% and 11% in the healthy group, while they were 25% and 75%, respectively, in the patient group. Polymorphism frequencies were not in HWE in both the positive and negative groups (P < 0.05). Logistic regression analysis showed that the AA genotype differed in co-dominant and recessive inheritance models, with odds ratios (OR) of OR = 4.06 (1.10 - 15.00) and OR = 4.21 (1.16 - 15.24). The ACE2 rs2285666 AA or A genotype was strongly associated with an increased risk of COVID-19 in this study.

Conclusions:

A significant difference in the distribution of the AA and A genotypes was observed in COVID-19 patients. Further studies involving larger and more diverse populations are necessary to explore the impact of this polymorphism on the susceptibility to and severity of COVID-19.

1. Background

Rapid person-to-person transmission and the increasing number of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) cases have significantly impacted the healthcare system, with millions of confirmed cases and deaths worldwide (1). Severe acute respiratory syndrome coronavirus-2 belongs to the genus Beta coronavirus within the family Coronaviridae. It is a helical, enveloped virus composed of an unsegmented, single-stranded, positive-sense RNA genome of approximately 30 kb, which is capped at the 5'-end and polyadenylated at the 3'-end. Severe acute respiratory syndrome coronavirus-2 primarily transmits through respiratory routes, and the disease it causes is known as Coronavirus disease 2019 (COVID-19) (2). The clinical presentation of COVID-19 can vary widely, ranging from asymptomatic infection to mild disease (81%), or progressing to severe disease (14%) that may result in multi-organ failure and death (5%) (3). Despite extensive research on SARS-CoV-2, the virus's pathophysiology remains incompletely understood and requires further investigation. The pathophysiology of acute respiratory distress syndrome (ARDS) is a multifactorial process, with recent evidence suggesting that host factors, including genetic variations, along with viral and environmental factors, may be involved in the disease’s pathogenesis and outcomes (4-6).

The renin-angiotensin system (RAS) is a hormone system composed of various angiotensin peptides and specific receptors. It is a critical regulator of the cardiovascular, renal, and pulmonary systems and plays a key role in the pathogenesis of many diseases. Severe acute respiratory syndrome coronavirus-2 targets human respiratory cells, including nasal and bronchial epithelial cells and pneumocytes, by binding its viral spike glycoprotein to the ACE2 receptor (7). ACE2 is a 120 kDa single-span transmembrane protein and a zinc-containing metalloenzyme that removes a carboxy-terminal residue from peptide hormones, such as angiotensin II, to produce angiotensin (1-7).

It functions as a negative regulator of the RAS (7-9). From the N-terminal to the C-terminal, ACE2 has four functional and structural domains: The peptidase domain (PD, aa 18 - 615), the collectrin-like domain (CLD, aa 616 - 740), the transmembrane domain (TM, aa 741 - 761), and the intracellular domain (ICD, aa 762 - 805). Additionally, the N-terminal contains a signal peptide (SigP) of 17 amino acids (7-9). Variations in the ACE2 gene can impact the renin-angiotensin pathway by increasing the level of ACE2 in the blood or at specific tissue sites. This can potentially damage the vascular endothelium and lung epithelium, thus increasing susceptibility to viral infections, such as SARS-CoV-2 (10-12).

Several in silico studies have described a correlation between ACE2 polymorphisms and COVID-19 outcomes (13-15). The fourth base of intron 3 in the ACE2 gene contains the rs2285666 polymorphism, also known as G8790A. This polymorphism can alter the process of messenger RNA (mRNA) alternative splicing, which, in turn, affects the level of ACE2 gene expression (13). Recently, researchers have revealed that the ACE2 rs2285666 polymorphism does not influence susceptibility to SARS-CoV-2 infection or the clinical outcome of the disease (16). However, some studies have found an association between ACE2 rs2285666 and COVID-19, suggesting that the GG genotype or the G-allele significantly increases the risk of infection and the severity of COVID-19 (13, 17-19). In contrast, other studies have found no association of this SNP with the disease or its severity (20). The variation in these results may be attributed to differences in ethnicity and the sample populations used in these studies. Therefore, analyzing this polymorphism and its different alleles in a diverse population of various ethnicities in Iran could provide more reliable conclusions.

2. Objectives

The current study aimed to investigate the ACE2 rs2285666 polymorphism in a diverse population of Iranian COVID-19 patients to determine its association with the disease. Additionally, we explored the relationship between the polymorphism and clinical-laboratory data from COVID-19 patients.

3. Methods

3.1. Population of the Study

The clinical samples consisted of 372 subjects who agreed to participate in this multicentric clinical study in Lorestan and Golestan provinces from May to October 2019. We enrolled 201 individuals who developed COVID-19 as contacts. Additionally, a family physician evaluated 172 Iranian medical staff members who showed no clinical signs of COVID-19, and an RT-qPCR test was conducted for confirmation. COVID-19 confirmation was based on the detection of the SARS-CoV-2 genome using RT-qPCR, following World Health Organization (WHO) guidelines. Cases with underlying diseases were excluded from the study population. A total of 5 mL of venous blood was collected from all participants in EDTA-containing tubes for further analysis. Hematological parameters, including erythrocyte sedimentation rate (ESR), hemoglobin (Hb), hematocrit (HCT), platelet (PLT), white blood cell (WBC), red blood cell (RBC), neutrophils, and lymphocytes, were assessed according to the manufacturer’s instructions. In this study, although there was no bias in sample selection, some confounding factors, such as co-morbidities, socioeconomic status, and treatments received, were not considered, which could be limitations of the study.

3.2. DNA Extraction and Genotyping

Genomic DNA was extracted from 350 μL of blood using a DNA isolation kit (Gene All, South Korea). To investigate the presence of the rs2285666 SNP in both groups, a PCR-based restriction fragment length polymorphism (PCR-RFLP) assay was performed. Briefly, 20 μL reaction mixtures were prepared containing 10 μL of Amplicon-PCR Master Mix (10× PCR buffer, 3.0 mM MgCl₂, 0.25 mM dNTPs, 1.5 units of Taq polymerase), 0.5 µmol/L of each primer (F: 5′-CATGTGGTCAAAAGGATATCT-3′, R: 5′-AAAGTAAGGTTGGCAGACAT-3′), 3 µL of purified genomic DNA, and 7 μL of distilled water. Thermal cycling conditions were as follows: One cycle at 95°C for 5 min, followed by 40 cycles consisting of three steps: 95°C for 15 seconds, 60°C for 30 seconds, and 72°C for 30 seconds, with a final extension at 72°C for 10 minutes. PCR product digestion was performed by adding 10 U/µL AluI restriction endonuclease (Fermentas, Latvia) (1.5 U), 2 μL of Tango buffer (10×), and 5 μL of PCR products. Distilled water was added to bring the total volume to 20 μL, and the mixture was incubated for 16 hours at 37°C.

3.3. Statistical Analysis

The demographic and clinical data were analyzed using SPSS version 26 software (IBM, Chicago, IL). Hardy-Weinberg equilibrium (HWE) and differences in genotypes were assessed using SNPSTATS online software. Odds ratios (OR) and 95% confidence intervals (CI) were calculated. The chi-square test was applied to determine intergroup differences, while the t-test and Mann-Whitney test were used to compare the means of the two studied groups. A P-value ≤ 0.05 was considered statistically significant for all tests and data analyses.

4. Results

4.1. Genotyping

The PCR-RFLP method was used to detect the ACE2 rs2285666 polymorphism. The PCR products were digested using the AluI enzyme. After digestion, the PCR products were electrophoresed on a 2% TBE agarose gel, along with a 100 bp DNA marker, and visualized. A 466 bp fragment was observed in individuals with the GG genotype, three fragments of 185, 281, and 466 bp indicated the GA genotype, and two fragments of 185 and 281 bp corresponded to the AA genotype (Figure 1).

The genotyping of the ACE2 rs2285666 polymorphism was performed using PCR-based restriction fragment length polymorphism (PCR-RFLP). After digestion, fragments of 281 and 185 bp identified the A allele, and a 466 bp band identified the G allele. After electrophoresis, the bands were visualized on Safe-stained agarose gels, and all genotypes were confirmed.
The genotyping of the ACE2 rs2285666 polymorphism was performed using PCR-based restriction fragment length polymorphism (PCR-RFLP). After digestion, fragments of 281 and 185 bp identified the A allele, and a 466 bp band identified the G allele. After electrophoresis, the bands were visualized on Safe-stained agarose gels, and all genotypes were confirmed.

4.2. Characteristics of the Study Population

The present study included 372 participants, of which 201 were COVID-19 confirmed cases, and 171 were negative or non-diagnostic results based on the RT-qPCR test. Among the samples, 173 (46.6%) were women, and 199 (53.4%) were men. The patients’ ages ranged from 0 to 96 years, with an average age of 56 years (Table 1).

Table 1.

Demographic, Clinical Manifestations, and Underlying Diseases of the Patients (n = 372) a, b

VariablesAll SamplesPositive Group (n = 201)Negative Group (n = 171)
Male199 (53)9499
Female173 (46)9677
Mean age (y)51 (13)5646
Lur198 (53)100 (50)98 (57)
Fars105 (28)76 (38)29 (17)
Turkmen56 (15)21 (10)35 (21)
Balooch12 (3)4 (2)8 (5)

4.3. Distribution of rs2285666 Genotypes

The genotype frequencies in SARS-CoV-2 positive and negative patients are shown in Table 2. Among females, the frequencies of the GG, GA, and AA genotypes were 61%, 35%, and 4% in the control group and 56%, 29%, and 15% in COVID-19 positive cases, respectively. The frequency of G and A genotypes was 89% and 11% in the control group, and 25% and 75% in the patient group, respectively. The comparison of genotype distributions between patients and the control group across different ethnicities is presented in Table 2. The genotype frequencies differed significantly among Fars participants, as well as across the entire population, for both males (P = 0.029 and P = 0.025, respectively) and females (P = 0.050 and P = 0.050, respectively).

Table 2.

Distribution of rs2285666 Genotypes in Male and Female Individuals Among Different Ethnicities a, b

Ethnicity and GenotypesNegative Group Positive GroupP-Value c
Female
Lur0.933
GG22 ( 56.4)30 (55.6)
GA15 (38.5)22 (40.7)
AA2 (5.1)2 (3.7)
Fars0.050
GG11 (73.3)17 (54.8)
GA4 (28.6)3 (9.3)
AA0 (0.0)11 (34.4)
Turkmen0.186
GG11 (68.8)7 (63.6)
GA5 (31.3)2 (18.2)
AA0 (0.0)2 (18.2)
Balooch-
GG3 (42.9)-
GA4 (57.1)
AA-
Total0.050
GG47 (61)54 (56.2)
GA27 (35.1)28 (29.2)
AA3 (3.9)14 (14.6)
Male
Lur0.594
G50 (84.7)37 (80.4)
A9 (15.3)9 (19.6)
Fars0.029
G23 (100.0)25 (69.4)
A0 (0.0)11 (30.6)
Turkmen0.578
G17 (89.5)8 (72.7)
A2 (12.5)3 (27.3)
Balooch0.250
G4 (100)0 (0.0)
A0 (0.0)1 (100.0)
Total0.025
G94 (89.5)70 (74.5)
A11 (10.5)24 (25.5)

The allele frequencies in the control group, patients, and total population were 0.83, 0.72, and 0.77 for the G allele, and 0.17, 0.28, and 0.23 for the A allele, respectively. The population was found not to be in HWE in both the positive and negative groups (P ≤ 0.05) (Table 3). Considering that ACE2 is an X-linked gene, females were studied to determine the inheritance pattern of the polymorphism. For both co-dominant (P = 0.046) and recessive (P = 0.014) inheritance types, there was a significantly higher likelihood of COVID-19 infection among individuals with the AA genotype, with OR of OR = 4.06 (1.10 - 15.00) and OR = 4.21 (1.16 - 15.24), respectively (Table 4). Additionally, there was a significantly higher frequency of the A allele in COVID-19 patients (P = 0.007).

Table 3.

Allele Frequencies of rs2285666 among Participants a

AllelesFrequencyHardy-Weinberg (P-Value)
ControlPatientTotalControlPatientTotal
G0.830.720.7710.00570.029
A0.170.280.23
Table 4.

Association of rs2285666 Genotypes with COVID-19 a, b

Model and GenotypesSNP Association with Response Group (n = 173, Crude Analysis)
Control GroupPatient GroupOR (95% CI)P-ValueAICBIC
Co-dominant0.046237.6247.1
G/G47 (61)54 (56.2)1.00
G/A27 (35.1)28 (29.2)0.90 (0.47 - 1.74)
A/A3 (3.9)14 (14.6)4.06 (1.10 - 15.00)
Dominant0.52241.3247.6
G/G47 (61)54 (56.2)1.00
G/A-A/A30 (39)42 (43.8)1.22 (0.66 - 2.24)
Recessive0.014235.7242
G/G-G/A74 (96.1)82 (85.4)1.00
A/A3 (3.9)14 (14.6)4.21 (1.16 - 15.24)
Overdominant0.41241.1247.4
G/G-A/A50 (64.9)68 (70.8)1.00
G/A27 (35.1)28 (29.2)0.76 (0.40 - 1.45)
Log-additive0.13239.4245.7
---1.43 (0.90 - 2.27)

4.4. Association of the rs2285666 Polymorphism with Clinical Symptoms and Laboratory variables

The study of clinical symptoms showed that patients exhibited fever (48.2%), cough (42.3%), muscle pain (17.4%), respiratory distress (52.2%), diabetes (12.9%), heart disease (8.4%), and high blood pressure (6.5%). The results showed that female patients with the GG genotype had significantly higher frequencies of fever (P = 0.041), respiratory distress (P = 0.015), and high blood pressure (P = 0.037). There were also significantly higher frequencies of respiratory distress among all female participants (P = 0.003). Among males with the G genotype and all male participants, the patient group showed significantly higher frequencies of fever (P = 0.008 and P = 0.004, respectively) and respiratory distress (P = 0.003 and P < 0.001, respectively) (Table 5).

Table 5.

Association Between the rs2285666 Polymorphism and Clinical Symptoms a, b

VariablesFeverCoughMuscle painRespiratory DistressDiabetesHeart DiseaseBlood Pressure
ACE2 G8790ANoYesNoYesNoYesNoYesNoYesNoYesNoYes
Female
GG
Control32 (68.1)15 (31.9)31 (65.9)16 (34.1)40 (85.1)7 (14.9)32 (68.1)15 (31.9)41 (87.2)6 (12.8)44 (93.6)3 (6.4)47 (100.0)0 (0.0)
Case26 (48.1)28 (51.9)32 (59.2)22 (40.8)42 (77.8)12 (22.2)24 (44.4)30 (55.6)46 (85.2)8 (14.8)48 (88.9)6 (11.1)48 (88.9)6 (11.1)
P- value c0.0410.5320.3160.0151.0000.5110.037
GA
Control15 (55.6)12 (44.4)17 (62.9)10 (37.1)26 (96.3)1 (3.7)20 (74.1)7 (25.9)26 (96.3)1 (3.7)25 (92.6)2 (7.4)25 (92.6)2 (7.4)
Case16 (57.1)12 (42.9)16 (57.1)12 (42.9)23 (82.1)5 (17.9)17 (60.7)11 (39.3)22 (78.6)6 (21.4)26 (92.9)2 (7.1)26 (92.9)2 (7.1)
P- value c1.0000.7790.1920.3720.1001.0001.000
AA
Control2 (66.6)1 (43.4)2 (66.6)1 (43.4)1 (43.4)2 (66.6)1 (43.4)2 (66.6)3 (100)0 (0.0)3 (100.0)0 (0.0)3 (100.0)0 (0.0)
Case5 (35.7)9 (64.3)8 (57.1)6 (42.9)12 (85.7)2 (14.3)5 (35.7)9 (64.3)11 (78.6)3 (21.4)13 (92.9)1 (17.1)13 (92.9)1 (17.1)
P- value c1.0001.0000.3311.0001.0001.0001.000
Total
Control49 (63.6)28 (36.4)50 (64.9)27 (35.1)67 (87.1)10 (11.9)53 (68.8)24 (31.2)70 (90.9)7 (9.1)72 (93.5)5 (6.5)75 (97.4)2 (2.6)
Case47 (48.9)49 (51.1)56 (58.3)40 (41.7)77 (80.2)19 (19.8)46 (47.9)50 (52.1)79 (82.3)17 (17.7)87 (90.6)9 (9.4)87 (90.6)9 (9.4)
P- value c0.0590.5190.2060.0030.1700.5930.203
Male
G
Control67 (71.3)27 (28.7)62 (65.9)32 (44.1)81 (86.2)13 (13.8)63 (67.1)31 (32.9)93 (98.9)1 (1.1)86 (91.5)8 (8.5)88 (93.6)6 (6.4)
Case34 (48.6)36 (51.4)42 (60)28 (40)57 (81.4)13 (18.6)29 (41.4)41 (58.6)64 (91.4)6 (8.6)65 (92.9)5 (7.1)66 (94.3)4 (5.7)
P- value c0.0080.5940.4860.0030.1150.7641.000
A
Control8 (72.7)3 (17.3)4 (36.4)7 (63.6)11 (100.0)0 (0.0)9 (81.8)2 (18.210 (90.9)1 (9.1)11 (100.0)0 (0.0)11 (100.0)0 (0.0)
Case12 (50.0)12 (50.0)7 (29.2)17 (70.8)21 (87.5)3 (12.5)10 (41.7)14 (58.3)21 (87.5)3 (12.5)21 (87.5)3 (12.5)24 (100.0)0 (0.0)
P- value c0.4460.6830.5340.0601.0000.534-
Total
Control75 (71.4)30 (28.6)66 (62.9)39 (37.1)92 (87.6)13 (12.4)73 (69.5)32 (30.5)103 (98.1)2 (1.9)97 (92.4)8 (7.6)99 (94.3)6 (5.7)
Case46 (48.9)48 (51.1)49 (52.1)45 (47.9)78 (82.9)16 (17.1)39 (41.5)55 (58.5)85 (90.4)9 (9.6)86 (91.5)8 (8.5)90 (95.7)4 (4.3)
P-value c0.0040.2120.287< 0.0010.0651.0001.000

In female patients, there was a significant increase in the ESR and a significant decrease in lymphocyte levels among both GG genotype females (P = 0.02 and P < 0.001, respectively) and all female participants (P = 0.031 and P = 0.006, respectively). In male patients with the A genotype, there was a significant increase in neutrophils (P = 0.007) and a significant decrease in lymphocyte levels (P = 0.004). In male patients with the G genotype, the Hb level was significantly decreased (P = 0.016). Among all male patients, there was a significant increase in ESR (P = 0.008) and neutrophils (P = 0.013), along with a significant decrease in RBC count (P = 0.029), Hb (P = 0.001), and HCT (P = 0.004) (Table 6). It should be noted that the results of this study could be influenced by limitations, such as the small number of individuals with the AA and A genotypes, as well as potential confounding factors that may affect laboratory values (e.g., age, medication, severity of COVID-19).

Table 6.

Association Between the rs2285666 Polymorphism and Laboratory Variables a, b

ACE2 PolymorphismESRWBCRBCHbHematocritPlateletNeutrophilLymphocytes
Female
GG
Control20.923 (20.2791)7.927 (4.8733)4.573 (1.1016)13.277 (2.9285)38.650 (7.9791)191.23 (77.563)60.43 (23.714)28.29 (12.669)
Case39.385 (23.3685)8.865 (4.4673)4.177 (.7337)12.437 (2.3279)36.223 (6.0680)203.58 (85.313)68.21 (28.567)16.77 (10.685)
P-value c0.0200.4470.1460.2200.1990.5760.292< 0.001
GA
Control33.889 (21.7913)7.433 (3.7397)4.533 (.4850)12.913 (2.9895)38.853 (4.3536)214.40 (83.494)53.93 (30.115)26.80 (10.605)
Case44.769 (32.2727)6.446 (1.9509)4.375 (.5252)13.546 (1.2079)38.058 (3.9715)212.46 (56.971)65.33 (19.482)25.42 (13.571)
P-value c0.3900.3560.3520.4460.5610.9320.1580.739
AA
Control33.500 (41.7193)7.950 (2.6163)3.750 (.7778)11.100 (3.8184)31.350 (6.7175)285.50 (122.329)62.00 (15.556)25.50 (9.192)
Case37.000 (27.8388)6.120 (2.8700)4.38 (.3962)13.340 (2.0428)38.500 (4.6707)154.20 (92.356)68.40 (23.681)25.60 (21.102)
P-value c0.9150.4730.1900.3330.1580.1740.7030.993
Total
Control26.833 (22.3270)7.738 (4.3031)4.515 (.8978)13.026 (2.9415)38.354 (6.7886)204.97 (82.296)57.95 (25.780)27.55 (11.505)
Case40.881 (26.1331)7.825 (3.8392)4.267 (.6412)12.888 (2.0333)37.089 (5.2866)203.09 (77.335)67.21 (25.071)20.42 (13.137)
P-value c0.0310.9150.1370.7970.2940.9050.0740.006
Male
G
Control27.893 (21.1351)9.649 (11.8762)4.772 (.8821)15.026 (2.9615)41.681 (7.2753)205.70 (83.349)66.46 (23.240)21.83 (13.124)
Case37.906 (24.1546)7.769 (4.3846)4.533 (.6895)13.758 (1.8831)39.079 (5.2343)195.53 (94.043)74.09 (13.042)20.04 (11.404)
P-value c0.0950.3210.1620.0160.0570.5880.0570.492
A
Control18.400 (24.2961)6.188 (2.1490)5.150 (.7597)15.013 (2.1735)41.725 (4.6879)232.88 (50.309)55.13 (23.314)29.25 (8.763)
Case52.000 (40.3044)9.500 (5.4299)4.442 (.7821)13.250 (2.2857)37.950 (5.6843)227.92 (131.524)77.58 (8.878)16.50 (8.274)
P-value c0.1130.1200.0600.1020.1370.9210.0070.004
Total
Control25.730 (21.3758)9.200 (10.6388)4.835 (.8433)15.070 (2.7676)41.845 (6.7401)210.15 (76.883)65.44 (22.684)22.56 (12.735)
Case41.511 (29.1194)8.636 (5.1439)4.514 (.7103)13.568 (2.0372)38.278 (6.3821)206.14 (101.863)74.18 (14.051)20.10 (12.899)
P-value c0.0080.7100.0290.0010.0040.8080.0130.295

5. Discussion

The clinical outcomes of SARS-CoV-2 infection vary widely among individuals, ranging from mild respiratory symptoms to multi-organ failure, which can lead to mortality (21). ACE2 serves as the functional receptor for viral attachment and penetration into host cells. In females, the frequencies of GG, GA, and AA genotypes were 61%, 35%, and 4% in the control group, and 56%, 29%, and 15% in COVID-19 positive patients, respectively. Among male participants, the prevalence of G and A genotypes was 89% and 11% in the control group, and 25% and 75% in the patient group, respectively. The polymorphism frequencies were not in HWE in both the positive and negative groups (P < 0.05). Logistic regression analysis revealed a significant difference for the AA genotype in both co-dominant and recessive inheritance models, with OR of OR = 4.06 (1.10 - 15.00) and OR = 4.21 (1.16 - 15.24), respectively. This study identified a strong association between the ACE2 rs2285666 AA or A genotypes and the likelihood of contracting COVID-19.

Several studies have shown that different SNPs in the ACE2 gene may contribute to the severity of COVID-19 and other diseases, such as hypertension, type 2 diabetes mellitus, and cardiovascular diseases (7). Numerous clinical studies have investigated these SNPs, including rs4830542, rs4240157, rs4646155, rs233575, rs2074192, rs2285666, and rs2106809 (7, 22, 23). Among these characterized SNPs, the G8790A (rs2285666, G > A) polymorphism is the most well-known and has been associated with susceptibility to COVID-19 and its severity (24). Most studies have been limited to in silico analyses, and there are few reports exploring the relationship between clinical outcomes and the rs2285666 polymorphism (12-16). In a recent study, Karakas Celik et al. analyzed ACE I/D, ACE2 rs2106809, and rs2285666 polymorphisms in 155 hospitalized Turkish COVID-19 patients, who were divided into three groups (mild, moderate, and severe). They found no association between these polymorphisms and COVID-19 severity (16).

In our study, we analyzed samples from 372 Iranian participants for the rs2285666 polymorphism. The results revealed differences in the genotype and allele frequencies between SARS-CoV-2-positive individuals and the control group. In contrast to the findings of Karakas Celik et al., we found a significant association between the ACE2 rs2285666 AA genotype or A-allele and the risk of developing COVID-19. Although a significant association was observed, it should be noted that some limitations, such as the small sample size of AA genotype cases, may affect the study. These findings provide promising perspectives for future research.

Srivastava et al. conducted a haplotype analysis and reported significantly lower rates of infection and case-fatality rate (CFR) for the alternate A allele of rs2285666 among Indian populations (13). Additionally, they found that the frequency of this allele was significantly higher in Indian populations compared to European, American, or African populations. However, while their study provided haplotype analysis, it did not establish a direct clinical association between the severity of COVID-19 and the rs2285666 polymorphism (13). In contrast, a cohort study found that the GG genotype or G allele of rs2285666 was significantly associated with higher infection (two-fold) and mortality (three-fold) rates (18). The current study found a significant association between the rs2285666 AA and A genotypes and COVID-19, along with some clinical-laboratory features, which is consistent with previous studies.

Sienko et al. studied 188 SARS-CoV-2-positive patients in two groups: Mild and severe COVID-19 cases that required hospitalization. The results indicated an increased risk of severity associated with specific ACE2 gene polymorphisms. The highest correlation was found for the variants rs2285666 (AA allele, OR = 2.12, P = 0.0189), rs2074192 (TT allele, OR = 2.05, P = 0.0016), rs4646174 (GG allele, OR = 1.93, P = 0.0016), rs4646156 (TT allele, OR = 1.71, P = 0.008), and rs2158083 (TT allele, OR = 1.84, P = 0.0025). The results of this study are consistent with the present study, which suggests the role of AA and A genotypes in COVID-19 infection susceptibility (25).

Najafi and Mahdavi examined 146 hospitalized COVID-19 patients, including 102 with mild and 44 with severe cases. There was no statistically significant relationship between the studied groups and genotype frequencies (P = 0.8472). Although the ACE2 rs2285666 polymorphism was not associated with disease severity, the SNP was significantly associated with ALT, ESR, and P. The results of this study were not consistent with those of the present study (20).

Alimoradi et al. studied 129 samples, including negative controls, mild, and severe groups. The results indicated an association between the ACE2 rs2285666 genotype (GG) and susceptibility to COVID-19 infection (P = 0.008; OR 5.0, 95% CI: 1.4 - 17.8), but no association with infection severity. These results were not consistent with our findings, which showed an association between the AA and A genotypes with the disease (17). The rs2285666 variant regulates alternative splicing and increases the expression level of the ACE2 receptor gene by up to 50%, potentially playing a significant role in SARS-CoV-2 susceptibility. In a study by Li, serum levels of ACE2 were higher in individuals with the AA genotype compared to those with the GG and GA genotypes (26). The overexpression of ACE2 in individuals with the AA genotype may be inconsistent with the current study’s results, which found the AA genotype to be associated with COVID-19. Since the rs2285666 variant is considered to be a splicing enhancer site, measuring serum ACE2 protein could potentially help predict the progression of COVID-19.

In the current study, we analyzed a large number of samples to explore the role of the ACE2 rs2285666 polymorphism in SARS-CoV-2 viral infection and to determine the prevalence of ACE2 rs2285666 among Iranian populations with different ethnic backgrounds. It is important to note that our results reflect a general population with some limitations in sample selection, rather than being the outcome of a highly selective or restricted study. The exact role of the ACE2 rs2285666 polymorphism in COVID-19 disease outcomes remains controversial. This controversy may be attributed to various confounding factors such as disease status, sample collection methods (including underlying health conditions or exposure levels), demographics (e.g., age, sex), and ethnic factors. The sample sizes for some ethnic groups (especially the Balooch) were small, making it challenging to draw reliable conclusions about genotype distribution.

5.1. Conclusions

Our study provides evidence of a relationship between the ACE2 rs2285666 (G8790A) polymorphism and COVID-19 susceptibility. The data from this study can be extended to general populations of different ethnicities, contributing to the understanding of the pathogenicity of SARS-CoV-2. Our findings highlight the need for large-scale studies to identify variations and polymorphisms in the ACE2 gene and determine their role in the susceptibility and severity of COVID-19 across different genetic backgrounds. Some limitations exist, including confounding factors in sample selection and the low number of participants from certain ethnic groups. Despite these limitations, the promising results of this study could guide future research in this area. The authors recommend further studies involving larger and more diverse populations to investigate the impact of this polymorphism on COVID-19 susceptibility and severity.

Acknowledgements

References

  • 1.

    Wu D, Wu T, Liu Q, Yang Z. The SARS-CoV-2 outbreak: What we know. Int J Infect Dis. 2020;94:44-8. [PubMed ID: 32171952]. [PubMed Central ID: PMC7102543]. https://doi.org/10.1016/j.ijid.2020.03.004.

  • 2.

    Harapan H, Itoh N, Yufika A, Winardi W, Keam S, Te H, et al. Coronavirus disease 2019 (COVID-19): A literature review. J Infect Public Health. 2020;13(5):667-73. [PubMed ID: 32340833]. [PubMed Central ID: PMC7142680]. https://doi.org/10.1016/j.jiph.2020.03.019.

  • 3.

    Grasselli G, Zangrillo A, Zanella A, Antonelli M, Cabrini L, Castelli A, et al. Baseline Characteristics and Outcomes of 1591 Patients Infected With SARS-CoV-2 Admitted to ICUs of the Lombardy Region, Italy. JAMA. 2020;323(16):1574-81. [PubMed ID: 32250385]. [PubMed Central ID: PMC7136855]. https://doi.org/10.1001/jama.2020.5394.

  • 4.

    Liu Y, Yan LM, Wan L, Xiang TX, Le A, Liu JM, et al. Viral dynamics in mild and severe cases of COVID-19. Lancet Infect Dis. 2020;20(6):656-7. [PubMed ID: 32199493]. [PubMed Central ID: PMC7158902]. https://doi.org/10.1016/S1473-3099(20)30232-2.

  • 5.

    Zheng Z, Peng F, Xu B, Zhao J, Liu H, Peng J, et al. Risk factors of critical & mortal COVID-19 cases: A systematic literature review and meta-analysis. J Infect. 2020;81(2):e16-25. [PubMed ID: 32335169]. [PubMed Central ID: PMC7177098]. https://doi.org/10.1016/j.jinf.2020.04.021.

  • 6.

    Li X, Xu S, Yu M, Wang K, Tao Y, Zhou Y, et al. Risk factors for severity and mortality in adult COVID-19 inpatients in Wuhan. J Allergy Clin Immunol. 2020;146(1):110-8. [PubMed ID: 32294485]. [PubMed Central ID: PMC7152876]. https://doi.org/10.1016/j.jaci.2020.04.006.

  • 7.

    Badawi S, Ali BR. ACE2 Nascence, trafficking, and SARS-CoV-2 pathogenesis: the saga continues. Hum Genomics. 2021;15(1):8. [PubMed ID: 33514423]. [PubMed Central ID: PMC7844112]. https://doi.org/10.1186/s40246-021-00304-9.

  • 8.

    Imai Y, Kuba K, Ohto-Nakanishi T, Penninger JM. Angiotensin-converting enzyme 2 (ACE2) in disease pathogenesis. Circ J. 2010;74(3):405-10. [PubMed ID: 20134095]. https://doi.org/10.1253/circj.cj-10-0045.

  • 9.

    Onabajo OO, Banday AR, Stanifer ML, Yan W, Obajemu A, Santer DM, et al. Interferons and viruses induce a novel truncated ACE2 isoform and not the full-length SARS-CoV-2 receptor. Nat Genet. 2020;52(12):1283-93. [PubMed ID: 33077916]. [PubMed Central ID: PMC9377523]. https://doi.org/10.1038/s41588-020-00731-9.

  • 10.

    Sarangarajan R, Winn R, Kiebish MA, Bountra C, Granger E, Narain NR. Ethnic prevalence of angiotensin-converting enzyme deletion (D) polymorphism and COVID-19 risk: Rationale for use of angiotensin-converting enzyme inhibitors/angiotensin receptor blockers. J Racial Ethn Health Disparities. 2020;8(4):1-8. [PubMed ID: 32901433]. [PubMed Central ID: PMC7478439]. https://doi.org/10.1007/s40615-020-00853-0.

  • 11.

    Hashemi SMA, Thijssen M, Hosseini SY, Tabarraei A, Pourkarim MR, Sarvari J. Human gene polymorphisms and their possible impact on the clinical outcome of SARS-CoV-2 infection. Arch Virol. 2021;166(8):2089-108. [PubMed ID: 33934196]. [PubMed Central ID: PMC8088757]. https://doi.org/10.1007/s00705-021-05070-6.

  • 12.

    Cao Y, Li L, Feng Z, Wan S, Huang P, Sun X, et al. Comparative genetic analysis of the novel coronavirus (2019-nCoV/SARS-CoV-2) receptor ACE2 in different populations. Cell Discov. 2020;6:11. [PubMed ID: 32133153]. [PubMed Central ID: PMC7040011]. https://doi.org/10.1038/s41421-020-0147-1.

  • 13.

    Srivastava A, Bandopadhyay A, Das D, Pandey RK, Singh V, Khanam N, et al. Genetic Association of ACE2 rs2285666 Polymorphism With COVID-19 Spatial Distribution in India. Front Genet. 2020;11:564741. [PubMed ID: 33101387]. [PubMed Central ID: PMC7545580]. https://doi.org/10.3389/fgene.2020.564741.

  • 14.

    Wang J, Xu X, Zhou X, Chen P, Liang H, Li X, et al. Molecular simulation of SARS-CoV-2 spike protein binding to pangolin ACE2 or human ACE2 natural variants reveals altered susceptibility to infection. J Gen Virol. 2020;101(9):921-4. [PubMed ID: 32538738]. [PubMed Central ID: PMC7654750]. https://doi.org/10.1099/jgv.0.001452.

  • 15.

    Calcagnile M, Forgez P, Iannelli A, Bucci C, Alifano M, Alifano P. Molecular docking simulation reveals ACE2 polymorphisms that may increase the affinity of ACE2 with the SARS-CoV-2 Spike protein. Biochimie. 2021;180:143-8. [PubMed ID: 33181224]. [PubMed Central ID: PMC7834737]. https://doi.org/10.1016/j.biochi.2020.11.004.

  • 16.

    Karakas Celik S, Cakmak Genc G, Piskin N, Acikgoz B, Altinsoy B, Kurucu Issiz B, et al. Polymorphisms of ACE (I/D) and ACE2 receptor gene (Rs2106809, Rs2285666) are not related to the clinical course of COVID-19: A case study. J Med Virol. 2021;93(10):5947-52. [PubMed ID: 34170561]. [PubMed Central ID: PMC8426884]. https://doi.org/10.1002/jmv.27160.

  • 17.

    Alimoradi N, Sharqi M, Firouzabadi D, Sadeghi MM, Moezzi MI, Firouzabadi N. SNPs of ACE1 (rs4343) and ACE2 (rs2285666) genes are linked to SARS-CoV-2 infection but not with the severity of disease. Virol J. 2022;19(1):48. [PubMed ID: 35305693]. [PubMed Central ID: PMC8934128]. https://doi.org/10.1186/s12985-022-01782-6.

  • 18.

    Mohlendick B, Schonfelder K, Breuckmann K, Elsner C, Babel N, Balfanz P, et al. ACE2 polymorphism and susceptibility for SARS-CoV-2 infection and severity of COVID-19. Pharmacogenet Genomics. 2021;31(8):165-71. [PubMed ID: 34001841]. [PubMed Central ID: PMC8415730]. https://doi.org/10.1097/FPC.0000000000000436.

  • 19.

    Sheikhian F, Sadeghi Mofrad S, Tarashi S, Ghazanfari Jajin M, Sakhaee F, Ahmadi I, et al. The impact of ACE2 polymorphisms (rs1978124, rs2285666, and rs2074192) and ACE1 rs1799752 in the mortality rate of COVID-19 in different SARS-CoV-2 variants. Hum Genomics. 2023;17(1):54. [PubMed ID: 37328914]. [PubMed Central ID: PMC10273585]. https://doi.org/10.1186/s40246-023-00501-8.

  • 20.

    Najafi M, Mahdavi MR. Association investigations between ACE1 and ACE2 polymorphisms and severity of COVID-19 disease. Mol Genet Genomics. 2023;298(1):27-36. [PubMed ID: 36255490]. [PubMed Central ID: PMC9579601]. https://doi.org/10.1007/s00438-022-01953-8.

  • 21.

    Zhou P, Yang XL, Wang XG, Hu B, Zhang L, Zhang W, et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Natur J. 2020;579(7798):270-3. [PubMed ID: 32015507]. [PubMed Central ID: PMC7095418]. https://doi.org/10.1038/s41586-020-2012-7.

  • 22.

    Fan XH, Wang YB, Wang H, Sun K, Zhang WL, Song XD, et al. Polymorphisms of angiotensin-converting enzyme (ACE) and ACE2 are not associated with orthostatic blood pressure dysregulation in hypertensive patients. Acta Pharmacol Sin. 2009;30(9):1237-44. [PubMed ID: 19684612]. [PubMed Central ID: PMC4007186]. https://doi.org/10.1038/aps.2009.110.

  • 23.

    Huang W, Yang W, Wang Y, Zhao Q, Gu D, Chen R. Association study of angiotensin-converting enzyme 2 gene (ACE2) polymorphisms and essential hypertension in northern Han Chinese. J Hum Hypertens. 2006;20(12):968-71. [PubMed ID: 17024138]. https://doi.org/10.1038/sj.jhh.1002090.

  • 24.

    Pinheiro DS, Santos RS, Jardim P, Silva EG, Reis AAS, Pedrino GR, et al. The combination of ACE I/D and ACE2 G8790A polymorphisms revels susceptibility to hypertension: A genetic association study in Brazilian patients. PLoS One. 2019;14(8). e0221248. [PubMed ID: 31430320]. [PubMed Central ID: PMC6701835]. https://doi.org/10.1371/journal.pone.0221248.

  • 25.

    Sienko J, Marczak I, Kotowski M, Bogacz A, Tejchman K, Sienko M, et al. Association of ACE2 Gene Variants with the Severity of COVID-19 Disease-A Prospective Observational Study. Int J Environ Res Public Health. 2022;19(19). [PubMed ID: 36231922]. [PubMed Central ID: PMC9564490]. https://doi.org/10.3390/ijerph191912622.

  • 26.

    Li YY. Lack of association of ACE2 G8790A gene mutation with essential hypertension in the Chinese population: A meta-analysis involving 5260 subjects. Front Physiol J. 2012;3. https://doi.org/10.3389/fphys.2012.00364.