Keywords
covid-19, single nucleotide polymorphisms, Covid-19 respiratory syndrome, SERPINE1, IFNAR2
This article is included in the Genomics and Genetics gateway.
This article is included in the Coronavirus (COVID-19) collection.
Coronavirus Disease 2019 (COVID-19), caused by SARS-CoV-2, is a global pandemic. Severe cases may develop post-COVID-19 respiratory syndrome, characterized by symptoms persisting ≥4 weeks. In this study, we aim to identify the association of post-COVID-19 respiratory syndromes with variants in the SERPINE1 and IFNAR2 genes which are responsible for regulating fibrinolysis and immune response against viral infection, respectively.
We recalled 85 COVID-19 positive subjects who were originally admitted to Persahabatan Central General Hospital, Jakarta, in 2020—2021. We performed TaqMan genotyping on clinically relevant variants at rs1051393 and rs2229207 in IFNAR2 and rs6092 in SERPINE1 and conducted statistical analysis to identify the association of these variants with the clinical manifestations of post-COVID-19 respiratory syndrome.
We found that patients with post-COVID-19 respiratory syndrome had greater disease severity at admission (p = 0.02) and during hospitalization (p = 0.004), along with longer hospital stays (p = 0.0001), possibly due to poorer radiological findings at admission (p = 0.013) and discharge (p = 0.001). Variants in SERPINE1 and IFNAR2 showed significant associations: rs1051393 (IFNAR2) with fever (p = 0.037), and rs2229207 with chest pain (p = 0.047), higher potassium (p = 0.001), chloride (p = 0.047), and pCO2 (p = 0.030). Additionally, rs6092 (SERPINE1) was linked to increased basophils (p = 0.027) and eGFR (p = 0.036). After Benjamini–Hochberg correction, only disease severity during hospitalization (padj = 0.021), radiological findings at admission (padj = 0.052) and discharge (padj = 0.008), hospitalization duration (padj = 0.00016), and the association of rs2229207 with potassium levels (padj = 0.003) remained significant.
Our findings suggested the potential contribution of genetic predispositions in the SERPINE1 and IFNAR2 genes to the emergence of post-COVID-19 respiratory symptoms, providing additional consideration in treatment management to improve the outcomes of patients positive for COVID-19.
covid-19, single nucleotide polymorphisms, Covid-19 respiratory syndrome, SERPINE1, IFNAR2
In this revised version of the manuscript, we have addressed the comments and suggestions provided by the reviewers. The revisions include improvements in language and clarity, as well as the addition of further data analyses. Specifically, we have incorporated additional genetic model analyses, including dominant, codominant, overdominant, and recessive models.
Furthermore, we have conducted multivariate analyses across all study parameters to provide a more comprehensive evaluation. Revisions have also been made to the Introduction and Discussion sections to enhance clarity and strengthen the interpretation of the findings. In addition, the tables presenting statistical results have been revised.
Methodological improvements have also been implemented, including a clearer justification for the use of saliva samples, refinement of the statistical analysis procedures, and more detailed inclusion criteria for study participants to ensure greater comprehensiveness.
See the authors' detailed response to the review by Shrikant Verma
See the authors' detailed response to the review by Abdullah Al Saba
Coronavirus Disease 2019 (COVID-19) was a global pandemic of respiratory illness caused by SARS-CoV-2.1,2 According to the Ministry of Health of the Republic of Indonesia in March 2022, the number of confirmed COVID-19 cases reached 5,891,022, with 154,221 reported deaths, while 5,658,238 patients had recovered.3 However, individuals recovered from COVID-19 may develop post-COVID-19 respiratory syndrome, characterized by persistent respiratory symptoms lasting ≥4 weeks after the initial onset.4 Zhou et al. (2020) described that among 55 COVID-19 patients, 14.5% experienced shortness of breath that persisted for up to three months,5 while another study revealed that 45.5% of 183 patients experienced shortness of breath for up to four weeks.6 Another study reported that among patients with post-COVID-19 respiratory syndrome, 15.5% experienced persistent cough, 11.2% reported shortness of breath, and 3.8% experienced sore throat.7
The emergence of post-COVID-19 symptoms is often preceded by severe COVID-19 symptoms, such as cytokine storms and dysregulated coagulation that may result in excessive blood clotting.7–9 Previous studies have highlighted the association between dysregulated PAI-1, the main inhibitor of fibrinolysis enzymes important for removing clot formation, with severe and critical COVID-19 conditions.10–12 A report by Fricke-Galindo et al. (2022) disclosed that rs6092 in SERPINE1 was associated with higher levels of P-selectin and tPA in subjects carrying GA + AA genotypes in severe COVID-19 patients, which may indicate the contribution of genetic variants to the impaired fibrinolysis in severe COVID-19.11 Similarly, variants in the IFNAR2 have been reported to be associated with significantly high mortality risk in subjects with severe COVID-19, as higher soluble IFNAR2 levels were also observed among the surviving patients.13 However, whether the genetic predisposition in SERPINE and IFNAR2 is associated with post-COVID-19 respiratory symptoms is still unclear and requires further confirmation.
Here, we investigated the association of rs6092 in SERPINE1, and rs1051393 and rs2229207 in IFNAR2 with prolonged respiratory symptoms and blood markers in 85 subjects of COVID-19 in Indonesia. We indicated significant associations of variants in SERPINE1 and IFNAR2, in which rs1051393 in IFNAR2 was significantly associated with fever, while rs2229207 was associated with chest pain, higher potassium and chloride levels in blood analysis, as well as higher pCO2 level. We also note that rs6092 in SERPINE1 was associated with elevated levels of basophils and eGFR in blood analysis. Our findings suggested the potential contribution of genetic predispositions in the SERPINE1 and IFNAR2 genes to the emergence of post-COVID-19 respiratory symptoms, providing additional consideration in treatment management to improve the outcomes of patients with COVID-19.
The ethical approval of this study was issued by the Health Research Ethics Committee, Faculty of Medicine, Universitas Indonesia (KET-872/UN2.FI/ETIK/PPM.00.02/2023) and by the Health Research Ethics Committee of Persahabatan General Hospital (98/KEPK-RSUPP/06/2023). This was a retrospective study involving subjects who were admitted to Persahabatan General Hospital (RSUP) Jakarta between 2020 and 2021. The study population consisted of 42 COVID-19 patients without respiratory syndrome as controls and 43 cases of post-COVID-19 respiratory symptoms who consented to be recalled in 2024 to provide their saliva samples. Subjects classified as cases if at least one of the following criteria was met: hospitalization duration of ≥28 days; or < 28 days with pulmonary fibrosis; or < 28 days with no improvement of radiological findings and at least one persistent respiratory symptom (cough, dyspnea, malaise, chest pain, or sore throat) at the time the patient was discharged.4 The control population consists of subjects with none of these criteria.
The cases in our study were classified as mild, moderate, severe, and critical, in which mild was characterized by symptoms without evidence of viral pneumonia or hypoxia, fever, cough, fatigue, myalgia, or anosmia/ageusia, with SpO2 >95% on room air. Moderate cases were defined by pneumonia, fever, cough, dyspnea, or tachypnea, without features of severe pneumonia, with SpO2 >93% on room air. Severe patients were indicated by pneumonia accompanied by at least one of the following: respiratory rate > 30 breaths/minute, severe respiratory distress, or SpO2 <93% on room air. Critical subjects were those with acute respiratory distress syndrome (ARDS), sepsis or septic shock, or other conditions requiring life-support interventions such as mechanical ventilation or vasopressor therapy.
DNA was obtained from 2 mL of saliva to prevent the feeling of discomfort during sample retrieval. Saliva was collected in a sterile 10 mL tube filled with preservative solution, then immediately transported to the laboratory for DNA extraction using Zymo Miniprep Kit (Zymo Research, Irvine, CA) according to the manufacturer’s instructions with modification. Quantification of the DNA samples was performed using Qubit Fluorometer (Thermo Fisher Scientific, Waltham, MA), while the purity was evaluated with Varioskan LUX multimode microplate reader (Thermo Fisher Scientific, Waltham, MA).
The genotypes of rs6092, rs1051393, and rs2229207 were evaluated with the TaqMan Assay (Applied Biosystems, Waltham, MA) using normalized DNA samples at 5 ng/μL. About 4 μL of the DNA samples were mixed with 12.5 μL of 2x TaqPath ProAmp Master Mix, 1.25 μL of TaqMan SNP Genotyping Assay, and 7.25 μL of nuclease-free water (NFW) in a 96-well plate. The thermal cycling was done in Applied Biosystem 7500 Real Time PCR (Applied Biosystems, Waltham, MA). The genotype of each variant was assigned by the autocalling feature of the TaqMan Genotyper Software v1.7.1, and then the analysis results were exported in a table file.
Data in this study were analyzed using SPSS v25.0. Statistical analysis of the demographic data, study population, and disease severity classification was assessed using Pearson’s Chi-Square test, while the analysis of blood laboratory results, arterial blood gas, and clotting blood factors was performed using the Mann-Whitney U test, Kruskal-Wallis, or the Independent t-test. We conducted linear and logistic regression to evaluate the association of rs6092, rs1051393, and rs2229207 with post-COVID-19 respiratory symptoms and blood laboratory parameters. To account for multiple comparisons, p-values were adjusted using the Benjamin–Hochberg (BH) false discovery rate (FDR), in which the global FDR correction was applied across all SNP–phenotype tests. Clustered FDR correction was conducted within predefined clinical and laboratory domains, including symptoms, electrolytes, complete blood count (CBC) parameters, liver function tests, renal function markers, and blood gas analyses. An adjusted p-value of <0.05 was considered significant.
We examined the demographic and hospitalization history of patients with post-COVID-19 respiratory syndrome and control subjects ( Table 1). The average age of patients with post-COVID-19 respiratory syndrome was higher than control, although the difference was not significant (51.23 vs. 47.33, p = 0.178). Despite no significant differences in terms of sex, year of admission, or body mass index (BMI), we noted significant differences in disease severity upon admission, in which we observe that more COVID-19 patients with respiratory syndrome were classified as severe to critical conditions, while the majority of patients with no respiratory syndrome were mostly had moderate severity (p = 0.02) ( Table 1). However, this difference was no longer significant after adjusting for BH-correction ( Table 2).
| Case | Control | p-value | |
|---|---|---|---|
| Age | 51.23 ± 13.388 | 47.33 ± 13.094 | 0.178* |
| Sex | M: 34, F: 9, n = 43 | M:29, F: 13, n = 42 | 0.292** |
| Year of Admission | 2020: 21, 2021: 22, n = 43 | 2020: 18, 2021: 24, n = 42 | 0.580** |
| BMI | 26.30 (23.59-28.53) | 26.81 (24.05-34.19) | 0.401*** |
| Severity on Admission | |||
| Mild | 2 (4.7%) | 0 (0.0%) | 0.023 |
| Moderate | 17 (39.5%) | 30 (71.4%) | |
| Severe | 17 (39.5%) | 11 (26.2%) | |
| Critical | 6 (14.0%) | 1 (2.4%) | |
| Radiology on Admission | |||
| Normal or Near Normal | 2 (4.7%) | 6 (14.3%) | 0.013 |
| Reversible Lesion | 36 (83.7%) | 36 (85.7%) | |
| Irreversible Lesion (Fibrosis) | 5 (11.6%) | 0 (0.0%) | |
| Hospitalization Time | 16.00 (11.00–20.00) | 9.50 (6.75–12.00) | 0.0001 |
| The Worst Clinical Severity Throughout Hospitalization | |||
| Mild | 0 (0.0%) | 0 (0.0%) | 0.004 |
| Moderate | 16 (37.2%) | 28 (66.7%) | |
| Severe | 19 (44.2%) | 12 (28.6%) | |
| Critical | 8 (18.6%) | 2 (4.8%) | |
| End of Radiology Series | |||
| Normal or Improvement with Near Normal result | 1 (2.3%) | 10 (31.3%) | 0.001 |
| Improvement, not Normal | 33 (25.6%) | 22 (68.8%) | |
| Stationary | 15 (34.9%) | 0 (0.0%) | |
| Deteroriation | 16 (37.2%) | 0 (0.0%) | |
| Discharge Outcome | |||
| Cured | 40 (93%) | 41 (97.6%) | 0.32*** |
| Cured with Complication | 3 (7%) | 1 (2.4%) | |
| Symptoms | |||
| Fever | 32 (74.4%) | 31 (73.85) | 0.949** |
| Cough | 35 (81.4%) | 32 (76.2%) | 0.557** |
| Malaise | 6 (14.0%) | 6 (14.3%) | 0.965** |
| Sorethroat | 2 (4.7%) | 2 (4.8%) | 1^ |
| Dyspnoea | 36 (83.7%) | 29 (69%) | 0.111** |
| Chest Pain | 1 (2.3%) | 2 (4.8%) | 0.616^ |
Subjects with post-COVID-19 respiratory syndrome had the highest disease severity during the course of hospitalization, in which we noted that these patients were further deteriorating towards the more severe and critical status (p = 0.004; padj = 0.021). This data was also supported by other clinical findings, in which we noted that patients with post-COVID-19 respiratory syndrome had significantly worse radiological examination upon admission (p = 0.013; padj = 0.052) and at the end of hospitalization (p = 0.001; padj = 0.008). Moreover, we observed that patients with post-COVID19 respiratory syndrome had longer hospitalization compared to patients with no respiratory syndrome (p = 0.0001; padj = 0.00016), although these differences may have been contributed by the inherently different classification between the case and control groups. In addition, we found that most patients in both groups had fever, cough, and malaise, but no significant differences were observed ( Table 1).
The allele and genotype frequencies of rs6092, rs1051393, and rs2229207 are summarized in Table 3. We found that 17.60% of individuals carry GG and 82% have GA genotypes for rs6092 (G > A). Meanwhile, genotype frequencies were 16.7% TT, 34.7% TG, and 48.6% GG for rs1051393 (T > G), while for rs2229207 (T > C), TT, TC, and CC genotypes were observed in 74.10%, 21.20%, and 4.70% of participants, respectively. The most frequent genotypes were GA for rs6092, GG for rs1051393, and TT for rs2229207. Allele frequencies were 58.82% G and 41.18% A for rs6092; 34% T and 66% G for rs1051393; and 84.7% T and 15.30% C for rs2229207. Hardy–Weinberg equilibrium (HWE) analysis demonstrated that all SNPs were in equilibrium, with χ2 values of 0.055 for rs6092 (df = 1, critical value = 3.41), 0.223 for rs1051393 (df = 2, critical value = 5.99), and 1.354 for rs2229207 (df = 2, critical value = 5.99), all below their respective critical values. These results suggested that the genotype distributions were consistent with HWE, supporting the genetic stability and reliability of the study population.
Our data showed no significant differences in genotype distributions of rs6092, rs1051393, and rs2229207 between patients with and without post-COVID-19 respiratory syndrome (p > 0.05). Genetic association analyses across various inheritance models indicated that none of rs6092, rs1051393, or rs2229207 were significantly associated with post-COVID-19 respiratory syndrome after correction for multiple comparisons ( Table 5). Consistently, genotype distributions for rs6092, rs1051393, and rs2229207 did not differ between the two groups, with cut-off p-values <0.05 and 95% CI, indicating no association of genotypes and post-COVID-19 respiratory syndrome in our study ( Table 6).
| SNP | Minor/Major Allele | Genetic Model | Genotype Frequency (Cases) n (%) | Genotype Frequency (Controls) n (%) | OR (95% CI) | p-value |
|---|---|---|---|---|---|---|
| rs6092 | A/G | Dominant | - | - | 0.91 (0.28–2.96) | 0.869 |
| - G/G | 8 (18.2%) | 7 (16.7%) | 1.00 (Reference) | - | ||
| - A/G + A/A | 36 (81.8%) | 35 (83.3%) | 0.91 (0.28–2.96) | 0.869 | ||
| Codominant | - | - | - | 0.874‡ | ||
| - G/G | 8 (18.2%) | 7 (16.7%) | 1.00 (Reference) | - | ||
| - A/G | 36 (81.8%) | 35 (83.3%) | 0.91 (0.27–3.01) | 0.874 | ||
| - A/A | 0 (0%) | 0 (0%) | - | - | ||
| Recessive | - | - | - | - | ||
| - G/G + A/G | 44 (100%) | 42 (100%) | 1.00 (Reference) | - | ||
| - A/A | 0 (0%) | 0 (0%) | - | - | ||
| Overdominant | - | - | 1.12 (0.34–3.68) | 0.852 | ||
| - A/G | 36 (81.8%) | 35 (83.3%) | 1.12 (0.34–3.68) | 0.852 | ||
| - G/G + A/A | 8 (18.2%) | 7 (16.7%) | 1.00 (Reference) | - | ||
| rs1051393 | G/T | Dominant | - | - | 1.05 (0.29–3.85) | 0.941 |
| - T/T | 5 (14.7%) | 7 (17.9%) | 1.00 (Reference) | - | ||
| - T/G + G/G | 29 (85.3%) | 32 (82.1%) | 1.05 (0.29–3.85) | 0.941 | ||
| Codominant | - | - | - | 0.965‡ | ||
| - T/T | 5 (14.7%) | 7 (17.9%) | 1.00 (Reference) | - | ||
| - T/G | 11 (32.4%) | 14 (35.9%) | 0.98 (0.24–4.03) | 0.973 | ||
| - G/G | 18 (52.9%) | 18 (46.2%) | 1.12 (0.28–4.52) | 0.872 | ||
| Recessive | - | - | 0.76 (0.30–1.94) | 0.569 | ||
| - G/G | 18 (52.9%) | 18 (46.2%) | 0.76 (0.30–1.94) | 0.569 | ||
| - T/T + T/G | 16 (47.1%) | 21 (53.8%) | 1.00 (Reference) | - | ||
| Overdominant | - | - | 1.12 (0.42–3.02) | 0.818 | ||
| - T/G | 11 (32.4%) | 14 (35.9%) | 1.12 (0.42–3.02) | 0.818 | ||
| - T/T + G/G | 23 (67.6%) | 25 (64.1%) | 1.00 (Reference) | - | ||
| rs2229207 | C/T | Dominant | - | - | 1.82 (0.60–5.47) | 0.289 |
| - T/T | 31 (70.5%) | 33 (78.6%) | 1.00 (Reference) | - | ||
| - T/C + C/C | 13 (29.5%) | 9 (21.4%) | 1.82 (0.60–5.47) | 0.289 | ||
| Codominant | - | - | - | 0.584‡ | ||
| - T/T | 31 (70.5%) | 33 (78.6%) | 1.00 (Reference) | - | ||
| - T/C | 10 (22.7%) | 8 (19.0%) | 1.85 (0.58–5.89) | 0.300 | ||
| - C/C | 3 (6.8%) | 1 (2.4%) | 1.21 (0.07–22.02) | 0.898 | ||
| Recessive | - | - | 1.00 (0.06–17.33) | >0.999 | ||
| - C/C | 3 (6.8%) | 1 (2.4%) | 1.00 (0.06–17.33) | >0.999 | ||
| - T/T + T/C | 41 (93.2%) | 41 (97.6%) | 1.00 (Reference) | - | ||
| Overdominant | - | - | 0.53 (0.18–1.61) | 0.263 | ||
| - T/C | 10 (22.7%) | 8 (19.0%) | 0.53 (0.18–1.61) | 0.263 | ||
| - T/T + C/C | 34 (77.3%) | 34 (81.0%) | 1.00 (Reference) | - |
In this study, we performed statistical analyses to identify associations of rs6092, rs1051393, and rs2229207 with post-COVID-19 respiratory syndromes ( Table 7). Our analysis indicates an association of IFNAR2 rs1051393 with fever (p = 0.037), in which further logistic regression analysis showed that rs1051393 was significantly associated with fever in two genetic models. While the recessive model indicated that an individual who carries the recessive genotype has higher odds of fever (OR = 4.42, 95% CI: 1.30–15.07, p = 0.018), the over-dominant model showed lower odds of fever (OR = 0.144, 95% CI: 0.029–0.726, p = 0.019). In this case, the codominant model did not show a statistically relevant association for either genotype comparison (heterozygote vs wild-type: OR = 4.17, 95% CI: 0.568–30.597, p = 0.160; homozygote vs wild-type: OR = 0.502, 95% CI: 0.109–2.310, p = 0.376). Our analysis also indicates an association of rs2229207 IFNAR2 in patients with chest pain (p = 0.047), although no relevant association of this rs2229207 with chest pain was revealed in our subsequent logistic regression analysis.
| Symptom | rs6092 SERPINE1 | rs1051393 IFNAR2 | rs2229207 IFNAR2 | ||||||
|---|---|---|---|---|---|---|---|---|---|
| GG | GA | TT | TG | GG | TT | CT | CC | ||
| Fever | Yes | 10 | 53 | 9 | 23 | 22 | 46 | 14 | 3 |
| No | 5 | 17 | 3 | 2 | 13 | 17 | 4 | 1 | |
| p-value | 0.468* | 0.07 **/0.037* | 0.708** | ||||||
| Cough | Yes | 12 | 55 | 10 | 19 | 28 | 51 | 13 | 3 |
| No | 3 | 15 | 2 | 6 | 7 | 12 | 5 | 1 | |
| p-value | 1.000** | 0.864* | 0.430** | ||||||
| Malaise | Yes | 11 | 62 | 10 | 23 | 28 | 53 | 17 | 3 |
| No | 4 | 8 | 2 | 2 | 7 | 10 | 1 | 1 | |
| p-value | 0.212** | 0.427** | 0.509** | ||||||
| Sorethroat | Yes | 1 | 3 | 0 | 1 | 3 | 3 | 1 | 0 |
| No | 14 | 67 | 12 | 24 | 32 | 60 | 17 | 4 | |
| p-value | 0.547** | 0.238** | 0.924** | ||||||
| Dyspnoe | Yes | 11 | 54 | 10 | 18 | 26 | 48 | 13 | 4 |
| No | 4 | 16 | 2 | 7 | 9 | 15 | 5 | 0 | |
| p-value | 0.745** | 0.751* | 0.812** | ||||||
| Chest Pain | Yes | 1 | 2 | 0 | 1 | 2 | 2 | 0 | 1 |
| No | 14 | 68 | 12 | 24 | 33 | 61 | 18 | 3 | |
| p-value | 0.446** | 0.433** | 0.564 **/0.047* | ||||||
We identified significant differences in basophil levels (p = 0.027) and eGFR (p = 0.036) for rs6092 of the SERPINE1, in which rs6092 was associated with higher levels of basophil and lower eGFR in blood. Meanwhile, we also note that rs2229207 of the IFNAR2 was significantly associated with potassium (p = 0.001) and chloride levels (p = 0.047), in which we observed that individuals with the CC genotype had higher potassium and chloride in blood workup ( Table 8). Blood gas analyses also revealed significant differences in CO2 partial pressure (p = 0.030) in rs2229207, with individuals carrying TC and CC genotypes having higher pCO2 ( Table 9). However, many of these associations did not withstand the global Benjamini–Hochberg FDR correction. Only rs2229207, which was associated with potassium levels, persisted after clustered Benjamini–Hochberg FDR adjustment (padj = 0.012) ( Table 11). Furthermore, our analysis on coagulation and inflammatory factors revealed no relevant association with rs6092, rs1051393, and rs2229207 in all subjects ( Table 10).
| Parameters | rs6092 SERPINE1 | p-value | rs1051393 IFNAR2 | p-value | rs2229207 IFNAR2 | p-value | |||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| GG (n = 15) | GA (n = 70) | TT (n = 12) | TC (n = 25) | GG (n = 35) | TT (n = 63) | TC (n = 18) | CC (n = 4) | ||||
| Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | ||||
| Sodium (Na), mmol/L | 135.9 ± 3.3 | 133.9 ± 4.01 | 0.071* | 133.7 ± 4.1 | 134.7 ± 3.7 | 133.8 ± 4.6 | 0.689** | 134.3 ± 4.3 | 133.6 ± 2.8 | 136.8 ± 2.5 | 0.344** |
| Potassium (K), mmol/L | 3.9 (3.5–4.5) | 3.9 (3.5–4.3) | 0.560*** | 3.9 (3.5–4.2) | 3.8 (3.6–4.2) | 3.9 (3.7–4.4) | 0.476† | 3.8 (3.5–4.1) | 4.2 (3.9–4.6) | 4.6 (4.5–4.8) | 0.001 † |
| Chloride (Cl), mmol/L | 102.1 ± 4.6 | 99.8 ± 4.2 | 0.063* | 99.7 ± 3.8 | 99.4 ± 4.4 | 100.8 ± 5.1 | 0.457** | 100.4 ± 4.2 | 98.9 ± 4.4 | 104.8 ± 4.2 | 0.047 ** |
| Hemoglobin (Hb), g/dL | 13.9 (11.6–15.0) | 13.7 (12.5–14.4) | 0.493*** | 14.3 (13.5–15.3) | 12.9 (11.6–14.9) | 13.8 (12.8–14.4) | 0.257† | 13.7 (12.4–14.5) | 13.9 (12.7–14.3) | 14.4 (11.3–15.2) | 0.766† |
| Leukocytes (Leu), 10 3/μL | 8.80 (6.28–9.78) | 7.74 (5.97–10.1) | 0.397*** | 7.12 (4.97–9.43) | 7.74 (6.45–8.96) | 8.23 (5.73–9.99) | 0.369† | 7.78 (6.10–9.20) | 7.96 (5.67–10.2) | 9.79 (7.39–13.9) | 0.518† |
| Aspartate Transaminase (AST), U/L | 273 (226–286) | 230 (185–309) | 0.440*** | 205 (183–268) | 263 (185–351) | 253 (207–295) | 0.718† | 234 (187–298) | 267 (179–304) | 243 (235–291) | 0.849† |
| Neutrophils (Neu), % | 81.0 (67–84.1) | 78.6 (67.0–85.8) | 0.708*** | 78.6 (67.1–81.5) | 79.1 (74.2–83.9) | 77.9 (66.7–86.4) | 0.704† | 78.3 (67.0–84.2) | 81.1 (70.3–90.4) | 77.4 (60.3–86.0) | 0.521† |
| Eosinophils (Eos), % | 0.0 (0.0–0.3) | 0.0 (0.0–0.3) | 0.747*** | 0.0 (0.0–0.2) | 0.0 (0.0–0.2) | 0.0 (0.0–0.3) | 0.973† | 0.0 (0.0–0.4) | 0.0 (0.0–0.1) | 0.2 (0.1–2.6) | 0.115† |
| Basophils (Bas), % | 0.1 (0.0–0.1) | 0.1 (0.0–0.2) | 0.027 *** | 0.1 (0.0–0.2) | 0.1 (0.0–0.2) | 0.1 (0.1–0.2) | 0.807† | 0.1 (0.0–0.2) | 0.1 (0.1–0.2) | 0.2 (0.1–0.4) | 0.884† |
| Lymphocytes (Lym), % | 12.8 (8.6–17.4) | 13.4 (8.6–22.9) | 0.475*** | 13.2 (9.7–17.6) | 13.4 (10.8–17.9) | 11.4 (6.7–23.0) | 0.663† | 14.0 (9.1–22.9) | 13.0 (6.0–22.4) | 6.7 (3.8–12.1) | 0.156† |
| Monocytes (Mo), % | 6.7 ± 3.4 | 7.3 ± 5.5 | 0.546* | 9.2 ± 4.6 | 6.8 ± 2.9 | 7.0 ± 3.7 | 0.165** | 7.4 ± 3.6 | 6.4 ± 3.7 | 8.1 ± 3.7 | 0.483** |
| NLR | 6.18 (3.86–9.0) | 5.76 (3.01–10.15) | 0.936*** | 5.08 (4.04–8.44) | 5.90 (4.34–7.49) | 6.05 (2.85–10.8) | 0.837† | 5.42 (3.06–9.03) | 6.19 (3.14–15.1) | 6.83 (3.10–17.2) | 0.685† |
| OT, U/L | 50.0 (31.0–70.0) | 44.0 (27.0–71.0) | 0.583*** | 46.0 (32.0–70.0) | 40.0 (28.0–74.0) | 45.0 (27.0–69.0) | 0.981† | 44.0 (26.0–70.0) | 48.0 (35.0–74.0) | 40.0 (25.0–56.0) | 0.606† |
| PT, U/L | 53.0 (39.0–67.0) | 43.0 (28.0–88.0) | 0.870*** | 41.0 (33.0–71.0) | 40.0 (28.0–100) | 47.0 (33.0–84.0) | 0.976† | 41.0 (28.0–88.0) | 58.0 (38.0–87.0) | 62.0 (34.0–90.0) | 0.599† |
| Blood Urea Nitrogen (BUN), mg/dL | 15.4 (10.28–20.07) | 13.6 (10.28–17.52) | 0.556*** | 13.1 (9.1–15.2) | 13.1 (10.2–18.3) | 13.6 (11.2–17.3) | 0.476† | 13.8 (10.3–17.8) | 11.9 (10.3–16.8) | 23.6 (11.2–31.1) | 0.575† |
| Creatinine (Cr), mg/dL | 0.9 (0.7–0.9) | 0.9 (0.8–1.1) | 0.498*** | 0.9 (0.8–1.1) | 0.8 (0.7–1.1) | 0.9 (0.8–1.1) | 0.801† | 0.9 (0.8–1.0) | 0.9 (0.8–1.1) | 1.0 (0.9–1.2) | 0.654† |
| eGFR, mL/min/1.73 m 2 | 95.0 (92.0–108) | 87.0 (72.0–105) | 0.036 * | 92.0 (75.0–105) | 92.0 (76.0–108) | 91.0 (78.0–100) | 0.849† | 92.0 (77.0–107) | 91.0 (76.0–103) | 91.0 (68.0–93.0) | 0.666† |
| Parameters | rs6092 SERPINE1 | p-value | rs1051393 IFNAR2 | p-value | rs2229207 IFNAR2 | p-value | |||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| GG (n = 15) | GA (n = 70) | TT (n = 12) | TC (n = 25) | GG (n = 35) | TT (n = 63) | TC (n = 18) | CC (n = 4) | ||||
| Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | Mean ± SD/Median (Q1-Q3) | ||||
| Bicarbonate (HCO3), mmol/L | 21.8 ± 2.4 | 21.6 ± 3.4 | 0.806 | 22.5 ± 2.6 | 22.0 ± 2.9 | 21.1 ± 3.7 | 0.426** | 21.4 ± 3.1 | 22.5 ± 3.8 | 22.5 ± 3.8 | 0.378** |
| O2 Saturation (O2_Sat), % | 96.8 (91.2–99.4) | 98.3 (95.9–99.4) | 0.286*** | 98.0 (96.7–98.6) | 98.5 (96.1–99.5) | 97.7 (95.1–99.4) | 0.511† | 98.2 (96.0–99.4) | 97.8 (92.8–99.4) | 97.7 (79.0–99.7) | 0.726† |
| pCO2, mmHg | 34.4 ± 4.3 | 33.4 ± 6.2 | 0.549* | 34.9 ± 5.2 | 33.6 ± 5.1 | 32.9 ± 7.1 | 0.697** | 32.7 ± 5.1 | 36.7 ± 7.7 | 35.3 ± 6.3 | 0.030 ** |
| pH | 7.4 ± 0.0 | 7.4 ± 0.0 | 0.463* | 7.414 ± 0.040 | 7.421 ± 0.048 | 7.414 ± 0.046 | 0.796** | 7.420 ± 0.041 | 7.394 ± 0.049 | 7.408 ± 0.025 | 0.084** |
| pO2, mmHg | 92.9 (66.4–131.5) | 102.5 (79.1–144.5) | 0.500*** | 93.1 (81.0–114) | 102.5 (79.8–183) | 96.3 (73.6–126) | 0.424† | 100 (78.4–146) | 99.9 (77.8–133) | 118 (63.1–172) | 0.911† |
| Standard HCO3, mmol/L | 22.9 ± 2.15 | 23.0 ± 2.5 | 0.887* | 23.7 ± 2.0 | 23.4 ± 2.4 | 22.7 ± 2.7 | 0.478** | 23.0 ± 2.4 | 23.4 ± 2.9 | 23.3 ± 2.5 | 0.785** |
| Total CO2, mmol/L | 22.8 ± 2.5 | 22.9 ± 4.0 | 0.989* | 24.4 ± 4.1 | 23.0 ± 3.0 | 22.1 ± 3.9 | 0.262** | 22.7 ± 3.7 | 23.6 ± 4.0 | 23.6 ± 4.0 | 0.615** |
| Parameters | rs6092 SERPINE1 | p-value | rs1051393 IFNAR2 | p-value | rs2229207 IFNAR2 | p-value | |||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| GG (n = 15) | GA (n = 70) | TT (n = 12) | TC (n = 25) | GG (n = 35) | TT (n = 63) | TC (n = 18) | CC (n = 4) | ||||
| Median (Q1-Q3) | Median (Q1-Q3) | Median (Q1-Q3) | Median (Q1-Q3) | Median (Q1-Q3) | Median (Q1-Q3) | Median (Q1-Q3) | Median (Q1-Q3) | ||||
| D-dimer, ng/mL | 820 (420–2060) | 610 (360–1070) | 0.440* | 560 (315–1455) | 880 (430–1280) | 600 (420–880) | 0.317† | 580 (360–1280) | 650 (420–1040) | 1955 (1345–18630) | 0.059† |
| C-Reactive Protein (CRP), mg/L | 78.7 (60.1–94.8) | 75.5 (21.6–114.2) | 0.881* | 88.2 (28.0–131) | 81.5 (41.3–109) | 75.8 (13.3–108) | 0.733† | 73.4 (21.6–107) | 96.6 (54.1–159) | 73.9 (34.2–113) | 0.284† |
All models adjusted for year of admission, age on admission, sex, and severity on admission. TT is the reference for IFNAR2, and GG is the reference for SERPINE1.
Our subsequent multivariable linear regression analysis (Table 14) showed that rs2229207 genotype was significantly associated with potassium levels (B = 0.317, 95% CI: 0.122–0.596, p = 0.003) after adjusted for age, sex, BMI, and severity, in which we noted an interesting finding that female individuals was independently associated with lower potassium levels (B = −0.219, 95% CI: −0.611 to −0.015, p = 0.040). To clarify the mode of inheritance, we performed further analysis of rs6092 and rs2229207 genotypes under different inheritance models ( Table 15). Our analysis revealed that rs6092 was significantly associated with basophil, eGFR, potassium, chloride, pCO2, pH, and D-dimer in blood across dominant, codominant, and over-dominant models (all p < 0.05). While for rs2229207, potassium showed consistent significant associations across all genetic models (p < 0.05), with additional significant associations observed for pCO2, pH, chloride, and D-dimer in specific models ( Table 15).
| SNPs | Genetic Model | n | p-value | |||||||
|---|---|---|---|---|---|---|---|---|---|---|
| Basophil | eGFR | K | Cl | pCO2 | pH | D-dimer | ||||
| rs6092 | Dominant | GG | 15 | 0.027 * | 0.036 *** | 0.56*** | 0.63* | 0.549* | 0.463* | 0.44*** |
| GA/AA | 70 | |||||||||
| Codominant | GG | 15 | 0.027 * | 0.036 *** | 0.56*** | 0.63 | 0.549* | 0.463* | 0.44*** | |
| GA | 70 | |||||||||
| AA | - | |||||||||
| Recessive | GG/GA | 85 | - | - | - | - | - | - | - | |
| AA | 0 | |||||||||
| Overdominant | GA | 70 | 0.027 * | 0.036 *** | 0.56*** | 0.63 | 0.549* | 0.463* | 0.44*** | |
| GG/AA | 15 | |||||||||
| rs2229207 | Dominant | TT | 63 | 0.851*** | 0.443*** | 0.001 *** | 0.722* | 0.032 * | 0.031 * | 0.211*** |
| TC/CC | 22 | |||||||||
| Codominant | TT | 63 | 0.884† | 0.666† | 0.001 † | 0.047 ** | 0.793** | 0.084** | 0.059† | |
| TC | 18 | |||||||||
| CC | 4 | |||||||||
| Recessive | CC | 4 | 0.62*** | 0.467*** | 0.009 *** | 0.033 * | 0.561* | 0.785* | 0.019 *** | |
| TT/TC | 81 | |||||||||
| Overdominant | TC | 18 | 0.955*** | 0.655*** | 0.02 *** | 0.139* | 0.054* | 0.03 * | 0.902*** | |
| TT/CC | 67 | |||||||||
Multivariable linear regression analysis was performed to evaluate the associations of rs6092 and rs2229207 with clinical parameters after adjusting for age, sex, severity, and year of admission ( Table 16). Our data revealed that rs6092 was not significantly relevant to eGFR (B = −0.120; 95% CI: −25.665 to 7.118; p = 0.264) and basophil percentage (B = 0.111; 95% CI: −0.201 to 0.588; p = 0.332). In contrast, rs2229207 was significantly associated with increased pCO2 levels under the dominant model (B = 2.417; 95% CI: 0.614 to 6.347; p = 0.018) and showed significant effects on pH under both dominant (B = −0.228; 95% CI: −0.043 to −0.001; p = 0.040) and over-dominant models (B = 0.219; 95% CI: 0.00008 to 0.046; p = 0.049). Furthermore, rs2229207 was strongly associated with potassium levels across dominant (B = 0.280; 95% CI: 0.099 to 0.701; p = 0.010), codominant (B = 0.317; 95% CI: 0.122 to 0.596; p = 0.003), and recessive models (B = −0.255; 95% CI: −1.398 to −0.109; p = 0.022), as well as with chloride (recessive model: B = −0.249; 95% CI: −9.456 to −0.772; p = 0.022) and D-dimer levels (recessive model: B = −0.333; 95% CI: −12976.164 to −2869.870; p = 0.003). Among covariates, age was negatively associated with eGFR (B = −0.224; 95% CI: −0.972 to −0.024; p = 0.040), while sex was associated with potassium levels (B = −0.219; 95% CI: −0.611 to −0.015; p = 0.040), and disease severity showed a negative association with chloride (B = −0.265; 95% CI: −2.948 to −0.284; p = 0.018). Multicollinearity diagnostics indicated no significant multicollinearity with all VIF < 1.2.
| Predictor | eGFR B (95% CI), *p* | Basophil % B (95% CI), *p* | Chloride B (95% CI), *p* | pCO2 B (95% CI), *p* | pH B (95% CI), *p* | D-dimer B (95% CI), *p* | Potassium B (95% CI), *p* |
|---|---|---|---|---|---|---|---|
| rs6092 | |||||||
| Dominant | −0.120 (−25.665, 7.118), 0.264 | 0.111 (−0.201, 0.588), 0.332 | - | - | - | - | - |
| Codominant | −0.120 (−25.665, 7.118), 0.264 | 0.111 (−0.201, 0.588), 0.332 | - | - | - | - | - |
| Overdominant | −0.120 (−25.665, 7.118), 0.264 | 0.111 (−0.201, 0.588), 0.332 | - | - | - | - | - |
| rs2229207 | |||||||
| Dominant | - | - | - | 2.417 (0.614, 6.347), 0.018 | −0.228 (−0.043, −0.001), 0.040 | - | 0.280 (0.099, 0.701), 0.010 |
| Codominant | - | - | 0.107 (−0.836, 2.513), 0.322 | - | - | - | 0.317 (0.122, 0.596), 0.003 |
| Recessive | - | - | −0.249 (−9.456, −0.772), 0.022 | - | - | −0.333 (−12976.164, −2869.870), 0.003 | −0.255 (−1.398, −0.109), 0.022 |
| Overdominant | - | - | - | - | 0.219 (0.00008, 0.046), 0.049 | - | −0.175 (−0.601, 0.064), 0.112 |
| Covariates | |||||||
| Age on admission | −0.224 (−0.972, −0.024), 0.040 1 | −0.89 (−0.016, 0.007), 0.4382 | −0.015 (−0.074, 0.064), 0.8873 | −0.011 (−0.099, 0.090), 0.9214 | 0.155 (−0.0002, 0.001), 0.1554 | 0.150 (−21.824, 136.281), 0.1545 | 0.040 (−0.008, 0.012), 0.7006 |
| Sex | 0.186 (−1.760, 26.771), 0.0851 | −0.028 (−0.386, 0.301), 0.8052 | 0.119 (−0.927, 3.291), 0.2683 | −0.181 (−5.306, 0.430), 0.0954 | −0.036 (−0.025, 0.018), 0.7434 | −0.077 (−3297.891, 1537.075), 0.4715 | −0.219 (−0.611, −0.015), 0.040 6 |
| Severity on admission | 0.152 (−2.578, 15.079), 0.1631 | 0.048 (−0.167, 0.258), 0.6732 | −0.265 (−2.948, −0.284), 0.018 3 | 0.068 (−1.238, 2.363), 0.5364 | −0.004 (−0.014, 0.013), 0.9714 | 0.124 (−632.107, 2378.494), 0.2525 | 0.057 (−0.138, 0.238), 0.6006 |
| Year of admission | 0.040 (−9.410, 15.361), 0.6341 | −0.032 (−0.341, 0.256), 0.7782 | −0.203 (−3.596, 0.059), 0.0583 | 1.110 (−1.384, 3.604), 0.3784 | −0.114 (−0.028, 0.009), 0.2954 | 0.001 (−2073.143, 2099.927), 0.9905 | 0.040 (−0.208, 0.309), 0.6996 |
COVID-19 patients with higher disease severity have an increased risk of developing post-COVID-19 respiratory syndrome.7 In unadjusted analyses, our study identified a significant association between higher disease severity upon admission (p = 0.023) and during the course of hospitalization (p = 0.004) within subjects with post-COVID-19 respiratory syndromes. However, we noted that only the disease severity during the course of hospitalization remained significant after adjusting for multiple comparisons using the Benjamini-Hochberg method (padj = 0.021). The severe symptoms in COVID-19 patients are often associated with a cytokine storm. This storm occurs when SARS-CoV-2 infects the lung tissue via the ACE2 receptor and induces an immune response that releases pro-inflammatory cytokines, such as IL-6, GM-CSF, and IFN-γ.14 This abundant release led to the infiltration of macrophages and neutrophils into the lung tissue.15,16 The persistence of the virus in tissues and a dysfunctional immune response during the acute phase can impair the enzymatic function of the ACE2 receptor and worsen respiratory symptoms. Histopathological studies have demonstrated that SARS-CoV-2 can persist in tissues for up to 230 days, potentially triggering repeated immune responses, which lead to the emergence of prolonged symptoms.17,18
The post-COVID-19 respiratory syndrome refers to respiratory symptoms that persist for at least four weeks following COVID-19 infection.4 In our study cohort, we found a higher prevalence of poor pulmonary radiological examination upon admission (p = 0.013) and at the end of hospitalization (p = 0.001) in cases. Additionally, we noted that subjects with post-COVID-19 respiratory syndrome had significantly longer hospitalization compared to control subjects. This finding might be attributed to the poor radiological examination, as we observed a significantly higher prevalence of fibrosis in subjects with post-COVID-19 respiratory syndrome. This finding is consistent with previous reports that individuals infected with COVID-19 who developed acute respiratory distress syndrome (ARDS) and require intubation during the acute phase have a threefold increased risk of developing irreversible pulmonary fibrosis as the main irreversible lesion.19
The irreversible lesion in the lungs of COVID-19 patients is potentially the result of hyperinflammation due to the downregulation of ACE2 and upregulation of angiotensin II, which had a negative impact on cells. The hyperinflammation in patients with COVID-19 is reported due to high cytokine release and immunothrombosis which activates the blood-clotting factors such as PAI-1.17,20 Previous studies also described that the infiltration of macrophages and neutrophils led to the expression of transforming growth factor-beta (TGF-β) and fibroblast activation. These fibroblasts differentiate into myofibroblasts, which produce extracellular matrix (ECM) components, such as collagen and fibronectin.21,22 The increased formation of ECM impaired the repair of the alveolar epithelium, increased extracellular collagen matrix formation, and led to the loss of normal lung architecture.19 Additionally, mechanical ventilation may further aggravate this condition, leading to ventilator-induced lung injury (VILI), a known contributor to pulmonary fibrosis.23
As we described previously, the incidence of pulmonary fibrosis might be attributed to the dysfunctional expression of blood-clotting factors in the lung.24,25 The human SERPINE1 gene encodes the plasminogen activator inhibitor 1 (PAI-1), which inhibits tissue plasminogen activator (tPA) by converting plasminogen to plasmin.26 Previous studies have indicated an association between the rs6092 of the SERPINE1 gene and coagulation abnormalities in COVID-19 patients.11 Moreover, another study reported that the high level of D-dimer and the inhibition of plasmin prevent fibrinolysis, leading to the excessive clotting in COVID-19 patients.26 Despite no association of rs6092 with the clinical symptoms of post-COVID-19 respiratory syndrome, we observed a significantly higher level of basophils and lower eGFR.
The increase in blood basophils might be the response to combat the presence of SARS-CoV-2, as basophils are responsible for the innate and adaptive immunity by releasing IL2RA and IL2RG.27 Despite more relevant evidence being required to explain the correlation of higher basophils in the blood of subjects with post-COVID-19 respiratory syndrome and pulmonary fibrosis, basophils-rich blood has been reported to increase the likelihood of lower extremity deep venous thrombosis (LEDVT) by 2.7 times in patients with spontaneous intracerebral hemorrhage (sICH), potentially by increasing the level of the enhanced factor II plasma coagulant.28,29 rs6092 remains a potential candidate marker with severe clinical manifestations in COVID-19 patients, as this variant caused an amino acid change of Ala15Thr located within the central hydrophobic core of the PAI-1 signal peptide. The Ala15Thr substitution increases hydrophobicity and enhances the propensity for α-helix formation, both of which are critical for proper signal peptide function.30
The ability of the immune system to combat viral infection is also contributed to by the expression of interferons. The interferon alpha and beta receptor 2, which is crucial in the antiviral response, is encoded by the IFNAR2 gene. When SARS-CoV-2 infects the cells, lymphocytes release interferons that bind to IFNAR2 receptors. This binding activates the Janus Kinase (JAK) pathway, resulting in the phosphorylation of STAT1 and STAT2. These transcription factors, in combination with IRF9, form the ISGF3 complex, which enters the nucleus and activates genes involved in the antiviral response, including MX1, OAS, and IRF7.31,32 Previous studies have demonstrated a relationship between variations in IFNAR2 and severe COVID-19 symptoms.13 In our study, both of the examined variants in the IFNAR2 gene were associated with clinical parameters, in which rs1051393 was significantly associated with fever (p = 0.037), while the rs2229207 was associated with chest pain (p = 0.047) and higher potassium (p = 0.001), chloride (p = 0.047), and pCO2 (p = 0.030).
To the best of our knowledge, this is the first study to report the associations of rs1051393 with fever and rs2229207 with chest pain. Notably, our findings also suggest an increasing trend across TT, CT, and CC genotypes of the rs2229207 and blood potassium (B = 0.317, 95% CI: 0.122 to 0.596, p = 0.003). The abnormalities in potassium levels were also described in different cohorts of studies. Ashammari et al. reported that 34 out of 72 COVID-19 patients with heart failure had higher potassium levels, which was attributed to acute kidney injury and the use of angiotensin receptor-neprilysin inhibitors.33 Interestingly, another report by Amin et al. found that hyperkalemia was observed in only 12% of COVID-19 patients, but was significantly more common among older subjects and in those with chronic kidney disease or diabetes mellitus. This study also showed that patients with hyperkalemia were more likely to be admitted to the ICU and had a longer duration of hospitalization.34
Despite limited evidence of biological relevance to the severity of COVID-19, the association of IFNAR2 variants with COVID-19 cases has been described elsewhere. A study by Nhung et al. in a Vietnamese population also found an association between rs2229207 and the risk of contracting COVID-19.35 GWAS data have also demonstrated that the IFNAR2 locus is also known to play a role in determining severe clinical manifestations of COVID-19, with one such variant being rs2236757, which is located in a non-coding region.13,36,37 In addition to its association with COVID-19, the rs2229207 has also been linked to chronic hepatitis patients, in which an amino acid substitution leads to a reduced expression of the IFNAR2 protein on the cell membrane, potentially disrupting interferon signaling.38,39 Both rs1051393 and rs2229207 were missense variants in the IFNAR2 receptor. These structural changes may interfere with the antiviral mechanism and reduce the expression of the IFNAR2 protein on the cell membrane, potentially disrupting interferon signaling.40
In our study, we did not observe a significant association of genotype distribution with control or cases. We noted that the minor allele frequency (MAF) for rs6092 in our study was 41%. The value is higher compared to other studies, which reported to be at 10.88% from 786 individuals and 7—11% as documented in NCBI dbSNP. This discrepancy may be attributed to the relatively limited sample size of our study population, which may not be fully representative when compared to global population databases. While for the rs1051393, we noted the MAF for the G allele at 34.03%, consistent with the previous study by Fricke-Galindo et al., which reported the MAF at 28—32%.11 We also found that the MAF for the C allele of rs2229207 was 15.29%. This result aligns with a prior investigation by Malvestiti et al. (2026), which documented an MAF of 8–10% for rs2229207 in a sample of 487 individuals.13 When compared to the NCBI dbSNP, the allele frequencies for both rs1051393 and rs2229207 in our study are relatively consistent with those observed in other Asian populations, which range from 42–44% for rs1051393 and 16–18% for rs2229207. Although the frequency of the T allele for rs1051393 in our study is slightly lower, the overall allelic distribution pattern remains comparable, suggesting a shared genetic background with Asians.
Our data has demonstrated possible contributions of genetic predisposition on the SERPINE1 and IFNAR2 genes to the emergence of post-COVID-19 respiratory syndromes but are subject to several limitations. While we observed significant associations respiratory symptoms with disease classification in our cohort, only disease severity during hospitalization, radiological examination results upon admission and at the end of hospitalization, and duration of hospitalization remained significant after applying the BH-correction. Moreover, we also note that only rs2229207 withstand the BH-correction after adjusting for age, sex, BMI, and disease severity. Additionally, the modest sample size would limit the statistical power in this study and restrain our analysis to explore the contribution of genetic variants on disease outcomes.
Our data showed that subjects with post-COVID-19 respiratory syndrome had significantly higher disease severity upon admission and during hospitalization. These patients had a longer hospitalization, which might be due to poor radiological examination upon admission and at the end of hospitalization. We argue that these were predisposed by genetic factors, in which we reported that rs1051393 was associated with fever, while rs2229207 was associated with chest pain, higher blood potassium, chloride, and pCO2. We also suggested that the association of rs6092 with higher basophil levels in blood might potentially contribute to the risk of blood clotting in cases with pulmonary fibrosis. However, due to modest sample size, only disease severity during hospitalization, radiological examination results upon admission and at the end of hospitalization, duration of hospitalization, and the association of rs2229207 with blood potassium levels remained significant after applying the Benjamini-Hochberg correction.
All patients involved in this study have agreed, and the study has passed the ethical review, based on the Letter of Approval of the Research and Health Ethics Committee of the Faculty of Medicine, University of Indonesia and Dr. Cipto Mangunkusumo National Central General Hospital (FKUI-RSCM), Jakarta: KET-872/UN2. F1/ETIK/PPM.002.02/2023. This research was also supported by the Persahabatan National Respiratory Hospital, Jakarta (98/KEPK-RSUPP/06/2023).
Figshare: The Genetic Association of Single Nucleotide Polymorphisms of SERPINE1 (rs6092) and IFNAR2 (rs1051393, rs2229207) Genes Is Related to Post Covid-19 Respiratory Syndrome. https://doi.org/10.6084/m9.figshare. 30338875.41
This project contains the following underlying data:
• Long Covid-19_F1000_dataset.xlsx – This file contains research data including: demographic information, clinical data, laboratory data, and genotyping results.
The data are available under the terms of the Creative Commons Attribution 4.0 International (CC BY 4.0) license.
This research is supported by the funding from the National Research and Innovation Agency (BRIN) through the Riset Inovasi dan Indonesia Maju (RIIM) grant (Contract no. 36/IV/KS/06/2022) Funding Year 2025.
| Views | Downloads | |
|---|---|---|
| F1000Research | - | - |
|
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Competing Interests: No competing interests were disclosed.
Reviewer Expertise: Molecular Biology and Genetics
Is the work clearly and accurately presented and does it cite the current literature?
Partly
Is the study design appropriate and is the work technically sound?
No
Are sufficient details of methods and analysis provided to allow replication by others?
Partly
If applicable, is the statistical analysis and its interpretation appropriate?
No
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
No
Competing Interests: No competing interests were disclosed.
Reviewer Expertise: Genetics and Genomics, Bioinformatics, Immunology, Population Genetics
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Partly
Competing Interests: No competing interests were disclosed.
Reviewer Expertise: Molecular Biology and Genetics
Alongside their report, reviewers assign a status to the article:
| Invited Reviewers | ||
|---|---|---|
| 1 | 2 | |
|
Version 2 (revision) 09 May 26 |
read | |
|
Version 1 21 Nov 25 |
read | read |
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)