Preprint
Article

Innate Immune Gene Polymorphisms and COVID‐19 Prognosis

Altmetrics

Downloads

173

Views

46

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

07 August 2023

Posted:

08 August 2023

You are already at the latest version

Alerts
Abstract
COVID-19 is characterized by a heterogeneous clinical presentation and prognosis. Risk factors contributing to the development of severe disease include old age and the presence of comorbidities. However, the genetic background of the host has also been recognized as an important determinant of disease prognosis. Considering the pivotal role of innate immunity in the control of SARS-CoV-2 infection, we analyzed the possible contribution of several innate immune gene polymorphisms (including TLR2-rs5743708, TLR4-rs4986790, TLR4-rs4986791, CD14-rs2569190, CARD8-rs1834481, IL18-rs2043211 and CD40-rs1883832) in disease severity and prognosis. A total of 249 individuals were enrolled and further divided into five (5) groups, according to the clinical progression scale provided by the World Health Organization (WHO) (asymptomatic, mild, moderate, severe and critical). We identified that elderly patients with obesity and/or diabetes mellitus were more susceptible to developing pneumonia and respiratory distress syndrome after SARS-CoV-2 infection, while the IL18-rs1834481 polymorphism was an independent risk factor for developing pneumonia. Moreover, individuals carrying either the TLR2-rs5743708 or the TLR4-rs4986791 polymorphisms exhibited a 3.6- and 2.5-fold increased probability for developing pneumonia and a more severe disease, respectively. Our data supports the notion that the host’s genetic background can significantly affect COVID-19 clinical phenotype, also suggesting that the IL18-rs1834481, TLR2-rs5743708 and TLR4-rs4986791 polymorphisms may be used as molecular predictors of COVID-19 clinical phenotype.
Keywords: 
Subject: Biology and Life Sciences  -   Immunology and Microbiology

1. Introduction

COVID-19 disease, which is caused by the emerging coronavirus SARS-CoV-2, is characterized by a heterogeneous and wide spectrum of clinical symptoms and prognoses. While most of infected individuals experience asymptomatic or mild disease, a considerable number of patients with SARS-CoV-2 infection can progress to a more severe and/or critical disease, accompanied by potentially life-threatening manifestations and complications [1]. Previous studies have already demonstrated that risk factors contributing to the development of severe COVID-19 include old age, low baseline ratio of arterial partial pressure of oxygen to fraction of inspired oxygen (pO2/FiO2), nationality, and the presence of comorbidities such as obesity, diabetes mellitus and hypertension [1,2,3,4,5,6].
Moreover, the host’s genetic background represents one of the most crucial parameters that modulate COVID-19 susceptibility and severity [7,8,9,10]. Several studies have already pointed out that variations in selected genes of the host play a pivotal role in disease progression, affecting the different COVID-19 clinical phenotypes and outcomes [7,8,9]. In this context, we and others have also analyzed functional innate immune gene polymorphisms (such as MBL2, TLR7 etc.), considering that the innate immune system plays a critical role in COVID-19 prognosis and outcome [3,11,12]. It has been established with certainty that innate immunity can determine early control of the SARS-CoV-2 infection, and its dysregulation may result in hyperinflammation events, which could lead to more severe disease and worse prognosis [1,13,14]. Thus, we could speculate that additional functional polymorphisms of the innate immune system may represent important risk factors for COVID-19 phenotype and outcome.
In this study, we explored the contribution of additional known functional innate immune gene polymorphisms, namely TLR2-rs5743708, TLR4-rs4986790, TLR4-rs4986791, CD14-rs2569190, CARD8-rs1834481, IL18-rs2043211 and CD40-rs1883832, in COVID-19 phenotype and prognosis since data on these polymorphisms are either limited or missing.

2. Materials and Methods

2.1. Subjects

A total of 249 individuals (male/female: 172/77, median age: 38.0 years, 18.0-89.0) were enrolled in the study. Participants’ ethnicities included 166 patients of European origin (145 Greeks), 78 of Asian origin (62 Turks) and 5 of American origin (3 Americans and 2 Cubans); most patients of non-Greek origin were cruise ferry passengers and crew members who became infected with COVID-19 during March 2020 [15]. All patients were infected by the initial Wuhan-Hu-1 strain of SARS-CoV-2 between March and October 2020, as confirmed by RNA reverse transcriptase-polymerase chain reaction (RT-PCR) testing on nasopharyngeal swab specimens.
Patients were further divided into five (5) groups, according to the World Health Organization (WHO) [16] guidelines and clinical progression scale as follows: (a) asymptomatic (58, 23.3%); (b) mild (88, 35.3%), including patients who had symptoms consistent with COVID-19, such as cough, chest pain, runny nose, anosmia, ageusia, digestive symptoms, headache or myalgia; but not exhibiting signs of viral pneumonia or hypoxia; (c) moderate (26, 10.4%), including hospitalized symptomatic patients who had clinical pneumonia symptoms (fever, cough, dyspnea, tachypnea) but no signs of severe pneumonia, such as SpO2 <90% on room air necessitating oxygen supplementation by mask or nasal prong; (d) severe (34, 13.7%), including patients with clinical signs of pneumonia plus one of the following: breath rate >30/min, severe respiratory distress or SpO2 < 90% on room air; and (e) critical (43, 17.3%), including patients with complications such as respiratory failure, acute respiratory distress syndrome (ARDS), sepsis or septic shock, thromboembolism, and/or multi-organ failure, as well as patients who died. According to the aforementioned progression scale for COVID-19 severity, patients with pneumonia could be classified into either stage three which denoted moderate disease, or stage four indicating severe disease, depending on the patient’s clinical assessment. Similarly, patients experiencing respiratory distress could be categorized as either stage four, indicating severe disease or stage five, indicating critical disease, based on the severity of their medical condition, PaO2/FiO2 ratio, and the life support measures required during hospitalization. An overview of enrolled patients’ demographic and clinical characteristics is shown in Table 1.
Written informed consent was obtained from all participants. The study was carried out in accordance with the principles of the Helsinki Declaration and was approved by the ethical committee of the Faculty of Medicine, University of Thessaly (No. 2115/22.04.2020).

2.2. Molecular analyses

Genomic DNA was initially extracted from peripheral blood samples using the EXTRACTME GENOMIC DNA KIT (Blirt, Gdansk, Poland), according to manufacturer’s instructions. DNA samples were stored at -20oC until use. We further selected seven (7) functional single nucleotide polymorphisms (SNPs) in six (6) innate immune genes: TLR2 (rs5743708), TLR4 (rs4986790, rs4986791), CD14 (rs2569190), CARD8 (rs1834481), IL18 (rs2043211) and CD40 (rs1883832). SNPs were genotyped by PCR and restriction fragment length polymorphism (PCR-RLFP) assays, as follows: 100-200 ng of DNA was amplified in a 30 μL reaction mixture containing 62.5 μM of each deoxynucleoside triphosphate, 20 pmol of forward and reverse primers, 1.5 mM MgCl2, and 0.8 U of DFS-Taq DNA polymerase (Bioron GmbH, Romerberg, Germany) in a buffer provided by the manufacturer. Afterwards, PCR products were subjected to digestion with the appropriate restriction enzymes (New England Biolabs, Ipswich, UK) at 37oC overnight; the digestion products were separated on 2-3% agarose gels and visualized by ethidium bromide staining. PCR and RFLP assays were conducted in an Applied Biosystems thermal cycler (Foster City, CA, USA). The sequences of the utilized primers, conditions of PCR reactions and the restriction enzymes used for the detection of each SNP are shown in detail in Table 2. It is worth of note that the primers for the detection of TLR2-rs5743708, TLR4-rs4986790, TLR4-rs4986791, and CD40-rs1883832 polymorphisms were previously reported by us and others in the literature [17,18,19].
To confirm our PCR-RFLP results, randomly chosen wild-type, heterozygous and homozygous PCR products were purified using a PCR purification kit (Qiagen, Crawley, UK) and subsequently sequenced using an ABI Prism 310 genetic analyzer (Applied Biosystems) and a BigDye Terminator DNA sequencing kit (Applied Biosystems).

2.3. Statistical Analysis

Categorical variables are described with the use of frequency and relative frequency. Continuous variables are described with medians and interquartile range (IQR). Categorical data were analyzed with the use of Chi-square test after Yate’s correction or Fisher’s exact test. The analysis of continuous variables was conducted using the Mann-Whitney U test, as the assumption of normal distribution was violated. Data were checked for deviation from normal distribution using the Shapiro-Wilk normality test. Multivariate analysis was performed in the form of binary and ordinal logistic regression. A 5% significance level was set for all analyses. The analysis was carried out with SPSS version 29.0 (IBM Corp. Released 2021. IBM SPSS Statistics for Windows, Version 29.0. Armonk, NY: IBM Corp).

3. Results

The distribution of patients’ genotype and allele frequencies regarding the analyzed SNPs are presented in Table 3. All genotype frequencies were in Hardy-Weinberg equilibrium (Table 3). We identified that the rs1834481 polymorphism significantly affected clinical presentation of COVID-19, since rs1834481-G carriers displayed a significantly increased risk for developing pneumonia and severe disease (Table 4).
In addition, we observed that carriers of TLR2-rs5743708 and TLR4-rs4986791 SNPs displayed a 3.7- and 2.2-fold increased risk for developing pneumonia, respectively, also affecting the risk for a more severe disease; however, the statistical differences were not marginally significant (Table 4). Thus, bearing in mind that both receptors exhibit common properties, they are the best characterized pattern recognition receptors (PRRs) and shape pathogen specific innate immunity through distinct ligand binding [20]. As detailed in Table 4, we further demonstrated that carriers of either TLR2 or TLR4 polymorphism, as a whole group, displayed a significant risk for developing pneumonia and more severe disease.
Moreover, we identified that the CARD8-rs2043211 polymorphism also affected the COVID-19 phenotype, providing a protective effect against severe disease (Table 4). Conversely, the other polymorphisms did not have any impact on COVID-19 severity (p>0.05 in all cases; Table 4).
COVID-19 severity was also significantly affected by both age and the presence of comorbidities (Table 1). Indeed, elderly patients as well as patients with obesity, hypertension or diabetes mellitus were more likely to present severe COVID-19-associated clinical symptoms, such as pneumonia and respiratory distress syndrome (Table 1).
Therefore, we performed multivariate analysis to identify risk factors affecting the clinical presentation and outcome of COVID-19, including all parameters that were ascertained in univariate analyses. Elderly patients with obesity and/or diabetes mellitus were more vulnerable to developing pneumonia and respiratory distress syndrome, following SARS-CoV-2 infection (Table 5). Moreover, individuals with a medical history of hematological and/or chronic immune-mediated diseases were more susceptible to developing pneumonia and severe disease, while patients with hypertension displayed an approximately 3-fold increased probability of developing more severe disease (Table 5). Additionally, we observed that the presence of the IL18-rs1834481 polymorphism was an independent risk factor for developing pneumonia, also displaying a 1.7-fold increased risk for more severe disease; although the latter was not marginally reached to statistical significance (Table 5). Finally, carriers of either TLR2-rs5743708 or TLR4-rs4986791 polymorphism also displayed a significantly increased risk for developing pneumonia and more severe disease (Table 5).

4. Discussion

COVID-19 outcome is the result of a complex interplay between viral components and host proteins which are involved in the immune response to the virus [1,11,12]. However, soon after the COVID-19 outbreak it became clear that the disease heterogeneity may be partially explained by the host’s genetic background [1,2,3,4,5,6,7,8,9,10,11,12]. In this context, the main objective of our study was to examine whether functional polymorphisms in genes coding for crucial mediators of innate immune response – namely TLR2, TLR4, CD14, CARD8, IL-18 and CD40 – can influence COVID-19 clinical severity in a multinational population of patients upon SARS-CoV-2 infection.
One of the principal findings of our study was that the IL18-rs1834481 polymorphism appeared to have a strong impact on COVID-19 clinical outcome, since heterozygotes for the rs1834481 polymorphism were at a higher risk of developing pneumonia (both in univariate and multivariate analyses). IL-18 is a proinflammatory cytokine released downstream of the activation of innate immune pathways, resulting in the formation of inflammasomes; it has already been correlated with a poor COVID-19 prognosis, since increased IL-18 levels have been found in patients with severe disease [21,22].
The rs1834481 is an intronic SNP (located into intron 3) of the IL18 gene [23]. Previous studies have demonstrated that along with other SNPs (namely rs2115763 and rs5744256), rs1834481 is a functional polymorphism affecting IL18 expression, and the presence of G alleles have significantly been associated with lower IL-18 serum levels [24,25]. The frequency of rs1834481-G allele is approximately 22.2% globally (https://www.ncbi.nlm.nih.gov/snp/rs1834481#frequency_tab), and its presence appears to protect carriers from coronary artery disease [26], non-small cell lung cancer [27] or human papilloma virus (HPV) persistence [28]. However, it has been associated with a poor prognosis of colon cancer [29]. According to our results, the rs1834481-G allele which has been associated with lower IL-18 levels after an immune stimulation, predisposes to more severe disease following SARS-CoV-2 infection. This may suggest that it could be used as a molecular predictor for COVID-19 severity and outcome. Considering that increased IL-18 levels were found in patients with a more severe COVID-19 phenotype, we could speculate that an inappropriate initial response against SARS-CoV-2 in rs1834481-G allele carriers may predispose to virus persistence, and subsequently to dysregulated immune responses. This could result in subsequent hyperinflammation events, and finally to more severe disease and poor outcome.
To the best of our knowledge, this is the first study in the literature indicating a strong association between IL18-rs1834481 and clinical severity of COVID-19. After a broad literature review, we found only one additional study analyzing the contribution of rs1834481 in COVID-19 severity. Leal et al. analyzed a larger cohort of Brazilian patients infected by SARS-CoV-2, but they did not observe any significant association of rs1834481 with COVID-19 severity; however, they used a different methodological approach comparing only patients with severe and mild disease, while they did not use, at least, the WHO classification as in our study [30].
Additionally, we clearly demonstrated that carriers of either TLR2-rs5743708 or TLR4-rs4986791 polymorphisms, as a whole group, displayed a significant risk for pneumonia development and more severe disease (Table 4 and Table 5). The aforementioned polymorphisms have been associated with receptor hypo-responsiveness and susceptibility to several infectious diseases (reviewed by Schroder and Schumann) [31]. Accumulating evidence demonstrates that the viral S protein of SARS-CoV-2 may bind to and activate TLR4, while TLR2 was shown to interact with the envelope protein of the virus [32,33,34]. Therefore, our data further supports the results of two previous studies from Egypt, where a positive correlation between the aforementioned polymorphisms and COVID-19 severity was observed [35,36], suggesting that both SNPs could be used as molecular predictors of COVID-19 prognosis and outcome.
Among the other innate immune gene polymorphisms analyzed, the CARD8-rs2043211 polymorphism was found to exhibit a statistically significant association. CARD8 negatively regulates inflammasome activation via inhibiting caspase 1 and NF-κB-mediated signaling pathways [37]. Presence of the polymorphism results in a substitution of phenylalanine by isoleucine at position 102 in the classical transcript variant 1 of CARD8 gene (NM_001184900.3), although it could lead to a termination codon in other splicing variants of the gene (transcript variants 12, 14 and 17-19; https://www.ncbi.nlm.nih.gov/snp/rs2043211#variant_details). Similar to the study conducted by de Sá et al. [38], in our cohort we observed a protective effect of CARD8-rs2043211 polymorphism in COVID-19 severity (Table 4), however this finding was lost in multivariate analysis (Table 5).
CD14 acts as a coreceptor for both of the aforementioned TLRs, facilitating the binding of several ligands (pathogen associated molecular patterns-PAMPS and danger associated molecular patterns-DAMPS) to their receptors [39]. Although Pati et al. suggested that CD14-rs2569190 which is located in the promoter region of the CD14 gene, may be linked to SARS-CoV-2-assocociated mortality in the European population [40], our findings do not support a positive association of this polymorphism with COVID-19 severity and outcome.
Finally, CD40 is a pivotal immune checkpoint molecule expressed by professional antigen-presenting cells, acting as a bridge between innate and adaptive immune responses [41]. To the best of our knowledge, this is the first study analyzing the contribution of the CD40-rs1883832 polymorphism with clinical presentation and severity of SARS-CoV-2 infection. Our group has recently reported that upon vaccination, this polymorphism affects the IgA antibody responses against SARS-CoV-2 [19], while it has also been correlated with a susceptibility to sepsis in a Chinese population [42]. However, we did not find a significant association between the CD40-rs1883832 polymorphism and the clinical course of COVID-19 in our study.

5. Conclusions

In conclusion, our study provides clear evidence that the host’s genetic background can significantly affect COVID-19 clinical presentation and severity, underlining the role of immunogenetics in COVID-19 prognosis. The identification of molecular markers such as the IL18-rs1834481, TLR2-rs5743708 and TLR4-rs4986791 polymorphisms, and the incorporation of this knowledge into clinical practice may allow for better stratification of patients for an improved therapeutic intervention.

Author Contributions

Conceptualization, C.H. and M.S.; Methodology, M.S., E.B., I.V., V.P., N.G., G.G., E.A., S.A., E.B., V.L., O.T., I.M., M.T., S.S., V.A.M., F.K., S.M., and G.D.; Validation, M.S. and C.H.; Formal analysis, M.S.; C.H. and A.D.; Investigation, E.B., I.V., V.P., N.G., G.G., E.A., S.A., E.B., V.L., O.T., I.M., M.T. and S.S.; Data curation, M.S., C.H., and E.B.; Writing—original draft preparation, M.S.; Writing—review and editing, M.S., C.H., G.G., S.M. and G.D.; Funding acquisition, M.S.; Supervision, M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. It was financed by the standard budget provided by the Research Committee of the University of Thessaly.

Institutional Review Board Statement

The study was carried out in accordance with the principles of the Helsinki Declaration and was approved by the ethical committee of the Faculty of Medicine, University of Thessaly (No. 2115/22.04.2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data that support the findings of this study are available from the corresponding authors upon request.

Acknowledgments

The authors express their sincere gratitude to all participants for blood donation and Lemonia Anagnostopoulou for her invaluable assistance in providing English corrections.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jin, Y.; Yang, H.; Ji, W.; Wu, W.; Chen, S.; Zhang, W.; Duan, G. Virology, Epidemiology, Pathogenesis, and Control of COVID-19. Viruses 2020, 12, 372. [Google Scholar] [CrossRef] [PubMed]
  2. Kompaniyets, L.; Pennington, A. F.; Goodman, A. B.; Rosenblum, H. G.; Belay, B.; Ko, J. Y.; Chevinsky, J. R.; Schieber, L. Z.; Summers, A. D.; Lavery, A. M.; et al. Underlying Medical Conditions and Severe Illness Among 540,667 Adults Hospitalized With COVID-19, March 2020–March 2021. Prev Chronic Dis 2021, 18, E66. [Google Scholar] [CrossRef] [PubMed]
  3. Speletas, M.; Dadouli, K.; Syrakouli, A.; Gatselis, N.; Germanidis, G.; Mouchtouri, V. A.; Koulas, I.; Samakidou, A.; Nikolaidou, A.; Stefos, A.; et al. MBL Deficiency-Causing B Allele (Rs1800450) as a Risk Factor for Severe COVID-19. Immunobiology 2021, 226, 152136. [Google Scholar] [CrossRef] [PubMed]
  4. Dalekos, G. N.; Stefos, A.; Georgiadou, S.; Lygoura, V.; Michail, A.; Ntaios, G.; Samakidou, A.; Giannoulis, G.; Gabeta, S.; Vlychou, M.; et al. Lessons from Pathophysiology: Use of Individualized Combination Treatments with Immune Interventional Agents to Tackle Severe Respiratory Failure in Patients with COVID-19. Eur J Intern Med 2021, 88, 52–62. [Google Scholar] [CrossRef]
  5. Moon, A. M.; Webb, G. J.; García-Juárez, I.; Kulkarni, A. V.; Adali, G.; Wong, D. K.; Lusina, B.; Dalekos, G. N.; Masson, S.; Shore, B. M.; et al. SARS-CoV-2 Infections Among Patients With Liver Disease and Liver Transplantation Who Received COVID-19 Vaccination. Hepatol Commun 2021, 6, 889–897. [Google Scholar] [CrossRef]
  6. Gatselis, N. K.; Lygoura, V.; Lyberopoulou, A.; Giannoulis, G.; Samakidou, A.; Vaiou, A.; Vatidis, G.; Antoniou, K.; Stefos, A.; Georgiadou, S.; et al. Soluble IL-2R Levels at Baseline Predict the Development of Severe Respiratory Failure and Mortality in COVID-19 Patients. Viruses 2022, 14, 787. [Google Scholar] [CrossRef]
  7. Ovsyannikova, I. G.; Haralambieva, I. H.; Crooke, S. N.; Poland, G. A.; Kennedy, R. B. The Role of Host Genetics in the Immune Response to SARS-CoV-2 and COVID-19 Sus-ceptibility and Severity. Immunol Rev 2020, 296, 205–219. [Google Scholar] [CrossRef]
  8. Adli, A.; Rahimi, M.; Khodaie, R.; Hashemzaei, N.; Hosseini, S. M. Role of Genetic Variants and Host Polymorphisms on COVID-19: From Viral Entrance Mechanisms to Immunological Reactions. J Med Virol 2022, 94, 1846–1865. [Google Scholar] [CrossRef]
  9. Zepeda-Cervantes, J.; Martínez-Flores, D.; Ramírez-Jarquín, J. O.; Tecalco-Cruz, Á. C.; Alavez-Pérez, N. S.; Vaca, L.; Sarmiento-Silva, R. E. Implications of the Immune Polymorphisms of the Host and the Genetic Variability of SARS-CoV-2 in the Development of COVID-19. Viruses 2022, 14, 94. [Google Scholar] [CrossRef]
  10. Thair, S. A.; He, Y. D.; Hasin-Brumshtein, Y.; Sakaram, S.; Pandya, R.; Toh, J.; Rawling, D.; Remmel, M.; Coyle, S.; Dalekos, G. N.; et al. Transcriptomic Similarities and Differences in Host Response between SARS-CoV-2 and Other Viral Infections. iScience 2020, 24, 101947. [Google Scholar] [CrossRef]
  11. El-Hefnawy, S. M.; Eid, H. A.; Mostafa, R. G.; Soliman, S. S.; Omar, T. A.; Azmy, R. M. COVID-19 Susceptibility, Severity, Clinical Outcome and Toll-like Receptor (7) MRNA Expression Driven by TLR7 Gene Polymorphism (Rs3853839) in Middle-Aged Individuals without Previous Comorbidities. Gene Rep 2022, 27, 101612. [Google Scholar] [CrossRef] [PubMed]
  12. Pojero, F.; Candore, G.; Caruso, C.; Di Bona, D.; Groneberg, D. A.; Ligotti, M. E.; Accardi, G.; Aiello, A. The Role of Immunogenetics in COVID-19. Int J Mol Sci 2021, 22, 2636. [Google Scholar] [CrossRef]
  13. Diamond, M. S.; Kanneganti, T.-D. Innate Immunity: The First Line of Defense against SARS-CoV-2. Nat Immunol 2022, 23, 165–176. [Google Scholar] [CrossRef] [PubMed]
  14. Paludan, S.R.; Mogensen, T.H. Innate immunological pathways in COVID-19 pathogenesis. Sci Immunol 2022, 7, eabm5505. [Google Scholar] [CrossRef] [PubMed]
  15. Hatzianastasiou, S.; Mouchtouri, V. A.; Pavli, A.; Tseroni, M.; Sapounas, S.; Vasileiou, C.; Dadouli, K.; Kyritsi, M.; Koureas, M.; Prezerakos, P.; et al. COVID-19 Outbreak on a Passenger Ship and Assessment of Response Measures, Greece, 2020. Emerg Infect Dis 2021, 27, 1927–1930. [Google Scholar] [CrossRef] [PubMed]
  16. World Health Organization. (2022). Clinical management of COVID-19: living guideline, 15 September 2022. World Health Organization. https://apps.who.int/iris/handle/10665/362783. License: CC BY-NC-SA 3.0 IGO. 15 September.
  17. Lorenz, E.; Frees, K. L.; Schwartz, D. A. Determination of the TLR4 Genotype Using Allele-Specific PCR. BioTechniques 2001, 31, 22–24. [Google Scholar] [CrossRef]
  18. Liadaki, K.; Petinaki, E.; Skoulakis, C.; Tsirevelou, P.; Klapsa, D.; Germenis, AE.; Speletas, M. Toll-like receptor 4 gene (TLR4), but not TLR2, polymorphisms modify the risk of tonsillar disease due to Streptococcus pyogenes and Haemophilus influenzae. Clin Vaccine Immunol 2011, 18, 217–222. [Google Scholar] [CrossRef]
  19. Speletas, M.; Bakaros, E.; Peristeri, A.-M.; Voulgaridi, I.; Sarrou, S.; Paliatsa, V.; Nasika, A.; Tseroni, M.; Anagnostopoulos, L.; Theodoridou, K.; et al. The Rs1883832 Polymorphism (CD40-1C>T) Affects the Intensity of IgA Responses after BNT162b2 Vaccination. Int J Mol Sci 2022, 23, 14056. [Google Scholar] [CrossRef]
  20. Mukherjee, S.; Karmakar, S.; Babu, S. P. S. TLR2 and TLR4 Mediated Host Immune Responses in Major Infectious Diseases: A Review. Braz J Infect Dis 2016, 20, 193–204. [Google Scholar] [CrossRef]
  21. Rodrigues, T. S.; de Sá, K. S. G.; Ishimoto, A. Y.; Becerra, A.; Oliveira, S.; Almeida, L.; Gonçalves, A. V.; Perucello, D. B.; Andrade, W. A.; Castro, R.; et al. Inflammasomes Are Activated in Response to SARS-CoV-2 Infection and Are Associated with COVID-19 Severity in Patients. J Exp Med 2020, 218, e20201707. [Google Scholar] [CrossRef]
  22. Satış, H.; Özger, H. S.; Aysert Yıldız, P.; Hızel, K.; Gulbahar, Ö.; Erbaş, G.; Aygencel, G.; Guzel Tunccan, O.; Öztürk, M. A.; Dizbay, M.; et al. Prognostic Value of Interleukin-18 and Its Association with Other Inflammatory Markers and Disease Severity in COVID-19. Cytokine 2021, 137, 155302. [Google Scholar] [CrossRef] [PubMed]
  23. Nolan, K. F.; Greaves, D. R.; Waldmann, H. The Human Interleukin 18 GeneIL18Maps to 11q22.2–Q22.3, Closely Linked to the DRD2 Gene Locus and Distinct from Mapped IDDM Loci. Genomics 1998, 51, 161–163. [Google Scholar] [CrossRef] [PubMed]
  24. He, M.; Cornelis, M. C.; Kraft, P.; van Dam, R. M.; Sun, Q.; Laurie, C. C.; Mirel, D. B.; Chasman, D. I.; Ridker, P. M.; Hunter, D. J.; et al. Genome-Wide Association Study Identifies Variants at IL18-BCO2 Locus Associated with Interleukin-18 Levels. Arterioscler Thromb Vasc Biol 2010, 30, 885–890. [Google Scholar] [CrossRef]
  25. Vergara, C.; Thio, C.; Latanich, R.; Cox, A. L.; Kirk, G. D.; Mehta, S. H.; Busch, M.; Murphy, E. L.; Villacres, M. C.; Peters, M. G.; et al. Genetic Basis for Variation in Plasma IL-18 Levels in Persons with Chronic Hepatitis C Virus and Human Immunodeficiency Virus-1 Infections. Genes Immun 2017, 18, 82–87. [Google Scholar] [CrossRef]
  26. Koch, W.; Wolferstetter, H.; Schatke, A.; Schömig, A.; Kastrati, A. Interleukin 18 Gene Variation and Risk of Acute Myocardial Infarction. Cytokine 2011, 56, 786–791. [Google Scholar] [CrossRef]
  27. Zhang, W.; Dang, S.; Zhang, G.; He, H.; Wen, X. Genetic Polymorphisms of IL-10, IL-18 and IL12B Are Associated with Risk of Non-Small Cell Lung Cancer in a Chinese Han Population. International Immunopharmacology 2019, 77, 105938. [Google Scholar] [CrossRef] [PubMed]
  28. Pontillo, A.; Bricher, P.; Leal, V. n. c.; Lima, S.; Souza, P. r. e.; Crovella, S. Role of Inflammasome Genetics in Susceptibility to HPV Infection and Cervical Cancer Development. Journal of Medical Virology 2016, 88, 1646–1651. [Google Scholar] [CrossRef]
  29. Cambui, R. A. G.; do Espírito Santo, G. F.; Fernandes, F. P.; Leal, V. N. C.; Galera, B. B.; Fávaro, E. G. P.; Rizzo, L. A.; Elias, R. M.; Pontillo, A. Double-Edged Sword of Inflammasome Genetics in Colorectal Cancer Prognosis. Clin Immunol 2020, 213, 108373. [Google Scholar] [CrossRef]
  30. Leal, V. N. C.; Paulino, L. M.; Cambui, R. A. G.; Zupelli, T. G.; Yamada, S. M.; Oliveira, L. A. T.; Dutra, V. de F.; Bub, C. B.; Sakashita, A. M.; Yokoyama, A. P. H.; et al. A Common Variant Close to the “Tripwire” Linker Region of NLRP1 Contributes to Severe COVID-19. Inflamm Res 2022, 1–8. [Google Scholar] [CrossRef]
  31. Schröder, N. W.; Schumann, R. R. Single Nucleotide Polymorphisms of Toll-like Receptors and Susceptibility to Infectious Disease. The Lancet Infectious Diseases 2005, 5, 156–164. [Google Scholar] [CrossRef]
  32. Choudhury, A.; Mukherjee, S. In Silico Studies on the Comparative Characterization of the Interactions of SARS-CoV-2 Spike Glycoprotein with ACE-2 Receptor Homologs and Human TLRs. J Med Virol 2020, 92, 2105–2113. [Google Scholar] [CrossRef] [PubMed]
  33. Zhao, Y.; Kuang, M.; Li, J.; Zhu, L.; Jia, Z.; Guo, X.; Hu, Y.; Kong, J.; Yin, H.; Wang, X.; et al. Publisher Correction: SARS-CoV-2 Spike Protein Interacts with and Activates TLR4. Cell Res 2021, 31, 825. [Google Scholar] [CrossRef] [PubMed]
  34. Zheng, M.; Karki, R.; Williams, E. P.; Yang, D.; Fitzpatrick, E.; Vogel, P.; Jonsson, C. B.; Kanneganti, T.-D. TLR2 Senses the SARS-CoV-2 Envelope Protein to Produce Inflammatory Cytokines. Nat Immunol 2021, 22, 829–838. [Google Scholar] [CrossRef] [PubMed]
  35. Alhabibi, A. M.; Hassan, A. S.; Abd Elbaky, N. M.; Eid, H. A.; Khalifa, M. A. A. A.; Wahab, M. A.; Althoqapy, A. A.; Abdou, A. E.; Zakaria, D. M.; Nassef, E. M.; et al. Impact of Toll-Like Receptor 2 and 9 Gene Polymorphisms on COVID-19: Susceptibility, Severity, and Thrombosis. J Inflamm Res 2023, 16, 665–675. [Google Scholar] [CrossRef]
  36. Taha, S. I.; Shata, A. K.; Baioumy, S. A.; Fouad, S. H.; Anis, S. G.; Mossad, I. M.; Moustafa, N. M.; Abdou, D. M.; Youssef, M. K. Toll-Like Receptor 4 Polymorphisms (896A/G and 1196C/T) as an Indicator of COVID-19 Severity in a Convenience Sample of Egyptian Patients. J Inflamm Res 2021, 14, 6293–6303. [Google Scholar] [CrossRef]
  37. Razmara, M.; Srinivasula, S. M.; Wang, L.; Poyet, J.-L.; Geddes, B. J.; DiStefano, P. S.; Bertin, J.; Alnemri, E. S. CARD-8 Protein, a New CARD Family Member That Regulates Caspase-1 Activation and Apoptosis. Journal of Biological Chemistry 2002, 277, 13952–13958. [Google Scholar] [CrossRef]
  38. de Sá, N. B. R.; Neira-Goulart, M.; Ribeiro-Alves, M.; Perazzo, H.; Geraldo, K. M.; Ribeiro, M. P. D.; Cardoso, S. W.; Grinsztejn, B.; Veloso, V. G.; Capão, A.; et al. Inflammasome Genetic Variants Are Associated with Protection to Clinical Severity of COVID-19 among Patients from Rio de Janeiro, Brazil. Biomed Res Int 2022, 9082455. [Google Scholar] [CrossRef]
  39. Di Gioia, M.; Zanoni, I. Toll-like Receptor Co-Receptors as Master Regulators of the Immune Response. Mol Immunol 2015, 63, 143–152. [Google Scholar] [CrossRef]
  40. Pati, A.; Padhi, S.; Panda, D.; Suvankar, S.; Panda, A. K. A Cluster of Differentiation 14 (CD14) Polymorphism (C-159T Rs2569190) Is Associated with SARS-CoV-2 Infection and Mortality in the European Population. J Infect Dis 2021, jiab180. [Google Scholar] [CrossRef]
  41. Sun, L.; Wang, X.; Saredy, J.; Yuan, Z.; Yang, X.; Wang, H. Innate-Adaptive Immunity Interplay and Redox Regulation in Immune Response. Redox Biol 2020, 37, 101759. [Google Scholar] [CrossRef]
  42. Liu, Z.-L.; Hu, J.; Xiao, X.-F.; Peng, Y.; Zhao, S.-P.; Xiao, X.-Z.; Yang, M.-S. The CD40 Rs1883832 Polymorphism Affects Sepsis Susceptibility and SCD40L Levels. Biomed Res Int 2018, 7497314. [Google Scholar] [CrossRef] [PubMed]
Table 1. Demographic and clinical characteristics of the study patients. Patients are divided into five (5) categories according to guidelines and the clinical progression scale provided by the World Health Organization (WHO), as detailed in the text.
Table 1. Demographic and clinical characteristics of the study patients. Patients are divided into five (5) categories according to guidelines and the clinical progression scale provided by the World Health Organization (WHO), as detailed in the text.
Total Asymptomatic disease Mild disease Moderate disease Severe disease Critical disease p value #
No 249 58 88 26 34 43 -
Sex (male/female) 172/77 40/18 60/28 16/10 27/7 29/14 0.648 (C)
Age 42.8 ± 18.5 35.3 ± 14.3 31.3 ± 11.6 48.2 ± 17.3 58.1 ± 14.8 60.8 ± 14.6 <0.001 (K-W)
Comorbidities (n, %) 89, 35.7 6, 10.3 9, 10.2 15, 57.7 23, 67.6 36, 83.7 <0.001 (C)
Obesity (n, %) 28, 11.2 0, 0.0 2, 2.3 4, 15.4 4, 11.8 18, 41.9 <0.001 (C)
Hypertension (n, %) 45, 18.1 1, 1.7 2, 2.3 5, 19.2 12, 35.3 25, 58.1 <0.001 (C)
Dyslipidemia (n, %) 17, 6.8 1, 1.7 0, 0.0 5, 19.2 4, 11.8 7, 16.3 <0.001 (C)
Diabetes mellitus (n, %) 20, 8.0 0, 0.0 2, 2.3 1, 3.8 8, 23.5 9, 20.9 <0.001 (C)
Thyroid disease (n, %) 10, 4.0 0, 0.0 3, 3.4 3, 11.5 1, 2.9 3, 7.0 0.116 (C)
CHD (n, %) 11, 4.4 0, 0.0 0, 0.0 1, 3.8 4, 11.8 6, 14.0 <0.001 (C)
CVD (n, %) 8, 3.2 0, 0.0 0, 0.0 1, 3.8 0, 0.0 7, 16.3 <0.001 (C)
CRD (n, %) 18, 7.2 1, 1.7 1, 1.1 0, 0.0 6, 17.6 10, 23.3 <0.001 (C)
Chronic liver and/or renal disease (n, %) 6, 2.4 1, 1.7 0, 0.0 2, 7.7 1, 2.9 2, 4.7 0.176 (C)
Malignancies (n, %) 6, 2.4 0, 0.0 2, 2.3 0, 0.0 1, 2.9 3, 7.0 0.204 (C)
Chronic immune and/or hematologic diseases (n, %) ^ 6, 2.4 1, 1.7 0, 0.0 0, 0.0 3, 8.8 2, 4.7 0.044 (C)
Comorbidity-others (n, %) * 15, 6.0 3, 5.2 0, 0.0 1, 3.8 5, 14.7 6, 14.0 0.004 (C)
Abbreviations: CHD, chronic heart disease (atrial fibrillation, heart failure, arrhythmias, prosthetic valve); CRD, chronic respiratory diseases (chronic obstructive pulmonary disease, asthma, sleep apnea); CVD, cerebrovascular disease (including also chronic venous disease); SD, standard deviation. ^ Chronic immune and/or hematologic diseases included patients with systemic lupus erythematosus (1), autoinflammatory disease (1), IgA deficiency (1), beta-thalassemia, (1), sickle-cell disease (1), splenectomy (unknown etiology) (1); * Comorbidity-others included patients with history of transplantation (liver or kidney), history of severe infections (meningitis, flu due to H1N1, recurrent urinary infections), history of major surgery (partial gastrectomy, nephrectomy, partial colectomy), dementia, depression. # Statistical analysis: C, Chi-square; F. Fisher’s exact test; K-W = Kruskal Wallis test; M-W, Mann-Whitney U test.
Table 2. PCR-RFLP protocols for the detection of single nucleotide polymorphisms (SNP) of innate immune genes of the study.
Table 2. PCR-RFLP protocols for the detection of single nucleotide polymorphisms (SNP) of innate immune genes of the study.
Gene SNP Primer sequences ^ PCR conditions PCR product (bp) Digestion Enzyme Expected RFLP fragments (bp)*
TLR2 rs5743708 F: 5′-TATGGTCCAGGAGCTGGAGA-3′ 94oC 2 min, 38 cycles (94oC 30 sec, 64oC 30 sec, 72oC 45 sec), 72oC 5 min 430 SfcI 430 (wt/GG), 430+307+123 (heter/GA), 307+123 (homo/AA)
R: 5′-TGACATAAAGATCCCAACTAGACAA-3′
TLR4 rs4986790 F: 5’-GATTAGCATACTTAGACTACTACCTCCATG-3’ 94oC 2 min, 32 cycles (94oC 30 sec, 57oC 30 sec, 72oC 60 sec), 72oC 5 min 249 NcoI 249 (wt/AA), 249+218+31 (heter/AG), 218+31 (homo/GG)
R: 5’-GATCAACTTCTGAAAAAGCATTCCCAC-3’
rs4986791 F: 5’-GGTTGCTGTTCTCAAAGTGATTTTGGGAGAA-3’ 94oC 2 min, 32 cycles (94oC 30 sec, 57oC 30 sec, 72oC 60 sec), 72oC 5 min 404 HinfI 404 (wt/CC), 404+375+29 (heter/CT), 375+29 (homo/TT)
R: 5’-ACCTGAAGACTGGAGAGTGAGTTAAATGC-3’
CD14 rs2569190 F: 5’-GCCAACAGATGAGGTTCACA-3’ 94oC 2 min, 30 cycles (94oC 30 sec, 60oC 30 sec, 72oC 30 sec), 72oC 1 min 641 HaeIII 308+174 (wt/CC), 522+308+174 (heter/CG), 522 (homo/GG)
R: 5’-GTTCGACCCCAAGACCCTAC-3’
CARD8 rs2043211 F: 5’-AAGAGCAGCATATAGTGGGAAA-3’ 94oC 2 min, 30 cycles (94oC 30 sec, 57oC 30 sec, 72oC 30 sec), 72oC 1 min 475 MboII 151+324 (wt/AA), 475+324+151 (heter/AT), 475 (homo/TT)
R: 5’-TGAGTTCGATGAAAAACACCC-3’
IL18 rs1834481 F: 5’-TTCTGTTGTAGGGGGAGCGT-3’ 94oC 2 min, 30 cycles (94oC 30 sec, 60oC 30 sec, 72oC 30 sec), 72oC 1 min 522 HaeIII 308+174 (wt/CC), 522+308+174 (heter/CG), 522 (homo/GG)
R: 5’-CATCAGGTTGCTTCCAGCCA-3’
CD40 rs1883832 F: 5’-ATAGGTGGACCGCGATTGGT-3’ 95oC 10 min, 32 cycles (94oC 30 sec, 58oC 30 sec, 72oC 30 sec), 72oC 10 min 216 NcoI 126+90 (wt/CC), 216+126+90 (heter/CT), 216 (homo/TT)
R: 5’-TCCCAACTCCCGTCTGGT-3’
Abbreviations: bp, base pairs; F, Forward; RFLP, restriction fragment length polymorphism; R, Reverse; * Within brackets we present the type of genotype (wt, wild-type; het, heterozygote; hom, homozygote) and the corresponding nucleotides of each SNP.
Table 3. Genotype and allele frequencies of analyzed single nucleotide polymorphisms (SNPs) of the study.
Table 3. Genotype and allele frequencies of analyzed single nucleotide polymorphisms (SNPs) of the study.
Gene SNP Genotype Allele frequency (n, %) HWE
WT (n, %) Heter (n, %) Homo (n, %)
TLR2 rs5743708 (g.25877G>A) 239 (96.0) 10 (4.0) 0 (0.0) G (488, 98.0) A (10, 2.0) 0.75
TLR4 rs4986790 (g.13622A>G) 229 (92.0) 19 (7.6) 1 (0.4) A (477, 95.8) G (21, 4.2) 0.38
rs4986791 (g.13922C>T) 230 (92.4) 18 (7.2) 1 (0.4) C (478, 96.0) T (20, 4.0) 0.33
CD14 rs2569190 (g.5371T>C) 61 (24.5) 123 (49.4) 65 (26.1) T (245, 49.2) C (253, 50.8) 0.85
IL18 rs1834481 (g.16014C>G) 187 (75.1) 54 (21.7) 8 (31.2) C (428, 85.9) G (70, 14.1) 0.11
CARD8 rs2043211 (g.26821A>T) 99 (37.8) 119 (47.8) 31 (12.4) A (317, 63.7) T (181, 36.3) 0.60
CD40 rs1883832 (g.5077C>T) 108 (43.4) 119 (47.8) 22 (8.8) C (335, 67.3) T (163, 32.7) 0.18
Abbreviations: Heter, heterozygous; Homo, homozygous; HWE, Hardy-Weiberg equilibrium; WT, wild-type.
Table 4. Innate immune gene single nucleotide polymorphisms (SNPs) and clinical phenotype of COVID-19.
Table 4. Innate immune gene single nucleotide polymorphisms (SNPs) and clinical phenotype of COVID-19.
Gene SNP Genotype (N, %)^ Pneumonia (yes vs no) Respiratory distress (yes vs no) Disease severity (according to the WHO classification) *
p OR (95%CI) p OR (95%CI) p OR (95%CI)
TLR2 rs5743708 WT (239, 96.0) Heter+Homo (10, 4.0) 0.061 3.73 (0.94 – 14.78) 0.255 2.13 (0.58 – 7.81) 0.103 2.57 (0.83 – 8.03)
TLR4 rs4986790 WT (229, 92.0) Heter+Homo (20, 8.0) 0.329 1.75 (0.63 – 3.93) 0.789 0.75 (0.24 – 2.35) 0.305 1.53 (0.68 – 3.47)
rs4986791 WT (230, 92.4) Heter+Homo (19, 7.6) 0.093 2.22 (0.86 – 5.73) 0.745 1.11 (0.38 – 3.22) 0.058 2.14 (0.92 – 4.94)
CD14 rs2569190 WT (65, 26.1) Heter+Homo (184, 73.9) 0.804 1.08 (0.60 – 1.92) 0.486 0.80 (0.42 – 1.51) 0.416 0.81 (0.49 – 1.34)
IL18 rs1834481 WT (187, 75.1) Heter+Homo (62, 24.9) 0.012 2.09 (1.17 – 3.73) 0.101 1.70 (0.90 – 3.20) 0.008 2.01 (1.20 – 3.38)
CARD8 rs2043211 WT (99, 39.8) Heter+Homo (150, 60.2) 0.220 0.72 (0.43 – 1.21) 0.084 0.72 (0.43– 1.21) 0.031 0.60 (0.38 – 0.96)
CD40 rs1883832 WT (108, 43.4) Heter+Homo (141, 56.6) 0.665 1.14 (0.63 – 2.05) 0.442 1.22 (0.73 – 2.05) 0.368 1.23 (0.78 – 1.93)
TLR2 or/and TLR4 rs5743708 or/and rs4986790 WT (223, 89.6) Heter+Homo(26, 10.4#) 0.020 2.70 (1.17 – 6.22) 0.432 1.43 (0.59 – 3.47) 0.022 2.35 (1.13 – 4.88)
Abbreviations: CO, confidence intervals; Heter, heterozygous; Homo, homozygous; OR, odds ratio; n.a., not-available (due to the low number of such groups of patients); WT, wild-type. The bold text highlights the statistically significant associations. ^ WT genotype served as the reference one for statistical analyses; * Disease severity, as described in Material and Methods section; # Three (3) individuals were carriers of both polymorphisms (rs5743708 and rs4986790) in heterozygous state.
Table 5. Multivariate analysis of risk factors for COVID-19 phenotype and prognosis.
Table 5. Multivariate analysis of risk factors for COVID-19 phenotype and prognosis.
Parameter Pneumonia Respiratory distress Disease severity (according to the WHO classification)*
p OR (95%CI) p OR (95%CI) p OR (95%CI)
Age <0.001 1.10 (1.06 - 1.13) <0.001 1.06 (1.02 - 1.09) <0.001 1.04 (1.02 - 1.06)
Sex (female/male) 0.723 0.84 (0.31 - 2.24) 0.118 0.42 (0.14 - 1.24) 0.481 0.82 (0.47 - 1.43)
Obesity 0.001 20.35 (3.46 - 119.68) 0.002 7.38 (2.13 - 25.61) <0.001 8.18 (3.08 - 21.73)
Hypertension 0.909 0.91 (0.19 - 4.34) 0.085 2.89 (0.86 – 9.65) 0.021 3.05 (1.18 – 7.87)
Diabetes mellitus 0.029 12.15 (1.30 - 113.81) 0.013 6.59 (1.48 - 29.36) 0.039 3.32 (1.06 - 10.38)
Dyslipidemia 0.195 4.46 (0.46 – 42.86) 0.813 1.21 (0.25 – 5.82) 0.801 1.16 (0.37 - 3.66)
Thyroid disease 0.323 2.85 (0.36 - 22.79) 0.489 2.05 (0.27 – 15.70) 0.220 2.29 (0.61 - 8.59)
CHD 0.999 - 0.305 0.37 (0.05 - 2.49) 0.801 0.82 (0.17 - 4.01)
CVD 0.999 - 0.172 5.61 (0.47 – 66.54) 0.055 9.26 (0.96 – 89.84)
CRD 0.775 1.31 (0.20 – 8.56) 0.307 2.05 (0.52 - 8.09) 0.058 3.01 (0.96 - 9.39)
Chronic liver and/or renal disease 0.754 1.54 (0.11 – 22.85) 0.749 0.67 (0.06 – 7.58) 0.917 0.91 (0.17 - 4.98)
Chronic immune and/or hematologic diseases 0.002 57.89 (4.30 - 780.17) 0.216 4.63 (0.41 – 52.35) 0.004 11.92 (2.19 – 64.82)
Malignancies 0.663 0.59 (0.06 - 6.28) 0.199 4.72 (0.44 – 50.53) 0.448 2.00 (0.33 - 12.05)
Comorbidity-others^ 0.897 1.13 (0.18 - 6.28) 0.892 1.12 (0.23 - 5.53) 0.466 1.55 (0.48 - 4.99)
TLR2-rs5743708 and/or TLR4-rs4986791 0.033 3.62 (1.11 - 11.83) 0.897 1.09 (0.31 - 3.83) 0.026 2.54 (1.12 - 5.79)
IL18-rs1834481 0.026 2.88 (1.13 - 7.32) 0.607 1.28 (0.50 - 3.26) 0.073 1.69 (0.95 – 2.99)
CARD8-rs2043211 0.751 1.15 (0.49 - 2.72) 0.210 0.57 (0.24 - 1.37) 0.204 0.72 (0.44 - 1.19)
Abbreviations: CHD, chronic heart disease (as in Table 1); CO, confidence intervals; CRD, chronic respiratory diseases (as in Table 1); CVD, cerebrovascular disease (as in Table 1); Heter, heterozygous; Homo, homozygous; WT, wild-type. The bold text highlights the statistically significant associations. ^ as described in Table 1; * Disease severity, as described in the “Material and Methods” section.
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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated