Habib Ali
Abrar Muhammad
Nafiu Bala Sanda
Ying Huang1,2
Youming Hou- 1State Key Laboratory of Ecological Pest Control for Fujian and Taiwan Crops, Fujian Agriculture and Forestry University, Fuzhou, China
- 2Fujian Provincial Key Laboratory of Insect Ecology, College of Plant Protection, Fujian Agriculture and Forestry University, Fuzhou, China
Bacterial symbionts of insects affect a wide array of host traits including fitness and immunity. Octodonta nipae (Maulik), commonly known as hispid leaf beetle is a destructive palm pest around the world. Understanding the dynamics of microbiota is essential to unravel the complex interplay between O. nipae and its bacterial symbionts. In this study, bacterial 16S rRNA V3-V4 region was targeted to decipher the diversity and dynamics of bacterial symbionts across different life stages [eggs, larvae, pupae, and adult (male and female)] and reproductive organs (ovaries and testis) of O. nipae. Clustering analysis at ≥97% similarity threshold produced 3,959 operational taxonomic units (OTUs) that belonged to nine different phyla. Proteobacteria, Actinobacteria, and Firmicutes represented the bulk of taxa that underwent notable changes during metamorphosis. Enterobacteriaceae and Dermabacteraceae were the most abundant families in immature stages (eggs, larvae, and pupae), while Anaplasmataceae family was dominated in adults (male and female) and reproductive organs (ovaries and testis). The genus Serratia and Lactococcus were most abundant in eggs, whereas Pantoea and Brachybacterium represented the bulk of larvae and pupae microbiota. Interestingly the genus Wolbachia found positive to all tested samples and was recorded extremely high (>64%) in the adults and reproductive organs. The bacteria varied across the developmental stages and responsible for various metabolic activities. Selection choice exerted by the insect host as a result of its age or developmental stage could be the main reason to ascertain the shift in the bacteria populations. Maternally inherited Wolbachia was found to be an obligate endosymbiont infecting all tested life stages, body parts, and tissues. These outcomes foster our understanding of the intricate associations between bacteria and O. nipae and will incorporate in devising novel pest control strategies against this palm pest.
Introduction
Insects are the most plentiful and diverse group of animals perform numerous essential functions in our ecosystem. They have occupied almost every feasible nutritional niche and adapted to wide range of environmental conditions, diets, and habitats (Douglas, 2009, 2015; Sugio et al., 2015). The bacterial consortia of insects contributes greatly to their evolutionary success, adaptation, and diversification (Price et al., 2011; Basset et al., 2012; Visser et al., 2012; Salem et al., 2013). Diverse bacterial symbionts associated with insect play many essential roles such as digestion and provision of essential nutrients (Warnecke et al., 2007; Jones et al., 2013; Muhammad et al., 2017), provides protection against parasites and predators (Osborne et al., 2009; Koch and Schmid-Hempel, 2011), synthesis of vitamins and amino acids, detoxification of toxic compounds, enhancement of social interactions, and stimulation of the immune system (Breznak, 2000; Kwong et al., 2017), thus enhances host’s fitness and immunity (Dillon and Dillon, 2004; Genta et al., 2006; Blatch et al., 2010; Bansal et al., 2014). Experimental evidences have suggested that elimination or dysbiosis of bacterial symbionts implicates consequences on host biology. For example, in medfly, Ceratitis capitata the microbiota affects measurable physiological and behavioral parameters related to fitness (Behar et al., 2008). Buchnera endosymbiont in aphids and Candidatus symbiont in sharpshooters complement their hosts with essential amino acids, micronutrients, and vitamins (Douglas, 2009). Female tsetse flies are unable to reproduce in the absence of Wigglesworthia, and it compromises the longevity, vectorial capacity, immunity, and digestion of the host (Pais et al., 2008). Metagenome-wide association (MGWA) study have identified bacterial genes that affect the development and lipid storage in Drosophila (Chaston et al., 2014).
The cytoplasmic inherited Wolbachia is a well-known reproductive parasite that can manipulate host’s reproduction through a number of phenotypic behaviors such as cytoplasmic incompatibility (Field et al., 1999), feminization (Bouchon et al., 1998), male killing (Fialho and Stevens, 2000; Hurst et al., 2000), and induction of parthenogenesis (Pannebakker et al., 2004). Interestingly, the infection of Wolbachia is widely prevailed in phylum Arthropoda (Bouchon et al., 1998; Dedeine et al., 2001; McMeniman et al., 2009; Moreira et al., 2009; Osborne et al., 2009; Ali et al., 2016, 2018a,b,c). The ability to manipulate host’s reproduction makes Wolbachia an excellent potential bio-control agent for insect pest management (Field et al., 1999; Zabalou et al., 2004), however, Wolbachia based pest management is still at its dawn. Recently, Wolbachia infection from different life stages, body parts and tissues of Octodonta nipae has been confirmed. The higher infection density in the adult stage and reproductive organs indicates the reproductive potential of Wolbachia in O. nipae, however, it yet to be vindicated (Ali et al., 2018c).
Although less appreciated, the microbiotas of invasive insects are believed to augment their establishment in newly introduced regions, mainly through enhancing the pathogenicity of their insect host (Jiu et al., 2007) or by increasing host’s fitness and immunity (Lu et al., 2016; Zhou et al., 2016). Moreover, insight into the insect symbionts could offer an effective tool for the pest management of agriculture and medical importance (Crotti et al., 2012). In order to exploit the microbial symbionts in practical control applications, it is important to understand the active microbial diversity and their physiological role in the living system (Marzorati et al., 2008). It is also important to understand the bacterial species richness, composition dynamics, and functional organization of the microbiota to solve the practical problems using microbial resource management approach (Verstraete, 2007; Marzorati et al., 2008).
Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae), also known as nipa palm hispid beetle, is a notorious insect pest of palm cultivation in China and other palm growing countries. The pest is native to Malaysia, but due to international trade and exchange of infested goods it has disseminated many non-native regions (Sun et al., 2003; Hou and Weng, 2010; Hou et al., 2011; Peng et al., 2018). In China, it was first reported from Hainan province in 2001, feasting on California fan palm, Washingtonia filifera (Sun et al., 2003) and later in 2007 its attack was noticed from Fujian province on Phoenix canariensis and Phoenix hanceana, Trachycarpus fortune (Hou and Weng, 2010). The pest causes damage by feeding on unopened young fronds of palm render them in stunted growth, hinders plant development, and in severe cases can even cause tree death (Vassiliou et al., 2011; Zhang et al., 2015). So far, the beetle’s attack has been reported on 15 genera of the family Palamae (Sun et al., 2003; Hua et al., 2014), preferably W. filifera (Linden ex. Andre), Calamus manan (Miquel) (Steiner, 2001), Syagrus romanzoffiana (Chamisso) (Vassiliou et al., 2011), T. fortune (Hooker), and P. canariensis (Chabaud) (Hou and Weng, 2010; Hou et al., 2011). Considering its huge economic importance, much work has been done during the last few years targeting various biological (Hou and Weng, 2010; Hou et al., 2011; Li et al., 2016; Peng et al., 2018) morphological (Zhang et al., 2015; Peng et al., 2018) and immunological aspects (Meng et al., 2016; Zhang et al., 2017), however, nothing is known about their microbial profile and its impact on host physiology and development. Deep elucidation of O. nipae bacterial symbionts will not only provide a comprehensive insight to the mechanism of successful invasion in newly introduced regions but also add to the development of novel pest management tool.
Therefore, it is crucial to unravel the complex interaction between bacterial symbionts and their host insects, especially for the invasive insects such as O. nipae. This study aimed to determine the bacterial community dynamics at different life stages and reproductive organs with emphasis on Wolbachia using culture independent highthroughput sequencing of bacterial 16S rRNA gene fragment.
Materials and Methods
Insect Sampling and Maintenance
To disclose the bacterial symbionts of O. nipae, the different life stages (Supplementary Figure S1) were collected during July 2016 from infested palm trees in a palm nursery at the Fuqing Entry-Exit Inspection and Quarantine Bureau in Fuqing (25°43′42′′N, 119°20′35′E), Fujian, China. The live beetles were identified based on their morphological characters and transported in plastic boxes with the host plant tissues to the laboratory for further analysis. The insect population was maintained at controlled conditions of temperature (27 ± 2°C), humidity (70–75% RH), and a light: dark (12:12) photoperiod on the fresh leaves of Trachycarpus fortunei (Hook) as previously described (Hou et al., 2014).
Insect Dissection and DNA Extraction
To extract bacterial DNA, individuals of O. nipae were randomly picked from the F2 laboratory population. The adult males and females were dissected to remove the reproductive organs (male testis and female ovary) as described previously (Ali et al., 2018b). Three replicates of each developmental stage [≈50 eggs, mature larvae, pupae (3–5 days old)], adult [male and female (5 individuals/sample)], and reproductive organs (testis and ovaries of 5 individuals/sample) were processed for DNA extraction. To avoid contamination, samples handling and processing were carried out inside the laminar flow hood using sterilized dissecting tools and reagents. After surface sterilization with 75% ethanol and sterile double distilled water, the samples were homogenized in 200 μl buffer ATL (animal tissue lysis) in a 1.5 ml micro-centrifuge tube with the help of stainless steel beads using high-throughput homogenizer. DNA extraction was carried out using DNeasy Blood and Tissue Kit (Qiagen, Valencia, CA, United States) according to the manufacturer’s guidelines with some modification as described (Ali et al., 2018a,b). DNA concentration and quality was assessed by NanoDrop 1000 (Thermo Scientific) and by running on 1% agarose gel. PCRs were run using the universal bacterial primers 27F and 1492 R (5′-AGAGTTTGATCATGGCTCAG-3′, 5′-TACGGYTACCTTGTTACGACTT-3′) (Sangon Bio-Technology, Shanghai, China) as quality control of DNA (Ali et al., 2018b). The PCR amplifications were carried out in a total volume of 25 μl containing 50 ng template DNA, 2 μl of each forward and reverse primers, 12.5 μl of 2X Taq PCR Mastermix (Tiangen Biotechnology Beijing, China). Negative controls were included where the DNA was replaced by ddH2O. The amplification conditions were setup as follow: initial denaturation at 94°C for 3 min, 30 cycles each of denaturation at for 30 s at 94°C, annealing, at 55°C for 30 s, elongation at 72°C for 60 s and a final extension at 72°C for 7 min.
16S rRNA Amplification
The bacterial 16S rRNA gene V3-V4 region was amplified by PCR using gene specific primers 341F:CCTACGGGNGGCWGCAG; 806R:GGACTACHVGGGTATCTAAT. The PCR reactions were carried out in a total volume of 50 μl mixture containing 5 μl of 10 × KOD Buffer, 5 μl of 2.5 mM dNTPs, 1.5 μl of each primer (5 μM), 1 μl of KOD Polymerase, and 100 ng of template DNA. The thermal conditions set for PCR were as follow: initial denaturation at 95°C for 2 min, followed by 27 cycles at 98°C for 10 s, 62°C for 30 s, and 68°C for 30 s and a final extension at 68°C for 10 min. The amplified PCR products along with negative controls were mixed with the same volume of 1× buffer containing SYBR green and run on 2% agarose gel.
Library Preparation and Sequencing
Next, the PCR products of all samples were purified using AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, United States) according to the manufacturer’s instructions and quantified using QuantiFluor-ST (Promega, United States). Sequencing libraries were generated using Illumina NEB (New England Biolabs) Next® UltraTM DNA Library preparation Kit and sequenced using HiSeq 2500 platform after quality assessment on the Qubit@ 2.0 Fluorometer (Thermo Scientific) and Agilent Bio-analyzer 2100 system.
Paired-End (PE) Reads Assembly and Quality Control
The paired-end reads were merged through FLASH (V1.2.7) (Magoc and Salzberg, 2011). The raw tags were passed through quality filtering to obtain high-quality clean tags using QIIME (Quantitative Insight into Microbial Ecology) (V1.7.0) (Caporaso et al., 2010; Bokulich et al., 2013) and were compared with reference database (Gold database) by applying UCHIME algorithm. The effective tags used for OTUs analysis were obtained by removing the chimeric sequences, primer sequences, and barcode sequences (Edgar et al., 2011).
OTUs (Operational Taxonomic Units) Cluster and Species Annotation
Operational taxonomic units were clustered with Uprase (V 7.0.1001)1 at ≥97% (3% cut-off) similarity threshold (Edgar, 2013). Based on RDP classifier algorithm (Version 2.2)2 (Wang et al., 2014), representative sequences were annotated for taxonomic information using GreenGene Database3 according to the online protocol. Multiple sequence alignment was performed with MUSCLE (Version 3.8.31) for phylogenetic analysis (Edgar, 2004). Between the groups, Venn analysis was performed in R to identify unique and common OTUs in R (Project R 3.0.2) (Caporaso et al., 2010).
Alpha (α) and Beta (β) Diversity Analysis
To investigate the bacterial species richness and community diversity in the samples, we calculated six α diversity indices (Shannon, Simpson, Chao 1, ACE, Observed species, and Good coverage) using the Vegan package in R (Project R 3.0.24) (Dixon, 2003; Caporaso et al., 2010) as explained by Montagna et al. (2015a,b, 2016). To investigate the pattern of variation in the structure of microbiota across different life stages, we employed β diversity weighted and un-weighted Unifrac matrices using QIIME (Version 1.7.0) (Caporaso et al., 2010). Furthermore, ANOSIM (analysis of similarity) analysis was carried out to reveal the differences in the structure of microbiota across the life stages (Montagna et al., 2016).
Putative Functional Profiling
The putative functional profiling was carried out using Tax4Fun as described by Aßhauer et al. (2015). The functional community profiling was predicted based on the bacterial16S rRNA gene OTUs associated with different life stages and reproductive organs. The sequenced prokaryotic genomes of 16S rRNA gene sequences were linked to KEGG (Kyoto Encyclopedia of Genes and Genomes) orthology for functional annotating (Kanehisa et al., 2012).
Statistical Analysis
Statistical differences between the two groups were measured using independent t test and multiple comparison was done with Tukey’s HSD post hoc analysis. Further explanation on statistical analysis for alpha and beta diversity indices can be found in the references provided. Wolbachia infection density dynamics across the life stages was estimated using Oneway ANOVA, while between the male and female or their reproductive organs, Independent sample t-test was used to establish the statistical differences. The p value was set <0.05 and the analysis was done with IBM SPSS Statistics (V.22.0).
Results
Bacterial Diversity Estimation
High-throughput sequencing analysis yielded a total of 2,075,281 raw reads from the 21 samples of various developmental stages [eggs, larvae, pupae, and adults (male and female)] and reproductive organs (Ovaries and Testis) of O. nipae. After quality filtering with QIIME and removal of chimeric sequences using UCHIME algorithm a total of 1,971,022 high quality clean tags were obtained for subsequent analysis (Supplementary Table S1). Clustering analysis of the pre-processed tags at 97% similarity threshold (3% cut off) generated a total of 3,959 OTUs (Table 1). The rarefaction analysis unraveled the sequencing depth and microbial diversity covered for all the samples, showed a degree of saturation (Supplementary Figure S2). The alpha diversity indices were estimated to uncover the bacterial diversity (Simpson and Shannon) and species richness (Chao1 and ACE) (Table 1). Diversity analysis revealed significant differences in microbial communities associated with O. nipae (Figure 1A and Supplementary Table S2). Across the various developmental stages the species richness estimated by Chao1 and ACE indicated significant differences between the microbiota of adult and pupal stages (Chao1 p = 0.034, ACE p = 0.037). However, no significant differences were observed in bacterial species richness and diversity between the two sexes (male and female, p > 0.05) and reproductive organs (ovaries and testis, p > 0.05). Similarly, the Beta diversity of microbiota associated with different life stages of this beetle investigated by weighted and un-weighted Unifrac matrices explained significant differences (Weighted Unifrac p = 0.00007, Un-weighted Unifrac p = 0.004), however, it did not vary in the reproductive organs. Also, in principal component analysis (PCA), the first two components PC1 and PC2 accounted for 91.9% (PC1 = 76.3%, PC2 = 15.6%) of the variation in bacterial symbionts (Figure 1B). Moreover, the ANOSIM analysis revealed high OTU turnover and low nestedness (R = 0.69, p = 0.001) across the different life stages which further confirmed that structure of microbiota was different (Supplementary Figure S3). Overall, diversity estimation revealed higher bacterial species richness and community diversity in the eggs followed by adult female and male and their respective reproductive organs and the lowest in the pupal and larval stages, respectively.
Table 1. Operational taxonomic units (OTUs), species richness and diversity estimation of the bacterial communities associated with various developmental stages and reproductive organs of Octodonta nipae.
Figure 1. Comparison of bacterial diversity (A) Alpha (α) and (B) Beta (β) in the samples from different life stages and reproductive organs of Octodonta nipae. The student t test compared the alpha diversity between the two groups while Tukey’s HSD multiple comparisons was used for comparing the differences in diversity across the life stages (see Supplementary Table S2). The line inside the boxplots represent medians, the dots in the center are the means and the whiskered bars are maximal and minimal values. Principal component analysis (PCA) comparing the bacterial symbionts between different life stages and reproductive organs. The Eigen values were calculated based on Bray Curtis distance matrix at 97% similarity threshold.
An Overview of O. nipae Bacterial Symbionts
Regardless of the developmental stages, sex or reproductive organs, the bacterial symbionts of O. nipae were mainly consist of the phyla Proteobacteria, Actinobacteria, Firmicutes, Cyanobacteria, and Bacteroidetes (>1%). Other least abundant taxa belonged to phyla Fusobacteria, Chloroflexi, Saccharibacteria, and Planctomycetes (<1%). Within the Proteobacteria, the microbiota was dominated by the families Anaplasmataceae, Enterobacteriaceae, Alcaligenaceae, Pseudomonadaceae, and Rhodocyclaceae. Family Dermabacteraceae, Tsukamurellaceae, Pseudonocardiaceae, and Actinomycetaceae were represented by the phylum Actinobacteria while the phylum Firmicutes was mainly comprised of bacteria in the family Streptococcaceae (relative abundance >1%) (Figure 2A). At least 70 bacterial genera were detected in our samples, among which the top 10 genera were; Wolbachia, Pantoea, Brachybacterium, Serratia, Lactococcus, Advenella, Pseudomonas, Tsukamurella, Pseudonocardia, and Acinetobacter, together accounted for 63–78% of the entire bacteria at genus level (Figure 2B).
Figure 2. Bacterial community composition at (A) family and (B) genus level associated with various developmental stages and reproductive organs of O. nipae.
Ontogenetic Changes of Microbiota Associated With O. nipae
Assessment of the complexity and variability of microbial communities across various developmental stages (egg, larvae, pupae, and adult) exhibited obvious differences. Two major phyla, Proteobacteria (ranging from 65 to 82%) and Actinobacteria (ranging from 11 to 26%) represented the bulk of taxa from all developmental stages (ranging from 79 to 99%). The relative abundance of other phyla varied markedly particularly the phylum Firmicutes in the eggs and pupal stages, Cyanobacteria (ranging from 0.4 to 2.8%) and Bacteroidetes (ranging from 0.006 to 2.1%) in all life stages. Almost 97% of the total read were successfully assigned to more than 50 bacterial families. Among them, the top ten families include; Anaplasmataceae, Enterobacteriaceae, Dermabacteraceae, Alcaligenaceae, Streptococcaceae, Pseudomonadaceae, Tsukamurellaceae, Rhodocyclaceae, Pseudonocardiaceae, and Actinomycetaceae accounted for 82.1–95.59% of the total taxa with distinct changes in relative abundances (Figure 2A). Enterobacteriaceae was the most abundant family during the early developmental stages (eggs 45.14%, larvae 69.21%, and pupae 49.79%), however, family Anaplasmataceae dominated the adult stage (male 51.29%, female 58.73%). Also, abundance of the phylum Dermabacteraceae and Alcaligenaceae was significantly higher in the pupal stage. At least 70 bacterial genera were detected in O. nipae, among which the top 10 genera were; Wolbachia, Pantoea, Brachybacterium, Serratia, Lactococcus, Advenella, Pseudomonas, Tsukamurella, Pseudonocardia, and Acinetobacter accounted for 63–78% of the bacteria at genus level with marked variation in relative abundance (Figure 2B). Serratia (23.53%), Lactococcus (14.88%), and Pseudomonas (7.44%) were represented the most in the eggs followed by Wolbachia (3.83%), Brachybacterium (3.67%), Acinetobacter (2.75%) Tsukamurella (2.73%), Pantoea (2.1%), and Advenella (1.99%). Conversely, the genus Pantoea represented the most in the larval (41.89%) and pupal (29.81%) stages, respectively. In pupal stage, Brachybacterium prevalence (24.71%) outnumbered as compared with other life stages (1.22–11.46%). Interestingly, the relative abundance of Wolbachia dramatically increased in the adult stage (male 64.74%, female 64.88%) (Figure 2B). Other OTUs with low abundances were Pseudonocardia, Acinetobacter, Actinomyces, Brevibacterium, Escherichia, Shigella, Amaricoccus, Corynebacterium-1, Proteiniphilum, Sphingobacterium, Stenotrophomonas, Paracoccus, Apibacter, Mycobacterium, Segniliparus, Terrisporobacter, and Sphingomonas. Taken together, the structure and composition of microbiota varied greatly at different developmental stages (Supplementary Figure S4).
Microbiota Associated With the Males and Females of O. nipae
Next, we compared the microbiota associated with the adult male and female to determine the impact of sex on microbial profile of O. nipae. Together, the two sexes shared a substantial 223 OTUs and the species richness and community diversity estimated by alpha and beta diversity analysis indicated that bacterial communities remained fairly unchanged (Supplementary Table S2). Moreover, the relative abundance of the bacterial taxa also didn’t differ significantly suggesting that gender doesn’t influence the bacterial symbionts associated with O. nipae (Supplementary Table S2). In adults, the most abundant OTU was a Wolbachia accounted for more than 64% of the total microbiota. Nonetheless, other prominent OTUs represented by Tsukamurella (male = 5.30%, female = 3.63%), Pseudonocardia (male = 5.11%, female = 3.13%), (male = 5.30%, female = 3.63%), Pantoea (male = 5.30%, female = 3.63%), Advenella (male = 2.43%, female = 1.68%), Brachybacterium (male = 2.38%, female = 1.22%), Serratia (male = 1.49%, female = 2.79%), Pseudomonas (male = 1.49%, female = 2.79%), and Acinetobacter (male = 0.64%, female = 0.18%) colonized O. nipae male and female, respectively.
Microbiota Associated the Reproductive Organs of O. nipae
Further interrogation on the bacterial composition and diversity between the ovaries and testis of O. nipae revealed that the microbial profile of reproductive organs was conserved in term of species richness and community diversity. A substantial conservation of adult microbiota was reflected in the reproductive organs of male and female dominated by Wolbachia (Ovaries = 68.01%, Testis = 74.94%). Also, abundances of the genera Lactococcus (Ovaries = 6.88%, Testis = 0.04%), Pseudomonas (Ovaries = 3.14%, Testis = 1.47%), and Serratia (Ovaries = 1.9 %, Testis = 0.5%) were highly variable between the reproductive organs of male and female O. nipae adults. Overall, analysis of PCR amplified bacterial 16S rRNA gene fragments from the different developmental stages of O. nipae revealed relatively complex bacterial communities that differed according to the developmental stage, however, the species richness and diversity between the sexes and reproductive organs were relatively equal except for the relative abundances of the dominant taxon (Figure 2).
Transovarial Transmission and Dynamics of Wolbachia
The most prevalent OTU belonged to Wolbachia being one of the most abundant genera particularly in the adult stage and reproductive organs of O. nipae accounted for more than 68% of the total microbiota at the genus level (Figure 3). Previously we have confirmed Wolbachia infection at different life stages and reproductive organs of O. nipae using wsp (Wolbachia surface protein) gene based diagnostic PCR approach (Ali et al., 2018b). Considering this, we assumed the similar transmission mechanism of Wolbachia and evaluated its presence. Diagnostic PCR assays and highthroughput sequencing analysis confirmed the ubiquitous presence of Wolbachia in all tested life stages and reproductive tissues (Ali et al., 2018b; Figure 3). These results suggest that Wolbachia is a transovarially transmitted obligate endosymbiont throughout the developmental stages of O. nipae. The relative abundance of Wolbachia across the life stages differed significantly, assessed by one-way ANOVA and variations within means were compared using Tukey’s HSD (Honest significant difference) test at p < 0.05 [ANOVA; F(4,10) = 689.82, p < 0.001]. In life stages, the highest abundance was detected in the adult stage; however, in immature stages (eggs, larvae, and pupae) the relative abundance did not vary significantly (Figure 3). Furthermore, comparing the Wolbachia prevalence between male and female individuals, the results showed that the infection density is similar and not affected by the sex (t = -0.074, p = 0.94). Although, the reproductive organs were favorably infected by Wolbachia (Ovaries = 68.01%, Testis = 74.9%) the infection density remained the same (t = 1.65, p = 0.173) between male and female (Figure 3).
Figure 3. Dynamics of Wolbachia infection relative abundance at different life stages and reproductive organs of O. nipae. Alphabetic letters indicate statistical differences across the life stages is estimated by one-way ANOVA with Tukey’s HSD test for multiple comparisons [F(4,10) = 689.82, p < 0.001]. Between the two sexes (male and female) and reproductive organs (ovaries and testis), Wolbachia density is compared by Independent sample t-test. NS, non-significant p > 0.05.
Inference of Putative Functions of O. nipae Microbiota
The metabolic potential of microbiota associated with O. nipae was inferred based on the 16S rRNA gene sequence of the prokaryotic genome of the existing KEGG database. The results predicted important putative functions played by the host’s bacterial symbionts crucial for the host physiology and development. Regardless of the developmental stage at least 275 OTUs were found to be involved in various functional categories including environmental information processing (signal transduction), genetic information processing, nucleotide metabolism (purine metabolism, pyrimidine metabolism), energy metabolism, nitrogen metabolism, cellular processes (cell growth and death) amino acid metabolism (Arginine and proline metabolism), carbohydrate metabolism (amino sugar and nucleotide sugar metabolism), metabolism of terpenoids and polyketides, and other processes. The metabolic potential was in congruent with the observed differences in the diversity and composition of the microbiota at different life stages where distinct shifts were observed in the putative functions (Figure 4). It is of prime importance to investigate taxa specific function to unravel the complex interplay between O. nipae and its symbionts.
Figure 4. Heat-map showing the metabolic potential inferred from the bacterial 16S rRNA gene sequences associated with different life stages. The different color pattern indicates the relative abundance of bacterial operational taxonomic units (OTUs) involved in various biological functions.
Discussion
Like other holometabolous insects, O. nipae undergoes complete metamorphosis (egg, larvae, pupae, and adult; Supplementary Figure S1) where the different stages are indeed differentiated in shape, structure and function. The present study provides the first insight into the bacterial symbionts of O. nipae across different life stages (eggs, larvae, pupae, and adult), two sexes (male and female) and reproductive organs (ovaries and testis). Diversity estimation analysis revealed that the bacterial communities differed significantly across the life stages, however, minimal influence of the sex was noticed on microbiota diversity between male and female individuals. Also, the bacterial communities in the reproductive organs (ovaries and testis) were conserved in term of species richness and diversity, however, marked differences were observed in the relative abundances of the major OTUs. The composition of bacterial community of O. nipae was relatively complex and diverse, yielded 149–269 OTUs per sample (Total 3,959 OTUs) (Table 1). Overall, diversity estimation revealed higher bacterial species richness and community diversity in the eggs followed by adult female and male and their respective reproductive organs and the lowest in the pupal and larval stages, respectively (Figure 1 and Table 1). A number of factors have been identified to affect the structure and composition of microbiota including host genetics, age, sex, diet, and geography (Adams et al., 2010; Elijah Powell et al., 2018; Mariño et al., 2018). The higher number of OTUs found in the eggs is comparable with the microbial diversity of Spodoptera littoralis, which might have caused due to exposure to a wide range of environmental microorganisms (Meng et al., 2016). Reduction in bacterial diversity in the proceeding developmental stages despite higher diversity in the eggs points out the control of host over its microbial titer (Chen et al., 2016). These results suggested that the O. nipae microbiota profile is more influenced by the developmental stage, whereas the impact of gender and reproductive tissues to showcase the dissimilarity among the bacterial communities was not that pronounced.
Like many other insects, O. nipae bacterial symbionts were mainly dominated by bacteria in the phyla Proteobacteria and Actinobacteria followed by Firmicutes, Cyanobacteria, and Bacteroidetes (>1%). Despite the significant variation, a substantial conservation of microbiota was observed across all life stages which indicates that the major microbiota are consistently present and may have a confounding impact on host physiology and invasion success (Dillon and Dillon, 2004; Muhammad et al., 2017). Proteobacteria frequently dominate the microbiota in various invertebrates including red palm weevil, R. ferrugineus (Muhammad et al., 2017), ground beetles, Harpalus pensylvanicus and Anisodactylus sanctaecruis (Schmid et al., 2014), desert locust, Schistocerca gregaria (Dillon et al., 2010), Bacttocera dorsalis (Andongma et al., 2015) and Hypothenemus hampei (Mariño et al., 2018). Actinobacteria was the second major bacterial component in O. nipae which have been shown to play various metabolic and physiological functions including synthesis of extracellular enzymes and secondary metabolites (Schrempf, 2001). Similarly, the phylum Firmicutes represent a major component of B. minax (Wang et al., 2014), R. ferrugineus (Muhammad et al., 2017), Mexican fruit fly, Anastrepha ludens (Kuzina et al., 2001) and the ant, Solenopsis invicta (Ishak et al., 2011). Enterobacteriaceae was the most abundant family during the early developmental stages (eggs 45.14%, larvae 69.21%, and pupae 49.79%), however, family Anaplasmataceae dominated the adult stage (male 51.29%, female 58.73%). Also, abundance of the family Dermabacteraceae and Alcaligenaceae abundance was significantly higher in the pupal stage. During metamorphosis, in the transition of life stages the structure of microbiota changes drastically (Tchioffo et al., 2016). In line with our observations numerous studies have shown diverse microbiota colonizes many insects and varies across the life stages. In B. dorsalis, the bacterial communities dominated by Proteobacteria and Firmicutes differed at different developmental stages (Andongma et al., 2015). A recent report on microbial profile of another beetle, H. hampei whose microbiota is highly diverse and variable displays a high degree of similarity in bacterial structure and composition with our study insect (Mariño et al., 2018). On contrary, the bacterial community structure in life stages and two sexes of the Chinese white pine beetle, Dendroctonus armandi exhibited no significant differences (Hu et al., 2013). Studies have suggested that the apparent discrepancies between male and female microbiotas profile is the result of different nutritional requirements and energy allocation (Basset et al., 2012; Ridley et al., 2012). Moreover, the microbiota structure of the reproductive organs of male and female O. nipae were almost the same in richness and diversity but varied in relative abundance (Figure 3). In line with this notion, the bacterial communities colonized the reproductive organs of Bactocera minax were similar between the male and female, dominated by phyla Proteobacteria, Firmicutes, and Actinobacteria (Wang et al., 2014).
In nature, Wolbachia infection is widespread and a large number of insects including many coleopteran live in symbiotic relationship with this endosymbiont (Fialho and Stevens, 2000; Hurst et al., 2000; Duron et al., 2008; Ali et al., 2016, 2018a,b,c). The abundance of Wolbachia differed significantly across the life stages and was the highest in the adult stage. The exact trend of Wolbachia infection have been witnessed in this beetle and its closest sister species namely Brontispa longissima using RT-qPCR (Ali et al., 2018c). Variation in relative abundance across life stages was expected because symbionts differentially thrive in the proceeding life stages shaped by host nutrition and functionality. The highest abundance in the adult stage and reproductive organs is concordant with its ability to manipulate host reproduction (Dobson et al., 1999). Wolbachia is a well-known reproductive parasites in various insects (Duron et al., 2008). It has been shown to manipulate host’s reproduction in several ways including cytoplasmic incompatibility (Field et al., 1999), parthenogenesis (Pannebakker et al., 2004), feminization (Bouchon et al., 1998), and male killing (Fialho and Stevens, 2000; Hurst et al., 2000). Additionally, Wolbachia infection density can be higher in both adult and egg, as explained for Wolbachia social transmission pattern being maternally transmitted with egg cytoplasm to the next generations that eventually regulate host’s reproduction (Ali et al., 2018c). Other than its phenomenal role in reproductive manipulation, Wolbachia also play a crucial role in host’s evolutionary success by inducing rapid speciation (Zhou et al., 1998; Shoemaker et al., 2004) or extinction of species (Dobson et al., 2002). Similarly, in parasitic wasp, Asobara tabida, Wolbachia is necessary for eggs development and oogenesis (Dedeine et al., 2001). Wolbachia infection in non-reproductive body parts such as nervous and muscular tissues has been documented in the flies (Drosophila), where it exerts a negative impact on host’s longevity (Juchault et al., 1994). Despite the enormous potential of Wolbachia as a reproductive manipulator and a prospective bio-control agent, its interactions with this insect pest of economic importance remains unknown. Taken together, these results along with our previous findings confirmed the ubiquitous presences of Wolbachia at all life stages, body parts and tissues of O. nipae and suggested that it is a transovarially transmitted highly adapted obligate endosymbiont worth further studies to unravel its biological function and impact on host physiology.
According to the annotation of KEGG ontology (Kanehisa et al., 2012) analysis the metabolic potential of bacterial communities are predicted. Many OTUs regardless of their affiliation with any life stage demonstrated important metabolic potential (Figure 4), suggested that O. nipae microbiota plays critical role in host physiology. For example, the most abundant OTUs in the immature stages belonged to family Enterobacteriaceae, is in line with several other studies (Andongma et al., 2015; Montagna et al., 2015a; Muhammad et al., 2017; Zhukova et al., 2017) which was significantly decreased in the adults and replaced by Anaplasmataceae as the most abundant family. Members of this family are involved in many important functions such as nitrogen fixation to make it available to its host (Dixon and Kahn, 2004; Behar et al., 2005), provides protection against the pathogen and enhances host fitness (Dillon and Dillon, 2004), degradation of uric acid by the enzyme urease to transform it to ammonia (Lauzon et al., 2000). In R. ferrugineus, Serratia have exhibited antimicrobial activity against bacterial pathogen and could be a potential bio-control agent of palm pests (Scrascia et al., 2016). In yellow-spotted longicorn beetle, Lactococcus sp. is involved in making of lactic acid and digest polysaccharides (Mazza et al., 2014). Acinetobacter degrades pesticides and other large molecular compounds such as polychlorinated phenols and polycyclic aromatic hydrocarbons for its insect hosts (Hao et al., 2002). Pseudomonas isolated from Plutella xylostella can produce antifungal siderophore compounds like pyoverdine, however the Brachybacterium was unable to siderophore synthesis (Indiragandhi et al., 2007). Also, Pantoea from grasshopper possesses antagonistic activity and inhibit the growth and spore formation of fungi Metarhizium anisopliae (Dillon and Charnley, 1995). Bacterium Tsukamurella is a human opportunistic pathogen (Safaei et al., 2018), however, its function is largely unknown in insects.
Conclusion
For the first time we brought to light the bacterial symbionts of O. nipae and quantified the dynamics of bacterial communities in life stages, sexes, and reproductive organs. The results revealed that bacterial diversity varied significantly across the life stages but not between the two sexes or reproductive organs. The bacterial composition and its relative abundances also varied greatly according to the life stages, indicating its influence on microbiota structure. A substantial conservation of the major taxa was observed across all life stages and reproductive organs, however, differed significantly in term of its relative abundance. Changes in the microbiota structure during metamorphosis was also reflected on its metabolic activities indicating the functional significance of O. nipae microbiota. Maternally inherited Wolbachia is proven to be an obligate endosymbiont of this insects which has enormous potential to manipulate host’s reproductive and behavioral responses. Finally, our study indicates the significance of characterizing O. nipae bacterial symbionts which will provide the basis for subsequent research to unravel the cross talk between the host and its bacterial associates.
Author Contributions
HA conducted the sampling and performed the experiments. AM, YH, NS, and HA analyzed the data and wrote the manuscript. YMH supervised and critically reviewed the manuscript.
Funding
We are very grateful to the grant from the National Key R&D Program of China (2017YFC1200605) and Fujian Science and Technology Special Project (2017NZ0003-1-6).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.00466/full#supplementary-material
Footnotes
- ^ http://drive5.com/uparse
- ^ http://sourceforge.net/projects/rdp-classifier/
- ^ http://greengenes.lbl.gov/
- ^ http://cran.r-project.org/
References
Adams, A. S., Adams, S. M., Currie, C. R., Gillette, N. E., and Raffa, K. F. (2010). Geographic variation in bacterial communities associated with the red turpentine beetle (Coleoptera: Curculionidae). Environ. Entomol. 39, 406–414. doi: 10.1603/EN09221
Ali, H., Hou, Y., Tang, B., Huang, B., Muhammad, A., Sanda, N. B., et al. (2016). A way of reproductive manipulation and biology of Wolbachia pipientis. J. Exp. Biol. Agric. Sci. 4, 156–168. doi: 10.18006/2016.4(2).156.168
Ali, H., Muhammad, A., Bala, N. S., Wang, G., Chen, Z., Peng, Z., et al. (2018a). Genomic evaluations of Wolbachia and mtDNA in the population of coconut hispine beetle, Brontispa longissima (Coleoptera: Chrysomelidae). Mol. Phylogenet. Evol. 127, 1000–1009. doi: 10.1016/j.ympev.2018.07.003
Ali, H., Muhammad, A., and Hou, Y. (2018b). Infection density dynamics and phylogeny of Wolbachia associated with coconut hispine beetle, Brontispa longissima (Gestro) (Coleoptera: Chrysomelidae) by multilocus sequence type (MLST) genotyping. J. Microbiol. Biotechnol. 28, 796–808. doi: 10.4014/jmb.1712.12019
Ali, H., Muhammad, A., Islam, S. U., Islam, W., and Hou, Y. (2018c). A novel bacterial symbiont association in the hispid beetle, Octodonta nipae (Coleoptera: Chrysomelidae), their dynamics and phylogeny. Microb. Pathog. 118, 378–386. doi: 10.1016/j.micpath.2018.03.046
Andongma, A. A., Wan, L., Dong, Y. C., Li, P., Desneux, N., White, J. A., et al. (2015). Pyrosequencing reveals a shift in symbiotic bacteria populations across life stages of Bactrocera dorsalis. Sci. Rep. 5:9470. doi: 10.1038/srep09470
Aßhauer, K. P., Wemheuer, B., Daniel, R., and Meinicke, P. (2015). Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data. Bioinformatics 31, 2882–2884. doi: 10.1093/bioinformatics/btv287
Bansal, R., Michel, A. P., and Sabree, Z. L. (2014). The crypt-dwelling primary bacterial symbiont of the polyphagous pentatomid pest Halyomorpha halys (Hemiptera: pentatomidae). Environ. Entomol. 43, 617–625. doi: 10.1603/EN13341
Basset, Y., Cizek, L., Cuenoud, P., Didham, R. K., Guilhaumon, F., Missa, O., et al. (2012). Arthropod diversity in a tropical forest. Science 338, 1481–1484. doi: 10.1126/science.1226727
Behar, A., Yuval, B., and Jurkevitch, E. (2005). Enterobacteria-mediated nitrogen fixation in natural populations of the fruit fly Ceratitis capitata. Mol. Ecol. 14, 2637–2643. doi: 10.1111/j.1365-294X.2005.02615.x
Behar, A., Yuval, B., and Jurkevitch, E. (2008). Gut bacterial communities in the mediterranean fruit fly (Ceratitis capitata) and their impact on host longevity. J. Insect Physiol. 54, 1377–1383. doi: 10.1016/j.jinsphys.2008.07.011
Blatch, S. A., Meyer, K. W., and Harrison, J. F. (2010). Effects of dietary folic acid level and symbiotic folate production on fitness and development in the fruit fly Drosophila melanogaster. Fly 4, 312–319. doi: 10.4161/fly.4.4.13258
Bokulich, N. A., Subramanian, S., Faith, J. J., Gevers, D., Gordon, J. I., Knight, R., et al. (2013). Quality-filtering vastly improves diversity estimates from illumina amplicon sequencing. Nat. Methods 10:57. doi: 10.1038/nmeth.2276
Bouchon, D., Rigaud, T., and Juchault, P. (1998). Evidence for widespread Wolbachia infection in isopod crustaceans: molecular identification and host feminization. Proc. R. Soc. B Biol. Sci. 265, 1081–1090. doi: 10.1098/rspb.1998.0402
Breznak, J. A. (2000). “Ecology of prokaryotic microbes in the guts of wood and litter-feeding termites,” in Termites: Evolution, Sociality, Symbioses, Ecology, eds T. Abe, D. E. Bignell, and M. Higashi (Berlin: Springer Science & Business Media), 209–231. doi: 10.1007/978-94-017-3223-9_10
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows high throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth0510-335
Chaston, J. M., Newell, P. D., and Douglas, A. E. (2014). Metagenome-wide association of microbial determinants of host phenotype in Drosophila melanogaster. mBio 5, e01631–e01714. doi: 10.1128/mBio.01631-14
Chen, B., Teh, B. S., Sun, C., Hu, S., Lu, X., Boland, W., et al. (2016). Biodiversity and activity of the gut microbiota across the life history of the insect herbivore Spodoptera littoralis. Sci. Rep. 6:29505. doi: 10.1038/srep29505
Crotti, E., Balloi, A., Hamdi, C., Sansonno, L., Marzorati, M., Gonella, E., et al. (2012). Microbial symbionts: a resource for the management of insect-related problems. Microb. Biotechnol. 5, 307–317. doi: 10.1111/j.1751-7915.2011.00312.x
Dedeine, F., Vavre, F., Fleury, F., Loppin, B., Hochberg, M. E., and Bouletreau, M. (2001). Removing symbiotic Wolbachia bacteria specifically inhibits oogenesis in a parasitic wasp. Proc. Natl. Acad. Sci. 98, 6247–6252. doi: 10.1073/pnas.101304298
Dillon, R. J., and Charnley, A. K. (1995). Chemical barriers to gut infection in the desert locust: in vivo production of antimicrobial phenols associated with the bacterium Pantoea agglomerans. J. Invertebr. Pathol. 66, 72–75. doi: 10.1006/jipa.1995.1063
Dillon, R. J., and Dillon, V. M. (2004). The gut bacteria of insects: nonpathogenic interactions. Annu. Rev. Entomol. 49, 71–92. doi: 10.1146/annurev.ento.49.061802.123416
Dillon, R. J., Webster, G., Weightman, A. J., and Keith Charnley, A. (2010). Diversity of gut microbiota increases with aging and starvation in the desert locust. Antonie Van Leeuwenhoek 97, 69–77. doi: 10.1007/s10482-009-9389-5
Dixon, P. (2003). VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930. doi: 10.1111/j.1654-1103.2003.tb02228.x
Dixon, R., and Kahn, D. (2004). Genetic regulation of biological nitrogen fixation. Nat. Rev. Microbiol. 2, 621–631. doi: 10.1038/nrmicro954
Dobson, S. L., Bourtzis, K., Braig, H. R., Jones, B. F., Zhou, W., Rousset, F., et al. (1999). Wolbachia infections are distributed throughout insect somatic and germ line tissues. Insect Biochem. Mol. Biol. 29, 153–160. doi: 10.1016/S0965-1748(98)00119-2
Dobson, S. L., Fox, C. W., and Jiggins, F. M. (2002). The effect of Wolbachia-induced cytoplasmic incompatibility on host population size in natural and manipulated systems. Proc. R. Soc. London B Biol. Sci. 269, 437–445. doi: 10.1098/rspb.2001.1876
Douglas, A. E. (2009). The microbial dimension in insect nutritional ecology. Funct. Ecol. 23, 38–47. doi: 10.1111/j.1365-2435.2008.01442.x
Douglas, A. E. (2015). Multiorganismal insects: diversity and function of resident microorganisms. Annu. Rev. Entomol. 60, 17–34. doi: 10.1146/annurev-ento-010814-020822
Duron, O., Bouchon, D., Boutin, S., Bellamy, L., Zhou, L., Engelstädter, J., et al. (2008). The diversity of reproductive parasites among arthropods: Wolbachia do not walk alone. BMC Biol. 6:27. doi: 10.1186/1741-7007-6-27
Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi: 10.1093/nar/gkh340
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10:996. doi: 10.1038/nmeth.2604
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Elijah Powell, J., Eiri, D., Moran, N. A., and Rangel, J. (2018). Modulation of the honey bee queen microbiota: effects of early social contact. PLoS One 13:e0200527. doi: 10.1371/journal.pone.0200527
Fialho, R. F., and Stevens, L. (2000). Male-killing Wolbachia in a flour beetle. Proc. R. Soc. B Biol. Sci. 267, 1469–1473. doi: 10.1098/rspb.2000.1166
Field, L. M., James, A. A., Turelli, M., and Hoffmann, A. A. (1999). Microbe-induced cytoplasmic incompatibility as a mechanism for introducing transgenes into arthropod populations. Insect Mol. Biol. 8, 243–255. doi: 10.1046/j.1365-2583.1999.820243.x
Genta, F. A., Dillon, R. J., Terra, W. R., and Ferreira, C. (2006). Potential role for gut microbiota in cell wall digestion and glucoside detoxification in Tenebrio molitor larvae. J. Insect Physiol. 52, 593–601. doi: 10.1016/j.jinsphys.2006.02.007
Hao, O. J., Kim, M. H., Seagren, E. A., and Kim, H. (2002). Kinetics of phenol and chlorophenol utilization by Acinetobacter species. Chemosphere 46, 797–807. doi: 10.1016/S0045-6535(01)00182-5
Hou, Y., Miao, Y., and Zhang, Z. (2014). Leaf consumption capacity and damage projection of Octodonta nipae (Coleoptera: Chrysomelidae) on three palm species. Ann. Entomol. Soc. Am. 107, 1010–1017. doi: 10.1603/AN14061
Hou, Y., and Weng, Z. (2010). Temperature-dependent development and life table parameters of Octodonta nipae (Coleoptera: Chrysomelidae). Environ. Entomol. 39, 1676–1684. doi: 10.1603/EN10015
Hou, Y. M., Wu, Z. J., Wang, C. F., Xie, L. H., You, M. S., and Hou, Y. M. (2011). in The Status and Harm of Invasive Insects in Fujian, China. Biological Invasions: Problems and Countermeasures, eds L. H. Xie, M. S. You, and Y. M. Hou (Beijing: Science Press), 111–114.
Hu, X., Wang, C., Chen, H., and Ma, J. (2013). Differences in the structure of the gut bacteria communities in development stages of the chinese white pine beetle (Dendroctonus armandi). Int. J. Mol. Sci. 14, 21006–21020. doi: 10.3390/ijms141021006
Hua, R., Hou, Y., and Shi, Z. (2014). Changes in the contents of physiologically active substances in Octodonta nipae (Coleoptera: Chrysomelidae) after low temperature acclimation. Acta Entomol. Sin. 57, 265–273.
Hurst, G. D., Johnson, A. P., Schulenburg, J. H., and Fuyama, Y. (2000). Male-killing Wolbachia in Drosophila: a temperature-sensitive trait with a threshold bacterial density. Genetics 156, 699–709.
Indiragandhi, P., Anandham, R., Madhaiyan, M., Poonguzhali, S., Kim, G. H., Saravanan, V. S., et al. (2007). Cultivable bacteria associated with larval gut of prothiofos-resistant, prothiofos-susceptible and field-caught populations of diamondback moth, Plutella xylostella and their potential for, antagonism towards entomopathogenic fungi and host insect nutriti. J. Appl. Microbiol. 103, 2664–2675. doi: 10.1111/j.1365-2672.2007.03506.x
Ishak, H. D., Plowes, R., Sen, R., Kellner, K., Meyer, E., Estrada, D. A., et al. (2011). Bacterial diversity in Solenopsis invicta and Solenopsis geminata ant colonies characterized by 16S amplicon 454 pyrosequencing. Microb. Ecol. 61, 821–831. doi: 10.1007/s00248-010-9793-4
Jiu, M., Zhou, X. P., Tong, L., Xu, J., Yang, X., Wan, F. H., et al. (2007). Vector-virus mutualism accelerates population increase of an invasive whitefly. PLoS One 2:e182. doi: 10.1371/journal.pone.0000182
Jones, R. T., Sanchez, L. G., and Fierer, N. (2013). A cross-taxon analysis of insect-associated bacterial diversity. PLoS One 8:e61218. doi: 10.1371/journal.pone.0061218
Juchault, P., Frelon, M., Bouchon, D., and Rigaud, T. (1994). New evidence for feminizing bacteria in terrestrial isopods - evolutionary implications. Cr. Acad. Sci. Iii-Vie 317, 225–230. doi: 10.1017/CBO9781107415324.004
Kanehisa, M., Goto, S., Sato, Y., Furumichi, M., and Tanabe, M. (2012). KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40, D109–D114. doi: 10.1093/nar/gkr988
Koch, H., and Schmid-Hempel, P. (2011). Socially transmitted gut microbiota protect bumble bees against an intestinal parasite. Proc. Natl. Acad. Sci. 108, 19288–19292. doi: 10.1073/pnas.1110474108
Kuzina, L. V., Peloquin, J. J., Vacek, D. C., and Miller, T. A. (2001). Isolation and identification of bacteria associated with adult laboratory mexican fruit flies, Anastrepha ludens (Diptera: Tephritidae). Curr. Microbiol. 42, 290–294. doi: 10.1007/s002840110219
Kwong, W. K., Mancenido, A. L., and Moran, N. A. (2017). Immune system stimulation by the native gut microbiota of honey bees. R. Soc. Open Sci. 4:170003. doi: 10.1098/rsos.170003
Lauzon, C. R., Sjogren, R. E., and Prokopy, R. J. (2000). Enzymatic capabilities of bacteria associated with apple maggot flies: a postulated role in attraction. J. Chem. Ecol. 26, 953–967. doi: 10.1023/A:1005460225664
Li, J.-L. L., Tang, B.-Z. Z., Hou, Y.-M. M., and Xie, Y.-X. X. (2016). Molecular cloning and expression of the vitellogenin gene and its correlation with ovarian development in an invasive pest Octodonta nipae on two host plants. Bull. Entomol. Res. 106, 1–9. doi: 10.1017/S0007485316000353
Lu, M., Hulcr, J., and Sun, J. (2016). The role of symbiotic microbes in insect invasions. Annu. Rev. Ecol. Evol. Syst. 47, 487–505. doi: 10.1146/annurev-ecolsys-121415-132050
Magoc, T., and Salzberg, S. L. (2011). FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963. doi: 10.1093/bioinformatics/btr507
Mariño, Y. A., Ospina, O. E., Verle Rodrigues, J. C., and Bayman, P. (2018). High diversity and variability in the bacterial microbiota of the coffee berry borer (Coleoptera: Curculionidae), with emphasis on Wolbachia. J. Appl. Microbiol. 125, 528–543. doi: 10.1111/jam.13768
Marzorati, M., Wittebolle, L., Boon, N., Daffonchio, D., and Verstraete, W. (2008). How to get more out of molecular fingerprints: practical tools for microbial ecology. Environ. Microbiol. 10, 1571–1581. doi: 10.1111/j.1462-2920.2008.01572.x
Mazza, G., Chouaia, B., Lozzia, G. C., and Montagna, M. (2014). The bacterial community associated to an italian population of Psacothea hilaris: a preliminary study. Bull. Insectol. 67, 281–285.
McMeniman, C. J., Lane, R. V., Cass, B. N., Fong, A. W. C., Sidhu, M., Wang, Y. F., et al. (2009). Stable introduction of a life-shortening Wolbachia infection into the mosquito Aedes aegypti. Science 323, 141–144. doi: 10.1126/science.1165326
Meng, E., Tang, B., Hou, Y., Chen, X., Chen, J., and Yu, X.-Q. Q. (2016). Altered immune function of Octodonta nipae (maulik) to its pupal endoparasitoid, Tetrastichus brontispae Ferrière. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 198, 100–109. doi: 10.1016/j.cbpb.2016.04.001
Montagna, M., Chouaia, B., Mazza, G., Prosdocimi, E. M., Crotti, E., Mereghetti, V., et al. (2015a). Effects of the diet on the microbiota of the red palm weevil (Coleoptera: Dryophthoridae). PLoS One 10:e0117439. doi: 10.1371/journal.pone.0117439
Montagna, M., Gomez-Zurita, J., Giorgi, A., Epis, S., Lozzia, G., and Bandi, C. (2015b). Metamicrobiomics in herbivore beetles of the genus Cryptocephalus (Chrysomelidae): toward the understanding of ecological determinants in insect symbiosis. Insect Sci. 22, 340–352. doi: 10.1111/1744-7917.12143
Montagna, M., Mereghetti, V., Gargari, G., Guglielmetti, S., Faoro, F., Lozzia, G., et al. (2016). Evidence of a bacterial core in the stored products pest Plodia interpunctella: the influence of different diets. Environ. Microbiol. 18, 4961–4973. doi: 10.1111/1462-2920.13450
Moreira, L. A., Iturbe-Ormaetxe, I., Jeffery, J. A., Lu, G., Pyke, A. T., Hedges, L. M., et al. (2009). A Wolbachia symbiont in Aedes aegypti limits Infection with dengue, chikungunya, and plasmodium. Cell 139, 1268–1278. doi: 10.1016/j.cell.2009.11.042
Muhammad, A., Fang, Y., Hou, Y., and Shi, Z. (2017). The gut entomotype of red palm weevil Rhynchophorus ferrugineus Olivier (Coleoptera: Dryophthoridae) and their effect on host nutrition metabolism. Front. Microbiol. 8:2291. doi: 10.3389/fmicb.2017.02291
Osborne, S. E., San Leong, Y., O’Neill, S. L., Johnson, K. N., Leong, Y. S., O’Neill, S. L., et al. (2009). Variation in antiviral protection mediated by different Wolbachia strains in Drosophila simulans. PLoS Pathog. 5:e1000656. doi: 10.1371/journal.ppat.1000656
Pais, R., Lohs, C., Wu, Y., Wang, J., and Aksoy, S. (2008). The obligate mutualist Wigglesworthia glossinidia influences reproduction, digestion, and immunity processes of its host, the tsetse fly. Appl. Environ. Microbiol. 74, 5965–5974. doi: 10.1128/AEM.00741-08
Pannebakker, B. A., Pijnacker, L. P., Zwaan, B. J., and Beukeboom, L. W. (2004). Cytology of Wolbachia-induced parthenogenesis in Leptopilina clavipes (Hymenoptera: Figitidae). Genome 47, 299–303. doi: 10.1139/g03-137
Peng, L., Chen, L., Li, J., Hou, Y., and Chen, Y. (2018). Mate recognition and antennal morphology of Octodonta nipae (Coleoptera: Chrysomelidae) adults. J. Asia. Pac. Entomol. 21, 268–278. doi: 10.1016/j.aspen.2018.01.008
Price, P. W. P., Denno, R. R. F., Eubanks, M. D. M., Finke, D. D. L., and Kaplan, I. (2011). Insect Ecology: Behavior, Populations and Communities. Cambridge: Cambridge University. doi: 10.1017/CBO9780511975387
Ridley, E. V., Wong, A. C. N., Westmiller, S., and Douglas, A. E. (2012). Impact of the resident microbiota on the nutritional phenotype of Drosophila melanogaster. PLoS One 7:e36765. doi: 10.1371/journal.pone.0036765
Safaei, S., Fatahi-Bafghi, M., and Pouresmaeil, O. (2018). Role of Tsukamurella species in human infections: first literature review. New Microb. New Infect. 22, 6–12. doi: 10.1016/j.nmni.2017.10.002
Salem, H., Kreutzer, E., Sudakaran, S., and Kaltenpoth, M. (2013). Actinobacteria as essential symbionts in firebugs and cotton stainers (Hemiptera, Pyrrhocoridae). Environ. Microbiol. 15, 1956–1968. doi: 10.1111/1462-2920.12001
Schmid, R. B., Lehman, R. M., and Lundgren, J. G. (2014). Sex-specific interactions of microbial symbioses on cricket dietary selection. Environ. Entomol. 43, 896–902. doi: 10.1603/EN13311
Schrempf, H. (2001). Recognition and degradation of chitin by streptomycetes. Antonie van Leeuwenhoek 79, 285–289. doi: 10.1023/A:1012058205158
Scrascia, M., Pazzani, C., Valentini, F., Oliva, M., Russo, V., D’Addabbo, P., et al. (2016). Identification of pigmented Serratia marcescens symbiotically associated with Rhynchophorus ferrugineus olivier (Coleoptera: curculionidae). Microbiologyopen 5, 883–890. doi: 10.1002/mbo3.377
Shoemaker, D. D. W., Dyer, K. A., Ahrens, M., McAbee, K., and Jaenike, J. (2004). Decreased diversity but increased substitution rate in host mtDNA as a consequence of Wolbachia endosymbiont infection. Genetics 168, 2049–2058. doi: 10.1534/genetics.104.030890
Steiner, H. (2001). Rattan Und Insekten – Eine Vergleichende Studie Der Herbivorenfauna Von Calamus Mananund Anderen Palmenarten Auf Der Halbinsel Malaysia. PhD. Thesis. Frankfurt: Johann Wolfgang Goethe-Universität,Frankfurt am Main.
Sugio, A., Dubreuil, G., Giron, D., and Simon, J. C. (2015). Plant-insect interactions under bacterial influence: ecological implications and underlying mechanisms. J. Exp. Bot. 66, 467–478. doi: 10.1093/jxb/eru435
Sun, J., Yu, P., Zhang, Y., and Wang, X. (2003). A new invasive coconut pest in hainan province. Entomol. Knowl. 40, 286–287.
Tchioffo, M. T., Boissière, A., Abate, L., Nsango, S. E., Bayibéki, A. N., Awono-Ambéné, P. H., et al. (2016). Dynamics of bacterial community composition in the malaria mosquito’s epithelia. Front. Microbiol. 6:1500. doi: 10.3389/fmicb.2015.01500
Vassiliou, V. A., Kazantzis, E., and Melifronidou-Pantelidou, A. (2011). First report of the nipa palm hispid Octodonta nipae on queen palms in cyprus. Phytoparasitica 39, 51–54. doi: 10.1007/s12600-010-0128-2
Verstraete, W. (2007). Microbial ecology and environmental biotechnology. ISME J. 1, 4–8. doi: 10.1038/ismej.2007.7
Visser, A. A., Nobre, T., Currie, C. R., Aanen, D. K., and Poulsen, M. (2012). Exploring the potential for actinobacteria as defensive symbionts in fungus-growing termites. Microb. Ecol. 63, 975–985. doi: 10.1007/s00248-011-9987-4
Wang, A., Yao, Z., Zheng, W., and Zhang, H. (2014). Bacterial communities in the gut and reproductive organs of Bactrocera minax (Diptera: Tephritidae) based on 454 pyrosequencing. PLoS One 9:e106988. doi: 10.1371/journal.pone.0106988
Warnecke, F., Luginbühl, P., Ivanova, N., Ghassemian, M., Richardson, T. H., Stege, J. T., et al. (2007). Metagenomic and functional analysis of hindgut microbiota of a wood-feeding higher termite. Nature 450:560. doi: 10.1038/nature06269
Zabalou, S., Riegler, M., Theodorakopoulou, M., Stauffer, C., Savakis, C., and Bourtzis, K. (2004). Wolbachia-induced cytoplasmic incompatibility as a means for insect pest population control. Proc. Natl. Acad. Sci. 101, 15042–15045. doi: 10.1073/pnas.0403853101
Zhang, H. J., Tang, B. Z., Lin, Y. P., Chen, Z. M., Zhang, X. F., Ji, T. L., et al. (2017). Identification of three prophenoloxidase-activating factors (PPAFs) from an invasive beetle Octodonta nipae maulik (Coleoptera: Chrysomelidae) and their roles in the prophenoloxidase activation. Arch. Insect Biochem. Physiol. 96:e21425. doi: 10.1002/arch.21425
Zhang, Y. K., Ding, X. L., Rong, X., and Hong, X. Y. (2015). How do hosts react to endosymbionts? A new insight into the molecular mechanisms underlying the Wolbachia-host association. Insect Mol. Biol. 24, 1–12. doi: 10.1111/imb.12128
Zhou, F., Lou, Q., Wang, B., Xu, L., Cheng, C., Lu, M., et al. (2016). Altered carbohydrates allocation by associated bacteria-fungi interactions in a bark beetle-microbe symbiosis. Sci. Rep. 6:20135. doi: 10.1038/srep20135
Zhou, W., Rousset, F., and O’Neill, S. (1998). Phylogeny and PCR-based classification of Wolbachia strains using wsp gene sequences. Proc. R. Soc. B Biol. Sci. 265, 509–515. doi: 10.1098/rspb.1998.0324
Keywords: Octodonta nipae, microbiota, 16S rRNA gene pyrosequencing, Wolbachia, life stages, reproductive organs
Citation: Ali H, Muhammad A, Sanda NB, Huang Y and Hou Y (2019) Pyrosequencing Uncovers a Shift in Bacterial Communities Across Life Stages of Octodonta nipae (Coleoptera: Chrysomelidae). Front. Microbiol. 10:466. doi: 10.3389/fmicb.2019.00466
Received: 28 September 2018; Accepted: 21 February 2019;
Published: 12 March 2019.
Edited by:
Suhelen Egan, University of New South Wales, AustraliaReviewed by:
Morten Schiøtt, University of Copenhagen, DenmarkDavide Sassera, University of Pavia, Italy
Copyright © 2019 Ali, Muhammad, Sanda, Huang and Hou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Youming Hou, ymhou@fafu.edu.cn