Skip to main content

Algorithms to predict cerebral malaria in murine models using the SHIRPA protocol

Abstract

Background

Plasmodium berghei ANKA infection in C57Bl/6 mice induces cerebral malaria (CM), which reproduces, to a large extent, the pathological features of human CM. However, experimental CM incidence is variable (50-100%) and the period of incidence may present a range as wide as 6-12 days post-infection. The poor predictability of which and when infected mice will develop CM can make it difficult to determine the causal relationship of early pathological changes and outcome. With the purpose of contributing to solving these problems, algorithms for CM prediction were built.

Methods

Seventy-eight P. berghei-infected mice were daily evaluated using the primary SHIRPA protocol. Mice were classified as CM+ or CM- according to development of neurological signs on days 6-12 post-infection. Logistic regression was used to build predictive models for CM based on the results of SHIRPA tests and parasitaemia.

Results

The overall CM incidence was 54% occurring on days 6-10. Some algorithms had a very good performance in predicting CM, with the area under the receiver operator characteristic (auROC) curve ≥ 80% and positive predictive values (PV+) ≥ 95, and correctly predicted time of death due to CM between 24 and 72 hours before development of the neurological syndrome (auROC = 77-93%; PV+ = 100% using high cut off values). Inclusion of parasitaemia data slightly improved algorithm performance.

Conclusion

These algorithms work with data from a simple, inexpensive, reproducible and fast protocol. Most importantly, they can predict CM development very early, estimate time of death, and might be a valuable tool for research using CM murine models.

Background

Cerebral malaria (CM) is one of the most important complications of Plasmodium falciparum infection. It is estimated that one in each four survivors may develop prolonged neurological and cognitive deficits [1]. CM is a multi-factorial disease with a poorly understood pathophysiology, but a number of host, parasite and epidemiological factors have been shown to play a role or influence pathogenesis [2–9]. Although no experimental model can fully reproduce all the aspects of the human disease, the murine models of CM have been instrumental in the study of CM physiopathogenesis, the most widely used model being the infection of different susceptible mouse strains by Plasmodium berghei ANKA [10], which shares many features with human CM [11]. Plasmodium berghei ANKA-infected susceptible mice develop a lethal neurological syndrome 6-12 days after infection with a cumulative incidence of 50-100% [12–14].

A number of approaches have been used to understand CM pathogenesis in the murine model, including interventional approaches, such as the blockage or addition of a given mediator [15], comparisons between CM-susceptible and resistant mouse strains [16, 17] or between CM-inducing and non-inducing parasite strains [18, 19]. Most approaches make it possible to demonstrate whether a factor is needed or not for the pathogenic process, but the precise role of each mediator and, mainly, a clear definition of the sequence of events during the course of the disease are difficult to determine. In particular, the association of physiopathological changes occurring early during infection with outcome, which is crucial for the understanding of CM pathogenesis, is curtailed by a number of factors related to the experimental model. First, the incidence of CM is variable among experiments, and can be as low as 50% [20]. Second, CM incidence can be spread over several days, most deaths occurring usually between days 6 and 9, but it can happen as early as day 5 and as late as day 12 of infection [12–14]. Therefore, the day of infection is usually a poor indicator of disease stage and, unless strict clinical assessment is performed, analyses based on the time of infection may give origin to bias derived from the different stages of the disease in different animals. Third, neurological signs (ataxia, convulsions, roll over, paralysis, coma) develop only few hours before death. Moreover, different host-parasite factors such as genetic background, age, inoculum size, course of parasitaemia and clonal variations of the parasite also interfere with the incidence of CM in mice [13, 17, 20]. These variations can, therefore, be a drawback for some experimental designs, for instance, in sequential studies aiming to identify early pathological changes related to CM development since, without a tool for predicting the outcome (development or not of the neurological syndrome), it may be complicated to determine which early changes are related to a given outcome. A way to solve or minimize this problem is to use factors or tests to build predictive models that can timely discriminate susceptible P. berghei ANKA-infected mice that will (CM+) or will not (CM-) develop CM.

It has been previously shown that a temperature below 30°C [21], the course of parasitaemia [22], and changes in locomotor activity [23] and in the TCR-β repertoire [10] are associated with death of P. berghei ANKA-infected mice. Lackner and colleagues [24], using a well-established protocol for behavioural evaluations in mice, the SHIRPA protocol, showed that CM+ mice present behavioural alterations strongly associated with a poor outcome 36 h prior to death. However, the study does not provide a tool for predicting CM development. The aim of the present study was to develop algorithms using logistic regression based on mouse behaviour data obtained with the SHIRPA primary screen and the levels of parasitaemia for early prediction of CM development in C57Bl/6 mice infected with P. berghei ANKA.

Methods

Mice

A total of 78 six- to eight-week-old female C57Bl/6 mice were obtained from CECAL/Fiocruz (Rio de Janeiro, Brazil) and used in five separate experiments. All mice were bred under specific pathogen-free (SPF) conditions and housed in groups of five (maximum), in plastic cages, with autoclaved food (Nuvital, Brazil) and filtered and autoclaved water ad libitum. All experimental protocols were reviewed and approved by the Fiocruz Ethics Committee on Animal Use (license number L-079/08).

Parasite and infection

Plasmodium berghei ANKA strain was stored in Alsever's solution in liquid nitrogen until use. A sample was thawed and 100 ÎĽ L were injected intraperitoneally (i.p.) into a donor mouse of the same age as the experimental group mice. Three days later, blood was collected and 1 Ă— 106 parasitized red blood cells (pRBC) were inoculated i.p. into each animal of the experimental group. Thin blood smears were daily made with a blood drop collected from the tip of the tail, stained according to the Panoptic method (Laborclin, Brazil) and examined under a light microscope (BH2, Olympus: Melville, New York, USA) with an oil immersion lens (1,000Ă— final magnification). Parasitaemia levels were measured by counting the percentage of pRBC in at least 2,000 RBC. Mice were observed three times a day and CM diagnosis was made based on presentation of neurological clinical signs (ataxia, disorientation, paraplegia, roll-over and coma) on days 5 - 12 post-infection.

Behavioural analyses

In order to provide a behavioural and functional profile of mice every 24 hours, from day 4 to day 10 of infection, SHIRPA primary screen was used as described elsewhere [25]. This protocol is composed of 40 separate parameters that provide a wide range of behavioural, neurological, and physiological measures. Briefly, each mouse was put under a clear cylinder with a diameter of 15 cm and 11 cm of height, which was located on a steel grid of about 1-cm mesh size. For 5 min, its behaviour was observed, the amount of rearing, grooming, and faecal pellets were counted and urination was scored as present or not. Subsequently, the mouse was transferred to a clear Perspex arena (55 Ă— 33 cm) for evaluating the transfer arousal and observation of motor behaviour. Afterwards, a sequence of manipulations using tail suspension was performed and grip strength, body tone, and reflexes were recorded. To complete the assessment, the animal was restrained in a supine position to record autonomic behaviours prior to measurement of the righting reflex. Throughout this procedure, vocalization, urination, fear, irritability and aggression were recorded. All behaviours were scored to provide a quantitative assessment of each parameter as proposed by Hatcher et al [26]. Startle response, body length and temperature were not assessed. The values of the variables tremor, palpebral closure, piloerection, gait, positional passivity, wire manoeuvre, lacrimation, salivation, provoked biting, righting reflex, negative geotaxis, and fear were reversed to get a more reasonable physiological behaviour; i.e., the best and worst scores valued maximum and minimum points, respectively [24]. A detailed description of the tests and scoring system is provided as supplemental material (See Additional file 1: SHIRPA Protocol - Primary screen). A total score and scores measuring different neurological functions - reflex and sensory function, neuropsychiatric state, motor behaviour, autonomous function, and muscle tone and strength - were generated grouping individual scores obtained in each evaluation, as previously described [24, 26]. However, different from Lackner et al [24], parameters were not normalized to a percentage of their respective maximum value before being summed up.

Statistical analyses

Results of parasitaemia levels were expressed as means and standard errors of the mean (SEM). The initial statistical analyses of the SHIRPA individual scores were performed by nonparametric tests. Based on the scale of the individual parameter, Fisher's exact test (≤ 3 parameter values) or Wilcoxon rank-sum test (if > 3 parameter values) was used to compare CM+ and CM- mice in the same day of infection. Wilcoxon matched-pairs signed-ranks test was also used to evaluate the equality of matched pairs of observations from CM- mice on days 4 - 6 of infection. A p value of less than 0.05 was considered to indicate statistical significance. Kaplan-Meier curves were used to evaluate survival pattern of mice.

Multiple logistic regression analyses were used to identify which SHIRPA individual items were the most important prognostic factors for the occurrence of CM. In multiple logistic regressions, a forward stepwise selection procedure was performed using a p-value of 0.3 for the association with CM as the criterion for entering variables in the model. A backward stepwise elimination procedure was then performed, with a p-value of 0.05 as the criterion for remaining in the model.

Predictive models based on a scoring system, with points allocated to each prognostic factor, were created from the final logistic regression model run in the test sample for days 4-6 of infection. The scoring system for each day was generated by dividing the value of the regression coefficient of each variable by the smallest coefficient and rounding the quotients to the closest integer [27]. As an example, tail elevation on day 4 had a regression coefficient of 1.115, which was then divided by the coefficient of faecal pellets (β = 0.2136), resulting in an approximate value of 5 (Item4 score, Figure 1a). The final score was obtained through the sum of each SHIRPA individual test score for those variables that remained in the final model multiplied by the coefficient generated. This approach resulted in the selection of individual SHIRPA tests that, when grouped, provided the best possible predictive power for CM development when applied on day 4 (Item4), day 5 (Item5) or day 6 (I tem6). For example, to obtain the final Item4 score for a mouse that had the following parameters on day 4: grooming in jar = 1; faecal pellets = 5; tail elevation = 3; limb tone = 1; trunk curl = 1; visual placing = 3; grip strength = 2; body tone = 1; pinna reflex = 1; corneal reflex = 1; contact righting reflex = 1; and aggression = 0, we used the formula shown in Figure 1a, as follows:

Figure 1
figure 1

Scores designed for prediction of CM using logistic regression. Item and Item+P scores were designed for prediction of CM using a stepwise multiple logistic regression analysis to identify the most important prognostic factors of SHIRPA individual scores for the occurrence of CM. Itemx and Item+Px scores are based on the day of infection and made to select animals that will develop CM. Itemy and Item+Py scores were made to predict if mice will develop CM after 24, 48 or 72 hours. RIJ = Rearing in Jar; GIJ = Grooming in Jar; U = Urination; FP = Faecal Pellets; BP = Body Position; SA = Spontaneous Activity; RRa = Respiration Rate; T = Tremor; TA = Transfer arousal; LA = Locomotor Activity; RIA = Rearing in Arena; PC = Palpebral Closure; Pi = Piloerection; G = Gait; PE = Pelvic Elevation; TE = Tail Elevation; TEs = Touch Escape; PP = Positional Passivity; TC = Trunk Curl; LG = Limb Grasping; VP = Visual Placing; GS = Grip Strength; BT = Body Tone; PR = Pinna Reflex; CR = Corneal Reflex; TP = Toe Pinch; WM = Wire Manoeuvre; SC = Skin colour; HR = Heart Rate; LT = Limb Tone; AT = Abdominal Tone; L = Lacrimation; S = Salivation; PB = Provoked Biting; RR = Righting Reflex; CRR = Contact Righting Reflex; NG = Negative Geotaxis; F = Fear; I = Irritability; A = Aggression; V = Vocalization.

Item4 score = 1 + 5 + (3 * 5) + (1 * (-4)) + (1 * (-10)) + (3 * (-10)) + (2 * (-10)) + (1 * (-10)) + (1 * (-10)) + (1 * (10)) + (1 * (-10)) + (0 * 10), which would lead to a score of -63 points for a mouse presenting those characteristics.

In the case-control study, cases were considered as mice that died in a certain time following infection (e.g., day 6) and controls were those CM- mice at that time. This selection procedure was repeated for each day of infection in which CM+ mice died. Based on this sampling scheme, separate data sets were built for predicting death with CM 24, 48, 72 and 96 hours before. For example, in the "24 h-data set", the SHIRPA individual scores and parasitaemia levels 24 hours before death were recorded for cases and 24 h before that selected day for controls. The same procedure was used for the 48-, 72-, and 96-hour data sets. A new variable indicating the day of infection was created and was always considered in the models. Then, predictive scores were built in a similar way than those generated for each day of infection.

A cross-validation procedure was performed by randomly deleting one observation at a time, then predicting each removed individual from a model fit from the remaining data. This procedure was repeated many times and sensitivity and specificity of which were calculated with predictions from the models. The data that did not include the observations were used to fit the model. Cross-validated goodness-of-fit tests were also calculated [28]. The area under the receiver operator characteristic (auROC) curve, sensitivity, specificity, and predictive values were used to evaluate models' performances. Predictive scores built in the case-control study were also evaluated in the original data set using the same cross-validation method to test models' performances in practice. All statistical analyses were performed using Stata 9.0.

Results

CM development

Inoculation of P. berghei ANKA to C57Bl/6 mice resulted in a lethal infection with a 100%-cumulative mortality on day 17. The cumulative incidence of CM was 54.3%, occurring on days 6 to 10 (Figure 2a). Mice that showed CM neurological signs (CM+ mice) died shortly thereafter, usually within a few hours. The remaining mice not affected by CM (CM- mice) died later on days 12-17 by other malaria-related pathologies (hyperparasitaemia, severe anemia) without showing any neurological sign.

Figure 2
figure 2

Survival curve (A) and course of parasitaemia (B) of mice infected with Plasmodium berghei ANKA. C57Bl/6 animals were infected with 1 × 106 Plasmodium berghei ANKA-parasitized red blood cells. CM+ = mice that developed cerebral malaria on days 6 - 8 after infection; CM- = mice that did not develop cerebral malaria; Late CM+ = mice that developed cerebral malaria on days 9 and 10 after infection. Results are expressed as means in (A) and means ± SEM in (B).

Course of parasitaemia

Parasitaemia was first detected 2 days after inoculation and grew rapidly for all infected mice up to day 4 of infection. As we have previously shown for other strains of mice [22], parasitaemia levels on subsequent days were critical for the outcome in terms of CM development. Indeed, parasitaemia levels in CM+ mice were significantly higher as compared to CM- mice on days 4 - 7 of infection (p < 0.05 in all time points, Figure 2b). Mice showing CM signs on days 9 - 10 had a different course of parasitaemia resembling CM- animals, being considered as "Late CM+" (Figure 2b).

Logistic regression models were used to assess whether parasitaemia levels on days 2-9 of infection predict the risk of CM development. Parasitaemia level of each mouse on days 4-7 was associated with the risk of CM development. Results from these logistic regression models on days 4-6 are shown in Table 1. On day 7, p-value was 0.015 and the odds ratio (OR) value was 1.068 (95% CI = 1.013-1.126). Parasitaemia level on other days was not associated with the risk of CM (p > 0.05).

Table 1 Association between parasitaemia level, total and functional category scores and the risk of CM development.

One example of the associative strength is the OR value for parasitaemia level observed on day 5 (OR = 1.080), indicating that an increase in the parasitaemia level by one unit raised the odds of CM development by 8% (Table 1). As CM+ mice started to show neurological signs and died on day 6, these data show that parasitaemia levels might be used to predict CM development at least 48 hours previous to death.

SHIRPA assessment

To evaluate whether parameters from the SHIRPA primary screen were associated with CM development on days 4-6 of infection, individual scores obtained in each evaluation from CM+ mice were compared with those obtained from CM- mice on the same day (See Additional file 2: Supplementary Table). Total and functional category scores between CM+ and CM- mice were also compared. On day 4, none of the individual tests were significantly associated with CM development (p > 0.05 for all tests) and concordantly total and functional scores were also not associated (p > 0.05 for all scores). On day 5, nine tests were associated with CM (p < 0.05 for rearing in jar, locomotor activity, rearing in arena, touch escape, limb grasping, wire manoeuvre, lacrimation, righting reflex and vocalization). The majority of these tests measure motor behaviour (locomotor activity, limb grasping, and wire manoeuvre) and neuropsychiatric state (touch escape and vocalization), explaining the association found among these functional scores and CM (p < 0.001 and p = 0.027 respectively) on day 5. On day 6, twenty five tests were significantly associated with CM development (p < 0.05 for rearing in jar, grooming in jar, body position, tremor, transfer arousal, locomotor activity, rearing in arena, palpebral closure, gait, pelvic elevation, tail elevation, touch escape, trunk curl, limb grasping, body tone, pinna reflex, toe pinch, wire manoeuvre, skin colour, limb tone, provoked biting, righting reflex, contact righting reflex, negative geotaxis and vocalization). Again the majority of tests significantly associated with CM measure motor behaviour (10 tests), but all functional scores included tests significantly associated with CM on day 6 (reflex and sensory function - 3 tests; neuropsychiatric state - 3 tests; autonomous function - 2 tests; muscle tone and strength - 2 tests). On day 6, only the autonomous function score was not associated with CM (p > 0.05). Two tests carried out on day 5 (rearing in jar and rearing in arena) and four tests on day 6 (rearing in jar, grooming in jar, rearing in arena and contact righting reflex) were significantly associated with CM in spite of not being considered for the building of functional category and total scores by Hatcher et al [26].

All comparisons were made considering CM+ and CM- mice on the same day of infection because CM- mice performance in 11 tests significantly changed during the course of infection (faecal pellets, p = 0.025; visual placing, p = 0.046; grip strength, p = 0.015; corneal reflex, p = 0.025; toe pinch, p = 0.003; wire manoeuvre, p = 0.026; righting reflex, p = 0.026; rearing in jar, p = 0.017; urination, p = 0.045; spontaneous activity, p = 0.047; piloerection, p = 0.025).

Results from logistic regression models used to assess whether total and functional category scores on days 4-6 of infection predict the risk of CM development are shown in Table 1. Total, motor behaviour and neuropsychiatric state scores of each mouse on days 5 - 6 were associated with the risk of CM development showing that these scores can be used to predict CM development at least 24 hours previous to death. Reflex and sensory function and muscle tone and strength scores were also associated with CM, but only later on (day 6) (Table 1).

Building and testing the predictive performance of scores

Since individual SHIRPA tests were not individually associated with CM development on day 4 of infection, different tests were combined in order to acquire predictive power. In addition, functional category scores indicate whether certain brain areas or functions are compromised in CM+ mice during the course of infection and are conceived to physiopathological purposes rather than to prediction. This could explain why these scores did not have a good predictive power when obtained on day 4 of infection. Therefore, new scores for days 4 - 6 were built assigning weights to each functional category score or making different combinations of specific SHIRPA tests (Itemx scores, Figure 1) aiming at improving early prediction of CM. As parasitaemia levels showed to be an early but not very powerful predictor of CM, we also built scores combining parasitemia with functional and SHIRPA test scores (Item+Px scores, Figure 1). Generated scores with the highest performance for each day are shown in Figures 1a, b, c, d, e and 1f.

As hypothesized, auROC curves increased when different weights were given to each functional score (from 38% to 66% on day 4, 72% to 78% on day 5 and 85% to 91% on day 6). auROC curves increased when parasitaemia levels were considered in the model on day 4 and 5 (from 66% to 73% on day 4 and from 78% to 82% on day 5), but this phenomenon did not happen on day 6. When only the most predictive tests for each day together with parasitaemia levels were considered, auROC curves increased even more (auROC curves = 80%, 96% and 95% on days 4 - 6, respectively). As expected, auROC curves increased with the proximity of the neurological syndrome manifestation, being more powerful on days 5 and 6 post-infection (72 - 96% and 85 - 95%, respectively) than on day 4 (38 - 80%). Although generated scores were generally more sensitive than specific, a relatively high prevalence of CM in the experiments may explain why it was observed positive predictive values (PV+) higher than negative predictive values (PV-) for most analysed scores. High PV+ values (≥ 92%) obtained for Item+Px scores on days 4-6 indicate that it is possible to early select mice with a high probability of developing CM.

Case-control study

Although it is possible to early select CM+ mice using the Item+Px scores generated as described above, this kind of approach does not allow knowing on which day post-infection mice will develop CM. Susceptible mice can develop CM between day 6 and day 10, a long period of at least 120 hours, showing that CM+ mice selected by the generated scores are in different pathophysiological states. This can be a problem because some factors necessary or involved in CM pathogenesis could be present only in the middle or late in the course of infection, not being identified when the entire CM+ selected population was analysed. A way to solve this problem is to build scores to predict on which day mice will develop CM. To accomplish this, a case-control study controlled by the day of inoculation was designed using the same experiments performed for construction of the previous scores.

Scores were generated to predict whether a mouse would develop CM in 24, 48, 72 or 96 hours later. Again scores were generated assigning different weights to each functional category score or making different combinations of specific SHIRPA tests (Itemy scores) together or not with parasitaemia (Item+Py scores). All generated scores for a 24 h-prediction showed an excellent performance (auROC ≥ 84%) with the Item+P24 score showing the best result (auROC = 95%). Scores built for a 48 h-prediction had a good performance (auROC between 78 - 94%) and the Item+P48 score had an equivalent performance when compared to the Item+P24 score (auROC = 94%). Only Item+P72 score had a good performance for a 72 h-prediction (auROC = 85%), the remaining scores having only a fair predictive power (auROC between 65 - 77%). All predictive scores for 96 h had a poor performance (auROC ≤ 72%). Again, auROC curves, sensitivity, and specificity increased with the proximity of the manifestation of the neurological signs, being more powerful within 24 hours of death. In general, specificity was greater than sensitivity for the generated scores, but a higher PV- obtained can be explained by the low prevalence of CM generated by the case-control study design. The high specificity (≥ 90%) of Item+P24 and Item+P48 scores indicates that they will possibly have higher PV+ when applied in practice and that they can be used to correctly classify CM+ mice according to the day of CM development. Generated scores with the highest performance are shown in Figures 1g, h, i, j, k and 1l.

Algorithms to predict CM development

Whether Item+Px scores would be used in situations where a minimum error of 10% is acceptable, any of them could be chosen to correctly identify CM+ mice. However, for research purposes, an error of 10% is not acceptable in many situations. It is possible to solve this problem using a cut-off that improves specificity and PV+ values. These analyses are shown in Table 2. As hypothesized, increases in the cut-off values also increases specificity and PV+ until values ≥ 95% for all Itemx and Item+Px scores.

Table 2 Best predictive performance of different logistic regression models.

To verify Itemy and Item+Py scores performance in practice and determine cut-offs that improve their specificities and PV+ values, these scores were cross-validated in the original data set used for the development of Itemx and Item+Px scores. The performance of the best predictive scores on the best days is also shown in Table 2. In general auROC curves were 5 - 20% worse on days 4 - 6 than when obtained in the case-control study. However, auROC curves remained above 79% for Item+Py and above 77% for Itemy scores on some specific days (Table 2). It was also possible to establish cut-offs with very high specificity and PV+ values for these scores on specific days (Table 2). These data show that at any time between days 4 and 6 it is possible to predict with high level of confidence which mice will develop CM (Figures 3, 4 and 5). Figures 3, 4 and 5 indeed provide researchers with a simple, fast and powerful screening/prediction tool that may prove useful in a wide range of applications. Furthermore, in some specific situations it is also possible to predict not only if but also when CM will occur: i) on day 4 it is possible to define which mice will develop CM on day 7 using Item72 and Item+P72; ii) on day 5 it is possible to define which mice will develop CM on day 6 using Item24 and Item+P24; iii) on day 6 it is possible to define which mice will develop CM on day 8 using Item48 and Item+P48 (Figures 3, 4 and 5). The cross-validation process revealed that combinations other than these do not provide high predictive powers.

Figure 3
figure 3

A step by step protocol applied on day 4 post infection for predicting CM. Item+P4 score should be used to estimate the probability of CM development on any time between days 5 - 10 post infection. Item72 and Item+P72 should be used to estimate the probability of CM development on day 7 post infection (72 hours after day 4 post infection). SHIRPA tests and its respective weights compounding each score are shown in the third line. Values obtained for each SHIRPA test should be multiplied by its respective weight and summed up to get a final score. A decision rule based on a specific final value for each score is shown. A final score above or equal (≥) this value is considered positive and is associated with a high probability of cerebral malaria (CM) development. On the other hand, a final score below this value (<) is considered negative and is associated with a low probability of CM. RIJ = Rearing in Jar; GIJ = Grooming in Jar; FP = Faecal Pellets; SA = Spontaneous Activity; RRa = Respiration Rate; TE = Tail Elevation; TC = Trunk Curl; LG = Limb Grasping; VP = Visual Placing; GS = Grip Strength; BT = Body Tone; PR = Pinna Reflex; CR = Corneal Reflex; LT = Limb Tone; RR = Righting Reflex; CRR = Contact Righting Reflex; A = Aggression; P = Parasitemia.

Figure 4
figure 4

A step by step protocol applied on day 5 post infection for predicting CM. Item5 and Item+P5 scores should be used estimate the probability of CM development on any time between days 6 - 10 post infection. Item24 and Item+P24 should be used to estimate the probability of CM development on day 6 post infection (24 hours after day 5 post infection). SHIRPA tests and its respective weights compounding each score are shown in the third line. Values obtained for each SHIRPA test should be multiplied by its respective weight and summed up to get a final score. A final score above or equal (≥) this value is considered positive and is associated with a high probability of cerebral malaria (CM) development. On the other hand, a final score under this value (<) is considered negative and is associated with a low probability of CM. RIJ = Rearing in Jar; U = Urination; BP = Body Position; RRa = Respiration Rate; T = Tremor; RIA = Rearing in Arena; G = Gait; PP = Positional Passivity; TC = Trunk Curl; LG = Limb Grasping; GS = Grip Strength; TP = Toe Pinch; WM = Wire Manoeuvre; SC = Skin colour; AT = Abdominal Tone; RR = Righting Reflex; NG = Negative Geotaxis; P = Parasitemia.

Figure 5
figure 5

A step by step protocol applied on day 6 post infection for predicting CM. Item6 and Item+P6 scores should be used estimate the probability of CM development on any time between days 7 - 10 post infection. Item48 and Item+P48 should be used to estimate the probability of CM development on day 8 post infection (48 hours after day 6 post infection). SHIRPA tests and its respective weights compounding each score are shown in the third line. Values obtained for each SHIRPA test should be multiplied by its respective weight and summed up to get a final score. A final score above or equal (≥) this value is considered positive and is associated with a high probability of cerebral malaria (CM) development. On the other hand, a final score under this value (<) is considered negative and is associated with a low probability of CM. U = Urination; RRa = Respiration Rate; TA = Transfer arousal; LA = Locomotor Activity; RIA = Rearing in Arena; PC = Palpebral Closure; Pi = Piloerection; PE = Pelvic Elevation; TE = Tail Elevation; PP = Positional Passivity; TC = Trunk Curl; LG = Limb Grasping; GS = Grip Strength; BT = Body Tone; CR = Corneal Reflex; TP = Toe Pinch; WM = Wire Manoeuvre; SC = Skin colour; HR = Heart Rate; LT = Limb Tone; AT = Abdominal Tone; PB = Provoked Biting; RR = Righting Reflex; CRR = Contact Righting Reflex; NG = Negative Geotaxis; F = Fear; A = Aggression; V = Vocalization; Cons = constant; P = Parasitemia.

Table 2 also show the score value that corresponds to the respective cut-off. This value represents the minimum value that a mouse must get in the respective score to be considered as a CM+ mouse. For example, a mouse with an Item+P4 score greater than -57 (-56, -55...) on day 4 will be considered as positive. For this cut-off point, the probability of CM development among those considered positive is 95% (Figure 3).

Discussion

Previous studies showed that parameters as temperature and parasitaemia levels are associated with CM development and death in P. berghei ANKA-infected mice [21–24], but only one study from our group [22] measured the degree of this association (using hazard ratio) and none of them showed estimates of sensitivity and specificity for CM development. Collete et al [10] showed that it is possible to discriminate CM+ from CM- B10.D2 mice infected with P. berghei ANKA by analysing the TCRβ repertoire using a throughput CDR3 spectratyping method. However, for being sophisticated, expensive, and time-consuming, the use of this method for purposes of predicting CM does not seem promising.

The SHIRPA protocol is reproducible across time and laboratories [29, 30] and it is a reliable tool for the evaluation of the neurological function in murine models of human disease [26, 31–36]. Lackner et al [24] showed that many scores obtained from the SHIRPA primary screen are associated with CM development. However, the purpose of the study was to provide researchers with a tool for improved performance of clinical assessment of CM mice, not as a tool for predicting CM development. Our analyses confirm their data that motor behaviour and neuropsychiatric state functions are early compromised in the course of infection (day 5) and that scores measuring muscle tone and strength, reflex and sensory function are only compromised late in CM+ mice. On the other hand, the performance of CM- mice also varied in the course of infection at least in eleven tests on days 4-6, showing that these mice have some degree of neurological impairment during P. berghei ANKA infection. These data agree with findings from Carvalho et al [37] showing that, although CM incidence was 70% overall in their series, 100% of P. berghei ANKA-infected CBA mice analysed presented some degree of histological alteration in the brain on days 6-8 after infection; CM- mice had microhaemorrhages, although in a much lower frequency than CM+ animals.

Animal models are important in guiding the discovery of mechanisms of disease and potential targets for prophylactic or therapeutic interventions, and they are as well a critical step in testing new developed curative or preventive interventions [38]. This is not different for CM [15]. As predictive scores that correctly select CM+ mice allow the study of these animals before the development of CM, they have multiple applications in CM studies.

A major application for a predictive score for CM murine models is on studies of pathophysiology. A sine qua non condition to attribute causality is that the factor must precede the outcome it is assumed to affect, however temporal direction might be difficult to establish if initial forms of disease, as happen in CM, were difficult to observe [39]. As Itemx and Item+Px scores can be used to correctly select CM+ mice before the development of CM on days 4 - 6, pathologic factors present in these mice that are not present in CM- mice have a high probability to be part of the causative mechanism of CM and not a consequence of it. Furthermore, Itemy and Item+Py scores can also help to determine whether a given mediator is present 24, 48 or 72 hours before death in CM+ mice and whether the blockage or addition of this mediator at the corresponding time interfere with the outcome of the disease. This approach is useful to analyse the temporal sequence of appearance and time of exposition needed so that different factors can cause CM in this model. These predictive scores can also be a suitable tool for preventive therapy and treatment studies allowing the determination of how much time before death a drug or intervention can change the outcome of the disease.

To be applied in practice, a score should be obtained from a simple and fast process. Employing proposed scores in practice is simple and direct using the scores values showed in Table 2 and Figures 3, 4 and 5. In addition, it is possible to calculate it very fast using a computer algorithm. Per mouse, it takes up to 30 minutes to carry out the entire SHIRPA protocol and to calculate parasitaemia levels by microscopy. This time may be too long for many experimental approaches, making Pond and Pond+P scores not suitable for practical use. However Item scores are composed of 8 to 21 tests, taking between 6-11 minutes per mouse to be applied. Unfortunately, parasitaemia levels alone had a poor performance in predicting CM (auROC < 72%). However, if we add this index to Item scores, time of testing increases but remains reasonable, and can still be shortened if other methods to estimate it, such as flow cytometry, are used [40, 41]. In addition, the performance of Item scores was very good and similar to that of Item+P scores to select CM+ mice (Table 2), making Item scores applicable in practice.

Conclusion

The present study provides consistent data supporting the idea that scores based on both the primary SHIRPA and from measuring parasitaemia levels are associated with and, then, can be used to predict the development of CM induced by a P. berghei ANKA infection in C57Bl/6 mice. The present study is the first to apply a commonly accepted behavioural protocol to build predictive scores of outcome in a well-established model of cerebral infectious disease.

As "Item" and "Item+P" scores are derived from a simple, highly reproducible, fast and cheap protocol, and since our results showed that they can reliably identify (or predict) CM+ mice and their time of death, we think that these scores might be a valuable tool to delineate the mechanisms underlying CM and for research with CM murine models. Neurophysiological, neuroanatomical, neurochemical and neuropharmacological research with CM murine models are some areas in which these scores can be applied.

References

  1. John CC, Bangirana P, Byarugaba J, Opoka RO, Idro R, Jurek AM, Wu B, Boivin MJ: Cerebral malaria in children is associated with long-term cognitive impairment. Pediatrics. 2008, 122: e92-99. 10.1542/peds.2007-3709.

    Article  PubMed Central  PubMed  Google Scholar 

  2. McKenzie FE, Smith DL, O'Meara WP, Riley EM: Strain theory of malaria: the first 50 years. Adv Parasitol. 2008, 66: 1-46. full_text.

    Article  PubMed Central  PubMed  Google Scholar 

  3. Idro R, Jenkins NE, Newton CR: Pathogenesis, clinical features, and neurological outcome of cerebral malaria. Lancet Neurol. 2005, 4: 827-840. 10.1016/S1474-4422(05)70247-7.

    Article  PubMed  Google Scholar 

  4. Hunt NH, Golenser J, Chan-Ling T, Parekh S, Rae C, Potter S, Medana IM, Miu J, Ball HJ: Immunopathogenesis of cerebral malaria. Int J Parasitol. 2006, 36: 569-582. 10.1016/j.ijpara.2006.02.016.

    Article  CAS  PubMed  Google Scholar 

  5. Alves FP, Durlacher RR, Menezes MJ, Krieger H, Silva LH, Camargo EP: High prevalence of asymptomatic Plasmodium vivax and Plasmodium falciparum infections in native Amazonian populations. Am J Trop Med Hyg. 2002, 66: 641-648.

    PubMed  Google Scholar 

  6. Fischer PR, Boone P: Short report: severe malaria associated with blood group. Am J Trop Med Hyg. 1998, 58: 122-123.

    CAS  PubMed  Google Scholar 

  7. Imperato PJ: Malaria parasitemia in healthy Africans in North Mara, Tanzania. J Community Health. 1986, 11: 92-97. 10.1007/BF01321510.

    Article  CAS  PubMed  Google Scholar 

  8. Lackner P, Beer R, Helbok R, Broessner G, Engelhardt K, Brenneis C, Schmutzhard E, Pfaller K: Scanning electron microscopy of the neuropathology of murine cerebral malaria. Malar J. 2006, 5: 116-10.1186/1475-2875-5-116.

    Article  PubMed Central  PubMed  Google Scholar 

  9. Dietmann A, Helbok R, Lackner P, Issifou S, Lell B, Matsiegui PB, Reindl M, Schmutzhard E, Kremsner PG: Matrix metalloproteinases and their tissue inhibitors (TIMPs) in Plasmodium falciparum malaria: serum levels of TIMP-1 are associated with disease severity. J Infect Dis. 2008, 197: 1614-1620. 10.1086/587943.

    Article  CAS  PubMed  Google Scholar 

  10. Collette A, Bagot S, Ferrandiz ME, Cazenave PA, Six A, Pied S: A profound alteration of blood TCRB repertoire allows prediction of cerebral malaria. J Immunol. 2004, 173: 4568-4575.

    Article  CAS  PubMed  Google Scholar 

  11. Hunt NH, Grau GE: Cytokines: accelerators and brakes in the pathogenesis of cerebral malaria. Trends Immunol. 2003, 24: 491-499. 10.1016/S1471-4906(03)00229-1.

    Article  CAS  PubMed  Google Scholar 

  12. Grau GE, Piguet PF, Engers HD, Louis JA, Vassalli P, Lambert PH: L3T4+ T lymphocytes play a major role in the pathogenesis of murine cerebral malaria. J Immunol. 1986, 137: 2348-2354.

    CAS  PubMed  Google Scholar 

  13. Neill AL, Hunt NH: Pathology of fatal and resolving Plasmodium berghei cerebral malaria in mice. Parasitology. 1992, 105: 165-175. 10.1017/S0031182000074072.

    Article  PubMed  Google Scholar 

  14. Rest JR: Cerebral malaria in inbred mice. I. A new model and its pathology. Trans R Soc Trop Med Hyg. 1982, 76: 410-415. 10.1016/0035-9203(82)90203-6.

    Article  CAS  PubMed  Google Scholar 

  15. de Souza JB, Riley EM: Cerebral malaria: the contribution of studies in animal models to our understanding of immunopathogenesis. Microbes Infect. 2002, 4: 291-300. 10.1016/S1286-4579(02)01541-1.

    Article  PubMed  Google Scholar 

  16. Lou J, Lucas R, Grau GE: Pathogenesis of cerebral malaria: recent experimental data and possible applications for humans. Clin Microbiol Rev. 2001, 14: 810-820. 10.1128/CMR.14.4.810-820.2001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Bagot S, Campino S, Penha-Goncalves C, Pied S, Cazenave PA, Holmberg D: Identification of two cerebral malaria resistance loci using an inbred wild-derived mouse strain. Proc Natl Acad Sci USA. 2002, 99: 9919-9923. 10.1073/pnas.152215199.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Beghdadi W, Porcherie A, Schneider BS, Dubayle D, Peronet R, Huerre M, Watanabe T, Ohtsu H, Louis J, Mecheri S: Inhibition of histamine-mediated signaling confers significant protection against severe malaria in mouse models of disease. J Exp Med. 2008, 205: 395-408. 10.1084/jem.20071548.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Steen Van den PE, Van Aelst I, Starckx S, Maskos K, Opdenakker G, Pagenstecher A: Matrix metalloproteinases, tissue inhibitors of MMPs and TACE in experimental cerebral malaria. Lab Invest. 2006, 86: 873-888. 10.1038/labinvest.3700454.

    Article  PubMed  Google Scholar 

  20. Amani V, Boubou MI, Pied S, Marussig M, Walliker D, Mazier D, Renia L: Cloned lines of Plasmodium berghei ANKA differ in their abilities to induce experimental cerebral malaria. Infect Immun. 1998, 66: 4093-4099.

    PubMed Central  CAS  PubMed  Google Scholar 

  21. Curfs JH, Schetters TP, Hermsen CC, Jerusalem CR, van Zon AA, Eling WM: Immunological aspects of cerebral lesions in murine malaria. Clin Exp Immunol. 1989, 75: 136-140.

    PubMed Central  CAS  PubMed  Google Scholar 

  22. Martins YC, Smith MJ, Pelajo-Machado M, Werneck GL, Lenzi HL, Daniel-Ribeiro CT, Carvalho LJ: Characterization of cerebral malaria in the outbred Swiss Webster mouse infected by Plasmodium berghei ANKA. Int J Exp Pathol. 2009, 90: 119-130. 10.1111/j.1365-2613.2008.00622.x.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Ball HJ, MacDougall HG, McGregor IS, Hunt NH: Cyclooxygenase-2 in the pathogenesis of murine cerebral malaria. J Infect Dis. 2004, 189: 751-758. 10.1086/381503.

    Article  CAS  PubMed  Google Scholar 

  24. Lackner P, Beer R, Heussler V, Goebel G, Rudzki D, Helbok R, Tannich E, Schmutzhard E: Behavioural and histopathological alterations in mice with cerebral malaria. Neuropathol Appl Neurobiol. 2006, 32: 177-188. 10.1111/j.1365-2990.2006.00706.x.

    Article  CAS  PubMed  Google Scholar 

  25. Rogers DC, Fisher EM, Brown SD, Peters J, Hunter AJ, Martin JE: Behavioral and functional analysis of mouse phenotype: SHIRPA, a proposed protocol for comprehensive phenotype assessment. Mamm Genome. 1997, 8: 711-713. 10.1007/s003359900551.

    Article  CAS  PubMed  Google Scholar 

  26. Hatcher JP, Jones DN, Rogers DC, Hatcher PD, Reavill C, Hagan JJ, Hunter AJ: Development of SHIRPA to characterise the phenotype of gene-targeted mice. Behav Brain Res. 2001, 125: 43-47. 10.1016/S0166-4328(01)00275-3.

    Article  CAS  PubMed  Google Scholar 

  27. Barquet N, Domingo P, Cayla JA, Gonzalez J, Rodrigo C, Fernandez-Viladrich P, Moraga-Llop FA, Marco F, Vazquez J, Saez-Nieto JA, Casal J, Canela J, Foz M: Prognostic factors in meningococcal disease. Development of a bedside predictive model and scoring system. Barcelona Meningococcal Disease Surveillance Group. JAMA. 1997, 278: 491-496. 10.1001/jama.278.6.491.

    Article  CAS  PubMed  Google Scholar 

  28. Crossval. A macro for Stata. [http://www.biostat.jhsph.edu/~ejohnson/regression/crossval.ado]

  29. Mandillo S, Tucci V, Holter SM, Meziane H, Banchaabouchi MA, Kallnik M, Lad HV, Nolan PM, Ouagazzal AM, Coghill EL, Gale K, Golini E, Jacquot S, Krezel W, Parker A, Riet F, Schneider I, Marazziti D, Auwerx J, Brown SD, Chambon P, Rosenthal N, Tocchini-Valentini G, Wurst W: Reliability, robustness, and reproducibility in mouse behavioral phenotyping: a cross-laboratory study. Physiol Genomics. 2008, 34: 243-255. 10.1152/physiolgenomics.90207.2008.

    Article  PubMed Central  PubMed  Google Scholar 

  30. Rogers DC, Peters J, Martin JE, Ball S, Nicholson SJ, Witherden AS, Hafezparast M, Latcham J, Robinson TL, Quilter CA, Fisher EM: SHIRPA, a protocol for behavioral assessment: validation for longitudinal study of neurological dysfunction in mice. Neurosci Lett. 2001, 306: 89-92. 10.1016/S0304-3940(01)01885-7.

    Article  CAS  PubMed  Google Scholar 

  31. Hatcher JP, Virley D, Hadingham SJ, Roberts J, Hunter AJ, Parsons AA: The behavioural effect of middle cerebral artery occlusion on apolipoprotein-E deficient mice. Behav Brain Res. 2002, 131: 139-149. 10.1016/S0166-4328(01)00377-1.

    Article  CAS  PubMed  Google Scholar 

  32. Hunter AJ, Hatcher J, Virley D, Nelson P, Irving E, Hadingham SJ, Parsons AA: Functional assessments in mice and rats after focal stroke. Neuropharmacology. 2000, 39: 806-816. 10.1016/S0028-3908(99)00262-2.

    Article  CAS  PubMed  Google Scholar 

  33. Harms LR, Eyles DW, McGrath JJ, Mackay-Sim A, Burne TH: Developmental vitamin D deficiency alters adult behaviour in 129/SvJ and C57BL/6J mice. Behav Brain Res. 2008, 187: 343-350. 10.1016/j.bbr.2007.09.032.

    Article  CAS  PubMed  Google Scholar 

  34. Dadachova E, Moadel T, Schweitzer AD, Bryan RA, Zhang T, Mints L, Revskaya E, Huang X, Ortiz G, Nosanchuk JS, Nosanchuk JD, Casadevall A: Radiolabeled melanin-binding peptides are safe and effective in treatment of human pigmented melanoma in a mouse model of disease. Cancer Biother Radiopharm. 2006, 21: 117-129. 10.1089/cbr.2006.21.117.

    Article  CAS  PubMed  Google Scholar 

  35. Norreel JC, Jamon M, Riviere G, Passage E, Fontes M, Clarac F: Behavioural profiling of a murine Charcot-Marie-Tooth disease type 1A model. Eur J Neurosci. 2001, 13: 1625-1634. 10.1046/j.0953-816x.2001.01535.x.

    Article  CAS  PubMed  Google Scholar 

  36. Rafael JA, Nitta Y, Peters J, Davies KE: Testing of SHIRPA, a mouse phenotypic assessment protocol, on Dmd(mdx) and Dmd(mdx3cv) dystrophin-deficient mice. Mamm Genome. 2000, 11: 725-728. 10.1007/s003350010149.

    Article  CAS  PubMed  Google Scholar 

  37. Carvalho LJ, Lenzi HL, Pelajo-Machado M, Oliveira DN, Daniel-Ribeiro CT, Ferreira-da-Cruz MF: Plasmodium berghei: cerebral malaria in CBA mice is not clearly related to plasma TNF levels or intensity of histopathological changes. Exp Parasitol. 2000, 95: 1-7. 10.1006/expr.2000.4508.

    Article  CAS  PubMed  Google Scholar 

  38. Chia R, Achilli F, Festing MF, Fisher EM: The origins and uses of mouse outbred stocks. Nat Genet. 2005, 37: 1181-1186. 10.1038/ng1665.

    Article  CAS  PubMed  Google Scholar 

  39. Hofler M: The Bradford Hill considerations on causality: a counterfactual perspective. Emerg Themes Epidemiol. 2005, 2: 11-10.1186/1742-7622-2-11.

    Article  PubMed Central  PubMed  Google Scholar 

  40. Hein-Kristensen L, Wiese L, Kurtzhals JA, Staalsoe T: In-depth validation of acridine orange staining for flow cytometric parasite and reticulocyte enumeration in an experimental model using Plasmodium berghei. Exp Parasitol. 2009, 123: 152-157. 10.1016/j.exppara.2009.06.010.

    Article  CAS  PubMed  Google Scholar 

  41. Cabrales P, Zanini GM, Meays D, Frangos JA, Carvalho LJ: Murine Cerebral Malaria Is Associated with a Vasospasm-Like Microcirculatory Dysfunction, and Survival upon Rescue Treatment Is Markedly Increased by Nimodipine. Am J Pathol. 2010, 176: 1306-1315. 10.2353/ajpath.2010.090691.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

To Marilene Damazio dos Santos and Luiz Carlos da Silva from the PavilhĂŁo LeĂ´nidas Deane Animal House for caring mice. To Belmira Santos from the Fiocruz Animal House for providing the C57Bl/6 mice. To Gabriel Farah for helping with SHIRPA experiments. Frederico Meister de Freitas helped to format the manuscript. This work was supported with funds of the Laboratory of Malaria Research (POM/PEF, IOC-Fiocruz) and of the Conselho Nacional de Desenvolvimento TecnolĂłgico (CNPq) - Brazil. YCM and CTDR are recipients of a Doctoral (CAPES) and a Research Produtivity (CNPq) fellowships, respectively, and GLW was partially supported by CNPq (Process 308889/2007-0 and 410590/2006-1). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yuri C Martins.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

YCM carried out the SHIRPA experiments, participated in the statistical analysis and drafted the manuscript. GLW participated in the design of the study, performed the statistical analysis and helped to draft the manuscript. LJMC participated in the design of the study and helped to draft the manuscript. BPTS and BGA participated in the SHIRPA experiments. TMS trained YCM in the SHIRPA tests and helped to draft the manuscript. DOG participated in the design of the study and hosted some trainees in the SHIRPA tests in his laboratory. CTDR conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: A protocol showing how the primary screen of SHIRPA tests were performed and scored. (PDF 60 KB)

12936_2010_1166_MOESM2_ESM.DOC

Additional file 2: Values of individual SHIRPA tests and functional category scores of infected non-CM (CM-) mice on days 4 - 6 post-infection were compared with cerebral malaria (CM+) mice in this Table. (DOC 88 KB)

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Martins, Y.C., Werneck, G.L., Carvalho, L.J. et al. Algorithms to predict cerebral malaria in murine models using the SHIRPA protocol. Malar J 9, 85 (2010). https://doi.org/10.1186/1475-2875-9-85

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2875-9-85

Keywords