Abstract
The complex gut bacterial communities have a major impact on organismal health. However, knowledge of the effects of habitat change on the gut microbiota of wild birds is limited. In this study, we characterized the gut microbiota of two different subspecies of the Silver-eared Mesia (Leiothrix argentauris), the native subspecies (L. a. rubrogularis) and immigrant subspecies (L. a. vernayi), using 16S rRNA gene high-throughput sequencing. These two subspecies live in a trace metal-contaminated area, and L. a. vernayi was trafficked. They are an excellent system for studying how the gut microbiome of wild animal changes when they move to new habitats. We hypothesized that the immigrant subspecies would develop the same adaptations as the native subspecies in response to habitat changes. The results showed that there were no significant differences in the composition, diversity, or functional metabolism of gut microbiota between native and immigrant subspecies under the combined action of similar influencing factors (the p values of all analyses of variance >0.05). In addition, the composition and functional metabolism of gut microbiota in two subspecies showed adaptation against trace metal damage. Linear discriminant analysis effect size (LEfSe) analysis revealed that Massilia in the intestinal microbiota of immigrant subspecies was significantly higher than that of native subspecies, suggesting that immigrant subspecies suffered habitat change. Finally, we found that these two subspecies living in the mining area had an extremely high proportion of pathogenic bacteria in their gut microbiota (about 90%), much higher than in other species (about 50%) living in wild environment. Our results revealed the adaptation of intestinal microbiota of immigrant Silver-eared Mesias under heavy metals stress, which would provide guidance for biodiversity conservation and pollution management in mining area.
Introduction
Vertebrates exhibit a highly complex symbiotic relationship with their gut microbiota (McFall-Ngai et al., 2013). The gut microbiota is conducive to metabolism, immunity, behavior, and development of the organisms, with positive implications on host health (McFall-Ngai et al., 2013; Youngblut et al., 2019). For instance, the gut microbiome can regulate nutrient absorption efficiency, influence external temperature, and control metabolic rate (; ). Besides, some specific microorganisms can enhance the resistance of the host to toxicity (). In turn, the hosts provide a nutrient-rich, stable habitat for the gut microbiome (). Hosts and their gut microbiota have shaped a complex symbiotic relationship.
The animal’s gut microbiome is derived primarily from the environment (). As a result, the gut microbiome is an unstable element, constantly varying in response to external environment changes (). For example, animals can acquire new gut microbes from the environment (). The relative abundance of gut microbes can be reconfigured according to various environmental factors (). The hosts’ metabolism (), nutrition (Moszak et al., 2020), immunology (), behavior (Morais et al., 2021), morphology () and development () would be affected by gut microbiome changes. Understanding how environmental change dictates the microbiota in the intestine will facilitate the conservation and management of wildlife (; Rosenberg and Zilber-Rosenberg, 2018).
Many factors cause animals to alter their habitat. These include natural migrations, such as the migration of birds and fish, as well as forced migrations, such as captivity, global warming, and animal trade. Generally, intestinal microbiota changes caused by natural migration are usually weak - most bacterial taxa are significantly unaffected (Risely et al., 2017, 2018). Conversely, wildlife will lose their native microbiome and reduce the alpha diversity of gut microbiome when moving to captivity (; ). Forced migrations caused by climate warming often result in gut microbiota variations due to diet changes (,). The illegal trade of animals usually ends in captivities or the release of wildlife to new habitats (Ni et al., 2020). However, previous research always focused on the gut microbiota of animals in captivities (; ). The gut microbiome of illegally traded wildlife released into new natural habitats is poorly known.
The Silver-eared Mesia (Leiothrix argentauris), of least concern threat-status (IUCN), is a resident bird and distributed mainly in forests around China and India, as well as in Sumatra1. This bird does not show distinct sexual dimorphism. Its feathers are brightly colored, varying in different subspecies. Silver-eared Mesia is ubiquitous in international live-bird trade because of its unique characteristics (Li and Jiang, 2014; ). This eventually leads to the traded Silver-eared Mesias being raised in captivity or introduced into natural habitats where other native populations live. In addition, the immigrant and native population of Silver-eared Mesia may live in the same flock because of their gregarious habits. Hence, Silver-eared Mesia is an excellent model to study how the gut microbiota of wildlife would change after moving to a new habitat by international trade.
In this study, we reported the composition of gut microbiota in the two subspecies of Silver-eared Mesia, including native (L. a. rubrogularis) and immigrant subspecies (L. a. vernayi), at a mine area in Southern China. The overall goal was to compare the composition, diversity, and functional characteristics of the gut microbiota of both native (L. a. rubrogularis) and immigrant (L. a. vernayi) populations. In addition, we also analyzed the intestinal microbiome under trace element stress to explore whether the immigrant population enhanced adaptability to the new environment for the host.
Materials and methods
Ethical guidelines
All samples of feathers and cloacal swabs of Silver-eared Mesias were obtained under the permission of Chongzuo Forestry Department (2018. 1), and Animal Ethics Committee, Guangxi University (GXU2018-039), China, and procedures followed the laws of the People’s Republic of China.
Study area overview
We conducted this study in a mine tailing (22°58′15″N, 107°17′28″E) in the Daxin County of Chongzuo Prefecture in Guangxi, Zhuang Autonomous Region, Southern China. The French began mining here more than 100 years ago. However, due to environmental pollution and the depletion of mineral resources, the mine was closed in 2001 (). The region has a subtropical Marine monsoon climate, with an average annual precipitation of 1348.8 mm. Due to water leaching and surface runoff, the downstream rivers and soil are seriously polluted by trace elements. Cadmium levels in the soil were 11.3 times higher than the recommended limit set in 2000 (Pan et al., 2020).
Study species
There are two subspecies of the Silver-eared Mesia in this area. One is a native subspecies, L. a. rubrogularis, which is naturally distributed in the studied area (). Another subspecies, L. a. vernayi, which is distributed initially mainly in Myanmar and Southwestern China, has recently been observed in field monitoring. The plumage characteristics of the two subspecies are similar. The only apparent difference is that L. a. rubrogularis has a red collar and base coverings of upper tail, while L. a. vernayi is orange. The diet of Silver-eared Mesia is uniformly omnivorous, feeding on insects and their larvae, fruits, and seeds.
Sample collection
We captured the two subspecies living in the forests around the tailing pond using mist nets in January 2022. We used passive methods (no birdsong playback), placing nets from dawn to dusk and patrolling them at least once per hour. We eventually captured 19 Silver-eared Mesias (5 L. a. rubrogularis and 14 L. a. vernayi). Once birds are captured, we cleaned the outside of the cloaca with an alcohol pad, inserted a sterile flocking swab (Puritan 25-3,316-U Ultra Flocked Swab, United States) fully into the cloaca, turned for 3–5 s, and preserved the swab in RNAlater (Qiagen, Hilden, Germany). Then, samples were immediately placed in sterile vials, kept in a cool box in the field, and later stored at −20°C. Upon return to the laboratory, all samples were stored at −80°C until processed. Before color-mapping and release, we measured the body index of each bird (including body weight, body length, length of wing and length of tarsometatarsus). Besides, we collected each bird’s primary, secondary, chest and tail feathers to measure trace element levels. The measurement method was described in detail in the Supplementary materal.
DNA extraction, PCR amplification, and amplicon sequencing
Total genome DNA from samples was extracted using cetyltrimethylammonium bromide (CTAB) method. DNA concentration and purity were monitored on 1% agarose gel. According to the concentration, DNA was diluted to 1 ng/μL using sterile water.
PCR amplification of the V3–V4 hypervariable region of bacterial 16S rRNA genes was performed using the bacterial-specific forward primer 341F (5′-CCTAYGGGRBGCASCAG-3′) and reverse primer 805R (5′-GGACTACNNGGGTATCTAAT-3′) with the barcode. All PCR reactions were carried out with 15 μl of Phusion® High-Fidelity PCR Master Mix (New England Biolabs); 2 μM of forward and reverse primers, and about 10 ng template DNA. Thermal cycling consisted of initial denaturation at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, and elongation at 72°C for 30 s, followed by a final extension of 5 min at 72°C. Mix same volume of 1X TAE buffer with PCR products and operate electrophoresis on 2% agarose gel for detection. PCR products was mixed in equidensity ratios. Then, mixture PCR products was purified with Qiagen Gel Extraction Kit (Qiagen, Germany).
Sequencing libraries were generated using TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, United States) following manufacturer’s recommendations and index codes were added. The library quality was assessed on the Qubit@ 2.0 Fluorometer (Thermo Scientific). At last, the library was sequenced on an Illumina NovaSeq platform and 250 bp paired-end reads were generated (a total of 19 samples, 5 L. a. rubrogularis, and 14 L. a. vernayi). Raw sequences obtained in this study are available through the National Center for Biotechnology Information (NCBI) database (accession number PRJNA853520).
16S rRNA gene amplicon sequencing analysis and processing
Microbiome bioinformatics analysis were performed with the QIIME 2 (2022.2) process (). With reference to the review by , we made slightly modified based on the official tutorial2. In short, paired-end reads were merged using ‘vsearch join-pairs’ and quality filtered using ‘quality-filter q-score-joined’ () within QIIME2 after trimming the barcode and primer sequences from the reads. Next, sequences were quality filtered and denoised using the Deblur workflow (). All amplicon sequence variants (ASVs) were aligned with mafft (). Taxonomy was assigned to ASVs using ‘feature-classifier classify-sklearn’ plugin against the pre-trained Naive Bayes classifier () [based on SILVA v138 database (Quast et al., 2012)]. Then, the ASVs containing less than 0.1% of total sequences across all samples, any contaminating mitochondrial and chloroplast sequences and sequences not assigned to phylum were filtered out from the denoising results. Next, ‘fragment-insertion sepp’ (Mirarab et al., 2012) is used to generate our phylogenetic tree by inserting short molecular sequences into an existing phylogenetic tree and aligning each short molecular sequence to the aligment of the full-length sequences, then that alignment is used to find the optimal location in the phylogenetic tree for querying sequences. Finally, to minimize the difference of sequencing depth across samples, data was rarefied to a sampling depth of 17,400 reads per sample for the downstream diversity analysis.
Metagenomic pathway prediction by PICRUSt2
The 16S rRNA gene sequence of the sample was mapped to KEGG3 database for annotation, and the abundance of metabolic pathways was predicted based on Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) (). R software (v4.1.2) was utilized for statistical analyses and visualization of the identified pathways. Correlations between samples were visualized using the pheatmap package. Linear discriminant analysis effect size (LEfSe) analysis was used to detect the differentially abundant metabolic pathways between L. a. rubrogularis and L. a. vernayi.
Statistical analysis
Statistical analysis was completed in R v4.1.2 (R Core Team, 2021). Alpha diversity analysis was carried out using the “picante” package (v1.8.2) (). Differences in Alpha diversity between L. a. rubrogularis and L. a. vernayi were evaluated by one-way ANOVA (), and all p values were corrected via the false discovery rate (FDR) correction. Rarefaction curves were generated by “ggplot2” package (v3.3.5) (Wickham, 2016). The ASV-level ranked abundance curves were generated to compare the richness and evenness of ASVs among samples. Venn diagram was generated to visualize the shared and unique ASVs among groups using the R package “VennDiagram” (v1.7.3) (), based on the occurrence of ASVs across samples/groups regardless of their relative abundance. Beta diversity analysis was performed to investigate the structural variation in microbial communities across samples using weighted and unweighted UniFrac distance metrics (Lozupone et al., 2011) and visualized via principal coordinate analysis (PCoA) (Ramette, 2007). Differences in the UniFrac distances for pairwise comparisons among groups were determined by Anosim test with 999 permutations using “vegan” package (v2.5–7) (). The taxonomic composition of each sample was visualized as a stacked bar plot at the phylum, family and genus level with the “ggplot2” package (v3.3.5) (Wickham, 2016). The circles plot of L. a. rubrogularis and L. a. vernayi were generated by Circos4 (Krzywinski et al., 2009). Correlations between samples were visualized using the pheatmap package (v1.0.12) (Raivo, 2019). Linear discriminant analysis effect size (LEfSe) was performed to detect differentially abundant taxa across groups using the default parameters (Segata et al., 2011).
Results
Sequencing quality and ASVs distribution
After strict quality filtering, a total of 482,768 16S rRNA gene reads were obtained, with an average of 25,409 ± 3,882 sequences per sample (median: 25785), from which 501 ASVs were identified to be from gut microbiotas of 19 Silver-eared Mesias (L. a. rubrogularis = 462, L. a. vernayi = 501). The rank abundance and rarefaction curves supported the adequacy of the sequencing depth, indicating that almost all the bacterial species were found in intestinal samples (Figures 1A–D). A total of 92.22% of the ASVs were shared among the two groups, with unique ASVs entirely being present in the L. a. vernayi (7.78%). All ASVs can be detected in L. a. vernayi (Supplementary Figure S1).
Figure 1
Characteristics of gut microbial diversity of L. a. rubrogularis and L. a. vernayi
To evaluate the differences in community richness and diversity among the subspecies, the effective sequences were aligned to calculate the Good’s coverage, Richness, Shannon, and PD_whole_tree indices. Good’s coverage estimates of each sample ranged from 99.82 to 99.97%, suggesting excellent coverage (Figure 2A). Measurements of gut community alpha diversity did not vary significantly between L. a. rubrogularis and L. a. vernayi (based on one-way ANOVA, Richness index, value of p = 0.4350; Shannon index, value of p = 0.876; PD_whole_tree, value of p = 0.329). Thus, the diversity of intestinal microbiota between L. a. rubrogularis and L. a. vernayi was similar (Figures 2B–D).
Figure 2
Weighted and unweighted Unifrac distances are used to generate PCoA plots reflecting the gut microbial beta-diversity between individuals. No significant difference between L. a. rubrogularis and L. a. vernayi was detected no matter what weighted Unifrac distance or unweighted Unifrac distance was considered (Figures 3A,B, Anosim test: value of p = 0.251 and 0.217, R = 0.0945 and 0.1084, respectively).
Figure 3
Characteristics of gut microbiota composition of L. a. rubrogularis and L. a. vernayi
Nearly all reads were assignable to 13 phyla, 17 classes, 54 orders, 97 families, and 162 genera. At the phylum level, Proteobacteria were the most abundant bacterial phylum (88.94, 93.46%) among L. a. rubrogularis and L. a. vernayi gut communities followed by Firmicutes (6.05, 3.12%), Bacteroidetes (2.51, 1.24%), and Actinobacteria (1.65, 1.61%) (Figure 4A). The relative abundance of the remaining others phyla was low. All the phyla of gut microbiota in L. a. rubrogularis and L. a. vernayi were shown in Figure 4B. The relative abundance of each phylum across samples fluctuated little.
Figure 4
At the family level, Pseudomonadaceae predominated (66.04, 63.87%) among L. a. rubrogularis and L. a. vernayi gut communities followed by Enterobacteriaceae (8.80, 9.72%), Vibrionaceae (6.71, 7.48%), Comamonadaceae (1.96, 1.80%), Lachnospiraceae (1.27, 1.38%), Sphingomonadaceae (1.22, 1.25%), Rhodocyclaceae (0.87, 0.84%), Halomonadaceae (0.85, 0.94%), Aeromonadaceae (0.73, 0.63%) and Burkholderiaceae (0.66, 0.66%) (Figure 4C). The top 10 families of gut microbiota in L. a. rubrogularis and L. a. vernayi were shown in Figure 4D. We observed that the abundance of Enterobacteriaceae in two samples of L. a. vernayi was higher than that in other samples, and even exceeded that of Pseudomonadaceae.
At the genus level, Pseudomonas predominated (66.04, 63.87%) among L. a. rubrogularis and L. a. vernayi gut communities followed by Shigella (7.18, 7.90%), Vibrio (6.56, 7.31%), Halomonas (0.85, 0.94), Novosphingobium (0.84, 0.86), Methyloversatilis (0.81, 0.78%), Aeromonas (0.71, 0.61), Hydrogenophaga (0.65, 0.63%), Streptococcus (0.64, 0.70%), and hgcl_clade (0.63, 0.65%) (Figure 4E). The top 10 genera of gut microbiota in L. a. rubrogularis and L. a. vernayi were shown in Figure 4F. The abundance of Shigella was higher in two samples of L. a. vernayi than in the other samples.
Difference analysis of the gut microflora between L. a. rubrogularis and L. a. vernayi
To further investigate whether there were differences in intestinal microbial community composition between L. a. rubrogularis and L. a. vernayi. The top 40 genera with average abundance were used to construct the heatmap of the genus-level hierarchical clustering. Although the clustering was divided into two categories, L. a. rubrogularis and L. a. vernayi were not, respectively, assigned to the two categories, and the distribution of the two groups was very dispersed (Figure 5A). LEfSe analysis was used to determine the microbial communities with significant differences in abundance between the two groups across the phylum, class, order, family, and genus levels. Interestingly, we identified only one genus taxon, Massilia, as having a significantly higher abundance in L. a. vernayi than in L. a. rubrogularis (LDA > 2, value of p < 0.05) (Figure 5B). Lefse analysis confirmed the results of genus clustering and showed that the composition of intestinal flora in L. a. rubrogularis and L. a. vernayi was basically consistent.
Figure 5
Prediction and comparisons in the gut microbial metabolism pathways
According to the prediction, the KEGG Pathway Database classified metabolic pathways into six categories, including cellular processes, environmental information processing, genetic information processing, human diseases, metabolism, and organismal systems. Each metabolic pathway was further divided into several levels. The second level included 31 metabolic pathway sub-functions, and a total of 159 the third level functional pathways was predicted. Figure 6A demonstrates the average relative abundance of bacteria in all individuals of Silver-eared Mesia mapped to the secondary functional pathway of the KEGG database. The KEGG pathways showed that ‘metabolism’ was the most important pathway of gut microbiota.
Figure 6
Wilcoxon rank-sum test showed no significant difference in secondary metabolic pathways between L. a. rubrogularis and L. a. vernayi (Figure 6B) (secondary metabolic pathways with relative abundance of less than 1% in all samples were excluded and the p-values were corrected by Bonferroni). Clustering heat maps of the abundances of the top 40 third level functional pathways showed differences in functional abundances across samples, but the individuals of L. a. rubrogularis and L. a. vernayi were not clustered into two categories, respectively (Figure 6C). The results demonstrated no significant difference in the metabolic function of the intestinal microbiota between L. a. rubrogularis and L. a. vernayi. The LEfSe analysis was applied to all levels of KEGG metabolic pathways and also did not detect any variation between L. a. rubrogularis and L. a. vernayi (Supplementary Figure S2).
Discussion
This is the first study to show the gut microbiota of Silver-eared Mesia whose habitat has been altered by trade. Previous studies have focused on comparing native habitats with captivity (Wienemann et al., 2011; Sun et al., 2019; Oliveira et al., 2020; San Juan et al., 2021; ; Zhang et al., 2022). However, under captive conditions, the host’s living environment and diet are artificially controlled, so it is difficult to reflect the impact of the real field environment on the intestinal microbiota of immigrant species. In addition, under natural conditions, the host’s habitat environment and food sources are variable, such as the significant differences in diet composition between individuals within a species (), which is difficult to overcome for cross-species and cross-regional comparison. In this study, native (L. a. rubrogularis) and immigrant (L. a. vernayi) subspecies lived in the same group, so their diet, environment, and other external factors were the same, and there were only slight phylogenetic differences, which can supply a unique opportunity to compare the gut microbiota of immigrant and native subspecies. Wildlife trade is a common route for introducing invasive species (). Traded migrant species may adversely affect native ecosystems, economic activity, and human well-being (Simberloff et al., 2013). Our findings can help to assess the risk of invasion of this non-native species for better understanding and management.
The gut microbiota of immigrant subspecies coincides with the native subspecies
On the whole, our results demonstrated no significant differences in the composition, diversity, or functional metabolism of the gut microbiota between native and immigrant subspecies, which might be related to a similar diet, environment, and other external factors. The convergent evolution of gut microbiome, driven by similar dietary and environmental factors, has been found in other vertebrates as well. For example, the bamboo-eating Ailuropoda melanoleuca and Ailurus styani share more similarities in their gut microbiota structure and function with each other than their carnivorous relatives (). Similarly, this phenomenon has been observed in studies of immigrants of different races. The gut microbiota of immigrants in the United States would become more and more similar to that of native people over time (Vangay et al., 2018; Peters et al., 2020; ). To sum up, the intestinal microbiota of immigrant populations developed the same adaptive characteristics as that of native populations, and the generation of these adaptive characteristics may help immigrant populations adapt quickly to habitat changes ().
Characteristics of the gut microbiota of the Silver-eared Mesia living in mining area
Trace metals in the soil can be concentrated in the plants, which insects then feed on and carry trace elements into their bodies (Zhang et al., 2017). Thus, the extensive diet of Silver-eared Mesia will enrich a mass of trace elements into the body, making its intestinal microbiota under trace element stress environment. As a result, we also investigated the characteristics of its intestinal microbiota under trace elements stress. From the community composition of the gut microbes of native and immigrant subspecies, it was found that at the phylum level, Proteobacteria dominated the intestinal flora of the two subspecies (about 90%), followed by Firmicutes, Bacteroidetes, and Actinobacteria. Previous studies have found that Firmicutes are the main microbiome in the guts of mammals, chickens and wild birds (Waite and Taylor, 2014; ). Members of Firmicutes play an important role in the metabolism, digestion and absorption of proteins and other substances nutrients, and participate in the synthesis of digestive enzymes to assist the host in the digestion and absorption of nutrients (). Proteobacteria, however, are mostly pathogenic bacteria, usually associated with intestinal ecological imbalances, metabolic and immune disorders (). Some studies of other wild birds have found Firmicutes to be the dominant species in the gut microbiota (Wang et al., 2017; Oliveira et al., 2020; Wu et al., 2021). In several other studies on wild birds, nevertheless, Proteobacteria were found to have a higher abundance than Firmicutes (Zhou et al., 2020; Zhang et al., 2022). Even so, compared with our study in which Proteobacteria accounted for almost 90% of the intestinal flora, Proteobacteria in their study were only slightly higher than Firmicutes. Given that the samples in these studies came from captivity or normal wild environments, it may be suggested that adaptation of intestinal flora to trace element stress may lead to Proteobacteria dominance. This is consistent with the characteristics of intestinal microbiota in other species under trace element stress (Šrut et al., 2019; ; Wu et al., 2020).
The role of intestinal microbiota in reducing heavy metal toxicity
Previous studies have shown that the gut microbiome is essential in protecting the host from the toxic effects of heavy metal exposure (; ; Ninkov et al., 2015). In our study, we also found that the dominant group in the intestinal microbiome showed strong resistance to heavy metal toxicity. For instance, the Pseudomonas was dominant among almost all samples. Pseudomonas can degrade cellulose () and produce multitudinous antibiotic compounds, which can effectively stem the diseases caused by pathogenic bacteria and fungi (). Particularly important, Pseudomonas showed a sky-high degree of resistance to heavy metals, such as Cd, Pb and As (; Liu et al., 2021; Pramanik et al., 2021; Sun et al., 2021). In addition, several genera with high abundance, such as Shigella, Vibrio and Halomonas, also have high trace element resistance (Mukherjee et al., 2019; ; ). These bacteria may play an important role in resisting trace element stress and maintaining normal growth, development and life activities of Silver-eared Mesia.
We also predicted the metabolic function of the gut microbiota of the Silver-eared Mesia. The metabolic pathways predicted by PICRUSt2 were consistent with those predicted by other wild species (Figure 6A; Zhou et al., 2020; ; Zhang et al., 2022). Microbial metabolic pathways (such as xenobiotics biodegradation metabolism, carbohydrate metabolism, amino acid metabolism, and metabolism of cofactors and vitamins) accounted for more than 60% of the 31 predicted most abundant pathways. This suggests that the gut microbiome of Silver-eared Mesia may be involved in high levels of metabolic activity that may help to resist damage caused by trace elements.
Speculation about when L. a. vernayi immigrated
We cannot evaluate the specific time of immigration of immigranted L. a. vernayi. According to our field survey results and some data () in recent years, we can infer that its earliest migration time was around 2019, because it has not been found in the field survey before 2019. Interestingly, the Lefse analysis of our taxa may support this inference. Among all taxa of L. a. rubrogularis and L. a. vernayi, lefse analysis detected a genus, Massilia, with a significantly higher relative abundance in L. a. vernayi than that in L. a. rubrogularis. The Massilia are important resistant microbes in trace element stress environments, as it has been identified in mining soil, farmland soils, beach, and sludge polluted by trace elements (; ; Lee et al., 2017; Wang et al., 2020). The Massilia can secrete a large amount of cyclodextrin, whose special physical structure can contain trace element ions (Santos and Barbosa-Tessmann, 2019), thus reducing the content of trace element in the intestinal environment and reducing its toxicity to Silver-eared Mesia (; Tajuddin Sikder et al., 2014). Another study on earthworms (Eisenia fetida) also found that Massilia showed a higher abundance after cadmium contamination, and that the abundance of Massilia decreased in the middle stage of cadmium pollution and increased to a stable level in the later stage (Zhou et al., 2021). This finding suggests that L. a. vernayi did not migrate for very long. We speculate that in the future, the abundance of Massilia in the intestinal tract of L. a. vernayi will fluctuate until it becomes similar to that of L. a. rubrogularis.
Limitations
It should be noted that this study has several limitations. First, we do not understand the microbiome of the immigrant form, prior to our sampling. Secondly, we did not observe it in non-mining areas perhaps due to its rarity in China (second class national protected animal), which led to the inability to observe the adaptation process of the intestinal microbiota of Silver-eared Mesia under trace element stress. In addition, there are some differences in the sample numbers of L. a. rubrogularis and L. a. vernayi. Although the intestinal microbiota of L. a. rubrogularis and L. a. vernayi are very similar in our study, the differences in sample numbers may also cause some deviations in the results. From another perspective, however, the different sample numbers of the two subspecies may reflect their proportions in the wild-living taxa. There are far more immigrant subspecies than native ones, which could be detrimental to our biodiversity conservation. We suggest that some controlled experiments in non-mining areas and the habitat of L. a. vernayi should be carried out in future studies, which will be of significance for the protection of wild animals.
Conclusion
In conclusion, the composition, diversity, and function of the intestinal microbiota of immigrant subspecies (L. a. vernayi) were basically the same as those of the native subspecies (L. a. rubrogularis). The gut microbiota of the Silver-eared Mesia living in the mining area developed adaptive characteristics to trace element stress. The significance of this study lies in: (1) the results shed some light on the adaptation of the gut microbiota of immigrant species to cope with habitat changes after “migration”. It supplies new insights into the adaptation mechanism of gut microbiota in response to environmental changes in wild animals, thus providing important implications for wildlife conservation and biodiversity conservation. (2) The results also showed the adaptation mechanism of the intestinal microbiota of birds under trace element environmental stress (how to change the composition and function of the microbiota to resist the toxic effects of trace elements). It is helpful to understand the toxicity of trace elements to birds and the strategies of birds to adapt to trace element pollution, which has a guiding role in pollution control and biodiversity protection in mining areas.
Funding
This work was funded by the National Natural Science Foundation of China (31870370) and the Key Grant of Guangxi Nature and Science Foundation (2018GXNSFDA281016).
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, PRJNA853520.
Ethics statement
The animal study was reviewed and approved by Chongzuo Forestry Bureau and Ethics Committee of Guangxi University.
Author contributions
TZ, SL, and AJ contributed intellectual input and assistance to this study. TZ and AJ designed the research. TZ and SL did the sampling and laboratory work. TZ conducted data analysis and wrote the first draft of the manuscript. AJ contributed substantially to revisions. All authors contributed to the article and approved the submitted version.
Acknowledgments
We thank Madhava Meegaskumbura for his help in the manuscript revision process. We are grateful to the Chongzuo Forestry Bureau for granting permission for this study. We thank Wekemo Tech Group Co., Ltd. Shenzhen China for its sequencing services.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1076523/full#supplementary-material
References
1
Al-AnsariM. M.BenabdelkamelH.AlMalkiR. H.Abdel RahmanA. M.AlnahmiE.MasoodA.et al. (2021). Effective removal of heavy metals from industrial effluent wastewater by a multi metal and drug resistant Pseudomonas aeruginosa strain RA-14 using integrated sequencing batch reactor. Environ. Res.199:111240. doi: 10.1016/j.envres.2021.111240
2
AlberdiA.AizpuruaO.BohmannK.Zepeda-MendozaM. L.GilbertM. T. P. (2016). Do vertebrate gut metagenomes confer rapid ecological adaptation?Trends Ecol. Evol.31, 689–699. doi: 10.1016/j.tree.2016.06.008
3
AlberdiA.Martin BidegurenG.AizpuruaO. (2021). Diversity and compositional changes in the gut microbiota of wild and captive vertebrates: a meta-analysis. Sci. Rep.11:22660. doi: 10.1038/s41598-021-02015-6
4
AmirA.McDonaldD.Navas-MolinaJ. A.KopylovaE.MortonJ. T.XuZ. Z.et al. (2017). Deblur rapidly resolves single-nucleotide community sequence patterns. mSystems2:e00191-16. doi: 10.1128/mSystems.00191-16
5
BadruddozaA. Z. M.ShawonZ. B. Z.TayW. J. D.HidajatK.UddinM. S. (2013). Fe3O4/cyclodextrin polymer nanocomposites for selective heavy metals removal from industrial wastewater. Carbohydr. Polym.91, 322–332. doi: 10.1016/j.carbpol.2012.08.030
6
BodawattaK. H.HirdS. M.GrondK.PoulsenM.JønssonK. A. (2022a). Avian gut microbiomes taking flight. Trends Microbiol.30, 268–280. doi: 10.1016/j.tim.2021.07.003
7
BodawattaK. H.KlečkováI.KlečkaJ.PužejováK.KoaneB.PoulsenM.et al. (2022b). Specific gut bacterial responses to natural diets of tropical birds. Sci. Rep.12:713. doi: 10.1038/s41598-022-04808-9
8
BokulichN. A.KaehlerB. D.RideoutJ. R.DillonM.BolyenE.KnightR.et al. (2018). Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome6:90. doi: 10.1186/s40168-018-0470-z
9
BokulichN. A.SubramanianS.FaithJ. J.GeversD.GordonJ. I.KnightR.et al. (2013). Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods10, 57–59. doi: 10.1038/nmeth.2276
10
BolyenE.RideoutJ. R.DillonM. R.BokulichN. A.AbnetC. C.Al-GhalithG. A.et al. (2019). Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol.37, 852–857. doi: 10.1038/s41587-019-0209-9
11
BombaywalaS.PurohitH. J.DafaleN. A. (2021). Mobility of antibiotic resistance and its co-occurrence with metal resistance in pathogens under oxidative stress. J. Environ. Manag.297:113315. doi: 10.1016/j.jenvman.2021.113315
12
BretonJ.DanielC.DewulfJ.PothionS.FrouxN.SautyM.et al. (2013). Gut microbiota limits heavy metals burden caused by chronic oral exposure. Toxicol. Lett.222, 132–138. doi: 10.1016/j.toxlet.2013.07.021
13
BroderickN. A.BuchonN.LemaitreB.McFall-NgaiM. J. (2014). Microbiota-induced changes in Drosophila melanogaster host gene expression and gut morphology. mBio5, e01117–e01114. doi: 10.1128/mBio.01117-14
14
CandelaM.BiagiE.MaccaferriS.TurroniS.BrigidiP. (2012). Intestinal microbiota is a plastic factor responding to environmental changes. Trends Microbiol.20, 385–391. doi: 10.1016/j.tim.2012.05.003
15
CardadorL.TellaJ. L.AnadónJ. D.AbellánP.CarreteM. (2019). The European trade ban on wild birds reduced invasion risks. Conserv. Lett.12:e12631. doi: 10.1111/conl.12631
16
ChenS.HolyoakM.LiuH.BaoH.MaY.DouH.et al. (2022a). Effects of spatially heterogeneous warming on gut microbiota, nutrition and gene flow of a heat-sensitive ungulate population. Sci. Total Environ.806:150537. doi: 10.1016/j.scitotenv.2021.150537
17
ChenS.HolyoakM.LiuH.BaoH.MaY.DouH.et al. (2022b). Global warming responses of gut microbiota in moose (Alces alces) populations with different dispersal patterns. J. Zool.318, 63–73. doi: 10.1111/jzo.12998
18
ChevalierC.StojanovićO.ColinD. J.Suarez-ZamoranoN.TaralloV.Veyrat-DurebexC.et al. (2015). Gut microbiota orchestrates energy homeostasis during cold. Cells163, 1360–1374. doi: 10.1016/j.cell.2015.11.004
19
ColstonT. J.JacksonC. R. (2016). Microbiome evolution along divergent branches of the vertebrate tree of life: what is known and unknown. Mol. Ecol.25, 3776–3800. doi: 10.1111/mec.13730
20
CopelandJ. K.ChaoG.VanderhoutS.ActonE.WangP. W.BenchimolE. I.et al. (2021). The impact of migration on the gut metagenome of south Asian Canadians. Gut Microbes13, 1–29. doi: 10.1080/19490976.2021.1902705
21
DallasJ. W.WarneR. W. (2022). Captivity and animal microbiomes: potential roles of microbiota for influencing animal conservation. Microb. Ecol. doi: 10.1007/s00248-022-01991-0
22
DouglasG. M.MaffeiV. J.ZaneveldJ. R.YurgelS. N.BrownJ. R.TaylorC. M.et al. (2020). PICRUSt2 for prediction of metagenome functions. Nat. Biotechnol.38, 685–688. doi: 10.1038/s41587-020-0548-6
23
EatonJ.ShepherdC.RheindtF.HarrisJ. B. C.van BalenB.WilcoveD.et al. (2015). Trade-driven extinctions and near-extinctions of avian taxa in Sundaic Indonesia. Forktail31, 1–12.
24
FanY.PedersenO. (2021). Gut microbiota in human metabolic health and disease. Nat. Rev. Microbiol.19, 55–71. doi: 10.1038/s41579-020-0433-9
25
FengG.-D.YangS.-Z.LiH.-P.ZhuH.-H. (2016). Massilia putida sp. nov., a dimethyl disulfide-producing bacterium isolated from wolfram mine tailing. Int. J. Syst. Evol. Microbiol.66, 50–55. doi: 10.1099/ijsem.0.000670
26
GilbertS. F.BoschT. C. G.Ledón-RettigC. (2015). Eco-Evo-devo: developmental symbiosis and developmental plasticity as evolutionary agents. Nat. Rev. Genet.16, 611–622. doi: 10.1038/nrg3982
27
GrievesL. A.BottiniC. L. J.GloorG. B.MacDougall-ShackletonE. A. (2022). Uropygial gland microbiota differ between free-living and captive songbirds. Sci. Rep.12:18283. doi: 10.1038/s41598-022-22425-4
28
GrondK.SandercockB. K.JumpponenA.ZeglinL. H. (2018). The avian gut microbiota: community, physiology and function in wild birds. J. Avian Biol.49:e01788. doi: 10.1111/jav.01788
29
HaasD.DéfagoG. (2005). Biological control of soil-borne pathogens by fluorescent pseudomonads. Nat. Rev. Microbiol.3, 307–319. doi: 10.1038/nrmicro1129
30
HalttunenT.ColladoM. C.El-NezamiH.MeriluotoJ.SalminenS. (2008). Combining strains of lactic acid bacteria may reduce their toxin and heavy metal removal efficiency from aqueous solution. Lett. Appl. Microbiol.46, 160–165. doi: 10.1111/j.1472-765X.2007.02276.x
31
HanboC. (2022). VennDiagram: Generate High-Resolution Venn and Euler Plots. R package version 1.7.3. URL: https://CRAN.R-project.org/package=VennDiagram
32
HeC.SuT.LiuS.JiangA.GoodaleE.QiuG. (2020). Heavy metal, arsenic, and selenium concentrations in bird feathers from a region in southern China impacted by intensive mining of nonferrous metals. Environ. Toxicol. Chem.39, 371–380. doi: 10.1002/etc.4622
33
HehemannJ.-H.CorrecG.BarbeyronT.HelbertW.CzjzekM.MichelG. (2010). Transfer of carbohydrate-active enzymes from marine bacteria to Japanese gut microbiota. Nature464, 908–912. doi: 10.1038/nature08937
34
HuangG.WangX.HuY.WuQ.NieY.DongJ.et al. (2021). Diet drives convergent evolution of gut microbiomes in bamboo-eating species. Sci. China Life Sci.64, 88–95. doi: 10.1007/s11427-020-1750-7
35
JariO.BlanchetF. G.MichaelF.RoelandK.PierreL.DanM.et al. (2020). Vegan: Community Ecology Package. R package version 2.5-7. URL: https://CRAN.R-project.org/package=vegan
36
JiangF.GaoH.QinW.SongP.WangH.ZhangJ.et al. (2021). Marked seasonal variation in structure and function of gut microbiota in Forest and alpine musk deer. Front. Microbiol.12:699797. doi: 10.3389/fmicb.2021.699797
37
JiangA.QinC.HuangL. (2021). Field Guide to the Birds of Guangxi. Guangxi: Guangxi Science and Technology Press.
38
JiménezD. J.Dini-AndreoteF.van ElsasJ. D. (2014). Metataxonomic profiling and prediction of functional behaviour of wheat straw degrading microbial consortia. Biotechnol. Biofuels Bioprod.7:92. doi: 10.1186/1754-6834-7-92
39
JoS.ShinC.ShinY.KimP. H.ParkJ. I.KimM.et al. (2020). Heavy metal and antibiotic co-resistance in Vibrio parahaemolyticus isolated from shellfish. Mar. Pollut. Bull.156:111246. doi: 10.1016/j.marpolbul.2020.111246
40
JobsonJ. D. (1991). “Analysis of variance and experimental design” in Applied Multivariate Data Analysis: Regression and Experimental Design. ed. JobsonJ. D. (New York, NY: Springer New York), 399–544.
41
JumpertzR.LeD. S.TurnbaughP. J.TrinidadC.BogardusC.GordonJ. I.et al. (2011). Energy-balance studies reveal associations between gut microbes, caloric load, and nutrient absorption in humans. Am. J. Clin. Nutr.94, 58–65. doi: 10.3945/ajcn.110.010132
42
KakadeA.SalamaE.-S.PengyaF.LiuP.LiX. (2020). Long-term exposure of high concentration heavy metals induced toxicity, fatality, and gut microbial dysbiosis in common carp, Cyprinus carpio. Environ. Pollut266:115293. doi: 10.1016/j.envpol.2020.115293
43
KatohK.MisawaK.KumaK. I.MiyataT. (2002). MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res.30, 3059–3066. doi: 10.1093/nar/gkf436
44
KembelS. W.CowanP. D.HelmusM. R.CornwellW. K.MorlonH.AckerlyD. D.et al. (2010). Picante: R tools for integrating phylogenies and ecology. Bioinformatics26, 1463–1464. doi: 10.1093/bioinformatics/btq166
45
KnightR.VrbanacA.TaylorB. C.AksenovA.CallewaertC.DebeliusJ.et al. (2018). Best practices for analysing microbiomes. Nat. Rev. Microbiol.16, 410–422. doi: 10.1038/s41579-018-0029-9
46
KohlK. D.WeissR. B.CoxJ.DaleC.Denise DearingM. (2014). Gut microbes of mammalian herbivores facilitate intake of plant toxins. Ecol. Lett.17, 1238–1246. doi: 10.1111/ele.12329
47
KrishnamoorthyR.KimK.SubramanianP.SenthilkumarM.AnandhamR.SaT. (2016). Arbuscular mycorrhizal fungi and associated bacteria isolated from salt-affected soil enhances the tolerance of maize to salinity in coastal reclamation soil. Agric. Ecosyst. Environ.231, 233–239. doi: 10.1016/j.agee.2016.05.037
48
KrzywinskiM. I.ScheinJ. E.BirolI.ConnorsJ.GascoyneR.HorsmanD.et al. (2009). Circos: an information aesthetic for comparative genomics. Genome Res.19, 1639–1645. doi: 10.1101/gr.092759.109
49
LeeH.KimD.-U.ParkS.YoonJ.-H.KaJ.-O. (2017). Massilia chloroacetimidivorans sp. nov., a chloroacetamide herbicide-degrading bacterium isolated from soil. Antonie Van Leeuwenhoek110, 751–758. doi: 10.1007/s10482-017-0845-3
50
LiL.JiangZ. (2014). International trade of CITES listed bird species in China. PLoS One9:e85012. doi: 10.1371/journal.pone.0085012
51
LiuC.LinH.LiB.DongY.Gueret Yadiberet MenzembereE. R. (2021). Endophyte pseudomonas putida enhanced Trifolium repens L. growth and heavy metal uptake: a promising in-situ non-soil cover phytoremediation method of nonferrous metallic tailing. Chemosphere272:129816. doi: 10.1016/j.chemosphere.2021.129816
52
LozuponeC.LladserM. E.KnightsD.StombaughJ.KnightR. (2011). UniFrac: an effective distance metric for microbial community comparison. ISME J.5, 169–172. doi: 10.1038/ismej.2010.133
53
McFall-NgaiM.HadfieldM. G.BoschT. C. G.CareyH. V.Domazet-LošoT.DouglasA. E.et al. (2013). Animals in a bacterial world, a new imperative for the life sciences. Proc. Nat1. Acad. Sci. U.S.A.110, 3229–3236. doi: 10.1073/pnas.1218525110
54
MirarabS.NguyenN.WarnowT. (2012). SEPP: SATé-enabled phylogenetic placement. Pac Symp Biocomput, 247–258. doi: 10.1142/9789814366496_0024
55
MoraisL. H.SchreiberH. L.MazmanianS. K. (2021). The gut microbiota–brain axis in behaviour and brain disorders. Nat. Rev. Microbiol.19, 241–255. doi: 10.1038/s41579-020-00460-0
56
MoszakM.SzulińskaM.BogdańskiP. (2020). You are what you eat—the relationship between diet, microbiota, and metabolic disorders—a review. Nutrients12:1096. doi: 10.3390/nu12041096
57
MukherjeeP.MitraA.RoyM. (2019). Halomonas Rhizobacteria of Avicennia marina of Indian Sundarbans promote rice growth under saline and heavy metal stresses through exopolysaccharide production. Front. Microbiol.10:1207. doi: 10.3389/fmicb.2019.01207
58
NiQ.HeX.ZengB.MengX.XuH.LiY.et al. (2020). Variation in gut microbiota of captive Bengal slow Lorises. Curr. Microbiol.77, 2623–2632. doi: 10.1007/s00284-020-02035-x
59
NinkovM.Popov AleksandrovA.DemeneskuJ.MirkovI.MileusnicD.PetrovicA.et al. (2015). Toxicity of oral cadmium intake: impact on gut immunity. Toxicol. Lett.237, 89–99. doi: 10.1016/j.toxlet.2015.06.002
60
OliveiraB. C. M.MurrayM.TsengF.WidmerG. (2020). The fecal microbiota of wild and captive raptors. Anim Microbiome.2:15. doi: 10.1186/s42523-020-00035-7
61
PanJ. Q.LuoQ. W.HuangZ. J.QianP. R.HuangJ. L.WeiQ. S.et al. (2020). Analysis of cadmium pollution in Xijiang River basin of Guangxi (in Chinese with English abstract). Technol. Innov. Appl.2, 61–63.
62
PetersB. A.YiS. S.BeasleyJ. M.CobbsE. N.ChoiH. S.BeggsD. B.et al. (2020). US nativity and dietary acculturation impact the gut microbiome in a diverse US population. ISME J.14, 1639–1650. doi: 10.1038/s41396-020-0630-6
63
PramanikK.MandalS.BanerjeeS.GhoshA.MaitiT. K.MandalN. C. (2021). Unraveling the heavy metal resistance and biocontrol potential of pseudomonas sp. K32 strain facilitating rice seedling growth under cd stress. Chemosphere274:129819. doi: 10.1016/j.chemosphere.2021.129819
64
QuastC.PruesseE.YilmazP.GerkenJ.SchweerT.YarzaP.et al. (2012). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res.41, D590–D596. doi: 10.1093/nar/gks1219
65
R Core Team (2021). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/
66
RaivoK. (2019). Pheatmap: Pretty Heatmaps. R package version 1.0.12. https://CRAN.R-project.org/package=pheatmap
67
RametteA. (2007). Multivariate analyses in microbial ecology. FEMS Microbiol. Ecol.62, 142–160. doi: 10.1111/j.1574-6941.2007.00375.x
68
RiselyA.WaiteD. W.UjvariB.HoyeB. J.KlaassenM. (2018). Active migration is associated with specific and consistent changes to gut microbiota in Calidris shorebirds. J. Anim. Ecol.87, 428–437. doi: 10.1111/1365-2656.12784
69
RiselyA.WaiteD.UjvariB.KlaassenM.HoyeB. (2017). Gut microbiota of a long-distance migrant demonstrates resistance against environmental microbe incursions. Mol. Ecol.26, 5842–5854. doi: 10.1111/mec.14326
70
RosenbergE.Zilber-RosenbergI. (2018). The hologenome concept of evolution after 10 years. Microbiome6:78. doi: 10.1186/s40168-018-0457-9
71
San JuanP. A.CastroI.DhamiM. K. (2021). Captivity reduces diversity and shifts composition of the Brown kiwi microbiome. Anim Microbiome.3:48. doi: 10.1186/s42523-021-00109-0
72
SantosF. C. D.Barbosa-TessmannI. P. (2019). Recombinant expression, purification, and characterization of a cyclodextrinase from Massilia timonae. Protein Expr. Purif.154, 74–84. doi: 10.1016/j.pep.2018.08.013
73
SegataN.IzardJ.WaldronL.GeversD.MiropolskyL.GarrettW. S.et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol.12:R60. doi: 10.1186/gb-2011-12-6-r60
74
SimberloffD.MartinJ.-L.GenovesiP.MarisV.WardleD. A.AronsonJ.et al. (2013). Impacts of biological invasions: what's what and the way forward. Trends Ecol. Evol.28, 58–66. doi: 10.1016/j.tree.2012.07.013
75
ŠrutM.MenkeS.HöcknerM.SommerS. (2019). Earthworms and cadmium – heavy metal resistant gut bacteria as indicators for heavy metal pollution in soils?Ecotoxicol. Environ. Saf.171, 843–853. doi: 10.1016/j.ecoenv.2018.12.102
76
SunC.-H.LiuH.-Y.ZhangY.LuC.-H. (2019). Comparative analysis of the gut microbiota of hornbill and toucan in captivity. MicrobiologyOpen8:e00786:e786. doi: 10.1002/mbo3.786
77
SunW.ZhuB.YangF.DaiM.SeharS.PengC.et al. (2021). Optimization of biosurfactant production from pseudomonas sp. CQ2 and its application for remediation of heavy metal contaminated soil. Chemosphere265:129090. doi: 10.1016/j.chemosphere.2020.129090
78
Tajuddin SikderM.TanakaS.SaitoT.KurasakiM. (2014). Application of zerovalent iron impregnated chitosan-caboxymethyl-β-cyclodextrin composite beads as arsenic sorbent. J. Environ. Chem. Eng.2, 370–376. doi: 10.1016/j.jece.2014.01.009
79
VangayP.JohnsonA. J.WardT. L.Al-GhalithG. A.Shields-CutlerR. R.HillmannB. M.et al. (2018). US immigration westernizes the human gut microbiome. Cells175, 962–972.e10. doi: 10.1016/j.cell.2018.10.029
80
WaiteD. W.TaylorM. W. (2014). Characterizing the avian gut microbiota: membership, driving influences, and potential function. Front. Microbiol.5:223. doi: 10.3389/fmicb.2014.00223
81
WangL.WangL. A.ZhanX.HuangY.WangJ.WangX. (2020). Response mechanism of microbial community to the environmental stress caused by the different mercury concentration in soils. Ecotoxicol. Environ. Saf.188:109906. doi: 10.1016/j.ecoenv.2019.109906
82
WangW.ZhengS.SharshovK.SunH.YangF.WangX.et al. (2017). Metagenomic profiling of gut microbial communities in both wild and artificially reared bar-headed goose (Anser indicus). MicrobiologyOpen6:e00429. doi: 10.1002/mbo3.429
83
WickhamH. (2016). ggplot2: Elegant Graphics for Data Analysis. New York, NY: Springer. https://ggplot2.tidyverse.org
84
WienemannT.Schmitt-WagnerD.MeuserK.SegelbacherG.SchinkB.BruneA.et al. (2011). The bacterial microbiota in the ceca of Capercaillie (Tetrao urogallus) differs between wild and captive birds. Syst. Appl. Microbiol.34, 542–551. doi: 10.1016/j.syapm.2011.06.003
85
WuN.WangX.XuX.CaiR.XieS. (2020). Effects of heavy metals on the bioaccumulation, excretion and gut microbiome of black soldier fly larvae (Hermetia illucens). Ecotoxicol. Environ. Saf.192:110323. doi: 10.1016/j.ecoenv.2020.110323
86
WuH.WuF.-T.ZhouQ.-H.ZhaoD.-P. (2021). Comparative analysis of gut microbiota in captive and wild oriental white storks: implications for conservation biology. Front. Microbiol.12:649466. doi: 10.3389/fmicb.2021.649466
87
YoungblutN. D.ReischerG. H.WaltersW.SchusterN.WalzerC.StalderG.et al. (2019). Host diet and evolutionary history explain different aspects of gut microbiome diversity among vertebrate clades. Nat. Commun.10:2200. doi: 10.1038/s41467-019-10191-3
88
ZhangC.WangX.AshrafU.QiuB.AliS. (2017). Transfer of lead (Pb) in the soil-plant-mealybug-ladybird beetle food chain, a comparison between two host plants. Ecotoxicol. Environ. Saf.143, 289–295. doi: 10.1016/j.ecoenv.2017.05.032
89
ZhangK.WangX.GongX.SuiJ. (2022). Gut microbiome differences in rescued common kestrels (Falco tinnunculus) before and after captivity. Front. Microbiol.13:858592. doi: 10.3389/fmicb.2022.858592
90
ZhouL.HuoX.LiuB.WuH.FengJ. (2020). Comparative analysis of the gut microbial communities of the Eurasian kestrel (Falco tinnunculus) at different developmental stages. Front. Microbiol.11:592539. doi: 10.3389/fmicb.2020.592539
91
ZhouD.LiangX.WangJ.WangS.LiX.NingY. (2021). Study on the regulatory mechanism of the earthworm microbial community in vitro and in vivo under cadmium stress. Environ. Pollut.279:116891. doi: 10.1016/j.envpol.2021.116891
Summary
Keywords
gut microbiota, Leiothrix argentauris, trace element, immigrant subspecies, native subspecies
Citation
Zhou T, Liu S and Jiang A (2023) Comparison of gut microbiota between immigrant and native populations of the Silver-eared Mesia (Leiothrix argentauris) living in mining area. Front. Microbiol. 14:1076523. doi: 10.3389/fmicb.2023.1076523
Received
21 October 2022
Accepted
06 January 2023
Published
24 January 2023
Volume
14 - 2023
Edited by
Mohamed T. El-Saadony, Zagazig University, Egypt
Reviewed by
Mohsen Almutary, Imam Abdulrahman Bin Faisal University, Saudi Arabia; Ahmed M. Saad, Zagazig University, Egypt
Updates
Copyright
© 2023 Zhou, Liu and Jiang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Aiwu Jiang, aiwuu@gxu.edu.cn
This article was submitted to Microorganisms in Vertebrate Digestive Systems, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.