Evolution of insecticide resistance and its mechanisms in Anopheles stephensi in the WHO Eastern Mediterranean Region

Background: While Iran is on the path to eliminating malaria, the disease with 4.9 million estimated cases and 9300 estimated deaths in 2018 remains a serious health problem in the World Health Organization (WHO) Eastern Mediterranean Region. Anopheles stephensi is the main malaria vector in Iran and its range extends from Iraq to western China. Recently, the vector invaded new territories in Sri Lanka and countries in the Horn of Africa. Insecticide resistance in An. stephensi is a potential issue in controlling the spread of this vector. Methods: Data were collated from national and international databases, including PubMed, Google Scholar, Scopus, ScienceDirect, SID, and IranMedex using appropriate search terms. Results: Indoor residual spaying (IRS) with DDT was piloted in Iran in 1945 and subsequently used in the malaria eradication programme. Resistance to DDT in An. stephensi was detected in Iran, Iraq, Pakistan, and Saudi Arabia in the late 1960s. Malathion was used for malaria control in Iran in 1967, then propoxur in 1978, followed by pirimiphosmethyl from 1992 to 1994. The pyrethroid insecticide lambda-cyhalothrin was used from 1994 to 2003 followed by deltamethrin IRS and long-lasting insecticidal nets (LLINs). Some of these insecticides with the same sequence were used in other malaria-endemic countries of the region. Pyrethroid resistance was detected in An. stephensi in Afghanistan in 2010, in 2011 in India and in 2012 in Iran. The newly invaded population of An. stephensi in Ethiopia was resistant to insecticides of all four major insecticide classes. Different mechanisms of insecticide resistance, including metabolic and insecticide target site insensitivity, have been developed in An. stephensi. Resistance to DDT was initially glutathione S-transferase based. Target site knockdown resistance was later selected by pyrethroids. Esterases and altered acetylcholinesterase are the underlying cause of organophosphate resistance and cytochrome p450s were involved in pyrethroid metabolic resistance. Conclusions: Anopheles stephensi is a major malaria vector in Iran and many countries in the region and beyond. The species is leading in terms of development of insecticide resistance as well as developing a variety of resistance mechanisms. Knowledge of the evolution of insecticide resistance and their underlying mechanisms, in particular, are important to Iran, considering the final steps the country is taking towards malaria elimination, but also to other countries in the region for their battle against malaria. This systematic review may also be of value to countries and territories newly invaded by this species, especially in the Horn of Africa, where the malaria situation is already dire

Evidence of metabolic mechanisms playing a role in multiple insecticides resistance in Anopheles stephensi populations from Afghanistan.

Background Malaria is endemic in most parts of Afghanistan and insecticide-based vector control measures are central in controlling the disease. Insecticide resistance in the main malaria vector Anopheles stephensi from Afghanistan is increasing and attempts should be made to determine the underlying resistance mechanisms for its adequate management. Methods The contents of cytochrome P450s, esterases, glutathione S-transferases (GSTs) and acetylcholine esterase (AChE) activities were measured in the Kunar and Nangarhar populations of An. stephensi from Afghanistan and the results were compared with those of the susceptible Beech strain using the World Health Organization approved biochemical assay methods for adult mosquitoes. Results The cytochrome P450s enzyme ratios were 2.23- and 2.54-fold in the Kunar and Nangarhar populations compared with the susceptible Beech strain. The enzyme ratios for esterases with alpha-naphthyl acetate were 1.45 and 2.11 and with beta-naphthyl acetate were 1.62 and 1.85 in the Kunar and Nangarhar populations respectively compared with the susceptible Beech strain. Esterase ratios with para-nitrophenyl acetate (pNPA) were 1.61 and 1.75 in the Kunar and Nangarhar populations compared with the susceptible Beech strain. The GSTs enzyme ratios were 1.33 and 1.8 in the Kunar and Nangarhar populations compared with the susceptible Beech strain. The inhibition of AChE was 70.9 in the susceptible Beech strain, and 56.7 and 51.5 in the Kunar and Nangarhar populations. The differences between all values of the enzymes activities/contents and AChE inhibition rates in the Kunar and Nangarhar populations were statistically significant when compared with those of the susceptible Beech strain. Conclusions Based on the results, the reported resistance to pyrethroid and organophosphate insecticides, and tolerance to bendiocarb in the Kunar and Nangarhar populations of An. stephensi from Afghanistan are likely to be caused by a range of metabolic mechanisms, including esterases, P450s and GSTs combined with target site insensitivity in AChE

Status of insecticide resistance and its biochemical and molecular mechanisms in Anopheles stephensi (Diptera: Culicidae) from Afghanistan

Background Insecticide resistance of Anopheles stephensi, the main malaria vector in eastern Afghanistan, has been reported previously. This study describes the biochemical and molecular mechanisms of resistance to facilitate effective vector control and insecticide resistance management. Methods Mosquito larvae were collected from the provinces of Kunar, Laghman and Nangarhar from 2014 to 2017. The susceptibility of the reared 3–4 days old adults was tested with deltamethrin 0.05%, bendiocarb 0.1%, malathion 5%, permethrin 0.75% and DDT 4%. Cytochrome P450 content and general esterase, glutathione S-transferase (GST) and acetylcholinesterase (AChE) activities were measured in the three field populations and the results were compared with those of the laboratory susceptible An. stephensi Beech strain. Two separate allele-specific PCR assays were used to identify L1014, L1014F and L1014S mutations in the voltage gated sodium channel gene of An. stephensi. Probit analysis, ANOVA and Hardy–Weinberg equilibrium were used to analyse bioassay, biochemical assay and gene frequency data respectively. Results The population of An. stephensi from Kunar was susceptible to bendiocarb, apart from this, all populations were resistant to all the other insecticides tested. The differences between all values for cytochrome P450s, general esterases, GSTs and AChE inhibition rates in the Kunar, Laghman and Nangarhar populations were statistically significant when compared to the Beech strain, excluding GST activities between Kunar and Beech due to the high standard deviation in Kunar. The three different sodium channel alleles [L1014 (wild type), L1014F (kdr west) and L1014S (kdr east)] were all segregated in the Afghan populations. The frequencies of kdr east mutation were 22.9%, 32.7% and 35% in Kunar, Laghman and Nangarhar populations respectively. Kdr west was at the lowest frequency of 4.44%. Conclusions Resistance to different groups of insecticides in the field populations of An. stephensi from Kunar, Laghman and Nangarhar Provinces of Afghanistan is caused by a range of metabolic and site insensitivity mechanisms, including esterases, cytochrome P450s and GSTs combined with AChE and sodium channel target site insensitivity. The intensity and frequency of these mechanisms are increasing in these populations, calling for urgent reorientation of vector control programmes and implementation of insecticide resistance management strategies.d

Global water quality changes posing threat of increasing infectious diseases, a case study on malaria vector Anopheles stephensi coping with the water pollutants using age-stage, two-sex life table method

Background: Water pollution due to uncontrolled release of chemical pollutants is an important global problem. Its effect on medically important insects, especially mosquitoes, is a critical issue in the epidemiology of mosquito-borne diseases. Methods: In order to understand the effect of water pollutants on the demography of Anopheles stephensi, colonies were reared in clean, moderately and highly polluted water for three consecutive generations at 27 °C, 75% RH, and a photoperiod of 12:12 h (L:D). The demographic data of the 4th generation of An. stephensi were collected and analysed using the age-stage, two-sex life table. Results: The intrinsic rate of increase (r), finite rate of increase (λ), mean fecundity (F) and net reproductive rate (R0) of An. stephensi in clean water were 0.2568 d−1, 1.2927 d−1, 251.72 eggs, and 109.08 offspring, respectively. These values were significantly higher than those obtained in moderately polluted water (r = 0.2302 d−1, λ = 1.2589 d−1, 196.04 eggs, and R0 = 65.35 offspring) and highly polluted water (r = 0.2282 d−1, λ = 1.2564 d−1, 182.45 eggs, and R0 = 62.03 offspring). Female adult longevity in moderately polluted (9.38 days) and highly polluted water (9.88 days) were significantly shorter than those reared in clean water (12.43 days), while no significant difference in the male adult longevity was observed among treatments. Conclusions: The results of this study showed that An. stephensi can partially adapt to water pollution and this may be sufficient to extend the range of mosquito-borne diseases

Kdr genotyping and the first report of V410L and V1016I kdr mutations in voltage-gated sodium channel gene in Aedes aegypti (Diptera: Culicidae) from Iran

Background: Aedes aegypti is the main vector of arboviral diseases worldwide. The species invaded and became established in southern Iran in 2020. Insecticide-based interventions are primarily used for its control. With insecticide resistance widespread, knowledge of resistance mechanisms is vital for informed deployment of insecticidal interventions, but information from Iranian Ae. aegypti is lacking. Methods: Fifty-six Ae. aegypti specimens were collected from the port city of Bandar Lengeh in Hormozgan Province in the South of Iran in 2020 and screened for kdr mutations. The most common kdr mutations in Latin America and Asia (V410L, S989P, V1016G/I and F1534C), especially when present in combinations, are highly predictive of DDT and pyrethroid resistance were detected. Phylogenetic analyses based on the diversity of S989P and V1016G/I mutations were undertaken to assess the phylogeography of these kdr mutations. Results: Genotyping all four kdr positions of V410L, S989P, V1016G/I and F1534C revealed that only 16 out of the 56 (28.57%) specimens were homozygous wild type for all kdr mutation sites. Six haplotypes including VSVF (0.537), VSVC (0.107), LSVF (0.016), LSIF (0.071), VPGC (0.257) and LPGC (0.011) were detected in this study. For the first time, 11 specimens harbouring the V410L mutation, and 8 samples with V1016I mutation were found. V410L and V1016I were coincided in 8 specimens. Also, six specimens contained 1016G/I double mutation which was not reported before. Conclusions: The relatively high frequency of these kdr mutations in Iranian Ae. aegypti indicates a population exhibiting substantial resistance to pyrethroid insecticides, which are used widely in control operations and household formulations. The detection of the 410L/1016I kdr mutant haplotype in Iranian Ae. aegypti suggests possible convergence of invasive populations from West Africa or Latin America. However, as Iran has very limited maritime/air connections with those African countries, a Latin American origin for the invasive Ae. aegypti in Iran is more plausible

Measuring progress and projecting attainment on the basis of past trends of the health-related Sustainable Development Goals in 188 countries: an analysis from the Global Burden of Disease Study 2016

The UN’s Sustainable Development Goals (SDGs) are grounded in the global ambition of “leaving no one behind”. Understanding today’s gains and gaps for the health-related SDGs is essential for decision makers as they aim to improve the health of populations. As part of the Global Burden of Diseases, Injuries, and Risk Factors Study 2016 (GBD 2016), we measured 37 of the 50 health-related SDG indicators over the period 1990–2016 for 188 countries, and then on the basis of these past trends, we projected indicators to 2030

Global, regional, and national incidence, prevalence, and years lived with disability for 328 diseases and injuries for 195 countries, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016

As mortality rates decline, life expectancy increases, and populations age, non-fatal outcomes of diseases and injuries are becoming a larger component of the global burden of disease. The Global Burden of Diseases, Injuries, and Risk Factors Study 2016 (GBD 2016) provides a comprehensive assessment of prevalence, incidence, and years lived with disability (YLDs) for 328 causes in 195 countries and territories from 1990 to 2016

Global, regional, and national under-5 mortality, adult mortality, age-specific mortality, and life expectancy, 1970–2016: a systematic analysis for the Global Burden of Disease Study 2016

BACKGROUND: Detailed assessments of mortality patterns, particularly age-specific mortality, represent a crucial input that enables health systems to target interventions to specific populations. Understanding how all-cause mortality has changed with respect to development status can identify exemplars for best practice. To accomplish this, the Global Burden of Diseases, Injuries, and Risk Factors Study 2016 (GBD 2016) estimated age-specific and sex-specific all-cause mortality between 1970 and 2016 for 195 countries and territories and at the subnational level for the five countries with a population greater than 200 million in 2016. METHODS: We have evaluated how well civil registration systems captured deaths using a set of demographic methods called death distribution methods for adults and from consideration of survey and census data for children younger than 5 years. We generated an overall assessment of completeness of registration of deaths by dividing registered deaths in each location-year by our estimate of all-age deaths generated from our overall estimation process. For 163 locations, including subnational units in countries with a population greater than 200 million with complete vital registration (VR) systems, our estimates were largely driven by the observed data, with corrections for small fluctuations in numbers and estimation for recent years where there were lags in data reporting (lags were variable by location, generally between 1 year and 6 years). For other locations, we took advantage of different data sources available to measure under-5 mortality rates (U5MR) using complete birth histories, summary birth histories, and incomplete VR with adjustments; we measured adult mortality rate (the probability of death in individuals aged 15-60 years) using adjusted incomplete VR, sibling histories, and household death recall. We used the U5MR and adult mortality rate, together with crude death rate due to HIV in the GBD model life table system, to estimate age-specific and sex-specific death rates for each location-year. Using various international databases, we identified fatal discontinuities, which we defined as increases in the death rate of more than one death per million, resulting from conflict and terrorism, natural disasters, major transport or technological accidents, and a subset of epidemic infectious diseases; these were added to estimates in the relevant years. In 47 countries with an identified peak adult prevalence for HIV/AIDS of more than 0·5% and where VR systems were less than 65% complete, we informed our estimates of age-sex-specific mortality using the Estimation and Projection Package (EPP)-Spectrum model fitted to national HIV/AIDS prevalence surveys and antenatal clinic serosurveillance systems. We estimated stillbirths, early neonatal, late neonatal, and childhood mortality using both survey and VR data in spatiotemporal Gaussian process regression models. We estimated abridged life tables for all location-years using age-specific death rates. We grouped locations into development quintiles based on the Socio-demographic Index (SDI) and analysed mortality trends by quintile. Using spline regression, we estimated the expected mortality rate for each age-sex group as a function of SDI. We identified countries with higher life expectancy than expected by comparing observed life expectancy to anticipated life expectancy on the basis of development status alone. FINDINGS: Completeness in the registration of deaths increased from 28% in 1970 to a peak of 45% in 2013; completeness was lower after 2013 because of lags in reporting. Total deaths in children younger than 5 years decreased from 1970 to 2016, and slower decreases occurred at ages 5-24 years. By contrast, numbers of adult deaths increased in each 5-year age bracket above the age of 25 years. The distribution of annualised rates of change in age-specific mortality rate differed over the period 2000 to 2016 compared with earlier decades: increasing annualised rates of change were less frequent, although rising annualised rates of change still occurred in some locations, particularly for adolescent and younger adult age groups. Rates of stillbirths and under-5 mortality both decreased globally from 1970. Evidence for global convergence of death rates was mixed; although the absolute difference between age-standardised death rates narrowed between countries at the lowest and highest levels of SDI, the ratio of these death rates-a measure of relative inequality-increased slightly. There was a strong shift between 1970 and 2016 toward higher life expectancy, most noticeably at higher levels of SDI. Among countries with populations greater than 1 million in 2016, life expectancy at birth was highest for women in Japan, at 86·9 years (95% UI 86·7-87·2), and for men in Singapore, at 81·3 years (78·8-83·7) in 2016. Male life expectancy was generally lower than female life expectancy between 1970 and 2016, an

Population and fertility by age and sex for 195 countries and territories, 1950–2017: a systematic analysis for the Global Burden of Disease Study 2017

Background: Population estimates underpin demographic and epidemiological research and are used to track progress on numerous international indicators of health and development. To date, internationally available estimates of population and fertility, although useful, have not been produced with transparent and replicable methods and do not use standardised estimates of mortality. We present single-calendar year and single-year of age estimates of fertility and population by sex with standardised and replicable methods. Methods: We estimated population in 195 locations by single year of age and single calendar year from 1950 to 2017 with standardised and replicable methods. We based the estimates on the demographic balancing equation, with inputs of fertility, mortality, population, and migration data. Fertility data came from 7817 location-years of vital registration data, 429 surveys reporting complete birth histories, and 977 surveys and censuses reporting summary birth histories. We estimated age-specific fertility rates (ASFRs; the annual number of livebirths to women of a specified age group per 1000 women in that age group) by use of spatiotemporal Gaussian process regression and used the ASFRs to estimate total fertility rates (TFRs; the average number of children a woman would bear if she survived through the end of the reproductive age span [age 10–54 years] and experienced at each age a particular set of ASFRs observed in the year of interest). Because of sparse data, fertility at ages 10–14 years and 50–54 years was estimated from data on fertility in women aged 15–19 years and 45–49 years, through use of linear regression. Age-specific mortality data came from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2017 estimates. Data on population came from 1257 censuses and 761 population registry location-years and were adjusted for underenumeration and age misreporting with standard demographic methods. Migration was estimated with the GBD Bayesian demographic balancing model, after incorporating information about refugee migration into the model prior. Final population estimates used the cohort-component method of population projection, with inputs of fertility, mortality, and migration data. Population uncertainty was estimated by use of out-of-sample predictive validity testing. With these data, we estimated the trends in population by age and sex and in fertility by age between 1950 and 2017 in 195 countries and territories. Findings: From 1950 to 2017, TFRs decreased by 49\ub74% (95% uncertainty interval [UI] 46\ub74–52\ub70). The TFR decreased from 4\ub77 livebirths (4\ub75–4\ub79) to 2\ub74 livebirths (2\ub72–2\ub75), and the ASFR of mothers aged 10–19 years decreased from 37 livebirths (34–40) to 22 livebirths (19–24) per 1000 women. Despite reductions in the TFR, the global population has been increasing by an average of 83\ub78 million people per year since 1985. The global population increased by 197\ub72% (193\ub73–200\ub78) since 1950, from 2\ub76 billion (2\ub75–2\ub76) to 7\ub76 billion (7\ub74–7\ub79) people in 2017; much of this increase was in the proportion of the global population in south Asia and sub-Saharan Africa. The global annual rate of population growth increased between 1950 and 1964, when it peaked at 2\ub70%; this rate then remained nearly constant until 1970 and then decreased to 1\ub71% in 2017. Population growth rates in the southeast Asia, east Asia, and Oceania GBD super-region decreased from 2\ub75% in 1963 to 0\ub77% in 2017, whereas in sub-Saharan Africa, population growth rates were almost at the highest reported levels ever in 2017, when they were at 2\ub77%. The global average age increased from 26\ub76 years in 1950 to 32\ub71 years in 2017, and the proportion of the population that is of working age (age 15–64 years) increased from 59\ub79% to 65\ub73%. At the national level, the TFR decreased in all countries and territories between 1950 and 2017; in 2017, TFRs ranged from a low of 1\ub70 livebirths (95% UI 0\ub79–1\ub72) in Cyprus to a high of 7\ub71 livebirths (6\ub78–7\ub74) in Niger. The TFR under age 25 years (TFU25; number of livebirths expected by age 25 years for a hypothetical woman who survived the age group and was exposed to current ASFRs) in 2017 ranged from 0\ub708 livebirths (0\ub707–0\ub709) in South Korea to 2\ub74 livebirths (2\ub72–2\ub76) in Niger, and the TFR over age 30 years (TFO30; number of livebirths expected for a hypothetical woman ageing from 30 to 54 years who survived the age group and was exposed to current ASFRs) ranged from a low of 0\ub73 livebirths (0\ub73–0\ub74) in Puerto Rico to a high of 3\ub71 livebirths (3\ub70–3\ub72) in Niger. TFO30 was higher than TFU25 in 145 countries and territories in 2017. 33 countries had a negative population growth rate from 2010 to 2017, most of which were located in central, eastern, and western Europe, whereas population growth rates of more than 2\ub70% were seen in 33 of 46 countries in sub-Saharan Africa. In 2017, less than 65% of the national population was of working age in 12 of 34 high-income countries, and less than 50% of the national population was of working age in Mali, Chad, and Niger. Interpretation: Population trends create demographic dividends and headwinds (ie, economic benefits and detriments) that affect national economies and determine national planning needs. Although TFRs are decreasing, the global population continues to grow as mortality declines, with diverse patterns at the national level and across age groups. To our knowledge, this is the first study to provide transparent and replicable estimates of population and fertility, which can be used to inform decision making and to monitor progress. Funding: Bill & Melinda Gates Foundation

Global, regional, and national disability-adjusted life-years (DALYs) for 333 diseases and injuries and healthy life expectancy (HALE) for 195 countries and territories, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016

BACKGROUND: Measurement of changes in health across locations is useful to compare and contrast changing epidemiological patterns against health system performance and identify specific needs for resource allocation in research, policy development, and programme decision making. Using the Global Burden of Diseases, Injuries, and Risk Factors Study 2016, we drew from two widely used summary measures to monitor such changes in population health: disability-adjusted life-years (DALYs) and healthy life expectancy (HALE). We used these measures to track trends and benchmark progress compared with expected trends on the basis of the Socio-demographic Index (SDI). METHODS: We used results from the Global Burden of Diseases, Injuries, and Risk Factors Study 2016 for all-cause mortality, cause-specific mortality, and non-fatal disease burden to derive HALE and DALYs by sex for 195 countries and territories from 1990 to 2016. We calculated DALYs by summing years of life lost and years of life lived with disability for each location, age group, sex, and year. We estimated HALE using age-specific death rates and years of life lived with disability per capita. We explored how DALYs and HALE differed from expected trends when compared with the SDI: the geometric mean of income per person, educational attainment in the population older than age 15 years, and total fertility rate. FINDINGS: The highest globally observed HALE at birth for both women and men was in Singapore, at 75·2 years (95% uncertainty interval 71·9-78·6) for females and 72·0 years (68·8-75·1) for males. The lowest for females was in the Central African Republic (45·6 years [42·0-49·5]) and for males was in Lesotho (41·5 years [39·0-44·0]). From 1990 to 2016, global HALE increased by an average of 6·24 years (5·97-6·48) for both sexes combined. Global HALE increased by 6·04 years (5·74-6·27) for males and 6·49 years (6·08-6·77) for females, whereas HALE at age 65 years increased by 1·78 years (1·61-1·93) for males and 1·96 years (1·69-2·13) for females. Total global DALYs remained largely unchanged from 1990 to 2016 (-2·3% [-5·9 to 0·9]), with decreases in communicable, maternal, neonatal, and nutritional (CMNN) disease DALYs offset by increased DALYs due to non-communicable diseases (NCDs). The exemplars, calculated as the five lowest ratios of observed to expected age-standardised DALY rates in 2016, were Nicaragua, Costa Rica, the Maldives, Peru, and Israel. The leading three causes of DALYs globally were ischaemic heart disease, cerebrovascular disease, and lower respiratory infections, comprising 16·1% of all DALYs. Total DALYs and age-standardised DALY rates due to most CMNN causes decreased from 1990 to 2016. Conversely, the total DALY burden rose for most NCDs; however, age-standardised DALY rates due to NCDs declined globally. INTERPRETATION: At a global level, DALYs and HALE continue to show improvements. At the same time, we observe that many populations are facing growing functional health loss. Rising SDI was associated with increases in cumulative years of life lived with disability and decreases in CMNN DALYs offset by increased NCD DALYs. Relative compression of morbidity highlights the importance of continued health interventions, which has changed in most locations in pace with the gross domestic product per person, education, and family planning. The analysis of DALYs and HALE and their relationship to SDI represents a robust framework with which to benchmark location-specific health performance. Country-specific drivers of disease burden, particularly for causes with higher-than-expected DALYs, should inform health policies, health system improvement initiatives, targeted prevention efforts, and development assistance for health, including financial and research investments for all countries, regardless of their level of sociodemographic development. The presence of countries that substantially outperform others suggests the need for increased scrutiny for proven examples of best practices, which can help to extend gains, whereas the presence of underperforming countries suggests the need for devotion of extra attention to health systems that need more robust support. FUNDING: Bill & Melinda Gates Foundation