Next Article in Journal
Ecological and Hydrological Indicators of Climate Change Observed by Dryland Communities of Malipati in Chiredzi, Zimbabwe
Previous Article in Journal
It Is a Wild World in the City: Urban Wildlife Conservation and Communication in the Age of COVID-19
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thrips Microbiome Study in Commercial Avocado (Persea americana Mill.) from Northwest Colombian Andes (Antioquia, Colombia) Shows the Presence of Wolbachia, Ehrlichia, Enterobacter

by
Daniela Cano-Calle
1,2,3,*,
Luisa Maria Montoya-Porras
1,
Sebastian Ochoa-Giraldo
2,
Howard Junca
4,
Erika Garcia-Bonilla
4,
Clara Saldamando-Benjumea
2,3,
Claudia Ximena Moreno-Herrera
1,* and
Rafael E. Arango-Isaza
2,3,*
1
Microbiodiversity and Bioprospection Research Group, Laboratorio de Biología Celular 19A-310, Molecular, Facultad de Ciencias, Universidad Nacional de Colombia, Medellín 050034, Colombia
2
Plant Biotechnology UNALMED-CIB Group, Laboratorio de Ecología y Evolución de Insectos 16-223, Facultad de Ciencias, Universidad Nacional de Colombia, Medellín 050034, Colombia
3
Plant Biotechnology UNALMED-CIB Group, Corporación para Investigaciones Biológicas (CIB), Medellín 050034, Colombia
4
RG Microbial Ecology: Metabolism, Genomics & Evolution, Div. Ecogenomics & Holobionts, Microbiomas Foundation, LT11A, Chía 250008, Colombia
*
Authors to whom correspondence should be addressed.
Diversity 2022, 14(7), 540; https://doi.org/10.3390/d14070540
Submission received: 24 May 2022 / Revised: 23 June 2022 / Accepted: 29 June 2022 / Published: 5 July 2022
(This article belongs to the Topic Arthropod Biodiversity: Ecological and Functional Aspects)

Abstract

:
Microbiota associated with insects play several important roles in their host, including protection against pathogens, provision of nutrition, and survival in hostile environments. The aim of this work was to identify the bacterial community found in avocado thrips from Northwestern Colombia (Antioquia department) in order to find isolates for potential biocontrol purposes. Culture-dependent methods based on 16S rRNA and gyrase B gene sequencing in 42 bacterial isolates allowed the identification of the genera Bacillus, Serratia, Moraxella, Pantoea, and Sphingomonas. Microbial diversity detected with the temperature gradient gel electrophoresis (TGGE) technique on three morphotypes of thrips, named brown (Scirtothrips hansoni), black (Frankliniella panamensis), and pale (Frankliniella sp.), showed a low bacterial community density (Shannon–Wiener index = 1480, p > 0.05) with significant differences among morphotypes (R = 0.7877, p = 0.0004). Results obtained with Illumina sequencing on the V1–V2 hypervariable region of the subunit 16S rRNA showed a predominant sequence in the brown morphotype (Scirtothrips hansoni) that belongs to the genus Wolbachia. The 16S amplicon analyses were extended to more samples and higher resolution using the V4–V5 hypervariable region. The results showed six additional bacteria phyla, confirming the previous observation for the dominant bacterial groups made in S. hansoni and the detection of the alternation of highly predominant genera among these thrips. Our results demonstrate that endosymbiont such as Wolbachia sp. are part of the microbiota of these pests, thereby indicating the possibility of employing this type of bacterium to improve the management of avocado thrips globally.

1. Introduction

Avocado (Persea americana Mill.) is an evergreen, fruit-yielding tree. The production of this fruit (about six million tons) has significantly increased, and approximately 62% of the global crop is harvested in five main countries, namely, Mexico (33%), the Dominican Republic (10.5%), Peru (7.8%), Indonesia (5.7%) and Colombia (5.1%) [1] However, this crop is severely affected by different pests, including thrips, which are a group of tiny insects that are part of a complex of cryptic species that feed from the avocado epidermis through stylet penetration, thereby sucking out the cell nutrients [2,3]. Thrips are characterized by exhibiting body-color polymorphisms, including yellow, brown, and black patterns. Together with other morphological traits, colors are used as diagnostic characters for identification at the species level [4,5].
In general, insects are favored by the conditions and resources provided by their interactions with viruses, bacteria, and fungi [6]. It is known that bacterial communities play important roles in insects, including assisting in growth, fitness, nutrition, development, protection against pathogens and parasitoids, virus transmission, and survival in hostile environments. Many of these interactions are symbiotic and some are achieved through toxin metabolization, and degradation of phytotoxins and pesticides by the microbiota [7,8,9]. Most insects are believed to harbor heritable bacterial symbionts [10] that can be pathogenic, mutualistic, or commensal; some of these are required for survival [11]. In addition, insect–microbiota interactions have been related to insects’ neurophysiological regulation, thereby influencing their social behavior [12]. Insect microbiota composition and diversity may be influenced by environmental factors such as soil pH, amount of organic carbon, and total nitrogen, in addition to climatic factors [13].
Previous work undertaken on the microbiota of thrips that attack different crops includes investigation of the presence of Erwinia in Frankliniella occidentalis [14,15,16,17,18] and Thrips tabaci (collected in onion) [19]. Moreover, Pantoea ananatis and Pantoea agglomerans were detected in the Thrips tabaci gut by Dutta et al. [20]. Other studies have demonstrated the presence of Wolbachia in Franklinothrips vespiformis [21], Scirtothrips perseae [22], Echinothrips americanus, Suocerathrips linguis, Gynaikothrips ficoru, Hercinothrips femoralis and Parthenothrips dracaenae [23], Thrips palmi [24], and Hoplothrips carpathicus [25]. Finding the genus Wolbachia in a pest species is important because it increases the potential for the pest’s biological control.
Most microbiota studies of thrips that attack crops of economic importance have applied culture-dependent analyses. Microbial identification has been improved with the availability of next-generation sequencing technology because it detects microbes that are not cultivable and are present in very low amounts in the samples. This facilitates community analyses in the microbes’ natural habitat, with an accurate taxonomic identification and estimation of their relative abundances and diversities [11,26,27]. Given the limited information available on the microbiota of avocado thrips and the importance of the microbiota associated with thrips, this study aimed to combine both culture-dependent and culture-independent analyses to study the diversity of bacterial communities present in avocado thrips in crops located in Neotropical Northwestern Andes (Department of Antioquia, Colombia). This work constitutes the first report of research on the bacterial community in avocado thrips, and provides the first step towards biological and genetic control strategies for the improvement in the integrated pest management of thrips.

2. Materials and Methods

2.1. Insect Collection

Adult thrips were collected in eastern Antioquia from several commercial avocado crops. Specimen identification was performed according to the thrips’ morphotypes, including the coloration and various morphological features found in the taxonomic keys provided by Mound and Hoddle [4] and Mound and Marullo [5].
Adult specimens of leaves and flowers of two varieties of avocado were collected: Hass (n = 304 samples) and Fuerte (n = 77 samples). The sampling was carried out in twelve different farms: Fruty Green (FG) (6°3′ N 75°30′ W), La Paz (LP) (5°59′18″ N 75°30′32″ W), Luisines (LS) (6°4′25.262″ N 75°29′50.2″ W), El Vaticano (EV) (6°8′53″ N 75°20′32″ W), La Felisa (LF) (6°8′1.6″ N 75°20′01.8″ W), Villa Fátima (VF) (6°17′32.57″ N 75°24′58.61″ W), La Vega (LV) (6°3.73′0″ N 75°24.97′0″ W), Los Molinos (LM) (6°00′20.6″ N 75°23′32.4″ W), San Sebastián (AASS) (5°53′18″ N 75°25′29″ W), Barahonda (BH) (6°5′9″ N 75°26′40″ W), El Paraíso (EP) (5°57′6″ N 75°31′25″ W), and Aquella (AQ) (6°3′46.3″ N 75°20′54.4″ W). These farms are located in different municipalities of Antioquia, including El Retiro, Marinilla, Guarne, La Ceja, Abejorral, Río Negro, Montebello, and Carmen de Viboral (Figure 1). The collections were carried out with oral aspirators [28] during the months of March and August 2017. The insects were transported live to the Cellular and Molecular Biology laboratory at Universidad Nacional de Colombia, headquartered in Medellin, and stored at 4 °C for later processing. A total of 731 thrips were collected, and 381 thrips were used for culture-dependent analysis and 350 thrips were used for culture-independent analysis.

2.2. Culture-Dependent Assays

2.2.1. Isolation and Purification of Bacteria from Thrips

Aerobic and cultivable bacteria were isolated from 381 adult specimens that were immediately processed after collection (Table S1). The samples were washed with 0.5% hypochlorite for 10 s then immersed in 70% ethanol for 10 s and washed three times in sterile water for 10 s to remove excess dust and exogenous bacteria. Each sample included 1 or 25 thrips that were homogenized to produce a pool used for serial dilutions in a 1 X Phosphate Saline Buffer (PBS) and further plated for bacterial growth on Luria–Bertani (LB) agar (BD Difco, East Rutherford, NJ, USA) and medium R2A agar (Merck KGaA, Darmstadt, Germany; [29,30]). Petri dishes were aerobically incubated at 25 °C for 24 to 48 h for LB and 7 days for R2A. LB and R2A media were chosen as non-selective media to promote the growth of a greater diversity of cultivable microorganisms. The total number of colony-forming units (CFUs) was estimated for each homogenate, and different colonies were selected based on the macroscopic (colony size, pigmentation, elevation) and microscopic characteristics by Gram stain. These were isolated from the different media plates and sub-cultured to obtain pure bacterial cultures. All isolates were cryopreserved in 20% glycerol at –80 °C. Comparisons between the mean number of CFUs between the Hass and Fuerte varieties were made using a t-test using the Past 3.0 software [31]. The t-test was estimated after performing a Kolmogorov–Smirnov normality test and a variance homogeneity test.

2.2.2. Extraction of Genomic DNA and PCR of the Transcribed Internal Spacer-ITS, 16S rRNA and gyrB Genes

For the amplification of ribosomal genes, bacterial genomic DNA was extracted from the selected bacteria using a heat lysis protocol [32] from freshly cultured cells. This DNA was later used as a template for PCR assays and was diluted at 1:10 in sterile water. Template DNA from all isolates was initially evaluated by amplifying the inter-ribosomal region (ITS) using primers L1 (5′-CAAGGCATCCACCGT-3′) and G1 (5′-GAAGTCGTAACAAGG-3′) [33], as previously described by Moreno et al. [34]. The amplified bands were run in polyacrylamide gels (8%) and then the software GelCompar II (Applied Biosystems Maths, Belgium) was used to construct a dendrogram, based on the obtained band patterns, using Pearson correlation distance matrices and complete-linkage grouping methods [35]. Fifty percent or more similarity values between ITS standards were established as criteria for subsequent molecular identification.
The DNA isolated from selected colonies was used to amplify the 16S rRNA gene (1.5 kbp), using the primers Eubac 27 F (5′-AGAGTTTGATCCTGGCTCAG-3′) and 1492R (5′-GGTTACCTTGTTACGACTT-3′) according to Moreno et al. [34]. This gene was further sequenced for the identification of each isolate affiliation. Additionally, bacterial DNA was extracted for amplification of the gyrB gene using the UltraClean Microbial DNA Isolation Kit (MOBIO) following the manufacturer’s instructions. Primers UP1 (5′-AGC AGG GTACGGATGTGCGAGCCRTCNACRTCNGCRTCNGTCAT-3′) and UP2R (5′-GAAGTC ATCATGACCGTTCTGCAYGCNGGNGGNAARTTYGA-3′) [36] were used in a PCR reaction to amplify a conserved region of approximately 1.26 kbp of the gyrB gene [37]. This locus was also used to corroborate the identification of each affiliation of the isolates. For this gene, the following thermal profile was used: 1 cycle of 94 °C for 5 min; 30 cycles of: 94 °C for 1 min, 53 °C for 1 min, 72 °C for 90 s; and a final cycle of 72 °C for 7 min. DNA positive (DNA from a pure culture of Bacillus cereus) and negative (ultrapure water) controls were routinely included in all PCR reactions. Amplification products were verified by visualization in 1% agarose gel and purified using Illustra GFX PCR DNA and a Gel Band Purification Kit (GE Healthcare, UK), and the double-stranded DNA was sequenced in both directions using the ABI PRISM 3700 DNA analyzer service of Applied Biosystems.

2.3. Culture-Independent Assays

After thrips collection, groups of 25 individuals were used for DNA extraction. These individuals were differentiated based on body coloration in three main morphotypes: S. hansoni (brown morphotype), Frankliniella (pale morphotype), and F. panamensis (dark morphotype) [4,5]. These samples were labeled according to sampling site as follows: six samples of S. hansoni (brown morphotype) collected in the farms La Vega (LV, LV2), Villa Fátima (VFm and VFp), and San Sebastian (AASS1 and AASS2); six samples belonging to Frankliniella (pale morphotype) collected on Barahonda (BH1), Luisines (LS1), El Paraíso (EP1, EP2), Los Molinos (LM), and La Vega (LV3 and LV5); and one sample of F. panamensis (dark morphotype) collected in La Vega (LVn). The pooled material was homogenized according to the three morphotypes and split into two aliquots, one for PCR-TGGE analysis and the other for NGS analysis using the amplification of two hypervariable regions V1–V2 and V4–V5, as explained below. Total DNA was extracted using chemical lysis with 50 mM Tris-HC1 buffer (pH 8.5), 1 mM EDTA, 0.5% SDS, 200 µg/mL, and proteinase K [38] followed by purification with a plus “FastDNA™ SPIN KIT” (MP Biomedicals).

2.3.1. Temperature Gradient Gel Electrophoresis (TGGE)

The TGGE protocol was used to monitor the dynamic changes in the microbe population in the three morphotypes. Total DNA samples were quantified using a ND-100 nanodrop Thermo Scientific spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). Subsequently, these samples were subjected to amplification of the 16S rRNA gene (variable regions V3–V6). These samples were amplified with forward primer 341 F (5′-CCTGCAGGAGGCAGCAG-3′), with an additional GC termination at the 5′ end and the reverse primer 907R (5′-CCCTGACGTGTTATTCAATTCY-3′). PCR reaction conditions were implemented as previously described for 16S rRNA gene amplification [39].
Gels were stained with 1X GelGreen® Nucleic Acid Gel Stain for 30 min and the obtained image was digitized using Gel-Compare II 6.6 Software (http://www.applied-maths.com accessed on 8 April 2019). A cluster analysis was constructed through Pearson correlation and complete-linkage clustering [35], to detect differences in diversity and abundance profiles amongst samples. Bands of interest were cut, and the DNA was eluted in 30 µL of ultrapure water. Five microliters of the DNA elution were re-amplified using the same primers, 341 F (without GC tail) and 907R. Additionally, some bands were selected to be cloned with the CloneJET® PCR Cloning Kit (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions. Bands and clones that were successfully re-amplified were sent for sequencing using the methods and procedures described above.

2.3.2. Sequencing and Phylogenetic Analysis

Sequences from PCR products of the 16S rRNA gene, gyrB gene of bacterial isolates and TGGE eluted bands were edited and cured using CLC Genomics Workbench software (Version 7.0.3) (CLC Bio, Aarhus, Denmark). Consensus sequences were contrasted with the GenBank database using the BLASTn algorithm. 16S rRNA sequences were also compared with Ribosomal Database Project (RDP) [40,41]. The homologous sequences from the Genbank database were aligned together with sequences obtained here using the MUSCLE function [42] found in MEGA X [43]. The sequences obtained here were further sent to the GenBank database.
The phylogenetic affiliations amongst sequences 16S rRNA and gyrase were analyzed using a Bayesian tree constructed in MrBayes (Version 3.2.0) [44] with an MCMC chain using 10 million generations and four chains for every 1,000,000 generations. Verification and validation of recombination events and the presence of chimeras were carried out with RDP2 software [41]. Methanococcus aeolicus (DQ136171) and Bacteroides caccae (NZ_CP022412) were used as outgroups for phylogenetic analysis of the 16S rRNA and gyrB genes, respectively. For the 16S rDNA partial sequences obtained by TGGE, a dendrogram was constructed with the neighbor-joining (NJ) and Kimura 2-parameter methods as evolutionary distances; finally, bootstrap analysis was undertaken using 5000 iterations with MEGA X [43].

2.3.3. Analysis of the Bacterial Community by PCR-TGGE

A binary matrix was constructed from the presence/absence of bands evidenced in the TGGE gel in the program PAST (Version 3.23). Principal component analysis (PCA) was used to determine differences amongst samples, using 10000 bootstraps and confidence intervals for the execution that represent 95% of the data [30]. Moreover, a grouping was performed with the unweighted pair group algorithm (UPGMA) and the Dice similarity index, and with 10,000 bootstraps. With the same binary matrix, a unidirectional analysis of similarities (ANOSIM) was performed using the Dice index to estimate the similarity of values between bacterial morphotypes of thrips of the bands obtained by TGGE [45].

2.3.4. 16S Amplicon Next Generation Sequencing (NGS)

The quality of the DNA of each group of thrips separated by morphotype (n = 10), obtained as previously described, was verified by spectrophotometry in a Nanodrop® 2000 (Thermo Fisher Scientific, USA). PCR amplification was performed using primers 27F (5′-AGAGTTTGATCMTGGCTCA-3′) and 338R (5′-GCTGCCTCCCGTAGGAGT-3′) [46] targeting hypervariable region V1 and V2 of the 16S rRNA gene, initially used to analyze three samples belonging to S. hansoni (brown morphotype) coming from two farms: LV (La Vega) and VFm. Further on, as is standard in studies of environmental microbiomes of the Earth Microbiome Project (EMP https://earthmicrobiome.org/protocols-and-standards/16s/ accessed on 17 April 2019), subsequently, the primers 515F (5′-GTGYCAGCMGCCGCGGTAA-3′) [47] and 806R (5′-GGACTACNVGGGTWTCTAAT-3′) [48] targeting the hypervariable region V4–V5 of 16S rRNA genes were used in the following samples as they could provide additional resolution inside the taxonomic groups found in the initial samples and could also detect Archaea. The samples were: three samples of S. hansoni (brown morphotype) were analyzed and taken from samples of farms: LV2 (La Vega), AASS1 (San Sebastián 1 = sampling period 1) and AASS2 (San Sebastián, sampling period 2); three Frankliniella samples (pale morphotype) from the farms: BH1 (Barahonda), EP1 (El Paraiso), and LS1 (Luisines); and a sample of F. panamensis (dark morphotype) from the farm LVn (La Vega). The PCR products were used for library construction using Illumina MiSeq platform 2 × 250 bp. The obtained sequences were paired and processed following the QIIME 1.9.1 pipeline for the V1–V2 region and QIIME 2 for the V4–V5 region [49]. Chimeras were removed and raw sequences were quality filtered, clustered into operational taxonomic units (OTUs) and taxonomically classified using the “Naive Bayesian Classifier” of the Ribosomal Database Project (RDP) classifier against the SILVA 132 database [50].
Sample comparisons were performed in terms of relative abundance, and alpha and beta diversity, using R scripts developed by Pérez-Jaramillo et al. [51] following the procedure of Castañeda-Monsalve et al. [27]. Results were plotted using the Phyloseq R package (V 1.24.2) [52]. Comparisons of bacterial diversities based on the bacterial genera dataset were performed with a chi-square test performed with R and significant p-values were considered after a Bonferroni correction (α = 0.05/28 = 0.0018). Observed OTUs, Shannon, Simpson, and Chao1 diversity metrics were also obtained. One-way ANOVA and Tukey HSD were used to test for significance. Additionally, a Venn diagram was calculated using the DISPLAYR web platform (https://www.displayr.com/ accessed on 14 January 2022) to determine the number of specific and shared genera among thrips morphotypes. These analyses were performed on the set of sequences of hypervariable region V4–V5 16S rRNA gene amplicons. For beta diversity calculations, OTU table normalizations were undertaken as described by Pérez-Jaramillo et al. [51]. Differences were visualized using principal coordinate analysis (PCoA), and constrained analysis of principal coordinates (CAP) [53]. Permutational multivariate analysis of variance (PERMANOVA) was used for hypothesis testing; differences between factors were considered to be significant if the likelihood of the observed statistics was <5% (p < 0.005).

3. Results

3.1. Bacterial Diversity through Culture-Dependent Assays

3.1.1. Isolation and Culture of Microbial Isolates

Bacterial growth was followed in 381 adult thrips that were processed in groups of 15 insect pools containing from 1 to 25 specimens (Table S1). Mean bacterial count ranged from 2.6 × 104 to 3.64 × 104 colony forming units (CFU/mL) for Hass and Fuerte varieties, respectively. No significant differences were found in these counts for the two avocado varieties studied (Hass vs. Fuerte) (Table S1).
After examining the macro and microscopic characteristics with Gram staining, 42 isolates were selected and further used for molecular characterizations based on an analysis pattern of the transcribed internal spacer region (ITS) (Figure S1). Profiles of ITS region showed significant similarities between isolates, with common bands between 700–900 bp and 400–600 bp. The dendrogram differentiated 29 clusters with a similarity percentage of 50%. Representative isolates of each cluster were further selected for 16S rRNA gene sequencing. No differences were found in bacterial isolates between Hass and Fuerte varieties, meaning that both share similar bacterial isolates (Table S1, Figures S1 and S2).

3.1.2. Identification of Bacterial Isolates Using 16S rDNA and gyrB Sequencing

Two different predominant phyla were identified: Proteobacteria and Firmicutes, with high percentages of similarity (99–100%; see Table S2). Furthermore, genera Serratia, Pantoea, Sphingomona, and Moraxella were the most abundant groups observed. Ten different bacterial species were identified (Figure S2, Table S2); five were Gram-negative: Serratia liquefaciens, Pantoea cypripedii, Pantoea agglomerans, Sphingomonas olei, and Moraxella osloensis, and five were Gram-positive, corresponding to the Bacillus genus.

3.1.3. Phylogenetic Analysis with Sequences of the 16S rRNA and gyrB Genes

Sequence-based molecular identification of the 16S rRNA and gyrB genes showed a high level of correspondence with sequences deposited in the NCBI at the species level and between different taxonomic divisions (Table S2). The Bayesian tree obtained with 16S rDNA sequences (Figure 2a) and gyrB (Figure 2b) produced two clusters; the first corresponded to the Bacillus genus and the second to Sphingomonas, Moraxella, Pantoea, and Serratia genera. Similarly, phylogenetic results showed that the Erwinia aphidicola (BFo1) strain was found as a symbiont in F. occidentalis (Chanbusarakum and Ullman, [14]; Facey et al. [18]) and this species has 97% identity with our sample T8H6 (Figure 3a). Furthermore, the sample T3H1_1 exhibits a 99% identity percentage with Pantoea agglomerans also reported in other thrips populations and thrips species such as F. occidentalis (De Vries et al. [54]) and T. tabaci (De Vries et al. [19]).

3.2. Bacterial Diversity through Culture-Independent Methods

3.2.1. PCR-TGGE Analysis

In total, 350 thrips were processed to obtain 14 pools of DNA that were further used for the PCR-TGGE analysis (Table S3). This technique allowed the differentiation of fragment profiles of the gene 16S rRNA of approximately 585 bp (Figure S3), assuming that each band (with a different pattern) represented a different taxonomic unit (OTU). From this gel, 68 unique and common DNA bands were selected based on their intensity and were then sequenced and identified according to their relationship with other reported sequences from GenBank. TGGE patterns observed for each sample were similar between replicates (Figure 3 and Figure S3).
After analyzing band patterns for each profile with GelCompar II® software, two groups were observed with similarity rates of 50% and 60% found in the UPGMA dendrogram. The groups of profiles obtained here clustered according to the three morphotypes of thrips studied, as shown in Figure 3 (Morse and Hoddle, 2006). The first cluster included the samples VFm, VFp, LV2, LV, AASS1, and AASS2 of S. hansoni (brown morphotype), CL of F. panamensis (dark morphotype), and LS1 of Frankliniella sp. (pale morphotype); the second cluster included the samples BH, LV3, LV5, EP1, EP2, and LM (pale morphotype) that corresponded to Frankliniella sp. (Figure S4).
Based on the migration and band intensity obtained from the TGGE profile, the bacterial diversity and richness amongst the samples were estimated. The Shannon–Wiener index was calculated, producing values of 1480 and 15.6404, respectively. Both estimators were low for the thrips studied.
Results obtained with ANOSIM analysis showed significant differences between the band patterns found in the three morphotypes of thrips (R = 0.7877, p = 0.0004) (Figure 3 and Figure S3). Wolbachia and non-cultivable Enterobacteriaceae bacteria prevailed in thrips of the brown morphotype (S. hansoni), non-cultivable Firmicutes bacteria were abundant in the pale morphotype (Frankliniella), and Pantoea and Anaplasma were also abundant in dark morphotype (F. panamensis).

3.2.2. Phylogenetic Affiliation of the Sequences Obtained from the TGGE Bands

For phylogenetic analysis, several bands were selected from profiles representing a bacterial operating taxonomic unit (OTU) (Figure 4). Forty-eight bands and seven clones from a total of sixty-eight bands produced good quality sequences that were then assembled, aligned, and identified with high percentages of similarity (98–100%). These bands correspond to specific bacterial groups sequenced with the 16S rDNA based on public databases (Table S4). The affiliations found showed that the phylum Proteobacteria was dominant, with Wolbachia being the most dominant genus (99%) in all samples; Burkholderia, Pantoea, and Ralstonia were also found (Table S4). In the phylum Firmicutes, the Bacillus genus and non-cultivable Firmicutes were found, and, for the phylum Actinobacteria, the genus Cutibacterium was found (Table S4).

3.2.3. Quality Analysis of 16S rDNA Sequence and Bacterial Identity Assignment by High-Throughput Sequencing-NGS

In general, all samples had C Good’s coverage greater than 98% indicating a high and comparable representation. The V1–V2 region allowed broad detection of 16S rRNA genes of bacterial origin. The samples analyzed with this hypervariable region showed the following characteristics: samples of the brown morphotype (S. hansoni) VFm, VFp, and LV exhibited 25.780, 36.381, and 39.812 reads and 941, 1.203, and 1.275 OTUs, respectively. To extend the simultaneous detection of bacteria and archaea, we also used the V4–V5 hypervariable region of the 16S rRNA gene, showing that samples that belong to the brown morphotype (S. hansoni), AASS1, AASS2, and LV2, had a total of 13.341, 23.504, and 42.623 reads and 366, 423, and 270 OTUs, respectively; the pale morphotype samples (Frankliniella), LS1, EP1, and BH1, had a total of 23.595, 67.764, and 99.843 reads and 352, 429, and 456 OTUs, respectively. Finally, for the only sample that belonged to the dark morphotype (F. panamensis, LVn), a total of 74.017 reads and 311 OTUs were obtained.
OTUs with readings representing less than 0.0001% of the total depth of the samples were filtered to eliminate so-called singletons as possible technical artifacts or low-frequency ribotypes that do not contribute significantly to the composition of the community. After filtering through rarefaction analysis, 24.773 reads remained for the V1–V2 region (Figure S5), assigned to 363 OTUs, and, for the V4–V5 region, 13.000 reads (Figure S6) were assigned to 627 OTUs.

3.2.4. Microbiota Composition

In general, for the V1–V2 region, only one phylum with a relative abundance >1% prevailed in S. hansoni (brown morphotype) and corresponded to Proteobacteria (98 to 99%) (Figure 5a). For the V4–V5 region, six phyla were identified with a relative abundance >1% in S. hansoni, (brown morphotype), F. panamensis (dark morphotype), and Frankliniella (pale morphotype). Within this group, for the S. hansoni samples (brown morphotype), the following groups were identified: sample AASS1 with five phyla: Proteobacteria (87.65%), Firmicutes (6.4%), Cyanobacteria (2.74%), Bacteroidetes (1.51%), and Actinobacteria (1.31%); sample AASS2, six phyla: Proteobacteria (79.97%), Firmicutes (7.48%), Deinococcus-Thermus (4.42%), Cyanobacteria (2.09%), Bacteroidetes (2.11%), and Actinobacteria (3.23%); and the sample LV2 showed five phyla: Proteobacteria (99.03%), Firmicutes (0.17%), Cyanobacteria (0.09%), Bacteroidetes (0.013%), and Actinobacteria (0.50%). Six phyla were found for the F. panamensis sample (dark morphotype): LVn: Proteobacteria (98.89%), Firmicutes (0.19%), Deinococcus-Thermus (0.02%), Cyanobacteria (0.78%), Bacteroidetes (0.06%), and Actinobacteria (0.06%). Finally, for the Frankliniella samples (pale morphotype): LS1 sample, five phyla were found: Proteobacteria (99.85%), Firmicutes (0.03%), Cyanobacteria (0.08%), Bacteroidetes (0.03%), and Actinobacteria (0.02%); sample BH1, six phyla: Proteobacteria (59.65%), Firmicutes (1.31%), Deinococcus-Thermus (0.23%), Cyanobacteria (37.81%), Bacteroidetes (0.17%), and Actinobacteria (0.59%); and EP1 showed five phyla: Proteobacteria (97.58%), Firmicutes (0.03%), Cyanobacteria (0.08%), Bacteroidetes (0.29%), and Actinobacteria (0.05%) (Figure 5b). The most common phyla found in all morphotypes were Proteobacteria, followed by Cyanobacteria, which were more common in Frankliniella (pale morphotype), and Firmicutes, which was also abundant in S. hansoni (brown morphotype). It is important to note that Euryarchaeota were detected in both pale and brown morphotypes.
With the V1–V2 region, two genera with a threshold >2% were identified in a sample of S. hansoni (brown morphotype). These genera were Wolbachia (93.98%) and Rosenbergiella (1.84%), while in the other samples of the same morphotypes Wolbachia prevailed with abundances of 95.89% and 92.85%, respectively (Figure 5c). Our results indicate a prevalence of the Wolbachia genus in samples of the S. hansoni species.
With the V4–V5 region, at the genus level, 28 genera were identified with a threshold >2%. This result suggests a great variation in the microbial communities amongst the three morphotypes (Figure 5d, Table S5). The most abundant and common genera found were: Wolbachia (except for a pale EP1 sample), Ralstonia, Phaseolibacter, and Enterobacter (mainly in the pale morphotype), and Streptococcus and Sphingomonas (mainly in the pale and brown morphotypes). Furthermore, the Halorubrum genus (Archaea) was present in both brown and pale morphotypes with abundances of 0.02%. Additionally, a species-specific association was found between S. hansoni (brown morphotype) and Wolbachia genus, in addition to Ehrlichia in F. panamensis (dark morphotype) and Enterobacter and Rickettsia genera in Frankliniella (pale morphotype).
Statistical chi-square analysis showed that Wolbachia, Sphingomonas, Parvibaculum, and Deinococcus-Thermus genera were significantly more abundant in S. hansoni (brown morphotype); Thalassospira, Enterobacter, and other bacteria in Frankliniella (pale morphotype); and, finally, Ehrlichia genus in F. panamensis (dark morphotype) (Table S5). Additionally, in the Venn diagram, it was found that three OTUs were shared for the three morphotypes of thrips belonging to the Enterobacter, Ralstonia, and Pseudomonas genera (Figure 5e, Table S5).

3.2.5. Bacterial Diversity

For the region V1–V2, alpha diversity analysis was performed after rarefaction, taking all samples to a depth of 24,773 sequences/sample; rarefaction curves reached a plateau between 22,000 to 25,000 reads. For the region V4–V5, alpha diversity analysis was performed after rarefaction, taking all samples to a depth of 13,000 sequences/sample; rarefaction curves reached a plateau between 12,000 to 13,000 reads.
Shannon and Simpson indices for the V1–V2 region were estimated to assess bacterial diversity at the morphotype level in S. hansoni (brown morphotype) (Figure 6a). Low diversity was observed, and it was not possible to perform an ANOVA since the three samples evaluated belonged to the same brown morphotype (S. hansoni) (single level). The estimated diversities for S. hansoni (brown morphotype) samples from LV, VFm, and VFp with the V1–V2 region were as follows: Shannon = 0.74 +/−0.11, Simpson = 0.19 +/−0.03. These indices were estimated in two farms: La Vega (LV) located in the municipality of La Ceja and Villa Fátima (VFm and VFp) located in the municipality of Guarne (Figure 6a).
Similarly, Shannon and Simpson indices for the V4–V5 region were evaluated to deepen the assessment of the bacterial diversity amongst S. hansoni (brown morphotype) and to expand the analyses in Frankliniella (pale morphotype) and F. panamensis (dark morphotype). Statistical differences were observed amongst them as the diversity was the highest in S. hansoni (brown morphotype) and the lowest in the dark morphotype. Alpha diversity indices estimated for the V4–V5 region showed that Shannon and Simpson indices were the highest for the brown morphotype (Shannon mean = 1.79 +/−1.05, Simpson mean = 0.56 +/−0.23) and the lowest for the dark morphotype (Shannon mean = 0.27, Simpson mean = 0.08) (Figure 6b). It is important to mention that, for the dark morphotype, only one sample was considered. This may reflect the estimates of low bacterial communities obtained for this specimen. ANOVA results based on the differences in these three indices between the three morphotypes were not significant (Table S6). This result may be due to the normalization of the samples since the AASS2 sample presented 13.000 reads (sample with the lowest reads) and this sample represented the threshold used for this analysis.
Principal coordinate analysis (PCoA) obtained with Bray–Curtis distances was only performed to compare the bacterial communities identified on the three morphotypes with the V4–V5 region. This analysis showed that the variable “Morphotype” explained 49.7% of the β-diversity of the microbial communities associated with avocado trips (PERMANOVA, p = 0.05, Table S7) (Figure 6c).

4. Discussion

In this research, we reveal the composition of the thrips microbiota using a culture-dependent approach with a taxonomic identification carried out with sequencing of the 16S rRNA and gyrase genes combined with culture-independent techniques supported by a next-generation sequencing analysis. A low diversity was observed both at the culture-dependent and culture-independent level and differences were observed between the three morphotypes evaluated and identified as S. hansoni (brown morphotype), F. panamensis (dark morphotype), and Frankliniella (pale morphotype).
With culture-dependent approaches, Proteobacteria and Firmicutes were identified as the predominant phyla. Although the culture-independent techniques showed the presence of the previously found phyla Proteobacteria and Firmicutes, plus Actinobacteria by TGGE, high-throughput sequencing of the 16S rDNA extended the detection of the microbial community to six phyla, mainly composed of Proteobacteria, with a minor fraction shared between representatives of Cyanobacteria, Actinobacteria, Deinococcus-Thermus, Bacteroidetes, and Firmicutes. With all methodological strategies used, Proteobacteria appeared as the dominant phylum in the thrips morphotypes. Although analyses of the microbiota of the order Thysanoptera are limited, especially in those that attack avocado, recent studies confirmed that Proteobacteria is the most abundant in the bacterial communities of other thrips species, such as S. dorsalis, Thrips palmi, and Hoplothrips carpathicus. Proteobacteria comprise ~46 to ~60% of the bacteria detected in thrips [7,11,25,55]. In contrast, Firmicutes have been found in low abundance and comprise no more than 10% [25].
Through culture-dependent methods we obtained 42 bacterial isolates according to their in vitro morphological characterization, which was selected for 16S rDNA sequencing. Additionally, 25 strains from this group were used for the analysis of the gyrB gene sequence to corroborate the identity of bacteria with both genes (Table S1, Figure 2). The results revealed the presence of Bacillus, Serratia, Moraxella, Pantoea, and Sphingomonas genera (Table S2, Figure 3a,b). The Gram-negative Pantoea and Serratia genera have been previously identified in Frankliniella occidentalis (Pergande) and Thrips tabaci (Lindeman) associated with other crops, such as flowers and tobacco [14,16,19]. The main Gram-negative bacteria previously reported in these two species of thrips were the genera Agrobacterium, Erwinia, Methylobacterium, Pseudomonas, and Serratia [7,14,17,19].
In more recent studies, symbiotic relationships between thrips and Erwinia have been described [14,19]. Erwinia was detected in several countries and in different crops, such as beans, cucumbers, and chrysanthemums, suggesting that this genus is transmitted horizontally. Specifically, two different types of bacteria, BFo-1 and BFo-2, were identified in two populations of F. occidentalis [14]. Genome sequencing revealed that BFo-1 is more closely related to species of the Erwinia genus and is a close relative of Erwinia aphidicola [18]. Furthermore, the BFo-2 type bacterium is an enterobacterial microorganism and is not phylogenetically close to other currently known bacterial species. The BFo-2 strain detected here shares a common ancestor with Pantoea and Erwinia genera [18]. The results obtained in the present work showed a value of 97% of BFo-1 with our sample T8H6 (Figure 2a), whereas the sample T3H1_1 had a 99% identity percentage with Pantoea agglomerans. Moreover, bands obtained by culture-independent methods showed the presence of non-cultivable bacteria sequences that are closely related to BFo-2 (Figure 5), suggesting that avocado thrips harbor both bacteria also identified in F. occidentalis [14,18]. All these data suggest that these species share some common microbiota regardless of the species or the crop that they live on.
Pantoea is an interesting Gram-negative genus detected in the present study. This genus has been extensively analyzed in various insects [56] and can degrade herbicides [57] and produce antimicrobials [58]. Most importantly, Pantoea is known to be a symbiote of various insects, as it helps to degrade cellulose in plants (Suen et al., 2010). Other Gram-negative species that were identified in this study include Sphingomonas olei, M. osloensis, and S. liquefaciens. These bacteria are capable of degrading aliphatic hydrocarbons [59], molluscicides [60], and insecticides [61].
Regarding the Gram-positive bacteria found here, several Bacillus strains were identified. Identification of Bacillus spp. using the 16S rRNA gene did not allow adequate identification at species level due to the high percentage of gene sequence similarities amongst closely related species [62]. However, identification using the gyrase gene gave an approximation of the possible Bacillus species, suggesting the presence of B. safensis, B. cereus, B. thuringiensis, B. velezensis, and B. amyloliquefaciens. These bacteria have the potential to be biocontrollers of plant pathogens, since they can produce lipopeptides [63]. B. amyloliquefaciens, B. cereus, B. pumilus, B. velezensis, and B. thurigiensis, amongst other species, can produce a significant reduction in the incidence of diseases in various host plants [63,64,65].
Regarding the culture-independent analysis, differences in microbiota composition and diversity were further evaluated in three different thrips. This analysis involved the extraction of biological material from the entire insect, which means that allochthonous and indigenous microbiota could be included in the DNA samples. Data showed that the three morphotypes shared phylotypes, but with stark differences between them (Figure 4 and Figure 5). This was additionally evidenced by a principal component analysis from the binary matrix of the banding profiles of the TGGE, which also showed significant differences between the three morphotypes (Figure 4) due to results obtained with the similarity indexes and the grouping UPGMA (Figure S3) and one-way ANOSIM. In culture-independent assays, and TGGE and metataxonomic methods, Proteobacteria were found most abundantly, particularly in S. hansoni (brown morphotype), a prevalence of the Wolbachia genus, and non-cultivable Enterobacteriaceae. Non-cultivable Firmicutes bacteria were found in Frankliniella (pale morphotype), and Pantoea and Anaplasma genera in F. panamensis (dark morphotype). Other genera, such as Burkholderia, Ralstonia, Bacillus, Cutibacterium, and Heliothrips haemorrhoidalis symbiotes were mainly present in Frankliniella (pale morphotype). These kinds of symbiotes have been previously described in other insects such as H. haemorrhoidalis, Telenomus tridentatus, and Aedes spidentatus [66,67,68].
Results obtained by high-throughput sequencing showed that Proteobacteria was the most abundant phylum found in the three morphotypes with a relative abundance >95 for the V1–V2 and V4–V5 regions of the 16S rRNA gene. This result corroborated previous studies of other species of thrips [7,11,25] and other insects [69,70], demonstrating the prevalence of this phylum in insects. Similarly, the bacterial communities detected with the V4–V5 region showed the presence of Firmicutes, Actinobacteria, Cyanobacteria, Bacteroidetes, and Deinococcus-Thermus: these phyla have also been previously reported in other species of thrips [7,11,25,71].
Furthermore, the relative abundance of the different genera detected here showed the prevalence of Wolbachia in S. hansoni (brown morphotype) of >92% with the V1–V2 region and >65% with the V4–V5 region, whereas the prevalence of this endosymbiont was lower for Frankliniella (pale morphotype), ranging from 0 to 2.43%, and F. panamensis (dark morphotype), equal to 0.16%.
Wolbachia has been previously identified in other species of thrips. For example, a study by Arakaki et al. [21] found Wolbachia supergroup B in Franklinothrips vespiformis microbiota. Moreover, Dickey et al. [7] detected this endosymbiont in S. dorsalis. Furthermore, Saurav et al. [24], observed that the female and male proportions found in the natural populations of T. palmi were 3:1 and identified Wolbachia in the abdomen of the insect. They suggested that differences in the proportions of the sexes were due to an incompatibility cytoplasmatic effect of this bacteria between genera in this insect. Another study by Kaczmarczyk et al. [25] found that, in H. carpathicus, the Wolbachia genus was present in 70% of the pupae and in 56% of adults. Ambika and Rajagopal [72] also found that Wolbachia was associated with Plicothrips apicalis in the abdomen of their larvae and adults. Additionally, observations made in this study revealed a higher proportion of females relative to males in nature, particularly for Frankliniella, whereas, for S. hansoni, no males were collected in the field. The presence of Wolbachia may be related to the differences in male and female proportions in these two species. Moreover, reproduction by parthenogenesis has been reported in thrips [21] and this may be due to an effect of Wolbachia. Detection of this bacteria is relevant in this study as further investigations should be undertaken to determine its prevalence in avocado thrips and also whether it can be used for biological control.
In the case of the Frankliniella sp. (pale morphotype), the most abundant genus was Enterobacter sp. To date, this bacterium has been reported as part of the S. dorsalis microbiota [7] but not in amounts as high as those found in this work. This bacterium has been associated with probiotic activities in the Mediterranean fly Ceratitis capitata, improving mating competitiveness, flight capacity, and longevity under starvation [73,74]. Another bacterial genus found here was Ehrlichia in the F. panamensis (dark morphotype). The presence of Ehrlichia suggests that these bacteria may be endosymbiont in pests. To date, no studies have shown the existence of this genus in thrips; however, Ehrlichia chaffeensis has been reported in other insects such as ticks, particularly in Amblyomma americanum.
Chao indices obtained from the V1–V2 region of the S. hansoni (brown morphotype) suggest that the richness of the bacterial communities is low. Furthermore, Shannon and Simpson indices were lower than Chao1, suggesting that there is community uniformity and OTU dominance in the S. hansoni samples used in this study. The evaluation of the microbiota diversity in the three morphotypes of thrips with the sequencing of the V4–V5 region of the 16S rRNA gene showed that there are no significant differences in alpha diversity indices among S. hansoni (brown morphotype), Frankliniella (pale morphotype), and F. panamensis (dark morphotype); however, regarding beta diversity and composition, the differences were very evident and strong. Similar studies have been conducted in T. tabaci suggesting that bacterial communities differ by the host [11]. Another study by Kaczmarczyk et al. [25] in H. carpathicus showed that the microbial communities varied according to the development of the life cycle, since they mostly found Wolbachia in the pupa stages.
Overall, the results obtained by culture-independent studies corroborate the findings from the culture-dependent analysis, where the presence of the Firmicutes and Proteobacteria phyla was evidenced, in addition to that of the genera Bacillus, Pantoea, Erwinia, and Sphingomonas. It is noteworthy that, due to the depth of the NGS sequencing, more genera were observed than those evidenced by culture-dependent analysis. In addition, the genus Wolbachia was found, which is an intracellular bacterium that generates symbiosis with the host and does not grow in culture medium. However, in the present study, anaerobic bacteria were not analyzed.

5. Conclusions

This work describes bacterial diversity present in avocado thrips species from Northeastern Antioquia, in Colombia. Our results complement the microbiota information found in various thrips species, including Frankliniella sp., F. panamensis, S. hansoni; F. occidentalis, T. tabaci, F. tritici, and H. haemorrhoidalis from avocado, bean, onion, nightshade, and wheat crops. Using molecular and traditional microbiological methods, a low bacterial diversity was found at the phylum and genus level. The application of this approach resulted in a prevalence of endosymbionts, particularly the Wolbachia, Pantoea, and Sphingomonas genera.
Additionally, this work describes the microbiota associated with three morphotypes of thrips: brown morphotype (S. hansoni), pale morphotype (Frankliniella), and dark morphotype (F. panamensis), found in avocado in Antioquia using next-generation sequencing approaches. Our results corroborate the presence of Wolbachia in the three morphotypes of thrips with prevalence in S. hansoni (brown morphotype). Furthermore, Enterobacter was found in Frankliniella sp. (pale morphotype), and Ehrlichia was found in F. panamensis (dark morphotype). Moreover, a higher α-diversity was estimated in S. hansoni.
The results obtained in this study showed that the bacterial community diversity present in adults of avocado thrips is low, having a high predominance of a few bacterial components. In addition, the microbiota in the morphotypes of thrips according to morphotype showed a large variation in the specimens analyzed. A strong predominance of the endosymbiont Wolbachia was also found here, which has been studied globally as a biological control in insects [75]. Therefore, these results are relevant to complementing the knowledge about the microbiota, as this information may be useful for biotechnological processes, including biological control methods based on symbiont employment or selective targeting.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/d14070540/s1, Figure S1: Dendrogram of Ribosomal Intergenic Spacer analysis (RISA) using complete linkage and Pearson correlation distance metric of bacterial isolates found in thrips; Figure S2: Morphology of bacterial isolates found from avocado thrips; Figure S3: Bands patterns on TGGE gels of the amplified product of the 16S rRNA gene from thrips DNA samples; Figure S4: Similarity clustering dendrogram using Dice-UPGMA of thrips samples; Figure S5: Rarefaction curve obtained for V1–V2 sequencing of the 16S rRNA gene in the brown morphotype collected from Antioquian avocado; Figure S6: Rarefaction curve obtained for V4–V5 sequencing of the 16S rRNA gene in the brown morphotype collected from Antioquian avocado; Table S1: Information on bacterial thrips samples, colony forming units (CFU), municipalities where the samples were taken and variety of avocado used for the tests; Table S2: Identification of bacteria isolated from adult thrips, according to their similarity with the sequences of the 16S RNAr gene. The identification of the species was considered positive when the % similarity was greater than 98%; Table S3: Information on used samples of thrips for independent cultivation, morphotypes, farm and avocado variety; Table S4: Taxonomic classification of bacteria and clones found in adults of thrips using temporal temperature gradient electrophoresis (TGGE) analysis, according to their similarity to the 16S rRNA gene sequences recorded in the GenBank and RDP II databases; Table S5: List of genera and percentages of relative abundance of bacteria identified in the three morphotypes with the V4–V5 region; Table S6: ANOVA test performed for the diversity indices found for the three morphotypes with V4–V5 region; Table S7: PERMAMOVA test used to compare the bacterial community between three morphotypes based on the Bray Curtis distance.

Author Contributions

All authors contributed to the revision of the manuscript. S.O.-G. and L.M.M.-P. worked on the culture-dependent experiments, D.C.-C. contributed to the development of all experiments in the laboratory, the collection trips, the data analyses, and the writing of the manuscript. E.G.-B. worked on the NGS analyses. H.J., C.X.M.-H., C.S.-B. and R.E.A.-I. designed the study, analyzed the data, supervised the experiments, and contributed to writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The research reported in this publication was supported by the Universidad Nacional de Colombia under Project Code 202010013471, contract number FP44842-132-2017. This material is based on work supported by the Ministry of science, technology and innovation (Minciencias) of Colombia Graduate Research Fellowship Program for Daniela Cano-Calle No. 727 in 2015 under Grant No. 201010020475. Funding was also obtained from the Foundation for the Promotion of Research and Technology, Bank of the Republic of Colombia Project No. 201423.

Institutional Review Board Statement

The insect collections were made between April and August 2017 through the certificate issued by Universidad Nacional de Colombia, Sede Medellín with the permission to collect specimens of wild species of biological diversity with scientific research purposes produced by the National Authority for Environmental Licenses (ANLA) to Universidad Nacional de Colombia through resolution No. 0255 of 14 March 2014 (article 3). Thrips were collected on private property and permission was received from the landowners before sampling procedures.

Data Availability Statement

The datasets supporting the conclusions of this article are included within the article and its additional file. The newly generated sequences are submitted to the GenBank database under accession numbers MZ222250-MZ222267 for 16S rRNA sequence, MZ233678-MZ233724 TGGE band, MZ230003-MZ230010 Clones and ON012588-ON012603 for gyrase sequence. The source data of the raw Illumina sequence reads of 16S amplicons from thrips samples can be found in the Sequence Read Archive (SRA) under BioProject ID: PRJNA837112 available at http://www.ncbi.nlm.nih.gov/bioproject/837112 accessed on 30 June 2022.

Acknowledgments

We are thankful to Roberto Aguilar and to the laboratorio de biología celular y molecular of the Universidad Nacional de Colombia, Sede Medellín. We thank the owners of the farms who permitted for the collection of biological material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Salazar-López, N.J.; Domínguez-Avila, J.A.; Yahia, E.M.; Belmonte-Herrera, B.H.; Wall-Medrano, A.; Montalvo-González, E.; González-Aguilar, G.A. Avocado fruit and by-products as potential sources of bioactive compounds. Food Res. Int. 2020, 138, 109774. [Google Scholar] [CrossRef] [PubMed]
  2. Cambero-campos, J.; Johansen-Naime, R.; García-Martínez, O.; Cantu-Sifuentes, M.; Cerna-Chavez, E.; Renata-Salazar, A. Especies Depredadoras de Trips (Thysanoptera) asociadas a huertas de aguacate en Nayarit, México. Acta Zoológica Mex. 2011, 27, 115–121. [Google Scholar] [CrossRef]
  3. Monje, B.; Delgadillo, D.; Gómez, J.; Varón, E. Manejo de Neohydatothrips signifer Priesner (Thysanoptera: Thripidae) en maracuyá (Passiflora edulis f. flavicarpa Degener) en el departamento del Huila (Colombia). Cienc. Y Tecnol. Agropecu. 2012, 13, 21–30. [Google Scholar] [CrossRef] [Green Version]
  4. Mound, L.; Hoddle, M. The Scirtothrips perseae species-group (Thysanoptera), with one new species from avocado, Persea americana. Zootaxa 2016, 4079, 388–392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Mound, L.A.; Marullo, R. The thrips of central and south America: An introduction (Insecta: Thysanoptera). Florida Entomol. 1996, 79, 270–271. [Google Scholar] [CrossRef]
  6. Douglas, A.E. Multiorganismal insects: Diversity and function of resident microorganisms. Annu. Rev. Entomol. 2015, 60, 17–34. [Google Scholar] [CrossRef] [Green Version]
  7. Dickey, A.M.; Trease, A.J.; Jara-Cavieres, A.; Kumar, V.; Christenson, M.K.; Potluri, L.P.; Morgan, J.K.; Shatters, R.G., Jr.; McKenzie, C.L.; Davis, P.H.; et al. Estimating bacterial diversity in Scirtothrips dorsalis (Thysanoptera: Thripidae) via next generation sequencing. Florida Entomol. 2014, 97, 362–366. [Google Scholar] [CrossRef]
  8. Douglas, A.E. The microbial dimension in insect nutritional ecology. Funct. Ecol. 2009, 23, 38–47. [Google Scholar] [CrossRef]
  9. Engel, P.; Moran, N.A. The gut microbiota of insects—Diversity in structure and function. FEMS Microbiol. Rev. 2013, 37, 699–735. [Google Scholar] [CrossRef]
  10. Duron, O.; Hurst, G.D. Arthropods and inherited bacteria: From counting the symbionts to understanding how symbionts count. BMC Biol. 2013, 11, 45. [Google Scholar] [CrossRef] [Green Version]
  11. Gawande, S.J.; Anandhan, S.; Ingle, A.; Roylawar, P.; Khandagale, K.; Gawai, T.; Jacobson, A.; Asokan, R.; Singh, M. Microbiome profiling of the onion thrips, Thrips tabaci Lindeman (Thysanoptera: Thripidae). PLoS ONE 2019, 14, e0223281. [Google Scholar] [CrossRef] [PubMed]
  12. Liberti, J.; Engel, P. The gut microbiota—Brain axis of insects. Curr. Opin. Insect Sci. 2020, 39, 6–13. [Google Scholar] [CrossRef] [PubMed]
  13. Huang, S.; Zhang, H. The Impact of environmental heterogeneity and life stage on the hindgut microbiota of Holotrichia Parallela larvae (Coleoptera: Scarabaeidae). PLoS ONE 2013, 8, e57169. [Google Scholar] [CrossRef]
  14. Chanbusarakum, L.; Ullman, D. Characterization of bacterial symbionts in Frankliniella occidentalis (Pergande), Western flower thrips. J. Invertebr. Pathol. 2008, 99, 318–325. [Google Scholar] [CrossRef]
  15. de Vries, E.J.; Jacobs, G.; Sabelis, M.W.; Menken, S.B.J.; Breeuwer, J.A.J. Diet-dependent effects of gut bacteria on their insect host: The symbiosis of Erwinia sp. and western flower thrips. Proc. R. Soc. Lond. B 2004, 271, 2171–2178. [Google Scholar] [CrossRef] [Green Version]
  16. De Vries, E.; Breeuwer, J.; Jacobs, G.; Mollema, C. The association of Western flower thrips, Frankliniella occidentalis, with a near Erwinia species gut bacteria: Transient or permanent? J. Invertebr. Pathol. 2001, 77, 120–128. [Google Scholar] [CrossRef]
  17. de Vries, E.J.; Jacobs, G.; Breeuwer, J.A.J. Growth and transmission of gut bacteria in the western flower thrips, Frankliniella occidentalis. J. Invertebr. Pathol. 2001, 77, 129–137. [Google Scholar] [CrossRef] [PubMed]
  18. Facey, P.D.; Méric, G.; Hitchings, M.D.; Pachebat, J.A.; Hegarty, M.J.; Chen, X.; Morgan, L.V.A.; Hoeppner, J.E.; Whitten, M.M.A.; Kirk, W.D.J.; et al. Draft genomes, phylogenetic reconstruction, and comparative genomics of two novel cohabiting bacterial symbionts isolated from Frankliniella occidentalis. Genome Biol. Evol. 2015, 7, 2188–2202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. de Vries, E.J.; Van Der Wurff, A.G.; Jacobs, G.; Breeuwer, J. Onion thrips, Thrips tabaci, have gut bacteria that are closely related to the symbionts of the western flower thrips, Frankliniella occidentalis. J. Insect Sci. 2008, 8, 23. [Google Scholar] [CrossRef] [Green Version]
  20. Dutta, B.; Barman, A.K.; Srinivasan, R.; Avci, U.; Ullman, D.; Langston, D.B., Jr.; Gitaitis, R. Transmission of Pantoea ananatis and Pantoea agglomerans causal agents of center rot of onion (Allium cepa L.) by onion thrips (Thrips tabaci Lindeman). Phytopathology 2014, 104, 812–819. [Google Scholar] [CrossRef] [Green Version]
  21. Arakaki, N.; Miyoshi, T.; Noda, H. Wolbachia-mediated parthenogenesis in the predatory thrips Franklinothrips vespiformis (Thysanoptera: Insecta). Proc. R. Soc. B Biol. Sci. 2001, 268, 1011–1016. [Google Scholar] [CrossRef] [Green Version]
  22. Weeks, A.R.; Velten, R.; Stouthamer, R. Incidence of a new sex-ratio-distorting endosymbiotic bacterium among arthropods. Proc. R. Soc. B Biol. Sci. 2003, 270, 1857–1865. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kumm, S.; Moritz, G. First detection of Wolbachia in arrhenotokous populations of thrips species (Thysanoptera: Thripidae and Phlaeothripidae) and its role in reproduction. Environ. Entomol. 2008, 37, 1422–1428. [Google Scholar] [CrossRef] [PubMed]
  24. Saurav, G.; Daimei, G.; Singh, V.R.; Popli, S.; Rajagopal, R. Detection and localization of Wolbachia in Thrips palmi Karny (Thysanoptera: Thripidae). Indian J. Microbiol. 2016, 56, 167–171. [Google Scholar] [CrossRef] [Green Version]
  25. Kaczmarczyk, A.; Kucharczyk, H.; Kucharczyk, M.; Kapusta, P.; Sell, J.; Zielińska, S. First insight into microbiome profile of fungivorous thrips Hoplothrips carpathicus (Insecta: Thysanoptera) at different developmental stages: Molecular evidence of Wolbachia endosymbiosis. Sci. Rep. 2018, 8, 14376. [Google Scholar] [CrossRef] [PubMed]
  26. Cao, Y.; Fanning, S.; Proos, S.; Jordan, K.; Srikumar, S. A review on the applications of next generation sequencing technologies as applied to food-related microbiome studies. Front. Microbiol. 2017, 8, 1829. [Google Scholar] [CrossRef]
  27. Castañeda-Monsalve, V.A.; Junca, H.; García-Bonilla, E.; Montoya-Campuzano, O.I.; Moreno-Herrera, C.X. Characterization of the gastrointestinal bacterial microbiome of farmed juvenile and adult white Cachama (Piaractus brachypomus). Aquaculture 2019, 512, 734325. [Google Scholar] [CrossRef]
  28. Martin, J.E. comp. Part 1. Collecting, preparing, and preserving insects, mites and spiders. In The Insects and Arachnids of Canada; Agriculture Canada: Ottawa, ON, Canada, 1976; pp. 9–175. ISBN 0660016508. [Google Scholar]
  29. Pham, V.H.T.; Kim, J. Improvement for isolation of soil bacteria by using common culture media. J. Pure Appl. Microbiol. 2016, 10, 49. [Google Scholar]
  30. Silva-Bedoya, L.M.; Sánchez-Pinzón, M.S.; Cadavid-Restrepo, G.E.; Moreno-Herrera, C.X. Bacterial community analysis of an industrial wastewater treatment plant in Colombia with screening for lipid-degrading microorganisms. Microbiol. Res. 2016, 192, 313–325. [Google Scholar] [CrossRef]
  31. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 2001, 4, 9. [Google Scholar]
  32. Sanders, E.R.; Miller, J.H. I, Microbiologist: A Discovery-Based Course in Microbial Ecology and Molecular Evolution; ASM Press: Los Angeles, CA, USA, 2010; ISBN 978-1-55581-470-0. [Google Scholar]
  33. Jensen, M.A.; Webster, J.A.; Straus, N. Rapid identification of bacteria on the basis of polymerase chain reaction-amplified ribosomal DNA spacer polymorphisms. Appl. Environ. Microbiol. 1993, 59, 945–952. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Moreno, C.; Romero, J.; Espejo, R.T. Polymorphism in repeated 16S PRNA genes is a common property of type strains and environmental isolates of the genus Vibrio. Microbiology 2002, 148, 1233–1239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Yim, O.; Ramdeen, K.T. Hierarchical cluster analysis: Comparison of three linkage measures and application to psychological data. Quant. Methods Psychol. 2015, 11, 8–21. [Google Scholar] [CrossRef] [Green Version]
  36. Yamamoto, S.; Harayama, S. PCR Amplification and direct sequencing of gyrB genes with universal primers and their application to the detection and taxonomic analysis of pseudomonas putida strains. Appl. Environ. Microbiol. 1995, 61, 1104–1109. [Google Scholar] [CrossRef] [Green Version]
  37. Yin, H.; Cao, L.; Qiu, G.; Wang, D.; Kellogg, L.; Zhou, J.; Liu, X.; Dai, Z.; Ding, J.; Liu, X. Molecular diversity of 16S rRNA and gyrB genes in copper mines. Arch. Microbiol. 2008, 189, 101–110. [Google Scholar] [CrossRef]
  38. Goldenberger, D.; Perschil, I.; Ritzler, M.; Altwegg, M. A Simple “Universal” DNA extraction procedure using SDS and proteinase K is compatible with direct PCR AMPLIFICATION. PCR Methods Appl. 1995, 4, 368–370. [Google Scholar] [CrossRef] [Green Version]
  39. Moreno, C.X.; Moy, F.; Daniels, T.J.; Godfrey, H.P.; Cabello, F.C. Molecular analysis of microbial communities identified in different developmental stages of Ixodes scapularis ticks from Westchester and Dutchess Counties, New York. Environ. Microbiol. 2006, 8, 761–772. [Google Scholar] [CrossRef]
  40. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Oxford Univ. Press 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [Green Version]
  41. Cole, J.R.; Wang, Q.; Cardenas, E.; Fish, J.; Chai, B.; Farris, R.J.; Kulam-Syed-Mohideen, A.S.; McGarrell, D.M.; Marsh, T.; Garrity, G.M.; et al. The Ribosomal Database Project: Improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009, 37, 141–145. [Google Scholar] [CrossRef] [Green Version]
  42. Edgar, R.C. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004, 5, 113. [Google Scholar] [CrossRef] [Green Version]
  43. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  44. Ronquist, F.; Teslenko, M.; Van Der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. Mrbayes 3.2: Efficient bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Vivero, R.J.; Gil Jaramillo, N.; Cadavid-Restrepo, G.; Uribe Soto, S.; Moreno Herrera, C.X. Structural differences in gut bacteria communities in developmental stages of natural populations of Lutzomyia evansi from Colombia’s Caribbean coast. Parasit. Vectors 2016, 9, 496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Camarinha-Silva, A.; Jáuregui, R.; Chaves-Moreno, D.; Oxley, A.P.A.; Schaumburg, F.; Becker, K.; Wos-Oxley, M.L.; Pieper, D.H. Comparing the anterior nare bacterial community of two discrete human populations using Illumina amplicon sequencing. Environ. Microbiol. 2014, 16, 2939–2952. [Google Scholar] [CrossRef] [PubMed]
  47. Parada, A.E.; Needham, D.M.; Fuhrman, J.A. Every base matters: Assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 2016, 18, 1403–1414. [Google Scholar] [CrossRef] [PubMed]
  48. Apprill, A.; Mcnally, S.; Parsons, R.; Weber, L. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat. Microb. Ecol. 2015, 75, 129–137. [Google Scholar] [CrossRef] [Green Version]
  49. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.; Costello, E.; Fierer, N.; Gonzalez, A.; Goodrich, J.; Gordon, J.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  50. Wang, Q.; Garrity, G.M.; Tiedje, J.M.; Cole, J.R. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 2007, 73, 5261–5267. [Google Scholar] [CrossRef] [Green Version]
  51. Pérez-Jaramillo, J.E.; Carrión, V.J.; Bosse, M.; Ferrão, L.F.V.; De Hollander, M.; Garcia, A.A.F.; Ramírez, C.A.; Mendes, R.; Raaijmakers, J.M. Linking rhizosphere microbiome composition of wild and domesticated Phaseolus vulgaris to genotypic and root phenotypic traits. ISME J. 2017, 11, 2244–2257. [Google Scholar] [CrossRef] [Green Version]
  52. McMurdie, P.J.; Holmes, S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef] [Green Version]
  53. Oksanen, J.; Blanchet, F.G.; Kindt, R.; Legendre, P.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; Stevens, M.H.M.; Wagner, H. Vegan: Community Ecology Package. Available online: http://vegan.r-forge.r-project.org/ (accessed on 12 February 2020).
  54. De Vries, E.; Van de Wetering, F.; Van der Hoek, M.; Jacobs, G.; Breeuwer, J. Symbiotic bacteria (Erwinia sp.) in the gut of Frankliniella occidentalis (Thysanoptera: Thripidae) do not affect its ability to transmit Tospovirus. Eur. J. Entomol. 2012, 109, 261–266. [Google Scholar] [CrossRef] [Green Version]
  55. Yun, J.; Roh, W.; Whon, W.; Jung, M.; Kim, M.; Park, D.; Yoon, C.; Nam, Y.-D.; Kim, Y.-J.; Choi, J.-H.; et al. Insect gut bacterial diversity determined by environmental habitat, diet, developmental stage, and phylogeny of host. Appl. Environ. Microbiol. 2014, 80, 5254–5264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Walterson, A.M.; Stavrinides, J. Pantoea: Insights into a highly versatile and diverse genus within the Enterobacteriaceae. FEMS Microbiol. Rev. 2015, 39, 968–984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Pileggi, M.; Alvim, S.; Pileggi, V.; Olchanheski, R.; Garbugio, P.; Munoz, A.M.; Koskinen, W.C.; Barber, B.; Sadowsky, M. Isolation of mesotrione—Degrading bacteria from aquatic environments in Brazil. Chemosphere 2012, 86, 1127–1132. [Google Scholar] [CrossRef] [Green Version]
  58. Walterson, A.M.; Smith, D.D.N.; Stavrinides, J. Identification of a Pantoea biosynthetic cluster that directs the synthesis of an antimicrobial natural product. PLoS ONE 2014, 9, e96208. [Google Scholar] [CrossRef] [Green Version]
  59. Kumar, D.; Kim, J. Sphingomonas olei sp. nov., with the ability to degrade aliphatic hydrocarbons, isolated from oil-contaminated soil. Int. J. Syst. Evol. Microbiol. 2017, 67, 2731–2738. [Google Scholar] [CrossRef]
  60. An, R.; Sreevatsan, S.; Grewal, P.S. Moraxella osloensis gene expression in the slug host Deroceras reticulatum. BMC Microbiol. 2008, 8, 19. [Google Scholar] [CrossRef] [Green Version]
  61. Kaibara, F.; Liyama, K.; Chieda, Y.; Lee, J.M.; Kusakabe, T.; Yasunaga-aoki, C.; Shimizu, S. Construction of serralysin-like metalloprotease-deficient mutants of Serratia liquefaciens and their virulence in the silkworm, Bombyx mori. J. Insect Biotechnol. Sericologynsect Biotechnol. Sericology 2012, 81, 55–61. [Google Scholar] [CrossRef]
  62. Wang, L.; Lee, F.; Tai, C.; Kasai, H. Comparison of gyrB gene sequences, 16S rRNA gene sequences and DNA–DNA hybridization in the Bacillus subtilis group. Int. J. Syst. Evol. Microbiol. 2007, 57, 1846–1850. [Google Scholar] [CrossRef] [Green Version]
  63. Shafi, J.; Tian, H.; Ji, M. Bacillus species as versatile weapons for plant pathogens: A review. Biotechnol. Biotechnol. Equip. 2017, 31, 446–459. [Google Scholar] [CrossRef] [Green Version]
  64. Cawoy, H.; Bettiol, W.; Fickers, P.; Ongena, M. Bacillus-based biological control of plant diseases. In Pesticides in the Modern World—Pests Control and Pesticides Exposure; Stoytcheva, M., Ed.; InTech: Rijeka, Croatia, 2011; pp. 273–302. ISBN 978-953-307-459-7. [Google Scholar]
  65. Choudhary, D.K.; Johri, B.N. Interactions of Bacillus spp. and plants—With special reference to induced systemic resistance (ISR). Microbiol. Res. 2009, 164, 493–513. [Google Scholar] [CrossRef] [PubMed]
  66. Nguyen, D.T.; Spooner-hart, R.N.; Riegler, M. Polyploidy versus endosymbionts in obligately thelytokous thrips. BMC Evol. Biol. 2015, 15, 23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. de Ramírez-Ahuja, M.L.; Gómez-Govea, M.A.; Lugo-Trampe, A.; Borrego-Soto, G.; Delgado-Enciso, I.; Ponce-Garcia, G.; Martínez-Fierro, M.L.; Ramírez-Valles, E.G.; Treviño, V.; Flores-Suarezs, A.E.; et al. Microbiota of Telenomus tridentatus (Platygastroidea: Scelionidae ): An unwanted parasitoid. Appl. Entomol. 2019, 143, 834–841. [Google Scholar] [CrossRef]
  68. Shelomi, M. Journal of Asia-Paci fi c Entomology Bacterial and eukaryote microbiomes of mosquito habitats in dengue- endemic southern Taiwan. J. Asia. Pac. Entomol. 2019, 22, 471–480. [Google Scholar] [CrossRef]
  69. Barak, H.; Kumar, P.; Zaritsky, A.; Mendel, Z.; Ment, D.; Kushmaro, A.; Ben-dov, E. Diversity of Bacterial Biota in Capnodis tenebrionis (Coleoptera: Buprestidae) Larvae. Pathogens 2019, 8, 4. [Google Scholar] [CrossRef] [Green Version]
  70. Chua, K.; Song, S.; Yong, H.; See-too, W.-S.; Yin, W.; Chan, K.-G. Microbial community composition reveals spatial variation and distinctive core microbiome of the weaver ant Oecophylla smaragdina in Malaysia. Sci. Rep. 2018, 8, 10777. [Google Scholar] [CrossRef] [Green Version]
  71. Powell, C.M.; Lopez Montiel, A.; Beddingfield, B.; Hanson, J.D.; Bextine, B.R. Comparison of bacterial communities of flower thrips (Frankliniella tritici) and potato psyllid (Bactericera cockerelli). Southwest. Entomol. 2015, 40, 765–773. [Google Scholar] [CrossRef]
  72. Ambika, S.; Rajagopal, R. Lumen anatomy and localization of Wolbachia sp. in the thrips, Plicothrips apicalis (Bagnall). Curr. Sci. 2018, 115, 1297–1304. [Google Scholar] [CrossRef]
  73. Augustinos, A.A.; Kyritsis, G.A.; Papadopoulos, N.T.; Abd-Alla, A.M.M.; Cáceres, C.; Bourtzis, K. Exploitation of the medfly gut microbiota for the enhancement of sterile insect technique: Use of Enterobacter sp. in larval diet-based probiotic applications. PLoS ONE 2015, 10, e0136459. [Google Scholar] [CrossRef] [Green Version]
  74. Kyritsis, G.A.; Augustinos, A.A.; Ntougias, S.; Papadopoulos, N.T.; Bourtzis, K.; Cáceres, C. Enterobacter sp. AA26 gut symbiont as a protein source for Mediterranean fruit fly mass-rearing and sterile insect technique applications. BMC Microbiol. 2019, 19, 288. [Google Scholar] [CrossRef] [Green Version]
  75. Caragata, E.P.; Rocha, M.N.; Pereira, T.N.; Mansur, S.B.; Heverton, L.; Id, C.D.; Id, L.A.M. Pathogen blocking in Wolbachia infected Aedes aegypti is not affected by Zika and dengue virus co-infection. PLoS Negl. Trop. Dis. 2019, 13, e0007443. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Map of the location of the studied area in Colombia, South America.
Figure 1. Map of the location of the studied area in Colombia, South America.
Diversity 14 00540 g001
Figure 2. Bayesian dendrograms of bacterial isolates from thrips: (a) dendrogram with the 16S rRNA gene, (b) dendrogram with the gyrB gene. The dendrograms were generated using the MrBayes program. The bacterial isolates are represented with codes, while the sequences obtained from GenBank have accessions with their respective identifications. The numbers in the nodes correspond to the posterior probability.
Figure 2. Bayesian dendrograms of bacterial isolates from thrips: (a) dendrogram with the 16S rRNA gene, (b) dendrogram with the gyrB gene. The dendrograms were generated using the MrBayes program. The bacterial isolates are represented with codes, while the sequences obtained from GenBank have accessions with their respective identifications. The numbers in the nodes correspond to the posterior probability.
Diversity 14 00540 g002
Figure 3. Information from independent-culture samples. (a) Clustering analysis by the Pearson and complete-linkage method of the PCR-TGGE profiles of the 16S rDNA amplification products (V3–V6 variable regions) using the GelCompar II software. Each band is assumed to represent a taxonomic unit and therefore the profiles reflect the bacterial diversity obtained. Nomenclature: Group 1: isolates of S. hansoni (brown morphotype) that came from the farms La Vega (LV), Villa Fátima (VFm and VFp) farms San Sebastian (AASS1 and AASS2), La Vega (LV2), F. panamensis (morphotype dark) La Vega farm (LVn), and Frankliniella sp. (pale morphotype) Luisines farm (LS1); Group 2: isolates of Frankliniella sp. (pale morphotype); the samples were Barahonda (BH1), Luisines (LS1), El Paraíso (EP1, EP2), Los Molinos (LM), and La Vega (LV3 and LV5). The samples that have the symbol ♦ were the selected samples used for the metataxonomic analysis. (b) Analysis of the main component of the quantitative TGGE matrix of the three thrips morphotypes. S. hansoni (brown morphotype); LV4, VFm, VFp, AASS1, AASS2 and LV2 (□), Frankliniella sp. (pale morphotype); BH, LS, LM, EP1, EP2, LV3 and LV5 (●), and F. panamensis (dark morphotype); CL (◊).
Figure 3. Information from independent-culture samples. (a) Clustering analysis by the Pearson and complete-linkage method of the PCR-TGGE profiles of the 16S rDNA amplification products (V3–V6 variable regions) using the GelCompar II software. Each band is assumed to represent a taxonomic unit and therefore the profiles reflect the bacterial diversity obtained. Nomenclature: Group 1: isolates of S. hansoni (brown morphotype) that came from the farms La Vega (LV), Villa Fátima (VFm and VFp) farms San Sebastian (AASS1 and AASS2), La Vega (LV2), F. panamensis (morphotype dark) La Vega farm (LVn), and Frankliniella sp. (pale morphotype) Luisines farm (LS1); Group 2: isolates of Frankliniella sp. (pale morphotype); the samples were Barahonda (BH1), Luisines (LS1), El Paraíso (EP1, EP2), Los Molinos (LM), and La Vega (LV3 and LV5). The samples that have the symbol ♦ were the selected samples used for the metataxonomic analysis. (b) Analysis of the main component of the quantitative TGGE matrix of the three thrips morphotypes. S. hansoni (brown morphotype); LV4, VFm, VFp, AASS1, AASS2 and LV2 (□), Frankliniella sp. (pale morphotype); BH, LS, LM, EP1, EP2, LV3 and LV5 (●), and F. panamensis (dark morphotype); CL (◊).
Diversity 14 00540 g003
Figure 4. Dendrogram of sequences of the 16S rDNA bands (~400 bp) obtained from the profiles (PCR-TGGE), together with sequences from GenBank. The dendrogram was constructed using the neighbor-joining method and the Kimura 2-parameter model. The numbers on the nodes represent bootstrap values. The robustness of the phylogeny was tested by bootstrap analysis using 1000 iterations.
Figure 4. Dendrogram of sequences of the 16S rDNA bands (~400 bp) obtained from the profiles (PCR-TGGE), together with sequences from GenBank. The dendrogram was constructed using the neighbor-joining method and the Kimura 2-parameter model. The numbers on the nodes represent bootstrap values. The robustness of the phylogeny was tested by bootstrap analysis using 1000 iterations.
Diversity 14 00540 g004
Figure 5. Relative abundance of bacterial OTUs at the phylum level found in (a) samples obtained from the sequencing of the V1–V2 region of S. hansoni. (b) Sequencing of the V4–V5 region of the samples of S. hansoni (brown), F. panamensis (dark morphotype), and Frankliniella (pale morphotype). Data are presented in percentages (%) for all sequences classified at this level of taxonomy. Relative abundance of bacterial OTUs at the genus level found for (c) samples obtained from the sequencing of the S. hansoni V1–V2 region. (d) Sequencing of the V4–V5 region of the samples of S. hansoni (brown), F. panamensis (dark morphotype), and Frankliniella (pale morphotype). Data are presented in percentages (%) for all sequences classified at this level of the taxonomy. (e) Venn diagram shows the “core” at the genus level. Relative abundance of genera based on the data obtained using the metataxonomic approach of the Antioquia thrips samples. Venn diagram showing overlays of OTUs, where 3 OTUs were common for the three thrips morphotypes (dark, pale, and brown) evaluated. Nomenclature: LV, LV2 and LVn correspond to the La Vega farm, VFm and VFp to the Villa Fátima farm, AASS1 and AASS2 to San Sebastian, BH1 to Barahonda, LS1 to Luisines, and EP1 to the farm El Paraíso.
Figure 5. Relative abundance of bacterial OTUs at the phylum level found in (a) samples obtained from the sequencing of the V1–V2 region of S. hansoni. (b) Sequencing of the V4–V5 region of the samples of S. hansoni (brown), F. panamensis (dark morphotype), and Frankliniella (pale morphotype). Data are presented in percentages (%) for all sequences classified at this level of taxonomy. Relative abundance of bacterial OTUs at the genus level found for (c) samples obtained from the sequencing of the S. hansoni V1–V2 region. (d) Sequencing of the V4–V5 region of the samples of S. hansoni (brown), F. panamensis (dark morphotype), and Frankliniella (pale morphotype). Data are presented in percentages (%) for all sequences classified at this level of the taxonomy. (e) Venn diagram shows the “core” at the genus level. Relative abundance of genera based on the data obtained using the metataxonomic approach of the Antioquia thrips samples. Venn diagram showing overlays of OTUs, where 3 OTUs were common for the three thrips morphotypes (dark, pale, and brown) evaluated. Nomenclature: LV, LV2 and LVn correspond to the La Vega farm, VFm and VFp to the Villa Fátima farm, AASS1 and AASS2 to San Sebastian, BH1 to Barahonda, LS1 to Luisines, and EP1 to the farm El Paraíso.
Diversity 14 00540 g005
Figure 6. Estimated alpha diversity indices for (a) S. hansoni (brown morphotype). The diversity of panel (a) was estimated with sequence data from the V1–V2 region; (b) S. hansoni (brown morphotype), F. panamensis (dark morphotype), and Frankliniella (pale morphotype). Panel (b) diversity was estimated with sequence data from the V4–V5 region. (c) Structure of the bacterial community of thrips. (a) Principal Coordinate Analysis (PCoA) with the 16S rRNA gene. Nomenclature: LV, LV2 and LVn correspond to the La Vega farm, VFm and VFp correspond to the Villa Fátima farm, AASS1 and AASS2 to San Sebastian, BH1 to Barahonda, LS1 to Luisines, and EP1 to farm El Paraíso.
Figure 6. Estimated alpha diversity indices for (a) S. hansoni (brown morphotype). The diversity of panel (a) was estimated with sequence data from the V1–V2 region; (b) S. hansoni (brown morphotype), F. panamensis (dark morphotype), and Frankliniella (pale morphotype). Panel (b) diversity was estimated with sequence data from the V4–V5 region. (c) Structure of the bacterial community of thrips. (a) Principal Coordinate Analysis (PCoA) with the 16S rRNA gene. Nomenclature: LV, LV2 and LVn correspond to the La Vega farm, VFm and VFp correspond to the Villa Fátima farm, AASS1 and AASS2 to San Sebastian, BH1 to Barahonda, LS1 to Luisines, and EP1 to farm El Paraíso.
Diversity 14 00540 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cano-Calle, D.; Montoya-Porras, L.M.; Ochoa-Giraldo, S.; Junca, H.; Garcia-Bonilla, E.; Saldamando-Benjumea, C.; Moreno-Herrera, C.X.; Arango-Isaza, R.E. Thrips Microbiome Study in Commercial Avocado (Persea americana Mill.) from Northwest Colombian Andes (Antioquia, Colombia) Shows the Presence of Wolbachia, Ehrlichia, Enterobacter. Diversity 2022, 14, 540. https://doi.org/10.3390/d14070540

AMA Style

Cano-Calle D, Montoya-Porras LM, Ochoa-Giraldo S, Junca H, Garcia-Bonilla E, Saldamando-Benjumea C, Moreno-Herrera CX, Arango-Isaza RE. Thrips Microbiome Study in Commercial Avocado (Persea americana Mill.) from Northwest Colombian Andes (Antioquia, Colombia) Shows the Presence of Wolbachia, Ehrlichia, Enterobacter. Diversity. 2022; 14(7):540. https://doi.org/10.3390/d14070540

Chicago/Turabian Style

Cano-Calle, Daniela, Luisa Maria Montoya-Porras, Sebastian Ochoa-Giraldo, Howard Junca, Erika Garcia-Bonilla, Clara Saldamando-Benjumea, Claudia Ximena Moreno-Herrera, and Rafael E. Arango-Isaza. 2022. "Thrips Microbiome Study in Commercial Avocado (Persea americana Mill.) from Northwest Colombian Andes (Antioquia, Colombia) Shows the Presence of Wolbachia, Ehrlichia, Enterobacter" Diversity 14, no. 7: 540. https://doi.org/10.3390/d14070540

APA Style

Cano-Calle, D., Montoya-Porras, L. M., Ochoa-Giraldo, S., Junca, H., Garcia-Bonilla, E., Saldamando-Benjumea, C., Moreno-Herrera, C. X., & Arango-Isaza, R. E. (2022). Thrips Microbiome Study in Commercial Avocado (Persea americana Mill.) from Northwest Colombian Andes (Antioquia, Colombia) Shows the Presence of Wolbachia, Ehrlichia, Enterobacter. Diversity, 14(7), 540. https://doi.org/10.3390/d14070540

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

Article Metrics

Back to TopTop