Open Peer Review
Risk factors for death among children 0–59 months of age with moderate-to-severe diarrhea in Manhiça district, southern Mozambique
BMC Infectious Diseasesvolume 19, Article number: 322 (2019)
Despite major improvements in child survival rates, the number of deaths due to diarrhea remains unacceptably high. We aimed to describe diarrhea-associated mortality and evaluate risk factors for death among Mozambican children with moderate-to-severe diarrhea (MSD).
Between December 2007 and November 2012, children under-five with MSD were enrolled in Manhiça district, as part of the Global Enteric Multicenter study (GEMS). Clinical, epidemiological, and socio-demographic characteristics were collected. Anthropometric measurements were performed and stool samples collected upon recruitment. A follow-up visit ~ 60 days post-enrolment was conducted and verbal autopsies performed in all death cases.
Of the 916 MSD-cases analyzed; 90% (821/916) completed 60 days follow-up and 69 patients died. The case fatality rate at follow-up was 8% (69/821), and the mortality rate 10.2 (95%CI: 7.75–13.59) deaths per 1000 persons-week at risk. Nearly half of the deaths 48% (33/69) among study participants clustered within 2 weeks of the onset of diarrhea.
Typical enteropathogenic Escherichia coli (typical EPEC) and Cryptosporidium were the two pathogens associated to an increased risk of death in the univariate analysis with (HR = 4.16, p = 0.0461) and (H = 2.84, p = 0.0001) respectively. Conversely, Rotavirus infection was associated to a decreased risk of death (HR = 0.52, p = 0.0198).
According to the multivariate analysis, risk factors for death included co-morbidities such as malnutrition (HR = 4.13, p < 0.0001), pneumonia/lower respiratory infection (HR = 3.51, p < 0.0001) or invasive bacterial disease (IBD) (HR = 6.80, p = 0.0009), presenting on arrival with lethargy or overt unconsciousness (HR = 1.73, p = 0.0302) or wrinkled skin (HR = 1.71, p = 0.0393), and cryptosporidium infection (HR = 2.14, p = 0.0038).
When restricting the analysis to those with available HIV results (n = 191, 22% of the total study sample), HIV was shown to be a significant risk factor for death (HR = 5.05, p = 0.0009).
Verbal autopsies were conducted in 100% of study deaths, and highlighted diarrhea as the main underlying cause of death 39%, (27/69); followed by HIV/AIDS related deaths 29.0% (20/69) and sepsis 11.6% (8/69).
Preventive strategies targeting Cryptosporidium, malnutrition and early identification and treatment of associated co-morbidities could contribute to the prevention of the majority of diarrhea associated deaths in Mozambican children.
In spite of very significant decreasing trends in the last three decades, childhood mortality remains unacceptably high globally, with 5.6 million children dying every year before reaching their 5th birthday, the majority of which in low and middle-income countries (LMIC) . Diarrheal diseases still represent a major cause of morbidity and mortality in childhood, and are believed to account for 499,000–526,000 annual child deaths, nearly 9% of all under five global mortality [2, 3], in spite of the good intake of life-saving interventions such as oral rehydration solution (ORS)  and the rotavirus vaccine . Deaths in sub-Saharan Africa and Southeast Asia account for ~ 78% of deaths due to diarrhea worldwide [5, 6], underscoring the inequities related to this particular syndrome.
The Global Enteric Multi-center Study (GEMS) reported that most episodes of moderate-to-severe diarrhea (MSD) among children under 5 years old were primarily due to four pathogens: Rotavirus, Cryptosporidium, Enterotoxigenic Escherichia coli (ST-ETEC), and Shigella . The risk of dying from diarrheal disease was reported to be higher among children younger than 2 years of age, albeit with relatively different rates from one region to another . The advent of the HIV/AIDS pandemic has also changed the incidence , clinical presentation and outcome of diarrheal diseases, as the immunosuppression derived from the infection favors a higher incidence of gastrointestinal infections, not only from “classical” diarrhea pathogens, but also from more aggressive opportunistic infections, typical of the immunocompromised host. This has resulted in significant changes in the last decades in the spectrum, clinical presentation, duration and prognosis of diarrheal episodes in those countries where HIV is highly prevalent [9, 10].
Malnutrition is another public health problem in children in LMIC, often associated with an increased risk of mortality and commonly a cause and a consequence of diarrheal disease. Indeed, diarrhea episodes among malnourished children seem to be more complicated and prolonged than among non-malnourished children, and convey a worse prognosis [11,12,13].
In Mozambique, a paradigmatic example of a resource-constrained sub-Saharan African country, it has been estimated that diarrheal diseases are the third leading cause of death among the pediatric population (0–14 years), accounting for one out of every 10 deaths in this age group . Data from the Manhiça district, in the Southern part of the country, confirm diarrheal diseases as the third leading cause of hospital admission among children 0–14 years and the fourth leading cause of death among children 12–59 months of age [15,16,17,18]. Significantly, in this same district, both HIV infection and malnutrition appear to be major public health problems, with HIV prevalence rates among adults at the community level peaking at 40% , and malnutrition (of any degree) affecting 47% of all children seen as outpatients .
Despite the significant diarrheal mortality burden, few studies have assessed the specific risk factors associated to a fatal outcome in LMIC during the diarrheal episodes. We propose in this analysis, to determine the possible risk factors associated to death among children aged 0–59 months presenting with MSD in Manhiça, a rural district in southern Mozambique with a high prevalence of infectious diseases.
Location of the study
The Manhiça Health Research Centre (CISM) coordinated this study, conducted in six health units in the District of Manhiça, southern Mozambique, including five peripheral health posts and the Manhiça District hospital, which acts as the referral hospital for the entire Manhiça District. This district, which is part of the wider Maputo province, is about 80 km north of Maputo, the countries’ capital. CISM established there in the year 1996 a research institution, which since then has been running in the district different research platforms, including a demographic surveillance system (DSS), covering when the study was conducted (2007–2012) an estimated population of around 92,000 inhabitants, equaling about half of the district’s total population . Previous reports have detailed the idiosyncrasies of the area .
In parallel to the DSS, CISM also established a morbidity surveillance platform at the Manhiça District Hospital (MDH) and other surrounding smaller health posts. Such platform was designed to passively capture all pediatric (< 15 years of age) outpatient visits and all ward hospitalizations . A permanent identification number (“perm-ID”), unique to each individual in the study area, and issued by the DSS, allows the round-the-clock linkage between the demographic and morbidity platforms. [23, 24].
The MDH is a busy district hospital which sees around 3500 pediatric admissions annually, and which receives over 50,000 outpatient visits. It is a 110-bed hospital, which dedicates 36 of its beds to child hospitalizations. It also includes a malnutrition specialized hospitalization unit. The MDH is supported in terms of staffing and resources by CISM. The morbidity platform uses standardized questionnaires for routine data collection, which collect demographic, clinical, laboratory and outcome data. Anthropometric measures collected include weight for all patients, and height/length among hospitalized patients. A basic microbiology surveillance protocol has also been ongoing since the year 2003, whereby for all admitted children under the age of 24 months, a blood culture is collected. For older children, this is also done in case of hyperpyrexia (documented temperature > 39°c), or based on the presence of severity criteria, or to the admitting clinician’s judgement [25, 26].
Study population and clinical procedures
The Global Enteric Multicenter Study (GEMS) of diarrheal disease in infants and young children was conducted in seven diverse, high-burden sites in Asia and Africa, including Manhiça, in southern Mozambique. In Manhiça, GEMS recruited patients during a five-year period (Dec. 2007 to Nov. 2012). GEMS was a case-control study, stratified by age. It aimed to unravel the epidemiology and causes of pediatric diarrheal disease globally and at each of the participating sites. Three age groups were targeted by the study: infants (0- < 12 months), toddlers (12- < 24 months), and older children (24–- < 60 months). During the 5 years during which the study was ongoing, study staff screened for diarrhea any child under the age of five registered in the center’s DSS, who actively attended any of the participating health units. The definition for a “diarrhea episode” included the occurrence of a minimum of three loose stools in the preceding twenty-four hours.
Diarrheal episodes needed to fulfil a set of inclusion criteria to be eligible for enrolment as cases, namely that the episode was new (starting at least after a minimum of seven days in which the child had not had diarrhea), acute (shorter than a week since its beginning), or showed any of the following signs of severity, so as to fulfil the definition of “moderate-to-severe diarrhea” (MSD):
Variables associated with dehydration: sunken eyes, decrease in skin turgor (defined by a long [≤2 s] or very long [> 2 s] abdominal skin pinch recoil), perceived need or confirmed administration of parenteral re-hydration,
dysentery (stools with discernable blood); or
hospitalization because of either diarrhea or dysentery.
As this was a case-control study, 1 to 3 healthy children (with no story of diarrhea in the preceding 7 days) were randomly selected using the DSS databases from the Manhiça district matched to the index case by neighborhood, age and gender. Controls were recruited no later than 2 weeks after the recruitment of each index case.
The study recorded for all enrolled individuals a variety of clinical (including anthropometric), epidemiological, and socio-demographic information. Fecal specimens were obtained for pathogen screening and characterization, following previously published methodologies . As Manhiça is a highly HIV-endemic, all children were actively screened for HIV infection, after appropriate HIV counseling. This was however only started in May 2010, when the study had already recruited about a half of the required patients and continued for all study participants subsequently. Screening procedures, which followed Mozambican recommendations, utilized a two-step approach, which included first a serological screening with the Determine HIV-1/2 Rapid Test (Abbott Laboratories, Abbott Park, IL), with confirmation for all initially positive cases with a second serological test (Uni-Gold Rapid Test; Trinity Biotech Co, Wicklow, Ireland). Molecular confirmation with DNA PCR was required to confirm HIV positivity in those children under the age of 18 months with positive HIV serology. This confirmation was not required for children above the age of 18 months which were positive with serology.
The GEMS protocol included a follow-up visit at the families’ household for each enrolled participant, approximately two months after recruitment (acceptable range: 49–90 days), with the objective of verifying the health status and outcome of the diarrheal episode. Up to a maximum of 3 different contacts were attempted to trace the children who could not be found during their 60-day visit, after which, if still not found, they were considered losses to follow-up [7, 27].
In the first three months after study recruitment, any death occurring among study participants was investigated using the verbal autopsy methodology. The INDEPTH-modified standard WHO VA questionnaire was utilized. All VA questionnaires were individually reviewed by two clinicians, and discrepant diagnoses were adjudicated by an expert cause of death reviewer (QB). Of note, each clinician could attribute up to a maximum of two plausible causes of death. The final CoD diagnosis was reached when there was agreement in either one or both of the diagnoses between the two reviewer’s and/or the adjudicator. Cases in which at least two of the reviewing clinicians considered the CoD diagnosis to be “unknown” were considered undetermined. When two different final diagnoses were assigned for a given death, each of these was individually mapped onto ICD-10 in order to calculate cause-specific rates.
Data management and statistical analysis
All data were collected using specific questionnaires, which were reviewed locally prior to being scanned and transferred to a centralized Data Collection Centre (DCC), in the US. Queries were issued centrally, but resolved locally, based on the revision of the original source documents stored on site. Stata/SE software (version 14.0,  and the package coxphf from R version 3.1.1 were utilized for data analyses [29, 30]. In order to estimate measures of disease burden such as diarrhea-associated mortality Rates in children, population weights based on sampling weights were calculated.
Mortality analyses were performed for MSD cases. Any death occurring during the period from recruitment to the day 60 follow-up visit was considered a study death. Time to death (or censoring) was calculated as the time span between the beginning of the diarrhea episode (first day of diarrhea, defined as the date of recruitment minus the duration of diarrhea at that moment) and the date of death (or the 60-day visit). Kaplan-Meier survival curves were estimated, and the difference between age groups was assessed by weighted-Cox regression estimation. Mortality Rates and appropriate centiles of survival time were also described (centiles in accordance to estimated survival curves).
Cox regression with Firth’s Penalized Likelihood adjusted for age [31, 32] was used to estimate the associations between the different covariates (including clinical characteristics, nutrition practices, breastfeeding and pathogens detected) and time-to-death. Multivariable models were estimated by a forward-stepwise procedure, with p-to-enter and p-to-leave both equal to 0.05 according to the Wald test. A specific multivariate analysis was conducted as a sensitivity analysis only among children with HIV results available.
In order to get a better understanding of the total community burden of diarrhea, Health Care Utilization and Attitudes Surveys (HUAS) were repeatedly conducted during the study period to investigate the what was in the area the proportion of cases usually seeking care at the different health units in the seven days following the onset of the MSD case (called r). Such surveys were conducted on random samples of children, twice a year, coinciding with each round of the DSS. Data from such surveys were pooled and weighted using the sampling weights (inverse of sampling fraction), based on the number of children in each age-sex stratum according to information derived from CISM’s DSS at the time of the round . After that, the values of “r” and its variances were estimated for each age stratum, using Kaplan Meier analysis.
A weighted analysis, as previously described by Sommerfelt et al , was utilized for the extrapolation of mortality rates to the population levels. Such method proposes to choose weights so as to make the proportion of cases in the analysis similar to that of the population from which the cases are drawn. For the calculation of weights applicable to controls, the inverse of sampling fraction within each age stratum was used, according to the formula: (N-cases in_pop)/controls, where N is the median population; cases_in pop is the number of cases in the population estimated by the number of eligible cases divided by r and controls implies the number of enrolled controls for that stratum. Weights for cases were constructed according to the sampling weights adjusted for fortnight quote (described above) divided by r.
This study, which was conducted under the GEMS framework, was coordinated by the Center for Vaccine Development (Baltimore, Maryland, USA) and led by CISM in Manhiça, in close collaboration and partnership with the Barcelona institute of Global Health (ISGlobal), in Spain. All GEMS-related protocols and informed consent documents were reviewed and approved prior to the initiation of the study by the National Bioethics Committee for Health of Mozambique (CNBS – IRB00002657), the ethics committee of the Hospital Clínic of Barcelona, and The Institutional Review Board at the University of Maryland. All parents/guardians of participating children were requested to sign a written informed consent after adequate information regarding the study was provided.
The study included data on childhood deaths among children recruited in Manhiça with MSD, between December 2007 and November 2012 as part of the GEMS. We enrolled a total of 916 children with MSD and 90% (821/916) of them completed the day 60 follow-up home visit (Fig. 1). Reasons for not completing the 60-day follow-up were essentially due to either changes of residence to other areas outside the DSS area, or temporary absences from their homes coinciding with the 60-day follow-up visits. In these cases, and according to the study protocol, children were declared lost to follow-up if 3 consecutive visits were conducted after day 60 with no success in finding the child.
As HIV counseling and testing was only introduced for study participants in May 2010, HIV data are missing for many recruited participants (100% for year 1 and 2; 73% for year 3, 36% for year 4 and 16% for year 5). However, no relevant differences were shown when comparing children with or without HIV available data (data not shown).
The proportion of children with MSD who were taken to the hospital within the first 7 days from onset of episode was 66, 59 and 41% among infants, toddlers and children, respectively (Fig. 2).
Sixty-nine deaths were detected among study participants during the study follow-up (recruitment-day up to last follow-up visit), providing a case fatality rate of 8% (69/821). The mortality rate after adequately weighting was 10.2 (95%CI: 7.75–13.59) deaths per 1000 persons-week at risk (PWAR). According to the age group the mortality rate was 10.9 (95%CI: 7.71–15.90) deaths per 1000 PWAR among infants, 11.9 (95%CI: 7.82–18.81) deaths per 1000 PWAR among toddlers and 5.9 (95%CI: 1.85–28.55) deaths per 1000 PWAR among children 24–59 months of age (p = 0.5357) (Fig. 3 and Table 1). No differences were found in terms of the distribution of deaths in each age group in relation to time, with 48% (33/69) of all deaths occurring within the first two weeks after the onset of the diarrhea episode, of which 25% (17/69) in the first 7 days, and 23% (16/69) after the first week and before completing 14 days of onset. In terms of distribution of deaths according to sex, no significant differences were observed [deaths among boys 57% (39/69) compared to 43% (30/69) among girls; hazard ratio of death [HR = 1.12 (95%CI: 0.79–1.80, p = 0.6199)].
Pathogens associated with an increased risk of death in the univariate analysis were Cryptosporidium (HR = 2.84 (95%CI: 1.72–4.58, p = 0.0001) and typical EPEC (HR = 4.16 (95%CI: 1.02–38.05, p = 0.0461), with Astrovirus infection (HR = 3.51, 95%CI: 0.97–8.93, p = 0.0545) bordering statistical significance. Conversely, cases infected with Rotavirus were less likely to die due to diarrhea (HR = 0.52, 95%CI: 0.38–0.90, p = 0.0198) (Table 2).
Univariate associations between clinical characteristics and time-to-death were also assessed (Table 3). Admission to hospital implied a higher risk of death (HR 3.53, 95%CI 1.95–7.52, p < 0.0001). The median time of hospitalization was significantly longer for those who died (6.5 days) than for those who survived (4.0 days, p < 0.0001). Clinical variables on arrival associated to mortality included, among others, lethargy or impaired consciousness, fever, signs of dehydration such as loss of skin turgor and skin pinch, chest indrawing and breathing difficulties, signs associated to malnutrition (“flaky paint” skin appearance, abnormal hair) and rectal prolapse. Only a small fraction of recruited children 22% (191/821) had available HIV results, and positivity rate among these was high 21% (40/191). Among these, HIV infection was strongly correlated with the risk of death (HR = 6.04, 95%CI 2.38–16.09, p = 0.0002). Invasive bacterial disease (bacteremia in an admitted child) (HR = 17.44, 95%CI 6.43–38.15, p < 0.0001) and pneumonia/lower respiratory infections (HR = 3.51, 95%CI 2.16–5.65, p < 0.0001) diagnosed during hospital admission were co-morbidities strongly associated with a lethal outcome. Malnourished children had also a significantly increased risk of dying (HR = 5.66, 95%CI 3.43–9.14, p < 0.0001). After adjusting for anthropometric characteristics at enrollment, breastfed children had a lower mortality risk (HR = 0.28, 95%CI 0.13–0.69 for partially breastfeed children; HR = 0.27, 95%CI 0.11–0.64 for exclusively breastfed children; p = 0.0036)]. Conversely, an admission diagnosis of dysentery was associated with a 70% reduction in the risk of dying (HR = 0.30, 95%CI 0.08–0.78, p = 0.0112).
Independent risk factors for death in the multivariate analysis included co-morbidities such as malnutrition (HR = 4.13, 95%CI: 2.42–6.96; p < 0.0001), pneumonia/lower respiratory infection (HR = 3.50, 95%CI 2.15–5.67, p < 0.0001) or invasive bacterial disease (IBD) (HR = 6.80, 95%CI 2.43–15.50, p = 0.0009), presenting on arrival with lethargy or overt unconsciousness (HR = 1.73, 95%CI 1.05–2.90, p = 0.0302) or wrinkled skin (HR = 1.71, 95%CI 1.02–2.84, p = 0.0393), and a diagnosis of cryptosporidium infection (HR = 2.14, 95%CI 1.28–3.58, p = 0.0038) (Table 4).
Verbal autopsy cause of death diagnoses
Verbal autopsies were conducted to all 69 deaths documented in study participants (Additional file 1: Table S1). Final cause of death (CoD) ascertained by this method highlighted diarrhea as the main underlying CoD 39% (27/69), followed by HIV/AIDS related deaths 29% (20/69) and sepsis 12% (8/69).
The unprecedented decreases in child mortality that the world has experienced in the last years is the result of multiple factors, including significant reductions in certain major causes of under 5 mortality, such as childhood diarrhea. Indeed, the last three decades have seen diarrhea deaths dropping down from the 1.748 million deaths in the year 1990 , to the estimated half a million deaths in 2015 . However, understanding that good preventive and curative strategies are now widely implemented, diarrhea associated mortality remains unacceptably high. This study has highlighted the associated impact of diarrheal disease in Manhiça, a rural area of Southern Mozambique, whereby up to 8.4% of the children recruited during a 5-year long period died within the first 2 months after the episode, yielding a population weighted mortality rate of 10.2 deaths/1000 persons-week at risk. The fact that nearly half of all deaths clustered within the first two weeks after the initiation of the diarrhea episode is strongly supportive of the significance of this pathology starting the chain of events leading to death, rather than other causes. This is also supported by the cause of death diagnoses provided by the verbal autopsies conducted among all study deaths, whereby 39% of all deaths were directly attributed to diarrhea. More work is therefore needed to put in place better preventive and management tools to specifically address diarrhea as a syndrome, and also all its secondary consequences.
This study allows a better understanding of the risk factors and associated problems that require more attention in order to address diarrheal mortality in this particular setting. Independent risk factors for death identified through this analysis are multifactorial and can be subdivided in those specifically related to pathogens, clinical signs and symptoms, or to concurrent co-morbidities.
Cryptosporidium was the only pathogen independently associated with death, even after adjusting for covariates. This pathogen had already been identified in the GEMS study as one of the major culprits of excess mortality among young children with MSD, being Manhiça the GEMS site where the highest overall prevalence was found , whereby cryptosporidium represented the second most frequent pathogens associated with MSD . Infection from this protozoan, that can affect both immunocompetent but also immunocompromised hosts, often derives from contact with potentially contaminated water or food sources. In the absence of a Cryptosporidium vaccine and considering the ineffectiveness of current treatment strategies in the youngest populations and the diagnostic challenges posed by this pathogen, minimizing the overall environmental burden (contaminated water/food) needs prioritization. Additionally, and given the high prevalence of HIV infection in the district of Manhiça , prompt initiation of antiretroviral therapy among those children confirmed to be HIV-infected may contribute to a reduction in the burden of Cryptosporidium infections, understanding that this is another well characterized risk factor for cryptosporidiosis .
The presence of an Astrovirus infection appeared as an independent risk factor for death in the multivariate analysis restricted to children with available HIV data, the association being borderline significant (HR 3.5, p = 0.0545) in the univariate analysis of the entire population. Although numbers are relatively small, more attention should be drawn to this particular pathogen, that has already demonstrated a substantial disease burden in neighboring countries , and for which neither curative nor preventive strategies have been described.
Although rotavirus-associated mortality has clearly been established globally , in this series rotavirus infections appeared to protect against death in the univariate analysis (HR = 0.5, p = 0.0198), highlighting overall good course secondary to this pathogen if adequate hydration measures are promptly started at the community of health system level.
Being admitted to hospital appeared as a risk factor for a poorer outcome. This is hardly surprising, considering that those patients with MSD requiring admission are those with higher severity, and with a major risk of co-morbidities , and that health care utilization in cases of MSD in this study was high. In terms of signs and symptoms detected in patients upon admission associated with a higher risk of death, those related to an altered consciousness and dehydration were most relevant, together with those suggestive of respiratory distress and/or nutritional deficiencies. These signs/symptoms reflect both the real mechanisms by which a diarrheal episode can become lethal and the typical complications of MSD. Indeed, guaranteeing an adequate intake of liquids and replacing excessively lost water and electrolytes during episodes of diarrhea remains the cornerstone strategy to prevent diarrheal associated deaths, and was at the basis of the WHO and UNICEF’s recommendation of promoting Oral Rehydration Solution (ORS) usage back in the 1970’s . However, diarrheal episodes can trigger both secondary infections or the exacerbation of pre-existing conditions, such as malnutrition, as historically and extensively demonstrated in the literature  and also as confirmed in this study, where malnutrition was one of the strongest risk factors associated with death. Breaking the vicious circle that can become established between malnutrition and diarrheal diseases will therefore be critical to address a significant proportion of all MSD deaths in settings such as Manhiça.
Breast-feeding, especially if this is the only source of nutrition, has been shown to protect children against diarrhea in Africa as elsewhere in the developing world . Previous analyses from this same study showed that in Manhiça, when controlled by other confounders, the particular breastfeeding pattern did not seem to protect against pediatric MSD  . Interestingly, and in accordance with the published evidence , this analysis has shown that breastfeeding (either partial or exclusive) during the first months of life has a significant protective effect against diarrheal attributable mortality. Our findings reinforce the current WHO recommendations that encourage exclusive breastfeeding during the first 6 months of life as a vital strategy to prevent infant mortality [44, 45].
A concomitant diagnosis of pneumonia (HR = 3.50, P < 0.0001) or invasive bacterial disease (HR = 6.79, P = 0.0009), were both clear and independent risk factors for death, signaling the importance of co-morbidities as triggers of a fatal outcome, and of concomitant antibiotic treatment as a potential justifiable strategy in patients with MSD which do not improve rapidly with adequate rehydration. Importantly, an admission diagnosis of dysentery was shown in this analysis to be protective against death in the univariate analysis (HR = 0.3, P = 0.0112). Typical pathogens causing dysentery such as Shigella or Entamoeba histolytica have been shown to cause longer and more severe diarrheal episodes, with an overall increased risk of associated mortality [46, 47]. Our somehow counterintuitively finding may be explained by the early and aggressive administration of antimicrobials in cases of diarrhea with blood in our setting, which may have attenuated the impact of this particular type of diarrheal disease.
In this series, MSD episodes among HIV-infected individuals carried a much worse prognosis. Although not even a fourth of all patients had HIV screening results available, due to the late initiation of this practice during the study, the prevalence of HIV-infected individuals among recruited patients was high (21%), and this infection was clearly recognized as an independent risk factor in the multivariate analysis restricted to those with HIV data available. Indeed, HIV infection has been shown to play a significant role not only in the burden and associated mortality of MSD, but also in the etiology and epidemiology of diarrheal disease among children in Manhiça . In this respect, in settings such as Manhiça with a high HIV prevalence, aggressive and prompt initiation of antiretrovirals among newly diagnosed HIV cases should be the pillar upon which basing any strategy destined to reducing the burden of diarrheal deaths. Similarly, the suitability of exclusive breastfeeding practices for this particular vulnerable group should also be strengthened .
This study has several limitations, starting with the aforementioned lack of HIV data for all study participants, which would have allowed a more precise characterization of the risk factors for MSD mortality. This happened because according to the national guidelines, HIV testing was not mandatory and required an additional patient consent; leading to the high proportion of participants without HIV results. However, the comparison of patients with or without HIV results available showed no differences (data not shown). This study failed to evaluate electrolyte or zinc levels among MSD cases and controls, which could have contributed to a better understanding of dehydration patterns or to the enrichment of the evidence supporting zinc supplementation in children with diarrhea. However, several previous studies have demonstrated the protective effect of this simple intervention against diarrhea mortality [50, 51] and its use remains warranted, something that was not in place in Mozambique during the study period. In addition, there was also a potential for recall bias during the administration of verbal autopsy forms in this study. Causes of death based on verbal autopsies were determined by physicians, although this method also has inherent limitations for its subjectivity and low reliability as described elsewere [52,53,54,55]. Finally, this analysis was conducted with GEM’s original set of microbiological data, and more recent re-analyses of the study samples, using more advanced molecular methods, have evidenced a different distribution of the pathogens associated to MSD, something that could have potentially altered our own analysis.
By characterizing diarrhea-associated deaths in a large group of Mozambican children, this study contributed to enlarging the pool of evidence regarding risk factors for mortality among pediatric patients with MSD. The knowledge of these factors can better define those children at higher risk of dying as a result of a diarrheal episode, and thus allow a more targeted and specific approach to prevent those deaths. In this setting, targeted measures against Cryptosporidium infection and aggressive management of dehydration, malnutrition and suspected co-morbidities, together with the promotion of continuous breastfeeding may help prevent most deaths due to diarrhea.
- atypical EPEC:
atypical enteropathogenic Escherichia coli
Cause of death
Data Collection Centre in the United States of America
Demographic Surveillance System
Enteroaggregative Escherichia coli
Global Enteric Multi-center Study
The International Network for the Demographic Evaluation of Populations and Their Health
Low and middle-income countries
Enterotoxigenic Escherichia coli producing heat-labile toxin
Oral rehydration solution
Enterotoxigenic Escherichia coli producing heat-stable toxin
- typical EPEC:
Typical enteropathogenic Escherichia coli
Liu L, Oza S, Hogan D, Chu Y, Perin J, Zhu J, et al. Global, regional, and national causes of under-5 mortality in 2000–15: an updated systematic analysis with implications for the sustainable development goals. Lancet. 2016;388(10063):3027-35.
Global Burden of Disease. Estimates of global, regional, and national morbidity, mortality, and aetiologies of diarrhoeal diseases. a systematic analysis for the Global Burden of Disease Study 2015. Lancet Infect Dis. 2017;17:909–48.
Munos MK, Fischer Walker CL, Black RE. The effect of oral rehydration solution and recommended home fluids on diarrhoea mortality. Int J Epidemiol. 2010;39(Suppl 1):i75-87.
Burnett E, Jonesteller CL, Tate JE, Yen C, Parashar UD. Global impact of rotavirus vaccination on childhood hospitalizations and mortality from diarrhea. J Infect Dis. 2017;215(11):1666-72.
Humphreys EH, Smith NA, Azman H, McLeod D, Rutherford GW. Prevention of diarrhoea in children with HIV infection or exposure to maternal HIV infection. Cochrane Database Syst Rev. 2010;(6):CD008563.
Boschi-Pinto C, Velebit L, Shibuya K. Estimating child mortality due to diarrhoea in developing countries. Bull World Health Organ. 2008;86(9):710-7.
Kotloff KL, Nataro JP, Blackwelder WC, Nasrin D, Farag TH, Panchalingam S, et al. Burden and aetiology of diarrhoeal disease in infants and young children in developing countries (the global enteric multicenter study, GEMS): a prospective, case-control study. Lancet. 2013;382:209–22.
van Eijk AM, Brooks JT, Adcock PM, Garrett V, Eberhard M, Rosen DH, et al. Diarrhea in children less than two years of age with known HIV status in Kisumu. Kenya Int J Infect Dis. 2010;14:e220–5. https://doi.org/10.1016/j.ijid.2009.06.001.
Groome MJ, Madhi SA. Five-year cohort study on the burden of hospitalisation for acute diarrhoeal disease in African HIV-infected and HIV-uninfected children: potential benefits of rotavirus vaccine. Vaccine. 2012;30(Suppl 1):A173-8.
Chhagan MK, Kauchali S. Comorbidities and mortality among children hospitalized with diarrheal disease in an area of high prevalence of human immunodeficiency virus infection. Pediatr Infect Dis J. 2006;25(4):333-8.
Nhampossa T, Sigaúque B, MacHevo S, MacEte E, Alonso P, Bassat Q, et al. Severe malnutrition among children under the age of 5 years admitted to a rural district hospital in southern Mozambique. Public Health Nutr. 2013;16(9):1565-74.
Talbert A, Thuo N, Karisa J, Chesaro C, Ohuma E, Ignas J, et al. Diarrhoea complicating severe acute malnutrition in kenyan children: a prospective descriptive study of risk factors and outcome. PLoS One. 2012;7:1–8.
el Samani EF, Willett WC, Ware JH. Association of malnutrition and diarrhea in children aged under five years. A prospective follow-up study in a rural Sudanese community. Am J Epidemiol. 1988;128(1):93-105.
Dgedge M, Novoa A, Macassa G, Sacarlal J, Black J, Michaud C, et al. The burden of disease in Maputo City, Mozambique: registered and autopsied deaths in 1994. Bull World Health Organ. 2001;79:546–52.
Sigaúque B, Roca A, Mandomando I, Morais L, Quintó L, Sacarlal J, et al. Community-acquired bacteremia among children admitted to a rural hospital in Mozambique. Pediatr Infect Dis J. 2009;28:108–13.
Mandomando I, MacEte E, Sigaúque B, Morais L, Quintó L, Sacarlal J, et al. Invasive non-typhoidal Salmonella in Mozambican children. Trop Med Int Health. 2009;14(12):1467-74.
Mandomando IM, Macete EV, Ruiz J, Sanz S, Abacassamo F, Vallès X, et al. Etiology of diarrhea in children younger than 5 years of age admitted in a rural hospital of southern Mozambique. Am J Trop Med Hyg. 2007;76(3):522-7.
Nhampossa T, Mandomando I, Acacio S, Quintó L, Vubil D, Ruiz J, et al. Diarrheal disease in rural Mozambique: burden, risk factors and etiology of diarrheal disease among children aged 0-59 months seeking care at health facilities. PLoS One. 2015;10(5):e0119824.
González R, Augusto OJ, Munguambe K, Pierrat C, Pedro EN, Sacoor C, et al. HIV incidence and spatial clustering in a rural area of southern Mozambique. PLoS One. 2015;10(7):e0132053.
Sigaúque B, Vubil D, Sozinho A, Quintó L, Morais L, Sacoor C, et al. Haemophilus influenzae type b disease among children in rural Mozambique: impact of vaccine introduction. J Pediatr. 2013;163(1 Suppl):S19-24.
Sacoor C, Nhacolo A, Nhalungo D, Aponte JJ, Bassat Q, Augusto O, et al. Profile: Manhica Health Research Centre (Manhica HDSS). Int J Epidemiol. 2013;42:1309–18. https://doi.org/10.1093/ije/dyt148.
Roca A, Quintó L, Abacassamo F, Morais L, Vallès X, Espasa M, et al. Invasive Haemophilus influenzae disease in children less than 5 years of age in Manhiça, a rural area of southern Mozambique. Trop Med Int Health. 2008;13(6):818-26
Guinovart C, Bassat Q, Sigaúque B, Aide P, Sacarlal J, Nhampossa T, et al. Malaria in rural Mozambique. Part I: children attending the outpatient clinic. Malar J. 2008;7:1–9.
Bassat Q, Machevo S, O’Callaghan-Gordo C, Sigaúque B, Morais L, Díez-Padrisa N, et al. Distinguishing malaria from severe pneumonia among hospitalized children who fulfilled integrated management of childhood illness criteria for both diseases: a hospital-based study in Mozambique. Am J Trop Med Hyg. 2011;85:626–34.
Sigaúque B, Roca A, Mandomando I, Morais L, Quintó L, Sacarlal J, et al. Community-acquired bacteremia among children admitted to a rural hospital in Mozambique. Pediatr Infect Dis J. 2009;28(2):108-13.
Roca A, Bassat Q, Morais L, Machevo S, Sigauque B, O’Callaghan C, et al. Surveillance of acute bacterial meningitis among children admitted to a district hospital in rural Mozambique. Clin Infect Dis. 2009;48(Suppl 2):S172-80.
Kotloff KL, Blackwelder WC, Nasrin D, Nataro JP, Farag TH, van Eijk A, et al. The global enteric multicenter study (GEMS) of diarrheal disease in infants and young children in developing countries: epidemiologic and clinical methods of the case/control study. Clin Infect Dis. 2012;55(Suppl 4):S232-45.
StataCorp. Stata Statistical Software: Release 13.1. College Station, TX: StataCorp LP. 2013.
Ploner M, Heinze G. Coxphf:Cox regression with Firth’s penalized likehood. R package version; 2015.
Team RC. R: A linguage and enviroment for statistical computing. Vienna: R Foundation for Statistical computing; 2015.
Firth D. Bias reduction of maximum likelihood estimates. Biometrika. 1993;80:27–38.
Heinze G, Schemper M. A solution to the problem of separation in logistic regression. Stat Med. 2002;21:2409–19.
Sommerfelt H, Steinsland H, van der Merwe L, Blackwelder WC, Nasrin D, Farag TH, et al. Case/control studies with follow-up: constructing the source population to estimate effects of risk factors on development, disease, and survival. Clin Infect Dis. 2012;55(Suppl 4):S262-70.
Lozano R, Naghavi M, Foreman K, Lim S, Shibuya K, Aboyans V, et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: A systematic analysis for the global burden of disease study 2010. Lancet. 2012;380(9859):2095-128.
Sow SO, Muhsen K, Nasrin D, Blackwelder WC, Wu Y, Farag TH, et al. The burden of cryptosporidium diarrheal disease among children < 24 months of age in moderate/high mortality regions of sub-Saharan Africa and South Asia, utilizing data from the global enteric multicenter study (GEMS). PLoS Negl Trop Dis. 2016;10:1–20.
Wang ZD, Liu Q, Liu HH, Li S, Zhang L, Zhao YK, et al. Prevalence of cryptosporidium, microsporidia and Isospora infection in HIV-infected people: a global systematic review and meta-analysis. Parasit Vectors. 2018;11(1):28.
Shioda K, Cosmas L, Audi A, Gregoricus N, Vinjé J, Parashar UD, et al. Population-based incidence rates of diarrheal disease associated with norovirus, sapovirus, and astrovirus in Kenya. PLoS One. 2016;11(4):e0145943.
Tate JE, Burton AH, Boschi-Pinto C. Parashar UD. Global, Regional, and National Estimates of Rotavirus Mortality in Children <5 Years of Age, 2000–2013. Clin Infect Dis. 2016;62(Suppl 2):S96-S105.
Nhampossa T, Mandomando I, Acacio S, Nhalungo D, Sacoor C, Nhacolo A, et al. Health care utilization and attitudes survey in cases of moderate-to-severe diarrhea among children ages 0-59 months in the district of manhica, southern Mozambique. Am J Trop Med Hyg. 2013;89(1 Suppl):41-48.
WHO. Pocket Book of Hospital Care for Children. Guidelines for the Management of Common Childhood Illnesses; 2013. https://doi.org/10.1016/j.cardfail.2011.02.010.
Ashworth A. International differences in child mortality and the impact of malnutrition. Hum Nutr Clin Nutr. 1982;36(4):279-88.
Lamberti LM, Fischer Walker CL, Noiman A, Victora C, Black RE. Breastfeeding and the risk for diarrhea morbidity and mortality. BMC Public Health. 2011;11(Suppl 3):S15.
Sankar MJ, Sinha B, Chowdhury R, Bhandari N, Taneja S, Martines J, et al. Optimal breastfeeding practices and infant and child mortality: a systematic review and meta-analysis. Acta Paediatr. 2015;104:3–13. https://doi.org/10.1111/apa.13147.
Habicht JP. Expert consultation on the optimal duration of exclusive breastfeeding - the process, recommendations, and challenges for the future. Adv Exp Med Biol. 2004;554:79-87.
World Health Organization, United Nations Children’s Fund. Guideline: updates on HIV and infant feeding: the duration of breastfeeding, and support from health services to improve feeding practices among mothers living with HIV. Geneva: World Health Organization; 2016.
Kotloff KL, Riddle MS, Platts-Mills JA, Pavlinac P, Zaidi AKM. Shigellosis. Lancet. 2018;391(10122):801-12.
Pires SM, Fischer-Walker CL, Lanata CF, Devleesschauwer B, Hall AJ, Kirk MD, et al. Aetiology-specific estimates of the global and regional incidence and mortality of diarrhoeal diseases commonly transmitted through food. PLoS One. 2015;10(12):e0142927.
Acácio S, Nhampossa T, Quinto L, Vubil D, Sacoor C, Kotloff K, et al. The role of HIV infection in the etiology and epidemiology of diarrheal disease among children aged 0-59 months in Manhica district. Int J Infect Dis. 2018;73:10-17.
WHO (World Health Organization), UNICEF, UNAIDS, UNFPA. Guidelines on HIV and Infant Feeding: Principles and recommendations for infant feeding in the context of HIV and summary of evidence. 2010.
Sazawal S, Black RE, Ramsan M, Chwaya HM, Dutta A, Dhingra U, et al. Effect of zinc supplementation on mortality in children aged 1-48 months: a community-based randomised placebo-controlled trial. Lancet. 2007;369(9565):927-34.
Yakoob MY, Theodoratou E, Jabeen A, Imdad A, Eisele TP, Ferguson J, et al. Preventive zinc supplementation in developing countries: impact on mortality and morbidity due to diarrhea, pneumonia and malaria. BMC Public Health. 2011;11(Suppl 3):S23.
WHO, Department of Measurement and Health Information Systems evidence and information for policy. WHO Technical Consultation on Verbal Autopsy Tools: Review of the literature and currently-used verbal autopsy tools; 2004. p. 1–45. http://www.who.int/healthinfo/statistics/mort_verbalautopsy.pdf.
Butler D. Verbal autopsy methods questioned. Nature. 2010;467(7319):1015.
Sacoor C, Payne B, Augusto O, Vilanculo F, Nhacolo A, Vidler M, et al. Health and socio-demographic profile of women of reproductive age in rural communities of southern Mozambique. PLoS One. 2018;13(2):e0184249.
Liu J, Platts-Mills JA, Juma J, Kabir F, Nkeze J, Okoi C, et al. Use of quantitative molecular diagnostic methods to identify causes of diarrhoea in children: a reanalysis of the GEMS case-control study. Lancet. 2016;388(10051):1291-301.
The authors would like to thank the participants and their parents in this study for allowing collecting data. We thank also Drs. Abel Nhama, Kizito Gondo, Sónia Machevo and all medical officers for clinical management, field workers and other staff from CISM and Manhiça District Hospital for collaborating on data collection and processing. Special thanks to João Campos Mucasse and Juvêncio Joaquim for excellent data management.
The study was supported by the Bill & Melinda Gates Foundation (GEMS: Project OPP 38874). ISGlobal is a member of the CERCA Programme, Generalitat de Catalunya (http://cerca.cat/en/suma/). CISM is supported by the Government of Mozambique and the Spanish Agency for International Development (AECID). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The datasets generated during and/or analyzed during the current study are not publicly available due to confidentiality reasons but are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study, conducted under the framework of the GEMS project, took place in seven LMICs and was coordinated by the Center for Vaccine Development (Baltimore, Maryland, USA) in collaboration with the Barcelona institute of Global Health (ISGlobal), in Spain. The GEMS protocol and informed consent were both approved by the National Bioethics Committee for Health of Mozambique (CNBS – IRB00002657), the ethics committee of the Hospital Clínic of Barcelona, and The Institutional Review Board at the University of Maryland. After informing the objectives and characteristics of the study a written informed consent was obtained from the child’s primary caretaker, and a witnessed consent was obtained for illiterate parents/caretakers.
Consent for publication
The authors declare that they have no Competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Main cause of death diagnoses among study participants as determined by verbal autopsy. (DOCX 14 kb)