PT  - JOURNAL ARTICLE
AU  - Shrestha, Himal
AU  - McCulloch, Karen
AU  - Chisholm, Rebecca H
AU  - Armoo, Samuel
AU  - Vierigh, Francis
AU  - Sirwani, Neha
AU  - Crawford, Katie E
AU  - Osei-Atweneboana, Mike
AU  - Grant, Warwick N
AU  - Hedtke, Shannon M
TI  - Synthesising environmental, epidemiological, and genetic data to assist decision making for onchocerciasis elimination
AID  - 10.1101/2023.02.14.23285937
DP  - 2023 Jan 01
TA  - medRxiv
PG  - 2023.02.14.23285937
4099  - http://medrxiv.org/content/early/2023/02/22/2023.02.14.23285937.short
4100  - http://medrxiv.org/content/early/2023/02/22/2023.02.14.23285937.full
AB  - Background Population genetics is crucial for understanding the transmission dynamics of diseases like onchocerciasis. Landscape genetics identifies the ecological features that impact genetic variation between sampling sites. Here, we have used a landscape genetics framework to understand the relationship between environmental features and gene flow of the filarial parasite Onchocerca volvulus and of its intermediate host and vector, blackflies in the genus Simulium. We analysed samples from the ecological transition region separating the savannah and forest ecological regions of Ghana, where the transmission of O. volvulus has persisted despite almost half a century of onchocerciasis control efforts.Methods We generated a baseline microfilarial prevalence map from the point estimates of pre-ivermectin microfilarial prevalence from 47 locations in the study area. We analysed mitochondrial data from 164 parasites and 93 blackflies collected from 15 communities and four breeding sites, respectively. We estimated population genetic diversity and identified correlations with environmental variables. Finally, we compared baseline prevalence maps to movement suitability maps that were based on significant environmental variables.Results We found that the resistance surfaces derived from elevation (r = 0.793, p = 0.005) and soil moisture (r = 0.507, p = 0.002) were significantly associated with genetic distance between parasite sampling locations. Similarly, for the vector populations, the resistance surfaces derived from soil moisture (r = 0.788, p = 0.0417) and precipitation (r = 0.835, p = 0.0417) were significant. The correlation between the baseline parasite prevalence map and the parasite resistance surface map was stronger than the correlation between baseline prevalence and the vector resistance surface map. The central parts of the transition region which were conducive for both the parasite and the vector gene flow were most strongly associated with high baseline onchocerciasis prevalence.Conclusions We present a framework for incorporating environmental, genetic, and prevalence data for identifying when ecological conditions are favourable for onchocerciasis transmission between communities. We identified areas with higher suitability for parasite and vector gene flow, which ultimately might help us gain deeper insights into defining transmission zones for onchocerciasis. Furthermore, this framework is translatable to other onchocerciasis endemic areas and to other vector-borne diseases.Competing Interest StatementThe authors have declared no competing interest.Funding StatementThis work was supported by funding from UNICEF/UNDP/World Bank/WHO Special Programme for Research and Training in Tropical Diseases (TDR) to SMH (P21-00481). HS was supported by Australian Government Research Training Program Scholarship, a La Trobe Graduate Research Scholarship and a La Trobe University Full-Fee Research Scholarship.Author DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesThe details of the IRB/oversight body that provided approval or exemption for the research described are given below:The study used source data for onchocerciasis prevalence that are openly available at the ESPEN data portal (https://espen.afro.who.int/tools-resources/download-data).I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).YesI have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesThe parasite sequence data are available at NCBI (Accession #: PRJNA560089), and the blackfly sequence data have been uploaded to EMBL-EBI (Accession #: PRJEB57094). The onchocerciasis prevalence data were obtained from the ESPEN data portal, and the sources for the environmental data are provided in the supplementary information. The scripts for the analysis pipeline are uploaded to the GitHub repository. https://github.com/himal2007/landscape_genetics_ghana https://espen.afro.who.int/tools-resources/download-data https://www.ebi.ac.uk/ https://www.ncbi.nlm.nih.gov/ BCIBayesian credible intervalBICBayesian information criteriaDAPCDiscriminant analysis of principal componentsDNADeoxyribonucleic AcidEBIEuropean Bioinformatics InstituteEMBLEuropean Molecular Biology LaboratoryENAEuropean Nucleotide ArchiveESPENExpanded Special Project for Elimination of Neglected Tropical DiseaseINLAIntegrated nested Laplace approximationsMDAiMass drug administration with ivermectinMLPEMaximum likelihood population effectsMMRRMixed matrix regression with randomisationNCBINational centre for biotechnology informationOCPOnchocerciasis Control ProgrammePCAPrincipal component analysisSDStandard deviationSEStandard errorSNPSingle nucleotide polymorphismSORSub-optimal response