Boosting pathogen genomics and bioinformatics workforce in Africa

Next-generation sequencing (NGS) of genomic data has established its fundamental value in public health surveillance, research and development, and precision medicine. In Africa, severe shortages of competent experts in genomics and bioinformatics, few opportunities for research, and inadequate genomic infrastructure have had a knock-on effect on the use of NGS technologies for research and public health practice. Several reasons-ranging from poor funding, inadequate infrastructure for training and practice, to brain drain-might partly account for the scarcity of genomics and bioinformatics expertise in the region. In recognition of these shortcomings and the importance of NGS genomic data, which was amplified during the COVID-19 pandemic in mid-2021, the Africa Centres for Disease Control and Prevention (Africa CDC) through the Africa Pathogen Genomics Initiative began building and expanding Africa's workforce in pathogen surveillance. By the end of 2022, the Africa CDC in collaboration with its partners and centres of excellence had trained 413 personnel, mostly from public health institutions, in 53 (96%) of 55 African Union Member States. Although this training has increased genomics, bioinformatics, and genomic epidemiology literacy, and genomic-informed pathogen surveillance, there is still a need for a strategic and sustainable public health workforce development in Africa.

Population structure and genetic connectivity of Plasmodium falciparum in pre-elimination settings of Southern Africa

To accelerate malaria elimination in the Southern African region by 2030, it is essential to prevent cross-border malaria transmission. However, countries within the region are highly interconnected due to human migration that aids in the movement of the parasite across geographical borders. It is therefore important to better understand Plasmodium falciparum transmission dynamics in the region, and identify major parasite source and sink populations, as well as cross-border blocks of high parasite connectivity. We performed a meta-analysis using collated parasite allelic data generated by microsatellite genotyping of malaria parasites from Namibia, Eswatini, South Africa, and Mozambique (N = 5,314). The overall number of unique alleles was significantly higher (P ≤ 0.01) in Namibia (mean A = 17.3 ± 1.46) compared to South Africa (mean A = 12.2 ± 1.22) and Eswatini (mean A = 13.3 ± 1.27, P ≤ 0.05), whilst the level of heterozygosity was not significantly different between countries. The proportion of polyclonal infections was highest for Namibia (77%), and lowest for Mozambique (64%). A was significant population structure was detected between parasites from the four countries, and patterns of gene flow showed that Mozambique was the major source area and Eswatini the major sink area of parasites between the countries. This study showed strong signals of parasite population structure and genetic connectivity between malaria parasite populations across national borders. This calls for strengthening the harmonization of malaria control and elimination efforts between countries in the southern African region. This data also proves its potential utility as an additional surveillance tool for malaria surveillance on both a national and regional level for the identification of imported cases and/or outbreaks, as well as monitoring for the potential spread of anti-malarial drug resistance as countries work towards malaria elimination.

Unlocking the Potential of Genomic Data to Inform Typhoid Fever Control Policy: Supportive Resources for Genomic Data Generation, Analysis, and Visualization

The global response to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic demonstrated the value of timely and open sharing of genomic data with standardized metadata to facilitate monitoring of the emergence and spread of new variants. Here, we make the case for the value of Salmonella Typhi (S. Typhi) genomic data and demonstrate the utility of freely available platforms and services that support the generation, analysis, and visualization of S. Typhi genomic data on the African continent and more broadly by introducing the Africa Centres for Disease Control and Prevention's Pathogen Genomics Initiative, SEQAFRICA, Typhi Pathogenwatch, TyphiNET, and the Global Typhoid Genomics Consortium.

SARS-CoV-2 Omicron variant of concern in the Seychelles: Introduction and spread

Background: The emergence of the Omicron variant of concern in late 2021 led to a resurgence of SARS-CoV-2 infections globally. By September 2022, Seychelles had experienced two major surges of SARS-CoV-2 infections driven by the Omicron variant. Here, we examine the genomic epidemiology of Omicron in the Seychelles between November 2021 and September 2022. Methods: We analysed 618 SARS-CoV-2 Omicron genomes identified in the Seychelles between November 2021 and September 2022 to infer virus introductions and local transmission patterns using phylogenetics and the ancestral state reconstruction approach. We then evaluated the impact of government coronavirus 2019 (COVID-19) countermeasures on the estimated number of viral introductions during the study period. Results: The genomes classified into 43 distinct Pango lineages. The first surge in Omicron cases (beginning November 2021 and peaking in January 2022) was predominated by the BA.1.1 lineage (59%) co-circulating with 11 other Omicron lineages. In the second surge (between April and June 2022), four lineages (BA.2, BA.2.10, BA.2.65 and BA.2.9) co-circulated and these were swiftly replaced by BA.5 subvariants in July 2022, which remained predominant through to September 2022. In the latter period, sporadic detections of BA.5 subvariants BQ.1, BE and BF were observed. We estimated 109 independent Omicron importations into Seychelles over the 11-month period, most of which occurred between December 2021 and March 2022 when strict government restrictions (SI>50%) were still in force. The districts Anse Royale, and Baie St. Anne Praslin appeared to be the major dispersal points fuelling local transmission. Conclusions: Our results suggest that the waves of Omicron infections in the Seychelles were driven by multiple lineages and multiple virus introductions. The introductions were followed by substantial local spread and successive lineage displacement that mirrored the global patterns.

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance

Investment in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequencing in Africa over the past year has led to a major increase in the number of sequences that have been generated and used to track the pandemic on the continent, a number that now exceeds 100,000 genomes. Our results show an increase in the number of African countries that are able to sequence domestically and highlight that local sequencing enables faster turnaround times and more-regular routine surveillance. Despite limitations of low testing proportions, findings from this genomic surveillance study underscore the heterogeneous nature of the pandemic and illuminate the distinct dispersal dynamics of variants of concern-particularly Alpha, Beta, Delta, and Omicron-on the continent. Sustained investment for diagnostics and genomic surveillance in Africa is needed as the virus continues to evolve while the continent faces many emerging and reemerging infectious disease threats. These investments are crucial for pandemic preparedness and response and will serve the health of the continent well into the 21st century.

Urgent need for a non-discriminatory and non-stigmatizing nomenclature for monkeypox virus

The current nomenclature for monkeypox virus is stigmatising and misleading. This Perspective article proposes a practical and neutral system of nomenclature that will allow efficient communication without the risk of further misconceptions, discrimination and stigmatisation.

A year of genomic surveillance reveals how the SARS-CoV-2 pandemic unfolded in Africa

The progression of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic in Africa has so far been heterogeneous, and the full impact is not yet well understood. In this study, we describe the genomic epidemiology using a dataset of 8746 genomes from 33 African countries and two overseas territories. We show that the epidemics in most countries were initiated by importations predominantly from Europe, which diminished after the early introduction of international travel restrictions. As the pandemic progressed, ongoing transmission in many countries and increasing mobility led to the emergence and spread within the continent of many variants of concern and interest, such as B.1.351, B.1.525, A.23.1, and C.1.1. Although distorted by low sampling numbers and blind spots, the findings highlight that Africa must not be left behind in the global pandemic response, otherwise it could become a source for new variants.

A review of the frequencies of Plasmodium falciparum Kelch 13 artemisinin resistance mutations in Africa

Artemisinin resistance (AR) emerged in South East Asia 13 years ago and the identification of the resistance conferring molecular marker, Plasmodium falciparum Kelch 13 (Pfk13), 7 years ago has provided an invaluable tool for monitoring AR in malaria endemic countries. Molecular Pfk13 surveillance revealed the resistance foci in the Greater Mekong Subregion, an independent emergence in Guyana, South America, and a low frequency of mutations in Africa. The recent identification of the R561H Pfk13 AR associated mutation in Tanzania, Uganda and in Rwanda, where it has been associated with delayed parasite clearance, should be a concern for the continent. In this review, we provide a summary of Pfk13 resistance associated propeller domain mutation frequencies across Africa from 2012 to 2020, to examine how many other countries have identified these mutations. Only four African countries reported a recent identification of the M476I, P553L, R561H, P574L, C580Y and A675V Pfk13 mutations at low frequencies and with no reports of clinical treatment failure, except for Rwanda. These mutations present a threat to malaria control across the continent, since the greatest burden of malaria remains in Africa. A rise in the frequency of these mutations and their spread would reverse the gains made in the reduction of malaria over the last 20 years, given the lack of new antimalarial treatments in the event artemisinin-based combination therapies fail. The review highlights the frequency of Pfk13 propeller domain mutations across Africa, providing an up-to-date perspective of Pfk13 mutations, and appeals for an urgent and concerted effort to monitoring antimalarial resistance markers in Africa and the efficacy of antimalarials by re-establishing sentinel surveillance systems.

Optimization of whole-genome sequencing of Plasmodium falciparum from low-density dried blood spot samples

Background

Whole-genome sequencing (WGS) is becoming increasingly useful to study the biology, epidemiology, and ecology of malaria parasites. Despite ease of sampling, DNA extracted from dried blood spots (DBS) has a high ratio of human DNA compared to parasite DNA, which poses a challenge for downstream genetic analyses. The effects of multiple methods for DNA extraction, digestion of methylated DNA, and amplification were evaluated on the quality and fidelity of WGS data recovered from DBS.

Methods

Low parasite density mock DBS samples were created, extracted either with Tween-Chelex or QIAamp, treated with or without McrBC, and amplified with one of three different amplification techniques (two sWGA primer sets and one rWGA). Extraction conditions were evaluated on performance of sequencing depth, percentiles of coverage, and expected SNP concordance.

Results

At 100 parasites/μL, Chelex-Tween-McrBC samples had higher coverage (5 × depth = 93% genome) than QIAamp extracted samples (5 × depth = 76% genome). The two evaluated sWGA primer sets showed minor differences in overall genome coverage and SNP concordance, with a newly proposed combination of 20 primers showing a modest improvement in coverage over those previously published.

Conclusions

Overall, Tween-Chelex extracted samples that were treated with McrBC digestion and are amplified using 6A10AD sWGA conditions had minimal dropout rate, higher percentages of coverage at higher depth, and more accurate SNP concordance than QiaAMP extracted samples. These findings extend the results of previously reported methods, making whole genome sequencing accessible to a larger number of low density samples that are commonly encountered in cross-sectional surveys.

Parasite genetic diversity reflects continued residual malaria transmission in Vhembe District, a hotspot in the Limpopo Province of South Africa

BACKGROUND

South Africa aims to eliminate malaria transmission by 2023. However, despite sustained vector control efforts and case management interventions, the Vhembe District remains a malaria transmission hotspot. To better understand Plasmodium falciparum transmission dynamics in the area, this study characterized the genetic diversity of parasites circulating within the Vhembe District.

METHODS

A total of 1153 falciparum-positive rapid diagnostic tests (RDTs) were randomly collected from seven clinics within the district, over three consecutive years (2016, 2017 and 2018) during the wet and dry malaria transmission seasons. Using 26 neutral microsatellite markers, differences in genetic diversity were described using a multiparameter scale of multiplicity of infection (MOI), inbreeding metric (Fws), number of unique alleles (A), expected heterozygosity (He), multilocus linkage disequilibrium (LD) and genetic differentiation, and were associated with temporal and geospatial variances.

RESULTS

A total of 747 (65%) samples were successfully genotyped. Moderate to high genetic diversity (mean He = 0.74 ± 0.03) was observed in the parasite population. This was ascribed to high allelic richness (mean A = 12.2 ± 1.2). The majority of samples (99%) had unique multi-locus genotypes, indicating high genetic diversity in the sample set. Complex infections were observed in 66% of samples (mean MOI = 2.13 ± 0.04), with 33% of infections showing high within-host diversity as described by the Fws metric. Low, but significant LD (standardised index of association, ISA = 0.08, P < 0.001) was observed that indicates recombination of distinct clones. Limited impact of temporal (F_ST range - 0.00005 to 0.0003) and spatial (F_ST = - 0.028 to 0.023) variation on genetic diversity existed during the sampling timeframe and study sites respectively.

CONCLUSIONS

Consistent with the Vhembe District's classification as a 'high' transmission setting within South Africa, P. falciparum diversity in the area was moderate to high and complex. This study showed that genetic diversity within the parasite population reflects the continued residual transmission observed in the Vhembe District. This data can be used as a reference point for the assessment of the effectiveness of on-going interventions over time, the identification of imported cases and/or outbreaks, as well as monitoring for the potential spread of anti-malarial drug resistance.

Understanding COVID-19 in Africa

In this Comment article, Sofonias Tessema and John Nkengasong provide an overview of the current state of the COVID-19 pandemic in Africa and the challenges posed by the triple burden of emerging, endemic and non-communicable diseases.

Africa Needs a New Public Health Order to Tackle Infectious Disease Threats

The SARS-CoV-2 pandemic has revealed that Africa needs a new public health order to be resilient, to adapt, and to cope with 21st-century disease threats. The new order will need strengthened continental and national public health institutions; local manufacturing of vaccines, therapeutics, and diagnostics; attraction, training, and retention of a public health workforce; and fostering of respectful local and international partnerships.

Sensitive, highly multiplexed sequencing of microhaplotypes from the Plasmodium falciparum heterozygome

Background

Targeted next-generation sequencing offers the potential for consistent, deep coverage of information-rich genomic regions to characterize polyclonal Plasmodium falciparum infections. However, methods to identify and sequence these genomic regions are currently limited.

Methods

A bioinformatic pipeline and multiplex methods were developed to identify and simultaneously sequence 100 targets and applied to dried blood spot (DBS) controls and field isolates from Mozambique. For comparison, whole-genome sequencing data were generated for the same controls.

Results

Using publicly available genomes, 4465 high-diversity genomic regions suited for targeted sequencing were identified, representing the P. falciparum heterozygome. For this study, 93 microhaplotypes with high diversity (median expected heterozygosity = 0.7) were selected along with 7 drug resistance loci. The sequencing method achieved very high coverage (median 99%), specificity (99.8%), and sensitivity (90% for haplotypes with 5% within sample frequency in dried blood spots with 100 parasites/µL). In silico analyses revealed that microhaplotypes provided much higher resolution to discriminate related from unrelated polyclonal infections than biallelic single-nucleotide polymorphism barcodes.

Conclusions

The bioinformatic and laboratory methods outlined here provide a flexible tool for efficient, low-cost, high-throughput interrogation of the P. falciparum genome, and can be tailored to simultaneously address multiple questions of interest in various epidemiological settings.

Confirmation of the absence of local transmission and geographic assignment of imported falciparum malaria cases to China using microsatellite panel

Background

Current methods to classify local and imported malaria infections depend primarily on patient travel history, which can have limited accuracy. Genotyping has been investigated as a complementary approach to track the spread of malaria and identify the origin of imported infections.

Methods

An extended panel of 26 microsatellites (16 new microsatellites) for Plasmodium falciparum was evaluated in 602 imported infections from 26 sub-Saharan African countries to the Jiangsu Province of People’s Republic of China. The potential of the 26 microsatellite markers to assign imported parasites to their geographic origin was assessed using a Bayesian method with Markov Chain Monte Carlo (MCMC) as implemented in the program Smoothed and Continuous Assignments (SCAT) with a modification to incorporate haploid genotype data.

Results

The newly designed microsatellites were polymorphic and are not in linkage disequilibrium with the existing microsatellites, supporting previous findings of high rate of recombination in sub-Saharan Africa. Consistent with epidemiology inferred from patients’ travel history, no evidence for local transmission was found; nearly all genetically related infections were identified in people who travelled to the same country near the same time. The smoothing assignment method assigned imported cases to their likely geographic origin with an accuracy (Angola: 59%; Nigeria: 51%; Equatorial Guinea: 40%) higher than would be achieved at random, reaching statistical significance for Angola and Equatorial Guinea.

Conclusions

Genotyping using an extended microsatellite panel is valuable for malaria case classification and programme evaluation in an elimination setting. A Bayesian method for assigning geographic origin of mammals based on genetic data was adapted for malaria and showed potential for identification of the origin of imported infections.

High Genetic Diversity of Plasmodium falciparum in the Low-Transmission Setting of the Kingdom of Eswatini

BACKGROUND

To better understand transmission dynamics, we characterized Plasmodium falciparum genetic diversity in Eswatini, where transmission is low and sustained by importation.

METHODS

Twenty-six P. falciparum microsatellites were genotyped in 66% of confirmed cases (2014-2016; N = 582). Population and within-host diversity were used to characterize differences between imported and locally acquired infections. Logistic regression was used to assess the added value of diversity metrics to classify imported and local infections beyond epidemiology data alone.

RESULTS

Parasite population in Eswatini was highly diverse (expected heterozygosity [HE] = 0.75) and complex: 67% polyclonal infections, mean multiplicity of infection (MOI) 2.2, and mean within-host infection fixation index (FWS) 0.84. Imported cases had comparable diversity to local cases but exhibited higher MOI (2.4 vs 2.0; P = .004) and lower mean FWS (0.82 vs 0.85; P = .03). Addition of MOI and FWS to multivariate analyses did not increase discrimination between imported and local infections.

CONCLUSIONS

In contrast to the common perception that P. falciparum diversity declines with decreasing transmission intensity, Eswatini isolates exhibited high parasite diversity consistent with high rates of malaria importation and limited local transmission. Estimates of malaria transmission intensity from genetic data need to consider the effect of importation, especially as countries near elimination.

Protective Immunity against Severe Malaria in Children Is Associated with a Limited Repertoire of Antibodies to Conserved PfEMP1 Variants

Extreme diversity of the major Plasmodium falciparum antigen, PfEMP1, poses a barrier to identifying targets of immunity to malaria. Here, we used protein microarrays containing hundreds of variants of the DBLα domain of PfEMP1 to cover the diversity of Papua New Guinean (PNG) parasites. Probing the plasma of a longitudinal cohort of malaria-exposed PNG children showed that group 2 DBLα antibodies were moderately associated with a lower risk of uncomplicated malaria, whereas individual variants were only weakly associated with clinical immunity. In contrast, antibodies to 85 individual group 1 and 2 DBLα variants were associated with a 70%-100% reduction in severe malaria. Of these, 17 variants were strong predictors of severe malaria. Analysis of full-length PfEMP1 sequences from PNG samples shows that these 17 variants are linked to pathogenic CIDR domains. This suggests that whereas immunity to uncomplicated malaria requires a broad repertoire of antibodies, immunity to severe malaria targets a subset of conserved variants. These findings provide insights into antimalarial immunity and potential antibody biomarkers for disease risk.

Applying next-generation sequencing to track falciparum malaria in sub-Saharan Africa

Next-generation sequencing (NGS) technologies are increasingly being used to address a diverse range of biological and epidemiological questions. The current understanding of malaria transmission dynamics and parasite movement mainly relies on the analyses of epidemiologic data, e.g. case counts and self-reported travel history data. However, travel history data are often not routinely collected or are incomplete, lacking the necessary level of accuracy. Although genetic data from routinely collected field samples provides an unprecedented opportunity to track the spread of malaria parasites, it remains an underutilized resource for surveillance due to lack of local awareness and capacity, limited access to sensitive laboratory methods and associated computational tools and difficulty in interpreting genetic epidemiology data. In this review, the potential roles of NGS in better understanding of transmission patterns, accurately tracking parasite movement and addressing the emerging challenges of imported malaria in low transmission settings of sub-Saharan Africa are discussed. Furthermore, this review highlights the insights gained from malaria genomic research and challenges associated with integrating malaria genomics into existing surveillance tools to inform control and elimination strategies.

Drug-Resistance and Population Structure of Plasmodium falciparum Across the Democratic Republic of Congo Using High-Throughput Molecular Inversion Probes

A better understanding of the drivers of the spread of malaria parasites and drug resistance across space and time is needed. These drivers can be elucidated using genetic tools. Here, a novel molecular inversion probe (MIP) panel targeting all major drug-resistance mutations and a set of microsatellites was used to genotype Plasmodium falciparum infections of 552 children from the 2013-2014 Demographic and Health Survey conducted in the Democratic Republic of the Congo (DRC). Microsatellite-based analysis of population structure suggests that parasites within the DRC form a homogeneous population. In contrast, sulfadoxine-resistance markers in dihydropteroate synthase show marked spatial structure with ongoing spread of double and triple mutants compared with 2007. These findings suggest that parasites in the DRC remain panmictic despite rapidly spreading antimalarial-resistance mutations. Moreover, highly multiplexed targeted sequencing using MIPs emerges as a cost-effective method for elucidating pathogen genetics in complex infections in large cohorts.

Declining trend of Plasmodium falciparum dihydrofolate reductase (dhfr) and dihydropteroate synthase (dhps) mutant alleles after the withdrawal of Sulfadoxine-Pyrimethamine in North Western Ethiopia

Antimalarial drug resistance is one of the major challenges in global efforts of malaria control and elimination. In 1998, chloroquine was abandoned and replaced with sulfadoxine/pyrimethamine, which in turn was replaced with artemether/lumefantrine for the treatment of uncomplicated falciparum malaria in 2004. Sulfadoxine/pyrimethamine resistance is associated with mutations in dihydrofolate reductase (Pfdhfr) and dihydropteroate synthase (Pfdhps) genes. The prevalence of mutation in Pfdhfr and Pfdhps genes were evaluated and compared for a total of 159 isolates collected in two different time points, 2005 and 2007/08, from Pawe hospital, in North Western Ethiopia. The frequency of triple Pfdhfr mutation decreased significantly from 50.8% (32/63) to 15.9% (10/63) (P<0.001), while Pfdhps double mutation remained high and changed only marginally from 69.2% (45/65) to 55.4% (40/65) (P = 0.08). The combined Pfdhfr/Pfdhps quintuple mutation, which is strongly associated with sulfadoxine/pyrimethamine resistance, was significantly decreased from 40.7% (24/59) to 13.6% (8/59) (P<0.0001). On the whole, significant decline in mutant alleles and re-emergence of wild type alleles were observed. The change in the frequency is explained by the reduction of residual drug-resistant parasites caused by the strong drug pressure imposed when sulfadoxine/pyrimethamine was the first-line drug, followed by lower fitness of these resistant parasites in the absence of drug pressure. Despite the decrease in the frequency of mutant alleles, higher percentages of mutation remain prevalent in the study area in 2007/08 in both Pfdhfr and Pfdhps genes. Therefore, further multi-centered studies in different parts of the country will be required to assess the re-emergence of sulfadoxine/pyrimethamine sensitive parasites and to monitor and prevent the establishment of multi drug resistant parasites in this region.

Phylogeography of var gene repertoires reveals fine-scale geospatial clustering of Plasmodium falciparum populations in a highly endemic area

Plasmodium falciparum malaria is a major global health problem that is being targeted for progressive elimination. Knowledge of local disease transmission patterns in endemic countries is critical to these elimination efforts. To investigate fine-scale patterns of malaria transmission, we have compared repertoires of rapidly evolving var genes in a highly endemic area. A total of 3680 high-quality DBLα-sequences were obtained from 68 P. falciparum isolates from ten villages spread over two distinct catchment areas on the north coast of Papua New Guinea (PNG). Modelling of the extent of var gene diversity in the two parasite populations predicts more than twice as many var gene alleles circulating within each catchment (Mugil = 906; Wosera = 1094) than previously recognized in PNG (Amele = 369). In addition, there were limited levels of var gene sharing between populations, consistent with local parasite population structure. Phylogeographic analyses demonstrate that while neutrally evolving microsatellite markers identified population structure only at the catchment level, var gene repertoires reveal further fine-scale geospatial clustering of parasite isolates. The clustering of parasite isolates by village in Mugil, but not in Wosera was consistent with the physical and cultural isolation of the human populations in the two catchments. The study highlights the microheterogeneity of P. falciparum transmission in highly endemic areas and demonstrates the potential of var genes as markers of local patterns of parasite population structure.