- 1Postgraduate Program in Parasitic Biology in The Amazon (PPGBPA), State University of Pará (UEPA), Belém, Brazil
- 2Nursing Department, University of the Amazon (UNAMA), Ananindeua, Brazil
- 3Postgraduate Program in Epidemiology and Health Surveillance (PPGEVS), Evandro Chagas Institute (IEC), Ananindeua, Brazil
- 4Bacteriology and Mycology Section, Evandro Chagas Institute (IEC), Ananindeua, Brazil
Coronavirus disease 2019 (COVID-19) is a multisystemic disease caused by SARS-CoV-2 that can lead to several pulmonary illnesses according to the immunological contexts of the individual. Haplotypes consist of single-nucleotide polymorphisms (SNPs) within candidate genes for diseases. TLR2 and TLR3 are genes located on human chromosome 4 (chr:4) and composite a haplotype that influence immune signaling and inflammatory pathways. The purpose of this article was to genetically analyze in silico a cohort of professionals from Belém-PA during the first wave of the pandemic using SNPs rs3804100, rs3775290, and rs3775291 on the human chr:4. This is a computational genomic design using bioinformatic software and machine-learning technologies on epidemiological data of Sanger sequencing data. Regarding the findings, none of the alleles formed by the haplotype showed statistical significance for symptomatology or disease severity. The haplotype block was not significant between the SNPs analyzed despite a high permutation rate of alleles at the beginning of the variance of the individual genomic data. Then, the TLR2–TLR3 haplotype (SNPs rs3804100, rs3775290, and rs3775291) showed little determination in the clinic of individuals with COVID-19 in Belém (Northern Brazil), which may indicate differences in collective genetic patterns and/or epigenetic influences compared to other more affected populations that have the same haplotype pattern.
1 Introduction
Coronavirus disease 2019 (COVID-19) was declared a pandemic by the World Health Organization (WHO) in March 2020 due to its high spread. The epidemic came on abruptly and left the Public Health System under severe pressure and not knowing how to manage the cases that arose, inducing infected people into a critical state and consequently many deaths (Chen et al., 2020; Choudhury and Mukherjee, 2020). As it is a highly transmissible pathogen, it was recommended that people should be restricted from circulating to reduce contamination rates (Bastos et al., 2021). However, health professionals did exactly the opposite of what was recommended for most of the population, which made them more susceptible to acquiring the disease (Coelho et al., 2022).
There are different clinical manifestations, depending on the host’s immune status, the genetic variability of the virus circulating at the time, and other genetic and even sociodemographic aspects; in this sense, the symptomatology of COVID-19 ranges from asymptomatic or with mild manifestations similar to other severe acute respiratory syndromes (SARS)—even acute respiratory distress, sepsis, septic shock, acute thrombosis, and/or multiple organ dysfunction. These complications often require hospitalizations for several days in ICUs, with many patients not surviving, and those who recovered would have to live with the sequelae left behind (post-COVID-19) (Cham et al., 2021; Fronteira et al., 2024).
The COVID-19 pandemic has affected the world in a big way, with the first pandemic wave characterized by rapid spread and high mortality. According to studies, from 1 January 2020 to 31 December 2021, there were approximately 5.94 million deaths worldwide. In Brazil, the first wave hit the country hard, showing the weaknesses of the health system and placing health professionals on the front line. Professionals faced long working hours, a lack of personal and collective protective equipment (PPE and CPE), and the anguish of dealing with a high number of critically ill patients. The emotional and physical impacts were significant, with many professionals suffering from anxiety and other mental health problems (Agrawal et al., 2021). According to the Brazilian Institute of Geography and Statistics (IBGE), approximately 1.5 million deaths occurred in 2020 in Brazil because of COVID-19. The front line of the pandemic in Brazil was marked by the dedication of health professionals, who fought to save lives in the midst of a troubled scenario (Almeida et al., 2020; COVID-19 Excess Mortality Collaborators, 2022).
Immunology explores the complex defense system of the body against pathogens. Innate immunity, the first line of defense, uses Toll-like receptors (TLRs) to recognize molecular patterns associated with pathogens (PAMPs). In particular, TLR2 detects lipopeptides and bacterial peptidoglycans, whereas TLR3 recognizes viral double-stranded RNA. Both TLRs play a crucial role in the detection of viruses, including SARS-CoV-2. Studies indicate that SARS-CoV-2 can be recognized using TLRs, including TLR2 and TLR3, through different mechanisms, with the interaction of viral proteins with these receptors, leading to the activation of the innate immune response, thus preventing excessive inflammatory response (Dhangadamajhi and Rout, 2021; Choudhary et al., 2022; Machado et al., 2023).
Haplotypes are combinations of alleles at different markers along the same chromosome that are inherited as a unit. In research, their main purpose is to identify the causes of diseases, marking recombination events in family and population-based studies, and clinically, they apply to the greater or lesser correlation in predicting the severity of diseases. Single-nucleotide polymorphisms (SNPs) can originate from point mutations in DNA, which are intragenic and cause important biological effects such as association with complex diseases and reactions in responses to treatments in humans (Bafica et al., 2005; Chen et al., 2017; Bortolotti et al., 2021).
This study of the TLR2–TLR3 haplotype, which includes two SNPs (rs3804100 in TLR2 and rs3775290 and rs3775291 in TLR3), may aid in the identification of genetic differences that influence COVID-19 susceptibility and clinical progression, hence facilitating the development of individualized disease containment and treatment methods. These SNPs (rs3804100 in TLR2 and rs3775290 and rs3775291 in TLR3) were selected because they are in nearby functional regions, previously associated in the literature with viral immune response and respiratory diseases (such as tuberculosis, influenza, and other viral infections) (Bucciol et al., 2022; Córdova-Dávalos et al., 2022; Chen et al., 2023; Silva et al., 2023a). Thus, investigating the haplotypes of SNP rs3804100 in the TLR2 gene and SNPs rs3775290 and rs3775291 in the TLR3 gene is a viable approach to understanding the immunogenetic mechanisms driving COVID-19 and its clinical symptoms (Chen et al., 2011; Vega et al., 2017; Sharif et al., 2020; Djordjevic et al., 2022).
Considering the impact of the pandemic on the workforce of healthcare professionals and the need to identify the association of haplotypes for the prognosis of COVID-19 in these professionals (Zangoue et al., 2021), along with the relevance of immunogenetics in this context, in this study, we aimed to analyze the haplotypes of the TLR2 and TLR3 genes on the human chromosome, in the genomic area near SNP rs3804100 of TLR2 and SNPs rs3775290 and rs3775291 of TLR3 in the labor cohort of Belém-PA, Brazil, from the first pandemic wave.
2 Materials and methods
2.1 Study design and ethical considerations
This study is characterized as a case–control type, with convenience sampling, and it is part of a larger scientific project. It followed the recommendations of the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) (Von Elm et al., 2014). This study protocol was approved by local ethics committees, and all subjects gave their written informed consent (Term of Free and Informed Consent—TCLE). This work was approved by the Research Ethics Committee of the State University of Pará—UEPA (CAAE: 38113620.1.0000.5174), with opinion number: 6.124.862. This research was carried out in accordance with the Helsinki Declaration (Helsinki Declaration, 1964) and Resolution No. 466/2012 of the Brazilian National Health Council (Brasil et al., 2012).
2.2 Settings and participants
This study was carried out in 10 health centers and institutions located in the city of Belém-PA, Brazil (in the Amazon Region—Northern Brazil). During the initial wave of the pandemic, individuals infected with COVID-19 were admitted to all collaborating institutions in Brazil (López-Juárez et al., 2021). Based on the acquisition of contacts that were appropriate for this research group, these institutions were selected at random. Epidemiological analyses were carried out using a questionnaire applied in person at health centers and analyses of genetic data sequenced from blood biological samples using molecular biology methods.
Professionals (in health, administrative, and general service areas) who actively worked during the first wave of the pandemic in the period between 1 April 2021 and 30 June 2020 (Bastos et al., 2021), who were directly and daily exposed to and infected by the virus in Belém, and who agreed to participate in the research by signing the consent form were included in this study. This cohort of professionals (convenience sampling methodology, institutions, and health centers selected) has already been portrayed by our research group in earlier studies involving other epidemiological and genetic analyses (Silva et al., 2023c; Silva et al., 2023b; Marcelino et al., 2025). All employees at healthcare facilities during this time may have been exposed to SARS-CoV-2. Because safety regulations were still developing in Brazil, medical facilities were overcrowded, masks were scarce, and N95 masks were not worn (Quigley et al., 2021).
2.3 Division of groups, variables, participants, sampling, and laboratory procedures
The study population questionnaires were evaluated to collect epidemiological data and potential association with genetic information of the participants, which were evidenced in demographic variables, comorbidities, signs and symptoms, professional category, and records of illnesses of intra- and extra-household contacts of the consanguineous and non-consanguineous types. Data included signs and symptoms associated with the first wave period, sex, age, employment, and comorbidities. Based on the patient’s self-report during the questionnaire, chest computed tomographic (CCT) data were included.
Sample collection was carried out between 1 June 2021 and 30 March 2022. Laboratory procedures were conducted in the Molecular Biology Laboratory—LABIMOL, Bacteriology and Mycology Section (SABMI) of the Evandro Chagas Institute (IEC). DNA extractions were performed using the DNeasy Blood and Tissue Kit (QIAGEN®, Venlo, Netherlands), following the manufacturer’s instructions. For all these professionals, SNP rs3804100 (T>C) of TLR2 and SNPs rs3775291 (C/T) and rs3775290 (C/T) of TLR3 were evaluated to correlate with the susceptibility and severity of COVID-19. The information on the SNPs (SNP ID) of TLR2 and TLR3 was retrieved from the National Center for Biotechnology Information—NCBI, dbSNP (http://www.ncbi.nlm.nih.gov/snp/; accessed on 1 January 2022) (Sayers et al., 2023). Between 1 November 2021 and 1 November 2022, the sample collections and questionnaires were performed. Polymerase chain reaction (PCR) was performed with Platinum Taq DNA Polymerase, DNA-free (Invitrogen®, Thermo Fisher Scientific Corporation, Waltham, Massachusetts, United States) using the thermocycler. The amplified products of PCR were subjected to electrophoresis in a 2% agarose gel containing 3.0 µL of Sybr Safe (Invitrogen®, Thermo Fisher Scientific Corporation, Waltham, MA, United States) to visualize the amplified DNA fragments in a photodocumenting device. TLR2 SNP rs3804100 and TLR3 SNPs rs3775291 and rs3775290 were selected together for the analysis of the present article due to their proximity in the genomic region of the same human chromosome.
The EasyPure PCR Purification Kit (TransGen Biotech Co.®, Beijing, China) was used to purify the PCR products in accordance with the manufacturer’s instructions. The BigDye X-Terminator Kit instructions were followed to submit the previously purified amplified products to the ABI 3130 Genetic Analyzer sequencer (Applied Biosystems®, Life Technologies, Thermo Fisher Scientific Corporation, Waltham, MA, United States) for the sequencing reaction. The Bioedit program version 7.2.5 was used to visualize and analyze the SNP areas of interest, and BLAST was then performed on the NCBI website (https://blast.ncbi.nlm.nih.gov/Blast.cgi/), obtained on 15 August 2023 (Hall, 1999).
The criteria and forms of grouping based on symptoms and severity of COVID-19 [AS, asymptomatic group—n = 91; SI, symptomatic group—n = 121 (divided into SCP, symptomatic with pulmonary involvement, with n = 34, and SSP, symptomatic without pulmonary involvement, with n = 5, based on the participants’ CCT data)], criteria for tomographic data evaluation of severity of COVID-19 through CCT (≥10% in CCT data, according to the guidelines provided by the WHO and the Ministry of Health) (Choi et al., 2020; UNA-SUS, 2020), sample size (N = 212) and power of sampling (p = 0.94), variables analyzed using the questionnaires, method of defining COVID-19 cases during the first wave in Brazil (based on Brazilian Ministry of Health recommendations), inclusion and exclusion criteria for participants, epidemiological data from the cohort, blood collection method, method of extracting DNA from samples, laboratory procedures, primers used, PCR and Sanger sequencing methods, and statistical investigation of agreement of each SNP with the Hardy–Weinberg equilibrium (HWE, p > 0.05) and comparative and associative analyses of genotyping results in a genotypic and allelic form of SNPs rs3804100, rs3775291, and rs3775290 with COVID-19 in this population evaluated (individually for each SNP) have already been reported in our previous studies, however, without the haplotypic investigation of this present work (Silva et al., 2023c; Silva et al., 2023b; Marcelino et al., 2025).
2.4 Bioinformatic and statistical analyses
Only subjects who had all SNPs correctly genotyped were included in the haplotype analysis, which was constructed from genotype data using the computer tool Haploview v4.2. At p < 0.05, statistical significance was established. Furthermore, the Haploview 4.2 tool was used to perform linkage disequilibrium (LD) blocks and haplotype association risk studies (Barrett et al., 2005). Thesias v3.1 software, which is based on the stochastic expectation–maximization technique for haplotype analysis, was used to estimate haplotype frequencies, HWE concordance, and their impacts. G*Power v3.1.9.7 software was used to measure sample size power using a Chi-square goodness-of-fit test (Kang, 2021). The power of the sample size, the normality of the variables, the Hardy–Weinberg equilibrium, and the epidemiological characteristics of the individuals in this cohort of professionals from Belém have already been established and calculated by previous studies carried out by our research group (Silva et al., 2023c; Silva et al., 2023b).
The level of connection between the examined haplotypes and COVID-19 (symptomatology, As and SI; and severity, AS, SCP, and SSP) was measured using the odds ratio (OR) and its 95% confidence interval (CI). The impact of each haplotype on COVID-19 is shown using the Thesias program as predicted phenotypic averages for a single dosage of each haplotype, together with the associated 95% CI in relation to the reference haplotype. The most prevalent haplotype in the study groups is the reference haplotype, which was determined using the Thesias program (Djordjevic et al., 2022). The statistical software program SPSS v26 (IBM Corp.®, Armonk, NY, United States) was used for all statistical analyses. The 95% CIs and ORs show the relative risk.
To contextualize the allele frequencies found in our local cohort within the spectrum of global genetic diversity, an in silico comparative analysis was performed using large-scale public data. The allele frequencies of SNPs rs3804100, rs3775290, and rs3775291, obtained through Sanger sequencing of our sample, were compared with the frequencies reported for the same variants in the Genome Aggregation Database (gnomAD), version 4.0.0. In particular, data from our Belém population were contrasted with the frequencies of the global aggregate and the “Latino/Admixed American” superpopulation available in gnomAD. The gnomAD database was accessed on 23 August 2025, and frequencies were extracted directly from the specific entries for each SNP in the portal. This analysis allowed us to assess how representative or divergent the genetic frequencies of our Amazonian cohort are in relation to broader reference populations (Gudmundsson et al., 2022).
3 Results
3.1 COVID-19 symptomatology parameter and haplotype analyses for this cohort
The chi-square goodness-of-fit test was used to calculate sample size power for the cohort population of symptomatic (N = 121) and asymptomatic (N = 91) COVID-19 individuals, with an error probability α of 0.05 and an effect size of 0.3. The actual real power (1-β error probability) is 0.94, which is higher than 0.80, and this is statistically acceptable. The results of the data sequenced and then analyzed revealed a pattern of allele frequency and haplotypic background without statistical significance in the influence of the investigated phenotype of COVID-19 symptomatology (Tables 1, 2).

Table 1. Allelic distribution and associative comparison of frequencies in the haplotypes for these groups.

Table 2. Distributive and association analyses of each haplotype in the evaluated groups to COVID-19, regarding the reference haplotype.
To further examine the association of the candidate SNPs between the AS × SI groups, we estimated the LD and haplotype blocks; SNPs rs3775290 and rs3804100, SNPs rs3804100 and rs3775291, and SNPs rs3775290 and rs3775291 were in low LD with each other (r2 = 0.08, r2 < 0.08, and r2 < 0.08), showing low r2 values, indicating low co-inheritance between the variants in the studied group and their limited impact on the haplotype (Figure 1).

Figure 1. Linkage disequilibrium (LD) blocks and R2 values between the analyzed polymorphisms in the TLR2 (rs3804100) and TLR3 (rs3775290 and rs3775291) genes for the total cohort sample (N = 212).
Permutation tests were performed using Haploview software, running 1,000 iterations for individual SNPs to analyze the possible change in the phenotype of the symptoms related to the changes in the alleles (Table 3), and this statistical profile of permutations is shown in Figure 2.

Table 3. Results of significance tests for 1,000 permutations of SNPs regarding the change in the symptomatology phenotype of COVID-19 in this cohort (N = 212).

Figure 2. Permutation test results of SNPs rs3804100 (TLR2), and rs3775290 and rs3775291 (TLR3) showing Chi-square distributions for COVID-19 symptomatology in the Belém cohort.
3.2 COVID-19 severity parameter and haplotype analyses for this cohort
In silico analysis for the paired evaluation of the SNPs of this haplotype was carried out on this cohort of professionals evaluated in relation to the severity of the disease (Table 4).

Table 4. Comparative analysis of the association of each haplotype with the profiles of COVID-19 severity groups (regarding reference haplotype: TCC).
3.3 Comparative in silico analysis of SNP allele frequencies
To contextualize our findings within the global genetic diversity, we compared the minor allele frequencies (MAFs) of the studied SNPs in our cohort (Belém) with those from the Genome Aggregation Database (gnomAD v4.0.0) (Table 5). As shown in Table 5, the allele frequencies in our population markedly differ from global aggregates, particularly for rs3804100 (MAFBelém: 4.7% vs MAFgnomAD Global: 24.5%). Notably, our cohort’s frequencies are more similar, yet still they are distinct from those of the gnomAD Latino/admixed American superpopulation (e.g., rs3804100 MAF: 4.7% vs 10.2%).

Table 5. Comparative minor allele frequency (MAF) analysis of SNPs rs3804100, rs3775291, and rs3775290 between the Belém cohort (this study) and the gnomAD v4.0.0 database.
3.4 Main findings of the article
The results demonstrated that the analyzed SNPs (rs3804100 in TLR2 and rs3775290 and rs3775291 in TLR3) presented specific allele frequencies in the studied cohort, with the most frequent minority variant being the T allele of rs3775290 (MAF = 25%) and the variants rs3775291 and rs3804100 presenting lower frequencies (MAF <5%). Analysis of the haplotypes derived from these SNPs showed a predominance of the TCC reference haplotype (∼70%), with other haplotypes distributed at lower frequencies.
No statistically significant association was found between these haplotypes/alleles and the presence or absence of COVID-19 symptoms (p > 0.05), either between the asymptomatic and symptomatic groups or between the different degrees of severity assessed through CCT (with lung involvement). The permutation test confirmed the absence of a significant association between the variants and the phenotypic variability of symptoms.
The results of the in silico comparative analysis of allele frequencies reveal a distinct genetic profile in the Belém cohort compared to reference population databases. The frequency of the minority C allele of the TLR2 rs3804100 SNP in our sample (4.7%) was drastically lower than the frequency reported for the global population in gnomAD (24.5%). More importantly, even when compared to the most genetically similar group—the Latino/admixed American cluster (10.2%)—the frequency in Belém remains less than half. This pattern is repeated, less pronounced, but still evident, for the TLR3 rs3775290 SNP, whose T allele has a frequency of 25.0% in Belém, compared to 28.4% in Latinos and 32.7% globally. SNP rs3775291, already rare globally (1.1%), was even rarer in the Belém sample (0.7%).
4 Discussion
COVID-19 continues to be a public health problem with numerous endemic foci around the world, and variants of concern (VOCs) continue to emerge in the face of parasite–host interactions and high morbidity in the elderly and people with comorbidities (Cui et al., 2025). Immunogenetic studies with haplotypes, such as this study, are important for revealing genetic patterns of ethnicities (populations), along with epigenetic and/or environmental aspects, and associating them with phenotypes/clinical–epidemiological conditions such as COVID-19 (Contes and Liu, 2025). Candidate genes, especially those encoding immune receptors that are essential in host defense, are preponderant in determining associative genetic biomarkers. In this sense, TLR genes (such as TLR2 and TLR3) are potential typical elements for consideration in immunogenetic analyses (Bochud et al., 2007; CHEN et al., 2021).
Previous studies have demonstrated associations between TLR2 and TLR3 polymorphisms and COVID-19 severity and susceptibility, underscoring the importance of these receptors in viral recognition and immune activation (Alseoudy et al., 2022; Alhabibi et al., 2023). Considering that both tuberculosis and COVID-19 affect the respiratory system and that TLR2 is crucial in IFN-γ-mediated cellular immunity, the analysis of SNP rs3804100 is relevant to understanding the heterogeneity of the host’s immune response to SARS-CoV-2 infection in terms of more frequent symptoms typical of lung involvement (Bafica et al., 2005).
Our work extends these insights by focusing on haplotype analyses within a Brazilian cohort, providing valuable data relevant to the genetic diversity and epidemiology of COVID-19 in this Amazon region (Secolin et al., 2021). The investigation of the TLR2–TLR3 haplotype that includes these three SNPs may help identify other genetic variations that influence susceptibility and clinical progression of COVID-19, aiding in the development of personalized containment and treatment strategies for the disease. Thus, the study of the haplotype of SNP rs3804100 in the TLR2 gene and SNP rs3775291 in the TLR3 gene represents a promising approach to unravel the immunogenetic mechanisms underlying COVID-19 and its clinical manifestations.
Given TLR2’s critical role in the innate immune response to viral infections, especially SARS-CoV-2, the haplotype including SNP rs3804100 in the TLR2 gene is a key topic of investigation in COVID-19 immunogenetics. The TLR2 gene activates the release of pro-inflammatory cytokines, including TNF-α and interferon-gamma, to defend against respiratory infections (Bafica et al., 2005; Bochud et al., 2007; Khanmohammadi and Rezaei, 2021).
SNP rs3804100, located in the exon 3 of the TLR2 gene, characterized by the substitution of thymine (T) with cytosine (C), may alter the stability of the TLR2 protein; however, studies conducted among healthcare professionals during the first wave of the pandemic indicated that this variant is benign and showed no significant association with the severity or symptoms of COVID-19 in that population (Cheng et al., 2015; Silva et al., 2023b). In a systematic review study by Schurz et al. (2015), more severe tuberculosis cases were associated in all genetic predictive models, with a higher risk of susceptibility due to the presence of the mutant allele (C) of the TLR2 rs3804100 SNP in Asian populations (Schurz et al., 2015).
TLR3 recognizes double-stranded viral RNA (dsRNA), activating signaling pathways that lead to the production of interferons and other antiviral cytokines, which are essential for controlling viral replication and modulating the inflammatory response (Bortolotti et al., 2021). SNP rs3775290, located in the exon 4 of the TLR3 gene, is a genetic variant that has been studied for its influence on susceptibility to COVID-19 due to the crucial role of TLR3 in the innate immune response against viruses, including SARS-CoV-2. TLR3 recognizes viral RNA and activates signaling pathways that promote the production of interferons and other antiviral cytokines, which are essential for controlling the infection. Studies have investigated the association of rs3775290 with the severity and susceptibility to COVID-19, seeking to understand how this variant may alter the efficiency of the host’s antiviral response. In particular, Alseoudy et al. (2022) found that the T allele of rs3775290 was associated with a higher risk of severe COVID-19 pneumonia (Alseoudy et al., 2022; Silva et al., 2023c).
The rs3775291 variant (located on exon 4 of the TLR3 gene) is a nucleotide substitution (C>T) that can cause structural and functional changes in the TLR3 protein, potentially impacting the effectiveness of the antiviral immune response. It results in a proline-to-serine change (P554S) and has also been associated with diminished TLR3 signaling and IFN-β production (Bucciol et al., 2022). Studies have investigated the association of this variation with COVID-19 susceptibility and severity; however, the results are inconsistent. Some studies found no significant association between SNP rs3775291 and COVID-19 severity or symptoms in specific populations (Dhangadamajhi and Rout, 2021; Silva et al., 2023c). However, research by Silva et al. (2023b) indicated that this variant may be associated with an increased risk of viral infections in Asian and American populations, suggesting a potential high-risk genetic pattern for COVID-19 in these groups (Xie and Zhu, 2020; Dhangadamajhi and Rout, 2021; Silva et al., 2023c; Silva et al., 2023c).
SNP rs3804100 was chosen because it (T > C) has already been widely studied in other respiratory diseases, such as tuberculosis, which, like COVID-19, affects the lungs (Flores-González et al., 2024). This suggests that variations in this SNP may influence the immune response in lung infections, causing instability in the formation of the TLR2 protein. More directly, Flores-González et al. (2024) demonstrated that the rs3804100 variant induces a distinct profile of TLR2 expression and cytokine release (e.g., IL-6 and TNF-α) in response to the SARS-CoV-2 spike protein, providing a direct functional link between this SNP and the COVID-19 immune response (Flores-González et al., 2024). The rs3775291 SNP has already been associated with susceptibility to various viral infections, and variations in this SNP in a non-synonymous way, that is, that alter the structure of the TLR3 protein, can impact its immune function and, consequently, the host’s response to COVID-19 infection (Andrade et al., 2024). In the case of SNP rs3775290, this biomarker was chosen for analysis in this study because of its location in the TLR3 gene, which encodes a key receptor for detecting RNA viruses and activating the innate immune response (Wang et al., 2015).
The TLR2 and TLR3 genes encode receptors essential for recognizing pathogen molecular patterns and activating essential innate immune responses. The rs3804100 variant in TLR2 has been linked to alterations in protein stability and modulation of inflammatory activation, whereas the rs3775290 and rs3775291 variants in TLR3 may influence viral RNA detection and subsequent interferon production (Dos Santos et al., 2025). The lack of a significant association suggests that in this group of professionals exposed and infected during the first wave of the pandemic in the Amazon region, these TLR2 and TLR3 variants do not play a determining role in modulating susceptibility or clinical severity of COVID-19. This may indicate that other genetic, epigenetic, environmental, or immunological factors have a greater influence on the host response, or that the genetic heterogeneity of the Amazonian population confers different risk patterns (Evdokimov et al., 2021; Minashkin et al., 2022).
The low LD observed between the SNPs evaluated (r2 < 0.08) indicates that these variants segregate virtually independently, which may explain the absence of strong haplotypic blocks in the study population. This genetic independence confirms that the combined effect of SNPs on the same chromosome 4 is limited, reinforcing the lack of a clearly defined risk haplotype associated with COVID-19 in the group evaluated (Salamaikina et al., 2023; Teräsjärvi et al., 2024).
From a clinical perspective, the lack of association between the tested SNPs and the severity or symptoms of COVID-19 in healthcare workers may reflect local genetic and environmental particularities. Previous studies have indicated that the immune response to viral infections is multifactorial, involving multiple genes, cellular pathways, and external influences beyond SNPs (Feng et al., 2022).
The comparison of the genetic diversity in our present study provides crucial context, showing that the allele frequencies in our Brazilian Amazonian cohort are significantly different from those observed in European populations but are more closely aligned, although still distinct, with those of other Latino/admixed American populations. This genetic distinctiveness strongly supports our conclusion that population-specific genetic architecture is a key factor in explaining the lack of association found in our study, which might be present in other, more well-studied populations. In addition, this comparative analysis underscores the unique genetic background of the Amazonian population studied and reinforces that the lack of association observed may be attributable to population-specific genetic patterns, limiting the generalizability of findings from other ethnic groups (Rodrigues de Moura et al., 2015; Silva et al., 2023c; Silva et al., 2023b).
This disparity has two critical implications, which provide a possible fundamental explanation for this study’s negative findings. First, it validates the unique genetic background of the Amazonian population studied, which is not adequately represented even in large databases that include Latinos (Lins et al., 2010). Second, and more crucial for the interpretation of our null results, it demonstrates that potentially risky or protective alleles, previously associated with COVID-19 in European and Asian populations where they are more common, are simply very rare or absent in this specific cohort. Consequently, the statistical power to detect any significant effect of these specific variants was inevitably low. Therefore, the lack of association does not necessarily invalidate findings from other studies but rather highlights that the genetic determinants of COVID-19 are highly dependent on the population context. Our negative findings are, in fact, a positive finding that highlights the importance of conducting genetic studies in diverse and underrepresented populations, as risk factors are not universal (Pena et al., 2020).
Despite the lack of statistical association, the present study contributes to mapping genetic variations and limitations of the role of these SNPs in the analyzed population. The findings reinforce the complexity of the COVID-19 genetic profile and the importance of studying diverse populations, especially those in the Amazon, whose genetic variability and environmental factors are poorly represented in large genomic databases. The lack of association between the analyzed SNPs highlights the need to address a broader panel of immunological variants and molecular mechanisms to understand how clinical response is modulated in COVID-19. These data are useful for directing future research, indicating that studies with larger sample sizes and integrative analyses involving genomics, transcriptomics, and function are necessary to elucidate relevant immunogenetic mechanisms in the Amazon population (Secolin et al., 2021).
This present study has some important limitations that should be considered when interpreting its results. A central limitation is the relatively small sample size (N = 212) and the use of a convenience sampling methodology, focused on healthcare professionals from a specific region (Belém, Pará, in the Brazilian Amazon Region). Although the calculated sample power (p = 0.94) is adequate for more common variants, this limited sample substantially reduces the statistical power to detect significant associations with low allele frequency variants (MAF <5%), such as SNPs rs3775291 and rs3804100. Consequently, the results for these rarer polymorphisms should be interpreted with extreme caution as the absence of a statistical association may reflect an inability to detect a true effect (type II error) rather than definitive evidence of its absence.
Another important consideration relates to the definition of COVID-19 cases during the early stages of the pandemic in Brazil. At that time, the criteria established by the Ministry of Health were still evolving and relied predominantly on clinical–epidemiological suspicion, with limited access to confirmatory laboratory testing (Hakim et al., 2021; Rodrigues et al., 2024). This context increased the probability of misclassification bias, with some true cases being missed (false negatives) and other respiratory infections being incorrectly categorized as COVID-19 (false positives). Additionally, the evolving understanding of COVID-19 symptomatology—ranging from asymptomatic to atypical systemic manifestations—further complicated the consistency of case definition. The lack of detailed longitudinal data also limits the analysis of the genetic influence on post-COVID-19 syndrome. Therefore, the findings should be considered hypothesis-generating for this unique population, highlighting the need for future studies with larger sample sizes, more comprehensive genetic panels, and a design that allows for more robust control of environmental and epigenetic variables (Sharpe, 2015).
The understanding of individual variations and associated genetic patterns in genomic blocks can aid in characterizing the immune response to SARS-CoV-2 and may, in the future, help in developing personalized treatment and prevention strategies for the disease (Caratachea, 2007; Aguiar et al., 2021). However, further investigations are necessary to clarify their exact role in the pathogenesis and clinical evolution of COVID-19.
5 Conclusion
Therefore, the TLR2–TLR3 haplotype (SNPs rs3804100, rs3775290, and rs3775291) showed little determination in the clinic of individuals with COVID-19 in Belém (Northern Brazil), which may indicate differences in collective genetic patterns and/or epigenetic influences compared to those in other more affected populations that have the same haplotype pattern worldwide. Moreover, the results obtained in this epidemiological study provide important insights into the Brazilian Amazon population genomics about COVID-19 involvement; however, the complexity of the factors involved and the identified methodological limitations indicate the need for additional studies to strengthen the prediction and confirmation of the data presented here.
Then, it is recommended to conduct complementary research, preferably with more robust analytical designs, greater control of confounding variables, and representative samples, to validate the raised hypotheses and improve the understanding of the causal mechanisms involved. In this way, it will be possible to increase the precision of the estimates and the validity of the inferences, contributing to the development of more effective preventive strategies and interventions in the studied population.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by the Research Ethics Committee—CEP of the State University of Pará (UEPA), Brazil, under opinion number 6.124.862. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
MS: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Resources, Software, Supervision, Validation, Visualization, Writing – original draft. LF: Conceptualization, Data curation, Formal Analysis, Investigation, Resources, Visualization, Writing – original draft. DS: Data curation, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft. EC: Data curation, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft. NA: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Writing – original draft. SB: Data curation, Formal Analysis, Investigation, Methodology, Resources, Writing – original draft. TA: Data curation, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft. ED: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Validation, Writing – original draft. AG: Formal Analysis, Investigation, Methodology, Resources, Validation, Visualization, Writing – original draft, Writing – review and editing. LL: Conceptualization, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing.
Funding
The author(s) declare that no financial support was received for the research and/or publication of this article.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Agrawal, V., Yadav, S. K., Agarwal, P., and Sharma, D. (2021). Strategies for optimizing the use of PPE during surgery in COVID-19 pandemic: rapid scoping review of guidelines. Indian J. Surg. 83, 17–27. doi:10.1007/s12262-020-02713-x
Aguiar, V. R. C., Augusto, D. G., Castelli, E. C., Hollenbach, J. A., Meyer, D., Nunes, K., et al. (2021). An immunogenetic view of COVID-19. Genet. Mol. Biol. 44, e20210036. doi:10.1590/1678-4685-gmb-2021-0036
Alhabibi, A. M., Hassan, A. S., Abd Elbaky, N. M., Eid, H. A., Khalifa, M. A. A. A., Wahab, M. A., et al. (2023). Impact of toll-like receptor 2 and 9 gene Polymorphisms on COVID-19: Susceptibility, Severity, and Thrombosis. J. Inflamm. Res. 16, 665–675. doi:10.2147/JIR.S394927
Almeida, C. A. R. P. N. D., Almeida, G. A. R. P. N. D., Carvalho, M. R. C. T. D., and Marcolino, A. B. D. L. (2020). Aspectos relacionados à saúde mental dos profissionais de saúde durante a pandemia do Covid-19: uma revisão integrativa da literatura/Mental health aspects of health professionals during the Covid-19 pandemic: an integrative literature review. BJHR 3, 19481–19491. doi:10.34119/bjhrv3n6-327
Alseoudy, M. M., Elgamal, M., Abdelghany, D. A., Borg, A. M., El-Mesery, A., Elzeiny, D., et al. (2022). Prognostic impact of toll-like receptors gene polymorphism on outcome of COVID-19 pneumonia: a case-control study. Clin. Immunol. 235, 108929. doi:10.1016/j.clim.2022.108929
Andrade, A. S., Bentes, A. A., Diniz, L. M., Carvalho, S. H., Kroon, E. G., and Campos, M. A. (2024). Association between Single-Nucleotide Polymorphisms in toll-like receptor 3 (tlr3), tlr7, tlr8 and tirap genes with severe symptoms in Children Presenting COVID-19. Viruses 17, 35. doi:10.3390/v17010035
Bafica, A., Scanga, C. A., Feng, C. G., Leifer, C., Cheever, A., and Sher, A. (2005). TLR9 regulates Th1 responses and cooperates with TLR2 in mediating optimal resistance to Mycobacterium tuberculosis. J. Exp. Med. 202, 1715–1724. doi:10.1084/jem.20051782
Barrett, J. C., Fry, B., Maller, J., and Daly, M. J. (2005). Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263–265. doi:10.1093/bioinformatics/bth457
Bastos, L. S., Ranzani, O. T., Souza, T. M. L., Hamacher, S., and Bozza, F. A. (2021). COVID-19 hospital admissions: brazil’s first and second waves compared. Lancet Respir. Med. 9, e82–e83. doi:10.1016/S2213-2600(21)00287-3
Bochud, P.-Y., Magaret, A. S., Koelle, D. M., Aderem, A., and Wald, A. (2007). Polymorphisms in TLR2 are associated with increased viral shedding and lesional rate in patients with genital Herpes Simplex virus type 2 infection. J. Infect. Dis. 196, 505–509. doi:10.1086/519693
Bortolotti, D., Gentili, V., Rizzo, S., Schiuma, G., Beltrami, S., Strazzabosco, G., et al. (2021). TLR3 and TLR7 RNA sensor activation during SARS-CoV-2 infection. Microorganisms 9, 1820. doi:10.3390/microorganisms9091820
Brasil, S., da, M., and Saúde, C. N. de (2012). Resolução n. 466/2012, de 12 de dezembro de 2012. Dispõe sobre diretrizes e normas regulamentadoras sobre pesquisas envolvendo seres humanos., (CNS Brasília).
Bucciol, G., Desmet, L., Moens, L., Delafontaine, S., Corveleyn, A., Meyts, I., et al. (2022). Pathogenic P554S variant in TLR3 in a patient with severe Influenza pneumonia. J. Clin. Immunol. 42, 430–432. doi:10.1007/s10875-021-01172-7
Caratachea, M. A. C. (2007). Polimorfismos genéticos: importancia y aplicaciones. Rev. Inst. Nal Enf. Resp. Mex. 20, 213–221.
Cham, L. B., Pahus, M. H., Grønhøj, K., Olesen, R., Ngo, H., Monrad, I., et al. (2021). Effect of Age on innate and adaptive immunity in hospitalized COVID-19 patients. J. Clin. Med. 10, 4798. doi:10.3390/jcm10204798
Chen, K., Gu, W., Zeng, L., Jiang, D., Zhang, L., Zhou, J., et al. (2011). Identification of haplotype tag SNPs within the entire TLR2 gene and their clinical relevance in patients with major trauma. Shock 35, 35–41. doi:10.1097/SHK.0b013e3181eb45b3
Chen, D., Xie, W., Lu, Y., Su, S., Nong, L., Jia, Y., et al. (2017). Gene polymorphisms of TLR2 and TLR3 in HBV clearance and HBV-related hepatocellular carcinoma in a Chinese male population. Int. J. Biol. Markers 32, e195–e201. doi:10.5301/jbm.5000238
Chen, N., Zhou, M., Dong, X., Qu, J., Gong, F., Han, Y., et al. (2020). Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. Lancet 395, 507–513. doi:10.1016/S0140-6736(20)30211-7
Chen, Y., Lin, J., Zhao, Y., Ma, X., and Yi, H. (2021). Toll-like receptor 3 (TLR3) regulation mechanisms and roles in antiviral innate immune responses. J. Zhejiang Univ. Sci. B 22, 609–632. doi:10.1631/jzus.B2000808
Chen, R., Wang, X., Li, Z., Dai, Y., Du, W., and Wu, L. (2023). Human Toll-like receptor 2 genetic polymorphisms with tuberculosis susceptibility: a systematic review and meta-analysis. Cytokine 172, 156405. doi:10.1016/j.cyto.2023.156405
Cheng, Y., Zhu, Y., Huang, X., Zhang, W., Han, Z., and Liu, S. (2015). Association between TLR2 and TLR4 gene polymorphisms and the susceptibility to inflammatory bowel disease: a meta-analysis. PLoS One 10, e0126803. doi:10.1371/journal.pone.0126803
Choi, H., Qi, X., Yoon, S. H., Park, S. J., Lee, K. H., Kim, J. Y., et al. (2020). Extension of coronavirus disease 2019 on chest CT and implications for chest radiographic interpretation. Radiol. Cardiothorac. Imaging 2, e200107. doi:10.1148/ryct.2020200107
Choudhary, M. L., Chaudhary, U., Salve, M., Shinde, P., Padbidri, V., Sangle, S. A., et al. (2022). Functional single-nucleotide polymorphisms in the MBL2 and TLR3 genes influence disease severity in influenza A (H1N1)pdm09 virus-infected patients from Maharashtra, India. Viral Immunol. 35, 303–309. doi:10.1089/vim.2021.0179
Choudhury, A., and Mukherjee, S. (2020). 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. Virology 92, 2105–2113. doi:10.1002/jmv.25987
Coelho, M., de, M. F., Cavalcante, V. M. V., Cabral, R. L., Oliveira, R. M., Nogueira, P. S. F., et al. (2022). Contexto de trabalho e manifestações clínicas da COVID-19 em profissionais de saúde. Acta Paul. Enferm. 35, eAPE0163345. doi:10.37689/acta-ape/2022ao0163345
Contes, K. M., and Liu, B. M. (2025). Epidemiology, clinical significance, and diagnosis of respiratory viruses and their Co-Infections in the Post-COVID era. Pathogens 14, 262. doi:10.3390/pathogens14030262
Córdova-Dávalos, L. E., Hernández-Mercado, A., Barrón-García, C. B., Rojas-Martínez, A., Jiménez, M., Salinas, E., et al. (2022). Impact of genetic polymorphisms related to innate immune response on respiratory syncytial virus infection in children. Virus Genes 58, 501–514. doi:10.1007/s11262-022-01932-6
COVID-19 Excess Mortality Collaborators (2022). Estimating excess mortality due to the COVID-19 pandemic: a systematic analysis of COVID-19-related mortality, 2020-21. Lancet 399, 1513–1536. doi:10.1016/S0140-6736(21)02796-3
Cui, S., Han, Q., Zhang, R., Zeng, S., Shao, Y., Li, Y., et al. (2025). Integration of metabolomics methodologies for the development of predictive models for mortality risk in elderly patients with severe COVID-19. BMC Infect. Dis. 25, 10. doi:10.1186/s12879-024-10402-3
Dhangadamajhi, G., and Rout, R. (2021). Association of TLR3 functional variant (rs3775291) with COVID-19 susceptibility and death: a population-scale study. Hum. Cell 34, 1025–1027. doi:10.1007/s13577-021-00510-6
Djordjevic, A., Zivkovic, M., Koncar, I., Stankovic, A., Kuveljic, J., and Djuric, T. (2022). Tag variants of LGALS-3 containing haplotype block in advanced carotid atherosclerosis. J. Stroke Cerebrovasc. Dis. 31, 106212. doi:10.1016/j.jstrokecerebrovasdis.2021.106212
Dos Santos, B. R. C., dos Santos, L. K. C., Ferreira, J. M., dos Santos, A. C. M., Sortica, V. A., and de Souza Figueiredo, E. V. M. (2025). Toll-like receptors polymorphisms and COVID-19: a systematic review. Mol. Cell Biochem. 480, 2677–2688. doi:10.1007/s11010-024-05137-3
Evdokimov, A. V., Suslova, T. A., Belyaeva, S. V., Burmistrova, A. L., and Stashkevich, D. S. (2021). Polymorphism of TLR genes and the course of COVID-19 bilateral pneumonia. Med. Acad. J. 21, 57–66. doi:10.17816/MAJ90324
Feng, S., Song, F., Guo, W., Tan, J., Zhang, X., Qiao, F., et al. (2022). Potential Genes associated with COVID-19 and comorbidity. Int. J. Med. Sci. 19, 402–415. doi:10.7150/ijms.67815
Flores-González, J., Monroy-Rodríguez, Z., Falfán-Valencia, R., Buendía-Roldán, I., Fricke-Galindo, I., Hernández-Zenteno, R., et al. (2024). Variants rs3804099 and rs3804100 in the TLR2 gene induce different profiles of TLR-2 expression and cytokines in response to spike of SARS-CoV-2. Int. J. Mol. Sci. 25, 11063. doi:10.3390/ijms252011063
Fronteira, I., Mathews, V., dos Santos, R. L. B., Matsumoto, K., Amde, W., Pereira, A., et al. (2024). Impacts for health and care workers of Covid-19 and other public health emergencies of international concern: living systematic review, meta-analysis and policy recommendations. Hum. Resour. Health 22, 10. doi:10.1186/s12960-024-00892-2
Gudmundsson, S., Singer-Berk, M., Watts, N. A., Phu, W., Goodrich, J. K., Solomonson, M., et al. (2022). Variant interpretation using population databases: lessons from gnomAD. Hum. Mutat. 43, 1012–1030. doi:10.1002/humu.24309
Hakim, M., Khattak, F. A., Muhammad, S., Ismail, M., Ullah, N., Atiq Orakzai, M., et al. (2021). Access and use experience of personal protective equipment among frontline healthcare workers in Pakistan during the COVID-19 emergency: a cross-sectional study. Health Secur. 19, 140–149. doi:10.1089/hs.2020.0142
Hall, T. A. (1999). “BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT,” in Nucleic acids symposium series (Oxford, United Kingdom: Oxford University Press), 95–98. [London].
Helsinki Declaration (1964). World Medical Association’s Declaration of Helsinki (with later amendments): ethical principles for medical research involving human subjects. World Med. Assoc. Hels.
Kang, H. (2021). Sample size determination and power analysis using the G*Power software. J. Educ. Eval. Health Prof. 18, 17. doi:10.3352/jeehp.2021.18.17
Khanmohammadi, S., and Rezaei, N. (2021). Role of toll-like receptors in the pathogenesis of COVID-19. J. Med. Virol. 10, 2735–2739. doi:10.1002/jmv.26826
Lins, T. C., Vieira, R. G., Abreu, B. S., Grattapaglia, D., and Pereira, R. W. (2010). Genetic composition of Brazilian population samples based on a set of twenty-eight ancestry informative SNPs. Am. J. Hum. Biol. 22, 187–192. doi:10.1002/ajhb.20976
López-Juárez, P., Serrano-Oviedo, L., Pérez-Ortiz, J. M., García-Jabalera, I., Bejarano-Ramírez, N., Gómez-Romero, F. J., et al. (2021). Comparative study of the COVID-19 admissions between first and second wave in a cohort of 1,235 patients. Rev. Esp. Quimioter. 34, 387–389. doi:10.37201/req/005.2021
Machado, M. H., Teixeira, E. G., Freire, N. P., Pereira, E. J., and Minayo, M. C. de S. (2023). Deaths of doctors and nursing staff by COVID-19 in Brazil: a sociological approach. Ciência and Saúde Coletiva 28, 405–419. doi:10.1590/1413-81232023282.05942022
Marcelino, B. dos R., Vieira, M. C. da S., Silva, M. J. A., Silva, L. C. S. S. da, Gurrão, E. P. da C., Santos, E. C. dos, et al. (2025). Study of TNF-α, IFN-γ, IL-10, TGF-β and IL-6 gene polymorphisms in a cohort of professionals who worked in the first pandemic wave in the Brazilian Amazon. CRI 45 45, 39–61. doi:10.1615/CritRevImmunol.2024055001
Minashkin, M. M., Grigortsevich, N. Y., Kamaeva, A. S., Barzanova, V. V., Traspov, A. A., Godkov, M. A., et al. (2022). The role of genetic factors in the development of acute respiratory viral infection COVID-19: predicting severe course and outcomes. Biomedicines 10, 549. doi:10.3390/biomedicines10030549
Pena, S. D. J., Santos, F. R., and Tarazona-Santos, E. (2020). Genetic admixture in Brazil. Am. J. Med. Genet. Part C Seminars Med. Genet. 184, 928–938. doi:10.1002/ajmg.c.31853
Quigley, A. L., Stone, H., Nguyen, P. Y., Chughtai, A. A., and MacIntyre, C. R. (2021). Estimating the burden of COVID-19 on the Australian healthcare workers and health system during the first six months of the pandemic. Int. J. Nurs. Stud. 114, 103811. doi:10.1016/j.ijnurstu.2020.103811
Rodrigues de Moura, R., Coelho, A. V. C., de Queiroz Balbino, V., Crovella, S., and Brandão, L. A. C. (2015). Meta-analysis of Brazilian genetic admixture and comparison with other Latin America countries. Am. J. Hum. Biol. 27, 674–680. doi:10.1002/ajhb.22714
Rodrigues, M. M., Trindade, J. Q., Lima, M. B. M., Silva, M. J. A., Kerr, L. R. F. S., Kendall, B. C., et al. (2024). Lived reality and impacts of COVID-19 on nurses who provided care in the first wave on Brazilian Amazon. JOJNHC 13. doi:10.19080/JOJNHC.2024.13.555869
Salamaikina, S., Karnaushkina, M., Korchagin, V., Litvinova, M., Mironov, K., and Akimkin, V. (2023). TLRs gene polymorphisms associated with pneumonia before and during COVID-19 pandemic. Diagnostics 13, 121. doi:10.3390/diagnostics13010121
Sayers, E. W., Bolton, E. E., Brister, J. R., Canese, K., Chan, J., Comeau, D. C., et al. (2023). Database resources of the National Center for biotechnology information in 2023. Nucleic acids Res. 51, D29–D38. doi:10.1093/nar/gkac1032
Schurz, H., Daya, M., Möller, M., Hoal, E. G., and Salie, M. (2015). TLR1, 2, 4, 6 and 9 variants associated with tuberculosis susceptibility: a systematic review and meta-analysis. PLOS ONE 10, e0139711. doi:10.1371/journal.pone.0139711
Secolin, R., De Araujo, T. K., Gonsales, M. C., Rocha, C. S., Naslavsky, M., Marco, L. D., et al. (2021). Genetic variability in COVID-19-related genes in the Brazilian population. Hum. Genome Var. 8, 15. doi:10.1038/s41439-021-00146-w
Sharif, E., Al-Wakeel, M., Mohamed, A., kerkadi, A., and Rizk, N. (2020). TLR4 receptor D299G/T399I haplotype polymorphism is associated with insulin resistance in Obese female subjects. Genes (Basel) 11, 814. doi:10.3390/genes11070814
Sharpe, D. (2015). Chi-square test is statistically significant: now what? Pract. Assess. Res. Eval. 20, 8.
Silva, M. J. A., Silva, C. S., Da Silva Vieira, M. C., Dos Santos, P. A. S., Frota, C. C., Lima, K. V. B., et al. (2023a). The relationship between TLR3 rs3775291 polymorphism and infectious diseases: a meta-analysis of case-control studies. Genes 14, 1311. doi:10.3390/genes14071311
Silva, M. J. A., Silva, C. S., Marinho, R. L., Cabral, J. G., Gurrão, E. P. D. C., Dos Santos, P. A. S., et al. (2023b). Analysis of epidemiological factors and SNP rs3804100 of TLR2 for COVID-19 in a cohort of professionals who worked in the first pandemic wave in Belém-PA, Brazil. Genes 14, 1907. doi:10.3390/genes14101907
Silva, M. J. A., Vieira, M. C. da S., Souza, A. B., dos Santos, E. C., Marcelino, B., dos, R., et al. (2023c). Analysis of associations between the TLR3 SNPs rs3775291 and rs3775290 and COVID-19 in a cohort of professionals of Belém-PA, Brazil. Front. Cell Infect. Microbiol. 13, 1320701. doi:10.3389/fcimb.2023.1320701
Teräsjärvi, J., Kainulainen, L., Peltola, V., Mertsola, J., Hakanen, A., and He, Q. (2024). Genetic polymorphisms of TLR1, TLR2, TLR3 and TLR4 in patients with recurrent or severe infections. Int. J. Immunogenetics 51, 242–251. doi:10.1111/iji.12676
UNA-SUS (2020). Organização Mundial de Saúde declara pandemia do novo Coronavírus. Available online at: https://www.unasus.gov.br/noticia/organizacao-mundial-de-saude-declara-pandemia-de-coronavirus (Accessed September 6, 2022).
Vega, J. A., Villegas-Ospina, S., Aguilar-Jiménez, W., Rugeles, M. T., Bedoya, G., and Zapata, W. (2017). Haplotypes in CCR5-CCR2, CCL3 and CCL5 are associated with natural resistance to HIV-1 infection in a Colombian cohort, biomedica 37. 267, 273. doi:10.7705/biomedica.v37i3.3237
Von Elm, E., Altman, D. G., Egger, M., Pocock, S. J., Gøtzsche, P. C., Vandenbroucke, J. P., et al. (2014). The strengthening the reporting of observational Studies in epidemiology (STROBE) statement: guidelines for reporting observational studies. Int. J. Surg. 12, 1495–1499. doi:10.1016/j.ijsu.2014.07.013
Wang, B.-G., Yi, D.-H., and Liu, Y.-F. (2015). TLR3 gene polymorphisms in cancer: a systematic review and meta-analysis. Chin. J. Cancer 34, 272–284. doi:10.1186/s40880-015-0020-z
Xie, J., and Zhu, Y. (2020). Association between ambient temperature and COVID-19 infection in 122 cities from China. Sci. Total Environ. 724, 138201. doi:10.1016/j.scitotenv.2020.138201
Keywords: COVID-19, TLR2, TLR3, Brazil, single-nucleotide polymorphism
Citation: Silva MJA, Figueira LRT, Sardinha DM, Costa da Cruz E, Andrade NCO, Bispo SKDS, Augusto Ferreira Dos Anjos T, Dos Santos EC, Garcia AJP and Lima LNGC (2025) Analyses of haplotypes of TLR2 and TLR3 genes for COVID-19 prognosis in a cohort of professionals who worked in the first pandemic wave in Belém-PA, Brazil. Front. Genet. 16:1659269. doi: 10.3389/fgene.2025.1659269
Received: 03 July 2025; Accepted: 15 September 2025;
Published: 02 October 2025.
Edited by:
Anteneh Getachew, University of Alabama at Birmingham, United StatesReviewed by:
Abhishek Kumar, Washington University in St. Louis, United StatesRakesh Kumar Sharma, University of Alabama at Birmingham, United States
Copyright © 2025 Silva, Figueira, Sardinha, Costa da Cruz, Andrade, Bispo, Augusto Ferreira Dos Anjos, Dos Santos, Garcia and Lima. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Marcos Jessé Abrahão Silva, amVzc2VhYnJhaGFvMTBAZ21haWwuY29t; Luana Nepomuceno Gondim Costa Lima, bHVhbmFsaW1hQGllYy5nb3YuYnI=