Background
Laboratory adaptation process used in sterile insect technique (SIT) programs can exert a significant impact on the insect-gut microbiome relationship, which may negatively impact the quality and performance of the fly. In the present study, changes in the gut microbiota that occur through laboratory adaptation of two Ceratitis capitata populations were investigated: Vienna 8 genetic sexing strain (GSS), a long-established control line, and a wild population recently introduced to laboratory conditions. The bacterial profiles were studied for both strains using amplicon sequencing of the 16S rRNA V3-V4 hypervariable region in larvae and in the gastrointestinal tract of teneral (1 day) and adults (5 and 15 days) reared under laboratory conditions for 14 generations (F0–F13). Findings demonstrated the development of distinct bacterial communities across the generations with differences in the bacterial composition, suggesting a strong impact of laboratory adaptation on the fly bacteriome. Moreover, different bacterial profiles were observed between wild and Vienna 8 FD-GSS displaying different patterns between the developmental stages. Proteobacteria, mainly members of the Enterobacteriaceae family, represented the major component of the bacterial community followed by Firmicutes (mainly in Vienna 8 FD-GSS adults) and Chlamydiae. The distribution of these communities is dynamic across the generations and seems to be strain- and age-specific. In the Vienna 8 FD-GSS population, Providencia exhibited high relative abundance in the first three generations and decreased significantly later, while Klebsiella was relatively stable. In the wild population, Klebsiella was dominant across most of the generations, indicating that the wild population was more resistant to artificial rearing conditions compared with the Vienna 8 FD-GSS colony. Analysis of the core bacteriome revealed the presence of nine shared taxa between most of the examined medfly samples including Klebsiella, Providencia, Pantoea, and Pseudomonas. In addition, the operational taxonomic unit co-occurrence and mutual exclusion networks of the wild population indicated that most of the interactions were classified as co-presence, while in the Vienna 8 FD-GSS population, the number of mutual exclusions and co-presence interactions was equally distributed. Obtained results provided a thorough study of the dynamics of gut-associated bacteria during the laboratory adaptation of different Ceratitis capitata populations, serving as guidance for the design of colonization protocols, improving the effectiveness of artificial rearing and the SIT application.