Association of a rapidly selected 4.3kb transposon-containing structural variation with a P450-based resistance to pyrethroids in the African malaria vector Anopheles funestus.

Deciphering the evolutionary forces controlling insecticide resistance in malaria vectors remains a prerequisite to designing molecular tools to detect and assess resistance impact on control tools. Here, we demonstrate that a 4.3kb transposon-containing structural variation is associated with pyrethroid resistance in central/eastern African populations of the malaria vector Anopheles funestus. In this study, we analysed Pooled template sequencing data and direct sequencing to identify an insertion of 4.3kb containing a putative retro-transposon in the intergenic region of two P450s CYP6P5-CYP6P9b in mosquitoes of the malaria vector Anopheles funestus from Uganda. We then designed a PCR assay to track its spread temporally and regionally and decipher its role in insecticide resistance. The insertion originates in or near Uganda in East Africa, where it is fixed and has spread to high frequencies in the Central African nation of Cameroon but is still at low frequency in West Africa and absent in Southern Africa. A marked and rapid selection was observed with the 4.3kb-SV frequency increasing from 3% in 2014 to 98% in 2021 in Cameroon. A strong association was established between this SV and pyrethroid resistance in field populations and is reducing pyrethroid-only nets' efficacy. Genetic crosses and qRT-PCR revealed that this SV enhances the expression of CYP6P9a/b but not CYP6P5. Within this structural variant (SV), we identified putative binding sites for transcription factors associated with the regulation of detoxification genes. An inverse correlation was observed between the 4.3kb SV and malaria parasite infection, indicating that mosquitoes lacking the 4.3kb SV were more frequently infected compared to those possessing it. Our findings highlight the underexplored role and rapid spread of SVs in the evolution of insecticide resistance and provide additional tools for molecular surveillance of insecticide resistance.

Genetic variation present in the CYP3A4 gene in Ni-Vanuatu and Kenyan populations in malaria endemicity.

Cytochrome P450 3A4 (CYP3A4) enzyme is involved in the metabolism of about 30 % of clinically used drugs, including the antimalarials artemether and lumefantrine. CYP3A4 polymorphisms yield enzymatic variants that contribute to inter-individual variation in drug metabolism. Here, we examined CYP3A4 polymorphisms in populations from malaria-endemic islands in Lake Victoria, Kenya, and Vanuatu, to expand on the limited data sets. We used archived dried blood spots collected from 142 Kenyan and 263 ni-Vanuatu adults during cross-sectional malaria surveys in 2013 and 2005-13, respectively, to detect CYP3A4 variation by polymerase chain reaction (PCR) and sequencing. In Kenya, we identified 14 CYP3A4 single nucleotide polymorphisms (SNPs), including the 4713G (CYP3A4∗1B; allele frequency 83.9 %) and 19382A (CYP3A4∗15; 0.7 %) variants that were previously linked to altered metabolism of antimalarials. In Vanuatu, we detected 15 SNPs, including the 4713A (CYP3A4∗1A; 88.6 %) and 25183C (CYP3A4∗18; 0.6 %) variants. Additionally, we detected a rare and novel SNP C4614T (0.8 %) in the 5' untranslated region. A higher proportion of CYP3A4 genetic variance was found among ni-Vanuatu populations (16 %) than among Lake Victoria Kenyan populations (8 %). Our work augments the scarce data sets and contributes to improved precision medicine approaches, particularly to anti-malarial chemotherapy, in East African and Pacific Islander populations.

Transcriptomic analyses of differentially expressed human genes, micro RNAs and long-non-coding RNAs in severe, symptomatic and asymptomatic malaria infection.

Malaria transmission and endemicity in Africa remains hugely disproportionate compared to the rest of the world. The complex life cycle of P. falciparum (Pf) between the vertebrate human host and the anopheline vector results in differential expression of genes within and between hosts. An in-depth understanding of Pf interaction with various human genes through regulatory elements will pave way for identification of newer tools in the arsenal for malaria control. Therefore, the regulatory elements (REs) involved in the over- or under-expression of various host immune genes hold the key to elucidating alternative control measures that can be applied for disease surveillance, prompt diagnosis and treatment. We carried out an RNAseq analysis to identify differentially expressed genes and network elucidation of non-coding RNAs and target genes associated with immune response in individuals with different clinical outcomes. Raw RNAseq datasets, retrieved for analyses include individuals with severe (Gambia-20), symptomatic (Burkina Faso-15), asymptomatic (Mali-16) malaria as well as uninfected controls (Tanzania-20; Mali-36). Of the total 107 datasets retrieved, we identified 5534 differentially expressed genes (DEGs) among disease and control groups. A peculiar pattern of DEGs was observed, with individuals presenting with severe/symptomatic malaria having the highest and most diverse upregulated genes, while a reverse phenomenon was recorded among asymptomatic and uninfected individuals. In addition, we identified 141 differentially expressed micro RNA (miRNA), of which 78 and 63 were upregulated and downregulated respectively. Interactome analysis revealed a moderate interaction between DEGs and miRNAs. Of all identified miRNA, five were unique (hsa-mir-32, hsa-mir-25, hsa-mir-221, hsa-mir-29 and hsa-mir-148) because of their connectivity to several genes, including hsa-mir-221 connected to 16 genes. Six-hundred and eight differentially expressed long non coding RNA (lncRNA) were also identified, including SLC7A11, LINC01524 among the upregulated ones. Our study provides important insight into host immune genes undergoing differential expression under different malaria conditions. It also identified unique miRNAs and lncRNAs that modify and/or regulate the expression of various immune genes. These regulatory elements we surmise, have the potential to serve a diagnostic purpose in discriminating between individuals with severe/symptomatic malaria and those with asymptomatic infection or uninfected, following further clinical validation from field isolates.

Entire expressed peripheral blood transcriptome in pediatric severe malarial anemia.

This study on severe malarial anemia (SMA: Hb < 6.0 g/dL), a leading global cause of childhood morbidity and mortality, compares the entire expressed whole blood host transcriptome between Kenyan children (3-48 mos.) with non-SMA (Hb ≥ 6.0 g/dL, n = 39) and SMA (n = 18). Differential expression analyses reveal 1403 up-regulated and 279 down-regulated transcripts in SMA, signifying impairments in host inflammasome activation, cell death, and innate immune and cellular stress responses. Immune cell profiling shows decreased memory responses, antigen presentation, and immediate pathogen clearance, suggesting an immature/improperly regulated immune response in SMA. Module repertoire analysis of blood-specific gene signatures identifies up-regulation of erythroid genes, enhanced neutrophil activation, and impaired inflammatory responses in SMA. Enrichment analyses converge on disruptions in cellular homeostasis and regulatory pathways for the ubiquitin-proteasome system, autophagy, and heme metabolism. Pathway analyses highlight activation in response to hypoxic conditions [Hypoxia Inducible Factor (HIF)-1 target and Reactive Oxygen Species (ROS) signaling] as a central theme in SMA. These signaling pathways are also top-ranking in protein abundance measures and a Ugandan SMA cohort with available transcriptomic data. Targeted RNA-Seq validation shows strong concordance with our entire expressed transcriptome data. These findings identify key molecular themes in SMA pathogenesis, offering potential targets for new malaria therapies.

Y chromosome shredding in Anopheles gambiae: Insight into the cellular dynamics of a novel synthetic sex ratio distorter.

Despite efforts to explore the genome of the malaria vector Anopheles gambiae, the Y chromosome of this species remains enigmatic. The large number of repetitive and heterochromatic DNA sequences makes the Y chromosome exceptionally difficult to fully assemble, hampering the progress of gene editing techniques and functional studies for this chromosome. In this study, we made use of a bioinformatic platform to identify Y-specific repetitive DNA sequences that served as a target site for a CRISPR/Cas9 system. The activity of Cas9 in the reproductive organs of males caused damage to Y-bearing sperm without affecting their fertility, leading to a strong female bias in the progeny. Cytological investigation allowed us to identify meiotic defects and investigate sperm selection in this new synthetic sex ratio distorter system. In addition, alternative promoters enable us to target the Y chromosome in specific tissues and developmental stages of male mosquitoes, enabling studies that shed light on the role of this chromosome in male gametogenesis. This work paves the way for further insight into the poorly characterised Y chromosome of Anopheles gambiae. Moreover, the sex distorter strain we have generated promises to be a valuable tool for the advancement of studies in the field of developmental biology, with the potential to support the progress of genetic strategies aimed at controlling malaria mosquitoes and other pest species.

SNP-slice resolves mixed infections: simultaneously unveiling strain haplotypes and linking them to hosts.

MOTIVATION: Multi-strain infection is a common yet under-investigated phenomenon of many pathogens. Currently, biologists analyzing SNP information sometimes have to discard mixed infection samples as many downstream analyses require monogenomic inputs. Such a protocol impedes our understanding of the underlying genetic diversity, co-infection patterns, and genomic relatedness of pathogens. A scalable tool to learn and resolve the SNP-haplotypes from polygenomic data is an urgent need in molecular epidemiology. RESULTS: We develop a slice sampling Markov Chain Monte Carlo algorithm, named SNP-Slice, to learn not only the SNP-haplotypes of all strains in the populations but also which strains infect which hosts. Our method reconstructs SNP-haplotypes and individual heterozygosities accurately without reference panels and outperforms the state-of-the-art methods at estimating the multiplicity of infections and allele frequencies. Thus, SNP-Slice introduces a novel approach to address polygenomic data and opens a new avenue for resolving complex infection patterns in molecular surveillance. We illustrate the performance of SNP-Slice on empirical malaria and HIV datasets and provide recommendations for using our method on empirical datasets. AVAILABILITY AND IMPLEMENTATION: The implementation of the SNP-Slice algorithm, as well as scripts to analyze SNP-Slice outputs, are available at https://github.com/nianqiaoju/snp-slice.

Novel Ion Channel Genes in Malaria Parasites.

Ion channels serve many cellular functions including ion homeostasis, volume regulation, signaling, nutrient acquisition, and developmental progression. Although the complex life cycles of malaria parasites necessitate ion and solute flux across membranes, the whole-genome sequencing of the human pathogen Plasmodium falciparum revealed remarkably few orthologs of known ion channel genes. Contrasting with this, biochemical studies have implicated the channel-mediated flux of ions and nutritive solutes across several membranes in infected erythrocytes. Here, I review advances in the cellular and molecular biology of ion channels in malaria parasites. These studies have implicated novel parasite genes in the formation of at least two ion channels, with additional ion channels likely present in various membranes and parasite stages. Computational approaches that rely on homology to known channel genes from higher organisms will not be very helpful in identifying the molecular determinants of these activities. Given their unusual properties, novel molecular and structural features, and essential roles in pathogen survival and development, parasite channels should be promising targets for therapy development.

From Genome-wide Association Studies to Functional Variants: ARL14 Cis-regulatory Variants Are Associated With Severe Malaria.

BACKGROUND: Genome-wide association studies have identified several nonfunctional tag single-nucleotide polymorphisms (SNPs) associated with severe malaria. We hypothesized that causal SNPs could play a significant role in severe malaria by altering promoter or enhancer activity. Here, we sought to identify such regulatory SNPs. METHODS: SNPs in linkage disequilibrium with tagSNPs associated with severe malaria were identified and were further annotated using FUMA. Then, SNPs were prioritized using the integrative weighted scoring method to identify regulatory ones. Gene reporter assays were performed to assess the regulatory effect of a region containing candidates. The association between SNPs and severe malaria was assessed using logistic regression models in a Senegalese cohort. RESULTS: Among 418 SNPs, the best candidates were rs116525449 and rs79644959, which were in full disequilibrium between them, and located within the ARL14 promoter. Our gene reporter assay results revealed that the region containing the SNPs exhibited cell-specific promoter or enhancer activity, while the SNPs influenced promoter activity. We detected an association between severe malaria and those 2 SNPs using the overdominance model and we replicated the association of severe malaria with the tagSNP rs116423146. CONCLUSIONS: We suggest that these SNPs regulate ARL14 expression in immune cells and the presentation of antigens to T lymphocytes, thus influencing severe malaria development.

An elevated level of interleukin-17A in a Senegalese malaria cohort is associated with rs8193038 IL-17A genetic variant.

Malaria infection is a multifactorial disease partly modulated by host immuno-genetic factors. Recent evidence has demonstrated the importance of Interleukin-17 family proinflammatory cytokines and their genetic variants in host immunity. However, limited knowledge exists about their role in parasitic infections such as malaria. We aimed to investigate IL-17A serum levels in patients with severe and uncomplicated malaria and gene polymorphism's influence on the IL-17A serum levels. In this research, 125 severe (SM) and uncomplicated (UM) malaria patients and 48 free malaria controls were enrolled. IL-17A serum levels were measured with ELISA. PCR and DNA sequencing were used to assess host genetic polymorphisms in IL-17A. We performed a multivariate regression to estimate the impact of human IL-17A variants on IL-17A serum levels and malaria outcomes. Elevated serum IL-17A levels accompanied by increased parasitemia were found in SM patients compared to UM and controls (P < 0.0001). Also, the IL-17A levels were lower in SM patients who were deceased than in those who survived. In addition, the minor allele frequencies (MAF) of two IL-17A polymorphisms (rs3819024 and rs3748067) were more prevalent in SM patients than UM patients, indicating an essential role in SM. Interestingly, the heterozygous rs8193038 AG genotype was significantly associated with higher levels of IL-17A than the homozygous wild type (AA). According to our results, it can be concluded that the IL-17A gene rs8193038 polymorphism significantly affects IL-17A gene expression. Our results fill a gap in the implication of IL-17A gene polymorphisms on the cytokine level in a malaria cohort. IL-17A gene polymorphisms also may influence cytokine production in response to Plasmodium infections and may contribute to the hyperinflammatory responses during severe malaria outcomes.

On the feasibility of malaria hypothesis.

In 1954, Allison proposed that hemoglobin S (HbS) gene causes protection against fatal malaria. This would explain the high HbS gene frequency observed in certain regions hyperendemic for malaria, so-called "malaria hypothesis". This in silico study was conducted to examine the feasibility of the hypothesis under more realistic initial conditions, where a mutant gene with heterozygous advantage against malaria (e.g., HbS) was introduced in a group of Neolithic hunter-gatherers who decided to start agriculture nearby water where malaria killed a proportion of population. The tribe population size, number of children born to each woman in each generation, mortality from malaria and sickle cell disease, the protection factor provided by the gene carriers against malaria, the probability of mating between the members of the parent and offspring populations, population growth, and increased fertility in women heterozygous for HbS, were also considered. For effectively confer protection against malaria within the shortest possible period, the mutation needs to be happened in a small population. For a large population, the process would take around 100 generations (~ 2500 years) or more to provide an effective protection. Even then, the probability that the new gene could survive and propagate to future generations is about 35%. Conventional population genetics equations with differential or difference equations, give totally incorrect estimates of the gene frequency in small populations; discrete mathematics should be used, instead. After introduction of the advantageous mutation, the gene frequency increased until a steady state value. This value is far less than the gene frequency reported in certain tribes of Africa. It seems that the malaria hypothesis, per se, could not explain such a high observed gene frequency, unless HbS is associated with lower mortality from other causes too.

A genome-wide association study of neutrophil count in individuals associated to an African continental ancestry group facilitates studies of malaria pathogenesis.

BACKGROUND: 'Benign ethnic neutropenia' (BEN) is a heritable condition characterized by lower neutrophil counts, predominantly observed in individuals of African ancestry, and the genetic basis of BEN remains a subject of extensive research. In this study, we aimed to dissect the genetic architecture underlying neutrophil count variation through a linear-mixed model genome-wide association study (GWAS) in a population of African ancestry (N = 5976). Malaria caused by P. falciparum imposes a tremendous public health burden on people living in sub-Saharan Africa. Individuals living in malaria endemic regions often have a reduced circulating neutrophil count due to BEN, raising the possibility that reduced neutrophil counts modulate severity of malaria in susceptible populations. As a follow-up, we tested this hypothesis by conducting a Mendelian randomization (MR) analysis of neutrophil counts on severe malaria (MalariaGEN, N = 17,056). RESULTS: We carried out a GWAS of neutrophil count in individuals associated to an African continental ancestry group within UK Biobank, identifying 73 loci (r2 = 0.1) and 10 index SNPs (GCTA-COJO loci) associated with neutrophil count, including previously unknown rare loci regulating neutrophil count in a non-European population. BOLT-LMM was reliable when conducted in a non-European population, and additional covariates added to the model did not largely alter the results of the top loci or index SNPs. The two-sample bi-directional MR analysis between neutrophil count and severe malaria showed the greatest evidence for an effect between neutrophil count and severe anaemia, although the confidence intervals crossed the null. CONCLUSION: Our GWAS of neutrophil count revealed unique loci present in individuals of African ancestry. We note that a small sample-size reduced our power to identify variants with low allele frequencies and/or low effect sizes in our GWAS. Our work highlights the need for conducting large-scale biobank studies in Africa and for further exploring the link between neutrophils and severe malaria.

Polymorphisms in the human angiotensin converting enzyme gene (ACE) linked to susceptibility of COVID-19 and malaria infections in the Ghanaian population.

Genetic variations in the human angiotensin converting enzyme gene (ACE) influence ACE enzyme expression levels in humans and subsequently influence both communicable and non-communicable disease outcomes. More recently, polymorphisms in this gene have been linked to susceptibility and outcomes of infectious diseases such as the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and malaria infections. This study is the first to investigate the genetic diversity of ACE and ACE2 polymorphisms in the Ghanaian population. Archived filter blood blot samples from malaria patients aged ≤9 years were used. Molecular analysis for the detection of ACE rs4646994 (I/D), ACE2 rs2106809 (C/T) and rs2285666 (G/A) alleles as well as ACE2 exons 1-4 polymorphisms was conducted on 300 samples. The D allele (54%,162/300) was the most dominant polymorphism observed in the ACE rs4646994 gene whilst the G (68%, 204/300) and T alleles (59.3%,178/300) were the most frequent ACE2 rs2285666 and rs2106809 polymorphisms observed. For the 300 samples sequenced for ACE2 exons 1-4, analyses were done on 268, 282 and 137 quality sequences for exons 1, 2 and 3-4 respectively. For exon 1, the mutation D38N (2.2%; 6/268) was the most prevalent. The S19P and E37K mutations previously reported to influence COVID-19 infections were observed at low frequencies (0.4%, 1/268 each). No mutations were observed in exon 2. The N121K/T variants were the most seen in exons 3-4 at frequencies of 5.1% (K121, 7/137) and 2.9% (T121, 4/137) respectively. Most of the variants observed in the exons were novel compared to those reported in other populations in the world. This is the first study to investigate the genetic diversity of ACE and ACE2 genes in Ghanaians. The observation of novel mutations in the ACE2 gene is suggesting selection pressure. The importance of the mutations for communicable and non-communicable diseases (malaria and COVID-19) are further discussed.

Plasmodium falciparum transmission in the highlands of Ethiopia is driven by closely related and clonal parasites.

Malaria cases are frequently recorded in the Ethiopian highlands even at altitudes above 2000 m. The epidemiology of malaria in the Ethiopian highlands, and, in particular, the role of importation by human migration from the highly endemic lowlands is not well understood. We sequenced 187 Plasmodium falciparum samples from two sites in the Ethiopian highlands, Gondar (n = 159) and Ziway (n = 28), using a multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug resistance loci. Here, we characterize the parasite population structure and genetic relatedness. We identify moderate parasite diversity (mean HE : 0.54) and low infection complexity (74.9% monoclonal). A significant percentage of infections share microhaplotypes, even across transmission seasons and sites, indicating persistent local transmission. We identify multiple clusters of clonal or near-clonal infections, highlighting high genetic relatedness. Only 6.3% of individuals diagnosed with P. falciparum reported recent travel. Yet, in clonal or near-clonal clusters, infections of travellers were frequently observed first in time, suggesting that parasites may have been imported and then transmitted locally. 31.1% of infections are pfhrp2-deleted and 84.4% pfhrp3-deleted, and 28.7% have pfhrp2/3 double deletions. Parasites with pfhrp2/3 deletions and wild-type parasites are genetically distinct. Mutations associated with resistance to sulphadoxine-pyrimethamine or suggested to reduce sensitivity to lumefantrine are observed at near-fixation. In conclusion, genomic data corroborate local transmission and the importance of intensified control in the Ethiopian highlands.

Limited threat of Plasmodium falciparum pfhrp2 and pfhrp3 gene deletion to the utility of HRP2-based malaria RDTs in Northern Uganda.

BACKGROUND: Rapid diagnostic tests (RDTs) that detect Plasmodium falciparum histidine-rich protein-2 (PfHRP2) are exclusively deployed in Uganda, but deletion of the pfhrp2/3 target gene threatens their usefulness as malaria diagnosis and surveillance tools. METHODS: A cross-sectional survey was conducted at 40 sites across four regions of Uganda in Acholi, Lango, W. Nile and Karamoja from March 2021 to June 2023. Symptomatic malaria suspected patients were recruited and screened with both HRP2 and pan lactate dehydrogenase (pLDH) detecting RDTs. Dried blood spots (DBS) were collected from all patients and a random subset were used for genomic analysis to confirm parasite species and pfhrp2 and pfhrp3 gene status. Plasmodium species was determined using a conventional multiplex PCR while pfhrp2 and pfhrp3 gene deletions were determined using a real-time multiplex qPCR. Expression of the HRP2 protein antigen in a subset of samples was further assessed using a ELISA. RESULTS: Out of 2435 symptomatic patients tested for malaria, 1504 (61.8%) were positive on pLDH RDT. Overall, qPCR confirmed single pfhrp2 gene deletion in 1 out of 416 (0.2%) randomly selected samples that were confirmed of P. falciparum mono-infections. CONCLUSION: These findings show limited threat of pfhrp2/3 gene deletions in the survey areas suggesting that HRP2 RDTs are still useful diagnostic tools for surveillance and diagnosis of P. falciparum malaria infections in symptomatic patients in this setting. Periodic genomic surveillance is warranted to monitor the frequency and trend of gene deletions and its effect on RDTs.

Large DNA fragment knock-in and sequential gene editing in Plasmodium falciparum: a preliminary study using suicide-rescue-based CRISPR/Cas9 system.

CRISPR/Cas9 technology applied to Plasmodium falciparum offers the potential to greatly improve gene editing, but such expectations including large DNA fragment knock-ins and sequential gene editing have remained unfulfilled. Here, we achieved a major advance in addressing this challenge, especially for creating large DNA fragment knock-ins and sequential editing, by modifying our suicide-rescue-based system that has already been demonstrated to be highly efficient for conventional gene editing. This improved approach was confirmed to mediate efficient knock-ins of DNA fragments up to 6.3 kb, to produce "marker-free" genetically engineered parasites and to show potential for sequential gene editing. This represents an important advancement in establishing platforms for large-scale genome editing, which might gain a better understanding of gene function for the most lethal cause of malaria and contribute to adjusting synthetic biology strategies to live parasite malaria vaccine development. Site-directed knock-in of large DNA fragments is highly efficient using suicide-rescue-based CRISPR/Cas9 system, and sequential gene insertion is feasible but further confirmation is still needed.

Geographical emergence of sulfadoxine-pyrimethamine drug resistance-associated P. falciparum and P. malariae alleles in co-existing Anopheles mosquito and asymptomatic human populations across Cameroon.

Malaria molecular surveillance remains critical in detecting and tracking emerging parasite resistance to anti-malarial drugs. The current study employed molecular techniques to determine Plasmodium species prevalence and characterize the genetic diversity of Plasmodium falciparum and Plasmodium malariae molecular markers of sulfadoxine-pyrimethamine resistance in humans and wild Anopheles mosquito populations in Cameroon. Anopheles mosquito collections and parasitological survey were conducted in villages to determine Plasmodium species infection, and genomic phenotyping of anti-folate resistance was accomplished by sequencing the dihydrofolate-reductase (dhfr) and dihydropteroate-synthase (dhps) genes of naturally circulating P. falciparum and P. malariae isolates. The malaria prevalence in Elende was 73.5% with the 5-15 years age group harboring significant P. falciparum (27%) and P. falciparum + P. malariae (19%) infections. The polymorphism breadth of the pyrimethamine-associated Pfdhfr marker revealed a near fixation (94%) of the triple-mutant -A16I51R59N108I164. The Pfdhps backbone mediating sulfadoxine resistance reveals a high frequency of the V431A436G437K540A581A613 alleles (20.8%). Similarly, the Pmdhfr N50K55L57R58S59S114F168I170 haplotype (78.4%) was predominantly detected in the asexual blood stage. In contrast, the Pmdhps- S436A437occured at 37.2% frequency. The combined quadruple N50K55L57R58S59S114F168I170_ S436G437K540A581A613 (31.9%) was the major circulating haplotype with similar frequency in humans and mosquitoes. This study highlights the increasing frequency of the P. malariae parasite mostly common in asymptomatic individuals with apparent P. falciparum infection. Interventions directed at reducing malaria transmission such as the scaling-up of SP are favoring the emergence and spread of multiple drug-resistant alleles between the human and mosquito host systems.

Anopheles gambiae on remote islands in the Indian Ocean: origins and prospects for malaria elimination by genetic modification of extant populations.

The mosquito Anopheles gambiae s.s. is a primary malaria vector throughout sub-Saharan Africa including the islands of the Comoros archipelago (Anjouan, Grande Comore, Mayotte and Mohéli). These islands are located at the northern end of the Mozambique Channel in eastern Africa. Previous studies have shown a relatively high degree of genetic isolation between the Comoros islands and mainland populations of A. gambiae, but the origin of the island populations remains unclear. Here, we analyzed phylogenetic relationships among island and mainland populations using complete mitochondrial genome sequences of individual A. gambiae specimens. This work augments earlier studies based on analysis of the nuclear genome. We investigated the source population of A. gambiae for each island, estimated the number of introductions, when they occurred and explored evidence for contemporary gene flow between island and mainland populations. These studies are relevant to understanding historical patterns in the dispersal of this important malaria vector and provide information critical to assessing their potential for the exploration of genetic-based vector control methods to eliminate this disease. Phylogenetic analysis and haplotype networks were constructed from mitogenome sequences of 258 A. gambiae from the four islands. In addition, 112 individuals from seven countries across sub-Saharan Africa and Madagascar were included to identify potential source populations. Our results suggest that introduction events of A. gambiae into the Comoros archipelago were rare and recent events and support earlier claims that gene flow between the mainland and these islands is limited. This study is concordant with earlier work suggesting the suitability of these oceanic islands as appropriate sites for conducting field trial releases of genetically engineered mosquitoes (GEMs).

The first genome sequence of Anopheles squamous from Macha, Zambia.

Despite efforts to minimize the impacts of malaria and reduce the number of primary vectors, malaria has yet to be eliminated in Zambia. Understudied vector species may perpetuate malaria transmission in pre-elimination settings. Anopheles squamosus is one of the most abundantly caught mosquito species in southern Zambia and has previously been found with Plasmodium falciparum sporozoites, a causal agent of human malaria. This species may be a critical vector of malaria transmission, however, there is a lack of genetic information available for An. squamosus. We report the first genome data and the first complete mitogenome (Mt) sequence of An. squamosus. The sequence was extracted from one individual mosquito from the Chidakwa area in Macha, Zambia. The raw reads were obtained using Illumina Novaseq 6000 and assembled through NOVOplasty alignment with related species. The length of the An. squamosus Mt was 15,351 bp, with 77.9 % AT content. The closest match to the whole mitochondrial genome in the phylogenetic tree is the African malaria mosquito, Anopheles gambiae. Its genome data is available through National Center for Biotechnology Information (NCBI) Sequencing Reads Archive (SRA) with accession number SRR22114392. The mitochondrial genome was deposited in NCBI GenBank with the accession number OP776919. The ITS2 containing contig sequence was deposited in GenBank with the accession number OQ241725. Mitogenome annotation and a phylogenetic tree with related Anopheles mosquito species are provided.