31
views
0
recommends
+1 Recommend
0 collections
    0
    shares
      • Record: found
      • Abstract: found
      • Article: found
      Is Open Access

      Promising approach to reducing Malaria transmission by ivermectin: Sporontocidal effect against Plasmodium vivax in the South American vectors Anopheles aquasalis and Anopheles darlingi

      research-article

      Read this article at

      Bookmark
          There is no author summary for this article yet. Authors can add summaries to their articles on ScienceOpen to make them more accessible to a non-specialist audience.

          Abstract

          Background

          The mosquito resistance to the insecticides threatens malaria control efforts, potentially becoming a major public health issue. Alternative methods like ivermectin (IVM) administration to humans has been suggested as a possible vector control to reduce Plasmodium transmission. Anopheles aquasalis and Anopheles darlingi are competent vectors for Plasmodium vivax, and they have been responsible for various malaria outbreaks in the coast of Brazil and the Amazon Region of South America.

          Methods

          To determine the IVM susceptibility against P. vivax in An. aquasalis and An. darlingi, ivermectin were mixed in P. vivax infected blood: ( 1) Powdered IVM at four concentrations (0, 5, 10, 20 or 40 ng/mL). ( 2) Plasma (0 hours, 4 hours, 1 day, 5, 10 and 14 days) was collected from healthy volunteers after to administer a single oral dose of IVM (200 μg/kg) ( 3) Mosquitoes infected with P. vivax and after 4 days was provided with IVM plasma collected 4 hours post-treatment ( 4) P. vivax-infected patients were treated with various combinations of IVM, chloroquine, and primaquine and plasma or whole blood was collected at 4 hours. Seven days after the infective blood meal, mosquitoes were dissected to evaluate oocyst presence. Additionally, the ex vivo effects of IVM against asexual blood-stage P. vivax was evaluated.

          Results

          IVM significantly reduced the prevalence of An. aquasalis that developed oocysts in 10 to 40 ng/mL pIVM concentrations and plasma 4 hours, 1 day and 5 days. In An. darlingi to 4 hours and 1 day. The An. aquasalis mortality was expressively increased in pIVM (40ng/mL) and plasma 4 hours, 1, 5 10 and 14 days post-intake drug and in An. darlingi only to 4 hours and 1 day. The double fed meal with mIVM by the mosquitoes has a considerable impact on the proportion of infected mosquitoes for 7 days post-feeding. The oocyst infection prevalence and intensity were notably reduced when mosquitoes ingested blood from P. vivax patients that ingested IVM+CQ, PQ+CQ and IVM+PQ+CQ. P. vivax asexual development was considerably inhibited by mIVM at four-fold dilutions.

          Conclusion

          In conclusion, whole blood spiked with IVM reduced the infection rate of P. vivax in An. aquasalis and An. darlingi, and increased the mortality of mosquitoes. Plasma from healthy volunteers after IVM administration affect asexual P. vivax development. These findings support that ivermectin may be used to decrease P. vivax transmission.

          Author summary

          Malaria is one of the most important infectious diseases in the world with hundreds of millions of new cases every year. The disease is caused by parasites of the genus Plasmodium where Plasmodium vivax represent most of the cases in the Americas. Current strategies to combat malaria transmission are being implemented; however, widespread insecticide resistance in vectors threatens the effectiveness of vector control programs. Ivermectin (IVM) has arisen as a new potential tool to be added to these programs as it has mosquito-lethal and sporontocidal properties making it a promising transmission reduction drug. Plasmodium vivax was drawn from patients, mixed with powdered IVM and metabolized IVM in plasma collected from healthy volunteers receiving IVM, and fed to mosquitoes via membrane feeding. Powdered and metabolized IVM interrupt P. vivax transmission, reducing oocyst infection and intensity rate of two South American malaria vectors An. aquasalis and An. darlingi. We also demonstrate the effect of IVM on asexual stages development of P. vivax, providing evidence that IVM may affect different parasite life cycle stages. Our findings place IVM as a strong candidate for malaria transmission reducing interventions.

          Related collections

          Most cited references61

          • Record: found
          • Abstract: found
          • Article: not found

          Natural microbe-mediated refractoriness to Plasmodium infection in Anopheles gambiae.

          Malaria parasite transmission depends on the successful transition of Plasmodium through discrete developmental stages in the lumen of the mosquito midgut. Like the human intestinal tract, the mosquito midgut contains a diverse microbial flora, which may compromise the ability of Plasmodium to establish infection. We have identified an Enterobacter bacterium isolated from wild mosquito populations in Zambia that renders the mosquito resistant to infection with the human malaria parasite Plasmodium falciparum by interfering with parasite development before invasion of the midgut epithelium. Phenotypic analyses showed that the anti-Plasmodium mechanism requires small populations of replicating bacteria and is mediated through a mosquito-independent interaction with the malaria parasite. We show that this anti-Plasmodium effect is largely caused by bacterial generation of reactive oxygen species.
            Bookmark
            • Record: found
            • Abstract: found
            • Article: found
            Is Open Access

            Midgut Microbiota of the Malaria Mosquito Vector Anopheles gambiae and Interactions with Plasmodium falciparum Infection

            Introduction Understanding how Plasmodium-Anopheles interactions contribute to the mosquito vector competence has received great attention lately, and the increasing knowledge promises to contribute to the development of new malaria control strategies [1]–[3]. Malaria still remains a serious health problem in developing African countries, causing more than 1 million deaths annually. Almost all these deaths are caused by the parasite Plasmodium falciparum whose major vector in Africa is Anopheles gambiae, which is widely distributed throughout the afro-tropical belt. A. gambiae s.s. is divided into two morphologically indistinguishable molecular forms, known as M and S, which are regarded as incipient species [4]–[6]. The M and S molecular forms exhibit ecological preferences [7], [8], but their respective epidemiological importance in malaria transmission has been poorly documented so far [9], [10]. The susceptibility of Anopheles mosquitoes to Plasmodium infection is under genetic control [11]–[13], but the large variability in oocyst number among closely related mosquitoes indicates that environmental factors also play a role. Multiple lines of evidence suggest that mosquito bacterial communities influence vector competence [14]–[17]. A protective role of Anopheles midgut bacteria against malaria infections was demonstrated by using antibiotic treatment to clear the gut microbiota, which resulted in enhanced Plasmodium infections [15], [18]. Consistently, coinfections of bacteria with Plasmodium reduced the number of developing oocysts in the mosquito midgut, both in laboratory and field conditions [15], [19], [20]–[24]. Interestingly, Cirimotich et al. [24] recently described an Enterobacter bacterium isolated from wild mosquitoes in Zambia that confers refractoriness to P. falciparum infection. Mechanisms mediating this refractory phenotype remain elusive. Instead of eliciting immune responses leading to reduced levels of parasite burden, the experiments conducted by Cirimotich et al. [24] revealed that the inhibition of Plasmodium development by commensal microbiota occurs through production of reactive oxygen species by the Enterobacter bacteria that directly target Plasmodium parasites in the midgut [17], [24]. Bacterial diversity in the Anopheles species is thought to be particularly low at the adult stage because of gut renewal during metamorphosis from pupae to adults. Nevertheless several bacterial species have been identified in the adult mosquito midgut using different conventional culture-mediated techniques [16], [23], [24]. These bacteria were acquired from the aquatic environment during immature stage development [25], [26], although vertical transmission (from mother to offspring) also has been documented [14], [16], [25], [27]. Generally, knowledge on mosquito midgut bacterial communities remains largely unknown, mostly because of the limitations of isolating techniques based on culturing and to the low resolution of fingerprint analysis. However, the recent deployment of next generation DNA sequencing technologies has provided new opportunities to explore microbial diversity of complex environments [28]–[31] as well as to further investigate disease susceptibilities and host-bacteria-pathogen interactions [32]–[34]. In this study, we performed a meta-taxogenomic analysis of microbial communities in the midguts of adult mosquitoes originating from natural larval habitats in Cameroon. We further investigated correlations between midgut microbiota and the mosquito malaria infection status. Previous investigations of bacteria-Plasmodium interactions in the mosquito vector have considered laboratory-reared mosquitoes challenged with cultured bacteria and infected with a cultured P. falciparum line. Here, we challenged wild female mosquitoes with a natural isolate of P. falciparum, thereby offering an opportunity to examine natural bacteria/parasite associations. Our analysis revealed that the midgut bacterial diversity represents an important force shaping the mosquito vector competence, where bacteria of the Enterobacteriaceae genera benefit P. falciparum development. Results Mosquito susceptibility to P. falciparum and bacterial flora A total of 92 A. gambiae mosquitoes collected at the larval stage and reared to adults in the insectary were successfully fed through membrane feeders on gametocyte containing blood from a single individual. The origin of mosquitoes and their genetic characteristics (molecular form and infection status) are summarized in Table S1. Mosquitoes of M molecular form were significantly more infected than those of S form (48.5% versus 27.1%; OR 0.39; 95% CI: 0.15–1.06; P = 0.044). However, the comparison of infection prevalence (number of infected mosquitoes) between the different localities revealed a “sampling site effect” (Fisher's exact test P 2% of the overall), the S1 domain reaches a lower percentage, indicating that this 16S library capable of identifying a greater number of minor clades (Figure 1). In addition, the S1 domain had better resolution, allowing more precise assigning of sequence tags (see Figure 1, mosquito NKD97). We then performed further analyses on the S1 domain, for all 30 samples that were successful for pyrosequencing. 10.1371/journal.ppat.1002742.g001 Figure 1 Comparison of bacterial diversity for the three 16S libraries at the genus level. Tag abundance was compared for three 16S libraries, and the graph shows the bacterial flora of six mosquito midguts. The three 16S libraries were obtained using primer sets targeting different 16S domains, as described in the Materials and Methods section. Only the most abundant categories (>2%) were considered. The S1 library only reached 95%, showing this domain allowed identifications for a greater number of minor clades. For mosquito NKD97, S2 and S3 primer sets only allowed the identification at the Enterobacteriaceae family level, whereas S1 reached the assignment at the genus level, Serratia. The bacterial communities of the mosquito midgut belonged to 26 different phyla, among which, 5 represented more than 99% of the total microbiota: Proteobacteria, Bacteroidetes, Actinobacteria, Firmicutes, and Fusobacteria. We examined the relative abundance of the major classes, that is, detected in more than 30% of the samples and having an average abundance of >0.1% (Figure 2). The gut microbiota presented a large inter-individual variability and was dominated by few taxa. The first striking result came from laboratory-reared mosquitoes that exhibited a drastically different composition of midgut bacteria from field mosquitoes. More than 96% of tags corresponding to bacteria in the midguts of mosquitoes from the Ngousso colony were assigned to Flavobacteria, although this class accounted for only 0.38% (±0.24) of the sequence tags in field mosquitoes. Similarity searches against the SSU SILVA database (release 108) identified Flavobacteria tags as belonging to Elizabethkingia sp., which already has been isolated from A. gambiae midguts from different insectaries [15], [35], [36]. The guts of Ngousso mosquitoes also contained, to a lower extent, Gammaproteobacteria (Pseudomonas sp.) and Alphaproteobacteria (Asaia sp.). 10.1371/journal.ppat.1002742.g002 Figure 2 Relative abundance of the different bacterial classes within each mosquito midgut sample. Mvan and NKD (Nkolondon) indicate the geographical origin of mosquitoes, NG being mosquitoes of the laboratory colony Ngousso. Pf+ and Pf− designate the P. falciparum infection status of the challenged mosquitoes, positive and negative, respectively. Only class rank groups that represented >1% of the total reads, and identified in at least 30% of mosquitoes, are shown. “Unclass” represents tags that could not be assigned to the class level, and were grouped into a higher taxonomical rank. In mosquitoes from natural habitats, midguts were mainly colonized by Proteobacteria (94%), and the most prominent classes were Gammaproteobacteria, Alphaproteobacteria and Betaproteobacteria (Figure 2). Figure 1 clearly shows the difference in bacterial composition between mosquitoes from the 2 breeding sites used in this study. In mosquitoes originating from Nkolondom, the intestinal bacterial flora is dominated by Alphaproteobacteria (68.65±7.38%), mainly of the genus Asaia sp. By contrast, the three major classes are almost equally represented in the midgut of mosquitoes from Mvan, although with large individual variability. The main taxa in this locality were Asaia, Sphingomonas, Burkholderia, Ralstonia, and Enterobacteriaceae. Enterobacteriaceae could not be assigned to more precise taxonomic ranks. Midgut bacteria were unevenly distributed among individual mosquitoes and between the different sampled localities. Several genera were found in all, or at least in a large majority, of the mosquitoes possibly representing the “mosquito midgut core microbiota” (Table S3). They included members of the genera Asaia, Burkholderia, Serratia, Ralstonia, Acinetobacter, Pseudomonas, Sphingomonas, Staphylococcus, Streptococcus, and Escherichia/Shigella. Of note, Asaia sp. was found in all samples, and its relative abundance showed great variation from one midgut to another, ranging from 1.49 to 98.95% in mosquitoes from Nkolondom and from 0.04 to 49.66% in those from Mvan. The group of unassigned Enterobacteriaceae also was identified in all field samples, with relative abundance ranging from 0.01 to 1.03% and from 0.04 to 71.51% in Nkolondom and Mvan, respectively. Other specific members of Enterobacteriaceae were frequent and represent a large proportion of the midgut bacterial communities: Serratia spp. accounted for 96.93% of the midgut bacteria in a mosquito (NKD97) from Nkolondom, and Cedecea spp. encompassed 12% of the gut bacterial content in 2 mosquitoes from Mvan. Escherichia/Shigella was found in more than 85% of the mosquitoes, at low densities. The sequence of the Esp_Z Enterobacter (JF690924), despite its presence in our reference database, was absent from the analyzed reads. Of note, the midgut bacterial flora was mainly composed of Gram-negative communities. No Gram-positive bacterium was identified in the laboratory mosquitoes, whereas they represented 5% of the total microbiota of the field mosquitoes. Gram-positive bacteria belonged to the classes Bacilli and Actinobacteria. To determine whether all phylotypes present in the mosquito midgut microbiota were detected in this study, we performed a rarefaction analysis for each sample on tags from the S1 domain; rarefaction curves are shown in Figure S1. The rarefaction curves decrease rapidly at approximately 2,000 sequences per sample and reach saturation at 3,000, indicating that our sequencing effort was sufficient to catch the overall bacterial diversity in the mosquito midgut. The rarefaction curves show the large variability in bacterial complexity among samples, varying from 13 to 340 operational taxonomic units (OTUs). In addition, they illustrate the paucity of clusters in the midgut of laboratory mosquitoes. They also revealed greater bacterial diversity in the samples from Mvan compared with those from Nkolondom (185±51 and 110±30, respectively; t-test t = 2.385, P = 0.025). Microbial diversity in the mosquito intestinal microbiota The Chao1 estimator, which gauges the number of unseen “species,” predicted that we covered, on average, 81% of the species diversity across all samples. To confirm this result, we computed the ACE and Jackknife estimator indexes; both had higher values than the observed richness, indicating an underestimation of the gut microbial diversity (see Table S2). We characterized the species diversity in our set of mosquito midguts using the species richness, the Shannon diversity index (H) and the Simpson's diversity index (D); data are shown in Table S2. No significant differences in the richness index were found comparing mosquito locality and/or P. falciparum prevalence. Significant differences of the diversity indexes were found when comparing mosquito locality (Shannon, P = 0.0091 and Simpson, P = 0.0097) but none when comparing the infection prevalence of mosquitoes. Thus, at the genus level, the microbiota of mosquitoes from Mvan was more diverse than that from Nkolondom, but Plasmodium-infected and non-infected mosquitoes did not exhibit differences in their microbial diversity. The relationship between the class taxonomic rank of bacteria and the origin of the mosquitoes (locality) was evaluated using redundancy analyses (RDA) (Figure 3). The Eigen values of the first four axes were recorded at 0.304, 0.214, 0.331, and 0.144, respectively. The first two constrained axes explained around 50% of the total variance in the bacterial community and 100% of the species-environment relationship. The unrestricted Monte Carlo permutation test (n = 499) indicated that all environmental variables were significant (Nkolondom variable: F = 10.99, P = 0.002; Mvan: F = 13.26, P = 0.004; Mvan and laboratory variables fit collinearly). Thus the different classes of bacteria were not randomly distributed but linked to the breeding site where mosquitoes grew up. As seen in Figure 2, Flavobacteria were related to laboratory mosquitoes, whereas Alphaproteobacteria were less diverse and related to the Nkolondom locality. The remaining classes cluster along the Mvan locality. These results confirm the higher diversity of bacterial taxa in mosquitoes collected in Mvan as compared with Nkolondom, and the paucity of the gut microbiota in laboratory-reared mosquitoes as compared with mosquitoes from the wild. 10.1371/journal.ppat.1002742.g003 Figure 3 Redundancy analysis for gut bacterial communities (taxonomic rank = class) in field and laboratory mosquitoes. The length of arrows indicates the strength of correlation between the variable and the ordination scores. Blue arrow: bacterial classes, green arrow: environmental variables. The Monte Carlo permutation test was used to test the statistical significance of the relationship between environmental variables and the bacterial classes. The “Flavo” (Flavobacteriaceae) segregates with “labo” environmental variable, “Alpha” (Alphabacteriaceae) with the “NKD” environmental variable (P 90%) felt into few taxa and only 21 bacterial families, and 28 genera had an abundance of >1% in at least one mosquito midgut. Thus, the mosquito midgut is colonized by few dominant bacteria species, probably involved in metabolic functions. We used three different pairs of primers (S1, S2, and S3) to amplify and analyze each midgut sample. The comparison of the bacterial diversity for the three libraries revealed similar patterns; however, the S1 library allowed more detailed identification of the mosquito microbiota. These data strengthen the previous observation that the SSU rRNA gene clone libraries are biased by the choice of the set of primers used for amplification and thereby distort the revealed biodiversity [45]. To our knowledge, this is the first 454 sequencing analysis where different couples of primers have been used to identify the bacterial diversity in biological samples. The analysis ensured that the primer sets used produced the most accurate view possible of the bacterial content of the mosquito midgut. Proteobacteria represented more than 90% of the bacterial gut content in the mosquitoes from the wild, whereas in laboratory-reared mosquitoes, more than 95% of sequence tags belonged to the Flavobacteria Elizabethkingia spp. The remaining tags from laboratory mosquitoes were assigned to the members of Gammaproteobacteria (Acinetobacter, Pseudomonas), Firmicutes (Staphylococcus, Streptococcus), and the Alphaproteobacterium Asaia sp. Bacterial richness and diversity seem to be particularly poor in the laboratory mosquitoes. We identified the Elizabethkingia spp. in 68% (19/28) of the field-collected mosquitoes, at low densities ( 50%. In insects, the gut microbiota differs according to the food source, and in blood-sucking insects, bacterial content is higher after a blood meal [15], [23], [59], [60]. Here, because adult mosquitoes were fed the same diet, the high variability of taxa abundance results from individual variation, and the most abundant lineages represent the mosquito “core gut microbiota.” The existence of a core gut microbiota, by which different bacteria species are sharing metabolic functions and maintain the gut homeostasis, is now emerging [34], [61], [62]. Because alteration of the microbiota composition has been related to the development of diseases or health disorders, the next challenge is to define members of the microbial community and/or the metabolic interdependencies essential to preserve optimal gut homeostasis [34], [63]. The characterization of the mosquito core microbiota during the time course of Plasmodium infection will be a next step toward understanding the impact of gut bacteria on parasite development within the mosquito midgut. P. falciparum traverses the intestinal epithelium within 24 h after blood meal, at the peak of the digestion process; and whether parasites take advantage of intense competitive interactions for nutrient resources between bacteria to thwart the immune surveillance has to be investigated. Indeed, the gut microbiota is known to play an important role in protecting the host from potentially pathogenic microbes [64], [65]. Protection occurs through different processes: stimulation of the mosquito immune response, competition for binding sites or nutrients and production of toxins [65]–[67]. However, despite the beneficial role of the microbiota, pathogens, such as helminthes and viruses have developed strategies for exploiting the gut microbiota to promote their transmission [68], [69], [70]. For the mosquito vector, our understanding is still at an early stage for how the natural resident microbial flora of the mosquito midgut contributes to its resistance to the Plasmodium [15], [17], [18], [24]. In this study, we found that the abundance of Enterobacteriaceae is higher in P. falciparum-infected mosquitoes, suggesting that some microbe-parasite interactions may contribute to the successful development of the malaria parasite. However, whether Enterobacteriaceae have an effect on parasite survival or whether the increased level of Enterobacteriaceae is a consequence of Plasmodium development remains elusive. Alternatively, genetic factors, such as allelic polymorphism of immune genes, could regulate the variable levels of permissiveness of the mosquitoes as has been previously shown [71]. In contrast to our findings, previous studies reported the deleterious effect of bacterial infections on Plasmodium development in the mosquito [19], [23]–[25]. Of interest in this context, several Enterobacteriaceae strains were able to inhibit the development of Plasmodium species in the mosquito midgut, among them Cedecea spp., Serratia spp., and Enterobacter spp. isolated from A. albimanus, A. stephensi, or A. arabiensis [19], [22]–[24]. The Esp_Z Enterobacter strain isolated from A. arabiensis caught in Zambia [24] was not identified in any of the reads we analyzed. The possibility that this Enterobacter strain would have been absent from the PCR products because of competition with a different clade is unlikely as we used three different sets of primers. Therefore, we expect that in the gut of A. gambiae mosquitoes in Cameroon, the Esp_Z Enterobacter strain was below the 0.1% abundance threshold or absent. This Enterobacter strain was isolated in Zambia from wild-caught A. arabiensis mosquitoes, and differences in the mosquito species, as well as differences between the study areas, may explain why we did not find this bacterium in our material. Cirimotich et al. [24] recovered the Esp_Z on LB media, and culturing methods can lead to the artificial amplification of a bacterial strain present in minute amounts in an environmental sample. Therefore, it would be of interest to examine the presence/abundance of Esp_Z in wild-caught Zambian A. arabiensis using the methodologies described here. In our study, we analyzed the gut resident microbiota and revealed a positive correlation between commensal Enterobacteriaceae and Plasmodium infection, indicating that the P. falciparum infection phenotype under natural conditions results from more complex interactions than previously thought. Our data suggest a possible protective role of the Enterobacteriaceae on natural P. falciparum infection. Interestingly, it has been shown that commensal Enterobacteriaceae may promote intestinal homeostasis by enhancing immune receptors in the human colon [72]. For the mosquito, as described for the insect model Drosophila [73], gut homeostasis could be maintained through the renewal of the intestinal epithelial layer that can be altered upon bacterial killing or through immune regulation. A major challenge now will be to correlate our data with quantitative phenotyping of the immune system of the gut. Previous reports on the susceptibility of the M and S molecular forms to P. falciparum infections relay contrasting findings [9], [10]. In Cameroon, mosquitoes of the two molecular forms collected in a sympatric area exhibited similar susceptibility to P. falciparum infection [10], whereas in Senegal, mosquitoes of the S form, derived from progenies of field-collected individuals, were more susceptible to P. falciparum than those of the M form [9]. In the present study with the mosquitoes collected in natural breeding sites and infected on the same blood donor, we found that the M form was more infected than the S form. However, a marked difference in the P. falciparum prevalence was observed according to the sampling site, and larger sample sizes of sympatric M and S populations of A. gambiae will be needed to further explore any difference of Plasmodium susceptibility between the two cryptic species. We propose that the composition of the gut microbiota may influence parasite transmission, which would explain the difference in infection levels between mosquito populations from diverse environments [9], [74]. The mosquito susceptibility to Plasmodium infection is under host genetic control, and several candidate genes have a recognized role in the establishment of the pathogen in the mosquito midgut [11]–[13]. However, how the mosquito gut microbiota influences Plasmodium transmission has to be unraveled. The mosquito gut ecosystem remains poorly understood, and elucidating the precise role of the symbiotic and commensal flora on the regulation of the insect immune response and on the infection course of pathogens, such as Plasmodium parasites, will be of great interest. Pathogens and microbes likely depend on similar mechanisms for interacting with their hosts, and a better knowledge of the mosquito-microbiota interactions would open new avenues for vector disease control through manipulation of gut microbial communities. Furthermore, unraveling these strategies mounted by the parasites to cooperate with the resident microbiota will allow a better understanding of co-evolution of host-pathogen interactions. Materials and Methods Ethics statement All procedures involving human subjects used in this study were approved by the Cameroonian national ethical committee (statement 099/CNE/SE/09). The gametocyte carrier used in this study was enrolled as a volunteer after his parents had signed a written informed consent. Mosquito collection and sample characterization A. gambiae mosquitoes were sampled in aquatic habitats at the L4 and pupae stages in four localities in Cameroon using standard dipping technique [75]. In each locality, breeding sites were inspected visually for presence of larval stages. At each breeding site, 10 dips were taken with a standard dipper (300 ml) and kept in a 5-liter container for transportation to the insectary at OCEAC. Anopheline larvae were identified morphologically; non-anopheline larvae and predators were removed. Larvae were kept in their original habitat water in a 3-liter plastic bucket and resulting pupae were collected daily for 2 days. Pupae were transferred to a plastic cup containing 20 ml of water from the breeding site, and the cup was placed in a 30×30 cm cage for emergence. The remaining larval collection was discarded after 2 days to avoid bias because of putative modifications of the biotic content of the aquatic habitats. Adult mosquitoes were maintained in standard insectary conditions (27±2°C, 85±5% RH, and 12 h light/dark) and provided with 8% sterile sucrose solution. Female mosquitoes were fed on a single P. falciparum gametocyte carrier to avoid infection rate variability because of the blood donor. Infectious feeding was performed as previously described [71], [76]. Females, 3 to 5 days old, were starved for 24 h and allowed to feed on the P. falciparum gametocyte containing blood for 35 minutes through membrane feeders. Unfed and partially fed mosquitoes were removed by aspiration and discarded. Fully engorged females were kept in the insectary until dissections 8 days after the infectious blood meal. Mosquitoes were surface sterilized in 70% ethanol for 5 minutes, then rinsed twice in sterile PBS solution, and midguts were dissected and stored individually at −20°C until processing. DNA was extracted using the DNeasy Blood &Tissue Kit from Qiagen (Valencia, CA) and quantified (Nanodrop ND-1000, NanoDrop Technologies, Montchanin, DE, USA). A 20-ng aliquot of DNAs was subjected to whole-genome amplification using the GenomiPhi V2 DNA Amplification Kit (GE HealthCare, Uppsala, Sweden), and the GenomiPhi templates served to characterize molecular forms of A. gambiae and the P. falciparum infection status. Molecular forms were determined according to Fanello et al. [77] and the identification of mosquitoes that successfully developed malaria infection using a P. falciparum specific PCR amplifying a Cox gene fragment [78]. Sample selection for the metagenomic analysis A total of 32 individual midguts were subjected to the 454-sequencing analysis. We included 2 samples of midguts dissected from mosquitoes of our local colony of A. gambiae, Ngousso. The Ngousso colony was established in January 2006 from larvae collected in breeding sites of an urbanized district of Yaounde, “Ngousso.” Larval collections were conducted during a 2-month period; mosquitoes were blood fed for oviposition and then PCR screened for molecular form of A. gambiae. Ngousso mosquitoes belong to the M molecular and Forest chromosomal forms. Since then, the colony has been routinely maintained at the OCEAC insectary. The Ngousso samples served to provide an overview of the bacterial content of laboratory mosquitoes reared under standard insectary conditions and as an experimental control in this study. The 30 remaining samples were chosen among field mosquitoes fed on blood from the same gametocyte carrier. We selected both P. falciparum positive and P. falciparum negative midguts to assess putative differences of microbiota between non-infected and infected individuals in our P. falciparum-challenged mosquitoes. 454 sequencing For each individual midgut DNA sample, we generated three PCR amplicon libraries. We targeted 3 different hypervariable regions of the 16S ribosomal RNA to allow accurate detection of the bacterial communities in each sample. Indeed, previous analyses showed that the set of primers used for amplification can have a strong impact on the biodiversity revealed; some abundant clades in a given sample could be foreseen, depending on the primers used [45]. The S1 library targeting the V4 hypervariable region was obtained using the forward primer 535F (5′-GTGCCAGCAGCCGCGGTAATA-3′) and the reverse primer 789R (5′-GCGTGGACTACCAGGGTATCT-3′), the S2 library for the V5–6 region using the 326F (5′-CAAACAGGATTAGATACCCTG-3′) and the 1082R (5′-CGTTRCGGGACTTAACCCAACA-3′) primers, and the S3 library targeting the V5–6 region with the 1065F (5′-CAGGTGCTGCATGGCYGTCGT-3′) and the 1336R (5′-CGATTACTAGCGATTCC-3′) primers. Amplified DNA was purified and quantified using Picogreen fluorescent dye (Molecular Probes, Eurogen, OR). Individual libraries were processed for 454 sequencing by ligating the 454 adapters coupled to MID tags, allowing the multiplexing of samples. The amplicon libraries were pooled in two separate batches and sequenced. The MIDs and 454 linkers were ligated after the PCR amplification and the pyrotags were sequenced unidirectionally. Pyrosequencing was performed at Genoscreen (Lille, France) using a Genome Sequencer FLX Titanium (GS-FLX) system (Roche, Basel, Switzerland). In total, we recovered 663,651 sequence reads (tags) that were subjected to quality controls. All 454 sequences were deposited in Genebank (SRS281724.1 and SRS281725.1). Data processing and taxonomic assignment Tags were extracted only if they contained the combination linker-MID-primer and the complement primer sequence at the 3′end. Tags were sorted in appropriate files according to their MID barcode and converted to the forward strand when necessary. A strict dereplication step was then applied that discarded tags with unidentified nucleotides (Ns) and those longer than 350 bp or shorter than 200 bp. Dereplicated tags were sorted by decreasing number of occurrences and clustered at k = 3 number of differences as described in Stoeck et al. [79]; this pipeline resulted in determining unique sequences. We next processed taxonomic assignments by implementing a new approach that clearly optimizes the successive assignments. We first extracted from the reference sequences of SILVA (release 106) domains corresponding to the various possible couples of the primers. Extraction was first performed requiring a perfect match between each primer and a sequence and, when no match was found, 1, 2 and 3 differences between each primer and a sequence were used successively. This pipeline then gave three reference databases, one per amplified 16S rDNA region, containing all reference amplicons putatively matching our tags. The S1, S2 and S3 tagged databases contained 424.634, 359.198 and 394.370 reference amplicons, respectively. In a second step, all unique tags were assigned a taxon using a global alignment method. Each amplicon of the reference database was considered if it had at least 70% similarity with a tag. The list of reference amplicons was sorted by decreasing percentage of similarity and rounded to an integer. For taxonomic assignments, the reference sequence with the highest percentage was used, and taxonomy to a given level was obtained by the consensus of these taxonomies when more than one result emerged. For example, a tag with 98% similarity to the class Gammaproteobacteria and Alphaproteobacteria was only assigned to the phylum Proteobacteria. When similarity was <80%, sequences were not assigned. Tags were clustered into OTUs according to their consensus taxonomy. For each mosquito sample and each amplified 16S rDNA region, OTU abundances represent relative abundances, the number of reads for the given OTU divided by the total of tags in the SSU region of that sample. Rarefaction curves were produced by plotting the number of unique sequence tags as a function of the number of randomly sampled tags. To generate rarefaction curves, we retained OTUs containing at least 2 sequence tags and encompassing the abundance threshold of 0.04% because rare sequences likely represent random sequencing errors and overestimate the overall diversity. Ecological indexes and statistical analysis Ecological indexes such as richness and diversity indexes (Simpson, Shannon), were computed using the Vegan [80] and BiodiversityR [81] packages under the R software (available at http://www.R-project.org) [82]. Chao1, ACE1 and Jackknife richness estimators were calculated using the SPADE software [83]. Indexes were calculated using values from the genus taxonomic rank, the lowest rank obtained with the 454 technology. Association between microbiota and environmental variables was assessed using a multivariate ordination test. We defined the different taxa present in the data set as “species variables” and the origin of the mosquitoes and the P. falciparum infection status as “environmental variables” for each individual. A detrended canonical correspondence analysis (DCA) was performed to determine the ordination method suitable for our data. The longest gradient we obtained was shorter than 3.0, indicating that the constrained form of linear ordination method, the RDA, was the most appropriate test [84]. Redundancy analysis was performed using Canoco v. 4.5 Software [85]. The environmental variables were set as dummy variables (0 or 1 values). RDA and associated Monte Carlo permutation tests (n = 499) were used to identify the measured environmental variable that contributed most significantly to the variation in the bacterial community data. The Monte Carlo test returns a p value associated with the effect of the environmental variable on the microbiota composition of the samples. Results were visualized on a biplot ordination diagram using CanoDraw extension. Supporting Information Figure S1 Rarefaction analyses for each mosquito midgut sample. Saturation curves were generated by plotting the number of unique sequence tags as a function of the number of randomly sampled tags. Tags were clustered at k = 3 differences, and OTUs were set when containing at least 2 sequence tags and encompassing the abundance threshold of 0.04%. (TIF) Click here for additional data file. Figure S2 Redundancy analysis for gut bacterial communities (taxonomic rank = family) in field mosquitoes. The length of arrows indicates the strength of correlation between the variable and the ordination scores. Blue arrow: bacterial classes, green arrow: environmental variables. The Monte Carlo permutation test was used to test the statistical significance of the relationship between environmental variables and the bacterial classes. “Family” and “Locality” variables segregate along the first axis, but “Entero” (Enterobacteriaceae) and “Infection” gathers along the second axis (P<0.05). (TIF) Click here for additional data file. Table S1 Molecular form identification and infection status of female A. gambiae mosquitoes collected at larval stage in different localities and challenged after emergence to a single P. falciparum donor. (DOC) Click here for additional data file. Table S2 Mosquito characteristics and diversity indexes for each individual analyzed upon 454 sequencing of the S1 library. (DOC) Click here for additional data file. Table S3 Main genera that composed the natural A. gambiae gut microbiota. (DOC) Click here for additional data file.
              Bookmark
              • Record: found
              • Abstract: found
              • Article: found
              Is Open Access

              Malaria in Brazil: an overview

              Malaria is still a major public health problem in Brazil, with approximately 306 000 registered cases in 2009, but it is estimated that in the early 1940s, around six million cases of malaria occurred each year. As a result of the fight against the disease, the number of malaria cases decreased over the years and the smallest numbers of cases to-date were recorded in the 1960s. From the mid-1960s onwards, Brazil underwent a rapid and disorganized settlement process in the Amazon and this migratory movement led to a progressive increase in the number of reported cases. Although the main mosquito vector (Anopheles darlingi) is present in about 80% of the country, currently the incidence of malaria in Brazil is almost exclusively (99,8% of the cases) restricted to the region of the Amazon Basin, where a number of combined factors favors disease transmission and impair the use of standard control procedures. Plasmodium vivax accounts for 83,7% of registered cases, while Plasmodium falciparum is responsible for 16,3% and Plasmodium malariae is seldom observed. Although vivax malaria is thought to cause little mortality, compared to falciparum malaria, it accounts for much of the morbidity and for huge burdens on the prosperity of endemic communities. However, in the last few years a pattern of unusual clinical complications with fatal cases associated with P. vivax have been reported in Brazil and this is a matter of concern for Brazilian malariologists. In addition, the emergence of P. vivax strains resistant to chloroquine in some reports needs to be further investigated. In contrast, asymptomatic infection by P. falciparum and P. vivax has been detected in epidemiological studies in the states of Rondonia and Amazonas, indicating probably a pattern of clinical immunity in both autochthonous and migrant populations. Seropidemiological studies investigating the type of immune responses elicited in naturally-exposed populations to several malaria vaccine candidates in Brazilian populations have also been providing important information on whether immune responses specific to these antigens are generated in natural infections and their immunogenic potential as vaccine candidates. The present difficulties in reducing economic and social risk factors that determine the incidence of malaria in the Amazon Region render impracticable its elimination in the region. As a result, a malaria-integrated control effort - as a joint action on the part of the government and the population - directed towards the elimination or reduction of the risks of death or illness, is the direction adopted by the Brazilian government in the fight against the disease.
                Bookmark

                Author and article information

                Contributors
                Role: ConceptualizationRole: Data curationRole: Formal analysisRole: Funding acquisitionRole: InvestigationRole: MethodologyRole: Writing – original draftRole: Writing – review & editing
                Role: ConceptualizationRole: Data curationRole: Formal analysisRole: Funding acquisitionRole: InvestigationRole: MethodologyRole: Project administrationRole: ResourcesRole: Writing – original draftRole: Writing – review & editing
                Role: Data curationRole: Formal analysisRole: InvestigationRole: Project administrationRole: Writing – original draftRole: Writing – review & editing
                Role: InvestigationRole: Writing – review & editing
                Role: Data curationRole: InvestigationRole: ResourcesRole: Writing – review & editing
                Role: InvestigationRole: Writing – review & editing
                Role: InvestigationRole: Writing – review & editing
                Role: Funding acquisitionRole: ResourcesRole: Writing – review & editing
                Role: ConceptualizationRole: Formal analysisRole: Funding acquisitionRole: InvestigationRole: Project administrationRole: ResourcesRole: Writing – review & editing
                Role: Formal analysisRole: InvestigationRole: Writing – review & editing
                Role: Data curationRole: InvestigationRole: Writing – review & editing
                Role: InvestigationRole: ResourcesRole: Writing – review & editing
                Role: ResourcesRole: Writing – review & editing
                Role: InvestigationRole: Writing – review & editing
                Role: Data curationRole: Formal analysisRole: InvestigationRole: ResourcesRole: Writing – review & editing
                Role: Formal analysisRole: InvestigationRole: ResourcesRole: ValidationRole: Writing – review & editing
                Role: ConceptualizationRole: Formal analysisRole: Funding acquisitionRole: InvestigationRole: MethodologyRole: Project administrationRole: ResourcesRole: SupervisionRole: ValidationRole: Writing – original draftRole: Writing – review & editing
                Role: ConceptualizationRole: Formal analysisRole: Funding acquisitionRole: InvestigationRole: MethodologyRole: Project administrationRole: ResourcesRole: SupervisionRole: ValidationRole: Writing – original draftRole: Writing – review & editing
                Role: Editor
                Journal
                PLoS Negl Trop Dis
                PLoS Negl Trop Dis
                plos
                plosntds
                PLoS Neglected Tropical Diseases
                Public Library of Science (San Francisco, CA USA )
                1935-2727
                1935-2735
                14 February 2018
                February 2018
                : 12
                : 2
                : e0006221
                Affiliations
                [1 ] International Center for Clinical Research in Malaria, Fundação de Medicina Tropical Dr. Heitor Vieira Dourado, Manaus, Brazil
                [2 ] Escola Superior de Ciências da Saúde, Universidade do Estado do Amazonas, Manaus, Brazil
                [3 ] Instituto de Pesquisas Leônidas & Maria Deane, Fundação Oswaldo Cruz, Manaus, Brazil
                [4 ] Centro de Pesquisas René Rachou, Fundação Oswaldo Cruz, Belo Horizonte, Brazil
                [5 ] Department of Entomology, Armed Forces Research Institute of Medical Sciences, Bangkok, Thailand
                [6 ] Department of Entomology, U.S. Naval Medical Research Unit No. 6 (NAMRU-6), Callao, Peru
                [7 ] Mahidol-Oxford Tropical Medicine Research Unit, Faculty of Tropical Medicine, Mahidol University, Bangkok, Thailand
                [8 ] Centre for Tropical Medicine and Global Health, Nuffield Department of Clinical Medicine, University of Oxford, Oxford, United Kingdom
                Institut Pasteur, FRANCE
                Author notes

                The authors have declared that no competing interests exist.

                Author information
                http://orcid.org/0000-0002-0848-1940
                Article
                PNTD-D-17-01546
                10.1371/journal.pntd.0006221
                5828505
                29444080
                0480930e-a74e-4694-aa9e-2b1e0a8a7086

                This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

                History
                : 16 October 2017
                : 8 January 2018
                Page count
                Figures: 7, Tables: 1, Pages: 23
                Funding
                This work was supported by National Counsel of Technological and Scientific development (CNPq), Coordination for the Improvement of Higher Education Personnel (CAPES) and Research Support Foundation of Amazonas (FAPEAM) through PPSUS and PAPAC projects supported this study. Bill & Melinda Gates Foundation has also funded this study through TransEpi Project. PFPP and MVGL are level 1 fellows from CNPq. VSS and YTP have fellowships from CAPES. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
                Categories
                Research Article
                Biology and Life Sciences
                Parasitology
                Parasite Groups
                Apicomplexa
                Plasmodium
                Biology and Life Sciences
                Anatomy
                Body Fluids
                Blood
                Medicine and Health Sciences
                Anatomy
                Body Fluids
                Blood
                Biology and Life Sciences
                Physiology
                Body Fluids
                Blood
                Medicine and Health Sciences
                Physiology
                Body Fluids
                Blood
                Biology and Life Sciences
                Parasitology
                Parasite Groups
                Apicomplexa
                Oocysts
                Medicine and Health Sciences
                Infectious Diseases
                Disease Vectors
                Insect Vectors
                Mosquitoes
                Biology and Life Sciences
                Species Interactions
                Disease Vectors
                Insect Vectors
                Mosquitoes
                Biology and Life Sciences
                Organisms
                Eukaryota
                Animals
                Invertebrates
                Arthropoda
                Insects
                Mosquitoes
                Biology and Life Sciences
                Anatomy
                Body Fluids
                Blood
                Blood Plasma
                Medicine and Health Sciences
                Anatomy
                Body Fluids
                Blood
                Blood Plasma
                Biology and Life Sciences
                Physiology
                Body Fluids
                Blood
                Blood Plasma
                Medicine and Health Sciences
                Physiology
                Body Fluids
                Blood
                Blood Plasma
                Medicine and Health Sciences
                Parasitic Diseases
                Malaria
                Medicine and Health Sciences
                Tropical Diseases
                Malaria
                Medicine and Health Sciences
                Pharmacology
                Pharmacokinetics
                Drug Metabolism
                Biology and Life Sciences
                Population Biology
                Population Metrics
                Death Rates
                Custom metadata
                vor-update-to-uncorrected-proof
                2018-02-27
                All relevant data are within the paper and its Supporting Information files.

                Infectious disease & Microbiology
                Infectious disease & Microbiology

                Comments

                Comment on this article