Study design and setting
Kenya is characterised by a diverse malaria ecology, ranging from semi-arid, sporadic transmission, epidemic prone highland malaria, coastal intermediary intensity transmission and high transmission areas surrounding Lake Victoria (Alegana et al., 2021). This study leveraged previous work that surveyed the prevalence and distribution of malaria infection among school children in Kenya from March to November 2022 across 45 of the 47 counties in Kenya (Alegana et al., 2021). For the current antimalarial drug-resistance surveillance study, we employed an opportunistic sampling strategy alongside this previous work whereby we included only 15 counties across malaria-endemic regions: Bungoma, Busia, Homa Bay, Kakamega, Kisumu, Migori, Siaya, and Vihiga, (lakeside), Kisii (highland), Kilifi and Kwale (Coastal), Kirinyaga (low risk), Turkana, Tana River and West Pokot (semi-arid). 121 schools were included in the current study as shown in Table 1. Dried blood spots (DBS) were collected from all participants from the 15 counties, but only the rapid diagnostic test (RDT) positive samples were processed for sequencing.
Table 1. All samples collected between March 2022 and November 2022 used for plasmodium parasite genotyping.
Study Site
|
Time point
|
Schools sampled
|
Samples collected
|
mRDT No. [%]
|
DNA extraction
|
18s Pf pos,
No. [%]
|
No. of PCR amplicons
|
Bungoma
|
Mar/Apr
|
10
|
987
|
74 [8]
|
74
|
54[73]
|
54
|
Busia
|
Mar/Apr
|
10
|
989
|
322 [33]
|
322
|
276[85]
|
181
|
Homabay
|
Mar/Apr
|
10
|
979
|
174 [18]
|
174
|
86[49]
|
55
|
Kakamega
|
Mar/Apr
|
12
|
1195
|
145 [12]
|
145
|
106[73]
|
69
|
Siaya
|
Mar/Apr
|
10
|
993
|
337 [34]
|
337
|
277[82]
|
201
|
Vihiga
|
Mar/Apr
|
10
|
996
|
68 [7]
|
68
|
51[75]
|
51
|
Kisumu
|
Mar/Apr
|
10
|
966
|
161 [17]
|
161
|
126[78]
|
92
|
Migori
|
Mar/Apr
|
10
|
981
|
275 [28]
|
275
|
206[75]
|
137
|
Kisii
|
May/Jun
|
4
|
182
|
12 [7]
|
12
|
12[100]
|
10
|
Kwale
|
May/Jun
|
6
|
373
|
18 [5]
|
18
|
8[44]
|
8
|
Kilifi
|
May/Jun
|
8
|
541
|
1 [0.2]
|
1
|
1[100]
|
0
|
Kirinyaga
|
May/Jun
|
2
|
135
|
0
|
0
|
0
|
0
|
Tana River
|
May/Jun
|
4
|
246
|
0
|
0
|
0
|
0
|
Turkana
|
Oct/Nov
|
11
|
1024
|
92 [9]
|
92
|
42[46]
|
42
|
West Pokot
|
Oct/Nov
|
4
|
286
|
32 [11]
|
32
|
20[63]
|
20
|
Total
|
|
121
|
10873
|
1711
|
1711
|
1260
|
920
|
School-based surveys were undertaken using sampling methods described elsewhere (Osoti et al., 2022). In brief, schools were randomly selected from each county proportional to schools within each county and 100 children aged 5-14 years randomly selected from the school register on the survey day. An initial meeting was held with the head teacher, school committee, and parents to explain the study, and before enrolment, the student’s parents or guardians provided informed consent while each child assented. Instead of written opt-in consent, parental consent was based on passive, opt-out consent (Ellickson & Hawes, 1989). A finger prick was taken to obtain blood for malaria rapid diagnostic test (mRDT) (CareStartTM), and a Dried Blood Spot Sample (DBS) was also collected on a Whatman CF12 filter paper (Cat No. 10535097, Cytiva, USA) for molecular diagnosis. After allowing the DBS samples to air dry for at least 1 hour, they were individually packed in zip-lock bags with a desiccant and shipped to the KEMRI-Wellcome Trust Research Programme laboratories. Children who tested mRDT positive for malaria were treated following the national malaria treatment guidelines with artemether-lumefantrine (AL) as a co-formulated tablet of 20mg of artemether and 120mg of lumefantrine. The study was approved by the KEMRI and National Ethics Review Committee (number KEMRI/SERU/ESACIPAC/11/3822). Additional approval was provided by the county health and education authorities.
DNA extraction, and amplicon deep sequencing (AmpSeq)
To extract parasite DNA from the RDT-positive DBS a previously published protocol (Osoti et al., 2022) was used. Briefly, 4 punches 6mm each were punched from two locations (at the center and periphery) of the DBS and placed into a 1.5ml Eppendorf tube using sterile tweezers. DNA extraction was done using the Chelex saponin method (Baidjoe et al., 2013). Parasite DNA was amplified using 18S rRNA Plasmodium falciparum qPCR assay (Hermsen et al., 2001). Samples from Kakamega, Busia, Homa Bay, Kisumu, Migori, and Siaya with a median cycling threshold (Ct) above 31.5 were excluded from sequencing due to the large sample size. In contrast, for samples from Bungoma, Kwale, Kisii, Vihiga, Turkana, and West Pokot, a minimum Ct threshold of 38.5 was set, as the sample sizes from these regions were smaller.
Using a previously published nested PCR approach (Osoti et al., 2022) amplicons were generated for the following genes, using molecular identifiers (MID) labelled primers (Osoti et al., 2022)for P. falciparum dihydrofolate reductase (Pfdhfr: PF3D7_0417200), P. falciparum dihydropteroate synthase (Pfdhps: PF3D7_0810800), P. falciparum kelch13 (Pfk13: PF3D7_1347700) and P. falciparum multidrug resistance 1 (Pfmdr1: PF3D7_0523000). Slight modifications were made to the initial protocol, generating PCR amplicons, and sequencing them in singlets rather than duplicates. The Pfk13 gene was split and amplified in two separate fragments. One fragment covered codon 469, while the other covered codon 675. The amplification used external PCR primers Pfk13ext675F (5’-GAAGCCTTGTTGAAAGAAGC-3’) and Pfk13ext675R (5’-CGGAGTGACCAAATCTGG-3’), and internal PCR primers Pfk13int675F (5’-GGGGGATATGATGGCTCTTCT-3’) and Pfk13int675R (5’-ACTAATAAAGATGGGCCAAGC-3’).
The PCR products were individually purified using the AMPure XP beads (Beckman Coulter, Inc.) as per the manufacturer’s instructions. Thereafter, the purified DNA was quantified using a Qubit double-strand DNA (dsDNA) high-sensitivity (HS) assay kit, according to the manufacturer’s instructions. Library preparation was done using the KAPA kit, while a size selection clean-up was done using 0.8X AMPure XP beads. The adapter-ligated libraries were amplified using Illumina primers and cleaned with 0.8X AMPure. The libraries were quantified using a Qubit dsDNA HS kit and sizes were verified by the DNA 1000 assay kit using the 2100 Bioanalyzer (Agilent). The libraries were mixed in equimolar concentrations, denatured, spiked with 8% PhiX DNA, and finally AmpSeq was done using a MiSeq reagent kit v3 (Illumina). The entire procedure has been reported by Osoti et al (2022).
Sequence data analysis.
SeekDeep v3.0.1 initially demultiplexed the sequences based on the MIDs. The paired consensus reads were trimmed and clustered to estimate the frequency of DNA clusters (referred to as microhaplotypes from here henceforth). Microhaplotypes were discarded if their relative frequency was <5%. Chimeric reads were considered PCR artefacts and discarded. The SNP and microhaplotype frequencies in the population were calculated as the number of samples containing the SNP or microhaplotype over the total number of samples genotyped to determine the number of individuals harbouring a SNP or haplotype. All statistical analysis were performed in R v4.0.3 (R Core Team, 2022).