Introduction

Rigorous public health efforts resulted in a significant decline in mortality due to diarrheal diseases (DD) over the past 20 years. However, these diseases continue to cause a major global disease burden, especially in children under 5 years of age. In addition, the incidence of childhood diarrhea in low-income countries has not declined as rapidly as mortality (GBD 2015 Risk Factors Collaborators 2016). In the most recent Global Burden of Disease (GBD) study, DD was the fourth-leading cause of death among children under 5, responsible for 499,000 deaths (95% UI: 447,000–558,000), representing 8.6% of all deaths in this age group (GBD 2015 Risk Factors Collaborators 2016). For those who survive these illnesses and suffer from repeated infections by enteric pathogens during the critical early years of life; DD can lead to serious, lifelong health consequences such as environmental enteric dysfunction (EED), growth faltering, impaired cognitive development, and reduced immune response to infection and vaccinations (Guerrant et al. 2013).

DD pathogen etiologic contribution may vary depending on the study’s geographic location, duration, or the population sampled (Lindsay et al. 2015). These infections are believed to be different in the developing world compared to the developed world with regard to a number of features, including earlier age of onset, multiple repeated exposures, greater diversity of pathogens, nutritional status of the host, and a number of others, such as co-infection, diet, and genetics (Heidt et al. 2014).

The Eastern Mediterranean Region (EMR) is home to more than 500 million people, representing a diverse group of 22 countries. EMR countries have diverse historical backgrounds, political and social contexts, and fiscal and cultural influences on their health care systems. The region has wide variation in per capita gross national product (GNP) (The World Bank 2016), which has a major influence on overall health spending and results in substantial health inequities both within and across countries. During recent years DD prevention efforts that focus on vaccines in the short term and improvements in water, sanitation, and hygiene in the long term have been impeded by warfare and political unrest in the region. These conflicts and wars also resulted in a huge problem of internal displacement and refugees (Mokdad et al. 2016).

Many countries in the EMR achieved important successes in the fight against DD in the 1970s and 1980s with the support of United Nations International Children’s Emergency Fund (UNICEF) and the World Health Organization (WHO) through the National Control of Diarrheal Diseases Project (NCDDP) (WHO 1992; Enzley et al. 1997). For example, Egypt’s program, which spanned from 1981 to 1991, was credited with significantly improving diarrheal case management (National Control Of Diarrheal Diseases Project 1988; El-Rafie et al. 1990; Cobb et al. 1996). However, over the last two decades, momentum has slowed (Forsberg et al. 2007). As of the date of this report, rotavirus vaccines have been introduced through National Immunization Programs in only nine countries in the region: Djibouti, Jordan, Libya, Morocco, Qatar, Saudi Arabia, Sudan, United Arab Emirates, and Yemen (PATH 2017).

In this report, we are updating our previous burden estimates (Khalil et al. 2016), pathogen distribution, and risk factors for diarrhea in children and adults in the EMR for 1990–2015.

Methods

The Global Burden of Disease Study (GBD 2015) is a systematic, comprehensive effort to quantify health loss from more than 300 diseases and injuries, including diarrheal diseases and associated risk factors. The GBD estimation strategy, including for diarrheal diseases, has been described in detail elsewhere (Foreman et al. 2012; Flaxman et al. 2015). The burden of diarrheal diseases is measured in deaths, incidence, and disability-adjusted life-years (DALYs), which are the sum of years of life lost (YLLs) and years lived with disability (YLDs). The etiological burden was also estimated for 13 pathogens associated with diarrhea.

All estimates are produced by year and by age, for both sexes, and for all countries. In accordance with the guidelines for accurate and transparent health estimates reporting (GATHER), code for each step of the estimation process is available online on GitHub (http://www.ghdx.healthdata.org/gbd-2015-code) (Institute for Health Metrics and Evaluation). The methods of each of these steps are summarized below.

Study region

The EMR countries were grouped according to per capita gross national income (GNI) into low-income countries (LICs) [Islamic Republic of Afghanistan (Afghanistan), Djibouti, Somalia, Republic of Yemen (Yemen)]; middle-income countries (MICs) [Arab Republic of Egypt (Egypt), Islamic republic of Iran (Iran), Iraq, Jordan, Lebanon, Libya, Morocco, Pakistan, Palestine, Sudan, Syrian Arab Republic (Syria), Tunisia]; and high-income countries (HICs) [Bahrain, Kuwait, Oman, Qatar, Saudi Arabia, and the United Arab Emirates (UAE)]. We defined LICs as those having a per capita GNI of $1045 or less, MICs as those with a per capita GNI between $1046 and $12,735, and HICs as countries with per capita GNI of $12,736 or greater (The World Bank 2016).

Mortality

Cause-specific mortality estimates for diarrheal diseases were modeled using a Bayesian ensemble modeling process (GBD 2015 Mortality and Causes of Death Collaborators 2016). Diarrhea mortality data included vital registration and verbal autopsy sources. The modeling process estimated the mortality rate due to diarrhea for both sexes from 1990 to 2015 for all age groups in every country and subnational regions in select countries. We considered the following covariates: education, lag-dependent income, underweight, latitude, population density, improved water and sanitation sources, diarrhea risk factors summary, Socio-demographic Index (SDI), and rotavirus vaccine coverage. The ensemble model approach allows for a suite of models, weighted by out-of-sample predictive validity, to inform the final estimates.

Morbidity

Diarrheal cases were defined as three or more loose stools in a 24-h period. As with mortality, morbidity was modeled at every year, sex, age, and geographic location in GBD 2015. The morbidity model used DisMod-MR 2.1, a Bayesian meta-analytic, age-integrating, mixed-effects model which is available online on EpiViz (http://www.vizhub/healthdata.org/epi) (Kassebaum et al. 2016). Diarrhea prevalence and incidence data from a systematic literature review, population-representative surveys, and hospital and health care utilization data informed the non-fatal model.

Etiologies

Diarrhea cases and deaths were attributed to pathogens using a counterfactual approach that accounts for exposure to pathogens and for the association between each pathogen and diarrhea. A systematic literature review on the proportion of diarrhea cases that test positive for a set of pathogens was updated for GBD 2015. These data were used in the DisMod-MR framework to estimate the age, sex, year, and geographic distribution of pathogens in diarrheal episodes. The population attributable fraction (PAF) was used to identify the fraction of diarrhea cases and deaths due to each pathogen. The PAF was calculated as: (GBD 2015 Mortality and Causes of Death Collaborators 2016).

$${\text{PAF}} = {\text{Proportion}} \times \left( {1 - \frac{1}{\text{OR}}} \right),$$

where Proportion is the proportion of cases positive for a pathogen and the odds ratio (OR) is the odds of diarrhea given pathogen detection. The odds ratios were from a systematic reanalysis of the Global Enteric Multicenter Study (GEMS) (Kotloff et al. 2013), a multi-site case–control study of moderate-to-severe diarrhea in children under 5 that systematically tested nearly half of the original GEMS samples using a molecular quantitative polymerase chain reaction (qPCR) diagnostic (Liu et al. 2016). A mixed-effects conditional logistic regression model estimated the odds ratios for diarrhea including random effects on site to account for geographic variation.

Since the odds of diarrhea given pathogen presence were calculated using the qPCR diagnostic, we adjusted our proportion estimates to be comparable to qPCR-based estimates. The sensitivity and specificity of the non-molecular diagnostic techniques from GEMS were evaluated compared to the qPCR diagnostic, and these values were used to make this adjustment with this formula (Wickham 2009):

$${\text{Proportion}}_{\text{True}} = \frac{{({\text{Proportion}}_{\text{Observed}} \times {\text{Specificity}} - 1)}}{{({\text{Sensitivity}} + {\text{Specificity}} - 1)}}.$$

Vibrio cholerae and Clostridium difficile were estimated separately from the other pathogens in GBD. Cholera cases were estimated using proportion data from published studies to calculate an expected number of annual cases for each country and year, and those estimates were compared to the World Health Organization case notification data (WHO 2016) to estimate underreporting of cholera. Cholera deaths were estimated using case fatality data in DisMod-MR. Since C. difficile is frequently associated with hospital and health care utilization, hospital incidence data were modeled in DisMod-MR using hospital and health care utilization data with ICD codes for C. difficile.

All etiologies were estimated independently and for each year, geography, age, and sex.

Risk factors

We also assessed diarrheal DALYs, YLLs, and YLDs attributable to childhood stunting, suboptimal breastfeeding, vitamin A deficiency, zinc deficiency, and water, sanitation, and hygiene (WASH). Risk factor attribution follows a general counterfactual approach where the exposure and relative risk of diarrhea were used to estimate the burden of the aforementioned risk factors (GBD 2015 Risk Factors Collaborators 2016).

Socio-demographic Index

We evaluate associations between diarrhea and socio-demographic status using the Socio-demographic Index (SDI). The SDI is a composite measure developed for GBD 2015 that accounts for fertility rate, lag-dependent income per capita, and education (GBD 2015 Mortality and Causes of Death Collaborators 2016). To capture the average relationships for each age–sex group, we applied a simple least squares spline regression of the diarrhea mortality rate on SDI. The predicted diarrhea mortality rates from this regression were used as expected mortality rates based on SDI. The SDI is scaled from 0 to 1 where 0 represents the lowest possible observed SDI and 1 is the highest. SDI in 2015 in the EMR ranged from 0.27 in Somalia to 0.83 in Kuwait.

Results

There were 103,692 diarrhea deaths in the EMR in 2015 (95% uncertainty interval (UI): 87,018–124,692), and the mortality rate was 16.0 deaths per 100,000 persons (95% UI: 13.4–19.2). The majority of these deaths occurred in children under 5 years old (63.3%) (65,670 deaths, 95% UI: 53,640–79,486). Although the greatest number of diarrhea deaths occurred in children under 5, diarrhea mortality was also high in the 70+ year age group. In fact, of the 22 countries in the EMR, the mortality rate was higher in the 70+ age group than the under-5 age group in 17 countries. Diarrhea mortality in the elderly was highest in Somalia (1695 per 100,000, 95% UI: 709–2896), followed by Djibouti and Pakistan (Table 1). Diarrhea mortality decreased over time in the EMR (Table 2; Fig. 1). Between 1990 and 2015, the number of diarrhea deaths decreased by 54% among all ages (95% UI: 43–62%) and 65% among children under 5 (95% UI: 55–73%). The fastest rate of decrease among children under 5 occurred in Iran (97%, 95% UI: 93–99%) and Syria (97%, 95% UI: 94–99%) and the slowest occurred in Somalia (11%, 95% UI: −40 to 45%) and Qatar (33%, 95% UI: −65 to 72%).

Table 1 The number of diarrhea deaths and mortality rate (per 100,000) in 2015 in the Eastern Mediterranean Region for each age group and sex (Global Burden of Disease 2015 study, Eastern Mediterranean Region, 2015)
Table 2 The number of deaths, disability-adjusted life-years (DALYs), and episodes of diarrhea in the Eastern Mediterranean Region and the diarrhea mortality (per 100,000) and incidence (per person-year) rates among under-5 and all ages in 1990, 2000, and 2015 (Global Burden of Disease 2015 study, Eastern Mediterranean Region, 1990–2015)
Fig. 1
figure 1

Diarrhea mortality rate per 100,000 population, by age, 1990–2015. The mortality rate in the Eastern Mediterranean Region is shown for six age groups over time from 1990 to 2015. Points represent mean estimates and error bars represent 95% uncertainty intervals (Global Burden of Disease 2015 study, Eastern Mediterranean Region, 1990–2015)

Diarrhea mortality among children under 5 varied by time and country. The under-5 diarrhea mortality rate was highest in Somalia, followed by Pakistan, Sudan, and Afghanistan (Figs. 2, 3). Due to its high population, the greatest number of under-5 deaths occurred in Pakistan (55,500 deaths, 95% UI: 43,258–70,027). Diarrhea mortality was associated with trends in the Socio-demographic Index in most geographies (Figs. 2b, 3), but the mortality rate was much lower in Palestine and Iraq than would be expected based on SDI alone. The ratio of the observed mortality rate to the expected mortality rate based on SDI alone was lowest in Syria and Palestine, where mortality was 3% of the expected value, and highest in Bahrain, where the mortality was 65% higher than expected based on SDI alone. Syria and Palestine have moderate SDI and very low diarrhea mortality rates.

Fig. 2
figure 2

Diarrhea mortality among children under 5 in the Eastern Mediterranean Region, 2015. a The diarrhea mortality rate per 100,000 among children under 5 in 2015. b Ratio of the observed under 5 diarrhea mortality rate in 2015 to the expected mortality rate based on the Socio-demographic Index only. Values below 1 indicate that the diarrhea mortality rate is lower than would be expected based on the global relationship between mortality and SDI, and values above 1 indicate higher mortality rates than would be expected (Global Burden of Disease 2015 study, Eastern Mediterranean Countries, 2015)

Fig. 3
figure 3

The diarrhea mortality rate per 100,000 for the Eastern Mediterranean Region by country among children under 5 in 1990, 2000, and 2015. Countries are ordered from top to bottom from lowest Socio-demographic Index (Somalia) to highest SDI (Qatar) based on SDI values in 2015. Error bars illustrate the 95% uncertainty interval for the estimates (Global Burden of Disease 2015 study, Eastern Mediterranean Countries, 2015)

Diarrhea was also responsible for a large number of illness episodes. In 2015, diarrhea incidence was 2.1 (95% UI: 1.9–2.4) per child-year, totaling more than 171 million episodes across the EMR. The incidence was much lower in adults, including those 70+ years old (0.7 per person-year, 95% UI: 0.69–0.75). The case fatality, expressed as the number of deaths over the number of cases, was 0.1% in those 70+ years old, significantly greater than in children under 5 years (0.038%). Case fatality increased non-linearly with incidence (e-Figure 1). Despite the large observed reductions in diarrhea mortality between 1990 and 2015, diarrhea incidence decreased much more marginally (20.7%).

Among diarrheal etiologies, rotavirus was the leading cause of death in 2015 among children under 5 (13,180 deaths, 95% UI: 9807–17,738), followed by Shigella (10,964 deaths) and enterotoxigenic E. coli (ETEC) (6885 deaths) (Table 3 and e-Figure 1). Among all ages, Shigella was the leading cause of diarrheal deaths (19,450 deaths, 95% UI: 10,026–33,996), followed by Aeromonas and rotavirus. Overall, 99% of diarrheal deaths in 2015 were attributed to at least one etiology (e-Figure 2).

Table 3 The number of diarrhea deaths and the population attributable fraction (PAF) for 13 diarrheal etiologies among under-5 and all ages in the Eastern Mediterranean Region, 2015 (Global Burden of Disease 2015 study, Eastern Mediterranean Region, 2015)

Unsafe water, sanitation, and hygiene were responsible for over 95% of diarrhea DALYs in the EMR in 2015 (95.1%, 95% UI: 89.2–98.1%), which is comparable with the global total in 2015. The proportion of diarrhea DALYs due to unsafe WASH was lowest in Jordan (87.0%) and highest in Sudan (97.0%) (e-Figure 3). Childhood undernutrition was the second-leading risk factor for diarrheal DALYs, responsible for 74.0% (95% UI: 67.4–79.1%), which is significantly higher than the global total (59.8%, 95% UI: 56.2–63.0%). Suboptimal breastfeeding was responsible for 39.5% of diarrhea DALYs among children under 5 (95% UI: 27.1–52.3%). Vitamin A and zinc deficiency were responsible for less than 10% of diarrhea DALYs among children under 5 (8 and 5%, respectively) and were lower than the global average for these risk factors (12.9 and 6.5%, respectively). Childhood undernutrition PAF ranged from 41.1% in the United Arab Emirates to 86.2% in Sudan. Diarrhea standardized exposure variables (SEV), a measure of the prevalence-weighted risk of diarrhea, was highest in Sudan, Afghanistan, and Yemen and lowest in the UAE, Jordan, and Kuwait.

Discussion

Our study is the most comprehensive assessment of diarrheal disease burden and the contributions of specific pathogens and risk factors in the EMR to date. In 2015, the estimated diarrhea-associated deaths and DALYs showed a slight decline from the GBD 2013 estimates of 125,000 and nearly 10 million, respectively. We also continue to find significant variation within the region, with LICs and MICs experiencing social unrest bearing the vast majority of diarrheal burden. Our data clearly illustrate the gross health inequity in the region: the HICs experienced a nominal diarrhea burden compared to the substantial burden in all LICs and some MICs.

There were no major methodological changes in diarrhea mortality modeling between GBD 2013 and GBD 2015 but there were updates to the cause of death and non-fatal data. A major difference between GBD 2013 and 2015 was the introduction of the molecular diagnostic case definition for diarrheal etiologies (Liu et al. 2014; GBD 2015 Mortality and Causes of Death Collaborators 2016; GBD 2015 Risk Factors Collaborators 2016; Kassebaum et al. 2016). In the current round of GBD, we used a systematic reanalysis of the Global Enteric Multicenter Study (GEMS) that retested roughly half of the original 22,000 moderate-to-severe stool samples with quantitative polymerase chain reaction (qPCR) (Liu et al. 2016). In transitioning to a molecular diagnostic case definition, we attribute a much larger proportion of diarrhea cases and deaths to etiologies compared to previous rounds of GBD that used non-molecular diagnostics.

The Socio-demographic Index was designed to be a measure of the socio-demographic status of a country and is well correlated with the Human Development Index (GBD 2015 Mortality and Causes of Death Collaborators 2016). Still, dramatic changes in development such as those that have occurred since the start of political unrest in some of the EMR countries, such as the Syrian civil war, may be missed or at least result in a lag in the components that make up the SDI including income, fertility, and education. We have previously shown that the SDI is a good marker of disease burden, but it may not entirely capture major sociopolitical upheavals (GBD 2015 Mortality and Causes of Death Collaborators 2016).

The introduction of vaccines against diarrheal pathogens may exacerbate inequalities in diarrhea burden. For example, although rotavirus infection was the largest contributor to the diarrheal burden of disease, in some countries in the region, rotavirus vaccine is only available in the private market. This means that wealthier families, who have less need for the vaccine, will gain the primary benefit from its availability. This is troubling because economic analyses of rotavirus vaccine introduction among a number of EMR countries have uniformly suggested that vaccine introduction would be cost-beneficial from a societal perspective (Ortega et al. 2009; Connolly et al. 2012; Javanbakht et al. 2015). One study in Somalia (the only LIC country) suggested that introduction of rotavirus as a special immunization program during a complex humanitarian emergency would meet WHO cost-effectiveness benchmarks (Gargano et al. 2015).

One unique contribution of this analysis is the inclusion of all age groups. Due to the high disease burden in young children, nearly all diarrhea interventions and most diarrhea burden studies are limited to those under 5 years of age (Boschi-Pinto et al. 2008; Fischer Walker et al. 2012; Walker et al. 2013). However, the burden among those over age 70 is substantial, with diarrheal disease-associated deaths totaling nearly one sixth of the number among those under 5. The elderly may face increased diarrhea risk due to immunosenescence and comorbidities, which may also necessitate special consideration in their treatment (Trinh and Prabhakar 2007). The increasing cholera burden in five EMR countries (Jordan, Palestine, Syria, Bahrain, and Oman) is a cause for concern, as neighboring countries remain at high risk of transmission due to the presence and movement of refugee populations among them.

We found no systematic difference in under-5 diarrhea deaths or DALYs when comparing females to males. A hypothesis that there may exists evidence of differential diarrhea mortality burden by sex is supported by a previous analysis of global demographic and health survey (DHS) data. This analysis suggested that reported that girls 1–4 years old, particularly in the Middle Eastern crescent, are at a mortality disadvantage compared to boys (Hill and Upchurch 1995). This could be explained by differences in health care access and nutritional status. Furthermore, an Egyptian study found some evidence that, even when parents sought care for their daughters with diarrhea, regional health care providers provided biased treatment in favor of boys (Yount 2003). In contrast, other studies, such as a 2009 verbal autopsy study in Iraq, found no difference in under-5 mortality by sex (Awqati et al. 2009).

Public health consequences of emergencies in the region

The EMR is facing numerous health challenges as a result of previous wars, and recent revolutions and the wars that followed (Mokdad et al. 2016; Charara et al. 2017) resulted in a huge refugee problem with millions of refugees and massive consequences for the health and well-being of millions of displaced people. The EMR now carries the largest burden of displaced populations globally.

Out of a total of 50 million refugees and internally displaced persons (IDPs) worldwide, more than 29 million (9 million refugees and 20 million IDPs) came from the region (WHO EMRO 2015). Syria is currently the world’s largest source of refugees and IDPs, with more than 40% of the population now displaced both inside the country and in neighboring states. Afghanistan and Somalia face two of the longest-spanning refugee situations, with Afghanis constituting the second-largest refugee group in the world, and Somalia facing one of the world’s most complex refugee situations (WHO EMRO 2015). Over the past few years, the region saw massive internal displacement in Iraq, with more than 3 million people fleeing their homes since June 2014, and in Yemen, where more than 2.3 million people have been internally displaced since March 2015 (ACAPS 2015).

The impact of these emergencies on public health is profound and enduring, affecting both the displaced populations and host communities. The risk of DD is increased in these settings due to limited access to safe water, rotavirus vaccines, and oral rehydration salts, along with other factors like infectious diseases, acute malnutrition, and inappropriate infant and young child feeding. The lack of safe drinking water as well as adequate sanitation and hygiene is especially concerning as a major risk factor for DD. Provision of a safe water supply, sanitation, and hygiene is a crucial priority, even in these emergency situations.

Study limitations

Our estimates of diarrhea mortality, morbidity, and etiology attribution are limited by data availability, and although our modeling process seeks to make use of all available data, the number of relevant publications in the region is limited and unbalanced between countries. The hierarchical modeling approach allows us to “borrow” strength across time and geography to generate the best possible estimates. Second, since we only account for the acute phase of diarrhea in our YLD estimates, the resulting DALYs severely underestimate diarrhea-associated morbidity. In future GBD updates, we expect to include long-term sequelae such as stunting and cognitive impairment, (Moore et al. 2010; MacIntyre et al. 2014; Colombara et al. 2016) which will better estimate the true burden of disease.

Despite these limitations, this analysis also has several strengths. GBD methodology ensures internal consistency so that morbidity and mortality cannot be simultaneously ascribed to competing causes and allows for comparability between countries and across regions.

Conclusion

Increased momentum of public health efforts is needed to reduce the burden of diarrhea in the EMR, especially in lower-income countries and countries experiencing political and social unrest.

Health inequities revealed in our study show that a coordinated approach that involves prevention and treatment is needed to address the multiple causes of diarrheal diseases. Regional health systems need to be strengthened to achieve the widespread availability and use of oral rehydration salts, improved nutrition, better sanitation and hygiene, and increased coverage of rotavirus immunization. Improved rates of breastfeeding should be strongly emphasized as one of the effective tools to prevent DD in countries’ public health programs.

Due to the high relevance of rotavirus infections in the EMR (Malek et al. 2010), there is also an urgent need to accelerate the rollout of rotavirus vaccine in the region through government immunization programs that would ensure access for the children who are most in need. In addition, regional governments should begin deliberation to integrate Shigella, ETEC, and other diarrheal vaccines that are currently in preclinical and clinical trials into their Expanded Programs on Immunization (EPI) as soon as they are approved and licensed.