Patterns of Microbiome Composition Vary Across Spatial Scales in a Specialist Insect

Microbial communities associated with animals vary based on both intrinsic and extrinsic factors. Of many possible determinants affecting microbiome composition, host phylogeny, host diet, and local environment are the most important. How these factors interact across spatial scales is not well understood. Here, we seek to identify the main influences on microbiome composition in a specialist insect, the western corn rootworm (WCR; Diabrotica virgifera virgifera), by analyzing the bacterial communities of adults collected from their obligate host plant, corn (Zea mays), across several geographic locations and comparing the patterns in communities to its congeneric species, the northern corn rootworm (NCR; Diabrotica barberi). We found that bacterial communities of WCR and NCR shared a portion of their bacterial communities even when collected from disparate locations. However, within each species, the location of collection significantly influenced the composition of their microbiome. Correlations of geographic distance between sites with WCR bacterial community composition revealed different patterns at different spatial scales. Community similarity decreased with increased geographic distance at smaller spatial scales (~25 km between the nearest sites). At broad spatial scales (>200 km), community composition was not correlated with distances between sites, but instead reflected the historical invasion path of WCR across the United States. These results suggest bacterial communities are structured directly by dispersal dynamics at small, regional spatial scales, while landscape-level genetic or environmental differences may drive community composition across broad spatial scales in this specialist insect.

Introduction

Animals have evolved while in constant contact with microorganisms. Associations between hosts and microorganisms exist on a continuum from beneficial to detrimental. In insects, bacterial communities can improve host fitness by enhancing nutrition (Ben-Yosef et al., 2014), disrupting plant defenses (Chu et al., 2013), and protecting against disease (Miller et al., 2021), but are also capable of inducing mortality in certain situations (Caccia et al., 2016; Mason et al., 2019). Variation in bacterial communities within hosts may result in differential survival with direct implications for design and implementation of conservation and pest management strategies (Paddock et al., 2021).

Inter- and intra-species variation in the microbiome can be influenced by numerous factors. Host species identity can dictate what microbes survive within the host. Insect guts vary in morphology, pH, and immune response (Caccia et al., 2019), which serve as filters for specific microbes resulting in communities that vary between closely related species (Adair et al., 2020). Host diet partially determines the local species pool with which the host interacts. Different feeding modalities (e.g., chewing vs. sucking mouthparts) constrain access to food sources, which can affect which microbes are associated with insects (Huang et al., 2021). Furthermore, certain food substrates can be digested by microbes within insects, and thus can alter communities through resource limitation or niche partitioning (Mason et al., 2020; Brochet et al., 2021). The host’s external environment can also generate variation between insect microbiomes by affecting the local source pool of microbes. Such factors as temperature (Wang et al., 2020), landscape context (Park et al., 2019), and plant diversity (Cohen et al., 2020) may be associated with variation between insect microbiomes in different local habitats.

Host microbiomes across a biogeographic space can best be described through a metacommunity framework (Adair and Douglas, 2017; Miller et al., 2018), where host-associated communities in local environments are subsets of the larger environmental metacommunity and linked through dispersal. Spatial limits on microbial dispersal, or “dispersal limitation,” can result in patterns of geospatial correlation in which microbial community similarity decays with increasing geographic distance (Finkel et al., 2012; Moeller et al., 2017). At different spatial scales, the determinants of microbiome composition may change and result in different communities. Both landscape composition and configuration can dictate the dispersal capacity of microorganisms within the metacommunity (Parajuli et al., 2020). For example, in humans, small-scale dispersal events may be disrupted by barriers such as human-made structures or vegetation (Parajuli et al., 2020). At a continental scale, dispersal may be very limited and other factors such as lifestyle, diet, age, and genetics may account for most variation between hosts (Yatsunenko et al., 2012). Which microbes are associated with hosts is determined by interactions of both deterministic and stochastic processes (Adair and Douglas, 2017), but how these processes interact across spatial scales is not well understood.

Corn rootworms (genus Diabrotica; Coleoptera, Chrysomelidae) represent a useful system to investigate the influence of biogeographical arrangement on host associated microbiome composition. Many studies have investigated determinants of microbiome composition in generalist feeding insects (Adair et al., 2018; Jones et al., 2019; Cohen et al., 2020). Studies with generalist insects have limitations on distinguishing between the influence of geographic location and the influence of diet on microbiome composition. In the United States, two rootworm species predominate east of the Rocky Mountains: the western corn rootworm (Diabrotica virgifera virgifera LeConte; WCR) and the northern corn rootworm (Diabrotica barberi Smith & Lawrence; NCR). Both species overlap in distribution, phenology, and host plant usage (Krysan and Miller, 1986), which allows us to better examine environmental factors influencing microbiome composition. In addition, they comprise the most damaging groups of corn pests in the United States, with management costs and yield losses combining for over $2 billion annually (Wechsler and Smith, 2018). Management continues to become more complicated as both species have evolved resistance to crop rotation and transgenic crops producing toxins derived from the bacterium Bacillus thuringiensis (Bt; Levine et al., 1992, 2002; Calles-Torrez et al., 2019). In WCR, the microbiome may help beetles overcome plant defenses, and changes in larval microbiome composition are linked to resistance to Bt (Chu et al., 2013; Paddock et al., 2021). The geospatial consistency of bacterial community composition within WCR has yet to be investigated.

Here, we characterize the bacterial communities in two sister Diabrotica species and compare the patterns of assembly at different spatial scales. We collected WCR beetles from corn fields across their ranges in the United States using two sampling scales, regional scale (~12–50 km between sites along linear transects) and continental scale (>200 km). We examined correlations between bacterial community dissimilarity and distance between collection sites for WCR. We hypothesize that environment contributes to variation in microbial communities, and thus, differentially influences the similarity of microbial communities across biogeographical space. Specifically, we predict that microbiome sequence similarity decreases with geographic distance (distance decay), and that the effect is stronger at the landscape level where dispersal limitation may be higher compared to local levels. Increased understanding of the factors influencing microbiome composition in insects can provide insight into how microbe-mediated effects on plant-insect interactions may have evolved.

Materials and Methods

Insects

WCR and NCR are univoltine species that emerge in late July into August. In this study, adult beetles were collected from corn fields between July and August in 2016, 2019, and 2020. Information including date of collection and location of the closest city or populated place to the collection site can be found in Supplementary Information. WCR were collected at two spatial scales (Figure 1). A small, regional-scale sampling scheme consisted of field sites located ~25 km apart along two transects in eastern Colorado and western Kansas. On average, WCR can disperse ~17 m a day. Distances between any two sites ranged from 11.46–276.73 km. Beetles were collected in ethanol and stored at −20°C at the Corn Insects and Crop Genetics Research Unit (CICGRU) in Ames, Iowa until DNA extraction. For the broad, continental-scale sampling scheme, field sites were ~250 km apart, with paired distances ranging from 248 to 3,122 km. Beetles were stored in 95% ethanol and shipped to Columbia, Missouri, where they remained until DNA was extracted. A total of 24 sites across 11 states were sampled for WCR in 2016, 2018, and 2020. Within the small regional-scale sampling scheme, 14 city level sample sites were collected in 2016. Eight beetles were collected at each site for a total of 192 WCR beetles. NCR were collected only at a continental scale and were processed in the same manner as WCR. A total of four sites were sampled for NCR in 2020. Within Missouri, two city level samples were collected. Eight beetles were collected at each site for a total of 32 NCR beetles.

FIGURE 1

Figure 1. Sample locations of wild, adult Diabrotica virgifera virgifera (black dots) and D. barberi (gray dots). Beetles were collected from cornfields during late July to early August. Colorado is zoomed out with small circles showing close-proximity sampling locations (note one location is across the border in Kansas). Each dot consists of one sample site where eight individual beetles were collected.

DNA Extraction and 16S rRNA Gene Amplification

For each sample location, DNA was extracted from individual beetles (n = 8). Before the day of extraction, beetles in ethanol were removed and placed in individual 2-ml tubes to dry overnight. Beetles were transferred to beaded tubes (MP Biomedicals, Santa Ana, CA, catalog no. 116913500), flash frozen in liquid nitrogen for 20 s, immediately placed in a single-tube bead beater, and shaken for 20 s to pulverize the sample. Bacterial DNA was extracted from pulverized beetles using PowerFecal® DNA Isolation Kit (QIAGEN, catalog no. 12830-50) in accordance with the manufacturer’s protocols.¹ Frozen beetles processed at CICGRU were ground by mortar and pestle and then transferred to sterile 1.5-ml microcentrifuge tubes. Bacterial DNA was extracted using DNA Blood and Tissue kit (QIAGEN, catalog no. 69504) in accordance with the manufacturer’s protocols.² The concentration of extracted DNA was quantified using a Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA). DNA was stored at −80°C until further downstream processing. Prior to amplification, DNA concentration was standardized to 3.51 ng/μl. The V4 hypervariable region of the 16S rRNA gene was amplified using single indexed universal primers (U515F/806R) with Illumina standard adapter sequences. Dual-indexed forward and reverse primers were used in all reactions. PCR reaction steps were as follows: 98°C^(3:00) + 25 cycles of [98°C^(0:15) + 50°C^(0:30) + 72°C^(0:30)]. The resulting amplicons (5 μl) were pooled before sequencing on the Illumina MiSeq 2 × 250 bp platform (Paddock et al., 2021). The construction and sequencing of 16S sequencing amplicon libraries were completed at the University of Missouri (MU) DNA Core facility.

16S rRNA Community Analysis

Sequence assembly and annotation were performed by the MU Informatics Research Core Facility. Primers were trimmed using Cutadapt in two rounds,³ first removing forward primers with an error rate of 0.11 mismatches and minimum length of 19 bp, followed by a second round of trimming from the 3′ end executed with an error rate of 0.1 mismatches and minimum length of 20 bp. A minimal overlap of 3 with the 3′ end of the primer sequence was required for removal. Untrimmed contigs were discarded between rounds of trimming. Using the Qiime2 plugin, DADA2 (Callahan et al., 2016; version 1.10.0), forward and reverse reads were truncated to 150 bp and discarded if the number of expected errors was >2.0. Chimeras were detected using the “consensus” method and removed. Resulting sequences were filtered to retain only sequences 249–257 bases long. Taxonomy was assigned to amplicon sequence variants (ASVs) using the Silva.v132 database with the “sklearn” classifier in Qiime2. ASVs were compiled into biom tables for data analysis. ASV tables with metadata and associated taxonomy were imported in to RStudio version 3.5.2 for downstream analysis. ASVs matching chloroplast, mitochondria, and archaea sequences were filtered and removed using phyloseq::filter_taxa in RStudio (McMurdie and Holmes, 2013). Taxa labeled “uncharacterized” at the phylum level were also removed. An extraction blank was used to remove contaminant sequences based on prevalence using decontam in RStudio. The resulting data set was rarefied to a read depth of 750 and subsequently used for analyses of alpha and beta diversity in RStudio.

Statistical Analysis

To examine differences in microbial communities within host species, we removed Wolbachia from the data. Wolbachia can account for over 90% of the relative abundance in microbial communities (Ludwick et al., 2019; Paddock et al., 2021), and its association is nearly ubiquitous (Giordano et al., 1997). Conversely, Wolbachia infection in NCR populations exhibits strong geographic partitioning (Roehrdanz and Levine, 2007). Analyses conducted with Wolbachia can be found in the Supplementary Table S4 and Supplementary Figures S2, S3 Alpha diversity indices (Chao-1 and inverse Simpson’s D) for both host species across sites were generated using estimate_richness in the phyloseq package. Analyses of beta-diversity were conducted using a permutational multivariate analysis of variance (PERMANOVA) based on Bray–Curtis and Jaccard distances at the ASV level. Comparisons of differentially abundant taxa present between species were conducted using DESeq2 in RStudio using unrarefied data (Love et al., 2014).

We asked whether the state where a sample was collected from affected microbiome composition within the WCR and NCR. To do this, we analyzed alpha- and beta-diversity separately for each species. For WCR, we also compared microbiome composition across cities within the smaller spatial scale in Colorado and Kansas. Alpha diversity indices (Chao-1 and inverse Simpson’s D) were log-transformed before analysis of variance to correct for non-normal distributions, and significance considered as p < 0.05. Analyses of beta-diversity were conducted using a permutational multivariate analysis of variance PERMANOVA based on Bray-Curtis and Jaccard distances at the ASV level. Year of collection was treated as a random factor by restricting permutations within year. Beta dispersion as measured by the average distance to centroid of each group was calculated using vegan::betadisper and compared using a permutational test with vegan::permutest. Pairwise comparisons were made for significant differences observed in PERMANOVA using pairwiseAdonis::pairwiseadonis2 at corrected p < 0.05 (Martinez Arbizu, 2020).

To better understand the patterns of community composition at different scales in WCR, we examined the correlation between community dissimilarity and geospatial arrangement. First, we calculated microbiome dissimilarity between individual insects using both Bray-Curtis and Jaccard distances at the ASV level. Geographic distance between sample locations was calculated using the Haversine formula for distance (Robusto, 1957). Resulting distance matrices were analyzed by mantel test using Spearman’s rank-ordered correlation permuted 9,999 times using vegan::mantel. Secondly, values for PC1 for the WCR, individual-species PCoA were extracted using the scores function and then used as the dependent variable in a linear model examining the relationship to geographic distance from a historically relevant origin source region defined as the sample site closest to Mexico along the invasion path (Arizona). To visualize differences between ASVs in WCR collected from different locations, we generated a heat map using the log10 abundance of 125 most abundant ASVs.

Data Availability

Raw sequences can be found on the NCBI SRA database under the project accession number PRJNA785968. Raw data and accompanying metadata can be found at FigShare at doi: 10.6084/m9.figshare.17130035. Code used in statistical analyses can be found at FigShare at doi: 10.6084/m9.figshare.17130455.

Results

WCR and NCR Bacterial Community Composition

Beetle microbiomes were mainly composed of bacteria from the phyla Proteobacteria, Firmicutes, Actinobacteria, and Bacteriodetes, in order of decreasing relative abundance (Figure 2A). A complete list of all ASVs found within WCR and NCR is located in Supplementary Table 1. Overall, 230 ASVs (~12% of the total ASVs identified) were shared in at least one WCR beetle and one NCR beetle (Figure 2B). These shared ASVs composed 88.2% (±1.5%) of reads in WCR and 91.2% (±3.0%) of NCR reads on average. Shared ASVs can be found in Supplementary Table 2. Seven ASVs were found in each species under a prevalence threshold of 50% presence across all samples per species (Supplementary Table 3). Four of those ASVs were shared between WCR and NCR: two from the family Enterobacteriaceae, and one each from the genera Pantoea and Lactococcus. ASVs from an unclassified genus in Enterobacteriaceae had the highest relative abundance in both WCR (54.3%) and NCR (40.4%). Other taxa with high relative abundance in WCR were in the genera Lactococcus, Acinetobacter, Pseudomonas, Serratia, Pantoea, Microbacterium, Massilia, Sphingomonas, and Exiguobacterium in order of decreasing relative abundance. In NCR, Lactococcus, Pantoea, Serratia, Acinetobacter, Sphingomonas, Pseudomonas, Chryseobacterium, Stenotrophomonas, and Microbacterium were found in high relative abundance (in decreasing order). Nevertheless, several taxa differed in relative abundance between the two rootworm species (Supplementary Figure 1). Alpha diversity measurements for WCR and NCR were averaged across sites (Table 1).

FIGURE 2

Figure 2. (A) Phylum level stacked bar chart of average relative abundance of bacterial communities from adult Diabrotica virgifera virgifera (western corn rootworm, WCR) and Diabrotica barberi (northern corn rootworm, NCR) collected from across the United States arranged in order of decreasing longitude (west to east). Each bar represents an individual beetle. (B) Venn diagram comparing amplicon sequence variant (ASV) overlap in WCR and NCR microbiomes using rarefied data. Data presented with Wolbachia removed from communities.

TABLE 1

Table 1. Average alpha diversity metrics of western corn rootworm (WCR) and northern corn rootworm (NCR) collected from different locations across the United States.

Biogeographical Impact on Bacterial Community Composition

Local environments have an impact on the composition of bacterial communities within WCR and NCR (Figure 3). We found bacterial communities differed at the city and state level in both host species based on the results of the PERMANOVA with Bray-Curtis and Jaccard distances (Table 2). In pairwise comparisons, WCR collected from New York, Pennsylvania, and Michigan consisted of unique bacterial communities different from all other states. Arizona was significantly different from all other locations except Colorado and Illinois. Texas was also different from Kansas and South Dakota. For samples collected from smaller regional-scale sampling scheme in Colorado and Kansas (n = 14), overall variation between communities was less than at the landscape level (Table 2). Wild NCR collected from Minnesota harbored unique bacterial communities compared to Missouri and North Dakota. We found no differences between North Dakota and Missouri at the city or state level (Supplementary Table 5). Patterns of differential abundance of ASVs across sites were observed (Figure 4).

FIGURE 3

Figure 3. Principal coordinates analysis of bacterial communities in Diabrotica virgifera virgifera (WCR) and Diabrotica barberi (NCR) collected from their natural host plant corn in the wild based on Bray-Curtis distances. Data presented with Wolbachia removed from communities.

TABLE 2

Table 2. Beta diversity metric comparisons between collection site within host species.

FIGURE 4

Figure 4. Heat map of log¹⁰ abundance of the top 125 most abundant ASVs in western corn rootworm (WCR) across sampling locations in the United States. Each row represents a single ASV and each column represents one beetle microbiome collected from one site within the state. Colorado consisted of multiple sampling sites. Family level grouping of ASV are provided on the right of the figure. Data presented with Wolbachia removed from communities.

Pairwise comparisons do not account for geospatial correlation; however, we used mantel tests to investigate correlations between geographic distance and microbiome community dissimilarity between collection locations. The full WCR data set encompassing the continental scale showed no significant correlation between geographic distance and microbiome dissimilarity for Bray-Curtis (p = 0.15) or Jaccard (p = 0.16) distances. At a smaller geospatial scale (~25 km apart), microbiome dissimilarity increased with increasing distance between sample locations (Figure 5A; Bray-Curtis, p = 0.03; Jaccard, p = 0.03). We suspected the historical geographic pattern of invasion may be correlated to microbiome similarity. Specifically, we examined the relationship between distance from Arizona, the most ancestral population sampled (Lombaert et al., 2018) and microbiome dissimilarity. The values from PC1 of the PCoA of WCR based on Jaccard distances revealed a significant correlation (Figure 5B; p = 0.009). This correlation explained 3% of the variation captured by PC1 (Figure 5B; R² = 0.03).

FIGURE 5

Figure 5. Correlations at various spatial scales between (A) Haversine distance and microbiome community dissimilarity (Jaccard) and (B) Haversine distance from ancestral population (Arizona) of Diabrotica virgifera virgifera (WCR) and values from axis 1 of PCoA of microbiome community similarity based on Jaccard distance.

Discussion

In this study, we found bacterial communities of two insects, WCR and NCR, are significantly impacted by location of collection. However, the geospatial patterning of bacterial communities differed at different spatial scales. The factors governing organization of microbial communities associated with hosts vary in strength across time and space. Populations of WCR can disperse between locations, encountering and mixing microbial populations from isolated patches. If selective environmental forces are similar between locations, the host-associated microbiome may be more affected by dispersal and drift. Thus, community similarity would be expected to decrease over increasing geographical space (Finkel et al., 2012). Here, we found, at smaller spatial scales, bacterial community similarity was correlated with geographic distance between WCR populations. Genetic analyses of the same WCR populations studied here revealed a weak but significant pattern of isolation by distance in Colorado (Kim and Sappington, unpublished data). Our results appear to follow this same pattern of decay by distance (Finkel et al., 2012). We conclude that the spatial turnover observed at the local level in WCR microbiomes is a result of dispersal limitation and ecological drift. However, at broad spatial scales, differences in landscape diversity or host genetics may overshadow dispersal limitation and account for most of the differences in observed microbiome communities separated by long distances.

In Drosophila, rapid genomic evolution can occur in response to the microbiome and, at broad spatial scales, differences in microbiome community composition could be driven by local adaptation (Rudman et al., 2019). Previous genetic analysis found a lack of genetic structuring in WCR across most of the United States, presumably because of a lingering lack of genetic equilibrium after the eastward range expansion out of the western Great Plains beginning in the mid-twentieth century (Kim and Sappington, 2005; Flagel et al., 2014). This eastward expansion of WCR (Gray et al., 2009) may also partially explain the lack of distance decay in microbial community composition at the continental scale observed in this study. Populations in New York, Pennsylvania, and Michigan relatively near the eastern front of WCR expansion may exhibit residual effects of founder populations as evidenced by their disparate microbiome communities. Cropping patterns in agricultural landscapes in the eastern United States are more diverse and include less continuously-planted corn than in the central United States, increasing isolation between populations of insects and microbes by hindering dispersal between locations (Onstad et al., 2003). Correlation with distance from a historically relevant origin region (Arizona) for WCR and microbiome composition may reflect an interaction between host genetic structure and isolation by distance. Alternatively, because the agricultural landscapes of the eastern United States differ markedly from those in the Corn Belt, local environmental selection of the microbial species pool could result in significantly different microbial communities found within the insect. It is not well understood what influence the soil microbiome has on adult corn rootworm microbiome composition or which microbes persist through metamorphosis (Ludwick et al., 2019). Feedback between host genetics, the microbial species pool, and dispersal limitation could account for the patterns found here. Further investigation is warranted to elucidate environmental effects (i.e., temperature, landscape composition, and altitude) on host associated microbiomes in rootworms and other insects.

It is possible that differences in microbiome composition derive from management tactics used within a crop field. While the beetles in this study were collected from corn, we do not know what hybrids or management practices were used in the fields of collection. Macro-level changes in diet (i.e., different host plant species) can have pronounced effects on microbiomes across insect species (Jones et al., 2019; Mason et al., 2020), but differences seem to diminish at finer levels of host plant taxonomic resolution (i.e., plant genotypes; Mason et al., 2021). Several corn genotypes harbor unique taxa, but it is unclear how they are affected by their local environment (Favela et al., 2021). Microbiomes of WCR larvae also respond to ingestion of transgenic crops producing B. thuringiensis (Bt). Larvae resistant to Bt host a unique microbial community that is less taxonomically rich than susceptible insects (Paddock et al., 2021). We do not know if any populations that we sampled are resistant to Bt. To that end, no studies have investigated whether adult WCR microbiomes also respond to Bt. The functionality of microbial communities in adult WCR are documented to involve adaptation to plant defenses and to influence oviposition sites (Lance, 1992; Chu et al., 2013). However, estimation of individual species functionality in the microbiome of WCR is lacking. Resolution provided by 16S rRNA sequencing does not allow for accurate prediction of functionality and requires further study. While we observed differences between Colorado and other sites, we cannot distinguish the biological signal from a possible kit artifact. In addition, we did not examine microbial strain-level diversity in this study. The increased resolution may illuminate differences in community composition that could not be detected by 16S rRNA sequencing.

Differences in the ecological properties and requirements of microbes may generate variation in the species pool at different locations, which can then affect assembly of host-associated microbiomes. Sample location affects microbiome composition in both species of corn rootworm (Table 1). The local species pool of microbes available to colonize the host likely varies with collection location. For instance, communities within a state are more similar than across states (Table 1), which suggests a local environmental or geospatial effect on bacterial assembly. Microbes associated with a host’s food comprise most of the local species pool for that host (Delsuc et al., 2014; Deb et al., 2019). Thus, specialist insects might have a more predictable microbiome (Gomes et al., 2020). Certain components of corn have a predictable and heritable microbiome (Walters et al., 2018). In our case, WCR adults are likely closely associated with the microorganisms colonizing their diet, sharing taxa found both in corn silk and leaves (Wagner et al., 2020; Khalaf et al., 2021). The most common family found in our samples, Enterobacteriaceae, is generally found in other insects feeding on corn plants as well (Jones et al., 2019), suggesting a close association with corn. Other microbial taxa may be more responsive to abiotic and biotic factors associated with the local environment or may vary in dispersal capacity. Future studies examining the site-specific environmental pool of bacteria would improve the understanding of the impact of food source on microbiome composition in insects.

We found WCR and NCR harbor a small set of commonly occurring bacteria regardless of location. Commonly occurring bacteria (i.e., those represented in 50% of samples) totaled seven for both D. virgifera virgifera and D. barberi, with four being shared between the two species. This is consistent with other insects which host facultative microbes with high variation between individuals (Blankenchip et al., 2018; Hammer et al., 2020; Silver et al., 2021). The number of shared taxa between WCR and NCR was relatively small (~12%) but accounted for a high amount of the total sequences within the communities (~90%). Specific filtering by host species can structure distinct microbial communities (Jones et al., 2019; Adair et al., 2020), and while we documented distinct taxa within each host species, we cannot be certain they are due to host filtering. A more controlled study is necessary to distinguish environmental and host effects. Still, host species identity does not account for all the interindividual variation observed in host associated microbiomes (Colman et al., 2012; Yun et al., 2014). In our study, the commonly occurring bacteria may capture a large amount of the variation between sites. Consequently, the correlations between locations we observed may be signatures of smaller parts of the bacterial communities found within WCR and NCR, given the effect sizes (Figures 5A,B). Generalist insects may exhibit weaker distance decay patterns due to variation in diet overshadowing local environmental influence.

Integrating landscape ecology, biogeography, and metacommunity theory into host-associated microbiome studies has been complicated. What microbes an animal encounters in its life are influenced in part by the host diet and the environment in which the host interacts (Littleford-Colquhoun et al., 2022). Both the diet and environment are impacted by selective and neutral forces. Host factors such as species identity or life stage also influence microbial community composition (Suárez-Moo et al., 2020). Specialist insects provide a unique opportunity to investigate factors affecting microbiome composition across geographical space. Dispersal limitation between geographically isolated environments results in unique microbial communities. This isolation can also impact host genetic factors that results in variation between microbiomes. Through a metacommunity framework, WCR can be thought of as “islands within islands,” where microorganisms can disperse and establish (or not) on sessile corn plants whereby they interact with another layer of dispersing organisms (insect) that impart their own selective filtering and dispersal barriers onto the microbial communities. This nested structure accounts for variation in dispersal capacity across scales. Understanding how all determinants interact to shape microbiomes within animals is important to better leverage the beneficial effects microbiomes can provide.

Data Availability Statement

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/Supplementary Material.

Author Contributions

All authors conceived of the study and designed the experiments. KP performed the experiments, conducted statistical analysis, and wrote the first draft of the manuscript. KP and DF interpreted the results. All authors contributed to the article and approved the submitted version.

Conflict of Interest

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.

Publisher’s Note

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.

Acknowledgments

The authors would like to thank everyone who collected and contributed beetle samples to this work. This includes Chris Difonzo, Elson Shields, Tony Testa, Joseph Spencer, John Tooker, Billy Fuller, Brad McManus, Lance Meinke, Andy Michel, Lucinda Wallace, Patrick Porter, Sarah Zukoff, Anthony Zukoff, Dan Mollenbeck, Janet Knodel, and Veronica Calles-Torres. We thank Randy Ritland and Anthony Bennett of USDA-ARS in Ames, Iowa, who helped with sample collection in Colorado and DNA extraction. We would like to acknowledge MU Metagenomics Center for their assistance in preparing samples for sequencing, MU DNA Core for performing the sequencing, and MU Informatics Core for their help with upstream processing of sequencing results. Funding for this project was provided by the USDA-ARS and the University of Missouri. Mention of a proprietary product does not constitute an endorsement or recommendation for its use by the USDA or the University of Missouri. USDA-ARS and the University of Missouri are equal opportunity employers.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.898744/full#supplementary-material

Footnotes

1. ^https://www.qiagen.com/us/resources/resourcedetail?id=00e4513c-597b-4bd5-a600-9259e6d62d07&lang=en

2. ^https://www.qiagen.com/us/resources/resourcedetail?id=68f29296-5a9f-40fa-8b3d-1c148d0b3030&lang=en

3. ^http://journal.embnet.org/index.php/embnetjournal/article/view/200/479

References

Adair, K. L., Bost, A., Bueno, E., Kaunisto, S., Kortet, R., Peters-Schulze, G., et al. (2020). Host determinants of among-species variation in microbiome composition in drosophilid flies. ISME J. 14, 217–229. doi: 10.1038/s41396-019-0532-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Adair, K. L., and Douglas, A. E. (2017). Making a microbiome: the many determinants of host-associated microbial community composition. Curr. Opin. Microbiol. 35, 23–29. doi: 10.1016/j.mib.2016.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Adair, K. L., Wilson, M., Bost, A., and Douglas, A. E. (2018). Microbial community assembly in wild populations of the fruit fly Drosophila melanogaster. ISME J. 12, 959–972. doi: 10.1038/s41396-017-0020-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ben-Yosef, M., Pasternak, Z., Jurkevitch, E., and Yuval, B. (2014). Symbiotic bacteria enable olive flies (B actrocera oleae) to exploit intractable sources of nitrogen. J Evol. Biol. 27, 2695–2705. doi: 10.1111/jeb.12527

PubMed Abstract | CrossRef Full Text | Google Scholar

Blankenchip, C. L., Michels, D. E., Elizabeth Braker, H., and Goffredi, S. K. (2018). Diet breadth and exploitation of exotic plants shift the core microbiome of Cephaloleia, a group of tropical herbivorous beetles. PeerJ 6:e4793. doi: 10.7717/peerj.4793

PubMed Abstract | CrossRef Full Text | Google Scholar

Brochet, S., Quinn, A., Mars, R. A. T., Neuschwander, N., Sauer, U., and Engel, P. (2021). Niche partitioning facilitates coexistence of closely related gut bacteria. elife 10:e68583. doi: 10.7554/eLife.68583

PubMed Abstract | CrossRef Full Text | Google Scholar

Caccia, S., Casartelli, M., and Tettamanti, G. (2019). The amazing complexity of insect midgut cells: types, peculiarities, and functions. Cell Tissue Res. 377, 505–525. doi: 10.1007/s00441-019-03076-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Caccia, S., Di Lelio, I., La Storia, A., Marinelli, A., Varricchio, P., Franzetti, E., et al. (2016). Midgut microbiota and host immunocompetence underlie Bacillus thuringiensis killing mechanism. Proc. Natl. Acad. Sci. U. S. A. 113, 9486–9491. doi: 10.1073/pnas.1521741113

PubMed Abstract | CrossRef Full Text | Google Scholar

Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869

PubMed Abstract | CrossRef Full Text | Google Scholar

Calles-Torrez, V., Knodel, J. J., Boetel, M. A., French, B. W., Fuller, B. W., and Ransom, J. K. (2019). Field-evolved resistance of northern and western corn rootworm (Coleoptera: Chrysomelidae) populations to corn hybrids expressing single and pyramided Cry3Bb1 and Cry34/35Ab1 Bt proteins in North Dakota. J. Econ. Entomol. 112, 1875–1886. doi: 10.1093/jee/toz111

PubMed Abstract | CrossRef Full Text | Google Scholar

Chu, C. C., Spencer, J. L., Curzi, M. J., Zavala, J. A., and Seufferheld, M. J. (2013). Gut bacteria facilitate adaptation to crop rotation in the western corn rootworm. Proc. Natl. Acad. Sci. U. S. A. 110, 11917–11922. doi: 10.1073/pnas.1301886110

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohen, H., McFrederick, Q. S., and Philpott, S. M. (2020). Environment shapes the microbiome of the blue orchard bee, Osmia lignaria. Microb. Ecol. 80, 897–907. doi: 10.1007/s00248-020-01549-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Colman, D. R., Toolson, E. C., and Takacs-Vesbach, C. D. (2012). Do diet and taxonomy influence insect gut bacterial communities? Mol. Ecol. 21, 5124–5137. doi: 10.1111/j.1365-294X.2012.05752.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Deb, R., Nair, A., and Agashe, D. (2019). Host dietary specialization and neutral assembly shape gut bacterial communities of wild dragonflies. PeerJ 7:e8058. doi: 10.7717/peerj.8058

PubMed Abstract | CrossRef Full Text | Google Scholar

Delsuc, F., Metcalf, J. L., Wegener Parfrey, L., Song, S. J., González, A., and Knight, R. (2014). Convergence of gut microbiomes in myrmecophagous mammals. Mol. Ecol. 23, 1301–1317. doi: 10.1111/mec.12501

PubMed Abstract | CrossRef Full Text | Google Scholar

Favela, A., Bohn, M. O., and Kent, A. D. (2021). Maize germplasm chronosequence shows crop breeding history impacts recruitment of the rhizosphere microbiome. ISME J. 15, 2454–2464. doi: 10.1038/s41396-021-00923-z

CrossRef Full Text | Google Scholar

Finkel, O. M., Burch, A. Y., Elad, T., Huse, S. M., Lindow, S. E., Post, A. F., et al. (2012). Distance-decay relationships partially determine diversity patterns of phyllosphere bacteria on Tamrix trees across the Sonoran Desert. Appl. Environ. Microbiol. 78, 6187–6193. doi: 10.1128/AEM.00888-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Flagel, L. E., Bansal, R., Kerstetter, R. A., Chen, M., Carroll, M., Flannagan, R., et al. (2014). Western corn rootworm (Diabrotica virgifera virgifera) transcriptome assembly and genomic analysis of population structure. BMC Genomics 15, 1–13. doi: 10.1186/1471-2164-15-195

CrossRef Full Text | Google Scholar

Giordano, R., Jackson, J. J., and Robertson, H. M. (1997). The role of Wolbachia bacteria in reproductive incompatibilities and hybrid zones of Diabrotica beetles and Gryllus crickets. Proc. Natl. Acad. Sci. U. S. A. 94, 11439–11444. doi: 10.1073/pnas.94.21.11439

PubMed Abstract | CrossRef Full Text | Google Scholar

Gomes, S. I., Kielak, A. M., Hannula, S. E., Heinen, R., Jongen, R., Keesmaat, I., et al. (2020). Microbiomes of a specialist caterpillar are consistent across different habitats but also resemble the local soil microbial communities. Anim. Microbiome 2, 1–12. doi: 10.1186/s42523-020-00055-3

CrossRef Full Text | Google Scholar

Gray, M. E., Sappington, T. W., Miller, N. J., Moeser, J., and Bohn, M. O. (2009). Adaptation and invasiveness of western corn rootworm: intensifying research on a worsening pest. Annu. Rev. Entomol. 54, 303–321. doi: 10.1146/annurev.ento.54.110807.090434

PubMed Abstract | CrossRef Full Text | Google Scholar

Hammer, T. J., Dickerson, J. C., McMillan, W. O., and Fierer, N. (2020). Heliconius butterflies host characteristic and phylogenetically structured adult-stage microbiomes. Appl. Environ. Microbiol. 86, e02007–e02020. doi: 10.1128/AEM.02007-20

CrossRef Full Text | Google Scholar

Huang, K., Wang, J., Huang, J., Zhang, S., Vogler, A. P., Liu, Q., et al. (2021). Host phylogeny and diet shape gut microbial communities within bamboo-feeding insects. Front. Microbiol. 12:633075. doi: 10.3389/fmicb.2021.633075

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, A. G., Mason, C. J., Felton, G. W., and Hoover, K. (2019). Host plant and population source drive diversity of microbial gut communities in two polyphagous insects. Sci. Rep. 9, 1–11. doi: 10.1038/s41598-019-39163-9

CrossRef Full Text | Google Scholar

Khalaf, E. M., Shrestha, A., Rinne, J., Lynch, M. D. J., Shearer, C. R., Limay-Rios, V., et al. (2021). Transmitting silks of maize have a complex and dynamic microbiome. Sci. Rep. 11, 1–17. doi: 10.1038/s41598-021-92648-4

CrossRef Full Text | Google Scholar

Kim, K. S., and Sappington, T. W. (2005). Genetic structuring of western corn rootworm (Coleoptera: Chrysomelidae) populations in the United States based on microsatellite loci analysis. Environ. Entomol. 34, 494–503. doi: 10.1603/0046-225X-34.2.494

CrossRef Full Text | Google Scholar

Krysan, J. L., and Miller, T. A., (eds). (1986). Methods for the Study of Pest Diabrotica. New York: Springer US.

Google Scholar

Lance, D. R. (1992). Odors influence choice of oviposition sites by Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae). J. Chem. Ecol. 18, 1227–1237. doi: 10.1007/BF00980076

PubMed Abstract | CrossRef Full Text | Google Scholar

Levine, E., Oloumi-Sadeghi, H., and Fisher, J. R. (1992). Discovery of multiyear diapause in Illinois and South Dakota northern corn rootworm (Coleoptera: Chrysomelidae) eggs and incidence of the prolonged diapause trait in Illinois. J. Econ. Entomol. 85, 262–267. doi: 10.1093/jee/85.1.262

CrossRef Full Text | Google Scholar

Levine, E., Spencer, J. L., Isard, S. A., Onstad, D. W., and Gray, M. E. (2002). Adaptation of the western corn rootworm to crop rotation: evolution of a new strain in response to a management practice. Am. Entomol. 48, 94–107. doi: 10.1093/ae/48.2.94

CrossRef Full Text | Google Scholar

Littleford-Colquhoun, B. L., Weyrich, L. S., Hohwieler, K., Cristescu, R., and Frère, C. H. (2022). How microbiomes can help inform conservation: landscape characteristics of gut microbiota helps shed light on additional population structure in specialist folivore. Anim. Microbiome 4:12. doi: 10.1186/s42523-021-00122-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Lombaert, E., Ciosi, M., Miller, N. J., Sappington, T. W., Blin, A., and Guillemaud, T. (2018). Colonization history of the western corn rootworm (Diabrotica virgifera virgifera) in North America: insights from random forest ABC using microsatellite data. Biol. Invasions 20, 665–677. doi: 10.1007/s10530-017-1566-2

CrossRef Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15:550. doi: 10.1186/s13059-014-0550-8

CrossRef Full Text | Google Scholar

Ludwick, D. C., Ericsson, A. C., Meihls, L. N., Gregory, M. L. J., Finke, D. L., Coudron, T. A., et al. (2019). Survey of bacteria associated with western corn rootworm life stages reveals no difference between insects reared in different soils. Sci. Rep. 9, 1–11. doi: 10.1038/s41598-019-51870-x

CrossRef Full Text | Google Scholar

Martinez Arbizu, P. (2020). pairwise Adonis: Pairwise multilevel comparison using adonis.

Google Scholar

Mason, C. J., Clair, A. S., Peiffer, M., Gomez, E., Jones, A. G., Felton, G. W., et al. (2020). Diet influences proliferation and stability of gut bacterial populations in herbivorous lepidopteran larvae. PLoS One 15:e0229848. doi: 10.1371/journal.pone.0229848

PubMed Abstract | CrossRef Full Text | Google Scholar

Mason, C. J., Hoover, K., and Felton, G. W. (2021). Effects of maize (Zea mays) genotypes and microbial sources in shaping fall armyworm (Spodoptera frugiperda) gut bacterial communities. Sci. Rep. 11, 1–10. doi: 10.1038/s41598-021-83497-2

CrossRef Full Text | Google Scholar

Mason, C. J., Ray, S., Shikano, I., Peiffer, M., Jones, A. G., Luthe, D. S., et al. (2019). Plant defenses interact with insect enteric bacteria by initiating a leaky gut syndrome. Proc. Natl. Acad. Sci. U. S. A. 116, 15991–15996. doi: 10.1073/pnas.1908748116

PubMed Abstract | CrossRef Full Text | Google Scholar

McMurdie, P. J., and Holmes, S. (2013). Phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, D. L., Smith, E. A., and Newton, I. L. G. (2021). A bacterial symbiont protects honey bees from fungal disease. mBio 12:e0050321. doi: 10.1128/mBio.00503-21

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, E. T., Svanbäck, R., and Bohannan, B. J. M. (2018). Microbiomes as metacommunities: Understanding host-associated microbes through metacommunity ecology. Trends Ecol. Evol. 33, 926–935. doi: 10.1016/j.tree.2018.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Moeller, A. H., Suzuki, T. A., Lin, D., Lacey, E. A., Wasser, S. K., and Nachman, M. W. (2017). Dispersal limitation promotes the diversification of the mammalian gut microbiota. Proc. Natl. Acad. Sci. U. S. A. 114, 13768–13773. doi: 10.1073/pnas.1700122114

PubMed Abstract | CrossRef Full Text | Google Scholar

Onstad, D. W., Crowder, D. W., Isard, S. A., Levine, E., Spencer, J. L., O’neal, M. E., et al. (2003). Does landscape diversity slow the spread of rotation-resistant western corn rootworm (Coleoptera: Chrysomelidae)? Environ. Entomol. 32, 992–1001. doi: 10.1603/0046-225X-32.5.992

CrossRef Full Text | Google Scholar

Paddock, K. J., Pereira, A. E., Finke, D. L., Ericsson, A. C., Hibbard, B. E., and Shelby, K. S. (2021). Host resistance to Bacillus thuringiensis is linked to altered bacterial community within a specialist insect herbivore. Mol. Ecol. 30, 5438–5453. doi: 10.1111/mec.15875

PubMed Abstract | CrossRef Full Text | Google Scholar

Parajuli, A., Hui, N., Puhakka, R., Oikarinen, S., Grönroos, M., Selonen, V. A. O., et al. (2020). Yard vegetation is associated with gut microbiota composition. Sci. Total Environ. 713:136707. doi: 10.1016/j.scitotenv.2020.136707

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, R., Dzialo, M. C., Spaepen, S., Nsabimana, D., Gielens, K., Devriese, H., et al. (2019). Microbial communities of the house fly Musca domestica vary with geographical location and habitat. Microbiome 7, 1–12. doi: 10.1186/s40168-019-0748-9

CrossRef Full Text | Google Scholar

Robusto, C. C. (1957). The Cosine-Haversine formula. Am. Math. Mon. 64, 38–40. doi: 10.2307/2309088

CrossRef Full Text | Google Scholar

Roehrdanz, R. L., and Levine, E. (2007). Wolbachia bacterial infections linked to mitochondrial DNA reproductive isolation among populations of northern corn rootworm (Coleoptera: Chrysomelidae). Ann. Entomol. Soc. Am. 100, 522–531. doi: 10.1603/0013-8746(2007)100[522:WBILTM]2.0.CO;2

CrossRef Full Text | Google Scholar

Rudman, S. M., Greenblum, S., Hughes, R. C., Rajpurohit, S., Kiratli, O., Lowder, D. B., et al. (2019). Microbiome composition shapes rapid genomic adaptation of Drosophila melanogaster. Proc. Natl. Acad. Sci. U. S. A. 116, 20025–20032. doi: 10.1073/pnas.1907787116

PubMed Abstract | CrossRef Full Text | Google Scholar

Silver, A., Perez, S., Gee, M., Xu, B., Garg, S., Will, K., et al. (2021). Persistence of the ground beetle (Coleoptera: Carabidae) microbiome to diet manipulation. PLoS One 16:e0241529. doi: 10.1371/journal.pone.0241529

PubMed Abstract | CrossRef Full Text | Google Scholar

Suárez-Moo, P., Cruz-Rosales, M., Ibarra-Laclette, E., Desgarennes, D., Huerta, C., and Lamelas, A. (2020). Diversity and composition of the gut microbiota in the developmental stages of the dung beetle Copris incertus Say (Coleoptera, Scarabaeidae). Front. Microbiol. 11:1698. doi: 10.3389/fmicb.2020.01698

PubMed Abstract | CrossRef Full Text | Google Scholar

Wagner, M. R., Busby, P. E., and Balint-Kurti, P. (2020). Analysis of leaf microbiome composition of near-isogenic maize lines differing in broad-spectrum disease resistance. New Phytol. 225, 2152–2165. doi: 10.1111/nph.16284

PubMed Abstract | CrossRef Full Text | Google Scholar

Walters, W. A., Jin, Z., Youngblut, N., Wallace, J. G., Sutter, J., Zhang, W., et al. (2018). Large-scale replicated field study of maize rhizosphere identifies heritable microbes. Proc. Natl. Acad. Sci. U. S. A. 115, 7368–7373. doi: 10.1073/pnas.1800918115

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Kapun, M., Waidele, L., Kuenzel, S., Bergland, A. O., and Staubach, F. (2020). Common structuring principles of the Drosophila melanogaster microbiome on a continental scale and between host and substrate. Environ. Microbiol. Rep. 12, 220–228. doi: 10.1111/1758-2229.12826

PubMed Abstract | CrossRef Full Text | Google Scholar

Wechsler, S., and Smith, D. (2018). Has resistance taken root in US corn fields? Demand for insect control. Am. J. Agric. Econ. 100, 1136–1150. doi: 10.1093/ajae/aay016

CrossRef Full Text | Google Scholar

Yatsunenko, T., Rey, F. E., Manary, M. J., Trehan, I., Dominguez-Bello, M. G., Contreras, M., et al. (2012). Human gut microbiome viewed across age and geography. Nature 486, 222–227. doi: 10.1038/nature11053

PubMed Abstract | CrossRef Full Text | Google Scholar

Yun, J. H., Roh, S. W., Whon, T. W., Jung, M. J., Kim, M. S., Park, D. S., et al. (2014). Insect gut bacterial diversity determined by environmental habitat, diet, developmental stage, and phylogeny of host. Appl. Environ. Microbiol. 80, 5254–5264. doi: 10.1128/AEM.01226-14

PubMed Abstract | CrossRef Full Text | Google Scholar