Skip to main content
Advertisement
  • Loading metrics

Potential mitigating role of ivermectin on the spread of Chlamydia trachomatis by Musca sorbens

  • Richard Selby ,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Supervision, Visualization, Writing – original draft

    rselby@sightsavers.org

    Affiliation Sightsavers international, 35 Perrymount Road, Haywards Heath, West Sussex, United Kingdom

  • Anita Jeyam,

    Roles Data curation, Formal analysis, Methodology, Validation, Writing – review & editing

    Affiliation Sightsavers international, 35 Perrymount Road, Haywards Heath, West Sussex, United Kingdom

  • Andrew Tate,

    Roles Data curation, Formal analysis, Investigation, Methodology, Resources, Writing – review & editing

    Affiliation Sightsavers international, 35 Perrymount Road, Haywards Heath, West Sussex, United Kingdom

  • Fikreab Kebede,

    Roles Data curation, Visualization, Writing – review & editing

    Affiliation Ministry of Health, Addis Ababa, Ethiopia

  • Philip Downs

    Roles Conceptualization, Formal analysis, Methodology, Writing – review & editing

    Affiliation Sightsavers international, 35 Perrymount Road, Haywards Heath, West Sussex, United Kingdom

Abstract

Trachoma is the world’s most frequent cause of blindness from an infectious agent. The disease caused by infection is associated with lack of access to sanitation and low hygiene standards. Trachoma is controlled through the Surgery, Antibiotics, Facial cleanliness, and Environmental improvement (SAFE) strategy, which delivers azithromycin (AZM) mass drug administration (MDA) in endemic areas. The putative vector Musca sorbens principally reproduce in human faecal matter left in the environment due to open defecation. Ivermectin (IVM) is on the WHO’s essential medicines list and is administered as preventative chemotherapy against two neglected tropical diseases (NTDs)—onchocerciasis, as an annual or bi-annual treatment, and lymphatic filariasis, as an annual treatment in combination with albendazole. Ivermectin has a known inhibitive effect on insects that reproduce in dung. To assess if IVM could be a viable vector control tool against M. sorbens, this study evaluates existing data from trachoma, onchocerciasis and lymphatic filariasis mass drug administration (MDA) operations in Ethiopia. Persistent and recrudescent trachoma in evaluation units (EUs) were examined for whether AZM MDA in EUs was accompanied by IVM MDA, and whether co-administration was associated with greater likelihood of trachoma control. Results show an association suggesting that EUs that received both IVM and AZM MDA benefit from improved control of trachoma in persistent or recrudescent areas, when compared to EUs that received AZM MDA. This initial investigation supports the potential for ivermectin’s use to support SAFE. Findings warrant further work to validate ivermectin’s impact on M. sorbens reproduction through controlled lab and field-based studies.

Author summary

Trachoma is the world’s most frequent cause of blindness from an infectious agent. Trachoma is controlled through the Surgery, Antibiotics, Facial cleanliness, and Environmental improvement (SAFE) strategy, which delivers azithromycin (AZM) mass drug administration (MDA) in endemic areas. The putative vector Musca sorbens principally reproduce in human faecal matter left in the environment due to open defecation. Ivermectin (IVM) is on the WHO’s essential medicines list, administered against two neglected tropical diseases (NTDs)–onchocerciasis and lymphatic filariasis. Ivermectin has an inhibitive effect on insects reproducing in dung.

To assess if IVM could be a viable vector control tool against M. sorbens, this study evaluates existing data from trachoma, onchocerciasis and lymphatic filariasis mass drug administration (MDA) operations in Ethiopia. Persistent and recrudescent trachoma in evaluation units (EUs) were examined for whether AZM MDA in EUs was accompanied by IVM MDA, and whether co-administration could have controlled the putative vector, unintentionally enhancing trachoma control.

Results suggest that EUs receiving both IVM and AZM MDA had improved control of trachoma, when compared to EUs that have only received AZM MDA.

This initial investigation supports conducting further work into ivermectin’s impact on M. sorbens reproduction through controlled lab and field-based studies.

Introduction

Trachoma is the world’s leading cause of blindness resulting from an infectious agent [1]. If left untreated trachoma leads to damage to the infected person’s eyes. This damage culminates in blindness, which is avoidable with intervention [2]. Trachoma is strongly associated with poverty along with a lack of access to sanitation [3,4]. The disease is prevalent throughout the tropics, with areas of Africa, Asia, the Middle East, Western Pacific and Latin America affected [5]. The causative bacterium Chlamydia trachomatis can spread between people by direct contact, fomites or by face seeking flies which mechanically transmit bacteria [6].

Of the face seeking flies, Musca sorbens is the putative principal vector responsible for C. trachomatis mechanical transmission [6]. M. sorbens are dung breeding flies that most commonly reproduce in human faeces that have been defecated in the open environment [7]. M. sorbens feed on sweat, mucosal secretions (particularly ocular and nasal), lesions, wounds and ulcers of the skin [8]. Female M. sorbens obtain the protein needed for egg production directly from human facial secretions [9]. These factors combine to form a prime mechanism to transmit eye infections, including trachoma. Studies conducted in Gambia demonstrated that reducing fly populations through insecticide spraying reduces trachoma prevalence by 56% [6].

Blindness from trachoma can be avoided through surgical intervention on infected individuals to correct the blinding stage (trachomatous trichiasis) [2]. Control efforts encompass mass drug administration (MDA) with azithromycin (AZM) [10] coupled with improving water, sanitation and hygiene (WASH) [11,12]. The combination of surgery, MDA and WASH contribute toward the SAFE strategy (Surgery, Antibiotics, Facial cleanliness, and Environmental improvement) to eliminate trachoma as a public health problem, and the transmission of C. trachomatis [13].

Ivermectin’s impact on dung breeding flies (Musca nevilli) has been investigated through treating cattle with sub-cutaneous injections of ivermectin (IVM) and letting flies oviposit on their dung [14]. Results showed reduced fly emergence from dung of IVM treated animals, lasting up to seven weeks post treatment. There was reduced fecundity amongst flies that successfully developed in dung of IVM treated animals [14]. Reproductive impact of IVM has been recorded for flies along with dung breeding beetles [15].

In addition, research into the acceptability of IVM MDA in Ghana, highlighted that communities had witnessed and appreciated IVM’s impact on dung breeding flies. Qualitative data collection revealed that communities had noted “Maggots in the public toilet died and that further convinced us of the efficacy of the drug” following IVM MDA. Existing work documents that IVM can be co-administered with other drugs including AZM in the context of NTD control [16,17] including in Ethiopia [18].

This led to initial questioning of whether IVM could provide a mechanism for M. sorbens control by effectively toxifying environmental human faeces? Would any IVM mediated impact upon M. sorbens populations affect trachoma, trachomatous inflammation–follicular (TF) persistence or recrudescence?

To investigate any potential impact of IVM on trachoma and TF, we conducted the following study focusing on Ethiopia, which has the highest burden of active trachoma and greatest recorded incidence of persistent and recrudescent TF globally [1,19]. Along with endemic areas for both onchocerciasis (OV) and lymphatic filariasis (LF), where MDA with IVM is conducted annually treating all individuals over 5 years of age [20,21]. As of 2020, sixteen percent of the Ethiopia population practice open defecation [22].

This document describes analysis conducted to investigate any previously unexplored association between IVM MDA and trachoma prevalence in Ethiopia. Authors highlight this study has been exploratory only, achieved by combining different data sources post-hoc to examine associations. This study aims to build a base of knowledge to guide further investigation.

Methods

For this secondary analysis, data were extracted for IVM and AZM administration throughout Ethiopia. Data concerning onchocerciasis and LF were obtained from the Expanded Special Project for Elimination of Neglected Tropical Diseases (ESPEN) platform between 2014 and 2020. Data for trachoma were also obtained through ESPEN, and using the GeoConnect administrative coding system. Trachoma data were evaluated between 2007 and 2020. Date ranges were selected to ensure completion and accuracy of data. From these datasets we identified locations based on two classifications:

1) Where onchocerciasis or lymphatic filariasis were co-endemic with trachoma (and received multiple rounds of IVM and AZM).

2) Where only trachoma is known, and AZM has been administered without IVM.

Trachoma data were summarized to one record per Evaluation Unit (EU). For each EU, information was gathered from the baseline survey up to the most recent impact (or surveillance) survey that categorised the EUs as either persistent or recrudescent.

To conduct this analysis we have used the following working definitions, as approved by WHO [23]:

A EU with persistent trachoma is classified as an EU with at least two impact surveys at which trachomatous inflammation-follicular (TF) prevalence in one to nine year olds (TF1–9) is ≥ 5%, without ever having had a TF1–9 < 5%.

A EU with recrudescent trachoma has at least one surveillance survey at which TF1–9 returns to ≥ 5%, and the current TF1–9 prevalence is ≥ 5%.

To evaluate trachoma, the data set was examined using a predefined sequential criterion, see Fig 1, using one record per EU defined by a unique GeoConnect identifier (ID). These records contain information on survey types, years and category of TF1–9 prevalence during each survey.

thumbnail
Fig 1. Exclusion and inclusion process of Ethiopian EUs for evaluating TF persistence and recrudescence against IVM and AZM MDA.

Note, 22 of the retained EUs matched criteria used for both recrudescence and persistence analysis.

https://doi.org/10.1371/journal.pntd.0011662.g001

Using the outlined definitions to explore persistence and recrudescence of trachoma [23]:

Recrudescence was examined among all EUs that had at least one surveillance survey.

Persistence was examined among the EUs that had 2 or more impact surveys.

Data for trachoma was merged with available onchocerciasis and LF data (recorded by ESPEN IDs), using a key Masterfile, matching ESPEN IDs and GeoConnect IDs.

Within the onchocerciasis and LF data, EUs were examined for if they had received any IVM MDA (either on its own or in association with another drug).

There are known variables that can impact M. sorbens reproduction and survival [24]. To examine the influence of these variables on this investigation available geographic data have been collated for elevation [2527], precipitation [28], landcover [29], population density [30], sanitation and water access [31]. These data were overlaid with the plots of EUs within our datasets for examination. Mean values for precipitation, elevation and population density per EU were calculated using zonal statistics in ArcGIS. For landcover, maximum coverage by land cover type (1km spatial resolution) per EU was utilised. Estimates for the proportion of people with access to improved sanitation and water were obtained from the ESPEN platform.

Statistical analyses were conducted at the EU level, using R v4.2.2. [32], using a significance level of 5%, no imputation of missing data was performed. Chi-square and fisher exact tests were used to examine the association between IVM and persistence/recrudescence. Within-sample quartiles were used to describe precipitation, elevation. Univariate logistic models were used to explore associations between potential individual confounders and persistence of trachoma. A multivariable model with simultaneous variable entry was used to assess the association between IVM and persistence adjusting for all potential confounders explored in this study. We supplemented this with a data-driven variable selection approach, using a backward selection procedure based on the AIC criterion to define an optimal parsimonious model.

Results

During data collation, EUs that could not be aligned from ESPEN into GeoConnect were discounted from our dataset prior to assessment commencing. This assessment commenced with 883 records for trachoma, one record mapped per EU, defined by a unique GeoConnect ID that contained information on the survey types, years of surveys and category of TF1–9 prevalence.

EUs were sequentially excluded due to lacking a baseline survey, retaining only EUs with baseline surveys. These EUs were then assessed to retain only those with a prevalence of TF1–9 ≥5% at baseline survey. Finally, to explore persistence and recrudescence, we only retained the 209 EUs with either 2 or more impact surveys or at least one impact survey and one surveillance survey.

We examined recrudescence among the 45 EUs that had at least one surveillance survey. Persistence was examined in the 186 EUs that were retained. Steps taken during the process are depicted in Fig 1.

As an initial exploration, simple cross-tabulations are presented between trachoma recrudescence and IVM MDA (Table 1) along with trachoma persistence and IVM distribution (Table 2).

thumbnail
Table 1. Cross tabulation of trachoma recrudescent EUs and mass IVM administration.

https://doi.org/10.1371/journal.pntd.0011662.t001

thumbnail
Table 2. Cross tabulation of trachoma persistent EUs and mass IVM administration.

https://doi.org/10.1371/journal.pntd.0011662.t002

Forty-four EUs that have received AZM MDA with at least one surveillance survey registered, were examined for trachoma recrudescence. Thirty-five of these EUs have not received IVM MDA, of which 13 did not fit the WHO criteria as being trachoma recrudescent and 21 meet classification as trachoma recrudescent. Among the ten EUs that have received both IVM MDA and AZM MDA, none of which met classification as TF recrudescence. There was a significant association between IVM MDA and recrudescence status (p<0.01).

In total 186 EUs had two or more impact surveys and were examined for trachoma persistence. In total 164 EUs had received only AZM MDA, of these 116 EUs were defined as trachoma persistent and 48 were not trachoma persistent. Twenty-two EUs had received both AZM and IVM. Of which 13 were classified as trachoma persistent and 9 were not trachoma persistent. There was a significant association between IVM MDA and persistence status (p<0.01).

Location of Ethiopian EUs are shown in Fig 2. Mapping EUs with and without IVM co-administration shows EUs with persistent and recrudescent trachoma have minimal overlap with IVM MDA. EUs that have received both AZM and IVM are less likely to have persistent or recrudescent trachoma among their population.

thumbnail
Fig 2. Locations of Ivermectin MDA related to Trachoma persistence and recrudescence status in Ethiopia.

Shapefile: https://espen.afro.who.int/tools-resources/cartography-database.

https://doi.org/10.1371/journal.pntd.0011662.g002

Confounders

Variables that could influence trachoma transmission rates are shown in maps that make up Fig 3, with one panel per variable. Maps are overlaid with the trachoma data and IVM MDA areas (for OV and LF) as presented in Fig 2.

thumbnail
Fig 3. Maps of influencing factors throughout EUs included in this study.

Clockwise from top left, elevation [27], precipitation, population density, access to improved sanitation, access to improved water and landcover. Shapefile: https://espen.afro.who.int/tools-resources/cartography-database.

https://doi.org/10.1371/journal.pntd.0011662.g003

Average precipitation, compiled by Fick and Hijmans [28] from August (peak rainfall in Ethiopia), indicates rainfall patterns are not vastly different between EUs receiving IVM MDA compared to trachoma associated EUs which met the criteria for inclusion in this study. Plotting landcover shows that EUs which are endemic for OV/LF where IVM MDA occurred were primarily classified as being forested with trees with some areas classified as crop, whereas the trachoma EUs are primarily classified as cropland and rangeland (grass and scrub). Population density throughout the OV/LF and trachoma EUs have a similar range, but it is notable that more IVM MDA EUs have lower categories of population density than the AZM MDA EUs. Comparison of access to improved water throughout the EUs included in this study does not reveal stark differences between trachoma EUs and OV/LF EUs. Access to improved sanitation is similar between the two classifications of EUs, with little improvement in the trachoma EUs that have also received IVM MDA.

EU characteristics by trachoma—persistence and recrudescence status

Table 3 details the EU confounder characteristics by trachoma status—recrudescence and persistence.

thumbnail
Table 3. EU characteristics by trachoma recrudescence and persistence status—descriptive statistics.

https://doi.org/10.1371/journal.pntd.0011662.t003

Among EUs with recrudescent trachoma, there was a higher proportion of EUs with lower baseline TF prevalence when compared to non-recrudescent trachoma EUs. The number of AZM administration years was lower. Mean population density was lower. The mean EU elevation was higher. The proportion of EUs with mean proportion of people with improved water access being in the highest quartile was lower, similarly for improved access to sanitation and the proportion of EUs with majority cropland was higher.

Compared to non-persistent EUs, among persistent EUs, the proportion with highest category of baseline TF prevalence was higher, as was the proportion with highest number of AZM years, the mean population density was lower, as was the mean elevation. The proportion of EUs with mean proportion of people with improved water access in the highest quartile was lower, similarly for improved sanitation, whereas the proportion of EUs in the highest precipitation quartile was lower. Finally, the proportion of EUs with majority cropland was lower while those with majority rangeland was higher.

Statistical modelling was only conducted on trachoma-persistence since the sample size for recrudescence was small (n = 45) and there was quasi-complete separation in the data with none of the recrudescent EUs having received IVM, the outcomes are given in Table 4.

Univariate models show that EUs that had received IVM were less likely to present persistent trachoma (OR = 0.29 [0.11, 0.71]). As were those with a lower number of AZM years (OR = 0.07 [0.01, 0.30] for 0–5 AZM years compared to more than 10 years). Those with highest mean population density (OR = 0.22 [0.08,0.53] for those with mean 318 or more people/km2 versus those with mean<152), those with higher elevation levels (OR = 0.35 [0.14, 0.85] for those with mean elevation 1740m/100m2 versus those with <1060m/100m2). Finally, the EUs with maximum land cover being rangeland were more likely to present persistent trachoma than those with maximum cropland (OR = 4.10 [1.86, 10.06]). There was no evidence of statistically significant associations between trachoma persistence and baseline TF prevalence levels, water or sanitation access nor precipitation levels.

In the full multivariable model adjusting for all potential confounders available, the association between IVM and persistence status is not statistically significant (p = 0.15). However, the optimal model resulting from the backward selection process includes IVM, number of AZM years, population density and precipitation levels. Results from this parsimonious model show that EUs with IVM, with a lower number of AZM years, higher population density and higher levels of precipitation were less likely to present persistent trachoma.

Discussion

This initial investigation suggests that MDA of IVM could be associated with lower odds of trachoma persistence and recrudescence. We hypothesise that IVM could be inadvertently delivering M. sorbens control by toxifying environmental human faeces, in turn this IVM mediated impact upon M. sorbens populations has supported SAFE strategy’s control of trachoma (TF). These findings warrant expansion of research to investigate the impact of IVM on M. sorbens, trachoma transmission and any environmental impact of IVM contaminated open defecation by humans. Such research would deliver an evaluation for repurposing IVM to support trachoma campaigns in difficult foci.

For the analysis presented here, we did not have access to dated historical data on IVM for OV and LF pre-2014. However, the earliest year of last trachoma survey considered for EUs was 2014 and IVM is not anticipated to have long-term residual impact capable of affecting M. sorbens reproduction. The plasma half-life of IVM in humans is approximately 18 hours following oral administration. Orally administered IVM and its metabolites are excreted in human faeces for 12 days post treatment [33].

The locations where IVM has been administered do not overlap with ongoing persistent or recrudescent trachoma. We have made efforts to investigate whether the spatial segregation is due to other influencing factors. Evaluating elevation, precipitation, land cover, population density, water and sanitation access throughout the EUs for trachoma, onchocerciasis and lymphatic filariasis transmission. Of note, authors made efforts to plot findings of this study alongside datasets for socio economic variables, but viable resolution for Ethiopia was not attainable for assessment.

It was expected that OV and LF endemic EUs would likely receive higher rainfall patterns than non-endemic EUs, as transmission of both depend on vectors which develop in aquatic environments. However, Fig 3 shows similar variations of rainfall throughout the investigated areas. This is relevant as rainfall accelerates the breakdown of faeces in the environment, removing the reproductive substrate for M. sorbens. which could have generated spatial differences in trachoma control. Plotting other factors that could influence trachoma transmissions indicates no stand-out difference in population density, access to improved sanitation or improved water access as presented. There is a viable difference in the landcover classifications throughout the OV and LF areas when compared to the trachoma endemic areas. But when considering EUs that had co-administration of AZM and IVM, the landcover includes both forest and cropland classifications present in these EUs that received IVM MDA, and trachoma reduced.

Variations in heavy rainfall could have driven differences in the abundance of environmental faecal matter and the speed of its decay. Thus, rainfall could impact M. sorbens reproduction between EUs with and without IVM MDA. While Fig 3 indicates this not to be the case. Authors acknowledge that Fig 3 uses the peak rainfall of August and other months may present variations that could influence M. sorbens reproduction. We have attempted to account for potential confounders that could affect the observed association between IVM MDA and enhanced trachoma control, using multivariable models. The model results are encouraging and highlight the need for future research to better understand the mechanisms at play behind the relationship between IVM, EU characteristics and trachoma so they are accurately reflected in the model building.

This research provides initial evidence that community wide IVM MDA may have a role in reducing trachoma transmission by delivering vector control. Future research into the impact on laboratory controlled and directly observed M. sorbens reproduction is vital. Laboratory research should be accompanied by case control investigations based in Ethiopia. Any field research should consider the introduction of IVM MDA in selected case EUs with both recrudescent or persistent trachoma with continuation of SAFE strategy including AZM MDA as currently managed.

Limitations

Authors highlight and acknowledge limitations within the study presented in this paper, which has been conducted as an assessment of the initial concept’s viability.

The analysed data is on an EU level only, repurposing data that had not been collected for the purpose of these analyses. We lack the granularity of individual-level data that would allow us to investigate socio-demographic and behavioural factors. However, population level effect on prevalence of TF is the principal outcome of interest while evaluating the potential of IVM for M. sorbens vector control. Authors highlight that some covariates used in our analysis have a degree of uncertainty, being themselves estimates derived from modelling, such as the WASH data for example.

As we have mentioned above and reiterate here, our analyses do not allow to infer a causal relationship between IVM MDA and lower odds of trachoma persistence or recrudescence and more research is warranted to understand the mechanisms driving this association.

Conclusion

Our analysis of existing data suggests an association between IVM MDA and improved control of trachoma in persistent or recrudescent areas. This finding warrants further research to ascertain the potential of IVM MDA implementation for trachoma control.

Following this initial assessment authors will expand the investigation by commencing insectary-based evaluation of M. sorbens reproduction in IVM contaminated human faeces. Concurrently, impact of IVM on adult flies will be investigated, to assess if adult fly mortality increases after imbibing facial secretions from an IVM treated person, as observed in heamatophagic ectoparasites [34,35].

Acknowledgments

Authors wish to express their gratitude to the Ethiopian Ministry of Health, International Trachoma Initiative (ITI), ESPEN and GeoConnect platforms for making the data used in this study openly available. We thank Andrew Balchin for support with graphical design of Fig 1 for use in the paper.

References

  1. 1. Mariotti SP, Pascolini D, Rose-Nussbaumer J. Trachoma: global magnitude of a preventable cause of blindness. British Journal of Ophthalmology. 2009;93(5):563–8. pmid:19098034
  2. 2. Kasi PM, Gilani AI, Ahmad K, Janjua NZ. Blinding trachoma: a disease of poverty. PLoS Med. 2004;1(2):e44. pmid:15578111
  3. 3. Garn JV, Boisson S, Willis R, Bakhtiari A, al-Khatib T, Amer K, et al. Sanitation and water supply coverage thresholds associated with active trachoma: Modeling cross-sectional data from 13 countries. PLOS Neglected Tropical Diseases. 2018;12(1):e0006110. pmid:29357365
  4. 4. Sullivan KM, Harding-Esch EM, Keil AP, Freeman MC, Batcho WE, Bio Issifou AA, et al. Exploring water, sanitation, and hygiene coverage targets for reaching and sustaining trachoma elimination: G-computation analysis. PLOS Neglected Tropical Diseases. 2023;17(2):e0011103. pmid:36780437
  5. 5. Polack S, Brooker S, Kuper H, Mariotti S, Mabey D, Foster A. Mapping the global distribution of trachoma. Bulletin of the World Health Organization. 2005;83(12):913–9. pmid:16462983
  6. 6. Emerson PM, Lindsay SW, Alexander N, Bah M, Dibba S-M, Faal HB, et al. Role of flies and provision of latrines in trachoma control: cluster-randomised controlled trial. The Lancet. 2004;363(9415):1093–8. pmid:15064026
  7. 7. Emerson PM, Bailey RL, Walraven GE, Lindsay SW. Human and other faeces as breeding media of the trachoma vector Musca sorbens. Med Vet Entomol. 2001;15(3):314–20. pmid:11583450
  8. 8. Robinson A, Bickford-Smith J, Abdurahman Shafi O, Abraham Aga M, Shuka G, Debela D, et al. Towards an odour-baited trap to control Musca sorbens, the putative vector of trachoma. Scientific Reports. 2021;11(1):14209. pmid:34244535
  9. 9. Emerson PM, Bailey RL, Mahdi OS, Walraven GE, Lindsay SW. Transmission ecology of the fly Musca sorbens, a putative vector of trachoma. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2000;94(1):28–32. pmid:10748893
  10. 10. West SK. Azithromycin for control of trachoma. Community Eye Health. 1999;12(32):55–6. pmid:17492006
  11. 11. Aragie S, Lietman TM, Keenan JD. Trachoma control using water, sanitation, and hygiene 2013; Authors’ reply. The Lancet Global Health. 2022;10(4):e480.
  12. 12. Stocks ME, Ogden S, Haddad D, Addiss DG, McGuire C, Freeman MC. Effect of Water, Sanitation, and Hygiene on the Prevention of Trachoma: A Systematic Review and Meta-Analysis. PLOS Medicine. 2014;11(2):e1001605. pmid:24586120
  13. 13. Jimenez V, Gelderblom HC, Mann Flueckiger R, Emerson PM, Haddad D. Mass drug administration for trachoma: how long is not long enough? PLoS Negl Trop Dis. 2015;9(3):e0003610. pmid:25799168
  14. 14. Krüger K, Scholtz CH. The effect of ivermectin on the development and reproduction of the dung-breeding fly Musca nevilli Kleynhans (Diptera, Muscidae). Agriculture, Ecosystems & Environment. 1995;53(1):13–8.
  15. 15. Strong L, Wall R, Woolford A, Djeddour D. The effect of faecally excreted ivermectin and fenbendazole on the insect colonisation of cattle dung following the oral administration of sustained-release boluses. Veterinary Parasitology. 1996;62(3):253–66. pmid:8686172
  16. 16. John LN, Gonzalez-Beiras C, Vall-Mayans M, Kolmau R, Houinei W, Wangi J, et al. Safety of mass drug coadministration with ivermectin, diethylcarbamazine, albendazole, and azithromycin for the integrated treatment of neglected tropical diseases: a cluster randomized community trial. The Lancet Regional Health–Western Pacific. 2022;18.
  17. 17. Romani L, Marks M, Sokana O, Nasi T, Kamoriki B, Wand H, et al. Feasibility and safety of mass drug coadministration with azithromycin and ivermectin for the control of neglected tropical diseases: a single-arm intervention trial. Lancet Glob Health. 2018;6(10):e1132–e8. pmid:30223985
  18. 18. McPherson S, Tafese G, Tafese T, Behaksra SW, Solomon H, Oljira B, et al. Safety of integrated mass drug administration of azithromycin, albendazole and ivermectin versus standard treatment regimens: a cluster-randomised trial in Ethiopia. EClinicalMedicine. 2023;59:101984. pmid:37152362
  19. 19. Kassim K, Kassim J, Aman R, Abduku M, Tegegne M, Sahiledengle B. Prevalence of active trachoma and associated risk factors among children of the pastoralist population in Madda Walabu rural district, Southeast Ethiopia: a community-based cross-sectional study. BMC Infectious Diseases. 2019;19(1):353. pmid:31035947
  20. 20. Gebrezgabiher G, Yewhalaw D, Hailu A, Mekonnen Z. Evaluation of mass treatment with ivermectin program reach and survey coverage for onchocerciasis elimination in selected endemic areas of Ethiopia. PLoS One. 2022;17(7):e0271518. pmid:35901026
  21. 21. Teshome A, Asfaw MA, Churko C, Yihune M, Chisha Y, Getachew B, et al. Coverage Validation Survey for Lymphatic Filariasis Treatment in Itang Special District of Gambella Regional State of Ethiopia: A Cross-Sectional Study. Infect Drug Resist. 2021;14:1537–43. pmid:33911881
  22. 22. Group WB. People practicing open defecation (% of population)—Ethiopia 2022 [Available from: https://data.worldbank.org/indicator/SH.STA.ODFC.ZS?locations=ET&view=chart.
  23. 23. World Health Organization. Informal consultation on end-game challenges for trachoma elimination, Task Force for Global Health, Decatur, United States of America, 7–9 December 2021. Geneva; 2022.
  24. 24. Ramesh A, Bristow J, Kovats S, Lindsay SW, Haslam D, Schmidt E, et al. The impact of climate on the abundance of Musca sorbens, the vector of trachoma. Parasit Vectors. 2016;9:48. pmid:26817815
  25. 25. Reuter HI, Nelson A, Jarvis A. An evaluation of void-filling interpolation methods for SRTM data. International Journal of Geographical Information Science. 2007;21(9):983–1008.
  26. 26. Alemayehu W, Melese M, Fredlander E, Worku A, Courtright P. Active trachoma in children in central Ethiopia: association with altitude. Trans R Soc Trop Med Hyg. 2005;99(11):840–3. pmid:16107273
  27. 27. Jarvis A, Reuter H, Nelson A, Guevara E. Hole-filled seamless SRTM data V4. Tech. rep. In: (CIAT) ICfTA, editor. Cali, Columbia2008.
  28. 28. Fick SE, Hijmans RJ. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. International Journal of Climatology. 2017;37(12):4302–15.
  29. 29. Karra K, Kontgis C, Statman-Weil Z, Mazzariello JC, Mathis M, Brumby SP, editors. Global land use / land cover with Sentinel 2 and deep learning. 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS; 2021 11–16 July 2021.
  30. 30. The spatial distribution of population density in 2020 based on country total adjusted to match corresponding UNDP estimate, Ethiopia [Internet]. 2020. Available from: https://hub.worldpop.org/geodata/summary?id=46094.
  31. 31. Deshpande A, Miller-Petrie MK, Lindstedt PA, Baumann MM, Johnson KB, Blacker BF, et al. Mapping geographical inequalities in access to drinking water and sanitation facilities in low-income and middle-income countries, 2000–17. The Lancet Global Health. 2020;8(9):e1162–e85. pmid:32827479
  32. 32. Team RC. R: A Language and Environment for Statistical Computing. In: Computing RFfS, editor. Vienna, Austria: https://www.R-project.org/; 2019.
  33. 33. Drugs.com. Ivermectin Tablets Auckland, New Zealand 2022 [Available from: https://www.drugs.com/pro/ivermectin-tablets.html.
  34. 34. McNeil D. Pill could join arsenal against bedbugs. The New York Times Dec. 2012;31:2012.
  35. 35. Martin RJ, Robertson AP, Choudhary S. Ivermectin: An Anthelmintic, an Insecticide, and Much More. Trends Parasitol. 2021;37(1):48–64. pmid:33189582