Skip to main content

The Interleukin-6 gene variants may protect against SARS-CoV-2 infection and the severity of COVID-19: a case-control study in a Moroccan population

Abstract

The symptoms of SARS-CoV-2 infection vary widely, ranging from asymptomatic cases to severe forms marked by acute respiratory distress syndrome, multi-organ damage, and fatalities. Studies indicate a correlation between specific genes and susceptibility to SARS-CoV-2 infection and disease severity, particularly involving variants in genes linked to inflammation and immune responses. The objective of this study is to investigate the association between rs1800795 (− 174 G > C) and rs1800797 (− 597 A > G) variants in the interleukin-6 (IL-6) promoter region and susceptibility to SARS-CoV-2 infection. Additionally, we aim to explore their correlation with COVID-19 severity in a Moroccan population. In this case-control study, we enrolled 270 unvaccinated COVID-19 patients, consisting of 132 with severe COVID-19 and 138 with asymptomatic-moderate COVID-19. Additionally, we included 339 SARS-CoV-2-negative group. Genotyping of rs1800795 and rs1800797 polymorphisms of the IL-6 gene was performed using predesigned TaqMan SNP genotyping. The median age of SARS-CoV-2-negative controls was 50 years, while severe COVID-19 cases exhibited a median age of 61 years. Additionally, individuals with asymptomatic to moderate COVID-19 had a median age of 36 years. We observed a significant age difference between severe and mild COVID-19 patients (p < 0.0001), and an association was noted between gender and the severity of COVID-19 (p = 0.011). The allele and genotype frequencies of the IL-6 − 597G > A and − 174G > C variants did not show significant associations with susceptibility to SARS-CoV-2 infection (p > 0.05). However, further analysis revealed that the linkage disequilibrium between rs1800797 and rs1800795 indicated that individuals with the GC* haplotype (OR = 0.04, 95% CI 0.01–0.30, p = 0.001) and AG* haplotype (OR = 0.11, 95% CI 0.03–0.46, p = 0.002) were significantly associated with protection against SARS-CoV-2 infection. Moreover, in the overdominant model, the IL-6 − 174 G/C genotype was found to be protective against the development of severe disease compared to those with the G/G-C/C genotypes (p = 0.03; OR = 0.41, 95% CI 0.18–0.96). However, correlations between complete blood count markers, hematological markers, D-dimer, C-reactive protein, and ferritin levels according to − 597 A > G and − 174G > C genotypes showed no significant differences (all p > 0.05). Our findings provide valuable insights into the pathogenesis of COVID-19, suggesting that genetic variations at the IL-6 gene may contribute to the susceptibility to severe SARS-CoV-2 infection within the Moroccan population.

Peer Review reports

Introduction

The coronavirus disease 2019 (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) rapidly emerged as a global pandemic round the world. As of December 24, 2023, there have been 773,119,173 million confirmed cases and more than 6,990,067 million deaths worldwide due to COVID-19, posing an unprecedented challenge to humanity (https://data.who.int/dashboards/covid19/deaths?n=c).

The clinical picture of patients infected with SARS-CoV-2 ranges from asymptomatic to severe courses of COVID-19 with acute respiratory distress syndrome (ARDS), multiorgan involvement and death. The initial symptoms of COVID-19 mainly include fever, cough, headache, fatigue, breathing difficulties, anosmia and ageusia [1]. Previous studies suggested that SARS-CoV-2 infection triggers an immune dysregulation known as a cytokine storm (CS) characterized by the high-level of proinflammatory cytokines and other inflammatory markers and has been suggested to be associated with COVID-19 severity [2, 3]. This hyperinflammatory syndrome triggers the severe inflammatory response observed in SARS-CoV-2 infection and causes coagulopathies, oxidative stress, organ damage, and death [4]. Interleukin 6 (IL-6) is one of the most widely recognised cytokine markers of inflammation, providing clinicians with valuable information for the early identification of patients with severe COVID-19 during the course of the disease [5].

SARS-CoV-2 triggers the activation of pathogenic Th1 cells, leading to the secretion of proinflammatory cytokines, including granulocyte-macrophage colony-stimulating factor (GM-CSF) and IL-6. GM-CSF, in turn, stimulates CD14 + CD16 + inflammatory monocytes to generate substantial amounts of IL-6, tumor necrosis factor-α (TNF-α), and other cytokines [4]. Furthermore, IL-6 binds to the soluble form of the IL-6 receptor (sIL-6R) and the membrane-bound IL-6 receptor (mIL-6R) through gp130, inducing various proinflammatory cytokines and chemokines. This pleiotropic effect extends to acquired and innate immune cells, respectively, ultimately leading to cytokine storms [4, 6]. Interestingly, a previous study found that tocilizumab (TCZ), an IL-6 receptor antagonist capable of inhibiting cytokine storms by blocking the IL-6 signal transduction pathway, could decrease the likelihood of progression to the composite outcome of mechanical ventilation or death. However, it did not demonstrate an improvement in overall survival [7].

The human IL-6 gene is located on chromosome 7p21 and harbour many single-nucleotide polymorphisms (SNPs) in the coding and non-coding regions. The IL-6 gene possesses a 303 bp promoter, where the adenine (A) to guanine (G) substitution at the − 597 position (− 597 A > G; rs1800797) and the guanine (G) to cytosine (C) replacement at the − 174 position (− 174G/C; rs1800795) are notable [8]. Indeed, previous data have highlighted that genetic polymorphisms in the promoter region exert an influence on IL-6 transcription [9].

Genetic factors within the host organism can play a significant role in the development and severity of SARS-CoV-2 infection. This includes immune-related genes, which can ultimately influence the clinical outcome of the disease [1, 10]. Investigating the role of genetic factors in the susceptibility and severity of COVID-19 can enhance our understanding of the pathogenetic mechanisms underlying complications and fatalities associated with this disease. It has been flagged that research into human protein-coding genes is disproportionately skewed towards a comparably small set of genes, and that genes that are identified by genome-wide datasets, and hence likely to have biological significance in the context of COVID-19, are at a risk of remaining ignored by researchers [11, 12]. Meta-analysis of genome-wide overlap (MAGO), or meta-analysis of over-lapping SNVs from genome-wide approaches, provides a relevant means of aggregating evidence for causal association, namely through results from meta-analysis by information content (MAIC) algorithm [13]. Notably, this study focusses the IL-6 gene which is amongst those within the MAIC top 1000 scores (https://baillielab.net/maic/covid19/, version April 14, 2021).

Therefore, using a gene candidate approach, we aimed to analyse the impact of IL-6 variants (rs1800797 and rs1800795) on SARS-CoV-2 infection and the clinical progression of COVID-19 in a Moroccan population.

Materials and methods

Patients and study design

The protocol for patient and control recruitment received approval from the ethical committee of the Faculty of Medicine of Rabat (S/22). All procedures conducted in this study adhered to the ethical standards established by the institutional and national research committee, following the principles outlined in the 1964 Helsinki Declaration and its subsequent amendments or comparable ethical standards. Written consent was obtained from all participants.

A case-control study was conducted using blood samples from unvaccinated COVID-19 patients and SARS-CoV-2-unaffected controls. The sample size was estimated using the sample size calculator server, aiming for 80% power and an expected size of 2, as previously described [14]. We enrolled 270 unvaccinated COVID-19 patients. The diagnosis of SARS-CoV-2 was established through positive real-time reverse-transcriptase-polymerase chain reaction (qRT-PCR), conducted on oro- and nasopharyngeal swab specimens, complemented by SARS-CoV-2 serology. SARS-CoV-2 infection can be grouped into the following severity of illness categories: 132 with severe COVID-19 (hospitalized for SpO2 < 94%, PaO2/FiO2 < 300 mm Hg, a respiratory rate > 30 breaths/min, or lung infiltrates > 50%), and 138 with asymptomatic to moderate COVID-19 (mild cases), (non-hospitalized), displayed fewer respiratory symptoms and did not require hospitalization for COVID-19.

Additionally, 339 SARS-CoV-2-negative group participated in routine check-up visits between December 20, 2020, and December 31, 2021, at the Pasteur Institute of Morocco, Moulay Youssef Hospital, Casablanca, and Morizgo Laboratory. The negative group consisted of participants who were PCR negative for SARS-CoV-2 RNA and negative of SARS-CoV-2 antibodies against the SARS-CoV-2 nucleocapsid protein at the time of recruitment.

Complete blood Count (CBC)

Samples for the complete blood count were collected aseptically from each participant using a 4-mL EDTA tube (BD Microtainer, Franklin Lakes, New Jersey). The analysis was conducted using a Mindray BC-6800 analyzer (Mindray, China).

D-dimer, C-reactive protein, and ferritin levels

D-dimer levels were assessed using the ELFA (Enzyme-Linked Fluorescent Assay) technique on a VIDAS® instrument (BioMérieux, France). C-reactive protein and ferritin were collected in a 0.5 ml serum separator tube (BD Microtainer, Franklin Lakes, New Jersey) and analyzed using a clinical chemistry analyzer (Architect ci4100, Abbott Laboratories, Abbott Park, Illinois, USA).

SARS-CoV-2 serology

A serological assay was performed to ascertain the presence of SARS-CoV-2 antibodies against the SARS-CoV-2 nucleocapsid protein using the Architect™ SARS-CoV-2 chemiluminescence microparticles (CMIA) immunoassay (Abbott Park, Illinois, USA). The interpretation of results followed the manufacturer’s criteria: negative when the cutoff index (S/C) was < 1.4 and positive when the cutoff index (S/C) was ≥ 1.4. The SARS-CoV-2 IgG II Quant assay was utilized for the quantitative determination of IgG antibodies specific to the receptor-binding domain (RBD). This assay was conducted using the Abbott Architect i2000SR instrument manufactured by Abbott Laboratories, located in Abbott Park, Illinois.

DNA isolation and SNPs genotyping

Genomic DNA extraction from blood was performed using the QIAamp DNA Blood Mini Kit (Qiagen, Venlo, The Netherlands). The concentration and purity of DNA were assessed using the Qubit Fluorometer (Thermo Fisher Scientific, MA, USA), following the manufacturer’s manual.

Genotyping of − 174G/C (rs1800795; assay ID C___1839697_20) and − 597 A > G (rs1800797; assay ID C___1839695_20) were conducted using predesigned TaqMan SNP genotyping assays (Thermo Fisher Scientific, MA, USA). Approximately 10 ng of genomic DNA was amplified in a 10 µL reaction mixture in a 96-well plate, containing 1x TaqMan SNP genotyping master mix (Thermo Fisher Scientific, MA, USA), SNP genotyping assay, and completed with DNase-free water. Negative controls were included in every run. The thermal profile of quantitative polymerase chain reaction (qPCR) included preincubation at 60 °C for 1 min and then 95 °C for 10 min, followed by 45 cycles of denaturation at 95 °C for 15 s, and annealing/extension at 60 °C for 1 min. This was followed by postincubation at 60 °C for 1 min and a final cooling at 40 °C for 30 s. Genotype analysis was conducted using LightCycler 480 software (Roche Diagnostics Ltd, Ferrenstrasse, Rotkreuz, Switzerland).

Statistical analysis

For descriptive statistics, categorical variables are expressed as frequencies (%), while continuous variables are presented as the median with the minimum and maximum values. Continuous variables were compared using the Mann-Whitney U test or the Kruskal-Wallis test. Departures from Hardy-Weinberg equilibrium were assessed by comparing observed genotype frequencies with expected genotype frequencies in all groups, calculated using observed allele frequencies by the chi-square G test (‘Goodness of Fit’) with 1 degree of freedom. Odds ratios (OR) and 95% confidence intervals (CI) were estimated using logistic regression analysis and adjusted for significant covariates to examine the association between wild-type or variant genotypes among cases and controls, considering co-dominant, dominant, recessive, and overdominant inheritance models. Haplotypes frequencies were estimated by the expectation-maximization algorithm with unknown phase. The association of haplotypes with susceptibility and severity was investigated through a logistic regression model, with the most common haplotype serving as the reference. Odds ratios were adjusted for significant covariates. We conducted univariate and multivariate logistic regression analyses to explore the independent early predictors and risk factors associated with the disease severity of COVID-19. Statistical analyses were performed using R software for Windows, GraphPad Prism 6e (GraphPad Software, San Diego, CA, USA) and SNPStats (https://www.snpstats.net/start.htm) [15]. A p-value < 0.05 was considered statistically significant.

Results

Demographic characteristics of participants included in the study

A total of 270 unvaccinated patients with a confirmed diagnosis of COVID-19 were included in the study, consisting of 106 women and 164 men, with a median age of 52 years (range: 18–91 years). Among these, 132 patients experienced severe COVID-19 (62 women, 70 men) with a median age of 61 years (range: 21–91 years), while 138 patients had asymptomatic to moderate COVID-19 (44 women, 94 men) with a median age of 36 years (range: 18–63 years). A significant age difference was observed between severe and mild COVID-19 patients (p < 0.0001), with severe cases being older. Gender was also associated with the severity of COVID-19 (p = 0.011). Additionally, 339 SARS-CoV-2-negative group (224 women, 115 men) with a median age of 50 years (range: 19–93 years) were included.

Analysis of single nucleotide polymorphisms in the IL-6 gene and susceptibility to SARS-CoV-2 infection

The genotypic and allele frequencies of two SNPs in the IL-6 gene (− 597 A > G, − 174G > C) were analyzed in 270 SARS-CoV-2-infected patients and 339 uninfected SARS-CoV-2 group, adjusting for age and sex through multivariate logistic regression analyses (Table 1). Genotypic distributions of both variants were in Hardy-Weinberg equilibrium for SARS-CoV-2-negative group (p = 0.28 for rs1800795 and p = 0.25 for rs1800797) and cases (p = 0.13 for rs1800795 and p = 0.20 for rs1800797). No significant associations were found between the allele and genotype frequencies of the IL-6− 597 A > G and − 174G > C variants and susceptibility to SARS-CoV-2 infection (p > 0.05).

Table 1 Genotype and allelic frequencies of two IL–6 promoter polymorphisms among SARS-CoV-2 cases and uninfected- controls and associations with risk of COVID-19

However, a subsequent analysis revealed significant linkage disequilibrium between rs1800797 − 597 A > G) and rs1800795 (− 174G > C), suggesting that individuals with the GC* haplotype (OR = 0.04, 95%CI 0.01–0.30, p = 0.001 and AG* haplotype (OR = 0.11, 95%CI 0.03–0.46, p = 0.002) were significantly associated with protection against SARS-CoV-2 infection (Table 2).

Table 2 Association of IL–6 promoter haplotypes with risk of SARS-CoV-2 infection

Analysis of single nucleotide polymorphisms in the IL-6 gene and susceptibility to severe COVID-19

To investigate the influence of variants in the promoter region of the IL-6 gene on the severity of COVID-19, we genotyped the − 597 A > G and − 174G > C polymorphisms in 132 patients with severe COVID-19 and 138 asymptomatic-moderate patients (mild cases) (Table 3). The genotype distributions for both SNPs were in Hardy-Weinberg equilibrium in the mild group (p = 1.00 for both rs1800795 and rs1800797 SNPs) but not in the severe COVID-19 cases (p = 0.01 for rs1800795 and p = 0.02 for rs1800797).

Table 3 Allele frequencies of two IL-6 polymorphisms among COVID-19 patients with severe and mild symptoms

In the overdominant model, the frequency of the G/C genotype of the rs1800795 polymorphism (− 174G > C) was significantly higher in mild cases than in severe patients (p = 0.03). Individuals with the G/C genotype were protected against progression to severe disease compared with those with the G/G or C/C genotypes (p = 0.03; OR = 0.41, 95%CI 0.18–0.96) (Table 3). After the Bonferroni correction (adjusted alpha = 0.005), we observed no significant difference in the rs1800795 polymorphism (− 174G > C) between severe and mild cases.

Regarding the rs1800797 (− 597 A > G) polymorphism, there was no significant statistical association between different genotypes and severity under various genetic models (p > 0.05). Haplotype analysis of the two IL-6 gene polymorphism loci (− 597 A > G and − 174G > C) was performed, as shown in Table 4. The probability of developing a severe form of COVID-19 was not linked to any haplotype (p > 0.05).

Table 4 Association between COVID-19 severity and major haplotypes of rs1800797 and rs1800795

Multivariate analysis indicated that the age of patients (OR = 1.06; 95% CI, 1.04–1.08, < 0.001), D-dimer (OR = 1.40; 95% CI, 1.06–1.56, p = 0.002), and CRP levels on admission (OR = 1.12; 95% CI, 1.06–1.26, p = 0.003) were risk factors associated with severe COVID-19 (Supplementary Table 1).

IL-6 gene polymorphisms and biological markers

We conducted an analysis to explore the potential impact of IL-6 polymorphisms, specifically rs1800795 and rs1800797, on complete blood count (CBC) and D-dimer levels in COVID-19 patients. The results are depicted in Figs. 1 and 2 for rs1800795 and rs1800797, respectively. Our findings revealed no significant differences between genotypes and the levels of hematological markers, C-reactive protein, and ferritin levels, based on the − 597 A > G and − 174G > C genotypes.

Fig. 1
figure 1

Biological markers in the COVID-19 cases according to IL-6 rs1800795 variant. A) IL-6 rs1800795 genotypes and leukocytes count. B) IL-6 rs1800795 genotypes and polynuclear neutrophils count. C) IL-6 rs1800795 genotypes and polynuclear eosinophils count. D) IL-6 rs1800795 genotypes and polynuclear basophils count. E) IL-6 rs1800795 genotypes and lymphocytes count. F) IL-6 rs1800795 genotypes and monocytes count. G) IL-6 rs1800795 genotypes and red blood cells count. H) IL-6 rs1800795 genotypes and platelets count. I) IL-6 rs1800795 genotypes and D-dimers count. J) IL-6 rs1800795 genotypes and C-reactive protein levels. K) IL-6 rs1800795 genotypes and ferritin levels. Data are presented as a scatter plot with median. Statistical tests were performed using the Mann Whitney U test

Fig. 2
figure 2

Biological markers in the COVID-19 cases according to IL-6 rs180097 polymorphism. A) IL-6 rs180097 genotypes and leukocytes count. B) IL-6 rs180097 genotypes and polynuclear neutrophils count. C) IL-6 rs180097 genotypes and polynuclear eosinophils count. D) IL-6 rs180097 genotypes and polynuclear basophils count. E) IL-6 rs180097 genotypes and lymphocytes count. F) IL-6 rs180097 genotypes and monocytes count. G) IL-6 rs180097 genotypes and red blood cells count. H) IL-6 rs180097 genotypes and platelets count. I) IL-6 rs180097 genotypes and D-dimers count. J) IL-6 rs180097 genotypes and C-reactive protein levels. K) IL-6 rs180097 genotypes and ferritin levels. Data are presented as a scatter plot with median. Statistical tests were performed using the Mann Whitney U test

Discussion

The clinical presentation of COVID-19 varies among patients. According to early Moroccan data, 64.4% of patients developed non-severe symptoms, while 33.6% were admitted to the intensive care unit [16]. An individual’s genetic background may influence the outcome of SARS-CoV-2 infection. Variability in susceptibility to SARS-CoV-2 and the severity of COVID-19 involves multiple factors, encompassing both viral and host elements, with genetic factors playing a role in affecting proinflammatory cascades [17]. IL-6, a pro-inflammatory cytokine, plays a crucial role in the immune response and inflammation. Previous suggestions have implicated its involvement in the severity of COVID-19, with elevated IL-6 levels being associated with a more severe disease course [18, 19]. Understanding the genetic factors that influence the severity of COVID-19 disease has clinical implications. Identifying individuals with potentially protective genotypes can contribute to risk stratification and the development of personalised treatment strategies. In our study, a significant age difference was observed between severe and mild COVID-19 patients, supporting the well-established finding that age plays an important role in disease severity [20,21,22,23]. Furthermore, the association between male gender and the severity of COVID-19 is consistent with previous observations. These findings underscore that sex hormones, angiotensin-converting enzyme levels, and severe inflammation in men may contribute to the progression from mild to severe cases [20, 24,25,26].

Earlier data revealed a systemic pro-inflammatory signature, in particular elevated levels of IL-6, which were associated with both severity of organ failure and 60-day mortality [27]. In this study, we explored the association between IL-6 polymorphisms (− 174G > C and − 597 A > G) and susceptibility to SARS-CoV-2 infection. However, allele and genotype frequencies revealed no significant association. These results are consistent with those of a study conducted in the Chinese population, which also reported non-significant associations in allele and genotype frequencies between healthy controls and COVID-19 patients [28].

Interestingly, our study showed interesting results when haplotypes were taken into account. Linkage disequilibrium analysis showed that individuals with specific haplotypes (GC* and AG*) were significantly associated with protection against SARS-CoV-2 infection. This finding suggests that certain combinations of alleles may confer a protective effect, potentially influencing the host response to the virus. This observation is in line with previous in vitro studies indicating that genetic polymorphisms in the IL-6 gene promoter influence its transcription by a complex haplotype-determined interaction, rather than by a simple additive mechanism [9]. Polymorphisms, such as the − 597 A > G and − 174G > C SNPs, located adjacent to the negative regulatory domain (nucleotides − 173 to − 151 relative to the transcription start site) and in close proximity to an NF-IL6 binding motif, appear to contribute to interindividual variation in IL-6 transcription and expression [28]. Additionally, in an ex vivo model, IL-6 promoter haplotypes defined by the − 597 and − 174 SNPs are associated with IL-6 production stimulated by endotoxin [28]. An early Turkish study found that the GG genotype of the − 174G > C polymorphism was significantly higher in COVID-19 patients with macrophagic activation syndrome. The G allele may serve as a risk factor for elevated serum IL-6 levels and progression to macrophagic activation syndrome [29]. An investigation conducted in Iraq supports the concept that the IL-6− 597 A > G variant is associated with susceptibility to HBV and HCV infections [30]. However, a Chinese study found no association between HBV resistance and the − 174G > C and − 597 A > G variants [31]. Interestingly, among HBV patients, the AG* haplotype was significantly associated with an increased risk of developing HBV infection (OR = 5.45; 95% CI: 3.63–8.18; p = 0.000), whereas a protective effect of the GC* haplotype against the progression of HBV infection was suggested. For HCV, the AC* haplotype was significantly associated with the risk of HCV infection (OR = 1.96, 95% CI: 1.31–2.93; p = 0.001) [30].

Under the overdominant model analysis, we demonstrated that the IL-6 − 174 G/C genotype confers protection against the development of severe COVID-19 after adjusting for sex and age. Consistent with a previous Chinese study, it was identified that rs1800796 C/C and rs1524107 T/T protect against severe COVID-19 [28]. Consistent with a previous meta-analysis, it was also shown that the C allele of the − 174G > C SNP is linked to higher IL-6 production and more severe COVID-19 in the Caucasian population [32]. Furthermore, a previous study demonstrated a high correlation between the C allele and elevated expression of IL-6 compared to the ancestral allele (G) [28]. In contrast to our study, an Indian report showed that individuals with the G/C genotype of the IL-6 (–174G > C) polymorphism have a significantly higher risk of severity (adjusted OR = 3.86, p < 0.001) [33].

For the rs1800797 (− 597 A > G) polymorphism, we did not find any association with the severity of COVID-19. Previous data from an Indian study [33] and the Iranian population [34] align with our findings and showed no association of COVID-19 severity with the rs1800797 (− 597 A > G) SNP [33].

Infection with SARS-CoV-2 induces an inflammatory response, resulting in a significant release of pro-inflammatory cytokines, leading to a cytokine storm. Numerous studies have investigated the cytokine profiles of individuals with COVID-19, revealing a strong correlation between the cytokine storm and a poor prognosis in severe cases of COVID-19 [4].

Furthermore, the probability of developing a severe form of COVID-19 was not linked to any haplotype. However, an Asian study showed that the common CTT* haplotype, represented by the rs1800796, rs1524107 and rs2066992 alleles, contributed to a better outcome from SARS-CoV-2 infection and was associated with a reduction in IL-6 expression [28].

These differences in results could be attributed to ethnic differences, as the viruses are also genetically distinct. The underlying immunopathogenesis of viral infections is thought to involve different mechanisms, influencing the outcome [35].

The absence of a statistically significant association between these IL-6 genotypes and the levels of complete blood count, D-dimer, CRP, and ferritin levels is consistent with findings from previous studies in other inflammatory contexts [36,37,38,39].

It’s important to acknowledge the limitations of the study, such as the relatively small sample size and potential confounding factors. Future research involving genome-wide association study (GWAS) with larger and more diverse cohorts and ex vivo functional studies could further validate these findings. Additionally, more attention should be paid to evaluating IL-6 serum levels in patients with GC* and AG* haplotypes in order to establish the genotype-to-phenotype relationship.

In conclusion, the GC* and AG* haplotypes of rs1800795 (− 174 G > C) and rs1800797 (− 597 A > G) variants in the IL-6 gene may serve as useful genetic markers to screen for high-risk SARS-CoV-2 infection in the Moroccan population, providing valuable insights into the genetic determinants of COVID-19 severity. The study’s results may have implications for public health strategies, such as identifying high-risk populations based on genetic factors and tailoring preventive measures accordingly. Further study is warranted to confirm these findings in larger cohorts and diverse populations, paving the way for more targeted approaches in managing and predicting severe COVID-19 cases.

Data availability

The data presented in this study are available upon reseanable request from the corresponding author. The data are not publicly available according to the ethical committee decision on the conduct of this study.

References

  1. Velavan TP, Pallerla SR, Rüter J, Augustin Y, Kremsner PG, Krishna S, Meyer CG. Host genetic factors determining COVID-19 susceptibility and severity. EBioMedicine. 2021;72:103629.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, Zhang L, Fan G, Xu J, Gu X, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395(10223):497–506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Mehta P, McAuley DF, Brown M, Sanchez E, Tattersall RS, Manson JJ. HLH across Speciality Collaboration U: COVID-19: consider cytokine storm syndromes and immunosuppression. Lancet. 2020;395(10229):1033–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Hu B, Huang S, Yin L. The cytokine storm and COVID-19. J Med Virol. 2021;93(1):250–6.

    Article  CAS  PubMed  Google Scholar 

  5. Aziz M, Fatima R, Assaly R. Elevated interleukin-6 and severe COVID-19: a meta-analysis. J Med Virol. 2020;92(11):2283–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Murakami M, Kamimura D, Hirano T. Pleiotropy and specificity: insights from the interleukin 6 family of cytokines. Immunity. 2019;50(4):812–31.

    Article  CAS  PubMed  Google Scholar 

  7. Salama C, Han J, Yau L, Reiss WG, Kramer B, Neidhart JD, Criner GJ, Kaplan-Lewis E, Baden R, Pandit L, et al. Tocilizumab in patients hospitalized with Covid-19 pneumonia. N Engl J Med. 2021;384(1):20–30.

    Article  CAS  PubMed  Google Scholar 

  8. Fishman D, Faulds G, Jeffery R, Mohamed-Ali V, Yudkin JS, Humphries S, Woo P. The effect of novel polymorphisms in the interleukin-6 (IL-6) gene on IL-6 transcription and plasma IL-6 levels, and an association with systemic-onset juvenile chronic arthritis. J Clin Invest. 1998;102(7):1369–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Terry CF, Loukaci V, Green FR. Cooperative influence of genetic polymorphisms on interleukin 6 transcriptional regulation. J Biol Chem. 2000;275(24):18138–44.

    Article  CAS  PubMed  Google Scholar 

  10. Vakil MK, Mansoori Y, Al-Awsi GRL, Hosseinipour A, Ahsant S, Ahmadi S, Ekrahi M, Montaseri Z, Pezeshki B, Mohaghegh P, et al. Individual genetic variability mainly of proinflammatory cytokines, cytokine receptors, and toll-like receptors dictates pathophysiology of COVID-19 disease. J Med Virol. 2022;94(9):4088–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. David S, Dorado G, Duarte EL, David-Bosne S, Trigueiro-Louro J, Rebelo-de-Andrade H. COVID-19: impact on Public Health and hypothesis-driven investigations on genetic susceptibility and severity. Immunogenetics. 2022;74(4):381–407.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Stoeger T, Nunes Amaral LA. COVID-19 research risks ignoring important host genes due to pre-established research patterns. Elife. 2020;9:e61981.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Parkinson N, Rodgers N, Head Fourman M, Wang B, Zechner M, Swets MC, Millar JE, Law A, Russell CD, Baillie JK, et al. Dynamic data-driven meta-analysis for prioritisation of host genes implicated in COVID-19. Sci Rep. 2020;10(1):22303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Noordzij M, Tripepi G, Dekker FW, Zoccali C, Tanck MW, Jager KJ. Sample size calculations: basic principles and common pitfalls. Nephrol Dial Transpl. 2010;25(5):1388–93.

    Article  Google Scholar 

  15. Solé X, Guinó E, Valls J, Iniesta R, Moreno V. SNPStats: a web tool for the analysis of association studies. Bioinformatics. 2006;22(15):1928–9.

    Article  PubMed  Google Scholar 

  16. El Aidaoui K, Haoudar A, Khalis M, Kantri A, Ziati J, El Ghanmi A, Bennis G, El Yamani K, Dini N, El Kettani C. Predictors of severity in Covid-19 patients in Casablanca, Morocco. Cureus. 2020;12(9):e10716.

    PubMed  PubMed Central  Google Scholar 

  17. Kaidashev I, Shlykova O, Izmailova O, Torubara O, Yushchenko Y, Tyshkovska T, Kyslyi V, Belyaeva A, Maryniak D. Host gene variability and SARS-CoV-2 infection: a review article. Heliyon. 2021;7(8):e07863.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Nikkhoo B, Mohammadi M, Hasani S, Sigari N, Borhani A, Ramezani C, Charajoo A, Badri S, Rostami F, Etemadi M, et al. Elevated interleukin (IL)-6 as a predictor of disease severity among Covid-19 patients: a prospective cohort study. BMC Infect Dis. 2023;23(1):311.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Guirao JJ, Cabrera CM, Jiménez N, Rincón L, Urra JM. High serum IL-6 values increase the risk of mortality and the severity of pneumonia in patients diagnosed with COVID-19. Mol Immunol. 2020;128:64–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Jin JM, Bai P, He W, Wu F, Liu XF, Han DM, Liu S, Yang JK. Gender differences in patients with COVID-19: focus on severity and mortality. Front Public Health. 2020;8:152.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Zimmermann P, Curtis N. Why does the severity of COVID-19 Differ with Age? Understanding the mechanisms underlying the Age Gradient in Outcome following SARS-CoV-2 infection. Pediatr Infect Dis J. 2022;41(2):e36–45.

    Article  PubMed  Google Scholar 

  22. Liu Y, Mao B, Liang S, Yang JW, Lu HW, Chai YH, Wang L, Zhang L, Li QH, Zhao L, et al. Association between age and clinical characteristics and outcomes of COVID-19. Eur Respir J. 2020;55(5):2001112.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Herrera-Esposito D, de Los Campos G. Age-specific rate of severe and critical SARS-CoV-2 infections estimated with multi-country seroprevalence studies. BMC Infect Dis. 2022;22(1):311.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Twitchell DK, Christensen MB, Hackett G, Morgentaler A, Saad F, Pastuszak AW. Examining male predominance of severe COVID-19 outcomes: a systematic review. Androg Clin Res Ther. 2022;3(1):41–53.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Acheampong DO, Barffour IK, Boye A, Aninagyei E, Ocansey S, Morna MT. Male predisposition to severe COVID-19: review of evidence and potential therapeutic prospects. Biomed Pharmacother. 2020;131:110748.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Bwire GM. Coronavirus: why men are more vulnerable to Covid-19 Than women? SN Compr Clin Med. 2020;2(7):874–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Lavillegrand JR, Garnier M, Spaeth A, Mario N, Hariri G, Pilon A, Berti E, Fieux F, Thietart S, Urbina T, et al. Elevated plasma IL-6 and CRP levels are associated with adverse clinical outcomes and death in critically ill SARS-CoV-2 patients: inflammatory response of SARS-CoV-2 patients. Ann Intensive Care. 2021;11(1):9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Chen T, Lin YX, Zha Y, Sun Y, Tian J, Yang Z, Lin SW, Yu F, Chen ZS, Kuang BH et al. A Low-Producing Haplotype of Interleukin-6 Disrupting CTCF Binding Is Protective against Severe COVID-19. mBio 2021, 12(5):e0137221.

  29. Kerget F, Kerget B. Frequency of Interleukin-6 rs1800795 (-174G/C) and rs1800797 (-597G/A) polymorphisms in COVID-19 patients in Turkey who develop macrophage activation syndrome. Jpn J Infect Dis. 2021;74(6):543–8.

    Article  CAS  PubMed  Google Scholar 

  30. Al-Saffar OB, Ad’hiah AH. Genetic variants in IL4RA, IL6, and IL12B genes and susceptibility to hepatitis B and C virus infections among Iraqi patients. J Med Virol. 2020;92(12):3448–58.

    Article  CAS  PubMed  Google Scholar 

  31. Dai Y, Liu XL, Chai QB, Gai YH, Li Y, Liu BG. Relationship between the-572G/C polymorphism in IL-6 gene promoter and Hepatitis virus B infection in Chinese Han population. Wor Chin J Dig 2009:1522–6.

  32. Ulhaq ZS, Soraya GV. Anti-IL-6 receptor antibody treatment for severe COVID-19 and the potential implication of IL-6 gene polymorphisms in novel coronavirus pneumonia. Med Clin (Engl Ed). 2020;155(12):548–56.

    CAS  PubMed  Google Scholar 

  33. Verma S, Verma S, Khan FH, Siddiqi Z, Raza ST, Abbas M, Mahdi F. Genetic polymorphisms of IL6 gene – 174G > C and – 597G > A are associated with the risk of COVID-19 severity. Int J Immunogenet. 2023;50(1):5–11.

    Article  CAS  PubMed  Google Scholar 

  34. Falahi S, Zamanian MH, Feizollahi P, Rezaiemanesh A, Salari F, Mahmoudi Z, Gorgin Karaji A. Evaluation of the relationship between IL-6 gene single nucleotide polymorphisms and the severity of COVID-19 in an Iranian population. Cytokine 2022, 154(155889).

  35. Rouse BT, Sehrawat S. Immunity and immunopathology to viruses: what decides the outcome? Nat Rev Immunol. 2010;10(7):514–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Brull DJ, Montgomery HE, Sanders J, Dhamrait S, Luong L, Rumley A, Lowe GD, Humphries SE. Interleukin-6 gene – 174g > c and – 572g > c promoter polymorphisms are strong predictors of plasma interleukin-6 levels after coronary artery bypass surgery. Arterioscler Thromb Vasc Biol. 2001;21(9):1458–63.

    Article  CAS  PubMed  Google Scholar 

  37. Kilpinen S, Hulkkonen J, Wang XY, Hurme M. The promoter polymorphism of the interleukin-6 gene regulates interleukin-6 production in neonates but not in adults. Eur Cytokine Netw. 2001;12(1):62–8.

    CAS  PubMed  Google Scholar 

  38. Humphries SE, Luong LA, Ogg MS, Hawe E, Miller GJ. The interleukin-6 -174 G/C promoter polymorphism is associated with risk of coronary heart disease and systolic blood pressure in healthy men. Eur Heart J. 2001;22(24):2243–52.

    Article  CAS  PubMed  Google Scholar 

  39. Basso F, Lowe GD, Rumley A, McMahon AD, Humphries SE. Interleukin-6 -174G > C polymorphism and risk of coronary heart disease in West of Scotland coronary prevention study (WOSCOPS). Arterioscler Thromb Vasc Biol. 2002;22(4):599–604.

    Article  CAS  PubMed  Google Scholar 

Download references

Funding

This study was supported by a ACIP (Action Concertées Interpasteuriennes) grant (ACIP n°505) from Institut Pasteur (Paris) and by the Taskforce COVID-19 of Institut Pasteur.

Author information

Authors and Affiliations

Authors

Contributions

Sayeh Ezzikouri contributed to the conception of the study and the acquisition of funding. Rachid Noureddine drafted the manuscript. Rachid Noureddine, Hanâ Baba, Safaa Aqillouch, Karima Abounouh, Oumaima Laazaazia, Mohcine Elmessaoudi-Idrissi, Fatima Zohra Bahmani, Ikram Allah Tanouti, and Ahd Ouladlahsen contributed to the sample collection, performed experiments in the laboratory, conducted statistical analyses, and interpreted the data. Hind Dehbi contributed to supervision. Sayeh Ezzikouri and M’hammed Sarih reviewed and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Sayeh Ezzikouri.

Ethics declarations

Ethics approval and consent to participate

This study was conducted in accordance with the Declaration of Helsinki. The Ethics Committee of the Faculty of Medicine of Rabat (S/22) approved the study. Informed consent was obtained from all subjects involved in the study.

Consent for publication

Not Applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Noureddine, R., Baba, H., Aqillouch, S. et al. The Interleukin-6 gene variants may protect against SARS-CoV-2 infection and the severity of COVID-19: a case-control study in a Moroccan population. BMC Med Genomics 17, 139 (2024). https://doi.org/10.1186/s12920-024-01911-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-024-01911-w

Keywords