Direct sequencing of Leishmania donovani from patients in Garissa County, Northern Kenya, reveals a newly emerging intra-specific hybrid genotype

Introduction: Leishmaniasis is endemic in many countries, including Kenya. Despite the rising frequency of visceral leishmaniasis (VL) outbreaks in Garissa County, northeastern Kenya, there is limited data on the genetic structure and epidemiology of Leishmania parasites originating from this region. This study used molecular methods to characterize the Leishmania parasites collected at Garissa County Referral Hospital during the 2019 - 2022 VL outbreak.Methods: 286 blood samples, collected from patients suspected of having VL at Garissa County Referral Hospital between 2019 and 2022, were used. Leishmania parasites were screened at genus level by a quantitative real-time PCR assay targeting the arginine permease gene AAP3 (AAP3-qRT-PCR). Species identification and targeted gene sequencing were made on Illumina MiSeq using PCR amplicons of Hsp70 gene and ITS regions. Whole genome sequencing (WGS) was performed directly on eight selected clinical samples using a target enrichment method after which data was analyzed using phylogenomic tools.Results: By AAP3-qRT-PCR, 128/286 (45%) blood specimens were determined to have Leishmania parasites. We obtained 86 Hsp70 and 79 ITS sequences that phylogenetically clustered with the L. donovani species complex. By WGS, the 8 selected samples belonged to L. donovani s.s., and clustered in two separate groups: one similar to the previously reported L. donovani group 5 and the other constituted a new and intra-specific hybrid genetic variant not reported previously. A scan of genomic signatures of resistance to anti-leishmanaia compounds identified SNPs that code for resistance to Antimony, Amphotericine B and Miltefosine in all the 8 Kenyan samples.Conclusion: This study reveals the complex nature of Leishmania epidemiology in Kenya and sheds light on the genomic polymorphism of L. donovani in the region, which in turn, may explain the evolving threat of VL in the region. As caveat, the genomic signatures of drug resistance genes that were identified should be interpreted with caution until their functional implication is clarified in future studies.