- Research article
- Open Access
The gut microbiota of larvae of Rhynchophorus ferrugineus Oliver (Coleoptera: Curculionidae)
BMC Microbiology volume 14, Article number: 136 (2014)
The red palm weevil (RPW) Rhynchophorus ferrugineus Olivier (Coleoptera: Curculionidae) is one of the major pests of palms. The larvae bore into the palm trunk and feed on the palm tender tissues and sap, leading the host tree to death. The gut microbiota of insects plays a remarkable role in the host life and understanding the relationship dynamics between insects and their microbiota may improve the biological control of insect pests. The purpose of this study was to analyse the diversity of the gut microbiota of field-caught RPW larvae sampled in Sicily (Italy).
The 16S rRNA gene-based Temporal Thermal Gradient Gel Electrophoresis (TTGE) of the gut microbiota of RPW field-trapped larvae revealed low bacterial diversity and stability of the community over seasons and among pools of larvae from different host trees. Pyrosequencing of the 16S rRNA gene V3 region confirmed low complexity and assigned 98% of the 75,564 reads to only three phyla: Proteobacteria (64.7%) Bacteroidetes (23.6%) and Firmicutes (9.6%) and three main families [Enterobacteriaceae (61.5%), Porphyromonadaceae (22.1%) and Streptococcaceae (8.9%)]. More than half of the reads could be classified at the genus level and eight bacterial genera were detected in the larval RPW gut at an abundance ≥1%: Dysgonomonas (21.8%), Lactococcus (8.9%), Salmonella (6.8%), Enterobacter (3.8%), Budvicia (2.8%), Entomoplasma (1.4%), Bacteroides (1.3%) and Comamonas (1%). High abundance of Enterobacteriaceae was also detected by culturing under aerobic conditions. Unexpectedly, acetic acid bacteria (AAB), that are known to establish symbiotic associations with insects relying on sugar-based diets, were not detected.
The RPW gut microbiota is composed mainly of facultative and obligate anaerobic bacteria with a fermentative metabolism. These bacteria are supposedly responsible for palm tissue fermentation in the tunnels where RPW larvae thrive and might have a key role in the insect nutrition, and other functions that need to be investigated.
The red palm weevil (RPW) Rhynchophorus ferrugineus Olivier (Coleoptera: Curculionidae) is widely considered the most damaging insect pest of palms in the world, even in all the countries where it has been accidentally introduced . RPW larvae feed within the apical growing point of the palms, producing a wet fermenting frass inside the tunnels , creating extensive damage to palm tissues and weakening the structure of the palm trunk; the resulting damage is often only visible long after infestation, when palms are close to death [3–5] (Additional file 1).
Insect intestinal tracts harbour rich communities of non-pathogenic microorganisms [6, 7] and a single gut can harbour 105–109 prokaryotic cells  that have been affiliated to twenty-six phyla, at least for the insects studied to date . It is increasingly evident that the microbiota of animals (humans included) plays a remarkable role in the host life. The genetic wealth of the microbiota affects all aspects of the holobiont’s (host plus all of its associated microorganisms) fitness such as adaptation, survival, development, growth, reproduction and evolution . When not strictly essential for survival, the insect gut microbiota affects many aspects of host phenotype; it can increase the digestive efficiency of soluble plant polysaccharides [10, 11] and can mediate interactions between the host and potential pathogens . Recent work suggests that the gut microbiota not only provide nutrients, but is also involved in the development and maintenance of the host immune system. However, the complexity, dynamics and types of interactions between the insect hosts and their gut microbiota are far from being well understood . Understanding the relationship dynamics between insects and their microbiota can improve the biocontrol of insect pests, which is a focus of much insect gut microbiology studies.
Despite the economic and environmental damages caused by the RPW in all the areas where it is endemic and where it has been accidentally introduced, little is known about its gut microbiota. The bacterial community that is embedded in the frass produced inside the tunnels of the palm Phoenix canariensis Chabaud by the RPW larvae is dominated by Enterobacteriaceae with a facultative fermentative metabolism .
The purpose of this study was to analyse the diversity of the gut microbiota of the R. ferrugineus larvae, that represent the development stage responsible for damages to palms. Field-caught larvae were sampled from its favourite host P. canariensis in different seasons and sites in Sicily (Italy), and analysed for the diversity of their gut microbiota. The analysis of the bacterial community was carried out by culture-independent methods using temporal thermal gradient gel electrophoresis (TTGE) and FLX454 pyrosequencing of PCR-generated amplicons from the 16S rRNA gene.
Total diversity of the gut microbiota of field caught RPW larvae
Bacterial TTGE profiles were generated using PCR-amplified bacterial 16S rRNA gene fragments from the content of pooled RPW larval guts collected from the trunks of infested P. canariensis palms in three different seasons and two areas in Sicily (Italy). TTGE band profiles indicate the presence of an average of 25 bands per sample, that correspond to putative bacterial phylotypes in RPW larval guts. An example of TTGE gel is shown in Figure 1, where three different pooled guts collected in December 2010 and April 2011 in Palermo (lanes 1 and 2, respectively), and in April 2011 in San Vito lo Capo (Trapani, lane 3) were analysed. All samples shared 16 bands, while 4, 2 and 4 bands were unique for samples 1, 2, 3, respectively. Similar profiles were obtained from larvae collected in October both in Palermo and Trapani (data not shown). Random sequencing of TTGE bands identified the presence of uncultured Gammaproteobacteria (of the genera Pantoea and Enterobacter) and Firmicutes (of genera Megasphaera and Clostridium) (Figure 1).
Pyrosequencing of the 16S rRNA gene amplicons was carried out on three pooled RPW larval guts sampled in Palermo in April 2011 (indicated as lot A). The analysis produced a total of 79,204 reads with an average length of 320.6 nucleotides that became, after quality filtering and clustering (needed for Ribosomal Database Project analysis), 75,564 for 97%, 76,724 for 95%, and 73,579 for 90% of similarity (Additional file 2). Reads were assigned to 41 operational taxonomic units (OTUs) at 90% of sequence identity threshold, and to 45 OTUs at 95% and 97% identity threshold, respectively, in order to perform rarefaction analysis. The total number of clusters obtained after filtering was of 2,107 (1,756 singletons) for 97%, 910 (530 singletons) for 95%, and 244 (124 singletons) for 90% of similarity, respectively. The rarefaction curves tended towards saturation at similar numbers of clusters at 97%, 95% and 90% pairwise ID thresholds (Figure 2). Subsequent analysis was, therefore, conducted at 97% ID.
Only three phyla account for 98% of the reads: these are Proteobacteria (64.7%), Bacteroidetes (23.6%) and Firmicutes (9.6%); the remaining 2% is represented by Tenericutes (1.4%) Fusobacteria (0.4%) and other Bacteria (0.2%) (Figure 3a). Proteobacteria are mainly represented by Gammaproteobacteria (96.7%) followed by Betaproteobacteria (2.71%) (Figure 3b). More than 98% of the reads were classified at the family level, with Enterobacteriaceae representing the 61.5% of the assemblage, followed by Porphyromonadaceae (22.1%) and Streptococcaceae (8.9%) (Additional file 3). More than half of the reads (52.7%) could be classified at the genus level and eight bacterial genera were detected in the larval RPW gut at an abundance ≥1% (Figure 4a). Dysgonomonas sequences account for the 21.8% of the whole sequences and this is the most represented genus in the gut of RPW larvae, followed by Lactococcus (8.9%) Salmonella (6.8%), Enterobacter (3.8%), Budvicia (2.8%), Entomoplasma (1.4%) Bacteroides (1.3%) and Comamonas (1%). Other twelve genera are represented at a value between 1% and 0.1% (Figure 4b). The phylogenetic tree of 16S rRNA gene amplicons clustered at 97% consensus is shown in the Additional file 4.
Diversity of cultivable bacteria
Bacterial isolation under aerobic conditions was carried out on three lots of three pooled RPW larval guts (lots A, B, C), all sampled in April 2011. The dilution plate counts on NA gave an average of 1.5 × 107 CFU gut-1, without differences among the three pools. Forty-four phenotypically different isolates, obtained from the dilution plates of the three pooled guts independently, were grouped into ten OTUs (Table 1) by ARDRA analysis with the restriction enzymes AfaI and AluI, that gave five and seven different restriction profiles, respectively (data not shown). 16S rRNA gene sequencing of representative isolates assigned the cultivable bacteria to the families Enterobacteriaceae (68.2%), Bacillaceae (20.5%), Comamonadaceae (9%) and Xanthomonadaceae (2.7%) (Table 1). The genus Citrobacter is the most abundant among the isolates (29.55%), followed by the genera Klebsiella (20.45%), Bacillus (20.45%) and Budvicia (11.36%).
Most of the sequences having homology with those of RPW isolates are from bacteria isolated from animals’ gut or from plants (endophytes), as well as from wastewater or bioremediation treatment plants and anaerobic marine sediments. Some of the Citrobacter and Klebsiella 16S rRNA sequences are almost identical to those from bacteria previously isolated from the frass produced by RPW larvae in the tunnels of palm trees (Additional file 5) .
Several attempts were made to surface-sterilize the larvae using different protocols; nevertheless the control plates, obtained by streaking on Nutrient Agar the cuticle of sterilized larvae, showed the growth of some colonies. Seven of these colonies were purified and analysed by ARDRA as described above. One representative isolate for each of the three OTUs obtained was chosen for 16S rRNA gene sequencing and all the three isolates were affiliated to the genus Bacillus, with their closest relatives in B. thuringiensis  and B. aerophilus  (Additional file 5).
Enrichment cultures were set for the isolation of acetic acid bacteria (AAB). AAB are known to establish symbiotic associations with the midgut of insects relying on sugar-based diets, such as nectars, fruit sugars, or phloem sap . At the end of the incubation period, four CaCO3 dissolving colonies were isolated from the enrichment cultures and identified by 16S rDNA sequencing. Unexpectedly, all the isolates that were able to use sorbitol and to dissolve CaCO3 in the agar plates were assigned to the genus Klebsiella (Additional file 5).
In this study, the diversity of the gut microbiota of Rhynchophorus ferrugineus (RPW), collected on infested palm trees Phoenix canariensis, was first analysed by TTGE of the PCR-amplified bacterial 16S rRNA gene fragments. The TTGE profiles obtained from different lots of larvae, sampled in different seasons and geographical sites, show relatively low complexity (average of 25 OTUs) and high similarities regardless the site of sampling and season, suggesting that the composition of the RPW microbiota is stable over time and among pools of larvae from different host trees.
In order to identify the gut bacterial community of RPW larvae, the variable region 2 (V2) of the bacterial 16S rRNA gene, already successfully employed in the analysis of several microbial communities [17–19], was analysed by pyrosequencing. The analysis confirmed that the bacterial community of the RPW larvae has low diversity although, as expected, more OTUs were identified in respect to TTGE analysis. Contrasting results are reported for bacterial diversity of gut microbiota of other coleopterans; high diversity and complexity was observed among tree xylophagous beetles that rely on the microbiota for efficient lignocellulose metabolism and thus survival , while low diversity was recorded in the gut of the red turpentine beetle . The RPW larvae are the major responsible for the palm damages because they live throughout their development inside the palm stem, feeding exclusively on palm tissues. This peculiar lifestyle may account for the low diversity detected in the gut of field sampled larvae of R. ferrugineus, regardless the investigation methods. There is strong evidence that mainly taxonomy and diet of the host can affect an organism’s gut microbial community [8, 21]. RPW larvae feed on nutrient-poor palm tissues and sap that contain mainly sucrose and glucose  but are poor of nitrogen [20, 23, 24]; an excess of sugars is known to reduce the complexity of the gut microbiota [25, 26]. Conversely, complex substrates, such as lignocellulose-derived materials, select complex gut bacterial communities even in highly divergent insect groups . Feeding on live trees may also expose the insect to tree physiological responses, that may play a selection pressure on its microbiota. Moreover, the antimicrobial activity of larvae  may further shape the gut community of RPW.
The gut of RPW larvae is dominated by three phyla, Proteobacteria, Bacteroidetes and Firmicutes, that account for 98% of the assemblage. These same phyla were also found in the sugarcane weevil Sphenophorus levis Vaurie , which belong to the Dryophthorinae subfamily as R. ferrugineus, and that is the only weevil, to our best knowledge, that has been characterized in its microbiota. Proteobacteria and Firmicutes represent also the predominant bacterial phyla in bark beetles  and, in general, in all insect guts studied so far, while Bacteroidetes are more prevalent in termites, detritivorous insects and, among Coleoptera, in the root feeding Melolontha melolontha L. (Coleoptera: Melolonthidae) .
The genus Dysgonomonas is, unexpectedly, the most represented in the gut of RPW larvae. Dysgonomonas (phylum Bacteroidetes) are facultative anaerobes with a fermentative metabolism producing acids and no gas, that were first recovered from a human infected gall bladder . Dysgonomonas is described as an opportunistic human pathogen but its habitat is unknown. Members of this genus were recently detected in microbial fuel cells (MFC) anode biofilms , in the gut of house flies (Musca domestica L.)  and in eight separate Drosophila populations where its presence is not restricted to any one locality, species, or diet type . Its presence in such a high number in the insect gut, and in RPW gut in particular, deserves to be further investigated because it might play an important role in the insect biology.
Salmonella, Enterobacter, Budvicia and other Enterobacteriaceae are highly represented in the 454 assemblage; as in other insects, they could play a beneficial role in nutrition, in the degradation of plant polymers and fermentation of sap sugars. Members of Enterobacteriaceae were also identified as intracellular symbionts of grain weevils Sitophilus spp. (Curculionidae)  and some isolates are able to fix nitrogen, thus contributing to a supplementary nitrogen source [20, 23].
Entomoplasma is the sixth genus to be represented in terms of abundance in the RPW gut (3%). Entomoplasma is a glucose fermenting non-helical mollicutes and its presence in the RPW gut is consistent with what is presently known of its habitat. This genus could be considered a marker of the Coleopteran microbiota. All five currently described Entomoplasma species, in fact, were isolated from the gut or haemolymph of various firefly beetles (Coleoptera: Lampyridae) and green tiger beetles (Coleoptera: Cicindelidae) .
In spite of being affiliated to three different phyla, all the first six dominating bacterial genera of the RPW gut are facultative or obligate anaerobes with a fermentative metabolism. Dysgonomonas is a facultative anaerobe that produces acids and no gas; Lactococcus is a omofermenting lactic acid bacterium and some species can produce acetoin. The other three dominating genera belong to the Enterobacteriaceae characterized by mixed acid fermentation with production of lactic, acetic, succinic acid and ethanol (Salmonella), or 2,3-butanediol fermentation, producing butanediol, ethanol, CO2 and H2 (Enterobacter and Budvicia). Entomoplasma is also a glucose fermenting bacterium. These results suggest that the peculiar life-style of RPW larva and its gut exert a strong selective pressure towards those microbial species that are specialised to grow in a high sugar environment and that these species probably have a competitive advantage on those that cannot tolerate organic acids.
Interestingly, two genera of Enterobacteriaceae, Pantoea and Rahnella, which had previously been isolated from frass, were not detected in the gut. Rahnella isolates from frass have their closest relatives in components of the microbiota of the red turpentine beetle Dendroctonus valens LeConte (Coleoptera: Scolytidae)  and of the larvae of the lepidopteran Hepialus gonggaensis Fu & Huang (Lepidoptera: Hepialidae) ; Pantoea from frass are close to bacteria of the fungus garden of the leaf-cutter ant Atta colombica Guérin-Méneville (Hymenoptera: Formicidae), where they contribute to external plant biomass degradation and nitrogen fixation  (Additional file 5). High identities of RPW gut isolates with frass isolates and with other beneficial insect-associated bacteria suggest that the RPW gut microbiota cooperates, in a continuum with the frass microbiota, to the fitness of the larva inside the palm. Thus, while a unique midgut-associated microbiota can be distinguished from the environmental bacterial community in some insects , the peculiar lifestyle of RPW larvae makes such discrimination difficult or probably meaningless. In fact, RPW larvae feed in a very confined environment, consisting of tunnels burrowed in the palm trunk, where they continuously ingest both fresh palm tissues and frass, composed of chewed and/or digested plant tissue, so that re-acquisition by ingestion of bacteria from the environment is highly probable to occur.
Beyond nutritional aspects, the gut and frass fermentation products, such as acetoin and organic acid derivatives, ethyl esters, act as insect aggregation pheromones playing a role of attraction to other insects and promoting new oviposition events on the same tree .
Acidification caused by bacterial fermentation could also confer other advantages to the insect host, as some microbial toxins of Lepidoptera, such as Bacillus thuringiensis toxins, are activated by alkaline conditions. Thus, the RPW microbiota might help protect this insect from B. thuringiensis toxin by decreasing the midgut pH . Moreover, together with that of fermenting yeasts, the bacterial metabolic activity increases the temperature inside the palm tissues, helping weevil overwintering .
It is well known that culture-based methods fail in describing the total diversity of natural bacterial communities so that, as expected, culture-independent techniques yielded a higher diversity of bacteria. Culturing under aerobic conditions led to the detection of nine bacterial genera in the RPW larval gut. Both pyrosequencing and culturing revealed that Enterobacteriaceae is the most represented bacterial family in the gut of RPW larvae.
In this work, the culture-based approach helped in obtaining a better description of some members of Enterobacteriaceae as the complete sequence of the 16S rRNA gene could be obtained from the isolated bacteria. The pyrosequencing approach, that relies upon a short 16S rRNA gene fragment, did not detect sequences of the genus Klebsiella, that was instead abundantly isolated by culturing. Failing of its detection could be due to low variability of the V2 region between Klebsiella and Enterobacter[12, 40] and the sequences of Klebsiella might have been included in the genus Enterobacter by the RDP Classifier software.
Another genus detected by cultivation but absent in the 454 assemblage was Bacillus that might be present at very low levels in the RPW gut, so that its detection might be impaired by PCR biases. Bacilli isolated from the gut are close to B. muralis and B. simplex, and cluster separately from palm endophyte bacilli and frass bacilli previously isolated, that are related to the B. cereus/thuringiensis group.
Cuticle Bacillus isolates, that survived sterilization procedures, form a separate cluster from gut bacilli and are closer to the Bacillus isolates previously obtained from frass and from healthy palms as endophytes  (Additional file 5). This suggest that they belong to a bacterial community external to the larvae, that might contribute to the fitness of larvae inside the plant tissues. The cuticle aerobic spore-forming bacteria might produce antimicrobial molecules that could negatively affect the sensitivity of the larvae to entomopathogenic fungi and bacteria .
A low bacterial diversity and the presence of a prevailing sugar-fermenting microbiota suggest that the digestion of plant polymers (cellulose, hemicellulose) is not a primary function of the RPW larvae. However, cellulolytic and hemicellulolytic bacteria were previously isolated by enrichment cultures from the gut of RPW larvae and were mainly affiliated to the Gamma and Alphaproteobacteria of the genera Pseudomonas, Enterobacter Microbacterium and Paenibacillus . The presence of these genera in the RPW gut was confirmed by pyrosequencing (Additional file 6). Matching the 454-reads with the 16S rRNA gene sequences of the gut cellulolytic isolates, we obtained up to 99% identity of cluster_3902 (3 sequences) with the cellulolytic isolate Pseudomonas sp. R-8 (Genbank accession JN167546) and 98% identity of five different clusters (for a total of 159 sequences) with the cellulolytic RPW gut isolate Enterobacter sp. R-10 (Genbank accession JN167548) (Additional file 6). Both isolates were positive to the Congo Red test and able to grow on xylan in pure culture  but their hydrolytic activity on plant polymers in situ has to be demonstrated (as, for example, it might be inhibited by sap sugars).
The gut of insects that rely on sugar-based diets, particularly those belonging to the orders Diptera, Hymenoptera and Hemiptera, are often dominated by acetic acid bacteria (AAB), . Although the larval RPW diet is almost exclusively based on sugars, we were unable to detect AAB using a consolidated method based on the enrichment culture technique . Moreover, the absence of AAB in the RPW gut was confirmed by deep sequencing, where only two sequences were affiliated to the genus Acidisoma (Acetobacteriaceae) (Additional file 2). AAB are common in sugary acidic and alcoholic habitats, but are usually limited by nutrients other that their primary carbon source. AAB are common in fruit-feeding Drosophila species but are absent in flower-feeding flies . Their absence in the RPW larvae could be explained by microbial interactions occurring inside the gut. The enrichment cultures set to specifically isolate AAB led, instead, to the isolation of Klebsiella strains that could outcompete AABs and that could fulfil also the nitrogen fixation function [20, 43], allowing the insect to live on a substrate with a high C/N ratio.
The RPW microbiota is composed mainly of facultative and obligate anaerobic bacteria with a fermentative metabolism. These bacteria might have a key role in the insect nutrition, and other functions that need to be investigated. Further research, focusing on the functional traits of the bacteria inhabiting the gut of R. ferrugineus, is critically important to establish if some bacteria may exert an essential role for the insect or might represent an obstacle for the optimization and promotion of the use of entomopathogenic fungi and bacilli in an integrated pest management approach.
Sampling of RPW larvae and gut extraction
Field caught RPW late instar larvae (hereafter called larvae) were collected in Winter and Spring from infested palms of the species Phoenix canariensis Chabaud, located in the urban and peri-urban area of Palermo, and in San Vito Lo Capo (Trapani), (Italy) (Additional file 1). The palms were cut down following phytosanitary measures for the control and eradication of R. ferrugineus (Regional Decree 6 March 2007). The palms were not treated by chemical or biological pesticides.
The temperature was measured in 6 healthy and 6 infested palm trees during sampling at April 2011. Temperature was measured using a Bi-metal control digital thermometer (Wika - 360A005A4HS) by burrowing a small hole in the trunks, where the probe was inserted inside the palm trees. The average temperature of infested palm trees was 32.13°C ± 0.83, while the average temperature calculated at the same time for healthy palm trees was 25.95°C ± 0.71 and the ambient atmosphere was 26.2°C.
All sampled larvae were maintained in a plastic box with their own frass, taken from tunnels, and immediately transported to the laboratory for analysis. Each specimen was weighed, placed at -80°C for 30 min and surface sterilized with sodium hypochlorite and ethanol as described elsewhere [2, 44]. Late-instar larvae (average weight = 3.5 g ± 0.7 g, body length 3 cm ± 0.6 head-capsule 6.0 mm ± 0.8), corresponding in general to the 7th instar, were used. Larvae sterilization control was performed by streaking each intact larva on the surface of a Nutrient Agar (NA, Difco) plate. Larvae were dissected, the whole gut was aseptically removed and used for DNA extraction and bacterial isolation. Each sample consisted of the content of three pooled guts extracted from three larvae of the same weight and caught at the same time in the same palm tree.
Total bacterial diversity was assessed by Temporal Thermal Gradient gel Electrophoresis (TTGE) of 16S rDNA PCR products. DNA extraction form guts was carried out using the QIAamp DNA Stool Mini Kit, QIAGEN® (Qiagen, Hilden, Germany) according to the manufacture’s protocol and performing a lysis step at 95°C in order to obtain better lysis of Gram positive bacteria. A DNA region of approximately 200 base pairs was PCR-amplified from total DNAs. PCR was carried out using universal eubacterial oligonucleotide primers 341f-GC (5′-CGCCCGCCGCGCGCGGCGGGCGGGGCGGGGGCACGGGGGGCCTACGGGAGGCAGCAG-3′) and 534r (5′-ATTACCGCGGCTGCTGG-3′) targeting the variable V3 region of the 16S rRNA gene .
PCR were carried out using Phire Hot Start II DNA Polymerase (Thermo Scientific), 1X PCR buffer, 500 nM each primer, 0.20 mM dNTP and. 100 ng of DNA in a final volume of 25 μl. Cycling conditions were: 98°C for 30 sec, followed by 35 cycles of 98°C for 10 sec, 58°C for 10 sec and 72°C for 15 sec, followed by a final extension at 72°C for 2 min.
PCR products were fractionated on polyacrylamide gel (polyacrylamide:bis 29:1) 8%, Urea 7 M, Formamide 10% v/v, TAE 1.5X, at 70 V for 21 h in DCode (Bio-Rad) apparatus with a starting temperature of 57°C and a temperature ramp rate of 0.4°C h-1. Gels were stained with SYBRGold nucleic acid gel stain (Molecular Probes, Invitrogen) for 30 min and viewed under UV light.
Random bands were excised with a sterile scalpel immediately after visualisation, rinsed in 100 μl of distilled water and incubated in 30–50 μl of water, depending on band intensity, to elute DNA. DNA was re-amplified using the PCR-DGGE primers and products checked by agarose gel electrophoresis. The PCR products were purified using the QIAGEN PCR purification kit (Qiagen Hilden, Germany) and sequenced using the 534r primer. Partial bacterial 16S rRNA gene sequences (approximately 160 bp) were subjected to a NCBI nucleotide BLAST search (http://blast.ncbi.nlm.nih.gov/Blast.cgi) to identify sequences of the highest similarity. TTGE band sequences longer than 150 bp were deposited in Genbank under accession number KC763481 to KC763483.
The variable region 2 (V2) of the bacterial 16S rRNA gene was amplified with the primers 27 F (5′-AGAGTTTGATCMTGGCTCAG-3′) and 338R (5′-TGCTGCCTCCCGTAGGAGT-3′) , modified with Adaptor A (CGTATCGCCTCCCTCGCGCCATCAG) and Adaptor B (CTATGCGCCTTGCCAGCCCGCTCAG), separated by the four nucleotides in italics, respectively, for pyrosequencing (Roche). The analysis was performed on DNAs extracted from a set of three larvae sampled in April 2011 (lot A) in the urban area of Palermo, Italy. PCRs for the biological samples and reagent control were carried out in five replicates with 0.6 U Platinum® Taq DNAPolymerase high fidelity (Invitrogen) in 1X PCR buffer, 2 mM MgCl2, 300 nM each primer, 0.24 mM dNTP and 100 ng of DNA in a final volume of 25 μl. Cycling conditions were: 94°C for 5 min, followed by 35 cycles of 94°C for 20 sec, 56°C for 30 sec and 68°C for 40 sec, followed by a final extension at 68°C for 5 min.
Equal volumes of the five reaction products were pooled and purified using the QiAquick Gel Extraction Kit (QIAGEN®). A further purification step was carried out using the Agencourt Ampure XP (Beckman Coulter Genomics), in order to obtain the required pyrosequencing-grade purity, that was assessed by loading a sample in a High Sensitivity DNA chip Agilent 2100 Bioanalyser.
PCR products were mixed for emulsion PCR at one copy per bead using only ‘A’ beads for unidirectional sequencing. Beads were subjected to sequencing on the Roche 454 GS FLX Titanium platform (Roche, Switzerland). Sequences obtained were directly clustered (no trimming was required) with CD-HIT 454 software  using three different similarity threshold: 90%, 95%, and 97%. This software was also used to extract representative cluster consensus sequences. After they were filtered and annotated using the Ribosomal Database Project (RDP) classifier software . Filtering consisted of deleting sequences shorter than 100 bp or containing a number of unknown nucleotides (N) greater than five. Finally, all sequences (clustered plus singletons) were annotated with RDP classifier using default parameters and then parsed to obtain a readable text file in output.
The most abundant unique sequence of each OTU cluster (family or, when possible, species) was selected as representative, then aligned by SINA , mounted in ARB  and subjected to chimera check (before submission in GenBank) by Pintail v. 1.1 software . Rarefaction curves were generated from families of clustered OTUs using EcoSim v.1.2d , separately for each percentage of similarity.
The 97% similarity clustered consensus sequences were deposited in Genbank under accession numbers KC896717-KC896758; raw reads were deposited in NCBI Sequence Read Archive with accession number SRR837401 (reference: BioProject PRJNA196888).
Bacterial isolation and identification
The aerobic cultivable bacterial fraction was analysed on three lots of larvae sampled in April 2011 in Palermo (lots A, B, C). Each pool consisted of three larval guts and their total average weight was 3.68 g (SD: 0.18).
RPW guts were aseptically extracted from each larva, then the content of three guts was pooled, serially diluted in sterile physiological solution, and plated on NA. The plates were incubated for 72 h at 28°C. At the end of the incubation period, colonies were counted and single colonies were streaked to purity on the same fresh medium. The isolates were grouped into OTUs by ARDRA analysis. The whole 16S gene was amplified by colony PCR using the bacterial universal primers fD1 and rD1 , as described elsewhere , and the amplicons were digested using the restriction enzymes AluI and AfaI. Representative isolates of each OTU were randomly chosen for bidirectional sequencing of the 16S rRNA gene. Colonies growing on sterilization control plates were streaked to purity and analysed by ARDRA and 16SrRNA gene partial sequencing. In the same time enrichment cultures in a sorbitol-containing medium at pH 3.5 were set as described by Yamada et al. , for the isolation of acetic acid bacteria (AAB). When microbial growth occurred, the microorganisms were streaked on CaCO3 agar plates and colonies capable of causing clearing of the CaCO3 were selected and identified by partial sequencing of PCR-amplified 16SrRNA gene. Sequences were subjected to NCBI nucleotide BLAST search as described above. Amplified sequences and close relatives were aligned using SILVA alignment tool . Alignment was merged with SSUref_108_Silva_NR database and manually checked with ARB . After alignment, the neighbour-joining algorithm of ARB package was used to generate the phylogenetic trees based on distance analysis for 16S rRNA genes. The robustness of inferred topologies was tested by bootstrap re-sampling using the same distance model (1000 replicates). 16S rRNA gene sequences were deposited in Genbank under accessions number KC584753 to KC584772 (gut isolates), KC763479-80 (cuticle isolates) and KC763478 (AAB enrichment culture isolate).
Recently, just before this manuscript was submitted to this journal, a study on the seasonal variation of the intestinal metagenomes of R. ferrugineus larvae and adults from date palms was published . This study reports that, at the phylum level, Proteobacteria dominate the gut metagenomes of date palm larvae, followed by Tenericutes or Firmicutes depending on the season. The authors identify Klebsiella pneumoniae and Lactococcus lactis as the dominant species of the microbiota. Bacteroidetes are found at negligible levels and the genus Dysgonomonas is not detected. Differences between larvae from date palm and those from Canary palm may be attributed to the host plant species.
The metagenomic analysis carried out by Jia et al.  describes a few ORFs involved in the hydrolisis of cellulolose and high abundance of ORFs related to utilization of disaccharides and simple sugars that are abundantly available in the palm sap . Sap sugars are presumably the main C and energy source for the RPW larvae and its microbiota, that is dominated by fermenting bacteria to obtain several metabolites including lactate and acetate.
Cox ML: Red palm weevil, Rhynchophorus ferrugineus in Egypt. FAO Pl Prot Bul. 1993, 41: 30-31.
Butera G, Ferraro C, Colazza S, Alonzo G, Quatrini P: The culturable bacterial community of frass produced by larvae of Rhynchophorus ferrugineus Olivier (Coleoptera: Curculionidae) in the Canary island date palm. Lett Appl Microbiol. 2012, 54: 530-536. 10.1111/j.1472-765X.2012.03238.x.
Dembilio Ó, Jacas JA: Basic bio-ecological parameters of the invasive Red Palm Weevil, Rhynchophorus ferrugineus (Coleoptera: Curculionidae), in Phoenix canariensis under Mediterranean climate. Bull Entomol Res. 2011, 101: 153-163. 10.1017/S0007485310000283.
Faleiro JR: A review on the issues and management of red palm weevil Rhynchophorus ferrugineus (Coleoptera: Rhynchophoridae) in coconut and date palm during the last one hundred years. Internat J Trop Insect Sci. 2006, 26: 135-154.
Ferry M, Gomez S: The Red Palm Weevil in the Mediterranean area. Palms. 2002, 46: 172-178.
Engel P, Moran NA: The gut microbiota of insects - diversity in structure and function. FEMS Microbiol Rev. 2013, 37: 699-735. 10.1111/1574-6976.12025.
Hongoh Y: Diversity and genomes of uncultured microbial symbionts in the termite gut. Biosci Biotechnol Biochem. 2010, 74: 1145-1151. 10.1271/bbb.100094.
Colman DR, Toolson EC, Takacs-Vesbach CD: Do diet and taxonomy influence insect gut bacterial communities?. Mol Ecol. 2012, 21: 5124-5137. 10.1111/j.1365-294X.2012.05752.x.
Rosenberg E, Zilber-Rosenberg I: Symbiosis and development: the hologenome concept. Birth Defects Res (Part C). 2011, 93: 56-66. 10.1002/bdrc.20196.
Dillon RJ, Dillon VM: The gut bacteria of insects: non-pathogenic interactions. Annu Rev Entomol. 2004, 49: 71-92. 10.1146/annurev.ento.49.061802.123416.
Kaufman MG, Klug MJ: The contribution of hindgut bacteria to dietary carbohydrate utilization by crickets (Orthoptera, Gryllidae). Comp Biochem Physiol A Physiol. 1991, 98: 117-123. 10.1016/0300-9629(91)90588-4.
Chakravorty S, Helb D, Burday M, Connell N, Alland D: A detailed analysis of 16S ribosomal RNA gene segments for the diagnosis of pathogenic bacteria. J Microbiol Methods. 2007, 69: 330-339. 10.1016/j.mimet.2007.02.005.
Tang X, Freitak D, Vogel H, Ping L, Shao Y, Cordero EA, Andersen G, Westermann M, Heckel DG, Boland W: Complexity and variability of gut commensal microbiota in polyphagous lepidopteran larvae. PLoS One. 2012, 7 (7): e36978-10.1371/journal.pone.0036978. doi:10.1371/journal.pone.0036978
Shivaji S, Chaturvedi P, Suresh K, Reddy GS, Dutt CB, Wainwright M, Narlikar JV, Bhargava PM: Bacillus aerius sp. nov., Bacillus aerophilus sp. nov., Bacillus stratosphericus sp. nov. and Bacillus altitudinis sp. nov., isolated from cryogenic tubes used for collecting air samples from high altitudes. Int J Syst Evol Microbiol. 2006, 56: 1465-1473. 10.1099/ijs.0.64029-0.
Suzuki T, Yamasato K: Phylogeny of spore-forming lactic acid bacteria based on 16S rRNA gene sequences. FEMS Microbiol Lett. 1994, 115: 13-17. 10.1111/j.1574-6968.1994.tb06607.x.
Crotti E, Rizzi A, Chouaia B, Ricci I, Favia G, Alma A, Sacchi L, Bourtzis K, Mandrioli M, Cherif A, Bandi C, Daffonchio D: Acetic acid bacteria, newly emerging symbionts of insects. Appl Environ Microbiol. 2010, 76: 6963-6970. 10.1128/AEM.01336-10.
Claesson MJ, Wang QO, O’Sullivan O, Greene-Diniz R, Cole JR, Ross RP, O’Toole PW: Comparison of two next-generation sequencing technologies for resolving highly complex microbiota composition using tandem variable 16S rRNA gene regions. Nucleic Acids Res. 2010, 38 (22): e200-10.1093/nar/gkq873. doi:10.1093/nar/gkq873
Milani C, Hevia A, Foroni E, Duranti S, Turroni F, Lugli A, Sanchez B, Martìn R, Gueimonde M, van Sinderen D, Margolles A, Ventura M: Assessing the fecal microbiota: an optimized ion torrent 16S rRNA gene-based analysis protocol. PLoS One. 2013, 8 (7): e68739-10.1371/journal.pone.0068739. doi:10.1371/journal.pone.0068739
Youssef N, Sheik CS, Krumholz LR, Najar FZ, Roe BA, Elshahed MS: Comparison of species richness estimates obtained using nearly complete fragments and simulated pyrosequencing-generated fragments in 16S rRNA gene-based environmental surveys. Appl Environ Microbiol. 2009, 75: 5227-5236. 10.1128/AEM.00592-09.
Morales-Jiménez J, Zuniga G, Villa-Tanaca L, Hernandez-Rodriguez C: Bacterial community and nitrogen fixation in the red turpentine beetle, Dendroctonus valens LeConte (Coleoptera: Curculionidae: Scolytinae). Microb Ecol. 2009, 58: 879-891. 10.1007/s00248-009-9548-2.
Chandler JA, Morgan Lang J, Bhatnagar S, Eisen JA, Kopp : A bacterial communities of diverse drosophila species: ecological context of a host–microbe model system. PLoS Genet. 2011, 7 (9): e1002272-10.1371/journal.pgen.1002272. doi:10.1371/journal.pgen.1002272
Nagnan P, Cain AH, Rochat D: Extraction et identification des composés volatils de la sève de palmier à huile fermentée (vin de palme) attractifs potentiels pour le charancon du palmier. Oléagineux. 1992, 47: 135-142.
Behar A, Yuval B, Jurkevitch E: Enterobacteria-mediated nitrogen fixation in natural populations of the fruit fly Ceratitis capitata. Mol Ecol. 2005, 14: 2637-2643. 10.1111/j.1365-294X.2005.02615.x.
Chandrasekhar K, Sreevani S, Seshapani P, Pramodhakumari J: A Review on palm wine. J Res Biol Sci. 2012, 2 (1): 33-38.
Santo Domingo JW, Kaufman MG, Klug MJ, Tiedje JM: Characterization of the cricket hindgut microbiota with fluorescently labeled rRNA-targeted oligonucleotide probes. Appl Environ Microbiol. 1998, 64: 752-755.
Wong CN, Ng P, Douglas AE: Low-diversity bacterial community in the gut of the fruitfly Drosophila melanogaster. Environ Microbiol. 2011, 13: 1889-1900. 10.1111/j.1462-2920.2011.02511.x.
Mazza G, Arizza V, Baracchi D, Barzanti GP, Benvenuti C, Francardi V, Frandi A, Gherardi F, Longo S, Manachini B, Perito B, Rumine P, Schillaci D, Turillazzi S, Cervo R: Antimicrobial activity of the red palm weevil Rhynchophorus ferrugineus. Bull Insectol. 2011, 64: 33-41.
Rinke R, Costa AS, Fonseca FP, Almeida LC, Delalibera Júnior I, Henrique-Silva F: Microbial diversity in the larval gut of field and laboratory populations of the sugarcane weevil Sphenophorus levis (Coleoptera, Curculionidae). Genet Mol Res. 2011, 10: 2679-2691. 10.4238/2011.November.1.1.
Hofstad T, Olsen I, Eribe ER, Falsen E, Collins MD, Lawson PA: Dysgonomonas gen. nov. to accommodate Dysgonomonas gadei sp. nov., an organism isolated from a human gall bladder, and Dysgonomonas capnocytophagoides (formerly CDC group DF-3). Int J Syst Evol Microbiol. 2000, 50: 2189-2195. 10.1099/00207713-50-6-2189.
Watanabe K, Miyahara M, Shimoyama T, Hashimoto K: Population dynamics and current-generation mechanisms in cassette-electrode microbial fuel cells. Appl Microbiol Biotechnol. 2011, 92: 1307-1314. 10.1007/s00253-011-3598-3.
Gupta AK, Nayduch D, Verma P, Shah B, Ghate HV, Patole MS, Shouche YS: Phylogenetic characterization of bacteria in the gut of house flies (Musca domestica L.). FEMS Microbiol Ecol. 2012, 79: 581-593. 10.1111/j.1574-6941.2011.01248.x.
Campbell BC, Bragg TS, Turner CE: Phylogeny of symbiotic bacteria of four weevil species (Coleoptera:Curculionidae) based on analysis of 16S ribosomal DNA. Insect Biochem Molec Biol. 1992, 22: 415-421. 10.1016/0965-1748(92)90136-3.
Tully JG, Whitcomb RF, Hackett KJ, Williamson DL, Laigret F, Carle P, Bové JM, Henegar RB, Ellis NM, Dodge DE, Adams J: Entomoplasma freundtii sp. nov., a new species from a green tiger beetle (Coleoptera: Cicindelidae). Int J Syst Bacteriol. 1998, 48: 1197-1204. 10.1099/00207713-48-4-1197.
Yu H, Wang Z, Liu L, Xia Y, Cao Y, Yin Y: Analysis of the intestinal microflora in Hepialus gonggaensis larvae using 16S rRNA sequences. Curr Microbiol. 2008, 56: 391-396. 10.1007/s00284-007-9078-4.
Suen G, Scott JJ, Aylward FO, Adams SM, Tringe SG, Pinto-Tomás AA, Foster CE, Pauly M, Weimer PJ, Barry KW, Goodwin LA, Bouffard P, Li L, Osterberger J, Harkins TT, Slater SC, Donohue TJ, Currie CR: An insect herbivore microbiome with high plant biomass-degrading capacity. PLoS Genet. 2010, 6 (9): e1001129-10.1371/journal.pgen.1001129. doi:10.1371/journal.pgen.1001129
Paoletti MG, Mazzon L, Martinez-Sañudo I, Simonato M, Beggio M, Dreon AL, Pamio A, Brilli M, Dorigo L, Engel AS, Tondello A, Baldan B, Concheri G, Squartini A: A unique midgut-associated bacterial community hosted by the cave beetle Cansiliella servadeii (Coleoptera: Leptodirini) reveals parallel phylogenetic divergences from universal gut-specific ancestors. BMC Microbiol. 2013, 13: 129-10.1186/1471-2180-13-129.
Guarino S, Lo Bue P, Peri E, Colazza S: Responses of Rhynchophorus ferrugineus adults to selected synthetic palm esters: electroantennographic studies and trap catches in an urban environment. Pest Manag Sci. 2011, 67: 77-81. 10.1002/ps.2035.
Broderick NA, Goodman RM, Handelsman J, Raffa KF: Effect of host diet and insect source on synergy of gypsy moth (Lepidoptera: Lymantriidae) mortality to Bacillus thuringiensis subsp. kurstaki by zwittermicin A. Environ Entomol. 2003, 32: 387-391. 10.1603/0046-225X-32.2.387.
Abe F, Ohkusu M, Kawamoto S, Sone K, Hata K: Isolation of yeasts from palm tissues damaged by the red palm weevil and their possible effect on the weevil overwintering. Mycoscience. 2010, 51: 215-223. 10.1007/S10267-009-0029-8.
Hassan AA, Akineden O, Kress C, Estuningsih S, Schneider E, Usleber E: Characterization of the gene encoding the 16S rRNA of Enterobacter sakazakii and development of a species-specific PCR method. Int J Food Microbiol. 2007, 116: 214-220. 10.1016/j.ijfoodmicro.2006.12.011.
Toledo AV, Alippi AM, de Remes Lenicov AM: Growth inhibition of Beauveria bassiana by bacteria isolated from the cuticular surface of the corn leafhopper, Dalbulus maidis and the planthopper, Delphacodes kuscheli, two important vectors of maize pathogens. J Insect Sci. 2011, 11: 1-13.
Yamada Y, Katsura K, Kawasaki H, Widyastuti Y, Saono S, Seki T, Uchimura T, Komagata K: Asaia bogorensis gen. nov., sp. nov., an unusual acetic acid bacterium in the alpha-Proteobacteria. Int J Syst Evol Microbiol. 2000, 50: 823-829. 10.1099/00207713-50-2-823.
Doolittle M, Raina A, Lax A, Boopathy R: Presence of nitrogen fixing Klebsiella pneumoniae in the gut of the Formosan subterranean termite (Coptotermes formosanus). Bioresour Technol. 2008, 99: 3297-3300. 10.1016/j.biortech.2007.07.013.
Pittman GW, Brumbley SM, Allsopp PG, O’Neill SL: “Endomicrobia” and other bacteria associated with the hindgut of Dermolepida albohirtum larvae. Appl Environ Microbiol. 2008, 74: 762-767. 10.1128/AEM.01831-07.
Muyzer G, DeWaal EC, Uitterlinden AG: Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes coding for 16S rRNA. Appl Environ Microbiol. 1993, 59: 695-700.
Hamady M, Walker JJ, Harris JH, Gold NJ, Knight R: Error-correcting barcoded primers allow hundreds of samples to be pyrosequenced in multiplex. Nat Methods. 2008, 5: 235-237. 10.1038/nmeth.1184.
Niu B, Fu L, Sun S, Li W: Artificial and natural duplicates in pyrosequencing reads of metagenomic data. BMC Bioinforma. 2010, 13: 11-187.
Cole JR, Wang Q, Fish JA, Chai B, McGarrell DM, Sun Y, Brown CT, Porras-Alfaro A, Kuske CR, Tiedje JM: Ribosomal database project: data and tools for high throughput rRNA analysis. Nucl Acids Res. 2014, 41 (Database issue): D633-D642. doi:10.1093/nar/gkt1244
Pruesse E, Peplies J, Glöckner FO: SINA: accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics. 2012, 28: 1823-1829. 10.1093/bioinformatics/bts252.
Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar , Buchner A, Lai T, Steppi S, Jobb G, Förster W, Brettske I, Gerber S, Ginhart AW, Gross O, Grumann S, Hermann S, Jost R, König A, Liss T, Lüssmann R, May M, Nonhoff B, Reichel B, Strehlow R, Stamatakis A, Stuckmann N, Vilbig A, Lenke M, Ludwig T, Bode A, Schleifer KH: ARB: a software environment for sequence data. Nucleic Acids Res. 2004, 32: 1363-1371. 10.1093/nar/gkh293.
Ashelford KE, Chuzhanova NA, Fry JC, Jones AJ, Weightman AJ: At least 1 in 20 16S rRNA sequence records currently held in public repositories is estimated to contain substantial anomalies. Appl Environ Microbiol. 2005, 71: 7724-7736. 10.1128/AEM.71.12.7724-7736.2005.
Entsminger GL: EcoSim Professional: Null Modelling Software for Ecologists, Version 1. 2012, Acquired Intelligence Inc., Kesey-Bear, & Pinyon Publishing, http://garyentsminger.com/ecosim/index.htm,
Weisburg WG, Barns SM, Pelletier DA, Lane DJ: 16S Ribosomal DNA amplification for phylogenetic study. J Bacteriol. 1991, 173: 697-703.
Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glöckner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007, 35: 7188-7196. 10.1093/nar/gkm864.
Jia S, Zhang X, Zhang G, Yin A, Zhang S, Li F, Wang L, Zhao D, Yun Q, Tala , Wang J, Sun G, Baabdullah M, Yu X, Hu S, Al-Mssallem IS, Yu J: Seasonally variable intestinal metagenomes of the red palm weevil (Rhynchophorus ferrugineus). Environ Microbiol. 2013, 15: 3020-3029.
The authors thank Maria Grazia Cusimano, Rosella Maggio and Flavia Contino for technical assistance in bacterial isolation, DNA extraction and amplification, and control of DNA quality for pyrosequencing, respectively. This work was partially financed by a grant from the Italian Ministry of Education (PRIN Program 2008 prot. 200847CA28-002) and by the University of Palermo (project FFR 2012 ATE-0322 N.2785).
The authors declare that they have no competing interests.
MT projected and carried out the microbiological and molecular analyses, EM performed the bioinformatic analyses, BM identified and collected the insects in the field and manipulated them for the gut microbiota analyses, SC constructed the phylogeny trees and helped to draft the manuscript, PQ conceived and coordinated the study and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Phoenix canariensis infested by Rhynchophorous ferrugineus (A and B); different infested palms cut in the higher part are shown. Larvae of the red palm weevil (RPW) Rhynchophorus ferrugineus, found inside the body of the infested palm (C). Female adult specimen of Rhynchophorus ferrugineus Olivier (Coleoptera, Curculionidae, Rhynchophorinae) (D). (PDF 171 KB)
Additional file 4: Phylogenetic tree of 16S rRNA gene amplicons clustered at 97% consensus. The tree was constructed by neighbour-joining method and Jukes Cantor distance matrix using the arb software. Bootstraps were calculated over 1000 random repetitions: values >60 and < =75 are shown as open circles, whereas values >75 are shown as filled circles. Sequences obtained in this study are indicated in bold. The scale bar represents 10% sequence divergence. (PDF 231 KB)
Additional file 5: Phylogenetic tree of 16S rDNA sequences of RPW gut isolates and related sequences, as determined by distance Jukes-Cantor analysis. One thousand boostrap analyses were conducted and values greater than 60% are reported. Two Archaea sequences of Methanopirus kandleri and Korarchaeum cryptophilum were used as outgroup. The scale bar represents the expected number of changes per nucleotide position. Colors indicate the isolation site or the isolation procedure described in this work and in . Blue: RPW gut isolates on NA; Red: frass bacteria; Green: palm bacterial endophytes; Fuchsia: gut isolates obtained from enrichment cultures on CMC; Yellow: larval cuticle bacteria isolated from sterilization control plates. Isolate_RPWenrichAAB* was isolated from the RPW larval gut from enrichment cultures set for for Acetic Acid Bacteria isolation . (PDF 271 KB)
Additional file 6: Identities of 16S rRNA gene sequences of cellulolytic bacterial strains isolated from the RPW gut by enrichment cultures  with the RDP 454 consensus sequences clustered at 97% ID of the RPW microbiota. Only identities above 90% are shown. (PDF 25 KB)
Authors’ original submitted files for images
About this article
Cite this article
Tagliavia, M., Messina, E., Manachini, B. et al. The gut microbiota of larvae of Rhynchophorus ferrugineus Oliver (Coleoptera: Curculionidae). BMC Microbiol 14, 136 (2014). https://doi.org/10.1186/1471-2180-14-136