- Open Access
Gut symbiont enhances insecticide resistance in a significant pest, the oriental fruit fly Bactrocera dorsalis (Hendel)
Microbiomevolume 5, Article number: 13 (2017)
Symbiotic bacteria affect insect physiology and ecology. They may also mediate insecticide resistance within their hosts and thereby impact pest and vector control practices. Here, we document a novel mechanism of insecticide resistance in which a gut symbiont of the tephritid pest fruit fly Bactrocera dorsalis enhances resistance to the organophosphate insecticide trichlorphon.
We demonstrated that the gut symbiont Citrobacter sp. (CF-BD) plays a key role in the degradation of trichlorphon. Based on a comparative genomics analysis with other Citrobacter species, phosphatase hydrolase genes were identified in CF-BD. These CF-BD genes had higher expression when trichlorphon was present. Bactrocera dorsalis inoculated with isolated CF-BD obtained higher trichlorphon resistance, while antibiotic-treated flies were less resistant confirming the key role of CF-BD in insecticide resistance.
Our findings suggest that symbiont-mediated insecticide resistance can readily develop in B. dorsalis and may represent a more widely relevant insecticide resistance mechanism than previously recognized.
Insects can possess symbiotic microorganisms in their gut lumen, in specialized organs, or within cells [1–4]. In general, such microbial partners can contribute to the nutrition of various insect groups , defense against natural enemies , reproductive traits , and other physiological and ecological properties of insects [8–12]. Some symbiotic bacteria also mediate detoxification of insect diets [13–16] and, similarly, of insecticides, therefore conferring insecticide resistance to their hosts as it has originally been discovered for the apple maggot Rhagoletis pomonella  and more recently demonstrated for stinkbugs [18, 19].
Chemical insecticides have been widely used to control insect pests and vectors ; however, many insect pests and vectors have evolved strong resistance to a diverse range of insecticides. The mechanisms underlying insecticide resistance vary across pesticides and include changes of drug target sites, increased expression of degrading enzymes, and enhanced drug excretion [21, 22]. The frequent failure of chemical control has globally drawn major research attention to resistance mechanisms and management. For example, it has been determined that certain bacteria also possess the ability to degrade pesticides [23, 24], suggesting that symbiotic bacteria of insects may also contribute to insecticide resistance. However, besides the examples of R. pomonella  and stinkbugs , it is not known whether bacterially facilitated insecticide resistance also occurs in other insect pest taxa of economic significance and, further, what the general mechanisms of symbiont-facilitated insecticide resistance are.
Previous studies have found that intensive insecticide application can accelerate insecticide biodegradation in the environment [25, 26], including by bacteria that are capable of degrading organophosphorus compounds . Studies have found that the biochemistry of organophosphorus compound degradation is identical in most bacteria. The functional enzyme in this process, organophosphate hydrolase or phosphotriesterase, is an organophosphate hydrolase encoded by the opd (organophosphate-degrading) gene, which has been isolated from taxonomically different bacterial species and from various geographical regions .
The oriental fruit fly Bactrocera dorsalis (Hendel) (Diptera: Tephritidae) is a significant pest species that damages a wide range of fruit and other horticultural products [29, 30], causing major financial losses to horticulture . Trichlorphon [dimethyl (2, 2, 2-trichloro-1-hydroxyethyl) phosphate] is a moderately toxic organophosphate insecticide that has been widely used to control this pest because of its low toxicity to humans and its high efficacy; however, resistance to this pesticide in B. dorsalis has been increasing [32, 33], thus threatening the effective management of the oriental fruit fly.
Several studies have been performed to elucidate the mechanism of B. dorsalis resistance to trichlorphon [33, 34]. The functional proteins involved in the resistance response to trichlorphon were identified based on proteomic analyses of B. dorsalis treated with trichlorphon . Independent from this, some pesticide-degrading bacteria have previously been isolated from trichlorphon-contaminated soil . We therefore hypothesized that symbiotic bacteria of B. dorsalis could increase its resistance to chemical insecticides by degrading trichlorphon. Here, we compared the diversity and abundance of the gut bacteria of B. dorsalis strains with different resistance levels in order to identify any changes in their bacterial community composition. We then isolated and cultivated the gut bacteria that were more prevalent in resistant B. dorsalis lines and obtained one bacterium, Citrobacter freundii BD (CF-BD) that was able to degrade trichlorphon. We then manipulated insecticide resistance by adding or removing CF-BD and provided further evidence for a close relationship between this bacterium and the pesticide resistance in B. dorsalis. Based on the annotated genome sequence of the isolated bacterium, together with a comparative genomics analysis, new phosphatase genes were identified in CF-BD when compared with other Citrobacter species and their expression studied when exposed to trichlorphon.
Trichlorphon-resistant and trichlorphon-susceptible fly strains exhibit different gut symbiotic bacterial communities
To confirm the resistance levels of the B. dorsalis susceptible (SS) and resistant (RS) strains, the toxicity of ingested trichlorphon was examined; the dose-response curves are shown in Fig. 1a. The median lethal concentration (LC50) of trichlorphon was 3.77 mg/L for SS and 123.57 mg/L for RS.
To assess the diversity of the gut symbiotic bacteria of SS and RS, the variable region of the 16S rDNA was sequenced via high-throughput amplicon sequencing. By filtering sequence reads, tags were generated for the six samples analyzed. As shown in Additional file 1: Table S1, the greatest number of tags was observed in RS1, and the fewest tags were found in SS3. OTUs were successfully generated for the six samples; as shown in Additional file 2: Table S2, Additional file 3: Table S3 produced the greatest number of OTUs. OTU abundance profiles were obtained for all of the samples by combining the OTU species annotation information and abundance information across different samples.
Specifically, more than 81.25% of the sequence reads were successfully annotated to the genus level, although less than 1.42% of the reads were annotated to the species level. The structure of bacterial communities in RS differed markedly from those in SS (Fig. 1b, c). A stack map at the genus level for the six samples demonstrates that the number of reads assigned to Citrobacter was significantly greater among RS than SS (independent-sample t test, t = 3.354, df = 4, p = 0.028; Fig. 1d, e).
Citrobacter sp. isolation and identification
In order to isolate Citrobacter sp., subculturing was performed on BHI agar flat plates, and the bacteria were identified by 16S rDNA sequencing and analyses of their physiological and biochemical characteristics. After 24 h, Citrobacter sp. colonies grown from the SS and RS isolates were identified by a white color, raised center, clean margin, and smooth, wet surface (Fig. 2a). 16S rDNA amplification and sequencing yielded a fragment of 1521 bp. Based on a BLAST search against GenBank, the 16S rDNA sequence exhibited 99% identity with Citrobacter freundii, and we named the strain CF-BD (Fig. 2b). The sequence was deposited in GenBank (accession number KR002082). Identification using the GYZ-15 eV system also indicated that CF-BD was most similar with Citrobacter freundii (Table 1). The Gram stain test was negative, and morphological analysis revealed that the bacterium was a short rod of approximately 0.5–0.9 μm × 1.2–2.0 μm with blunt, round ends (Fig. 2c).
In vivo localization of CF-BD symbionts
The gut of B. dorsalis was dissected and subjected to fluorescence in situ hybridization (FISH) targeting 16S rRNA of CF-BD symbionts. CF-BD signals were consistently localized in the midgut crypts of the whole gut of B. dorsalis (Fig. 2d, e).
Detection of CF-BD symbionts from diverse populations
Based on the genome sequence information, a primer pair for recA PCR amplification was designed to screen for CF-BD infections in flies. Flies from wild populations were subjected to diagnostic and quantitative PCR. As a result, CF-BD was detected in all 78 individuals from 13 populations, with 100% infection. However, individuals from different populations differed significantly in their CF-BD densities. Flies from Zhengzhou, Henan, had the highest density of CF-BD, and flies from Zhanjiang, Guangdong, had the lowest density of CF-BD (F = 6.38, df = 12, p < 0.01, Fig. 3).
The insecticide resistance of RS flies is decreased by CF-BD deprivation
Antibiotic sensitivity tests revealed that CF-BD is highly sensitive to streptomycin, sulfisoxazole, sulfamethoxazole, chloramphenicol, ciprofloxacin, tetracycline, and nalidixic acid. However, CF-BD was resistant to amikacin, ampicillin, cefazolin, and amoxicillin (Additional file 4: Table S4). Thus, we used streptomycin to clear CF-BD from the guts of RS flies. After feeding with streptomycin, the number of CF-BD colonies obtained from RS fly gut samples was significantly reduced compared with the RS gut samples from flies that were fed sterile water (independent-sample t test, t = 15.496, df = 8, p < 0.01; Fig. 4a). This indicated that the abundance of CF-BD in the fly gut was significantly decreased after feeding with an antibiotic solution. Moreover, the RS flies that were fed the streptomycin solution were more sensitive to trichlorphon than the RS flies (two-way ANOVA F = 66.905, df = 1, p < 0.0001; LC50 = 89.78 mg/L for RS+ streptomycin, LC50 = 133.31 mg/L for RS; Fig. 4b), indicating that CF-BD was required for enhanced resistance to trichlorphon.
The insecticide resistance of SS flies is increased after CF-BD supplementation
We also attempted to increase the abundance of CF-BD in the guts of SS flies by feeding them with CF-BD inoculum. After feeding SS flies CF-BD, the amount of gut CF-BD was significantly higher than in the control flies that were given water (independent-sample t test, t = 7.67, df = 8, p < 0.01; Fig. 4c). Moreover, the SS flies that were fed the CF-BD inoculum were more resistant to trichlorphon than the control flies (two-way ANOVA F = 7.848, df = 1, p = 0.011; LC50 = 2.35 mg/L for SS+ CF-BD inocula, LC50 = 1.28 mg/L for SS; Fig. 4d), which indicated that increases in CF-BD in the guts of flies increased their resistance to trichlorphon.
Trichlorphon tolerance and degradation ability of CF-BD
To investigate the function of CF-BD in relation to trichlorphon resistance, the trichlorphon tolerance and degradation ability of CF-BD was tested by adding CF-BD to trichlorphon-enriched BHI agar plates. As the growth rates and diameters of the colonies did not differ between the trichlorphon-enriched and non-enriched plates (independent-sample t test, t = 1.247, df = 10, p = 0.241), CF-BD was not inhibited by trichlorphon (Fig. 4e, f). Gas chromatography analysis revealed a greater amount of trichlorphon in the filtrate purified from the trichlorphon-enriched medium without CF-BD (Fig. 4g, independent-sample t test, t = 4.217, df = 4, p = 0.014), which suggests that trichlorphon was degraded by the bacterium. Furthermore, GC-MS analysis revealed that trichlorphon was degraded into chloral hydrate and dimethyl phosphite, which are significantly less toxic than trichlorphon [36, 37] (Fig. 4h).
Predicting trichlorphon degradation ability through a genome analysis of CF-BD
To predict the CF-BD trichlorphon degradation ability, the CF-BD genome was further analyzed. The draft genome of CF-BD contains one scaffold with 5,098,016 bp, 4762 predicted coding sequences (CDS), and 25 rRNA and 85 tRNA genes. Putative functions were assigned to 4689 of the predicted CDS. The average G + C content was 52.06%. Among the function-annotated genes, 55 phosphatase genes were identified. Comparative genomics that employed the sequences of the fully available genomes of C. rodentium, C. koseri, C. freundii, and CF-BD revealed large clusters that were only found in CF-BD (Fig. 5a, c). Also, the genomic components of CF-BD were different from those of C. rodentium, C. koseri, and C. freundii (Fig. 5b). In particular, the clusters highlighted at positions I, II, and III contained five phosphatase genes (gene ID 0086, 1012, 1747, 2752. and 4498) that were specifically present in CF-BD (Fig. 5c). Analysis of sequence similarities suggested that these five phosphatase genes were more similar to the organophosphorus hydrolase genes (OPH) of other bacteria (Fig. 5d), which assigned potential function to degrade the organophosphorus to the five genes.
In order to verify the function of the five identified OPH-like genes, the expression of these genes was measured with the CF-BD exposed to trichlorphon at different concentrations. The results showed the expression of genes 0086, 1012, 1747, and 4498 were significantly higher when exposed to trichlorphon (0086 F = 13.995, df = 34, p < 0.001; 1012 F = 4.36, df = 34, p = 0.011; 1747 F = 9.291, df = 34, p < 0.001; 4498 F = 21.092, df = 34, p < 0.001; Fig. 6).
We report the first isolation of Citrobacter freundii strain (CF-BD) associated with trichlorphon resistance from the digestive tract of the significant pest fruit fly B. dorsalis. Using 16S rDNA and genome sequencing, Gram staining, physiological and biochemical identification procedures, and fluorescence in situ hybridization, we confirmed that some of the bacteria residing in the midgut crypts of the fly are similar to C. freundii and enhance the fly’s insecticide resistance. With the advent of sequencing technology, many bacterial communities have been identified in the gut of various insects [38–40], including of tephritid fruit flies [11, 41, 42]. Citrobacter is a commonly found symbiotic taxon in insects and fruit flies specifically  and belongs to the Gammaproteobacteria, a class that includes dominant symbiotic bacteria in many insect lineages [44, 45]. We determined that the dominant symbiotic bacteria in the gut of oriental fruit flies were in the family of Enterobacteriaceae and, in particular, in the genus Klebsiella (Fig. 1d), suggesting that these bacterial symbionts may play a role in the biology of this fly species. CF-BD was also present in all flies at variable and lower densities. Citrobacter freundii includes pathogenic isolates that can cause respiratory and urinary tract infections in humans ; strains of this species can also be pathogenic to fish . In contrast, oriental fruit flies do not exhibit pathology associated with CF-BD, and an increased titer of this bacterium was found in the insecticide resistant strain (RS), indicating a different function of this bacterium in the gut of resistant versus susceptible flies. Moreover, several genomic differences were found between CF-BD and other C. freundii strains, suggesting that CF-BD is not necessarily pathogenic to humans or fish.
Growth on trichlorphon-enriched plates and the degradation of trichlorphon in the mineral medium demonstrated that CF-BD can degrade trichlorphon (Fig. 4e–h). So far, many bacterial species have been shown to possess this ability, although these bacteria were largely isolated from environmental samples such as from soil [28, 48]. We determined that CF-BD isolated from fruit flies was able to degrade trichlorphon and enhance the resistance of the flies to trichlorphon. Similar results have been obtained for a bacterium isolated from R. pomonella  and stinkbugs, in particular, the bean bug Riptortus pedestris . In R. pedestris, resistance has been observed with both oral and percutaneous applications of the insecticide. In this stinkbug species, gut symbiotic bacteria of the genus Burkholderia can be acquired by nymphs from soil. While our study demonstrated that CF-BD in B. dorsalis enhances resistance to trichlorphon, it remains unknown how CF-BD is transmitted through life stages and in fly populations. Our field survey detected CF-BD at low densities in all tested adult flies; however, it is not clear whether CF-BD is maternally inherited, horizontally transmitted between larvae or adults, or taken up from the environment.
Microbial symbionts of insect species can have important functions in their hosts. In bark beetle, symbiotic bacteria can supply essential nutrients and assist with digestion and detoxification of plant compounds , and this was also seen in the coffee berry borer . In aphids, Hamiltonella defensa  can increase resistance to hymenopteran parasitoids [6, 51] via toxin-encoding bacteriophages  that kill parasitoid eggs or early larvae. Some microbial symbionts are able to assist the host in escaping predation  and increase host resistance to Bacillus thuringiensis . Furthermore, microorganisms present in the gut can provide additional functions, including the digestion of plant polymers and the detoxification of plant-produced toxins . In bean bugs, Burkholderia was rapidly established within a single insect generation after environmental acquisition and shared horizontally between different individuals . Although we have demonstrated that CF-BD can be obtained by oral feeding, additional evidence is needed to determine whether it can be transmitted horizontally in a manner similar to Burkholderia.
Many organophosphate pesticide-degrading bacteria have been reported in recent years [48, 56, 57]. Moreover, many Citrobacter species can degrade phenols and phosphates [58–60], and their degradation pathways are widely recognized. Myo-inositol hexakisphosphate can be degraded by a phytase from Citrobacter braakii , and the degradation of m-cresol via the ortho cleavage pathway by Citrobacter farmeri SC01 has also been reported . In our study, we isolated CF-BD from the gut of B. dorsalis, although its organophosphate-degrading ability has previously not been reported. Genome sequence analysis of CF-BD revealed 55 phosphatase genes, and within these, five new genes were unique to CF-BD. Sequence similarities of the five genes revealed that they were OPH-like, and their expression was stimulated by the presence of trichlorphon. We also found that trichlorphon was degraded into chloral hydrate and dimethyl phosphite by CF-BD. Although the characteristics of these metabolites are extremely similar to results obtained with C. braakii, additional evidence is needed to determine whether this degradation pathway is the same as that utilized by C. braakii and whether phosphatase is the critical enzyme in this pathway. We only sequenced the C. freundii genome in RS flies, which were collected in 2008 and then experienced serial selection for insecticide resistance in the laboratory. Whether the CF-BD genome has changed between the field collection and after selection for insecticide resistance remains unknown. Citrobacter freundii in SS flies may have lower insecticide degradation capacity than CF-BD of RS flies, and future research should address this, in combination with comparative genomics approaches. Until this has been tested further, we conclude that the enhanced insecticide resistance in B. dorsalis in our study is due to increased abundance of CF-BD in the gut.
In this study, we isolated a bacterium with an important role in the resistance of B. dorsalis to trichlorphon, and this symbiont was widespread in the wild population of the flies (Fig. 3); thus, regulation of this symbiont might be useful for managing trichlorphon resistance. Moreover, our findings also demonstrate the necessity of considering the gut symbiotic bacteria of insects during the implementation of novel pest control measures . Because the generation times of bacteria are considerably shorter than those of the host insects, the evolution of insecticide resistance in insects may proceed much more rapidly via symbiont-mediated processes.
We document a novel mechanism of insecticide resistance in which a gut bacterium (Citrobacter freundii) of Bactrocera dorsalis enhances the fruit flies’ resistance to trichlorphon, and we experimentally demonstrated that C. freundii degrades trichlorphon into chloral hydrate and dimethyl phosphite. These results indicate that symbiont-mediated trichlorphon resistance might develop in oriental fruit flies.
The susceptible strain (SS) of B. dorsalis was collected from a carambola (Averrhoa carambola) orchard in Guangzhou, Guangdong Province, in April 2008, and was maintained in the laboratory for approximately 80 generations under the following conditions: 25 ± 1 °C; 16:8 h light:dark cycle; 70–80% relative humidity (RH); a maize-based artificial larval diet containing 150 g corn flour, 150 g banana, 0.6 g sodium benzoate, 30 g yeast, 30 g sucrose, 30 g paper towel, 1.2 mL hydrochloric acid and 300 mL water, and an adult diet consisting of water, yeast hydrolysate, and sugar. Pesticide exposure was avoided during rearing.
A resistant strain (RS) of B. dorsalis was obtained by selection after adult exposure to a trichlorphon treated surface over the course of 33 generations. Trichlorphon was diluted with acetone to LC50 for each generation and evenly coated onto the inside of a 250-mL conical flask by shaking. After the acetone had volatilized, 30 3–5-day-old fly pairs were placed in the flask for 24 h. The surviving flies were selected for breeding of the next generation . The resistance levels were calculated based on Benson’s formula .
Fly toxicity assay
A total of 15 3–5-day-old mixed (female and male) adult flies were placed in plastic cups that had six intruding 200-μL pipette tips on the top and ventilation holes on the sides (Additional file 5: Figure S1). The tips contained aqueous trichlorphon solutions at different concentrations: 0 mg/L, 0.625 × 10− mg/L, 1.25 × 10−3 mg/L, 2.5 × 10−3 mg/L, 5 × 10−3 mg/L, 10 × 10−3 mg/L, and 20 × 10−3 mg/L for SS; 0 mg/L, 25 × 10−3 mg/L, 50 × 10−3 mg/L, 100 × 10−3 mg/L, 200 × 10−3 mg/L, 400 × 10−3 mg/L, and 800 × 10−3 mg/L for RS. Three replicate cups with six tips each were prepared for each concentration. All of the treatments were maintained at a temperature of 25 ± 1 °C under a 16:8 h light:dark cycle. Mortality was recorded after 24 h.
Extraction of DNA from SS and RS fly guts
Three-to-five-day-old mixed (female and male) adult flies were selected and soaked in absolute ethanol for 3 min. The soaked flies were dissected under a stereo-microscope, and the guts were transferred into centrifuge tubes containing DNA extraction buffer. For each sample, 15 flies were dissected, and three samples were collected for both the SS (SS1, SS2, SS3) and the RS (RS1, RS2, RS3) strains. Total DNA of the dissected gut samples was extracted using a DNA extraction kit (Tiangen, Beijing, China) following the manufacturer’s instructions.
V3 + V4 region of 16S rDNA amplification and sequencing
Approximately 465 bp of the V3 + V4 region of the bacterial 16S rDNA gene was amplified by PCR according to a standard protocol (Additional file 6: Supplementary methods). And the amplified DNA was sequenced using the Illumina sequencing kit and the Illumina MiSeq sequencer (Illumina, San Diego, CA, USA).
Bioinformatic analysis of sequencing results
After quantity control (Additional file 6: Supplementary methods), sequence reads were subjected to redundancy treatment with Mothur software  to count the number of identical tags. For species annotation, an RDP classifier  was used with naïve Bayesian settings; the confidence threshold was set to 0.5. To obtain additional information regarding species diversity composition, we subjected the tags to OTU abundance analysis (Additional file 6: Supplementary methods).
Citrobacter sp. isolation and culture
Six 3–5-day-old adult fruit flies (three flies from each strain) were collected and immediately soaked in 70% ethanol for 3 min to remove surface bacteria. The guts of the flies were dissected and collected in six sterile centrifuge tubes to which 20 μL of sterile water was added. The guts were then ground with sterile grinding pestles, and the fluid was streaked and cultivated for 24 h at 30 °C on brain heart infusion (BHI) agar flat plates, which is the specific medium for cultivation of Citrobacter sp. Colonies with the same morphology were selected for subculturing. The pure cultures were inoculated into BHI medium, and the liquid cultures were stored in 25% glycerol solution under −80 °C.
Bacteria identification by 16S rDNA amplification
16S rDNA of the cultivated bacteria was amplified sequenced (Additional file 6: Supplementary methods). The sequences were subjected to a BLAST search against the NCBI database for sequence homology analysis.
Physiological and biochemical identification of the bacterium and Gram staining
The physiological and biochemical characteristics of the Citrobacter isolates were tested using GYZ-15 eV biochemical detection kits produced by Guangzhou Huankai Microbial Sci. & Tech. Co., Ltd. (Guangzhou, China) according to the manufacturer’s instructions. Gram staining kits produced by Guangzhou Huankai Microbial Sci. & Tech. Co., Ltd. were used to determine bacterial morphology under a light microscope.
Antibiotic and trichlorphon sensitivity testing
To test bacterial sensitivities to antibiotics, 10-μL inocula were streaked onto BHI agar flat plates to which drug susceptibility test papers were subsequently attached (Table 1). The plates were incubated for 24 h at 30 °C, and the diameters of the inhibition zones were then measured. According to CLSI standards, the cultures were classified as sensitive (S), medium (I), or resistant (R). Plates containing 1 × 10−3 mg/L trichlorphon were also prepared to test the sensitivity of the bacteria to trichlorphon with a 10-μL inoculum and an incubation period of 24 h. As a control, plates without trichlorphon were prepared and incubated with inocula. The diameters of colonies were measured after 24 h.
Resistance of flies inoculated with CF-BD
The bacterial inoculum for flies was prepared by selecting and incubating a colony of CF-BD in BHI medium at 30 °C until an OD600 of 0.8. Next, >300 newly emerged SS flies were fed sugar and bacterial inoculum instead of water for 3 days. After 3 days, the guts of five flies were dissected, ground, and diluted with 200 μL sterile water in 1.5-mL centrifuge tubes. Next, 10 μL of the fluid was plated and cultivated for 12 h, and the numbers of colonies were recorded and compared with the SS flies. Then, cohorts of 15 CF-BD treated 3-day old SS flies each were collected and fed a trichlorphon solution to measure the LC50. The tested trichlorphon concentrations were 0 mg/L, 0.625 × 10−3 mg/L, 1.25 × 10−3 mg/L, 2.5 × 10−3 mg/L, 5 × 10−3 mg/L, and 10 × 10−3 mg/L. After 24 h, the number of dead flies was recorded. As control, newly emerged SS flies were fed pure water and after 3 days, provided trichlorphon solutions in cohorts of 15 flies each. For each trichlorphon concentration, three replicates were performed.
Resistance of flies cleared of CF-BD
More than 300 newly emerged RS flies were fed a streptomycin solution for 2 days to clear them of bacteria. Afterwards, the guts of five flies were dissected, ground, and diluted with 200 μL of sterile water; the sample was used to cultivate bacteria on plates; the number of colonies was recorded and compared with that of the control. Then, cohorts of 15 RS flies each were fed with a trichlorphon solution, and the number of dead flies was counted after 24 h. The trichlorphon concentrations were 0 mg/L, 25 × 10−3 mg/L, 50 × 10−3 mg/L, 100 × 10−3 mg/L, and 200 × 10−3 mg/L. As control, newly emerged RS flies were fed pure water and, after 3 days, provided trichlorphon in cohorts of 15 flies each. For each trichlorphon concentration, three replicates were performed.
Testing of the trichlorphon degradation characteristics of purified bacteria
Mineral medium (MM) was prepared containing the following salts (g/L): NaCl 1.00, (NH4)2SO4 1.00, K2HPO4 1.50, KH2PO4 0.50, and MgSO4·7H2O 0.50 (pH: 7.0–7.5). The enrichment medium was prepared by adding trichlorphon to MM at a concentration of 100 mg/L. To measure degradation by the isolated bacteria, a 2-mL inoculum was added to MM (100 mL); 2 mL pure water added to the MM (100 mL) was used as the control. Three replicates were prepared for each treatment. The inoculated media were shaken and cultivated at 30 °C for 24 h. Then, 2 mL of the cultures was collected for quantitative analyses of trichlorphon. To purify trichlorphon, 4 mL acetone was added to the culture, which was then shaken for 1 h. NaCl was added until saturation; the solution was shaken for 1 min and allowed to stand for 5 min at room temperature before the acetone layer was collected and filtered through a 45-μm bacterial membrane filter. The filtrate was dried with anhydrous sodium sulfate and used for gas chromatography. The detection conditions were as follows: air flow 5 mL/min, hydrogen flow 5.5 mL/min, carrier gas (nitrogen) flow 8 mL/min, detector temperature 230 °C, injection port temperature 200 °C, chromatographic column temperature 150 °C, and injection volume 0.5 μL. Acetone was used as the solvent, and the peak area was used to quantify trichlorphon. We further examined the filtrate by GC-MS to identify the degradation products.
Fluorescence in situ hybridization
To locate the CF-BD in the B. dorsalis gut, FISH was conducted for RS flies. The oligonucleotide probe called the CFBD-probe (5′-AATGGCGTACACAAAGAG-3′) which was used for in situ hybridization, was labeled with Cy3 at the 5′ end. The probe specifically targeted different regions in the 16S rRNA of the CF-BD symbiont, and it was simultaneously used to enhance the hybridization signals. The dissected gut samples were incubated in a hybridization buffer [20 mM Tris-HCl (pH 8.0), 0.9 M NaCl, 0.01% sodium dodecyl sulfate (SDS), and 30% formamide] containing 50 nM probe. Following overnight incubation, the samples were thoroughly washed in phosphate-buffered saline (PBS) and mounted in SlowFade antifade solution (Molecular Probes, Chuo-ku, Japan). The samples were observed under an epifluorescent microscope (Axiophot, Carl Zeiss, Shinjuku-ku, Japan).
CF-BD genome sequencing and annotation
To investigate the function of CF-BD at the genomic level, the genome of CF-BD was sequenced. After genome assembly and comparative genome analyses, potential organophosphorus hydrolase genes (OPH) were identified and compared with other similar genes from other bacteria (Additional file 6: Supplementary methods).
Expression analysis of OPH-like genes by real-time qPCR
To identify the potential function of the OPH-like genes, total RNA of CF-BD cultivated for 20 h in the BHI medium containing 0, 12.5, 25, and 50 mg trichlorphon were extracted to analyze the expression of the OPH-like genes. Complementary DNA (cDNA) was reverse-transcribed from 2 μg total RNA using MMLV reverse transcriptase (Promega). The recA gene was used as the reference gene . Real-time quantitative PCR amplification was performed using Mx3000P spectrofluorometric thermal cycler (Agilent Technologies, Santa Clara, CA, USA) and Real Master Mix (SYBR Green) kit (Tiangen), starting with a 2 min incubation at 95 °C, followed by 40 cycles of 95 °C, 20 s; 55 °C, 1 min; and 72 °C, 30 s. Primer information for the genes was described in Additional file 3: Table S3.
Detection of CF-BD in different B. dorsalis populations
To detect CF-BD infection in different B. dorsalis populations, 78 flies collected from 13 populations in 2015 (Additional file 7: Table S5) were subjected to real-time quantitative PCR with CF-BD-specific primers in order to obtain the cycle threshold (Ct) value for each fly. A 282-bp region of the recA gene in CF-BD was amplified by the specific primers (Additional file 3: Table S3) under a temperature profile of 95 °C for 5 min followed by 35 cycles at 95 °C for 30 s, 55 °C for 30 s and 72 °C for 30 s, and ending at 72 °C for 5 min.
Mortality data for the SS and RS flies were corrected using Abbott’s formula , and LC50 values were examined with a probit analysis conducted with SPSS (Statistical Package for the Social Sciences) 16.0 software. Variation in the OTU profile similarities was visualized via 2D, non-metric multidimensional scaling (nMDS) plots and was statistically evaluated using nonparametric MANOVA (NPMANOVA). Differences between the treatments and controls were compared with independent-sample t tests. Differences were considered significant when the p values were <0.05. The data were analyzed using SPSS.
- LC50 :
The median lethal concentration
Bright M, Bulgheresi S. A complex journey: transmission of microbial symbionts. Nat Rev Micro. 2010;8(3):218–30.
Koga R, Moran NA. Swapping symbionts in spittlebugs: evolutionary replacement of a reduced genome symbiont. ISME J. 2014;8(6):1237–46.
Cariveau DP, Powell JE, Koch H, Winfree R, Moran NA. Variation in gut microbial communities and its association with pathogen infection in wild bumble bees (Bombus). ISME J. 2014;8(12):2369–79.
Stursova M, Snajdr J, Cajthaml T, Barta J, Santruckova H, Baldrian P. When the forest dies: the response of forest soil fungi to a bark beetle-induced tree dieback. ISME J. 2014;8(9):1920–31.
Douglas AE. The microbial dimension in insect nutritional ecology. Funct Ecol. 2009;23(1):38–47.
Oliver KM. Variation in resistance to parasitism in aphids is due to symbionts not host genotype. Proc Natl Acad Sci U S A. 2005;102:12795–800.
Sharon G, Segal D, Ringo JM, Hefetz A, Zilber-Rosenberg I, Rosenberg E. Commensal bacteria play a role in mating preference of Drosophila melanogaster. Proc Natl Acad Sci U S A. 2010;107(46):20051–6.
Desai MS, Brune A. Bacteroidales ectosymbionts of gut flagellates shape the nitrogen-fixing community in dry-wood termites. ISME J. 2012;6(7):1302–13.
Senderovich Y, Halpern M. The protective role of endogenous bacterial communities in chironomid egg masses and larvae. ISME J. 2013;7(11):2147–58.
Tsuchida T, Koga R, Fukatsu T. Host plant specialization governed by facultative symbiont. Science. 2004;303(5666):1989.
Morrow JL, Frommer M, Shearman DC, Riegler M. The microbiome of field-caught and laboratory-adapted Australian tephritid fruit Fly species with different host plant use and specialisation. Microb Ecol. 2015;70(2):498–508.
Dunbar HE, Wilson AC, Ferguson NR, Moran NA. Aphid thermal tolerance is governed by a point mutation in bacterial symbionts. PLoS Biol. 2007;5(5):e96.
Genta FA, Dillon RJ, Terra WR, Ferreira C. Potential role for gut microbiota in cell wall digestion and glucoside detoxification in Tenebrio molitor larvae. J Insect Physiol. 2006;52(6):593-601.
Lauzon CR, Potter S, Prokopy RJ. Degradation and detoxification of the dihydrochalcone phloridzin by Enterobacter agglomerans, a bacterium associated with the apple pest, Rhagoletis pomonella (Walsh) (Diptera : Tephritidae). Environ Entomol. 2003;32(5):953–62.
Ben-Yosef M, Pasternak Z, Jurkevitch E, Yuval B. Symbiotic bacteria enable olive fly larvae to overcome host defences. R Soc Open Sci. 2015;2(7):150170.
Cejanavarro JA, Vega FE, Karaoz U, Zhao H, Jenkins S, Lim HC, Kosina P, Infante F, Northen TR, Brodie EL. Gut microbiota mediate caffeine detoxification in the primary insect pest of coffee. Nat Commun. 2015;6:7618.
Boush MG, Matsumura F. Insecticidal degradation by Pseudomonas melophthora, the bacterial symbiote of the apple maggot. J Econ Entomol. 1967;60:918–920(913).
Kikuchi Y, Hayatsu M, Hosokawa T, Nagayama A, Tago K, Fukatsu T. Symbiont-mediated insecticide resistance. Proc Natl Acad Sci U S A. 2012;109(22):8618–22.
Kikuchi Y, Hosokawa T, Fukatsu T. An ancient but promiscuous host—symbiont association between Burkholderia gut symbionts and their heteropteran hosts. ISME J. 2011;5(3):446–60.
Perry AS, Yamamoto I, Ishaaya I, Perry R. Toxicology of Insecticides. In: Insecticides in agriculture and environment. Berlin Heidelberg: Springer; 1998. p. 11.
Roush RT, McKenzie JA. Ecological genetics of insecticide and acaricide resistance. Annu Rev Entomol. 1987;32:361–80.
Devine G. Global Pesticide Resistance in Arthropods - By M. E. Whalon, D. Mota-Sanchez & R. M. Hollingworth. Entomol Exp Appl. 2009;131(1):106-106.
Robertson BK, Alexander M. Growth-linked and cometabolic biodegradation: possible reason for occurrence or absence of accelerated pesticide biodegradation. Pestic Sci. 1994;41(4):311–8.
Sonia Rodríguez-Cruz M, Jones JE, Bending GD. Field-scale study of the variability in pesticide biodegradation with soil depth and its relationship with soil characteristics. Soil Biol Biochem. 2006;38(9):2910–8.
Felsot AS. Enhanced biodegradation of insecticides in soil: implications for agroecosystems. Annu Rev Entomol. 1989;34(1):453–76.
Arbeli Z, Fuentes CL. Accelerated biodegradation of pesticides: an overview of the phenomenon, its basis and possible solutions; and a discussion on the tropical dimension. Crop Prot. 2007;26(12):1733–46.
Ortiz-Hernández ML, Quintero-Ramírez R, Nava-Ocampo AA, Bello-Ramírez AM. Study of the mechanism of Flavobacterium sp. for hydrolyzing organophosphate pesticides. Fundam Clin Pharmacol. 2003;17(6):717–23.
Singh BK, Walker A. Microbial degradation of organophosphorus compounds. FEMS Microbiol Rev. 2006;30(3):428–71.
Yee WL, Goughnour RB. Host plant use by and new host records of apple maggot, western cherry fruit fly, and other Rhagoletis species (Diptera: Tephritidae) in western Washington state. Pan Pac Entomol. 2008;84(3):179–93.
Hollingsworth RG, Drew RAI, Allwood AJ, Romig M, Vagalo M, Tsatsia F. Host plants and relative abundance of fruit fly (Diptera: Tephritidae) species in the Solomon Islands. Aust J Entomol. 2003;42(2):95–108.
Clarke AR, Armstrong KF, Carmichael AE, Milne JR, Raghu S, Roderick GK, Yeates DK. Invasive phytophagous pests arising through a recent tropical evolutionary radiation: the Bactrocera dorsalis complex of fruit flies. Annu Rev Entomol. 2005;50(1):293–319.
Zhang YP, Zeng L, Yong-Yue LU, Liang GW. Genetic analysis of Bactrocera dorsalis resistance to trichlorphon. J South China Agric Univ. 2008;29:39–43.
Zhang YP, Zeng L, Yong-Yue LU, Liang GW. Resistance stability and re-growth in adults of the oriental fruit fly, Bactrocera dorsalis(Diptera:Tephritidae) to trichlorphon. Acta Entomol Sin. 2008;51(10):1044–9.
Jin T, Zeng L, Lin YY, Lu YY, Liang GW. Characteristics of protein variants in trichlorphon-resistant Bactrocera dorsalis (Diptera; Tephritidae) larvae. Genet Mol Res. 2012;11(3):2608–19.
Singh BK. Organophosphorus-degrading bacteria: ecology and industrial applications. Nat Rev Microbiol. 2009;7(2):156–64.
Kuz’minov BP, Kokot VR, Sharova TI, Zhuk S, Vus MM. Experimental data on the hygienic regulation of dimethyl phosphite in the air of a work area. Gig Tr Prof Zabol. 1992;4:22–3.
Sing K, Erickson T, Amitai Y, Hryhorczuk D. Chloral hydrate toxicity from oral and intravenous administration. J Toxicol Clin Toxicol. 1996;34(1):101–6.
Hu X, Wang C, Chen H, Ma J. Differences in the structure of the Gut bacteria communities in development stages of the Chinese white pine beetle (dendroctonus armandi). Int J Mol Sci. 2013;14(10):21006–20.
Fang W, Fang Z, Liu Z, Yuan J, Zhang X, Peng H, Hong Y, Xiao Y. Phylogenetic analysis of bacterial community in the gut of American cockroach (Periplaneta americana). Acta Microbiol Sin. 2013;53(9):984–94.
Dillon RJ, Webster G, Weightman AJ, Dillon VM, Blanford S, Charnley AK. Composition of Acridid gut bacterial communities as revealed by 16S rRNA gene analysis. J Invertebr Pathol. 2008;97(3):265–72.
Aharon Y, Pasternak Z, Ben Yosef M, Behar A, Lauzon C, Yuval B, Jurkevitch E. Phylogenetic, metabolic, and taxonomic diversities shape mediterranean fruit fly microbiotas during ontogeny. Appl Environ Microbiol. 2013;79(1):303–13.
Andongma AA, Wan L, Dong YC, Li P, Desneux N, White JA, Niu CY. Pyrosequencing reveals a shift in symbiotic bacteria populations across life stages of Bactrocera dorsalis. Sci Rep. 2015;5:9470.
Behar A, Jurkevitch E, Yuval B. Bringing back the fruit into fruit fly—bacteria interactions. Mol Ecol. 2008;17(5):1375–86.
Dunn AK, Stabb EV. Culture-independent characterization of the microbiota of the ant lion Myrmeleon mobilis (Neuroptera: Myrmeleontidae). Appl Environ Microbiol. 2006;71:8784–94.
Babendreier D, Joller D, Romeis J, Bigler F, Widmer F. Bacterial community structures in honeybee intestines and their response to two insecticidal proteins. FEMS Microbiol Ecol. 2007;59(3):600–10.
Long-Hua HU, Xue-Fei HU, Xiong JQ, Jia KR, Zhang BL, Xiao-Yan HU. Clinical distribution and characteristics of drug resistance of Citrobacter freundii. Chin J Microbiol Immunol. 2010;22:724–9.
Joh SJ, Ahn EH, Lee HJ, Shin GW, Kwon JH, Park CG. Bacterial pathogens and flora isolated from farm-cultured eels (Anguilla japonica) and their environmental waters in Korean eel farms. Vet Microbiol. 2013;163(1–2):190–5.
Zhang YH, Xu D, Liu JQ, Zhao XH. Enhanced degradation of five organophosphorus pesticides in skimmed milk by lactic acid bacteria and its potential relationship with phosphatase production. Food Chem. 2014;164:173–8.
Adams AS, Currie CR, Cardoza Y, Klepzig KD, Raffa KF. Effects of symbiotic bacteria and tree chemistry on the growth and reproduction of bark beetle fungal symbionts. Can J For Res. 2009;39(6):1133–47.
Moran NA, Russell JA, Koga R, Fukatsu T. Evolutionary relationships of three new species of Enterobacteriaceae living as symbionts of aphids and other insects. Appl Environ Microbiol. 2005;71(6):3302–10.
Ferrari J, Darby AC, Daniell TJ, Godfray HCJ, Douglas AE. Linking the bacterial community in pea aphids with host-plant use and natural enemy resistance. Ecol Entomol. 2004;29(1):60–5.
Degnan PH, Moran NA. Diverse phage-encoded toxins in a protective insect endosymbiont. Appl Environ Microbiol. 2008;74(21):6782–91.
Kaltenpoth M. Actinobacteria as mutualists: general healthcare for insects? Trends Microbiol. 2009;17(12):529–35.
Patil C, Borase H, Salunke B, Patil S. Alteration in Bacillus thuringiensis toxicity by curing gut flora: novel approach for mosquito resistance management. Parasitol Res. 2013;112(9):3283–8.
Hansen AK, Moran NA. The impact of microbial symbionts on host plant utilization by herbivorous insects. Mol Ecol. 2014;23(6):1473–96.
Siripattanakul-Ratpukdi S, Vangnai A, Sangthean P, Singkibut S. Profenofos insecticide degradation by novel microbial consortium and isolates enriched from contaminated chili farm soil. Environ Sci Pollut Res. 2015;22(1):320–8.
Akbar S, Sultan S, Kertesz M. Determination of cypermethrin degradation potential of soil bacteria along with plant growth-promoting characteristics. Curr Microbiol. 2015;70(1):75–84.
Pontoppidan K, Glitsoe V, Guggenbuhl P, Quintana AP, Nunes CS, Pettersson D, Sandberg AS. In vitro and in vivo degradation of myo-inositol hexakisphosphate by a phytase from Citrobacter braakii. Arch Anim Nutr. 2012;66(6):431–44.
Ren Y, Peng L, Zhao G, Wei C. Degradation of m-cresol via the ortho cleavage pathway by Citrobacter farmeri SC01. Biochem Eng J. 2014;88:108–14.
Ma Y, Wang L, Liu L, Zhang X. Biodegradation of tylosin residue in pharmaceutical solid waste by a novel Citrobacter amalonaticus strain. Environ Prog Sustain Energy. 2015;34(1):99–104.
Engel P, Moran NA. The gut microbiota of insects—diversity in structure and function. FEMS Microbiol Rev. 2013;37(5):699–735.
Jin T, Zeng L, Lin Y, Lu Y, Liang G. Insecticide resistance of the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), in mainland China. Pest Manag Sci. 2011;67(3):370–6.
Benson RL. On the necessity of controlling the level of insecticide resistance in insect populations. Bioscience. 1971;21(23):1160–5.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.
Huse SM, Welch DM, Morrison HG, Sogin ML. Ironing out the wrinkles in the rare biosphere through improved OTU clustering. Environ Microbiol. 2010;12(7):1889–98.
Takle GW, Toth IK, Brurberg MB. Evaluation of reference genes for real-time RT-PCR expression studies in the plant pathogen Pectobacterium atrosepticum. BMC Plant Biol. 2007;7(1):50.
Abbott WS. A method of computing the effectiveness of an insecticide. J Econ Entomol. 1925;18:265–76.
We thank Gary W. Felton at Penn State University and Tao Jin at Chinese Academy of Tropical Agricultural Sciences for their constructive comments regarding this manuscript.
This study was supported by the Science Foundation for the Excellent Youth Scholars of Guangdong Province (No. Yq2013031) and the National Key Research and Development Project (2016YC1201200). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
Sequence data has been deposited at NCBI under Bioproject PRJNA357667.
YX conceived and designed the experiments. DC and ZG analyzed the data. GL and YX contributed the reagents/materials/analysis tools. DC, MR, ZX, and YX wrote the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Primer information of the genes. (DOCX 17 kb)
Prevalence of CF-BD in individual field flies. (DOCX 18 kb)
Tag data. (DOCX 17 kb)
Numbers of OTUs in different samples. (DOCX 16 kb)
Device used for testing the toxicity of trichlorphon to flies. (DOCX 103 kb)
Supplementary methods. (DOCX 27 kb)
Antibiotic sensitivity of CF-BD. (DOCX 17 kb)
Gene information used for phylogenetic tree construction. (DOCX 20 kb)