1. Introduction
The COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has been a global challenge since its emergence [1]. In Brazil, the number of confirmed cases surpassed thirty-eight million, making it the fifth most affected country in the world. Furthermore, with over 712,000 registered deaths, Brazil holds the second position in mortality rate, according to the Health Ministry (
Additionally, high levels of IL-6, as observed in various COVID-19 studies, correlate with predictive mortality and association with respiratory failure, shock, and severe symptoms in hospitalized patients [16,17,18]. Similarly, elevated TNF-α negatively impacts the lower respiratory system, resulting in airway narrowing, bronchial hyperreactivity, damage to the respiratory epithelium, and pulmonary fibrosis [19,20,21]. Several studies reviewed by Qudus et al. (2023) [22] have associated genetic variants with cytokine expression and the development of ARDS through cytokine storm mechanisms. Genetic variants in the promoter region of the IL1B gene have been associated with negative outcomes in inflammatory conditions, such as rs4848306 in reactive septic arthritis and rs16944 in chronic obstructive pulmonary disease (COPD). A study showed that the T allele of IL1B +3953C > T provides protection against severe COVID-19, while the CT genotype increases the risk of severity [23]. For the IL6 gene, intronic variants have been related to respiratory diseases, including ARDS and COPD, with the GG genotype of rs1800795 associated with severe COVID-19 cases among Kurdish patients [24]. Variants in the promoter region of the TNF gene are associated with a range of clinical conditions and severe outcomes. Additionally, the A allele of the rs1800629 variant has been correlated with increased severity and higher mortality risk associated with COVID-19, suggesting its potential as a genetic marker for predicting adverse outcomes [25]. Other TNF variants such as rs361525G, rs1799724, and rs1800630A have also been associated with inflammatory syndromes, like Thrombocytopenia Syndrome (SFTS) and ARDS.
These findings illustrate how genetic variants, particularly in regulatory regions of cytokine genes, may influence host susceptibility to severe disease outcomes. The immune response to SARS-CoV-2 is shaped by a combination of environmental and genetic factors, and the variability in COVID-19 severity—from mild symptoms to fatal complications—highlights the importance of exploring genetic contributions. Despite growing evidence, data on the role of specific immunogenetic variants remain limited, especially regarding their impact on genetically diverse populations. Understanding these associations is essential to clarifying individual variability in immune response and disease progression.
The thirteen SNPs selected in this study are located in promoter or enhancer regions of the IL1B, IL6, and TNF genes and were chosen based on evidence from the HaploReg [26] and ENCODE databases [27], indicating functional regulatory activity. These variants have been associated with transcriptional modulation, altered binding motifs, and epigenetic marks—e.g., histone modifications— in immune-relevant tissues, supporting their role in inflammatory and respiratory responses. Although self-reported ancestry data were not available, the study population represents a genetically admixed cohort from southern Brazil, predominantly of European, African, and Native American backgrounds. To date, few studies have evaluated the association of functionally relevant variants in cytokine genes with COVID-19 severity in Brazilian populations—particularly in the South, where genetic profiles and epidemiological patterns differ from other regions [28]. This study aims to address this gap by investigating the association between thirteen functional variants in IL1B, IL6, and TNF and severe COVID-19 outcomes in 500 unvaccinated individuals from southern Brazil.
2. Materials and Methods
2.1. Study Design and Data Collection Methodology
The sample consisted of 500 patients who were diagnosed with COVID-19 through the real-time reverse transcription polymerase chain reaction (RT-PCR) test for the detection of SARS-CoV-2 viral RNA, using oro-nasopharyngeal swabs. These samples were collected from April 2020 to March 2021. All the samples were provided by the Central Laboratory of Rio Grande do Sul (LACEN-RS), one of the main laboratories responsible for COVID-19 testing during the pandemic in Rio Grande do Sul, the most southern state of Brazil. Demographic and clinical data were collected from the Brazilian notification system SIVEP-GRIPE (Influenza Epidemiological Surveillance Information System—Severe Acute Respiratory Syndrome Notifications). The clinical and demographic information collected comprised gender, age, permanent place of residence, and comorbidities (including heart diseases, lung diseases, obesity, hematological diseases, Down Syndrome, liver diseases, asthma, diabetes mellitus, immunodeficiencies, kidney diseases, or neurological diseases). Regarding COVID-19, the date of symptom onset, date of the RT-PCR test, date of hospital admission, and symptoms at the time of COVID-19 diagnosis, such as fever, cough, sore throat, difficulty breathing, oxygen saturation rate < 95%, fatigue, and loss of smell or taste, were collected from the form of COVID-19 notification. For COVID-19 cases that required hospitalization, data related to oxygen therapy, intensive care unit care, and outcome (death or recovery) were also recorded from medical records. Following processing and filtering of clinical data, participants were classified into two main groups: hospitalized and non-hospitalized. Additionally, patients were further categorized based on their clinical condition into the following groups: (1) requiring intensive care unit (ICU) admission versus not requiring ICU admission and (2) deceased versus discharged from the hospital. Given that participants were classified according to disease outcomes, this study follows a retrospective case–control design. Clinical outcomes were classified using standardized criteria based on national COVID-19 guidelines [29]. Clinical severity levels were determined using a hierarchical logic as follows:
0 (Asymptomatic): Case explicitly recorded as “asymptomatic” in the system;
1 (Mild): No ICU admission, no ventilatory support, and no death;
2 (Moderate): Received ventilatory support (non-invasive or invasive), but no ICU admission or death;
3 (Severe): ICU admission and/or death.
For statistical analyses, individuals classified as moderate (2) and severe (3) were grouped as hospitalized patients. Apart from a confirmed COVID-19 diagnosis, no further exclusion criteria were applied. After curating the clinical data, cases with incomplete or inconclusive information were removed from the final analysis.
2.2. DNA Extraction and Genotyping Methods
Genomic DNA was extracted from nasopharyngeal swab samples using the DNA extraction methodology developed by Sbruzzi et al., 2024 [30]. Subsequently, the samples were genotyped for 13 functional variants in the TNF, IL1B, and IL6 genes (Table 1), which were selected based on their known roles in immune response modulation and associations with inflammatory diseases.
Genotyping was performed using TaqMan® assays (Thermo Fisher Scientific, Waltham, MA, USA). Real-time PCR amplification was carried out on the StepOnePlus™ Real-Time PCR System, following standard protocols. Quality control criteria included call rates > 98% for all SNPs and >99% concordance in duplicate samples. All genotyped SNPs were evaluated for minor allele frequency (MAF) and tested for deviation from the Hardy–Weinberg equilibrium (HWE). SNPs showing significant deviation from the HWE were further investigated by stratified testing in cases and controls. In addition, all discrepant genotypes were re-genotyped using validated internal controls to confirm the accuracy of the calls. Only SNPs with MAF > 0.01 and appropriate genotyping quality were retained for downstream analysis.
2.3. Genetic Association Analysis
To investigate the relationship between cytokine gene variants and COVID-19 severity—including hospitalization, ICU admission, and death—both univariate and multivariate analyses were conducted. Allele and genotype frequencies were compared using Pearson’s chi-square test or Fisher’s exact test. Three genetic inheritance models—additive, dominant, and recessive—were tested for each SNP. Univariate p-values were adjusted using the Bonferroni correction to control for multiple testing across the 13 SNPs. Variants with Bonferroni-adjusted p-values ≤ 0.100 were retained for multivariate analysis. Multivariate logistic regression models were performed using the Generalized Linear Models (GENLIN) procedure in IBM SPSS Statistics v.18.0.0 (IBM, Armonk, NY, USA;
3. Results
3.1. Symptoms and Comorbidities Associations
The demographic and clinical characteristics of the study sample, stratified by hospitalization status, are presented in Table 2.
Samples with DNA concentrations below 1.5 ng/µL and incomplete clinical data were excluded. The final study included 231 individuals who required hospitalization (46.2%) and 235 subjects who did not require hospital admission (47.0%) (Table 2). Non-hospitalized individuals displayed a significantly lower median age of 38 years (IQR 29.0; 53.0) compared to the hospitalized group (median 65 years; IQR 51.0; 77.0; (p < 0.01). There was a higher proportion of females (56.5%) in the non-hospitalized group, while males were more prevalent (56.4%) among the hospitalized cases (p = 0.006). Higher prevalence rates of comorbidities were observed among hospitalized individuals, including obesity (81.8%), heart diseases (81.3%), diabetes mellitus (80.5%), and chronic lung diseases (79.3%). These findings showed a strong association between these comorbidities and the need for hospitalization. Regarding the severity indicators, statistical analysis of the hospitalized cases revealed significant p-values (<0.01), further supporting the findings. These indicators include the need for hospitalization, oxygen therapy, intensive care unit (ICU) admission, and the outcome (death). Specifically, 51.8% (n = 108) of patients were admitted to critical care, 93.8% (n = 213) required ventilatory support, and 54.6% (n = 125) experienced mortality (Table 2).
3.2. Significant Associations Between IL1B and TNF Genetic Variants and COVID-19 Severity
Among the thirteen genotyped variants, six SNPs in the IL1B and TNF genes demonstrated the strongest and most consistent associations with COVID-19 severity, particularly in relation to hospitalization, ICU admission, and mortality. We evaluated the HWE and MAF for all 13 genotyped SNPs. The variants IL1B rs4848306 and TNF rs1800630 showed significant HWE deviations in the control group (IL1B rs4848306: p ≈ 0.04; TNF rs1800630: p < 0.001). Re-genotyping and comparison with internal samples and population-based data from the ABraOM database confirmed the consistency of the genotyping results. Thus, we maintain both of them in further analyses.
Analysis of genotypic distribution, including dominant and recessive models, was performed and shown in Supplementary Table S1. The results indicated that the AG genotype of the IL1B rs16944 variant was significantly associated with hospitalization, with an odds ratio (OR) of 1.98 (95% CI: 1.22–3.23, p = 0.006) (Table 3).
Similarly, the comparison between the AA + AG genotypes and GG also demonstrated an association with hospitalization, with an OR of 1.62 (95% CI: 1.04–2.55, p = 0.036). For the IL1B rs1143627 variant, the CT genotype was associated with hospitalization compared to the CC genotype, with an OR of 1.82 (95% CI: 1.02–3.27, p = 0.043). In the case of TNF gene variants, the CT genotype for rs1799964 showed a significant association with hospitalization (OR = 1.97; 95% CI: 1.22–3.22, p = 0.006). The comparison between the CC + CT and TT genotypes also revealed a significant association, with an OR of 1.93 (95% CI: 1.23–3.06, p = 0.005). For rs1800630, the AA genotype was significantly associated with hospitalization, with an OR of 2.37 (95% CI: 1.08–5.33, p = 0.034). The comparison between the AA + AC and CC genotypes was also significant, with an OR of 1.76 (95% CI: 1.09–2.84, p = 0.021) (Table 3).
When analyzing susceptibility to severity associated with genetic variants, the genes IL1B, IL6, and TNF exhibited statistically significant allelic and genotypic distributions (see Supplementary Tables S2 and S3). Regarding the need for ICU admission, the results revealed that for the IL1B rs16944 variant, individuals with the AG genotype presented an 83.4% higher risk of ICU admission compared to those with the GG genotype (OR = 1.83; 95% CI = 1.09–3.10; p = 0.022) (Table 4).
Additionally, for the IL1B rs1143627 variant, individuals with the CT genotype have a 108.2% higher risk of ICU admission compared to those with the CC genotype (OR = 2.08; 95% CI = 1.12–3.85; p = 0.019). Regarding mortality, the CC homozygote for rs1799964 presented an OR of 3.73 (95% CI: 1.21–14.37), indicating that individuals with this genotype have a significantly higher likelihood of death compared to those with the TT genotype (p = 0.034) (Table 4). Furthermore, the comparison between the CC + CT and TT genotypes yielded an OR of 1.66 (95% CI: 1.01–2.77), reinforcing the association of the C allele with an increased risk of death relative to the T allele (p = 0.049) (Table 5).
Taken together, the results indicate that four functional variants—IL1B rs16944 and rs1143627, and TNF rs1799964 and rs1800630—were consistently associated with worse clinical outcomes across different genetic models. These associations were observed for hospitalization, ICU admission, and mortality, reinforcing the role of host genetic variation in inflammatory mediators as a potential determinant of COVID-19 severity.
4. Discussion
The primary objective of this study was to investigate the correlation between genetic variations in IL1B, IL6, TNF, and critical outcomes of COVID-19. The heterozygous genotype of the rs16944 (IL1B) variant and the presence of the A allele showed a significant association with increased susceptibility to severe COVID-19 outcomes, including hospitalization and ICU admission, suggesting a possible contributory role in disease severity. This finding is consistent with the data reported by Balzanelli et al. (2022) [31], who identified the AG genotype as a potential risk factor for COVID-19. The A allele has also been linked to other severe inflammatory conditions, such as chronic obstructive pulmonary disease (COPD) and septic shock, suggesting that it may predispose individuals to exaggerated inflammatory responses. These immune reactions, particularly cytokine storms, are also observed in severe COVID-19 cases and may help contextualize a potential role for IL1B variants in modulating disease severity. Specifically, the A allele is known to regulate IL-1β expression, which leads to the overproduction of this pro-inflammatory cytokine in various conditions, including lung cancer [32]. Moreover, the A allele’s association with increased COPD risk in East Asian populations [33] and a higher likelihood of septic shock and death in patients with systemic inflammatory response syndrome (SIRS) [34] raises the possibility that it may promote a stronger inflammatory response, potentially contributing to more severe COVID-19 progression. The rs1143627 variant, in its heterozygous genotype, was significantly associated with an increased risk of hospitalization and ICU admission. This variant has been implicated in various clinical outcomes related to inflammatory and infectious diseases. The C allele is considered a potential independent risk factor for ICU admission in COVID-19 cases [35], and the CC genotype is associated with an increased risk of severe COVID-19 [36]. A key point is that the presence of the C allele disrupts the TATA box in the promoter region of the gene, reducing IL1B expression [32,37], as corroborated by studies showing a lower expression rate of this gene [38]. In contrast, the T allele is associated with increased IL-1β expression, which heightens the risk of tuberculosis [39].
The heterozygous genotype of rs1799964 (TNF) was significantly associated with COVID-19 hospitalizations, while individuals with the CC genotype had an increased risk of death. Compared to TT, the CC and CT genotypes were associated with higher odds of hospitalization and mortality. The rs1799964 T > C variant of the TNF gene has been associated with various clinical outcomes in different inflammatory and infectious conditions. The increased TNF expression observed in the CC genotype [40] may exacerbate the inflammatory response, worsening COVID-19. In studies on SARS-CoV, Wang et al. (2008) [41] reported that the CC genotype was associated with severe complications, while the CT genotype was identified as a potential protective factor against SARS-CoV infection. However, Li et al. (2022) [42] found an association between the C allele and a mild COVID-19 phenotype, suggesting that gene expression may play a complex role, with significant variations among different coronavirus infections. The rs1800630 C > A variant was significantly associated with a higher susceptibility to severe COVID-19 outcomes, such as hospitalization, with the AA homozygous showing an elevated risk and AA + AC compared to CC also suggesting a potential genetic predisposition to hospitalization. The A allele has been consistently linked to increased TNF-α expression, as reviewed by Smith and Humphries (2009) [43] and Bank et al. (2014) [37]. Additionally, the A allele was associated with a higher promoter affinity for NF-κB, which may influence gene expression [44]. Clinically, the A allele is related to an increased risk of acute respiratory distress syndrome (ARDS) [45]. In contrast, the CC genotype has been associated with severe outcomes in patients with Severe Fever with Thrombocytopenia Syndrome (SFTS) [46]. Non-hospitalized individuals were notably younger compared to hospitalized ones. This may be owing to the fact that older individuals are more susceptible to both infection and severe COVID-19 cases due to a higher likelihood of having multiple comorbidities [47], in addition to the already observed increase in exhausted T-cell response and a decrease in naive T-helper lymphocytes, as well as a reduction in pro-inflammatory cytokines associated with age [48]. Male individuals were the majority in the hospitalized group, as male sex tends to be associated with greater severity and mortality from COVID-19, while young women exhibited an overrepresentation of pathways related to T- and B-cell activation. This suggests that women may have a more robust immune activation profile. Additionally, males tend to exhibit higher levels of inflammatory cytokines and a lower percentage of lymphocytes during SARS-CoV-2 infection and recovery [49]. The most common comorbidities found in this sample are those already been identified as associated with severe COVID-19 cases, including hypertension, diabetes, cardiovascular diseases, obesity, and, less commonly, respiratory diseases, liver diseases, kidney diseases, and neoplasms [47,50,51,52,53]. In this context, the association between genetic variants in the IL1B, IL6, and TNF genes and the severity of COVID-19 was explored, with particular emphasis on outcomes such as hospitalization, the need for ICU admission, and mortality.
This study has limitations that should be considered. The relatively small sample size may have limited the statistical power to detect modest associations, and the lack of longitudinal follow-up precludes assessment of late clinical outcomes such as reinfections or post-acute complications. These factors may limit the generalizability of the findings. Furthermore, serum cytokine levels were not measured, hindering direct biochemical genotype–phenotype correlation. Although prior experimental studies support the functional relevance of the investigated variants [32,37,38,39,40,41,42,43,44], future research should integrate genetic analysis with cytokine profiling (e.g., ELISA or bead-based multiplex assays) to elucidate the mechanistic link between host genetics, immune modulation, and disease severity. While the lack of cytokine profiling is a limitation in this study, it is important to consider that genotype–phenotype correlations may be observed only in particular situations. Functional variants in genes such as IL6, IL1B, or TNF may only exert measurable effects when their downstream pathways are actively engaged—for instance, during an inflammatory response [42,43,44,45,46,47,48,49,50,51,52,53]. In individuals with low baseline cytokine levels, the same risk alleles may not manifest any detectable phenotypic impact. The context-dependent nature of genotype effects means that they are not fixed but continue to be modulated as a result of immunological state or environmental factors [54]. In the case of COVID-19, in which dysregulated cytokine release is a hallmark of severe disease, genetic variants that promote pro-inflammatory signaling might intensify this response, contributing to worse outcomes [54]. Without protein-level data at the moment of the infection, however, it is not possible to determine whether the observed genetic associations reflect true functional consequences or are conditional upon unmeasured immune activation. Future studies should incorporate cytokine quantification to disentangle constitutive versus context-specific genetic effects. Additionally, common genetic variants, such as those we analyzed in this study, could have a small effect size, and this could make it difficult to identify this association [53].
Data on self-reported ancestry and genomic population structure were not available. This study was conducted in a cohort from southern Brazil, characterized by a genetically admixed background [28]. Population stratification and genetic ancestry are important considerations in genomic studies since differences in allele frequencies across ancestral groups can confound associations if not properly accounted for [55]. Including individuals from diverse ancestral backgrounds increases the chances of identifying new associations and can reveal population-specific effects and motivate future replication efforts in those groups [55,56]. Importantly, this case–control study did not adjust for population stratification, which is critical in admixed populations, such as those in Latin America. In such settings, linkage disequilibrium (LD) patterns vary as a function of ancestral contributions: populations with higher Native American ancestry exhibit slower LD decay, whereas African-derived populations show more rapid decay [57]. This heterogeneity, shaped by historical admixture (e.g., the trans-Atlantic slave trade and post-colonial migrations), influences both haplotype structure and allele frequency distributions, potentially leading to spurious associations in the absence of appropriate correction. Future studies should incorporate ancestry-informative markers or genome-wide principal components to account for population substructure and increase the robustness and replicability of genetic findings.
Additionally, many candidate gene associations—especially those identified in small, single-center cohorts—have failed replication in broader or more diverse populations [58,59]. This often happens due to several factors, including differences in LD structure, environmental exposures, phenotype definitions, or analytical approaches across studies [60,61]. Brazil is a high admixture country [28], and LD structure and allelic frequencies, as well as gene–environment interactions, might make different such associations. Failure to replicate the results does not necessarily negate a true biological signal but may reflect context-dependent effects or unmeasured confounding. This is particularly relevant in populations such as Brazil, where unique population-specific genetic variants—such as loss-of-function variants identified by the DNA do Brasil initiative—may act as unmeasured confounders influencing disease susceptibility and outcomes [28]. Therefore, validation in larger, ancestrally diverse cohorts with harmonized protocols is essential to better understand genetic associations. Key strengths of this study include the cohort of 500 unvaccinated individuals from the pre-vaccine era, the use of national epidemiological registry data (SIVEP-Gripe), and the focus on functional variants in central inflammatory pathways. Multivariate models adjusted for major confounders further strengthen the associations found.
5. Conclusions
This study identified significant associations between genetic variants in the IL1B, IL6, and TNF genes and the severity of clinical outcomes in COVID-19 patients. Specific genotypes, such as the AG of rs16944 in the IL1B gene and the CT of rs1799964 in the TNF gene, were associated with an increased risk of hospitalization, ICU admission, and mortality. These findings suggest that these genetic variants may help to develop strategies of future markers for identifying individuals at higher risk of developing severe forms of the disease. Beyond their relevance to COVID-19, these genetic variants have implications for other viral infections, such as SARS-CoV and severe acute respiratory syndrome, with which they have been associated with severe inflammatory outcomes. The presence of these variants might also influence the immune response to vaccines, affecting both the efficacy and the durability of the protection conferred. As a case–control study conducted in an admixed population with a limited sample size, our findings should be interpreted as hypothesis generating. Replication in larger and more genetically diverse cohorts is required to confirm these associations and to deepen the understanding of their biological implications. In addition, future studies should include extended follow-up periods to better assess the progression of conditions related to genetic variants and the impact of vaccination. Additionally, further studies should explore the interaction of these variants with other genetic and environmental factors to provide a more comprehensive understanding of the genetic influence on the severity of viral infections and vaccine efficacy.
Conceptualization, M.F.F., R.C.S., T.W.K., and F.S.L.V.; methodology, M.F.F., R.C.S., M.F.M.-F., V.C.G., N.A.C., J.A.B.C., T.W.K., and F.S.L.V.; validation, M.F.F., R.C.S., and T.W.K.; formal analysis, M.F.F., R.C.S., T.W.K., and F.S.L.V.; investigation, M.F.F., R.C.S., T.W.K., and F.S.L.V.; resources, M.F.F., R.C.S., T.S.G., L.G.M., J.A.B.C., T.W.K., and F.S.L.V.; data curation, M.F.F., T.S.G., L.G.M., T.W.K., and F.S.L.V.; writing—original draft preparation, M.F.F., T.W.K., and F.S.L.V.; writing—review and editing, M.F.F., R.C.S., M.F.M.-F., V.C.G., T.S.G., L.G.M., N.A.C., J.A.B.C., T.W.K., and F.S.L.V.; supervision, F.S.L.V.; project administration, F.S.L.V.; funding acquisition, F.S.L.V. All authors have read and agreed to the published version of the manuscript.
The study was conducted according to the guidelines of the Brazilian regulations and approved by the Research Ethics Committee of the Hospital de Clínicas de Porto Alegre (HCPA) (protocol code CAAE 30797220.9.0000.5327), and it was approved on 29 April 2020.
The Ethical Review Board waived patient consent on 29 April 2020 due to the risks of infection and the inability to obtain informed consent from hospitalized patients, many of whom died shortly after the sample collection.
The data that support the findings of this study are available from the corresponding author upon reasonable request. Due to ethical and legal restrictions, the data are not publicly available.
The authors would like to thank all the healthcare professionals and scientists who worked tirelessly to address this pandemic threat.
The authors declare no conflicts of interest. The sponsors had no role in the design, execution, interpretation, or writing of the study.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Variants of the IL1B, IL6, and TNF genes.
Gene | rsID | Positions * | Chr | Type | Probe | Allele | MAF | DNA Do Brasil ** | HaploReg and ENCODE |
---|---|---|---|---|---|---|---|---|---|
IL1B | rs4848306 | −3737 | 2 | Intergenic | C__11725735_10 | G/A | A (0.36) | 0.36 | Involved in inflammatory response; proximal enhancer; maximum activity near the IL1B gene. |
rs1143623 | −1473 | 2 | Intergenic | C___1839941_10 | C/G | C (0.29) | 0.28 | Proximal enhancer; involved in the regulation of IL1B expression. | |
rs16944 | −511 | 2 | Intergenic | C___1839943_10 | G/A | G (0.49) | 0.47 | Potential inflammatory regulation; located near the IL1B gene, with potential regulatory impact. | |
rs1143627 | −31 | 2 | Intronic | C___1839944_10 | C/T | T (0.47) | 0.49 | Potential regulator of inflammatory expression; promoter; maximum activity near IL1B gene. | |
rs1800795 | −174 | 7 | Intronic | C___1839697_20 | G/C | C (0.14) | 0.21 | Inflammatory response; promoter; involved in the regulation of IL6 expression. | |
IL6 | rs2069832 | +615 | 7 | Intronic | C__15957646_10 | G/A | A (0.14) | 0.21 | Inflammation modulation; promoter; related to increased expression of IL6. |
rs2069840 | +3437 | 7 | Intronic | C__15804104_10 | C/G | G (0.19) | 0.30 | Involvement in inflammatory response; intronic; associated with increased expression of IL6. | |
rs2069845 | +3331 | 7 | Intronic | C___1839699_10 | G/A | G (0.25) | 0.34 | Involvement in inflammatory response; distal enhancer; related to IL6 expression. | |
TNF | rs1799964 | −1031 | 6 | Intronic | C___7514871_10 | T/C | C (0.22) | 0.16 | Potential regulator of inflammatory diseases; proximal enhancer; involved in regulation of TNF expression. |
rs1800630 | −863 | 6 | Intronic | Customized | C/A | A (0.15) | 0.12 | Regulation of immune response; proximal enhancer; highly relevant for the regulation of TNF expression. | |
rs1799724 | −857 | 6 | Intronic | C__11918223_10 | C/T | T (0.10) | 0.07 | Essential for immune response; enhancer; involved in the regulation of TNF. | |
rs1800629 | −308 | 6 | Intronic | C___7514879_10 | G/A | A (0.09) | 0.07 | Regulation of inflammatory response; intronic; involved in the modulation of TNF expression. | |
rs361525 | −238 | 6 | Intronic | C___2215707_10 | G/A | A (0.06) | 0.04 | Promoter; involved in the regulation of TNF; essential for immune response. |
Chr = chromosome; MAF = minor allele frequency; ** DNA do Brasil Project Variant Browser/ABraOM: Brazilian genomic variants (MAF). * Nucleotide positions are relative to the first ATG codon at exon 1, where the adenine is assigned position +1.
Demographic and clinical characteristics of COVID-19 patients (hospitalization, ICU, and death) [ n (%)].
Non-Hospitalized (n = 235) | Hospitalized | p | Non-ICU | ICU | p | Discharged from the Hospital | Death | p | |
---|---|---|---|---|---|---|---|---|---|
Age | |||||||||
n | 235 | 231 | 342 | 108 | 337 | 125 | |||
Years | 38 | 65 | 47 | 64.5 | 44 | 72 | |||
Total | 52 [36.0; 69.0] | <0.01 * | 51.50 (36.75; 63.50) | <0.01 * | 52.50 [29.0; 61.75] | <0.01 * | |||
Sex | |||||||||
Female | 139 (56.5) | 107 (43.5) | 186 (78.2) | 52 (21.8) | 188 (77.0) | 56 (23.0) | |||
Male | 96 (43.6) | 124 (56.4) | 156 (73.6) | 56 (26.4) | 149 (68.3) | 69 (31.7) | |||
Total | 235 (50.4) | 231 (49.6) | 0.01 ‡ | 342 (76.0) | 108 (24.0) | 0.26 ‡ | 337 (72.9) | 125 (27.1) | 0.04 ‡ |
Comorbidities | |||||||||
Heart conditions | |||||||||
No | 158 (57.0) | 119 (43.0) | 214 (79.3) | 56 (20.7) | 218 (79.3) | 57 (20.7) | |||
Yes | 24 (18.7) | 104 (81.3) | 79 (62.7) | 47 (37.3) | 67 (52.3) | 61 (47.7) | |||
Total | 182 (44.9) | 231 (55.1) | <0.01 ‡ | 293 (74.0) | 103 (26.0) | <0.01 ‡ | 285 (70.7) | 118 (29.3) | <0.01 ‡ |
Chronic lung diseases | |||||||||
No | 175 (47.3) | 195 (52.7) | 270 (74.8) | 91 (25.2) | 266 (72.3) | 102 (27.7) | |||
Yes | 6 (20.7) | 23 (79.3) | 19 (65.5) | 10 (34.5) | 15 (51.7) | 14 (48.3) | |||
Total | 181 (45.4) | 218 (54.6) | 0.01 ‡ | 289 (74.1) | 101 (25.9) | 0.27 ‡ | 281 (70.8) | 116 (29.2) | 0.02 ‡ |
Obesity | |||||||||
No | 173 (49.3) | 178 (50.7) | 266 (77.6) | 77 (22.4) | 255 (73.1) | 94 (26.9) | |||
Yes | 8 (18.2) | 36 (81.8) | 19 (44.2) | 24 (55.8) | 23 (52.3) | 21 (47.7) | |||
Total | 181 (45.8) | 214 (54.2) | <0.01 ‡ | 285 (73.8) | 101 (26.2) | <0.01 ‡ | 278 (70.7) | 115 (29.3) | 0.004 ‡ |
Diabetes mellitus | |||||||||
No | 168 (51.4) | 159 (48.6) | 252 (78.5) | 69 (21.5) | 248 (76.1) | 78 (23.9) | |||
Yes | 15 (19.5) | 62 (80.5) | 39 (53.4) | 34 (46.6) | 36 (47.4) | 40 (52.6) | |||
Total | 183 (45.3) | 221 (54.7) | <0.01 ‡ | 291 (73.9) | 103 (26.1) | <0.01 ‡ | 284 (70.6) | 118 (29.4) | <0.01 ‡ |
Severity § | |||||||||
Admission to critical care (n = 414) | |||||||||
No | 116 (51.8) | ||||||||
Yes | 108 (48.2) | ||||||||
Total | 224 (49.8) | <0.01 ‡ | |||||||
Ventilatory support required (n = 453) | |||||||||
No | 14 (6.2) | ||||||||
Yes | 213 (93.8) | ||||||||
Total | 227 (50.1) | <0.01 ‡ | |||||||
Death (n = 427) | |||||||||
No | 104 (45.4) | ||||||||
Yes | 125 (54.6) | ||||||||
Total | 229 (49.6) | <0.01 ‡ |
[IQR−–IQR+]; * Mann–Whitney test; ‡ chi-square test. § The present analysis was restricted to hospitalized individuals.
Multivariate logistic regression analysis for variants of IL1B, IL6, and TNF genes—hospitalization *.
Gene | Variant | Genotype | Odds Ratio (OR) | 95% CI | p |
---|---|---|---|---|---|
IL1B | rs16944 | AA | 1.015 | 0.548–1.879 | 0.961 |
AG | 1.976 | 1.218–3.234 | 0.006 | ||
AA + AG vs. GG | 1.620 | 1.035–2.552 | 0.036 | ||
GG + AG vs. AA | 0.697 | 0.400–1.208 | 0.200 | ||
rs1143627 | TT | 0.958 | 0.526–1.745 | 0.887 | |
CT | 1.820 | 1.021–3.269 | 0.043 | ||
TT + CT vs. CC | 1.363 | 0.801–2.330 | 0.255 | ||
CC + CT vs. TT | 0.639 | 0.405–1.001 | 0.052 | ||
IL6 | rs2069845 | GG | 1.006 | 0.489–2.080 | 0.987 |
AG | 0.665 | 0.420–1.050 | 0.081 | ||
GG + AG vs. AA | 0.724 | 0.468–1.115 | 0.143 | ||
AA + AG vs. GG | 1.229 | 0.620–2.455 | 0.555 | ||
TNF | rs1799964 | CC | 1.720 | 0.684–4.422 | 0.252 |
CT | 1.970 | 1.220–3.215 | 0.006 | ||
CC + CT vs. TT | 1.927 | 1.226–3.060 | 0.005 | ||
TT + CT vs. CC | 1.364 | 0.554–3.426 | 0.501 | ||
rs1800630 | AA | 2.367 | 1.079–5.326 | 0.034 | |
AC | 1.584 | 0.919–2.681 | 0.101 | ||
AA + AC vs. CC | 1.755 | 1.094–2.842 | 0.021 | ||
CC + AC vs. AA | 2.112 | 0.976–4.687 | 0.061 | ||
rs1799724 | TT | 2.852 | 0.562–21.694 | 0.242 | |
CT | 1.411 | 0.823–2.434 | 0.212 | ||
TT + CT vs. CC | 1.491 | 0.886–2.529 | 0.134 | ||
CC + CT vs. TT | 2.654 | 0.524–20.149 | 0.275 | ||
rs1800629 | AA | 0.350 | 0.027–8.621 | 0.433 | |
AG | 1.054 | 0.606–1.833 | 0.852 | ||
AA + AG vs. GG | 1.020 | 0.590–1.762 | 0.944 | ||
GG + AG vs. AA | 0.347 | 0.027–8.530 | 0.429 |
* All models were adjusted for age and sex. Bold numbers indicate statistical significance.
Multivariate logistic regression analysis for variants of IL1B, IL6, and TNF genes—intensive care unit (ICU) *.
Gene | Variant | Genotype | Odds Ratio (OR) | 95% CI | p |
---|---|---|---|---|---|
IL1B | rs16944 | AA | 0.835 | 0.452–1.559 | 0.567 |
AG | 1.834 | 1.092–3.101 | 0.022 | ||
AA + AG vs. GG | 1.425 | 0.891–2.272 | 0.138 | ||
GG + AG vs. AA | 0.608 | 0.349–1.077 | 0.083 | ||
rs1143627 | TT | 1.085 | 0.587–1.983 | 0.793 | |
CT | 2.082 | 1.123–3.848 | 0.019 | ||
TT + CT vs. CC | 1.526 | 0.873–2.626 | 0.131 | ||
CC + CT vs. TT | 0.669 | 0.418–1.073 | 0.094 | ||
IL6 | rs2069845 | GG | 1.590 | 0.715–3.852 | 0.276 |
AG | 0.653 | 0.400–1.061 | 0.086 | ||
GG + AG vs. AA | 0.776 | 0.485–1.233 | 0.285 | ||
AA + AG vs. GG | 1.979 | 0.934–4.630 | 0.091 | ||
TNF | rs1799964 | CC | 1.844 | 0.687–5.908 | 0.256 |
CT | 1.252 | 0.762–2.087 | 0.381 | ||
CC + CT vs. TT | 1.328 | 0.827–2.158 | 0.245 | ||
TT + CT vs. CC | 1.707 | 0.646–5.404 | 0.315 | ||
rs1800630 | AA | 1.474 | 0.651–3.689 | 0.376 | |
AC | 1.083 | 0.633–1.893 | 0.773 | ||
AA + AC vs. CC | 1.175 | 0.722–1.940 | 0.521 | ||
CC + AC vs. AA | 1.442 | 0.646–3.572 | 0.396 | ||
rs1799724 | TT | … | … | … | |
CT | … | … | … | ||
TT + CT vs. CC | 1.218 | 0.702–2.179 | 0.493 | ||
CC + CT vs. TT | … | … | … | ||
rs1800629 | AA | 0.568 | 0.025–6.342 | 0.654 | |
AG | 0.541 | 0.262–1.037 | 0.077 | ||
AA + AG vs. GG | 0.542 | 0.269–1.024 | 0.710 | ||
GG + AG vs. AA | 0.621 | 0.028–6.931 | 0.706 |
* All models were adjusted for age and sex. Note : “…” indicates that estimates could not be calculated due to an insufficient number of observations. Bold numbers indicate statistical significance.
Multivariate logistic regression analysis for variants of IL1B, IL6, and TNF genes—death *.
Gene | Variant | Genotype | Odds Ratio (OR) | 95% CI | p |
---|---|---|---|---|---|
IL1B | rs16944 | AA | 0.704 | 0.365–1.362 | 0.295 |
AG | 1.274 | 0.750–2.168 | 0.370 | ||
AA + AG vs. GG | 1.065 | 0.652–1.731 | 0.799 | ||
GG + AG vs. AA | 0.617 | 0.342–1.123 | 0.110 | ||
rs1143627 | TT | 1.357 | 0.709–2.590 | 0.354 | |
CT | 1.615 | 0.860–3.028 | 0.134 | ||
TT + CT vs. CC | 1.493 | 0.833–2.653 | 0.173 | ||
CC + CT vs. TT | 0.979 | 0.601–1.607 | 0.933 | ||
IL6 | rs2069845 | GG | 1.266 | 0.579–2.882 | 0.563 |
AG | 0.930 | 0.562–1.537 | 0.776 | ||
GG + AG vs. AA | 0.990 | 0.614–1.594 | 0.966 | ||
AA + AG vs. GG | 1.313 | 0.627–2.877 | 0.481 | ||
TNF | rs1799964 | CC | 3.728 | 1.208–14.366 | 0.034 |
CT | 1.452 | 0.863–2.479 | 0.165 | ||
CC + CT vs. TT | 1.658 | 1.008–2.766 | 0.049 | ||
TT + CT vs. CC | 3.279 | 1.081–12.472 | 0.052 | ||
rs1800630 | AA | 1.498 | 0.652–3.668 | 0.355 | |
AC | 1.619 | 0.904–2.973 | 0.112 | ||
AA + AC vs. CC | 1.583 | 0.945–2.700 | 0.086 | ||
CC + AC vs. AA | 1.328 | 0.587–3.205 | 0.509 | ||
rs1799724 | TT | … | … | … | |
CT | … | … | … | ||
TT + CT vs. CC | 1.682 | 0.931–3.141 | 0.092 | ||
CC + CT vs. TT | … | … | … | ||
rs1800629 | AA | 1.386 | 0.112–33.456 | 0.805 | |
AG | 1.033 | 0.550–1.895 | 0.918 | ||
AA + AG vs. GG | 1.047 | 0.566–1.898 | 0.880 | ||
GG + AG vs. AA | 1.379 | 0.111–33.228 | 0.808 |
* All models were adjusted for age and sex. Note : “…” indicates that estimates could not be calculated due to an insufficient number of observations. Bold numbers indicate statistical significance.
Supplementary Materials
The following supporting information can be downloaded at
1. Negi, K.; Agarwal, M.; Pahuja, I.; Bhardwaj, B.; Rawat, M.; Bhaskar, A.; Dwivedi, V.P. Combating the Challenges of COVID-19 Pandemic: Insights into Molecular Mechanisms, Immune Responses and Therapeutics against SARS-CoV-2. Oxf. Open Immunol.; 2023; 4, iqad001. [DOI: https://dx.doi.org/10.1093/oxfimm/iqad001]
2. Tort, S.; Ciapponi, A. What Is the Accuracy of World Health Organization (WHO)-Specified and Related COVID-19 Symptoms for the Diagnosis of COVID-19?. Cochrane Clin. Answ.; 2022; Available online: https://www.cochranelibrary.com/cca/doi/10.1002/cca.3998/full# (accessed on 23 October 2024). [DOI: https://dx.doi.org/10.1002/cca.3998]
3. Struyf, T.; Deeks, J.J.; Dinnes, J.; Takwoingi, Y.; Davenport, C.; Leeflang, M.M.; Spijker, R.; Hooft, L.; Emperador, D.; Domen, J.
4. Zhang, Q.; Matuozzo, D.; Le Pen, J.; Lee, D.; Moens, L.; Asano, T.; Bohlen, J.; Liu, Z.; Moncada-Velez, M.; Kendir-Demirkol, Y.
5. Costela-Ruiz, V.J.; Illescas-Montes, R.; Puerta-Puerta, J.M.; Ruiz, C.; Melguizo-Rodríguez, L. SARS-CoV-2 Infection: The Role of Cytokines in COVID-19 Disease. Cytokine Growth Factor. Rev.; 2020; 54, pp. 62-75. [DOI: https://dx.doi.org/10.1016/j.cytogfr.2020.06.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32513566]
6. Fara, A.; Mitrev, Z.; Rosalia, R.A.; Assas, B.M. Cytokine Storm and COVID-19: A Chronicle of pro-Inflammatory Cytokines. Open Biol.; 2020; 10, 200160. [DOI: https://dx.doi.org/10.1098/rsob.200160] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32961074]
7. Hu, B.; Huang, S.; Yin, L. The Cytokine Storm and COVID-19. J. Med. Virol.; 2021; 93, pp. 250-256. [DOI: https://dx.doi.org/10.1002/jmv.26232]
8. Ben Moftah, M.; Eswayah, A. Intricate Relationship between SARS-CoV-2–Induced Shedding and Cytokine Storm Generation: A Signaling Inflammatory Pathway Augmenting COVID-19. Health Sci. Rev.; 2022; 2, 100011. [DOI: https://dx.doi.org/10.1016/j.hsr.2021.100011]
9. Manik, M.; Singh, R.K. Role of Toll-like Receptors in Modulation of Cytokine Storm Signaling in SARS-CoV-2-induced COVID-19. J. Med. Virol.; 2022; 94, pp. 869-877. [DOI: https://dx.doi.org/10.1002/jmv.27405]
10. Wong, L.-Y.R.; Perlman, S. Immune Dysregulation and Immunopathology Induced by SARS-CoV-2 and Related Coronaviruses —Are We Our Own Worst Enemy?. Nat. Rev. Immunol.; 2022; 22, pp. 47-56. [DOI: https://dx.doi.org/10.1038/s41577-021-00656-2]
11. Gomes, S.M.R.; Brito, A.C.D.S.; Manfro, W.F.P.; Ribeiro-Alves, M.; Ribeiro, R.S.D.A.; Da Cal, M.S.; Lisboa, V.D.C.; Abreu, D.P.B.D.; Castilho, L.D.R.; Porto, L.C.D.M.S.
12. Lechowicz, K.; Drożdżal, S.; Machaj, F.; Rosik, J.; Szostak, B.; Zegan-Barańska, M.; Biernawska, J.; Dabrowski, W.; Rotter, I.; Kotfis, K. COVID-19: The Potential Treatment of Pulmonary Fibrosis Associated with SARS-CoV-2 Infection. JCM; 2020; 9, 1917. [DOI: https://dx.doi.org/10.3390/jcm9061917]
13. Bell, L.C.K.; Meydan, C.; Kim, J.; Foox, J.; Butler, D.; Mason, C.E.; Shapira, S.D.; Noursadeghi, M.; Pollara, G. Transcriptional Response Modules Characterize IL-1β and IL-6 Activity in COVID-19. iScience; 2021; 24, 101896. [DOI: https://dx.doi.org/10.1016/j.isci.2020.101896] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33319166]
14. Hu, H.; Pan, H.; Li, R.; He, K.; Zhang, H.; Liu, L. Increased Circulating Cytokines Have a Role in COVID-19 Severity and Death With a More Pronounced Effect in Males: A Systematic Review and Meta-Analysis. Front. Pharmacol.; 2022; 13, 802228. [DOI: https://dx.doi.org/10.3389/fphar.2022.802228]
15. Santos, A.A.C.D.; Rodrigues, L.E.; Alecrim-Zeza, A.L.; De Araújo Ferreira, L.; Trettel, C.D.S.; Gimenes, G.M.; Silva, A.F.D.; Sousa-Filho, C.P.B.; Serdan, T.D.A.; Levada-Pires, A.C.
16. Aziz, M.; Fatima, R.; Assaly, R. Elevated Interleukin-6 and Severe COVID-19: A Meta-analysis. J. Med. Virol.; 2020; 92, pp. 2283-2285. [DOI: https://dx.doi.org/10.1002/jmv.25948]
17. Galván-Román, J.M.; Rodríguez-García, S.C.; Roy-Vallejo, E.; Marcos-Jiménez, A.; Sánchez-Alonso, S.; Fernández-Díaz, C.; Alcaraz-Serna, A.; Mateu-Albero, T.; Rodríguez-Cortes, P.; Sánchez-Cerrillo, I.
18. Yin, J.-X.; Agbana, Y.L.; Sun, Z.-S.; Fei, S.-W.; Zhao, H.-Q.; Zhou, X.-N.; Chen, J.-H.; Kassegne, K. Increased Interleukin-6 Is Associated with Long COVID-19: A Systematic Review and Meta-Analysis. Infect. Dis. Poverty; 2023; 12, 43. [DOI: https://dx.doi.org/10.1186/s40249-023-01086-z]
19. Lucas, C.; Wong, P.; Klein, J.; Castro, T.B.R.; Silva, J.; Sundaram, M.; Ellingson, M.K.; Mao, T.; Oh, J.E.; Israelow, B.
20. Mortaz, E.; Tabarsi, P.; Jamaati, H.; Dalil Roofchayee, N.; Dezfuli, N.K.; Hashemian, S.M.; Moniri, A.; Marjani, M.; Malekmohammad, M.; Mansouri, D.
21. Ablamunits, V.; Lepsy, C. Blocking TNF Signaling May Save Lives in COVID-19 Infection. Mol. Biol. Rep.; 2022; 49, pp. 2303-2309. [DOI: https://dx.doi.org/10.1007/s11033-022-07166-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35076845]
22. Qudus, M.S.; Tian, M.; Sirajuddin, S.; Liu, S.; Afaq, U.; Wali, M.; Liu, J.; Pan, P.; Luo, Z.; Zhang, Q.
23. Verma, S.; Verma, S.; Khan, F.H.; Siddiqi, Z.; Raza, S.T.; 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, pp. 5-11. [DOI: https://dx.doi.org/10.1111/iji.12605] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36323530]
24. Rashid, P.M.A.; Salih, G.F. Genetic Polymorphism of Interleukin-6 in Asymptomatic and ICU-Admitted COVID-19 Patients in Sulaymaniyah Province, Kurdistan Region of Iraq. Genet. Test. Mol. Biomark.; 2024; 28, pp. 297-303. [DOI: https://dx.doi.org/10.1089/gtmb.2023.0304]
25. Deng, X.; Tang, K.; Wang, Z.; He, S.; Luo, Z. Impacts of Inflammatory Cytokines Variants on Systemic Inflammatory Profile and COVID-19 Severity. J. Epidemiol. Glob. Health; 2024; 14, pp. 363-378. [DOI: https://dx.doi.org/10.1007/s44197-024-00204-w]
26. Ward, L.D.; Kellis, M. HaploReg: A Resource for Exploring Chromatin States, Conservation, and Regulatory Motif Alterations within Sets of Genetically Linked Variants. Nucleic Acids Res.; 2012; 40, pp. D930-D934. [DOI: https://dx.doi.org/10.1093/nar/gkr917]
27. The ENCODE Project Consortium. An Integrated Encyclopedia of DNA Elements in the Human Genome. Nature; 2012; 489, pp. 57-74. [DOI: https://dx.doi.org/10.1038/nature11247]
28. Nunes, K.; Kimura, L.; de Sena, L.S.; Krieger, J.E.; Varela, M.C.; Camargo, L.M.A.; Ferreira, R.G.M.; Krieger, H.; Bortolini, M.C.; Mill, J.G.
29. Secretaria de Atenção Primária à Saúde (SAPS). Protocolo de Manejo Clínico do Coronavírus (COVID-19) na Atenção Primária à Saúde (v9); Ministério da Saúde: Brasília, DF, Brasil, 2020.
30. Sbruzzi, R.C.; Feira, M.F.; Cadore, N.A.; Giudicelli, G.C.; Kowalski, T.W.; Gregianini, T.S.; Chies, J.A.B.; Vianna, F.S.L. An Efficient Extraction Method Allowing the Genetic Evaluation of Host DNA from Samples Collected for Virus Infection Diagnosis in Viral Transport Medium. Biopreservation Biobank.; 2024; 22, pp. 166-173. [DOI: https://dx.doi.org/10.1089/bio.2022.0188]
31. Balzanelli, M.; Distratis, P.; Lazzaro, R.; Pham, V.; Tran, T.; Dipalma, G.; Bianco, A.; Serlenga, E.; Aityan, S.; Pierangeli, V.
32. Behzadi, P.; Sameer, A.S.; Nissar, S.; Banday, M.Z.; Gajdács, M.; García-Perdomo, H.A.; Akhtar, K.; Pinheiro, M.; Magnusson, P.; Sarshar, M.
33. Xie, Z.-K.; Huang, Q.-P.; Huang, J.; Xie, Z.-F. Association between the IL1B, IL1RN Polymorphisms and COPD Risk: A Meta-Analysis. Sci. Rep.; 2014; 4, 6202. [DOI: https://dx.doi.org/10.1038/srep06202] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25174605]
34. Jiménez-Sousa, M.Á.; Medrano, L.M.; Liu, P.; Almansa, R.; Fernández-Rodríguez, A.; Gómez-Sánchez, E.; Rico, L.; Heredia-Rodríguez, M.; Gómez-Pesquera, E.; Tamayo, E.
35. Galán-Huerta, K.A.; Zamora-Márquez, M.A.; Flores-Pérez, R.O.; Bocanegra-Ibarias, P.; Salas-Treviño, D.; Rivas-Estilla, A.M.G.; Flores-Treviño, S.; Lozano-Sepúlveda, S.A.; Martínez-Acuña, N.; Camacho-Ortiz, A.
36. Tanimine, N.; Takei, D.; Tsukiyama, N.; Yoshinaka, H.; Takemoto, Y.; Tanaka, Y.; Kobayashi, T.; Tanabe, K.; Ishikawa, N.; Kitahara, Y.
37. Bank, S.; Andersen, P.S.; Burisch, J.; Pedersen, N.; Roug, S.; Galsgaard, J.; Turino, S.Y.; Brodersen, J.B.; Rashid, S.; Rasmussen, B.K.
38. Chen, H.; Wilkins, L.M.; Aziz, N.; Cannings, C.; Wyllie, D.H.; Bingle, C.; Rogus, J.; Beck, J.D.; Offenbacher, S.; Cork, M.J.
39. Zhang, G.; Zhou, B.; Li, S.; Yue, J.; Yang, H.; Wen, Y.; Zhan, S.; Wang, W.; Liao, M.; Zhang, M.
40. Nourian, M.; Chaleshi, V.; Pishkar, L.; Azimzadeh, P.; Baradaran Ghavami, S.; Balaii, H.; Alinaghi, S.; Shahrokh, S.; Asadzadeh Aghdaei, H.; Zali, M.R. Evaluation of Tumor Necrosis Factor (TNF)-α mRNA Expression Level and the Rs1799964 Polymorphism of the TNF-α Gene in Peripheral Mononuclear Cells of Patients with Inflammatory Bowel Diseases. Biomed. Rep.; 2017; 6, pp. 698-702. [DOI: https://dx.doi.org/10.3892/br.2017.908]
41. Wang, S.; Wei, M.; Han, Y.; Zhang, K.; He, L.; Yang, Z.; Su, B.; Zhang, Z.; Hu, Y.; Hui, W. Roles of TNF-α Gene Polymorphisms in the Occurrence and Progress of SARS-Cov Infection: A Case-Control Study. BMC Infect. Dis.; 2008; 8, 27. [DOI: https://dx.doi.org/10.1186/1471-2334-8-27]
42. Li, P.; Ke, Y.; Shen, W.; Shi, S.; Wang, Y.; Lin, K.; Guo, X.; Wang, C.; Zhang, Y.; Zhao, Z. Targeted Screening of Genetic Associations with COVID-19 Susceptibility and Severity. Front. Genet.; 2022; 13, 1073880. [DOI: https://dx.doi.org/10.3389/fgene.2022.1073880] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36531218]
43. Smith, A.J.P.; Humphries, S.E. Cytokine and Cytokine Receptor Gene Polymorphisms and Their Functionality. Cytokine Growth Factor. Rev.; 2009; 20, pp. 43-59. [DOI: https://dx.doi.org/10.1016/j.cytogfr.2008.11.006]
44. Skoog, T.; Hamsten, A.; Eriksson, P. Allele-Specific Chromatin Remodeling of the Tumor Necrosis Factor-α Promoter. Biochem. Biophys. Res. Commun.; 2006; 351, pp. 777-783. [DOI: https://dx.doi.org/10.1016/j.bbrc.2006.10.114]
45. Azevedo, Z.M.; Moore, D.B.; Lima, F.C.; Cardoso, C.C.; Bougleux, R.; Matos, G.I.; Luz, R.A.; Xavier-Elsas, P.; Sampaio, E.P.; Gaspar-Elsas, M.I.
46. Xing, B.; Li, X.-K.; Zhang, S.-F.; Lu, Q.-B.; Du, J.; Zhang, P.-H.; Yang, Z.-D.; Cui, N.; Guo, C.-T.; Cao, W.-C.
47. Waqar, M.A.; Riaz, T.; Munir, M.; Ashfaq, A.; Zulfiqar, H.; Sandhu, M.; Mahmood, A.; Zahir Durrani, F.; Faheem, S.; Tayyab, S. COVID-19 and Various Comorbidities: A Review: COVID-19 and Various Comorbidities. Pak. J. Health Sci.; 2023; 4, pp. 2-7. [DOI: https://dx.doi.org/10.54393/pjhs.v4i01.485]
48. Aragon, L.; Iribarren-López, A.; Alberro, A.; Iparraguirre, L.; Von Wichmann, M.; Marimon, J.M.; Saiz-Calderon, N.; Agudo, J.; Gálvez, M.I.; Cipitria, M.C.
49. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X.
50. Chatterjee, S.; Nalla, L.V.; Sharma, M.; Sharma, N.; Singh, A.A.; Malim, F.M.; Ghatage, M.; Mukarram, M.; Pawar, A.; Parihar, N.
51. Fedorchenko, Y.; Zimba, O. Comorbidities in the COVID-19 Pandemic: Scopus-Based Bibliometric Analysis. J. Korean Med. Sci.; 2023; 38, e93. [DOI: https://dx.doi.org/10.3346/jkms.2023.38.e93]
52. Jonnalagadda, S.; Prathyusha, Y.; Kumar, V.S.; Aparna, Y. Systematic Review and Meta Analysis on Symptoms and Comorbidities Associated with COVID-19 Disease. J. Pure Appl. Microbiol.; 2022; 16, pp. 2337-2351. [DOI: https://dx.doi.org/10.22207/JPAM.16.4.50]
53. Cobat, A.; Zhang, Q. COVID Human Genetic Effort Abel, L.; Casanova, J.L.; Fellay, J. Human Genomics of COVID-19 Pneumonia: Contributions of Rare and Common Variants. Annu. Rev. Biomed. Data Sci.; 2023; 6, pp. 465-486. [DOI: https://dx.doi.org/10.1146/annurev-biodatasci-020222-021705] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37196358]
54. Hiti, L.; Markovič, T.; Lainscak, M.; Farkaš Lainščak, J.; Pal, E.; Mlinarič-Raščan, I. The Immunopathogenesis of a Cytokine Storm: The Key Mechanisms Underlying Severe COVID-19. Cytokine Growth Factor. Rev.; 2025; 82, pp. 1-17. [DOI: https://dx.doi.org/10.1016/j.cytogfr.2024.12.003]
55. Caliebe, A.; Tekola-Ayele, F.; Darst, B.F.; Wang, X.; Song, Y.E.; Gui, J.; Sebro, R.A.; Balding, D.J.; Saad, M.; Dubé, M.
56. Taylor, D.J.; Chhetri, S.B.; Tassia, M.G.; Biddanda, A.; Yan, S.M.; Wojcik, G.L.; Battle, A.; McCoy, R.C. Sources of Gene Expression Variation in a Globally Diverse Human Cohort. Nature; 2024; 632, pp. 122-130. [DOI: https://dx.doi.org/10.1038/s41586-024-07708-2]
57. Da Cruz, P.R.S.; Ananina, G.; Secolin, R.; Gil-da-Silva-Lopes, V.L.; Lima, C.S.P.; De França, P.H.C.; Donatti, A.; Lourenço, G.J.; De Araujo, T.K.; Simioni, M.
58. Roberts, G.H.L.; Partha, R.; Rhead, B.; Knight, S.C.; Park, D.S.; Coignet, M.V.; Zhang, M.; Berkowitz, N.; Turrisini, D.A.; Gaddis, M.
59. Povysil, G.; Butler-Laporte, G.; Shang, N.; Wang, C.; Khan, A.; Alaamery, M.; Nakanishi, T.; Zhou, S.; Forgetta, V.; Eveleigh, R.J.M.
60. Tkachenko, A.A.; Changalidis, A.I.; Maksiutenko, E.M.; Nasykhova, Y.A.; Barbitoff, Y.A.; Glotov, A.S. Replication of Known and Identification of Novel Associations in Biobank-Scale Datasets: A Survey Using UK Biobank and FinnGen. Genes; 2024; 15, 931. [DOI: https://dx.doi.org/10.3390/genes15070931]
61. Zou, J.; Zhou, J.; Faller, S.; Brown, R.P.; Sankararaman, S.S.; Eskin, E. Accurate Modeling of Replication Rates in Genome-Wide Association Studies by Accounting for Winner’s Curse and Study-Specific Heterogeneity. G3 Genes|Genomes|Genet.; 2022; 12, jkac261. [DOI: https://dx.doi.org/10.1093/g3journal/jkac261]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details






1 Postgraduate Program in Genetics and Molecular Biology (PPGBM), Department of Genetics, Institute of Biosciences, Universidade Federal do Rio Grande do Sul, Porto Alegre 91501-970, RS, Brazil; [email protected] (M.F.F.); [email protected] (R.C.S.); [email protected] (M.F.M.-F.); [email protected] (N.A.C.); [email protected] (J.A.B.C.); [email protected] (T.W.K.), Laboratory of Genomic Medicine, Center of Experimental Research, Hospital de Clínicas de Porto Alegre (HCPA), Porto Alegre 90035-903, RS, Brazil; [email protected], Laboratory of Immunobiology and Immunogenetics, Department of Genetics, Institute of Biosciences, Universidade Federal do Rio Grande do Sul, Porto Alegre 91501-970, RS, Brazil
2 Laboratory of Genomic Medicine, Center of Experimental Research, Hospital de Clínicas de Porto Alegre (HCPA), Porto Alegre 90035-903, RS, Brazil; [email protected], Laboratory of Immunobiology and Immunogenetics, Department of Genetics, Institute of Biosciences, Universidade Federal do Rio Grande do Sul, Porto Alegre 91501-970, RS, Brazil
3 Laboratório Central do Centro Estadual de Vigilância em Saúde da do Rio Grande do Sul (LACEN), Secretaria da Saúde do estado do Rio Grande do Sul (LACEN/CEVS/SES-RS), Porto Alegre 90620-000, RS, Brazil; [email protected]
4 Centro Estadual de Vigilância em Saúde da Secretaria de Saúde do estado do Rio Grande do Sul (CEVS/SES-RS), Porto Alegre 90620-000, RS, Brazil; [email protected]
5 Postgraduate Program in Genetics and Molecular Biology (PPGBM), Department of Genetics, Institute of Biosciences, Universidade Federal do Rio Grande do Sul, Porto Alegre 91501-970, RS, Brazil; [email protected] (M.F.F.); [email protected] (R.C.S.); [email protected] (M.F.M.-F.); [email protected] (N.A.C.); [email protected] (J.A.B.C.); [email protected] (T.W.K.), Laboratory of Immunobiology and Immunogenetics, Department of Genetics, Institute of Biosciences, Universidade Federal do Rio Grande do Sul, Porto Alegre 91501-970, RS, Brazil
6 Postgraduate Program in Genetics and Molecular Biology (PPGBM), Department of Genetics, Institute of Biosciences, Universidade Federal do Rio Grande do Sul, Porto Alegre 91501-970, RS, Brazil; [email protected] (M.F.F.); [email protected] (R.C.S.); [email protected] (M.F.M.-F.); [email protected] (N.A.C.); [email protected] (J.A.B.C.); [email protected] (T.W.K.), Laboratory of Genomic Medicine, Center of Experimental Research, Hospital de Clínicas de Porto Alegre (HCPA), Porto Alegre 90035-903, RS, Brazil; [email protected], Laboratory Genetics Unit, Medical Genetics Service, Hospital de Clínicas de Porto Alegre (HCPA), Porto Alegre 90035-903, RS, Brazil