Article What Can the Bacterial Community of sexdens (Linnaeus, 1758) Tell Us about the Habitats in Which This Evolves?

Manuela de Oliveira Ramalho 1,2,*, Cintia Martins 3, Maria Santina Castro Morini 4 and Odair Correa Bueno 1

1 Universidade Estadual Paulista, UNESP, Centro de Estudos de Insetos Sociais—CEIS, Instituto de Biociências, Campus Rio Claro, Avenida 24A, 1515, Bela Vista, Rio Claro, SP 13506-900, Brazil; [email protected] 2 Cornell University, Department of Entomology, 129 Garden Ave, Ithaca, NY, 14850, USA 3 Universidade Federal do Piauí - Campus Ministro Reis Velloso, Av. São Sebastião, 2819, Parnaíba 64.202-020, Piauí, Brazil; [email protected] 4 Universidade de Mogi das Cruzes, Núcleo de Ciências Ambientais, Av. Dr. Cândido Xavier de Almeida e Souza, 200, Centro Cívico, Mogi das Cruzes, SP 08780-911, Brazil; [email protected] * Correspondence: [email protected]

Received: 5 March 2020; Accepted: 22 May 2020; Published: 28 May 2020

Abstract: Studies of bacterial communities can reveal the evolutionary significance of symbiotic interactions between hosts and their associated , as well as identify environmental factors that may influence host biology. Atta sexdens is an ant species native to Brazil that can act as an agricultural pest due to its intense behavior of cutting plants. Despite being extensively studied, certain aspects of the general biology of this species remain unclear, such as the evolutionary implications of the symbiotic relationships it forms with bacteria. Using high-throughput amplicon sequencing of 16S rRNA genes, we compared for the first time the bacterial community of A. sexdens (whole ant workers) populations according to the habitat (natural versus agricultural) and geographical location. Our results revealed that the bacterial community associated with A. sexdens is mainly influenced by the geographical location, and secondarily by the differences in habitat. Also, the bacterial community associated with citrus differed significantly from the other communities due to the presence of Tsukamurella. In conclusion, our study suggests that environmental shifts may influence the bacterial diversity found in A. sexdens.

Keywords: environment; extensive cultivation; bacterial diversity; sustainable development; conservation

1. Introduction Several studies have focused on understanding the factors that influence bacterial communities and their implications for the host, ranging from [1] to [2–8]. Previous studies on vertebrate and invertebrate , including , have shown that bacterial communities are influenced by different factors such as host phylogeny [5,9,10], stages of development [4,11] and host environment [3]. In addition, the diet [2,12] and age [13] are other important host attributes that should be considered in studies on associated bacterial communities. However, in order to fully understand the ecological patterns and mechanisms behind these naturally established symbioses, more studies exploring new hosts are needed.

Insects 2020, 11, 332; doi:10.3390/insects11060332 www.mdpi.com/journal/insects Insects 2020, 11, 332 2 of 21

The composition and ecological drivers of host bacterial communities have been extensively studied in ants of the genera Camponotus and Cephalotes [2,4,9,10,14,15]. These symbioses can be so intimate and fundamental to the participating organisms that in many cases, hosts cannot survive without their symbionts, in this case called the primary. Examples include Blochmannia associated with Camponotus species [3,4,16–18], and Rhizobiales with Cephalotes species [2]. However, because the Fomicidae family is very diverse on Earth and includes a large number of species with very different lifestyles [19], it is important to extend the study of the associated microbiota to other species [20]. Exclusively distributed in the Neotropical and Neartic regions, the Atta has 17 valid species and is most diverse in Brazil, where these ants are popularly known as leaf-cutting ants, with nine species recorded [19,21]. The evolutionary success of the genus is largely attributed to the mutualistic relationship with the Leucogaricus gongylophorus and other associated organisms such as bacteria [22–27]. Leaf-cutter ant workers cut leaves and other pieces of plant material and use them as a substrate for the fungus symbiont, causing great economic damage mainly in extensive agricultural habitats [26,28]. There are estimates that the economic damage caused by leaf-cutting ants reach a scale of billions of dollars in the global economy [29]. Most control techniques have temporary effects and, in some cases, otherwise negatively impact the environment, as is the case for the chemical method that is most frequently utilized by farmers [30,31]. Among the leaf-cutting ant species, Atta sexdens (Linnaeus, 1758) (Formicidae: ) stands out in terms of abundance and wide distribution in Brazil [32] and the intimate association of leaf-cutting ants with human establishes this genus as an important group of ants to be studied. However, there are only a few studies about the ecological factors that can affect host bacterial communities in any ant genus [1,3,4,10,33,34], and even fewer that address these questions in Atta [23,35,36]. São Paulo state (Brazil) hosts a great diversity of natural and agricultural habitats in close proximity to each other, which makes it the ideal place to study the effect of these different habitats on the bacterial community of Atta sexdens. The Cerrado (Brazilian Savannah) and Atlantic Forest are the dominant natural habitats in the São Paulo state, but they are also heavily affected by agricultural practices. The Cerrado is characterized by trees with contorted trunks ranging from evergreen to deciduous in the dry season and large shrubs of about 2–8 meters in height, and a grassy ground [37]. The remaining natural coverage of the Cerrado in the São Paulo state is about 13% of its original (which corresponded to 33% of the state), highlighting the intense anthropogenic impact on its coverage in the state [38]. The Atlantic Forest is characterized by being an Ombrophilous Dense Forest [39] with several endemic species of flora and fauna. Historically, this habitat covered an extensive part of the Brazilian coast, and today it is reduced to just 12% of the original size [40]. In São Paulo state, the Atlantic Forest is reduced and now its remnants cover around 16.3% of the state [41]. The main agricultural habitats in the São Paulo state are sugar cane, followed in smaller scale by citrus and eucalyptus reforestation [42]. These activities use large quantities of pesticides, and fipronil and sulfluramide are specifically used for the control of leaf-cutting ants [28]. With the shift from natural to agricultural habitat, simplification of the landscape can occur due to floral simplification and application of insecticides. These changes can be beneficial to some insects and harmful to others [43]. In the case of pollinators, it was shown that the loss of environmental diversity caused by the agricultural process impacts on the fitness and success of the insects’ survival [44]. This fitness loss in agricultural habitats can be attributed to changes in the microbiota that can compromise the health of the insect, as has been proven in some studies of bees [45,46]. Therefore, associated bacterial communities can be modified with environmental shifts [47], and anthropogenic modifications can have an impact on the composition of the microbiota [48]. However, a study has never been conducted regarding the impact of changes in habitat directly affecting the microbiota associated with Atta sexdens. We hypothesize that differences in habitat (natural and agricultural) can influence the diversity of the bacterial community associated with leaf-cutting ants. Furthermore, analyzing the factors that alter the bacterial community can provide insights into understanding the general biology of the host. Thus, in this study we tested whether 1) natural (Cerrado and Atlantic Forest) and agricultural (sugar

Insects 2020, 11, 332 3 of 21 cane, citrus and eucalyptus) habitats, and 2) geographical location (West, Center, East, South and North, from São Paulo state in Brazil) can influence the A. sexdens microbiota. We expect all these factors to influence and explain the diversity of bacterial communities associated with this leaf- cutting ant.

2. Materials and Methods

2.1. Sample Collection Samples were collected in natural habitats (n = 4 areas of Cerrado, n = 5 areas of Atlantic Forest) and agricultural habitats (n = 5 sugar cane; n = 5 citrus; n = 6 eucalyptus), located in São Paulo state – Brazil (Figure S1). In each habitat, a linear transect was traversed along which foraging workers (with average size and weight between 10–20 mg and head width between 2.0 and 2.8 mm) were collected. The distance between colonies (n = 3 for each location) was at least 200 meters. We collected 12 workers in the Cerrado and 15 in the Atlantic Forest, 15 each in the sugar cane and citrus agricultural habitats, and 18 in the eucalyptus agricultural habitats (Figure 1; more details in Table S1). It is important to highlight that all agricultural habitats from this study are managed with agrochemicals: insecticides fungicides and herbicides for weed control. Synthetic were also applied as part of cultural practices [42]. The taxonomic identifications were determined following morphological characters for species identification available in Mackay and Mackay (1986) [49] and vouchers were deposited in the collection de Mirmecologia do Alto Tietê. Samples used were collected in 95% ethanol in the field and kept at −20 °C until DNA extraction.

Figure 1. Samples of Atta sexdens collected from different natural and agricultural habitats.

Insects 2020, 11, 332 4 of 21

2.2. DNA Extraction and Bacterial DNA Sequencing In total, we used 75 whole ant workers for DNA extraction using the DNeasy Power Kit (Qiagen, Germantown, MD, USA) following recommended methods from Moreau [50] to avoid contamination. Furthermore, three negative controls were added to the extraction methods and to subsequent 16S rRNA amplifications and library sequencing. These negative controls are blank samples (corresponding to water instead of DNA in the polymerase chain reaction (PCR)) that were used to remove contaminants that may have been added in the process of obtaining the data. Amplification of the V4 region of 16S rRNA (515F/806R primers) was performed in triplicate following Earth Microbiome Project (EMP) protocols (http://www.earthmicrobiome.org/emp- standard-protocols/), also described in Caporaso [51]. PCR reactions followed recommendations from the Taq Platinum Kit (Invitrogen, Schwerte, Germany) with 2.5 μL Buffer 10×, 1.0 μL MgCl2 50 mM, 1.0 μL dNTPs 10 Mm, 1.0 μL BSA (bovine serum albumin) 1 mg/mL, 0.5 μL 515F primer 10 uM, 0.5 μL 806R primer 10 uM, 0.5 μL Taq Platinum 5U/μL, 5 μL template DNA (10 ng/μL), 13.0 μL PCR Water (Certified DNA-free). The reaction was submitted to the thermal cycler under the following conditions: 95 °C for 3 min, 30 cycles at 95 °C for 30 s, 60 °C for 30 s and 72 °C for 30 s, with final of 72 °C for 10 min and 4 °C hold. Purification was performed with the Ampure X beads and adapter connections (NEXTERA XT) were also performed following the manufacturer's recommendations. The samples were quantified via Nanodrop (Thermo Fisher Scientific, Waltham, MA, USA). Pooling and normalization were performed so that the samples had the same concentration. Then, the DNA pool was diluted down to 2 nM, denatured, and diluted to a final concentration of 16 pM with a 20% PhiX for sequencing on the Illumina MiSeq with run 2× 250 pb at the Centro de Genômica (ESALQ, Piracicaba, SP, Brazil). All raw sequence data are publicly available in NCBI SRA accession number PRJNA556881 and study SAMN12373152.

2.3. Bioinformatic Analysis Data were imported into Qiime2 [52] following the Casava 1.8 paired-end demultiplexed fastq protocol. The Dada2 plugin was used to create a feature table and quality control [53,54]. Paired-end sequence reads were trimmed (forward reads trim 30 and trunc 245, reverse reads trim 30 and trunc 220) for better results and SILVA_132_QIIME database [55,56] with 99% similarity was used to access ASVs (amplicon sequence variants) for the taxonomic information from the data sequenced. We trained our specific classifier based on our sample preparation and sequencing, using the “feature- classifier classify-sklearn" command [57]. Negative controls were filtered by the prevalence method using the Decontam package in R software (R Development Core Team, Vienna, Austria, 2019) following default settings [58]. Once our feature table was filtered, an extra filtering step was applied in Qiime2 in order to remove mitochondrial, chloroplast, and low abundance ASVs. Reads were aligned with the “align-to-tree- mafft-fasttree” command [59] and the phylogeny was inferred before alpha and beta diversity analyses [60,61]. Samples were normalized at 2000 reads in subsequent analyses. The categories used in order to test for potential influences in leaf-cutting ants were conducted in two classes: 1) natural habitats (Cerrado and Atlantic Forest) and main agricultural habitats (sugar cane, citrus and eucalyptus), and 2) geographical location where the colonies were collected. To test the influence of geographical location within São Paulo State (Brazil), we divided the data according to their cardinal points of collection: West, Center, East, South, and North of the state and samples were assigned to those locations independent of their status belonging to natural or agricultural habitats (Figure S1). Permutational multivariate analysis of variance (permanova) tests with the “p- pairwise” parameter were conducted with Unifrac distance, a qualitative measure that includes phylogenetic relationships between features [62]. Analysis of non-metric multidimensional scaling (NMDS), principle coordinates analysis (PCoA), and related statistics were implemented to better visualize our results of relationships between ecological communities and in the PAST3 software package [63–65] and Qiime2 in Emperor software [61]. Similarity percentage (SIMPER) analysis [63] was used to verify the ASVs’ contribution to structuring the bacterial community found in the study.

Insects 2020, 11, 332 5 of 21

Network Analysis and heatmap plots of bacterial communities were conducted in R software [66] with the packages phyloseq [67] and ggplot2 [68].

2.4. Data Accessibility All raw sequence data are publicly available in NCBI SRA accession number PRJNA556881 and study SAMN12373152.

3. Results Samples were normalized to a minimum of 2000 reads per sample with sequencing reaching the expected depth (Figure S2). After this cutoff, 75 samples remained and represented a total of 6,383,699 reads distributed in 273 different ASVs. Sampling ranged from a maximum of 177,158 reads to a minimum of 2244 reads, confirming the successful sequencing of the bacterial community associated with A. sexdens. Subsequent analyses were performed at the family taxonomic level of the bacterial communities or at the high hierarchy taxonomic level when available in the database. The ASV table and the overview of relative abundance can be seen in Table S2 and Figure S3, respectively. The main bacterial community recovered revealed Mesoplasma as the predominant bacterial genus (7%), followed by Tsukamurella (6%), Aeromicrobium (Nocardioidaceae) (5%) and others in smaller quantities. Overall when we include all natural and agricultural habitats, the bacterial community of A. sexdens consists of several taxa per sample (Figure 2), which makes it difficult to detect any bacterial structure or community patterns.

Insects 2020, 11, 332 6 of 21

Figure 2. Summary of main microbiome amplicon sequence variants (ASVs) found in Atta sexdens with 16S rRNA amplicon sequencing. Alpha diversity of microbiota in each sample of A. sexdens. Samples were grouped according to the habitat (natural or agricultural) to which they belong. Bar graphs for each library (one column = community from a sample) show the abundance of sequence reads classified to selected 99% similarity. Each color represents a distinct bacterium. ** indicates samples in natural habitat.

3.1. Beta Diversity and the Influence of Natural and Agricultural Habitats The results of the beta diversity analysis suggest that the bacterial community of leaf-cutting ants do not differ between natural habitats and agricultural habitats (permanova weighted unifrac distance: p-value = 0.177, pseudo F = 1.329, Figure 3A; permanova unweighted unifrac distance: p- value = 0.423, pseudo F = 1.014, Figure 3B). When beta diversity analysis was conducted only with samples from natural habitats, we also did not find statistically significant differences (permanova weighted unifrac distance: p-value = 0.599, Pseudo F = 0.858, Figure 3C; permanova unweighted unifrac distance p-value = 0.086, Pseudo F = 1.186, Figure 3D). Overall, there is no difference in the bacterial community of A. sexdens associated with the Cerrado or Atlantic Forest. When analyzing only leaf-cutting ant bacteria associated with agricultural habitats (cane, citrus and eucalyptus) the beta diversity results indicated an influence in the weighted bacterial community

Insects 2020, 11, 332 7 of 21

(permanova weighted unifrac distance: p-value = 0.04, Pseudo F = 1,802) (Figure 3E). Unweighted analysis differences were not significant (permanova unweighted unifrac distance: p-value = 0.151, Pseudo F = 1.095) (Figure 3F), confirming that the abundance of bacteria recovered in weighted unifrac distance are important in helping explain associations with agricultural habitats, but not the bacterial composition.

Figure 3. Non-metric multidimensional scaling (NMDS) plot illustrating bacterial community structure among Atta sexdens workers of natural vs. agricultural habitats. (A). Weighted unifrac distance matrix for samples from natural and agricultural habitats. (B). Unweighted unifrac distance matrix for samples from natural and agricultural habitats. (C). Weighted unifrac distance matrix for samples from Atlantic forest and Cerrado. (D). Unweighted unifrac distance matrix for samples from Atlantic Forest and Cerrado. (E). Weighted unifrac distance matrix for samples from agricultural habitats. (F). Unweighted unifrac distance matrix for samples from agricultural habitats. Dots were colored according to the different groups. The distance between bacterial communities represents their underlying distance in the multivariate space. Statistics highlighted in red show p values greater than 0.05. Besides one significant p value in the weighted unifrac distance for the agricultural habitat, note that there is not a clear structuring of the different groupings.

Insects 2020, 11, 332 8 of 21

To explore this result more deeply, we performed the paired permanova test (weighted unifrac distance) that allows all samples associated with the agricultural habitat to be analyzed in pairs (see Figure 4). The only pair that obtained a significant result was of bacterial communities associated with citrus and eucalyptus agricultural habitats (p-value = 0.04 and Pseudo F = 3.045), contributing to the difference found in Figure 3. In addition, through the results of SIMPER analysis, the main bacteria (ASVs) associated with citrus and eucalyptus agriculture are identified in Table 1 (for all groupings see Table S3). Note that Tsukamurella is responsible for 12.5% of the difference found in this pair (citrus and eucalyptus), followed by Aeromicrobium with 7.07% and others in smaller abundance. Our heatmap analysis gives a better view of Tsukamurella's association with A. sexdens, especially citrus (Figure 5).

Figure 4. Pairwise permutational multivariate analysis of variance (permanova) results of agricultural sites. Results highlighted in red indicate p value greater than 0.05. Bacterial communities associated with citrus and eucalyptus agriculture contributes to the difference found in Figure 3.

Table 1. Similarity percentage (SIMPER) analysis of main bacteria in citrus and eucalyptus group. This test indicates the contribution of specific ASVs to the observed differences in community structure among citrus and eucalyptus. Note that for this group, Tsukamurella is the most influential bacteria.

Overall Average Most Influential Percent Contribuition Habitats Cumulative % Dissimilarity ASV/ to Difference Tsukamurella 12.05 12.05 Aeromicrobium 7.07 19.12 Anthococcus 5.89 25.01 Mesoplasma 5.11 30.12 Micrococcales 4.3 34.42 Citrus vs. 85.00% Alphaproteobacteria 3.91 38.32 Eucalyptus Enterobacteriaceae 2.84 41.17 Nocardioides 2.66 43.83 Pseudomonas 2.53 46.36 Gammaproteobacteria 2.09 48.45 Burkholderiaceae 2.04 50.49

Insects 2020, 11, 332 9 of 21

Figure 5. Main bacterial ASVs found in Atta sexdens are visualized in this heatmap. The shading in the heatmap indicate variation in the relative abundance of different bacteria in all samples from this study, ranging from 0% (white) to 100% (black). Highlighted in blue is the Tsukamurella bacteria, which although found in different habitats, is mainly associated with citrus (blue boxes).

Insects 2020, 11, 332 10 of 21

3.2. Beta Diversity and the Influence of Geographic Location Another factor that may be responsible for influencing the diversity of bacterial communities is the locality where the host was collected. Our results show the influence of the geographic regions of the state of São Paulo (Brazil) (West, Center, East, South, and North) on the sampled bacterial communities in terms of abundance (permanova weighted unifrac distance: p-value = 0.026, Pseudo F = 1.534) and presence/absence (permanova unweighted unifrac distance: p-value = 0.001, Pseudo F = 1.438). That is, there is a difference not only in the composition but also in the abundance of bacterial communities associated with leaf-cutting ants when considering the geographic regions of the state of São Paulo (Brazil) (Figure 6A and 6B). The network analysis shows the complexity of the bacterial community associated with A. sexdens (Figure 6C).

Figure 6. Plot illustrating bacterial community structure among Atta sexdens workers of different regions (West, Center, East, South, and North), and habitats (natural and agricultural). (A). Principal coordinates analysis (PCoA) with weighted unifrac distance matrix. (B). PCoA with unweighted unifrac distance matrix. (C). Network analyses of bacteria present in Atta sexdens workers in different habitats. Note the complexity of the bacterial communities associated with leaf-cutting ants. Colors represent different habitats, and shapes represent different geographical regions.

Insects 2020, 11, 332 11 of 21

We also tested for any significant differences between the habitats (natural: Cerrado and Atlantic Forest; and agricultural: sugar cane, citrus and eucalyptus) present within each geographic region. For the Center region we found significant results in both the abundance analysis and composition analysis (which tests for presence or absence) (Table 2). For the other regions analyzed separately, our results show that there is a significant difference in bacterial composition (presence and absence) associated with different habitats, but not in abundance. When we conducted analyses between natural and agricultural habitats, although we were able to reveal interesting results regarding the Tsukamurella influence mainly associated with citrus, some data were scrambled and diluted in the vast amount of information generated by this technique. When first investigating the influence of geographic location (West, Center, East, South and North), we determined that it does have a large impact on the weighted and unweighted composition of the bacterial community associated with A. sexdens. However, when the analyses were conducted by region (separately), the results showed significant differences in the composition (unweighted) of bacterial communities associated with different habitats, i.e., natural (Cerrado and Atlantic Forest) and agricultural (sugar cane, citrus and eucalyptus).

Table 2. Abundance (permanova weighted unifrac distance) and composition (permanova unweighted unifrac distance) analysis of the microbiota associated with A. sexdens collected at different geographical locations. This table shows significant differences mainly in the composition between the habitat (natural: Cerrado and Atlantic Forest; and agricultural: sugarcane, citrus and eucalyptus) in each geographical location.

Permanova Weighted Unifrac Permanova Unweighted Unifrac Distance Regions Distance p-Value Pseudo F p-Value Pseudo F Center 0.025 1.940 0.001 1.162 East 0.254 1.127 0.018 1.212 North 0.494 0.977 0.001 1.192 West 0.075 1.43 0.014 1.280 South 0.376 1.092 0.016 1.198

To test which pairs of geographic regions contributed to the difference found in Figure 6, paired permanova tests verified abundance (weighted unifrac distance) and composition (unweighted unifrac distance) and found that several groupings showed significant differences (see Figure 7). In general, all regions (West, Center, East, South, and North) influenced the abundance and composition of bacterial communities associated with this leaf-cutting ant. Therefore, the pairs that obtained significant results were subjected to SIMPER analysis, where we obtained the main bacteria associated with these groupings (peer regions) (see Table 3). Our results are consistent, showing that Tsukamurella is one of the major bacteria responsible for differences in all groupings of geographic regions. These results were also confirmed by the analysis in Table 1, further supporting that this bacterium has a significant impact on the bacterial community of A. sexdens. In addition, Mesoplasma, the main bacterium found in the present study (see Figure 2), was also frequent and one of the major bacterial genera responsible for the differences found in almost all clusters, except Center vs. North.

Insects 2020, 11, 332 12 of 21

Figure 7. Pairwise permanova results of different regions. A. Weighted unifrac distance matrix B. Weighted unifrac distance matrix. Results highlighted in red indicate p values greater than 0.05. Bars represent the pairs of regions with significant differences.

Insects 2020, 11, 332 13 of 21

Table 3. SIMPER analysis of main bacteria contributing to differences between each region pair. This test indicates the contribution of specific ASVs to the observed differences in community structure among different groups. Note that Tsukamurella remains the most influential bacteria, present in all clusters.

Overall Average Most Influential Percent Contribution to Regions Cumulative % Dissimilarity ASV/Taxonomy Difference Alphaproteobacteria 10.63 10.63 Aeromicrobium 8.38 19.01 Gammaproteobacteria 7.64 26.65 Center vs. East 74.77% Mesoplasma 6.925 33.58 Pseudomonas 6.403 39.98 Tsukamurella 5.928 45.91 Micrococcales 5.739 51.65 Alphaproteobacteria 8.30 8.30 Gammaproteobacteria 8.14 16.44 Aeromicrobium 7.96 24.40 East vs. North 75.89% Tsukamurella 7.43 31.83 Mesoplasma 7.30 39.13 Nocardioides 6.23 45.36 Pseudomonas 4.92 50.28 Aeromicrobium 12.93 12.93 Tsukamurella 11.34 24.27 Micrococcales 6.81 31.08 Center vs. 73.01% Gammaproteobacteria 6.52 37.60 North Nocardioides 5.188 42.79 Alphaproteobacteria 4.634 47.42 Pseudomonas 4.315 51.74 Tsukamurella 11.01 11.01 Aeromicrobium 9.871 20.88 Enterobacteriaceae 7.052 27.93 Center vs. West 82.17% Mesoplasma 6.649 34.58 Micrococcales 6.587 41.17 Alphaproteobacteria 6.284 47.45 Bacteria 5.37 52.82 Aeromicrobium 8.51 8.51 Mesoplasma 8.15 16.66 Anthococcus 8.05 24.70 Center vs. Tsukamurella 7.87 32.57 79.77% South Micrococcales 5.74 38.31 Burkholderiaceae 5.10 43.41 Pseudomonas 4.73 48.14 Alphaproteobacteria 4.57 52.71 Mesoplasma 12.77 12.77 Alphaproteobacteria 10.09 22.86 Enterobacteriaceae 7.458 30.32 East vs. West 82.48% Tsukamurella 6.881 37.20 Gammaproteobacteria 6.798 44.00 Pseudomonas 5.869 49.87 Bacteria 4.881 54.75 Tsukamurella 11.83 11.83 Aeromicrobium 9.297 21.13 Mesoplasma 6.928 28.06 Nocardioides 5.837 33.89 North vs. West 81.97% Gammaproteobacteria 5.473 39.37 Enterobacteriaceae 5.438 44.8 Bacteria 4.399 49.2 Nocardiopsis 4.18 53.39 Anthococcus 9.13 9.13 Tsukamurella 8.84 17.98 Mesoplasma 8.41 26.39 North vs. South 80.48% Aeromicrobium 8.16 34.55 Nocardioides 5.86 40.41 Gammaproteobacteria 5.39 45.79 Burkholderiaceae 4.52 50.31 West vs. South 86.00% Mesoplasma 13.33 13.33

Insects 2020, 11, 332 14 of 21

Tsukamurella 8.571 21.90 Anthococcus 8.282 30.18 Enterobacteriaceae 5.651 35.83 Aeromicrobium 5.253 41.08 Burkholderiaceae 4.228 45.31 Bacteria 4.195 49.51 Pseudomonas 4.104 53.61

3.3. Alpha Diversity: Different Habitats and Geographical Locations Alpha diversity also presented distinct results when the analysis was conducted comparing only the different habitats and geographic regions with the Faith phylogeny index. According to a Kruskall–Wallis test, only geographic region influences the alpha diversity associated with A. sexdens (H = 16,795, p-value = 0.002). As this richness index is qualitative, it takes into account the composition of the bacterial community (presence and absence). However, the same analysis conducted with the Shannon index regardless of abundance did not yield any significant results (see Figure S4). Therefore, the different geographical regions may influence the composition of alpha diversity associated with A. sexdens, but do not appear to influence the abundance.

4. Discussion There are several factors that can predict and influence the diversity (composition and abundance) of microbial communities associated with a given host [8]. The microbiota of leaf-cutting ants has already been addressed in several studies [23,27,34–36,69–71], but there remains much to understand about the drivers of these symbioses. This lack of knowledge motivates the present study, which is the first to compare the bacterial communities of ants in natural (Cerrado and Atlantic Forest) to agricultural (sugar cane, citrus and eucalyptus) habitats, as well as consider geographical location (West, Center, East, South, and North of the state of São Paulo-Brazil), to assess the factors that may influence microbial diversity. We found differences in microbiota of the different samples from the same species in different habitats and geographic locations. While each strain of bacteria can be acquired in a different way [72], our results may reflect the fact that these bacteria are being acquired from different environmental conditions (e.g., landscape simplification, different plant availabilities, agrochemicals, and soil), through horizontal more than vertical transmission. That is, even with sophisticated defense systems and behavior utilized by A. sexdens [49,73,74], the bacterial communities of the present study were altered mainly due to differing geographical locations, and secondarily due to natural and agricultural habitats. There are several studies that show that the composition and abundance of the microbiota is impacted by horizontal acquisition [3,47,70,75–77]. Our study corroborates those that have found that the geographic location of the host influences the diversity of the bacterial community, whether in the case of mammals [75], birds [78,79], other insects [47,77,80], or also of ants [3,70,76,81]. For example, the bacterial diversity recovered from the mesosoma and head of Camponotus workers is influenced by the habitat where these ants were collected, which suggests that these bacteria are being acquired via horizontal transfer [3]. Another study that reinforces the impact of geographical location on the diversity of bacterial communities was conducted by Martins and Moreau [76] with Pheidole species collected around the world. The same pattern of acquisition of microbial communities was also reported in fourth-instar larvae of Solenopsis invicta ants [81] and the lower attine species Mycocepurus smithii [70]. workers also showed a diverse and unstable microbiota with Mollicutes being the main bacteria recovered [36]. The bacterial community recovered from A. sexdens is also very complex, and this study was able to identify several important bacterial taxa associated with this species. ASVs identified as Mesoplasma sp. (Mollicutes, Entomoplasmatales) is one of the most prevalent bacteria in our sample with 7% relative frequency. Our results corroborate a previous study [23] that reports evidence of this bacterium mainly associated with the post-pharyngeal gland (70%) and also in smaller quantities in the midgut of the queen and the whole body of the worker of A. sexdens. Meirelles et al. [35] describes the association of Mesoplasma sp. with infrabucal pellets carried

Insects 2020, 11, 332 15 of 21 by A. texana queens. Also in this study, through phylogenetic analysis, the authors show specificity of this bacterium associated with eight different genera of Attine, since it formed a monophyletic clade. Mesoplasma sp. associated to the army ants [82] was also grouped in another clade, also showing specificity to army ants [83]. The function of this ant-associated bacterium was investigated with the genome of two strains, and it was found that each strain was acquired independently from a different origin. These bacteria are capable of decomposing arginine and providing nitrogen-rich amino-acid [84]. In that same study, the authors also concluded that the relative abundance of Mesoplasma sp. is related to the type of substrate that the ant worker offers to the fungus garden. Our study found a high relative abundance of these bacteria, which suggests that it must have been facilitated by the substrate of A. sexdens. The next most abundant ASV recovered in the present study is Tsukamurella sp. (Actinobacteria), which despite occurring in other habitats was strongly related to the citrus agricultural habitat. Although its function is not fully understood, there are studies reporting that these bacteria may be an endophytic growth promoter [85] minimizing the growth of phytopathogenic fungi [86], and also may act as a nematicide [87] used in agriculture. The third most common ASV was identified as Aeromicrobium sp. (Actinobacteria), which has also been reported on species of ants as well [88] but is generally associated with the soil [89]. Other studies have reported the presence of ant-associated Actinobacteria, especially the genus Pseudonocardia sp., often associated with other Attines ants [71,90]. These bacteria are capable of producing compounds with properties that defend the fungal garden and the ant colony from disease [25,71,90]. The samples of A. sexdens from our study were also associated with these bacteria in smaller quantities compared to the literature, which suggests that perhaps A. sexdens is also relying on other measures to maintain the health of the fungus garden and the colony. Further studies are needed to address this issue. The phylum Proteobacteria has also been reported in A. cephalotes Linnaeus 1758 and A. colombica Guérin-Méneville 1844 [69,91]. Pseudomonas and Burkonderialles are bacteria belonging to this phylum, and were associated with the workers of leaf-cutting ants of the present study, although in smaller abundances. Some bacteria present in the insect microbiome may be able to metabolize agrochemicals, making the host pesticide-resistant [92–94], and it has been shown that certain strains of the genus Pseudomonas sp. are associated with insect resistance to insecticides [95]. The bacterial community is subject to several environmental stressors that may alter its composition [96]. As an example, certain strains of Burkholderia sp. are capable of degrading insecticide fenitrothion (organophosphate) [97]. Tago et al.[98] found an increase of this bacterium in the soil after the application of the pesticide in sugar cane agricultural habitats. In another study of Acromyrmex spp., besides Entomoplasmatales, the bacteria Wolbachia and Rhizobiales were highlighted in association with these ants. These latter two bacteria are highly associated with several genera of ants [2,8,36,99]. While the function of Wolbachia associated with ants is still unknown, it is transmitted vertically, although its diversity is influenced by geography [17,36,100]. On the other hand, Rhizobiales is also often associated with nitrogen fixation and recycling [2,69], which is extremely important especially in the case of herbivorous ants. Our results also recovered both bacteria associated with A. sexdens. The microbiota associated with a host can be composed of mutualists, commensals and pathogens, and many of the bacteria associated with A. sexdens are probably transient bacteria that have no associated role in ants. With the 16S rRNA amplicon technique, we cannot be sure of the microbe–host relationships or their functions. However, overall these bacteria may have been acquired horizontally and thus may reflect the environment (habitats). In addition, although the bacterial community of A. sexdens is unstable and highly variable, we do not know how this complexity varies in the different castes within the colony, and its implications for the host. In the case of Atta cephalotes it was found that the bacteria Mollicutes passes via social transmission of workers to the larvae [36]. Future studies that consider the different castes and location of this microbiota in the ant body may help to elucidate the mode of transmission and consequently the function of these bacteria associated with A. sexdens.

Insects 2020, 11, 332 16 of 21

5. Conclusions Genomic approaches have become increasingly accessible and can be used as excellent tools to explain several biological phenomena such as intra- and interspecific relationships, species ecology such as diet and symbioses, as well as molecular, physiological, behavioral and evolutionary features of the host. These approaches have the potential to generate data that can contribute to the understanding of symbiotic interactions within the field of myrmecology [2,101–105]. Recent next- generation sequencing techniques have made it possible to analyze the communities of bacteria present in several hosts. Specifically for ants, many studies have demonstrated the importance of these symbionts and how they can influence host attributes [2,6,8,10,20,106]. However, studies of microbial interactions are challenging because they involve dynamic systems that can be shaped by the sum of abiotic factors, as presented here. This is the first study that compared the bacterial communities of A. sexdens in natural habitats, such as Cerrado and Atlantic Forest, to those in agricultural habitats (sugar cane, citrus and eucalyptus) and that considered geographic location (West, Center, East, South, and North of the state of São Paulo-Brazil), showing that all these factors may influence microbial diversity. Knowledge about this insect-associated microbiota may help to develop strategies to increase insect pest resistance in agricultural habitats, as well identify bacteria that can act as bioremediation of pesticides and help recover degraded areas [107].

Supplementary Materials: The following are available online at www.mdpi.com/2075-4450/11/6/332/s1, Figure S1: Location of Atta sexdens collection areas in São Paulo State, Brazil, Figure S2: Rarefaction curves were used to estimate richness in the Shannon index for 16S rRNA. The vertical axis shows the Shannon index observed and the number of sequences per sample is shown on the horizontal axis, Figure S3: The overview of bacteria relative abundance recovered from this study of Atta sexdens, Figure S4: Alpha diversity associated with Atta sexdens, Table S1: Samples of Atta sexdens collected from different natural and agriculture habitats, Table S2: Amplicon Sequence variants (ASVs) table highlighting each bacterium recovered from each sample of Atta sexdens, Table S3: SIMPER analysis of main bacteria in all differents habitats. This test indicates the contribution of specific ASVs to the observed differences in community structure among different groups. Note that Tsukamurella remains the most influential bacteria associated with citrus crop,

Author Contributions: M.d.O.R., O.C.B., and M.S.C.M. designed research, M.d.O.R. performed research, M.d.O.R. and C.M. analyzed data, M.d.O.R., C.M., M.S.C.M. and O.C.B. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding: This research was funded by National Science Foundation (NSF DEB 1900357), and supported by Fundunesp and CNPq.

Acknowledgments: We thank Amanda Carlos, Alexandro Santana Vieira, Sergio Pascon, and Nathalia Sampaio da Silva for assistance in the field. We also thank all members of the Formigas Cortadeiras Lab from the Social Insects Study Center (CEIS, Unesp RC) and the Moreau Lab for their support, and especially Corrie S. Moreau, Peter J. Flynn, and Benjamin D. Blanchard for valuable comments and English corrections. M.d.O.R. thanks the U.S. National Science Foundation (NSF DEB 1900357) for financial support. This research was supported by Fundunesp and CNPq.

Conflicts of Interest: All authors have declared no conflict of interest.

References

1. Ley, R.E.; Hamady, M.; Lozupone, C.; Turnbaugh, P.J.; Ramey, R.R.; Bircher, J.S.; Schlegel, M.L.; Tucker, T.A.; Schrenzel, M.D.; Knight, R.; et al. of Mammals and Their Gut Microbes. Science 2008, 320, 1647–1651. 2. Hu, Y.; Sanders, J.G.; Łukasik, P.; D’Amelio, C.L.; Millar, J.S.; Vann, D.R.; Lan, Y.; Newton, J.A.; Schotanus, M.; Kronauer, D.J.C.; et al. Herbivorous turtle ants obtain essential nutrients from a conserved nitrogen- recycling gut microbiome. Nat. Commun. 2018, 9, 964. 3. Ramalho, M.O.; Moreau, C.S.; Bueno, O.C. The Potential Role of Environment in Structuring the Microbiota of Camponotus across Parts of the Body. Adv. Entomol. 2019, 7, 47–70. 4. Ramalho, M.O.; Bueno, O.C.; Moreau, C.S. Species-specific signatures of the microbiome from Camponotus and Colobopsis ants across developmental stages. PLoS ONE 2017, 12, e0187461.

Insects 2020, 11, 332 17 of 21

5. Ramalho, M.O.; Bueno, O.C.; Moreau, C.S. Microbial composition of spiny ants (: Formicidae: Polyrhachis) across their geographic range. BMC Evol. Biol. 2017, 17, 96. 6. Rubin, B.E.; Kautz, S.; Wray, B.D.; Moreau, C.S. Dietary specialization in mutualistic acacia-ants affects relative abundance but not identity of host-associated bacteria. Mol. Ecol. 2018, 28, 900–916. 7. Sanders, J.G.; Łukasik, P.; Frederickson, M.E.; Russell, J.A.; Koga, R.; Knight, R.; Pierce, N.E. Dramatic Differences in Gut Bacterial Densities Correlate with Diet and Habitat in Ants. Integr. Comp. Biol. 2017, 57, 705–722. 8. Ramalho, M.O.; Duplais, C.; Orivel, J.; Gibson, J.C.; Dejan, A.; Suarez, A.V.; Moreau, C.S. Development but not diet alters microbial communities in the Neotropical arboreal trap jaw ant Daceton armigerum: An exploratory study. Sci. Rep. 2020, 10, 7350. 9. Russell, J.A.; Moreau, C.S.; Goldman-Huertas, B.; Fujiwara, M.; Lohman, D.J.; Pierce, N.E. Bacterial gut symbionts are tightly linked with the evolution of herbivory in ants. Proc. Natl. Acad. Sci. USA 2009, 106, 21236–21241. 10. Sanders, J.G.; Powell, S.; Kronauer, D.J.C.; Vasconcelos, H.L.; Frederickson, M.E.; Pierce, N.E. Stability and phylogenetic correlation in gut microbiota: Lessons from ants and apes. Mol. Ecol. 2014, 23, 1268–1283. 11. Roeselers, G.; Mittge, E.K.; Stephens, W.Z.; Parichy, D.M.; Cavanaugh, C.M.; Guillemin, K.; Rawls, J.F. Evidence for a core gut microbiota in the zebrafish. ISME J. 2011, 5, 1595–1608. 12. Turnbaugh, P.J.; Ridaura, V.K.; Faith, J.J.; Rey, F.E.; Knight, R.; Gordon, J.I. The effect of diet on the human gut microbiome: A metagenomic analysis in humanized gnotobiotic mice. Sci. Transl. Med. 2009, 1, 6ra14. 13. Dillon, R.J.; Webster, G.; Weightman, A.J.; Keith Charnley, A. Diversity of gut microbiota increases with aging and starvation in the desert locust. Antonie Van Leeuwenhoek 2010, 97, 69–77. 14. Hu, Y.; Łukasik, P.; Moreau, C.S.; Russell, J.A. Correlates of gut community composition across an ant species (Cephalotes varians) elucidate causes and consequences of symbiotic variability. Mol. Ecol. 2014, 23, 1284–1300. 15. Lanan, M.C.; Augusto, P.; Rodrigues, P.; Agellon, A.; Jansma, P.; Wheeler, D.E. A bacterial filter protects and structures the gut microbiome of an insect. ISME J. 2015, 12. 16. Degnan, P.H.; Lazarus, A.B.; Wernegreen, J.J. Genome sequence of Blochmannia pennsylvanicus indicates parallel evolutionary trends among bacterial mutualists of insects. Genome Res. 2005, 15, 1023–33. 17. Ramalho, M.O.; Vieira, A.S.; Pereira, M.C.; Moreau, C.S.; Bueno, O.C. Transovarian Transmission of Blochmannia and Wolbachia Endosymbionts in the Neotropical Weaver Ant Camponotus textor (Hymenoptera, Formicidae). Curr. Microbiol. 2018, 75, 866–873. 18. Wernegreen, J.J.; Kauppinen, S.N.; Brady, S.G.; Ward, P.S. One nutritional symbiosis begat another: Phylogenetic evidence that the ant tribe Camponotini acquired Blochmannia by tending sap-feeding insects. BMC Evol. Biol. 2009, 9, 292. 19. AntWeb AntWeb. Available online: https://www.antweb.org (accessed on 24 February 2017). 20. Russell, J.A.; Sanders, J.G.; Moreau, C.S. Hotspots for symbiosis: Function, evolution, and specificity of ant- microbe associa-tions from trunk to tips of the ant phylogeny (Hymenoptera: Formicidae). Myrmecol. News 2017, 24, 43–69. 21. Bolton, B. An Online Catalog of the Ants of the World. Available online: http://www.antcat.org/(accessed on 20 October 2016). 22. Andersen, S.B.; Boye, M.; Nash, D.R.; Boomsma, J.J. Dynamic Wolbachia prevalence in Acromyrmex leaf- cutting ants: Potential for a nutritional symbiosis. J. Evol. Biol. 2012, 25, 1340–1350. 23. Vieira, A.S.; Ramalho, M.O.; Martins, C.; Martins, V.G.; Bueno, O.C. Microbial Communities in Different Tissues of Atta sexdens rubropilosa Leaf-cutting Ants. Curr. Microbiol. 2017, 74, 1216–1225. 24. Santos, A.V.; Dillon, R.J.; Dillon, V.M.; Reynolds, S.E.; Samuels, R.I. Ocurrence of the antibiotic producing bacterium Burkholderia sp. in colonies of the leaf-cutting ant Atta sexdens rubropilosa. FEMS Microbiol. Lett. 2004, 239, 319–323. 25. Poulsen, M.; Cafaro, M.J.; Erhardt, D.P.; Little, A.E.F.; Gerardo, N.M.; Tebbets, B.; Klein, B.S.; Currie, C.R. Variation in Pseudonocardia antibiotic defence helps govern parasite-induced morbidity in Acromyrmex leaf-cutting ants. Environ. Microbiol. Rep. 2009, 2, 534–540. 26. Della Lucia, T.M.; Gandra, L.C.; Guedes, R.N. Managing leaf-cutting ants: Peculiarities, trends and challenges. Pest Manag. Sci. 2014, 70, 14–23. 27. Mueller, U.G.; Dash, D.; Rabeling, C.; Rodrigues, A. Coevolution between Attine ants and Actinomycete bacteria: A reevaluation. Evolution (N. Y) 2008, 62, 2894–2912.

Insects 2020, 11, 332 18 of 21

28. Della Lucia, T. Formigas Cortadeiras: Da Bioecologia ao Manejo; UFV: Viçosa, Brazil, 2011. 29. Montoya-Lerma, J.; Giraldo-Echeverri, C.; Armbrecht, I.; Farji-Brener, A.; Calle, Z. Leaf-cutting ants revisited: Towards rational management and control. Int. J. Pest Manag. 2012, 58, 225–247. 30. Bueno, O.C.; Bueno, F.C. Manejo de Pragas Urbanas. In Controle de Formigas em Áreas Urbanas.; São FEALQ: Paulo, Brazil, 2007; pp. 67–77. 31. Ortiz, G.; Cintra-Socolowski, P.; Vieira, A.S.; Bueno, O.C. Ação de produtos utilizados no cultivo da cana- de-açúcar sobre as formigas-cortadeiras. In Cana-de-Açúcar e seus Impactos: Uma visão Acadêmica; Fontanetti, C.S., Bueno, O.C., Eds.; Canal6 editora: Bauru, Brazil, 2017; pp. 179–196. 32. Delabie, J.H.C.; Alves, H.R.; Reuss-Strenzel, G.M.; Carmo, A.D.; Nascimento, I.D. Distribuição das formigas cortadeiras dos gêneros Acromyrmex e Atta no Novo Mundo. In Formigas Cortadeiras: Da Bioecologia ao Manejo; Universidade Federal de Viçosa: Viçosa, Brazil, 2011; pp. 80–101. 33. Groussin, M.; Mazel, F.; Sanders, J.G.; Smillie, C.S.; Lavergne, S.; Thuiller, W.; Alm, E.J. Unraveling the processes shaping mammalian gut microbiomes over evolutionary time. Nat. Commun. 2017, 8, 14319. 34. Teseo, S.; van Zweden, J.S.; Pontieri, L.; Kooij, P.W.; Sørensen, S.J.; Wenseleers, T.; Poulsen, M.; Boomsma, J.J.; Sapountzis, P. The scent of symbiosis: Gut bacteria may affect social interactions in leaf-cutting ants. Anim. Behav. 2019, 150, 239–254. 35. Meirelles, L.A.; McFrederick, Q.S.; Rodrigues, A.; Mantovani, J.D.; de Melo Rodovalho, C.; Ferreira, H.; Bacci, M.; Mueller, U.G. Bacterial microbiomes from vertically transmitted fungal inocula of the leaf-cutting ant Atta texana. Environ. Microbiol. Rep. 2016, 8, 630–640. 36. Zhukova, M.; Sapountzis, P.; Schiøtt, M.; Boomsma, J.J. Diversity and transmission of gut bacteria in Atta and Acromyrmex leaf-cutting ants during development. Front. Microbiol. 2017, 8, 1–14. 37. Ratter, J. The Brazilian Cerrado Vegetation and Threats to its Biodiversity. Ann. Bot. 1997, 80, 223–230. 38. Sano, E.E.; Rosa, R.; Luís, J.; Laerte, S.B.; Ferreira, G. Documentos 190 Mapeamento de Cobertura Vegetal do Bioma Cerrado: Estratégias e resultados Embrapa Cerrados Planaltina, DF 2007 Empresa Brasileira de Pesquisa Agropecuária Embrapa Cerrados Ministério da Agricultura, Pecuária e Abastecimento. 2007. Available: https://www.infoteca.cnptia.embrapa.br/bitstream/doc/570887/1/doc190.pdf. Accessed at 03 November 2019. 39. Joly, C.A.; Metzger, J.P.; Tabarelli, M. Experiences from the Brazilian Atlantic Forest: Ecological findings and conservation initiatives. New Phytol. 2014, 204, 459–473. 40. Ribeiro, M.C.; Metzger, J.P.; Martensen, A.C.; Ponzoni, F.J.; Hirota, M.M. The Brazilian Atlantic Forest: How much is left, and how is the remaining forest distributed? Implications for conservation. Biol. Conserv. 2009, 142, 1141–1153. 41. SOS Mata Atlântica. Atlas dos Remanescentes Florestais de Mata Atlântica. Relatório Técnico; 2019. Available: https://www.sosma.org.br/wp-content/uploads/2019/06/Atlas-mata-atlanticaDIGITAL.pdf. Accessed on 3 November 2019. 42. Cepea PIB do Agronegócio de São Paulo-Centro de Estudos Avançados em Economia Aplicada-CEPEA- Esalq/USP. Available online: https://www.cepea.esalq.usp.br/br/pib-do-agronegocio-de-sao-paulo.aspx (accessed on 5 September 2019). 43. Power, E.F.; Kelly, D.L.; Stout, J.C. Organic Farming and Landscape Structure: Effects on Insect-Pollinated Plant Diversity in Intensively Managed Grasslands. PLoS ONE 2012, 7, e38073. 44. Kaluza, B.F.; Wallace, H.M.; Heard, T.A.; Minden, V.; Klein, A.; Leonhardt, S.D. Social bees are fitter in more biodiverse environments. Sci. Rep. 2018, 8, 1–10. 45. Dharampal, P.S.; Diaz-Garcia, L.; Haase, M.A.B.; Zalapa, J.; Currie, C.R.; Hittinger, C.T.; Steffan, S.A. Microbial Diversity Associated with the Pollen Stores of Captive-Bred Bumble Bee Colonies. Insects 2020, 11, 250. 46. Dharampal, P.S.; Carlson, C.; Currie, C.R.; Steffan, S.A. Pollen-borne microbes shape bee fitness. Proc. R. Soc. B Biol. Sci. 2019, 286, 20182894. 47. Hannula, S.E.; Zhu, F.; Heinen, R.; Bezemer, T.M. Foliar-feeding insects acquire microbiomes from the soil rather than the host plant. Nat. Commun. 2019, 10, 1–9. 48. Moran, N.A.; Ochman, H.; Hammer, T.J. Evolutionary and Ecological Consequences of Gut Microbial Communities. Annu. Rev. Ecol. Evol. Syst. 2019, 50, 451–475. 49. Mackay, W.; Mackay, E. Las hormigas de Colombia: Arrieras del Genero/Atta/(Hymenoptera-Formicidae). Revista Colombiana de Entomología. 1986, 12, 23-30.

Insects 2020, 11, 332 19 of 21

50. Moreau, C.S. A practical guide to DNA extraction, PCR, and gene-based DNA sequencing in insects. Halteres 2014, 5, 32–42. 51. Caporaso, J.G.; Lauber, C.L.; Walters, W.A.; Berg-Lyons, D.; Huntley, J.; Fierer, N.; Owens, S.M.; Betley, J.; Fraser, L.; Bauer, M.; et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012, 6, 1621–1624. 52. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. 53. McDonald, D.; Clemente, J.C.; Kuczynski, J.; Rideout, J.R.; Stombaugh, J.; Wendel, D.; Wilke, A.; Huse, S.; Hufnagle, J.; Meyer, F.; et al. The Biological Observation Matrix (BIOM) format or: How I learned to stop worrying and love the ome-ome. Gigascience 2012, 1, 7. 54. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High- resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. 55. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glockner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013, 41, D590–D596. 56. Yilmaz, P.; Parfrey, L.W.; Yarza, P.; Gerken, J.; Pruesse, E.; Quast, C.; Schweer, T.; Peplies, J.; Ludwig, W.; Glöckner, F.O. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucleic Acids Res. 2014, 42, D643–D648. 57. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Gregory Caporaso, J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2- feature-classifier plugin. Microbiome 2018, 6, 90. 58. Davis, N.M.; Proctor, D.; Holmes, S.P.; Relman, D.A.; Callahan, B.J. Simple statistical identification and removal of contaminant sequences in marker-gene and metagenomics data. Microbiome 2018, 6, 226. 59. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 2013, 30, 772–780. 60. Price, M.N.; Dehal, P.S.; Arkin, A.P. FastTree 2—Approximately Maximum-Likelihood Trees for Large Alignments. PLoS ONE 2010, 5, e9490. 61. Vázquez-Baeza, Y.; Pirrung, M.; Gonzalez, A.; Knight, R. EMPeror: A tool for visualizing high-throughput microbial community data. Gigascience 2013, 2, 16. 62. Anderson, M.J. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001, 26, 32–46. 63. Clarke, K.R. Non-parametric multivariate analyses of changes in community structure. Austral Ecol. 1993, 18, 117–143. 64. Hammer, Ø. Harper, D.A.T.; Ryan, P.D. Paleontological Statistics Software: Package for Education and Data Analysis. Palaeontol. Electron. 2001, 1–9. 65. McCune, B.; Grace, J.B.; Urban, D. Analysis of Ecological Communities: Gleneden Beach. Oregon, MJM Software Design.; 2002. 66. R Development Core Team. R: A Language and Environment for Statistical Computing; R Development Core Team: Vienna, Austria, 2019. 67. McMurdie, P.J.; Holmes, S. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. 68. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer: Cham, Switzerland, 2009; ISBN 9780387981406. 69. Sapountzis, P.; Zhukova, M.; Hansen, L.H.; Sørensen, S.J.; Schiøtt, M.; Boomsma, J.J. Acromyrmex leaf- cutting ants have simple gut microbiota with nitrogen-fixing potential. Appl. Environ. Microbiol. 2015, 81, 5527–5537. 70. Kellner, K.; Ishak, H.D.; Linksvayer, T.A.; Mueller, U.G. Bacterial community composition and diversity in an ancestral ant fungus symbiosis. FEMS Microbiol. Ecol. 2015, 91. 71. Woodhams, D.C.; Brucker, R.M. Disease defence through generations: Leaf-cutter ants and their symbiotic bacteria. Mol. Ecol. 2013, 22, 4141–4143. 72. Tseng, S.P.; Wetterer, J.K.; Suarez, A.V.; Lee, C.Y.; Yoshimura, T.; Shoemaker, D.W.; Yang, C.C.S. Genetic Diversity and Wolbachia Infection Patterns in a Globally Distributed Invasive Ant. Front. Genet. 2019, 10, 838.

Insects 2020, 11, 332 20 of 21

73. Schultz, T.R.; Brady, S.G. Major Evolutionary Transitions in Ant Agriculture; National Academy of Sciences: Washington, DC, USA, 2008; Volume 105, pp. 5435–5440. 74. Chevrette, M.G.; Bratburd, J.R.; Currie, C.R.; Stubbendieck, R.M. Experimental Microbiomes: Models Not to Scale. mSystems 2019, 4, e00175-19. 75. Linnenbrink, M.; Wang, J.; Hardouin, E.A.; Künzel, S.; Metzler, D.; Baines, J.F. The role of biogeography in shaping diversity of the intestinal microbiota in house mice. Mol. Ecol. 2013, 22, 1904–1916. 76. Martins, C.; Moreau, C.S. Influence of host phylogeny, geographical location and seed harvesting diet on the bacterial community of globally distributed Pheidole ants. PeerJ 2020, 2020, e8492. 77. Kikuchi, Y.; Hosokawa, T.; Fukatsu, T. Insect-microbe mutualism without vertical transmission: A stinkbug acquires a beneficial gut symbiont from the environment every generation. Appl. Environ. Microbiol. 2007, 73, 4308–4316. 78. Ruiz-Rodríguez, M.; Soler, J.J.; Martín-Vivaldi, M.; Martín-Platero, A.M.; Méndez, M.; Peralta-Sánchez, J.M.; Ananou, S.; Valdivia, E.; Martínez-Buenob, M. Environmental factors shape the community of symbionts in the hoopoe uropygial gland more than genetic factors. Appl. Environ. Microbiol. 2014, 80, 6714– 6723. 79. Song, S.J.; Sanders, J.G.; Delsuc, F.; Metcalf, J.; Amato, K.; Taylor, M.W.; Mazel, F.; Lutz, H.L.; Winker, K.; Graves, G.R.; et al. Comparative analyses of vertebrate gut microbiomes reveal convergence between birds and bats. MBio 2020, 10.1128/mBio.02901-19. 80. Zhao, Y.; Zhang, S.; Luo, J.Y.; Wang, C.Y.; Lv, L.M.; Cui, J.J. Bacterial communities of the cotton aphid Aphis gossypii associated with Bt cotton in northern China. Sci. Rep. 2016, 6, 1–8. 81. Lee, A.H.; Husseneder, C.; Hooper-Bùi, L. Culture-independent identification of gut bacteria in fourth- instar red imported fire ant, Solenopsis invicta Buren, larvae. J. Invertebr. Pathol. 2008, 98, 20–33. 82. Funaro, C.F.; Kronauer, D.J.C.; Moreau, C.S.; Goldman-Huertas, B.; Pierce, N.E.; Russell, J.A. Army ants harbor a host-specific clade of Entomoplasmatales bacteria. Appl. Environ. Microbiol. 2011, 77, 346–50. 83. Łukasik, P.; Newton, J.A.; Sanders, J.G.; Hu, Y.; Moreau, C.S.; Kronauer, D.J.C.; O’Donnell, S.; Koga, R.; Russell, J.A. The structured diversity of specialized gut symbionts of the New World army ants. Mol. Ecol. 2017, 26, 3808–3825. 84. Sapountzis, P.; Zhukova, M.; Shik, J.Z.; Schiott, M.; Boomsma, J.J. Reconstructing the functions of endosymbiotic mollicutes in fungus-growing ants. Elife 2018, 7, e39209 85. Araújo Teixeira, M.; Soares de Melo, I.; Faria Vieira, R.; Eduardo Carvalho Costa, F.; Harakava, R. Microrganismos Endofíticos de Mandioca de Áreas Comerciais e Etnovariedades em três Estados Brasileiros. Pesquisa Agripecuária Brasileira. 2007, 42, 42-49. 86. De Assumpção, L.C.; Lacava, P.T.; Dias, A.C.F.; de Azevedo, J.L.; Menten, J.O.M. Diversidade e potencial biotecnológico da comunidade bacteriana endofítica de sementes de soja. Pesqui. Agropecuária Bras. 2009, 44, 503–510. 87. Velazco, L.A.; Fernández, D.B. Study of Oral Acute Toxicity/Pathogenicity of the Tsukamurella Paurometabola Active Agent of the Biological Nematicide C-924 (Hebernem®). Adv. Pharmacoepidemiol. Drug Saf. 2012. 10.4172/apds.1000106 88. Ishak, H.D.; Plowes, R.; Sen, R.; Kellner, K.; Meyer, E.; Estrada, D.A.; Dowd, S.E.; Mueller, U.G. Bacterial Diversity in Solenopsis invicta and Solenopsis geminata Ant Colonies Characterized by 16S amplicon 454 Pyrosequencing. Microb. Ecol. 2011, 61, 821–831. 89. Reyes, R.D.H.; Cafaro, M.J. Paratrechina longicornis ants in a tropical dry forest harbor specific Actinobacteria diversity. J. Basic Microbiol. 2015, 55, 11–21. 90. Andersen, S.B.; Hansen, L.H.; Sapountzis, P.; Sørensen, S.J.; Boomsma, J.J. Specificity and stability of the Acromyrmex-Pseudonocardia symbiosis. Mol. Ecol. 2013, 22, 4307–4321. 91. Aylward, F.O.; Burnum, K.E.; Scott, J.J.; Suen, G.; Tringe, S.G.; Adams, S.M.; Barry, K.W.; Nicora, C.D.; Piehowski, P.D.; Purvine, S.O.; et al. Metagenomic and metaproteomic insights into bacterial communities in leaf-cutter ant fungus gardens. ISME J. 2012, 6, 1688–1701. 92. Bartha, R.; Lanzilotta, R.P.; Pramer, D. Stability and effects of some pesticides in soil. Appl. Microbiol. 1967, 15, 67–75. 93. Tu, C.M. Influence of five pyrethroid insecticides on microbial populations and activities in soil. Microb. Ecol. 1980, 5, 321–327. 94. van den Bosch, T.J.M.; Welte, C.U. Detoxifying symbionts in agriculturally important pest insects. Microb. Biotechnol. 2017, 10, 531–540.

Insects 2020, 11, 332 21 of 21

95. Saikia, N.; Das, S.K.; Patel, B.K.C.; Niwas, R.; Singh, A.; Gopal, M. Biodegradation of beta-cyfluthrin by Pseudomonas stutzeri strain S1. Biodegradation 2005, 16, 581–589. 96. Almeida, L.G.D.; Moraes, L.A.; Trigo, J.R.; Omoto, C.; Cônsoli, F.L. The gut microbiota of insecticide- resistant insects houses insecticide-degrading bacteria: A potential source for biotechnological exploitation. PLoS ONE 2017, 12, e0174754. 97. Kikuchi, Y.; Hayatsu, M.; Hosokawa, T.; Nagayama, A.; Tago, K.; Fukatsu, T. Symbiont-mediated insecticide resistance. Proc. Natl. Acad. Sci. USA 2012, 109, 8618–22. 98. Tago, K.; Okubo, T.; Itoh, H.; Kikuchi, Y.; Hori, T.; Sato, Y.; Nagayama, A.; Hayashi, K.; Ikeda, S.; Hayatsu, M. Insecticide-Degrading Burkholderia Symbionts of the Stinkbug Naturally Occupy Various Environments of Sugarcane Fields in a Southeast Island of Japan. Microbes Environ. 2015, 30, 29–36. 99. Reeves, D.D.; Price, S.L.; Ramalho, M.O.; Moreau, C.S. The Diversity and Distribution of Wolbachia, Rhizobiales, and Ophiocordyceps Within the Widespread Neotropical Turtle Ant, Cephalotes atratus (Hymenoptera: Formicidae). Neotrop. Entomol. 2020, 49, 52–60. 100. Kelly, M.; Price, S.L.; de Oliveira Ramalho, M.; Moreau, C.S. Diversity of Wolbachia Associated with the Giant Turtle Ant, Cephalotes atratus. Curr. Microbiol. 2019, 76, 1330–1337. 101. Gadagkar, R. The birth of ant genomics. Proc. Natl. Acad. Sci. USA 2011, 108, 5477–5478. 102. Gadau, J.; Helmkampf, M.; Nygaard, S.; Roux, J.; Simola, D.F.; Smith, C.R.; Suen, G.; Wurm, Y.; Smith, C.D. The genomic impact of 100 million years of social evolution in seven ant species. Trends Genet. 2012, 28, 14– 21. 103. Libbrecht, R.; Oxley, P.R.; Kronauer, D.J.; Keller, L. Ant genomics sheds light on the molecular regulation of social organization. Genome Biol. 2013, 14, 212. 104. Rubin, B.E.R.; Moreau, C.S. ARTICLE Comparative genomics reveals convergent rates of evolution in ant– plant mutualisms. Nat. Commun. 2016, 7. https://doi.org/10.1038/ncomms12679 105. Tsutsui, N.D. Dissecting ant recognition systems in the age of genomics. Biol. Lett. 2013, 9. 106. Russell, J.A.; Funaro, C.F.; Giraldo, Y.M.; Goldman-Huertas, B.; Suh, D.; Kronauer, D.J.C.; Moreau, C.S.; Pierce, N.E. A Veritable Menagerie of Heritable Bacteria from Ants, Butterflies, and Beyond: Broad Molecular Surveys and a Systematic Review. PLoS ONE 2012, 7, e51027. 107. Douglas, A.E. Strategies for Enhanced Crop Resistance to Insect Pests. Annu. Rev. Plant Biol. 2018, 69, 637– 660.

© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).