Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks, 1990-2015: a systematic analysis for the Global Burden of Disease Study 2015

SummaryBackground The Global Burden of Diseases, Injuries, and Risk Factors Study 2015 provides an up-to-date synthesis of the evidence for risk factor exposure and the attributable burden of disease. By providing national and subnational assessments spanning the past 25 years, this study can inform debates on the importance of addressing risks in context. Methods We used the comparative risk assessment framework developed for previous iterations of the Global Burden of Disease Study to estimate attributable deaths, disability-adjusted life-years (DALYs), and trends in exposure by age group, sex, year, and geography for 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks from 1990 to 2015. This study included 388 risk-outcome pairs that met World Cancer Research Fund-defined criteria for convincing or probable evidence. We extracted relative risk and exposure estimates from randomised controlled trials, cohorts, pooled cohorts, household surveys, census data, satellite data, and other sources. We used statistical models to pool data, adjust for bias, and incorporate covariates. We developed a metric that allows comparisons of exposure across risk factors—the summary exposure value. Using the counterfactual scenario of theoretical minimum risk level, we estimated the portion of deaths and DALYs that could be attributed to a given risk. We decomposed trends in attributable burden into contributions from population growth, population age structure, risk exposure, and risk-deleted cause-specific DALY rates. We characterised risk exposure in relation to a Socio-demographic Index (SDI). Findings Between 1990 and 2015, global exposure to unsafe sanitation, household air pollution, childhood underweight, childhood stunting, and smoking each decreased by more than 25%. Global exposure for several occupational risks, high body-mass index (BMI), and drug use increased by more than 25% over the same period. All risks jointly evaluated in 2015 accounted for 57·8% (95% CI 56·6–58·8) of global deaths and 41·2% (39·8–42·8) of DALYs. In 2015, the ten largest contributors to global DALYs among Level 3 risks were high systolic blood pressure (211·8 million [192·7 million to 231·1 million] global DALYs), smoking (148·6 million [134·2 million to 163·1 million]), high fasting plasma glucose (143·1 million [125·1 million to 163·5 million]), high BMI (120·1 million [83·8 million to 158·4 million]), childhood undernutrition (113·3 million [103·9 million to 123·4 million]), ambient particulate matter (103·1 million [90·8 million to 115·1 million]), high total cholesterol (88·7 million [74·6 million to 105·7 million]), household air pollution (85·6 million [66·7 million to 106·1 million]), alcohol use (85·0 million [77·2 million to 93·0 million]), and diets high in sodium (83·0 million [49·3 million to 127·5 million]). From 1990 to 2015, attributable DALYs declined for micronutrient deficiencies, childhood undernutrition, unsafe sanitation and water, and household air pollution; reductions in risk-deleted DALY rates rather than reductions in exposure drove these declines. Rising exposure contributed to notable increases in attributable DALYs from high BMI, high fasting plasma glucose, occupational carcinogens, and drug use. Environmental risks and childhood undernutrition declined steadily with SDI; low physical activity, high BMI, and high fasting plasma glucose increased with SDI. In 119 countries, metabolic risks, such as high BMI and fasting plasma glucose, contributed the most attributable DALYs in 2015. Regionally, smoking still ranked among the leading five risk factors for attributable DALYs in 109 countries; childhood underweight and unsafe sex remained primary drivers of early death and disability in much of sub-Saharan Africa. Interpretation Declines in some key environmental risks have contributed to declines in critical infectious diseases. Some risks appear to be invariant to SDI. Increasing risks, including high BMI, high fasting plasma glucose, drug use, and some occupational exposures, contribute to rising burden from some conditions, but also provide opportunities for intervention. Some highly preventable risks, such as smoking, remain major causes of attributable DALYs, even as exposure is declining. Public policy makers need to pay attention to the risks that are increasingly major contributors to global burden. Funding Bill & Melinda Gates Foundation

Trends in HIV/AIDS morbidity and mortality in Eastern Mediterranean countries, 1990–2015: findings from the Global Burden of Disease 2015 study

OBJECTIVES: We used the results of the Global Burden of Disease 2015 study to estimate trends of HIV/AIDS burden in Eastern Mediterranean Region (EMR) countries between 1990 and 2015. METHODS: Tailored estimation methods were used to produce final estimates of mortality. Years of life lost (YLLs) were calculated by multiplying the mortality rate by population by age-specific life expectancy. Years lived with disability (YLDs) were computed as the prevalence of a sequela multiplied by its disability weight. RESULTS: In 2015, the rate of HIV/AIDS deaths in the EMR was 1.8 (1.4–2.5) per 100,000 population, a 43% increase from 1990 (0.3; 0.2–0.8). Consequently, the rate of YLLs due to HIV/AIDS increased from 15.3 (7.6–36.2) per 100,000 in 1990 to 81.9 (65.3–114.4) in 2015. The rate of YLDs increased from 1.3 (0.6–3.1) in 1990 to 4.4 (2.7–6.6) in 2015. CONCLUSIONS: HIV/AIDS morbidity and mortality increased in the EMR since 1990. To reverse this trend and achieve epidemic control, EMR countries should strengthen HIV surveillance, and scale up HIV antiretroviral therapy and comprehensive prevention services

The global burden of tuberculosis: results from the Global Burden of Disease Study 2015

Background: An understanding of the trends in tuberculosis incidence, prevalence, and mortality is crucial to tracking of the success of tuberculosis control programmes and identification of remaining challenges. We assessed trends in the fatal and non-fatal burden of tuberculosis over the past 25 years for 195 countries and territories. Methods: We analysed 10 691 site-years of vital registration data, 768 site-years of verbal autopsy data, and 361 site-years of mortality surveillance data using the Cause of Death Ensemble model to estimate tuberculosis mortality rates. We analysed all available age-specific and sex-specific data sources, including annual case notifications, prevalence surveys, and estimated cause-specific mortality, to generate internally consistent estimates of incidence, prevalence, and mortality using DisMod-MR 2.1, a Bayesian meta-regression tool. We assessed how observed tuberculosis incidence, prevalence, and mortality differed from expected trends as predicted by the Socio-demographic Index (SDI), a composite indicator based on income per capita, average years of schooling, and total fertility rate. We also estimated tuberculosis mortality and disability-adjusted life-years attributable to the independent effects of risk factors including smoking, alcohol use, and diabetes. Findings: Globally, in 2015, the number of tuberculosis incident cases (including new and relapse cases) was 10·2 million (95% uncertainty interval 9·2 million to 11·5 million), the number of prevalent cases was 10·1 million (9·2 million to 11·1 million), and the number of deaths was 1·3 million (1·1 million to 1·6 million). Among individuals who were HIV negative, the number of incident cases was 8·8 million (8·0 million to 9·9 million), the number of prevalent cases was 8·9 million (8·1 million to 9·7 million), and the number of deaths was 1·1 million (0·9 million to 1·4 million). Annualised rates of change from 2005 to 2015 showed a faster decline in mortality (–4·1% [–5·0 to –3·4]) than in incidence (–1·6% [–1·9 to –1·2]) and prevalence (–0·7% [–1·0 to –0·5]) among HIV-negative individuals. The SDI was inversely associated with HIV-negative mortality rates but did not show a clear gradient for incidence and prevalence. Most of Asia, eastern Europe, and sub-Saharan Africa had higher rates of HIV-negative tuberculosis burden than expected given their SDI. Alcohol use accounted for 11·4% (9·3–13·0) of global tuberculosis deaths among HIV-negative individuals in 2015, diabetes accounted for 10·6% (6·8–14·8), and smoking accounted for 7·8% (3·8–12·0). Interpretation: Despite a concerted global effort to reduce the burden of tuberculosis, it still causes a large disease burden globally. Strengthening of health systems for early detection of tuberculosis and improvement of the quality of tuberculosis care, including prompt and accurate diagnosis, early initiation of treatment, and regular follow-up, are priorities. Countries with higher than expected tuberculosis rates for their level of sociodemographic development should investigate the reasons for lagging behind and take remedial action. Efforts to prevent smoking, alcohol use, and diabetes could also substantially reduce the burden of tuberculosis

Restoration of IFNγR Subunit Assembly, IFNγ Signaling and Parasite Clearance in Leishmania donovani Infected Macrophages: Role of Membrane Cholesterol

Despite the presence of significant levels of systemic Interferon gamma (IFNγ), the host protective cytokine, Kala-azar patients display high parasite load with downregulated IFNγ signaling in Leishmania donovani (LD) infected macrophages (LD-MØs); the cause of such aberrant phenomenon is unknown. Here we reveal for the first time the mechanistic basis of impaired IFNγ signaling in parasitized murine macrophages. Our study clearly shows that in LD-MØs IFNγ receptor (IFNγR) expression and their ligand-affinity remained unaltered. The intracellular parasites did not pose any generalized defect in LD-MØs as IL-10 mediated signal transducer and activator of transcription 3 (STAT3) phosphorylation remained unaltered with respect to normal. Previously, we showed that LD-MØs are more fluid than normal MØs due to quenching of membrane cholesterol. The decreased rigidity in LD-MØs was not due to parasite derived lipophosphoglycan (LPG) because purified LPG failed to alter fluidity in normal MØs. IFNγR subunit 1 (IFNγR1) and subunit 2 (IFNγR2) colocalize in raft upon IFNγ stimulation of normal MØs, but this was absent in LD-MØs. Oddly enough, such association of IFNγR1 and IFNγR2 could be restored upon liposomal delivery of cholesterol as evident from the fluorescence resonance energy transfer (FRET) experiment and co-immunoprecipitation studies. Furthermore, liposomal cholesterol treatment together with IFNγ allowed reassociation of signaling assembly (phospho-JAK1, JAK2 and STAT1) in LD-MØs, appropriate signaling, and subsequent parasite killing. This effect was cholesterol specific because cholesterol analogue 4-cholestene-3-one failed to restore the response. The presence of cholesterol binding motifs [(L/V)-X1–5-Y-X1–5-(R/K)] in the transmembrane domain of IFNγR1 was also noted. The interaction of peptides representing this motif of IFNγR1 was studied with cholesterol-liposome and analogue-liposome with difference of two orders of magnitude in respective affinity (KD: 4.27×10−9 M versus 2.69×10−7 M). These observations reinforce the importance of cholesterol in the regulation of function of IFNγR1 proteins. This study clearly demonstrates that during its intracellular life-cycle LD perturbs IFNγR1 and IFNγR2 assembly and subsequent ligand driven signaling by quenching MØ membrane cholesterol

Estimates of the global, regional, and national morbidity, mortality, and aetiologies of lower respiratory infections in 195 countries, 1990-2016: a systematic analysis for the Global Burden of Disease Study 2016.

BACKGROUND: Lower respiratory infections are a leading cause of morbidity and mortality around the world. The Global Burden of Diseases, Injuries, and Risk Factors (GBD) Study 2016, provides an up-to-date analysis of the burden of lower respiratory infections in 195 countries. This study assesses cases, deaths, and aetiologies spanning the past 26 years and shows how the burden of lower respiratory infection has changed in people of all ages. METHODS: We used three separate modelling strategies for lower respiratory infections in GBD 2016: a Bayesian hierarchical ensemble modelling platform (Cause of Death Ensemble model), which uses vital registration, verbal autopsy data, and surveillance system data to predict mortality due to lower respiratory infections; a compartmental meta-regression tool (DisMod-MR), which uses scientific literature, population representative surveys, and health-care data to predict incidence, prevalence, and mortality; and modelling of counterfactual estimates of the population attributable fraction of lower respiratory infection episodes due to Streptococcus pneumoniae, Haemophilus influenzae type b, influenza, and respiratory syncytial virus. We calculated each modelled estimate for each age, sex, year, and location. We modelled the exposure level in a population for a given risk factor using DisMod-MR and a spatio-temporal Gaussian process regression, and assessed the effectiveness of targeted interventions for each risk factor in children younger than 5 years. We also did a decomposition analysis of the change in LRI deaths from 2000-16 using the risk factors associated with LRI in GBD 2016. FINDINGS: In 2016, lower respiratory infections caused 652 572 deaths (95% uncertainty interval [UI] 586 475-720 612) in children younger than 5 years (under-5s), 1 080 958 deaths (943 749-1 170 638) in adults older than 70 years, and 2 377 697 deaths (2 145 584-2 512 809) in people of all ages, worldwide. Streptococcus pneumoniae was the leading cause of lower respiratory infection morbidity and mortality globally, contributing to more deaths than all other aetiologies combined in 2016 (1 189 937 deaths, 95% UI 690 445-1 770 660). Childhood wasting remains the leading risk factor for lower respiratory infection mortality among children younger than 5 years, responsible for 61·4% of lower respiratory infection deaths in 2016 (95% UI 45·7-69·6). Interventions to improve wasting, household air pollution, ambient particulate matter pollution, and expanded antibiotic use could avert one under-5 death due to lower respiratory infection for every 4000 children treated in the countries with the highest lower respiratory infection burden. INTERPRETATION: Our findings show substantial progress in the reduction of lower respiratory infection burden, but this progress has not been equal across locations, has been driven by decreases in several primary risk factors, and might require more effort among elderly adults. By highlighting regions and populations with the highest burden, and the risk factors that could have the greatest effect, funders, policy makers, and programme implementers can more effectively reduce lower respiratory infections among the world's most susceptible populations. FUNDING: Bill & Melinda Gates Foundation

The global distribution of lymphatic filariasis, 2000–18: a geospatial analysis

Background Lymphatic filariasis is a neglected tropical disease that can cause permanent disability through disruption of the lymphatic system. This disease is caused by parasitic filarial worms that are transmitted by mosquitos. Mass drug administration (MDA) of antihelmintics is recommended by WHO to eliminate lymphatic filariasis as a public health problem. This study aims to produce the first geospatial estimates of the global prevalence of lymphatic filariasis infection over time, to quantify progress towards elimination, and to identify geographical variation in distribution of infection. Methods A global dataset of georeferenced surveyed locations was used to model annual 2000–18 lymphatic filariasis prevalence for 73 current or previously endemic countries. We applied Bayesian model-based geostatistics and time series methods to generate spatially continuous estimates of global all-age 2000–18 prevalence of lymphatic filariasis infection mapped at a resolution of 5 km2 and aggregated to estimate total number of individuals infected. Findings We used 14 927 datapoints to fit the geospatial models. An estimated 199 million total individuals (95% uncertainty interval 174–234 million) worldwide were infected with lymphatic filariasis in 2000, with totals for WHO regions ranging from 3·1 million (1·6–5·7 million) in the region of the Americas to 107 million (91–134 million) in the South-East Asia region. By 2018, an estimated 51 million individuals (43–63 million) were infected. Broad declines in prevalence are observed globally, but focal areas in Africa and southeast Asia remain less likely to have attained infection prevalence thresholds proposed to achieve local elimination. Interpretation Although the prevalence of lymphatic filariasis infection has declined since 2000, MDA is still necessary across large populations in Africa and Asia. Our mapped estimates can be used to identify areas where the probability of meeting infection thresholds is low, and when coupled with large uncertainty in the predictions, indicate additional data collection or intervention might be warranted before MDA programmes cease

Global, regional, and national incidence, prevalence, and mortality of HIV, 1980–2017, and forecasts to 2030, for 195 countries and territories: a systematic analysis for the Global Burden of Diseases, Injuries, and Risk Factors Study 2017

Background Understanding the patterns of HIV/AIDS epidemics is crucial to tracking and monitoring the progress of prevention and control efforts in countries. We provide a comprehensive assessment of the levels and trends of HIV/AIDS incidence, prevalence, mortality, and coverage of antiretroviral therapy (ART) for 1980–2017 and forecast these estimates to 2030 for 195 countries and territories. Methods We determined a modelling strategy for each country on the basis of the availability and quality of data. For countries and territories with data from population-based seroprevalence surveys or antenatal care clinics, we estimated prevalence and incidence using an open-source version of the Estimation and Projection Package—a natural history model originally developed by the UNAIDS Reference Group on Estimates, Modelling, and Projections. For countries with cause-specific vital registration data, we corrected data for garbage coding (ie, deaths coded to an intermediate, immediate, or poorly defined cause) and HIV misclassification. We developed a process of cohort incidence bias adjustment to use information on survival and deaths recorded in vital registration to back-calculate HIV incidence. For countries without any representative data on HIV, we produced incidence estimates by pulling information from observed bias in the geographical region. We used a re-coded version of the Spectrum model (a cohort component model that uses rates of disease progression and HIV mortality on and off ART) to produce age-sex-specific incidence, prevalence, and mortality, and treatment coverage results for all countries, and forecast these measures to 2030 using Spectrum with inputs that were extended on the basis of past trends in treatment scale-up and new infections. Findings Global HIV mortality peaked in 2006 with 1·95 million deaths (95% uncertainty interval 1·87–2·04) and has since decreased to 0·95 million deaths (0·91–1·01) in 2017. New cases of HIV globally peaked in 1999 (3·16 million, 2·79–3·67) and since then have gradually decreased to 1·94 million (1·63–2·29) in 2017. These trends, along with ART scale-up, have globally resulted in increased prevalence, with 36·8 million (34·8–39·2) people living with HIV in 2017. Prevalence of HIV was highest in southern sub-Saharan Africa in 2017, and countries in the region had ART coverage ranging from 65·7% in Lesotho to 85·7% in eSwatini. Our forecasts showed that 54 countries will meet the UNAIDS target of 81% ART coverage by 2020 and 12 countries are on track to meet 90% ART coverage by 2030. Forecasted results estimate that few countries will meet the UNAIDS 2020 and 2030 mortality and incidence targets. Interpretation Despite progress in reducing HIV-related mortality over the past decade, slow decreases in incidence, combined with the current context of stagnated funding for related interventions, mean that many countries are not on track to reach the 2020 and 2030 global targets for reduction in incidence and mortality. With a growing population of people living with HIV, it will continue to be a major threat to public health for years to come. The pace of progress needs to be hastened by continuing to expand access to ART and increasing investments in proven HIV prevention initiatives that can be scaled up to have population-level impact

Global, regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990-2015:a systematic analysis for the Global Burden of Disease Study 2015

Background Non-fatal outcomes of disease and injury increasingly detract from the ability of the world's population to live in full health, a trend largely attributable to an epidemiological transition in many countries from causes affecting children, to non-communicable diseases (NCDs) more common in adults. For the Global Burden of Diseases, Injuries, and Risk Factors Study 2015 (GBD 2015), we estimated the incidence, prevalence, and years lived with disability for diseases and injuries at the global, regional, and national scale over the period of 1990 to 2015.Methods We estimated incidence and prevalence by age, sex, cause, year, and geography with a wide range of updated and standardised analytical procedures. Improvements from GBD 2013 included the addition of new data sources, updates to literature reviews for 85 causes, and the identification and inclusion of additional studies published up to November, 2015, to expand the database used for estimation of non-fatal outcomes to 60 900 unique data sources. Prevalence and incidence by cause and sequelae were determined with DisMod-MR 2.1, an improved version of the DisMod-MR Bayesian meta-regression tool first developed for GBD 2010 and GBD 2013. For some causes, we used alternative modelling strategies where the complexity of the disease was not suited to DisMod-MR 2.1 or where incidence and prevalence needed to be determined from other data. For GBD 2015 we created a summary indicator that combines measures of income per capita, educational attainment, and fertility (the Socio-demographic Index [SDI]) and used it to compare observed patterns of health loss to the expected pattern for countries or locations with similar SDI scores.Findings We generated 9.3 billion estimates from the various combinations of prevalence, incidence, and YLDs for causes, sequelae, and impairments by age, sex, geography, and year. In 2015, two causes had acute incidences in excess of 1 billion: upper respiratory infections (17.2 billion, 95% uncertainty interval [UI] 15.4-19.2 billion) and diarrhoeal diseases (2.39 billion, 2.30-2.50 billion). Eight causes of chronic disease and injury each affected more than 10% of the world's population in 2015: permanent caries, tension-type headache, iron-deficiency anaemia, age-related and other hearing loss, migraine, genital herpes, refraction and accommodation disorders, and ascariasis. The impairment that affected the greatest number of people in 2015 was anaemia, with 2.36 billion (2.35-2.37 billion) individuals affected. The second and third leading impairments by number of individuals affected were hearing loss and vision loss, respectively. Between 2005 and 2015, there was little change in the leading causes of years lived with disability (YLDs) on a global basis. NCDs accounted for 18 of the leading 20 causes of age-standardised YLDs on a global scale. Where rates were decreasing, the rate of decrease for YLDs was slower than that of years of life lost (YLLs) for nearly every cause included in our analysis. For low SDI geographies, Group 1 causes typically accounted for 20-30% of total disability, largely attributable to nutritional deficiencies, malaria, neglected tropical diseases, HIV/AIDS, and tuberculosis. Lower back and neck pain was the leading global cause of disability in 2015 in most countries. The leading cause was sense organ disorders in 22 countries in Asia and Africa and one in central Latin America; diabetes in four countries in Oceania; HIV/AIDS in three southern sub-Saharan African countries; collective violence and legal intervention in two north African and Middle Eastern countries; iron-deficiency anaemia in Somalia and Venezuela; depression in Uganda; onchoceriasis in Liberia; and other neglected tropical diseases in the Democratic Republic of the Congo.Interpretation Ageing of the world's population is increasing the number of people living with sequelae of diseases and injuries. Shifts in the epidemiological profile driven by socioeconomic change also contribute to the continued increase in years lived with disability (YLDs) as well as the rate of increase in YLDs. Despite limitations imposed by gaps in data availability and the variable quality of the data available, the standardised and comprehensive approach of the GBD study provides opportunities to examine broad trends, compare those trends between countries or subnational geographies, benchmark against locations at similar stages of development, and gauge the strength or weakness of the estimates available. Copyright (C) The Author(s). Published by Elsevier Ltd.</p

Mapping geographical inequalities in oral rehydration therapy coverage in low-income and middle-income countries, 2000-17

Background Oral rehydration solution (ORS) is a form of oral rehydration therapy (ORT) for diarrhoea that has the potential to drastically reduce child mortality; yet, according to UNICEF estimates, less than half of children younger than 5 years with diarrhoea in low-income and middle-income countries (LMICs) received ORS in 2016. A variety of recommended home fluids (RHF) exist as alternative forms of ORT; however, it is unclear whether RHF prevent child mortality. Previous studies have shown considerable variation between countries in ORS and RHF use, but subnational variation is unknown. This study aims to produce high-resolution geospatial estimates of relative and absolute coverage of ORS, RHF, and ORT (use of either ORS or RHF) in LMICs. Methods We used a Bayesian geostatistical model including 15 spatial covariates and data from 385 household surveys across 94 LMICs to estimate annual proportions of children younger than 5 years of age with diarrhoea who received ORS or RHF (or both) on continuous continent-wide surfaces in 2000-17, and aggregated results to policy-relevant administrative units. Additionally, we analysed geographical inequality in coverage across administrative units and estimated the number of diarrhoeal deaths averted by increased coverage over the study period. Uncertainty in the mean coverage estimates was calculated by taking 250 draws from the posterior joint distribution of the model and creating uncertainty intervals (UIs) with the 2 center dot 5th and 97 center dot 5th percentiles of those 250 draws. Findings While ORS use among children with diarrhoea increased in some countries from 2000 to 2017, coverage remained below 50% in the majority (62 center dot 6%; 12 417 of 19 823) of second administrative-level units and an estimated 6 519 000 children (95% UI 5 254 000-7 733 000) with diarrhoea were not treated with any form of ORT in 2017. Increases in ORS use corresponded with declines in RHF in many locations, resulting in relatively constant overall ORT coverage from 2000 to 2017. Although ORS was uniformly distributed subnationally in some countries, within-country geographical inequalities persisted in others; 11 countries had at least a 50% difference in one of their units compared with the country mean. Increases in ORS use over time were correlated with declines in RHF use and in diarrhoeal mortality in many locations, and an estimated 52 230 diarrhoeal deaths (36 910-68 860) were averted by scaling up of ORS coverage between 2000 and 2017. Finally, we identified key subnational areas in Colombia, Nigeria, and Sudan as examples of where diarrhoeal mortality remains higher than average, while ORS coverage remains lower than average. Interpretation To our knowledge, this study is the first to produce and map subnational estimates of ORS, RHF, and ORT coverage and attributable child diarrhoeal deaths across LMICs from 2000 to 2017, allowing for tracking progress over time. Our novel results, combined with detailed subnational estimates of diarrhoeal morbidity and mortality, can support subnational needs assessments aimed at furthering policy makers' understanding of within-country disparities. Over 50 years after the discovery that led to this simple, cheap, and life-saving therapy, large gains in reducing mortality could still be made by reducing geographical inequalities in ORS coverage. Copyright (c) 2020 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license.Peer reviewe