Deciphering the dynamics and evolution of insecticide resistance in malaria vectors is crucial for successful vector control. This study reports an increase of resistance intensity and a rise of multiple insecticide resistance in Anopheles funestus in Malawi leading to reduced bed net efficacy. Anopheles funestus group mosquitoes were collected in southern Malawi and the species composition, Plasmodium infection rate, susceptibility to insecticides and molecular bases of the resistance were analysed. Mosquito collection revealed a predominance of An. funestus group mosquitoes with a high hybrid rate (12.2 %) suggesting extensive species hybridization. An. funestus sensu stricto was the main Plasmodium vector (4.8 % infection). Consistently high levels of resistance to pyrethroid and carbamate insecticides were recorded and had increased between 2009 and 2014. Furthermore, the 2014 collection exhibited multiple insecticide resistance, notably to DDT, contrary to 2009. Increased pyrethroid resistance correlates with reduced efficacy of bed nets (<5 % mortality by Olyset® net), which can compromise control efforts. This change in resistance dynamics is mirrored by prevalent resistance mechanisms, firstly with increased over-expression of key pyrethroid resistance genes (CYP6Pa/b and CYP6M7) in 2014 and secondly, detection of the A296S-RDL dieldrin resistance mutation for the first time. However, the L119F-GSTe2 and kdr mutations were absent. Such increased resistance levels and rise of multiple resistance highlight the need to rapidly implement resistance management strategies to preserve the effectiveness of existing insecticide-based control interventions.

Insecticide resistance in malaria vectors threatens to reverse recent gains in malaria control. Deciphering patterns of gene flow and resistance evolution in malaria vectors is crucial to improve control strategies and preventing malaria resurgence. A genome-wide survey of Anopheles funestus genetic diversity Africa-wide revealed evidences of a major division between southern Africa and elsewhere, associated with different population histories. Three genomic regions exhibited strong signatures of selective sweep, each spanning major resistance loci (CYP6P9a/b, GSTe2 and CYP9K1). However, a sharp regional contrast was observed between populations correlating with gene flow barriers. Signatures of complex molecular evolution of resistance were detected with evidence of copy number variation, transposon insertion and a gene conversion between CYP6P9a/b paralog genes. Temporal analyses of samples before and after bed net scale up suggest that these genomic changes are driven by this control intervention. Multiple independent selective sweeps at the same locus in different parts of Africa suggests that local evolution of resistance in malaria vectors may be a greater threat than trans-regional spread of resistance haplotypes. Malaria control currently relies heavily on insecticide-based vector control interventions. Unfortunately, resistance to insecticides is threatening their continued effectiveness. Metabolic resistance has the greatest operational significance, yet it remains unclear how resistant mosquito populations, evolutionarily respond to the massive selection pressure from control interventions including insecticide-treated nets. Deciphering patterns of gene flow between populations of major malaria vectors such as Anopheles funestus and elucidating genomic signature of resistance evolution are crucial to design resistance management strategies and preventing malaria resurgence. Here, we performed a genome-wide survey of An. funestus genetic diversity from across its continental range using reduced-genome representation (ddRADseq) and whole genome (PoolSeq) approaches revealing evidences of significant barriers to gene flow impacting the spread of insecticide resistance alleles. This study detected signatures of strong selective sweeps occurring in genomic regions controlling cytochrome P450-based and glutathione s-transferase metabolic resistance to insecticides in this species. Fine-scale analysis of the major pyrethroid resistance-associated genomic regions revealed complex molecular evolution with evidence of copy number variation, transposon insertion and a gene conversion highlighting the risk that if this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised.

The phylogenetic relationships of representative species of the superfamily Lepocreadioidea were assessed using partial lsrDNA and nad1 sequences. Forty-two members of the family Lepocreadiidae, six putative members of the Enenteridae, six gyliauchenid species and one Gorgocephalidae, were studied along with 22 species representing 8 families. The Lepocreadioidea is found to be monophyletic, except for the two species of the putative enenterid genus Cadenatella, which are found to be only distantly related to the lepocreadioids. The Lepocreadioidea is formed of five clades in a polytomy, the Gorgocephalidae, a clade containing the Enenteridae and Gyliauchenidae, a small clade of atypical lepocreadiines and the deep-sea lepidapedine lepocreadiids, a small clade consisting of a freshwater form and a group of shallow-water putative lepidapedines and the final clade includes the remaining lepocreadiids. Thus, the generally accepted concept of the Lepocreadiidae is polyphyletic. The Enenteridae (minus Cadenatella) and the Gyliauchenidae are jointly and individually monophyletic, and are sister groups. The nad1 gene on its own places a deep-sea lepocreadiine with the deep-sea lepidapedines, whereas lsrDNA, combined sequences and morphology place this deep-sea lepocreadiine within a group of typical lepocreadiids. It could not be demonstrated that a significant proportion of sites in the nad1 gene evolved under positive selection; this anomalous relationship therefore remains unexplained. Most deep-sea species are in a monophyletic group, a few of which also occur in shallow waters, retaining some characters of the deep-sea clade. Many lepocreadioid species infect herbivorous fish, and it may be that the recently discovered life-cycle involving a bivalve first intermediate host and metacercariae encysted on vegetation is a common life-cycle pattern. The host relationships show no indication of co-speciation, although the host-spectrums exhibited are not random, with related worms tending to utilize related hosts. There are, however, many exceptions. Morphology is found to be of limited value in indicating higher level relationships. For example, even with the benefit of hindsight the gyliauchenids show little morphological similarity to their sister group, the Enenteridae.

Background: Several eukaryotic parasites form cysts that transmit infection. The process is found in diverse organisms such as Toxoplasma, Giardia, and nematodes. In Entamoeba histolytica this process cannot be induced in vitro, making it difficult to study. In Entamoeba invadens, stage conversion can be induced, but its utility as a model system to study developmental biology has been limited by a lack of genomic resources. We carried out genome and transcriptome sequencing of E. invadens to identify molecular processes involved in stage conversion. Results: We report the sequencing and assembly of the E. invadens genome and use whole transcriptome sequencing to characterize changes in gene expression during encystation and excystation. The E. invadens genome is larger than that of E. histolytica, apparently largely due to expansion of intergenic regions; overall gene number and the machinery for gene regulation are conserved between the species. Over half the genes are regulated during the switch between morphological forms and a key signaling molecule, phospholipase D, appears to regulate encystation. We provide evidence for the occurrence of meiosis during encystation, suggesting that stage conversion may play a key role in recombination between strains. Conclusions: Our analysis demonstrates that a number of core processes are common to encystation between distantly related parasites, including meiosis, lipid signaling and RNA modification. These data provide a foundation for understanding the developmental cascade in the important human pathogen E. histolytica and highlight conserved processes more widely relevant in enteric pathogens.

Insecticide resistance in mosquito populations threatens recent successes in malaria prevention. Elucidating patterns of genetic structure in malaria vectors to predict the speed and direction of the spread of resistance is essential to get ahead of the ‘resistance curve’ and to avert a public health catastrophe. Here, applying a combination of microsatellite analysis, whole genome sequencing and targeted sequencing of a resistance locus, we elucidated the continent-wide population structure of a major African malaria vector, Anopheles funestus. We identified a major selective sweep in a genomic region controlling cytochrome P450-based metabolic resistance conferring high resistance to pyrethroids. This selective sweep occurred since 2002, likely as a direct consequence of scaled up vector control as revealed by whole genome and fine-scale sequencing of pre- and post-intervention populations. Fine-scaled analysis of the pyrethroid resistance locus revealed that a resistance-associated allele of the cytochrome P450 monooxygenase CYP6P9a has swept through southern Africa to near fixation, in contrast to high polymorphism levels before interventions, conferring high levels of pyrethroid resistance linked to control failure. Population structure analysis revealed a barrier to gene flow between southern Africa and other areas, which may prevent or slow the spread of the southern mechanism of pyrethroid resistance to other regions. By identifying a genetic signature of pyrethroid-based interventions, we have demonstrated the intense selective pressure that control interventions exert on mosquito populations. If this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised. Malaria control currently relies heavily on insecticide-based vector control interventions. Unfortunately, resistance to insecticides threatens the continued effectiveness of these measures. Metabolic resistance, caused by increased detoxification of insecticides, presents the greatest threat to vector control, yet it remains unclear how these mechanisms are linked to underlying genetic changes driven by the massive selection pressure from these interventions, such as the widespread use of Long Lasting Insecticide Nets (LLINs) across Africa. Therefore, understanding the direction and speed at which this operationally important form of resistance spreads through mosquito populations is essential if we are to get ahead of the ‘resistance curve’ and avert a public health catastrophe. Here, using microsatellite markers, whole genome sequencing and fine-scale sequencing at a major resistance locus, we elucidated the Africa-wide population structure of Anopheles funestus, a major African malaria vector, and detected a strong selective sweep occurring in a genomic region controlling cytochrome P450-based metabolic pyrethroid resistance in this species. Furthermore, we demonstrated that this selective sweep is driven by the scale-up of insecticide-based malaria control in Africa, highlighting the risk that if this level of selection and spread of resistance continues unabated, our ability to control malaria with current interventions will be compromised.