Skip to main content
Download PDF

The landscape of human genes involved in the immune response to parasitic worms

BMC Evolutionary Biology volume 10, Article number: 264 (2010) Cite this article

Abstract

Background

More than 2 billion individuals worldwide suffer from helminth infections. The highest parasite burdens occur in children and helminth infection during pregnancy is a risk factor for preterm delivery and reduced birth weight. Therefore, helminth infections can be regarded as a strong selective pressure.

Results

Here we propose that candidate susceptibility genes for parasitic worm infections can be identified by searching for SNPs that display a strong correlation with the diversity of helminth species/genera transmitted in different geographic areas. By a genome-wide search we identified 3478 variants that correlate with helminth diversity. These SNPs map to 810 distinct human genes including loci involved in regulatory T cell function and in macrophage activation, as well as leukocyte integrins and co-inhibitory molecules. Analysis of functional relationships among these genes identified complex interaction networks centred around Th2 cytokines. Finally, several genes carrying candidate targets for helminth-driven selective pressure also harbour susceptibility alleles for asthma/allergy or are involved in airway hyper-responsiveness, therefore expanding the known parallelism between these conditions and parasitic infections.

Conclusions

Our data provide a landscape of human genes that modulate susceptibility to helminths and indicate parasitic worms as one of the major selective forces in humans.

Background

Helminth infections are estimated to infect about 2 billion individuals worldwide (reviewed in [1]). Although rarely fatal, these parasites cause high rates of morbidity by establishing chronic infections. In particular, the highest parasite burdens are observed in pre-school and school-aged children, often resulting in anemia, undernourishment and growth stunting (reviewed in [1]). During pregnancy helminth infection is a risk factor for preterm delivery, reduced birth weight and maternal mortality (reviewed in [1]). Moreover, by chronically infecting their host, parasitic worms increase the susceptibility to other pathogens such as viruses, bacteria and protozoa [1]. Previous works have indicated that the intensity of helminth infection is a heritable trait, although measures of heritability vary among studies and depend on the parasite beanalyzed [2].

These observations suggest that helminth infections have represented a very strong selective pressure for humans, a selective pressure that is very likely to also be remarkably constant over time. Indeed, most vertebrates have been hosting a wide range of parasitic worms for million years and humans have had their share before emerging as a species (reviewed in [3]). We have previously addressed the role of helminths as selective agents in human evolution analyzing a large set of human genes encoding interleukins and their receptors; we demonstrated that the pressure imposed by parasitic worms on these genes has been stronger than the one due to viral and microbial agents [4]. The reasons for this observation likely lie in the long-term relationship between humans and helminths, in the relatively slow evolutionary rates of these parasites and in their geographic distribution being considerably stable. Here we aimed at exploiting the selection signatures left by these pathogens on human genes to identify, at the genome-wide level, candidate genes and variants that may have been subjected to helminth-driven selective pressure.

Results

Helminth diversity and prevalence correlate across geographic locations

In order to search for candidate variants subjected to helminth-driven selection, an estimate of the selective pressure exerted by these pathogens needs to be defined. It has been previously suggested [4–6] that the selective pressure exerted by infectious diseases in different geographic areas can be estimated by counting the number of different pathogen species/genera that are transmitted in these regions. In particular, in the case of helminthiases, we consider parasite diversity to be a better estimate of helminth-driven selective pressure than prevalence for different reasons. First of all, comprehensive data on prevalence are impossible to retrieve for many parasite species/genera or countries and even when available, prevalence data may vary considerably within the same country depending on the surveyed regions and the population surveyed (e.g. city dwellers rather than farmers/bushmen/nomads, children rather than adults). Moreover, the prevalence of specific helminthiases might have changed greatly over recent years as a result of eradication campaigns, and historical prevalence data are rarely available. Also, it should be considered that prevalence data are difficult to combine since in endemic regions polyparasitism is common [1]; indeed, in these regions subjects infected with multiple helminths tend to harbor the most intense infections, possibly due to an additive and/or multiplicative impact on nutrition and organ pathology (reviewed in [7]).

These same observations also apply to other possible measures such as parasite burden or infection pathogenicity, which are very difficult to quantify and that may have varied considerably along human evolutionary history. Conversely, the diversity of helminth species transmitted in different geographic locations has been shown to depend upon climatic variables [8] which, in turn, may be considered as relatively stable over time, suggesting that diversity may better describe long-term evolutionary pressures. Thus, we calculated helminth diversity from the Global Infectious Disease and Epidemiology Network database (Gideon); as described in the method section, all entries for single helminthiases were manually inspected and all species/genera transmitted in a given location were counted as present irrespective of their prevalence. The number of different helminth species/genera per country is reported in Additional file 1 (Table S1).

We next wished to verify whether the prevalence of the most common helminth infections (as reported in [1]) correlates with helminth diversity across the 52 populations genotyped in the HGDP-CEPH panel (Additional file 1, Table S1). Again, we retrieved prevalence and diversity data from Gideon by manually inspecting single and retrieving all prevalence surveys. The prevalence of single helminthiases per country was obtained by averaging all surveys; data obtained on HIV-seropositive individuals were not included in the analyses because of the known modulation of HIV on the susceptibility on helminth-based infections [1, 9, 10]. Finally, countries with no survey data were not included in the analysis for the specific helminth. Parasites were then grouped into five major classes (following [1]). As shown in table 1, the prevalence of all parasite groups correlated with helminth diversity.

Table 1 Correlation between the prevalence of all parasite groups and helminth diversity
Full size table

SNPs associated with helminth diversity are over-represented in immune response genes

Previous analysis have indicated that several genes coding for interleukins and interleukin receptors are subjected to helminth-driven selective pressure [4]. Interleukins are central mediators of immunity and inflammation and, in general, we might expect genes with a role in immune response to be preferential targets of helminth-driven selective pressure. Specifically, we expect variants within these genes to be more frequently associated with helminth diversity than observed for randomly sampled loci. We verified this prediction by analyzing the ImmPort list which contains 2,287 genes involved in immune response and covered in the HGDP-CEPH panel, this latter containing data for more than 660,000 single nucleotide polymorphisms genotyped in almost 950 individuals sampled throughout the world (Additional file 1, Table S1) [11].

We calculated Kendall's τ rank correlation coefficient between allele frequencies of SNPs in ImmPort genes and helminth diversity. A SNP was defined as being significantly associated with helminth-diversity if it displayed a significant correlation (p value after Bonferroni correction < 0.01; uncorrected p value < 1.94 × 10-7) and a τ value higher than the 95th percentile in the distribution of correlation coefficients calculated over all SNPs having minor allele frequency (MAF) similar (in the 1% range) to that of the SNP being analyzed. This latter requirement stems from the need to account for the influence of non-selective events, as a few HGDP-CEPH population result from recent or ancient admixture [11] and population history, migratory events and genetic drift have affected human genetic variability [12, 13]. Among 2,287 ImmPort genes, 246 contained at least one SNP significantly associated with helminth diversity (Additional file 2, Table S2). The likelihood to obtain an equal or higher number of genes carrying significantly associated SNPs was assessed by a re-sampling approach. Specifically, we divided all genes with at least one SNP typed in the HGDP-CEPH panel in 24 intervals based on the number of typed SNPs; 10,000 re-samplings were then performed by selecting for each ImmPort gene, a randomly selected gene with a similar number of SNPs (i.e. a gene located in the same interval) (see methods). The number of SNPs in ImmPort genes did not differ significantly from the average number in the re-samplings (p = 0.22) and the empirical probability of obtaining 246 genes with at least one significant SNP resulted equal to 0.041, indicating that immune response genes more frequently display variants correlating with helminth diversity compared to randomly chosen loci.

When the same analysis was performed using the prevalence of soil-transmitted nematodes, filarial nematodes, schistosomes and food-borne trematodes no significant enrichment for ImmPort genes was observed (empirical p values = 0.83, 0.96, 0.45, and 0.39, respectively).

In agreement with previous suggestions [4], these data indicate that helminth diversity may be a reliable estimator helminth-driven selective pressure.

Genome-wide search for variants subjected to helminth-driven selective pressure

Given these results, we wished to identify SNPs significantly associated with helminth diversity on a genome-wide base. We therefore calculated Kendall's rank correlations between allele frequency and helminth diversity for all SNPs (n = 660,832) typed in the HGDP-CEPH panel. We next searched for instances which withstood Bonferroni correction with α = 0.05 (i.e. uncorrected p value < 7.6 × 10-8) and displayed a τ percentile rank higher than the 95th among MAF-matched SNPs. A total of 3,478 SNPs mapping to 810 distinct genes satisfied both requirements (Additional file 3, Table S3). We next verified whether climatic variables could be responsible for the correlations detected between these SNPs and helminth diversity. Hence, for all countries where HGDP-CEPH populations are located we obtained the following parameters: average annual minimum and maximum temperature, short wave radiation flux and precipitation rate (annual maximum and mean). None of these SNPs withstood Bonferroni corrections in these analyses.

Previous works have reported an enrichment of selection signatures within or in close proximity to human genes [12, 14–17]. In line with these data we verified that helminth-associated SNPs are more frequently located within gene regions compared to a control set of MAF-matched variants (χ2 test, p = 0.014).

A full list of the 3,478 SNPs that showed a significant correlation with helminth diversity is available as Additional file 3 (Table S3). Table 2 shows the 20 strongest correlations together with a short comment on the possible role of candidate genes in immune response or helminth resistance.

Table 2 Top 20 SNP (or SNP clusters) associated with helminth diversity.
Full size table

Among the 810 genes subjected to helminth-driven selective pressure, we identified ectodysplasin A (EDA) and its receptor (EDAR). EDAR has been subjected to strong positive selection in populations of Asian descent [17] and is responsible for hair thickness in these populations [18]. One coding variant (rs3827760, 370Val/Ala) in EDAR is thought to be the selection target but it has not been genotyped in the HGDP-CEPH panel. We therefore wished to verify whether the variants we found to correlate with helminth diversity are located on the same haplotype as the selected 370Ala allele. We used the Sweep software [19] to analyze the haplotype structure in the genomic region encompassing EDAR; we selected a core containing the 370Ala/Val variant and used HapMap data from Asian individuals (Chinese and Japanese). The results indicated that most chromosomes carrying the putatively selected 370Ala variant also harbor four alleles we found to be significantly correlated with helminth diversity (Figure 1).

Figure 1
figure1

Haplotype analysis of the genomic region encompassing EDAR. Haplotype bifurcation diagram for the genomic region previously identified as being subjected to a selective sweep (see text). The location of known genes in the region is also shown. As described in the text, we selected a core containing the putatively selected 370Val/Ala variant (dark blue circle). The four SNPs that correlate with helminth diversity are represented as light blue circles.

Full size image

Functional characterization of genes subjected to helminth-driven selective pressure

We investigated the role and functional relationship among helminth diversity-associated genes using the Ingenuity Pathway Analysis and the PANTHER classification system [20–22]. For these analyses, a SNP was ascribed to a given gene if it was located within the transcribed region or in the 25 kb upstream of the transcription start site.

Unsupervised IPA analysis retrieved five high-scoring networks (p < 10-12) (Figure 2 and Additional file 4, Figure S1) and two additional networks with lower scores (p < 10-9). The two highest scoring networks were merged, as well as networks 3 and 5. As shown in figure 2, 37 and 32 genes in merged networks 1 and 2 correlated with helminth diversity, respectively, corresponding to almost 60% of network nodes.

Figure 2
figure2

Network analysis of genes associated with helminth diversity. Genes are represented as nodes, edges indicate known interactions between proteins (sold lines depicts direct and dashed lines depict indirect interaction). Genes are color coded as follows: green, genes with at least one SNP significantly associated with helminth diversity; gray, genes covered by at least one SNP in the HGDP-CEPH panel; white, genes with no SNPs in the panel. Panel A and B represents the networks obtained by merging networks 1 with 2 and 2 with 5, respectively.

Full size image

We next investigated the over-representation of PANTHER classification categories among genes significantly associated with helminth diversity. Table 3 shows the 5 significantly over-represented PANTHER pathways with the contributing genes, as well as the most significantly over-represented molecular functions and biological processes. Notably, genes involved in cytokine-mediated inflammation and integrin signaling accounted for two significantly over-represented pathways. While the multiple functions of integrins and cytokines on the immune system are established, the role of glutamate and α adrenergic receptor signaling in immune related processes are less understood. Recent evidence has indicated that ionotropic glutamate receptors are expressed by T lymphocytes (reviewed in [23]) and glutamate can exert different actions on these cells including triggering adhesion, proliferation and chemotaxis (reviewed in [23]). Similarly, α adrenergic receptors are expressed by different immunocompetent cells (reviewed in [24]); in particular, expression of ADRA1A is induced in human monocytes by inflammatory cytokines [25]. Finally, it is interesting to notice that components of the treleasing hormone (TRH) receptor signaling pathway are also over-represented among genes subjected to helminth driven selective pressure. Tacts on specific receptors within the pituitary gland to stimulate the release of thyroid stimulating hormone and prolactin. The immunomodulatory role of thyroid hormones on immune functions are clearlyestablished (reviewed in [26]). Moreover, experiments in mice have indicated that thyroxine plays a role in the establishment of Schistosoma mansoniinfection [27, 28] and animals treated with the hormone display increased parasite numbers and development of giant worms [28]. Notably, among variants correlating with helminth diversity at the genome-wide level, we also found one SNP relatively close to the gene encoding prolactin (rs13198653, τ = 0.64, p = 4.4 × 10-9) and one variant within PRLR (prolactin receptor, rs4235652, τ = 0.62, p = 1.8 × 10-9) (Additional file 3, Table S3). PRL was shown to be an immunomodulator and acts as a cytokine on many different immune cells; indeed, PRLR is expressed by B, T and NK cells, as well as macrophages (reviewed in [26]) and PRL expression in T cells is regulated by IL-2, IL-4 and IL-1β. These data therefore suggest a role for both thyroid hormones and prolactin in the resistance to parasitic worms.

Table 3 Panther over-represented categories among genes showing correlation with helminth diversity
Full size table

Helminth-driven selection and susceptibility to allergy and asthma

We next wished to analyze the relationship between variants/genes associated with helminth diversity and the genetic susceptibility to asthma and allergy. We searched among published genome-wide association studies (GWAS) for SNPs that have been associated with allergy, asthma or related traits (serum IgE levels and plasma eosinophil count). Only 12 SNPs were retrieved, 9 of them genotyped in the HGDP-CEPH panel. One of these SNPs (rs12619285) displayed a significant correlation (p = 5.8 × 10-8) with helminth diversity (Table 4) and the allele associated with high eosinophil counts [29] positively correlated with the diversity of parasitic worms. In order to gain further insight into this issue, we focused on genes rather than variants, in line with the view that the gene rather than the allele should be regarded as the replication unit. Given the low consistency that often plagues association studies, only robust allergy/asthma susceptibility genes were considered (see methods). As shown in Table 4, we observed that 12 allergy/atopy genes displayed at least one SNP (either genic or intergenic) significantly associated, at the genome-wide level, with helminth diversity. Among these genes, one SNP in IL4 that has been associated with asthma (rs2070874, +33C/T) has also been genotyped in the HGDP-CEPH panel; again, the allele that positively correlates with helminth diversity (T) is associated with asthma [30].

Table 4 SNPs that correlate with helminth diversity in asthma/allergy genes
Full size table

It is worth noting that we found several SNPs located upstream the transcription start site of CLTA4 and subjected to helminth driven selective pressure. CLTA4 is located on the long arm of chromosome 2, telomeric to CD28. This raises the possibility that the observed allele associations at these two genes (see Figure 3) derive from linkage to a single selected allele. Yet, analysis of linkage disequilibrium (LD) (Additional file 5, Figure S2) indicates that LD is not extensive across the genomic region (also due to the presence of a recombination hot-spot in between the two genes), suggesting that CD28 and the promoter region of CTLA4 are independent selection targets.

Figure 3
figure3

Analysis of SNP located in the genomic region encompassing CD28 and CTLA4. SNPs significantly associated with helminth diversity are shown in red while the region covered by CD28 and CTLA4 are shown in blue. The location of DNAse hypersensitive sites in CD4+ T cells is shown in gray while recombination hotspots are in black. The image was generated by using the "add custom track" utility available through the UCSC Genome Browser.

Full size image

Discussion

Here we propose that candidate susceptibility genes for parasitic worm infections can be identified by searching for SNPs that display a strong correlation with the diversity of helminth species/genera transmitted in different geographic areas. Our approach relies on the assumption that helminth-driven selective pressure has affected the spatial distribution of these variants and it may suffer from a few limits and caveats. First, we used the diversity of parasitic worms as a measure of helminth-driven selective pressure for the reasons reported above. Yet, helminth diversity may be affected by report biases (e.g. under-reporting in developing countries due to limited research and clinical facilities) and it weights equally all helminth species irrespective of their prevalence and disease burden. As an example, rare helminth species that are reported in Gideon and included in our diversity measure may have exerted a very limited selective force (although prevalence may have changed over time and it is difficult to infer selective pressure for single species). In order to evaluate the impact of rare species, we recalculated helminth diversity by taking into account only parasites that are common in at least one country (see methods and Additional file 1, Table S1). This diversity measure strongly correlated (τ = 0.75, p < 10-5) with parasite diversity calculated over all species, suggesting that the inclusion of rare helminths should not largely affect our results. Another possible confounding factor is accounted for by the co-variation of helminth diversity with other environmental variables (e.g. climate and other infectious agents). We verified that variants associated with helminth diversity are significantly enriched within immune response genes and that the SNPs we identified do not correlate with climatic factors, suggesting that climate does not act as an important confounding factor. Yet, we cannot rule out the possibility that other infectious agents have affected the spatial distribution of the variants we identified as the diversity of helminths correlates with that of other human pathogens [4, 8] across geographic locations.

Finally, although we searched for SNPs strongly associated with helminth diversity (uncorrected p value < 7.6 × 10-8) and we applied a correction based on the MAF-matched distribution of Kendall's rank correlation coefficients, we cannot exclude that the spatial distribution of a fraction of the variants we identified is due to population demography, migration history and drift, therefore representing false positives.

Despite these limitations we were able to identify several genes that can be regarded as good candidates as modulators of susceptibility to helminth infection as testified by our interaction network and PANTHER analyses. Mammalian hosts respond to parasitic worms in a relatively uniform manner by producing specific cytokines (mainly IL-4, IL-10, IL-5 and IL-13) and IgE, as well as through the activation of effector cells such as eosinophils, basophils and mast cells [31]. Overall, the response to helminth infection is Th2-dominated and serves to both oppose the parasite and to contain tissue-damage. In line with this concept, IL4 and IL10, as well as IL13 are hub genes in two interaction networks we identified (Figure 2). Nonetheless, the role of other immune components during helminth infections is becoming increasingly clear. In addition to Th1 cells that mediate host response in some stages of Schistosoma and Brugia malayi infection (reviewed in [31, 32]), the role of regulatory T cells (Treg) is now recognized (reviewed in [31, 32]). An increase in Tregs has been observed in different experimental mouse models of helminthiases and in infected humans (reviewed in [31]). Notably, among the 20 genes more strongly associated with helminth diversity (Table 2) we identified two loci, namely CD200R1L [33] and STIM2 [34] that are involved in the development and function of Treg cells. Additionally, we found several SNPs located upstream the transcription start site of CTLA4 to strongly correlate with helminth diversity (Table 4). The gene encodes a co-inhibitory lymphocyte molecule that is preferentially expressed by Treg cells [35] and is thought to be at least partially responsible for the hypo-responsive phenotype of Th2 effector cells (referred to as "conditioned Th2") which is often observed in helminth infections (reviewed in [31]). As an example, blockade of CTLA-4 during Nippostrongylus brasiliensis infection results in higher Th2 cytokine production and decreased parasite numbers [36]. CTLA-4 competes with CD28 for binding to CD86 (and CD80) (reviewed in [37]). Both CD28 and CD86 carry variants significantly associated with helminth diversity (Figure 3 and Additional file 3, Table S3) and their binding provides a co-stimulatory signal for naive T cells; however, binding of CTLA-4 to CD86/CD80 on dendritic and T cells leads to functional inhibition (reviewed in [37]). Similarly to CTLA4, both CD86 and CD28 have been implicated in the immune response against helminths. Specifically, previous studies have indicated that anti-CD86 treatment blocks immune response to Schistosoma mansoni and Heligmosomoides polygyrus and cd28-/- mice display increased susceptibility to S. mansoni [38]. Interestingly, we also found CD247, another co-inhibitory molecule to correlate with helminth diversity (Additional file 3, Table S3). A recent report [39] has shown that Schistosoma induces anergy of CD4+ and CD8+ T cells by up-regulation of CD247 expression on macrophages. These latter cells are considered an important component of anti-helminth response and an alternative form of macrophages has been described in subjects infected by parasitic worms. These cells up-regulate arginase instead of iNOS and express specific molecules including RETNLB (in humans the entry corresponding to FIZZ1/retnla has been discontinued and replaced with RETNLB) and CHIA (acidic chitinase or AMCase) (reviewed in [31, 32]). Notably, we found one SNP in RETNLB to significantly correlate with helminth diversity (Additional file 3, Table S3); with respect to CHIA, we noticed that one variant (rs10494133) displayed a strong correlation with helminth diversity although it did not withstand Bonferroni correction at the genome-wide level (τ = 0.49, p = 3.3 × 10-6).

We found several integrins and adhesion molecules to correlate with helminth diversity (Table 3). Among these, ITGAM, ITGAL and ITGAX (Figure 2) encode integrin α chains that combine with the β 2 chain (encoded by ITGB2) to form leukocyte-specific heterodimeric integrins. These molecules regulate lymphocyte adhesion and transendothelial migration, playing therefore a central role in inflammatory processes. ITGAL and ITGAM are bound by neutrophil inhibitory factor (NIF), an antiadhesive glycoprotein isolated from the canine hookworm Ancylostoma caninum [40, 41], suggesting that leukocyte integrins are relevant to the immune response to helminth infections. Interestingly, more recent evidence [42] has revealed that the genome of the human parasite Necator americanus encodes at least 9 genes with similarity to NIF, suggesting that their products might play a similar role in establishing an immunocompromised niche for the parasite. The ITGAM/ITGB2 and ITGAX/ITGB2 integrins bind iC3b (Figure 2), a cleavage product of complement component 3 (C3). Previous studies have shown that iC3b is deposited on N. brasiliensis larvae [43, 44] and c3 deficient mice carry high lung larval burdens [43]. In line with these findings, C3 is essential for killing Strongyloides stercoralis larvae in mice [45] and c3-/- mice do not develop an effective Th2 response after infection with S.mansoni and cannot clear the parasite after chemotherapy [46].

As far as the second network is concerned (Figure 2), it is worth mentioning that NOTCH1 and ETS1 have been implicated in multiple immune functions. The NOTCH signaling pathway is involved in the intrathymic differentiation of T cells, as well as in Th cell development in the periphery [47]. In line with the role of Treg cells in helminth infection, NOTCH1 has recently been shown to be involved in Treg function [48] and to cooperate with TGFb for regulation of the FOXP3 promoter [49]. With respect to ETS1, it functions as a transcriptional regulator of several cytokine genes including IL5, IL2 and GMCSF [50–52]. It also regulates expression of CD226 [53] and it is known to induce of Th1 mediated inflammation [54].

Also, network B (Figure 2) contains two genes coding for ryanodine receptors (RYR1 and RYR2); we also found RYR3 to correlate significantly with helminth diversity (Additional file 3, Table S3). The function of these molecules in the immune system is poorly understood yet, both RYR1 and RYR3 have been involved in calcium signaling in T cells [55, 56]. Moreover, recent evidences have indicated that dendritic cells express RYR1 and activation of the receptor causes a rapid increase in the expression of MHCII molecules on the surface of these cells [57].

Finally, it is interesting to notice that among the genes subjected to helminth-driven selective pressure in network B we found SYK, encoding a tyrosine kinase that interacts with the high affinity IgE receptor and mediates IgE signaling in mast cells and basophils [58, 59]. Similarly, CD226 and FYN (both in network A, Figure 2) have been involved in mast cell activation mediated by the high affinity IgE receptor [60], suggesting a role for these genes in allergic inflammation. Indeed, syk has been shown to mediate airway hyper-responsiveness in an experimental mouse model [61] and most genes discussed above have been involved in the elicitation of allergic phenomena. In addition genes reported in Table 4, the interaction of CD28 with CD86 is central to induction of allergic airway inflammation in mice [62] and CD86 antisense oligonucleotides suppress airway hyper-responsiveness in allergic animals [63]. Variants in C3 have been associated with asthma [64] and mice deficient in C3 exhibit diminished airway hyper-responsiveness and lung eosinophilia when challenged with allergen [65]; also, NOTCH1 is involved CD8+ T cell-mediated development of airway hyper-responsiveness and inflammation [66], while ITGAL/ITGB2 mediates altered responsiveness of atopic asthmatic airway smooth muscle in rabbits [67]. Finally, Ets-1 induces tenascin expression in bronchial fibroblasts [68]. In this respect it is worth mentioning that, although subjects genotyped in the HGDP-CEPH panel are supposed to be healthy, a proportion of them may suffer from relatively mild diseases including asthma, atopy and related disorders; this may be especially true in some areas such as Latin America, for example, where urban centres have the highest reported prevalence of asthma worldwide (reviewed in [69]). While this possibility does not affect the results we reported herein, it highlights the fact that the epidemiology of these disorders is rapidly changing, and several reports have revealed a general increase in prevalence with urbanization, leading to the suggestion that environmental factors (including helminth infections) may play a central role in modulating the susceptibility to these diseases (reviewed in [69]). The relationship between asthma/allergy susceptibility and parasitic worms is though to be complex (reviewed in [70]). On one hand helminth-driven selective pressure is expected to favor individuals carrying alleles that allow a strong Th2 response and, therefore to promote the transmission and spread of asthma-susceptibility variants. On the other hand, lack of parasites in developed countries has likely removed the immunomodulatory role of these organisms, eventually leading to the increased incidence of atopic conditions. The current knowledge of asthma/allergy susceptibility alleles (12 alleles identified by GWAS) is too limited to warrant extensive speculation on the first issue. Still, our data indicate that many genes we identified carry variants associated with asthma/allergy or have been involved in the elicitation of airway hyper-responsiveness. Therefore, our results expand the previously noticed parallelism between genes involved in the development of asthma/allergy and those responsible for responding to parasitic worms, suggesting that the evolutionary scenario underlying the increase in asthma, allergy and related phenotype envisages a relevant role for these long-standing parasites.

Among the genes subjected to helminth-driven selective pressure we identified EDAR and EDA, its ligand. Binding of ectodysplasin to EDAR activates the NF-kB pathway through the NEMO protein. The EDA/EDAR pair mediates signals needed for the development of ectodermal appendages and mutations in both genes result in hypohidrotic ectodermal dysplasia. Many studies [15, 17, 71, 72] have indicated that EDAR has been subjected to a strong selective pressure resulting in the rapid spread of the putatively selected 370Ala allele in Asian populations. This allele is responsible for the hair phenotype of these populations but the selective pressure underlying the selective sweep is unknown. Hypotheses have been proposed that increased hair thickness might be protective against cold climates or be favored through sexual selection [18, 73]. We found that in Asian populations, most chromosomes carrying the selected allele also carry four SNPs subjected to helminth-driven selective pressure (Figure 1). Both EDA and EDAR are expressed in human lymphocytes and dendritic cells (see methods), suggesting that they may function as NF-kB activators in these cell types, as well. It is therefore tempting to speculate that helminths represent the selective pressure underlying the spread of a selected allele in Asia. This idea is consistent with the concept whereby infectious agents have represented one of the major selective forces for human populations.

Conclusions

In summary, our data are consistent with the notion whereby parasitic worms have acted as a powerful selective force on human populations and have contributed to shape nucleotide variability at a number of genes involved in immune responses. We also show that several genes associated with helminth diversity are involved in the pathogenesis of atopic conditions or in airway hyper-responsiveness.

Methods

Data retrieval and statistical analysis

Helminth absence/presence matrices for the 21 countries where HGDP-CEPH populations are located were derived from the Gideon database. Information in Gideon is weekly updated and derives from World Health Organization reports, National Health Ministries, PubMed searches and epidemiology meetings. The Gideon Epidemiology module follows the status of known infectious diseases globally, as well as in individual countries, with specific notes indicating the disease's history, incidence and distribution per country. We manually curated helminth absence/presence matrices by extracting information from single Gideon entries. Following previous suggestions [4–6], we recorded only helminths that are transmitted in the 21 countries, meaning that cases of transmission due to tourism and immigration were not taken into account. A total of 60 helminth species were identified in at least one country (Additional file 6, Table S4). Prevalence data for single helminth infections were similarly obtained from Gideon, as described in the text. In order to calculate parasite diversity for species that are common in at least one country, we inspected Gideon entries for survey data or prevalence notes; helminth infections reported as "rare in humans" were discarded; similarly, parasites with no prevalence estimates or notes were considered as rare; therefore, this diversity measure should be regarded as an approximate estimate.

The annual minimum and maximum temperature were retrieved from the NCEP/NCAR database (Legates and Willmott Average, re-gridded dataset) using the geographic coordinates reported by HGDP-CEPH website for each population. Similarly, net short wave radiation flux data were obtained from NCEP/NCAR (Reanalysis 1: Surface Flux); these data were read using Grid Analysis and Display System (GrADS).

Since helminth diversity, due to data organization in Gideon, can only be calculated per country (rather than per population), the same procedure was applied to climatic variables. Therefore the values of annual temperature, radiation flux and precipitation rate were averaged for populations located in the same country. This assures that a similar number of ties is maintained in all correlation analyses.

Data concerning the HGDP-CEPH panel derive from a previous work [11]. Atypical or duplicated samples and pairs of close relatives were removed [74]. Following previous indications [4, 5], Bantu individuals (South Africa) were considered as one population.

A SNP was ascribed to a specific gene if it was located within the transcribed region or no farther than 500 bp upstream the transcription start site. MAF for any single SNP was calculated as the average over all populations. The list of immune response genes was derived from the Immunology Database and Analysis Portal (ImmPort). Expression data were obtained from SymAtlas. SNPs identified in GWAS and associated with allergy, asthma or related traits (serum IgE levels and plasma eosinophil count) were derived from the A Catalog of Published Genome-Wide Association Studies. The list of allergy/asthma susceptibility genes was obtained from a previous review [9] or from the Online Medelian Inheritance in Man website (MIM: 600807).

All correlations were calculated by Kendall's rank correlation coefficient (τ), a non-parametric statistic used to measure the degree of correspondence between two rankings. The reason for using this test is that even in the presence of ties, the sampling distribution of τ satisfactorily converges to a normal distribution for values of n larger than 10 [75].

In order to estimate the probability of obtaining 246 genes carrying at least one significantly associated SNP out of a group of 2,287 genes (the number of ImmPort genes), we applied a re-sampling approach after dividing genes on the basis of the number of SNPs typed in the HGDP-CEPH panel. In particular, all genes covered by at least one SNP in the HGDP-CEPH panel (number of genes = 15,280) were divided in 24 intervals based on the distribution of typed SNPs per gene (Additional file 7, Table S5). Samples of 2,287 genes were randomly extracted from a list of all genes covered by at least one SNP in the HGDP-CEPH panel by applying the criterion that for each ImmPort gene, a control gene was selected from the same interval. For each sample the number of genes with at least one significant SNP were counted. The empirical probability of obtaining 246 genes was then calculated from the distribution of counts deriving from 10,000 random samples. Similarly, the number of SNPs in ImmPort genes was compared to the distribution of SNPs in the 10,000 re-samplings.

Analysis of PANTHER over-represented functional categories and pathways was performed using the "Compare Classifications of Lists" tool available at the PANTHER classification system website. Briefly, gene lists are compared to the reference list using the binomial test [22] for each molecular function, biological process, or pathway term in PANTHER. All p values were Bonferroni corrected. All calculation were performed in the R environment [76]. For PANTHER analysis we widened the inclusion criteria in that SNPs located within the transcribed region or in the 25 kb upstream the transcription start site were ascribed to the gene.

eQTL data were derived from the eQTL Resource web site held at the University of Chicago (Prtichard Lab).

Network construction

Biological network analysis was performed with Ingenuity Pathways Analysis (IPA) software using an unsupervised analysis. IPA builds networks by querying the Ingenuity Pathways Knowledge Base for interactions between the identified genes and all other gene objects stored in the knowledge base; it then generates networks with a maximum network size of 35 genes/proteins. We used all genes showing at least one significantly associated SNP as the input set; in this case a SNP was ascribed to a gene if it was located within the transcribed region or in the 25 kb upstream. All network edges are supported by at least one published reference or from canonical information stored in the Ingenuity Pathways Knowledge Base. To determine the probability of the analyzed genes to be found together in a network from Ingenuity Pathways Knowledge Base due to random chance alone, IPA applies a Fisher's exact test. The network score represents the -log (p value).

LINKS

The Immunology Database and Analysis Portal, https://www.immport.org

Ingenuity Pathway Analysis, Ingenuity Systems, http://www.ingenuity.com

NCEP/NCAR, Surface flux, http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.surfaceflux.html

GrADS, http://www.iges.org/grads/

SymAtlas,: http://symatlas.gnf.org

Catalog of published GWAS, http://www.genome.gov

Panther,: http://www.pantherdb.org

eQTL resources @ the pritchard lab, http://eqtl.uchicago.edu/Home.html

UCSC Genome Browser, http://genome.ucsc.edu

HGDP-CEPH Panel, http://hagsc.org/hgdp/

Sweep software, http://www.broadinstitute.org/mpg/sweep/

Abbreviations

SNP:

single nucleotide polymorphism

Treg:

regulatory T cell

LD:

linkage disequilibrium

MAF:

minor allele frequency

NIF:

neutrophil inhibitory factor

TRH:

treleasing hormone

PRL:

prolactin.

References

  1. 1.

    Hotez PJ, Brindley PJ, Bethony JM, King CH, Pearce EJ, Jacobson J: Helminth infections: the great neglected tropical diseases. J Clin Invest. 2008, 118 (4): 1311-1321. 10.1172/JCI34261.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  2. 2.

    Quinnell RJ: Genetics of susceptibility to human helminth infection. Int J Parasitol. 2003, 33 (11): 1219-1231. 10.1016/S0020-7519(03)00175-9.

    Article  CAS  PubMed  Google Scholar 

  3. 3.

    Dunne DW, Cooke A: A worm's eye view of the immune system: consequences for evolution of human autoimmune disease. Nat Rev Immunol. 2005, 5 (5): 420-426. 10.1038/nri1601.

    Article  CAS  PubMed  Google Scholar 

  4. 4.

    Fumagalli M, Pozzoli U, Cagliani R, Comi GP, Riva S, Clerici M, Bresolin N, Sironi M: Parasites represent a major selective force for interleukin genes and shape the genetic predisposition to autoimmune conditions. J Exp Med. 2009, 206 (6): 1395-1408. 10.1084/jem.20082779.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  5. 5.

    Fumagalli M, Cagliani R, Pozzoli U, Riva S, Comi GP, Menozzi G, Bresolin N, Sironi M: Widespread balancing selection and pathogen-driven selection at blood group antigen genes. Genome Res. 2009, 19 (2): 199-212. 10.1101/gr.082768.108.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  6. 6.

    Prugnolle F, Manica A, Charpentier M, Guegan JF, Guernier V, Balloux F: Pathogen-driven selection and worldwide HLA class I diversity. Curr Biol. 2005, 15 (11): 1022-1027. 10.1016/j.cub.2005.04.050.

    Article  CAS  PubMed  Google Scholar 

  7. 7.

    Pullan R, Brooker S: The health impact of polyparasitism in humans: are we under-estimating the burden of parasitic diseases?. Parasitology. 2008, 135 (7): 783-794. 10.1017/S0031182008000346.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  8. 8.

    Guernier V, Hochberg ME, Guegan JF: Ecology drives the worldwide distribution of human diseases. PLoS Biol. 2004, 2 (6): e141-10.1371/journal.pbio.0020141.

    PubMed Central  Article  PubMed  Google Scholar 

  9. 9.

    Vercelli D: Discovering susceptibility genes for asthma and allergy. Nat Rev Immunol. 2008, 8 (3): 169-182. 10.1038/nri2257.

    Article  CAS  PubMed  Google Scholar 

  10. 10.

    Clerici M, Butto S, Lukwiya M, Saresella M, Declich S, Trabattoni D, Pastori C, Piconi S, Fracasso C, Fabiani M, Ferrante P, Rizzardini G, Lopalco L: Immune activation in africa is environmentally-driven and is associated with upregulation of CCR5. Italian-Ugandan AIDS Project. AIDS. 2000, 14 (14): 2083-2092. 10.1097/00002030-200009290-00003.

    Article  CAS  PubMed  Google Scholar 

  11. 11.

    Li JZ, Absher DM, Tang H, Southwick AM, Casto AM, Ramachandran S, Cann HM, Barsh GS, Feldman M, Cavalli-Sforza LL, Myers RM: Worldwide human relationships inferred from genome-wide patterns of variation. Science. 2008, 319 (5866): 1100-1104. 10.1126/science.1153717.

    Article  CAS  PubMed  Google Scholar 

  12. 12.

    Coop G, Pickrell JK, Novembre J, Kudaravalli S, Li J, Absher D, Myers RM, Cavalli-Sforza LL, Feldman MW, Pritchard JK: The role of geography in human adaptation. PLoS Genet. 2009, 5 (6): e1000500.-10.1371/journal.pgen.1000500.

    PubMed Central  Article  PubMed  Google Scholar 

  13. 13.

    Handley LJ, Manica A, Goudet J, Balloux F: Going the distance: human population genetics in a clinal world. Trends Genet. 2007, 23 (9): 432-439. 10.1016/j.tig.2007.07.002.

    Article  CAS  PubMed  Google Scholar 

  14. 14.

    Barreiro LB, Laval G, Quach H, Patin E, Quintana-Murci L: Natural selection has driven population differentiation in modern humans. Nat Genet. 2008, 40 (3): 340-345. 10.1038/ng.78.

    Article  CAS  PubMed  Google Scholar 

  15. 15.

    Voight BF, Kudaravalli S, Wen X, Pritchard JK: A map of recent positive selection in the human genome. PLoS Biol. 2006, 4 (3): e72-10.1371/journal.pbio.0040072.

    PubMed Central  Article  PubMed  Google Scholar 

  16. 16.

    International HapMap Consortium, Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal SM, Pasternak S, Wheeler DA, Willis TD, Yu F, Yang H, Zeng C, Gao Y, Hu H, Hu W, Li C, Lin W, Liu S, Pan H, Tang X, Wang J, Wang W, Yu J, Zhang B, Zhang Q, Zhao H, et al: A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007, 449 (7164): 851-861. 10.1038/nature06258.

    Article  Google Scholar 

  17. 17.

    Sabeti PC, Varilly P, Fry B, Lohmueller J, Hostetter E, Cotsapas C, Xie X, Byrne EH, McCarroll SA, Gaudet R, Schaffner SF, Lander ES, International HapMap Consortium, Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal SM, Pasternak S, Wheeler DA, Willis TD, Yu F, Yang H, Zeng C, Gao Y, Hu H, et al: Genome-wide detection and characterization of positive selection in human populations. Nature. 2007, 449 (7164): 913-918. 10.1038/nature06250.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  18. 18.

    Fujimoto A, Kimura R, Ohashi J, Omi K, Yuliwulandari R, Batubara L, Mustofa MS, Samakkarn U, Settheetham-Ishida W, Ishida T, Morishita Y, Furusawa T, Nakazawa M, Ohtsuka R, Tokunaga K: A scan for genetic determinants of human hair morphology: EDAR is associated with Asian hair thickness. Hum Mol Genet. 2008, 17 (6): 835-843. 10.1093/hmg/ddm355.

    Article  CAS  PubMed  Google Scholar 

  19. 19.

    Sabeti PC, Reich DE, Higgins JM, Levine HZ, Richter DJ, Schaffner SF, Gabriel SB, Platko JV, Patterson NJ, McDonald GJ, Ackerman HC, Campbell SJ, Altshuler D, Cooper R, Kwiatkowski D, Ward R, Lander ES: Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002, 419 (6909): 832-837. 10.1038/nature01140.

    Article  CAS  PubMed  Google Scholar 

  20. 20.

    Thomas PD, Kejariwal A, Guo N, Mi H, Campbell MJ, Muruganujan A, Lazareva-Ulitsky B: Applications for protein sequence-function evolution data: mRNA/protein expression analysis and coding SNP scoring tools. Nucleic Acids Res. 2006, 34 (Web Server): W645-50. 10.1093/nar/gkl229.

    PubMed Central  Article  PubMed  Google Scholar 

  21. 21.

    Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, Diemer K, Muruganujan A, Narechania A: PANTHER: a library of protein families and subfamilies indexed by function. Genome Res. 2003, 13 (9): 2129-2141. 10.1101/gr.772403.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  22. 22.

    Cho RJ, Campbell MJ: Transcription, genomes, function. Trends Genet. 2000, 16 (9): 409-415. 10.1016/S0168-9525(00)02065-5.

    Article  CAS  PubMed  Google Scholar 

  23. 23.

    Levite M: Neurotransmitters activate T-cells and elicit crucial functions via neurotransmitter receptors. Curr Opin Pharmacol. 2008, 8 (4): 460-471. 10.1016/j.coph.2008.05.001.

    Article  CAS  PubMed  Google Scholar 

  24. 24.

    Kavelaars A: Regulated expression of alpha-1 adrenergic receptors in the immune system. Brain Behav Immun. 2002, 16 (6): 799-807. 10.1016/S0889-1591(02)00033-8.

    Article  CAS  PubMed  Google Scholar 

  25. 25.

    Heijnen CJ, Rouppe van der Voort C, van de Pol M, Kavelaars A: Cytokines regulate alpha(1)-adrenergic receptor mRNA expression in human monocytic cells and endothelial cells. J Neuroimmunol. 2002, 125 (1-2): 66-72. 10.1016/S0165-5728(02)00034-6.

    Article  CAS  PubMed  Google Scholar 

  26. 26.

    Kelley KW, Weigent DA, Kooijman R: Protein hormones and immunity. Brain Behav Immun. 2007, 21 (4): 384-392. 10.1016/j.bbi.2006.11.010.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  27. 27.

    Wahab MF, Warren KS, Levy RP: Function of the thyroid and the host-parasite relation in murine schistosomiasis mansoni. J Infect Dis. 1971, 124 (2): 161-171.

    Article  CAS  PubMed  Google Scholar 

  28. 28.

    Saule P, Adriaenssens E, Delacre M, Chassande O, Bossu M, Auriault C, Wolowczuk I: Early variations of host thyroxine and interleukin-7 favor Schistosoma mansoni development. J Parasitol. 2002, 88 (5): 849-855.

    Article  CAS  PubMed  Google Scholar 

  29. 29.

    Gudbjartsson DF, Bjornsdottir US, Halapi E, Helgadottir A, Sulem P, Jonsdottir GM, Thorleifsson G, Helgadottir H, Steinthorsdottir V, Stefansson H, Williams C, Hui J, Beilby J, Warrington NM, James A, Palmer LJ, Koppelman GH, Heinzmann A, Krueger M, Boezen HM, Wheatley A, Altmuller J, Shin HD, Uh ST, Cheong HS, Jonsdottir B, Gislason D, Park CS, Rasmussen LM, Porsbjerg C, et al: Sequence variants affecting eosinophil numbers associate with asthma and myocardial infarction. Nat Genet. 2009, 41 (3): 342-347. 10.1038/ng.323.

    Article  CAS  PubMed  Google Scholar 

  30. 30.

    Kabesch M, Tzotcheva I, Carr D, Hofler C, Weiland SK, Fritzsch C, von Mutius E, Martinez FD: A complete screening of the IL4 gene: novel polymorphisms and their association with asthma and IgE in childhood. J Allergy Clin Immunol. 2003, 112 (5): 893-898. 10.1016/j.jaci.2003.08.033.

    Article  CAS  PubMed  Google Scholar 

  31. 31.

    Maizels RM, Balic A, Gomez-Escobar N, Nair M, Taylor MD, Allen JE: Helminth parasites--masters of regulation. Immunol Rev. 2004, 201: 89-116. 10.1111/j.0105-2896.2004.00191.x.

    Article  CAS  PubMed  Google Scholar 

  32. 32.

    Maizels RM, Yazdanbakhsh M: Immune regulation by helminth parasites: cellular and molecular mechanisms. Nat Rev Immunol. 2003, 3 (9): 733-744. 10.1038/nri1183.

    Article  CAS  PubMed  Google Scholar 

  33. 33.

    Gorczynski R, Khatri I, Lee L, Boudakov I: An interaction between CD200 and monoclonal antibody agonists to CD200R2 in development of dendritic cells that preferentially induce populations of CD4+CD25+ T regulatory cells. J Immunol. 2008, 180 (9): 5946-5955.

    Article  CAS  PubMed  Google Scholar 

  34. 34.

    Oh-Hora M, Yamashita M, Hogan PG, Sharma S, Lamperti E, Chung W, Prakriya M, Feske S, Rao A: Dual functions for the endoplasmic reticulum calcium sensors STIM1 and STIM2 in T cell activation and tolerance. Nat Immunol. 2008, 9 (4): 432-443. 10.1038/ni1574.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  35. 35.

    Takahashi T, Tagami T, Yamazaki S, Uede T, Shimizu J, Sakaguchi N, Mak TW, Sakaguchi S: Immunologic self-tolerance maintained by CD25(+)CD4(+) regulatory T cells constitutively expressing cytotoxic T lymphocyte-associated antigen 4. J Exp Med. 2000, 192 (2): 303-310. 10.1084/jem.192.2.303.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  36. 36.

    McCoy K, Camberis M, Gros GL: Protective immunity to nematode infection is induced by CTLA-4 blockade. J Exp Med. 1997, 186 (2): 183-187. 10.1084/jem.186.2.183.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  37. 37.

    Chen L: Co-inhibitory molecules of the B7-CD28 family in the control of T-cell immunity. Nat Rev Immunol. 2004, 4 (5): 336-347. 10.1038/nri1349.

    Article  CAS  PubMed  Google Scholar 

  38. 38.

    King CL, Xianli J, June CH, Abe R, Lee KP: CD28-deficient mice generate an impaired Th2 response to Schistosoma mansoni infection. Eur J Immunol. 1996, 26 (10): 2448-2455. 10.1002/eji.1830261027.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

    Smith P, Walsh CM, Mangan NE, Fallon RE, Sayers JR, McKenzie AN, Fallon PG: Schistosoma mansoni worms induce anergy of T cells via selective up-regulation of programmed death ligand 1 on macrophages. J Immunol. 2004, 173 (2): 1240-1248.

    Article  CAS  PubMed  Google Scholar 

  40. 40.

    Lo SK, Rahman A, Xu N, Zhou MY, Nagpala P, Jaffe HA, Malik AB: Neutrophil inhibitory factor abrogates neutrophil adhesion by blockade of CD11a and CD11b beta(2) integrins. Mol Pharmacol. 1999, 56 (5): 926-932.

    CAS  PubMed  Google Scholar 

  41. 41.

    Moyle M, Foster DL, McGrath DE, Brown SM, Laroche Y, De Meutter J, Stanssens P, Bogowitz CA, Fried VA, Ely JA: A hookworm glycoprotein that inhibits neutrophil function is a ligand of the integrin CD11b/CD18. J Biol Chem. 1994, 269 (13): 10008-10015.

    CAS  PubMed  Google Scholar 

  42. 42.

    Daub J, Loukas A, Pritchard DI, Blaxter M: A survey of genes expressed in adults of the human hookworm, Necator americanus. Parasitology. 2000, 120 (Pt 2): 171-184. 10.1017/S0031182099005375.

    Article  CAS  PubMed  Google Scholar 

  43. 43.

    Giacomin PR, Wang H, Gordon DL, Botto M, Dent LA: Loss of complement activation and leukocyte adherence as Nippostrongylus brasiliensis develops within the murine host. Infect Immun. 2005, 73 (11): 7442-7449. 10.1128/IAI.73.11.7442-7449.2005.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  44. 44.

    Giacomin PR, Gordon DL, Botto M, Daha MR, Sanderson SD, Taylor SM, Dent LA: The role of complement in innate, adaptive and eosinophil-dependent immunity to the nematode Nippostrongylus brasiliensis. Mol Immunol. 2008, 45 (2): 446-455. 10.1016/j.molimm.2007.05.029.

    Article  CAS  PubMed  Google Scholar 

  45. 45.

    Kerepesi LA, Hess JA, Nolan TJ, Schad GA, Abraham D: Complement component C3 is required for protective innate and adaptive immunity to larval strongyloides stercoralis in mice. J Immunol. 2006, 176 (7): 4315-4322.

    Article  CAS  PubMed  Google Scholar 

  46. 46.

    La Flamme AC, MacDonald AS, Huxtable CR, Carroll M, Pearce EJ: Lack of C3 affects Th2 response development and the sequelae of chemotherapy in schistosomiasis. J Immunol. 2003, 170 (1): 470-476.

    Article  CAS  PubMed  Google Scholar 

  47. 47.

    Osborne BA, Minter LM: Notch signalling during peripheral T-cell activation and differentiation. Nat Rev Immunol. 2007, 7 (1): 64-75. 10.1038/nri1998.

    Article  CAS  PubMed  Google Scholar 

  48. 48.

    Asano N, Watanabe T, Kitani A, Fuss IJ, Strober W: Notch1 signaling and regulatory T cell function. J Immunol. 2008, 180 (5): 2796-2804.

    Article  CAS  PubMed  Google Scholar 

  49. 49.

    Samon JB, Champhekar A, Minter LM, Telfer JC, Miele L, Fauq A, Das P, Golde TE, Osborne BA: Notch1 and TGFbeta1 cooperatively regulate Foxp3 expression and the maintenance of peripheral regulatory T cells. Blood. 2008, 112 (5): 1813-1821. 10.1182/blood-2008-03-144980.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  50. 50.

    Blumenthal SG, Aichele G, Wirth T, Czernilofsky AP, Nordheim A, Dittmer J: Regulation of the human interleukin-5 promoter by Ets transcription factors. Ets1 and Ets2, but not Elf-1, cooperate with GATA3 and HTLV-I Tax1. J Biol Chem. 1999, 274 (18): 12910-12916. 10.1074/jbc.274.18.12910.

    Article  CAS  PubMed  Google Scholar 

  51. 51.

    Romano-Spica V, Georgiou P, Suzuki H, Papas TS, Bhat NK: Role of ETS1 in IL-2 gene expression. J Immunol. 1995, 154 (6): 2724-2732.

    CAS  PubMed  Google Scholar 

  52. 52.

    Thomas RS, Tymms MJ, McKinlay LH, Shannon MF, Seth A, Kola I: ETS1, NFkappaB and AP1 synergistically transactivate the human GM-CSF promoter. Oncogene. 1997, 14 (23): 2845-2855. 10.1038/sj.onc.1201125.

    Article  CAS  PubMed  Google Scholar 

  53. 53.

    Jian JL, Zhu CS, Xu ZW, Ouyang WM, Ma DC, Zhang Y, Chen LJ, Yang AG, Jin BQ: Identification and characterization of the CD226 gene promoter. J Biol Chem. 2006, 281 (39): 28731-28736. 10.1074/jbc.M601786200.

    Article  CAS  PubMed  Google Scholar 

  54. 54.

    Grenningloh R, Kang BY, Ho IC: Ets-1, a functional cofactor of T-bet, is essential for Th1 inflammatory responses. J Exp Med. 2005, 201 (4): 615-626. 10.1084/jem.20041330.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  55. 55.

    Schwarzmann N, Kunerth S, Weber K, Mayr GW, Guse AH: Knock-down of the type 3 ryanodine receptor impairs sustained Ca2+ signaling via the T cell receptor/CD3 complex. J Biol Chem. 2002, 277 (52): 50636-50642. 10.1074/jbc.M209061200.

    Article  CAS  PubMed  Google Scholar 

  56. 56.

    Dammermann W, Zhang B, Nebel M, Cordiglieri C, Odoardi F, Kirchberger T, Kawakami N, Dowden J, Schmid F, Dornmair K, Hohenegger M, Flugel A, Guse AH, Potter BV: NAADP-mediated Ca2+ signaling via type 1 ryanodine receptor in T cells revealed by a synthetic NAADP antagonist. Proc Natl Acad Sci USA. 2009, 106 (26): 10678-10683.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  57. 57.

    Vukcevic M, Spagnoli GC, Iezzi G, Zorzato F, Treves S: Ryanodine receptor activation by Ca v 1.2 is involved in dendritic cell major histocompatibility complex class II surface expression. J Biol Chem. 2008, 283 (50): 34913-34922. 10.1074/jbc.M804472200.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  58. 58.

    Shiue L, Green J, Green OM, Karas JL, Morgenstern JP, Ram MK, Taylor MK, Zoller MJ, Zydowsky LD, Bolen JB: Interaction of p72syk with the gamma and beta subunits of the high-affinity receptor for immunoglobulin E, Fc epsilon RI. Mol Cell Biol. 1995, 15 (1): 272-281.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  59. 59.

    Zhang J, Berenstein EH, Evans RL, Siraganian RP: Transfection of Syk protein tyrosine kinase reconstitutes high affinity IgE receptor-mediated degranulation in a Syk-negative variant of rat basophilic leukemia RBL-2H3 cells. J Exp Med. 1996, 184 (1): 71-79. 10.1084/jem.184.1.71.

    Article  CAS  PubMed  Google Scholar 

  60. 60.

    Bachelet I, Munitz A, Mankutad D, Levi-Schaffer F: Mast cell costimulation by CD226/CD112 (DNAM-1/Nectin-2): a novel interface in the allergic process. J Biol Chem. 2006, 281 (37): 27190-27196. 10.1074/jbc.M602359200.

    Article  CAS  PubMed  Google Scholar 

  61. 61.

    Matsubara S, Li G, Takeda K, Loader JE, Pine P, Masuda ES, Miyahara N, Miyahara S, Lucas JJ, Dakhama A, Gelfand EW: Inhibition of spleen tyrosine kinase prevents mast cell activation and airway hyperresponsiveness. Am J Respir Crit Care Med. 2006, 173 (1): 56-63. 10.1164/rccm.200503-361OC.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  62. 62.

    Mathur M, Herrmann K, Qin Y, Gulmen F, Li X, Krimins R, Weinstock J, Elliott D, Bluestone JA, Padrid P: CD28 interactions with either CD80 or CD86 are sufficient to induce allergic airway inflammation in mice. Am J Respir Cell Mol Biol. 1999, 21 (4): 498-509.

    Article  CAS  PubMed  Google Scholar 

  63. 63.

    Crosby JR, Guha M, Tung D, Miller DA, Bender B, Condon TP, York-DeFalco C, Geary RS, Monia BP, Karras JG, Gregory SA: Inhaled CD86 antisense oligonucleotide suppresses pulmonary inflammation and airway hyper-responsiveness in allergic mice. J Pharmacol Exp Ther. 2007, 321 (3): 938-946. 10.1124/jpet.106.119214.

    Article  CAS  PubMed  Google Scholar 

  64. 64.

    Barnes KC, Grant AV, Baltadzhieva D, Zhang S, Berg T, Shao L, Zambelli-Weiner A, Anderson W, Nelsen A, Pillai S, Yarnall DP, Dienger K, Ingersoll RG, Scott AF, Fallin MD, Mathias RA, Beaty TH, Garcia JG, Wills-Karp M: Variants in the gene encoding C3 are associated with asthma and related phenotypes among African Caribbean families. Genes Immun. 2006, 7 (1): 27-35. 10.1038/sj.gene.6364267.

    Article  CAS  PubMed  Google Scholar 

  65. 65.

    Drouin SM, Corry DB, Kildsgaard J, Wetsel RA: Cutting edge: the absence of C3 demonstrates a role for complement in Th2 effector functions in a murine model of pulmonary allergy. J Immunol. 2001, 167 (8): 4141-4145.

    Article  CAS  PubMed  Google Scholar 

  66. 66.

    Guseh JS, Bores SA, Stanger BZ, Zhou Q, Anderson WJ, Melton DA, Rajagopal J: Notch signaling promotes airway mucous metaplasia and inhibits alveolar development. Development. 2009, 136 (10): 1751-1759. 10.1242/dev.029249.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  67. 67.

    Grunstein MM, Hakonarson H, Maskeri N, Kim C, Chuang S: Intrinsic ICAM-1/LFA-1 activation mediates altered responsiveness of atopic asthmatic airway smooth muscle. Am J Physiol Lung Cell Mol Physiol. 2000, 278 (6): L1154-63.

    CAS  PubMed  Google Scholar 

  68. 68.

    Nakamura Y, Esnault S, Maeda T, Kelly EA, Malter JS, Jarjour NN: Ets-1 regulates TNF-alpha-induced matrix metalloproteinase-9 and tenascin expression in primary bronchial fibroblasts. J Immunol. 2004, 172 (3): 1945-1952.

    Article  CAS  PubMed  Google Scholar 

  69. 69.

    Cooper PJ: Interactions between helminth parasites and allergy. Curr Opin Allergy Clin Immunol. 2009, 9 (1): 29-37. 10.1097/ACI.0b013e32831f44a6.

    PubMed Central  Article  PubMed  Google Scholar 

  70. 70.

    Hopkin J: Immune and genetic aspects of asthma, allergy and parasitic worm infections: evolutionary links. Parasite Immunol. 2009, 31 (5): 267-273. 10.1111/j.1365-3024.2009.01104.x.

    Article  CAS  PubMed  Google Scholar 

  71. 71.

    Williamson SH, Hubisz MJ, Clark AG, Payseur BA, Bustamante CD, Nielsen R: Localizing recent adaptive evolution in the human genome. PLoS Genet. 2007, 3 (6): e90-10.1371/journal.pgen.0030090.

    PubMed Central  Article  PubMed  Google Scholar 

  72. 72.

    Tang K, Thornton KR, Stoneking M: A new approach for using genome scans to detect recent positive selection in the human genome. PLoS Biol. 2007, 5 (7): e171-10.1371/journal.pbio.0050171.

    PubMed Central  Article  PubMed  Google Scholar 

  73. 73.

    Bryk J, Hardouin E, Pugach I, Hughes D, Strotmann R, Stoneking M, Myles S: Positive selection in East Asians for an EDAR allele that enhances NF-kappaB activation. PLoS One. 2008, 3 (5): e2209-10.1371/journal.pone.0002209.

    PubMed Central  Article  PubMed  Google Scholar 

  74. 74.

    Rosenberg NA: Standardized subsets of the HGDP-CEPH Human Genome Diversity Cell Line Panel, accounting for atypical and duplicated samples and pairs of close relatives. Ann Hum Genet. 2006, 70 (Pt 6): 841-847. 10.1111/j.1469-1809.2006.00285.x.

    Article  CAS  PubMed  Google Scholar 

  75. 75.

    Salkind NJ: Encyclopedia of measurement and statistics. 2007, Thousand Oaks, CA: Sage Publications

    Google Scholar 

  76. 76.

    R Development Core Team: R: A Language and Environment for Statistical Computing. Vienna, Austria. 2008

    Google Scholar 

  77. 77.

    Kraus DM, Elliott GS, Chute H, Horan T, Pfenninger KH, Sanford SD, Foster S, Scully S, Welcher AA, Holers VM: CSMD1 is a novel multiple domain complement-regulatory protein highly expressed in the central nervous system and epithelial tissues. J Immunol. 2006, 176: 4419-4430.

    Article  CAS  PubMed  Google Scholar 

  78. 78.

    Park CG, Lee SY, Kandala G, Lee SY, Choi Y: A novel gene product that couples TCR signaling to Fas(CD95) expression in activation-induced cell death. Immunity. 1996, 4: 583-591. 10.1016/S1074-7613(00)80484-7.

    Article  CAS  PubMed  Google Scholar 

  79. 79.

    Besecker B, Bao S, Bohacova B, Papp A, Sadee W, Knoell DL: The human zinc transporter SLC39A8 (Zip8) is critical in zinc-mediated cytoprotection in lung epithelia. Am J Physiol Lung Cell Mol Physiol. 2008, 294: 1127-1136. 10.1152/ajplung.00057.2008.

    Article  Google Scholar 

  80. 80.

    Begum NA, Kobayashi M, Moriwaki Y, Matsumoto M, Toyoshima K, Seya T: Mycobacterium bovis BCG cell wall and lipopolysaccharide induce a novel gene, BIGM103, encoding a 7-TM protein: identification of a new protein family having Zn-transporter and Zn-metalloprotease signatures. Genomics. 2002, 80: 630-645. 10.1006/geno.2002.7000.

    Article  CAS  PubMed  Google Scholar 

  81. 81.

    Blanchet F, Cardona A, Letimier FA, Hershfield MS, Acuto O: CD28 costimulatory signal induces protein arginine methylation in T cells. J Exp Med. 2005, 202: 371-377. 10.1084/jem.20050176.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  82. 82.

    Li H, Kim JH, Koh SS, Stallcup MR: Synergistic effects of coactivators GRIP1 and beta-catenin on gene activation: cross-talk between androgen receptor and Wnt signaling pathways. J Biol Chem. 2003, 279: 4212-4220. 10.1074/jbc.M311374200.

    Article  PubMed  Google Scholar 

  83. 83.

    Kwon HJ, Breese EH, Vig-Varga E, Luo Y, Lee Y, Goebl MG, Harrington MA: Tumor necrosis factor alpha induction of NF-kappaB requires the novel coactivator SIMPL. Mol Cell Biol. 2004, 24: 9317-9326. 10.1128/MCB.24.21.9317-9326.2004.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  84. 84.

    van Es MA, van Vught PW, Blauw HM, Franke L, Saris CG, Van den Bosch L, de Jong SW, de Jong V, Baas F, van't Slot R, Lemmens R, Schelhaas HJ, Birve A, Sleegers K, Van Broeckhoven C, Schymick JC, Traynor BJ, Wokke JH, Wijmenga C, Robberecht W, Andersen PM, Veldink JH, Ophoff RA, van den Berg LH: Genetic variation in DPP6 is associated with susceptibility to amyotrophic lateral sclerosis. Nat Genet. 2008, 40: 29-31. 10.1038/ng.2007.52.

    Article  CAS  PubMed  Google Scholar 

  85. 85.

    Wang J, Kudoh J, Takayanagi A, Shimizu N: Novel human BTB/POZ domain-containing zinc finger protein ZNF295 is directly associated with ZFP161. Biochem Biophys Res Commun. 2005, 327: 615-627. 10.1016/j.bbrc.2004.12.048.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

M.S. is part of the Doctorate School in Molecular Medicine, University of Milan.

This study was supported by grants from Istituto Superiore di Sanita' "Programma Nazionale di Ricerca sull' AIDS", the EMPRO and AVIP EC WP6 Projects, the nGIN EC WP7 Project, the Japan Health Science Foundation, 2008 Ricerca Finalizzata [Italian Ministry of Health], 2008 Ricerca Corrente [Italian Ministry of Health], Progetto FIRB RETI: Rete Italiana Chimica Farmaceutica CHEM-PROFARMA-NET [RBPR05NWWC], and Fondazione CARIPLO.

Author information

Author notes

    Affiliations

    1. Scientific Institute IRCCS E. Medea, Bioinformatic Lab, Via don L. Monza 20, 23842, Bosisio Parini (LC), Italy

      Matteo Fumagalli, Uberto Pozzoli, Rachele Cagliani, Nereo Bresolin & Manuela Sironi

    2. Bioengineering Department, Politecnico di Milano, P.zza L. da Vinci, 32, 20133, Milan, Italy

      Matteo Fumagalli

    3. Dino Ferrari Centre, Department of Neurological Sciences, University of Milan, IRCCS Ospedale Maggiore Policlinico, Mangiagalli and Regina Elena Foundation, Via F. Sforza 35, 20100, Milan, Italy

      Giacomo P Comi & Nereo Bresolin

    4. Department of Biomedical Sciences and Technologies LITA Segrate, University of Milan, Via F.lli Cervi 93, 20090, Milan, Italy

      Mario Clerici

    5. Don C. Gnocchi ONLUS Foundation IRCCS, Via Capecelatro 66, 20148, Milan, Italy

      Mario Clerici

    Authors
    1. Matteo Fumagalli
      View author publications

      You can also search for this author in PubMed Google Scholar

    2. Uberto Pozzoli
      View author publications

      You can also search for this author in PubMed Google Scholar

    3. Rachele Cagliani
      View author publications

      You can also search for this author in PubMed Google Scholar

    4. Giacomo P Comi
      View author publications

      You can also search for this author in PubMed Google Scholar

    5. Nereo Bresolin
      View author publications

      You can also search for this author in PubMed Google Scholar

    6. Mario Clerici
      View author publications

      You can also search for this author in PubMed Google Scholar

    7. Manuela Sironi
      View author publications

      You can also search for this author in PubMed Google Scholar

    Corresponding author

    Correspondence to Manuela Sironi.

    Additional information

    Authors' contributions

    MF, UP, RC and MS performed the analyzes, analyzed and interpreted the data; GPC and NB participated in the study coordination; MC and MS conceived the study and wrote the paper.

    All authors read and approved the final manuscript.

    Matteo Fumagalli, Uberto Pozzoli contributed equally to this work.

    Electronic supplementary material

    Additional file 1: Table S1. Populations in the HGDP-CEPH panel and helminth diversity estimates. (XLS 12 KB)

    Additional file 2: Table S2. Genes in the ImmPort list that display at least one SNP significantly associated with helminth diversity. For each gene the SNP showing the strongest correlation is reported. SNPs are ranked according to the value of τ. (XLS 28 KB)

    Additional file 3: Table S3. SNPs significantly associated with helminth diversity. The table reports all SNPs that withstood Bonferroni correction at the genome-wide level (with alpha = 0.5) and displayed a tau percentile rank higher than the 95th among MAF-matched SNPs, as described in the main text. SNPs are ranked according to the value of τ. If the SNP is located within a genic region (or in the 500 upstream nucleotides) the gene symbol is reported. Alternatively, the gene closest to the SNP and its distance (in bp) are indicated. (XLS 445 KB)

    Figure S1

    Additional file 4: . In addition to the two merged networks in the main text, IPA identified three additional networks (A-C) with p < 10-9. Genes are represented as nodes, edges indicate known interactions between proteins (sold lines depicts direct and dashed lines depict indirect interaction). Genes are color coded as follows: green, genes with at least one SNP significantly associated with helminth diversity; gray, genes covered by at least one SNP in the HGDP-CEPH panel; white, genes with no SNPs in the panel. (TIFF 280 KB)

    Figure S2

    Additional file 5: . Analysis of LD in the genomic region encompassing CD28 and CTLA4. SNPs significantly associated with helminth diversity are shown in red while the region covered by CD28 and CTLA4 are shown in blue. The location of DNAse hypersensitive sites in CD4+ T cells is shown in gray while recombination hot-spots are in black. LD plots (r2) are shown for Yoruba (YRI), Europeans (CEU) and Asians (JPT+CHB). The image was generated by using the "add custom track" utility available through the UCSC Genome Browser. (TIFF 2 MB)

    Additional file 6: Table S4. Helminth species/genera transmitted in at least one country and that are common in at least one country. (DOC 18 KB)

    Additional file 7: Table S5. Gene subdivision on the basis of SNP number. Genes were divided in 24 intervals according to the number of SNPs typed in the HGDP-CEPH panel. (PDF 27 KB)

    Authors’ original submitted files for images

    Below are the links to the authors’ original submitted files for images.

    Authors’ original file for figure 1

    Authors’ original file for figure 2

    Authors’ original file for figure 3

    Rights and permissions

    Reprints and Permissions

    About this article

    Cite this article

    Fumagalli, M., Pozzoli, U., Cagliani, R. et al. The landscape of human genes involved in the immune response to parasitic worms. BMC Evol Biol 10, 264 (2010). https://doi.org/10.1186/1471-2148-10-264

    Download citation

    Keywords

    • Treg Cell
    • Ingenuity Pathway Analysis
    • Helminth Infection
    • Immune Response Gene
    • Helminth Species

    BMC Ecology and Evolution

    ISSN: 2730-7182

    Contact us