- 1Laboratorio de Gerociencias, Dirección General, Medicina de Rehabilitación, Laboratorio de Infectología, Departamento de Reconstrucción Articular, Laboratorio de Medicina Genómica, Laboratorio Facilitador, Instituto Nacional de Rehabilitación Luis Guillermo Ibarra Ibarra, Secretaría de Salud, Mexico City, Mexico
- 2Postgrado en Biología Experimental, Dirección de Ciencias Biológicas y de la Salud (DCBS), Universidad Autónoma Metropolitana Iztapalapa, Mexico City, Mexico
- 3Departamento de Nutrición Animal, Departamento de Inmunogenética, Instituto Nacional de Ciencias Médicas y Nutrición Salvador Zubirán, Secretaría de Salud, Mexico City, Mexico
- 4Centro de Innovación Médica Aplicada, Subdirección de Epidemiología e Infectología, Hospital General Dr. Manuel Gea González, Mexico City, Mexico
- 5Departamento de Biología Molecular y Unidad de Cuidados Intensivos, Instituto Nacional de Cardiología Ignacio Chavez, Mexico City, Mexico
- 6Nuevo Hospital General Delegación Regional Sur de la Ciudad de México ISSSTE, Mexico City, Mexico
- 7Laboratorio de Biología Celular y Tisular, Laboratorio de Embriología, Escuela Militar de Medicina, Universidad del Ejército y Fuerza Aérea, Ciudad de México, Mexico
- 8Servicio de Cirugía General, Hospital Central Norte Petróleos Mexicanos (PEMEX), Mexico City, Mexico
- 9Programa de Investigadoras e investigadores por México del Consejo Nacional de Ciencia y Tecnología (CONACYT), Instituto Nacional de Cancerología, Mexico City, Mexico
- 10Instituto Nacional de Perinatología Isidro Espinosa de los Reyes, Mexico City, Mexico
- 11Independent Consultant, San Antonio, TX, United States
- 12Instituto de Investigación en Ciencias Biomédicas, Centro Universitario de Ciencias de la Salud, Universidad de Guadalajara, Guadalajara, Mexico
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the current coronavirus disease 2019 (COVID-19) pandemic, affecting more than 219 countries and causing the death of more than 5 million people worldwide. The genetic background represents a factor that predisposes the way the host responds to SARS-CoV-2 infection. In this sense, genetic variants of ACE and ACE2 could explain the observed interindividual variability to COVID-19 outcomes. In order to improve the understanding of how genetic variants of ACE and ACE2 are involved in the severity of COVID-19, we included a total of 481 individuals who showed clinical manifestations of COVID-19 and were diagnosed by reverse transcription PCR (RT-PCR). Genomic DNA was extracted from peripheral blood and saliva samples. ACE insertion/deletion polymorphism was evaluated by the high-resolution melting method; ACE single-nucleotide polymorphism (SNP) (rs4344) and ACE2 SNPs (rs2285666 and rs2074192) were genotyped using TaqMan probes. We assessed the association of ACE and ACE2 polymorphisms with disease severity using logistic regression analysis adjusted by age, sex, hypertension, type 2 diabetes, and obesity. The severity of the illness in our study population was divided as 31% mild, 26% severe, and 43% critical illness; additionally, 18% of individuals died, of whom 54% were male. Our results showed in the codominant model a contribution of ACE2 gene rs2285666 T/T genotype to critical outcome [odds ratio (OR) = 1.83; 95%CI = 1.01–3.29; p = 0.04] and to require oxygen supplementation (OR = 1.76; 95%CI = 1.01–3.04; p = 0.04), in addition to a strong association of the T allele of this variant to develop critical illness in male individuals (OR = 1.81; 95%CI = 1.10–2.98; p = 0.02). We suggest that the T allele of rs2285666 represents a risk factor for severe and critical outcomes of COVID-19, especially for men, regardless of age, hypertension, obesity, and type 2 diabetes.
Introduction
The ongoing coronavirus disease 2019 (COVID-19) pandemic has affected around 336 million individuals, causing the death of nearly 5,560,718 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infected subjects around the world (1). Current evidence shows that COVID-19 patients experience clinical manifestations ranging from asymptomatic to severe pneumonia with multiple organ failure (2–4). The severity of COVID-19 and its clinical manifestations and outcomes are related to the internalization mechanism of the virus into the host cell, host genetics variants, advanced age, gender, and comorbidities (5). In this sense, the homeostasis of the renin–angiotensin system is another risk factor underlying the pathogenesis of COVID-19 because angiotensin-converting enzyme 2 (ACE2) is the predominant receptor by which the SARS-CoV-2 virus enters and infects cells (6, 7). An altered ACE/ACE2 expression ratio could contribute to severe outcomes in COVID-19 patients (8), as it does for cardiovascular diseases (9).
Recent studies have suggested the association of genetic polymorphisms of ACE and ACE2 with the case rate of COVID-19. However, the clinical implication of ACE genetic variants in the severity and prognosis of COVID-19 remains unclear. The insertion/deletion (I/D) polymorphisms play a pivotal role in cardiovascular and respiratory diseases; for instance, the D/D genotype has been associated with SARS progression (10, 11). Moreover, it has been associated with poor clinical outcomes of acute respiratory distress syndrome (ARDS), where individuals with the D/D genotype had significantly higher mortality than those who carry the I/I genotype (12). In addition, different reports have proposed that ACE2 polymorphisms rs2285666 (G870A) could modulate the susceptibility to SARS-CoV-2 infection by contributing to higher expression of ACE2 receptor (13, 14).
Furthermore, several reports have shown that ACE2 gene polymorphism is related to acute lung injury, making COVID-19 patients significantly prone to develop ARDS (12, 14, 15). To improve the understanding of how genetic variations and biological mechanisms are involved in COVID-19 severity, we consider that the association of case fatality rate and genetic variants of ACE/ACE polymorphism could represent a strategy to identify possible SNPs as susceptibility and prognostic markers in the Mexican population.
Material and Methods
Setting and Participants
From June 2020 to March 2021, we performed a multicenter cross-sectional study in which the following institutions were included: the Instituto Nacional de Rehabilitación Luis Guillermo Ibarra-Ibarra, Instituto Nacional de Cardiología Ignacio Chávez, Hospital Central Militar, Hospital Central Norte Pemex, and Hospital General Dr. Manuel Gea González. A non-probability sampling method was conducted, as the participants were patients at the COVID-19 triage facilities.
All subjects were evaluated with the following comprehensive clinical procedures: from each patient, oxygen saturation levels and ferritin, D-dimer, lactate dehydrogenase, and C-reactive protein (CRP) were recorded. Ferritin and D-dimer levels were measured in plasma with automatic analyzers (DxI 800 Beckman Coulter, COBAS pro e-801, COBAS integra 400 plus, and ACL TOP 550 CTS Werfem, COBAS integra 400 plus, respectively), while lactate dehydrogenase and C-reactive protein (CRP) were measured in serum using Olympus AU 800, COBAS pro c-503, COBAS integra 400 plus, and DxC 700 AU Beckman Coulter.
Ethical Statement
This study was conducted following the good clinical practice and the Declaration of Helsinki. Informed consent was obtained from each participant before entering the study. The study was approved by the ethics committee of the Instituto Nacional de Rehabilitación LGII (INR-LGII: 17/20). The inclusion criteria were age ≥18 years and a positive test for SARS-CoV-2 infection by RT-PCR. The exclusion criteria were incomplete clinical history and related individuals.
Outcomes
According to the criteria of Gandhi et al. (16), we defined each COVID-19 patient group as mild, severe, and critical. Mild group included those ambulatory subjects with oxygen saturation level ≥94%, and with symptoms such as fever, headache, fatigue, odynophagia, cough, rhinorrhea, diarrhea, anosmia, or dysgeusia, with or without dyspnea or pneumonia, and not requiring hospitalization. Severe status was defined as those hospitalized individuals with saturation levels <94%, and any of the following symptoms: tachypnea (FR > 30 bpm), pulmonary infiltrate >50%, and dyspnea for small efforts. Finally, the critical group were those patients requiring invasive mechanical ventilation who could course with shock and multiorgan failure. In addition, we stratified patients by their requirement of oxygen therapy in two groups: those who did not need any type of oxygen therapy and those who had oxygen requirement by nasal tips, mask, or intubation.
Genotyping
Genomic DNA was obtained from peripheral blood and saliva using a commercial kit (QIAamp DNA blood Mini kit, Qiagen, Germany), and its quality was verified in agarose gels stained with SYBR® Green I nucleic acid gel stain (Invitrogen, CA, USA). Then, the DNA concentration was measured and adjusted to 20 ng/µl with a NanoDrop spectrophotometer (Thermo Fisher Scientific, CA, USA).
A systematic review of ACE and ACE2 gene SNPs was performed, from which we selected the SNPs with a minor allele frequency (MAF) ≥20% according to 1000 Genomes Project or Hapmap Project in the Mexican population (MXL) with previous reports showing significant association with other diseases.
The I/D polymorphism was genotyped by real-time polymerase chain reaction, coupled to a high-resolution melting curve (HRM) analysis allowing the identification of the dissociation temperature for the deletion and insertion amplicon. The DNA was denatured at 95°C for 5 min, followed by 40 cycles of denaturation at 95°C for 25 s, annealing at 67°C for 60 s, and extension at 72°C for 75 s, melt ramp from 77°C to 95°C. This protocol was implemented in the RotorGen Q thermal cycler (Qiagen, Germany). The genotype for this variant was assigned by HRM analysis and confirmed by the melting curve. Additionally, the rs4344 of ACE and rs2285666 and rs2074192 of ACE2 were analyzed. These SNPs were determined with TaqMan genotyping assays on a StepOne Real-Time PCR equipment (Thermo Fisher Scientific, CA, USA).
Statistical Analysis
The normality of the distribution of the variables was evaluated. Kruskal–Wallis was used for comparing non-parametric continuous variables between studied groups, and the results were described using the median and the interquartile range (IQR). For the categorical variables, the chi-squared test was performed. For all tests, a value of p <0.05 was considered statistically significant. Hardy–Weinberg expectations (HWEs) were assessed for all polymorphisms in the mild group, and the linkage disequilibrium (LD) between variants of the same gene using HaploView software V4.2 (17). A logistic regression analysis was used to evaluate the association between genetic variants and COVID-19 outcomes in the three main inheritance models, i.e., codominant, dominant, and recessive, adjusted by age, sex, hypertension, type 2 diabetes, and obesity. The final models were evaluated using the Hosmer–Lemeshow goodness-of-fit test. The association between SNPs and clinical features was assessed by comparing their distribution among alleles and genotypes by Kruskal–Wallis test and stratified by disease outcome. The analysis was performed using the STATA v.13 statistical package (StataCorp, TX, USA).
Results
Patients
From June 2020 to March 2021, 489 individuals with COVID-19 were selected. Nonetheless, eight patients did not have complete clinical data. Therefore, a total of 481 patients were included in this study. We classified the study population as 31% (149) mild, 26% (125) severe, and 43% (207) critical, of which 86 subjects died. The median age was 51 (IQR = 43–63), and a total of 60% (290) were male subjects. In the overall population, the common symptoms were cough, headache, myalgia, and dyspnea. Most individuals had at least one coexisting illness such as hypertension and type 2 diabetes. We observed significant differences in the majority of the subjects’ clinical characteristics depending on the group they belong to (Table 1). However, no significant differences were observed for abdominal pain (p > 0.05).
Association Between ACE and ACE2 Polymorphisms With COVID-19 Outcomes
Allelic and genotype frequencies were calculated (Table 2). For ACE gene, I/D alleles in the total sample were 61% for I and 39% for D allele; in addition, for rs4344 polymorphism, the MAF was 39% (A allele). Only 42% of patients with mild illness carried the II genotype, while 33% of severe patients and 39% of critical patients showed this genotype. These frequencies were also observed for the rs4344 variant and were not significantly different between study groups.
When we compared the genotypes frequencies of the evaluated SNPs between mild, severe, and critical groups, we found significant differences in rs2285666 polymorphism and marginally difference in rs2074192 polymorphism both of ACE2 gene (p = 0.009, p = 0.05, respectively). When comparing mild illness and the requirement for oxygen therapy (severe or critical illness), we found significant differences in rs2285666 and rs2074192 ACE2 variants genotypic frequencies (p = 0.002 and p = 0.01, respectively). Among those patients that needed oxygen therapy, 36% carried the T/T genotype of the rs2285666 gene variants in comparison to 23% of the patients who did not need it, while for the variant of the rs2074192, the T/T genotype showed a 35% in comparison to 28% of the patients who did not require it. For the alleles frequencies, the MAF of rs2285666 variant was 41% (T allele), while the MAF of rs2071192 was 42% (T allele). We did not detect any significant difference in the allelic frequencies of both polymorphisms between outcomes.
The genotypes of ACE I/D and rs4344 variants were in Hardy–Weinberg equilibrium. However, the genotypes of ACE2 were not in equilibrium. We observed a strong LD between both ACE gene variants, showing a D′ of 0.97. Regarding the ACE2 gene variants, the LD was slightly lower, having a D′ of 0.88.
Multiple Logistic Regression Analysis
In the logistic regression analysis, we did not find a statistically significant association of ACE variants with COVID-19 outcomes. Notwithstanding, for the C/T genotype of the rs2285666 polymorphism, we found an OR = 0.52 (95%CI = 0.29–0.94), and for T/T genotype, an OR = 1.66 (95%CI = 1.01–2.73) in the codominant model with critical outcome (Supplementary Material). After adjusting by age, hypertension, sex, type 2 diabetes, and obesity, we found a significant positive association between the T/T genotype in the codominant model of rs2285666 polymorphism and critical outcome (OR = 1.83; 95%CI = 1.01–3.29), and for oxygen requirement (OR = 1.76; 95%CI = 1.01–3.04). Moreover, for this same ACE2 genetic variant but under a recessive model, we observed an OR = 1.89 (95%CI = 1.06–3.35) for critical outcome and an OR = 1.80 (95%CI = 1.06–3.05) for oxygen requirement. This positive association with critical outcome and oxygen requirement found with the T/T genotype of ACE2 rs2285666 polymorphism was also observed with the T allele (OR = 1.58, 95%CI = 1.09–2.30; and OR = 1.52, 95%CI = 1.08–2.14). Finally, the T allele of this ACE2 variant was marginally associated with a higher risk of severe outcome (OR = 1.45; 95%CI = 0.99–2.13) (Table 3).
Given that ACE2 gene polymorphisms are in the X chromosome, an allele analysis stratified by sex was performed. We observed that the positive association found with the rs2285666 T allele and severe illness was maintained and even had a higher magnitude and significance among men than the one observed in the whole study population (OR = 1.72; 95%CI = 1.02–2.89; p =0.03). The same happened with the association seen between ACE2 rs228566 T allele and critical illness, and with oxygen requirement, which for men showed an OR = 1.81 (95%CI = 1.10–2.98) and OR = 1.77 (95%CI, 1.12–2.83), respectively (Table 4).
We further evaluated the potential to identify patients who needed oxygen therapy using the predicted values of the multivariate logistic regression model of the associated polymorphism of ACE2 by receiver operating characteristic (ROC) curve analysis. We observed that with this model, we can categorize patients that required oxygen therapy from those who did not (Figure 1).
Figure 1 ROC curve for the association model of rs2285666 ACE2 polymorphism with oxygen requirement.
Association of ACE and ACE2 Genes Polymorphisms With COVID-19 Clinical Biomarkers and Oxygen Saturation Levels
We performed a refined approach to determine the behavior of the main COVID-19 clinical biomarkers and oxygen saturation level in our study population. In this sense, we could observe a significant reduction in oxygen saturation with increasing disease severity. Furthermore, statistically significant differences were observed among outcomes for COVID-19 biomarkers (p < 0.001) (Supplementary Table S2). D-Dimer and ferritin plasma concentrations showed significantly higher values in critical COVID-19 subjects (679.5 ng/ml for D-dimer and 619 ng/ml for ferritin) than in mild disease patients. We also evaluated lactate dehydrogenase and CRP in serum. A clear tendency for increased levels of these biomarkers was observed among outcomes.
We evaluated the distribution of clinical features and biomarkers among the four genotypes polymorphisms and COVID-19 outcomes; nevertheless, we only found significant differences with the ACE2 rs2285666 polymorphism. In the mild group, we found a higher amount of D-dimer in T allele carriers (p50 = 323 ng/ml, IQR = 217–501 ng/ml) in comparison to C allele carriers, and for C/T genotype (p50 = 360 ng/ml, IQR = 214.5–501 ng/ml) and T/T genotype (p50 = 285, IQR = 217,430 ng/ml) carriers in comparison to C/C patients. Regarding ferritin levels, differences among genotypes were also observed among mild, severe, and critical patients. In this context, the lowest level of ferritin was observed in the carriers of the heterozygote genotype. Moreover, CPR was significantly increased in T allele carriers of the mild group (p50 = 4.95 mg/ml, IQR = 2–13.4 mg/ml), and in the critical patients with the T/T genotype (p50 = 26.55 mg/ml, IQR = 16.61–113.19 mg/ml). In the severe group, we only found significant differences in lactate dehydrogenase levels among ACE2 rs2285666 genotypes, showing the highest level those patients carrying the T/T genotype (p50 = 322.8 U/L, IQR = 213.3–439 U/L) (Supplementary Table S3).
Discussion
A total of 481 individuals with COVID-19 were included in this study, 22% of them were in a severe stage and 86 subjects died with a median age of 64.5 years. These findings agree with those reported by Zhang et al., since he observed that patients older than 60 years were more prone to develop complications with fatal outcomes (18). Additionally, COVID-19 has been reported to be more severe in male individuals (19), which is in agreement with our results. The 67% of individuals who developed a critical outcome were male, and 28 (67%) of the total deaths were attributable to them. This difference in gender susceptibility has been reported in previous studies performed in Italy (20), Spain, China, Germany, the UK, and South Korea (21).
A recent meta-analysis showed that fever, cough, fatigue, and dyspnea were the most prevalent symptoms of COVID-19 (22). On the other hand, in the present work, the majority of symptoms were significantly different among the outcomes (p < 0.001); cough and dyspnea were the highlights in critical disease conditions. A previous report from the WHO–China Joint Mission on COVID-19 showed that 87.9% of confirmed COVID-19 cases presented fever, while in our study, only 32% had fever. Interestingly, 33% of mild patients in our study presented fever, while 41% of those who developed a critical illness had fever. Cough is another common symptom with a previously reported prevalence of 67.7% (4). In our study, we observed a 67% of frequency; however, in critical cases, it was major than previously reported (76%).
Frequent comorbidities observed in individuals with COVID-19 include hypertension, type 2 diabetes, and cardiovascular disease, being serious risk factors in severely ill individuals compared with non-severely ill patients (23). In our study, critically ill adults with COVID-19 presented type 2 diabetes (28%), obesity (27%), and hypertension (23%). These findings are consistent with the published data by Wu et al. (24) and Grasselli et al. (25).
Some risk factors have been described for COVID-19 outcomes; however, understanding the role of host genetic variants for risk or protective effect could provide some insights into COVID-19 outcomes. Differences in the severity of COVID-19 due to sex have already been reported where men are the most vulnerable to severe or fatal COVID-19 outcomes (26, 27). Possible explanations include exposure to environmental factors such as smoking, diet, and physical activity, and genetic factors associated with the mechanisms of infection by SARS-CoV-2 including the genes of the ACE/ACE2 pathways.
Although there are several studies stating the association of different polymorphisms of ACE and ACE2 with hypertension and other cardiovascular diseases, there is no consistent evidence describing the role of rs2285666 and rs2074192 of ACE2 gene or of I/D and rs4344 polymorphisms of ACE gene in COVID-19 outcomes.
Regarding ACE I/D polymorphisms, the frequencies observed in our study were similar to those reported in the Mexican Mestizo population by Vargas-Alarcón et al. They indicated that the I/D polymorphism distinguishes the Amerindian population from other populations, suggesting that the I/D polymorphism could be a distinctive genetic susceptibility marker for some diseases in Amerindian population (28). Sarangarajan et al. reported the prevalence of I/D in different populations including South Americans with I/D prevalence similar to our findings; it is important to state that I/D genotypes prevalence were different among geographical locations (29). This evidence points out the consistency of our results regarding ACE I/D variant frequencies.
Some studies of ACE gene in the Mexican population have reported the I/D polymorphism as a genetic marker of susceptibility for cardiovascular disease, diabetes, and hypertension diseases (28, 30). This evidence suggests that the I/D polymorphism of ACE could impact COVID-19 outcomes (31, 32). A recent study by Gómez et al. associated the D/D genotype of ACE only with hypertensive individuals who presented a severe COVID-19 outcome (12). Conversely, our results showed no significant association of this variant with COVID-19 outcomes, but we found a significant association with hypertension that could be an indirect effect of COVID-19 complications (see Supplementary Material). In agreement with our results, Karakas et al. reported that I/D was not associated with the severity of COVID-19 infection in Turkey populations (33). Nevertheless, it is important to point out that the D/D genotype of ACE gene has a role on the renin–angiotensin system, by increasing ACE levels and angiotensin-inactivating AT-1 receptor, and downstream pathophysiological effects. Therefore, future studies should not be discarded. In the same way, rs4344 polymorphism has been associated with hypertension (34) and angiotensin-converting enzyme inhibitors-related cough (15). Here, we observed the same results; nonetheless, no association with critical COVID-19 outcomes was identified.
The ACE2 rs2285666 variant falls in the third intron of the gene, affecting gene expression by an alternative splicing mechanism (7). Asselta and colleagues reported that substitution of C for T increased the strength of the splice site in 9.2%, having as result a higher expression of the ACE2 protein (35). Gemmati and Tisato reported that the ACE2 gene is a “first genetic gateway” involved in infection, severity, and COVID-19 (26); however, differences among populations and geographic regions could exist, as reported by Shikov et al. (36). Moreover, ACE2 has been associated with hypertension, heart failure, and diabetes. Falahi and Kenarkoohi have already suggested that more studies focused on the role of ACE2 gene in the pathogenesis and outcomes of COVID-19 are needed due to the existing controversies whether ACE2 genotypes might explain the differences in the infection severity (37). Our results showed that the T allele of rs2285666 is a risk factor of severe and critical outcomes in COVID-19, especially in men. Choudhary et al. described the possible physiological roles of ACE2 in COVID-19, suggesting that the study of ACE2 variants will help understand the pathophysiology of this disease (14). Moreover, Benetti et al. reported in a cohort of COVID-19 patients that genetic variants of ACE2 have an impact on its protein function, which reinforces the hypothesis that at least some variants identified, or the cumulative effect of them, might confer different susceptibility to infection and progression of COVID-19 (38). Our results indicate an association between the T/T genotype of ACE2 rs2285666 (G8790A) and an increased risk for severe–critical outcomes (oxygen requirement) in individuals with COVID-19. Srivastava and colleagues stated that genetic variations of ACE2 affect the susceptibility to COVID-19 and found a lower infection rate in the carriers of the T allele in Indian populations. However, it is important to mention that there are many factors involved in the SARS-CoV2 infection, and ACE2 might not be the only gene involved (39). For the rs2074192 ACE2 polymorphism, Shikov et al. (36) reported similar frequencies between controls, mild, and severe outcome of COVID-19 with the ones presented here. However, a recent study by Cafiero et al. (40) reported in 104 subjects of Italy that even though genotypes frequencies were similar to our results (frequencies CC>TT>CT), they suggested that rs2074192 could predict the clinical outcome of COVID-19. We did not find a significant association between this ACE2 variant and disease outcome; nevertheless, this could be explained by the difference in sample sizes and the difference in the genetic background among both populations.
Regarding risk differences by sex, Gemmati et al. reported that an unbalance of ACE/ACE2 might show marked differences in the outcomes of COVID-19 in both sexes, given that this unbalance could induce a higher inflammatory mediators/receptors expression; hence, men might show a worse clinical scenario than women, as women could activate a mosaic advantage due to their X-heterozygosity (27). After all, many genes associated with the regulation of the immune system are located in the X chromosome. In that sense, ACE2 gene locus is Xp22.2, and it has been reported that the rs2285666 T/T genotype of ACE2 decreased the gene expression level by up to 50% compared with the C/C genotype in Italian and other European populations (27, 38). Supporting this notion, we found that men had a higher risk of severe and critical COVID-19 disease than women carrying the T allele of this ACE2 polymorphism.
Among the clinical features evaluated in the present study that have been used as biomarkers for COVID-19, we could observe that ACE2 genetic variant had an impact on D-dimer, ferritin, lactate dehydrogenase, and CRP levels. Even though the direct effect of ACE2 on these biomarkers levels is not quite clear, the difference showed in these COVID-19-associated biomarkers supports the notion that ACE2 rs228566 polymorphism could be consider as a genetic susceptibility marker for COVID-19 outcome.
Another factor important to mention is that there are no previous studies of ACE2 genetic variants in the Mexican population with COVID-19, although the frequency difference of these variants among different populations could affect the association of ACE2 with COVID-19 illness. However, not only the genetic background might explain the differences in the outcomes of COVID-19, but also environmental risk factors such as smoking, drinking habits, and personal hygiene could be implicated.
The strengths of this study include a representative and multicentric hospital-based from Mexico City sample, obtaining clinical information, and the diagnosis of COVID-19 by RT-PCR test and blood sample for DNA extraction. Additionally, the patient outcomes were classified by Gandhi criteria. This adequate classification allowed us to stratify patients and to identify the association of genetic variants with the outcomes of individuals with COVID-19. However, the limitations of this study include that we could not have lifestyle information such as smoking habits, which has been reported to increase the severity of COVID-19 outcomes in men. Another limitation worth mentioning is that we could not evaluate other genes involved in ACE and ACE2 pathways that could be contributing to the increased risk for severe COVID-19 outcome seen for the T allele of ACE2 rs228566 polymorphism.
In conclusion, our results show that the T allele of rs2285666 is associated with more severe outcomes of COVID-19, particularly in men, regardless of their age or the presence of hypertension and type 2 diabetes. This genetic variant could be useful as a prediction and susceptibility marker in the prognosis of COVID-19, which will also help to personalize the COVID-19 treatment.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics Statement
This study was conducted following the good clinical practice and the Declaration of Helsinki. Informed consent was obtained from each participant before entering the study. The study was approved by the ethics committee of (INR-LGII: 17/20). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
ALR, GAMN, and CP: study conception and design. LEMG, CSA, GGS, GVZ, and PVC: literature search. RFC, LLJ, RVJ, JRH, RVV, YRS, RBD, FMR, MMM, MCR, DZA, EBG, LBS, DZA, and ICZ: included patients. LEMG, SHD, DDS, LRT, AMC, and VLT: acquisition of data. LEMG, CCR, OHG, JMRP, SOP, BHL, JMF, and GRV: performed experiments. LMG, GMN, and JJM: data analysis and interpretation. CMA: figures. LEMG, CMA, and ALR: drafting the article. GVA, GJ, FMV, GMN, and ALR: critical revision of the article and final approval of the version to be published. All authors contributed to the article and approved the submitted version.
Funding
This study was funded by the Consejo Nacional de Ciencia y Tecnología; CONACYT 312513 SARS-COV 2.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We gratefully thank the extraordinary effort of healthcare workers who sacrificed their lives while saving patients. In memorial to Pedro Alejandro Lomelí-Mejía, a researcher at the National Institute of Rehabilitation LGII.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.812940/full#supplementary-material
References
1. WHO Coronavirus (COVID-19) Dashboard. Available at: https://covid19.who.int (Accessed January 21, 2021).
2. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical Features of Patients Infected With 2019 Novel Coronavirus in Wuhan, China. Lancet (2020) 395(10223):497–506. doi: 10.1016/S0140-6736(20)30183-5
3. Baj J, Karakula-Juchnowicz H, Teresinski G, Buszewicz G, Ciesielka M, Sitarz E, et al. COVID-19: Specific and Non-Specific Clinical Manifestations and Symptoms: The Current State of Knowledge. J Clin Med (2020) 9(6):1753. doi: 10.3390/jcm9061753
4. Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, et al. Clinical Characteristics of Coronavirus Disease 2019 in China. N Engl J Med (2020) 382(18):1708–20. doi: 10.1056/NEJMoa2002032
5. Monticelli M, Mele BH, Andreotti G, Cubellis MV, Riccio G. Why Does SARS-CoV-2 Hit in Different Ways? Host Genetic Factors can Influence the Acquisition or the Course of COVID-19. Eur J Med Genet (2021) 64(6):104227. doi: 10.1016/j.ejmg.2021.104227
6. Choudhary S, Malik YS, Tomar S. Identification of SARS-CoV-2 Cell Entry Inhibitors by Drug Repurposing Using in Silico Structure-Based Virtual Screening Approach. Front Immunol (2020) 11:1664. doi: 10.3389/fimmu.2020.01664
7. Gheblawi M, Wang K, Viveiros A, Nguyen Q, Zhong JC, Turner AJ, et al. Angiotensin-Converting Enzyme 2: SARS-CoV-2 Receptor and Regulator of the Renin-Angiotensin System: Celebrating the 20th Anniversary of the Discovery of ACE2. Circ Res (2020) 126(10):1456–74. doi: 10.1161/CIRCRESAHA.120.317015
8. Bakhshandeh B, Sorboni SG, Javanmard AR, Mottaghi SS, Mehrabi MR, Sorouri F, et al. Variants in ACE2; Potential Influences on Virus Infection and COVID-19 Severity. Infect Genet Evol: J Mol Epidemiol Evol Genet Infect Dis (2021) 90:104773. doi: 10.1016/j.meegid.2021.104773
9. Pagliaro P, Penna C. ACE/ACE2 Ratio: A Key Also in 2019 Coronavirus Disease (Covid-19)? Front Med (Lausanne) (2020) 7:335. doi: 10.3389/fmed.2020.00335
10. Verma S, Abbas M, Verma S, Khan FH, Raza ST, Siddiqi Z, et al. Impact of I/D Polymorphism of Angiotensin-Converting Enzyme 1 (ACE1) Gene on the Severity of COVID-19 Patients. Infect Genet Evol: J Mol Epidemiol Evol Genet Infect Dis (2021) 91:104801. doi: 10.1016/j.meegid.2021.104801
11. Vaduganathan M, Vardeny O, Michel T, McMurray JJV, Pfeffer MA, Solomon SD. Renin-Angiotensin-Aldosterone System Inhibitors in Patients With Covid-19. N Engl J Med (2020) 382(17):1653–9. doi: 10.1056/NEJMsr2005760
12. Gómez J, Albaiceta GM, García-Clemente M, López-Larrea C, Amado-Rodríguez L, Lopez-Alonso I, et al. Angiotensin-Converting Enzymes (ACE, ACE2) Gene Variants and COVID-19 Outcome. Gene (2020) 762:145102. doi: 10.1016/j.gene.2020.145102
13. Saengsiwaritt W, Jittikoon J, Chaikledkaew U, Udomsinprasert W. Genetic Polymorphisms of ACE1, ACE2, and TMPRSS2 Associated With COVID-19 Severity: A Systematic Review With Meta-Analysis. Rev Med Virol (2022) 8:e2323. doi: 10.1002/rmv.2323
14. Choudhary S, Sreenivasulu K, Mitra P, Misra S, Sharma P. Role of Genetic Variants and Gene Expression in the Susceptibility and Severity of COVID-19. Ann Lab Med (2021) 41(2):129–38. doi: 10.3343/alm.2021.41.2.129
15. Grilo A, Saez-Rosas MP, Santos-Morano J, Sanchez E, Moreno-Rey C, Real LM, et al. Identification of Genetic Factors Associated With Susceptibility to Angiotensin-Converting Enzyme Inhibitors-Induced Cough. Pharmacogenet Genomics (2011) 21(1):10–7. doi: 10.1097/FPC.0b013e328341041c
16. Gandhi RT, Lynch JB, Del Rio C. Mild or Moderate Covid-19. N Engl J Med (2020) 383(18):1757–66. doi: 10.1056/NEJMcp2009249
17. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: Analysis and Visualization of LD and Haplotype Maps. Bioinformatics (2005) 21(2):263–5. doi: 10.1093/bioinformatics/bth457
18. Zhang J, Wang X, Jia X, Li J, Hu K, Chen G, et al. Risk Factors for Disease Severity, Unimprovement, and Mortality in COVID-19 Patients in Wuhan, China. Clin Microbiol Infect (2020) 26(6):767–72. doi: 10.1016/j.cmi.2020.04.012
19. Pradhan A, Olsson PE. Sex Differences in Severity and Mortality From COVID-19: Are Males More Vulnerable? Biol Sex Differ (2020) 11(1):53. doi: 10.1186/s13293-020-00330-7
20. Onder G, Rezza G, Brusaferro S. Case-Fatality Rate and Characteristics of Patients Dying in Relation to COVID-19 in Italy. JAMA (2020) 323(18):1775–6. doi: 10.1001/jama.2020.4683
21. The Sex, Gender and COVID-19 Project. Available at: https://globalhealth5050.org/the-sex-gender-and-covid-19-project/ (Accessed November 3, 2021).
22. Yang J, Zheng Y, Gou X, Pu K, Chen Z, Guo Q, et al. Prevalence of Comorbidities and Its Effects in Patients Infected With SARS-CoV-2: A Systematic Review and Meta-Analysis. Int J Infect Dis: IJID: Off Publ Int Soc Infect Dis (2020) 94:91–5. doi: 10.1016/j.ijid.2020.03.017
23. Gallo Marin B, Aghagoli G, Lavine K, Yang L, Siff EJ, Chiang SS, et al. Predictors of COVID-19 Severity: A Literature Review. Rev Med Virol (2021) 31(1):1–10. doi: 10.1002/rmv.2146
24. Wu J, Liu J, Zhao X, Liu C, Wang W, Wang D, et al. Clinical Characteristics of Imported Cases of Coronavirus Disease 2019 (COVID-19) in Jiangsu Province: A Multicenter Descriptive Study. Clin Infect Dis (2020) 71(15):706–12. doi: 10.1093/cid/ciaa199
25. Grasselli G, Greco M, Zanella A, Albano G, Antonelli M, Bellani G, et al. Risk Factors Associated With Mortality Among Patients With COVID-19 in Intensive Care Units in Lombardy, Italy. JAMA Intern Med (2020) 180(10):1345–55. doi: 10.1001/jamainternmed.2020.3539
26. Gemmati D, Tisato V. Genetic Hypothesis and Pharmacogenetics Side of Renin-Angiotensin-System in COVID-19. Genes (2020) 11(9):1044. doi: 10.3390/genes11091044
27. Gemmati D, Bramanti B, Serino ML, Secchiero P, Zauli G, Tisato V. COVID-19 and Individual Genetic Susceptibility/Receptivity: Role of ACE1/ACE2 Genes, Immunity, Inflammation and Coagulation. Might the Double X-Chromosome in Females Be Protective Against SARS-CoV-2 Compared to the Single X-Chromosome in Males? Int J Mol Sci (2020) 21(10):3474. doi: 10.3390/ijms21103474
28. Vargas-Alarcón G, Hernández-Pacheco G, Rodríguez-Pérez JM, Pérez-Hernández N, Pavón Z, Fragoso JM, et al. Angiotensin-Converting Enzyme Gene (ACE) Insertion/Deletion Polymorphism in Mexican Populations. Hum Biol (2003) 75(6):889–96. doi: 10.1353/hub.2004.0012
29. Sarangarajan R, Winn R, Kiebish MA, Bountra C, Granger E, NRJJoR N, et al. 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 (2021) 8(4):973–80. doi: 10.1007/s40615-020-00853-0
30. Palomo-Pinon S, Gutierrez-Rodriguez ME, Diaz-Flores M, Sanchez-Barrera R, Valladares-Salgado A, Utrera-Barillas D, et al. DD Genotype of Angiotensin-Converting Enzyme in Type 2 Diabetes Mellitus With Renal Disease in Mexican Mestizos. Nephrol (Carlton) (2009) 14(2):235–9. doi: 10.1111/j.1440-1797.2008.01034.x
31. Pabalan N, Tharabenjasin P, Suntornsaratoon P, Jarjanazi H, Muanprasat C. Ethnic and Age-Specific Acute Lung Injury/Acute Respiratory Distress Syndrome Risk Associated With Angiotensin-Converting Enzyme Insertion/Deletion Polymorphisms, Implications for COVID-19: A Meta-Analysis. Infect Genet Evol: J Mol Epidemiol Evol Genet Infect Dis (2021) 88:104682. doi: 10.1016/j.meegid.2020.104682
32. Pati A, Mahto H, Padhi S, Panda AK. ACE Deletion Allele Is Associated With Susceptibility to SARS-CoV-2 Infection and Mortality Rate: An Epidemiological Study in the Asian Population. Clin Chim Acta (2020) 510:455–8. doi: 10.1016/j.cca.2020.08.008
33. 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. doi: 10.1002/jmv.27160
34. Martínez-Rodríguez N, Posadas-Romero C, Villarreal-Molina T, Vallejo M, Del-Valle-Mondragón L, Ramírez-Bello J, et al. Single Nucleotide Polymorphisms of the Angiotensin-Converting Enzyme (ACE) Gene Are Associated With Essential Hypertension and Increased ACE Enzyme Levels in Mexican Individuals. PloS One (2013) 8(5):e65700. doi: 10.1371/journal.pone.0065700
35. Asselta R, Paraboschi EM, Mantovani A, Duga S. ACE2 and TMPRSS2 Variants and Expression as Candidates to Sex and Country Differences in COVID-19 Severity in Italy. Aging (Albany NY) (2020) 12(11):10087–98. doi: 10.18632/aging.103415
36. Shikov AE, Barbitoff YA, Glotov AS, Danilova MM, Tonyan ZN, Nasykhova YA, et al. Analysis of the Spectrum of ACE2 Variation Suggests a Possible Influence of Rare and Common Variants on Susceptibility to COVID-19 and Severity of Outcome. Front Genet (2020) 11:551220. doi: 10.3389/fgene.2020.551220
37. Falahi S, Kenarkoohi A. Sex and Gender Differences in the Outcome of Patients With COVID-19. J Med Virol (2021) 93(1):151–2. doi: 10.1002/jmv.26243
38. Benetti E, Tita R, Spiga O, Ciolfi A, Birolo G, Bruselles A, et al. ACE2 Gene Variants may Underlie Interindividual Variability and Susceptibility to COVID-19 in the Italian Population. Eur J Hum Genet (2020) 28(11):1602–14. doi: 10.1038/s41431-020-0691-z
39. 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. doi: 10.3389/fgene.2020.564741
Keywords: COVID-19, SARS-CoV-2, ACE2, ACE, genetic variants, polymorphism
Citation: Martínez-Gómez LE, Herrera-López B, Martinez-Armenta C, Ortega-Peña S, Camacho-Rea MdC, Suarez-Ahedo C, Vázquez-Cárdenas P, Vargas-Alarcón G, Rojas-Velasco G, Fragoso JM, Vidal-Vázquez P, Ramírez-Hinojosa JP, Rodríguez-Sánchez Y, Barrón-Díaz D, Moreno ML, Martínez-Ruiz FdJ, Zayago-Angeles DM, Mata-Miranda MM, Vázquez-Zapién GJ, Martínez-Cuazitl A, Barajas-Galicia E, Bustamante-Silva L, Zazueta-Arroyo D, Rodríguez-Pérez JM, Hernández-González O, Coronado-Zarco R, Lucas-Tenorio V, Franco-Cendejas R, López-Jácome LE, Vázquez-Juárez RC, Magaña JJ, Cruz-Ramos M, Granados J, Hernández-Doño S, Delgado-Saldivar D, Ramos-Tavera L, Coronado-Zarco I, Guajardo-Salinas G, Muñoz-Valle JF, Pineda C, Martínez-Nava GA and López-Reyes A (2022) ACE and ACE2 Gene Variants Are Associated With Severe Outcomes of COVID-19 in Men. Front. Immunol. 13:812940. doi: 10.3389/fimmu.2022.812940
Received: 10 November 2021; Accepted: 25 January 2022;
Published: 17 February 2022.
Edited by:
Yoshio Koyanagi, Kyoto University, JapanReviewed by:
Elena Ciaglia, University of Salerno, ItalyCarolinne De Sales Marques, Federal University of Alagoas, Brazil
Copyright © 2022 Martínez-Gómez, Herrera-López, Martinez-Armenta, Ortega-Peña, Camacho-Rea, Suarez-Ahedo, Vázquez-Cárdenas, Vargas-Alarcón, Rojas-Velasco, Fragoso, Vidal-Vázquez, Ramírez-Hinojosa, Rodríguez-Sánchez, Barrón-Díaz, Moreno, Martínez-Ruiz, Zayago-Angeles, Mata-Miranda, Vázquez-Zapién, Martínez-Cuazitl, Barajas-Galicia, Bustamante-Silva, Zazueta-Arroyo, Rodríguez-Pérez, Hernández-González, Coronado-Zarco, Lucas-Tenorio, Franco-Cendejas, López-Jácome, Vázquez-Juárez, Magaña, Cruz-Ramos, Granados, Hernández-Doño, Delgado-Saldivar, Ramos-Tavera, Coronado-Zarco, Guajardo-Salinas, Muñoz-Valle, Pineda, Martínez-Nava and López-Reyes. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alberto López-Reyes, YWxsb3JleUB5YWhvby5jb20=; Gabriela Angélica Martínez-Nava, YW1lcmlhLmp1c3RpY2VAZ21haWwuY29t