Involvement of IL10 rs1800872 and rs1800896 in the Acquisition and Progression of Dengue Virus Infection to Severe Forms of the Disease in Burkina Faso ()
1. Introduction
Dengue is the most important viral disease transmitted to humans by the bite of the Aedes aegypti mosquito in most tropical and subtropical countries of the world [1] [2]. There are an estimated 50 to 100 million cases of dengue fever worldwide every year [3]. Isolated for the first time in 1925 in Burkina Faso, the country recorded 79,867 suspected cases including 34,687 probable cases and 349 deaths in 2023 [4] [5]. Dengue virus (DENV) is a single-stranded RNA virus of positive polarity of 10.7 Kb, belonging to the genus Flavivirus [6].
There are 5 serotypes of the virus which are genetically related but antigenically distinct, classified as DENV 1 to 5 [7]. DENV 1, 2, 3 and 4 are each capable of causing self-limiting fevers to fatal states. Infection with any of the four viral serotypes confers protective immunity against reinfection with the same serotype only, while subsequent infection with other serotypes can result in severe dengue [6]. DENV 5, isolated in 2013, follows a purely sylvatic cycle [7]. The main target cells of DENV in vitro are monocytes, macrophages and dendritic cells [8]. Generally, dengue can be asymptomatic or cause mild illness during a primary infection. However, a second infection can be fatal [9] [10]. Environmental factors, the serotype/genotype of the dengue virus, the immune response and the genetic make-up of the host are considered to influence the development of the clinical manifestations of dengue, as well as the severity of the disease [11] [12]. And the host immune response has been highlighted as a genetic marker of the disease with the production and release of several pro-inflammatory, antiviral and immunoregulatory cytokines [11] [13]. The host immune response is primarily determined by the recognition of the virus by the host cells and the subsequent antiviral response. Th1 cells are typically associated with an acute phase of fever, while Th2 cells are expressed throughout the infection, leading to endothelial cell dysfunction and causing plasma leakage. The activation of T lymphocytes triggers a cytokine cascade. During DENV infection, these cytokines can be secreted by macrophages and monocytes [7].
Previous studies shown that cytokines have a significant impact on the severity of dengue fever, including alpha tumor necrosis factor (TNF-α), interleukin 6 (IL6), interleukin 10 (IL10), interleukin 17 (IL17) and gamma interferons (IFN-γ) [10] [12] [14]. In addition, single-nucleotide polymorphisms (SNPs) in the genes for these cytokines have contributed significantly to our understanding of pathophysiology and the role of host genetics in dengue infection [15] [16]. Serum concentrations of TNF-α, IL6, IL2, and INF-γ are elevated during the first three days of infection of DENV, while high levels of IL10, IL5, and IL4 appear later. The interaction of these cytokines, which can either increase or decrease their production, may be responsible for triggering the inflammatory process in infectious diseases [7].
Interleukin 10 (IL10) plays a pleiotropic immune and inflammatory regulatory role in infectious diseases and is produced by monocytes, macrophages, dendritic cells, regulatory T cells (Tregs) [13] [17]. IL10 contains three distinct SNPs in its promoter region, at nucleotide positions 1082, 819, 592 [17]. In the pathogenesis of DENV, IL-10 has immunomodulatory activity, resulting in persistent viral infection, promoting inflammation and disease severity [12]. Polymorphisms of the IL10 gene could influence the Antibody-Dependent Enhancement reaction (ADE) [3].
Given that cytokine production is influenced by genetic polymorphisms and that different genetic variants have been associated with dengue disease, it would be relevant to investigate the potential role of cytokine genetic polymorphisms, in the pathogenesis of dengue.
The aim of our study was to evaluate the role of interleukin 10 (IL-10) gene polymorphisms rs1800872 and rs1800896 in the pathogenesis of dengue fever and to assess their involvement in the progression of infection to severe forms of the disease in a Burkinabe population.
2. Materials and Methods
2.1. Ethics Approval and Consent to Participate
The study was approved by the Health Research Ethics Committee of the Burkina Faso Ministry of Health: No. 2022-02-034 and written informed consent was obtained from participants prior to blood sampling.
2.2. Type and Study Population
This case-control study, involved patients of all ages, including children and blood donors from various professions and social categories in Ouagadougou. Blood samples from dengue cases were collected in the laboratories of the Hopital Saint Camille de Ouagadougou (HOSCO) and the Pietro Annigoni Biomolecular Research Center (CERBA) from October 2021 to June 2022. Blood samples from subjects with no known history of dengue, received during the collection period, served as controls. Patients with well-documented pre-existing pathologies were excluded from this study.
2.3. Genomic DNA Extraction
Peripheral venous blood was collected from each individual in tubes containing EDTA anticoagulant. Genomic DNA was extracted using the FAVORGEN Mini kit, following the protocol provided by the manufacturer. It was then quantified and checked for purity using the BioDrop apparatus. The DNA was stored at −20˚C until analysis.
2.4. Genotyping
Two bi-allelic polymorphisms of the IL-10 promoter were detected by the PCR-RFLP method using primers that amplify a short fragment of DNA containing the polymorphisms. A pair of primers (sense: 5'-GGTGAGCACTACCTGACTAGC-3' and antisense: 5'-CCTAGGTCACAGTGACGTGG) was used to amplify the IL10-592 fragment (C/A; rs1800872) and for the IL10-1082 (G/A; rs1800896), we used the primer pair (sense: 5'-CCAGATATCTGAAGAAGTCCTG-3'; antisense: 5'-CTCTTACCTATCCCTACTTCC-3'). Amplification was performed in a 20 μL reaction volume containing 1 μL of each primer, 100 ng of template DNA and 10 μL of 5X FIREPOL® master mix and 7 μL of water. Polymerase chain reactions were run for 30 cycles according to the following program: initial denaturation for 5 minutes at 95˚C, denaturation for 30 seconds at 95˚C, hybridization for 30 seconds at 55˚C, extension for 30 seconds at 72˚C and final extension for 5 minutes at 72˚C.
2.5. Enzymatic Digestion
The amplification products, either 5 μL of each sample, were digested with the restriction enzymes Fast Digest RSaI (Thermo Scientific TM) and Fast Digest MnlI (Thermo Scientific TM) for the IL-10-592 and IL10-1082 genes respectively (see Table 1) at 37˚C for 30 minutes in a water bath and then subjected to electrophoresis on a 2% agarose gel and stained with ethidium bromide. The enzymatic digestion reaction of each sample was carried out in a total reaction volume of 25 µL containing the enzyme buffer, the enzyme, pure water, and the PCR product. The IL10 A-1082 allele also gave 2 fragments of 134 and 65 bp and the IL10 G-1082 allele gave 3 fragments of 112 bp, 65 bp and 22 bp (Figure 1). The IL10 A-592 allele gave 2 fragments of 236 and 176 bp and the IL10 C-592 allele gave a single fragment of 412 bp (Figure 2).
Table 1. PCR products from enzymatic digestion.
Genes |
Polymorphism |
Size (bp) |
Restriction enzyme |
Concentration of electrophoresis gel |
Expected digestion products |
IL10 |
IL10-592C/A; (rs1800872) |
412 |
RSAI |
2% |
AA (264 pb et 176 pb) AC (412 pb, 236 pb et 176 pb) CC (412 pb) |
|
IL10-1082G/A; (rs1800896) |
189 |
MnLI |
|
AA (134 pb et 65 pb) AG (134 pb~112 pb55 pb et 22pb) GG (112 pb 55 pbet 22 pb) |
Figure 1. Genotypic illustration of the IL10-1082 G/A polymorphism in agarose gel.
Figure 2. Illustration of IL10-592 C/A genotypic polymorphisms in agarose gel.
2.6. Statistical Analysis
The data were analysed using R software. The chi-square (χ2) test was used for frequency comparisons. Odds ratios (OR) and 95% confidence intervals (CI) were calculated. Results were considered statistically significant for p value less than 0.05.
3. Results
3.1. Study Population
A total of 246 peoples were recruited. These included 129 people with clinical signs suggestive of dengue fever, confirmed by diagnostic tests at the medical biology laboratory, and 117 peoples who, based on biological tests, had not been in contact with DENV (Table 2). The study population consisted of 43.09% men infected with DENV compared to 56.91% women. Table 2 also shows that the age group most affected by DENV infection was over 40, with a DF infection rate of 66.89%. The criteria used for cases of severe dengue fever included, plasma leakage characterized by a rise in hematocrit, severe thrombocytopenia, and symptoms noted at the time of consultation, such as abdominal pain with vomiting and splenomegaly. The majority of the population (68.70%, (169/246)) was over 40. This age group also had the highest proportion of cases of dengue fever (65.89%), compared with 12.40% in the least active age group.
Table 2. Socio-demographic data of the study population.
Variables |
DF n (%) |
DS n (%) |
DF + DS n (%) |
Controls n (%) |
Total n (%) |
OR |
95% CI |
p-value |
N |
98 |
31 |
129 |
117 |
246 |
|
|
|
Gender, n (%) |
|
|
|
|
|
|
|
|
Male |
41 (41.84) |
15 (48.39) |
56 (43.41) |
50 (42.74) |
106 (43.09) |
- |
- |
- |
Female |
57 (58.16) |
16 (51.61) |
73 (56.59) |
67 (57.26) |
140 (56.91) |
1.303 |
[0.579 - 2.932] |
0.5397 |
Year |
|
|
|
|
|
|
|
|
0 - 19 |
14 (14.29) |
2 (6.45) |
16 (12.40) |
11 (9.40) |
27 (10.98) |
- |
- |
- |
20 - 39 |
19 (19.39) |
9 (29.03) |
28 (21.71) |
22 (18.80) |
50 (20.33) |
0.309 |
[0.0283 - 1.8449] |
0.1621 |
≥40 |
65 (66.33) |
20 (64.52) |
85 (65.89) |
84 (71.79) |
169 (68.70) |
0.3016 |
[0.0562 - 1.6190) |
0.2775 |
3.2. Serological Characteristics of the Study Population
Table 3 shows the proportions of DF and DS according to the stage of infection. Primary DENV infection was 46.94% (46/98) and 12.90% (4/31), respectively. Among severe dengue cases, 45.16% (14/31) had at least one previous infection, with 22.58% in the acute phase. The proportion of the population (42.86% (42/98)) had antibody serology corresponding to a phase of infection secondary to another type of DENV or in the recovery phase.
Table 3. Serological data specific to DENV infection in the study population.
Serologicalmarkers DENV |
Controls |
Dengue fever (DF) |
Severe dengue (DS) |
Total |
N |
% |
N |
% |
N |
% |
N |
% |
AgNS1 |
0 |
0 |
61 |
62.24 |
17 |
54.84 |
78 |
31.71 |
Ac-IgM−/IgG− |
117 |
100 |
46 |
46.94 |
04 |
12.90 |
167 |
67.89 |
Ac-IgM−/IgG+ |
0 |
0 |
42 |
42.86 |
14 |
45.16 |
56 |
22.76 |
Ac-IgM+/IgG− |
0 |
0 |
01 |
01.02 |
07 |
22.58 |
08 |
03.25 |
Ac-IgM+/IgG+ |
0 |
0 |
09 |
09.18 |
06 |
19.35 |
15 |
06.10 |
TOTAL |
117 |
100 |
98 |
100 |
31 |
100 |
246 |
100 |
3.3. Relationship between Hematological Parameters and Dengue
Acquisition
The following figures show three blood cell count (CBC) parameters in a whisker plot to assess their involvement in DF. Figure 3 and Figure 4 show the median lymphocyte and hematocrit levels, which are almost identical in both groups. Figure 5, on the other hand, shows the platelet count in DF cases versus controls. From these graphs, we can see that the drop in platelet count is linked to dengue fever (p = 0.00).
Figure 3. Moustache plot of lymphocyte counts in the study population.
Figure 4. Graph of haematocrit levels in the study population.
Figure 5. Involvement of thrombocytopenia in dengue disease.
3.4. Genotypic Distribution of IL10-592 and IL10-1082 Genes
The frequencies of IL10-592 and IL10-1082 gene polymorphisms were assessed, and the genotypic distribution was in Hardy-Weinberg equilibrium in both patients and controls (p-value < 0.05) (see Table 4). Statistical analysis of IL10-592 and IL10-1082 polymorphisms indicates their involvement in DF infection. The AA genotype (48.36%, p < 0.05) of the IL10-592 SNP is associated with the risk of developing the DS form, while the CA heterozygote (21.31%) tends not to promote DF infection and even slows its progression to the DS form.
Table 4. IL10 genotype frequencies in different forms of dengue fever.
SNP |
Genotypes |
DF/Control n (%) |
DS/Controls n (%) |
Controls n (%) |
DS/DF n (%) |
IL-10-592
(rs1800872) |
CC |
19 (20.21) |
7 (25.00) |
49 (41.88) |
26 (21.31) |
CA |
31 (32.98) |
6 (21.43) |
41 (35.04) |
37 (30.33) |
AA |
44 (46.81) |
15 (53.57) |
27 (23.08) |
59 (48.36) |
HWE
p-value |
0.007019 |
0.009695 |
0.003803 |
0.0001697 |
IL10-1082(rs1800896) |
AA |
18 (19.35) |
6 (19.35) |
29 (25.66) |
24 (19.35) |
AG |
26 (27.96) |
7 (22.58) |
38 (33.63) |
33 (26.61) |
GG |
49 (52.69) |
18 (58.06) |
46 (40.71) |
67 (54.03) |
HWE
p-value |
0.00042 |
0.01144 |
0.001011 |
0.00001418 |
As for the SNP of IL10-1082, its GG genotype would be associated with a predisposition to progression from DF to DS, while its AG genotype would confer greater protection against progression to the severe form of dengue.
3.5. IL10 Genotype Frequencies According to Circulating Dengue
Serotypes in Burkina Faso
In the present study, DENV1 was the most prevalent, and DENV4 was the least frequent (Table 5). However, this difference was not statistically significant in terms of the relationship between the genotypes of the genes studied and the different DENV serotypes circulating in Burkina Faso.
Table 5. IL10 genotype frequencies by DENV serotype.
SNP |
Genotypes |
DENV 1 n (%) |
DENV 2 n (%) |
DENV 3 n (%) |
DENV 4 n (%) |
IL10-592(rs1800872) |
CC |
0 (0.00) |
2 (1.55) |
0 (0.00) |
0 (0.00) |
CA |
7 (24.14) |
4 (13.19) |
1 (3.45) |
0 (0.00) |
AA |
7 (24.14) |
5 (17.24) |
4 (13.19) |
1 (3.45) |
HWE p-value |
0.513 |
0.5377 |
1.00 |
NA |
IL10-1082(rs1800896) |
AA |
2 (6.90) |
0 (0.00) |
1 (3.45) |
0 (0.00) |
AG |
3 (10.34) |
4 (13.19) |
0 (0.00) |
0 (0.00) |
GG |
8 (6.20) |
5 (17.24) |
5 (17.24) |
1 (3.45) |
HWE p-value |
0.1652 |
1.00 |
0.09091 |
NA |
3.6. Genotypic and Allelic Polymorphisms of IL10 Genes and
Dengue Pathogenesis
The allelic distribution shows a higher frequency of the allele A of both genes in controls than in DS and DF patients (59.40% for IL10-592 and 42.48% for IL10-1082). Indeed, with p-value < 0.05, we could suggest that the allele A would confer less risk of infection in DF and would not favor progression from DF to the DS form. The C allele of the IL10-592 gene and the G allele of the IL10-1082 genotype had higher frequencies in patients with DF compared with DS (63.30% and 66.67%) (Table 6).
Table 6. Allelic frequency and involvement of IL10-592 and IL10-1082 genes in dengue pathogenesis.
SNP |
Genotypes and alleles |
DF n (%) |
DS n (%) |
Controls n (%) |
OR |
CI (95%) |
p-value |
IL10-592(rs1800872) |
CA vs CC |
31 (62.00) |
6 (46.15) |
41 (45.56) |
0.591 |
[0.305 - 1.132] |
0.108 |
AA vs CC |
44 (69.84) |
15 (68.18) |
27 (35.53) |
0.246 |
[0.125 - 0.471] |
0.000 |
CA + AA vs CC |
75 (79.79) |
21 (75.00) |
68 (58.12) |
0.378 |
[0.212 - 0.664] |
0.001 |
C (%) |
119 (63.30) |
36 (64.29) |
95 (40.60) |
- |
- |
- |
A (%) |
69 (36.70) |
20 (35.71) |
139 (59.40) |
2.523 |
[1.747 - 3.662] |
0.000 |
IL10-1082(rs1800896)G/A |
AG vs GG |
26 (34.67) |
7 (28.00) |
38 (45.24) |
1.671 |
[0.918 - 3.062] |
0.089 |
AA vs GG |
18 (26.87) |
6 (25.00) |
29 (38.67) |
1.752 |
[0.906 - 3.418] |
0.090 |
AG + AA vs GG |
44 (47.31) |
13 (41.94) |
67 (59.29) |
1.707 |
[1.020 - 2.873] |
0.040 |
G (%) |
124 (66.67) |
43 (69.36) |
130 (57.52) |
- |
- |
- |
A (%) |
62 (33.33) |
19 (30.64) |
96 (42.48) |
1.521 |
[1.046 - 2.215] |
0.027 |
4. Discussion
IL-10 is a major anti-inflammatory cytokine, essential in controlling the host immune response by regulating the production of several pro-inflammatory cytokines. It has been associated with several diseases and is considered an important immunoregulatory mediator produced by monocytes, dendritic cells and T and B lymphocytes [12] [14]. IL-10 is also likely to be regulated at transcriptional level by several polymorphisms, including 819 and 1082 [12].
The aim of our study was to evaluate the involvement of IL10 gene polymorphisms in the pathogenesis of dengue disease and its progression to severe forms in a Burkinabe population. In the present study, the different IL10 polymorphisms showed that IL10-592 and IL10-1082 were statistically associated with DF and even progression from DF to the severe DS form.
Examining the genotypic distribution of IL10-592 polymorphisms, we found that the AA genotype (46.81% of DF cases and 53.57% in DS cases, with a p-value < 0.001) could be associated with the risk of developing DS, while the CA heterozygote (32.98%s of DF cases and 21.43% in DS cases, with a p-value < 0.001) could confer resistance in the progression to the DS form; and homozygous CC (20.21% of DF cases and 25.00% in DS cases, with a p-value < 0.001) would tend not to favor DF infection.
The GG genotype of the IL10-1082 gene could be a predisposing factor for DF and even progression to DS (52.69% DF and 58.06% DS). Whereas the AA genotype (19.35% DS/DF with p-value < 0.001) could be a protective factor against progression to DS. Gyeltshen et al. obtained a frequency of GG genotype associated with higher plasma levels of IL10-1082 [16]. According to the virus serotypes circulating in the study population, we obtained equal frequencies of AA and CA genotypes (24.14%) in serotype DENV1, which was the most frequently encountered. Serotype DENV4 was also present in the Burkinabe population, and the AA genotype of IL10-592 and the GG genotype of IL10-1082 were present at 3.45%. The study by Féitosa et al. in Brazil also shows an association of high IL10 levels with DENV1 infections [18]. The association of AG + AA mutant genotypes (47.31% of DF cases and 41.94% of DS cases; OR = 1.707; CI [1.020 - 1.873] with p-value = 0.040) of the IL10-1082 gene could confer slight protection against DS.
The allele C allele of IL10-592 is proportionally predominant in dengue cases with all CA+AA mutated genotypes, compared with the wild-type AA genotype, which would be a factor slowing down progression to DS (79.79% of DF cases and 75.00% in DS cases; OR = 0.378; CI [0.212 - 0.664] with a p-value = 0.001). Also, the A allele of IL10-1082 (OR = 1.521 CI [1.046 - 2.215], p-value = 0.027) was associated with the development of dengue and would not be indifferent in the progression to the DS form; this may be in line with the study of cytokine polymorphism frequencies in a population from Western Europe, Africa, Asia, the Middle East and South America, which obtained high frequencies of the IL10-1082 A allele, and also with the study on a Bhutanese population which found a frequency of the IL10-592 C allele at 53.86% [16] [17].
Our results are in line with a large number of studies that have shown the association of the IL10 gene with dengue fever. Such as those obtained in Indonesia, which obtained higher IL10 levels in DF than in DS [10]; in Cuba, which found that SNPs-1082 (G/A), and -592 (C/A) were associated with DHF risk [14]; in Sri Lanka, which confirmed the association between the same genotypes and DHF, suggesting a risk factor for the development of DS, while alternative combinations were associated with protection against dengue hemorrhagic fever [11]; and also in Brazil, which obtained a distribution of IL10-1082 genotypes [18].
Knowing that one study had associated these genotypes with low IL10-expression we can suggest with other authors that dengue-associated genotypes could be responsible for the ineffective immune response and resulting severity [7] [11] [14].
5. Conclusion
The results of this study showed a statistically significant association of IL10-592, IL-1082 gene polymorphisms in dengue virus infection in the Burkinabe population, and also their implications in disease progression to severe forms of the disease. Despite the studies carried out to understand the pathogenesis of dengue disease, questions remain as to how many genes contribute to susceptibility to dengue, and to what extent these genes interact with each other to cause severe disease. Further studies with other genes and also with the involvement of other arboviruses, particularly in cases of co-infections, in different population groups are essential for a complete understanding of the genetic associations with dengue infections.
Acknowledgements
We would like to thanks The World Academy of Sciences and the Swedish International Development Cooperation Agency (Sida) for funding this research through the grant No. 17-403 RG/BIO/AF/AC_I–FR3240297757. We thank also the UNESCO chaire in “Génie génétique et Biologie Moléculaire” for the technical support.
Informed Consent Statement
See ethical approval.
Data Availability
The datasets used and/or analyzed during the current study are available from the corresponding authors on reasonable request.
Author Contributions
Conceptualization, data curation, formal analysis, methodology, writing—original draft, writing—review and editing by ASAT, FT, JS and FWD.
FWD, ASAT, FT, LT, SOTB, TWCO, RAO, PDI, MS, PSD and JS contributed for the samples analysis, validation, review and editing the manuscript.
FWD mobilized the funding for the study.
All authors read and approved the final version of the manuscript.