Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Midgut Microbial Community of Culex quinquefasciatus Mosquito Populations from India

  • Kshitij Chandel ,

    chandelkc@gmail.com

    Affiliation Vector Management Division, Defence Research and Development Establishment, Jhansi Road, Gwalior, India

  • Murlidhar J. Mendki,

    Affiliation Vector Management Division, Defence Research and Development Establishment, Jhansi Road, Gwalior, India

  • Rasesh Y. Parikh,

    Affiliation National Centre for Cell Sciences, Pune University Campus, Pune, Maharashtra, India

  • Girish Kulkarni,

    Affiliation National Centre for Cell Sciences, Pune University Campus, Pune, Maharashtra, India

  • Sachin N. Tikar,

    Affiliation Vector Management Division, Defence Research and Development Establishment, Jhansi Road, Gwalior, India

  • Devanathan Sukumaran,

    Affiliation Vector Management Division, Defence Research and Development Establishment, Jhansi Road, Gwalior, India

  • Shri Prakash,

    Affiliation Vector Management Division, Defence Research and Development Establishment, Jhansi Road, Gwalior, India

  • Brahma D. Parashar,

    Affiliation Vector Management Division, Defence Research and Development Establishment, Jhansi Road, Gwalior, India

  • Yogesh S. Shouche,

    Affiliation National Centre for Cell Sciences, Pune University Campus, Pune, Maharashtra, India

  • Vijay Veer

    Affiliation Vector Management Division, Defence Research and Development Establishment, Jhansi Road, Gwalior, India

Abstract

The mosquito Culex quinquefasciatus is a ubiquitous species that serves as a major vector for west nile virus and lymphatic filariasis. Ingestion of bloodmeal by females triggers a series of physiological processes in the midgut and also exposes them to infection by these pathogens. The bacteria normally harbored in the midgut are known to influence physiology and can also alter the response to various pathogens. The midgut bacteria in female Cx. quinquefasciatus mosquitoes collected over a large geographical area from India was studied. Examination of 16S ribosomal DNA amplicons from culturable microflora revealed the presence of 83 bacterial species belonging to 31 bacterial genera. All of these species belong to three phyla i.e. Proteobacteria, Firmicutes and Actinobacteria. Phylum Proteobacteria was the most dominant phylum (37 species), followed by Firmicutes (33 species) and Actinobacteria (13 species). Phylum Proteobacteria, was dominated by members of γ-proteobacteria class. The genus Staphylococcus was the largest genus represented by 11 species whereas Enterobacter was the most prevalent genus and recovered from all the field stations except Leh. Highest bacterial prevalence was observed from Bhuj (22 species) followed by Nagrota (18 species), Masimpur (18 species) and Hathigarh (16 species). Whereas, least species were observed from Leh (8 species). It has been observed that individual mosquito harbor extremely diverse gut bacteria and have very small overlap bacterial taxa in their gut. This variation in midgut microbiota may be one of the factors responsible for variation in disease transmission rates or vector competence within mosquito population. The present data strongly encourage further investigations to verify the potential role of the detected bacteria in mosquito for the transmission of lymphatic filariasis and west nile virus. To the best of our knowledge this is the first study on midgut microbiota of wild Cx. quinquefasciatus from over a large geographical area.

Introduction

Among the disease transmitting insects, the mosquitoes are the primary hosts for transmission of diseases like malaria, dengue, chikungunya, lymphatic filariasis, yellow fever etc., which together are responsible for several million deaths and hundreds of millions of cases every year. Culex quinquefasciatus Say is an important vector of west nile virus and filarial nematode, Wuchereria bancrofti. The later causes lymphatic filariasis in humans and presently over 120 million peoples are infected with filarial worm[1]. Isolation of Japanese encephalitis virus was also reported from field collected Cx. quinquefasciatus mosquitoes. [2].

Lymphatic filariasis is a major public health problem in India and Worldwide, it is estimated that 1.3 billion people from 83 countries are living at the risk of infection; however in India LF is endemic in 17 states and six union territories, and is responsible for one third of the global disease burden with about 554.2 million people at risk of infection, with 31 million parasite carriers and 23 million cases of symptomatic filariasis[1].

The bacteria colonizing midguts of insect vectors have drawn special attention for their interaction with both the insect hosts and pathogenic organisms[3]. Little is known about the midgut microflora of Culex mosquitoes and very few studies has been conducted to study the midgut microbiota of Culex mosquitoes [4], [5]. A comprehensive study has been conducted by us on the diversity of microbiota in the midgut lumen of Cx. quinquefasciatus in order to ascertain their potential role in disease transmission and for their exploitation in vector management. In the present study, wild populations of Cx. quinquefasciatus were collected from different locations from India to investigate geographical variation in their midgut bacterial community.

Materials and Methods

Sample collection

Study locations comprised 10 field stations situated in different climatic zones viz. coastal, arid, semi-arid, mountainous and subtropical (Fig. 1). Mosquitoes were collected from army cantonments from these field stations, with the permission of Director General Armed Forces Medical Services (DGAFMS), Ministry of Defence, Government of India. From each station, indoor resting adult Cx. quinquefasciatus females were collected from different sites at dawn and dusk. Collected mosquitoes were kept in pre-sterilized cages. Mosquitoes were anesthetized with chloroform and species was identified morphologically using standard taxonomic key. Adult mosquitoes were surface sterilized with 70% ethanol for 5 min followed by washing in phosphate buffered saline (PBS) (twice) before dissection for further processing of midgut isolates and bacterial cultivation according to Pidiyar et. al.[4]. Midguts were microscopically dissected out under sterile conditions and transferred individually to 100 µl of brain heart infusion broth (BHI broth) in 1.5 ml microcentrifuge tubes and enriched for 4 hrs at room temperature. After enrichment, an equal volume of 40% glycerol was added in each tube and stored in liquid nitrogen and transported to the laboratory for further study. From each station, 100 mosquitoes were dissected out and midgut samples were collected for bacterial diversity. To confirm the sterility of the procedure, two controls were also taken from each station which contained PBS from a mosquito's second wash.

thumbnail
Figure 1. Collection sites of Culex quinquefasciatus mosquito in India.

1: Jamnagar(JMN), 2: Bhuj(BHJ), 3: Barmer(BAR), 4: Jodhpur(JOD), 5: Bathinda(BHT), 6: Amritsar(ASR), 7: Nagrota(NAG), 8: Leh(LEH), 9: Hathigarh(HGR), 10: Masimpur(MSM).

https://doi.org/10.1371/journal.pone.0080453.g001

Isolation of bacterial flora

Midgut contents were ten-fold serially diluted up to 10−8 dilution in PBS and 100 µl of each dilution was spread on tryptose soya agar (TSA) supplemented with 5% sheep blood and incubated at 30°C for 48–72Hrs. The resulting bacterial colonies were grouped based on their colony morphology. Bacterial colonies that were morphologically distinct were selected and subcultured on TSA plates until a pure culture was obtained for further analysis.

DNA extraction and PCR amplification of 16S rDNA

Pure bacterial isolates from mosquito midguts were subcultured in 5 ml tryptose soya broth (TSB) at 30°C for 24 hrs. Cell pellets were suspended in distilled water and lysed using repeated cycles of freezing and thawing, lysozyme and proteinase K treatment. DNA was isolated by phenol∶chloroform∶isoamyl alcohol (24∶25∶1) extraction and isopropanol precipitation protocol[6]. Complete 16S rRNA gene (Approx. 1.5 kb size) was amplified from extracted DNA of isolates as described by Pidiyar et al., [4] using eubacteria specific primers 27F 5′– CCA GAG TTT GAT CMT GGC TCA G – 3′ and 1525 R 5′ –TTC TGC AGT CTA GAA GGA GGT GWT CCA GCC-3′. Amplification of the 16S rRNA gene was confirmed by gel electrophoresis using 1% agarose. PCR products were purified using PEG-NaCl precipitation[6] and then sequenced on an ABI 3730 (Applied Biosystems Inc. Foster City, CA) automated DNA sequencer using the Big Dye termination kit at the National Centre for Cell Sciences (NCCS), Pune, India. For sequencing, three internal primers were used for generating large overlapping sequences. These sequences were then assembled into full length amplicons for analysis. Sequences were submitted to GenBank under the accession numbers JN644482 to JN644627.

Sequence analysis

To identify the closest related sequence, obtained sequences were analyzed at the EzTaxon server (http://147.47.212.35:8080)[7]. For phylogenetic analysis, related sequences were retrieved from the database and multiple sequence alignment of sequences was carried out using the CLUSTAL W program available at http://www.ebi.ac.uk/Tools/msa/clustalw2/. Unaligned sequences at the beginning and end of the alignment file were trimmed using DAMBE Ver. 5.2.9 (http:/www.dambe.bio.uottawa.ca) and the alignment file was converted to MEGA format. Phylogenetic trees were constructed using >1200 base pair aligned sequences by the neighbor joining method using Kimura 2 distances parameter in MEGA 4.0 package[8]. One thousand bootstrap replicates were generated, and a consensus tree was derived.

Statistical analysis

To study the species richness and diversity the following diversity indices were estimated: Simpson Index [9], Shannon Index, Evenness and Sφrrensen Index [10]. Good's coverage was calculated using following formula[11]: “Percent coverage  = [1−(n/N)]×100” Where n = number of bacterial species represented by single isolate and N =  is the total number of bacterial isolates from particular location.

Results

Microbial diversity in the midgut of the Cx. quinquefasciatus mosquito was analyzed for the presence of aerobic bacteria. Adult female Cx. quinquefasciatus mosquitoes were collected during the post monsoon season where mosquito population density was quite high.

To study aerobic bacterial diversity, we used a 16S rRNA gene sequence based approach. Characterization of gut microbes in Cx. quinquefasciatus collected from all 10 locations using a culture dependent method led to the identification of 31 bacterial genera, including 83 species from 3 phyla. The 16S rDNA sequence of bacterial isolates were aligned with reference strains available in GenBank and used for construction of phylogenetic trees as shown in figures S1S3 (supplementary material), revealed the relatedness among the bacteria identified. The most prominent phylum was Proteobacteria with the greatest number of species isolated from mosquito midguts and included 37 bacterial species (45%) from 15 genera. Firmicutes was the second largest phylum with 33 bacterial species (40%) from 9 genera. The least frequently isolated phylum was Actinobacteria where only 13 (15%) species from 7 genera were identified.

Bhuj belongs to a coastal climatic zone. In this region, 65 samples (65%) of mosquito midguts were found to be positive for bacterial growth and 22 bacterial species from 14 genera were reported from this region (Table 1). Of the 22 bacterial species, 50% (11 species) belong to phylum Proteobacteria and 41% (9 species) and 9% (2 species) belong to phyla Firmicutes and Actinobacteria, respectively (Table 2). Bacillus thuringiensis was the most frequently isolated bacterial species and isolated from 36% of the samples analyzed (Table 2). Four species viz., Bacillus nealsonii, Staphylococcus caprae, Acinetobacter baumannii and Pantoea anthophila were isolated only once. Percent coverage of bacterial diversity from this region was 97.6% (Table 1).

thumbnail
Table 1. showing occurrence of bacterial taxa identified, Good's coverage and diversity indices at 10 different field stations.

https://doi.org/10.1371/journal.pone.0080453.t001

thumbnail
Table 2. of bacterial species identified from Cx. quinquefasciatus midgut and their occurrence in each location with their sequence accession number.

https://doi.org/10.1371/journal.pone.0080453.t002

Second location from coastal region is Jamnagar, from where 69 (out of 100) samples were found to be positive for bacterial growth (Table 1) and 12 bacterial species from 9 genera and 3 phyla were isolated. In Jamnagar, Firmicutes was the most dominating phylum with 6 species (belonging to 5 genera) accounting for 50% of the total species identified. Two and three species from 2 genera belonged to phyla Proteobacteria and Actinobacteria respectively (Table 2). Most frequently isolated species was Enterococcus faecalis which was recorded in 21% of the samples analyzed, followed by Acinetobacter soli which was isolated from 19% of the samples analyzed and Enterobacter cloacae which was isolated from 17% of the samples. Arthrobacter creatinolyticus, Microbacterium imperiale and Exiguobacterium profundum were represented by single individual (Table 2). Similar to Bhuj the percent coverage from Jamnagar was 97.4% (Table 1).

From arid climatic zone two stations namely Jodhpur and Barmer situated in thar desert were selected for study. From Jodhpur, only 25% of the samples were found to be positive for bacteria and 13 species (belonging to 9 genera) from 3 phyla were identified (Table 1). Firmicutes was the most prevalent phylum, represented by 62% of the total bacterial species identified. Three species (23%) from 3 genera belong to phylum Proteobacteria and the remaining 2 species (15%) from 2 genera belong to phylum Actinobacteria. The most frequently isolated species was Pseudomonas stutzeri which is isolated from 7% of the total samples analyzed (Table 2), and in Jodhpur the percent coverage was 95.5% (Table 1).

In Barmer region, percent positivity of mosquito midguts for bacteria was only 22%. Among these, a total of 11 species (belonging to 9 genera) were identified (Table 1). Phylum Proteobacteria was the largest phylum represented by 7 species (64%) belonging to 6 genera. Two species with 2 genera belonged to Firmicutes and 2 species from 1 genus belonged to phylum Actinobacteria (Table 2). The most frequently isolated species was Morganella morganii, that was isolated from 6% of the total samples analyzed followed by Microbacterium oxydans which was recovered from 5% of the samples. The percent coverage in Barmer region was calculated as 85.7% (Table 1).

Amritsar is a semiarid climatic zone showed a 25% positivity of mosquito midgut for presence of bacteria. A total of 13 species belonging to 8 genera and 2 phyla were identified (Table 1). Phylum Proteobacteria was the most dominant phylum, wherein 8 of 13 species (72%) belonged to this phylum. The remaining 5 species (38%) belonged to phylum Firmicutes (Table 2). Pantoea dispersa was the most frequently isolated bacterial species and was isolated from 5% of the samples. Five species, Bacillus subtilis, Enterococcus silesiacus, Enterobacter asburiae, Escherichia hermannii and Shigella flexneri were isolated only once. The Good's coverage from Amritsar was quite low and calculated 83.3% (Table 1).

In Bathinda, thirty two midguts (32%) were found to be positive for bacterial flora. A total of 15 species from 9 genera were identified from this region (Table 1). The most prevalent phylum was Firmicutes, followed by Proteobacteria, and then Actinobacteria. Eight species (53%) belonged to Firmicutes, six species (40%) belonged to Proteobacteria, and 1 species (7%) belonged to phylum Actinobacteria. Enterobacter cloacae was the most frequently isolated species and was isolated from 8% of the total samples. In this region, Bacillus circulans, Bacillus thuringiensis, Staphylococcus haemolyticus and Providencia alcalifaciens were represented by single isolate (Table 2). In comparison to Amritsar, the percent coverage from Bathinda was quite high and calculated as 90.1% (Table 1).

The mountainous climatic zones namely Nagrota and Leh were studied. In Nagrota, the percent positivity was 67%, with eighteen species (belonging to 14 genera) isolated (Table 1). In Nagrota, 50% of total species identified belonged to phylum Proteobacteria, the most dominant phylum. Five species from 4 genera belonged to phylum Firmicutes, and 4 species from 3 genera belonged to phylum Actinobacteria (Table 2). Bacillus subtilis was the most frequently isolated species and was isolated from 22% of the mosquito midguts analyzed. Two species, Kytococcus schroeteri and Aerococcus urinaeequi were represented by only one isolate. Percent coverage for bacterial diversity from Nagrota was 98.6% (Table 1).

Another field station from the mountainous climatic zone was Leh. From Leh only 20 mosquito midguts were found to be positive for bacteria and 8 species (belonging to 4 genera) from 3 phyla were isolated (Table 1). Phylum Firmicutes was the biggest phylum with 4 species (50%), followed by phylum Proteobacteria where only 3 species were isolated, and 1 species belonged to phylum Actinobacteria. Out of these, Staphylococcus epidermidis, was the most frequently isolated bacterium isolated from 8% of the total samples analyzed. Kocuria marina and Pseudomonas cuatrocienegasensis were represented by only one isolate(Table 2). The percent coverage of the bacterial diversity from Leh was 93.3% (Table 1).

Hathigarh and Masimpur represent the tropical climate of Assam. Eighty five percent of the samples from Hathigarh were found to be positive for bacteria and 16 species (belonging to 14 genera) were identified (Table 1). Phylum Proteobacteria was the most dominant phylum with 11 species (69%) from 7 genera represented. Three species (19%) from 3 genera belong to phylum Firmicutes. The least frequently isolated phylum in this region was Actinobacteria with 2 species (12%) belonging to 2 genera (Table 2). In Hathigarh, Aeromonas hydrophila was the most frequently isolated species identified in 31% of the total samples analyzed followed by Escherichia coli and Klebsiella pneumoniae which were recovered from 24 and 22% of the mosquito population. Three species viz., Lysinibacillus macroides, Acinetobacter schindleri and Aeromonas enteropelogenes were represented by single isolates. Bacterial diversity from Hathigarh was quite rich as Good's percent coverage was recorded 98.4% (Table 1).

The percent positivity was 82% with 18 species (belonging to 14 genera) from 2 phyla isolated from Masimpur (Table 1). In this region, 56% of the isolated bacterial species (10 species from 6 genera) belonged to phylum Proteobacteria, 33% (6 species from 5 genera) belonged to phylum Firmicutes and the remaining 11% (2 species from 2 genera) belonged to phylum Actinobacteria (Table 2). The most frequently isolated species was Staphylococcus saprophyticus which was isolated from 34% of the total samples studied followed by Klebsiella pneumoniae which was recovered from 31% samples. The least frequently isolated species were Acinetobacter beijerinckii and Shigella flexneri which were represented by single isolates. The Good's coverage at Masimpur was 99% (Table 1).

Simpson richness index, Shannon diversity index and Evenness diversity indices were estimated from all the locations and are summarized in Table 1. The values of Simpson diversity index ranged from 0.813 to 0.908, the lowest value was estimated for Leh and the highest value was for Nagrota region. Shannon diversity index is another widely used index for comparing diversity between various habitats [12]. In the present study, values of Shannon (H) index ranged from 1.832 to 2.615 (Table 1) with the highest value for Bhuj and lowest value for the Leh region. Evenness was used for estimating how well the species are evenly distributed among the individuals. Values of evenness were between 0.622 – 0.852 (Table 1). The highest evenness was recorded for Amritsar and the least was estimated for the Bhuj region. Sφrrenson coefficient is also very widely used similarity index. Pair wise matrix of Sφrrenson index and number of species shared by two populations are shown in Table 3. The values of Sφrrenson index ranged from 0 to 0.529 with the highest value calculated between Masimpur and Hathigarh. The lowest values of Sφrrenson index (0) were calculated between the communities where none of the species was shared.

thumbnail
Table 3. Pair wise matrix of common species (Upper half) and Sφrrensen similarity index (Lower half).

https://doi.org/10.1371/journal.pone.0080453.t003

Discussion

The bacteria inhabiting mosquito midgut has drawn special attention in recent past due to their interactions with both mosquito host as well as disease causing parasites. The present work generated detailed information about aerobic bacterial flora in the midgut lumen of field caught Cx. quinquefasciatus female mosquitoes. Richness and diversity of microbes associated with the field collected adult mosquitoes was found to be quite high in all the samples.

Midgut contents from any of the field stations have not shown presence of anaerobic bacteria or fungi (data not shown). Though Aspergillus sp. has been isolated from Anopheles stephensi larval gut but adult mosquitoes did not harbour any fungus[13]. Absence of fungus in adult gut indicates that the environment of midgut lumen of the adult mosquito is not conducive for the fungi [14]. Similarly anaerobic bacteria were also not recovered from midgut content from any of the field stations.

Conversely, a large number of aerobic bacteria were isolated from the Culex midgut contents from all the field stations. Our result demonstrates that the aerobic microbial flora of the adult mosquito midgut is a complex one and is dominated by gram negative gammaproteobacteria. The maximum number of the bacterial species was identified from Bhuj followed by Nagrota, Masimpur and Hathigarh. The least number of bacterial species were isolated from Leh region (8 species belonging to 4 genera). The climate of Leh is not favorable for mosquito breeding and mosquitoes can only breed in the short period where the temperature regime is suitable for their breeding. Mosquitoes mainly acquire bacteria either from their larval habitat or from the environment during the process of nectar feeding or blood feeding. Because of very low temperatures, environmental bacterial load is also very low which was also reflected in mosquito midgut bacterial diversity. High numbers of bacterial species were recorded from the areas where rainfall and RH was high such as Bhuj, Nagrota, Hathigarh and Masimpur . Though the rainfall and high humidity stimulates the bacterial growth in environment, in the present study no significant correlation was observed with rainfall, relative humidity and temperature. Since many biotic and abiotic factors account for the distribution of microbes in environment, it is possible that other unidentified and ecologically important factors, or interactions between these factors may account for the differences in richness observed across populations.

The extent of prevalence of culturable bacteria in different species of mosquitoes and in different population of the same species seems to be quite variable. In present study percent positivity for presence of bacteria ranged from 20% to 85% and maximum incidence was recorded in Cx. quinquefasciatus population from Hathigarh whereas least was found in Leh region. Most of the species were unevenly distributed among the host population. Overall bacterial prevalence is less than 50%. Similarly Lindh et. al. (2005) reported bacterial prevalence of culturable bacteria among 15% of the mosquito population[15]. However, Djadid et. al. (2011) observed very high bacterial prevalence in field collected An. stephensi from Iran but could not retrieve any culturable bacteria from midgut of Anopheles maculipennis[16]. It can be inferred from present study that either possibly only a part of mosquito population acquires bacteria or large number of bacteria are not stable midgut residents. Alternatively, it could also be possible that mosquitoes harbouring bacteria were not culturable under the culture conditions provided in present study.

The diversity index quantifies diversity in a community and describes its numerical structure. Simpson gave the probability of any two individuals drawn at random from an infinitely large community belonging to different species [9]. In the present study, the values of Simpson index ranged from 0.813 to 0.908 with the highest and the lowest between Nagrota and Leh, respectively. The value of Simpson index increases with diversity. The Shannon (H) index is another widely used index for comparing diversity between various habitats[12]. In the present study, the values of Shannon index ranged between 1.832 to 2.615 being the highest for Bhuj and lowest from the Leh region. Generally the values of Shannon index ranged between 1.5 to 3.5, where values above 3 indicates that the diversity of the habitat is rich and stable, whereas values under 1.0 indicate that there is pollution and degradation of habitat structure and diversity is not stable. The Shannon index is a representation of both species abundance and evenness, when either of these two factors increases, the diversity index also increases. The Shannon index has maximum values when all the species are equally abundant in community. Evenness was used for the estimating how well the species are evenly distributed among the individuals. The highest evenness was recorded for Amritsar station indicating that bacterial species at Amritsar are evenly distributed among the individuals as compared to other field stations. The lowest evenness was recorded from Bhuj region indicating that species are less evenly distributed and some species more dominant than the others. The Sφrrenson coefficient is a very widely used similarity index. This mainly depends upon the number of species common in two populations analyzed. Values of the Sφrrenson Index vary between 0 to 1 where 0 indicates none of the species is shared by two communities, whereas a value of 1 indicates that both the communities having exactly the same bacterial composition. Values of Sφrrenson index ranged from 0 to 0.529 in the present study. The maximum numbers of common species were found between Masimpur and Hathigarh, which resulted in the highest Sφrrenson value. Nine bacterial species were shared by these two stations and the value of Sφrrenson index is 0.529. Conversely, none of the species were shared between LEH-BHJ, LEH-JMN, LEH-BAR, LEH-ASR, LEH-NAG and LEH-MSM, hence the values of Sφrrenson index was also 0 for these pairs of stations.

There were a large number of soil and environmental bacteria isolated in this study, such as species of Acinetobacter, Microbacterium, Micrococcus, Stenotrophomonas and Bacillus sp., and these isolates showed regional variation. This suggests that the local soil and water environment plays an important role in colonization of the mosquito midgut with regional bacteria encountered at breeding sites or during nectar or blood feeding.

Most of the genera recovered in present study have already been reported from midgut of various mosquito vector species. Of the 31 bacterial genera 23 has been previously reported to inhabit mosquito gut. However, 8 genera Arthrobacter, Janibacter, Kytococcus, Leucobacter, Aerococcus, Sporosarcina, Vagococcus and Delftia are the first report as midgut microflora of Culex mosquito. Members of genera Acinetobacter, Aeromonas, Bacillus, Enterobacter, Enterococcus, Klebsiella, Pantoea, Pseudomonas, Serratia, Staphylococcus and Stenotrophomonas, have been frequently reported from mosquito gut in previous studies and our results are consistent with those of the earlier reports[4], [5], [15], [17][28]. This suggests that at least a fraction of mosquito midgut inhabitants could be common for different mosquito species inhabiting similar environments and may represent evolutionary conservation of association between bacteria and mosquito gut. Comparative analysis of bacterial diversity from adult Culex mosquito revealed the high prevalence of genus Enterobacter (14.4%), Bacillus (11.9%), Staphylococcus (10.0%), Enterococcus (9.5%), Acinetobacter (9.4%), Klebsiella (8.2%), Pseudomonas (5.9%) and Aeromonas (4.8%) in present study. Bacterial isolates belong to genus Enterobacter were recovered from all the field locations except Leh and comprise of a major part of midgut microbiota of Culex mosquitoes in the present study which is in accordance with earlier reports [4], [5], [18], [29]. It has been reported previously that species of Enterobacter are the most common bacteria isolated from insect gut[30].

The midgut bacterial infection in wild mosquito populations may influence parasite transmission and could contribute to understanding variation in vectorial capacity observed by same species in different locations because naturally existing microorganism in mosquito midgut have important roles to determine parasite survival and development. Mosquitoes are known to respond to infection by disease causing pathogen and elicit a specific immune response against them [31]. Same immune response gene are also expressed in response to midgut bacteria and this raises the possibility that the presence of specific bacteria in the midgut may alter the vectorial efficiency at which a pathogen is transmitted by a vector mosquito [24]. Thus, the midgut bacterial composition has a considerable effect on the survival of pathogens in the midgut environment. Earlier studies have indicated that the susceptibility of Cx. quinquefasciatus mosquito for JE virus increases when Pseudomonas sp. and Acinetobacter sp. were incorporated in the mosquito blood meal [32]. Similarly, in a study on occurrence of Klebsiella sp. and Pseudomonas sp. in mosquito midguts reported that mosquitoes with Pseudomonas in their midgut showed a higher prevalence of malaria sporozoites, whereas females infected with Klebsiella sp. could not support parasite development [33], [34]. Enterobacter isolated from wild mosquito populations from Zambia inhibits midgut epithelium invasion of Plasmodium falciparum through generation of reactive oxygen species (ROS) in Anopheles mosquitoes and resulted in inhibition of parasite development in mosquito gut [35]. Similarly, improved conversion of oocystes from ookinets was observed when mosquito fed with plasmodium infected blood containing antibodies raised against midgut bacteria[36].

The large numbers of bacteria present in the vector midgut are capable of producing factors that can kill parasites. Haemolysin produced by Enterobacter cloacae, Serratia marcescens, Escherichia coli, Enterococcus faecalis exerts activity against both prokaryotic and eukaryotic cells[3]. Other molecules produced by Serratia marcescens are protease and prodigiosin. Prodigiosin has been shown to have potent activity against Trypenosoma curzi and derivatives of prodigiosin also have shown marked activity against Plasmodium falciparum [3].

To the best of our knowledge this is the first attempt at comparative cataloguing of midgut microbiota of Cx. quinquefasciatus mosquitoes over a large geographical area in India. Identification and characterization of mosquito midgut flora is likely to contribute towards better understanding of mosquito biology including longevity, reproduction and mosquito microbe interaction that may be important to develop novel strategies for vector control. We embarked on this microbiological survey to identify bacteria that might be used in future efforts to develop paratransgenesis as a mechanism to block vectorial transmission by mosquito. Riehle et al., have mentioned essential requirement for suitable candidate for paratransgenesis[37] and some bacteria isolated in present study could be a suitable candidate for screening for paratransgenesis.

Supporting Information

Figure S1.

Dendrogram showing phylogenetic affiiliation of bacterial isolates belonging to phylum Actinobacteria. The tree was constructed using neighbor joining algorithm with Kimura 2 parameter distances. Number at the nodes indicate percent bootstrap values (1000 replicates). The bar indicates the Jukes-Cantor evolutionary distance.

https://doi.org/10.1371/journal.pone.0080453.s001

(TIF)

Figure S2.

Dendrogram showing phylogenetic affiiliation of bacterial isolates belonging to phylum Firmicutes. The tree was constructed using neighbor joining algorithm with Kimura 2 parameter distances. Number at the nodes indicate percent bootstrap values (1000 replicates). The bar indicates the Jukes-Cantor evolutionary distance.

https://doi.org/10.1371/journal.pone.0080453.s002

(TIF)

Figure S3.

Dendrogram showing phylogenetic affiiliation of bacterial isolates belonging to phylum Proteobacteria. The tree was constructed using neighbor joining algorithm with Kimura 2 parameter distances. Number at the nodes indicate percent bootstrap values (1000 replicates). The bar indicates the Jukes-Cantor evolutionary distance.

https://doi.org/10.1371/journal.pone.0080453.s003

(TIF)

Acknowledgments

The authors are thankful to the Prof. (Dr.) M. P. Kaushik, Director, DRDE, Gwalior for encouragement, interest and for providing all necessary facilities to conduct this research work. Sincere thanks also to the National Center for Cell Sciences (NCCS), Pune for providing DNA sequencing facilities.

Author Contributions

Conceived and designed the experiments: YSS SP BDP VV. Performed the experiments: KC RYP GK. Analyzed the data: KC MJM RYP. Contributed reagents/materials/analysis tools: MJM SNT DS. Wrote the paper: KC MJM.

References

  1. 1. Agrawal VK, Sashindran VK (2006) Lymphatic Filariasis in India : Problems, Challenges and New Initiatives. Medical Journal Armed Forces India 62: 359–362.
  2. 2. Mourya DT, Ilkal MA, Mishra AC, Jacob PG, Pant U, et al. (1989) Isolation of Japanese encephalitis virus from mosquitoes collected in Karnataka state, India from 1985 to 1987. Trans R Soc Trop Med Hyg 83: 550–552.
  3. 3. Azambuja P, Garcia ES, Ratcliffe NA (2005) Gut microbiota and parasite transmission by insect vectors. Trends Parasitol 21: 568–572.
  4. 4. Pidiyar VJ, Jangid K, Patole MS, Shouche YS (2004) Studies on cultured and uncultured microbiota of wild culex quinquefasciatus mosquito midgut based on 16s ribosomal RNA gene analysis. Am J Trop Med Hyg 70: 597–603.
  5. 5. Demaio J, Pumpuni CB, Kent M, Beier JC (1996) The midgut bacterial flora of wild Aedes triseriatus, Culex pipiens, and Psorophora columbiae mosquitoes. Am J Trop Med Hyg 54: 219–223.
  6. 6. Sambrook J, Russell DW (2001) Molecular Cloning : A Laboratory Manual; 3, editor: Cold Spring Harbor Laboratory Press.
  7. 7. Chun J, Lee JH, Jung Y, Kim M, Kim S, et al. (2007) EzTaxon: a web-based tool for the identification of prokaryotes based on 16S ribosomal RNA gene sequences. Int J Syst Evol Microbiol 57: 2259–2261.
  8. 8. Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 24: 1596–1599.
  9. 9. Simpson EH (1949) Measurement of Diversity. Nature 163: 1.
  10. 10. Sorensen TA (1948) A method of establishing groups of equal amplitude in plant sociology based on similarity of species content, and its application to analyses of the vegetation on Danish commons. Kongelige Danske Videnskabernes Selskabs Biologiske Skrifter 5: 34.
  11. 11. Good IJ (1953) The population frequencies of species and the estimation of population parameters. Biometrica 40: 237–264.
  12. 12. Clarke KR, M WR (2001) Change in marine communities: an approach to statistical analysis and interpretation; 2nd, editor. U K: Primer – E Ltd. Plymouth.
  13. 13. Tajedin L, Hashemi J, Abaei M, Hosseinpour L, Rafei F, et al. (2009) Study on Fungal Flora in the Midgut of the Larva and Adult of the Different Populations of the Malaria Vector Anopheles stephensi. Iran J Arthropod Borne Dis 3: 36–40.
  14. 14. Ignatova EA, Nagornaia SS, Povazhnaia TN, Ianishevskaia GS (1996) The yeast flora of blood-sucking mosquitoes. Mikrobiol Z 58: 12–15.
  15. 15. Lindh JM, Terenius O, Faye I (2005) 16S rRNA Gene-Based identification of midgut bacteria from field-cought Anopheles gambiae Sensu Lato and A. funestus mosquitoes reveals new species related to known insect symbionts. Appl Environ Microbiol 71: 7217–7223.
  16. 16. Djadid ND, Jazayeri H, Raz A, Favia G, Ricci I, et al. (2011) Identification of the midgut microbiota of An. stephensi and An. maculipennis for their appliation as a paratransgenic tool against malari. PLoS One 6: 1–7.
  17. 17. Boissiere A, Tchioffo MT, Bachar D, Abate L, Marie A, et al. (2012) Midgut microbiota of the malaria mosquito vector Anopheles gambiae and interactions with Plasmodium falciparum infection. PLoS Pathog 8: e1002742.
  18. 18. Chao J, Wistreich G (1959) Microbial isolation from the midgut of Culex tarsalis conquillet. Journal of Insect Pathology 1: 311–318.
  19. 19. Chavshin AR, Oshaghi MA, Vatandoost H, Pourmand MR, Raeisi A, et al. (2012) Identification of bacterial microflora in the midgut of the larvae and adult of wild caught Anopheles stephensi: a step toward finding suitable paratransgenesis candidates. Acta Trop 121: 129–134.
  20. 20. Ferguson MJ, Micks DW (1961) Micro organisms associated with mosquitoes, Bacteria isolate from adult Culex fatigans Wiedmann. Journal of Insect Pathology 3: 112–119.
  21. 21. Fouda MA, Hassan MI, Al-Daly AG, Hammad KM (2001) Effect of midgut bacteria of Culex pipiens L. on digestion and reproduction. Journal of the Egyptian Society of Parasitology 31: 767–780.
  22. 22. Gonzalez L, Santillan F, Rodríguez MH, Menolez D, Availa JE (2003) Bacteria in midguts of field-collected Anopheles albimanus block Plasmodium vivax sporogonic development. Journal of Medical Entomology 40: 371–374.
  23. 23. Gusamao DS, Santos AV, Marini DC, Bacci M, Berbert-Molina MA, et al. (2010) Culture-dependent and culture-independent characterization of microorganisms associated with Aedes aegypti (Diptera: Culicidae)(L.) and dynamics of bacterial colonization in the midgut. Acta Trop 115: 7.
  24. 24. Pumpuni CB, Demaio J, Kent M, Davis JR, Beier JC (1996) Bacterial population dynamics in three anopheline species: the impact on Plasmodium sporogonic development. Am J Trop Med Hyg 54: 214–218.
  25. 25. Rani A, Sharma A, Rajagopal R, Adak T, Bhatnagar RK (2009) Bacterial diversity analysis of larvae and adult midgut microflora using culture-dependent and culture-independent methods in lab-reared and field-collected Anopheles stephensi-an Asian malarial vector. BMC Microbiol 9: 96.
  26. 26. Straif SC, Mbogo CN, Toure AM, Walker ED, Kaufman M, et al. (1998) Midgut bacteria in Anopheles gambiae and An. funestus (Diptera: Culicidae) from Kenya and Mali. J Med Entomol 35: 222–226.
  27. 27. Terenius O, Oliveira CD, Pinheiro WD, Tadei WP, James AA, et al. (2008) 16S rRNA gene sequences fro bacteria associated with adult Anopheles darlingi (Diptera : Culicidae) mosquitoes. J Med Entomol 45: 172–175.
  28. 28. Wang Y, Gilbreath TM, Kukutla P, Yan G, Xu J (2011) Dynamic gut microbiome across life history of the malaria mosquito Anopheles gambiae in Kenya. PLoS One 6: e24767.
  29. 29. Beier MS, Pumpuni CB, Beier JC, Davis JR (1994) Effects of para-aminobenzoic acid, Insulin and gentamicin on Plasmodium falciparum development in anopheline mosquitoes (Diptera, Culicidae). J Med Entomol 31: 561–565.
  30. 30. Tanada Y, Kaya HK (1993) Insect Pathology. San Diego Academic Press.
  31. 31. Dimopoulos G, Richman A, Muller HM, Kafatos FC (1997) Molecular immune responses of the mosquito Anopheles gambiae to bacteria and malaria parasites. Proc Natl Acad Sci U S A 94: 11508–11513.
  32. 32. Mourya DT, Gokhale MD, Pidiyar V, Barde PV, Patole M, et al. (2002) Study of the effect of the midgut bacterial flora of Culex quinquefasciatus on the susceptibility of mosquitoes to Japanese encephalitis virus. Acta Virol 46: 257–260.
  33. 33. Jadin J (1965) Can purple photosynthetic bacteria play a role in the sporogony of plasmodia? Bull Acad Natl Med 149: 470–472.
  34. 34. Jadin J, Vincke IH, Dunjic A, Delville JP, Wery M, et al. (1966) Role of Pseudomonas in the sporogenesis of the hematozoon of malaria in the mosquito. Bull Soc Pathol Exot Filiales 59: 514–525.
  35. 35. Cirimotich CM, Dong Y, Clayton AM, Sandifor SL, Souza-Neto J, et al. (2011) Natural microbe-mediated refractoriness to Plasmodium infection in Anopheles gamibae. Science 332: 855–858.
  36. 36. Noden BH, Vaughan JA, Pumpuni CB, Beier JC (2011) Mosquito ingestion of antibodies against mosquito midgut microbiota improves conversion of ookinetes to oocysts for Plasmodium falciparum, but not P. yoelii. Parasitol Int 60: 440–446.
  37. 37. Riehle MA, Lorena MJ (2005) Using bacteria to express and display anti parasite molecules in mosquitoes, Current and future strategies. Insect Biochemistry and Molecular Biology 35: 699–707.