Incorporating genetics in identifying peanut allergy risk and tailoring allergen immunotherapy: A perspective on the genetic findings from the LEAP trial.

Examining the genetics of peanut allergy (PA) in the context of clinical trial interventions and outcomes provides an opportunity to not only understand gene-environment interactions for PA risk but to also understand the benefit of allergen immunotherapy. A consistent theme in the genetics of food allergy is that in keeping with the dual allergen exposure hypothesis, barrier- and immune-related genes are most commonly implicated in food allergy and tolerance. With a focus on PA, we review how genetic risk factors across 3 genes (FLG, MALT1, and HLA-DQA1) have helped delineate distinct allergic characteristics and outcomes in the context of environmental interventions in the Learning Early about Peanut Allergy (LEAP) study and other clinical trials. We specifically consider and present a framework for genetic risk prediction for the development of PA and discuss how genetics, age, and oral consumption intertwine to predict PA outcome. Although there is some promise in this proposed framework, a better understanding of the mechanistic pathways by which PA develops and persists is needed to develop targeted therapeutics for established disease. Only by understanding the mechanisms by which PA develops, persists, and resolves can we identify adjuvants to oral immunotherapy to make older children and adults immunologically similar to their younger, more malleable counterparts and thus more likely to achieve long-term tolerance.

Epidermal differentiation complex genetic variation in atopic dermatitis and peanut allergy.

BACKGROUND: Deleterious variation in the epidermal differentiation complex (EDC) on chromosome 1 is a well-known genetic determinant of atopic dermatitis (AD) and has been associated with risk of peanut allergy (PA) in population-based studies. OBJECTIVE: Our aim was to determine the effect of genetic variation in the EDC on AD trajectory and risk of PA in early life. METHODS: Genome sequencing was used to measure genetic variation in the EDC in the Learning Early about Peanut Allergy (LEAP) study participants. Association tests were done to identify gene- and variant-level predicted deleterious variation associated with AD severity by using the Scoring Atopic Dermatitis (SCORAD) tool (n = 559) at baseline and each follow-up visit, as well as PA and food allergy in peanut avoiders (n = 275). Predicted deleterious variants included missense variants that were frameshift insertions, frameshift deletions, stop-gain mutations, or stop-loss mutations. Associations between variant load, SCORAD score, and PA were tested by using linear and generalized linear regression models. RESULTS: The genes FLG, FLG2, HRNR, and TCHH1 harbored the most predicted deleterious variation (30, 6, 3, and 1 variant, respectively). FLG variants were associated with SCORAD score at all time points; 4 variants (R1798X, R501X, S126X, and S761fs) drove the association with SCORAD score at each time point, and higher variant load was associated with greater AD severity over time. There was an association between these variants and PA, which remained significant independent of baseline AD severity (odds ratio = 2.63 [95% CI = 1.11-6.01] [P = .02]). CONCLUSIONS: Variation in FLG predicted to be deleterious is associated with AD severity at baseline and longitudinally and has an association with PA independent of baseline severity.

HLA-associated outcomes in peanut oral immunotherapy trials identify mechanistic and clinical determinants of therapeutic success.

Rationale: Previous studies identified an interaction between HLA and oral peanut exposure. HLA-DQA1*01:02 had a protective role with the induction of Ara h 2 epitope-specific IgG4 associated with peanut consumption during the LEAP clinical trial for prevention of peanut allergy, while it was a risk allele for peanut allergy in the peanut avoidance group. We have now evaluated this gene-environment interaction in two subsequent peanut oral immunotherapy (OIT) trials - IMPACT and POISED - to better understand the potential for the HLA-DQA1*01:02 allele as an indicator of higher likelihood of desensitization, sustained unresponsiveness, and peanut allergy remission. Methods: We determined HLA-DQA1*01:02 carrier status using genome sequencing from POISED (N=118, age: 7-55yr) and IMPACT (N=126, age: 12-<48mo). We tested for association with remission, sustained unresponsiveness (SU), and desensitization in the OIT groups, as well as peanut component specific IgG4 (psIgG4) using generalized linear models and adjusting for relevant covariates and ancestry. Results: While not quite statistically significant, a higher proportion of HLA-DQA1*01:02 carriers receiving OIT in IMPACT were desensitized (93%) compared to non-carriers (78%); odds ratio (OR)=5.74 (p=0.06). In this sample we also observed that a higher proportion of carriers achieved remission (35%) compared to non-carriers (22%); OR=1.26 (p=0.80). In POISED, carriers more frequently attained continued desensitization (80% versus 61% among non-carriers; OR=1.28, p=0.86) and achieved SU (52% versus 31%; OR=2.32, p=0.19). psIgG4 associations with HLA-DQA1*01:02 in the OIT arm of IMPACT which included younger study subjects recapitulated patterns noted in LEAP, but no associations of note were observed in the older POISED study subjects. Conclusions: Findings across three clinical trials show a pattern of a gene environment interaction between HLA and oral peanut exposure. Age, and prior sensitization contribute additional determinants of outcomes, consistent with a mechanism of restricted antigen recognition fundamental to driving protective immune responses to OIT.

Secondary analyses for genome-wide association studies using expression quantitative trait loci.

Genome-wide association studies (GWAS) have successfully identified thousands of single nucleotide polymorphisms (SNPs) associated with complex traits; however, the identified SNPs account for a fraction of trait heritability, and identifying the functional elements through which genetic variants exert their effects remains a challenge. Recent evidence suggests that SNPs associated with complex traits are more likely to be expression quantitative trait loci (eQTL). Thus, incorporating eQTL information can potentially improve power to detect causal variants missed by traditional GWAS approaches. Using genomic, transcriptomic, and platelet phenotype data from the Genetic Study of Atherosclerosis Risk family-based study, we investigated the potential to detect novel genomic risk loci by incorporating information from eQTL in the relevant target tissues (i.e., platelets and megakaryocytes) using established statistical principles in a novel way. Permutation analyses were performed to obtain family-wise error rates for eQTL associations, substantially lowering the genome-wide significance threshold for SNP-phenotype associations. In addition to confirming the well known association between PEAR1 and platelet aggregation, our eQTL-focused approach identified a novel locus (rs1354034) and gene (ARHGEF3) not previously identified in a GWAS of platelet aggregation phenotypes. A colocalization analysis showed strong evidence for a functional role of this eQTL.

HLA alleles and sustained peanut consumption promote IgG4 responses in subjects protected from peanut allergy.

We investigated the interplay between genetics and oral peanut protein exposure in the determination of the immunological response to peanut using the targeted intervention in the LEAP clinical trial. We identified an association between peanut-specific IgG4 and HLA-DQA1*01:02 that was only observed in the presence of sustained oral peanut protein exposure. The association between IgG4 and HLA-DQA1*01:02 was driven by IgG4 specific for the Ara h 2 component. Once peanut consumption ceased, the association between IgG4-specific Ara h 2 and HLA-DQA1*01:02 was attenuated. The association was validated by observing expanded IgG4-specific epitopes in people who carried HLA-DQA1*01:02. Notably, we confirmed the previously reported associations with HLA-DQA1*01:02 and peanut allergy risk in the absence of oral peanut protein exposure. Interaction between HLA and presence or absence of exposure to peanut in an allergen- and epitope-specific manner implicates a mechanism of antigen recognition that is fundamental to driving immune responses related to allergy risk or protection.

Genome sequencing unveils a regulatory landscape of platelet reactivity.

Platelet aggregation at the site of atherosclerotic vascular injury is the underlying pathophysiology of myocardial infarction and stroke. To build upon prior GWAS, here we report on 16 loci identified through a whole genome sequencing (WGS) approach in 3,855 NHLBI Trans-Omics for Precision Medicine (TOPMed) participants deeply phenotyped for platelet aggregation. We identify the RGS18 locus, which encodes a myeloerythroid lineage-specific regulator of G-protein signaling that co-localizes with expression quantitative trait loci (eQTL) signatures for RGS18 expression in platelets. Gene-based approaches implicate the SVEP1 gene, a known contributor of coronary artery disease risk. Sentinel variants at RGS18 and PEAR1 are associated with thrombosis risk and increased gastrointestinal bleeding risk, respectively. Our WGS findings add to previously identified GWAS loci, provide insights regarding the mechanism(s) by which genetics may influence cardiovascular disease risk, and underscore the importance of rare variant and regulatory approaches to identifying loci contributing to complex phenotypes.

Gene and protein expression in human megakaryocytes derived from induced pluripotent stem cells.

BACKGROUND: There is interest in deriving megakaryocytes (MKs) from pluripotent stem cells (iPSC) for biological studies. We previously found that genomic structural integrity and genotype concordance is maintained in iPSC-derived MKs. OBJECTIVE: To establish a comprehensive dataset of genes and proteins expressed in iPSC-derived MKs. METHODS: iPSCs were reprogrammed from peripheral blood mononuclear cells (MNCs) and MKs were derived from the iPSCs in 194 healthy European American and African American subjects. mRNA was isolated and gene expression measured by RNA sequencing. Protein expression was measured in 62 of the subjects using mass spectrometry. RESULTS AND CONCLUSIONS: MKs expressed genes and proteins known to be important in MK and platelet function and demonstrated good agreement with previous studies in human MKs derived from CD34+ progenitor cells. The percent of cells expressing the MK markers CD41 and CD42a was consistent in biological replicates, but variable across subjects, suggesting that unidentified subject-specific factors determine differentiation of MKs from iPSCs. Gene and protein sets important in platelet function were associated with increasing expression of CD41/42a, while those related to more basic cellular functions were associated with lower CD41/42a expression. There was differential gene expression by the sex and race (but not age) of the subject. Numerous genes and proteins were highly expressed in MKs but not known to play a role in MK or platelet function; these represent excellent candidates for future study of hematopoiesis, platelet formation, and/or platelet function.

Genome-wide association study of asthma, total IgE, and lung function in a cohort of Peruvian children.

BACKGROUND: Genetic ancestry plays a role in asthma health disparities. OBJECTIVE: Our aim was to evaluate the impact of ancestry on and identify genetic variants associated with asthma, total serum IgE level, and lung function. METHODS: A total of 436 Peruvian children (aged 9-19 years) with asthma and 291 without asthma were genotyped by using the Illumina Multi-Ethnic Global Array. Genome-wide proportions of indigenous ancestry populations from continental America (NAT) and European ancestry from the Iberian populations in Spain (IBS) were estimated by using ADMIXTURE. We assessed the relationship between ancestry and the phenotypes and performed a genome-wide association study. RESULTS: The mean ancestry proportions were 84.7% NAT (case patients, 84.2%; controls, 85.4%) and 15.3% IBS (15.8%; 14.6%). With adjustment for asthma, NAT was associated with higher total serum IgE levels (P < .001) and IBS was associated with lower total serum IgE levels (P < .001). NAT was associated with higher FEV1 percent predicted values (P < .001), whereas IBS was associated with lower FEV1 values in the controls but not in the case patients. The HLA-DR/DQ region on chromosome 6 (Chr6) was strongly associated with total serum IgE (rs3135348; P = 3.438 × 10-10) and was independent of an association with the haplotype HLA-DQA1∼HLA-DQB1:04.01∼04.02 (P = 1.55 × 10-05). For lung function, we identified a locus (rs4410198; P = 5.536 × 10-11) mapping to Chr19, near a cluster of zinc finger interacting genes that colocalizes to the long noncoding RNA CTD-2537I9.5. This novel locus was replicated in an independent sample of pediatric case patients with asthma with similar admixture from Brazil (P = .005). CONCLUSION: This study confirms the role of HLA in atopy, and identifies a novel locus mapping to a long noncoding RNA for lung function that may be specific to children with NAT.

Current insights into the genetics of food allergy.

Food allergy (FA), a growing public health burden in the United States, and familial aggregation studies support strong roles for both genes and environment in FA risk. Deepening our understanding of the molecular and cellular mechanisms driving FAs is paramount to improving its prevention, diagnosis, and clinical management. In this review, we document lessons learned from the genetics of FA that have aided our understanding of these mechanisms. Although current genetic association studies suffer from low power, heterogeneity in definition of FA, and difficulty in our ability to truly disentangle FA from food sensitization (FS) and general atopy genetics, they reveal a set of genetic loci, genes, and variants that continue to implicate the importance of barrier and immune function genes across the atopic march, and FA in particular. The largest reported effects on FA are from MALT1 (odds ratio, 10.99), FLG (average odds ratio, ∼2.9), and HLA (average odds ratio, ∼2.03). The biggest challenge in the field of FA genetics is to elucidate the specific mechanism of action on FA risk and pathogenesis for these loci, and integrative approaches including genetics/genomics with transcriptomics, proteomics, and metabolomics will be critical next steps to translating these genetic insights into practice.

Transcriptional profile of platelets and iPSC-derived megakaryocytes from whole-genome and RNA sequencing.

Genome-wide association studies have identified common variants associated with platelet-related phenotypes, but because these variants are largely intronic or intergenic, their link to platelet biology is unclear. In 290 normal subjects from the GeneSTAR Research Study (110 African Americans [AAs] and 180 European Americans [EAs]), we generated whole-genome sequence data from whole blood and RNA sequence data from extracted nonribosomal RNA from 185 induced pluripotent stem cell-derived megakaryocyte (MK) cell lines (platelet precursor cells) and 290 blood platelet samples from these subjects. Using eigenMT software to select the peak single-nucleotide polymorphism (SNP) for each expressed gene, and meta-analyzing the results of AAs and EAs, we identify (q-value < 0.05) 946 cis-expression quantitative trait loci (eQTLs) in derived MKs and 1830 cis-eQTLs in blood platelets. Among the 57 eQTLs shared between the 2 tissues, the estimated directions of effect are very consistent (98.2% concordance). A high proportion of detected cis-eQTLs (74.9% in MKs and 84.3% in platelets) are unique to MKs and platelets compared with peak-associated SNP-expressed gene pairs of 48 other tissue types that are reported in version V7 of the Genotype-Tissue Expression Project. The locations of our identified eQTLs are significantly enriched for overlap with several annotation tracks highlighting genomic regions with specific functionality in MKs, including MK-specific DNAse hotspots, H3K27-acetylation marks, H3K4-methylation marks, enhancers, and superenhancers. These results offer insights into the regulatory signature of MKs and platelets, with significant overlap in genes expressed, eQTLs detected, and enrichment within known superenhancers relevant to platelet biology.

Advancing Food Allergy Through Omics Sciences.

Since the publication of the first draft of the human genome, there has been an explosion of new technologies with increasing power to interrogate the totality of biological molecules (eg, DNA, RNA, proteins, metabolites) and their modifications (eg, DNA methylation, histone modifications). These technologies, collectively called omics, have been widely applied in the last 2 decades to study biological systems to gain deeper insight into mechanisms driving the physiology and pathophysiology of human health and disease. Because of its complex, multifactorial nature, food allergy is especially well suited to be investigated using omics approaches. In this rostrum, we review how omic technologies have been applied to explore diverse aspects of food allergy, including adaptive and innate immune processes in food-allergic responses, the role of the microbiome in food allergy risk, metabolic changes in the gut and blood associated with food allergy, and the identification of biomarkers and potential therapeutic targets for the condition. We discuss the strengths and limitations of the studies performed thus far and the need to adopt systems biology approaches that integrate data from multiple omics to fully leverage the potential of these technologies to advance food allergy research and care.

Genomic integrity of human induced pluripotent stem cells across nine studies in the NHLBI NextGen program.

Human induced pluripotent stem cell (hiPSC) lines have previously been generated through the NHLBI sponsored NextGen program at nine individual study sites. Here, we examined the structural integrity of 506 hiPSC lines as determined by copy number variations (CNVs). We observed that 149 hiPSC lines acquired 258 CNVs relative to donor DNA. We identified six recurrent regions of CNVs on chromosomes 1, 2, 3, 16 and 20 that overlapped with cancer associated genes. Furthermore, the genes mapping to regions of acquired CNVs show an enrichment in cancer related biological processes (IL6 production) and signaling cascades (JNK cascade & NFκB cascade). The genomic region of instability on chr20 (chr20q11.2) includes transcriptomic signatures for cancer associated genes such as ID1, BCL2L1, TPX2, PDRG1 and HCK. Of these HCK shows statistically significant differential expression between carrier and non-carrier hiPSC lines. Overall, while a low level of genomic instability was observed in the NextGen generated hiPSC lines, the observation of structural instability in regions with known cancer associated genes substantiates the importance of systematic evaluation of genetic variations in hiPSCs before using them as disease/research models.

Tissue-specific impact of FADS cluster variants on FADS1 and FADS2 gene expression.

Omega-6 (n-6) and omega-3 (n-3) long (≥ 20 carbon) chain polyunsaturated fatty acids (LC-PUFAs) play a critical role in human health and disease. Biosynthesis of LC-PUFAs from dietary 18 carbon PUFAs in tissues such as the liver is highly associated with genetic variation within the fatty acid desaturase (FADS) gene cluster, containing FADS1 and FADS2 that encode the rate-limiting desaturation enzymes in the LC-PUFA biosynthesis pathway. However, the molecular mechanisms by which FADS genetic variants affect LC-PUFA biosynthesis, and in which tissues, are unclear. The current study examined associations between common single nucleotide polymorphisms (SNPs) within the FADS gene cluster and FADS1 and FADS2 gene expression in 44 different human tissues (sample sizes ranging 70-361) from the Genotype-Tissue Expression (GTEx) Project. FADS1 and FADS2 expression were detected in all 44 tissues. Significant cis-eQTLs (within 1 megabase of each gene, False Discovery Rate, FDR<0.05, as defined by GTEx) were identified in 12 tissues for FADS1 gene expression and 23 tissues for FADS2 gene expression. Six tissues had significant (FDR< 0.05) eQTLs associated with both FADS1 and FADS2 (including artery, esophagus, heart, muscle, nerve, and thyroid). Interestingly, the identified eQTLs were consistently found to be associated in opposite directions for FADS1 and FADS2 expression. Taken together, findings from this study suggest common SNPs within the FADS gene cluster impact the transcription of FADS1 and FADS2 in numerous tissues and raise important questions about how the inverse expression of these two genes impact intermediate molecular (such a LC-PUFA and LC-PUFA-containing glycerolipid levels) and ultimately clinical phenotypes associated with inflammatory diseases and brain health.

Interferon-γ (IFNG) microsatellite repeat and single nucleotide polymorphism haplotypes of IFN-α receptor (IFNAR1) associated with enhanced malaria susceptibility in Indian populations.

Pro-inflammatory cytokines IFNγ and IFNα function through their cellular receptors IFNγR1 and IFNαR1, respectively to mediate immune processes during malaria infection. A total of 21 SNPs, 2 ins/del polymorphisms and a microsatellite repeat, selected on the basis of their reported association with infectious diseases including malaria in world populations, were analysed for association with Plasmodium falciparum malaria susceptibility in a case-control study with adult patients and ethnically-matched controls drawn from a disease meso- to hyperendemic and a nonendemic region of India. Among the five IFNG SNPs tested, an intron 3 and a 3'UTR SNP associated with disease in the endemic region. In addition, large (CA)n repeats of IFNG intron 1 associated with protection from severe malaria in the endemic region (severe vs. control, odds ratio=0.21, 95% CI=0.08-0.52, P=1.3 × 10(-4)). The TA11CAG haplotype (rs2069705 T/C, rs2430561 A/T, rs3138557 (CA)n, rs2069718 T/C, rs2069727 A/G, rs2069728 G/A) carrying a short CA11 repeat also exhibited very strong association with severe malaria, particularly in the endemic region (severe vs. control, OR=14.56, 95% CI=3.39-85.81, P=3 × 10(-5)). One SNP each from the IFNA8 and IFNA17 of IFNA gene cluster had a protective effect in the non-endemic region but not in the endemic region. A promoter and an intron 2 SNP of IFNAR1 were risk factors for disease and the IFNAR1 haplotype GCCAGG (rs2843710 C/G, rs2850015 C/T, +6993 C/T, rs2243594 A/G, rs1012335 G/C, rs2257167 G/C) carrying both the risk alleles strikingly associated with disease manifestation in the endemic region (severe vs. control, OR=27.14, 95% CI=3.12-1254, P=2 × 10(-5); non-severe vs. control, OR=61.87, 95% CI=10.08-2521, P=1 × 10(-8)). The data indicates dissimilar contribution of cytokine and cytokine receptor variants to disease in populations residing in areas of differential malaria endemicity.

Deletion of the APOBEC3B gene strongly impacts susceptibility to falciparum malaria.

APOBEC3B, a gene involved in innate response, exhibits insertion-deletion polymorphism across world populations. We observed the insertion allele to be nearly fixed in malaria endemic regions of sub-Saharan Africa as well as populations with high malaria incidence in the past. This prompted us to investigate the possible association of the polymorphism with falciparum malaria. We studied the distribution of APOBEC3B, in 25 diverse Indian populations comprising of 500 samples and 176 severe or non-severe Plasmodium falciparum patients and 174 ethnically-matched uninfected individuals from a P. falciparum endemic and a non-endemic region of India. The deletion frequencies ranged from 0% to 43% in the Indian populations. The frequency of the insertion allele strikingly correlated with the endemicity map of P. falciparum malaria in India. A strong association of the deletion allele with susceptibility to falciparum malaria in the endemic region (non-severe vs. control, Odds ratio=4.96, P value=9.5E(-06); severe vs. control, OR=4.36, P value=5.76E(-05)) was observed. Although the frequency of deletion allele was higher in the non-endemic region, there was a significant association of the homozygous deletion genotype with malaria (OR=3.17, 95% CI=1.10-10.32, P value=0.0177). Our study also presents a case for malaria as a positive selection force for the APOBEC3B insertion and suggests a major role for this gene in innate immunity against malaria.

Negative epistasis between α+ thalassaemia and sickle cell trait can explain interpopulation variation in South Asia.

Recent studies in Kenya and Ghana have shown that individuals who inherit two malaria-protective genetic disorders of haemoglobin-α(+) thalassaemia and sickle cell trait-experience a much lower level of malaria protection than those who inherit sickle cell trait alone. We have previously demonstrated that this can limit the frequency of α(+) thalassaemia in a population in which sickle cell is present, which may account for the frequency of α(+) thalassaemia in sub-Saharan Africa not exceeding 50%. Here we consider the relationship between α(+) thalassaemia and sickle cell in South Asian populations, and show that very high levels of α(+) thalassaemia combined with varying levels of malaria selection can explain why sickle cell has penetrated certain South Asian populations but not others.

Distinct cytokine profiles define clinical immune response to falciparum malaria in regions of high or low disease transmission.

The immune effector response to Plasmodium falciparum infection involves a finely-tuned interplay between different cell types and cytokines. However, the processes by which they mediate the development of clinical immunity, in areas of different endemicity, are poorly understood. We analyzed circulating levels of pro-inflammatory (TNF, IFN-γ, IL-12, IL-16) and anti-inflammatory (IL-4, IL-10, IL-13) cytokines in control and patient groups drawn from a P. falciparum-endemic and a non-endemic region of India. The endemic region control population exhibited a lower pro- to anti-inflammatory cytokine ratio, indicating a shift towards a high basal Th2 response. Levels of IL-10 contributed most towards the region-specific difference in basal cytokine response. IL-10 was also the strongest predictor of disease in the endemic region, while IL-12, along with IL-10 and IL-6, contributed most to disease outcome in the non-endemic region. A low, mean IFN-γ/IL-10 ratio was associated with disease severity in the endemic region (p < 0.0001). In contrast, a low mean IL-12/IL-10 ratio correlated with disease outcome in the non-endemic region (p < 0.0001). In the endemic region, IL-13 correlated negatively with IFN-γ in severe patients (Spearman's ρ: -0.49; p : 0.013), while in the non-endemic region, IL-13 correlated negatively with IL-6 in severe malaria patients (Spearman's ρ: -0.485; p : 0.001). In conclusion, levels of pro- and anti-inflammatory cytokines and the relative balance between the Th1 and Th2 response, illustrates how populations residing in areas of varying disease endemicity may respond to P. falciparum-induced immune challenge.

Variations in host genes encoding adhesion molecules and susceptibility to falciparum malaria in India.

BACKGROUND: Host adhesion molecules play a significant role in the pathogenesis of Plasmodium falciparum malaria and changes in their structure or levels in individuals can influence the outcome of infection. The aim of this study was to investigate the association of SNPs of three adhesion molecule genes, ICAM1, PECAM1 and CD36, with severity of falciparum malaria in a malaria-endemic and a non-endemic region of India. METHODS: The frequency distribution of seven selected SNPs of ICAM1, PECAM1 and CD36 was determined in 552 individuals drawn from 24 populations across India. SNP-disease association was analysed in a case-control study format. Genotyping of the population panel was performed by Sequenom mass spectroscopy and patient/control samples were genotyped by SNaPshot method. Haplotypes and linkage disequilibrium (LD) plots were generated using PHASE and Haploview, respectively. Odds-ratio (OR) for risk assessment was estimated using EpiInfotrade mark version 3.4. RESULTS: Association of the ICAM1 rs5498 (exon 6) G allele and the CD36 exon 1a A allele with increased risk of severe malaria was observed (severe versus control, OR = 1.91 and 2.66, P = 0.02 and 0.0012, respectively). The CD36 rs1334512 (-53) T allele as well as the TT genotype associated with protection from severe disease (severe versus control, TT versus GG, OR = 0.37, P = 0.004). Interestingly, a SNP of the PECAM1 gene (rs668, exon 3, C/G) with low minor allele frequency in populations of the endemic region compared to the non-endemic region exhibited differential association with disease in these regions; the G allele was a risk factor for malaria in the endemic region, but exhibited significant association with protection from disease in the non-endemic region. CONCLUSION: The data highlights the significance of variations in the ICAM1, PECAM1 and CD36 genes in the manifestation of falciparum malaria in India. The PECAM1 exon 3 SNP exhibits altered association with disease in the endemic and non-endemic region.