Plasmodium falciparum population dynamics in East Africa and genomic surveillance along the Kenya-Uganda border.

East African countries accounted for ~ 10% of all malaria prevalence worldwide in 2022, with an estimated 23.8 million cases and > 53,000 deaths. Despite recent increases in malaria incidence, high-resolution genome-wide analyses of Plasmodium parasite populations are sparse in Kenya, Tanzania, and Uganda. The Kenyan-Ugandan border region is a particular concern, with Uganda confirming the emergence and spread of artemisinin resistant P. falciparum parasites. To establish genomic surveillance along the Kenyan-Ugandan border and analyse P. falciparum population dynamics within East Africa, we generated whole-genome sequencing (WGS) data for 38 parasites from Bungoma, Western Kenya. These sequences were integrated into a genomic analysis of available East African isolate data (n = 599) and revealed parasite subpopulations with distinct genetic structure and diverse ancestral origins. Ancestral admixture analysis of these subpopulations alongside isolates from across Africa (n = 365) suggested potential independent ancestral populations from other major African populations. Within isolates from Western Kenya, the prevalence of biomarkers associated with chloroquine resistance (e.g. Pfcrt K76T) were significantly reduced compared to wider East African populations and a single isolate contained the PfK13 V568I variant, potentially linked to reduced susceptibility to artemisinin. Overall, our work provides baseline WGS data and analysis for future malaria genomic surveillance in the region.

Knockdown resistance mutations in Phlebotomus argentipes sand flies in Bihar, India.

BACKGROUND: Vector control based on indoor residual spraying (IRS) is one of the main components of the visceral leishmaniasis (VL) elimination programme in India. Dichlorodiphenyltrichloroethane (DDT) was used for IRS until 2015 and was later replaced by the synthetic pyrethroid alpha-cypermethrin. Both classes of insecticides share the same target site, the voltage-gated sodium channel (Vgsc). As high levels of resistance to DDT have been documented in the local sand fly vector, Phlebotomus argentipes, it is possible that mutations in the Vgsc gene could provide resistance to alpha-cypermethrin, affecting current IRS pyrethroid-based vector control. METHODS: This study aimed to compare frequencies of knockdown resistance (kdr) mutations in Vgsc between two sprayed and two unsprayed villages in Bihar state, India, which had the highest VL burden of the four endemic states. Across four villages, 350 female P. argentipes were collected as part of a 2019 molecular xenomonitoring study. DNA was extracted and used for sequence analysis of the IIS6 fragment of the Vgsc gene to assess the presence of kdr mutations. RESULTS: Mutations were identified at various positions, most frequently at codon 1014, a common site known to be associated with insecticide resistance in mosquitoes and sand flies. Significant inter-village variation was observed, with sand flies from Dharampur, an unsprayed village, showing a significantly higher proportion of wild-type alleles (55.8%) compared with the three other villages (8.5-14.3%). The allele differences observed across the four villages may result from selection pressure caused by previous exposure to DDT. CONCLUSIONS: While DDT resistance has been reported in Bihar, P. argentipes is still susceptible to pyrethroids. However, the presence of kdr mutations in sand flies could present a threat to IRS used for VL control in endemic villages in India. Continuous surveillance of vector bionomics and insecticide resistance, using bioassays and target genotyping, is required to inform India's vector control strategies and to ensure the VL elimination target is reached and sustained.

Uncovering the genetic diversity in Aedes aegypti insecticide resistance genes through global comparative genomics.

Aedes aegypti is vector of many arboviruses including Zika, dengue, yellow fever, West Nile, and Chikungunya. Its control efforts are hampered by widespread insecticide resistance reported in the Americas and Asia, while data from Africa is more limited. Here we use publicly available 729 Ae. aegypti whole-genome sequencing samples from 15 countries, including nine in Africa, to investigate the genetic diversity in four insecticide resistance linked genes: ace-1, GSTe2, rdl and vgsc. Apart from vgsc, the other genes have been less investigated in Ae. aegypti, and almost no genetic diversity information is available. Among the four genes, we identified 1,829 genetic variants including 474 non-synonymous substitutions, some of which have been previously documented, as well as putative copy number variations in GSTe2 and vgsc. Global insecticide resistance phenotypic data demonstrated variable resistance in geographic areas with resistant genotypes. Overall, our work provides the first global catalogue and geographic distribution of known and new amino-acid mutations and duplications that can be used to guide the identification of resistance drivers in Ae. aegypti and thereby support monitoring efforts and strategies for vector control.

Multi-platform whole genome sequencing for tuberculosis clinical and surveillance applications.

Whole genome sequencing (WGS) of Mycobacterium tuberculosis offers valuable insights for tuberculosis (TB) control. High throughput platforms like Illumina and Oxford Nanopore Technology (ONT) are increasingly used globally, although ONT is known for higher error rates and is less established for genomic studies. Here we present a study comparing the sequencing outputs of both Illumina and ONT platforms, analysing DNA from 59 clinical isolates in highly endemic TB regions of Thailand. The resulting sequence data were used to profile the M. tuberculosis pairs for their lineage, drug resistance and presence in transmission chains, and were compared to publicly available WGS data from Thailand (n = 1456). Our results revealed isolates that are predominantly from lineages 1 and 2, with consistent drug resistance profiles, including six multidrug-resistant strains; however, analysis of ONT data showed longer phylogenetic branches, emphasising the technologies higher error rate. An analysis incorporating the larger dataset identified fifteen of our samples within six potential transmission clusters, including a significant clade of 41 multi-drug resistant isolates. ONT's extended sequences also revealed strain-specific structural variants in pe/ppe genes (e.g. ppe50), which are candidate loci for vaccine development. Despite some limitations, our results show that ONT sequencing is a promising approach for TB genomic research, supporting precision medicine and decision-making in areas with less developed infrastructure, which is crucial for tackling the disease's significant regional burden.

New reference genomes to distinguish the sympatric malaria parasites, Plasmodium ovale curtisi and Plasmodium ovale wallikeri.

Despite Plasmodium ovale curtisi (Poc) and wallikeri (Pow) being important human-infecting malaria parasites that are widespread across Africa and Asia, little is known about their genome diversity. Morphologically identical, Poc and Pow are indistinguishable and commonly misidentified. Recent rises in the incidence of Poc/Pow infections have renewed efforts to address fundamental knowledge gaps in their biology, and to develop diagnostic tools to understand their epidemiological dynamics and malaria burden. A major roadblock has been the incompleteness of available reference assemblies (PocGH01, PowCR01; ~ 33.5 Mbp). Here, we applied multiple sequencing platforms and advanced bioinformatics tools to generate new reference genomes, Poc221 (South Sudan; 36.0 Mbp) and Pow222 (Nigeria; 34.3 Mbp), with improved nuclear genome contiguity (> 4.2 Mbp), annotation and completeness (> 99% Plasmodium spp., single copy orthologs). Subsequent sequencing of 6 Poc and 15 Pow isolates from Africa revealed a total of 22,517 and 43,855 high-quality core genome SNPs, respectively. Genome-wide levels of nucleotide diversity were determined to be 2.98 × 10-4 (Poc) and 3.43 × 10-4 (Pow), comparable to estimates for other Plasmodium species. Overall, the new reference genomes provide a robust foundation for dissecting the biology of Poc/Pow, their population structure and evolution, and will contribute to uncovering the recombination barrier separating these species.

Genomic surveillance of Anopheles mosquitoes on the Bijagós Archipelago using custom targeted amplicon sequencing identifies mutations associated with insecticide resistance.

BACKGROUND: Insecticide resistance is reducing the efficacy of vector control interventions, consequently threatening efforts to control vector-borne diseases, including malaria. Investigating the prevalence of molecular markers of resistance is a useful tool for monitoring the spread of insecticide resistance in disease vectors. The Bijagós Archipelago (Bijagós) in Guinea-Bissau is a region of stable malaria transmission where insecticide-treated nets are the mainstay for malaria control. However, the prevalence of molecular markers of insecticide resistance in malaria vectors is not well understood. METHODS: A total of 214 Anopheles mosquitoes were analysed from 13 islands across the Bijagós. These mosquitoes were collected using CDC light traps in November 2019, during the peak malaria transmission season. High-throughput multiplex amplicon sequencing was used to investigate the prevalence of 17 different molecular markers associated with insecticide resistance in four genes: vgsc, rdl, ace1 and gste2. RESULTS: Of the 17 screened mutations, four were identified in mosquitoes from the Bijagós: vgsc L995F (12.2%), N1570Y (6.2%) and A1746S (0.7%) and rdl A269G (1.1%). This study is the first to report the L995F knock-down resistance (kdr)-west allele in Anopheles melas on the Archipelago. An additional eight non-synonymous single-nucleotide polymorphisms were identified across the four genes which have not been described previously. The prevalences of the vgsc L995F and N1570Y mutations were higher on Bubaque Island than on the other islands in this study; Bubaque is the most populous island in the archipelago, with the greatest population mobility and connection to continental Guinea-Bissau. CONCLUSIONS: This study provides the first surveillance data for genetic markers present in malaria vectors from islands across the Bijagós Archipelago. Overall prevalence of insecticide resistance mutations was found to be low. However, the identification of the vgsc L995F and N1570Y mutations associated with pyrethroid resistance warrants further monitoring. This is particularly important as the mainstay of malaria control on the islands is the use of pyrethroid insecticide-treated nets.

A computational method for predicting the most likely evolutionary trajectories in the stepwise accumulation of resistance mutations.

Pathogen evolution of drug resistance often occurs in a stepwise manner via the accumulation of multiple mutations that in combination have a non-additive impact on fitness, a phenomenon known as epistasis. The evolution of resistance via the accumulation of point mutations in the DHFR genes of Plasmodium falciparum (Pf) and Plasmodium vivax (Pv) has been studied extensively and multiple studies have shown epistatic interactions between these mutations determine the accessible evolutionary trajectories to highly resistant multiple mutations. Here, we simulated these evolutionary trajectories using a model of molecular evolution, parameterised using Rosetta Flex ddG predictions, where selection acts to reduce the target-drug binding affinity. We observe strong agreement with pathways determined using experimentally measured IC50 values of pyrimethamine binding, which suggests binding affinity is strongly predictive of resistance and epistasis in binding affinity strongly influences the order of fixation of resistance mutations. We also infer pathways directly from the frequency of mutations found in isolate data, and observe remarkable agreement with the most likely pathways predicted by our mechanistic model, as well as those determined experimentally. This suggests mutation frequency data can be used to intuitively infer evolutionary pathways, provided sufficient sampling of the population.

Rapid profiling of Plasmodium parasites from genome sequences to assist malaria control.

BACKGROUND: Malaria continues to be a major threat to global public health. Whole genome sequencing (WGS) of the underlying Plasmodium parasites has provided insights into the genomic epidemiology of malaria. Genome sequencing is rapidly gaining traction as a diagnostic and surveillance tool for clinical settings, where the profiling of co-infections, identification of imported malaria parasites, and detection of drug resistance are crucial for infection control and disease elimination. To support this informatically, we have developed the Malaria-Profiler tool, which rapidly (within minutes) predicts Plasmodium species, geographical source, and resistance to antimalarial drugs directly from WGS data. RESULTS: The online and command line versions of Malaria-Profiler detect ~ 250 markers from genome sequences covering Plasmodium speciation, likely geographical source, and resistance to chloroquine, sulfadoxine-pyrimethamine (SP), and other anti-malarial drugs for P. falciparum, but also providing mutations for orthologous resistance genes in other species. The predictive performance of the mutation library was assessed using 9321 clinical isolates with WGS and geographical data, with most being single-species infections (P. falciparum 7152/7462, P. vivax 1502/1661, P. knowlesi 143/151, P. malariae 18/18, P. ovale ssp. 5/5), but co-infections were identified (456/9321; 4.8%). The accuracy of the predicted geographical profiles was high to both continental (96.1%) and regional levels (94.6%). For P. falciparum, markers were identified for resistance to chloroquine (49.2%; regional range: 24.5% to 100%), sulfadoxine (83.3%; 35.4- 90.5%), pyrimethamine (85.4%; 80.0-100%) and combined SP (77.4%). Markers associated with the partial resistance of artemisinin were found in WGS from isolates sourced from Southeast Asia (30.6%). CONCLUSIONS: Malaria-Profiler is a user-friendly tool that can rapidly and accurately predict the geographical regional source and anti-malarial drug resistance profiles across large numbers of samples with WGS data. The software is flexible with modifiable bioinformatic pipelines. For example, it is possible to select the sequencing platform, display specific variants, and customise the format of outputs. With the increasing application of next-generation sequencing platforms on Plasmodium DNA, Malaria-Profiler has the potential to be integrated into point-of-care and surveillance settings, thereby assisting malaria control. Malaria-Profiler is available online (bioinformatics.lshtm.ac.uk/malaria-profiler) and as standalone software ( https://github.com/jodyphelan/malaria-profiler ).

Functional genetic variation in pe/ppe genes contributes to diversity in Mycobacterium tuberculosis lineages and potential interactions with the human host.

Introduction: Around 10% of the coding potential of Mycobacterium tuberculosisis constituted by two poorly understood gene families, the pe and ppe loci, thought to be involved in host-pathogen interactions. Their repetitive nature and high GC content have hindered sequence analysis, leading to exclusion from whole-genome studies. Understanding the genetic diversity of pe/ppe families is essential to facilitate their potential translation into tools for tuberculosis prevention and treatment. Methods: To investigate the genetic diversity of the 169 pe/ppe genes, we performed a sequence analysis across 73 long-read assemblies representing seven different lineages of M. tuberculosis and M. bovis BCG. Individual pe/ppe gene alignments were extracted and diversity and conservation across the different lineages studied. Results: The pe/ppe genes were classified into three groups based on the level of protein sequence conservation relative to H37Rv, finding that >50% were conserved, with indels in pe_pgrs and ppe_mptr sub-families being major drivers of structural variation. Gene rearrangements, such as duplications and gene fusions, were observed between pe and pe_pgrs genes. Inter-lineage diversity revealed lineage-specific SNPs and indels. Discussion: The high level of pe/ppe genes conservation, together with the lineage-specific findings, suggest their phylogenetic informativeness. However, structural variants and gene rearrangements differing from the reference were also identified, with potential implications for pathogenicity. Overall, improving our knowledge of these complex gene families may have insights into pathogenicity and inform the development of much-needed tools for tuberculosis control.

Dried Blood Spots (DBS): A suitable alternative to using whole blood samples for diagnostic testing of visceral leishmaniasis in the post-elimination era.

BACKGROUND: Serum or whole blood collection, processing, transport and storage still present significant challenges in low resource settings where mass surveillance is required to sustain disease elimination. Therefore, in this study, we explored the diagnostic efficacy of dried blood spots (DBS) as a minimally invasive and potentially cost-effective alternative sampling technique to whole blood sampling procedures for subsequent detection of Leishmania donovani antibodies or DNA. METHODOLOGY AND PRINCIPAL FINDINGS: Archived serum, DNA samples from whole blood of visceral leishmaniasis (VL) cases and healthy controls, and DBS from corresponding cases and controls, were used. Both molecular and serological assays were optimized to detect L. donovani antibodies or DNA in DBS elute and results were compared against those obtained with whole blood. Serological assays (both rK28 ELISA and rK39 ELISA) of DBS samples showed sensitivity and specificity of 100% and had excellent agreement with results from whole blood samples (kappa value ranged from 0.98-1). Bland-Altman analysis of OD values from rK28-ELISA with DBS elute and patients' serum showed an excellent agreement (ICC = 0.9) whereas a good agreement (ICC = 0.8) was observed in the case of rK39-ELISA. However, qPCR and RPA of DBS samples had a diminished sensitivity of 76% and 68%, respectively, and poor agreement was observed with the whole blood samples. CONCLUSION: Our results demonstrate that DBS offer excellent diagnostic efficiency for serological assays and represent a viable alternative to whole blood sampling procedures.

Detection of insecticide resistance markers in Anopheles funestus from the Democratic Republic of the Congo using a targeted amplicon sequencing panel.

Vector control strategies have been successful in reducing the number of malaria cases and deaths globally, but the spread of insecticide resistance represents a significant threat to disease control. Insecticide resistance has been reported across Anopheles (An.) vector populations, including species within the An. funestus group. These mosquitoes are responsible for intense malaria transmission across sub-Saharan Africa, including in the Democratic Republic of the Congo (DRC), a country contributing > 12% of global malaria infections and mortality events. To support the continuous efficacy of vector control strategies, it is essential to monitor insecticide resistance using molecular surveillance tools. In this study, we developed an amplicon sequencing ("Amp-seq") approach targeting An. funestus, and using multiplex PCR, dual index barcoding, and next-generation sequencing for high throughput and low-cost applications. Using our Amp-seq approach, we screened 80 An. funestus field isolates from the DRC across a panel of nine genes with mutations linked to insecticide resistance (ace-1, CYP6P4, CYP6P9a, GSTe2, vgsc, and rdl) and mosquito speciation (cox-1, mtND5, and ITS2). Amongst the 18 non-synonymous mutations detected, was N485I, in the ace-1 gene associated with carbamate resistance. Overall, our panel represents an extendable and much-needed method for the molecular surveillance of insecticide resistance in An. funestus populations.

Mixed infections in genotypic drug-resistant Mycobacterium tuberculosis.

Tuberculosis disease (TB), caused by Mycobacterium tuberculosis, is a major global public health problem, resulting in more than 1 million deaths each year. Drug resistance (DR), including multi-drug (MDR-TB), is making TB control difficult and accounts for 16% of new and 48% of previously treated cases. To further complicate treatment decision-making, many clinical studies have reported patients harbouring multiple distinct strains of M. tuberculosis across the main lineages (L1 to L4). The extent to which drug-resistant strains can be deconvoluted within mixed strain infection samples is understudied. Here, we analysed M. tuberculosis isolates with whole genome sequencing data (n = 50,723), which covered the main lineages (L1 9.1%, L2 27.6%, L3 11.8%, L4 48.3%), with genotypic resistance to isoniazid (HR-TB; n = 9546 (29.2%)), rifampicin (RR-TB; n = 7974 (24.4%)), and at least MDR-TB (n = 5385 (16.5%)). TB-Profiler software revealed 531 (1.0%) isolates with potential mixed sub-lineage infections, including some with DR mutations (RR-TB 21/531; HR-TB 59/531; at least MDR-TB 173/531). To assist with the deconvolution of such mixtures, we adopted and evaluated a statistical Gaussian Mixture model (GMM) approach. By simulating 240 artificial mixtures of different ratios from empirical data across L1 to L4, a GMM approach was able to accurately estimate the DR profile of each lineage, with a low error rate for the estimated mixing proportions (mean squared error 0.012) and high accuracy for the DR predictions (93.5%). Application of the GMM model to the clinical mixtures (n = 531), found that 33.3% (188/531) of samples consisted of DR and sensitive lineages, 20.2% (114/531) consisted of lineages with only DR mutations, and 40.6% (229/531) consisted of lineages with genotypic pan-susceptibility. Overall, our work demonstrates the utility of combined whole genome sequencing data and GMM statistical analysis approaches for providing insights into mono and mixed M. tuberculosis infections, thereby potentially assisting diagnosis, treatment decision-making, drug resistance and transmission mapping for infection control.

Comparison of collection methods for Phlebotomus argentipes sand flies to use in a molecular xenomonitoring system for the surveillance of visceral leishmaniasis.

BACKGROUND: The kala-azar elimination programme has resulted in a significant reduction in visceral leishmaniasis (VL) cases across the Indian Subcontinent. To detect any resurgence of transmission, a sensitive cost-effective surveillance system is required. Molecular xenomonitoring (MX), detection of pathogen DNA/RNA in vectors, provides a proxy of human infection in the lymphatic filariasis elimination programme. To determine whether MX can be used for VL surveillance in a low transmission setting, large numbers of the sand fly vector Phlebotomus argentipes are required. This study will determine the best method for capturing P. argentipes females for MX. METHODOLOGY/PRINCIPAL FINDINGS: The field study was performed in two programmatic and two non-programmatic villages in Bihar, India. A total of 48 households (12/village) were recruited. Centers for Disease Control and Prevention light traps (CDC-LTs) were compared with Improved Prokopack (PKP) and mechanical vacuum aspirators (MVA) using standardised methods. Four 12x12 Latin squares, 576 collections, were attempted (12/house, 144/village,192/method). Molecular analyses of collections were conducted to confirm identification of P. argentipes and to detect human and Leishmania DNA. Operational factors, such as time burden, acceptance to householders and RNA preservation, were also considered. A total of 562 collections (97.7%) were completed with 6,809 sand flies captured. Females comprised 49.0% of captures, of which 1,934 (57.9%) were identified as P. argentipes. CDC-LTs collected 4.04 times more P. argentipes females than MVA and 3.62 times more than PKP (p<0.0001 for each). Of 21,735 mosquitoes in the same collections, no significant differences between collection methods were observed. CDC-LTs took less time to install and collect than to perform aspirations and their greater yield compensated for increased sorting time. No significant differences in Leishmania RNA detection and quantitation between methods were observed in experimentally infected sand flies maintained in conditions simulating field conditions. CDC-LTs were favoured by householders. CONCLUSIONS/SIGNIFICANCE: CDC-LTs are the most useful collection tool of those tested for MX surveillance since they collected higher numbers of P. argentipes females without compromising mosquito captures or the preservation of RNA. However, capture rates are still low.

Large-scale reference-free analysis of flavivirus sequences in Aedes aegypti whole genome DNA sequencing data.

Flaviviruses are a diverse group of RNA viruses, which include the etiological agents of Zika, dengue and yellow fever that are transmitted by mosquitoes. Flaviviruses do not encode reverse transcriptase and cannot reverse transcribe into DNA, yet DNA sequences of flaviviruses are found both integrated in the chromosomes of Aedes aegypti mosquitoes and as extrachromosomal sequences. We have previously examined the Ae. aegypti reference genome to identify flavivirus integrations and analyzed conservation of these sequences among whole-genome data of 464 Ae. aegypti collected across 10 countries globally. Here, we extended this analysis by identifying flavivirus sequences in these samples independently of the Ae. aegypti reference assembly. Our aim was to identify the complete set of viral sequences, including those absent in the reference genome, and their geographical distribution. We compared the identified sequences using BLASTn and applied machine learning methods to identify clusters of similar sequences. Apart from clusters of sequences that correspond to the four viral integration events that we had previously described, we identified 19 smaller clusters. The only cluster with a strong geographic association consisted of Cell-fusing agent virus-like sequences specific to Thailand. The remaining clusters did not have a geographic association and mostly consisted of near identical short sequences without strong similarity to any known flaviviral genomes. The short read sequencing data did not permit us to determine whether identified sequences were extrachromosomal or integrated into Ae. aegypti chromosomes. Our results suggest that Liverpool strain and field Ae. aegypti mosquitoes have a similar variety of conserved flaviviral DNA, whose functional role should be investigated in follow-up studies.

Genomic and functional portrait of multidrug-resistant, hydrogen sulfide (H2S)-producing variants of Escherichia coli.

Atypical Escherichia coli forms exhibit unusual characteristics compared to typical strains. The H2S-producing variants of some atypical E. coli strains cause a wide range of illnesses in humans and animals. However, there are sparse reports on such strains worldwide. We performed whole-genome sequencing (WGS) and detailed characterization of four H2S-producing E. coli variants from poultry and human clinical sources in Dhaka, Bangladesh. All four isolates were confirmed as E. coli using biochemical tests and genomic analysis, and were multidrug-resistant (MDR). WGS analysis including an additional Chinese strain, revealed diverse STs among the five H2S-producing E. coli genomes, with clonal complex ST10 being detected in 2 out of 5 genomes. The predominant phylogroup detected was group A (n = 4/5). The blaTEM1B (n = 5/5) was the most predominant extended-spectrum beta-lactamase (ESBL) gene, followed by different alleles of blaCTX-M (blaCTX-M -55,-65,-123; n = 3/5). Multiple plasmid replicons were detected, with IncX being the most common. One E. coli strain was classified as enteropathogenic E. coli. The genomes of all five isolates harbored five primary and four secondary function genes related to H2S production. These findings suggest the potential of these isolates to cause disease and spread antibiotic resistance. Therefore, such atypical E. coli forms should be included in differential diagnosis to understand the pathogenicity, antimicrobial resistance and evolution of H2S-producing E. coli.

Case report: A successfully treated case of community-acquired urinary tract infection due to Klebsiella aerogenes in Bangladesh.

Klebsiella aerogenes, a nosocomial pathogen, is increasingly associated with extensive drug resistance and virulence profiles. It is responsible for high morbidity and mortality. This report describes the first successfully treated case of community-acquired urinary tract infection (UTI) caused by Klebsiella aerogenes in an elderly housewife with Type-2 diabetes (T2D) from Dhaka, Bangladesh. The patient was empirically treated with intravenous ceftriaxone (500 mg/8 h). However, she did not respond to the treatment. The urine culture and sensitivity tests, coupled with bacterial whole-genome sequencing (WGS) and analysis, revealed the bacteria to be K. aerogenes which was extensively drug-resistant but was susceptible to carbapenems and polymyxins. Based on these findings, meropenem (500 mg/8 h) was administered to the patient, who then responded to the treatment and recovered successfully without having a relapse. This case raises awareness of the importance of diagnosis of not-so-common etiological agents, correct identification of the pathogens, and targeted antibiotic therapy. In conclusion, correctly identifying etiological agents of UTI using WGS approaches that are otherwise difficult to diagnose could help improve the identification of infectious agents and improve the management of infectious diseases.

Comparison of in silico predicted Mycobacterium tuberculosis spoligotypes and lineages from whole genome sequencing data.

Bacterial strain-types in the Mycobacterium tuberculosis complex underlie tuberculosis disease, and have been associated with drug resistance, transmissibility, virulence, and host-pathogen interactions. Spoligotyping was developed as a molecular genotyping technique used to determine strain-types, though recent advances in whole genome sequencing (WGS) technology have led to their characterization using SNP-based sub-lineage nomenclature. Notwithstanding, spoligotyping remains an important tool and there is a need to study the congruence between spoligotyping-based and SNP-based sub-lineage assignation. To achieve this, an in silico spoligotype prediction method ("Spolpred2") was developed and integrated into TB-Profiler. Lineage and spoligotype predictions were generated for > 28 k isolates and the overlap between strain-types was characterized. Major spoligotype families detected were Beijing (25.6%), T (18.6%), LAM (13.1%), CAS (9.4%), and EAI (8.3%), and these broadly followed known geographic distributions. Most spoligotypes were perfectly correlated with the main MTBC lineages (L1-L7, plus animal). Conversely, at lower levels of the sub-lineage system, the relationship breaks down, with only 65% of spoligotypes being perfectly associated with a sub-lineage at the second or subsequent levels of the hierarchy. Our work supports the use of spoligotyping (membrane or WGS-based) for low-resolution surveillance, and WGS or SNP-based systems for higher-resolution studies.

Drug resistance profiling of asymptomatic and low-density Plasmodium falciparum malaria infections on Ngodhe island, Kenya, using custom dual-indexing next-generation sequencing.

Malaria control initiatives require rapid and reliable methods for the detection and monitoring of molecular markers associated with antimalarial drug resistance in Plasmodium falciparum parasites. Ngodhe island, Kenya, presents a unique malaria profile, with lower P. falciparum incidence rates than the surrounding region, and a high proportion of sub-microscopic and low-density infections. Here, using custom dual-indexing and Illumina next generation sequencing, we generate resistance profiles on seventy asymptomatic and low-density P. falciparum infections from a mass drug administration program implemented on Ngodhe island between 2015 and 2016. Our assay encompasses established molecular markers on the Pfcrt, Pfmdr1, Pfdhps, Pfdhfr, and Pfk13 genes. Resistance markers for sulfadoxine-pyrimethamine were identified at high frequencies, including a quintuple mutant haplotype (Pfdhfr/Pfdhps: N51I, C59R, S108N/A437G, K540E) identified in 62.2% of isolates. The Pfdhps K540E biomarker, used to inform decision making for intermittent preventative treatment in pregnancy, was identified in 79.2% of isolates. Several variants on Pfmdr1, associated with reduced susceptibility to quinolones and lumefantrine, were also identified (Y184F 47.1%; D1246Y 16.0%; N86 98%). Overall, we have presented a low-cost and extendable approach that can provide timely genetic profiles to inform clinical and surveillance activities, especially in settings with abundant low-density infections, seeking malaria elimination.

High throughput human genotyping for variants associated with malarial disease outcomes using custom targeted amplicon sequencing.

Malaria has exhibited the strongest known selective pressure on the human genome in recent history and is the evolutionary driving force behind genetic conditions, such as sickle-cell disease, glucose-6-phosphatase deficiency, and some other erythrocyte defects. Genomic studies (e.g., The 1000 Genomes project) have provided an invaluable baseline for human genetics, but with an estimated two thousand ethno-linguistic groups thought to exist across the African continent, our understanding of the genetic differences between indigenous populations and their implications on disease is still limited. Low-cost sequencing-based approaches make it possible to target specific molecular markers and genes of interest, leading to potential insights into genetic diversity. Here we demonstrate the versatility of custom dual-indexing technology and Illumina next generation sequencing to generate a genetic profile of human polymorphisms associated with malaria pathology. For 100 individuals diagnosed with severe malaria in Northeast Tanzania, variants were successfully characterised on the haemoglobin subunit beta (HBB), glucose-6-phosphate dehydrogenase (G6PD), atypical chemokine receptor 1 (ACKR1) genes, and the intergenic Dantu genetic blood variant, then validated using pre-existing genotyping data. High sequencing coverage was observed across all amplicon targets in HBB, G6PD, ACKR1, and the Dantu blood group, with variants identified at frequencies previously observed within this region of Tanzania. Sequencing data exhibited high concordance rates to pre-existing genotyping data (> 99.5%). Our work demonstrates the potential utility of amplicon sequencing for applications in human genetics, including to personalise medicine and understand the genetic diversity of loci linked to important host phenotypes, such as malaria susceptibility.

Genome dynamics of high-risk resistant and hypervirulent Klebsiella pneumoniae clones in Dhaka, Bangladesh.

Klebsiella pneumoniae is recognized as an urgent public health threat because of the emergence of difficult-to-treat (DTR) strains and hypervirulent clones, resulting in infections with high morbidity and mortality rates. Despite its prominence, little is known about the genomic epidemiology of K. pneumoniae in resource-limited settings like Bangladesh. We sequenced genomes of 32 K. pneumoniae strains isolated from patient samples at the International Center for Diarrhoeal Disease Research, Bangladesh (icddr,b). Genome sequences were examined for their diversity, population structure, resistome, virulome, MLST, O and K antigens and plasmids. Our results revealed the presence of two K. pneumoniae phylogroups, namely KpI (K. pneumoniae) (97%) and KpII (K. quasipneumoniae) (3%). The genomic characterization revealed that 25% (8/32) of isolates were associated with high-risk multidrug-resistant clones, including ST11, ST14, ST15, ST307, ST231 and ST147. The virulome analysis confirmed the presence of six (19%) hypervirulent K. pneumoniae (hvKp) and 26 (81%) classical K. pneumoniae (cKp) strains. The most common ESBL gene identified was blaCTX-M-15 (50%). Around 9% (3/32) isolates exhibited a difficult-to-treat phenotype, harboring carbapenem resistance genes (2 strains harbored blaNDM-5 plus blaOXA-232, one isolate blaOXA-181). The most prevalent O antigen was O1 (56%). The capsular polysaccharides K2, K20, K16 and K62 were enriched in the K. pneumoniae population. This study suggests the circulation of the major international high-risk multidrug-resistant and hypervirulent (hvKp) K. pneumoniae clones in Dhaka, Bangladesh. These findings warrant immediate appropriate interventions, which would otherwise lead to a high burden of untreatable life-threatening infections locally.