Differential Effects of Antibiotic Therapy on the Structure and Function of Human Gut Microbiota

The human intestinal microbiota performs many essential functions for the host. Antimicrobial agents, such as antibiotics (AB), are also known to disturb microbial community equilibrium, thereby having an impact on human physiology. While an increasing number of studies investigate the effects of AB usage on changes in human gut microbiota biodiversity, its functional effects are still poorly understood. We performed a follow-up study to explore the effect of ABs with different modes of action on human gut microbiota composition and function. Four individuals were treated with different antibiotics and samples were taken before, during and after the AB course for all of them. Changes in the total and in the active (growing) microbiota as well as the functional changes were addressed by 16S rRNA gene and metagenomic 454-based pyrosequencing approaches. We have found that the class of antibiotic, particularly its antimicrobial effect and mode of action, played an important role in modulating the gut microbiota composition and function. Furthermore, analysis of the resistome suggested that oscillatory dynamics are not only due to antibiotic-target resistance, but also to fluctuations in the surviving bacterial community. Our results indicated that the effect of AB on the human gut microbiota relates to the interaction of several factors, principally the properties of the antimicrobial agent, and the structure, functions and resistance genes of the microbial community.     

Microarray analysis and barcoded pyrosequencing provide consistent microbial profiles depending on the source of human intestinal samples

Large-scale and in-depth characterization of the intestinal microbiota necessitates application of high-throughput 16S rRNA gene-based technologies, such as barcoded pyrosequencing and phylogenetic microarray analysis. In this study, the two techniques were compared and contrasted for analysis of the bacterial composition in three fecal and three small intestinal samples from human individuals. As PCR remains a crucial step in sample preparation for both techniques, different forward primers were used for amplification to assess their impact on microbial profiling results. An average of 7,944 pyrosequences, spanning the V1 and V2 region of 16S rRNA genes, was obtained per sample. Although primer choice in barcoded pyrosequencing did not affect species richness and diversity estimates, detection of Actinobacteria strongly depended on the selected primer. Microbial profiles obtained by pyrosequencing and phylogenetic microarray analysis (HITChip) correlated strongly for fecal and ileal lumen samples but were less concordant for ileostomy effluent. Quantitative PCR was employed to investigate the deviations in profiling between pyrosequencing and HITChip analysis. Since cloning and sequencing of random 16S rRNA genes from ileostomy effluent confirmed the presence of novel intestinal phylotypes detected by pyrosequencing, especially those belonging to the Veillonella group, the divergence between pyrosequencing and the HITChip is likely due to the relatively low number of available 16S rRNA gene sequences of small intestinal origin in the DNA databases that were used for HITChip probe design. Overall, this study demonstrated that equivalent biological conclusions are obtained by high-throughput profiling of microbial communities, independent of technology or primer choice.     

Characterization of the fecal bacteria communities of forage-fed horses by pyrosequencing of 16S rRNA V4 gene amplicons

The diversity of the equine fecal bacterial community was evaluated using pyrosequencing of 16S rRNA gene amplicons. Fecal samples were obtained from horses fed cool-season grass hay. Fecal bacteria were characterized by amplifying the V4 region of bacterial 16S rRNA gene. Of 5898 mean unique sequences, a mean of 1510 operational taxonomic units were identified in the four fecal samples. Equine fecal bacterial richness was higher than that reported in humans, but lower than that reported in either cattle feces or soil. Bacterial classified sequences were assigned to 16 phyla, of which 10 were present in all samples. The largest number of reads belonged to Firmicutes (43.7% of total bacterial sequences), Verrucomicrobia (4.1%), Proteobacteria (3.8%), and Bacteroidetes (3.7%). The less abundant Actinobacteria, Cyanobacteria, and TM7 phyla presented here have not been previously described in the gut contents or feces of horses. Unclassified sequences represented 38.1% of total bacterial sequences; therefore, the equine fecal microbiome diversity is likely greater than that described. This is the first study to characterize the fecal bacterial community in horses by the use of 16S rRNA gene amplicon pyrosequencing, expanding our knowledge of the fecal microbiota of forage-fed horses.     

Characterization of the Fecal Microbiota Using High-Throughput Sequencing Reveals a Stable Microbial Community during Storage

The handling and treatment of biological samples is critical when characterizing the composition of the intestinal microbiota between different ecological niches or diseases. Specifically, exposure of fecal samples to room temperature or long term storage in deep freezing conditions may alter the composition of the microbiota. Thus, we stored fecal samples at room temperature and monitored the stability of the microbiota over twenty four hours. We also investigated the stability of the microbiota in fecal samples during a six month storage period at −80°C. As the stability of the fecal microbiota may be affected by intestinal disease, we analyzed two healthy controls and two patients with irritable bowel syndrome (IBS). We used high-throughput pyrosequencing of the 16S rRNA gene to characterize the microbiota in fecal samples stored at room temperature or −80°C at six and seven time points, respectively. The composition of microbial communities in IBS patients and healthy controls were determined and compared using the Quantitative Insights Into Microbial Ecology (QIIME) pipeline. The composition of the microbiota in fecal samples stored for different lengths of time at room temperature or −80°C clustered strongly based on the host each sample originated from. Our data demonstrates that fecal samples exposed to room or deep freezing temperatures for up to twenty four hours and six months, respectively, exhibit a microbial composition and diversity that shares more identity with its host of origin than any other sample.     

Bacterial Community Mapping of the Mouse Gastrointestinal Tract

Keeping mammalian gastrointestinal (GI) tract communities in balance is crucial for host health maintenance. However, our understanding of microbial communities in the GI tract is still very limited. In this study, samples taken from the GI tracts of C57BL/6 mice were subjected to 16S rRNA gene sequence-based analysis to examine the characteristic bacterial communities along the mouse GI tract, including those present in the stomach, duodenum, jejunum, ileum, cecum, colon and feces. Further analyses of the 283,234 valid sequences obtained from pyrosequencing revealed that the gastric, duodenal, large intestinal and fecal samples had higher phylogenetic diversity than the jejunum and ileum samples did. The microbial communities found in the small intestine and stomach were different from those seen in the large intestine and fecal samples. A greater proportion of Lactobacillaceae were found in the stomach and small intestine, while a larger proportion of anaerobes such as Bacteroidaceae, Prevotellaceae, Rikenellaceae, Lachnospiraceae, and Ruminococcaceae were found in the large intestine and feces. In addition, inter-mouse variations of microbiota were observed between the large intestinal and fecal samples, which were much smaller than those between the gastric and small intestinal samples. As far as we can ascertain, ours is the first study to systematically characterize bacterial communities from the GI tracts of C57BL/6 mice.     

Long-Term Temporal Analysis of the Human Fecal Microbiota Revealed a Stable Core of Dominant Bacterial Species

Next-generation sequencing has greatly contributed to an improved ecological understanding of the human gut microbiota. Nevertheless, questions remain regarding the characteristics of this ecosystem and the ecological processes that shape it, and controversy has arisen regarding the stability of the bacterial populations and the existence of a temporal core. In this study, we have characterized the fecal microbial communities of three human individuals over a one-year period by 454 pyrosequencing of 16S rRNA tags in order to investigate the temporal characteristics of the bacterial communities. The findings revealed a temporal core of 33 to 40 species-level Operational Taxonomic Units (OTUs) within subjects. Although these OTUs accounted only for around 12% of the total OTUs detected, they added up to >75% of the total sequences obtained for each individual. In order to determine the capacity of the sequencing and bioinformatic approaches applied during this study to accurately determine the proportion of a core microbiota, we analyzed the fecal microbiota of nine mice with a defined three-member community. This experiment revealed that the sequencing approach inflated the amount of rare OTUs, which introduced a significant degree of artificial variation across samples, and hence reduced the apparent fraction of shared OTUs. However, when assessing the data quantitatively by focusing on dominant lineages, the sequencing approaches deliver an accurate representation of the community. In conclusion, this study revealed that the human fecal microbiota is dominated by around 40 species that maintain persistent populations over the duration of one year. The findings allow conclusions about the ecological factors that shape the community and support the concept of a homeostatic ecosystem controlled largely by deterministic processes. Our analysis of a three-member community revealed that methodological artifacts of OTU-based approaches complicate core calculations, and these limitations have to be considered in the interpretation of microbiome studies.     

Changes in fecal microbiota of healthy dogs administered amoxicillin

The effect of oral amoxicillin treatment on fecal microbiota of seven healthy adult dogs was determined with a focus on the prevalence of bacterial antibiotic resistance and changes in predominant bacterial populations. After 4–7 days of exposure to amoxicillin, fecal Escherichia coli expressed resistance to multiple antibiotics when compared with the pre-exposure situation. Two weeks postexposure, the susceptibility pattern had returned to pre-exposure levels in most dogs. A shift in bacterial populations was confirmed by molecular fingerprinting of fecal bacterial populations using denaturing gradient gel electrophoresis (PCR-DGGE) of the 16S V3 rRNA gene region. Much of the variation in DGGE profiles could be attributed to dog-specific factors. However, permutation tests indicated that amoxicillin exposure significantly affected the DGGE profiles after controlling for the dog effect ( P =0.02), and pre-exposure samples were clearly separated from postexposure samples. Sequence analysis of DGGE bands and real-time PCR quantification indicated that amoxicillin exposure caused a shift in the intestinal ecological balance toward a Gram-negative microbiota including resistant species in the family Enterobacteriaceae .     

Microbial Diversity and Evidence of Novel Homoacetogens in the Gut of Both Geriatric and Adult Giant Pandas (Ailuropoda melanoleuca)

Recent studies have described the bacterial community residing in the guts of giant pandas, together with the presence of lignocellulolytic enzymes. However, a more comprehensive understanding of the intestinal microbial composition and its functional capacity in giant pandas remains a major goal. Here, we conducted a comparison of bacterial, fungal and homoacetogenic microbial communities from fecal samples taken from two geriatric and two adult captive giant pandas. 16S rDNA amplicon pyrosequencing revealed that Firmicutes and Proteobacteria are the most abundant microbiota in both geriatric and adult giant pandas. However, members of phylum Actinobacteria found in adult giant pandas were absent in their geriatric counterparts. Similarly, ITS1 amplicon pyrosequencing identified developmental changes in the most abundant fungal classes from Sordariomycetes in adult pandas to Saccharomycetes in geriatric pandas. Geriatric pandas exhibited significantly higher abundance of a potential probiotic fungus (Candida tropicalis) as compared to adult pandas, indicating their importance in the normal digestive physiology of aged pandas. Our study also reported the presence of a lignocellulolytic white-rot fungus, Perenniporia medulla-panis, and the evidence of novel homoacetogens residing in the guts of giant pandas.     

The Pervasive Effects of an Antibiotic on the Human Gut Microbiota,as Revealed by Deep 16S rRNA Sequencing

The human intestinal microbiota is essential to the health of the host and plays a role in nutrition, development, metabolism, pathogen resistance, and regulation of immune responses. Antibiotics may disrupt these coevolved interactions, leading to acute or chronic disease in some individuals. Our understanding of antibiotic-associated disturbance of the microbiota has been limited by the poor sensitivity, inadequate resolution, and significant cost of current research methods. The use of pyrosequencing technology to generate large numbers of 16S rDNA sequence tags circumvents these limitations and has been shown to reveal previously unexplored aspects of the “rare biosphere.” We investigated the distal gut bacterial communities of three healthy humans before and after treatment with ciprofloxacin, obtaining more than 7,000 full-length rRNA sequences and over 900,000 pyrosequencing reads from two hypervariable regions of the rRNA gene. A companion paper in PLoS Genetics (see Huse et al., doi: 10.1371/journal.pgen.1000255) shows that the taxonomic information obtained with these methods is concordant. Pyrosequencing of the V6 and V3 variable regions identified 3,300–5,700 taxa that collectively accounted for over 99% of the variable region sequence tags that could be obtained from these samples. Ciprofloxacin treatment influenced the abundance of about a third of the bacterial taxa in the gut, decreasing the taxonomic richness, diversity, and evenness of the community. However, the magnitude of this effect varied among individuals, and some taxa showed interindividual variation in the response to ciprofloxacin. While differences of community composition between individuals were the largest source of variability between samples, we found that two unrelated individuals shared a surprising degree of community similarity. In all three individuals, the taxonomic composition of the community closely resembled its pretreatment state by 4 weeks after the end of treatment, but several taxa failed to recover within 6 months. These pervasive effects of ciprofloxacin on community composition contrast with the reports by participants of normal intestinal function and with prior assumptions of only modest effects of ciprofloxacin on the intestinal microbiota. These observations support the hypothesis of functional redundancy in the human gut microbiota. The rapid return to the pretreatment community composition is indicative of factors promoting community resilience, the nature of which deserves future investigation.     

Bacteria,phages and pigs: the effects of in-feed antibiotics on the microbiome at different gut locations

Disturbance of the beneficial gut microbial community is a potential collateral effect of antibiotics, which have many uses in animal agriculture (disease treatment or prevention and feed efficiency improvement). Understanding antibiotic effects on bacterial communities at different intestinal locations is essential to realize the full benefits and consequences of in-feed antibiotics. In this study, we defined the lumenal and mucosal bacterial communities from the small intestine (ileum) and large intestine (cecum and colon) plus feces, and characterized the effects of in-feed antibiotics (chlortetracycline, sulfamethazine and penicillin (ASP250)) on these communities. 16S rRNA gene sequence and metagenomic analyses of bacterial membership and functions revealed dramatic differences between small and large intestinal locations, including enrichment of Firmicutes and phage-encoding genes in the ileum. The large intestinal microbiota encoded numerous genes to degrade plant cell wall components, and these genes were lacking in the ileum. The mucosa-associated ileal microbiota harbored greater bacterial diversity than the lumen but similar membership to the mucosa of the large intestine, suggesting that most gut microbes can associate with the mucosa and might serve as an inoculum for the lumen. The collateral effects on the microbiota of antibiotic-fed animals caused divergence from that of control animals, with notable changes being increases in Escherichia coli populations in the ileum, Lachnobacterium spp. in all gut locations, and resistance genes to antibiotics not administered. Characterizing the differential metabolic capacities and response to perturbation at distinct intestinal locations will inform strategies to improve gut health and food safety.     

Variation in Taxonomic Composition of the Fecal Microbiota in an Inbred Mouse Strain across Individuals and Time

Genetics, diet, and other environmental exposures are thought to be major factors in the development and composition of the intestinal microbiota of animals. However, the relative contributions of these factors in adult animals, as well as variation with time in a variety of important settings, are still not fully understood. We studied a population of inbred, female mice fed the same diet and housed under the same conditions. We collected fecal samples from 46 individual mice over two weeks, sampling four of these mice for periods as long as 236 days for a total of 190 samples, and determined the phylogenetic composition of their microbial communities after analyzing 1,849,990 high-quality pyrosequencing reads of the 16S rRNA gene V3 region. Even under these controlled conditions, we found significant inter-individual variation in community composition, as well as variation within an individual over time, including increases in alpha diversity during the first 2 months of co-habitation. Some variation was explained by mouse membership in different cage and vendor shipment groups. The differences among individual mice from the same shipment group and cage were still significant. Overall, we found that 23% of the variation in intestinal microbiota composition was explained by changes within the fecal microbiota of a mouse over time, 12% was explained by persistent differences among individual mice, 14% by cage, and 18% by shipment group. Our findings suggest that the microbiota of controlled populations of inbred laboratory animals may not be as uniform as previously thought, that animal rearing and handling may account for some variation, and that as yet unidentified factors may explain additional components of variation in the composition of the microbiota within populations and individuals over time. These findings have implications for the design and interpretation of experiments involving laboratory animals.     

Microbiome Analysis of Stool Samples from African Americans with Colon Polyps

Background

Colonic polyps are common tumors occurring in ~50% of Western populations with ~10% risk of malignant progression. Dietary agents have been considered the primary environmental exposure to promote colorectal cancer (CRC) development. However, the colonic mucosa is permanently in contact with the microbiota and its metabolic products including toxins that also have the potential to trigger oncogenic transformation.

Aim

To analyze fecal DNA for microbiota composition and functional potential in African Americans with pre-neoplastic lesions.

Materials & Methods

We analyzed the bacterial composition of stool samples from 6 healthy individuals and 6 patients with colon polyps using 16S ribosomal RNA-based phylogenetic microarray; the Human intestinal Tract Chip (HITChip) and 16S rRNA gene barcoded 454 pyrosequencing. The functional potential was determined by sequence-based metagenomics using 454 pyrosequencing.

Results

Fecal microbiota profiling of samples from the healthy and polyp patients using both a phylogenetic microarraying (HITChip) and barcoded 454 pyrosequencing generated similar results. A distinction between both sets of samples was only obtained when the analysis was performed at the sub-genus level. Most of the species leading to the dissociation were from the Bacteroides group. The metagenomic analysis did not reveal major differences in bacterial gene prevalence/abundances between the two groups even when the analysis and comparisons were restricted to available Bacteroides genomes.

Conclusion

This study reveals that at the pre-neoplastic stages, there is a trend showing microbiota changes between healthy and colon polyp patients at the sub-genus level. These differences were not reflected at the genome/functions levels. Bacteria and associated functions within the Bacteroides group need to be further analyzed and dissected to pinpoint potential actors in the early colon oncogenic transformation in a large sample size.     

Intestinal bacterial communities that produce active estrogen-like compounds enterodiol and enterolactone in humans

Lignans are dietary diphenolic compounds which require activation by intestinal bacteria to exert possible beneficial health effects. The intestinal ecosystem plays a crucial role in lignan metabolism, but the organisms involved are poorly described. To characterize the bacterial communities responsible for secoisolariciresinol (SECO) activation, i.e., the communities that produce the enterolignans enterodiol (ED) and enterolactone (EL), a study with 24 human subjects was undertaken. SECO activation was detected in all tested fecal samples. The intestinal bacteria involved in ED production were part of the dominant microbiota (6 x 10(8) CFU g(-1)), as revealed by most-probable-number enumerations. Conversely, organisms that catalyzed the formation of EL occurred at a mean concentration of approximately 3 x 10(5) CFU g(-1). Women tended to have higher concentrations of both ED- and EL-producing organisms than men. Significantly larger amounts of EL were produced by fecal dilutions from individuals with moderate to high concentrations of EL-producing bacteria. Two organisms able to demethylate and dehydroxylate SECO were isolated from human feces. Based on 16S rRNA gene sequence analyses, they were named Peptostreptococcus productus SECO-Mt75m3 and Eggerthella lenta SECO-Mt75m2. A new 16S rRNA-targeted oligonucleotide probe specific for P. productus and related species was designed and further used in fluorescent in situ hybridization experiments, along with five additional group-specific probes. Significantly higher proportions of P. productus and related species (P = 0.012), as well as bacteria belonging to the Atopobium group (P = 0.035), were typical of individuals with moderate to high concentrations of EL-producing communities.     

Terminal restriction fragment length polymorphism analysis for human fecal microbiota and its application for analysis of complex bifidobacterial communities

Terminal restriction fragment length polymorphism (T-RFLP) analysis was used to characterize and compare human fecal microbiota among individuals. T-RFLP patterns of fecal 16S ribosomal DNA (rDNA) PCR products from three adults revealed host-specific bacterial communities and were in good agreement with those reported in our previous study. In addition, we applied T-RFLP analysis for the analysis of complex bifidobacterial communities in human fecal samples. The developed method based on Bifidobacterium genus-specific PCR and T-RFLP could identify more than one bifidobacterial species. T-RFLP patterns of Bifidobacterium genus-specific PCR products from the fecal samples were host-specific as well as those of fecal 16S rDNA PCR products. These results were confirmed by PCR-denaturing gradient gel electrophoresis (DGGE) with primers specific for the genus Bifidobacterium and Bifidobacterium species- and group-specific PCR. Our study demonstrates that T-RFLP analysis is useful for assessment of the diversity of the human fecal microbiota and rapid comparison of the community structure among individuals, and that the applied method is useful for rapid and sensitive analysis of bifidobacterial community.     

The effect of the macrolide antibiotic tylosin on microbial diversity in the canine small intestine as demonstrated by massive parallel 16S rRNA gene sequencing

Background  

Recent studies have shown that the fecal microbiota is generally resilient to short-term antibiotic administration, but some bacterial taxa may remain depressed for several months. Limited information is available about the effect of antimicrobials on small intestinal microbiota, an important contributor to gastrointestinal health. The antibiotic tylosin is often successfully used for the treatment of chronic diarrhea in dogs, but its exact mode of action and its effect on the intestinal microbiota remain unknown. The aim of this study was to evaluate the effect of tylosin on canine jejunal microbiota. Tylosin was administered at 20 to 22 mg/kg q 24 hr for 14 days to five healthy dogs, each with a pre-existing jejunal fistula. Jejunal brush samples were collected through the fistula on days 0, 14, and 28 (14 days after withdrawal of tylosin). Bacterial diversity was characterized using massive parallel 16S rRNA gene pyrosequencing.     

Antibiotics conspicuously affect community profiles and richness, but not the density of bacterial cells associated with mucosa in the large and small intestines of mice

The influence of three antibiotics (bacitracin, enrofloxacin, and neomycin sulfate) on the mucosa-associated enteric microbiota and the intestines of mice was examined. Antibiotics caused conspicuous enlargement of ceca and an increase in overall length of the intestine. However, there were no pathologic changes associated with increased cecal size or length of the intestine. Conspicuous reductions in the richness of mucosa-associated bacteria and changes to community profiles within the small (duodenum, proximal jejunum, middle jejunum, distal jejunum, and ileum) and large (cecum, ascending colon, and descending colon) intestine occurred in mice administered antibiotics. Communities in antibiotic-treated mice were dominated by a limited number of Clostridium-like (i.e. clostridial cluster XIVa) and Bacteroides species. The richness of mucosa-associated communities within the small and large intestine increased during the 14-day recovery period. However, community profiles within the large intestine did not return to baseline (i.e. relative to the control). Although antibiotic administration greatly reduced bacterial richness, densities of mucosa-associated bacteria were not reduced correspondingly. These data showed that the antibiotics, bacitracin, enrofloxacin, and neomycin sulfate, administered for 21 days to mice did not sterilize the intestine, but did impart a tremendous and prolonged impact on mucosa-associated bacterial communities throughout the small and large intestine.     

Effects of feed additives on the development on the ileal bacterial community of the broiler chicken

Intensifying concerns about the use of antimicrobials in meat and poultry production has enhanced interest in the application of prebiotics, probiotics and enzymes to enhance growth and prevent disease in food animals. Growth-promoting antibiotics enhance growth of animals by reducing the load of bacteria in the intestine, by reducing colonization by intestinal pathogens or by enhancing the growth and/or metabolism of beneficial bacteria in the intestine. Recently, molecular ecology, utilizing DNA-sequence heterogeneity of the 16S rRNA gene, has revealed a surprising diversity of uncharacterized bacteria inhabiting this ecosystem. We used this approach to determine the effect of growth-promoting antibiotics on the development and composition of the ileal bacterial community. Pairwise comparisons, correspondence analysis and community diversity indices revealed significant differences among the treatments (bacitracin/virginiamycin or monensin) and controls. Antibiotics reduced the diversity of the ileal bacterial community and induced communities rich in Clostridia throughout the life of the broiler chicken. These results indicate that some bacterial species, such as lactobacilli, were suppressed and also suggest that many intestinal Clostridia may be non-pathogenic. Future studies should focus on characterizing the important bacterial species needed to stabilize the intestinal microbiota and identifying those commensals that stimulate and enhance development of intestinal function.     

Bacteriophage versus antibiotic therapy on gut bacterial communities of juvenile green turtle,Chelonia mydas

Green turtles are endangered marine herbivorous hindgut fermenters that contribute to a variety of marine ecosystems. Debilitated turtles are often rehabilitated in turtle hospitals. Since accurate diagnosis of disease is difficult, broad-spectrum antibiotics are routinely used as a general treatment, potentially causing collateral damage to the gut microbiome of the patient. Here, we evaluated the concept of the application of bacteriophage (phages) to eliminate targeted intestinal bacteria as an alternative to a broad-spectrum antibiotic (enrofloxacin) in clinically healthy, captive green turtles. Additionally, the impact of a broad-spectrum antibiotic (enrofloxacin) and phage therapy on the gut bacterial communities of green turtles was evaluated. Gut bacterial communities in faecal samples were analysed by sequencing the V1-V3 regions of the bacterial 16S rRNA. Bacteria-specific phage cocktails significantly (P < 0.05) reduced targeted Acinetobacter in phage-treated turtles during the therapy. Compared to control, no significant difference was observed in the bacterial diversity and compositions in phage-treated turtles. In contrast, bacterial diversity was significantly (P < 0.05) reduced in antibiotic-treated turtles at day 15 and throughout the trial. The alteration in the bacterial microbiota of antibiotic-treated turtles was largely due to an increase in abundance of Gram-positive Firmicutes and a concurrent decrease in Gram-negative Bacteroidetes, Proteobacteria and Verrucomicrobia. Additionally, we observed the relative abundance of several bacteria at lower taxonomic level was much less affected by phages than by antibiotics. These data offer the proof of concept of phage therapy to manipulate transient as well as indigenous bacterial flora in gut-related dysbiosis of turtles.