Introduction
The oriental fruit fly Bactrocera dorsalis (Hendel) is a significant invasive fruit fly pest in parts of Asia and Africa and can cause considerable economic losses (Clarke et al., Reference Clarke, Armstrong, Carmichael, Milne, Raghu, Roderick and Yeates2005; Schutze et al., Reference Schutze, Mahmood, Pavasvici, Bo, Newman, Clarke, Krosch and Cameron2015). The damage caused by this fruit fly can be attributed to its high fecundity and wide host range (Li & Ye, Reference Li and Ye2000). Much attention has been paid to the population genetics of B. dorsalis to uncover the origin and invasion pathways of this pest in China (Aketarawong et al., Reference Aketarawong, Bonizzoni, Thanaphum, Gomulski, Gasperi and Malacrida2007; Li et al., Reference Li, Yang, Tang, Zeng and Liang2007; Shi et al., Reference Shi, Kerdelhue and Ye2010). The invasions of B. dorsalis in China can be divided into two branches: one originated from South Asia, and the other originated from Southeast Asia. The Southeast Asia-originated fruit flies landed in China at Yunnan Province and then expanded to Sichuan, Chongqing, Guizhou, and Guangxi Provinces. The South Asia-originated flies may have landed at Guangdong and Guangxi Provinces and then expanded to south and east China, including Fujian, Zhejiang, and Jiangsu Provinces and Shanghai city (Li et al., Reference Li, Yang, Tang, Zeng and Liang2007, Reference Li, Wu, Chen, Wu and Li2012; Shi et al., Reference Shi, Kerdelhue and Ye2010).
Many studies have attempted to investigate the associated bacteria of fruit fly species (Tsiropoulos, Reference Tsiropoulos1983; Bextine et al., Reference Bextine, Lampe, Lauzon, Jackson and Miller2005; Capuzzo et al., Reference Capuzzo, Firrao, Mazzon, Squartini and Girolami2005; Konstantopoulou et al., Reference Konstantopoulou, Raptopoulos, Stavrakis and Mazomenos2005; Mazzon et al., Reference Mazzon, Piscedda, Simonato, Martinez-Sañudo, Squartini and Girolami2008, Reference Mazzon, Martinez-Sanudo, Simonato, Squartini, Savio and Girolami2010; Sacchetti et al., Reference Sacchetti, Granchietti, Landini, Viti, Giovannetti and Belcari2008; Kounatidis et al., Reference Kounatidis, Crotti, Sapountzis, Sacchi, Rizzi and Chouaia2009; Morrow et al., Reference Morrow, Frommer, Shearman and Riegler2015a , Reference Morrow, Frommer, Royer, Shearman and Riegler b ). Previous studies have shown that gut bacteria play an important role in the development and reproduction of their host fruit flies (Lauzon, Reference Lauzon, Bourtzis and Miller2003). Firstly, Tephritidae mainly feed on nitrogen-poor fruits and are unable to synthesize essential amino acid by themselves. Associated bacteria can participate in the carbon and nitrogen cycle of the host fruit flies and help synthesize essential amino acids and minerals (Miyazaki et al., Reference Miyazaki, Boush and Baerwald1968; Lauzon et al., Reference Lauzon, Sjogren and Prokopy2000, Reference Lauzon, Bussert, Sjogren and Prokopy2003). Associated Enterobacter agglomerans of the Mexican fruit fly Anastrepha ludens (Loew) can produce the volatile 2, 5-dimethyl pyrazine, which can attract the fly. This attractiveness could be the second function of associated bacteria (Robacker et al., Reference Robacker, Lauzon and He2004, Reference Robacker, Lauzon, Patt, Margara and Sacchetti2009). Thirdly, E. agglomerans isolated from apple maggot fly Rhagoletis pomonella (Walsh) can also help flies to degrade the toxic purine compounds of host plants (Lauzon et al., Reference Lauzon, Sjogren and Prokopy2000). Last but not least, intestinal associated bacteria can also prevent the colonization of pathogenic bacteria in Ceratitis capitata (Wiedemann) and improve its fitness (Behar et al., Reference Behar, Yuval and Jurkevitch2008b ). Besides associated bacteria, some tephritids have bacterial symbionts that are more tightly associated with their host, which means there exist potential co-evolutionary interactions between these bacteria and its host. The bacteria found in olive fly (Bactrocera oleae) and Tephritinae subfamily flies, developing in composite flower heads, is a good example (Mazzon et al., Reference Mazzon, Piscedda, Simonato, Martinez-Sañudo, Squartini and Girolami2008, Reference Mazzon, Martinez-Sanudo, Simonato, Squartini, Savio and Girolami2010).
Several reports have analyzed the associated bacteria of B. dorsalis. Volatile 2-butanone and amino compounds produced by associated Citrobacter freundii, Klebsiella pneumoniae and Enterobacter cloacae, can attract B. dorsalis to feed on the special host plant (Jang & Nishijima, Reference Jang and Nishijima1990). The presence of Wolbachia infection in B. dorsalis, even if with low infection rate, was detected by Sun et al. (Reference Sun, Cui and Li2007). The gut bacterial diversity from three populations of B. dorsalis, including laboratory-reared, laboratory sterile sugar-reared, and field-collected populations were studied and compared (Wang et al., Reference Wang, Jin and Zhang2011). In order to identify the cultivable bacteria for developing bacterial biocontrol agents, the attractive potential of the cultivable bacteria in the intestine of the three populations was also discussed (Wang et al., Reference Wang, Jin, Peng, Zhang, Chen and Hua2014). The bacteria associated with different developmental stages of B. dorsalis were also identified, and were assigned to 172 Operational Taxonomic Units belonging to six phyla (Andongma et al., Reference Andongma, Wan, Dong, Li, Desneux, White and Niu2015). All these reports provided useful information insights into biological control and the development of environmentally friendly pesticides (Wang et al., Reference Wang, Jin, Peng, Zhang, Chen and Hua2014). However, up to now no studies have compared the associated bacterial diversity between different field populations from different regions within the invasive range of B. dorsalis.
In this study, we characterized and compared the associated bacterial communities present in geographically diverse wild populations of B. dorsalis collected from different invasion areas of China through the construction of 16S rRNA gene libraries. Moreover, the putative function of associated bacteria of B. dorsalis was discussed. The knowledge of the microbial communities associated with wild populations of B. dorsalis from different invasion areas will help to understand the impact of the bacterial diversity on the ability of fruit fly to invade new areas.
Materials and methods
Fruit fly sampling
In this study, B. dorsalis samples from one laboratory colony and from 15 wild populations were analyzed. The laboratory population was provided by Guangdong Inspecting and Quarantine Technology Centre (IQTC) and reared using an artificial diet at 25°C with a relative humidity of 75%. The wild sampling sites included 14 sites in China and one site in Thailand. Samples were collected in each site at a single collecting time. The Chinese populations were sampled in eight representative south-eastern provinces to cover as much of the invasion area of B. dorsalis in China as possible. The only population sampled in Thailand was Sisaket, which is located in the northeast part of Thailand. Detailed information on sample collection is summarized in table 1 and fig. 1.
Fruit flies collection site. Geographical map of the Chinese provinces from where Bactrocera dorsalis samples were collected. The site numbers correspond with those shown in table 1.
Characteristics of 16 collection sites of Bactrocera dorsalis and bacterial 16S rRNA gene sequence clone libraries.
1 Total number of clones sequenced for each individual PCR libraries.
2 The OTU was determined based on a 3% distance level.
3 ID assigned for the recovered OTUs based on a 3% distance level.
4 Coverage was estimated as the ratio between the OTUs observed and the total number of OTUs estimated in the population (Chao1).
5 GenBank numbers of the submitted sequences based on a 3% distance level.
All of the wild samples were collected by trapping with methyl eugenol (ME). Approximately 3 h later insects were removed from traps and transported to the laboratory and identified by IQTC experts. As the ME is a male-specific attractant only males were collected. The samples of B. dorsalis were maintained at −80°C until DNA extraction.
Invasion regions
Flies were sampled from provinces in which they have had different time periods of invasive populations. According to these periods three different regions were considered: the early term invasion region (ET region), the medium term invasion region (MT region) and the late term invasion region (LT region). As shown in table 1, the ET region includes Yunnan, Guangdong, Guangxi and Guizhou provinces in which B. dorsalis has been found since 1970s (Bai & Song, Reference Bai and Song1997). The LT region includes the four sites in Shanghai city, where B. dorsalis had been established from 2001 (Zhou et al., Reference Zhou, Ye, Yuan and Pan2006), and the other collection sites were attributed to the MT region (Li et al., Reference Li, Wu, Chen, Wu and Li2012).
DNA extraction and 16S rRNA gene libraries construction
Before DNA extraction, the whole fruit flies were sterilized in 1% sodium hypochlorite (NaClO) and 97% ethanol for 5 min and then washed five times with sterile distilled water following the procedure described by Tang et al. (Reference Tang, Lv, Jing, Zhu and He2010).
The whole body of three adults from each B. dorsalis population were processed to construct the gene libraries. The bacterial DNA was extracted using the commercial Tissue/Cell DNA Mini Kit (Tiangen, China) following the manufacturer's protocol. The 16S rRNA gene was amplified using the universal primers 27F/1492R (Rani et al., Reference Rani, Sharma, Rajagopal, Adak and Bhatnagar2009). Each polymerase chain reaction (PCR) was conducted in a reaction volume of 50 µl containing 5 µl of 10 × reaction buffer (Mg2+), 1 µl of the extracted nucleic acid DNA (115.8 ng µl−1) as the template, 1U of Taq polymerase (Tiangen, China), 4 µl of the dNTP mixture (2.5 mM), and 1 µl of the forward and reverse primers (10 µM), respectively. The DNA amplification was performed on a Veriti TM 96-well Thermal Cycler (ABI, USA) with the following reaction conditions: 95°C for 3 min, 35 cycles of 94°C for 3 min, 55°C for 1 min, and 72°C for 1 min, and a final incubation at 72°C for 10 min. The PCR products were purified using a DNA Purification Kit (Tiangen, China). The purified PCR products were cloned into the PGMT-19 vector (TaKaRa) and transformed into Escherichia coli DH5α cells (TaKaRa). The transformation was verified using the RV-M/M13-47 universal primers. The amplification products were sequenced using the same primers.
Sequencing and phylogenetic analysis
For each gene library, at least 50 clones were chosen for sequencing. The sequence chromatograms were visually inspected, and the sequences were edited and aligned using MEGA 4.0 (Tamura et al., Reference Tamura, Dudley, Nei and Kumar2007). The chimeric sequences were removed using Mallard 1.02 (Ashelford et al., Reference Ashelford, Chuzhanova, Fry, Jones and Weightman2006), and the remaining sequences were blasted in The National Center for Biotechnology (NCBI; http://www.ncbi.nlm.nih.gov). The distance matrices were constructed using the DNADIST program of PHYLIP 3.69 (Felsenstein, Reference Felsenstein1989). The operational taxonomic units (OTUs), based on their identity and 3% distance level, were determined by DOTUR 1.53 (Schloss et al., Reference Schloss, Delalibera, Handelsman and Raffa2006). This program assigns sequences to OTUs based on sequence data by using values that are less than the cut-off level.
In order to estimate the adequacy of sampling, rarefaction curves were generated and the Chao1 richness estimator (Chao, Reference Chao1984) was calculated using the same software.
The Chao1 richness index estimates diversity of a community based on the number of singletons (OTUs represented by only one sequence) and doubletons (OTUs represented by two sequences) found in a sample (Bohannan & Hughes, Reference Bohannan and Hughes2003). The richness estimates are reported for 3% difference between sequences.
To study the phylogenetic relationships among the associated bacteria of B. dorsalis, we considered bacteria found in other fruit flies and the closely related bacterial sequences present in GenBank. Sequences of at least 1000 bp were considered for phylogenetic tree construction. Bacteria retrieved in C. capitata, B. oleae, Bactrocera zonata, four taxa of the subfamily Tephritinae (Campiglossa guttella, Dioxyna bidentis, Trupanea amoena, Tephritis matricariae and Noeeta genus) were selected as the reference sequences (Supplementary table 1). The phylogenetic relationships among these bacteria were estimated using the neighbour-joining method and the Jukes-Cantor model. For each 16S rRNA gene library, representative sequences for each OTU obtained were submitted to GenBank under the accession numbers JQ918001–JQ918075 (table 1).
Bacterial diversity comparison
We provided three kinds of bacterial diversity comparisons: male and female of the laboratory colony, among 15 wild populations, and among three different invasion regions. Clone percentage, clone frequencies across number of PCR libraries was used to describe the relative abundance of each genus. The bacterial diversity among three different invasion regions was statistically analyzed performing ANOVA analysis with Tukey-Kramer test using CoStat-Statistics Software v6.204 (http://www.cohort.com/costat.html).
Results
Bacterial 16S rRNA libraries and phylogenetic analysis
16S rRNA gene libraries were constructed for 15 wild populations and one laboratory population, and 789 sequences with a length of 1500 bp were obtained. All 789 sequences were grouped into 133 OTUs based on their identity (similarity = 100%) and 24 OTUs based on their 3% distance level (similarity >97%) (table 1). The rarefaction curves for most of the clone libraries tended to be saturated at 3% difference between sequences (Supplementary figure 1). Moreover, except for three populations (GZ, SJ and TH), the coverage was higher than 95% (table 1) suggesting that the number of clones sampled was enough to provide an adequate estimation of B. dorsalis bacterial diversity.
The 24 OTUs retrieved were affiliated to three classes, namely gamma -Proteobacteria (21 OTUs), delta -Proteobacteria (one OTU), and Bacilli (two OTUs), belonging to seven families (Xanthomonadaceae, Enterobacteriaceae, Pseudomonadaceae, Orbaceae, Desulfovibrionaceae, Streptococcaceae, and Enterococcaceae) and 13 genera (Stenotrophomonas, Pseudomonas, Enterobacter, Morganella, Kluyvera, Providencia, Klebsiella, Serratia, Citrobacter, Orbus, Erwinia, Lactococcus, and Vagococcus) (table 2). The results show that the Enterobacteriaceae family and Enterobacter genus were the most preponderant bacterial taxa of B. dorsalis with percentages of approximately 83 and 35% in the 789 clones, respectively (fig. 2).
Clone percentage of symbiotic bacteria of each genus. The percentage was calculated as the percentage of clones corresponding to the indicated genus of the total 789 clones.
Taxonomical assignment of the cloned 16S rRNA bacterial amplifications from Bactrocera dorsalis
1 The OTU was determined using the DOTUR program based on a 3% distance level.
2 Organism; GenBank Accession No. (BLAST using both the type and non-type strains in NCBI).
3 The clone percentages were calculated from a total of 789 clones in all of the 16 populations.
4 The bold values represented the clone percentage of the top three OTUs.
The 24 OTUs were distributed in seven families and all the seven families were well supported in seven clades, respectively (fig. 3). Sixteen OTUs were grouped in the Enterobacteriaceae clade together with bacteria retrieved in other fruit flies. However, within the Enterobacteriaceae, none of the OTUs found in B. dorsalis were present in the highly supported group, which was shared by the specific symbiotic bacteria of Tephritinae subfamily (Ca. Stammerula sp.) and B. oleae (Ca. Erwinia dacicola). Moreover, besides the clade of the family Orbaceae, that were represented by three OTUs (11, 14 and 16), the remaining family clades contained only one OTU (fig. 3).
Neighbour-joining (NJ) phylogenetic tree of bacteria associated with Bactrocera dorsalis based on the 16S rDNA sequences. The NJ tree was constructed using MEGA 4.0 with the NJ method and the Jukes-Cantor model and shows the phylogenetic relationship of symbiotic bacteria affiliated with 24 OTUs, the relative free-living bacteria and other bacteria found in other fruit flies. The grey box contains specific and unculturable symbiotic bacteria of Bactrocera oleae and Tephritinae subfamily. The scale bar corresponds to 0.02 substitutions per nucleotide position. The percentage bootstrap values above 50 (1000) replicates are indicated as notes.
Bacterial diversity associated with the laboratory colony
In the laboratory population, the bacteria found in female fruit flies belonged to two genera (Stenotrophomonas and Enterobacter), whereas the bacteria found in males belonged to three genera (Stenotrophomonas, Enterobacter and Klebsiella), of which Enterobacter was the most dominant genus (with a clone percentage of 64 and 88%, respectively) (fig. 4).
Relative abundance of the different genera in the associated bacteria from one laboratory and 15 wild populations.
Bacterial diversity associated with different wild populations
Among the total 24 OTUs based on 3% distance level, the Guangdong-Guangzhou (GZ) population had the largest number of OTUs (n = 8), whereas the Fujian-Xiamen (XM) and Shanghai-Qingpu (QP), had the least number of OTUs (n = 2) (table 1).
The genus-level bacterial diversity was analyzed for each population, and the detailed results are shown in fig. 4. An average of 3.9 genera of bacteria per population was present in the wild populations of B. dorsalis. Of the 15 wild populations, Guangdong-GZ was the population with the highest number of genera. Seven different genera were found, including Stenotrophomonas, Orbus, Citrobacter, Morganella, Providencia, Serratia, and the predominant Enterobacter (with a clone percentage of 69%) (fig. 4).
The populations of Shanghai Region (SJ, QP, BS, and JD) showed the fewest number of genera. Particularly in the Shanghai-JD population, the bacteria were affiliated with two genera, namely Klebsiella (25% clone percentage) and Pseudomonas (75% clone percentage), whereas three genera were present in the SJ, QP, and BS populations. Enterobacter was the most predominant genus in SJ (80% clone percentage), whereas Providencia and Citrobacter were the most predominant genera in QP (50% clone percentage) and BS (53% clone percentage), respectively.
We consider the genera that were recovered in more than five Chinese populations of B. dorsalis as the most common genera. Thus, the most frequent genus was Klebsiella, which was found in 13 populations, followed by Citrobacter (ten populations), Enterobacter (nine populations), Orbus (eight populations) and Providencia (5 populations).
Bacterial diversity in different invasion regions
As is shown in table 3, the bacteria obtained from the flies collected in the ET region were assigned to three classes and grouped into 10 genera: Stenotrophomonas, Enterobacter, Klebsiella, Orbus, Citrobacter, Providencia, Morganella, Unknown Delta-proteobacteria, Lactococcus, and Vagococcus. The bacteria of the flies collected in the MT region were attributed to two classes and eight genera. Interestingly, the bacteria found in B. dorsalis populations collected from the LT region were attributed only to the Gamma-Proteobacteria class and classified into six genera, four of which are affiliated with the family Enterobacteriaceae (table 3). Four genera, including Enterobacter, Klebsiella, Citrobacter and Providencia were found in all three invasion regions.
Bacteria found in different invasion regions of Bactrocera dorsalis
“+” indicates that the genus was detected in this region.
The result associated with three invasion regions showed that, considering OTUs based on a 3% distance level, the bacterial diversity in the ET region was statistically different from that in the LT regions (ANOVA, F = 5.15, df = 2,11, P = 0.036). However, there was no significant difference between the ET and MT regions or between the MT and LT regions (fig. 5).
Bacterial diversity (%) in the three invasion regions. ANOVA followed by a Tukey post-hoc test. Letters above the bars indicate significant differences at P < 0.05 (bars are means ± SE). Numbers inside bars represent sample sizes.
Discussion
Bacterial diversity of B. dorsalis
The 16S rRNA gene libraries from the single laboratory and 15 wild populations indicate that the associated bacteria of B. dorsalis are affiliated with the classes Gamma-proteobacteria, Delta-proteobacteria, and Bacilli, of which Gamma-proteobacteria is the most predominant class with a percentage of 90.68%. Gamma-proteobacteria include taxa of important symbiotic bacteria of many insects, such as other fruit flies, mealybugs, mosquitoes, and beetles (Thao et al., Reference Thao, Gullan and Baumann2002; Capuzzo et al., Reference Capuzzo, Firrao, Mazzon, Squartini and Girolami2005; Schloss et al., Reference Schloss, Delalibera, Handelsman and Raffa2006; Mazzon et al., Reference Mazzon, Piscedda, Simonato, Martinez-Sañudo, Squartini and Girolami2008; Rani et al., Reference Rani, Sharma, Rajagopal, Adak and Bhatnagar2009). Within the Gamma-proteobacteria, our results also confirmed that Enterobacteriaceae (in particular with the genera Enterobacter, Klebsiella, Citrobacter, and Providencia), are the predominant bacteria family as previously found in B. dorsalis (Jang & Nishijima, Reference Jang and Nishijima1990) and in other tephritid flies (i.e. Behar et al., Reference Behar, Yuval and Jurkevitch2008b ; Wang et al., Reference Wang, Jin, Peng, Zhang, Chen and Hua2014) . Similarly, in a single field-collected population of B. dorsalis, Wang et al. (Reference Wang, Jin and Zhang2011) reported a predominance of Gamma-proteobacteria followed by Firmicutes and few taxa of Delta-proteobacteria, Actinobacteria, Bacteroidetes, Flavobacteria, and Alpha-proteobacteria. Moreover, the analyses showed Klebsiella, Citrobacter, Enterobacter, Pectobacterium and Serratia the common bacterial phylotypes for all three libraries. On the contrary, Andongma et al. (Reference Andongma, Wan, Dong, Li, Desneux, White and Niu2015) compared the microbiota across the life stage of B. dorsalis and came to the conclusion that Enterococcaceae (Bacilli) was the most abundant family in the adult flies.
While the cloning method used in this study is optimal for identifying bacteria at the genus level it is limited for the quantitative analysis because of eventual biases of clone frequencies occurring at several levels during the process (e.g. primer bias, amplification efficiency bias, ligation bias etc.). Future studies considering next generation sequencing techniques will help to overcome these limitations.
No bacteria belonging to the genus Wolbachia were detected in any of the populations used for this study, even though our sample collecting locations cover almost all the distribution provinces of B. dorsalis in China. Similarly, another study reported the absence of Wolbachia in Chinese populations of B. dorsalis (Yunnan, Fujian and Wuhan population) (Augustinos et al., Reference Augustinos, Drosopoulou, Gariou-Papalexiou, Asimakis, Cáceres, Tsiamis, Bourtzis, Mavragani-Tsipidou and Zacharopoulou2015). In conflict with these results, Sun et al. (Reference Sun, Cui and Li2007), detected Wolbachia in four of five Chinese populations but with low infection rates (0.7–3%). As widely reported and discussed by Morrow et al. (Reference Morrow, Frommer, Shearman and Riegler2015a , Reference Morrow, Frommer, Royer, Shearman and Riegler b ) for Australian fruit flies, infections of Wolbachia at low frequency may be indicative of transient infections that result from spillover events from a yet unknown source, questioning their vertical inheritance.
In our study, we generally found more bacterial diversity, considering the OTUs based on a 3% distance level, in wild populations than in the laboratory population (table 1), and the bacterial diversity in ET regions was statistically higher than that in LT regions. These differences could be influenced by environmental factors such as the habitat and the type of diet (i.e. host plant) that these fruit flies are exposed to as proposed by several authors. Wang et al. (Reference Wang, Jin and Zhang2011) reported that the gut bacteria diversity of field collected B. dorsalis was higher than that of fruit flies reared on sterile artificial food. They came to a conclusion that different environmental conditions and food supply could influence the diversity of the harboured bacterial communities and increase community variations. Morrow et al. (Reference Morrow, Frommer, Shearman and Riegler2015a , Reference Morrow, Frommer, Royer, Shearman and Riegler b ) detected and compared the microbiome associated with field-caught and laboratory-adapted Australian Tephritid fruit fly species, and also came to the conclusion that the environment, species identity and ecology can influence the microbial composition.
Possible function of associated bacteria
The functions of fruit flies’ symbionts and free-living associated bacteria have been reported via four aspects, including participating in the carbon and nitrogen cycle of fruit flies to help them synthesize essential amino acids and minerals; attracting the fly to feed on host plant; degrading the toxic compounds from host plants and preventing the colonization of pathogenic bacteria (Robacker & Lauzon, Reference Robacker and Lauzon2002; Lauzon et al., Reference Lauzon, Bussert, Sjogren and Prokopy2003; Behar et al., Reference Behar, Yuval and Jurkevitch2005, Reference Behar, Jurkevitch and Yuval2008a , Reference Behar, Yuval and Jurkevitch b , Reference Behar, Yuval and Jurkevitch c ; Robacker et al., Reference Robacker, Lauzon, Patt, Margara and Sacchetti2009). Many bacteria retrieved in this study have been reported showing beneficial function to host flies such as, Lactococcus that has been reported to inhibit the growth of pathogenic bacteria in Anastrepha ludens (Kuzina et al., Reference Kuzina, Peloquin, Vacek and Miler2001). Moreover, it has been reported that some of the bacteria retrieved in this study (Klebsiella, Citrobacter, Providencia, and most Enterobacter) have been reported to play a key role in the biological fixation of nitrogen (Robacker & Lauzon, Reference Robacker and Lauzon2002; Lauzon et al., Reference Lauzon, Bussert, Sjogren and Prokopy2003). Almost all the bacteria we found in wild male populations belonged to the genus Enterobacter, which can provide advantages to host flies through contributing essential nutrients during the adult life stage (Ben-Yosef et al., Reference Ben-Yosef, Jurkevitch and Yuval2008). E. cloacae and C. freundii have been reported for their high attractive ability to host flies (Wang et al., Reference Wang, Jin, Peng, Zhang, Chen and Hua2014). The abundant bacteria retrieved in our study could be useful to find more bacteria with attracting ability, which may be important for the biocontrol management of B. dorsalis.
One previous report has suggested an important role of the gut microbiota of insects in the host range and host suitability. The Asian longhorn beetle (Anoplophora glabripennis), an invasive pest with a broad host range, presents a wide bacterial diversity, whereas the linden borer (Saperda vestita), a cerambycid with a more restricted host range, contains only a small subset of the same bacteria (Schloss et al., Reference Schloss, Delalibera, Handelsman and Raffa2006). Recently, collaborative invasion (co-invasion) between the invasive species and its niche related species has also been a hot topic of research, and the most typical examples include the co-invasion of white fly Bemisia tabaci (Gennadius) biotype with the gemini viruses and the co-invasion of red turpentine beetle Dendroctonus valens (LeConte) with its associated fungus Leptographium procerum, (Liu et al., Reference Liu, Zhao, Jiang, Zhou and Liu2009; Lu et al., Reference Lu, Zhou, De Beer, Wingfield and Sun2009; Himler et al., Reference Himler, Adachi-Hagimori, Bergen, Kozuch, Kelly and Tabashnik2011). Given the significance of microbial symbionts in the invasion of insects and the important role of associated bacteria (Feldhaar, Reference Feldhaar2011), the research on the microbiota may help us to uncover the invasion mechanism of invasive insect pests such as B. dorsalis. However, the speculated co-invasion between associated bacteria and host flies needs further investigation, including verifying the function of the microbiota through exposing flies to different bacteria. Meanwhile, the bacterial community associated with females from different wild populations and different invasion regions should also be investigated in the future, because the bacterial harboured by females can be transmitted to the larvae and may thus contribute to the invasive potential of the fly since they may determine the ability of larvae to utilize new hosts. Moreover, fruit flies from more instars, more distribution areas, and more genera should also be used in future studies to uncover the co-invasion mechanism of fruit flies and their associated bacteria.
Supplementary Material
The supplementary material for this article can be found at http://dx.doi.org/10.1017/S0007485316000390.
Acknowledgements
This work was supported by the National 973 Project (No. 2009CB119204) and Basic Research Program (No. 2015KC001). The authors would like to thank the experts in Guangdong Inspecting and Quarantine Technology Centre (IQTC) for species identification of oriental fruit fly and Guangchao Wu in Forest Pest Management and Quarantine Station of Shanghai for providing fruit fly samples.
