key: cord-0002342-lx8lowf7
authors: Lupiañez, Carmen B.; Villaescusa, María T.; Carvalho, Agostinho; Springer, Jan; Lackner, Michaela; Sánchez-Maldonado, José M.; Canet, Luz M.; Cunha, Cristina; Segura-Catena, Juana; Alcazar-Fuoli, Laura; Solano, Carlos; Fianchi, Luana; Pagano, Livio; Potenza, Leonardo; Aguado, José M.; Luppi, Mario; Cuenca-Estrella, Manuel; Lass-Flörl, Cornelia; Einsele, Hermann; Vázquez, Lourdes; Ríos-Tamayo, Rafael; Loeffler, Jurgen; Jurado, Manuel; Sainz, Juan
title: Common Genetic Polymorphisms within NFκB-Related Genes and the Risk of Developing Invasive Aspergillosis
date: 2016-08-12
journal: Front Microbiol
DOI: 10.3389/fmicb.2016.01243
sha: db3824a1c673ce4d1bb41c065065598089c89578
doc_id: 2342
cord_uid: lx8lowf7

Invasive Aspergillosis (IA) is an opportunistic infection caused by Aspergillus, a ubiquitously present airborne pathogenic mold. A growing number of studies suggest a major host genetic component in disease susceptibility. Here, we evaluated whether 14 single-nucleotide polymorphisms within NFκB1, NFκB2, RelA, RelB, Rel, and IRF4 genes influence the risk of IA in a population of 834 high-risk patients (157 IA and 677 non-IA) recruited through a collaborative effort involving the aspBIOmics consortium and four European clinical institutions. No significant overall associations between selected SNPs and the risk of IA were found in this large cohort. Although a hematopoietic stem cell transplantation (HSCT)-stratified analysis revealed that carriers of the IRF4(rs12203592T/T) genotype had a six-fold increased risk of developing the infection when compared with those carrying the C allele (OR(REC) = 6.24, 95%CI 1.25–31.2, P = 0.026), the association of this variant with IA risk did not reach significance at experiment-wide significant threshold. In addition, we found an association of the IRF4(AATC) and IRF4(GGTC) haplotypes (not including the IRF4(rs12203592T) risk allele) with a decreased risk of IA but the magnitude of the association was similar to the one observed in the single-SNP analysis, which indicated that the haplotypic effect on IA risk was likely due to the IRF4(rs12203592) SNP. Finally, no evidence of significant interactions among the genetic markers tested and the risk of IA was found. These results suggest that the SNPs on the studied genes do not have a clinically relevant impact on the risk of developing IA.

Invasive Aspergillosis (IA) is an opportunistic infection caused by Aspergillus, a ubiquitously present airborne pathogenic mold. A growing number of studies suggest a major host genetic component in disease susceptibility. Here, we evaluated whether 14 single-nucleotide polymorphisms within NFκB1, NFκB2, RelA, RelB, Rel, and IRF4 genes influence the risk of IA in a population of 834 high-risk patients (157 IA and 677 non-IA) recruited through a collaborative effort involving the aspBIOmics consortium and four European clinical institutions. No significant overall associations between selected SNPs and the risk of IA were found in this large cohort. Although a hematopoietic stem cell transplantation (HSCT)-stratified analysis revealed that carriers of the IRF4 rs12203592T/T genotype had a six-fold increased risk of developing the infection when compared with those carrying the C allele (OR REC = 6.24, 95%CI 1.25-31.2, P = 0.026), the association of this variant with IA risk did not reach significance at experiment-wide significant threshold. In addition, we found an association of the IRF4 AATC and IRF4 GGTC haplotypes (not including the IRF4 rs12203592T risk allele) with a decreased risk of IA but the magnitude of the association was similar to the one observed in the single-SNP analysis, which indicated that the haplotypic effect on IA risk was likely due to the IRF4 rs12203592 SNP. Finally, no evidence of significant interactions among the genetic markers tested and the risk of IA was found. These results suggest that the SNPs on the studied genes do not have a clinically relevant impact on the risk of developing IA.

Invasive Aspergillosis (IA) is an opportunistic infection often caused by species of Aspergillus, a common saprophytic filamentous fungus that is ubiquitously present in the environment. The LIFE initiative has estimated that around 30 million patients worldwide are at high risk of developing IA every year and over 200,000 patients develop the infection annually. Despite the substantial improvement in efficacy of newly developed anti-fungal drugs (mainly azoles alone or in combination with other antifungal drugs), IA continues to be a life-threatening infection in an increasing proportion of immunocompromised or critically ill subjects (Neofytos et al., 2009; Steinbach et al., 2012) . IA is frequently encountered in patients of intensive care (Meersseman et al., 2007) and solid organ transplantation units (Singh et al., 2013) but also among those patients who undergo allogeneic hematopoietic stem cell transplantation (allo-HSCT) or are diagnosed either with acute myeloid leukemia (AML) or acute lymphoid leukemia (ALL) and receive high-dose chemotherapy regimens (Steinbach et al., 2012; Neofytos et al., 2013a,b) .

Although major clinical risk factors for IA have been identified (Kousha et al., 2011) and the management of high-risk patients has been improved through the optimization of prevention strategies and early initiation of anti-fungal prophylaxis, mortality rates for IA remain still today unacceptably high (up to 30-60%; Neofytos et al., 2009; Steinbach et al., 2012; Karthaus and Buchheidt, 2013) . These observational findings suggest that additional factors may contribute to the risk of developing IA. In this regard, a growing number of studies have suggested that host genetic polymorphisms within or near immune-related genes may contribute to determine the risk of developing the infection (Kesh et al., 2005; Sainz et al., 2007a Sainz et al., ,b, 2008a Sainz et al., ,b, 2010 Sainz et al., , 2012 Bochud et al., 2008; Mezger et al., 2008; Ramaprakash et al., 2009; Cunha et al., 2010 Cunha et al., , 2011 Cunha et al., , 2014 Chai et al., 2011; Grube et al., 2013; Stappers et al., 2014; Wojtowicz et al., 2015a,b) . Interestingly, a substantial proportion of the susceptibility markers identified to date for IA are located in genes directly or indirectly implicated in the activation of the nuclear factor-kappa B (NFκB) signaling pathway, which suggests a relevant role of this biological route in determining the risk of developing this fungal infection. In particular, it has been reported that single nucleotide polymorphisms (SNPs) within toll-like receptors (TLRs; Mambula et al., 2002; Kesh et al., 2005; Bochud et al., 2008; Pamer, 2008; Ramaprakash et al., 2009; Carvalho et al., 2012; Grube et al., 2013) , C-type lectins (Cunha et al., 2010; Chai et al., 2011; Sainz et al., 2012) , PTX3 (Cunha et al., 2014) , and tumor necrosis factor receptors (TNFRs; Sainz et al., 2007b Sainz et al., , 2010 , which are pathogen recognition receptors (PRRs) that often culminate in the activation of NFκB pathway, may render patients more susceptible to develop IA.

Based on these findings but also those that have demonstrated that NFκB1 (p105/p50), NFκB2 (p100/p52), Rel A (p65), RelB, c-Rel genes may form homo-and hetero-dimers to regulate the activation of the canonical and non-canonical NFκB pathways (Moynagh, 2005; Gilmore, 2006; Hoffmann et al., 2006; Schlitzer et al., 2013; Bajaña et al., 2016) but also IRF4-dependent immune responses (Boddicker et al., 2015) , we hypothesized that the presence of common genetic polymorphisms within NFκB1, NFκB2, RelA, RelB, Rel, and IRF4 genes might influence the risk of developing IA in high-risk patients. Thus, the aim of the present study was to investigate the relationship between 14 single nucleotide polymorphisms within these genes and the risk of IA but also to determine whether these variants could interact with each other to modify the risk of developing the infection.

Eight hundred and thirty-four high-risk European Caucasian patients undergoing allo-HSCT or being diagnosed with acute myeloid leukemia (AML) or acute lymphoid leukemia (ALL) and receiving intensive remission-induction chemotherapy were recruited in this case-control population-based study (Lupiañez et al., 2015) . Three hundred and thirty-five patients were ascertained from the aspBIOmics consortium (http://www. aspbiomics.eu) whereas 341 patients were collected from two Spanish medical institutions (University Hospital of Salamanca and Clinic University Hospital of Valencia) and through a Spanish multicenter clinical trial (PCRAGA, EU clinical trial number: 2010-019406-17; Aguado et al., 2015) . In addition, 148 patients were recruited from two Italian medical institutions (Università Cattolica del S. Cuore, Rome; and University of Modena and Reggio Emilia, AOU Policlinico, Modena). Of those 834 patients, a total of 157 patients were diagnosed with proven or probable IA according to the updated EORTC/MSG criteria (De Pauw et al., 2008) whereas the remaining 677 patients did not show any sign of fungal infection.

Fourteen polymorphisms within the NFκB1, NFκB2, RelA, RelB, Rel, and IRF4 genes were selected to be genotyped in the whole population ( Table 1) . SNP selection was based on previously reported associations with cancer (Curtin et al., 2010; Do et al., 2010; Han et al., 2011; Slattery et al., 2011; Seufert et al., 2013; Wang et al., 2014) and immune-related diseases including infections (Chen et al., 2006 (Chen et al., , 2016 Gregersen et al., 2009; Trynka et al., 2009; Eyre et al., 2010; Varadé et al., 2011; Bowes et al., 2012; Ellinghaus et al., 2012; Ali et al., 2013; Leung et al., 2014; Pan et al., 2015) but also because their potential functionality according to the Haploreg (http://www.broadinstitute.org/ mammals/haploreg/haploreg.php) and ENCODE annotation data (https://genome.ucsc.edu/ENCODE/). The genotyping of the selected polymorphisms was carried out at GENYO (Center for Genomics and Oncological Research: Pfizer/University of Granada/Andalusian Regional Government, Granada, Spain) using KASPar R assays (LGC Genomics, Hoddesdon, UK) according to manufacturer's instructions. For internal quality control, 5% of samples were randomly selected and included as duplicates. Concordance between the original and the duplicate samples for the 14 SNPs was ≥99.0%. Call rates for all SNPs were ≥90.0% with the exception of the IRF4 rs872071 SNP that was excluded from further analysis.

The Hardy-Weinberg Equilibrium (HWE) test was performed in the control group (non-IA patients) by a standard observedexpected chi-square (χ 2 ). Logistic regression analysis adjusted for age, sex, country of origin, allo-HSCT, underlying disease and prophylaxis status was used to assess the main effects of the selected SNPs on IA risk. We also performed gene-HSCT interaction analyses to determine whether the association between SNPs and IA was of similar magnitude in HSCT (at highest risk) and non-HSCT patients. Although the selection of variables for adjustment was based on well-established risk factors for IA, the partial availability of data regarding prophylaxis status, type of chemotherapy, immunosuppressive drugs (HSCT), HLA mismatch (HSCT), or CMV status did not allow us to assess the impact of these factors on the genetic associations tested. Statistical power of the overall and HSCTstratified analyses was estimated using Quanto software (http:// hydra.usc.edu/gxe/). All tests were conducted using the statistical software SPSS (v.20) and STATA (v.12) for MAC. In order to account for multiple comparisons, we calculated an adjusted significance level using the Meff method (Nyholt, 2004) but also considering the number of inheritance models tested (codominant, dominant, recessive, and log-additive). Thus, the significant threshold used for the main effect analysis was 0.001 ([0.05]/13 independent genetic markers]/4 inheritance models).

Haplotype blocks were constructed from the genotyping data of the non-IA control group using the SNPtool (http://www. dkfz.de/de/molgen_epidemiology/tools/SNPtool.html; Chen et al., 2009) and the Haploview software (v.4.2). Selected polymorphisms within the same locus were not in linkage disequilibrium ensuring that these variants represented independent variability signals (Supplementary Figure 1) . In addition to the analysis based on a single SNP, we performed haplotype frequency estimation and haplotype association analysis using SNPStats (http://bioinfo.iconcologia. net/SNPstats). Haplotype frequencies were determined using the Expectation-maximization (EM) algorithm and the minimum haplotype frequency was set at 0.01, therefore assessing association for common haplotypes.

We also assessed whether selected polymorphisms correlated with mRNA expression in the publicly available eQTL IGV browser for primary cells (http://www.gtexportal.org/home/) or using the Haploreg data.

We also were interested in testing whether NFκB-related SNPs could interact with each other to modify the risk of IA using the multifactor dimensionality reduction (MDR) software. A detailed description of the MDR method has been reported elsewhere (Ritchie et al., 2003; Moore, 2004) . A 10-fold cross-validation analysis (exhaustive search) and permutation testing were used to confirm the best interaction models. MDR models were considered statistically significant at P < 0.05 (P sign). Statistical significance of each particular model was then re-evaluated using a 1.000-fold permutation test to compare observed testing balanced accuracies with those expected under the null hypothesis of no association (using the MDR permutation testing module 0.4.9 alpha). Interactions were visualized by constructing an interaction dendrogram according to the method described by Moore et al. (Moore, 2004) . MDR software and MDR permutation testing module are open-source and freely available from http://www.epistasis.org.

In this population-based case-control study a total of 834 hematological patients were recruited. Demographic and clinical characteristics of these patients are summarized in Table 2 . IA and non-IA patients had similar age but male patients were more prone to develop IA than females (male/female ratio = 1.86 vs. 1.18, P = 0.015). HSCT was common and equally distributed between IA and non-IA groups (45.85 vs. 45.49%) and the underlying disease (hematological disorder) was also uniform between both groups ( Table 2) . As expected, patients without prophylaxis were more prone to develop IA compared to those patients who take antifungal drugs (23.28 vs. 14.51%; Table 2 ). All SNPs were in Hardy-Weinberg in the control group with the exception of the REL rs13031237 that was excluded from further analyses (non-IA; P < 0.001). Allele and genotype frequencies of selected SNPs were in line with those reported in Hapmap. Logistic regression analysis adjusted for age, sex, country of origin, allo-HSCT, and underlying disease showed that none of the selected SNPs was significantly associated with the risk of IA according to codominant, dominant, recessive and log-additive models of inheritance (Table 3 ). In addition, although prophylaxis status was only partially available in our population (n = 618), no significant changes in association estimates were observed when this clinical variable was included as covariate for adjustment (data not shown). No significant associations between the NFκB1, NFκB2, cRel, RelB, and IRF4 polymorphisms and IA risk were also found in an allo-HSCTrestricted analysis considering donor genotypes and IA episodes occurred after transplantation ( Table 4) . Although our statistical power was limited (Supplementary Table 1) , we found a fairly weak association of the IRF4 rs12203592T/T genotype with the risk of developing the infection at the nominal significance level of P < 0.05. Carriers of the IRF4 rs12203592T/T genotype showed a six-fold increased risk of developing IA when compared with patients carrying the wild type C-allele (OR REC = 6.24, 95%CI 1.25-31.2, P = 0.026). Interestingly, we found that the presence of the minor allele of the rs12203592 (T) strongly correlated with IRF4 mRNA expression levels in whole blood samples and Epstein-Barr virus (EBV)-transformed lymphocytes (P = 6.0 • 10 −7 and P = 3.0 • 10 −7 , respectively), which suggested a possible functional role of this polymorphism. Based on these interesting results, we decided to explore the potential functional impact of this polymorphism using Haploreg and ENCODE annotation data. This analysis revealed that the IRF4 rs12203592 SNP resides on a strong enhancer and near of an active promoter that might constitute a regulatory element for IRF4. In addition, this analysis showed that this intronic polymorphism was predicted to change binding motifs for NFκB, HDAC2, and HMG-IY, which are proteins implicated in the transcriptional regulation of multiple innate and adaptive immune-related genes (Liu et al., 2001; Hayden et al., 2006; Falkenberg and Johnstone, 2014) . In line with the findings suggesting a functional role of this variant, we also found that the rs12203592 SNP mapped among promoter and enhancer histone marks in different primary T-cell subtypes (TCD8+, TCD4+, Treg, and Th17) of different origin (peripheral blood, hematopoietic stem cells, cord blood, etc.) but also in primary B-cells, natural killer cells and monocytes. Although at this point it was tempting to speculate that the IRF4 rs12203592 SNP might play a role in modulating the risk of developing IA likely through the regulation of IRF4 mRNA expression levels, the association of this SNP with IA risk did not remained significant after correction for multiple testing (P = 0.001; Table 4 ) and, therefore, requires further confirmation. In agreement with the single-SNP analysis, haplotype association analysis showed fairly weak associations of the IRF4 AATC and IRF4 GGTC haplotypes (not including the IRF4 rs12203592T risk allele) with a decreased risk of IA. However, these associations were of similar magnitude to the one observed for the IRF4 rs12203592 SNP in the single-SNP analysis (OR = 0.28, 95%CI 0.08-0.95, P = 0.042 and OR = 0.04, 95%CI 0.00-0.71; Table 5 ) suggesting that the haplotypic effect was due to this intronic variant. Finally, given the involvement of physiological complexes (p50/p65, p52/RelB, and p65-c-Rel) in the activation of the NFκB signaling pathway and the role of the p52/RelB complex in the transcriptional activation of the IRF4 (Boddicker et al., 2015) , we also decided to investigate whether NFκB-related SNPs might interact to each other to modify the risk of IA. Results of the MDR analysis evaluating all possible combinations among the NFκBrelated SNPs are shown in Table 6 and Figure 1 . The best model suggested an interaction between the REL rs842647 , NFκB2 rs1056890 , IRF4 rs7768807 , and REL rs13017599 SNPs to synergistically increase the risk of developing IA (TA = 0.5488, P Sign = 0.0010). Although this model was the best to predict IA and presented a high cross validation consistency (CVC) (9/10), it did not remain statistically significant following 1000-fold permutation test (P = 0.31). In addition, we found a significant 2-locus interaction model including the NFκB2 rs1056890 and IRF4 rs7768807 SNPs to increase the risk of getting the infection (TA = 0.5353, P Sign = 0.0010) but this interaction model also failed to retain statistical significance after 1000-fold permutation test (P = 0.47; Table 6 ). None of the best models included the IRF4 rs12203592 SNP. 

The NFκB pathway is implicated in fostering a wide variety of physiological processes such as immune cell turnover (Hayden et al., 2006) , inflammation (Karin and Greten, 2005) , T-cell differentiation (Th1, Th2, Th17, and Treg subsets; Oh and Ghosh, 2013) , DC maturation (Burkly et al., 1995; Abe et al., 2003) , cell apoptosis (Wang et al., 1998; Karin and Greten, 2005) , cell adhesion (Chen et al., 1995; Lockyer et al., 1998; Chiang et al., 2008) , and angiogenesis (Beinke and Ley, 2004) . In the presence of fungal pathogens such as Aspergillus fumigatus, the canonical NFκB pathway involving NEMO-IKKα-IKKβ and NFκB-p50/p65 complexes can be activated by membrane pattern recognition receptors (PRRs) such as toll-like receptors (TLRs; Mambula et al., 2002; Bellocchio et al., 2004; Roeder et al., 2004) , C-type lectin receptors (CLRs; Serrano-Gómez et al., 2004; LeibundGut-Landmann et al., 2007; Gringhuis et al., 2009; Rivera et al., 2011; Saijo and Iwakura, 2011) , cytokine receptors (IL1R, TNFRs) (Bozza et al., 2014) , or soluble PRRs like pentraxins (PTX3; Garlanda et al., 2002 ) that leads to the immediate secretion of pro-inflammatory cytokines (TNF, IL1α, IL1β, IL6, IL12; Nicholson et al., 1996; Roeder et al., 2004) , chemokines (CCL2, CCL3, CXCL2 and CXCL10; Shahan et al., 1998; Phadke and Mehrad, 2005; Reid et al., 2009) , and reactive oxygen species (ROS; superoxide, hydroxyl radical, nitric oxide; Philippe et al., 2003; Chiang et al., 2008) but also to the expression of specific immune receptors and cell surface adhesion molecules (VCAM-1, E-Selectin; Akira et al., 2001; Chiang et al., 2008; Gringhuis et al., 2009 ). In addition, binding of fungal compounds to C-type lectins (Dectin-1 and Dectin-2) and the activation of certain TNF family cytokines (CD40L, BAFF, or LT-β) observed during fungal infections leads to the activation of the non-canonical NFκB pathway (p52/RelB), which is implicated in B-cell and DC maturation (Hu et al., 2011; Sun, 2012) and γδTh17 cell development (Powolny-Budnicka et al., 2011) but also in the control of the expression of a wide range of immune-related genes (Geijtenbeek and Gringhuis, 2009; Vallabhapurapu and Cross-validation and permutation testing were used to identify the best models. All possible two-way SNP interactions were tested using 10-fold cross-validation and the exhaustive search. The model with the highest testing balanced accuracy (TA) and cross validation consistency (CVC) was selected as "best model." Statistical significance was evaluated by the Sign test and confirmed using a 1.000-fold permutation test to compare observed testing balanced accuracies with those expected under the null hypothesis of no association (using the MDR permutation testing module 0.4.9 alpha).*1000-fold permutation test (α = 0.001, TA = 0.6370; α = 0.01, TA = 0.5984; α = 0.05, TA = 0.5830; α = 0.10, TA = 0.5717).

FIGURE 1 | Interaction dendrogram. The interaction dendrogram reveals no significant interactions among the studied SNPs to modulate the risk of IA. Karin, 2009; Plato et al., 2013) . On the other hand, it has recently been suggested the existence of a CD30-p52/Relb-IRF4 loop to regulate the NFκB pathway and cell proliferation (Boddicker et al., 2015) and a role of IRF4 in modulating the differentiation of different DC (Bajaña et al., 2016) and Th17-mediated immune responses against A. fumigatus (Schlitzer et al., 2013) .

Considering the central role of the PRRs-NFκB pathways and IRF4 in the immune responses against A. fumigatus (Oh and Ghosh, 2013; Schlitzer et al., 2013; Williams et al., 2013) and those studies suggesting that genetic host factors may account for differences in susceptibility to IA (Ok et al., 2011; van der Velden et al., 2011) , we decided to investigate the link between genetic polymorphisms within NFκB1, NFκB2, RelA, RelB, Rel, and IRF4 genes and the risk of IA. In spite of the growing number of studies assessing the role of genetic polymorphisms within PRRs (TLR2, TLR4, TLR9, Dectin-1, Dectin-2, DC-SIGN, PTX3; Kesh et al., 2005; Bochud et al., 2008; Pamer, 2008; Cunha et al., 2010 Cunha et al., , 2011 Cunha et al., , 2014 Chai et al., 2011; Sainz et al., 2012; Grube et al., 2013; Wojtowicz et al., 2015b) , cytokines (IL1, IL10, IFNG; Sainz et al., 2007a Sainz et al., , 2008a Lupiañez et al., 2015) , and their receptors (IL4R, TNFR1, and TNFR2; Sainz et al., 2007b Sainz et al., , 2010 Lupiañez et al., 2015) in determining the susceptibility to invasive fungal infections, this is the first study that attempts to evaluate the impact of polymorphisms within NFκB1, NFκB2, RelA, RelB, REL, and IRF4 genes on the risk of developing IA. Our data showed no significant overall associations between selected SNPs and IA infection. The best result was the association of the IRF4 rs12203592 SNP with an increased risk of IA in HSCT patients that suggested a weak effect of this locus on the risk of IA that might become evident only in those patients with a more profound degree of immunosuppression. However, despite the potential interest of these results and the evidences supporting a functional role of the IRF4 rs12203592 variant in regulating IRF4 mRNA expression and thereby modulating the NFκB pathway (Boddicker et al., 2015) and IRF4-dependent immune responses (Schlitzer et al., 2013; Bajaña et al., 2016) , the association of this variant did not remain significant at the experiment-wide significance threshold, suggesting that the association of this polymorphism with IA could be due to chance. Haplotype analysis showed a small effect of common haplotypes on IA risk but the magnitude of these effects suggested that the observed association was likely due to the IRF4 rs12203592 SNP rather than the combined effect of the IRF4 SNPs. In addition, when we evaluated whether there was any significant interaction among the SNPs analyzed, we did not observe any consistent interaction model that could affect the risk of developing IA.

In conclusion, this case-control study does not provide strong evidence of a relationship between polymorphisms within NFκB1, NFκB2, cRel, RelB, and IRF4 genes and IA risk. Nonetheless, given the limited statistical power of the HSCTstratified analysis (80% to detect odds ratio of 2.1 at α = 0.001 for a SNP with a frequency of 0.25, dominant model) and the evidences suggesting a functional role of the IRF4 rs12203592 SNP, we cannot dismiss the possibility of a small but still real effect of this variant or its haplotypes on the risk of IA in HSCT patients. Future case-control population-based studies conducted in larger HSCT populations are now warranted to further evaluate whether the IRF4 locus may have a role in determining the susceptibility to IA.

The study protocol was approved by the local ethics review boards of all participating centers and written informed consent was obtained from each patient before inclusion in accordance with the Declaration of Helsinki. Ethical approval for this study was provided by the Comité de Ética e Investigación Clinica (CEIC) of the Virgen de las Nieves Hospital (Granada), University Hospital of Salamanca (Salamanca), Clinic University Hospital of Valencia (Valencia), and the Centro Nacional de Microbiologia (27_2012). The PCRAGA trial is registered with ClinicalTrials.gov (NCT01742026) and EudraCT (2010-019406-17) . Ethical approval was also provided by the Subcomissao de Etica para as Ciencias da Vida e Saude ( JS analysed the data. MJ and JS drafted the manuscript. All authors read and approved the final version of the manuscript.

This study was supported by grants PI12/02688 from Fondo de Investigaciones Sanitarias (Instituto de Salud Carlos III, Madrid, Spain), the ERA-NET PathoGenoMics (03159000A; Ministerio de Ciencia e Innovación PIM2010EPA-00756, Madrid, Spain), the Collaborative Research Center / Transregio 124 FungiNet, the Austrian Science Fundation (FWF I-656-B09), the Fundação para a Ciência e Tecnologia (FCT), cofunded by Programa Operacional Regional do Norte (ON.2-O Novo Norte), the Quadro de Referência Estratégico Nacional (QREN) through the Fundo Europeu de Desenvolvimento Regional (FEDER) and the Projeto Estratégico -LA 26 -2013-2014 (PEst-C/SAU/LA0026/2013). Agostinho Carvalho and Cristina Cunha were supported by the Fundação para a Ciência e Tecnologia (FCT), Portugal (IF/00735/2014 and SFRH/BPD/96176/2013, respectively). The PCRAGA trial was supported by an unrestricted grant from Pfizer, which had no involvement or control over the collection, analysis, and interpretation of data; the writing of the report; or the decision to submit the paper for publication. This study was also supported by Astellas Pharma Inc. and a donation from Consuelo González Moreno, an acute myeloid leukemia survivor.

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.

Supplementary Table 1 

Distinct contributions of TNF and LT cytokines to the development of dendritic cells in vitro and their recruitment in vivo

Serum galactomannan versus a combination of galactomannan and polymerase chain reaction-based Aspergillus DNA detection for early therapy of invasive aspergillosis in high-risk hematological patients: a randomized controlled trial

Toll-like receptors: critical proteins linking innate and acquired immunity

An investigation of rheumatoid arthritis loci in patients with early-onset psoriasis validates association of the REL gene

IRF4 and IRF8 Act in CD11c+ cells to regulate terminal differentiation of lung tissue dendritic cells

Functions of NF-κB1 and NF-κB2 in immune cell biology

TLRs govern neutrophil activity in aspergillosis

Toll-like receptor 4 polymorphisms and aspergillosis in stem-cell transplantation

The oncogenic transcription factor IRF4 is regulated by a novel CD30/NF-κB positive feedback loop in peripheral T-cell lymphoma

Comprehensive assessment of rheumatoid arthritis susceptibility loci in a large psoriatic arthritis cohort

PTX3 binds MD-2 and promotes TRIF-dependent immune protection in aspergillosis

Expression of relB is required for the development of thymic medulla and dendritic cells

TLR3 essentially promotes protective class I-restricted memory CD8(+) T-cell responses to Aspergillus fumigatus in hematopoietic transplanted patients

The Y238X stop codon polymorphism in the human β-glucan receptor dectin-1 and susceptibility to invasive aspergillosis

SNP_tools: A compact tool package for analysis and conversion of genotype data for MS-Excel

Selective inhibition of E-selectin, vascular cell adhesion molecule-1, and intercellular adhesion molecule-1 expression by inhibitors of I κ B-α phosphorylation

Genetic variation in the REL gene increases risk of Behcet's Disease in a Chinese han population but that of PRKCQ does not

Nasopharyngeal shedding of severe acute respiratory syndromeassociated coronavirus is associated with genetic polymorphisms

Aspergillus fumigatus stimulates leukocyte adhesion molecules and cytokine production by endothelial cells in vitro and during invasive pulmonary disease

Genetic PTX3 deficiency and aspergillosis in stem-cell transplantation

Dectin-1 Y238X polymorphism associates with susceptibility to invasive aspergillosis in hematopoietic transplantation through impairment of both recipient-and donor-dependent mechanisms of antifungal immunity

Genetically-determined hyperfunction of the S100B/RAGE axis is a risk factor for aspergillosis in stem cell transplant recipients

Exploring multilocus associations of inflammation genes and colorectal cancer risk using hapConstructor

Revised definitions of invasive fungal disease from the European Organization for Research and Treatment of Cancer/Invasive Fungal Infections Cooperative Group and the National Institute of Allergy and Infectious Diseases Mycoses Study Group (EORTC/MSG) Consensus Group

An intronic polymorphism of IRF4 gene influences gene transcription in vitro and shows a risk association with childhood acute lymphoblastic leukemia in males

Genome-wide meta-analysis of psoriatic arthritis identifies susceptibility locus at REL

Confirmation of association of the REL locus with rheumatoid arthritis susceptibility in the UK population

Histone deacetylases and their inhibitors in cancer, neurological diseases and immune disorders

Non-redundant role of the long pentraxin PTX3 in anti-fungal innate immune response

Signalling through C-type lectin receptors: shaping immune responses

Introduction to NF-κB: players, pathways, perspectives

REL, encoding a member of the NF-κB family of transcription factors, is a newly defined risk locus for rheumatoid arthritis

Dectin-1 directs T helper cell differentiation by controlling noncanonical NF-κB activation through Raf-1 and Syk

TLR5 stop codon polymorphism is associated with invasive aspergillosis after allogeneic stem cell transplantation

A germline variant in the interferon regulatory factor 4 gene as a novel skin cancer risk locus

NF-κB and the immune response

Transcriptional regulation via the NF-κB signaling module

Noncanonical NF-κB regulates inducible costimulator (ICOS) ligand expression and T follicular helper cell development

NF-κB: linking inflammation and immunity to cancer development and progression

Invasive aspergillosis: new insights into disease, diagnostic and treatment

TLR1 and TLR6 polymorphisms are associated with susceptibility to invasive aspergillosis after allogeneic stem cell transplantation

Pulmonary aspergillosis: a clinical review

Syk-and CARD9-dependent coupling of innate immunity to the induction of T helper cells that produce interleukin 17

Cytokine candidate genes predict the development of secondary lymphedema following breast cancer surgery

The HMG I proteins: dynamic roles in gene activation, development, and tumorigenesis

Inhibition of nuclear factor-κB-mediated adhesion molecule expression in human endothelial cells

Polymorphisms in host immunitymodulating genes and risk of invasive aspergillosis: results from the AspBIOmics Consortium

Toll-like receptor (TLR) signaling in response to Aspergillus fumigatus

Invasive aspergillosis in the intensive care unit

Polymorphisms in the chemokine (C-X-C motif) ligand 10 are associated with invasive aspergillosis after allogeneic stem-cell transplantation and influence CXCL10 expression in monocyte-derived dendritic cells

Computational analysis of gene-gene interactions using multifactor dimensionality reduction

The NF-κB pathway

Epidemiology and outcome of invasive fungal infection in adult hematopoietic stem cell transplant recipients: analysis of Multicenter Prospective Antifungal Therapy (PATH) Alliance registry

Epidemiology, outcomes, and risk factors of invasive fungal infections in adult patients with acute myelogenous leukemia after induction chemotherapy

Epidemiology, outcomes, and mortality predictors of invasive mold infections among transplant recipients: a 10-year, single-center experience

Inhibition of the transcription factors NF-κ B and AP-1 underlies loss of cytokine gene expression in rat alveolar macrophages treated with a diffusible product from the spores of Aspergillus fumigatus

A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other

NF-κB: roles and regulation in different CD4(+) T-cell subsets

Genetic susceptibility to Aspergillus fumigatus infections

TLR polymorphisms and the risk of invasive fungal infections

Identification of haplotype tag single nucleotide polymorphisms within the nuclear factor-κB family genes and their clinical relevance in patients with major trauma

Cytokines in host defense against Aspergillus: recent advances

Killing of Aspergillus fumigatus by alveolar macrophages is mediated by reactive oxidant intermediates

C-type lectin-like receptors of the dectin-1 cluster: ligands and signaling pathways

RelA and RelB transcription factors in distinct thymocyte populations control lymphotoxin-dependent interleukin-17 production in gammadelta T cells

Toll-like receptor 9 modulates immune responses to Aspergillus fumigatus conidia in immunodeficient and allergic mice

Pattern recognition: recent insights from Dectin-1

Power of multifactor dimensionality reduction for detecting gene-gene interactions in the presence of genotyping error, missing data, phenocopy, and genetic heterogeneity

Dectin-1 diversifies Aspergillus fumigatus-specific T cell responses by inhibiting T helper type 1 CD4 T cell differentiation

Toll-like receptors as key mediators in innate antifungal immunity

Dectin-1 and Dectin-2 in innate immunity against fungi

Interleukin-10 promoter polymorphism as risk factor to develop invasive pulmonary aspergillosis

Dectin-1 and DC-SIGN polymorphisms associated with invasive pulmonary Aspergillosis infection

IL1 gene cluster polymorphisms and its haplotypes may predict the risk to develop invasive pulmonary aspergillosis and modulate C-reactive protein level

Genetic variants of IL6 gene promoter influence on C-reactive protein levels but are not associated with susceptibility to invasive pulmonary aspergillosis in haematological patients

Variable number of tandem repeats of TNF receptor type 2 promoter as genetic biomarker of susceptibility to develop invasive pulmonary aspergillosis

TNFR1 mRNA expression level and TNFR1 gene polymorphisms are predictive markers for susceptibility to develop invasive pulmonary aspergillosis

IRF4 transcription factor-dependent CD11b+ dendritic cells in human and mouse control mucosal IL-17 cytokine responses

Dendritic cell-specific intercellular adhesion molecule 3-grabbing nonintegrin mediates binding and internalization of Aspergillus fumigatus conidia by dendritic cells and macrophages

IkappaBKbeta and NFκB1, NSAID use and risk of colorectal cancer in the Colon Cancer Family Registry

Concentration-and time-dependent upregulation and release of the cytokines MIP-2, KC, TNF, and MIP-1alpha in rat alveolar macrophages by fungal spores implicated in airway inflammation

Aspergillosis in solid organ transplantation

Interferon-signaling pathway: associations with colon and rectal cancer risk and subsequent survival

Polymorphisms in cytokine genes IL6, TNF, IL10, IL17A and IFNG influence susceptibility to complicated skin and skin structure infections

Clinical epidemiology of 960 patients with invasive aspergillosis from the PATH Alliance registry

The noncanonical NF-κB pathway

Coeliac disease-associated risk variants in TNFAIP3 and REL implicate altered NF-κB signalling

Regulation and function of NF-κB transcription factors in the immune system

Genetic variants and the risk for invasive mould disease in immunocompromised hematology patients

Analysis of the REL polymorphism rs13031237 in autoimmune diseases

NF-κB antiapoptosis: induction of TRAF1 and TRAF2 and c-IAP1 and c-IAP2 to suppress caspase-8 activation

Association of interferon regulatory factor 4 gene polymorphisms rs12203592 and rs872071 with skin cancer and haematological malignancies susceptibility: a meta-analysis of 19 case-control studies

Transcription factor IRF4 drives dendritic cells to promote Th2 differentiation

IL1B and DEFB1 polymorphisms increase susceptibility to invasive mold infection after solid-organ transplantation

PTX3 polymorphisms and invasive mold infections after solid organ transplant