A comparison of nonlethal sampling methods for amphibian gut microbiome analyses

Contrasting recovery of metagenome‑assembled genomes and derived bacterial communities and functional profiles from lizard fecal and cloacal samples

Genome-resolved metagenomics, based on shotgun sequencing, has become a powerful strategy for investigating animal-associated bacterial communities, due its heightened capability for delivering detailed taxonomic, phylogenetic, and functional insights compared to amplicon sequencing-based approaches. While genome-resolved metagenomics holds promise across various non-lethal sample types, their effectiveness in yielding high-quality metagenome-assembled genomes remains largely unexplored. Our investigation of fecal and cloacal microbiota of the mesquite lizards (Sceloporus grammicus) using genome-resolved metagenomics revealed that fecal samples contributed 97% of the 127 reconstructed bacterial genomes, whereas only 3% were recovered from cloacal swabs, which were largely enriched with host DNA. Taxonomic, phylogenetic and functional alpha bacterial diversity was greater in fecal samples than in cloacal swabs. We also observed significant differences in bacterial community composition between sampling methods, and higher inter-individual variation in cloacal swabs. Bacteroides, Phocaeicola and Parabacteroides (all Bacteroidota) were more abundant in the feces, whereas Hafnia and Salmonella (both Pseudomonadota) increased in the cloaca. Functional analyses showed that metabolic capacities of the microbiota to degrade polysaccharides, sugars and nitrogen compounds were enriched in fecal samples, likely reflecting the role of intestinal bacteria in nutrient metabolism. Overall, our results indicate that fecal samples outperform cloacal swabs in characterizing bacterial assemblages within lizards using genome-resolved metagenomics.

Short-term dynamics of fecal microbiome and antibiotic resistance in juvenile rainbow trout (Oncorhynchus mykiss) following antibiotic treatment and withdrawal

BACKGROUND

In aquaculture, the secretions of cultured organisms contribute to the development of aquatic antibiotic resistance. However, the antibiotic-induced changes in fish feces remain poorly understood. This study aimed to assess the short-term dynamics of fecal microbiome and antibiotic resistance in juvenile rainbow trout (Oncorhynchus mykiss) upon antibiotic treatment and withdrawal period.

METHODS

Fish were orally administered diets supplemented with oxytetracycline (OTC) or sulfadiazine/trimethoprim (SDZ/TMP) for 10 consecutive days, followed by a 25-day withdrawal period. Fecal samples were collected before antibiotic treatment (day 0), and at 1, 3, 7, and 10 days post antibiotic administration (dpa), as well as 1, 3, 7, 14, and 25 days post antibiotic cessation (dpc). The fecal microbiome community was profiled using both culture-dependent and -independent methods. The relative abundance of antibiotic resistance genes (ARGs) and the class 1 integron-integrase gene (intI1) in the feces were quantified using real-time PCR.

RESULTS

Antibiotic treatment disrupted the fecal microbial communities, and this alteration persisted even after antibiotic cessation. Moreover, OTC treatment increased the relative abundance of tet genes, while sul and dfr genes increased in the SDZ/TMP-treated group. Notably, Flavobacterium, Pseudomonas, and Streptococcus exhibited a significant correlation with the abundance of ARGs, suggesting their potential role as carriers for ARGs.

CONCLUSION

This study demonstrates the antibiotic-induced changes in the fecal microbiome and the increase of ARGs in rainbow trout feces. These findings provide novel insights into the dynamics of microbiome recovery post-antibiotic cessation and suggest that fish feces provide a non-invasive approach to predict changes in the fish gut microbiome and resistome.

Insights into the Gut Microbiome of the South American Leaf-Toed Gecko (Phylodactylus gerropygus) Inhabiting the Core of the Atacama Desert

Living in arid environments presents unique challenges to organisms, including limited food and water, extreme temperatures, and UV exposure. Reptiles, such as the South American leaf-toed gecko (Phyllodactylus gerrhopygus), have evolved remarkable adaptations to thrive in such harsh conditions. The gut microbiome plays a critical role in host adaptation and health, yet its composition remains poorly characterized in desert reptiles. This study aimed to characterize the composition and abundance of the gut microbiome in P. gerrhopygus inhabiting the hyperarid Atacama Desert, taking into account potential sex differences. Fecal samples from adult female and male geckos were analyzed by 16S rRNA gene amplicon sequencing. No significant differences in bacterial alpha diversity were observed between the sexes. However, the phylum Bacteroidota was more abundant in females, while males had a higher Firmicutes/Bacteroidota ratio. The core microbiome was dominated by the phyla Bacteroidota, Firmicutes, and Proteobacteria in both sexes. Analysis of bacterial composition revealed 481 amplicon sequence variants (ASVs) shared by female and male geckos. In addition, 108 unique ASVs were exclusive to females, while 244 ASVs were unique to males. Although the overall bacterial composition did not differ significantly between the sexes, certain taxa exhibited higher relative abundances in each sex group. This study provides insight into the taxonomic structure of the gut microbiome in a desert-adapted reptile and highlights potential sex-specific differences. Understanding these microbial communities is critical for elucidating the mechanisms underlying host resilience in Earth's most arid environments, and for informing conservation efforts in the face of ongoing climate change.

Host-associated helminth diversity and microbiome composition contribute to anti-pathogen defences in tropical frogs impacted by forest fragmentation

Habitat fragmentation can negatively impact wildlife populations by simplification of ecological interactions, but little is known about how these impacts extend to host-associated symbiotic communities. The symbiotic communities of amphibians play important roles in anti-pathogen defences, particularly against the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). In this study, we analyse the role of macroparasitic helminth communities in concert with microbial communities in defending the host against Bd infection within the context of forest fragmentation. We found that skin microbial and helminth communities are disrupted at fragmented habitats, while gut microbiomes appear more resilient to environmental change. We also detected potential protective roles of helminth diversity and anti-pathogen microbial function in limiting Bd infection. Microbial network analysis revealed strong patterns of structure in both skin and gut communities, with helminths playing central roles in these networks. We reveal consistent roles of microbial and helminth diversity in driving host-pathogen interactions and the potential implications of fragmentation on host fitness.

Microbial communities are thermally more sensitive in warm-climate lizards compared with their cold-climate counterparts

Environmental temperature affects the composition, structure, and function of the gut microbial communities in host animals. To elucidate the role of gut microbiota in thermal adaptation, we designed a 2 species × 3 temperatures experiment, whereby we acclimated adult males of two agamid lizard species (warm-climate Leiolepis reevesii and cold-climate Phrynocephalus przewalskii) to 20, 28, and 36°C for 2 weeks and then collected their fecal and small-intestinal samples to analyze and compare the microbiota using 16S rRNA gene amplicon sequencing technology. The fecal microbiota displayed more pronounced interspecific differences in microbial community than the small-intestinal microbiota in the two species occurring in thermally different regions. The response of fecal and small-intestinal microbiota to temperature increase or decrease differed between the two species, with more bacterial taxa affected by acclimation temperature in L. reevesii than in P. przewalskii. Both species, the warm-climate species in particular, could cope with temperature change by adjusting the relative abundance of functional categories associated with metabolism and environmental information processing. Functional genes associated with carbohydrate metabolism were enhanced in P. przewalskii, suggesting the contribution of the fecal microbiota to cold-climate adaptation in P. przewalskii. Taken together, our results validate the two hypotheses tested, of which one suggests that the gut microbiota should help lizards adapt to thermal environments in which they live, and the other suggests that microbial communities should be thermally more sensitive in warm-climate lizards than in cold-climate lizards.

Comparative study of the gut microbial communities collected by scraping and swabbing in a fish model: a comprehensive guide to promote non-lethal procedures for gut microbial studies

In the present study, we propose the use of swabs in non-lethal sampling procedures to collect the mucosa-adhered gut microbiota from the posterior intestine of fish, and therefore, we compare the bacterial communities collected by conventional scraping and by swabbing methods. For this purpose, samples of the posterior intestine of rainbow trout (Oncorhynchus mykiss) were collected first using the swabbing approach, and after fish euthanasia, by mucosa scraping. Finally, bacterial communities were compared by 16S rRNA gene Illumina sequencing. Results from the current study revealed that similar values of bacterial richness and diversity were found for both sampling procedures. Similarly, there were no differences between procedures when using qualitative metrics (Jaccard and unweighted UniFrac) for estimating inter-individual diversity, but the quantitative metrics (Bray-Curtis and weighted UniFrac) showed a higher dispersion when samples were obtained by swabbing compared to scraping. In terms of bacterial composition, there were differences in abundance for the phyla Firmicutes and Proteobacteria. The cause of these differential abundances may be the inability of the swab to access to certain areas, such as the basal region of the intestinal villi. Moreover, swabbing allowed a higher representation of low abundant taxa, which may also have an important role in host microbiome regardless of their low abundance. Overall, our results demonstrate that the sampling method is a factor to be considered in experimental design when studying gut bacterial communities to avoid potential biases in the interpretation or comparison of results from different studies. In addition, the advantages and disadvantages of each procedure (swabbing vs scraping) are discussed in detail, concluding that swabbing can be implemented as a reliable and non-lethal procedure for posterior gut microbiota studies, which is of particular interest for animal welfare and the 3Rs principle, and may offer a wide range of novel applications.

Comparative analysis of two nonlethal methods for the study of the gut bacterial communities in wild lizards

Fecal samples or cloacal swabs are preferred over lethal dissections to study vertebrate gut microbiota for ethical reasons, but it remains unclear which nonlethal methods provide more accurate information about gut microbiota. We compared the bacterial communities of three gastrointestinal tract (GIT) segments, that is, stomach, small intestine (midgut), and rectum (hindgut) with the bacterial communities of the cloaca and feces in the mesquite lizard Sceloporus grammicus. The hindgut had the highest taxonomic and functional alpha diversity, followed by midgut and feces, whereas the stomach and cloaca showed the lowest diversities. The taxonomic assemblages of the GIT segments at the phylum level were strongly correlated with those retrieved from feces and cloacal swabs (rs > 0.84 in all cases). The turnover ratio of Amplicon Sequence Variants (ASVs) between midgut and hindgut and the feces was lower than the ratio between these segments and the cloaca. More than half of the core-ASVs in the midgut (24 of 32) and hindgut (58 of 97) were also found in feces, while less than 5 were found in the cloaca. At the ASVs level, however, the structure of the bacterial communities of the midgut and hindgut were similar to those detected in feces and cloaca. Our findings suggest that fecal samples and cloacal swabs of spiny lizards provide a good approximation of the taxonomic assemblages and beta diversity of midgut and hindgut microbiota, while feces better represent the bacterial communities of the intestinal segments at a single nucleotide variation level than cloacal swabs.

Assessing the validity of fecal sampling for characterizing variation in threespine stickleback's gut microbiota

The gut microbiota is crucial for many aspects of their hosts' biology, and it has been characterized for many species across the animal kingdom. Yet, we still don't have a good understanding of whether non-lethal sampling can accurately capture the diversity of gut-associated bacterial communities, as estimated from lethal sampling of intestinal tissue. We further lack knowledge on whether non-lethal sampling methods are suitable for detecting gut microbiota shifts associated with changes in environmental factors (e.g., diet). We addressed these questions in threespine stickleback fish, a model system for evolutionary ecology, by comparing bacterial communities from intestinal tissue and feces. Despite some differences in community composition between the two sample types and considerable temporal variation among fecal samples, bacterial communities appear to largely overlap. Further, we detected consistent and significant changes of fecal bacterial communities associated with an experimental diet manipulation. This suggests that fecal sampling can represent an adequate non-lethal method to characterize the gut microbiota of threespine stickleback, but additional studies will be necessary before drawing general conclusions regarding the validity of fecal sampling for gut microbiota studies. To this end, we give recommendations to improve the characterization of the gut microbiota via fecal sampling. Fecal sampling allows studying temporal gut microbiota shifts associated with environmental change at the individual level, which increases opportunities for future experimental gut microbiota research.

The structural discrepancy between the small and large gut microbiota of Asiatic toad (Bufo gargarizans) during hibernation

Hibernating amphibians are suitable for the research on the adaptation of gut microbiota to long-term fasting and cold stresses. However, the previous studies mainly focus on the large or whole gut microbiota but not the small gut microbiota. To test the structural discrepancy between the small and large gut microbiota during hibernation, we performed two independent batches of 16S rRNA gene amplicon sequencing to profile the small and large gut microbiota of hibernating Asiatic toad (Bufo gargarizans) from two wild populations. Both batches of data revealed that Proteobacteria, Bacteroidetes, and Firmicutes were the three most dominant phyla in the small and large gut microbiota. Three core OTUs with 100% occurrence in all gut microbiotas were annotated as Pseudomonas. A significant structural discrepancy was detected between the small and large gut microbiota. For instance, Proteobacteria assembled in the small intestine with a higher proportion than it did in the large intestine, but Bacteroidetes and Firmicutes assembled in the large intestine with a higher proportion than they did in the small intestine. The large gut microbiota exhibited higher diversity than the small gut microbiota. Nevertheless, a severe batch effect existed in the structural analysis of the gut microbiotas. The large gut microbiota showed a better resistance to the batch effect than the small gut microbiota did. This study provides preliminary evidence that microbes assemble in the small and large intestines of amphibians with discrepant patterns during hibernation.

Screening Salamanders for Symbionts

Microbial symbionts are broadly categorized by their impacts on host fitness: commensals, pathogens, and mutualists. However, recent investigations into the physiological basis of these impacts have revealed nuanced microbial influences on a wide range of host developmental, immunological, and physiological processes, including regeneration. Exploring these impacts begins with knowing which microbes are present. This methodological pipeline contains both targeted assays using PCR and culturing, as well as culture-independent approaches, to survey host salamander tissues for common and unknown microbial symbionts.

Lessons from the diet: Captivity and sex shape the gut microbiota in an oviparous lizard (Calotes versicolor)

Studies have indicated that the abundance and community structure of gut microbiota are altered by diet. In this study, next-generation sequencing of the 16S rRNA gene amplicon was performed to evaluate variations in the gut microbiota of wild and captive individuals of both sexes of Calotes versicolor. The results showed that there was a significant sex difference in microbial community structure for wild C. versicolor, Bacteroide was the dominant genus in wild females (WF), whereas Ochrobactrum was the dominant genus in wild males (WM). Acinetobacter and Hymenobacter were the dominant genera in WF, while Clostridium was the dominant genus in captive females (CF). The results indicated that differences in diet between wild and captive C. versicolor also resulted in variations in gut microbiota. Thus, it was not surprising that captivity and sex shape the gut microbiota in C. versicolor. In summary, the fundamental information presented about the gut microbiota of both sexes of wild (and captive females) C. versicolor, indicates that the artificial environments are not suitable for the wild C. versicolor.

Diet diversity and environment determine the intestinal microbiome and bacterial pathogen load of fire salamanders

Diverse communities of symbiotic microbes inhabit the digestive systems of vertebrates and play a crucial role in animal health, and host diet plays a major role in shaping the composition and diversity of these communities. Here, we characterized diet and gut microbiome of fire salamander populations from three Belgian forests. We carried out DNA metabarcoding on fecal samples, targeting eukaryotic 18S rRNA of potential dietary prey items, and bacterial 16S rRNA of the concomitant gut microbiome. Our results demonstrated an abundance of soft-bodied prey in the diet of fire salamanders, and a significant difference in the diet composition between males and females. This sex-dependent effect on diet was also reflected in the gut microbiome diversity, which is higher in males than female animals. Proximity to human activities was associated with increased intestinal pathogen loads. Collectively, the data supports a relationship between diet, environment and intestinal microbiome in fire salamanders, with potential health implications.

Environment-Dependent Variation in Gut Microbiota of an Oviparous Lizard (Calotes versicolor)

Simple Summary

The different gut sections potentially provide different habitats for gut microbiota. We found that Bacteroidetes, Firmicutes, and Proteobacteria were the three primary phyla in gut microbiota of C. versicolor. The relative abundance of dominant phyla Bacteroidetes and Firmicutes exhibited an increasing trend from the small intestine to the large intestine, and there was a higher abundance of genus Bacteroides (Class: Bacteroidia), Coprobacillus and Eubacterium (Class: Erysipelotrichia), Parabacteroides (Family: Porphyromonadaceae) and Ruminococcus (Family: Lachnospiraceae), and Family Odoribacteraceae and Rikenellaceae in the hindgut, and some metabolic pathways were higher in the hindgut. Our results reveal the variations of gut microbiota composition and metabolic pathways in different parts of the lizards’ intestine.

Abstract

Vertebrates maintain complex symbiotic relationships with microbiota living within their gastrointestinal tracts which reflects the ecological and evolutionary relationship between hosts and their gut microbiota. However, this understanding is limited in lizards and the spatial heterogeneity and co-occurrence patterns of gut microbiota inside the gastrointestinal tracts of a host and variations of microbial community among samples remain poorly understood. To address this issue and provide a guide for gut microbiota sampling from lizards, we investigated the bacteria in three gut locations of the oriental garden lizard (Calotes versicolor) and the data were analyzed for bacterial composition by 16S ribosomal RNA (16S rRNA) gene amplicon sequencing. We found the relative abundance of the dominant phyla exhibited an increasing trend from the small intestine to the large intestine, and phyla Firmicutes, Bacteroidetes and Proteobacteria were the three primary phyla in the gut microbiota of C. versicolor. There were a higher abundance of genus Bacteroides (Class: Bacteroidia), Coprobacillus and Eubacterium (Class: Erysipelotrichia), Parabacteroides (Family: Porphyromonadaceae) and Ruminococcus (Family: Lachnospiraceae), and Family Odoribacteraceae and Rikenellaceae in the sample from the hindgut. The secondary bile acid biosynthesis, glycosaminoglycan degradation, sphingolipid metabolism and lysosome were significantly higher in the hindgut than that in the small intestine. Taken together our results indicate variations of gut microbiota composition and metabolic pathway in different parts of the oriental garden lizard.

Decisive Effects of Life Stage on the Gut Microbiota Discrepancy Between Two Wild Populations of Hibernating Asiatic Toads (Bufo gargarizans)

Until now, the effects of driving factors on the gut microbiota of amphibians are still mostly confounded. Due to a long-term fasting, hibernating amphibians are ideal experimental materials to explore this question. In this study, we characterized the small intestine microbiota of adult hibernating Asiatic toads (Bufo gargarizans) collected from two geographical populations using 16S rRNA amplicon sequencing technique and evaluated the effects of non-dietary factors (e.g., sex and host genetic background). Proteobacteria (0.9196 ± 0.0892) was characterized as the most dominant phylum in the small gut microbiota of hibernating Asiatic toads, among which five core OTUs were identified and three were classified into Pseudomonas. In view of the coincidence between the dominant KEGG pathways (such as the two-component system) and Pseudomonas, Pseudomonas appeared to be a key adaptor for small gut microbiota during hibernation. Furthermore, we detected a greater discrepancy of gut microbiota between geographical populations than between sexes. Both sex and host genetic background showed a minor effect on the gut microbiota variation. Finally, life stage was determined to be the decisive factor driving the gut microbiota discrepancy between populations. However, a large proportion of the gut microbiota variation (∼70%) could not be explained by the measured deterministic factors (i.e., sex, location, body length, and routine blood indices). Therefore, other factors and/or stochastic processes may play key roles in shaping gut bacterial community of hibernating amphibians.

Characterization of the microbiome of the invasive Asian toad in Madagascar across the expansion range and comparison with a native co-occurring species

Biological invasions are on the rise, with each invader carrying a plethora of associated microbes. These microbes play important, yet poorly understood, ecological roles that can include assisting the hosts in colonization and adaptation processes or as possible pathogens. Understanding how these communities differ in an invasion scenario may help to understand the host's resilience and adaptability. The Asian common toad, Duttaphrynus melanostictus is an invasive amphibian, which has recently established in Madagascar and is expected to pose numerous threats to the native ecosystems. We characterized the skin and gut bacterial communities of D. melanostictus in Toamasina (Eastern Madagascar), and compared them to those of a co-occurring native frog species, Ptychadena mascareniensis, at three sites where the toad arrived in different years. Microbial composition did not vary among sites, showing that D. melanostictus keeps a stable community across its expansion but significant differences were observed between these two amphibians. Moreover, D. melanostictus had richer and more diverse communities and also harboured a high percentage of total unique taxa (skin: 80%; gut: 52%). These differences may reflect the combination of multiple host-associated factors including microhabitat selection, skin features and dietary preferences.

Direct Comparison of Fecal and Gut Microbiota in the Blue Mussel (Mytilus edulis) Discourages Fecal Sampling as a Proxy for Resident Gut Community

Bivalves have ecological and economic importance but information regarding their associated microbiomes is lacking. As suspension feeders, bivalves capture and ingest a myriad of particles, and their digestive organs have a high throughput of particle-associated microbiota. To better understand the complement of transient and resident microbial communities, standard methods need to be developed. For example, fecal sampling could represent a convenient proxy for the gut microbiome and is simple, nondestructive, and allows for sampling of individuals through time. The goal of this study was to evaluate fecal sampling as a reliable proxy for gut microbiome assessment in the blue mussel (Mytilus edulis). Mussels were collected from the natural environment and placed into individual sterilized microcosms for 6 h to allow for fecal egestion. Feces and gut homogenates from the same individuals were sampled and subjected to 16S rRNA gene amplicon sequencing. Fecal communities of different mussels resembled each other but did not resemble gut communities. Fecal communities were significantly more diverse, in terms of amplicon sequence variant (ASV) richness and evenness, than gut communities. Results suggested a mostly transient nature for fecal microbiota. Nonetheless, mussels retained a distinct resident microbial community in their gut after fecal egestion that was dominated by ASVs belonging to Mycoplasma. The use of fecal sampling as a nondestructive substitute for direct sampling of the gut is strongly discouraged. Experiments that aim to study solely resident bivalve gut microbiota should employ an egestion period prior to gut sampling to allow time for voidance of transient microbes.