Two decades of molecular surveillance in Senegal reveal rapid changes in known drug resistance mutations over time.

BACKGROUND: Drug resistance in Plasmodium falciparum is a major threat to malaria control efforts. Pathogen genomic surveillance could be invaluable for monitoring current and emerging parasite drug resistance. METHODS: Data from two decades (2000-2020) of continuous molecular surveillance of P. falciparum parasites from Senegal were retrospectively examined to assess historical changes in malaria drug resistance mutations. Several known drug resistance markers and their surrounding haplotypes were profiled using a combination of single nucleotide polymorphism (SNP) molecular surveillance and whole genome sequence based population genomics. RESULTS: This dataset was used to track temporal changes in drug resistance markers whose timing correspond to historically significant events such as the withdrawal of chloroquine (CQ) and the introduction of sulfadoxine-pyrimethamine (SP) in 2003. Changes in the mutation frequency at Pfcrt K76T and Pfdhps A437G coinciding with the 2014 introduction of seasonal malaria chemoprevention (SMC) in Senegal were observed. In 2014, the frequency of Pfcrt K76T increased while the frequency of Pfdhps A437G declined. Haplotype-based analyses of Pfcrt K76T showed that this rapid increase was due to a recent selective sweep that started after 2014. DISCUSSION (CONCLUSION): The rapid increase in Pfcrt K76T is troubling and could be a sign of emerging amodiaquine (AQ) resistance in Senegal. Emerging AQ resistance may threaten the future clinical efficacy of artesunate-amodiaquine (ASAQ) and AQ-dependent SMC chemoprevention. These results highlight the potential of molecular surveillance for detecting rapid changes in parasite populations and stress the need to monitor the effectiveness of AQ as a partner drug for artemisinin-based combination therapy (ACT) and for chemoprevention.

Evaluating the performance of Plasmodium falciparum genetic metrics for inferring National Malaria Control Programme reported incidence in Senegal.

BACKGROUND: Genetic surveillance of the Plasmodium falciparum parasite shows great promise for helping National Malaria Control Programmes (NMCPs) assess parasite transmission. Genetic metrics such as the frequency of polygenomic (multiple strain) infections, genetic clones, and the complexity of infection (COI, number of strains per infection) are correlated with transmission intensity. However, despite these correlations, it is unclear whether genetic metrics alone are sufficient to estimate clinical incidence. METHODS: This study examined parasites from 3147 clinical infections sampled between the years 2012-2020 through passive case detection (PCD) across 16 clinic sites spread throughout Senegal. Samples were genotyped with a 24 single nucleotide polymorphism (SNP) molecular barcode that detects parasite strains, distinguishes polygenomic (multiple strain) from monogenomic (single strain) infections, and identifies clonal infections. To determine whether genetic signals can predict incidence, a series of Poisson generalized linear mixed-effects models were constructed to predict the incidence level at each clinical site from a set of genetic metrics designed to measure parasite clonality, superinfection, and co-transmission rates. RESULTS: Model-predicted incidence was compared with the reported standard incidence data determined by the NMCP for each clinic and found that parasite genetic metrics generally correlated with reported incidence, with departures from expected values at very low annual incidence (< 10/1000/annual [‰]). CONCLUSIONS: When transmission is greater than 10 cases per 1000 annual parasite incidence (annual incidence > 10‰), parasite genetics can be used to accurately infer incidence and is consistent with superinfection-based hypotheses of malaria transmission. When transmission was < 10‰, many of the correlations between parasite genetics and incidence were reversed, which may reflect the disproportionate impact of importation and focal transmission on parasite genetics when local transmission levels are low.

Investigating the etiologies of non-malarial febrile illness in Senegal using metagenomic sequencing.

The worldwide decline in malaria incidence is revealing the extensive burden of non-malarial febrile illness (NMFI), which remains poorly understood and difficult to diagnose. To characterize NMFI in Senegal, we collected venous blood and clinical metadata in a cross-sectional study of febrile patients and healthy controls in a low malaria burden area. Using 16S and untargeted sequencing, we detected viral, bacterial, or eukaryotic pathogens in 23% (38/163) of NMFI cases. Bacteria were the most common, with relapsing fever Borrelia and spotted fever Rickettsia found in 15.5% and 3.8% of cases, respectively. Four viral pathogens were found in a total of 7 febrile cases (3.5%). Sequencing also detected undiagnosed Plasmodium, including one putative P. ovale infection. We developed a logistic regression model that can distinguish Borrelia from NMFIs with similar presentation based on symptoms and vital signs (F1 score: 0.823). These results highlight the challenge and importance of improved diagnostics, especially for Borrelia, to support diagnosis and surveillance.

Evaluating the performance of Plasmodium falciparum genetics for inferring National Malaria Control Program reported incidence in Senegal.

Genetic surveillance of the Plasmodium falciparum parasite shows great promise for helping National Malaria Control Programs (NMCPs) assess parasite transmission. Genetic metrics such as the frequency of polygenomic (multiple strain) infections, genetic clones, and the complexity of infection (COI, number of strains per infection) are correlated with transmission intensity. However, despite these correlations, it is unclear whether genetic metrics alone are sufficient to estimate clinical incidence. Here, we examined parasites from 3,147 clinical infections sampled between the years 2012-2020 through passive case detection (PCD) across 16 clinic sites spread throughout Senegal. Samples were genotyped with a 24 single nucleotide polymorphism (SNP) molecular barcode that detects parasite strains, distinguishes polygenomic (multiple strain) from monogenomic (single strain) infections, and identifies clonal infections. To determine whether genetic signals can predict incidence, we constructed a series of Poisson generalized linear mixed-effects models to predict the incidence level at each clinical site from a set of genetic metrics designed to measure parasite clonality, superinfection, and co-transmission rates. We compared the model-predicted incidence with the reported standard incidence data determined by the NMCP for each clinic and found that parasite genetic metrics generally correlated with reported incidence, with departures from expected values at very low annual incidence (<10/1000/annual [‰]). When transmission is greater than 10 cases per 1000 annual parasite incidence (annual incidence >10 ‰), parasite genetics can be used to accurately infer incidence and is consistent with superinfection-based hypotheses of malaria transmission. When transmission was <10 ‰, we found that many of the correlations between parasite genetics and incidence were reversed, which we hypothesize reflects the disproportionate impact of importation and focal transmission on parasite genetics when local transmission levels are low.

Malaria surveillance reveals parasite relatedness, signatures of selection, and correlates of transmission across Senegal.

We here analyze data from the first year of an ongoing nationwide program of genetic surveillance of Plasmodium falciparum parasites in Senegal. The analysis is based on 1097 samples collected at health facilities during passive malaria case detection in 2019; it provides a baseline for analyzing parasite genetic metrics as they vary over time and geographic space. The study's goal was to identify genetic metrics that were informative about transmission intensity and other aspects of transmission dynamics, focusing on measures of genetic relatedness between parasites. We found the best genetic proxy for local malaria incidence to be the proportion of polygenomic infections (those with multiple genetically distinct parasites), although this relationship broke down at low incidence. The proportion of related parasites was less correlated with incidence while local genetic diversity was uninformative. The type of relatedness could discriminate local transmission patterns: two nearby areas had similarly high fractions of relatives, but one was dominated by clones and the other by outcrossed relatives. Throughout Senegal, 58% of related parasites belonged to a single network of relatives, within which parasites were enriched for shared haplotypes at known and suspected drug resistance loci and at one novel locus, reflective of ongoing selection pressure.

Improving diagnosis of non-malarial fevers in Senegal: Borrelia and the contribution of tick-borne bacteria.

The worldwide decline in malaria incidence is revealing the extensive burden of non-malarial febrile illness (NMFI), which remains poorly understood and difficult to diagnose. To characterize NMFI in Senegal, we collected venous blood and clinical metadata from febrile patients and healthy controls in a low malaria burden area. Using 16S and unbiased sequencing, we detected viral, bacterial, or eukaryotic pathogens in 29% of NMFI cases. Bacteria were the most common, with relapsing fever Borrelia and spotted fever Rickettsia found in 15% and 3.7% of cases, respectively. Four viral pathogens were found in a total of 7 febrile cases (3.5%). Sequencing also detected undiagnosed Plasmodium, including one putative P. ovale infection. We developed a logistic regression model to distinguish Borrelia from NMFIs with similar presentation based on symptoms and vital signs. These results highlight the challenge and importance of improved diagnostics, especially for Borrelia, to support diagnosis and surveillance.

Ex vivo RSA and pfkelch13 targeted-amplicon deep sequencing reveal parasites susceptibility to artemisinin in Senegal, 2017.

BACKGROUND: Malaria control is highly dependent on the effectiveness of artemisinin-based combination therapy (ACT), the current frontline malaria curative treatment. Unfortunately, the emergence and spread of parasites resistant to artemisinin (ART) derivatives in Southeast Asia and South America, and more recently in Rwanda and Uganda (East Africa), compromise their long-term use in sub-Saharan Africa, where most malaria deaths occur. METHODS: Here, ex vivo susceptibility to dihydroartemisinin (DHA) was evaluated from 38 Plasmodium falciparum isolates collected in 2017 in Thiès (Senegal) expressed in the Ring-stage Survival Assay (RSA). Both major and minor variants were explored in the three conserved-encoding domains of the pfkelch13 gene, the main determinant of ART resistance using a targeted-amplicon deep sequencing (TADS) approach. RESULTS: All samples tested in the ex vivo RSA were found to be susceptible to DHA (parasite survival rate < 1%). The non-synonymous mutations K189T and K248R in pfkelch13 were observed each in one isolate, as major (99%) or minor (5%) variants, respectively. CONCLUSION: The results suggest that ART is still fully effective in the Thiès region of Senegal in 2017. Investigations combining ex vivo RSA and TADS are a useful approach for monitoring ART resistance in Africa.

Malaria surveillance reveals parasite relatedness, signatures of selection, and correlates of transmission across Senegal.

Parasite genetic surveillance has the potential to play an important role in malaria control. We describe here an analysis of data from the first year of an ongoing, nationwide program of genetic surveillance of Plasmodium falciparum parasites in Senegal, intended to provide actionable information for malaria control efforts. Looking for a good proxy for local malaria incidence, we found that the best predictor was the proportion of polygenomic infections (those with multiple genetically distinct parasites), although that relationship broke down in very low incidence settings (r = 0.77 overall). The proportion of closely related parasites in a site was more weakly correlated ( r = -0.44) with incidence while the local genetic diversity was uninformative. Study of related parasites indicated their potential for discriminating local transmission patterns: two nearby study areas had similarly high fractions of relatives, but one area was dominated by clones and the other by outcrossed relatives. Throughout the country, 58% of related parasites proved to belong to a single network of relatives, within which parasites were enriched for shared haplotypes at known and suspected drug resistance loci as well as at one novel locus, reflective of ongoing selection pressure.

Ex vivo RSA and Pfkelch13 targeted-amplicon deep sequencing reveal parasites susceptibility to artemisinin in Senegal, 2017.

INTRODUCTION: Malaria control is highly dependent on the effectiveness of artemisinin-based combination therapies (ACTs), the current frontline malaria curative treatments. Unfortunately, the emergence and spread of parasites resistant to artemisinin (ART) derivatives in Southeast Asia and South America, and more recently in Rwanda and Uganda (East Africa), compromise their long-term use in Sub-Saharan Africa where most malaria deaths occur. METHODS: Here, we evaluated ex vivo susceptibility to dihydroartemisinin (DHA) from 38 P. falciparum isolates collected in 2017 in Thiès (Senegal) expressed with the Ring-stage Survival Assay (RSA). We explored major and minor variants in the full Pfkelch13 gene, the main determinant of ART resistance using a targeted-amplicon deep sequencing (TADS) approach. RESULTS: All samples tested in the ex vivo RSA were found to be susceptible to DHA. Both non-synonymous mutations K189T and K248R were observed each in one isolate, as major (99%) or minor (5%) variants, respectively. CONCLUSION: Altogether, investigations combining ex vivo RSA and TADS are a useful approach for monitoring ART resistance in Africa.

Genomic investigation of a dengue virus outbreak in Thiès, Senegal, in 2018.

Dengue virus is a major and rapidly growing public health concern in tropic and subtropic regions across the globe. In late 2018, Senegal experienced its largest dengue virus outbreak to date, covering several regions. However, little is known about the genetic diversity of dengue virus (DENV) in Senegal. Here we report complete viral genomes from 17 previously undetected DENV cases from the city of Thiès. In total we identified 19 cases of DENV in a cohort of 198 individuals with fever collected in October and November 2018. We detected 3 co-circulating serotypes; DENV 3 was the most frequent accounting for 11/17 sequences (65%), 4 (23%) were DENV2 and 2 (12%) were DENV1. Sequences were most similar to recent sequences from West Africa, suggesting ongoing local circulation of viral populations; however, detailed inference is limited by the scarcity of available genomic data. We did not find clear associations with reported clinical signs or symptoms, highlighting the importance of testing for diagnosing febrile diseases. Overall, these findings expand the known range of DENV in Senegal, and underscore the need for better genomic characterization of DENV in West Africa.