- Open Access
Expression of complement and toll-like receptor pathway genes is associated with malaria severity in Mali: a pilot case control study
Malaria Journal volume 15, Article number: 150 (2016)
The host response to infection by Plasmodium falciparum, the parasite most often responsible for severe malaria, ranges from asymptomatic parasitaemia to death. The clinical trajectory of malaria is influenced by host genetics and parasite load, but the factors determining why some infections produce uncomplicated malaria and some proceed to severe disease remain incompletely understood.
To identify molecular markers of severe falciparum malaria, human gene expression patterns were compared between children aged 6 months to 5 years with severe and uncomplicated malaria who were enrolled in a case–control study in Bandiagara, Mali. Microarrays were used to obtain expression data on severe cases and uncomplicated controls at the time of acute disease presentation (five uncomplicated and five severe), 1 week after presentation (three uncomplicated and three severe) and treatment initiation, and in the subsequent dry season (late convalescence, four uncomplicated and four severe). This is a pilot study for the first use of microarray technology in Mali.
Complement and toll-like receptor (TLR) pathways were differentially expressed, with severe cases showing higher expression of the C1q, TLR2, TLR4, TLR8, and CR1 genes. Other genes previously associated with malaria pathogenesis, GZMB, FOS and HSPA6, were also higher among severe cases. TLR2, TLR4, TLR8, CR1, GZMB, FOS, and HSPA6 genes were expressed at lower levels in severe cases at late convalescence.
Overexpression of genes previously associated with uncomplicated malaria was associated with severe disease. Low baseline expression of these genes may represent candidate markers for severe malaria. Despite the small sample size, results of this pilot study offer promising targets for follow-up analyses.
Malaria continues to be a significant cause of morbidity and mortality with an estimated 438,000 deaths and 214 million cases in 2015 . The response to infection by Plasmodium falciparum, the parasite responsible for most severe malaria, is varied and complex, ranging from asymptomatic parasitaemia to death. It has been estimated that out of 400 people who are bitten by P. falciparum infected Anopheles mosquitoes, 200 will be infected, 100 will develop the uncomplicated form of the disease, two will develop severe symptoms, and one will die . Numerous factors affect the clinical manifestation of malaria, including a patient’s age , level of P. falciparum transmission , parasite virulence , and gene expression levels in both the host and the parasite [6, 7], among others . Uncomplicated malaria is typically defined by fever accompanied by a blood smear positive for asexual P. falciparum. Severe malaria manifestations include severe anaemia, hypoglycaemia, acidosis, neurological symptoms from cerebral involvement, end-organ damage and/or respiratory distress . The biological bases for differences in clinical manifestations remain largely unexplained.
Selective pressure of P. falciparum on human evolution has helped to shape the human genome; in addition, the P. falciparum and human genomes have co-evolved because sub-clinical infection or uncomplicated disease is more advantageous to the parasite than severe disease . Examples of selection on the host genome include sickle cell trait (HbS), thalassaemias, haemoglobin C, and other erythrocyte variants . These variants have been argued to play a more prominent role than environmental factors in the outcome of malaria infection . Despite the importance attributed to genetically determined host responses, relatively few genetic variants have been identified that affect the clinical trajectory of infection. This is highlighted by the observation that sickle cell trait, the most prominent protective genetic mechanism, only explains 2 % of the total variation in response to P. falciparum infection . To address this, analyses of human genomic variation using genome-wide level data have been performed to identify additional factors , but because of linkage disequilibrium patterns in Africans, such genome-wide analyses are problematic . Alternative approaches that specifically target human genes in toto may be more fruitful. One technically feasible design is to assess genome-wide levels of expression to gauge the most important genes in response to infection.
Few studies to date have used microarray technology to identify host genes involved in response to P. falciparum infection. One study demonstrated a possible role for up-regulation of neutrophil and erythrocytic activity in 22 Kenyan children with fever and P. falciparum parasitaemia or other acute febrile illnesses . In another study, toll-like receptors were shown to be induced by malarial infection irrespective of prior exposure, demonstrating the same pattern of expression in 22 experimentally infected, naïve volunteers in the USA and 22 naturally infected Cameroonians . Although generally informative regarding innate and humoral immunity at a single time point, these studies provide data on temporally static measures against varied ethnic backgrounds, or only consider P. falciparum infection as a single clinical entity. Both of these factors may be critical in evaluating the results of expression studies.
To address the possible complicating factors of both the temporal differences in gene expression and host genetic diversity, a study in a high-transmission, ethnically homogeneous endemic region of Bandiagara, Mali was carried out. This provides a potentially powerful perspective from which to study host immune responses. By using a case–control study design, with cases of severe malaria compared to controls of uncomplicated malaria, the transcriptional profiles of these divergent host responses can be compared to detect differentially expressed gene clusters.
Children of both sexes aged between 6 months and 5 years, residing in Bandiagara, Mali and presenting to the Bandiagara District Hospital with acute clinical malaria were eligible for enrolment. Prior to enrolment, informed consent was obtained from parents or legal guardians. Other eligibility criteria included meeting definitions of either uncomplicated or severe falciparum malaria and availability for the entire six-month study period. World Health Organization (WHO) criteria for uncomplicated and severe malaria were used . The presence of either respiratory distress, anaemia (<15 % haematocrit or <5 g/dL haemoglobin), hypoglycaemia (<40 mg/dL), persistent vomiting, or the inability to eat or drink necessitating parenteral feeding also qualified the patient as having severe malaria . Exclusion criteria included the presence of any disease with clinical features similar to severe malaria (e.g., meningitis), a prior splenectomy or prolonged use of immunosuppressants, and any immunological intervention such as a vaccine or administration of immunoglobulins in the 6 months prior to enrolment in the study. Diagnoses were made during the August-December malaria transmission season.
Blood was collected from each study participant (2.5 mL) in a PAXgeneTM Blood RNA tube (PreAnalytix, Germany) that stabilizes mRNA , and 20 µL was used to prepare a Giemsa- stained blood smear to confirm P. falciparum infection. The PAXgeneTM tubes were kept at room temperature overnight in accordance with the manufacturer’s instructions (>2 h), then kept at −20 °C for 24 h and finally transferred to −80 °C for long-term storage. RNA was extracted in accordance with protocol using the PAXgeneTM Blood RNA kit. The quality of the RNA samples was evaluated with the use of Agilent RNA 6000 nano kit (Agilent Technologies, Germany) and the corresponding Agilent 2100 Bioanalyzer, using the RNA integrity number (RIN) metric . The A260/A280 ratio was also used to assess RNA purity. All microarray analyses were carried out at the Malaria Research and Training Centre in Bamako, Mali. The Affymetrix Human Genome U133 Plus 2.0 GeneChip platform and protocols (Santa Clara, USA) were used for the expression analysis. The Affymetrix GeneChip® Scanner 3000 was used to analyse the chips and Affymetrix GeneChip Operating Software version 1.2 (GCOS) was used to manage Affymetrix GeneChip data.
Data were analysed using R statistical software and the BioConductor ‘limma’ package . Preprocessing and normalization were carried out using the robust multi-array average (RMA) . Quality control measures were assessed by analysing GAPDH-6 and actin ratios as they compare to scaling factors, as well as evaluating the RNA degradation patterns for the various chips. Hierarchical clustering analyses were performed to check for any probe co-regulation patterns . In this pilot study representing the first use of microarray technology in Mali, the use of Bonferroni or false discovery rate multiple-testing adjustments was not appropriate for analysing this small dataset (Additional file 1: Figure S1). Instead the data were analysed using the criteria of greater than two-fold change [log2(2) = 1] and a significance level of p < 0.05 for testing the null hypothesis of whether the two sets are the same. To assess whether any specific biological pathways were over- or under-expressed the probes matching the above criteria were filtered using the Kyoto Encyclopedia of Genes and Genomes  pathway filter in the statistical software R. The KEGG analyses were carried out for the probes matching these criteria separately during the three time points of the study: acute illness, convalescence and late convalescence. The probes matching the criteria for the acute illness time point were also filtered using Ingenuity Pathway Analysis (Ingenuity Systems, Redwood City, CA, USA), and their expression was followed over the remaining two time points. Only the probes whose functions could be defined by Ingenuity Pathway Analysis were included in this analysis.
Real-time PCR was carried out on samples from seven individuals (four severe cases, three uncomplicated controls) for whom RNA remained after the microarray analysis (microarray analysis was carried out again on only those seven individuals for this comparison). Ten genes were analysed, eight matching the significance threshold of greater than two-fold change [log2(2) = 1] and a p value of <0.05 and two were housekeeping genes (beta-actin and GAPDH). The selected genes were: CYP1B1, C1QB, CD163, CCR2, CLEC4D, NAIP, TLR2, and IL18R1. Pre-designed primer and probe sets for all genes were obtained from Applied Biosystems . RT-PCR was carried out using the StepOnePlus™ Real-Time PCR System (Applied Biosystems), with TaqMan® reagents in accordance with the manufacturer’s recommendations. TaqMan® One-Step RT-PCR Master Mix Reagents Kit (PN 4304437) was used. Briefly, 50 ng of total RNA was used as a template for the one-step RT-PCR assay. The PCR conditions were as follows: 30 min at 48 °C, 15 s at 95 °C, 1 min at 60 °C, 40 cycles. Data were analysed using ABI software to generate delta Ct values, and the delta Ct method was used with GAPDH as Ref. . For each gene, at least one uncomplicated and one severe malaria patient had to have expression data for all of the tests done in triplicate to be included in the final analysis. For genes that had multiple probes in the Gene Chip, expression values were averaged for comparison with PCR data. The RT-PCR was carried out at Genomics/Affymetrix unit of the Malaria Research and Training Centre in Bamako, Mali.
This study was approved by the ethics committee of the Faculty of Medicine, Pharmacy and Odonto-Stomatology, University of Sciences, Techniques and Technology of Bamako, Mali.
Study population characteristics
Twenty-three uncomplicated malaria controls and 24 severe cases (with cerebral and other severe symptoms) were enrolled in this study. Five severe cases with cerebral symptoms and five uncomplicated malaria controls were chosen for pilot microarray studies of gene expression on the day of presentation at the clinic (acute illness time point). Three each of controls and cases returned seven days later to submit blood samples for a gene expression analysis a week later (convalescence time point). The clinical characteristics of the acute illness and convalescence time point patients in the study are summarized in Table 1 (Additional file 2: Table S1).
Due to loss to follow-up and development of other disease between the relevant time points, eight different patients were followed to study the late convalescence, low transmission dry season gene expression following either severe cerebral (four patients) or uncomplicated (four patients) malarial disease. These patients were closely followed to assure that no other clinical episodes occurred between their acute illness time point presentation and their late convalescence follow-up. The clinical data for this group are summarized in Table 2 (Additional file 3: Table S2).
Top differentially expressed genes during acute illness
A total of 215 probes representing 134 genes were differentially expressed between severe and uncomplicated cases and matched the inclusion criteria of >twofold change and a p value <0.05. The complete list of these probes is included in Additional file 4: Table S3. Furthermore, 37 probes corresponding to 33 genes showed differential expression at a p value <0.005 and a greater than twofold change (Table 3), and four probes, representing MB21D1, SAMSN1, PLBD1, and MKRN1 genes, had a p value <0.001. These top four genes have not been previously associated with malaria pathogenesis and/or severity. Of the remaining probes, four were in three genes, HSPA6, TLR8 and FOS, previously identified in experiments for malaria severity in human or murine systems [24–26]. Two probes interrogated the expression of the HSPA6 gene (p values of 0.0020 and 0.0040), with single association in TLR8 and FOS (p values of 0.0038 and 0.0037, respectively) (Table 3).
KEGG pathway analysis of the acute illness time point samples
A KEGG analysis based on the 215 significant probes showed that the malaria pathway had the second-highest odds ratio (Table 4). A hierarchical clustering heat map demonstrates the co-regulation of the 215 significant probes (Fig. 1a), and reveals the degree of separation between severe cases and uncomplicated ones. The odds ratios for the three top pathways, Staphylococcus aureus infection, malaria, and leishmaniasis were 15.2, 11.8 and 10.3, respectively, and all had p values of less than 0.001. Five of the 86 probed genes in the malaria pathway showed differential expression between the severe and uncomplicated case groups, with the expected value of 1. Both toll-like receptors two and four (TLR2, TLR4), part of the KEGG malaria pathway, had higher expression in the severe cases compared to the uncomplicated malaria group (2.07-fold higher for TLR2, 2.27- and 2.46-fold higher for the TLR4 probes). The other three significant genes associated with malaria in KEGG were the Complement Receptor 1 (2.05-fold higher expression in severe cases) along with glycophorin A and B (2.89- and 2.88-fold higher uncomplicated case expression, respectively). Since Staphylococcus aureus infection and leishmaniasis are infectious diseases, and each had a similar odds ratio to malaria, the top three pathways were assessed for any overlap in significant probes. The leishmaniasis pathway had three probes in common with malaria (TLR2, TLR4, CR1), while Staphylococcus aureus infection had one (CR1). Other significant KEGG pathways of relevance to malaria pathogenesis included complement and coagulation cascades and antigen processing and presentation with odds ratios of 6.4 and 4.5, respectively.
KEGG pathway analyses of the convalescence and late convalescence time point samples
Analyses of differentially expressed KEGG pathways during convalescence revealed that the ‘malaria’ pathway was no longer significant at the <0.05 level when comparing uncomplicated versus severe patients (Additional file 5: TableS4). KEGG assessment of expression differences in the two groups at the late convalescence time point showed that the malaria pathway was again statistically significant (Additional file 6: Table S5). Hierarchical clustering heat maps were created for both the convalescence (Fig. 1b) and the late convalescence (Fig. 1c) time points, and the separation of severe from uncomplicated cases observed during acute illness was maintained. For the three severe and three uncomplicated case patients followed from acute disease to convalescence, the list of KEGG pathways differentiated at the 0.05 level includes the antigen processing and presentation pathway. During acute illness, this pathway had an odds ratio of 4.5 and a p value of 0.034, and it is still significant at convalescence seven days later with an odds ratio of 8.5 and a p value of 0.002 (Additional file 5: Table S4). With only one gene expected to show significance at the given threshold, HLA-DQA1 (3.97-fold higher expression in severe cases), HLA-DQB1 (10.34-fold higher expression in severe cases), HLA-DPA1 (2.01-fold higher expression in uncomplicated cases), and KLRC3 (3.39-fold higher expression in severe cases) all matched the above stated criteria. The unrelated cases and controls from the low transmission, late convalescence time point list malaria as the 12th most significantly differentiated pathway with an odds ratio of 4.1 and a p value of 0.043 (Additional file 6: Table S5). Three genes were observed to be expressed above the predefined threshold, CD36 (4.38- and 3.81-fold higher expression in the uncomplicated cases for two probes), IL18R (2.30-fold higher expression in uncomplicated cases) and TLR4 (2.94- and 2.61-fold higher expression in uncomplicated cases), where only one was expected.
Ingenuity pathway analysis of the acute illness time point samples
The top canonical pathways produced by applying Ingenuity Pathway Analysis to significant probes at acute illness were the complement system, toll-like receptor signaling and role of pattern recognition receptors in recognition of bacteria and viruses. Figure 2 shows a stacked bar chart of these pathways and other significant pathways considered to be relevant to the pathogenesis of malaria: Cytotoxic T lymphocyte-mediated apoptosis of target cells (11th ranked pathway in terms of significance), dendritic cell maturation (13th), T cell receptor signaling (14th), T helper cell differentiation (18th), and IL-10 signaling (20th) . Of the significant probes that overlap with the above listed canonical pathways, the great majority shows lower expression in the patients who developed uncomplicated disease as compared to those who had severe symptoms. Differentially expressed probes form these canonical pathways that were previously shown to play a role in severe malaria pathogenesis were TLR2, TLR4, TLR8, HSPA6, CR1, C1qA, C1qB, C1qC, FOS, and GZMB. The TLR probes showed a higher expression in severe malaria cases during acute illness (2.07-fold for TLR2, 2.19-fold for TLR8, 2.27 and 2.46-fold for the TLR4 probes). HSPA6 levels were also higher in severe cases by 2.56-fold during this time point, as was CR1, which was expressed at a higher level in severe cases by 2.04. C1q is comprised of three chains, all of which were found to be expressed at higher levels in the severe patients (C1qA by 2.79-fold, C1qB by 3.12-fold and C1qC by 2.64-fold). The FOS and GZMB genes were also found to be expressed at elevated levels in severe malaria cases in this study (by 2.09- and 2.31-fold, respectively).
Expression of genes involved in acute illness phenotype differentiation during convalescence and late convalescence
Longitudinal expression dynamics of the probes relevant to severe versus uncomplicated malaria presentation at acute illness from the section above were evaluated during convalescence. Eleven probes corresponding to eight genes were analysed; these genes were previously found to be experimentally relevant in severe malaria. Figure 3 displays the evolution of the expression of these genes over the three time points of this study. Seven of the eight genes showed characteristics that make them reasonable targets as molecular markers of severe malaria.
TLR probes showed a higher expression in severe malaria cases during acute illness (2.07-fold for TLR2, 2.19-fold for TLR8, 2.27 and 2.46-fold for TLR4 probes). By convalescence, most of the probes still showed higher expression in severe cases, but by a smaller margin (1.03-fold for TLR2, 1.85- and 1.33-fold for the TLR4 probes), with the exception of the TLR8 probe which showed a higher expression in the uncomplicated cases (1.36-fold). This trend continued into late convalescence, where all probes showed higher expression in uncomplicated cases as compared to severe (1.94-fold for TLR2, 1.73-fold for TLR8, 2.61- and 2.94-fold for the TLR4 probes). HSPA6 levels were higher in severe cases by 2.56-fold during the acute illness time point, by convalescence uncomplicated cases were 1.17-fold higher and in late convalescence uncomplicated cases were expressed 1.90 times higher. During acute illness, CR1 was expressed at a higher level in severe cases by 2.04, 7 days later during convalescence it was 1.09 times higher is severe cases, and by late convalescence CR1 was expressed at higher levels in the uncomplicated case patients by 1.50-fold. FOS expression was higher in severe cases during acute illness as well as during convalescence, by 2.08 and 1.11, respectively, but by late convalescence, the uncomplicated cases that showed a 1.90-fold higher expression. GZMB had a 2.33-fold higher expression in severe cases during acute illness, 2.04-fold higher expression in severe cases during convalescence, but by the late convalescence time point uncomplicated cases expressed the gene 1.66-fold higher.
C1q was the only protein out of the eight candidates that did not follow a similar expression profile. Three probes corresponding to the chains of the C1q protein show a change from having higher expression in the severe patients during acute illness (C1qA by 2.79-fold, C1qB by 3.12-fold, C1qC by 2.64-fold) to showing higher expression in uncomplicated cases during convalescence (C1qA by 1.65-fold, C1qB by 1.54-fold, C1qC by 1.69-fold). By the late convalescence time point, the difference in these probes between uncomplicated and severe cases was seemingly indistinguishable between the groups (C1qA 1.10-fold higher expression in severe, C1qB no difference, C1qC 1.20-fold higher expression in severe).
The Spearman Rho for the correlation of the GeneChip and RT-PCR expression was 0.87. The Pearson’s correlation between the GeneChip expression data and RT-PCR values was 0.72, but there were not enough observations to assure a normal distribution (Additional file 7: Table S6). The correlations indicate consistency between the two approaches considering the low fold change inclusion criteria for the expression data.
Seven molecular marker candidate genes of severe malaria were identified in this study. Evaluating the baseline expression of these genes may identify individuals at risk prior to infection. When comparing gene expression between uncomplicated and severe malaria patients at acute illness, the list of differentially expressed genes was enriched for those consistently associated with malaria in prior studies. Previous studies reported differences in expression of toll-like receptors and the complement system in healthy as compared to uncomplicated P. falciparum malaria [7, 27, 28]. This study demonstrated that TLR and complement systems are not just implicated in human gene expression in response to uncomplicated symptoms of a P. falciparum infection, but also appear to play a role in the pathophysiology of severe cerebral malaria. These genes were expressed at significantly higher levels in patients presenting with severe symptoms during acute illness. All of the sub-units of C1q were elevated in severe cases, as was CR1, TLR2, TLR4, and TLR8. These results are consistent with prior studies showing dysregulation of innate immunity and inflammatory mechanisms in severe malaria morbidity [23, 25, 27, 28].
The relevance of both complement and TLR pathways is well established in response to malaria. The classical pathway of complement activation occurs via binding of antibodies in immune complexes initiated by C1q . CR1 has been shown to be protective against P. falciparum infections in low concentrations, but when over-expressed it can lead to rosetting of blood cells and subsequent severe malaria [8, 28, 29]. Single nucleotide polymorphisms (SNPs) in CR1 also associate with disease manifestation in India . The elevation of TLR probes in severe patients in this study also recapitulates previous findings . TLR4 has been associated with clinical malaria outcomes in human studies as its suppression in dendritic cells has been associated with uncomplicated P. falciparum infection in Malian children of the Dogon ethnicity . TLR8 is thought to prime and potentiate the inflammatory response to pathogens, thereby contributing to the pathogenesis of severe malaria . HSPA6 expression has also been previously implicated in malaria pathogenesis, although it has never been considered when studying the severe versus uncomplicated phenotype. Upon P. falciparum infection of nucleated erythroid precursor cells, the orthochromatic cells, the level of HSPA6 was shown to be significantly increased above what it was in controls .
This study also replicates in human P. falciparum malaria findings from mouse models for severe versus uncomplicated Plasmodium berghei malaria, which implicated genes involved in both complement and TLR, along with FOS and GZMB . In a murine model, a TLR 2/4 −/− genotype leads to increased resistance against developing experimentally induced P. berghei cerebral malaria . The c-fos protein was found to be expressed at elevated levels in experimental murine P. berghei-induced cerebral malaria . While the c-fos protein was also found to be slightly elevated in non-cerebral murine malaria cases, its level and extent of neurological involvement was far greater in the cerebral cases. This is the first study directly linking FOS expression levels to severe malaria in humans. Another mouse model found the GZMB −/− genotype to be protective of P. berghei cerebral malaria .
The use of three time points throughout the evolution of disease permitted study of the changes in expression as disease progressed from the acute illness episode, through convalescence and into the late convalescence, low transmission season. Although the patient population is not the same for the acute and late convalescence time points, the results from the late convalescence time point can stand alone, as these are children who experienced severe versus uncomplicated malaria in the previous rainy season; therefore, differences in gene expression from this time point have a prognostic value. While the various C1Q probes appeared to have similar expression patterns, the TLR probes, GZMB, FOS, HSPA6, and CR1 were all expressed at higher levels in the uncomplicated cases during the late convalescence, dry season time point. The higher expression of the TLR genes, GZMB, FOS, HSPA6, and CR1, in uncomplicated malaria cases during the late convalescence time point can be interpreted two ways. It is possible that prolonged suppression of these genes following severe infection is causal. Considering that the average time elapsed since the acute illness presentation was 150 days, and no intermittent malaria episodes occurred, this scenario is unlikely; although it cannot be ruled out. Obtaining baseline expression levels, prior to the development of malaria, would have been more informative, however there are substantial logistical obstacles on account of the relatively low prevalence of cerebral malaria. It can also be argued that the late convalescence levels represent the baseline expression values for these individuals. Under this hypothesis, the baseline expression levels of these genes are candidate markers for developing severe malarial disease upon infection.
The WHO definition of severe malaria allows for a broad spectrum of clinical symptoms, ranging from respiratory distress to cerebral complications, amongst others . While this is a clinically convenient method of clustering individuals and rationing appropriate care, it is problematic from the perspective of gene expression profiling, in that each clinical scenario of severe malaria spectrum may have a different underlying gene expression profile. All patients defined as severe cases at all time points of this study were classified based on cerebral complications; however, there is some heterogeneity even within this group (Tables 1, 2). Considering the broad spectrum of clinical features that fall into the severe malaria category, having all of the patients fall into the cerebral complications phenotype is a strength of this study.
In an experiment looking strictly into host gene expression, the differences in P. falciparum genetics and their interaction with host factors are a possible source of error . The importance of parasite genetics in determining the clinical trajectory of disease has been demonstrated in previous studies delineating the effects of var gene expression in governing clinical outcomes [33–35]. Larger studies are needed to examine parasite genetic polymorphism and gene expression differences in relation to human gene expression and clinical malaria syndromes. Host circadian rhythm expression differences could have also added to the variation between the severe and uncomplicated cases of malaria as exact blood collection times were not measured, therefore, cases and controls could not be matched on this basis. It is also possible that the findings are confounded by variable white cell counts in the patient population; however, the capacity to separate the cells was not yet available in Mali during this pilot study. Furthermore, since the objective was to find easily tested molecular markers, whole blood analyses were implemented.
The availability and increasing affordability of genome-wide expression arrays has intensified the search of human molecular markers for susceptibility to severe malaria. The low incidence of severe malaria limits the potential for genome-wide approaches from a sequencing/genome-wide association perspective because of the sample size necessary to attain significance in such studies. The presence of various sub-types of severe malaria further complicates such studies due to phenotyping considerations, again lowering the potential sample size for a well-designed study. The uncomplicated malaria phenotype has been shown to be problematic in genome-wide association studies (GWAS), where HBS did not attain genome-wide significance in the initial screen . Expression approaches such as the one described here might provide a better solution to unravelling severe malaria pathogenesis. Expression studies are more powered for finding relevant genes because they only carry out 22,000 independent tests, as opposed to the >1 million in GWAS, and gene expression is closer to the clinical phenotype than sequence variation. This is supported by the results of this study, as they were highly enriched for genes consistently associated with the phenotype in prior experiments, despite the small sample size. The results verify expression findings previously limited to the mouse model organism, or human data only previously associated with uncomplicated disease. The gene clusters highlighted by this study can guide future inquiries into unravelling the intricacies of the protective host immune response against severe malarial infection. Understanding the genetics of the host response, via gene expression, may accelerate the identification of clinically useful prognostic indicators, and ultimately, inform treatment and prevention guidelines for severe malaria.
This study suggests that overexpression of genes previously associated with uncomplicated malaria may also be associated with severe disease. Here, models of severe malaria were replicated in humans for the first time, and genes previously thought only to contribute to uncomplicated malaria were implicated in severe disease. In particular, genes in the complement and toll-like receptor pathways were expressed at higher levels in severe cases during acute illness. The longitudinal component of this study demonstrated that 7 days after presentation and treatment, the levels of these complement and toll-like receptor genes were indistinguishable in the uncomplicated and severe malaria patients. However, when expression of these genes was evaluated during the dry season following acute illness, patients with previous severe malaria showed lower levels of expression in these genes. These low baseline values should be evaluated as possible candidate biomarkers for severe malaria.
WHO. World Malaria Report. Geneva: World Health Organization; 2015.
Marsh K. Malaria–a neglected disease? Parasitology. 1992;104(Suppl):S53–69.
Djimdé AA, Doumbo OK, Traore O, Guindo AB, Kayentao K, Diourte Y, et al. Clearance of drug-resistant parasites as a model for protective immunity in Plasmodium falciparum malaria. Am J Trop Med Hyg. 2003;69:558–63.
Snow RW, Omumbo JA, Lowe B, Molyneux CS, Obiero JO, Palmer A, et al. Relation between severe malaria morbidity in children and level of Plasmodium falciparum transmission in Africa. Lancet. 1997;349:1650–4.
Gupta S, Hill AV, Kwiatkowski D, Greenwood AM, Greenwood BM, Day KP. Parasite virulence and disease patterns in Plasmodium falciparum malaria. Proc Natl Acad Sci USA. 1994;91:3715–9.
McGuire W, Hill AV, Allsopp CE, Greenwood BM, Kwiatkowski D. Variation in the TNF-alpha promoter region associated with susceptibility to cerebral malaria. Nature. 1994;371:508–10.
Idaghdour Y, Quinlan J, Goulet JP, Berghout J, Gbeha E, Bruat V, et al. Evidence for additive and interaction effects of host genotype and infection in malaria. Proc Natl Acad Sci USA. 2012;109:16786–93.
Miller LH, Baruch DI, Marsh K, Doumbo OK. The pathogenic basis of malaria. Nature. 2002;415:673–9.
World Health Organization. Communicable diseases cluster. severe falciparum malaria. Trans R Soc Trop Med Hyg. 2000;94(Suppl 1):S1–90.
Hill AV, Jepson A, Plebanski M, Gilbert SC. Genetic analysis of host-parasite coevolution in human malaria. Philos Trans R Soc Lond B Biol Sci. 1997;352:1317–25.
Kwiatkowski DP. How malaria has affected the human genome and what human genetics can teach us about malaria. Am J Hum Genet. 2005;77:171–92.
Ntoumi F, Kwiatkowski DP, Diakite M, Mutabingwa TK, Duffy PE. New interventions for malaria: mining the human and parasite genomes. Am J Trop Med Hyg. 2007;77:270–5.
Jallow M, Teo YY, Small KS, Rockett KA, Deloukas P, Clark TG, et al. Genome-wide and fine-resolution association analysis of malaria in West Africa. Nat Genet. 2009;41:657–65.
Teo YY, Small KS, Kwiatkowski DP. Methodological challenges of genome-wide association analysis in Africa. Nat Rev Genet. 2010;11:149–60.
Griffiths MJ, Shafi MJ, Popper SJ, Hemingway CA, Kortok MM, Wathen A, et al. Genomewide analysis of the host response to malaria in Kenyan children. J Infect Dis. 2005;91:1599–611.
Ockenhouse CF, Hu WC, Kester KE, Cummings JF, Stewart A, Heppner DG, et al. Common and divergent immune response signaling pathways discovered in peripheral blood mononuclear cell gene expression patterns in presymptomatic and clinically apparent malaria. Infect Immun. 2006;74:5561–73.
Rainen L, Oelmueller U, Jurgensen S, Wyrich R, Ballas C, Schram J, et al. Stabilization of mRNA expression in whole blood samples. Clin Chem. 2002;48:1883–90.
Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, Gassmann M, et al. The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 2006;7:3.
Tuimala J. DNA Microarray data analysis using Bioconductor. CSC—IT Center for Science Ltd 2008.
Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3:Article3.
Genomes KEoGa. http://www.genome.jp/kegg/. Accessed 19 Mar 2012.
Biosystems A. https://bioinfo.appliedbiosystems.com/genome-database/. Accessed 12 Nov 2012.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25:402–8.
Tamez PA, Liu H, Wickrema A, Haldar KP. falciparum modulates erythroblast cell gene expression in signaling and erythrocyte production pathways. PLoS ONE. 2011;6:e19307.
Franklin BS, Ishizaka ST, Lamphier M, Gusovsky F, Hansen H, Rose J, et al. Therapeutical targeting of nucleic acid-sensing Toll-like receptors prevents experimental cerebral malaria. Proc Natl Acad Sci USA. 2011;108:3689–94.
Ma N, Harding AJ, Pamphlett R, Chaudhri G, Hunt NH. Increased c-fos expression in the brain during experimental murine cerebral malaria: possible association with neurologic complications. J Infect Dis. 1997;175:1480–9.
Perkins DJ, Were T, Davenport GC, Kempaiah P, Hittner JB, Ong’echa JM. Severe malarial anemia: innate immunity and pathogenesis. Int J Biol Sci. 2011;7:1427–42.
Sinha S, Jha GN, Anand P, Qidwai T, Pati SS, Mohanty S, et al. CR1 levels and gene polymorphisms exhibit differential association with falciparum malaria in regions of varying disease endemicity. Hum Immunol. 2009;70:244–50.
Silver KL, Higgins SJ, McDonald CR, Kain KC. Complement driven innate immune response to malaria: fuelling severe malarial diseases. Cell Microbiol. 2010;12:1036–45.
Arama C, Giusti P, Bostrom S, Dara V, Traore B, Dolo A, et al. Interethnic differences in antigen-presenting cell activation and TLR responses in Malian children during Plasmodium falciparum malaria. PLoS ONE. 2011;6:e18319.
Kordes M, Matuschewski K, Hafalla JC. Caspase-1 activation of interleukin-1beta (IL-1beta) and IL-18 is dispensable for induction of experimental cerebral malaria. Infect Immun. 2011;79:3633–41.
Haque A, Best SE, Unosson K, Amante FH, de Labastida F, Anstey NM, et al. Granzyme B expression by CD8 + T cells is required for the development of experimental cerebral malaria. J Immunol. 2011;186:6148–56.
Kaestli M, Cockburn IA, Cortes A, Baea K, Rowe JA, Beck HP. Virulence of malaria is associated with differential expression of Plasmodium falciparum var gene subgroups in a case-control study. J Infect Dis. 2006;193:1567–74.
Rottmann M, Lavstsen T, Mugasa JP, Kaestli M, Jensen AT, Muller D, et al. Differential expression of var gene groups is associated with morbidity caused by Plasmodium falciparum infection in Tanzanian children. Infect Immun. 2006;74:3904–11.
Kyriacou HM, Stone GN, Challis RJ, Raza A, Lyke KE, Thera MA, et al. Differential var gene transcription in Plasmodium falciparum isolates from patients with cerebral malaria compared to hyperparasitaemia. Mol Biochem Parasitol. 2006;150:211–8.
RSS participated in the analytical design, statistical analyses, interpretation, writing the draft of the manuscript and revisions. AD contributed to the conception and design of the project, microarray analyses, interpretation, writing and revising the manuscript. JEM contributed to the conception and design, data collection, microarray analyses, writing and revising the manuscript. AN contributed to the conception and design, the microarray analyses, writing and revising the manuscript. JAB contributed to the statistical analyses, interpretation and revising the manuscript. AKK contributed to the design, writing and revising the manuscript. MAT contributed to the design, writing and revising the manuscript. AAD contributed to the design, writing and revising the manuscript. GV contributed to the microarray analyses, statistical analyses, and revising the manuscript. PL contributed to the microarray analyses, statistical analyses, and revising the manuscript. SMW contributed to the statistical analyses, interpretation, writing and revising the manuscript. CVP contributed to the analytical design, interpretation, writing, formatting and revising the manuscript. OKD contributed to the analytical design, interpretation, writing and revising the manuscript. All authors read and approved the final manuscript.
This work was conducted in Mali supported by a Christophe Mérieux grant to OKD. The Affymetrix platform and technology were donated, transferred and maintained at the MRTC by the Foundation Mérieux and BioMérieux. This work was also supported by the National Institutes of Health Office of the Director, Fogarty International Center, Office of AIDS Research, National Cancer Center, National Eye Institute, National Heart, Blood, and Lung Institute, National Institute of Dental and Craniofacial Research, National Institute On Drug Abuse, National Institute of Mental Health, National Institute of Allergy and Infectious Diseases, and National Institutes of Health Office of Women’s Health and Research through the Fogarty International Clinical Research Scholars and Fellows Programme at Vanderbilt University (R24 TW007988) and the American Relief and Recovery Act. This work was supported by Public Health Service award T32 GM07347 from the National Institute of General Medical Studies for the Vanderbilt Medical-Scientist Training Programme. SMW and RS were partially supported by NIH grant P20 GM103534. MAT, AK, AD, AAD, and AN were supported by grant D43TW001589 from the NIH Fogarty International Center. CVP was supported by the Doris Duke Charitable Foundation and JAB and CVP by the Howard Hughes Medical Institute.
The authors have declared that there are no competing interests.
Rafal S. Sobota, Antoine Dara and Jessica E. Manning contributed equally to this work
Additional file 1: Figure S1. Power as a function of sample size for an Affymetrix Human Genome U133 Plus 2.0 GeneChip with alpha 0.05, FDR 0.05, and a standard deviation of 0.6 under the assumption that 90% of the genes on the chip remain undifferentiated in their expression.
Additional file 2: Table S1. Clinical characteristics of enrolled patients for the acute illness and convalescence time points.
Additional file 3: Table S2. Clinical characteristics of enrolled patients for the Dry Season time point.
Additional file 4: Table S3. All statistically significant probes for testing gene expression difference between severe cases and uncomplicated controls during acute illness (function designated by Ingenuity Pathway Analysis).
Additional file 5: Table S4. Top differentially expressed KEGG pathways between severe cases and uncomplicated controls during convalescence.
Additional file 6: Table S5. Top differentially expressed KEGG pathways between severe cases and uncomplicated controls during late convalescence.
Additional file 7: Table S6. RT-PCR validation of GeneChip expression data with RT-PCR (probes for CYP1B1, IL18R1 and C1QB failed to consistently amplify, failing the inclusion criteria; therefore they were not included in the final validation analysis).
About this article
Cite this article
Sobota, R.S., Dara, A., Manning, J.E. et al. Expression of complement and toll-like receptor pathway genes is associated with malaria severity in Mali: a pilot case control study. Malar J 15, 150 (2016). https://doi.org/10.1186/s12936-016-1189-6
- Genetic Epidemiology