Abstract
Dengue is considered one of the most challenging public health threats in the world. Infection may be clinically asymptomatic but can result in severe forms. The indoleamine 2,3 dioxygenase (IDO) gene encodes one of first enzymes (IDO) of the kynurenine pathway. This study aimed to verify the association between G2431A IDO1 gene single nucleotide polymorphism (SNP) (rs3739319) and dengue fever development. We included 299 dengue-infected individuals in the study and 96 dengue-free controls. We collected clinical and diagnostic test data and divided the patients with dengue infection into three groups, based on World Health Organization (WHO) criteria: 131 Dengue without warning signs (DWOS), 143 Dengue with warning signs (DWS), and 25 severe dengue (SD). We genotyped 193 of the dengue cases using quantitative polymerase chain reaction to the SNP rs3739319. The other 106 dengue cases and 96 dengue-free controls had previously been genotyped using the Illumina Genotyping Kit. Genotyping of the infected patients revealed frequencies of 106 GG (35.4%), 126 GA (42.1%), and 67 AA (22.4%), whereas the nondengue exposed control group showed similar frequencies, 29 GG (30.2%), 52 GA (54.2%), and 15 AA (15.6%). Under risk analysis we found that AA genotype patients had a higher risk of developing SD in a codominant model (AA × GG; odds ratio [OR] = 11.5-fold in comparison to non-SD group -DWOS and -DWS patients; confidence interval [CI] = 0.02–0.32; Yates correction = 1.9e-05) and in a recessive model (AA × AG+GG; OR = 9.41; CI = 3.62–26.7; Yates correction = 4.8e-08). An allelic model reinforced the association between A allele and SD phenotype development that was found in the SD versus DWOS+DWS analysis (OR = 3.59; CI = 1.50–9.56; Yates correction = 0.0033). Our data show an association between the IDO G2431A variant and the risk for SD. This SNP may be relevant for further investigation into disease mechanisms and host factors in future genetic and pathophysiological studies.
Introduction
Dengue virus (DENV) is one of the biggest public health challenges in the world (5). Estimates point to 50–100 million cases annually worldwide, encompassing tropical and subtropical regions of the globe, and infection may cause a spectrum of disease: from asymptomatic to severe dengue (SD) forms (5,24). The DENV is enveloped, having an icosahedral capsid whose genetic material is an 11 Kb RNA single strand (12,21). It belongs to the Flavivirus genus, Flaviviridae family, and is transmitted by hematophagous arthropod vectors, usually from Aedes genus (especially Aedes aegypti and Aedes Albopictus) (20).
Dengue clinical pathogenesis is controversial (21). The World Health Organization (WHO) classification system, developed in 2009, groups cases on the basis of clinical features. Dengue without warning signs (DWOS) cases are characterized by headache, myalgia, arthralgia, nausea, vomiting, rash, and other common symptoms to other arboviruses. Dengue with warning signs (DWS) encompasses disease where patients develop mucosal bleeding, ascites, abdominal pain, and persistent vomiting in the defervescence phase of the disease. In SD patients presents bleeding, shock, and severe organ involvement (26,45).
Several important factors have been associated with disease progression, including environment, immune status, virus strain, and host genetic background (40). For the latter, it is relevant to mention the influence of single nucleotide polymorphisms (SNPs) on disease progression, beyond complex multiple signatures between SNPs (6,13,46). Some SNPs have been discovered as potential influencers of severe disease development, such as rs1799964 (TNFα), rs1143627 (IL1β), and rs1800871 (IL10), which are cytokine genes related to induction of fever during disease and increase of vascular permeability (6,37,38). Other SNPs from genes, such as rs222270 (VDR), rs4804803 (CD209), exon 1 and promoter SNPs (MBL2), rs3753394, and rs800292, belonging to the complement factor H gene (CFH), have also been described to contribute to the development of severe disease (1,7,15,31,43,44).
The gene indoleamine 2,3 dioxygenase (IDO1) is 14,9 Kb and encodes the IDO protein, an important regulator of the immune response (22). IDO is one of first enzymes of the kynurenine pathway that converts tryptophan to N-formylkynurenine by means of oxidation of its pyrrole ring (17). N-formylkynurenine is metabolized to kynurenine, which can be released or catabolized to generate new products, among them nicotinamide adenine dinucleotide (NAD+), the end product of the whole reaction.
In dengue, infected febrile individuals may have high IDO transcripts in blood cells, suggesting that this enzyme may be involved in the disease pathophysiology (9). Taking into account the connection between some innate immune genes (SNPs) and dengue development, and high IDO levels in dengue in acute dengue patients, this study aimed to verify the relationship between the SNP rs3739319 (IDO1 gene) and SD development. There are very few studies characterizing the relevance of the IDO SNP in dengue patients. We describe an association between the presence of the IDO1 G2431A variant and an increased risk for SD, highlighting the importance of host genetics to DENV pathophysiology. This SNP may be a useful target for further investigation in future genetic and pathophysiological studies.
Material and Methods
Patients and clinical definition
Patients with dengue-related symptoms ≥5 years of age, without any other disease that might interfere with diagnostic test results–such as Zika virus, another flavivirus—were screened from three hospitals in Recife, Brazil: Oswaldo Cruz Hospital, Santa Joana Hospital, Esperança Hospital and Instituto Materno Infantil de Pernambuco (IMIP). All enrolled patients gave their informed consent and the project was approved by the Ethics Committee of FIOCRUZ-PE (CEP/CPQAM no.11/11, C.A.A.E. 0009.0.095.000-11, IORG0001419). Clinical details of the patients obtained included age, sex, and serum anti-DENV IgG levels (Table 1).
Clinical Details of the Control Subgroup and Brazilian Dengue 3 Patients
All control patients are IgG negative.
DWOS, Dengue without warning signs; DWS, Dengue with warning signs; SD, Severe dengue.
Patients were examined from the first to the fourth day of the disease, with additional evaluations every 24–48 h. Each volunteer donated ∼20 mL of blood over the course of the first 30 days after symptom onset. DENV (serotype 3) was detected by reverse transcriptase–polymerase chain reaction (RT-PCR) according to Lanciotti et al. (27) and the virus isolation in C6/36 cell line was identified by immunofluorescence testing (19) with serotype-specific and anti-dengue monoclonal antibodies (Bio-Manguinhos, Fundação Oswaldo Cruz, Brazil).
Dengue cases were classified in accordance with the WHO guidelines, as described (45). After clinical and laboratorial evaluation, a subset of 299 well-characterized dengue cases was selected for genotyping of the SNP rs3739319 in IDO1 gene: 131 DWOS, 143 DWS, and 25 SD. The control group was formed of 96 healthy volunteers that were negative for dengue tests described above.
DNA extraction and patient genotyping
The IDO1 gene G2431A region of 96 healthy Brazilian blood donors and 299 dengue-positive selected was genotyped to identify the allele and genotype frequencies as described below. Genomic DNA was extracted from PBMCs of the patients by using the Qiagen Genomic Prep Mini Spin Kit following the manufacturer's protocol.
For 193 of the 299 dengue-infected cases, we used a TaqMan Applied Genotyping quantitative PCR (qPCR) Assay (C_9036149_10) to genotype the rs3739319 locus following the recommended protocol: 6.25 μL Master Mix (TaqMan Genotyping Master Mix; Applied Biosystems/Thermo Fisher), 0.625 μL probe, 3.125 μL sterile H2O, and 3 μL DNA to 13 μL of final volume. The amplification reaction was performed in Applied Biosystems 7500 machine, using the following cycling profile: 95°C/10 min, 40 cycles of: 92°C/15 sec and annealing/extension 60°C/1 min, besides the post-PCR at 60°C/30 sec.
We used the Illumina GoldenGate Genotyping Kit (Cat no. GT-901-1001) to genotype the other 106 dengue-infected cases and 96 controls included in the study. The Illumina system uses chips, with specific primers for regions to be studied, and genotypes for a locus of interest through a multiplex system associated with a scanning system. The BeadChip is read by a BeadArray Reader through two channels of high resolution lasers. This reading generates one file per channel that is interpreted by the iScan software to determine intensity values for each bead. Finally, Genome Studio analyzes the data in conjunction with bead pool map (* .bpm) or manifest file (* .bgx).
Statistical analysis
The risk analysis and plotting were carried out using the open source R statistical package version i386 3.5.0. To study the association between genotype groups and clinical dengue categories, Fisher's exact test was performed, using the R function fisher test, which was also used to find the associated odds ratio (OR), using a conditional maximum-likelihood method. Codominant model (GG vs. AA), dominant model (GG vs. GA+AA), recessive model (GG+GA vs. AA), overdominant model (GA vs. GG+AA), and allelic model (G vs. A) were used to investigate risk factors for SD with a hemorrhagic phenotype development. The Hardy–Weinberg equilibrium was evaluated using chi-square test and values were considered significant when p < 0.05. Yates correction test was applied to verify the possibility of risk maintenance in the population.
Results
Two hundred ninety-nine DNA samples from well-characterized dengue patients (DENV-3 genotype III, Srilankan–India strain) were genotyped for locus G2431A (rs3739319) in IDO1 gene. Dengue patients revealed genome frequencies of 106 GG (35.4%), 126 GA (42.1%), and 67 AA (22.4%), whereas the control group had frequencies of 29 GG (30.2%), 52 GA (54.2%), and 15 AA (15.6%) (Table 2).
IDO G2431A (rs3739319) Genotype/Allele Frequencies in Control Group and Brazilian Dengue 3 Patients
DEN, Dengue group; SNP, single nucleotide polymorphism.
An interesting fact, depicted in Table 2, is the high frequency of the AA genotype in SD patients (68%) compared with DWOS and DWS patients, with AA frequencies of 21.4% and 15.4%, respectively. This tendency was either observed for allele frequencies (Table 2), where A allele shows a frequency of 73%, much higher than dengue group (45%).
To determine if there was any association between severity of dengue disease and the 2431A locus (rs3739319 SNP), the five models cited in the methodology were used. We describe our results fully in Table 3. In the codominant model, AA genotype patients showed a risk of 11.5-fold to develop SD in comparison with DWOS/DWS group (OR = 11.5; confidence interval [CI] = 0.02–0.32; Yates correction = 1.9e-05). Corroborating with this result, the dominant model revealed that a GG genotype provides protection against disease progression to SD (OR = 0.22; CI = 0.04–0.78, Yates correction = 0.0194). In the same way, when recessive analyses were considered, the AA genotype exhibited an increased risk for SD (SD vs. control, SD vs. DWOS+DWS). One expressive result was seen when SD was compared with DWOS+DWS (OR = 9.41; CI = 3.62–26.7; Yates correction = 4.8e-08).
IDO G2431A (rs3739319) Allele/Genotype Risk Analysis (rs3739319) Between Control Group and Brazilian Dengue 3 Patients
The analyses were performed using the Fisher exact test. The results that presented statistical significance (P ≤ 0.05) are highlighted in bold, and these are also specified by the Yates correction.
NS, no significance; CI, confidence interval; OR, odds ratio.
The overdominant model showed that an AG genotype protects against SD (OR = 0.31; CI = 0.09–0.9; Yates correction = 0.033). This can be explained by the fact that only one copy of the G allele is sufficient to confer protection against SD, since this allele is dominant when compared with the A allele.
Finally, the risk analysis was done based on the allelic model. This reinforced the connection of the A allele with the development of SD, since there was a risk relation in the three analyzes performed (SD vs. control; SD vs. DWOS; SD vs. DWOS+DWS). The higher risk of A allele was found when SD versus DWOS+DWS analysis was performed (OR = 3.59; CI = 1.50–9.56; Yates correction = 0.0033).
Discussion
In this work we performed a study to understand the relationship between 2431A polymorphism (rs3739319) and severity of dengue disease. To do this, 299 dengue patients and 96 healthy controls were genotyped. It was found that patients with an AA genotype have a higher risk of developing SD, while the presence of a G allele is protective against this severe disease phenotype.
The IDO1 gene is located on chromosome 8 and encodes the enzyme indoleamine 2,3 dioxygenase (IDO), a very important immunoregulatory protein for the catabolism of the amino acid tryptophan (Trp). It is usually expressed by antigen-presenting cells under the stimulus of proinflammatory factors, such as tumor necrosis factor (TNF-α) and interferons (IFN) (8,10,23). In addition, its activity is associated with the functioning of various metabolic pathways and the maintenance of immune homeostasis (35,36).
The SNP G2431A is located in an intronic region between exons 9 and 10, near enhancer gene regulatory regions and promoter flanks. Its role in immunopathology of viral and nonviral diseases is still poorly understood, a coincident fact shared by other IDO1. Functional associations have been found between the SNP rs7820268 and a series of diseases, including systemic sclerosis, accompanied by a decrease in CD8+ Treg activity (42); aspirin-induced asthma (39) in a Japanese population; and, singly, with type 1 diabetes mellitus (DM1), with an OR of 1.69 (p = 0.014), whereas in the haplotype model (H1/H1 haplotype) the OR was 1.9 (p = 0.042) for disease development (30).
For viral disease, results from chronic hepatitis C patients indicated that the rs9657182 was related with moderate or severe clinical depressive behavior during treatment with IFN-α, affecting its efficacy (41). Moreover, this SNP and G2431A (rs3739319) were analyzed in patients vaccinated or infected by hepatitis B virus (HBV) to assess HBV anti-HBs antibody responses. Both SNPs were not associated with antibody response activity. However, the TT genotype of rs9657182 was related to increase in mortality in HBV-infected patients (18).
It is observed that IDO1 is overexpressed in viral diseases due to the presence of ISREs (IFN-stimulated response elements) in the promoter region in addition to other sites that are activated by IFN-γ (34). Under gene ablation or inhibition of the enzyme by a drug during a given infection, the organism can become reactive by developing a memory T cell response, as well as having high expression of IFN-γ by CD4+ and CD8+ T cells, as occurs during influenza virus infection (35). In HIV-positive patients, plasma levels of IDO have already been used as a biomarker for tuberculosis (3).
It has also been reported in viral diseases that IDO is involved in growth inhibition of the herpes virus, cytomegalovirus and parainfluenza virus type 3 (HPIV3) (2,11,32), whereas respiratory syncytial virus induces the expression and activity of IDO, which is possibly related to the development of allergic disease (4,33). In Puumala hantavirus infection, IDO activity is associated with a severe course of Nephropathia epidemica (form of hemorrhagic fever with renal syndrome) due to the suppressive capability of Treg cells (26). Another study has shown that IDO inhibition generates a more robust T cell response to influenza virus and the authors suggest this as a potential approach to increase the cell-mediated immune response to vaccination (16).
Some studies have evaluated the activity of the IDO enzyme in viral infections caused by flaviviruses and have shown its importance in the immunological response. In Japanese Encephalitis (JE), it has been observed that decreased IDO1 expression is associated with suppressed viral replication, which causes an increase in resistance to the development of the disease (25,47). In West Nile virus fever, the decrease in L-Trp caused by increased IDO1 expression attenuated viral replication (29). Moreover, overexpression of hepatic IDO before hepatitis C virus (HCV) infection markedly impaired HCV replication in hepatocytes, suggesting that IDO limits the spread of HCV within the liver beyond IFN to be stimulated. Furthermore, it also regulates host immune responses through an inhibitory effect on CD4+ T cell proliferation (28).
Regarding dengue, one study observed a synergistic effect between DENV and IFN-γ on induction of the IDO1 gene and a partial dependence of the activity of this enzyme on the antiviral response to DENV caused by IFN-γ. In this study, an increase in the IDO1 gene expression in response to infection by DENV in human umbilical vein endothelial cells (HUVECs) was found in vitro. In the serum of dengue-infected subjects, a decrease in tryptophan levels and an increase in kynurenine levels was observed, indicating an increase in tryptophan catabolism by IDO activity (14). Another study carried out with DENV-4-infected individuals showed an increase in IDO in blood monocytes from dengue patients, especially in nonclassical monocytes with high CD14+ and CD16+ expression. A correlation of IDO with IL-10, an anti-inflammatory cytokine, was also found. Based on these findings, it has been proposed that this enzyme may be exerting both antiviral effects and modulating the exacerbated immune response during dengue (14). However, no difference was found between the levels of IDO in pairs with classical dengue and SD (43).
While our results deliver novel understanding, we recognize that sample size is a limitation of our study. There is now a need to perform larger studies and to study other SNPs that might influence this pathway, such as SNPs in cytokine genes, or molecules that are associated with this response, such as IFN-γ.
We are presenting the first significant evidence of an association between the IDO1 G2431A variant and the risk for SD. We believe that this SNP can be a potential target to other genetic and pathophysiological studies that examine the role of host genetic factors in human dengue disease.
Footnotes
Acknowledgments
The authors are grateful to Dr. Suzannah Lant for all the support for this work.
Author Disclosure Statement
No competing financial interests exist.
