Next Article in Journal
Populations and Host/Non-Host Plants of Spittlebugs Nymphs in Olive Orchards from Northeastern Portugal
Previous Article in Journal
Bacterial Communities of Ixodes scapularis from Central Pennsylvania, USA
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Divergence in Gut Bacterial Community Among Life Stages of the Rainbow Stag Beetle Phalacrognathus muelleri (Coleptera: Lucanidae)

1
School of Resources and Environmental Engineering, Anhui University, Hefei 230601, China
2
Anhui Province Key Laboratory of Wetland Ecosystem Protection and Restoration, Hefei 230601, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this article.
Insects 2020, 11(10), 719; https://doi.org/10.3390/insects11100719
Submission received: 19 September 2020 / Revised: 16 October 2020 / Accepted: 17 October 2020 / Published: 21 October 2020

Abstract

:

Simple Summary

Phalacrognathus muelleri is naturally distributed in Queensland (Australia) and New Guinea, and this species can be successfully bred under artificial conditions. In this study, we compared gut bacterial community structure among different life stages. There were dramatic shifts in gut bacterial community structure between larvae and adults, which was probably shaped by their diet. The significant differences between early instar and final instars larvae suggested that certain life stages are associated with a defined gut bacterial community. Our results contribute to a better understanding of the potential role of gut microbiota in a host’s growth and development, and the data will benefit stag beetle conservation in artificial feeding conditions.

Abstract

Although stag beetles are popular saprophytic insects, there are few studies about their gut bacterial community. This study focused on the gut bacterial community structure of the rainbow stag beetle (i.e., Phalacrognathus muelleri) in its larvae (three instars) and adult stages, using high throughput sequencing (Illumina Miseq). Our aim was to compare the gut bacterial community structure among different life stages. The results revealed that bacterial alpha diversity increased from the 1st instar to the 3rd instar larvae. Adults showed the lowest gut bacterial alpha diversity. Bacterial community composition was significantly different between larvae and adults (p = 0.001), and 1st instar larvae (early instar) had significant differences with the 2nd (p = 0.007) and 3rd (p = 0.001) instar larvae (final instar). However, there was little difference in the bacterial community composition between the 2nd and 3rd instar larvae (p = 0.059). Our study demonstrated dramatic shifts in gut bacterial community structure between larvae and adults. Larvae fed on decaying wood and adults fed on beetle jelly, suggesting that diet is a crucial factor shaping the gut bacterial community structure. There were significant differences in bacterial community structure between early instar and final instars larvae, suggesting that certain life stages are associated with a defined gut bacterial community.

Graphical Abstract

1. Introduction

Insects are the most diverse and abundant class of animals, living in multiple habitats, and feeding on various substrates [1,2]. The guts of insects are colonized by diverse microorganisms that play integral roles in their hosts, including affecting the hosts’ metabolism, providing essential amino acids, vitamins and nitrogen for the host, promoting efficient digestion of nutrient-poor diets and recalcitrant foods, aiding the defense and detoxification ability and protecting hosts from potentially harmful microbes [3]. Associations between microorganisms and insect hosts are widespread in nature. Studies on Cerambycidae species revealed that their gut bacterial community is dominated by diet type, developmental stage and gut compartment [4]. Insect gut microorganisms are recognized to originate from their environment and diet [5]. Beetles (Coleoptera) are the most diverse group of insects [6]. Almost all beetles undergo holometabolism, including a pupal stage compared with incomplete metamorphosis. The study of Agrilus mali (Coleoptera: Buprestidae) and Popillia japonica (Coleoptera: Scarabaeidae) has suggested that the metabolic activity of the intestine decreases in the pupal stage and morphology changes during the insects’ metamorphosis influence the associated bacteria communities [3,7]. Moreover, larvae and adults feed on different foods, thus, occupying distinct ecological niches in the carbon cycle of forest ecosystems [8]. The guts of beetles are excellent models to study the variety of gut microorganisms that contribute to hosts’ digestion, detoxification, development, pathogen resistance and physiology [1,2,5,9].
Previous studies have focused on the gut bacterial community of coleopteran pests to advance progress in pest control [10,11], including invasive bark beetles (Scolytidae) and longhorn beetles (Cerambycidae), that mainly feed on tree xylem and phloem. These species are considered to be serious pests of forests, as they grow and develop while feeding on lignocellulose [12,13]. Many studies have investigated gut microbial species that are critical for the development and survival of the hosts [14,15]. The aim of these studies was to identify how to damage and alter the core microbiota to cause the death of pests, thereby developing new pest control strategies [16,17,18,19]. Further, gut bacteria in these pests contributed to understand the degrading functions by digesting cellulose and lignin (i.e., wood, litter, and humus) to promote host growth [10,11,12,20,21].
In addition, some beetle gut microbial studies have involved dung beetles (Scarabaeidae) and burying beetles (Silphidae). Adult dung beetles feed primarily on nutritionally rich dung particles, while larvae consume coarser dung particles with a higher C/N ratio [22]. There are studies comparing the gut microbiota differences between larvae and adults [22], suggesting that diet is a factor in the composition of gut microbiota. Individual studies have focused on maternal transmission [23] to explore the similarity of gut microbes between larvae and adult females or males [24]. Burying beetles use small vertebrate carcasses as a source of food for their developing larvae [25]. Studies of their gut microflora have focused on the functions of the gut microbiota, such as the digestion of carrion and the detoxification of microbial toxins [26,27].
There is little information about the gut bacterial community of stag beetles, although they are one of the most important beetles in forests. The larval stages live in and feed on decaying wood, and the adults utilize fermented tree sap and over–ripe fruits [28]. To date, studies of microorganisms associated with stag beetles involve fungi, particularly those transported externally on the exoskeleton or in specialized structures known as mycangia [29,30,31,32]. Fungi in the larval gut could contribute to forest matter cycling through turning a low nutrition food source into available nutrients [32,33]. Likewise, little work has been done on the gut bacteria of stag beetles using Illumina MiSeq high-throughput sequencing, with the exception of a study by Jiang [34] who analyzed the gut bacterial structure in Odontolabis fallaciosa adults. However, little work has focused on gut bacterial communities of Phalacrognathus muelleri (Coleptera: Lucanidae) during successive life stages with different diets.
This study represents the first comprehensive analysis of the gut bacterial community of the rainbow stag beetle P. muelleri during different life stages using high- throughput sequencing. This species is naturally distributed in Queensland (Australia) and New Guinea. The beetle has become a very popular pet due to the colorful body and interesting male mandibles and is successfully bred under artificial conditions in many countries. We investigated 1st, 2nd, and 3rd instar larvae and adults that were raised under standard rearing conditions (see Materials and Methods). Specifically, the main goal of this study was to explore the differences of gut bacterial community composition and diversity during successive life stages of P. muelleri with different diets.

2. Materials and Methods

2.1. Sample Collection and Rearing Conditions

All the samples of P. muelleri were obtained from the Mu–Ye Insect Company (Lishui, Zhejiang, China), including the 1st instar (L1, eight individuals), the 2nd instar (L2, twelve individuals), the 3rd instar (L3, nine males and nine females) and the adults (six males and ten females) (Table S1). These samples were artificially reared under constant conditions (temperature: 18–22 °C). Larvae were reared in pudding boxes containing decaying wood with a 50% water content; the larvae were frequently observed as to feeding status, and the feed was regularly replaced and continual1y rehydrated. Adults were reared on beetle jelly and their feed was regularly replaced.

2.2. Sample Dissection

Before dissection, all samples were disinfected for 3 min with 70% ethanol, and then rinsed with distilled water to ensure no contamination on the body surface [35,36]. Beetles were transported to a horizontal clean bench, and dissected in 10–fold diluted phosphate–buffered solution (PBS) (500 mL, NaCl 1.37 M, KCL 26.8 mM, Na2HPO4 81.0 mM, KH2PO4 17.6 mM, pH 7.2–7.4) [37] under a stereomicroscope. The midgut and hindgut were removed using sterile fine tip forceps and placed into 2 mL Lysing Matrix E under sterile conditions to avoid contamination. In addition, 1st instar, 2nd instar, 3rd instar and adults were dissected at five days, ten days, ten days, and ten days, respectively.

2.3. Microbial DNA Extraction and PCR Amplification

DNA extractions were performed on each gut sample using QIAamp FAST DNA Stool Mini Kits (Qiagen Inc. Valencia, CA, USA) according to the operating instructions. The extracted DNA was dissolved in 75 μL of elution buffer, quantified by NanoDrop ND–1000 (Thermo Scientific, Wilmington, DE, USA) and stored at −20 °C. The purified DNA of each sample served as an amplification template. The primer sets 515F (5′-GTGCCAGCMGCCGCGG-3′) and 907R (5′-GGACTACHVGGGT WTCTAAT-3′) [38], were used to amplify the V4 to V5 variable regions of the bacterial 16S rRNA gene fragments for the Illumina Mi–Seq platform (PE 250) at Majorbio (Shanghai, China). The PCR reactions were carried out in 20 μL reaction mixtures containing 4 μL 5 × FastPfu Buffer, 2 μL 2.5 mM dNTPs, 0.8 μL 5 uM of each primer, 0.4 μL DNA polymerase, 0.2 μL BSA, 10 ng template DNA and deionized μltrapure water (to 20 μL). The PCR conditions were as follow: initial at 95 °C (3 min) follows by 27 cycles at 95 °C (30 s), 55 °C (30 s), and 72 °C (45 s), and final extension at 72 °C for 10 min. To check for contamination, PCR negative controls were performed without added DNA template. Negative PCR controls did not contain detectable PCR product and were not processed for sequencing. Triplicate reaction mixtures per sample were pooled together and purified using an agarose gel DNA purification kit (TaKaRa, Bio Inc., Kusatsu, Japan). The PCR products were pooled in equimolar amounts (10 pg for each sample) before sequencing.

2.4. Processing of Sequence Data

The raw data were processed by the Quantitative Insights Into Microbial Ecology (QIIME v.1.9 software [39]. Low–quality sequences (below an average quality score of 30 and the length < 250 bp) were deleted. High–quality sequences were clustered into Operational Taxonomic Units (OTUs; 97% similarity; de novo approach). The chimeras were eliminated by the USEARCH (V.1.8.0). The ribosome database project classifier selects the most abundant sequence in each OTU as the representative sequence [40], which was aligned by PyNAST [39]. In order to perform similar sequencing and homogenization among samples, we used randomly selected subsets of 26,000 sequences (lowest sequence read depth; repetition with 20 times) per sample to compare bacterial community composition and diversity for all samples [41].

2.5. Statistical Analysis

The differences in alpha diversity and the relative abundance of dominant phyla among life stages were based on one-way analysis of variance (ANOVA) with Tukey Honestly Significant Difference (HSD) Post-Hoc testing (SPSS 20.0 for Windows, Chicago, IL, USA) [41]. Linear discriminant analysis (LDA) effect size (LEfSe) was used to identify intestinal bacterial taxa with significant differences among host stages. This method uses the non–parametric Kruskal–Wallis test with default settings (an alpha value of 0.05 and an effect size threshold of 2) in a rank sum test to identify biomarkers [42]. The differences in bacterial community composition between different stages were analyzed by non–metric multidimensional scaling (NMDS) and analysis of similarity (ANOSIM; permutations = 999) using the vegan package (Version 2.0–10 [43] in R v.2.8.1 [44,45]. The contribution of bacterial OTUs (operational taxonomic units) to the differences between stages was analyzed by SIMPER using the vegan package in R software. Indicator analysis was performed using the labdsv package in R software [46]. Bacterial diversity and the relative abundance of dominant phyla were analyzed by one-way ANOVA.

2.6. Data Availability Statement

The raw data were submitted to the Sequence Read Archive (SRA) of NCBI under the accession number SUB8060438.

3. Results

3.1. Intestinal Bacterial Alpha Diversity

In this study, a total of 1,997,394 quality–filtered bacterial sequences were retrieved from the 54 samples for the primer pair F515/R907, ranging from 26,163 to 61,214 sequences per sample (Table S2). A total of 1136 bacterial OTUs were detected (relative abundance (%) <0.01%), with all samples ranging from 200 to 5912 (97% similarity), and 12 OTUs were shared across all life stages. There were 553, 806, 913 and 229 OTUs in L1, 1st instar; L2, 2nd instar; L3, 3rd instar and Ad, Adults, respectively. In addition, there were unique OTUs within each stages, especially for adults, where the unique OTUs accounted for 72.4% with the 229 distinct OTUs observed, and the unique OTUs of 1st instar, 2nd instar, 3rd instar accounted for 0.7%, 3.7% and 11.9%, separately (Figure 1).
Bacterial alpha diversity indicators (i.e., OTU richness, Shannon index, evenness, and phylogenetic diversity) were calculated at a depth of 26,000 randomly selected sequences per sample. According to the index estimation results, the bacterial alpha diversity of the larvae was higher than that of the adults, and the bacterial diversity of larvae increased with successive instars (Figure 2). In different sexes alpha diversity analysis, with the exception of OTU richness, which was relatively higher in females than in males of L3, all other indicators showed little differences in L3. There were no significant differences in the alpha diversity of adults between females and males (Figure S1).
LEfSe analysis further identified specific bacterial taxa that were differentially abundant across all life stages. The results showed that bacteria in four classes (i.e., Spartobacteria, Solibacteres, Thermoplasmata, Methanopyri) and three orders (i.e., Chthoniobacterale, Solibacterales, Methanopyrales) were significantly more abundant in L1. Bacteria in two phyla (i.e., Planctomycetales, Acidobacteria) and one class (i.e., Planctomycetia) were significantly more abundant in L2. Bacteria from two classes (i.e., Bacilli, Methanobacteri) and five orders (i.e., Myxococcales, Thiotrichales, Rhodospirillales, Methanobacteriales, Legionellalesi) were significantly more abundant in the gut of L3 (Figure 3). SIMPER analysis revealed that OTU_46981 (Sporomusa; 2.91%), OTU_63703 (Ruminococcaceae; 2.14%) made primary contributions to community differences between L1 and L2; OTU_46981 (Sporomusa; 2.94%) and OTU_63703 (Ruminococcaceae; 2.17%) were the main OTUs responsible for the differences between L1 and L3; OTU_46981 (Sporomusa; 2.89%) and OTU_38714 (Dysgonomonas; 7.72%) contributed to differences in the bacterial community composition between L1 and Ad; OTU_46981 (Sporomusa; 0.80%) and OTU_9882 (Dysgonomonas; 1.74%) produced the significant differences in composition between L2 vs. L3; OTU_9882 (Dysgonomonas; 1.60%) and OTU_38714 (Dysgonomonas; 7.71%) were the principal OTUs contributing to differences between L2 and Ad, OTU_38714 (Dysgonomonas; 7.72%) and OTU_53348 (Lactobacillales; 4.04%), contributed to the differences in bacterial community composition between L3 and Ad (Table S3). Indicator analysis was used to identify bacterial OTUs that were specifically associated with various life stages. There were different indicator species in the L1, L2, L3 and Ad (Table S4), and some indicators were different in females and males (Table S5). We also compared the bacterial community composition across different life stages using NMDS analysis; significantly different communities were present in larvae and adults (ANOSIM: p = 0.001, Table 1; Figure 4). And 1st instar larvae (early instar) had significant differences with the 2nd (p = 0.007) and 3rd (p = 0.001) instar larvae (final instar) (Table 1; Figure 4).

3.2. Intestinal Bacterial Community Structure

The dominant intestinal bacterial phyla across all samples were Firmicutes (74.72%), Proteobacteria (12.67%), Bacteroidetes (10.81%) and Tenericutes (1.80%) (Figure 5). These dominant phyla had significant differences in relative abundance in larvae and adults except for Tenericutes (one–way ANOVA: p < 0.05), and there were no significant differences among the three larval instars (one–way ANOVA: p > 0.05). The relative abundance of Firmicutes showed a significant decrease, and the relative abundances of Proteobacteria, and Bacteroidetes showed significant increases in adults relative to larvae (one–way ANOVA: p < 0.05) (Figure 5). There were no significant differences in dominant phyla between females and males of L3. Proteobacteria showed significant differences between male and female adults, with a higher relative abundance in males (Figure S2).

4. Discussion

We investigated the gut bacterial community of P. muelleri across different life stages under artificial breeding conditions through Illumina MiSeq high-throughput sequencing. This is the first study of the gut bacteria of P. muelleri focusing on different life stages. Four phyla—Firmicutes, Proteobacteria, Bacteroidetes and Tenericutes—were predominant in the gut bacteria across all stages, which is consistent with prior studies in other beetles (i.e., bark beetles, longhorn beetles, herbivorous beetles and burying beetles) [12,13,26,47]. Firmicutes and Proteobacteria were the dominant phyla in Anoplophora glabripennis [12], Hylobius abietis [48], Monochamus alternatus and Psacothea hilaris [4], Nicrophorus vespilloides [26], and Popillia japonica [7]. Firmicutes, Proteobacteria and Bacteroidetes were also the dominant phyla in Holotrichia parallela [47], Dendroctonus valens [13], and Dendroctonus rhizophagus [10], and Tenericutes also was the dominant phylum in Odontolabis fallaciosa (Lucanidae) [34]. Prior research demonstrated that intestinal Firmicutes play significant role in the degradation of complex plant carbohydrates. Bacteroides species are able to degrade diverse plant polysaccharides [49], thereby improving the host’s ability to digest food. Proteobacteria contribute to nitrogen fixation and food metabolism to keep the host healthy [50]. The aggregation of these communities in insect gut is dependent on physicochemical conditions of gut compartments [5], such as the available oxygen and pH [4].
In our study, the larval alpha diversity was significantly higher than that of adults (Figure 2), in consistent with a previous study found in species of Melolontha hippocastani [51]. Intestinal bacteria communities are considered as originate from the environment and diet [3]. Thus, the high diversity displayed in larvae may due to the vast microbial colonized via specific nutritional complementation from different diet [50,52]. Moreover, adults rarely feed after emergence, even if enough beetle jelly is given. Among the three larval instars of P. muelleri, the alpha diversity increased in successive instars, a result in accordance with gut bacteria in the Melolontha hippocastani.
According to the NMDS analysis and ANOSIM analysis, gut bacterial community composition was significantly different between larvae and adults (ANOSIM: p = 0.001, Table 1, Figure 4). During the growth of P. muelleri, the dominant phyla of Firmicutes decreased its relative abundance, while the relative abundance of Bacteroidetes and Proteobacteria were increased in adults (Figure 5). Previous studies have demonstrated that gut bacteria communities can be influenced by host diet [14,21,35,53]. Among the three larval instars of P. muelleri, their gut bacterial structures reflected significant divergence. Gut bacterial community composition was significantly different and between L1 and L3 (ANOSIM: p = 0.001, Figure 4) and between L1 and L2 (ANOSIM: p = 0.007, Figure 4), while the difference was relatively small between L2 and L3 (ANOSIM: p = 0.059, Figure 4). We inferred that the compositional differences were related to food intake. Chen [54] also suggested that the physiological and biochemical conditions within the larval alimentary tract affect bacterial community structure. Different food intake and content among three instars could affect the intestinal microbial composition, as in scarab beetle larvae [55]. Relatively speaking, as hosts grew, the food intake increased.
Notably, understanding the function among intestinal bacteria is important to explore the complexity host development process. We observed Ruminococcaceae, Veillonellaceae, Christensenellaceae and Lachnospiraceae was dominated in the larval stages (Table S4), three families that were previously observed to co-occur in the bacterial gut microbiota of the longhorned beetle and bark beetle larvae [54]. They can help their host degrade the decaying wood and metabolize cellulose to promote larval digestion. Such as Ruminococcaceae are vital to degradation of lignocellulose [56]. Christensenellaceae and Lachnospiraceae are important for the degradation of plant material and cellulose catabolism [49,57], and Christensenellaceae is related to host health [58]. Interestingly, the family Veillonellaceae was considered as a probiotic providing benefits for energy balance and producing various volatile fatty acids to lower gut pH [59], which is conducive to host gut fitness. Enterobacteriaceae was dominanted in the adults (Table S5), in consistent with a previous study found in bark beetle [37]. Bacteria in this family can ferment glucose, reduce nitrates to nitrites, and degrade sucrose to provide nutrients for hosts [60].
Likewise, bacteria in the five genera Sporomusa, Candidatus, Coprococcus, Prevotella and Turicibacter were isolated from larvae (Table S4). These groups contribute to external immunity and food fermentation [5,61]. Members of genera Dysgonomonas, Enterococcus are abundant in the adult stages, similar in the Odontolabis fallacios adults [34], Lactococcus, Trabulsiella, Neorickettsia, casseliflavus, Selenomonas were also found in the termites gut [61]. They were primarily involved in lignocellulose decomposition, food fermentation and enhancement of metabolic capacity of the host [4,62,63]. This information illustrated that different diets cause differences in gut bacteria, a result that has been reported in the scarab beetle, Popillia japonica [7], in the fruit fly Drosophila suzukiiis and in the pine weevil Hylobius abietis [52,64].
We further considered the influence of gender on intestinal bacteria. In L3, there was a barely significant difference between males and females (ANOSIM: p = 0.047). Adults also differed little between males and females (ANOSIM: p = 0.578), supporting similar findings in bark beetles from earlier studies [13,21]. This suggested that male and female individuals subjected to standard rearing conditions harbored communities that were highly conserved in structure and membership. We suspect that this may be due to: (i) the beetles being bred in captivity, having a homogeneous diet, and subject to the same living conditions; (ii) the beetles were raised alone with no opportunistic communication. Both factors can affect colonization by gut microbiota [65]. In addition, P. muelleri in our study had not yet attained sexual maturity (this usually takes 30 to 60 days, while our sampling was done at 10 days after emergence), and this may have affected microbial colonization.

5. Conclusions

Our study revealed the structure of the gut bacterial community in different life stages of the rainbow stag beetle, P. muelleri. Results showed that bacterial community composition was significantly different between larvae and adults, and between early instar and final instar larvae. Diet and life stages can thus collectively influence the gut bacterial community composition. This research provides a basis for subsequent studies on the roles of these intestinal bacteria in stag beetle development and ecology. Results also contribute to a better understanding of the potential role of gut microbiota in a host’s growth and development, and the data will benefit stag beetle conservation in artificial feeding conditions. However, there were certain limitations in this research, as the level of biological replication was relatively low. In addition, we did not sequence the intestinal bacterial in pupae. These limitations should be addressed in future studies. In addition, future research will develop more complex bioinformatics tools on the basis of high-throughput sequencing to analyze the undiscovered microorganisms in the insect intestines, and gain insight into their functions for their hosts.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4450/11/10/719/s1, Figure S1. Intestinal bacterial alpha-diversity in different sexes, Figure S2. Relative abundance of bacterial taxa at the phylum level in the guts of different sexes of P. muelleri, Table S1. All sample materials of Phalacrognathus Muelleri, after removing outlier sample values, Table S2. Gut bacterial sequences across the samples, Table S3. SIMPER analysis showing the contribution of bacteria OTUs to the differences in bacterial community composition in the gut depending on different life stages, Table S4. Indicator species in different life stages, Table S5. Indicator species in different sexes, Table S6. The total OTU table.

Author Contributions

X.W. designed the project. M.W. and X.X. analyzed the data and wrote the paper. M.W. performed the experiments. All authors gave final approval for publication. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 31872276, 31801989, 31572311).

Acknowledgments

We thank Qing-Ye Sun and Yang Li for their assistance at the beginning of this study. We appreciate our lab colleagues Yu-Yan Cao, Le-le Jing, Yong-Ning Zhai for their useful discussion. We also thank Sai Huang for his kind help.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rajagopal, R. Beneficial interactions between insects and gut bacteria. Indian J. Microbiol. 2009, 49, 114–119. [Google Scholar] [CrossRef] [Green Version]
  2. Dillon, R.J.; Dillon, V.M. The gut bacteria of insects: Nonpathogenic interactions. Annu. Rev. Entomol. 2004, 49, 71–92. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, Z.; Jiao, S.; Li, X.; Li, M. Bacterial and fungal gut communities of Agrilus mali at different developmental stages and fed different diets. Sci. Rep. 2018, 8. [Google Scholar] [CrossRef] [Green Version]
  4. Kim, J.M.; Choi, M.Y.; Kim, J.W.; Lee, S.A.; Ahn, J.H.; Song, J.; Kim, S.H.; Weon, H.Y. Effects of diet type, developmental stage, and gut compartment in the gut bacterial communities of two Cerambycidae species (Coleoptera). J. Microbiol. 2017, 55, 21–30. [Google Scholar] [CrossRef] [PubMed]
  5. Engel, P.; Moran, N.A. The gut microbiota of insects-diversity in structure and function. FEMS Microbiol. Rev. 2013, 37, 699–735. [Google Scholar] [CrossRef] [PubMed]
  6. Lawrence, J.F.; Ślipiński, A.; Seago, A.E.; Thayer, M.K.; Newton, A.F.; Marvaldi, A.E. Phylogeny of the Coleoptera Based on Morphological Characters of Adults and Larvae. Ann. Zool. 2011, 61, 1–217. [Google Scholar] [CrossRef]
  7. Chouaia, B.; Goda, N.; Mazza, G.; Alali, S.; Florian, F.; Gionechetti, F.; Callegari, M.; Gonella, E.; Magoga, G.; Fusi, M.; et al. Developmental stages and gut microenvironments influence gut microbiota dynamics in the invasive beetle Popillia japonica Newman (Coleoptera: Scarabaeidae). Environ. Microbiol. 2019, 21, 4343–4359. [Google Scholar] [CrossRef]
  8. Price, P.W.; Denno, R.F.; Eubanks, M.D.; Finke, D.L.; Kaplan, I. Insect Ecology: Behavior, Populations and Communities; Cambridge University Press: New York, NY, USA, 2011; pp. 1–71. [Google Scholar] [CrossRef]
  9. Douglas, A.E. Multiorganismal insects: Diversity and function of resident microorganisms. Annu. Rev. Entomol. 2015, 60, 17–34. [Google Scholar] [CrossRef] [Green Version]
  10. Briones-Roblero, C.I.; Hernandez-Garcia, J.A.; Gonzalez-Escobedo, R.; Soto-Robles, L.V.; Rivera-Orduna, F.N.; Zuniga, G. Structure and dynamics of the gut bacterial microbiota of the bark beetle, Dendroctonus rhizophagus (Curculionidae: Scolytinae) across their life stages. PLoS ONE 2017, 12, e0175470. [Google Scholar] [CrossRef]
  11. Mohammed, W.S.; Ziganshina, E.E.; Shagimardanova, E.I.; Gogoleva, N.E.; Ziganshin, A.M. Comparison of intestinal bacterial and fungal communities across various xylophagous beetle larvae (Coleoptera: Cerambycidae). Sci. Rep. 2018, 8, 10073. [Google Scholar] [CrossRef]
  12. Geib, S.M.; Jimenez-Gasco Mdel, M.; Carlson, J.E.; Tien, M.; Jabbour, R.; Hoover, K. Microbial community profiling to investigate transmission of bacteria between life stages of the wood-boring beetle, Anoplophora glabripennis. Microb. Ecol. 2009, 58, 199–211. [Google Scholar] [CrossRef] [PubMed]
  13. Xu, L.; Lu, M.; Xu, D.; Chen, L.; Sun, J. Sexual variation of bacterial microbiota of Dendroctonus valens guts and frass in relation to verbenone production. J. Insect Physiol. 2016, 95, 110–117. [Google Scholar] [CrossRef] [PubMed]
  14. Hu, X.; Li, M.; Chen, H. Community structure of gut fungi during different developmental stages of the Chinese white pine beetle (Dendroctonus armandi). Sci. Rep. 2015, 5, 8411. [Google Scholar] [CrossRef] [PubMed]
  15. Oyedokun, A.V.; Adeniyi, D.O. Microbial diversity in the gut of cashew stem girdler, Analeptes trifasciata Fabricius (Coleoptera: Cerambycidae), in Ibadan, Nigeria. Int. J. Insect Sci. 2016, 8, 17–22. [Google Scholar] [CrossRef]
  16. Xu, L.-T.; Lu, M.; Sun, J.-H. Invasive bark beetle-associated microbes degrade a host defensive monoterpene. Insect Sci. 2016, 23, 183–190. [Google Scholar] [CrossRef]
  17. Xu, L.; Deng, J.; Zhou, F.; Cheng, C.; Zhang, L.; Zhang, J.; Lu, M. Gut microbiota in an invasive bark beetle infected by a pathogenic fungus accelerates beetle mortality. J. Pest Sci. 2018, 92, 343–351. [Google Scholar] [CrossRef]
  18. Mason, C.J.; Campbell, A.M.; Scully, E.D.; Hoover, K. Bacterial and fungal midgut community dynamics and transfer between mother and brood in the Asian longhorned beetle (Anoplophora glabripennis), an invasive xylophage. Microb. Ecol. 2019, 77, 230–242. [Google Scholar] [CrossRef]
  19. Soto-Robles, L.V.; Torres-Banda, V.; Rivera-Orduna, F.N.; Curiel-Quesada, E.; Hidalgo-Lara, M.E.; Zuniga, G. An overview of genes from Cyberlindnera americana, a symbiont yeast isolated from the gut of the bark beetle Dendroctonus rhizophagus (Curculionidae: Scolytinae), involved in the detoxification process using genome and transcriptome data. Front. Microbiol. 2019, 10, 2180. [Google Scholar] [CrossRef] [Green Version]
  20. Reid, N.M.; Addison, S.L.; Macdonald, L.J.; Lloyd-Jones, G. Biodiversity of active and inactive bacteria in the gut flora of wood-feeding huhu beetle larvae (Prionoplus reticularis). Appl. Environ. Microbiol. 2011, 77, 7000–7006. [Google Scholar] [CrossRef] [Green Version]
  21. 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, 21006–21020. [Google Scholar] [CrossRef]
  22. Shukla, S.P.; Sanders, J.G.; Byrne, M.J.; Pierce, N.E. Gut microbiota of dung beetles correspond to dietary specializations of adults and larvae. Mol. Ecol. 2016, 25, 6092–6106. [Google Scholar] [CrossRef] [PubMed]
  23. Schwab, D.B.; Riggs, H.E.; Newton, I.L.; Moczek, A.P. Developmental and ecological benefits of the maternally transmitted microbiota in a dung beetle. Am. Nat. 2016, 188, 679–692. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Parker, E.S.; Dury, G.J.; Moczek, A.P. Transgenerational developmental effects of species-specific, maternally transmitted microbiota in Onthophagus dung beetles. Ecol. Entomol. 2019, 44, 274–282. [Google Scholar] [CrossRef]
  25. Shukla, S.P.; Vogel, H.; Heckel, D.G.; Vilcinskas, A.; Kaltenpoth, M. Burying beetles regulate the microbiome of carcasses and use it to transmit a core microbiota to their offspring. Mol. Ecol. 2018, 27, 1980–1991. [Google Scholar] [CrossRef] [PubMed]
  26. Vogel, H.; Shukla, S.P.; Engl, T.; Weiss, B.; Fischer, R.; Steiger, S.; Heckel, D.G.; Kaltenpoth, M.; Vilcinskas, A. The digestive and defensive basis of carcass utilization by the burying beetle and its microbiota. Nat. Commun. 2017, 8, 15186. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Heise, P.; Liu, Y.; Degenkolb, T.; Vogel, H.; Schaberle, T.F.; Vilcinskas, A. Antibiotic-producing beneficial bacteria in the gut of the burying beetle Nicrophorus vespilloides. Front. Microbiol. 2019, 10, 1178. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Ulyshen, M.D.; Zachos, L.G.; Stireman, J.O.; Sheehan, T.N.; Garrick, R.C.; Basset, Y.; Keyghobadi, N. Insights into the ecology, genetics and distribution of Lucanus elaphus Fabricius (Coleoptera: Lucanidae), North America’s giant stag beetle. Insect Conserv. Divers. 2017, 10, 331–340. [Google Scholar] [CrossRef] [Green Version]
  29. Tanahashi, M.; Kubota, K.; Matsushita, N.; Togashi, K. Discovery of mycangia and the associated xylose-fermenting yeasts in stag beetles (Coleoptera: Lucanidae). Naturwissenschaften 2010, 97, 311–331. [Google Scholar] [CrossRef]
  30. Tanahashi, M.; Kubota, K. Utilization of the nutrients in the soluble and insoluble fractions of fungal mycelium by larvae of the stag beetle, Dorcus rectus (Coleoptera: Lucanidae). Eur. J. Entomol. 2013, 110, 611–615. [Google Scholar] [CrossRef] [Green Version]
  31. Miyashita, A.; Hirai, Y.; Sekimizu, K.; Kaito, C. Antibiotic-producing bacteria from stag beetle mycangia. Drug Discov. Ther. 2015, 9, 33–37. [Google Scholar] [CrossRef] [Green Version]
  32. Tanahashi, M.; Hawes, C.J. The presence of a mycangium in European Sinodendron cylindricum (Coleoptera: Lucanidae) and the associated yeast symbionts. J. Insect Sci. 2016, 16, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Tanahashi, M.; Ikeda, H.; Kubota, K. Elementary budget of stag beetle larvae associated with selective utilization of nitrogen in decaying wood. Naturwissenschaften 2018, 105, 33. [Google Scholar] [CrossRef] [PubMed]
  34. Jiang, Y.; Sun, B.; Cao, Y.; Zhai, Y.; Wan, X. Diversity of gut bacterial communities in male adults of Odontolabis fallaciosa (Coleoptera: Scarabaeoidea: Lucanidae) with different mandibular forms. Acta Entomol. Sin. 2018, 61, 322–330. [Google Scholar] [CrossRef]
  35. Wang, J.; Chen, H.; Tang, M. Community structure of gut bacteria of Dendroctonus armandi (Coleoptera: Curculionidae: Scolytinae) larvae during overwintering stage. Sci. Rep. 2017, 7, 14242. [Google Scholar] [CrossRef] [Green Version]
  36. Vasanthakumar, A.; Delalibera, I.; Handelsman, J.; Klepzig, K.D.; Schloss, P.D.; Raffa, K.F. Characterization of gut-associated bacteria in larvae and adults of the southern pine beetle, Dendroctonus frontalis Zimmermann. Environ. Entomol. 2006, 35, 1710–1717. [Google Scholar] [CrossRef]
  37. Morales-Jimenez, J.; Zuniga, G.; Ramirez-Saad, H.C.; Hernandez-Rodriguez, C. Gut-associated bacteria throughout the life cycle of the bark beetle Dendroctonus rhizophagus Thomas and Bright (Curculionidae: Scolytinae) and their cellulolytic activities. Microb. Ecol. 2012, 64, 268–278. [Google Scholar] [CrossRef]
  38. Yusoff, M.Z.; Hu, A.; Feng, C.; Maeda, T.; Shirai, Y.; Hassan, M.A.; Yu, C.P. Influence of pretreated activated sludge for electricity generation in microbial fuel cell application. Bioresour. Technol. 2013, 145, 90–96. [Google Scholar] [CrossRef]
  39. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Pena, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods. 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  40. Wang, Q.; Garrity, G.M.; Tiedje, J.M.; Cole, J.R. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 2007, 73, 5261–5267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Xiang, X.; He, D.; He, J.-S.; Myrold, D.D.; Chu, H. Ammonia-oxidizing bacteria rather than archaea respond to short-term urea amendment in an alpine grassland. Soil Biol. Biochem. 2017, 107, 218–225. [Google Scholar] [CrossRef]
  42. Xiang, X.; Zhang, F.; Fu, R.; Yan, S.; Zhou, L. Significant differences in bacterial and potentially pathogenic communities between sympatric hooded crane and greater white-fronted goose. Front. Microbiol. 2019, 10, 163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package. R Package Version 2.0-10. 2016. Available online: https://www.researchgate.net/publication/258996451 (accessed on 1 October 2020).
  44. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2012; Available online: http://www.gnu.org/copyleft/gpl.html (accessed on 1 October 2020).
  45. Xiang, X.; Gibbons, S.M.; Li, H.; Shen, H.; Fang, J.; Chu, H. Shrub encroachment is associated with changes in soil bacterial community composition in a temperate grassland ecosystem. Plant Soil 2018, 425, 539–551. [Google Scholar] [CrossRef]
  46. Xiang, X.; Liu, J.; Zhang, J.; Li, D.; Xu, C.; Kuzyakov, Y. Divergence in fungal abundance and community structure between soils under long-term mineral and organic fertilization. Soil Tillage Res. 2020, 196. [Google Scholar] [CrossRef]
  47. Huang, S.; Sheng, P.; Zhang, H. Isolation and identification of cellulolytic bacteria from the gut of Holotrichia parallela larvae (Coleoptera: Scarabaeidae). Int. J. Mol. Sci. 2012, 13, 2563–2577. [Google Scholar] [CrossRef]
  48. Berasategui, A.; Axelsson, K.; Nordlander, G.; Schmidt, A.; Borg-Karlson, A.K.; Gershenzon, J.; Terenius, O.; Kaltenpoth, M. The gut microbiota of the pine weevil is similar across Europe and resembles that of other conifer-feeding beetles. Mol. Ecol. 2016, 25, 4014–4031. [Google Scholar] [CrossRef]
  49. Flint, H.J.; Scott, K.P.; Duncan, S.H.; Louis, P.; Forano, E. Microbial degradation of complex carbohydrates in the gut. Gut Microbes 2012, 3, 289–306. [Google Scholar] [CrossRef] [Green Version]
  50. Yun, J.H.; Roh, S.W.; Whon, T.W.; Jung, M.J.; Kim, M.S.; Park, D.S.; Yoon, C.; Nam, Y.D.; Kim, Y.J.; Choi, J.H.; et al. Insect gut bacterial diversity determined by environmental habitat, diet, developmental stage, and phylogeny of host. Appl. Environ. Microbiol. 2014, 80, 5254–5264. [Google Scholar] [CrossRef] [Green Version]
  51. Arias-Cordero, E.; Ping, L.; Reichwald, K.; Delb, H.; Platzer, M.; Boland, W. Comparative evaluation of the gut microbiota associated with the below- and above-ground life stages (larvae and beetles) of the forest cockchafer, Melolontha hippocastani. PLoS ONE 2012, 7, e51557. [Google Scholar] [CrossRef] [Green Version]
  52. Vacchini, V.; Gonella, E.; Crotti, E.; Prosdocimi, E.M.; Mazzetto, F.; Chouaia, B.; Callegari, M.; Mapelli, F.; Mandrioli, M.; Alma, A.; et al. Bacterial diversity shift determined by different diets in the gut of the spotted wing fly Drosophila suzukiiis primarily reflected on acetic acid bacteria. Environ. Microbiol. Rep. 2017, 9, 91–103. [Google Scholar] [CrossRef] [Green Version]
  53. Colman, D.R.; Toolson, E.C.; Takacs-Vesbach, C.D. Do diet and taxonomy influence insect gut bacterial communities? Mol. Ecol. 2012, 21, 5124–5137. [Google Scholar] [CrossRef]
  54. Chen, B.; Teh, B.S.; Sun, C.; Hu, S.; Lu, X.; Boland, W.; Shao, Y. Biodiversity and activity of the gut microbiota across the life history of the insect herbivore Spodoptera littoralis. Sci. Rep. 2016, 6, 29505. [Google Scholar] [CrossRef] [PubMed]
  55. Hongoh, Y.; Ekpornprasit, L.; Inoue, T.; Moriya, S.; Trakulnaleamsai, S.; Ohkuma, M.; Noparatnaraporn, N.; Kudo, T. Intracolony variation of bacterial gut microbiota among castes and ages in the fungus-growing termite Macrotermes gilvus. Mol. Ecol. 2006, 15, 505–516. [Google Scholar] [CrossRef] [PubMed]
  56. Ziganshina, E.E.; Mohammed, W.S.; Shagimardanova, E.I.; Vankov, P.Y.; Gogoleva, N.E.; Ziganshin, A.M. Fungal, bacterial, and archaeal diversity in the digestive tract of several beetle larvae (Coleoptera). BioMed Res. Int. 2018, 2018, 6765438. [Google Scholar] [CrossRef] [PubMed]
  57. Biddle, J.F.; Fitz-Gibbon, S.; Schuster, S.C.; Brenchley, J.E.; House, C.H. Metagenomic signatures of the Peru Margin subseafloor biosphere show a genetically distinct environment. Proc. Natl. Acad. Sci. USA 2008, 105, 10583–10588. [Google Scholar] [CrossRef] [Green Version]
  58. Waters, J.L.; Ley, R.E. The human gut bacteria Christensenellaceae are widespread, heritable, and associated with health. BMC Biol. 2019, 17, 83. [Google Scholar] [CrossRef]
  59. Marchandin, H.; Jumas-Bilak, E. The Family Veillonellaceae. In The Prokaryotes: Firmicutes and Tenericutes; Springer: Berlin/Heidelberg, Germany, 2014; pp. 433–453. [Google Scholar] [CrossRef]
  60. Gutierrez, C.; Somoskovi, A. Human Pathogenic Mycobacteria. In Reference Module in Biomedical Sciences; Elsevier: Amsterdam, The Netherlands, 2014. [Google Scholar] [CrossRef]
  61. Scher, J.U.; Sczesnak, A.; Longman, R.S.; Segata, N.; Ubeda, C.; Bielski, C.; Rostron, T.; Cerundolo, V.; Pamer, E.G.; Abramson, S.B.; et al. Expansion of intestinal Prevotella copri correlates with enhanced susceptibility to arthritis. eLife 2013, 2, e01202. [Google Scholar] [CrossRef]
  62. Shil, R.K.; Mojumder, S.; Sadida, F.F.; Uddin, M.; Sikdar, D. Isolation and identification of cellulolytic bacteria from the gut of three phytophagus insect species. Braz. Arch. Biol. Technol. 2014, 57, 927–932. [Google Scholar] [CrossRef]
  63. Sun, X.; Yang, Y.; Zhang, N.; Shen, Y.; Ni, J. Draft genome sequence of Dysgonomonas macrotermitis strain JCM 19375T, isolated from the gut of a termite. Genome Announc. 2015, 3, e00963-15. [Google Scholar] [CrossRef] [Green Version]
  64. Berasategui, A.; Salem, H.; Paetz, C.; Santoro, M.; Gershenzon, J.; Kaltenpoth, M.; Schmidt, A. Gut microbiota of the pine weevil degrades conifer diterpenes and increases insect fitness. Mol. Ecol. 2017, 26, 4099–4110. [Google Scholar] [CrossRef]
  65. Bellinvia, S.; Johnston, P.R.; Reinhardt, K.; Otti, O. Bacterial communities of the reproductive organs of virgin and mated common bedbugs Cimex lectularius. Ecol. Entomol. 2019, 45, 142–154. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Venn diagram showing the co-occurence of the OTUs among samples from different stages. Numbers in parentheses indicate total OTUs in each stage group, and numbers inside the Venn diagram indicate unique and shared OTUs. OTU, operational taxonomic unit.
Figure 1. Venn diagram showing the co-occurence of the OTUs among samples from different stages. Numbers in parentheses indicate total OTUs in each stage group, and numbers inside the Venn diagram indicate unique and shared OTUs. OTU, operational taxonomic unit.
Insects 11 00719 g001
Figure 2. Intestinal bacterial alpha-diversity in different developmental stages. The bottom and top of the box denote the first and third quartiles respectively; the band inside the box denotes median; error bars denote standard deviations; different letters above bars represent significant differences from Tukey’s HSD comparisons (p < 0.05). L1, 1st instar; L2, 2nd instar; L3, 3rd instar; Ad, Adults.
Figure 2. Intestinal bacterial alpha-diversity in different developmental stages. The bottom and top of the box denote the first and third quartiles respectively; the band inside the box denotes median; error bars denote standard deviations; different letters above bars represent significant differences from Tukey’s HSD comparisons (p < 0.05). L1, 1st instar; L2, 2nd instar; L3, 3rd instar; Ad, Adults.
Insects 11 00719 g002
Figure 3. LEfSe analysis of intestinal bacterial biomarkers associated with host types. Identified phylotype biomarkers ranked by effect size, and the alpha value was < 0.05. Cladogram representing the taxonomic hierarchical structure of the phylotype biomarkers identified among four host types, red, phylotypes overrepresented in gut of L1; green, phylotypes statistically overrepresented in gut of L2; blue, phylotypes statistically overrepresented in gut of L3; purple, phylotypes statistically overrepresented in gut of Ad. L1, 1st instar; L2, 2nd instar; L3, 3rd instar; Ad, Adults.
Figure 3. LEfSe analysis of intestinal bacterial biomarkers associated with host types. Identified phylotype biomarkers ranked by effect size, and the alpha value was < 0.05. Cladogram representing the taxonomic hierarchical structure of the phylotype biomarkers identified among four host types, red, phylotypes overrepresented in gut of L1; green, phylotypes statistically overrepresented in gut of L2; blue, phylotypes statistically overrepresented in gut of L3; purple, phylotypes statistically overrepresented in gut of Ad. L1, 1st instar; L2, 2nd instar; L3, 3rd instar; Ad, Adults.
Insects 11 00719 g003
Figure 4. Non-metric multidimensional scaling (NMDS) plot showing bacterial community composition between larvae and adults (A), and contrasts between different larval instars (B). L1, 1st instar; L2, 2nd instar; L3, 3rd instar; Ad, Adults.
Figure 4. Non-metric multidimensional scaling (NMDS) plot showing bacterial community composition between larvae and adults (A), and contrasts between different larval instars (B). L1, 1st instar; L2, 2nd instar; L3, 3rd instar; Ad, Adults.
Insects 11 00719 g004
Figure 5. Relative abundance of bacterial taxa at phylum in the gut of different life stages of P. muelleri. Bars represent mean; error bars denote standard deviation, letters above bars represents significant differences from one-way ANOVA (p < 0.05). L1, 1st instar; L2, 2nd instar; L3, 3rd instar; Ad, Adults.
Figure 5. Relative abundance of bacterial taxa at phylum in the gut of different life stages of P. muelleri. Bars represent mean; error bars denote standard deviation, letters above bars represents significant differences from one-way ANOVA (p < 0.05). L1, 1st instar; L2, 2nd instar; L3, 3rd instar; Ad, Adults.
Insects 11 00719 g005
Table 1. Differences in bacteria community composition depending on different development stages and sexes by analyses of similarities (ANOSIM).
Table 1. Differences in bacteria community composition depending on different development stages and sexes by analyses of similarities (ANOSIM).
StagesANOSIM
Rp
L1 vs. L20.3020.007
L1 vs. L30.4970.001
L1 vs. Ad0.9990.001
L2 vs. L30.1090.059
L2 vs. Ad1.0000.001
L3 vs. Ad1.0000.001
L3♀ VS L3♂0.2000.047
Ad♀ VS Ad♂−0.0400.578
L1, 1st instar; L2, 2nd instar; L3, 3rd instar; Ad, Adults. ♀, female; ♂ male.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, M.; Xiang, X.; Wan, X. Divergence in Gut Bacterial Community Among Life Stages of the Rainbow Stag Beetle Phalacrognathus muelleri (Coleptera: Lucanidae). Insects 2020, 11, 719. https://doi.org/10.3390/insects11100719

AMA Style

Wang M, Xiang X, Wan X. Divergence in Gut Bacterial Community Among Life Stages of the Rainbow Stag Beetle Phalacrognathus muelleri (Coleptera: Lucanidae). Insects. 2020; 11(10):719. https://doi.org/10.3390/insects11100719

Chicago/Turabian Style

Wang, Miaomiao, Xingjia Xiang, and Xia Wan. 2020. "Divergence in Gut Bacterial Community Among Life Stages of the Rainbow Stag Beetle Phalacrognathus muelleri (Coleptera: Lucanidae)" Insects 11, no. 10: 719. https://doi.org/10.3390/insects11100719

APA Style

Wang, M., Xiang, X., & Wan, X. (2020). Divergence in Gut Bacterial Community Among Life Stages of the Rainbow Stag Beetle Phalacrognathus muelleri (Coleptera: Lucanidae). Insects, 11(10), 719. https://doi.org/10.3390/insects11100719

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop