<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2023.1199187</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Influence of native and exotic plant diet on the gut microbiome of the Gray&#x00027;s Malayan stick insect, <italic>Lonchodes brevipes</italic></article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Lim</surname> <given-names>Yan Zhen</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2324245/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Poh</surname> <given-names>Yan Hong</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Lee</surname> <given-names>Kevin C.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/377762/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Pointing</surname> <given-names>Stephen Brian</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/202495/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Wainwright</surname> <given-names>Benjamin J.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1037860/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Tan</surname> <given-names>Eunice Jingmei</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/459398/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Division of Science, Yale-NUS College</institution>, <addr-line>Singapore</addr-line>, <country>Singapore</country></aff>
<aff id="aff2"><sup>2</sup><institution>School of Science, Auckland University of Technology</institution>, <addr-line>Auckland</addr-line>, <country>New Zealand</country></aff>
<aff id="aff3"><sup>3</sup><institution>Department of Biological Sciences, National University of Singapore</institution>, <addr-line>Singapore</addr-line>, <country>Singapore</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Amit Kumar, Nanjing University of Information Science and Technology, China</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Irene Stefanini, University of Turin, Italy; Kara Tinker, National Energy Technology Laboratory (DOE), United States</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Eunice Jingmei Tan <email>ynctje&#x00040;nus.edu.sg</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>27</day>
<month>07</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1199187</elocation-id>
<history>
<date date-type="received">
<day>03</day>
<month>04</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>05</day>
<month>07</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2023 Lim, Poh, Lee, Pointing, Wainwright and Tan.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Lim, Poh, Lee, Pointing, Wainwright and Tan</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license> </permissions>
<abstract>
<p>Herbivorous insects require an active lignocellulolytic microbiome to process their diet. Stick insects (phasmids) are common in the tropics and display a cosmopolitan host plant feeding preference. The microbiomes of social insects are vertically transmitted to offspring, while for solitary species, such as phasmids, it has been assumed that microbiomes are acquired from their diet. This study reports the characterization of the gut microbiome for the Gray&#x00027;s Malayan stick insect, <italic>Lonchodes brevipes</italic>, reared on native and introduced species of host plants and compared to the microbiome of the host plant and surrounding soil to gain insight into possible sources of recruitment. Clear differences in the gut microbiome occurred between insects fed on native and exotic plant diets, and the native diet displayed a more species-rich fungal microbiome. While the findings suggest that phasmids may be capable of adapting their gut microbiome to changing diets, it is uncertain whether this may lead to any change in dietary efficiency or organismal fitness. Further insight in this regard may assist conservation and management decision-making.</p></abstract>
<kwd-group>
<kwd>insect microbiomes</kwd>
<kwd>phasmatodea</kwd>
<kwd><italic>Lonchodes brevipes</italic></kwd>
<kwd>host plant</kwd>
<kwd>herbivores</kwd>
</kwd-group>
<counts>
<fig-count count="5"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="49"/>
<page-count count="9"/>
<word-count count="5631"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Microbial Symbioses</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1. Introduction</title>
<p>Herbivorous insects can host microorganisms, such as bacteria, viruses, fungi, archaea, and protozoa in their digestive tract, referred to as their gut microbiome (Gurung et al., <xref ref-type="bibr" rid="B14">2019</xref>; Mu&#x000F1;oz-Benavent et al., <xref ref-type="bibr" rid="B32">2021</xref>). The gut microbiome of insects has been implicated in the transformation of food substrates and acting as a direct dietary source providing essential nutrients lacking in the plant diet (Dillon and Dillon, <xref ref-type="bibr" rid="B11">2003</xref>; Moran et al., <xref ref-type="bibr" rid="B31">2019</xref>). Evidence for an adapted lignocellulosic microbiome is strong for some insect groups, notably termites (K&#x000F6;nig et al., <xref ref-type="bibr" rid="B23">2013</xref>; Jang and Kikuchi, <xref ref-type="bibr" rid="B19">2020</xref>). Studies have shown that bacteria and fungi in insect guts possess enzymes breaking down insoluble plant polymers, including cellulose and hemicellulose, in addition to soluble compounds such as starch and pectin (Watanabe and Tokuda, <xref ref-type="bibr" rid="B45">2010</xref>; Jang and Kikuchi, <xref ref-type="bibr" rid="B19">2020</xref>). Additional roles for the gut microbiome have been inferred from studies of other animal systems and include protection against pathogen infection (Mendes and Raaijmakers, <xref ref-type="bibr" rid="B30">2015</xref>; Raymann and Moran, <xref ref-type="bibr" rid="B35">2018</xref>; Jang and Kikuchi, <xref ref-type="bibr" rid="B19">2020</xref>) and behavioral influences (Davidson et al., <xref ref-type="bibr" rid="B8">2018</xref>, <xref ref-type="bibr" rid="B9">2020</xref>; Jones et al., <xref ref-type="bibr" rid="B21">2018</xref>).</p>
<p>The recruitment of gut microbiomes in insects is not well constrained, with studies on social insects such as bees and termites suggesting direct transfer between individuals [e.g., (Koch and Schmid-Hempel, <xref ref-type="bibr" rid="B22">2011</xref>; Bourguignon et al., <xref ref-type="bibr" rid="B4">2018</xref>)], while other studies indicate that gut microbiomes may reflect recruitment from dietary substrates (Jones et al., <xref ref-type="bibr" rid="B20">2019</xref>). Studies on lepidopterans, for instance, demonstrated a correlation between diet and gut microbial diversity, where the host plants were the main contributors to the herbivorous insects&#x00027; gut microbiome (Mason et al., <xref ref-type="bibr" rid="B27">2020</xref>). The extent of the effect of diet on gut microbiome differed depending on both the lepidopteran host and host plant species (Jones et al., <xref ref-type="bibr" rid="B20">2019</xref>) and their exposure to environments with soil (Hannula et al., <xref ref-type="bibr" rid="B15">2019</xref>). In herbivorous insects, the gut microbiome has been shown to largely reflect the host insect phylogeny rather than their ecology (McLean et al., <xref ref-type="bibr" rid="B28">2019</xref>). The extent to which the insect gut acts as an environmental filter that shapes the diversity of its microbiome toward taxa that are beneficial symbionts for the insect host appears to depend upon the insect taxa (Moran et al., <xref ref-type="bibr" rid="B31">2019</xref>; Mason, <xref ref-type="bibr" rid="B25">2020</xref>).</p>
<p>Phasmids (stick insects) are non-social folivorous insects that feed almost exclusively on leaves but are able to adapt to foreign diets, that is, plant species that the insects are not usually exposed to or would consume, if the favored host is unavailable (Boucher and Varady-Szabo, <xref ref-type="bibr" rid="B3">2005</xref>; Bl&#x000FC;thgen et al., <xref ref-type="bibr" rid="B2">2006</xref>; Shelomi et al., <xref ref-type="bibr" rid="B37">2014a</xref>). The current understanding of gut microbiomes in phasmids is not well constrained. In one study, bacterial gut microbiomes showed no significant differences in composition across several phasmid families (Shelomi et al., <xref ref-type="bibr" rid="B37">2014a</xref>). By sequencing the gut contents, one study found that Proteobacteria was dominant, followed by Actinobacteria and Firmicutes, and <italic>Spiroplasma</italic> spp. and <italic>Sphingobium</italic> spp. were the most common microbes isolated (Shelomi et al., <xref ref-type="bibr" rid="B38">2013</xref>). In another study, cultured bacterial strains from phasmid guts affiliated with the Actinomycetales, Enterobacteriales, Lactobacillales, Pseudomonadales, and Xanthomonadales (Shelomi et al., <xref ref-type="bibr" rid="B39">2015</xref>). It has been suggested that the environment and diet can affect the gut microbiome of phasmids (Ignasiak and Maxwell, <xref ref-type="bibr" rid="B18">2017</xref>). Furthermore, gut microorganisms found in phasmids were shown to produce cell wall degrading enzymes including cellulases important to the digestion of plant matter (Shelomi et al., <xref ref-type="bibr" rid="B40">2014b</xref>), and it has also been postulated that some phasmids possess their own genes for breakdown of soluble plant carbohydrates, for example, pectin (Shelomi et al., <xref ref-type="bibr" rid="B36">2016</xref>). For species of phasmids that do not adhere their eggs to the leaf substrate, this stage in life history is unlikely to be a route for microbiome recruitment from the plant substrate.</p>
<p>Here, we set out to further interrogate the hypothesis that different plant diets in phasmids would lead to shifts in the taxonomic composition of the gut microbiome. This is of relevance to Singapore because of the high levels of urbanization with a mixture of native and exotic plant host species and the potential for a shift in plant hosts of herbivores. Specifically, we sought to test the following hypotheses: (a) Does phasmid gut microbiota vary with diet, and (b) is phasmid gut microbiota a subset of food (i.e., plant) or local soil microbiomes? We reared a local species, Gray&#x00027;s Malayan stick insect, <italic>Lonchodes brevipes</italic>, on two species of host plants. This phasmid is encountered in mangroves and rainforests (Seow-Choen, pers. obs.) and feeds on sea hibiscus (SH), <italic>Hibiscus tiliaceus</italic> in the wild, and does not adhere their eggs to the leaf substrate. In hobbyist cultures, this phasmid is commonly reared on a widespread non-native ornamental plant, golden penda (GP), <italic>Xanthostemon chrysanthus</italic>. We compared the gut microbiome of individuals reared on both plants, as well as the microbiome of the host plant leaves and soil because soil is the most microbially emissive component of the habitat.</p>
</sec>
<sec id="s2">
<title>2. Materials and methods</title>
<sec>
<title>2.1. Experimental subjects and environmental samples</title>
<p>We used laboratory-reared offspring from field-caught <italic>Lonchodes brevipes</italic> in this study. This ensured that no adaptations to domestication were introduced. Phasmids were housed individually in a temperature-controlled laboratory at 24&#x000B0;C in a 12-h:12-h light: dark cycle. All individuals were provided food plants <italic>ad libitum</italic> and monitored to be feeding and behaving normally prior to the trials. Individuals with missing or misshapen legs or body parts that may have affected their feeding behavior were not used in our experiments. We used a total of 20 nymphs between 12 and 16 weeks of age; 10 nymphs were raised since hatching on the leaves of the golden penda while another 10 nymphs were raised on sea hibiscus. From initial tests, frass collected from single days yielded insufficient DNA; therefore, frass was collected from individual nymphs daily over 7-day intervals and pooled for DNA extraction. In total, 10 samples of each of the host plant leaves and topsoil within a meter around the tree trunks were collected aseptically into sterile containers. As the vegetation on campus is regularly sprayed with or exposed to pesticides, all leaves were washed prior to feeding the phasmids. Therefore, for consistency, leaves used for sequencing were also washed prior to DNA extraction.</p>
</sec>
<sec>
<title>2.2. DNA extraction, PCR amplification, and sequencing</title>
<p>First, leaf and frass samples were flash frozen in liquid nitrogen and ground with the mortar and pestle. We extracted the DNA from leaf samples using FavorPrep<sup>TM</sup> Plant Genomic DNA Extraction Mini Kit (Favorgen), while DNA extraction from frass and soil samples was performed using QIAamp PowerFecal Pro DNA Kit (QIAGEN). To maximize yield, phenol:chloroform:isoamyl alcohol was added to the frass prior to homogenization. Extracted DNA samples were then cleaned with the DNeasy PowerClean Pro Cleanup Kit (QIAGEN).</p>
<p>PCR amplification of the V4 region of the 16S SSU rRNA gene was performed using the bacterial and archaeal primers 515F and 806R [515F&#x02014;GTG CCA GCM GCC GCG GTA A; 806R&#x02014;GGA CTA CHV GGG TWT CTA AT; Caporaso (<xref ref-type="bibr" rid="B6">2011</xref>)]. Forward and reverse primers were modified to include Illumina adapters, a linker, and a unique barcode (Caporaso, <xref ref-type="bibr" rid="B6">2011</xref>). Each reaction was performed in a total volume of 25 &#x003BC;l, containing 0.75 &#x003BC;l of each primer at 10 &#x003BC;M, 12.5 &#x003BC;l of buffer, 1.5 &#x003BC;l BSA at 10 mg/ml, 0.5 &#x003BC;l of MgCl<sub>2</sub>, 0.1 &#x003BC;l of enzyme (KAPA 3G, Kapa Biosystems, Inc, Wilmington, MA, USA), 1 &#x003BC;l of undiluted DNA, with H<sub>2</sub>O added to make up to 25 &#x003BC;l. PCR cycling condition was set at 94&#x000B0;C for 180 s, followed by 35 cycles of 94&#x000B0;C for 45 s, 50&#x000B0;C for 60 s, and 72&#x000B0;C for 90 s, and a final extension at 72&#x000B0;C for 10 min. Negative PCR controls were included to identify any potential laboratory reagent contamination.</p>
<p>Fungal DNA amplification was performed with the ITS1F [5&#x02032;-CTT GGT CAT TTA GAG GAA GTA A-3&#x02032;; Gardes and Bruns (<xref ref-type="bibr" rid="B13">1993</xref>)] and ITS2 [5&#x02032;-GCT GCG TTC TTC ATC GAT GC-3&#x02032;; White et al. (<xref ref-type="bibr" rid="B47">1990</xref>)] primer set. Primers were modified to include a unique barcode, linker, and Illumina adapter [see Smith and Peay (<xref ref-type="bibr" rid="B41">2014</xref>) for exact details]. Each reaction was performed in a total volume of 25 &#x003BC;l, containing 0.75 &#x003BC;l of each primer at 10 &#x003BC;M, 12.5 &#x003BC;l of buffer, 1.5 &#x003BC;l BSA at 10 mg/ml, 0.5 &#x003BC;l of MgCl<sub>2</sub>, 0.1 &#x003BC;l of enzyme (KAPA 3G, Kapa Biosystems, Inc, Wilmington, MA, USA), 1 &#x003BC;l of undiluted DNA, with H<sub>2</sub>O added to make up to 25 &#x003BC;l. PCR cycling condition was set at 95&#x000B0;C for 3 min, followed by 35 cycles of 95&#x000B0;C for 20 s, 53&#x000B0;C for 15 s, and 72&#x000B0;C for 20 s with a final extension at 72&#x000B0;C for 60 s. Negative PCR controls substituting water for DNA were included to aid in the detection of contamination.</p>
<p>Successful PCR amplification was confirmed on a 1% TAE buffer agarose gel. normalization and cleaning of PCR products were performed in sequalprep normalization plates (Invitrogen, Frederick, 129 Maryland, USA). Both libraries were sequenced independently on the Illumina MiSeq platform (600 130 cycles, V3 chemistry, 300-bp paired end reads) with a 30% PhiX spike (Macrogen, Korea).</p>
</sec>
<sec>
<title>2.3. Data analysis</title>
<p>The amplicon sequencing data were processed using the R package dada2 v.1.18.0 (Callahan et al., <xref ref-type="bibr" rid="B5">2016</xref>) to infer distinct amplicon sequence variants (ASVs). The workflow was as follows: primer sequences were first removed via cutadapt v3.3 (Martin, <xref ref-type="bibr" rid="B24">2011</xref>). The bacterial reads were uniformly trimmed to the length of 230/220 bp (forward/reverse) and filtered to remove reads exceeding 2/5 (forward/reverse) expected errors or reads shorter than 175 bp. The fungal reads were not trimmed, but reads exceeding 5/8 bp (forward/reverse) expected errors or reads shorter than 50 bp were removed. In both cases, any reads with N or those identified as possible PhiX sequences were removed. The ASVs were inferred with pseudo-pooling to increase sensitivity and preserve singletons for richness estimations. Potential chimeric sequences were removed, and the inferred ASVs were classified using the naive Bayesian classifier implemented in dada2 with the reference database UNITE v8.2 eukaryote (Nilsson et al., <xref ref-type="bibr" rid="B33">2018</xref>) for fungi and SILVA v138 (Quast et al., <xref ref-type="bibr" rid="B34">2013</xref>) for bacteria. The R package phyloseq v1.36.0 (McMurdie and Holmes, <xref ref-type="bibr" rid="B29">2013</xref>) was used to handle the resulting ASV abundance table for ecological analyses and visualizations. The R package decontam v1.12.0 (Davis et al., <xref ref-type="bibr" rid="B10">2018</xref>) was used to identify putative contaminant ASVs by removing ASVs identified by either the frequency (threshold = 0.1) or the prevalence (threshold = 0.5) method, based on the taxa in the control samples and the DNA concentrations of the initial extractants (Nanodrop, Thermo Fisher Scientific, USA). After removing ASVs associated with non-bacterial or non-fungal (genomic) sources (e.g., mitochondria and chloroplast), samples with &#x0003C; 1,000 reads were removed due to insufficient information and potential bias that they may introduce. Excluding control samples, the processing and filtering resulted in 69 fungi samples with a total of 3,417,410 counts over 3,452 ASVs, as well as 72 bacterial samples with a total of 2,683,111 counts over 8,314 ASVs. The taxonomy and taxa named are based on the SILVA reference database that we used.</p>
<p>General processing of the community data including the calculation of relative abundance and estimates of alpha diversity was conducted using the R package phyloseq (McMurdie and Holmes, <xref ref-type="bibr" rid="B29">2013</xref>) and visualized using ggplot2 (Wickham, <xref ref-type="bibr" rid="B48">2016</xref>). The proportional differences of bacterial and fungal families in the samples belonging to the two categories (host plant) were compared using the R package ANCOM-BC v.1.2.2; the differences (enriched in one group vs. another) were expressed as the <italic>W</italic> statistic, which reflects the effect size. A larger positive <italic>W</italic> indicates enrichment in sea hibiscus, and a more negative W indicates enrichment in golden penda. Significant differences were based on <italic>p</italic>-values adjusted by the Holm&#x02013;Bonferroni method (adj. <italic>p</italic> &#x0003C; 0.05). Therefore, the two extremes (positive or negative <italic>W</italic>) indicated the largest differences between the two groups of frass samples. Significance testing of the community structure in relation to sample type and the associated host plant species was conducted using &#x0201C;manyglm,&#x0201D; with a negative binomial distribution, from the R package mvabund v.4.2.1 (Wang et al., <xref ref-type="bibr" rid="B44">2012</xref>).</p>
</sec>
</sec>
<sec id="s3">
<title>3. Results and discussion</title>
<p>The high-throughput sequencing of bacterial and fungal diversity in phasmid frass was used as a proxy for gut microbiomes. This approach is congruent with that adopted for most other studies of microbiomes in animal and human systems, for example, Huttenhower et al. (<xref ref-type="bibr" rid="B17">2012</xref>). We used this approach to interrogate the central hypothesis that native and introduced plant food sources impacted the composition of phasmid gut microbiomes. Comparison of microbiota associated with food sources and surrounding soils provided triangulation of potential sources and influence on the composition of the gut microorganisms. The DNA yield from samples was used as a general proxy for overall (i.e., bacterial and fungal) microbial biomass, although we recognize that host DNA may also have been a confounding factor. This indicated that while the food substrate and frass had similar levels of biomass (<xref ref-type="fig" rid="F1">Figure 1A</xref>), the soil samples yielded significantly greater microbial biomass per gram of substrate than frass (<italic>p</italic> = 0.001) and leaves (<italic>p</italic> &#x0003C; 0.001). Of note was the markedly higher biomass in soil associated with native vs. exotic plants, and this may indicate a more developed rhizosphere microbiome for the native plants that reflects their natural history.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>The plots are organized into two columns with the bacterial assemblages on the left and the fungal assemblages on the right. <bold>(A)</bold> shows DNA yields from the various substrates (frass, leaf, soil). <bold>(B, C)</bold> show taxonomic richness using the Chao1 richness estimator. <bold>(D, E)</bold> show taxonomic diversity as Shannon&#x00027;s index (H&#x02032;). Each data point in a box plot represents a biological sample. The data points are grouped into sample types (frass, leaf, and soil) and further separated by associated host plant species.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1199187-g0001.tif"/>
</fig>
<p>Contrasting alpha diversity estimates were observed for bacteria and fungi across the sample types. Bacterial richness estimates (Chao1) were significantly higher in the frass of phasmids fed on native (sea hibiscus) vs. exotic (golden penda) leaf substrates (<xref ref-type="fig" rid="F1">Figure 1B</xref>; Mann&#x02013;Whitney <italic>U</italic> = 0, <italic>n</italic><sub>1</sub> = 16, <italic>n</italic><sub>2</sub> = 19, <italic>p</italic> &#x0003C; 0.001). Similarly, fungal richness was significantly higher in the frass of phasmids reared on native vs. exotic leaves (<xref ref-type="fig" rid="F1">Figure 1C</xref>; Mann&#x02013;Whitney <italic>U</italic> = 0, <italic>n</italic><sub>1</sub> = 16, <italic>n</italic><sub>2</sub> = 19, <italic>p</italic> &#x0003C; 0.001). The Shannon&#x00027;s index estimates that consider richness and evenness were not significantly different for bacteria between native and exotic substrates (<xref ref-type="fig" rid="F1">Figure 1D</xref>), whereas for fungi estimates, richness and evenness were significantly higher for the native substrate (<xref ref-type="fig" rid="F1">Figure 1E</xref>; Mann&#x02013;Whitney <italic>U</italic> = 2, <italic>n</italic><sub>1</sub> = 11, <italic>n</italic><sub>2</sub> = 19, <italic>p</italic> &#x0003C; 0.001). The leaves of native and exotic plants displayed similar bacterial richness to that of frass while the soil as expected supported a markedly greater diversity of bacteria (Thompson et al., <xref ref-type="bibr" rid="B42">2017</xref>). In contrast, both leaves and soil supported greater fungal richness than frass, and this reflected the importance of the phyllosphere as a reservoir of fungal diversity (Xu et al., <xref ref-type="bibr" rid="B49">2022</xref>).</p>
<p>Beta diversity patterns measured by Hellinger distances revealed that gut bacterial assemblages were more similar to each other than their associated leaf or soil assemblages regardless of the leaf diet, suggesting that bacteria in the gut had undergone environmental filtering and recruitment (<xref ref-type="fig" rid="F2">Figure 2A</xref>). In contrast, the fungal assemblages in frass from the different diets displayed less overlap with each other but more overlap with the corresponding leaf substrate, thus suggesting recruitment from the leaf substrate may be important for fungi (<xref ref-type="fig" rid="F2">Figure 2B</xref>). For both bacteria and fungi, the frass assemblages displayed no overlap in overall similarity with soil assemblages. This raises the interesting question of what the recruitment sources are for gut bacteria and fungi in phasmids. The leaves and soil were unable to fully explain the observed diversity in frass and so additional sources might include vertical transmission by phasmids, or other aeolian or environmental sources. Interestingly, the results also revealed that bacteria and fungi in soil respond differently to native and exotic plants, with bacteria assemblages remaining essentially unchanged while a significant shift in fungal assemblages occurred in soils planted with native vs. exotic plants (<xref ref-type="fig" rid="F2">Figures 2A</xref>, <xref ref-type="fig" rid="F2">B</xref>). The altered soil mycobiome with exotic plants may thus have undetermined impacts on the functioning of the soil and phyllosphere. This has potential implications for horticultural management and conservation.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>Ordination of similarities between <bold>(A)</bold> bacterial and <bold>(B)</bold> fungal communities via principal coordinate analysis of Hellinger distances of the rarefied ASVs. Each assemblage derived from a sample is shown as a data point, with the sample type indicated by its shape and the associated host plant by its color. The clustering of each group is highlighted by a 95% confidence interval ellipse assuming a multivariate t-distribution.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1199187-g0002.tif"/>
</fig>
<p>Taxonomic diversity estimation revealed that frass bacterial assemblages were enriched in taxa within Alpha- and Gamma-proteobacteria. The most frequently encountered bacterial families were Pseudomonadaceae and Enterobacteriaceae for the native plant-fed insects and Comamonadaceae and Enterobacteriaceae in the exotic plant-fed insects (<xref ref-type="fig" rid="F3">Figure 3A</xref>; <xref ref-type="supplementary-material" rid="SM1">Supplementary Figure 1</xref>). Interestingly, while Enterobacteriaceae is abundant in both groups of frass samples, in SH, the family consisted almost entirely of the genus <italic>Klebsiella</italic>, whereas in GP, <italic>Salmonella</italic> was also present. Patterns in diversity for fungi were less evident (<xref ref-type="fig" rid="F3">Figure 3B</xref>; <xref ref-type="supplementary-material" rid="SM2">Supplementary Figure 2</xref>) due to the high taxonomic diversity and large numbers of unclassified taxa. Many of the fungal genera (e.g., <italic>Aspergillus</italic> and <italic>Penicillium</italic>) include known cellulose-degrading taxa (V&#x000E1;zquez-Montoya et al., <xref ref-type="bibr" rid="B43">2020</xref>), although whether genetic signals represent vegetative cells rather than dormant spores requires these data to be interpreted with some caution. In order to further interrogate the likelihood that gut microbiomes of insects reared on native and exotic plant feed had distinct assembly patterns, we conducted differential abundance analysis to identify bacterial and fungal families that differed significantly (i.e., were enriched) between frass samples from phasmids fed on the two plant hosts (<xref ref-type="fig" rid="F4">Figure 4A</xref>). Notably, the strong differential abundance pattern for <italic>Cladosporium</italic> species (enriched in SH frass and virtually absent in GP frass) may indicate the effect of selection pressure due to the gut environment affected by the different diets (<xref ref-type="fig" rid="F4">Figure 4B</xref>). The core ASVs shared between the sample type (frass, leaf, or soil) and associated host plant (golden penda or sea hibiscus) varied greatly, with the greatest number of core ASVs shared between the leaves of the two host plant species (<xref ref-type="fig" rid="F5">Figure 5</xref>), while the majority of core ASVs were only found in the respective host plant leaves and frass of insects feeding on sea hibiscus, that is, not shared between communities in different sample types (<xref ref-type="fig" rid="F5">Figure 5</xref>).</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Overall community composition of <bold>(A)</bold> bacteria and <bold>(B)</bold> fungi samples at the class level, by sample type and source plant species, and expressed as relative abundances in each of the groups. Classes with mean relative abundance of &#x0003C;0.5% are pooled into the &#x0201C;Others&#x0201D; category.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1199187-g0003.tif"/>
</fig>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Differential abundance analysis of fungi. <bold>(A)</bold> Overall dominant fungal families in frass samples that were over &#x0003E;1% mean relative abundance. Effect sizes where W = &#x0003E;0 indicated the family was overall more abundant in frass from sea hibiscus leaf substrate, whereas none were more abundant in frass from golden penda leaf substrate (i.e., W = &#x0003C; 0). A teal-colored data point indicates the difference in abundance of the family is statistically significant (adj. <italic>p</italic> &#x0003C; 0.05), while a red-colored data point indicates the association was not supported. <bold>(B)</bold> Species within the fungal genus <italic>Cladosporium</italic> were strongly associated with frass produced from sea hibiscus and virtually absent in the golden penda frass.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1199187-g0004.tif"/>
</fig>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>The number of core ASVs shared between groups based on sample type (frass, leaf, soil) and associated host plant (golden penda or sea hibiscus) in the <bold>(A)</bold> bacteria and <bold>(B)</bold> fungi assemblages. A core ASV for a group of samples is here defined as being above 0.1% in relative abundance in over 25% of samples.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1199187-g0005.tif"/>
</fig>
<p>Overall, this study shows that community structure was driven by host selection instead of dispersal from the environment, as the environments (soil/leaf) were not as different as the frass associated with two source plant species. Nonetheless, the sample types (i.e., frass, leaf, and soil), the associated host plant species, and the interaction term of the two factors were all statistically significant in structuring the bacterial and fungal communities observed in this study (manyglm, all <italic>p</italic> &#x0003C; 0.001). Furthermore, the microbiomes varied significantly in response to a native or exotic plant diet, with fungi accounting for more variation than bacteria. While the gut communities of some herbivorous insects are characterized by transient and food-derived bacteria (Whitaker et al., <xref ref-type="bibr" rid="B46">2016</xref>; Hern&#x000E1;ndez-Garc&#x000ED;a et al., <xref ref-type="bibr" rid="B16">2017</xref>), other studies do not find a strong correlation between diet and the gut microbiome (Mason and Raffa, <xref ref-type="bibr" rid="B26">2014</xref>; Chaturvedi et al., <xref ref-type="bibr" rid="B7">2017</xref>). A previous study suggests that phasmids are well adapted to changes in diets and may not be reliant on their gut microbiome for specialized digestion of plants (Shelomi et al., <xref ref-type="bibr" rid="B40">2014b</xref>). Our study provides additional insight that indicates a shift in diet, resulting in measurable changes to the insect gut microbiome that includes known cellulolytic taxa. The data also indicate that exotic plants result in altered soil and herbivorous insect microbiomes and this is relevant to the conservation and management issues where factors that potentially affect insect herbivory are important (Bahrndorff et al., <xref ref-type="bibr" rid="B1">2016</xref>). Finally, the characterization of the <italic>Lonchodes brevipes</italic> microbiome also adds to the growing global inventory of animal microbiomes (Douglas, <xref ref-type="bibr" rid="B12">2019</xref>).</p>
</sec>
<sec sec-type="data-availability" id="s4">
<title>Data availability statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/<xref ref-type="sec" rid="s8">Supplementary material</xref>.</p>
</sec>
<sec sec-type="author-contributions" id="s5">
<title>Author contributions</title>
<p>ET and SP conceived the project, interpreted the findings, and wrote the manuscript. YL and YP performed the specimen collections and DNA extractions. BW performed PCR amplification and sequencing. KL processed the sequencing data. All authors read and agreed on the final manuscript.</p>
</sec>
</body>
<back>
<sec sec-type="funding-information" id="s6">
<title>Funding</title>
<p>This project was funded by a Yale-NUS College Instrumentation grant to ET, BW, and SP, a Yale-NUS College Summer Research Program to YL, and the Ministry of Education, Singapore, and Yale-NUS College Start-up Grant to ET.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s7">
<title>Publisher&#x00027;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec sec-type="supplementary-material" id="s8">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2023.1199187/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmicb.2023.1199187/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Image_1.tif" id="SM1" mimetype="image/tif" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Figure 1</label>
<caption><p>Distribution of bacteria in samples sorted according to PCoA1 (45.1%). <bold>(A)</bold> Relative abundance of bacterial families; <bold>(B, C)</bold> Relative abundance of bacterial genera for the two host plants.</p></caption> </supplementary-material>
<supplementary-material xlink:href="Image_2.tif" id="SM2" mimetype="image/tif" xmlns:xlink="http://www.w3.org/1999/xlink">
<label>Supplementary Figure 2</label>
<caption><p>Distribution of fungal samples sorted according to PCoA1 (21.6%). <bold>(A)</bold> Relative abundance of fungal families; <bold>(B, C)</bold> Relative abundance of fungal genera for the two host plants.</p></caption> </supplementary-material>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bahrndorff</surname> <given-names>S.</given-names></name> <name><surname>Alemu</surname> <given-names>T.</given-names></name> <name><surname>Alemneh</surname> <given-names>T.</given-names></name> <name><surname>Lund Nielsen</surname> <given-names>J.</given-names></name></person-group> (<year>2016</year>). <article-title>The microbiome of animals: implications for conservation biology</article-title>. <source>Int. J. Genomics</source> <volume>2016</volume>, <fpage>5304028</fpage>. <pub-id pub-id-type="doi">10.1155/2016/5304028</pub-id><pub-id pub-id-type="pmid">27195280</pub-id></citation></ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bl&#x000FC;thgen</surname> <given-names>N.</given-names></name> <name><surname>Metzner</surname> <given-names>A.</given-names></name> <name><surname>Ruf</surname> <given-names>D.</given-names></name></person-group> (<year>2006</year>). <article-title>Food plant selection by stick insects (Phasmida) in a bornean rain forest</article-title>. <source>J. Trop. Ecol.</source> <volume>22</volume>, <fpage>35</fpage>&#x02013;<lpage>40</lpage>. <pub-id pub-id-type="doi">10.1017/S0266467405002749</pub-id></citation>
</ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Boucher</surname> <given-names>S.</given-names></name> <name><surname>Varady-Szabo</surname> <given-names>H.</given-names></name></person-group> (<year>2005</year>). <article-title>Effects of different diets on the survival, longevity and growth rate of the annam stick insect, &#x0003C; span class= &#x0201C;genus-species&#x0201D; &#x0003E;medauroidea extradentata &#x0003C; /span&#x0003E; (Phasmatodea: Phasmatidae)</article-title>. <source>J. Orthoptera Res.</source> <volume>14</volume>, <fpage>115</fpage>&#x02013;<lpage>118</lpage>. <pub-id pub-id-type="doi">10.1665/1082-6467(2005)14(115:EODDOT)2.0.CO;2</pub-id></citation>
</ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bourguignon</surname> <given-names>T.</given-names></name> <name><surname>Lo</surname> <given-names>N.</given-names></name> <name><surname>Dietrich</surname> <given-names>C.</given-names></name> <name><surname>&#x00160;obotn&#x000ED;k</surname> <given-names>J.</given-names></name> <name><surname>Sidek</surname> <given-names>S.</given-names></name> <name><surname>Roisin</surname> <given-names>Y.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>Rampant host switching shaped the termite gut microbiome</article-title>. <source>Curr. Biol.</source> <volume>28</volume>, <fpage>649</fpage>&#x02013;<lpage>654</lpage>.e642. <pub-id pub-id-type="doi">10.1016/j.cub.2018.01.035</pub-id><pub-id pub-id-type="pmid">29429621</pub-id></citation></ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Callahan</surname> <given-names>B. J.</given-names></name> <name><surname>McMurdie</surname> <given-names>P. J.</given-names></name> <name><surname>Rosen</surname> <given-names>M. J.</given-names></name> <name><surname>Han</surname> <given-names>A. W.</given-names></name> <name><surname>Johnson</surname> <given-names>A. J. A.</given-names></name> <name><surname>Holmes</surname> <given-names>S. P.</given-names></name></person-group> (<year>2016</year>). <article-title>DADA2: high-resolution sample inference from Illumina amplicon data</article-title>. <source>Nat. Methods</source> <volume>13</volume>, <fpage>581</fpage>&#x02013;<lpage>583</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.3869</pub-id><pub-id pub-id-type="pmid">27214047</pub-id></citation></ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Caporaso</surname> <given-names>J. G.</given-names></name> <name><surname>Lauber</surname> <given-names>C. L.</given-names></name> <name><surname>Costello</surname> <given-names>E. K.</given-names></name> <name><surname>Berg-Lyons</surname> <given-names>D.</given-names></name> <name><surname>Gonzalez</surname> <given-names>A.</given-names></name> <name><surname>Stombaugh</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Moving pictures of the human microbiome</article-title>. <source>Genome Biol</source>. <volume>12</volume>, <fpage>R50</fpage>. <pub-id pub-id-type="doi">10.1186/gb-2011-12-5-r50</pub-id><pub-id pub-id-type="pmid">21624126</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chaturvedi</surname> <given-names>S.</given-names></name> <name><surname>Rego</surname> <given-names>A.</given-names></name> <name><surname>Lucas</surname> <given-names>L. K.</given-names></name> <name><surname>Gompert</surname> <given-names>Z.</given-names></name></person-group> (<year>2017</year>). <article-title>Sources of variation in the gut microbial community of lycaeides melissa caterpillars</article-title>. <source>Sci. Rep.</source> <volume>7</volume>, <fpage>11335</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-017-11781-1</pub-id><pub-id pub-id-type="pmid">28900218</pub-id></citation></ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Davidson</surname> <given-names>G. L.</given-names></name> <name><surname>Cooke</surname> <given-names>A. C.</given-names></name> <name><surname>Johnson</surname> <given-names>C. N.</given-names></name> <name><surname>Quinn</surname> <given-names>J. L.</given-names></name></person-group> (<year>2018</year>). <article-title>The gut microbiome as a driver of individual variation in cognition and functional behaviour</article-title>. <source>Philos. Trans. R. Soc. B Biol. Sci.</source> <volume>373</volume>, <fpage>20170286</fpage>. <pub-id pub-id-type="doi">10.1098/rstb.2017.0286</pub-id><pub-id pub-id-type="pmid">30104431</pub-id></citation></ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Davidson</surname> <given-names>G. L.</given-names></name> <name><surname>Raulo</surname> <given-names>A.</given-names></name> <name><surname>Knowles</surname> <given-names>S. C. L.</given-names></name></person-group> (<year>2020</year>). <article-title>Identifying microbiome-mediated behaviour in wild vertebrates</article-title>. <source>Trends Ecol. Evol.</source> <volume>35</volume>, <fpage>972</fpage>&#x02013;<lpage>980</lpage>. <pub-id pub-id-type="doi">10.1016/j.tree.2020.06.014</pub-id><pub-id pub-id-type="pmid">32736804</pub-id></citation></ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Davis</surname> <given-names>N. M.</given-names></name> <name><surname>Proctor</surname> <given-names>D. M.</given-names></name> <name><surname>Holmes</surname> <given-names>S. P.</given-names></name> <name><surname>Relman</surname> <given-names>D. A.</given-names></name> <name><surname>Callahan</surname> <given-names>B. J.</given-names></name></person-group> (<year>2018</year>). <article-title>Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data</article-title>. <source>Microbiome</source> <volume>6</volume>, <fpage>226</fpage>&#x02013;<lpage>226</lpage>. <pub-id pub-id-type="doi">10.1186/s40168-018-0605-2</pub-id><pub-id pub-id-type="pmid">30558668</pub-id></citation></ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dillon</surname> <given-names>R. J.</given-names></name> <name><surname>Dillon</surname> <given-names>V. M.</given-names></name></person-group> (<year>2003</year>). <article-title>The gut bacteria of insects: nonpathogenic interactions</article-title>. <source>Ann. Rev. Entomol.</source> <volume>49</volume>, <fpage>71</fpage>&#x02013;<lpage>92</lpage>. <pub-id pub-id-type="doi">10.1146/annurev.ento.49.061802.123416</pub-id><pub-id pub-id-type="pmid">14651457</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Douglas</surname> <given-names>A. E.</given-names></name></person-group> (<year>2019</year>). <article-title>Simple animal models for microbiome research</article-title>. <source>Nat. Rev. Microbiol.</source> <volume>17</volume>, <fpage>764</fpage>&#x02013;<lpage>775</lpage>. <pub-id pub-id-type="doi">10.1038/s41579-019-0242-1</pub-id><pub-id pub-id-type="pmid">31417197</pub-id></citation></ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gardes</surname> <given-names>M.</given-names></name> <name><surname>Bruns</surname> <given-names>T. D.</given-names></name></person-group> (<year>1993</year>). <article-title>ITS primers with enhanced specificity for basidiomycetes - application to the identification of mycorrhizae and rusts</article-title>. <source>Mol. Ecol</source>. <volume>2</volume>, <fpage>113</fpage>&#x02013;<lpage>118</lpage>. <pub-id pub-id-type="doi">10.1111/j.1365-294X.1993.tb00005.x</pub-id><pub-id pub-id-type="pmid">8180733</pub-id></citation></ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gurung</surname> <given-names>K.</given-names></name> <name><surname>Wertheim</surname> <given-names>B.</given-names></name> <name><surname>Falcao Salles</surname> <given-names>J.</given-names></name></person-group> (<year>2019</year>). <article-title>The microbiome of pest insects: it is not just bacteria</article-title>. <source>Entomol. Exp. Appl.</source> <volume>167</volume>, <fpage>156</fpage>&#x02013;<lpage>170</lpage>. <pub-id pub-id-type="doi">10.1111/eea.12768</pub-id></citation>
</ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hannula</surname> <given-names>S. E.</given-names></name> <name><surname>Zhu</surname> <given-names>F.</given-names></name> <name><surname>Heinen</surname> <given-names>R.</given-names></name> <name><surname>Bezemer</surname> <given-names>T. M.</given-names></name></person-group> (<year>2019</year>). <article-title>Foliar-feeding insects acquire microbiomes from the soil rather than the host plant</article-title>. <source>Nat. Commun</source>. <volume>10</volume>, <fpage>1254</fpage>. <pub-id pub-id-type="doi">10.1038/s41467-019-09284-w</pub-id><pub-id pub-id-type="pmid">30890706</pub-id></citation></ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hern&#x000E1;ndez-Garc&#x000ED;a</surname> <given-names>J. A.</given-names></name> <name><surname>Briones-Roblero</surname> <given-names>C. I.</given-names></name> <name><surname>Rivera-Ordu&#x000F1;a</surname> <given-names>F. N.</given-names></name> <name><surname>Z&#x000FA;&#x000F1;iga</surname> <given-names>G.</given-names></name></person-group> (<year>2017</year>). <article-title>Revealing the gut bacteriome of dendroctonus bark beetles (curculionidae: scolytinae): diversity, core members and co-evolutionary patterns</article-title>. <source>Sci. Rep.</source> <volume>7</volume>, <fpage>13864</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-017-14031-6</pub-id><pub-id pub-id-type="pmid">29066751</pub-id></citation></ref>
<ref id="B17">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Huttenhower</surname> <given-names>C.</given-names></name> <name><surname>Gevers</surname> <given-names>D.</given-names></name> <name><surname>Knight</surname> <given-names>R.</given-names></name> <name><surname>Abubucker</surname> <given-names>S.</given-names></name> <name><surname>Badger</surname> <given-names>J. H.</given-names></name> <name><surname>Chinwalla</surname> <given-names>A. T.</given-names></name> <etal/></person-group>. (<year>2012</year>). <article-title>Structure, function and diversity of the healthy human microbiome</article-title>. <source>Nature</source> <volume>486</volume>, <fpage>207</fpage>&#x02013;<lpage>214</lpage>. <pub-id pub-id-type="doi">10.1038/nature11234</pub-id><pub-id pub-id-type="pmid">22699609</pub-id></citation></ref>
<ref id="B18">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ignasiak</surname> <given-names>K.</given-names></name> <name><surname>Maxwell</surname> <given-names>A.</given-names></name></person-group> (<year>2017</year>). <article-title>Antibiotic-resistant bacteria in the guts of insects feeding on plants: prospects for discovering plant-derived antibiotics</article-title>. <source>BMC Microbiol.</source> <volume>17</volume>, <fpage>223</fpage>. <pub-id pub-id-type="doi">10.1186/s12866-017-1133-0</pub-id><pub-id pub-id-type="pmid">29191163</pub-id></citation></ref>
<ref id="B19">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jang</surname> <given-names>S.</given-names></name> <name><surname>Kikuchi</surname> <given-names>Y.</given-names></name></person-group> (<year>2020</year>). <article-title>Impact of the insect gut microbiota on ecology, evolution, and industry</article-title>. <source>Curr. Opin. Insect Sci.</source> <volume>41</volume>, <fpage>33</fpage>&#x02013;<lpage>39</lpage>. <pub-id pub-id-type="doi">10.1016/j.cois.2020.06.004</pub-id><pub-id pub-id-type="pmid">32634703</pub-id></citation></ref>
<ref id="B20">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jones</surname> <given-names>A. G.</given-names></name> <name><surname>Mason</surname> <given-names>C. J.</given-names></name> <name><surname>Felton</surname> <given-names>G. W.</given-names></name> <name><surname>Hoover</surname> <given-names>K.</given-names></name></person-group> (<year>2019</year>). <article-title>Host plant and population source drive diversity of microbial gut communities in two polyphagous insects</article-title>. <source>Sci. Rep.</source> <volume>9</volume>, <fpage>2792</fpage>. <pub-id pub-id-type="doi">10.1038/s41598-019-39163-9</pub-id><pub-id pub-id-type="pmid">30808905</pub-id></citation></ref>
<ref id="B21">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jones</surname> <given-names>J. C.</given-names></name> <name><surname>Fruciano</surname> <given-names>C.</given-names></name> <name><surname>Marchant</surname> <given-names>J.</given-names></name> <name><surname>Hildebrand</surname> <given-names>F.</given-names></name> <name><surname>Forslund</surname> <given-names>S.</given-names></name> <name><surname>Bork</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>The gut microbiome is associated with behavioural task in honey bees</article-title>. <source>Insectes Soc.</source> <volume>65</volume>, <fpage>419</fpage>&#x02013;<lpage>429</lpage>. <pub-id pub-id-type="doi">10.1007/s00040-018-0624-9</pub-id><pub-id pub-id-type="pmid">30100619</pub-id></citation></ref>
<ref id="B22">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Koch</surname> <given-names>H.</given-names></name> <name><surname>Schmid-Hempel</surname> <given-names>P.</given-names></name></person-group> (<year>2011</year>). <article-title>Socially transmitted gut microbiota protect bumble bees against an intestinal parasite</article-title>. <source>Proc. Nat. Acad. Sci.</source> <volume>108</volume>, <fpage>19288</fpage>&#x02013;<lpage>19292</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1110474108</pub-id><pub-id pub-id-type="pmid">22084077</pub-id></citation></ref>
<ref id="B23">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>K&#x000F6;nig</surname> <given-names>H.</given-names></name> <name><surname>Li</surname> <given-names>L.</given-names></name> <name><surname>Fr&#x000F6;hlich</surname> <given-names>J.</given-names></name></person-group> (<year>2013</year>). <article-title>The cellulolytic system of the termite gut</article-title>. <source>Appl. Microbiol. Biotechnol.</source> <volume>97</volume>, <fpage>7943</fpage>&#x02013;<lpage>7962</lpage>. <pub-id pub-id-type="doi">10.1007/s00253-013-5119-z</pub-id><pub-id pub-id-type="pmid">23900801</pub-id></citation></ref>
<ref id="B24">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Martin</surname> <given-names>M.</given-names></name></person-group> (<year>2011</year>). <article-title>Cutadapt removes adapter sequences from high-throughput sequencing reads</article-title>. <source>EMBnet. J.</source> <volume>17</volume>, <fpage>10</fpage>&#x02013;<lpage>12</lpage>. <pub-id pub-id-type="doi">10.14806/ej.17.1.200</pub-id><pub-id pub-id-type="pmid">28715235</pub-id></citation></ref>
<ref id="B25">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mason</surname> <given-names>C. J.</given-names></name></person-group> (<year>2020</year>). <article-title>Complex relationships at the intersection of insect gut microbiomes and plant defenses</article-title>. <source>J. Chem. Ecol.</source> <volume>46</volume>, <fpage>793</fpage>&#x02013;<lpage>807</lpage>. <pub-id pub-id-type="doi">10.1007/s10886-020-01187-1</pub-id><pub-id pub-id-type="pmid">32537721</pub-id></citation></ref>
<ref id="B26">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mason</surname> <given-names>C. J.</given-names></name> <name><surname>Raffa</surname> <given-names>K. F.</given-names></name></person-group> (<year>2014</year>). <article-title>Acquisition and structuring of midgut bacterial communities in gypsy moth (lepidoptera: erebidae) larvae</article-title>. <source>Environ. Entomol.</source> <volume>43</volume>, <fpage>595</fpage>&#x02013;<lpage>604</lpage>. <pub-id pub-id-type="doi">10.1603/EN14031</pub-id><pub-id pub-id-type="pmid">24780292</pub-id></citation></ref>
<ref id="B27">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mason</surname> <given-names>C. J.</given-names></name> <name><surname>St. Clair</surname> <given-names>A.</given-names></name> <name><surname>Peiffer</surname> <given-names>M.</given-names></name> <name><surname>Gomez</surname> <given-names>E.</given-names></name> <name><surname>Jones</surname> <given-names>A. G.</given-names></name> <name><surname>Felton</surname> <given-names>G. W.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Diet influences proliferation and stability of gut bacterial populations in herbivorous lepidopteran larvae</article-title>. <source>PLoS ONE</source> <volume>15</volume>, <fpage>e0229848</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0229848</pub-id><pub-id pub-id-type="pmid">32168341</pub-id></citation></ref>
<ref id="B28">
<citation citation-type="thesis"><person-group person-group-type="author"><name><surname>McLean</surname> <given-names>A. H. C.</given-names></name> <name><surname>Godfray</surname> <given-names>H. C. J.</given-names></name> <name><surname>Ellers</surname> <given-names>J.</given-names></name> <name><surname>Henry</surname> <given-names>L. M.</given-names></name></person-group> (<year>2019</year>). <article-title>Host relatedness influences the composition of aphid microbiomes</article-title>. <source>Environ. Microbiol. Rep.</source> <volume>11</volume>, <fpage>808</fpage>&#x02013;<lpage>816</lpage>. <pub-id pub-id-type="doi">10.1111/1758-2229.12795</pub-id><pub-id pub-id-type="pmid">31573138</pub-id></citation></ref>
<ref id="B29">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>McMurdie</surname> <given-names>P. J.</given-names></name> <name><surname>Holmes</surname> <given-names>S.</given-names></name></person-group> (<year>2013</year>). <article-title>phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data</article-title>. <source>PLoS ONE</source> <volume>8</volume>:<fpage>e61217</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0061217</pub-id><pub-id pub-id-type="pmid">23630581</pub-id></citation></ref>
<ref id="B30">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mendes</surname> <given-names>R.</given-names></name> <name><surname>Raaijmakers</surname> <given-names>J. M.</given-names></name></person-group> (<year>2015</year>). <article-title>Cross-kingdom similarities in microbiome functions</article-title>. <source>ISME J.</source> <volume>9</volume>, <fpage>1905</fpage>&#x02013;<lpage>1907</lpage>. <pub-id pub-id-type="doi">10.1038/ismej.2015.7</pub-id><pub-id pub-id-type="pmid">25647346</pub-id></citation></ref>
<ref id="B31">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Moran</surname> <given-names>N. A.</given-names></name> <name><surname>Ochman</surname> <given-names>H.</given-names></name> <name><surname>Hammer</surname> <given-names>T. J.</given-names></name></person-group> (<year>2019</year>). <article-title>Evolutionary and ecological consequences of gut microbial communities</article-title>. <source>Annu. Rev. Ecol. Evol. Syst.</source> <volume>50</volume>, <fpage>451</fpage>&#x02013;<lpage>475</lpage>. <pub-id pub-id-type="doi">10.1146/annurev-ecolsys-110617-062453</pub-id><pub-id pub-id-type="pmid">32733173</pub-id></citation></ref>
<ref id="B32">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mu&#x000F1;oz-Benavent</surname> <given-names>M.</given-names></name> <name><surname>P&#x000E9;rez-Cobas</surname> <given-names>A. E.</given-names></name> <name><surname>Garc&#x000ED;a-Ferris</surname> <given-names>C.</given-names></name> <name><surname>Moya</surname> <given-names>A.</given-names></name> <name><surname>Latorre</surname> <given-names>A.</given-names></name></person-group> (<year>2021</year>). <article-title>Insects&#x00027; potential: understanding the functional role of their gut microbiome</article-title>. <source>J. Pharm. Biomed. Anal.</source> <volume>194</volume>:<fpage>113787</fpage>. <pub-id pub-id-type="doi">10.1016/j.jpba.2020.113787</pub-id><pub-id pub-id-type="pmid">33272789</pub-id></citation></ref>
<ref id="B33">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nilsson</surname> <given-names>R. H.</given-names></name> <name><surname>Larsson</surname> <given-names>K. -H.</given-names></name> <name><surname>Taylor</surname> <given-names>A. F. S.</given-names></name> <name><surname>Bengtsson-Palme</surname> <given-names>J.</given-names></name> <name><surname>Jeppesen</surname> <given-names>T. S.</given-names></name> <name><surname>Schigel</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2018</year>). <article-title>The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications</article-title>. <source>Nucleic Acids Res</source>. <volume>47</volume>, <fpage>D259</fpage>&#x02013;<lpage>D264</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gky1022</pub-id><pub-id pub-id-type="pmid">30371820</pub-id></citation></ref>
<ref id="B34">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quast</surname> <given-names>C.</given-names></name> <name><surname>Pruesse</surname> <given-names>E.</given-names></name> <name><surname>Yilmaz</surname> <given-names>P.</given-names></name> <name><surname>Gerken</surname> <given-names>J.</given-names></name> <name><surname>Schweer</surname> <given-names>T.</given-names></name> <name><surname>Yarza</surname> <given-names>P.</given-names></name> <etal/></person-group>. (<year>2013</year>). <article-title>The SILVA ribosomal RNA gene database project: improved data processing and web-based tools</article-title>. <source>Nucleic Acids Res.</source> <volume>41</volume>, <fpage>D590</fpage>&#x02013;<lpage>D596</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gks1219</pub-id><pub-id pub-id-type="pmid">23193283</pub-id></citation></ref>
<ref id="B35">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Raymann</surname> <given-names>K.</given-names></name> <name><surname>Moran</surname> <given-names>N. A.</given-names></name></person-group> (<year>2018</year>). <article-title>The role of the gut microbiome in health and disease of adult honey bee workers</article-title>. <source>Curr. Opin. Insect Sci.</source> <volume>26</volume>, <fpage>97</fpage>&#x02013;<lpage>104</lpage>. <pub-id pub-id-type="doi">10.1016/j.cois.2018.02.012</pub-id><pub-id pub-id-type="pmid">29764668</pub-id></citation></ref>
<ref id="B36">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shelomi</surname> <given-names>M.</given-names></name> <name><surname>Danchin</surname> <given-names>E. G. J.</given-names></name> <name><surname>Heckel</surname> <given-names>D.</given-names></name> <name><surname>Wipfler</surname> <given-names>B.</given-names></name> <name><surname>Bradler</surname> <given-names>S.</given-names></name> <name><surname>Zhou</surname> <given-names>X.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Horizontal gene transfer of pectinases from bacteria preceded the diversification of stick and leaf insects</article-title>. <source>Sci. Rep.</source> <volume>6</volume>, <fpage>26388</fpage>. <pub-id pub-id-type="doi">10.1038/srep26388</pub-id><pub-id pub-id-type="pmid">34312432</pub-id></citation></ref>
<ref id="B37">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shelomi</surname> <given-names>M.</given-names></name> <name><surname>Jasper</surname> <given-names>W. C.</given-names></name> <name><surname>Atallah</surname> <given-names>J.</given-names></name> <name><surname>Kimsey</surname> <given-names>L. S.</given-names></name> <name><surname>Johnson</surname> <given-names>B. R.</given-names></name></person-group> (<year>2014a</year>). <article-title>Differential expression of endogenous plant cell wall degrading enzyme genes in the stick insect (Phasmatodea) midgut</article-title>. <source>BMC Genomics</source> <volume>15</volume>, <fpage>917</fpage>. <pub-id pub-id-type="doi">10.1186/1471-2164-15-917</pub-id><pub-id pub-id-type="pmid">25331961</pub-id></citation></ref>
<ref id="B38">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shelomi</surname> <given-names>M.</given-names></name> <name><surname>Lo</surname> <given-names>W. -S.</given-names></name> <name><surname>Kimsey</surname> <given-names>L. S.</given-names></name> <name><surname>Kuo</surname> <given-names>C. -H.</given-names></name></person-group> (<year>2013</year>). <article-title>Analysis of the gut microbiota of walking sticks (Phasmatodea)</article-title>. <source>BMC Res. Notes</source> <volume>6</volume>, <fpage>368</fpage>. <pub-id pub-id-type="doi">10.1186/1756-0500-6-368</pub-id><pub-id pub-id-type="pmid">24025149</pub-id></citation></ref>
<ref id="B39">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shelomi</surname> <given-names>M.</given-names></name> <name><surname>Sitepu</surname> <given-names>I. R.</given-names></name> <name><surname>Boundy-Mills</surname> <given-names>K. L.</given-names></name> <name><surname>Kimsey</surname> <given-names>L. S.</given-names></name></person-group> (<year>2015</year>). <article-title>Review of the gross anatomy and microbiology of the phasmatodea digestive tract</article-title>. <source>J. Orthoptera Res.</source> <volume>24</volume>, <fpage>29</fpage>&#x02013;<lpage>40</lpage>. <pub-id pub-id-type="doi">10.1665/034.024.0105</pub-id></citation>
</ref>
<ref id="B40">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shelomi</surname> <given-names>M.</given-names></name> <name><surname>Watanabe</surname> <given-names>H.</given-names></name> <name><surname>Arakawa</surname> <given-names>G.</given-names></name></person-group> (<year>2014b</year>). <article-title>Endogenous cellulase enzymes in the stick insect (Phasmatodea) gut</article-title>. <source>J. Insect Physiol.</source> <volume>60</volume>, <fpage>25</fpage>&#x02013;<lpage>30</lpage>. <pub-id pub-id-type="doi">10.1016/j.jinsphys.2013.10.007</pub-id><pub-id pub-id-type="pmid">24216471</pub-id></citation></ref>
<ref id="B41">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname> <given-names>D. P.</given-names></name> <name><surname>Peay</surname> <given-names>K. G.</given-names></name></person-group> (<year>2014</year>). <article-title>Sequence depth, not PCR replication, improves ecological inference from next generation DNA sequencing</article-title>. <source>PLoS ONE</source> <volume>9</volume>, <fpage>e90234</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0090234</pub-id><pub-id pub-id-type="pmid">24587293</pub-id></citation></ref>
<ref id="B42">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Thompson</surname> <given-names>L. R.</given-names></name> <name><surname>Sanders</surname> <given-names>J. G.</given-names></name> <name><surname>McDonald</surname> <given-names>D.</given-names></name> <name><surname>Amir</surname> <given-names>A.</given-names></name> <name><surname>Ladau</surname> <given-names>J.</given-names></name> <name><surname>Locey</surname> <given-names>K. J.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>A communal catalogue reveals Earth&#x00027;s multiscale microbial diversity</article-title>. <source>Nature</source> <volume>551</volume>, <fpage>457</fpage>&#x02013;<lpage>463</lpage>. <pub-id pub-id-type="doi">10.1038/nature24621</pub-id><pub-id pub-id-type="pmid">29088705</pub-id></citation></ref>
<ref id="B43">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>V&#x000E1;zquez-Montoya</surname> <given-names>E. L.</given-names></name> <name><surname>Castro-Ochoa</surname> <given-names>L. D.</given-names></name> <name><surname>Maldonado-Mendoza</surname> <given-names>I. E.</given-names></name> <name><surname>Luna-Su&#x000E1;rez</surname> <given-names>S.</given-names></name> <name><surname>Castro-Mart&#x000ED;nez</surname> <given-names>C.</given-names></name></person-group> (<year>2020</year>). <article-title>Moringa straw as cellulase production inducer and cellulolytic fungi source</article-title>. <source>Rev. Argent. Microbiol.</source> <volume>52</volume>, <fpage>4</fpage>&#x02013;<lpage>12</lpage>. <pub-id pub-id-type="doi">10.1016/j.ram.2019.02.005</pub-id><pub-id pub-id-type="pmid">31204058</pub-id></citation></ref>
<ref id="B44">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>Y.</given-names></name> <name><surname>Naumann</surname> <given-names>U.</given-names></name> <name><surname>Wright</surname> <given-names>S. T.</given-names></name> <name><surname>Warton</surname> <given-names>D. I.</given-names></name></person-group> (<year>2012</year>). <article-title>mvabund&#x02013; an R package for model-based analysis of multivariate abundance data</article-title>. <source>Methods Ecol. Evol.</source> <volume>3</volume>, <fpage>471</fpage>&#x02013;<lpage>474</lpage>. <pub-id pub-id-type="doi">10.1111/j.2041-210X.2012.00190.x</pub-id></citation>
</ref>
<ref id="B45">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Watanabe</surname> <given-names>H.</given-names></name> <name><surname>Tokuda</surname> <given-names>G.</given-names></name></person-group> (<year>2010</year>). <article-title>Cellulolytic systems in insects</article-title>. <source>Annu. Rev. Entomol.</source> <volume>55</volume>, <fpage>609</fpage>&#x02013;<lpage>632</lpage>. <pub-id pub-id-type="doi">10.1146/annurev-ento-112408-085319</pub-id><pub-id pub-id-type="pmid">19754245</pub-id></citation></ref>
<ref id="B46">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Whitaker</surname> <given-names>M. R. L.</given-names></name> <name><surname>Salzman</surname> <given-names>S.</given-names></name> <name><surname>Sanders</surname> <given-names>J.</given-names></name> <name><surname>Kaltenpoth</surname> <given-names>M.</given-names></name> <name><surname>Pierce</surname> <given-names>N. E.</given-names></name></person-group> (<year>2016</year>). <article-title>Microbial communities of lycaenid butterflies do not correlate with larval diet</article-title>. <volume>7</volume>, <fpage>1920</fpage>. <pub-id pub-id-type="doi">10.3389/fmicb.2016.01920</pub-id><pub-id pub-id-type="pmid">27965647</pub-id></citation></ref>
<ref id="B47">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>White</surname> <given-names>T. J.</given-names></name> <name><surname>Bruns</surname> <given-names>T.</given-names></name> <name><surname>Lee</surname> <given-names>S.</given-names></name> <name><surname>Taylor</surname> <given-names>J.</given-names></name></person-group> (<year>1990</year>). <article-title>&#x0201C;Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics,&#x0201D;</article-title> in <source>PCR Protocols, A Guide to Methods and Applications</source>, eds M. A. Innes, D. H. Gelfand, J. S. Sninsky, and T. J. White (<publisher-loc>San Diego, CA</publisher-loc>: <publisher-name>Academic Press</publisher-name>), <fpage>315</fpage>&#x02013;<lpage>322</lpage>.</citation>
</ref>
<ref id="B48">
<citation citation-type="book"><person-group person-group-type="author"><name><surname>Wickham</surname> <given-names>H.</given-names></name></person-group> (<year>2016</year>). <source>ggplot2: Elegant Graphics for Data Analysis</source>. <publisher-loc>New York, NY</publisher-loc>: <publisher-name>Springer-Verlag</publisher-name>.</citation>
</ref>
<ref id="B49">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xu</surname> <given-names>N.</given-names></name> <name><surname>Zhao</surname> <given-names>Q.</given-names></name> <name><surname>Zhang</surname> <given-names>Z.</given-names></name> <name><surname>Zhang</surname> <given-names>Q.</given-names></name> <name><surname>Wang</surname> <given-names>Y.</given-names></name> <name><surname>Qin</surname> <given-names>G.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Phyllosphere microorganisms: sources, drivers, and their interactions with plant hosts</article-title>. <source>J. Agric. Food Chem.</source> <volume>70</volume>, <fpage>4860</fpage>&#x02013;<lpage>4870</lpage>. <pub-id pub-id-type="doi">10.1021/acs.jafc.2c01113</pub-id><pub-id pub-id-type="pmid">35435673</pub-id></citation></ref>
</ref-list> 
</back>
</article> 