Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Combinatorial analysis of ACE and ACE2 polymorphisms reveals protection against COVID-19 worsening: A genetic association study in Brazilian patients

  • Romes Bittencourt Nogueira de Sousa,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Writing – original draft, Writing – review & editing

    Affiliation Department of Physiological Sciences, Institute of Biological Sciences, Federal University of Goiás (UFG), Goiânia, Goiás, Brazil

  • Lis Raquel Silva do Nascimento,

    Roles Data curation, Investigation, Methodology, Validation

    Affiliation Department of Physiological Sciences, Institute of Biological Sciences, Federal University of Goiás (UFG), Goiânia, Goiás, Brazil

  • Luiz Henrique Alves Costa,

    Roles Data curation, Formal analysis, Methodology, Writing – original draft

    Affiliation Department of Physiological Sciences, Institute of Biological Sciences, Federal University of Goiás (UFG), Goiânia, Goiás, Brazil

  • Vanessa Rafaela Milhomem Cruz Leite,

    Roles Data curation, Methodology

    Affiliation Laboratory of Molecular Biology, Institute of Biological Sciences, Federal University of Goiás (UFG), Goiânia, Goiás, Brazil

  • Clayton Luiz Borges,

    Roles Data curation, Methodology

    Affiliation Laboratory of Molecular Biology, Institute of Biological Sciences, Federal University of Goiás (UFG), Goiânia, Goiás, Brazil

  • José Miguel de Deus,

    Roles Conceptualization, Data curation, Investigation, Methodology, Supervision, Validation

    Affiliation Department of Gynecology and Obstetrics, Federal University of Goiás, Goiânia, Goiás, Brazil

  • Ana Cristina Silva Rebelo,

    Roles Project administration, Supervision

    Affiliation Department of Morphology, Institute of Biological Sciences, Federal University of Goiás (UFG), Goiânia, Goiás, Brazil

  • Denise da Silva Pinheiro ,

    Contributed equally to this work with: Denise da Silva Pinheiro, Gustavo Rodrigues Pedrino

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    denisepinheirobiomed@hotmail.com

    Affiliation Laboratory of Clinical Analysis and Health Education, Institute of Biological Sciences, Federal University of Goiás (UFG), Goiânia, Goiás, Brazil

  • Gustavo Rodrigues Pedrino

    Contributed equally to this work with: Denise da Silva Pinheiro, Gustavo Rodrigues Pedrino

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing

    Affiliation Department of Physiological Sciences, Institute of Biological Sciences, Federal University of Goiás (UFG), Goiânia, Goiás, Brazil

Abstract

Since angiotensin-converting enzyme 2, ACE2, was identified as the receptor for SARS-CoV-2 and considering the intense physiological interplay between the two angitensinases isoforms, ACE and ACE2, as counter-regulatory axis of the renin-angiotensin system, we proposed the evaluation of polymorphisms in these two key regulators in relation to COVID-19 severity. A genetic association study involving 621 COVID-19 hospitalized patients from Brazil was performed. All subjects had a confirmed diagnosis of COVID-19 via RT-PCR. Patients were categorized into two groups: the "mild" group (N = 296), composed of individuals hospitalized in ward beds who progressed to cure, and the "severe" group (N = 325), composed of individuals who required hospitalization in an intensive care unit (ICU), or who died. Blood samples were genotyped for ACE I/D polymorphism and ACE2 G8790A polymorphism by real-time PCR via TaqMan assay. The analysis of combined polymorphisms revealed a protective role for genotypic profile II/A_ (ORA = 0,26; p = 0,037) against the worsening of COVID-19 in women. The results indicate a protection profile to COVID-19 progression, in which the II/A_ carriers have almost four times less chance of a severe outcome. It is proposed that a decreased activity of ACE (deleterious effects) in conjunction with an increased ACE2 activity (protective effects), should be the underlying mechanism. The findings are unprecedented once other studies have not explored the genotypic combination analysis for ACE and ACE2 polymorphisms and bring perspectives and expectations for dealing with the COVID-19 pandemic based on definitions of genetically-based risk groups within the context of personalized medicine.

Introduction

Coronavirus 19 Disease (COVID-19) is a condition that has been associated with pathophysiological impairment of the Renin-Angiotensin System (RAS) as a result of infection with the viral agent SARS-CoV-2 [14]. Virus entry into the intracellular environment involves the sequestration of Angiotensin-Converting Enzyme 2 (ACE2), an angitensinase identified as the receptor of the new coronavirus, in the membrane of target cells [2, 57]. This process, which has made ACE2 one of the most studied proteins in recent times [8], can directly lead to RAS imbalance affecting several RAS components, which may be the pivotal link of the new coronavirus with acute lung injury and multiorgan failure [912]. The key role played by angitensinases, ACE2 and its homolog ACE, in the regulation of the RAS, has led researchers to launch a hypothesis regarding genetic susceptibility to COVID-19 involving polymorphisms that alter the activity and/or expression of these enzymes, which has been confirmed in some studies.

With the invasion of the coronavirus, ACE2 undergoes downregulation, a process involving viral attachment, internalization together with viral particles, and subsequent degradation, resulting in ACE2 depletion from the membrane of the host cells [68, 13]. In addition to the internalization of ACE2 upon viral attachment, activation of ADAM17 (a disintegrin and metalloproteinase domain-containing protein 17), also contributes to the decline in ACE2 expression on cell membranes, with the generation of circulating soluble ACE2 [14]. In this sense, the rise in Ang II in COVID infection reflects its reduced biotransformation to Ang 1–7 due to declining ACE2 expression on endothelial cell membranes, which triggers acute deterioration of the lung and other tissues due to alterations such as inflammation, oxidative stress, fibrosis, and thrombosis [5, 10, 11, 14, 15]. In this way, this may be one of the main mechanisms involved in the pathophysiology of coronaviruses, particularly SARS-CoV-2, which still demands better scientific elucidation [14, 16].

Therefore, physiological profiles of ACE and ACE2 must be associated with the outcome of patients affected by COVID-19 by the impairment of angiotensin II to angiotensin 1–7 balance [7, 12]. Corroborating this line of evidence, patients with previous decreased ACE2 activity, such as in the elderly, systemic arterial hypertension (SAH), diabetes mellitus (DM), and cardiovascular diseases (CVDs), are more subject to longer hospitalizations and death from COVID-19, that is, tend to worsen the disease, and such conditions are considered risk factors [1, 3, 1719].

The literature has demonstrated the association between polymorphic variations in the ACE (chromosome 17q23.3) and ACE2 (chromosome Xp22) genes and susceptibility to cardiometabolic diseases such as SAH, DM, and CVDs [2022]. In this context, ACE I/D (rs2106809) and ACE2 G8790A (rs2285666) polymorphisms can be highlighted as the most studied due to their effect on the functionality of the enzymes by impacting genetic mechanisms of alternative splicing (ACE2 G8790A polymorphism in intron 3) or even for presenting linkage disequilibrium with several other genetic polymorphisms (ACE I/D polymorphism in intron 16) [2329]. In this sense, it is reasonable to infer that those polymorphic changes could also affect the prognosis of patients with COVID-19, with the D allele of ACE, associated with increased enzyme activity, and the ACE2 G allele, associated with decreased enzyme activity, prone to deleterious effects in the organisms, as previously shown in a study of our research group related to hypertension susceptibility [22].

In this regard, considering the recent literature, Srivastava et al. (2020) [30] demonstrated the correlation of the A allele of ACE2 G8790A polymorphism with a lower rate of infection and lethality (r = −0.699; p = 0.002) in an ecological study involving Indian populations affected by the disease. Möhlendick et al., (2021) [31] have obtained an association between the G allele of ACE2 G8790A polymorphism and the risk for SARS-CoV-2 infection and the clinical course of COVID-19 in a case-control study conducted in the German population, while the ACE I/D polymorphism had no association. Karakaş Çelik et al. (2021) [32] did not find an association between both polymorphisms and severity in patients affected by COVID-19 in a case study performed in the Turkish population.

The divergent results among the studies cited imply the need for further research involving different populations, and it is noteworthy that the genotypic combination analysis of the two RAS components was not explored, as proposed in the present study, to achieve the provable interplay between them in disease outcome. The current study sought to assess the impact of ACE and ACE2 polymorphisms on the prognosis of COVID-19 patients with a confirmed diagnosis of SARS-CoV-2 infection being treated in a field hospital in Goiânia, GO, Brazil. The results indicated, for the first time in the literature, a genetic protection factor for the susceptibility to the worsening of COVID-19, which was mediated by the combined polymorphisms of ACE and ACE2. We should also highlight that this is the first study to investigate a genetic association in the outcome of COVID-19 carried out in the Brazilian population.

Material and methods

Ethical statement

This study was approved by the Ethics Committee of the Federal University of Goiás (4140331/2020). It was conducted in accordance with the ethical principles for medical research involving human subjects of the World Medical Association Declaration of Helsinki, and, after explaining the research, the volunteers or their families signed a free and informed consent form.

Subjects

A total of 621 individuals (348 men and 273 women, age = 61.45±15.90 years) were recruited during the period from August to December 2020, while SARS-CoV-2 were the prevalent circulating serotype of Coronavirus in Brazil and before the beginning of vaccination [33]. Patients were subdivided into the "moderate group" (n = 296, 161 men, 135 women, age = 57.49±15.78 years) and the "severe group" (n = 325, 187 men, 138 women, age = 65, 06±15.20 years). All subjects were diagnosed with COVID-19 by RT-PCR. The moderate group was composed of patients who received health care at the Municipal Field Hospital of Goiânia-Goiás, and who were admitted to a hospital bed in a ward, having evolved to cure. The severe group consisted of those who required admission to the Intensive Care Unit (ICU) in the same hospital, as previously informed, including individuals who died while hospitalized in ICU beds.

Sociodemographic and clinical data (gender, age, and comorbidities) were collected by consulting medical records, in which the daily evolution of each patient was sought to allow discrimination between the research groups. Comorbidities were grouped into the following subgroups, according to the etiology of the disease or disorder (only data of preceding diseases were collected): hypertension, diabetes, respiratory system diseases (RSDs), cardiovascular diseases (CVDs), obesity, neuropsychiatric diseases, and disorders (NDDs) and acute or chronic kidney diseases (ACKD). Data on smoking was also collected. All these data were confirmed and considered in the statistical analyses for comparison between groups in the regression modeling.

Peripheral blood samples were collected through venipuncture in 4 ml ETDA tubes and centrifuged at 5000 rpm for 20 minutes to allow separation of the leukocyte ring. This was collected by Pasteur pipettes and stored at -20°C in 1.5 ml microtubes, aiming at DNA extraction.

Genotyping of ACE and ACE2 polymorphisms

DNA extraction was performed using the DNA PureLink™Genomic DNA Mini Kit, according to the manufacturer’s suggestions. The extraction products were labeled and stored at -80o C for use in the patronized genotyping reaction. Molecular analysis of genetic polymorphisms was performed by a real-time polymerase chain reaction (qPCR) TaqMan assay.

In the genotyping of ACE I/D polymorphism, it was adopted the protocol developed by Koch et al. [34] adapted for use in the QuantStudio 5 (Thermofiher®) apparatus. Primers sequences were: ACE111: 5’-CCCATCCTTTCTCCCATTTCTC-3’; ACE112: 5’-AGCTGGAATAAAATTGGCGAAAC-3’; ACE113: 5’-CCTCCCAAAGTGCTGGGATTA-3’. Probes sequences were: VIC-5’-AGGCGTGATACAGTCA-3’ (I allele); FAM-5’-TGCTGCCTATACAGTCA-3’ (D allele).

Thermocycling conditions were: Pre-Hold—60°C for 30 seconds; Hold-Stage—95°C for 10 minutes; PCR-Stage (40 cycles) - 95°C for 15 seconds and 60°C for 1 minute; Post-Read Stage—60°C for 30 seconds.

The discrimination of genotypes II, ID, and DD was performed by analyzing the amplification curves at the end of thermocycling. It was considered: Allele 1- VIC–Insertion (ACE111/ACE113-ACE112); Allele 2- FAM–Deletion (AC111-AC112).

For the genotyping of ACE2 G8790A SNP polymorphism, a previously developed and validated TaqMan SNP allelic discrimination assay was adopted (ThermofisherScientific®, Assay ID: C__2551626_1; SNP ID: rs2285666) for use in the QuantStudio 5 device. The thermocycling conditions were the same as previously described for ACE I/D. Allelic discrimination, as well as the sequence of probes and primers, followed the protocol established by Thermofisher®. The notation A_ is used to refer to AA and GA genotypes, while G_ is used to designate GG and GA genotypes.

Statistical analysis

Sociodemographic, clinical, and genotyping data were tabulated in Microsoft Office Excel 2007 software and analyzed with BioEstat v. 5.3 (Mamiraua Institute, available at https://mamiraua.org.br) and IBM SPSS v. 22. For comparison between groups, continuous variables were expressed as mean and standard deviation and analyzed using the T-test, while qualitative variables, expressed as yes/no, were analyzed using the χ2 test or Fisher’s exact test when necessary. To analyze the susceptibility to COVID-19 outcomes (moderate and severe groups), the Odds Ratio (OR) test was performed, with a 95% confidence interval, and adjusted Odds Ratio (ORA) values were obtained by binary multifactorial logistic regression analysis, with control of confounding factors. The confounding factors used were age, hypertension, diabetes, respiratory system diseases, cardiovascular diseases, obesity, neuropsychiatric diseases and disorders, acute or chronic kidney diseases and smoking (Table 1). In relation to the genotypes, the fidelity to the Hardy-Weinberg Equilibrium was measured by the χ2 test. The p-value considered indicative of significance was p<0.05. As ACE2 gene is located on the X chromosome, the multinomial logistic regression was performed separately for the two sexes where it appears.

thumbnail
Table 1. Variables in the multifactorial logistic regression analysis for ACE I/D and ACE2 G8790A polymorphisms in the moderate and severe COVID-19 groups.

https://doi.org/10.1371/journal.pone.0288178.t001

Results

Clinical and sociodemographic variables for the studied groups are shown in Table 2. The distribution by sex, obesity and smoking variables did not show disagreement between the groups, pointing to homogeneity for the subsequent analysis. As expected, the variable age had increased values in the severe outcome group (p<0.001), once it was a considerable risk factor for the worsening of patients affected by COVID-19. The same occurred for all the clinical conditions and diseases evaluated, which were significantly more frequent in the severe group, corroborating the well-known role of comorbidities in disease aggravation.

thumbnail
Table 2. Characterization of clinical variables for the moderate and severe groups.

https://doi.org/10.1371/journal.pone.0288178.t002

Regarding the genotypic analysis (Table 3), in the moderate group, the frequencies were: ACE I/D: II: 19.6%; ID: 48.0%, DD: 32.4%. For the ACE2 G8790A polymorphism, in females, the frequencies were: 50.4% for GG, 38.5% for GA, and 11.1% for AA. The allelic frequencies for males were 76.4% (G) and 23.6% (A). In the severe group, the following values were observed for ACE I/D polymorphism: 20.3% for II; 53.5% for ID, and 26.8% for DD. For the ACE2 G8790A polymorphism, in females, the frequencies were: 58.7% for GG, 34.8% for GA, and 6.5% for AA. Allelic frequencies for males were: 78.1% (G) and 21.9% (A). No statistically significant differences were observed in the frequency distribution between groups (analyses by the chi-square test), and no association was verified in the susceptibility analysis for each genetic polymorphism individually.

thumbnail
Table 3. Genotypic frequency distribution of ACE I/D and ACE2 G8790A polymorphisms in the moderate and severe COVID-19 groups and risk analysis for the outcome in COVID-19.

https://doi.org/10.1371/journal.pone.0288178.t003

In the association analysis for the combined genotypes of the two genes (Table 4), there was a predominance of occurrence of the ACE/ACE2 genotypic profile II/A_ (p = 0.028) in the moderate group for females (analysis by the dominant model of ACE2 –grouping GA and AA individuals). The result was confirmed in the susceptibility analysis after adjustment by logistic regression, which revealed that it is a protective genotype against the worsening of COVID-19 (ORA = 0.26; p = 0.037) conferring almost four times less chance of having a worse prognosis in the evolution of the disease. For the other genotypic combination models, no statistically significant differences were observed.

thumbnail
Table 4. Genotypic frequency distribution for the combination of ACE I/D and ACE2 G8790A polymorphisms in moderate and severe COVID-19 groups and risk analysis for the outcome in COVID-19.

https://doi.org/10.1371/journal.pone.0288178.t004

The distribution of ACE and ACE2 genotypes (Table 5) for the sampled population was consistent with that expected by the Hardy-Weinberg Equilibrium, suggesting the absence of selective pressures involving the genetic locus evaluated and a representative sampling of the studied population in the study field.

thumbnail
Table 5. Hardy-Weinberg equilibrium test for ACE I/D and ACE2 G8790A polymorphisms in moderate and severe groups of COVID-19.

https://doi.org/10.1371/journal.pone.0288178.t005

Discussion

The possibility of a genetic predisposition to SARS-CoV-2 infection has been raised in the scientific literature, especially considering the role that polymorphisms of ACE and ACE2 could play in the infection process and disease progression, considering the activity of these enzymes in the pathophysiology of acute lung injury and in the invasion mechanism used by the new coronavirus [3538]. The results of the present study point to the existence, already known in the literature, of health conditions as risk factors, as well as the unprecedented existence of a genetic protection factor regarding the susceptibility to worsening of the patient with COVID-19 mediated by the combined polymorphisms of ACE and ACE2.

As observed, the elderly (p<0.001) and the presence of comorbidities such as systemic arterial hypertension (p = 0.0310), diabetes mellitus (p = 0.0016), respiratory diseases (p = 0.0018), cardiovascular diseases (p = 0.0003), neuropsychiatric diseases and disorders (p = 0.0001) and acute and chronic kidney diseases (p = 0.0002) were significantly more frequent in individuals with more severe outcomes compared to less severe ones when infected with SARS-CoV-2, being associated with the worsening of the disease as independent risk factors. These associations are already well recognized in the literature for the Brazilian and worldwide populations [1, 39, 40], as previously expected, due to the systemic pathophysiological changes caused by COVID-19, especially in serious situations [41].

The study of the role played by key enzymes that make up RAS is essential for understanding the endogenous contributors to the worsening of COVID-19. The role of increased ACE2 expression in the complications of the disease is debated because that would facilitate virus entry into the cells, but at the same time, expressing less, or having less activity of this enzyme would also contribute to other aggravating conditions imposed by the virus that is involved in lung injury and multiorgan failure [3, 6, 7, 12, 4244]. In this sense, the study of the impact of polymorphisms, epigenetic changes, and other molecular mechanisms in the modulation of ACE and ACE2 activity or expression is important to elucidate the worsening of the disease [16]. ACE I/D (rs2106809) and ACE2 G8790A (rs2285666) polymorphisms are studied in a great variety of conditions by their effect on the plasma and tissue levels of these key RAS enzymes, having been described as functional polymorphisms [45]. It´s well known the ACE D allele relation with increased ACE activity [24, 25, 45], while evidence points to the ACE2 A allele association with enhanced gene expression and the possibility of interference in the protein product by alternative splicing [30, 46, 47]. The molecular mechanisms underlying these effects still need to be better elucidated.

In the results of this study, the genotypic profiles of ACE and ACE2, individually, did not represent susceptibility or protection for the patient affected by COVID-19. However, it was found that the genotypic combination of the II genotype of ACE with the A allele of ACE2 (II/A_) was significantly more frequent in female individuals from the mild outcome group of COVID-19 (χ2 = 4.83; p = 0.028), presenting a protective role against the worsening of the disease (ORA = 0.26; p = 0.037).

The association results obtained mean that II/A_ carriers have almost four times (1 over 0.26) less chance of a severe outcome in females. This result corroborates gender-related differences in COVID-19 susceptibility in the sense of a male increased mortality rate, leading to reflections concerning a possible differential ACE2 expression in woman involving mechanisms of X chromosome inactivation escape, epigenetic modifications, and sex hormones regulation [48]. Epigenetic modifications during aging, marked by reduced DNA methylation and reversion of inactivated genes, would also reflect changes in the incidence and prognosis of the disease [49].

This unprecedented result is in line with Hatami et al. 2020 [50], who have pointed to a possible association between the increase in I/D ratio and the increase in COVID-19 recovery rate by meta-regression analysis involving prevalence data from 30 different countries (point estimate: 0.48, CI 95%: 0.05–0.91, p  =  0.027). Then, the presence of the I allele and, mainly, the II genotype of ACE may be considered a benefical factor in disease progression.

In previous studies, the D allele or DD genotype of ACE has been associated with increased enzyme activity and implied susceptibility to diseases, such as hypertension and CVDs, in which there would be greater activity of the enzyme [2022]. Thus, the I allele would be associated with lower ACE activity, resulting in lower Ang-II production, and, therefore, it would confer a more favorable physiological status to the success in coping with the pathology, leading to a lower pro-inflammatory cascade during the infection course, marked by ACE2 downregulation and, consequently, RAS imbalance.

In this line of evidence, Sieńko et al. (2020) [51], in a narrative review, pointed to the importance of the D allele in the risk of serious outcomes in COVID-19. Delanghe et al. (2020) [35] draw attention to an important geographical variation of ACEI/D polymorphism that began to be associated with COVID-19 prevalence and mortality data [36]. Then, an ecological study conducted by Yamamoto et al. (2020) [52] involving data from 26 countries in Europe and Asia reported a negative ACE II genotype correlation with the number of COVID-19 cases (R = − 0.847) and deaths (R = − 0.755).

Yamamoto et al. (2020) [52] highlight that ethnic-geographical influences can impact the findings of the studies, given that populations with different genetic backgrounds can have different outcomes. In this sense, the present work, despite not having found a profile of genetic susceptibility, is in line with the results indicated, once it was verified an association of protection for the II genotype of ACE when in combination with the A allele of ACE2, having to take into account that it sampled individuals from the Brazilian population, which is recognizably marked by great racial miscegenation in its origins and in its historical course.

In agreement with the findings of the present study, Manning & Fink (2020) [53] in an ecological study involving data on COVID-19 fatalities among 41 nations suggested that genotype II could have a protective role against worsening. The authors speak in favor of a supposed “oxygen mechanism”, according to which there would be a polymorphic contribution of the D allele to lung damage and reduced blood O2 saturation. In our study, it was not possible to verify an association of the D allele with aggravation in patients affected by the disease.

Similar results were verified by Aung et al. (2020) [54] in another ecological study contemplating 25 countries from diverse geographical regions of the world. The study pointed out that no association was found between the frequency of the DD genotype and COVID-19 mortality rates (OR = 4.3; p = 0.2). Despite this, the researchers verified that increased II genotype frequency was significantly associated with decreased COVID-19 mortality rates (IRR: 0.3; p = 0.03).

Regarding the G8790A polymorphism of the ACE2 gene, in an ecological study involving Indian populations, Srivastava et al. (2020) [30] pointed to a correlation of the A allele with a lower rate of infection (r = −0.571; p = 0.021) and lethality (r = −0.699; p = 0.002) in patients affected by the disease. The result agrees with our work in the sense that a protective factor was observed for the A allele in patients already affected by COVID-19.

Otherwise, Möhlendick et al., (2021) [31] have identified the ACE2 G allele and GG genotype as a genetic risk factor for SARS-CoV-2 infection and COVID-19 mortality, evaluating 297 SARS-CoV-2-positive and 253 SARS-CoV-2-negative tested patients in the German population. A lack of association for the ACE I/D polymorphism was observed. Although the cited association was not verified in our study, it corroborates the hypothesis of the functional aspects associated with ACE2 polymorphisms, especially that the G allele is associated with greater pulmonary and systemic injuries due to the imbalance in the RAS with the decrease in ACE2 activity.

In this study, we did not aim at estimating the influence of genetic polymorphisms on the prevalence of infection by the disease, but rather at associating them with the type of clinical outcome during hospitalization. However, regarding the allelic frequencies, it can be observed that the frequencies obtained in this study are close to those found by Pinheiro et al. (2019) [22], evaluating 123 controls and 117 hypertensive patients in the Brazilian population, considering the control group, for both polymorphisms evaluated in severe and moderate groups, with ACE I/D frequencies around 0.5 for each allele and, for ACE2, a proportion about 0.2 for A and 0.8 for G alleles. Thus, it is likely that there would be no change in the rate of infection by COVID-19 associated with these polymorphisms for the Brazilian population.

The initial results of non-association for ACE and ACE2 polymorphisms analyzed separately are in agreement with a case study conducted by Karakaş Çelik et al. (2021) [32] that evaluated 155 patients divided into three COVID-19 groups (mild, moderate, and severe); however, the referred research has not explored the genotypic combination analyses, which have brought interesting association findings to this work and have been an important tool to be employed on genetic association studies [22].

The synergistic effects of ACE and ACE2 are in accordance with their counterbalancing biological roles in RAS and constitute a finding that has already been highlighted in another study of our research group [22], which revealed an important association of the DD/G_ profile with hypertension susceptibility (ORA = 3.57; p = 0.026), comorbidity, which is also a COVID-19 risk factor. Interestingly, the results obtained in this study indicate protection for the complementary combination profile (II/A_), which can be considered along the same line of understanding, since if we have susceptibility on one side, we can expect protection on the other.

Similarly, if in that study we have pointed out that the increased activity of ACE (DD genotype) in conjugation with decreased activity of ACE2 (G allele) could be the underlining mechanism in hypertension susceptibility, here we hypothesize that the associated opposite effects on enzyme activity in the complementary genotype profile, that is, decreased ACE activity in conjunction with the increased ACE2 activity in patients with the protective genotype (II/A_), should be the protective physiological mechanism involved in COVID-19 outcome.

It’s interesting to highlight epigenetic mechanisms, which include DNA methylation, histone post-translational modifications, and microRNAs, that can impact ACE2 expression, have been associated with conditions of disease [49]. When considering the findings of this study, the epigenetic inactivation of the A allele (protective) in heterozygous women (GA) would have a prejudicial effect, mainly on ACE DD individuals. Through aging, the reversion of epigenetic silencing in female would have a positive impact on heterozygotes due to the possibility of A allele reactivation. This would help explain aspects of the differential sex-related susceptibility to COVID-19, once male subjects are hemizygous for ACE2 gene, and, therefore, cannot be beneficed by these epigenetic changes.

Considering the survey of the current literature, attention is drawn to the scarcity of studies addressing the G8790A polymorphism of ACE2 in the ontogeny of COVID-19 infection and entry of the virus into host cells, while most of the studies have been dedicated to exploring genetic variations in the coding region of the ACE2 gene [55, 56] with important results being obtained regarding variants that interfere with the virus-receptor interaction [16, 5557]. However, it is worth mentioning that polymorphisms that influence enzymatic activity, even in an intronic region, can bring interesting findings because they have the potential to impact the clinical course of COVID-19 by impairing RAS balance, as verified in this study and evidenced by the results of Khayat et al. (2020) [47].

So far, no other work has presented the results found in this study for the genetic association between ACE and ACE2 polymorphisms. The non-association for each polymorphism alone supports the hypothesis of genetic complexity behind outcomes for COVID-19, which explains the difficulty in the search for completely effective treatments to combat severe outcomes for the disease. It is reasonable to assume that a greater understanding of polymorphisms associated with the worsening of COVID-19 is highly relevant in decision-making in public policies, whether, for example, in the development of vaccination strategies by genotypic risk groups, or in personalized care for each risk profile.

As limitations and complements to this study, we highlight the expansion of the sample, both in number and in population, as well as the evaluation of the concentration of detached circulating enzymes in the plasma. The corroboration of these results by further research can substantially benefit personalized medicine, which favors the adoption of measures linked to evidence-based medicine given the need to continue reducing mortality and controlling the disease.

Conclusions

The study of genetic risk factors is one arm of personalized medicine development. The results obtained in this study indicate a protection profile for COVID-19 progression for the population under analysis, in which the II/A_ carriers have almost four times less chance to evolve to a severe outcome. This means that, not having this genetic combination can be a signal of poor prognosis for the patient.

The finding is unprecedented in the literature, as other studies have not explored the genotypic combination analysis for ACE and ACE2 polymorphisms. In line with the literature, it is reasonable to assume that the protective alleles identified are associated with greater homeostasis in pulmonary RAS, reducing the pro-inflammatory profile generated by the disease, and thus preventing worse outcomes.

The number of independent risk factors, including genetic contributors, associated with aggravation or protection against COVID-19 identified in this and other research bring to light the multifactorial nature of the disease, with the study of new genes being highly suggestible for a better understanding of COVID-19 pathophysiology and the development of new treatment strategies.

Supporting information

S1 File. Data from mild and severe COVID-19 groups.

https://doi.org/10.1371/journal.pone.0288178.s001

(DOCX)

Acknowledgments

The authors thank the Laboratory of Molecular Biology and the Laboratory of Clinical Analysis and Health Education, both at the Federal University of Goiás, for their great contributions in carrying out the experiments and storing samples for this study. Special acknowledgments are due to the Municipal Field Hospital of Goiânia/GO for the team’s availability in obtaining samples and clinical data. The authors thank the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) for supporting this work.

References

  1. 1. 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. The Lancet. 2020;395: 497–506. pmid:31986264
  2. 2. Lu R, Zhao X, Li J, Niu P, Yang B, Wu H, et al. Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. The Lancet. 2020;395: 565–574. pmid:32007145
  3. 3. Kumar M, al Khodor S. Pathophysiology and treatment strategies for COVID-19. J Transl Med. 2020;18: 1–9.
  4. 4. Kirtipal N, Bharadwaj S, Kang SG. From SARS to SARS-CoV-2, insights on structure, pathogenicity and immunity aspects of pandemic human coronaviruses. Infection, Genetics and Evolution. 2020;85: 104502. pmid:32798769
  5. 5. Wan Y, Shang J, Graham R, Baric RS, Li F. Receptor recognition by novel coronavirus from Wuhan: An analysis based on decade-long structural studies of SARS. J Virol. 2020;94. pmid:31996437
  6. 6. Bastolla U, Chambers P, Abia D, Garcia-Bermejo M-L, Fresno M. Is Covid-19 Severity Associated With ACE2 Degradation? Frontiers in Drug Discovery. 2022;1: 5.
  7. 7. Silhol F, Sarlon G, Deharo JC, Vaïsse B. Downregulation of ACE2 induces overstimulation of the renin–angiotensin system in COVID-19: should we block the renin–angiotensin system? Hypertension Research 2020 43:8. 2020;43: 854–856. pmid:32439915
  8. 8. Bader M, Turner AJ, Alenina N. ACE2, a multifunctional protein–from cardiovascular regulation to COVID-19. Clin Sci. 2020;134: 3229–3232. pmid:33305790
  9. 9. Hamming I, Timens W, Bulthuis MLC, Lely AT, Navis GJ, van Goor H. Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis. Journal of Pathology. 2004;203: 631–637. pmid:15141377
  10. 10. Imai Y, Kuba K, Penninger JM. [Lessons from SARS: a new potential therapy for acute respiratory distress syndrome (ARDS) with angiotensin converting enzyme 2 (ACE2)]. Masui. 2008;57: 302–10. Available: http://www.ncbi.nlm.nih.gov/pubmed/18340998 pmid:18340998
  11. 11. Imai Y, Kuba K, Penninger JM. The discovery of angiotensin-converting enzyme 2 and its role in acute lung injury in mice. Experimental Physiology. Blackwell Publishing Ltd; 2008. pp. 543–548. pmid:18448662
  12. 12. Verdecchia P, Cavallini C, Spanevello A, Angeli F. The pivotal link between ACE2 deficiency and SARS-CoV-2 infection. Eur J Intern Med. 2020;76: 14–20. pmid:32336612
  13. 13. Heyman SN, Kinaneh S, Abassi Z. The Duplicitous Nature of ACE2 in COVID-19 Disease. EBioMedicine. 2021;67. pmid:33910120
  14. 14. Abassi Z, Higazi AAR, Kinaneh S, Armaly Z, Skorecki K, Heyman SN. ACE2, COVID-19 Infection, Inflammation, and Coagulopathy: Missing Pieces in the Puzzle. Front Physiol. 2020;11: 1253. pmid:33123031
  15. 15. Kuba K, Imai Y, Penninger JM. Angiotensin-converting enzyme 2 in lung diseases. Current Opinion in Pharmacology. Elsevier; 2006. pp. 271–276. pmid:16581295
  16. 16. Lubbe L, Cozier GE, Oosthuizen D, Ravi Acharya K, Sturrock ED. ACE2 and ACE: structure-based insights into mechanism, regulation and receptor recognition by SARS-CoV. Clin Sci. 2020;134: 2851–2871. pmid:33146371
  17. 17. Guan W, Ni Z, Hu Y, Liang W, Ou C, He J, et al. Clinical Characteristics of Coronavirus Disease 2019 in China. New England Journal of Medicine. 2020; NEJMoa2002032. pmid:32109013
  18. 18. 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. Journal of Allergy and Clinical Immunology. 2020;146: 110–118. pmid:32294485
  19. 19. Rosenthal N, Cao Z, Gundrum J, Sianis J, Safo S. Risk Factors Associated With In-Hospital Mortality in a US National Sample of Patients With COVID-19. JAMA Netw Open. 2020;3: e2029058–e2029058. pmid:33301018
  20. 20. Nawaz SK, Hasnain S. Pleiotropic effects of ACE polymorphism. Biochem Med (Zagreb). 2009; 36–49.
  21. 21. Lambert DW, Clarke NE, Turner AJ. Not just angiotensinases: New roles for the angiotensin-converting enzymes. Cellular and Molecular Life Sciences. 2010. pmid:19763395
  22. 22. Pinheiro DS, Santos RS, Jardim PCBV, 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: e0221248. pmid:31430320
  23. 23. Rigat B, Hubert C, Corvol P, Soubrier F. PCR detection of the insertion/deletion polymorphism of the human angiotensin converting enzyme gene (DCP1) (dipeptidyl carboxypeptidase 1). Nucleic Acids Res. 1992;20: 1433. Available: http://www.ncbi.nlm.nih.gov/pubmed/1313972 pmid:1313972
  24. 24. Tiret L, Rigat B, Visvikis S, Breda C, Corvol P, Cambien F, et al. Evidence, from combined segregation and linkage analysis, that a variant of the angiotensin I-converting enzyme (ACE) gene controls plasma ACE levels. Am J Hum Genet. 1992;51: 197–205. Available: http://www.ncbi.nlm.nih.gov/pubmed/1319114 pmid:1319114
  25. 25. Rigat B, Hubert C, Alhenc-Gelas F, Cambien F, Corvol P, Soubrier F. An insertion/deletion polymorphism in the angiotensin I-converting enzyme gene accounting for half the variance of serum enzyme levels. J Clin Invest. 1990;86: 1343–6. pmid:1976655
  26. 26. Zhu X, Bouzekri N, Southam L, Cooper RS, Adeyemo A, McKenzie CA, et al. Linkage and Association Analysis of Angiotensin I–Converting Enzyme (ACE)–Gene Polymorphisms with ACE Concentration and Blood Pressure. The American Journal of Human Genetics. 2001;68: 1139–1148. pmid:11283791
  27. 27. Patel SK, Velkoska E, Freeman M, Wai B, Lancefield TF, Burrell LM. From gene to protein-experimental and clinical studies of ACE2 in blood pressure control and arterial hypertension. Front Physiol. 2014;5 JUN: 1–12. pmid:25009501
  28. 28. Burrell LM, Harrap SB, Velkoska E, Patel SK. The ACE2 gene: its potential as a functional candidate for cardiovascular disease. Clin Sci (Lond). 2013;124: 65–76. pmid:23013041
  29. 29. Rieder MJ, Nickerson DA. Hypertension and Single Nucleotide Polymorphisms. Current Hypertension Reports Current Science Inc ISSN. 2000;2: 44–49. Available: https://link.springer.com/content/pdf/10.1007%2Fs11906-000-0057-4.pdf pmid:10981126
  30. 30. 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: 1163. pmid:33101387
  31. 31. Möhlendick B, Schönfelder 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; 165–171. pmid:34001841
  32. 32. Karakaş Çelik S, Çakmak Genç G, Pişkin N, Açikgöz B, Altinsoy B, Kurucu İşsiz 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: 5947–5952. pmid:34170561
  33. 33. Alcantara LCJ, Nogueira E, Shuab G, Tosta S, Fristch H, Pimentel V, et al. SARS-CoV-2 epidemic in Brazil: how the displacement of variants has driven distinct epidemic waves. Virus Res. 2022;315. pmid:35461905
  34. 34. Koch W, Latz W, Eichinger M, Ganser C, Schömig A, Kastrati A. Genotyping of the Angiotensin I-Converting Enzyme Gene Insertion/Deletion Polymorphism by the TaqMan Method. Clin Chem. 2005;51: 1547–1549. pmid:16040856
  35. 35. Delanghe JR, Speeckaert MM, de Buyzere ML. The host’s angiotensin-converting enzyme polymorphism may explain epidemiological findings in COVID-19 infections. Clinica Chimica Acta. 2020 [cited 5 Apr 2020]. pmid:32220422
  36. 36. Delanghe JR, Speeckaert MM, De Buyzere ML. COVID-19 infections are also affected by human ACE1 D/I polymorphism. Clinical Chemistry and Laboratory Medicine. De Gruyter; 2020. pmid:32286246
  37. 37. Jacobs M, Lahousse L, van Eeckhoutte HP, Wijnant SRA, Delanghe JR, Brusselle GG, et al. Effect of ACE1 polymorphism rs1799752 on protein levels of ACE2, the SARS-CoV-2 entry receptor, in alveolar lung epithelium. ERJ Open Res. 2021;7: 00940–02020. pmid:33889638
  38. 38. Fang L, Karakiulakis G, Roth M. Are patients with hypertension and diabetes mellitus at increased risk for COVID-19 infection? Lancet Respir Med. 2020;8: e21. pmid:32171062
  39. 39. Souza FSH de, Hojo-Souza NS, Batista BD de O, Silva CM da, Guidoni DL. On the analysis of mortality risk factors for hospitalized COVID-19 patients: A data-driven study using the major Brazilian database. PLoS One. 2021;16: e0248580. pmid:33735272
  40. 40. Iftimie S, López-Azcona AF, Vicente-Miralles M, Descarrega-Reina R, Hernández-Aguilera A, Riu F, et al. Risk factors associated with mortality in hospitalized patients with SARS-CoV-2 infection. A prospective, longitudinal, unicenter study in Reus, Spain. Adrish M, editor. PLoS One. 2020;15: e0234452. pmid:32881860
  41. 41. dong Gao Y, Ding M, Dong X, jin Zhang J, Kursat Azkur A, Azkur D, et al. Risk factors for severe and critically ill COVID-19 patients: A review. Allergy. 2021;76: 428–455. pmid:33185910
  42. 42. Hofmann H, Geier M, Marzi A, Krumbiegel M, Peipp M, Fey GH, et al. Susceptibility to SARS coronavirus S protein-driven infection correlates with expression of angiotensin converting enzyme 2 and infection can be blocked by soluble receptor. Biochem Biophys Res Commun. 2004;319: 1216–1221. pmid:15194496
  43. 43. Hoffmann M, Kleine-Weber H, Schroeder S, Krüger N, Herrler T, Erichsen S, et al. SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor. Cell. 2020;181: 271–280.e8. pmid:32142651
  44. 44. Ferrario CM, Ahmad S, Groban L. Twenty years of progress in angiotensin converting enzyme 2 and its link to SARS-CoV-2 disease. Clin Sci. 2020;134: 2645–2664. pmid:33063823
  45. 45. Gintoni I, Adamopoulou M, Yapijakis C. The Impact of ACE and ACE2 Gene Polymorphisms in Pulmonary Diseases Including COVID-19. In Vivo (Brooklyn). 2022;36: 13. pmid:34972696
  46. 46. Wu YH, Li JY, Wang C, Zhang LM, Qiao H. The ACE2 G8790A Polymorphism: Involvement in Type 2 Diabetes Mellitus Combined with Cerebral Stroke. J Clin Lab Anal. 2017. pmid:27500554
  47. 47. Khayat AS, De Assumpção PP, Khayat BCM, Araújo TMT, Batista-Gomes JA, Imbiriba LC, et al. ACE2 polymorphisms as potential players in COVID-19 outcome. PLoS One. 2020;15: e0243887. pmid:33370311
  48. 48. Foresta C, Rocca MS, Di Nisio A. Gender susceptibility to COVID-19: a review of the putative role of sex hormones and X chromosome. J Endocrinol Invest. 2021;44: 951–956. pmid:32936429
  49. 49. Lima RS, Rocha LPC, Moreira PR. Genetic and epigenetic control of ACE2 expression and its possible role in COVID‐19. Cell Biochem Funct. 2021;39: 713. pmid:34075603
  50. 50. Hatami N, Ahi S, Sadeghinikoo A, Foroughian M, Javdani F, Kalani N, et al. Worldwide ACE (I/D) polymorphism may affect COVID-19 recovery rate: an ecological meta-regression. Endocrine. 2020;68: 479–484. pmid:32542429
  51. 51. Sieńko J, Kotowski M, Bogacz A, Lechowicz K, Drożdżal S, Rosik J, et al. COVID-19: The Influence of ACE Genotype and ACE-I and ARBs on the Course of SARS-CoV-2 Infection in Elderly Patients. Clin Interv Aging. 2020;15: 1231–1240. pmid:32764907
  52. 52. Yamamoto N, Ariumi Y, Nishida N, Yamamoto R, Bauer G, Gojobori T, et al. SARS-CoV-2 infections and COVID-19 mortalities strongly correlate with ACE1 I/D genotype. Gene. 2020;758: 144944. pmid:32628976
  53. 53. Manning JT, Fink B. Understanding COVID-19: A hypothesis regarding digit ratio (2D:4D), ACE I/D polymorphism, oxygen metabolism and national case fatality rates. Early Hum Dev. 2020;151: 105161. pmid:32905897
  54. 54. Aung AK, Aitken T, Teh BM, Yu C, Ofori-Asenso R, Chin KL, et al. Angiotensin converting enzyme genotypes and mortality from COVID-19: An ecological study. Journal of Infection. 2020;81: 961–965. pmid:33197472
  55. 55. Suryamohan K, Diwanji D, Stawiski EW, Gupta R, Miersch S, Liu J, et al. Human ACE2 receptor polymorphisms and altered susceptibility to SARS-CoV-2. Communications Biology 2021 4:1. 2021;4: 1–11. pmid:33846513
  56. 56. Li J, Wang Y, Liu Y, Zhang Z, Zhai Y, Dai Y, et al. Polymorphisms and mutations of ACE2 and TMPRSS2 genes are associated with COVID-19: a systematic review. Eur J Med Res. 2022;27: 1–10.
  57. 57. 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: 1129. pmid:33133145