- Research article
- Open Access
A comparative study of microbial community and dynamics of Asaia in the brown planthopper from susceptible and resistant rice varieties
BMC Microbiologyvolume 19, Article number: 139 (2019)
The brown planthopper (BPH) is likely the most destructive, piercing and sucking monophagous insect pest of rice that causes substantial economic losses to farmers. Although yeast-like symbionts (YLS) and virus transmission have been observed in the BPH, the bacterial population inhabiting the BPH has received minimal research attention. Labelling BPH-associated bacterial species may shed light on BPH biology and the interaction between the BPH and rice to provide novel approaches for the efficient control of this insect pest.
We examined RNA-seq results to identify bacterial populations present in different generations of BPHs maintained on susceptible or resistant rice varieties. Overall, 87 operational taxonomic units (OTUs) were determined from the BPH-F0, F6 and F16 generations. These OTUs had Shannon and Simpson index values of 0.37–0.6 and 0.56–1.19, respectively. The evenness values of 0.7–1.00 showed the vastness of the bacterial diversity recovered from the BPH samples. The results showed high species diversity in the BPHs collected from susceptible rice and a high number of members of unclassified bacteria in the BPHs isolated from resistant rice. We noticed that Proteobacteria OTUs were predominant across all samples. Furthermore, PCR data of Asaia species showed variable DNA amplification across the BPH samples collected from susceptible or resistant varieties. The identification of Asaia in BPH eggs and BPH-egg-infected rice revealed its influence on the interaction between the BPH egg and rice.
The BPHs had clear differences in their microbiomes and in their ability to feed on different rice hosts. These variations could have an essential impact on host adaptation and interaction. These results provide a better understanding of the bacterial diversity and interaction of the microbiome of different generations of BPHs. Furthermore, PCR data of Asaia sp. variation across the BPH samples (isolated from different host genotypes selected from the field and laboratory, including BPH eggs and egg-infected rice tissues), suggest that Asaia could be an important member of the insect microbiome involved in adaptation, its interaction with rice and, most importantly, as a paratransgenic tool for insect control.
Microbes and insects have co-evolved for several hundred million years, and these associations are often reflected by extensive symbiotic relationships between partners. Microbes have been found to be involved in various phenomena of their host insects [1,2,3,4]. Previous reports have suggested that two planthopper species, Laodelphax striatellus and Sogatella furcifera, harbour microbes that can be transferred horizontally between different species of insects and strongly influence their hosts’ immune response, sexual reproduction and cytoplasmic incompatibility [4,5,6]. The role of microbes in insects using various molecular techniques has unveiled the unexplored world of microbe-insect interactions. However, further investigations are required to derive fruitful conclusions.
The brown planthopper, BPH, Nilaparvata lugens (Stål) (Hemiptera: Delphacidae) is a major monophagous agricultural pest of rice crops in the temperate and tropical regions of Asia . The BPH has two winged forms (macropterous and brachypterous) at the adult stage. The macropterous (long-winged) BPH has long-distance migration behaviour  and infests rice fields, while the brachypterous (short-winged) BPH is known for laying many eggs and producing numerous offspring in rice hosts [9, 10] and can also quickly overcome resistance genes of its rice host by developing new virulence factors [11,12,13]. These characteristics make the BPH difficult to forecast and control. The BPH infests rice at all stages of plant growth. As a result of feeding by both nymphs and adults at the base of tillers, plants become yellow and ultimately die. This symptom of crop damage is known as “hopper burn” “(http://www.rkmp.co.in/content/brown-plant-hopper-bph)”.
Furthermore, recent reports revealed that endosymbionts of insects , which show an important role in uric acid metabolism , might also play a major role in insect-plant interactions [16, 17], thereby forming an integral part of BPH adaption to its rice host. However, while BPH-associated microbiomes and their modes of action are known, their stability is not clear. Here, we examine and classify the bacterial population of the BPH in an attempt to explain this interaction and gain important insights into this interaction to develop a new approach to manage this serious insect pest.
Taxonomic diversity index
The taxonomic diversity (normal distribution and alpha- and beta-diversity) of the bacterial population was analysed using RNA-seq data from different generations of adult female BPH samples. Normal distributions of the bacterial population within the BPH samples were calculated using the Shapiro-Wilk W and Anderson-Darling A test (see Additional file 1). The Shapiro-Wilk W index of normality ranged from 0.6 to 0.820, with a P-value of 0.0002 to 0.064. The observed Anderson-Darling A probability distribution ranged from 0.565 to 1.423, with a P-value of 0.0003 to 0.088 (see Additional file 1). Alpha-diversity indices (Shannon, Simpson, Evenness, Brillouin, Fisher alpha, Berger-Parker and Chao1) of the bacterial population were notably distinct within the BPH samples (Table 1). Shannon diversity values varied from 0.37 to 0.64, Simpson diversity values varied from 0.56 to 1.19, Evenness diversity values varied from 0.70 to 1.00 and Brillouin information diversity values varied from 0.34 to 0.99 (Table 1). The Fisher alpha and Berger-Parker species richness values ranged from 0.00 to 2.38 and 0.40 to 0.75, respectively (Table 1). The species richness (Chao1) estimator indicated that species richness ranged from 2 to 4, with the highest value in BPH-F16-TN1 and BPH-F16-IR36 among the BPH samples (Table 1). Notably, the BPH-F16-TN1 sample had more operational taxonomic units (OTUs) than the other BPH samples (Table 1). PAST (Palaeontological Statistics, version 3) software was used to analyse the beta diversity of the bacterial population between the BPH samples (as shown by PCoA). PCoA (principal coordinate analysis; similarity and distance) was calculated using eigenvalues and eigenvectors (coordinates). The similarity/distance eigenvalues were raised to the power of C (the transformation exponent) before Eigen analysis, and the standard value was C = 2. The “Eigenvalue scaling” option scale was used for each axis using the square root of the eigenvalue, and the minimal spanning tree option was based on the selected similarity or distance index. Similarity- and distance-PCoA estimated the presence or absence of bacterial species between the BPH (BPH-F0, BPH-F6-TN1, BHP-F16-TN1, BPH-F6-IR36, BPH-F16-IR36, BPH-F6-RH, and BPH-F16-RH) samples (Fig. 1). Similarity- (Fig. 1a) and distance-PCoA (Fig. 1b) plots revealed the bacterial diversity in the samples and noticeably separated the BPH-F0, BPH-F6-TN1, BHP-F16-TN1, BPH-F6-IR36, BPH-F16-IR36, BPH-F6-RH, and BPH-F16-RH samples. Evaluation of the PCoA data of pair-wise BPH sample comparisons was performed using Bonferroni correction (Table 2). Eight pair-wise BPH sample comparisons (BPH-F0 and BPH-F6-TN1, BPH-F0 and BPH-F6-IR36, BPH-F6-TN1 and BPH-F6-IR36, BPH-F6-TN1 and BPH-F6-RH, BPH-F16-TN1 and BPH-F6-IR36, BPH-F6-IR36 and BPH-F16-IR36, BPH-F6-IR36 and BPH-F6-RH, and BPH-F6-IR36 and BPH-F16-RH) had a Bonferroni corrected P-value of 1.0, and the remaining 34 pair-wise sample comparisons had Bonferroni corrected P-values in the range of 0.987 to 0.00 (Table 2). Furthermore, the Whittaker indices for similarities between bacterial communities of the BPH samples ranged from 0.143–0.666 (Table 3).
Distribution of operational taxonomic units (OTUs) and relative abundance analysis
The identified diverse bacterial species in the BPH samples were classified at different taxonomic levels, phyla, classes, orders, families, and genera (see Additional file 2). The diversity of OTUs at the phylum level across all BPH samples was determined as follows: Proteobacteria, Actinobacteria, Bacteroidetes, Firmicutes, Fusobacteria, Euryarchaeota, p_Unclassified, p_ Unclassified and p_Unclassified (see Additional file 2). Gene Cluster 3.0 and TreeView software was used to generate a heat map representation of the OTUs, and the results evidently revealed the bacterial diversity and the relative abundance of OTUs in the bacterial composition (phylum level) among the BPH samples (BPH-F0, BPH-F6-TN1, BPH-F16-TN1, BPH-F6-IR36, BPH-F16-IR36, BPH-F6-RH, and BPH-F16-RH) (Fig. 1c).
With a total of 51 bacterial OTUs, Proteobacteria was observed as the predominant phylum among the BPH samples. Proteobacteria alone shared 7, 4, 15, 1, 11, 4 and 9 of the total OTUs in the samples BPH-F0, BPH-F6-TN1, BPH-F16-TN1, BPH-F6-IR36, BPH-F16-IR36, BPH-F6-RH and BPH-F16-RH, respectively (Fig. 2a). However, bacterial OTUs identified as Actinobacteria were the second dominant phylum, and Actinobacteria shared 22 OTUs of the total OTUs classified. The remaining 8, 1, 1, 1 and 3 of the total OTUs represented the phylum Firmicutes, Bacteroidetes, Fusobacteria, Euryarchaeota and p_Unclassified, respectively (Fig. 2a, see Additional file 3). The predominant bacterial OTUs, belonging to Proteobacteria among the BPH samples (Fig. 2b), represented the genera Acinetobacter (10 OTUs); Corynebacterium (8 OTUs); Asaia and Cutibacterium (6 OTUs each); Staphylococcus and Brevundimonas (5 OTUs each); Moraxella (4 OTUs); Micrococcus, Sphingomonas and Escherichia (3 OTUs each); and Ralstonia, Massilia and Occidentia (2 OTUs each); Fusobacterium, Rhizobium, Erwinia, Serratia, Leucobacter, Actinomyces, Fluviicola, Rhodococcus, Alkalibacterium, Anoxybacillus, Delftia, Hydrogenophilus, Bartonella, Chryseobacterium, Cuniculiplasma, Propionibacterium, Lactobacillus, Rhizobacter, Neokamagataea, Thiobacillus, Rhodovastum, Vitreoscilla, and Uncultured Enhydrobacter (1 OTUs each).
Comparative analysis and identification of bacterial populations between BPH samples collected from susceptible and resistant rice varieties
To study the comparative abundance of the bacterial species across the BPH samples, we organized the different BPH (F6 and F16 generations) samples, isolated from susceptible (TN1) and resistant (IR36 and RH) rice varieties, into three groups: BPH-(F6 + F16)-TN1, BPH-(F6 + F16)-IR36 and BPH-(F6 + F16)-RH samples (Fig. 2c, see Additional file 4). We observed that the number of bacterial species was higher in the BPH-(F6 + F16)-TN1 sample and that the number of bacterial species was lower in the BPH-(F6 + F16)-RH sample (Fig. 2c). A total of 36, 20 and 18 bacterial species were detected in the BPH-(F6 + F16)-TN1, BPH-(F6 + F16)-IR36 and BPH-(F6 + F16)-RH samples, respectively (see Additional file 4). Comparing the bacterial populations of the BPH-(F6 + F16)-TN1 sample with the BPH-(F6 + F16)-IR36 sample, Micrococcus, Staphylococcus, Moraxella, Chryseobacterium, Cutibacterium, Sphingomonas, Asaia, Thiobacillus, Rhodovastum, Fusobacterium, Occidentia, Bartonella, Rhodococcus, Alkalibacterium, and Anoxybacillus were lower in the BPH-(F6 + F16)-IR36 sample (Fig. 2b, see Additional file 4), while upon comparison of the bacterial populations in the BPH-(F6 + F16)-TN1 sample with the BPH-(F6 + F16)-RH sample, the bacterial population of Micrococcus, Staphylococcus, Corynebacterium, Moraxella, Chryseobacterium, Ralstonia, Massilia, Cutibacterium, Brevundimonas, Sphingomonas, Asaia, Thiobacillus, Rhodovastum, Fusobacterium, Occidentia, Bartonella, Rhodococcus, Alkalibacterium, and Anoxybacillus were lower in the BPH-(F6 + F16)-RH sample (Fig. 2b, see Additional file 4). In contrast, Cuniculiplasma, Propionibacterium, Vitreoscilla, uncultured Enhydrobacter, and Lactobacillus were found to be higher in number in the BPH-(F6 + F16)-IR36 sample, while Acinetobacter, Rhizobacter, and Neokamagataea were found to be higher in the BPH-(F6 + F16)-RH sample (Fig. 2b, see Additional file 4). These decreases or increases in the number of bacterial species may affect the survival rates of BPHs on resistant (IR36 and RH) rice varieties compared with the susceptible (TN1) rice varieties.
Specific and shared bacterial population
The number of specific and shared OTU sequences present in the BPH samples were analysed using a Venn diagram (TBtools (v0.6653) and http://bioinformatics.psb.ugent.be/cgi-bin/liste/Venn/calculate_venn.htpl) (Fig. 3, see Additional file 5). Thirty-nine specific and 14 shared OTU sequences were identified across all samples. Acinetobacter and Asaia species were shared between five BPH samples, while Staphylococcus, Ralstonia, Massilia, Sphingomonas, Occidentia and uncultured bacteria were shared between only two BPH samples. The remaining 25 bacterial species were specific for each BPH sample (Fig. 3, see Additional file 5).
Dynamic variability in the occurrence of Asaia in different BPH samples collected from susceptible and resistant rice varieties
The RNA-seq results revealed that Asaia was the most common bacterial genus between the BPH samples isolated from different rice varieties. Asaia have all of the necessary ecological characteristics (discussed in the discussion section), which indicate its potential to be a member of the paratransgenic tool against the BPH. To study the presence of Asaia in the different BPH samples, the Asaia-specific 16S rRNA gene was amplified by PCR using an Asaia-specific primer (Fig. 4a-h). Amplification of the Asaia-specific 16S rRNA gene (400 bp) was detected in the long-winged adult (male and female) BPH (LWBM-HHZ and LWBF-HHZ) samples collected from the “HHZ” susceptible rice variety. However, we noticed a weak PCR signal for Asaia in the field-collected long-winged adult BPH (Field-LWBM and Field-LWBF) samples (Fig. 4a, b). Furthermore, amplification of the Asaia-specific 16S rRNA gene (400 bp) was observed in the long-winged male BPH (LWBM-IR36) sample, but a weak PCR signal for Asaia was observed in the long-winged female BPH (LWBF-IR36) sample collected from the resistant rice IR36 (Fig. 4c, d). Conversely, amplification of the Asaia-specific 16S rRNA gene (400 bp) was detected in the short-winged adult (male and female) BPH (SWBM-HHZ and SWBF-HHZ) samples from the “HHZ” susceptible rice (Fig. 4e-h), while we noticed a weak PCR signal for Asaia in the field-collected short-winged adult (male and female) BPH (Field-SWBM and Field-SWBF) samples (Fig. 4e, f). Furthermore, amplification of the Asaia-specific 16S rRNA gene (400 bp) was observed in the BPH (SWBM-IR36 and SWBF-IR36) samples (Fig. 4g, h). Each amplified PCR product, Asaia-specific 16S rRNA gene (400 bp), of the BPH samples was commercially sequenced using the Asaia-specific primers (discussed in the Methods section) at www.igebiotech.com, Guangzhou, China, and the obtained sequencing data were reconfirmed as Asaia sp. using the NCBI-BLAST 16S rRNA database (data not shown). Nevertheless, there is a need to study the significant role of Asaia within their insect host.
Identification of Asaia in the BPH eggs and BPH-egg-infected rice
To recognize the possible role of Asaia in the interaction between the BPH and its rice host, we carried out PCR amplification of the Asaia-specific 16S rRNA gene (400 bp) from BPH-egg, BPH-egg-infected and uninfected rice tissue (HHZ susceptible rice) samples. Amplification of the Asaia-specific gene (400 bp) was detected in the BPH eggs and BPH-egg infected rice samples, but there was no amplification of the Asaia-specific gene in the uninfected rice samples (Fig. 4i, j). These PCR products, Asaia-specific 16S rRNA gene (400 bp), were sequenced using the Asaia-specific primers (discussed in the methods section) at www.igebiotech.com, Guangzhou, China, and Sanger sequencing generated sequence data that reconfirmed the sequences as belonging to Asaia sp. using the NCBI-BLAST 16S rRNA database (data not shown). The detection of the PCR signal of the Asaia-specific gene in BPH-egg-infected rice but not in uninfected rice revealed that this bacterium might have an important role in the interaction between the BPH eggs and the host plant upon which they are laid. However, this requires further verification.
The insect-associated bacterial community is one of the key determinants of insect physiology, but there is limited research on the functionality of the BPH bacterial community during the interaction between BPH and rice. The small size of BPH eggs and poor artificial food make research on the BPH difficult . Therefore, this investigation was carried out to identify the bacterial population and to determine its influential role in the interaction between the BPH and rice [19,20,21,22]. Further, it would be worth studying whether BPHs feeding on susceptible and resistant rice genotypes can affect survival due to an increase or decrease in the number of favourable bacteria within them.
This project described various bacterial communities in seven different samples of BPHs. Proteobacteria OTUs were the predominant constituents in all BPH samples, possibly because the feeding stage of the BPH on rice possessed Proteobacteria, which could then be relocated from the eggs to adults . Furthermore, the BPH samples, depicting only the feeding stage, revealed distinct alpha-diversity (Table 2). The previous investigation characterized Proteobacteria as the major bacterial community structure in the habitat  due to their effective colonization [25, 26]. The Proteobacteria community has also been revealed as the major constituent in the Asian rice gall midge, a monophagous dipteran insect , Drosophila, a polyphagous  dipteran insect and blood feeder bugs i.e., Cimex lectularius , Aedes sp.  and ticks .
The bacterial microbiomes in the BPHs contained both unique and abundant OTUs. These unique species, which were unique to a respective sample, may be indicative of their important role in BPH physiology and may influence their cooperation with the rice plant. The survival rates of nymphs (BPH) have been reported to be significantly lower on resistant rice hosts than on susceptible rice hosts [31, 32]. Compared to the bacterial community noted in the BPH-(F6 + F16)-TN1 sample, the BPH-(F6 + F16)-IR36 and BPH-(F6 + F16)-RH samples from the resistant rice (IR36 and RH) hosts revealed a lower bacterial population from the genera Micrococcus, Staphylococcus, Corynebacterium, Moraxella, Chryseobacterium, Ralstonia, Massilia, Cutibacterium, Brevundimonas, Sphingomonas, Asaia, Thiobacillus, Rhodovastum, Fusobacterium, Occidentia, Bartonella, Rhodococcus, Alkalibacterium, and Anoxybacillus (see Additional file 4), although there was an enhanced bacterial population from the genera Cuniculiplasma, Propionibacterium, Vitreoscilla, uncultured Enhydrobacter, Lactobacillus, Acinetobacter, Rhizobacter, and Neokamagataea (see Additional file 4). This decrease in various bacterial species in the BPH-(F6 + F16)-IR36 and BPH-(F6 + F16)-RH samples may be the cause for BPH-(F6 + F16)-IR36 and BPH-(F6 + F16)-RH mortality because host invariability and demise are equated to change in the native bacterial population, as hypothesized by researchers .
Acinetobacter, Corynebacterium and Asaia are opportunistic bacteria that cause infections in humans [34,35,36]. These 3 bacterial species were identified in all BPH samples. The ten observed Acinetobacter in the samples had maximum sequence identity to A. vivianii, A. radioresistens, A. halotolerans, A. soli, A johnsonii and A. calcoaceticus (see Additional file 2). Furthermore, it has been reported that plasmids are key elements in prokaryotic evolution and their adaptation to variable environmental circumstances . The existence and broad diversity of plasmids in Acinetobacter are well known for natural transformation , metabolic activities or metal resistance [39,40,41] in their host. However, in the context of BPH, there is still no report regarding whether a similar function of Acinetobacter plasmids within the BPH shows any impact on insight into the BPH-microbe interaction. This is yet to be proven.
Insect cells cannot secrete all kinds of enzymes and thus require microbes to supply specific enzymes for nutrient metabolism . A recent study observed that enzymes are secreted by microbes into their insect host to balance the insect life cycle . In our study, Corynebacterium and Staphylococcus were present among the BPH samples, and we speculate that these microbes may be involved in manipulating the biological functions of their insect host [42, 44, 45]. Nevertheless, this is yet to be elucidated.
The presence of Asaia sp. in agricultural pests such as Nilaparvata lugens , Scaphoideus titanus , and Pieris rapae  has previously been demonstrated. A previous study reported that Asaia was observed only in front of the ovipositor sheath in the biotype 2 (Mudgo, carrying the resistant gene (Bph1) female BPH . In our study, all RNA-sequencing data showed that Asaia was present in all BPH (excluding BHP-F0 and BPH-F16-IR36) samples isolated from susceptible (TN1) and resistant (IR36 and RH) rice varieties, and further PCR signals reconfirmed the presence of Asaia sp. in all BPH samples collected from field, susceptible (Huang Hua Zhan) and resistant (IR36), rice varieties. In the case of Gluconacetobacter , Asaia may contribute to metabolic functions that influence the adaptation of the BPH to both susceptible and resistant rice varieties. However, more information is required to confirm this hypothesis.
Furthermore, Asaia have also been reported in multiple symbioses with the mosquito. They influence the interaction with both Plasmodium and Anopheles , possess strong adhesive properties , act as an immune-modulator activating antimicrobial peptide expression , accelerate larval development [51, 52], are transmitted from mother to offspring likely by a mechanism of egg smearing , are capable of cross-colonizing insects of phylogenetically distant genera and order , are proposed as a tool for the control of mosquito-borne diseases, specifically malaria [55, 56], and have positive markers for functional nitrogenise activity . Therefore, a similar role is likely to be possible in the BPH as well. In our study, Asaia was observed in BPH-egg and BPH-egg-infected rice stem tissues but was absent in uninfected rice stem tissues. Asaia may transfer from female BPHs to rice, likely mediated by an egg-smearing mechanism while BPHs are laying eggs inside the rice stem. The strong adhesive and antibacterial properties of Asaia may influence the attachment and survival of the eggs on rice. Taking this observation further, we speculate that Asaia may also play a role in the fitness of the BPH and that its cross-colonization property may support its use as a paratransgenic tool for the control of BPH in the infestation of rice. However, this requires further study.
This study demonstrated that members of Proteobacteria were predominant in all BPH samples. The BPH-associated bacterial diversity may be due to the recruitment of various microbiomes from the extrinsic habitat. This diverse bacterial population can play an important role in BPH-rice interactions and influence BPH survival in rice (susceptible and resistant) hosts. Furthermore, the thorough Asaia documentation thus confirmed and established its interaction with adult BPHs, BPH eggs and rice. Now, with the extensive bacterial population catalogue thus established, this study will facilitate future experiments including the interaction of bacterial species with the insect phenotype. Even more broadly, it is necessary to further investigate the bacterial effect on insect associations with the host and to use native bacteria as a paratransgenic tool for insect control.
Total RNA-sequence analysis for bacterial populations in different BPH samples
Laboratory strains of BPH were collected in September 2007 from Guangdong Academy of Agricultural Sciences (GAAS), Guangzhou, in Guangdong Province, China, and a colony has been reared continuously on susceptible and resistant rice varieties in our laboratory since then. Different generations (F0, F6, and F16) of adult female (each 15) BPH samples from susceptible (TN1) and resistant rice (IR36 carrying the bph2 resistance gene and RH carrying the Bph3 resistance gene)  varieties were used for the current study. Total RNA was isolated from adult BPH females using a total RNA extraction kit (Omega, Norcross, GA, USA) according to the manufacturer’s instructions. Library construction and sequencing were performed by Vazyme Biotech Co., Ltd. (Nanjing, China) using the Illumina-HiSeq 2500 platform and PE100 as a sequencing strategy (data not shown in this manuscript). A total of 54.61 GB high quality reads (data not shown in this manuscript) were obtained from the BPH samples. An average of 7.8 G reads from seven different BPH samples named BPH-F0 (F0 generation), BPH-F6-TN1 (F6 generation BPH from TN1), BHP-F16-TN1 (F16 generation BPH from TN1), BPH-F6-IR36 (F6 generation BPH from IR36), BPH-F16-IR36 (F16 generation BPH from IR36), BPH-F6-RH (F6 generation BPH from RH) and BPH-F16-RH (F16 generation BPH from RH) were analysed to identify the microbial community structure across the samples using the 16S ribosomal RNA Sequences database from National Center for Biotechnology Information (NCBI), USA (https://www.ncbi.nlm.nih.gov) (see Additional files 6, 7, 8, 9, 10, 11 and 12).
To study the dynamics of Asaia in the short- and long-winged (adults; male and female) BPHs, the BPH samples were randomly collected from the entire rice field in Shaoguan city in Guangdong Province, China and stored in 100% alcohol. The rest of the BPH samples used in this study were maintained in the greenhouse at State Key Laboratory of Biocontrol, Sun Yat-sen University, Guangzhou in Guangdong Province, China under standard conditions . The procedure described by Ojha et al.  was applied to grow different rice varieties. The BPHs (50 each of females and males) were released on 10–15-day-old plants, and on day 1, the plastic trays with plants were transferred to a humidity (70–90%) chamber for the process of egg incubation. The adult BPHs were collected from “Huang Huas Zhan” (HHZ, a Chinese susceptible rice variety) and IR36. Entomological needles were used to isolate intact eggs from “HHZ” rice within 48 h and preserved in alcohol (100%). Nearly 200 intact eggs were collected and used for the entire experiment. Thirty short- and long-winged BPH adults were also collected from field rice, both susceptible and resistant rice varieties. To study the Asaia identification on the rice host, 300 mg of BPH egg-infected and uninfected tissues of susceptible rice were collected.
Isolation of DNA from the BPH and rice samples
The collected BPH samples were frozen at − 20 °C for 5 min, immersed in 5% NaOCl (sodium hypochlorite), and rinsed with sterile distilled water five times to remove surface microorganisms. Total DNA was isolated from BPH adults and eggs using the Insect DNA kit (Omega Bio-tek, Inc., Norcross, GA, USA) according to the manufacturer’s instructions. Total DNA from infected and uninfected rice tissues was extracted using the Plant DNA Kit (Omega Bio-tek, Inc., Norcross, GA, USA) according to the manufacturer’s instructions. DNA concentrations were estimated using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, USA).
Primer design and PCR amplification of Asaia from different BPH and rice samples
The Asaia-specific primers were designed using the Asaia-specific 16S rRNA gene sequences belonging to the accession numbers:- NR_122089.1, NR_112880.1, NR_024810.1, NR_114144.1, NR_041564.1, NR_113849.1, NR_024728.1, NR_112879.1, NR_112953.1, NR_113845.1, and NR_024738.1 (see Additional file 13). PCR amplification of the ~ 400 bp 16S rRNA gene sequence of Asaia was performed using the Asaia primers (forward: 5′-GGCGCGTAGGCGGTTTACAC-3′ and reverse: 5′-TGCGCGTTGCTTCGAATTAAACCA-3′) to evaluate their presence in various samples of BPHs and rice. Actin-gene (NCBI-GenBank: EU179846.1) specific-primers (forward: 5′-TGCGTGACATCAAGGAGA-3′ and reverse: 5′-TGTTCCAGCCTTCCTT-3′) were used and PCR-amplified ~ 174 bp actin gene, a housekeeping gene used as an internal control for the normalization of PCR data. Ten ng DNA of short- and long-winged (adults), egg, egg-infected and uninfected rice tissues were used as the template for PCR, and the PCR volume was 25 μl. PCRs were performed with 200 μM dNTPs (New England Biolabs Inc.), 0.5 U Q5 Hot start high-fidelity Taq DNA polymerase (New England Biolabs Inc.), and 13 μM primers (each). The thermal cycler protocol was 98 °C for 5 min, 35 cycles at 98 °C for 30 s, 56 °C for 45 s, and 72 °C for 30 s, and a final 5 min extension at 72 °C. The PCR conditions for actin-gene amplification from each BPH sample were 94 °C for 5 min, 35 cycles at 94 °C for 30 s, 54 °C for 30 s, and 72 °C for 45 s, and a final 5 min extension at 72 °C. Amplicons were run on 1% (w/v) agarose gel to analyse band size. The EtBr-stained gels were viewed and photographed in our gel documentation system (GL212PRO Imager, Carestream Health Inc., USA), and captured-gel images were saved in JPEG format. Then, ImageJ software (https://imagej.nih.gov/ij/) was used to estimate the relative colour intensities of each amplicons.
The normal distribution and alpha diversity (Shannon and Simpson’s diversity indices, Evenness, Fisher alpha, Brillouin, Fisher alpha, Berger-Parker and Chao1) index were used to evaluate the differences in bacterial populations between BPH samples (obtained from both susceptible and resistant rice varieties) used in the current study. PCoA plots were plotted using PAST (Palaeontological Statistics, version 3) software. To show the relative abundance of microbial phyla for the respective BPH samples, heat maps were created using Gene Cluster 3.0 and TreeView software.
operational taxonomic units
polymerase chain reaction
Douglas AE. Nutritional interactions in insect-microbial symbiosis: aphids and their symbiotic bacteria Buchnera. Annu Rev Entomol. 1998;43:17–37.
Baumann P. Biology of bacteriocyte-associated endosymbionts of plant sap-sucking insects. Annu Rev Entomol. 2005;59:155–89.
Hurst GD, Jiggins FM. Male-killing bacteria in insects: mechanisms, incidence, and implications. Emerg Infect Dis. 2000;6(4):329–36.
Stouthamer R, Breeuwer JA, Hurst GD. Wolbachia pipientis: microbial manipulator of arthropod reproduction. Annu Rev Entomol. 1999;53:71–102.
Heath BD, Butcher RD, Whitfield WG, Hubbard SF. Horizontal transfer of Wolbachia between phylogenetically distant insect species by a naturally occurring mechanism. Curr Biol. 1999;9:313–6.
Siozios S, Sapountzis P, Ioannidis P, Bourtzis K. Wolbachia symbiosis and insect immune response. Insect Science. 2008;15:89–100.
Brown planthopper: threat to rice production in Asia. International Rice Research Institute, Los Banos, Laguna, Philippines, The brown planthopper 1979;p.3–17.
Cheng JA, Holt J. A systems analysis approach to brown planthopper control on rice in Zhejinag Province, China 1 simulation of outbreaks. J Appl Ecol. 1990;271:85–99.
Sun ZX, Zhai YF, Zhang JQ, Kang K, Cai JH, Fu Y, et al. The genetic basis of population fecundity prediction across multiple field populations of Nilaparvata lugens. Mol Ecol. 2015;24(4):771–84.
Noda H, Kawai S, Koizumi Y, Matsui K, Zhang Q, Furukawa S, et al. Annotated ESTs from various tissues of the brown planthopper Nilaparvata lugens: a genomic resource for studying agricultural pests. BMC Genomics. 2008;9:117.
Horgan F. Mechanisms of resistance: a major gap in understanding planthopper-rice interactions. Planthoppers: new threats to the sustainability of intensive rice production systems in Asia. Edited by Heong KL, Hardy B. Los Baños (Philippines): International Rice Research Institute 2009;281–302.
Sōgawa K. The rice brown planthopper: feeding physiology and host plant interactions. Annu Rev Entomol. 1982;27:49–73.
Cook AG, Denno RF. Planthopper/plant interactions: feeding behavior, plant nutrition, plant defense, and host plant specialization. In: Denno RF, Pefect TJ, editors. Planthoppers. Boston, MA: Springer; 1994. p. 114–39.
Chen CC, Cheng LL, Hou RF. Studies on the intracellular yeast like symbiote in the brown planthopper, Nilaparvata lugens Stal. II. Effects of antibiotics and elevated temperature on the symbiotes and their host. J Appl Entomol. 1981;92:440–9.
Sasaki T, Kawamura M, Ishikawa H. Nitrogen recycling in the brown planthopper, Nilaparvata lugens: involvement of yeast-like endosymbionts in uric acid metabolism. J Insect Physiol. 1996;42:125–9.
Lu ZX, Yu XP, Chen JM, Zheng XS, Xu HX, Zhang JF, Chen LZ. Dynamics of yeats-like symbiote and its relationship with the virulence of brown planthopper, Nilaparvata lugens Stal, to resistant rice varieties. J Asia Pac Entomol. 2004;7(3):317–23.
Chen YH, Bernal CC, Tan J, Horgan FJ, Fitzgerald MA. Planthopper ‘adaptation’ to resistant rice varieties: changes in amino acid composition over time. J Insect Physiol. 2011;57:1375–84.
Qiang F, Zhang Z, Cui H, Lai F, Sun Z. A chemically defined diet enables continuous rearing of the brown planthopper, Nilaparvata lugens (Stal) (Homoptera: Delphacidae). Appl Entomol Zool. 2001;36(1):111–6.
Ojha A, Sinha DK, Padmakumari AP, Bentur JS, Nair S. Bacterial community structure in the Asian rice gall midge reveals a varied microbiome rich in Proteobacteria. Sci Rep. 2017;7(1):9424.
Karakus M, Karabey B, Orçun Kalkan Ş, Özdemir G, Oğuz G, Erişöz Kasap Ö, et al. Midgut bacterial diversity of wild populations of Phlebotomus (P.) papatasi, the vector of zoonotic cutaneous Leishmaniasis (ZCL) in Turkey. Sci Rep. 2017;7(1):14812.
Krams IA, Kecko S, Jõers P, Trakimas G, Elferts D, Krams R, et al. Microbiome symbionts and diet diversity incur costs on the immune system of insect larvae. J Exp Biol. 2017;220:4204–12.
Yun JH, Roh SW, Whon TW, Jung MJ, Kim MS, Park DS, et al. Insect gut bacterial diversity determined by environmental habitat, diet, developmental stage, and phylogeny of host. Appl Environ Microbiol. 2014;80:5254–64.
Douglas AE. Multiorganismal insects: diversity and function of resident microorganisms. Annu Rev Entomol. 2015;60:17–34.
Galand PE, Casamayor EO, Kirchman DL, Lovejoy C. Ecology of the rare microbial biosphere of the Arctic Ocean. Proc Natl Acad Sci U S A. 2009;106:22427–32.
Rampelotto PH, de Siqueira Ferreira A, Barboza AD, Roesch LF. Changes in diversity, abundance, and structure of soil bacterial communities in Brazilian savanna under different land use systems. Microb Ecol. 2013;66:593–607.
Roesch LF, Fulthorpe RR, Riva A, Casella G, Hadwin AK, Kent AD, et al. Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 2007;1:283–90.
Broderick NA, Lemaitre B. Gut-associated microbes of Drosophilla melanogaster. Gut Microbes. 2012;3:307–21.
Meriweather M, Mathews S, Rio R, Baucom RS. A 454 survey reveals the community composition and core microbiome of the common bed bug (Cimex lectularius) across an urban landscape. PLoS One. 2013;8(4):e61465.
Zouache K, Raharimalala FN, Raquin V, Tran-Van V, Raveloson LH, Ravelonandro P, Mavingui P. Bacterial diversity of field-caught mosquitoes, Aedes albopictus and Aedes aegypti, from different geographic regions of Madagascar. FEMS Microbiol Ecol. 2011;75:377–89.
Narasimhan S, Fikrig E. Tick microbiome: the force within. Trends Parasitol. 2015;31:315–23.
Wei Z, Hu W, Lin Q, Cheng X, Tong M, Zhu L, et al. Understanding rice plant resistance to the brown planthopper (Nilaparvata lugens): a proteomic approach. Proteomics. 2009;9(10):2798–808.
Alagar M, Suresh S, Samiyappan R, Saravanakumar D. Reaction of resistant and susceptible rice genotypes against brown planthopper (Nilaparvata lugens). Phytoparasitica. 2007;35:346.
Blaser MJ, Webb GF. Host demise as a beneficial function of indigenous microbiota in human hosts. mBio. 2014;5(6):e02262–14.
Doughari HJ, Ndakidemi PA, Human IS, Benade S. The ecology, biology and pathogenesis of Acinetobacter spp.: an overview. Microbes Environ. 2011;26(2):101–12.
Valero-Guillen PL, Yague G, Segovia M. Characterization of acyl-phosphatidylinositol from the opportunistic pathogen Corynebacterium amycolatum. Chem Phys Lipids. 2005;133(1):17–26.
Alauzet C, Teyssier C, Jumas-Bilak K, Gouby A, Chiron R, Rabaud C, et al. Gluconobacter as well as Asaia species, newly emerging opportunistic human pathogens among acetic acid bacteria. J Clin Microbiol. 2010;48(11):3935–42.
Eberhard WG. Evolution in bacterial plasmids and levels of selection. Q Rev Biol. 1990;65(1):3–22.
Palmen R, Vosman B, Buijsman P, Breek CK, Hellingwerf KJ. Physiological characterization of natural transformation in Acinetobacter calcoaceticus. J Gen Microbiol. 1993;139(2):295–305.
Winstanley C, Taylor SC, Williams PA. pWW174: a large plasmid from Acinetobacter calcoaceticus encoding benzene catabolism by the beta-ketoadipate pathway. Mol Microbiol. 1987;1(2):219–27.
Schembri MA, Bayly RC, Davies JK. Cloning and analysis of the polyhydroxyalkanoic acid synthase gene from an Acinetobacter sp: evidence that the gene is both plasmid and chromosomally located. FEMS Microbiol Lett. 1994;118:145–52.
Vaneechoutte M, Nemec A, Musilek M, van der Reijden TJ, van der Barselaar M, Tjernberg I, et al. Description of Acinetobacter venetianus ex Di cello et al. 1997 sp. nov. Int J Syst Evol Microbiol. 2009;59:1376–81.
Engel P, Moran NA. The gut microbiota of insects-diversity in structure and function. FEMS Microbiol Rev. 2013;37(5):699–735.
Visôtto LE, Oliveira MG, Guedes RN, Ribon AO, Good-God PI. Contribution of gut bacteria to digestion and development of the velvetbean caterpillar, Anticarsia gemmatalis. J Insect Physiol. 2009;55(3):185–91.
Feng W, Wang XQ, Zhou W, Liu GY, Wan YJ. Isolation and characterization of lipase-producing bacteria in the intestine of the silkworm, Bombyx mori, reared on different forage. J Insect Sci. 2011;11(135).
Khurana S, Sanli G, Powers DB, Anderson S, Blaber M. Molecular modeling of substrate binding in wild-type and mutant Corynebacteria 2,5-diketo-D-gluconate reductases. Proteins. 2000;39(1):68–75.
Tang M, Lv L, Jing S, Zhu L, He G. Bacterial symbionts of the brown planthopper, Nilaparvata lugens (Homoptera: Delphacidae). Appl Environ Microbiol. 2010;76(6):1740–5.
Gonella E, Crotti E, Rizzi A, Mandrioli M, Favia G, Daffonchio D, Alma A. Horizontal transmission of the symbiotic bacterium Asaia sp. in the leafhopper Scaphoideus titanus ball (Hemiptera: Cicadellidae). BMC Microbiol. 2012;12(Suppl 1):S4.
Robinson CJ, Schloss P, Ramos Y, Raffa K, Handelsman J. Robustness of the bacterial community in the cabbage white butterfly larval midgut. Microb Ecol. 2010;59(2):199–211.
Capone A, Ricci I, Damiani C, Mosca M, Rossi P, Scuppa P, et al. Interactions between Asaia, Plasmodium and Anopheles: new insights into mosquito symbiosis and implications in malaria symbiotic control. Parasit Vectors. 2013;6(1):182.
Kruk T, Szczepanowicz K, Kregiel D, Szyk-Warszynska L, Warszynski P. Nanostructured multilayer polyelectrolyte films with silver nanoparticles as antibacterial coatings. Colloids SurfB Biointerfaces. 2015;137:158–66.
Mitraka E, Stathopoulos S, Siden-Kiamos I, Christophides GK, Louis C. Asaia accelerates larval development of Anopheles gambiae. Pathog Glob Health. 2013;107(6):305–11.
Chouaia B, Rossi P, Epis S, Mosca M, Ricci I, Damiani C, et al. Delayed larval development in Anopheles mosquitoes deprived of Asaia bacterial symbionts. BMC Microbiol. 2012;12(Suppl1):S2.
Damiani C, Ricci I, Crotti E, Rossi P, Rizzi A, Scuppa P, et al. Mosquito-bacteria symbiosis: the case of Anopheles gambiae and Asaia. Microb Ecol. 2010;60(3):644–54.
Crotti E, Damiani C, Pajoro M, Gonella E, Rizzi A, Ricci I, et al. Asaia, a versatile acetic acid bacterial symbiont, capable of cross-colonizing insects of phylogenetically distant genera and orders. Environ Microbiol. 2009;11(12):3252–64.
Epis S, Gaibani P, Ulissi U, Chouaia B, Ricci I, Damiani C, et al. Do mosquito-associated bacteria of the genus Asaia circulate in humans? Eur J Clin Microbiol Infect Dis. 2012;31(6):1137–40.
Bongio NJ, Lampe DJ. Inhibition of Plasmodium berghei development in mosquitoes by effector proteins secreted from Asaia sp. bacteria using a novel native secretion signal. PLoS One. 2015;10(12):e0143541.
Samaddar N, Paul A, Chakravorty S, Chakraborty W, Mukherjee J, Chowdhuri D, Gachhui R. Nitrogen fixation in Asaia sp. (family Acetobacteraceae). Curr Microbio. 2011;63:226–31.
Jena KK, Kim SM. Current status of brown planthopper (BPH) resistance and genetics. Rice. 2010;3:161–71.
Pang R, Qiu J, Li T, Yang P, Yue L, Pan Y, Zhang W. The regulation of lipid metabolism by a hypothetical P-loop NTPase and its impact on fecundity of the brown planthopper. Biochim Biophys Acta Gen Subj. 2017;1861(7):1750–8.
We thank Longgyu Yuan for their assistance with insect maintaining in green house at State Key Laboratory of Biocontrol, SunYat-sen University, China. The Authors thank Dr. Ruchi Agrawal, International Centre for Genetic Engineering and Biotechnology (ICGEB), New Delhi, for critical review of the manuscript.
This work was funded by the National Natural Science Foundation of China (U1401212) to WZ. AO thanks State Key Laboratory of Biocontrol visiting scholar foundation for a research grant (SKLBC15F02), SunYat-sen University, China. The funders had no role in study design, data collection and analysis, preparation of the manuscript, or decision to publish.
Ethics approval and consent to participate
Insect samples were either maintained at our greenhouse at State Key Laboratory of Biocontrol, Sun Yat-sen University or obtained from the entire rice field in the Shaoguan city in Guangdong Province, China - no permission required.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Normal distribution test of bacterial population in the BPH samples. (DOCX 14 kb)
Details of OTU ID and its NCBI-blast, 16S rRNA database, and details of the closest-matched accession numbers for BPH-F0, BPH-F6-TN1, BPH-F16-TN1, BPH-F6-IR36, BPH-F16-IR36, BPH-F6-RH, and BPH-F16-RH samples. (PDF 92 kb)
Overview of the bacterial population across the BPH samples. (PDF 85 kb)
Comparison of the bacterial populations across the BPH (F6 + F16 generation BPHs collected from susceptible, TN1, and resistant, IR36 and RH, rice varieties) samples at the genera level. Reduction and increase in the bacterial population across the BPH (F6 + F16 generation BPHs collected from the susceptible, TN1, rice variety) samples when compared with the BPH (F6 + F16 generation BPHs collected from resistant, IR36 and RH, rice variety) samples at the genus level. (DOCX 17 kb)
Identification of a unique bacterial population among the BPH samples. Total number of unique bacteria among the BPH samples (a), shared bacteria between the BPH samples (b), and (c) specific bacterial name among the BPH samples. (DOCX 16 kb)
Bacterial sequences of the F0 generation BPHs. (PDF 90 kb)
Bacterial sequences of the F6 generation BPHs from the susceptible, TN1, rice variety. (PDF 84 kb)
Bacterial sequences of the F16 generation BPHs from the TN1 rice variety. (PDF 86 kb)
Bacterial sequences of the F6 generation BPHs from the resistant, IR36, rice variety. (TXT 22 kb)
Bacterial sequences of the F16 generation BPHs from the IR36 rice variety. (XLSX 18 kb)
Bacterial sequences of the F6 generation BPHs from the resistant, RH, rice variety. (XLSX 10 kb)
Bacterial sequences of the F16 generation BPHs from the RH rice variety. (PDF 87 kb)
Asaia-specific primer design. Asaia-specific 16S rRNA gene sequences belonging to the NCBI Accession Numbers NR_122089.1, NR_112880.1, NR_024810.1, NR_114144.1, NR_041564.1, NR_113849.1, NR_024728.1, NR_112879.1, NR_112953.1, NR_113845.1, and NR_024738.1 were used. (PDF 85 kb)