Abstract

Genetic susceptibility factors, parasite strain, and an adequate modulation of the immune system seem to be crucial for disease progression after Trypanosoma cruzi infection. HLA-G and its murine functional homolog Qa2 have well-recognized immunomodulatory properties. We evaluated the HLA-G 3′ untranslated region (3′UTR) polymorphic sites (associated with mRNA stability and target for microRNA binding) and HLA-G tissue expression (heart, colon, and esophagus) in patients presenting Chagas disease, stratified according to the major clinical variants. Further, we investigated the transcriptional levels of Qa2 and other pro- and anti-inflammatory genes in affected mouse tissues during T. cruzi experimental acute and early chronic infection induced by the CL strain. Chagas disease patients exhibited differential HLA-G 3′UTR susceptibility allele/genotype/haplotype patterns, according to the major clinical variant (digestive/cardiac/mixed/indeterminate). HLA-G constitutive expression on cardiac muscle and colonic cells was decreased in Chagasic tissues; however, no difference was observed for Chagasic and non-Chagasic esophagus tissues. The transcriptional levels of Qa2 and other anti and proinflammatory (CTLA-4, PDCD1, IL-10, INF-γ, and NOS-2) genes were induced only during the acute T. cruzi infection in BALB/c and C57BL/6 mice. We present several lines of evidence indicating the role of immunomodulatory genes and molecules in human and experimental T. cruzi infection.

1. Introduction

Chagas disease or American trypanosomiasis is caused by the parasite Trypanosoma cruzi, which may produce acute and chronic manifestations. During the acute phase, most cases are indeterminate, but some patients may develop severe myocarditis or meningoencephalitis, which can be lethal. The chronic phase exhibits four major clinical variants: (i) the indeterminate form represents 50–70% of cases and develops without evident clinical and pathological signs [1]; (ii) the cardiac form encompasses 10–20% of cases and usually presents progressive congestive heart failure, various cardiac arrhythmias, thromboembolic events, and sudden death [1, 2]; (iii) the digestive form (8–10% of cases) is characterized by clinical signs of megaesophagus, megacolon, or both; (iv) the cardiodigestive or mixed form (8% of cases) comprises clinical and pathological signs of cardiac and digestive involvement. These clinical syndromes are considered to be a consequence of autonomic neuronal loss, microvascular derangements, and chronic organ inflammation directly dependent on parasite persistence or by immune system cells [3, 4].

In parallel with the chronic parasite infection, the immune response against parasite antigens is associated with the Chagas heart disease pathogenesis. Key features of the immunity in the Chagas disease include (i) predominance of partially activated T CD8 lymphocytes in cardiac inflammatory infiltrates [5], accompanied by high production of nitric oxide (NO), IL-12, monocyte chemoattractant protein-1 (MCP-1), and IFN-γ by infiltrating macrophages; (ii) preponderant IFN-γ production by cytolytic natural killer (NK) cells in the acute phase of the disease to control tissue and systemic parasite burden; and (iii) polyclonal B lymphocyte response [613].

Cytokine patterns in Chagas disease are characterized by T-helper (Th)1 polarization in the acute phase, with predominance of IFN-γ and TNF production. During the chronic phase, a coexistence of the Th1 (IL-12, IFN-γ, and TNF) and Th2 profiles (IL-4 and IL-10) may be observed, and the equilibrium between these profiles may be relevant to disease morbidity [3, 14, 15]; that is, a Th1 response aggravates and a Th2 response produces a better outcome in murine trypanosomiasis and human Chagas disease [1620]. Besides cytokines, other immunomodulatory molecules such as PD-1 and CTLA-4 are also involved in acute and chronic trypanosomiasis [2025]. In this context, immunomodulatory molecules may play an important role for disease evolution, controlling or exacerbating the immune response against T. cruzi itself or tissue modifications induced by parasites.

Among immunomodulatory molecules, the human leucocyte antigen- (HLA-) G is a nonclassical histocompatibility class Ib molecule, which has a well-recognized role in controlling several branches of the immune response, inhibiting T cell proliferation, NK, and cytotoxic T lymphocytes, inducing regulatory T cells and tolerizing dendritic cells [26]. These effects are primarily due to the preferential interaction of the HLA-G molecule with the Ig-like transcript (ILT)2, ILT4, and killer cell immunoglobulin-like (KIR2DL4) receptors, which induce inhibitory intracellular signals via tyrosine-based (ITIM) motifs [27].

HLA-G mRNA has restricted expression in nonpathological conditions, for organs such as the placenta, thymus, heart, intestines, brain, and skin [28], and the constitutive expression of the HLA-G molecule has been observed on placenta, thymus, cornea, pancreas, brain, erythroid cells, and blood cell surfaces [2934]. On the other hand, HLA-G neoexpression has been observed in several situations, including tumors, viral infections, engrafted tissues, and autoimmune and inflammatory diseases. Depending on the underlying condition, the expression may be advantageous when the blockade of the immune response is desirable, such as in autoimmune disorders and in allografting, whereas HLA-G expression in tumors and chronic viral infections may be harmful [35, 36].

At the coding region, the HLA-G gene presents limited polymorphisms compared to classical HLA class I (HLA-A/B/C) genes, exhibiting only 50 alleles (The International Immunogenetics Database (IMGT), v3.16.0) and coding 16 different membrane-bound full length proteins. The regulatory regions also present several polymorphisms that coincide with or are close to binding sites for transcription factors (promoter region) or are targets for microRNA binding (3′ untranslated region (3′UTR)). At least 29 variation sites have been described at the promoter region [37] and 16 variants have been identified at the HLA-G 3′UTR [38]. Considering that the regulatory regions are involved on the magnitude of HLA-G gene expression, the study of these regions is relevant. In the healthy Brazilian population, at least eight HLA-G 3′UTR variation sites have been described [37, 39], and some of them have been associated with plasma soluble HLA-G (sHLA-G) levels [40].

The HLA-G murine functional homolog Qa2, encoded by H2-Q7/Q9 gene, also presents restricted tissue expression, and Qa2 mRNA has been reported in thymus, liver, intestines, spleen, placenta, and brain. Qa2 also modulates the immune response by interacting with as yet unidentified NK cell receptors [4145].

Little information is available regarding HLA-G and the Qa2 murine functional homolog in acute and chronic parasitic diseases. Considering that (i) no information is available regarding the role of these immunoregulatory molecules in T. cruzi infection; (ii) Chagas disease is a chronic disorder, in which several mechanisms of immunomodulation have been described in association with its pathogenesis; (iii) clinical variants of Chagas disease may depend on parasite and host genetic factors; (iv) experimental trypanosomiasis may help in the understanding of the human disease counterpart, we designed a study encompassing the HLA-G tissue expression and HLA-G 3′UTR polymorphic site typing in patients presenting Chagas disease, stratified according to major clinical variants. In addition, we evaluated the transcriptional level of the H2-Q7/Q9 (Qa2), and other immunoregulatory H2-T23 (Qa1), CTLA-4, and PDCD1 (PD-1) genes in mouse (BALB/c and C57BL/6) affected tissues during experimental acute and early chronic infection caused by the CL strain of T. cruzi and correlated the expression of these genes with other mediators of inflammation, including INF-γ, inducible nitric oxide synthase (NOS-2), and IL-10.

2. Methods

2.1. Subjects

The protocol of the study was approved by the local Research Ethics Committee (Protocol number 11237/2009) and written informed consent was obtained from all participants. A total of 177 chronic Chagas disease patients exhibiting positive serology for T. cruzi antigens followed at the Divisions of Cardiology and Gastroenterology of the Department of Medicine of the School of Medicine of Ribeirão Preto, University of São Paulo, Brazil, and 155 healthy individuals from the same region of patients, exhibiting negative serology for T. cruzi infection, were studied. The Chagas disease patients were submitted to clinical and laboratory examination. Electro- and echocardiography, esophagus and colon contrast X-ray examinations, and esophagus electromanometry were performed to classify patients into cardiac (, exhibiting or not heart failure), digestive (, exhibiting megaesophagus, megacolon, or both), cardiodigestive or mixed (, presenting any combination of cardiac and digestive forms), and indeterminate () variants.

2.2. DNA Extraction and Genotyping

Genomic DNA was extracted from peripheral blood leucocytes using a standard salting-out procedure. Genotyping of the variation sites at HLA-G 3′UTR was performed by sequencing analyses as previously described [39].

2.3. Immunohistochemical Analysis of HLA-G

Fifty-four specimens of tissues obtained from Chagas disease deceased patients were analyzed, of which 30 were derived from heart, 13 from colon, and 11 from esophagus. In parallel, 20 tissue specimens (8 from heart, 6 from esophagus, and 6 from sigmoid colon) of deceased individuals without Chagas disease and exhibiting no histopathological tissue alteration were analyzed. Three specimens of trophoblast tissue (positive control for HLA-G expression) were also analyzed. The immunohistochemical analysis was performed using the primary antibody MEM/G9 (HLA-G1 and sHLA-G5 specific mouse IgG) (Exbio, Prague, Czech Republic), diluted 1 : 200, as previously described [46]. Briefly, after defining the range of brown color considered being positive, images were converted to 256 shades (8-bit) of gray. Then, the grayscale images were converted into a binary (black and white) variable to define the cutoff point. The threshold was adjusted, and the brown areas became black portions in the binary image [47]. Immunostained areas were evaluated using the ImageJ software (National Institutes of Health, Bethesda, MD).

2.4. Experimental Infections and Transcriptional Level Analysis

Male eight-week-old C57BL/6 and BALB/c mice were intraperitoneously (IP) injected with bloodstream CL strain forms of T. cruzi, using three infected and three noninfected animals in two independent experiments. The mice were cared for according to the institutional guidelines on ethics in animal experiments and all protocols were approved by the local Ethics Committee on Animal Care and Research (process number 172/2009).

For the acute infection, 103 bloodstream forms were IP injected. Mouse survival was verified daily for 30 days and the parasitemia was quantified microscopically by counting the parasites in 5 μL of citrated blood obtained from the tail lateral vein from day 7 until day 29 after infection. Parasite load was analyzed in the heart at day 24 after infection [48]. For the early chronic infection, 102 bloodstream forms of T. cruzi were IP injected and mouse survival was verified daily for 60 days. Parasite load was analyzed in the heart at day 60 after infection.

Gene expression was detected in the heart and esophagus at 24 and 60 days after infection (dpi), for acute and early chronic infections, respectively. Although Qa2, Qa1, and PD-1 are the encoded molecules of the H2-Q7/Q9, H2-T23, and PDCD1 genes, respectively, to facilitate understanding we will refer in this text to Qa2, Qa1, and PD-1 expression as synonyms of their respective gene expressions.

Total RNA was extracted from tissues using the TRIzol reagent (Invitrogen, Carlsbad, CA) and was reverse transcribed to obtain cDNA with the SuperScript III (Invitrogen) transcriptase reverse. SYBR Green Mix-based real-time quantitative PCR assays were performed using a StepOnePlus system (Applied Biosystems, Foster City, CA). Data were normalized according to the expression of the glyceraldehyde phosphate dehydrogenase (GADPH) housekeeping gene and relative expression of each mRNA was calculated with the method [49]. Primer sequences included (i) H2-Q7/Q9-forward 5′-ATG GCG ACC ATT GCT GTT GT-3′, H2-Q7/Q9-reverse 5′-TCC AAT GAT GGC CAC AGC T-3′; (ii) H2-T23-fwd 5′-GCA CAA GTC AGA GGC AGT CG-3′, H2-T23-rev 5′-TGC AGG TAT GCC CTC TGT TG-3′; (iii) CTLA-4-fwd 5′-ACC TCT GCA AGG TGG AAC TCA-3′, CTLA-4-rev 5′-CCA TGC CCA CAA AGT ATG GC-3′; (iv) PD-1-fwd 5′-TTC AGG TTT ACC ACA AGC TGG-3′, PD-1-rev 5′-TGA CAA TAG GAA ACC GGG AA-3′; (v) INF-γ-fwd 5′-GCA TCT TGG CTT TGC AGC T-3′, INF-γ-rev 5′-CCT TTT TCG CCT TGC TGT TG-3′; (vi) NOS-2-fwd 5′-CGA AAC GCT TCA CTT CCA A-3′, NOS-2-rev 5′-TGA GCC TAT ATT GCT GTG GCT-3′; (vii) IL-10-fwd 5′-TGG ACA ACA TAC TGC TAA CC-3′, IL-10-rev 5′-GGA TCA TTT CCG ATA AGG CT-3′; (viii) GAPDH-fwd 5′-TGC AGT GGC AAA GTG GAG AT-3′, GADPH-rev 5′-CGT GAG TGG AGT CAT ACT GGA A-3′.

2.5. Statistical Analysis

Allele and genotype frequencies were estimated by direct counting, and adherences of phenotypical proportions to expectations under the Hardy-Weinberg equilibrium (HWE) were tested by the complete enumeration method using the GENEPOP 3.4 software [50]. Linkage disequilibrium (LD) between variation sites at HLA-G 3′UTR was evaluated by means of a likelihood ratio test of linkage disequilibrium implemented at the ARLEQUIN software [51]. HLA-G 3′UTR haplotypes of each individual were inferred for each group of patients (stratified according to the four clinical variants) and for the whole group of patients. Two distinct computational methods that do not take any prior information into account were used for this purpose: (1) the EM algorithm [52] implemented with the PL-EM software [53] and (2) a coalescence-based method implemented in the PHASE v2 software [54]. Therefore, the haplotype pair of each subject was independently inferred by four approaches; whenever these four approaches resulted in nonconsensual inference for a given subject, he was removed from all procedures that used haplotype data as input. The frequency of each allele, genotype, or haplotype was compared between patients and controls by the two-sided Fisher exact test, with the aid of the GraphPad Instat 3.05 software, which was also used to estimate the Odds Ratio (OR) and its 95% Confidence Interval (CI).

For the HLA-G tissue expression, gene expression, and parasitism load, statistical analyses were performed using Student’s -test or Mann-Whitney tests and the relationship between immunomodulatory and mediators of inflammation genes was performed using the Spearman Correlation, both with the aid of the GraphPad Prism 5 v5.0b software. To facilitate data presentation, variables were expressed as mean ± SEM.

For all analyses, significance was defined at .

3. Results

3.1. Genotype and Allele Frequencies of Variation Sites at HLA-G 3′UTR

HLA-G 3′UTR variation sites were compared between (i) patients (considered as whole) with controls (healthy individuals); (ii) patients presenting clinically detectable disease (CDM) (cardiac (C) plus digestive (D) plus mixed (M) forms) with controls; (iii) patients presenting CDM with indeterminate patients (indeterminate form—I); and (iv) patients stratified according to the four clinical variants with each other and with controls. Association tests were performed for alleles and genotypes of each variation site.

The following HLA-G 3′UTR variation sites were observed in this study: 14-base pair insertion/deletion (I/D; rs66554220); +3001C/T (rs not available); +3003T/C (rs115689421); +3010C/G (rs116152775); +3027A/C (rs115810666); +3035C/T (rs115100128); +3142G/C (rs115928989); +3187A/G (rs114317070); and +3196C/G (rs115045214). Overall, these variation sites adhered to the Hardy-Weinberg Equilibrium, except +3010C/G in cardiac and whole group of patients and +3142G/C in digestive and whole group of patients. Table 1 shows the genotype and allele frequencies in patients with chronic Chagas disease and controls.

14-base pair I/D and +3001C/T: no significant differences were observed when the five types of comparisons were performed. The +3001T allele was observed in only one Chagas disease patient presenting exclusively the digestive form.

+3003T/C: the +3003CC genotype was not observed in Chagasic patients. Compared to controls, (i) the +3003TC genotype frequency was decreased in the whole group of patients () and in the digestive form (); (ii) the frequency of the +3003C allele was decreased in the whole group of patients () and in patients exhibiting the digestive form (); (iii) +3003TT genotype frequency was overrepresented in the whole group () and in digestive () patients; and (iv) +3003T allele frequency was increased in whole group () and in the digestive group () of patients.

+3010C/G: compared to controls, the +3010CG genotype was underrepresented in the whole group of patients () and in those presenting the cardiac variant ().

+3027A/C: a decreased frequency of the +3027AC genotype was observed in the digestive form compared to control () and to patients with the indeterminate (), cardiac (), and mixed () forms. The frequency of the +3027CC genotype was increased in patients exhibiting the digestive form compared to indeterminate (), cardiac (), mixed (), and control () groups. The +3027AA and +3027AC genotypes were not observed in patients with the digestive form. As consequence, the digestive group exhibited a decreased frequency of the +3027A allele compared to controls () and to patients with the indeterminate (), cardiac (), and mixed () forms. The +3027C allele presented the opposite association pattern.

+3035C/T: the +3035CC genotype was overrepresented in patients with the digestive form compared to control () and to patients with cardiac () and mixed () forms. A decreased frequency of the +3035CT genotype was observed in patients exhibiting the digestive form compared to mixed () and control () groups. The +3035C allele was overrepresented in the group of patients with the digestive form compared to control (), indeterminate (), cardiac (), and mixed () groups. The +3035T allele presented the opposite association pattern.

+3142G/C: the +3142GG genotype frequency was significantly lower in digestive patients compared to indeterminate (), cardiac (), and mixed () forms. The +3142GC genotype frequency was lower in cardiac patients compared to digestive () and control () groups. In addition, +3142GC genotype frequency was decreased in the whole group of patients group compared to the control group ().

+3187A/G: compared to controls, an increased frequency of the +3187GG genotype was observed in the whole group () and in the group of patients with the cardiac form ().

+3196C/G: patients exhibiting the mixed form presented a decreased frequency of the +3196CC genotype () and an increased frequency of +3196GC genotype (). In addition, +3196CG genotype was increased in the mixed group compared to the indeterminate variant ().

When CDM patients (excluding the indeterminate group) were compared to controls, the following results were observed: (i) the +3003CT genotype frequency was decreased in patients (); (ii) the frequency of the +3003TT genotype was increased in patients (); (iii) +3003C allele frequency was decreased in patients () and +3003T allele presented the opposite association pattern (); (iv) the +3187GG genotype frequency was overrepresented in patients (); (v) +3196CC genotype frequency was decreased () and +3196GC genotype frequency was increased () in patients.

When CDM patients were compared to the indeterminate form, only the +3196GC genotype was increased in CDM patients (). On the other hand, when indeterminate patients were compared to healthy controls, no significant results were observed.

Table 2 shows the Odds Ratio and 95% Confidence Interval values obtained for all significant comparisons.

3.2. Haplotype Frequency of Variation Sites at HLA-G 3′UTR

To further understand how the ensemble of variation sites participate in Chagas disease susceptibility and considering that these variation sites were in linkage disequilibrium (results not shown), we analyzed the frequency of the HLA-G 3′UTR haplotypes observed for patients and controls, and the results of the statistical analyses are shown in Tables 3(a) and 3(b). Eleven out of 177 patients have been removed from all analyses that used haplotype data as input, since results for the four approaches of haplotype inference resulted in nonconsensual results. The names of the 3′UTR haplotypes (UTR-1, UTR-2, etc.) follow the nomenclature previously described by our group [39], which has been extensively used in the literature.

Compared to controls, (i) UTR-4 was underrepresented in the whole group of patients (), CDM (), digestive (), and cardiac () forms; (ii) UTR-13 was overrepresented in the indeterminate group (); (iii) UTR-14 was overrepresented in cardiac patients (). A decreased frequency of the UTR-7 haplotype was observed in the digestive group compared to indeterminate (), mixed (), and control () groups. In addition, UTR-13 was also increased in patients presenting the indeterminate form when compared to CDM patients ().

3.3. HLA-G Expression in Specimens Obtained from Chronic Chagas Disease Patients

HLA-G molecule expression was evaluated in the major organs affected by the disease, including heart, colon, and esophagus. Compared to non-Chagasic tissues, HLA-G expression was significantly decreased in Chagasic heart (; Figures 1(a), 1(b), and 1(c)) and colon (, Figures 1(d), 1(e), and 1(f)). HLA-G expression in individuals without Chagas disease was primarily observed on cardiac muscle cells and no cellular infiltration was observed (Figure 1(a)), whereas specimens from Chagas patients exhibiting cardiomegaly showed lesser HLA-G expression on cardiac muscle cells together with an infiltration of mononuclear cells (lymphocytes and plasma cells) exhibiting HLA-G expression (Figure 1(b)). The HLA-G immunolabeling in esophagus of Chagas patients with esophagomegaly was closely similar to that observed for individuals without Chagas disease (Figures 1(g), 1(h), and 1(i)).

3.4. Qa2 Expression during Experimental Acute and Early Chronic Infections

The transcriptional level of Qa2 in acute and early chronic infections was studied in heart and esophagus specimens obtained from BALB/c and C57BL/6 infected mice, using the T. cruzi CL strain.

First, we characterized the acute and early chronic infections regarding the animal survival rate, tissue parasite load, and blood parasitism. Survival rates showed no significant differences for BALB/c and C57BL/6 mice for both acute (Figure 2(a)) and early chronic (Figure 2(b)) infections. During acute infection, the heart parasite load was 2-fold increased in C57BL/6 mice in relation to BALB/c group (; Figure 2(c)). No significant difference was observed between the two groups during the early chronic infection (Figure 2(d)). In addition, when we analyzed the bloodstream forms in the acute infection, BALB/c mice showed an increased parasitism level at days 21, 23, and 25 compared to C57BL/6 mice (; Figure 2(e)).

During acute infection, the heart and esophagus transcriptional levels of Qa2 were significantly increased for both mouse strains, when compared to noninfected mice. Compared to the control group, the heart expression of Qa2 in BALB/c and C57BL/6 was 28-fold and 25-fold increased, respectively (; Figure 3(a)), and the esophagus expression of Qa2 was 17-fold and 16-fold increased in BALB/c and C57BL/6, respectively (; Figure 3(b)).

Acutely infected BALB/c mice showed increased transcriptional levels of Qa2 in heart and esophagus, which were 7-fold and 14-fold higher than early chronically infected mice, respectively (; Figures 3(a) and 3(b)). The transcriptional level of Qa2 was 6-fold higher in heart of acutely infected C57BL/6 mice compared to early chronic infection (; Figure 3(a)). No significant difference in terms of Qa2 expression was observed between noninfected and early chronically infected mice in heart and esophagus of BALB/c and C57BL/6 mice (Figures 3(a) and 3(b)).

To verify whether other known immunomodulatory genes were concomitantly modulated during the infection by CL strain of T. cruzi, we analyzed the transcriptional levels of Qa1, CTLA-4 and PD-1 genes. Our results showed that all of them were significantly augmented in heart and esophagus of BALB/c strain when compared to noninfected and early chronically infected mice (; Figures 3(a) and 3(b)). C57BL/6 mice showed an increased transcriptional level of Qa1 in heart and increased levels of CTLA-4 in esophagus and PD-1 in both heart and esophagus, during the acute infection (), compared to controls and early chronically infected animals (Figures 3(a) and 3(b)). Moreover, during the acute infection, Qa1 was 4-fold higher in esophagus of BALB/c strain compared to C57BL/6 animals, and PD-1 expression was 2-fold higher in heart of C57BL/6 mice compared to BALB/c strain (; Figures 3(a) and 3(b)). No significant difference was observed between noninfected and early chronically infected mice in heart and esophagus specimens obtained from BALB/c and C57BL/6 mice (Figures 3(a) and 3(b)).

Considering that nonclassical histocompatibility genes may be influenced by the action of inflammatory mediators, we evaluated the transcriptional levels of INF-γ, NOS-2, and IL-10 genes in the heart to observe the relationship between these genes and the nonclassical genes.

INF-γ, NOS-2, and IL-10 expression was increased in both BALB/c and C57BL/6 strains during the acute infection compared to controls and early chronically infected animals (; Figure 3(c)). The comparisons between strains showed that INF-γ expression was 4-fold higher in C57BL/6 than in BALB/c (Figure 3(c)). During the early chronic infection, the transcriptional levels of these genes were similar to noninfected mice for both BALB/c and C57BL/6 strains (Figure 3(c)). When we studied the relationship between these genes, we observed a positive correlation between Qa2 and INF-γ (; ) and Qa2 and NOS-2 (; ) in heart specimens of both mouse strains during acute T. cruzi infection.

4. Discussion

The present study provides consistent data regarding the involvement of HLA-G and its Qa2 murine functional homolog in the T. cruzi infection. HLA-G gene polymorphism and HLA-G tissue expression have been evaluated in several diseases; however, little is known about the role of these molecules in the human or experimental T. cruzi infections.

The association of genetic markers with infectious disease has posed intriguing questions regarding the comparisons of gene frequencies between individuals with the disease with healthy individuals who had never been exposed to the specific infectious agent. Due to the myriad of chronic clinical forms observed in Chagas disease, part of these concerns may be circumvented by the possibility of comparing patients presenting well recognized clinical variants with patients who were infected by T. cruzi and have not developed clinically detected forms (indeterminate). When patients as a whole were compared to healthy individuals, some HLA-G 3′UTR alleles/genotypes were overrepresented (+3003T allele and +3003TT and +3187GG genotypes), while others were underrepresented (+3003C allele, and +3003CT, +3010GC, +3042GC genotypes) in Chagas disease patients. When symptomatic patients (CDM) were compared to healthy controls, similar results were obtained, indicating that these polymorphic sites were indeed associated with the clinically detected Chagas disease. In contrast, when indeterminate patients were compared to controls, the frequency of the HLA-G 3′UTR alleles/genotypes was closely similar. This result corroborates the role of the aforementioned HLA-G 3′UTR alleles and genotypes in the clinically recognized Chagas variants.

The stratification of Chagas patients into clinical variants showed that the +3027CC and +3035CC genotypes and the +3027C and +3035C alleles were strongly associated with the digestive form of Chagas disease. In addition, the +3142GC genotype was associated with a decreased risk to cardiac form development, when compared to the digestive form. These results indicate that genetic susceptibility to digestive or cardiac forms may be distinct. Indeed, previous studies conducted by our group evaluating classical histocompatibility (HLA-A, -B, -DRB1 and DQB1) antigens/alleles [55] and other immunoregulatory genes, such as CTLA-4 [25], also showed a differential pattern of susceptibility according to the Chagas clinical variant. No specific posttranscription regulation mechanism has been reported concerning the +3027A/C and +3035C/T polymorphic sites; however, in a previous in silico study, we reported several miRNAs that may target the +3027A/C and +3035C/T polymorphic sites [56].

When HLA-G 3′UTR polymorphic sites were considered as haplotypes, decreased UTR-4 and UTR-7 frequencies were associated with the clinically detectable Chagas disease and with the digestive form, respectively. On the other hand, UTR-14 was associated with the cardiac variant and UTR-13 with the indeterminate form. Noteworthy, particular posttranscription elements that are present in these UTRs also exhibited a similar differential behavior when their frequencies were individually evaluated, further corroborating the role of HLA-G 3′UTR sites in disease susceptibility. Our group demonstrated that specific HLA-G 3′UTR haplotypes are associated with differential soluble levels of HLA-G (sHLA-G), being UTR-1 associated with high sHLA-G and UTR-5 and UTR-7 with decreased sHLA-G levels [40]. However, there are no studies associating HLA-G 3′UTR haplotypes with tissue HLA-G expression.

Considering that cytokines, chemokine, and immunomodulatory molecules play an important role in the pathogenesis of many diseases, including Chagas disease, and gene polymorphisms may influence the expression of these molecules, different studies have reported the association between different clinical forms of Chagas disease and polymorphisms at TNF [5761], IL-1β [62], IL-10 [63], IL-4 [64], INF-γ [65], TGF-β1 [66], IL-12B [67], CXCL9 [68], CXCL10 [68], CCR5 [6870], CCL2 [71], and CTLA-4 [25] genes.

Regarding HLA-G expression on tissue cells, patients with the digestive form exhibited decreased expression in colon but not in esophagus specimens, contrasting with patients exhibiting the cardiac variant who presented a decreased expression on myocardium cells. Many segments of the digestive tract and normal myocardium fibers [28] may constitutively express HLA-G; however, the meaning of these findings has not been completely elucidated. The best characterized constitutive HLA-G expression is in the placenta, where it may provide protection to the fetus against the cytotoxic cell immune response against paternal antigens [72]. HLA-G constitutive expression in gut and heart may in fact present a similar protective effect. Indeed, increased plasma levels of sHLA-G and myocardium HLA-G expression have been associated with decreased cellular acute rejection and better graft survival in patients submitted to heart transplantation [7375]. Then, the loss of HLA-G expression on myocardium cells may facilitate the action of infiltrating lymphomononuclear cells after tissue injury caused by T. cruzi. On the other hand, myocardium infiltrating lymphomononuclear cells also exhibited HLA-G expression, which could decrease their functional properties. Since myocardium damage has been primarily attributed to fiber losses due to the action of immune cytotoxic cells, one may consider that HLA-G lymphomononuclear cell expression may not be the only regulatory surface cell molecule. In addition, it is important to emphasize that all these findings were observed in heart specimens of deceased patients, and the temporal HLA-G expression during the chronic phase has not been evaluated as yet.

Overall, the outcome of the experimental T. cruzi infection may depend on several factors, including (i) the inoculated strain, (ii) the amount of parasites, and (iii) the genetic background of the animal. Several of these combinations have been used to induce acute or early chronic infection. BALB/c mouse has been considered to be susceptible to acute T. cruzi infection induced by the Y and Tulahuén strains, due to the predominance of a Th2 immune response, characterized by high production of IL-4 and IL-10 [3, 76]. On the other hand, C57BL/6 mouse is considered to be resistant to acute T. cruzi infection induced by Y and Tulahuén strains, since the mouse produces a Th1 immune response with the production of INF-γ, IL-12, TNF, and NO [3, 14, 15]. In the present study, we used the less virulent T. cruzi CL strain that depending on the dose may produce severe acute or lead to the development of early chronic infection. In BALB/c mice, the acute infection by CL strain induced higher heart parasitism and parasitemia with increased production of proinflammatory mediators [77]. Since BALB/c is more susceptible to T. cruzi infection, chronic infection is difficult to be induced by more virulent strains. On the other hand, little is known about the acute and early chronic infection induced by CL strain in C57BL/6 mice. With the use of the CL strain we did not observe a differential susceptibility in BALB/c and C57BL/6 mice in terms of mouse survival; however, BALB/c showed higher parasitemia, whereas C57BL/6 presented higher heart parasitism. Therefore, the C57BL/6 controlled CL strain parasitemia but not heart tropism, and the inverse occurred with the BALB/c mouse.

Similarly to HLA-G, the Qa2 molecule may also control the immune response by inhibiting NK cells [45]. During acute T. cruzi infection, the transcriptional level of Qa2 was increased in heart and esophagus in both BALB/c and C57BL/6 mouse strains, and the increased expression was not related to resistance or susceptibility to the acute infection. Although little attention has been devoted to the role of immunomodulatory nonclassical MHC molecules in experimental infection, Qa2 expression may be induced to counterbalance the action of proinflammatory and other anti-inflammatory or immunomodulatory molecules. Indeed, the transcriptional levels of the Qa1, CTLA-4, PDCD1, INF-γ, NOS-2, and IL-10 genes were also induced in the acute T. cruzi infection, as shown in this study. Additional literature findings corroborate this idea. Tissue damage observed in the affected tissues is due to the inflammatory process generated by an exacerbated immune response triggered against the parasite, and the high production of CTLA-4, PD-1, and IL-10 is related to control inflammatory process during T. cruzi infection [3, 14, 15, 23, 24, 78]. Our data suggest that the augmentation of immunomodulatory gene expression during the acute phase, including Qa2 gene, may be associated with the control of tissue damage caused by inflammatory process, cytolysis, and fibrosis.

Regarding the early chronic infection, the expression of immunomodulatory, proinflammatory, and anti-inflammatory genes analyzed in this study was closely similar to noninfected mice in both BALB/c and C57BL/6 strains, and both animal groups were able to evolve into the early chronic infection, presenting only a little amount of T. cruzi in the heart. Considering that mice were exposed to decreased number of the less virulent CL strain, it is reasonable to conclude that both BALB/c and C57BL/6 mice had time and chance to mount an adequate immune response against the parasites.

5. Conclusion

We present several lines of evidence pointing to the participation of the immunomodulatory molecules HLA-G and its functional murine homolog Qa2 on human and experimental T. cruzi infection. In terms of Chagas disease susceptibility, HLA-G 3′UTR alleles/genotypes/haplotypes exhibited differential frequencies in infected/diseased patients when compared to only infected patients and healthy controls. The observation of a decreased HLA-G expression on cardiac muscle and colonic cells associated with the increased expression of HLA-G on myocardium infiltrating lymphomononuclear cells may reflect a lack of protection of these tissues that constitutively express HLA-G. In experimental infection, the mouse genetic background exerted only a mild influence on the overall response against the T. cruzi CL strain. The transcriptional levels of Qa2, Qa1, CTLA-4, PDCD1, INF-γ, NOS-2, and IL-10 genes were induced only during the acute T. cruzi infection in BALB/c and C57BL/6 mice, indicating a fine balance between pro- and anti-inflammatory genes.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Ana Maria da Rocha for technical assistance. This work was supported by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Grant nos. CAPES/COFECUB 653/09, PRODOC 288/05-5, and PNPD 02883/09-0), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq Science Without Borders Program, Grant no. 236754/2012-2), and Núcleo de Apoio a Pesquisa em Doenças Inflamatórias (NAP-DIN). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the paper.