Study site and population
Consented children aged 5 -15 years were screened following a cross-sectional study design from sub locations in Chulaimbo, Kisumu County (Figure 1), during the wet season from October to December 2019 and October to December 2020, and the dry season (January to March 2020). Chulaimbo is a rural site 19 km north-west of Kisumu City, located at 0.03572°S, 34.621°E, and an altitude range of 1328-1458 meters above sea level 21. The region has a mean annual temperature range of 120C - 35°C. This region experiences an average annual rainfall of 1352 mm and an average relative humidity between 66 and 83%. Malaria transmission in this area is endemic, with Plasmodium falciparum as the dominant parasite species in the area 22. Most residents are small-scale subsistence farmers.
Mosquitoes used for the study
Laboratory reared Anopheles gambiae female mosquitoes (Kisumu strain) between 3-5 days post-emergence were used for membrane feeding assays. This colony was selected and maintained at the Center of Excellence for Malaria Research in Homa Bay, Western Kenya. They were reared at temperatures of 27 - 290C, 69 - 80% relative humidity (RH), and a 12 hours light and 12 hours dark cycle. The colony was then constantly maintained on 10% sucrose 23 after the blood meal until the dissection day.
Identifying gametocyte carriers
Parasitological assessments to detect P. falciparum gametocyte carriers were conducted in school-aged children between the ages of 5 and 15 years, who had assented and had their guardians' consent to participate in the study. Blood samples were obtained from the children using finger pricks on well-labeled Whatman® 903 Protein Saver Cards (GE Healthcare WB100014) with the participants' information. A total of 50 µl of blood were collected and placed onto the cards, which were then air dried before being stored and stored at -20oC for further molecular analyses. As the blood was being collected on cards, thick and thin smears were also prepared for the same participant. The blood films were stained with 10% Giemsa and read after drying. Parasites were viewed under a compound microscope and Plasmodium species identified in thick smears Malaria parasites counts were read against 500 white blood cells. Gametocyte densities were determined in slides for all P. falciparum positive participants by counting the number of gametocytes per 500 leukocytes by microscopy and expressed as parasites per μl assuming a standard white blood cell (WBC) concentration of 8000/μl 24. Two trained microscopists took two readings per slide smear, and 20% of the slides were randomly selected for quality control verification by a senior external microscopist. Membrane feeding was limited to slides of gametocyte positive subjects only. Individuals who tested positive for malaria and had symptoms were referred to a local health center and treated according to Ministry of Health recommendations.25
Mosquito infections using membrane feeding assays
Individuals identified as carriers among screened volunteers, who tested positive for P. falciparum gametocytes donated blood for whole blood and serum replaced experiments in the laboratory. Blood was drawn intravenously by a professional phlebotomist using butterfly needles. Approximately 3ml of blood was collected by venipuncture in heparinized tubes for each volunteer. An aliquot of 1ml of blood was immediately placed into pre-warmed hemotek feeders (1ml capacity) at 37°C, while another 1.5ml was transferred into 2ml Eppendorf tubes and centrifuged at 2000 rounds per minute for 2 minutes before adding it to the hemotek feeders. The supernatant of serum was discarded and replaced with a naïve human serum type AB (Bio Whittaker, Cambrex Bio Science Walkersville, MD, USA). Replaced blood was then quickly transferred to the feeders to allow the starved mosquitoes to feed.
Participants’ blood positive for gametocytes was used to feed the insectary reared Anopheles gambiae mosquitoes using membrane feeding assays (serum replacement and whole blood experiments) 26. Aggressive 3-5 day old female Anopheles gambiae mosquitoes were starved for 6 to 8 hours prior to feeding on infected blood. Whole blood and serum replacement experiments were conducted for each participant. A total of 37 paired experiments were conducted. Each feeding cup contained 100 mosquitoes. The mosquitoes were allowed to feed from different feeders of the same infected blood for 15-30 minutes through a parafilm membrane. All membrane feeding procedures were conducted at 370C using the hemotek system. Only fully engorged blood-fed mosquitoes were selected and maintained at 27-290C temperatures and 69-80% relative humidity following a 12 hours light and 12 hours dark cycle. They were given 10% sucrose for 9 days post-feeding and the ones that survived were dissected for midgut oocysts enumeration. The unfed and partially fed mosquitoes were discarded by freezing them for 15 minutes at -200C. After membrane feeding, volunteers were treated with artemether-lumefantrine (Coartem®) according to the Ministry of Health guideline 25.
Oocysts counts
All fully engorged mosquitoes that survived on day 8 or 9 post-feeding were dissected under a dissecting microscope as described by Afrane et al 27. Briefly, each mosquito gut was carefully pulled out from the abdomen in 0.5% mercurochrome and allowed to stain for 10 minutes. The midguts were then examined for the presence of oocysts under a light microscope. The number of oocysts observed were counted and recorded per mosquito gut. The oocysts load was expressed as the number of oocysts per infected mosquito. Mosquito carcasses corresponding to their infected midguts were labeled and preserved for further molecular assays to determine TEP1 genotypes28. Briefly, TEP1 was genotyped using a nested PCR-RFLP targeting 892 base pairs for nest 1 and a final fragment length of 758 base pairs after nest 2. Both PCR reaction conditions were set as denaturation at 95 °C for 3 min, 35 cycles of 94°C for 30s, annealing at 55°C for 30s, extension at 72°C for 30s, and a final step at 72 °C for 6 min using Dream Taq Green Master Mix (Thermo Fisher Scientific). PCR products were then digested using restriction enzymes Bam HI, Hind III, or Bse NI (New England Biolabs Inc) according to the manufacturer’s instructions and the result analyzed on 2.5% agarose gel electrophoresis. The TEP1 allelic classes were determined by fragment size of restriction enzyme digestion. A subset was also randomly selected for the genotype confirmatory purposes by sequencing.
DNA extraction and parasite genotyping
The Chelex technique was used to obtain Plasmodium parasite DNA from the dried blood spots 29. As previously reported 30, a multiplex real-time PCR (RT-PCR) was utilized to identify Plasmodium species. Pfs47 was genotyped using PCR and Sanger sequencing, as previously published 30.
Ethical approval
The ethical review board of the Maseno University, Kenya (MSU/DRPI/MUERC/00456/17) reviewed and approved the protocol for screening of P. falciparum gametocyte carriers and subsequent intravenous blood drawing. A detailed written informed assent and consent to participate in the study was provided by all study volunteers and their parents or guardians. Feeding of mosquitoes was conducted in a secure, insect-proof room at the Chulaimbo health center. All experiments and methods were performed in accordance with the institution’s guidelines and regulations.
Statistical analysis
Data from the participants was tabulated in Microsoft Excel V16. Computing descriptive statistics (sum, mean, standard deviation, standard error, and 95% confidence interval) and comparing means were done using Graph Pad Prism v.8.0.1 and SPSS version 25 for Windows software. The Shapiro–Wilk normality test was used to check data normality before performing multiple mean comparisons and chi-square tests. Data were considered statistically significant at P< 0.05. The codoncode Aligner 11.0.1 (CodonCode Corp., Centerville, MA) was used to check the sequence quality and trim low-quality bases. Bio-Edit software was used to align the sequences and determine the nonsynonymous mutations and codon changes based on reference sequence (Pf3D7_1346800). MEGA software was used to construct the UPGMA (unweighted pair group method with arithmetic mean) tree based on the Kimura 2-parameter (K2P) distance model with 1,000 bootstrap replicates.