TLR9 polymorphisms in African populations: no association with severe malaria, but evidence of cis-variants acting on gene expression
© Campino et al; licensee BioMed Central Ltd. 2009
Received: 14 August 2008
Accepted: 13 March 2009
Published: 13 March 2009
During malaria infection the Toll-like receptor 9 (TLR9) is activated through induction with plasmodium DNA or another malaria motif not yet identified. Although TLR9 activation by malaria parasites is well reported, the implication to the susceptibility to severe malaria is not clear. The aim of this study was to assess the contribution of genetic variation at TLR9 to severe malaria.
This study explores the contribution of TLR9 genetic variants to severe malaria using two approaches. First, an association study of four common single nucleotide polymorphisms was performed on both family- and population-based studies from Malawian and Gambian populations (n>6000 individual). Subsequently, it was assessed whether TLR9 expression is affected by cis-acting variants and if these variants could be mapped. For this work, an allele specific expression (ASE) assay on a panel of HapMap cell lines was carried out.
No convincing association was found with polymorphisms in TLR9 for malaria severity, in either Gambian or Malawian populations, using both case-control and family based study designs. Using an allele specific expression assay it was observed that TLR9 expression is affected by cis-acting variants, these results were replicated in a second experiment using biological replicates.
By using the largest cohorts analysed to date, as well as a standardized phenotype definition and study design, no association of TLR9 genetic variants with severe malaria was found. This analysis considered all common variants in the region, but it is remains possible that there are rare variants with association signals. This report also shows that TLR9 expression is potentially modulated through cis-regulatory variants, which may lead to differential inflammatory responses to infection between individuals.
Exploring the idea that malaria motifs can induce the innate immune system and initiate an inflammatory response may lead to a greater understanding of resistance/susceptibility to severe malaria. Several reports have found that malaria parasites express molecules that are recognized by Toll-like receptors (TLR), which play a critical role in the early innate immune response to invading pathogens. For instance, the malaria pigment haemozoin has been shown to function as carrier for plasmodium DNA which stimulates TLR9 and induces inflammation [1–3]. Activation of TLR9 expressed on dendritic cell has also been proposed as a mechanism used by malaria parasites to trigger regulatory T cells and evade the immune system . Regulatory T cells were also shown to contribute to the pathogenesis of severe malaria by suppressing anti-malarial immunity during Plasmodium berghei ANKA infection . The correct balance between induction of an immune response against the malaria parasite and limitation of suppression of host immunity by the parasite seems essential to control infection.
To assess the clinical relevance of TLR9 polymorphisms and severe malaria two association studies were conducted in Ghana. While one of the studies did not observe any clear association with severe malaria, another found TLR9 polymorphisms to be associated with the clinical manifestation of malaria during pregnancy . Recently, the TLR-9-1486C/T variants was found to be associated with high parasitaemia in a cohort of patients with mild malaria from the Amazonian region of Brazil .
This report explores the contribution of TLR9 genetic variants to severe malaria using two approaches. At first, a large case-control and family-based association studies was performed in Malawian and Gambian populations using four common polymorphisms located in the TLR9 gene region. It is also known that distant SNPs can affect gene expression, and these can only be identified using functional in vitro assays [9–11]. Therefore, subsequently it was assessed whether TLR9 expression is affected by cis-acting variants and whether these variants could be mapped. Such polymorphism can be further investigated as candidates for possible association with severe malaria. Exploiting the fact that TLR9 is a major TLR expressed by human B cells , an allele specific expression (ASE) assay was performed using HapMap  B-lymphoblastoid cell lines (LCLs), for which genotypes on hundreds-of-thousands SNPs exist, permitting the correlation between expression data and genetic variability [11, 14–17] and thus facilitating the discovery of putative functional polymorphisms.
Patient samples were collected as part of ongoing epidemiological studies of severe malaria at the Royal Victoria Hospital, Banjul, The Gambia (885 severe malaria cases, 629 controls; 1,044 family trios) and the Queen Elizabeth Central Hospital, Blantyre, Malawi (712 severe malaria cases, 416 controls; 248 family trios). The set of nuclear family trios comprise a severe malaria affected child and its two (biological) parents, and were assessed as 'true trios' using the Nuclear software package. All DNA samples were collected and genotyped following approval from the relevant research ethics committees and informed consent from participants.
All cases were children admitted to hospital with evidence of Plasmodium falciparum on blood film and clinical features of severe malaria . Subjects were defined as having had 'cerebral malaria' if their Blantyre coma score was less than or equal to 3 on presentation or early during admission. A second phenotypic subset of 'severe malarial anaemia' was defined as those subjects having had a haemoglobin density of less than 5 g/dl or a haematocrit less than 15%. Participants with co-existing severe or chronic medical conditions (e.g. bacterial pneumonia, kwashiorkor) unrelated to a severe malarial infection were excluded. Controls were cord blood samples obtained from birth clinics in the same locality as the cases, and thought to approximate a random sample of the population thus, reflecting the true allele frequency. Cerebral malaria (CM) and severe malarial anaemia (SMA) were the most common symptoms defining severe malaria.
Sample preparation and genotyping
For the association studies (participants described below) genomic DNA samples underwent whole genome amplification through either Primer Extension Pre-amplification  or Multiple Displacement Amplification , before genotyping on a Sequenom MassArray genotyping platform. Only the four common SNPs (minor allele frequency (MAF) of at least 5% in the HapMap Yoruba population) in the TLR9 region were genotyped. The majority of polymorphisms in this region are monomorphic.
For the ASE assay 19 unrelated CEU and 20 YRI individuals were used, selected from the HapMap collection. Established lymphoblastoid cell lines (LCLs) from these individuals were obtained from the Centre d'Etude du Polymorphisme Humain (CEPH) collection (Coriell Institute for Medical Research). Cell lines were cultured has described previously .
For expression analysis, total RNA was extracted from cells lines from each individual using Tri reagent (Sigma). Poly A+ RNA was isolated from total RNA using the Dynabeads mRNA Purification Kit (Dynal). Synthesis of First-Strand cDNA was processed according to the StrataScript™ First-Strand Synthesis System (Stratagene). The Allelotype platform from Massarray (Sequenom) was utilized for accurate relative quantification of allele specific cDNA and gDNA transcripts . Primers were designed using the dedicated software Spectrodesigner (Sequenom). Allele specific expression quantification was performed according to the protocol of Forton et al  and Campino et a l . Briefly, for each cDNA and gDNA assays, nine technical replicates were performed. Biological replicates from independent cultures were assayed equally. For a given cDNA or gDNA assay the allelic transcript ratio was calculated on each of the nine technical replicates from the relative quantity of the two allele-specific transcripts. The mean allelic transcript ratio for the whole assay was then calculated, and normalized to the mean allelic transcript ratio for the genomic controls. An assay was accepted for further analysis if the standard error of the mean for the technical and biological replicates was less than 10%.
Haplotypes were downloaded from the HapMap Phase II database http://www.hapmap.org. For data analysis, 770,394 markers were used; all monomorphic markers were excluded from the analysis. Locations of the SNP markers were based on those of the human reference sequence http://genome.ucsc.edu of May 2004 (hg 17, build 35).
Statistical analysis of the case-control and family trio association studies
Genotypic deviations from Hardy-Weinberg equilibrium (HWE) were assessed using a chi-square statistical test. Case-control association analysis was undertaken by logistic regression and included the covariates: ethnic group, gender and the HbS polymorphism. Family-based association analysis was performed using the transmission disequilibrium test (TDT) , as well as a case-pseudo-control approach with conditional logistic regression . The TDT test is robust to the effects of population stratification or sub-structure. All analysis was performed using the R statistical package.
Statistical analysis of allelic expression
For the association analysis, log-transformed ratios of allele expression values for individuals with the heterozygous genotype at the exonic SNP (rs352140) were used. To map putative cis-variants a statistical analysis approach was used, based on the linear regression method (similar to that developed by Teare et al and Campino et al)[24, 27]. In particular, it was assessed whether there is evidence of an increasing or decreasing (linear) relationship between the ratio of expressions for the phase-known heterozygous type 1 (e.g. AG), homozygous types (e.g. AA and GG), and heterozygous type 2 (e.g. GA) genotypes. The rationale for the approach is that allelic imbalances would be observed only in individuals heterozygous for the cis-acting polymorphism. All analysis was performed using the R statistical package.
Results and discussion
Association analyses of common TLR9 SNPs and severe malaria
To assess the contribution TLR9 polymorphisms to severe malaria susceptibility/resistance an association analysis was performed using case-control and trio cohorts from Malawian (712 cases, 416 controls; 248 family trios) and The Gambian (885 cases, 629 controls; 1,044 family trios). Cerebral malaria and severe anaemia were the most common symptoms defining severe malaria. The majority of the Gambian (70%) and Malawian (91%) cohorts comprise cases of cerebral malaria.
Allelic frequencies of TLR9 polymorphisms
Severe malaria and TLR9 allelic-based association analysis
C vs T
C vs T
A vs G
T vs C
C vs T
C vs T
A vs G
T vs C
Severe malaria and TLR9 genotype-based association analysis
CT vs TT
CC v TT
CT vs TT
CC vs TT
AG vs GG
AA vs GG
TC vs CC
TT vs CC
The analysis of the Gambia case-control cohorts showed that there was no evidence of association between severe malaria and the polymorphisms studied. Similarly for the Malawi case-controls, there was no strong evidence of association, but there are some marginal effects consistent with the haplotype and allelic analyses above. Accounting for multiple testing would make us less likely to believe the marginal evidence (P ≈ 0.05) of associations for some of the polymorphisms and potential haplotype effects. Mantel-Haenszel pooled estimates of risk for each SNP across the populations and studies did not lead to significant results (p > 0.05). In conclusion, no convincing association was found with polymorphisms at TLR9 for malaria severity, in either Gambian or Malawian populations, using both case-control and family based study designs
Exploring the existence of functional polymorphism affecting TLR9 allelic expression
Candidate gene genetic studies are generally restricted to SNPs covering the gene area as potential regulatory SNPs localized away from a gene are difficult to identify. Therefore, as no association was found with SNPs in the TLR9 gene area, there was further exploration into whether other SNPs which can be far away from the gene could have an important role on the regulation of TLR9 expression.
Numerous studies have recently contributed to the identification of putative regulatory variants by correlating total gene expression or allele specific expression data with genetic variability[11, 16, 17, 28] Using a similar approach, in this report an allele specific expression (ASE) assay was performed on a panel of HapMap B-lymphoblastoid cell lines (LCL) for which genotyped data is available and facilitates the mapping of functional polymorphisms. These cell lines have been widely used to map cis/trans regulatory regions and recent publications have shown high consistence between cis-regulation results obtained using HapMap cell lines and primary tissue.
The strength of the ASE assay is the possibility to compare the relative expression of the two alleles in heterozygous individuals within the same cellular sample. Therefore, both alleles will be exposed to the same environmental, technical and genetic factors (e.g. trans-acting effect). If the expression of the alleles is not equal (allelic imbalance), it suggests that the expression of the gene is under cis- regulation.
Cell lines from individuals of northern and western Europe origin (CEU) and from Yoruba of Ibadan, Nigeria (YRI) were selected from the HapMap collection  and were genotyped for a exonic synonymous SNP present on the TLR9 gene. Since this SNP is located within a part of the gene that is transcribed, it is possible to discriminate between allelic transcripts in heterozygous individuals.
Twenty-one unrelated CEU and nineteen YRI individuals were found to be heterozygous for this SNP and were used for the ASE analysis. For these individuals, their genomic DNA was used to obtain the reference equi-molar ratio between the two alleles (allele1 and allele2), and this was then compared with their allelic RNA expression. Allelic expression imbalance (AEI) was determined if independent replicate assays in the RNA samples showed differential allele expression ratios that deviated from the genomic DNA. If there is no allelic imbalance then the ratio of the two alleles (allele1/allele2) will be equal to one. A ratio greater than one indicates an over-expression of allele1, while a ratio smaller than one indicates an over-expression of allele2.
Overall, while no functional polymorphism were mapped using the current HapMap genotyping data, this data shows that TLR9 expression is potentially modulated through cis-regulatory variants.
To determine if genetic variation in the TLR9 gene influences disease outcome in malaria, a genetic association study was conducted to investigate the possible effect of TLR9 polymorphisms on severe malaria in two African populations. This study uses a standardized phenotype definition and study design, is the largest of its kind to date for TLR9, and includes both case-control and family-based data. No association was found with TLR9 polymorphism and severe manifestation of malaria in the Gambian cohorts and only a weak effect was observed in the Malawi case-control study. In the latter, it was not possible to adjust the potentially confounding effect of ethnicity, and possible population structure. When the more robust TDT analysis was applied to the Malawian trios, no association was observed. This analysis considered all common variants in the TLR9 region, but it remains possible that there are rare variants with association signals.
It is also possible that there exist TLR9 regulatory variants that contribute to severe malaria resistance/susceptibility that have not been mapped. It is known that regulatory polymorphisms which control gene expression can be localized hundreds of kb from the gene they influence[10, 15]. Selection of candidate SNPs for association analysis close to the gene of interest is therefore likely to miss distant functional polymorphisms. One approach that will enrich potential functional polymorphism in association analysis is to use expression studies to identify SNPs that modulate gene expression. In an attempt to find putative cis-variants affecting TLR9 expression, an allele-specific expression assays was carried out in LCL from both CEU and YRI HapMap individuals. Interestingly, allelic expression imbalances were observed in both sample sets, which is reproducible in the same individuals, indicating strongly the existence of cis-regulation. This observation indicates clearly that different individuals might differ in the expression of TLR9, which can lead to different responses to the stimulus induced by malaria motifs or other pathogens. The identification of these functional genetic variants will be crucial to attest any involvement of TLR9 in severe malaria. The finding that TLR9 expression is under cis-regulation will strongly pave the path for further studies of infectious diseases, including the consequences of differential TLR9 activation.
Although evidence of cis-regulation of TLR9 expression was found, it was not possible to identify polymorphisms that correlated with the observed distortion. This is likely to be due to the incomplete polymorphism coverage of the HapMap and by the small sample sizes of the current available LCL. Increasing the number of available cell lines with abundant genotyping data or the application of new re-sequencing technologies to discover new and rare polymorphism will greatly facilitate the identification of functional variants. This work illustrates the need to upscale approaches and resources in expression analysis in the search for disease associated variants.
This work was in part funded by the Wellcome Trust and Medical Research Council UK. SC was funded by a Marie-Curie intra-European fellowship (FP6). JF and AF were funded by a Wellcome Trust Training Fellowship. TGC, SA, KR and DPK are in part funded by the Bill and Melinda Gates Foundation, Wellcome Trust, and Medical Research Council UK.
- Coban C, Ishii KJ, Kawai T, Hemmi H, Sato S, Uematsu S, Yamamoto M, Takeuchi O, Itagaki S, Kumar N, Horii T, Akira S: Toll-like receptor 9 mediates innate immune activation by the malaria pigment hemozoin. J Exp Med. 2005, 201: 19-25. 10.1084/jem.20041836.PubMed CentralView ArticlePubMedGoogle Scholar
- Pichyangkul S, Yongvanitchit K, Kum-arb U, Hemmi H, Akira S, Krieg AM, Heppner DG, Stewart VA, Hasegawa H, Looareesuwan S, Shanks GD, Miller RS: Malaria blood stage parasites activate human plasmacytoid dendritic cells and murine dendritic cells through a Toll-like receptor 9-dependent pathway. J Immunol. 2004, 172: 4926-4933.View ArticlePubMedGoogle Scholar
- Parroche P, Lauw FN, Goutagny N, Latz E, Monks BG, Visintin A, Halmen KA, Lamphier M, Olivier M, Bartholomeu DC, Gazzinelli RT, Golenbock DT: Malaria hemozoin is immunologically inert but radically enhances innate responses by presenting malaria DNA to Toll-like receptor 9. Proc Natl Acad Sci USA. 2007, 104: 1919-1924. 10.1073/pnas.0608745104.PubMed CentralView ArticlePubMedGoogle Scholar
- Hisaeda H, Tetsutani K, Imai T, Moriya C, Tu L, Hamano S, Duan X, Chou B, Ishida H, Aramaki A, Shen J, Ishii KJ, Coban C, Akira S, Takeda K, Yasutomo K, Torii M, Himeno K: Malaria Parasites Require TLR9 Signaling for Immune Evasion by Activating Regulatory T Cells. J Immunol. 2008, 180: 2496-2503.View ArticlePubMedGoogle Scholar
- Amante FH, Stanley AC, Randall LM, Zhou Y, Haque A, McSweeney K, Waters AP, Janse CJ, Good MF, Hill GR, Engwerda CR: A role for natural regulatory T cells in the pathogenesis of experimental cerebral malaria. Am J Pathol. 2007, 171: 548-559. 10.2353/ajpath.2007.061033.PubMed CentralView ArticlePubMedGoogle Scholar
- Mockenhaupt FP, Cramer JP, Hamann L, Stegemann MS, Eckert J, Oh NR, Otchwemah RN, Dietz E, Ehrhardt S, Schroder NW, Bienzle U, Schumann RR: Toll-like receptor (TLR) polymorphisms in African children: common TLR-4 variants predispose to severe malaria. J Commun Dis. 2006, 38: 230-245.PubMedGoogle Scholar
- Mockenhaupt FP, Hamann L, von Gaertner C, Bedu-Addo G, von Kleinsorgen C, Schumann RR, Bienzle U: Common polymorphisms of toll-like receptors 4 and 9 are associated with the clinical manifestation of malaria during pregnancy. J Infect Dis. 2006, 194: 184-188. 10.1086/505152.View ArticlePubMedGoogle Scholar
- Leoratti FM, Farias L, Alves FP, Suarez-Mutis MC, Coura JR, Kalil J, Camargo EP, Moraes SL, Ramasawmy R: Variants in the Toll-Like Receptor Signaling Pathway and Clinical Outcomes of Malaria. J Infect Dis. 2008Google Scholar
- Pastinen T, Hudson TJ: Cis-acting regulatory variation in the human genome. Science. 2004, 306: 647-650. 10.1126/science.1101659.View ArticlePubMedGoogle Scholar
- Stranger BE, Nica AC, Forrest MS, Dimas A, Bird CP, Beazley C, Ingle CE, Dunning M, Flicek P, Koller D, Montgomery S, Tavare S, Deloukas P, Dermitzakis ET: Population genomics of human gene expression. Nat Genet. 2007, 39: 1217-1224. 10.1038/ng2142.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheung VG, Spielman RS, Ewens KG, Weber TM, Morley M, Burdick JT: Mapping determinants of human gene expression by regional and genome-wide association. Nature. 2005, 437: 1365-1369. 10.1038/nature04244.PubMed CentralView ArticlePubMedGoogle Scholar
- Bernasconi NL, Onai N, Lanzavecchia A: A role for Toll-like receptors in acquired immunity: up-regulation of TLR9 by BCR triggering in naive B cells and constitutive expression in memory B cells. Blood. 2003, 101: 4500-4504. 10.1182/blood-2002-11-3569.View ArticlePubMedGoogle Scholar
- The International HapMap Project. Nature. 2003, 426: 789-796. 10.1038/nature02168.Google Scholar
- Hull J, Campino S, Rowlands K, Chan MS, Copley RR, Taylor MS, Rockett K, Elvidge G, Keating B, Knight J, Kwiatkowski D: Identification of common genetic variation that modulates alternative splicing. PLoS Genet. 2007, 3: e99-10.1371/journal.pgen.0030099.PubMed CentralView ArticlePubMedGoogle Scholar
- Forton JT, Udalova IA, Campino S, Rockett KA, Hull J, Kwiatkowski DP: Localization of a long-range cis-regulatory element of IL13 by allelic transcript ratio mapping. Genome Res. 2007, 17: 82-87. 10.1101/gr.5663007.PubMed CentralView ArticlePubMedGoogle Scholar
- Stranger BE, Forrest MS, Clark AG, Minichiello MJ, Deutsch S, Lyle R, Hunt S, Kahl B, Antonarakis SE, Tavare S, Deloukas P, Dermitzakis ET: Genome-wide associations of gene expression variation in humans. PLoS Genet. 2005, 1: e78-10.1371/journal.pgen.0010078.PubMed CentralView ArticlePubMedGoogle Scholar
- Pastinen T, Ge B, Gurd S, Gaudin T, Dore C, Lemire M, Lepage P, Harmsen E, Hudson TJ: Mapping common regulatory variants to human haplotypes. Hum Mol Genet. 2005, 14: 3963-3971. 10.1093/hmg/ddi420.View ArticlePubMedGoogle Scholar
- Teo YYFA, Sanjoaquin MA, Pederson B, Small KS, Rockett KA, Kwiatkowski DP, Clark TG: Assessing genuine parents-offspring trios for genetic association studies. Hum Heredity. 2008Google Scholar
- Marsh K, Forster D, Waruiru C, Mwangi I, Winstanley M, Marsh V, Newton C, Winstanley P, Warn P, Peshu N: Indicators of life-threatening malaria in African children. N Engl J Med. 1995, 332: 1399-1404. 10.1056/NEJM199505253322102.View ArticlePubMedGoogle Scholar
- Zhang L, Cui X, Schmitt K, Hubert R, Navidi W, Arnheim N: Whole genome amplification from a single cell: implications for genetic analysis. Proc Natl Acad Sci USA. 1992, 89: 5847-5851. 10.1073/pnas.89.13.5847.PubMed CentralView ArticlePubMedGoogle Scholar
- Gonzalez JM, Portillo MC, Saiz-Jimenez C: Multiple displacement amplification as a pre-polymerase chain reaction (pre-PCR) to process difficult to amplify samples and low copy number sequences from natural environments. Environ Microbiol. 2005, 7: 1024-1028. 10.1111/j.1462-2920.2005.00779.x.View ArticlePubMedGoogle Scholar
- Knight JC, Keating BJ, Kwiatkowski DP: Allele-specific repression of lymphotoxin-alpha by activated B cell factor-1. Nat Genet. 2004, 36: 394-399. 10.1038/ng1331.View ArticlePubMedGoogle Scholar
- Elvidge GP, Price TS, Glenny L, Ragoussis J: Development and evaluation of real competitive PCR for high-throughput quantitative applications. Anal Biochem. 2005, 339: 231-241. 10.1016/j.ab.2005.01.040.View ArticlePubMedGoogle Scholar
- Campino S, Forton J, Raj S, Mohr A, Auburn S, Fry A, Mangano V, Vandiedonck C, Richardson A, Rockett K, Clark T, Kwiatkowski D: Validating discovered cis-acting regulatory genetic variants: application of an allele specific expression approach to HapMap populations. PLoS ONE. 2008, 3 (12):Google Scholar
- Spielman RS, McGinnis RE, Ewens WJ: Transmission test for linkage disequilibrium: the insulin gene region and insulin-dependent diabetes mellitus (IDDM). Am J Hum Genet. 1993, 52: 506-516.PubMed CentralPubMedGoogle Scholar
- Cordell HJ: Properties of case/pseudocontrol analysis for genetic association studies: Effects of recombination, ascertainment, and multiple affected offspring. Genet Epidemiol. 2004, 26: 186-205. 10.1002/gepi.10306.View ArticlePubMedGoogle Scholar
- Teare MD, Heighway J, Santibanez Koref MF: An expectation-maximization algorithm for the analysis of allelic expression imbalance. Am J Hum Genet. 2006, 79: 539-543. 10.1086/506968.PubMed CentralView ArticlePubMedGoogle Scholar
- Morley M, Molony CM, Weber TM, Devlin JL, Ewens KG, Spielman RS, Cheung VG: Genetic analysis of genome-wide variation in human gene expression. Nature. 2004, 430: 743-747. 10.1038/nature02797.PubMed CentralView ArticlePubMedGoogle Scholar
- Goring HH, Curran JE, Johnson MP, Dyer TD, Charlesworth J, Cole SA, Jowett JB, Abraham LJ, Rainwater DL, Comuzzie AG, Mahaney MC, Almasy L, MacCluer JW, Kissebah AH, Collier GR, Moses EK, Blangero J: Discovery of expression QTLs using large-scale transcriptional profiling in human lymphocytes. Nat Genet. 2007, 39: 1208-1216. 10.1038/ng2119.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.