- 1College of Wildlife and Protected Area, Northeast Forestry University, Harbin, China
- 2Department of Genetics, College of Life Science, Northeast Forestry University, Harbin, China
Introduction: Animal intestines are extremely rich in microbial ecosystems. Numerous studies in different fields, such as epidemiology and histology, have revealed that gut microorganisms considerably mediate the survival and reproduction of animals. However, gut microbiology studies of homogeneously distributed wild cranes are still rare. This study aimed to understand the structural composition of the gut microbial community of wild cranes and elucidate the potential roles of the microorganisms.
Methods: We used high-throughput sequencing to analyze the gut microbial community structure of wild cranes in the Zhalong Nature Reserve.
Results: A total of 1,965,683 valid tags and 5248 OTUs were obtained from 32 fecal samples. Twenty-six bacteria phyla and 523 genera were annotated from the intestinal tract of the red-crowned crane. Twenty-five bacteria phyla and 625 genera were annotated from the intestine of the white-naped crane. Firmicutes, Proteobacteria, and Bacteroidetes are the dominant bacterial phyla in the intestinal tract of red-crowned cranes, while Catellicoccus, Lactobacillus, Neisseria, and Streptococcus were the dominant genera. The dominant bacterial phyla in the intestinal tract of white-naped cranes were Firmicutes, Proteobacteria, Bacteroidetes, Epsilonbacteraeota, Actinobacteria, and Fusobacteria. However, the dominant genera were Catellicoccus, Lactobacillus, Neisseria, Campylobacter, Streptococcus, Anaerobiospirillum, Romboutsia, Turicibacter, Haemophilus, and Lautropia. Firmicutes had significantly higher relative abundance in the intestine of the red-crowned than white-naped cranes (P < 0.05). However, the relative abundance of Actinobacteria and Bacteroidetes was significantly higher (P < 0.05) in the intestines of white-naped than red-crowned cranes. The diversity of the intestinal flora between the two crane species was significantly different (P < 0.05). Besides, the alpha diversity of the intestinal flora was higher for white-naped than red-crowned cranes. Eight of the 41 functional pathways differed in the gut of both crane species (P < 0.05).
Discussion: Both species live in the same area and have similar feeding and behavioral characteristics. Therefore, host differences are possibly the main factors influencing the structural and functional differences in the composition of the gut microbial community. This study provides important reference data for constructing a crane gut microbial assessment system. The findings have implications for studying deeper relationships between crane gut microbes and genetics, nutrition, immunity, and disease.
Introduction
The animal digestive tract is a complex and dynamic ecosystem with various microorganisms in the digestive tract (Hooper et al., 2001). Studies have shown that bacteria, archaea, and fungal microorganisms are the major groups of intestinal microorganisms (Muegge et al., 2011). Bacteria are the most common intestinal microorganisms (Chen et al., 2018) especially, the phyla that account for a large proportion of the total are Firmicutes, Bacteroidetes, and Proteobacteria (Sommer and Bäckhed, 2013). Numerous studies have proven that intestinal flora is crucial for host survival and reproduction (Hooper et al., 2001; McFall-Ngai et al., 2013). Intestinal flora controls energy homeostasis (Gabanyi et al., 2022), resistance disease and nutrition (Olvera-Rosales et al., 2021) and regulates the immune system (Grizotte-Lake et al., 2018).
Studies have shown that intestinal flora influences many physiological activities such as signaling, immune response, and nutrient metabolism in birds (Kohl, 2012; Grond et al., 2018). Bird gut flora can be influenced by various factors, such as habitat changes (Xie et al., 2016; Wu et al., 2018; Sun et al., 2019), weather changes (Wienemann et al., 2011; Zhang et al., 2020), host genotype (Trevelline et al., 2020; Gao et al., 2021), gender (Cui, 2019), age (Kohl et al., 2019), diet (Xie et al., 2016; Liu et al., 2020a), and health status (Ruiz-Rodríguez et al., 2009; Tong et al., 2018) among others. The dominant gut flora is key in maintaining the normal physiological functions of migratory birds (Grond et al., 2019). However, long-distance migratory birds may carry and transmit microorganisms across geographic areas (Canuti et al., 2019). Trevelline et al. (2019) and West et al. (2019) argue that understanding the microbial communities hosted by individual animals can improve conservation measures for rare wildlife and wild birds. Birds are important ecosystem indicator organisms, and the study of bird-gut microbes is ecologically important. Currently, there are approximately 10,000 bird species worldwide, nearly 2.5 times of mammals, yet research on bird gut microbes is only about 10% of the work on mammalian gut microbes (Grond et al., 2018).
Red-crowned (Grus japonensis) and white-naped (Grus vipio) cranes are large, wading, omnivorous birds. Adult cranes are monogamous and co-inhabit wetland environments. They are considered indicator species for wetland environmental quality (Zou and Wu, 2006; Wu, 2012). In February 2021, both crane species appeared in the National Class I Key Wildlife Protection of China’s National List of Key Wildlife Protection. Moreover, the IUCN Red List (2021) lists both the red-crowned crane and the white-naped crane as vulnerable (VU). The International Crane Foundation (ICF) estimates that 2,800–3,430 red-crowned cranes are in the wild globally. Furthermore, two major populations of red-crowned cranes exist worldwide: island resident populations and continental migratory populations. For white-naped cranes, the worldwide population in the wild is between 6,700–7,700 individuals. The western population winters in the middle and lower reaches of the Yangtze River in China, and the eastern population winters on the Korean Peninsula and Kyushu Island, Japan. Unfortunately, the populations of both crane species are declining in China.
Between March to November, wild red-crowned and white-naped cranes breed in some wetlands of Heilongjiang Province (Zou et al., 2018). The species of plants eaten by the cranes are similar during this breeding period, but the proportions vary. At courtship, the overlap of nutritional and ecological niches between the two populations was 91.67% (Wu, 2012). Although the two crane species are co-dominantly distributed, they differ in their selection and use of the microhabitats (Li, 2017).
Zhalong National Nature Reserve, located in the middle of the East Asia-Australasia bird migration corridor, has the largest reed wetland in Asia (Yang, 2015). The reserve acceded to the Convention on Wetlands of International Importance in 1992 as a waterfowl habitat. Zhalong National Nature Reserve focuses on protecting rare birds and wetland ecosystems, among them are birds including red-crowned and white-naped cranes, the oriental white stork (Ciconia boyciana), Eurasian Spoonbill (Platalea leucorodia), bustard (Otis tarda), and reed parrotbill (Paradoxornis heudei). The Zhalong Reserve is rich in crane resources and is an essential global crane research hotspot. Since the two crane species have similar breeding rhythms and reproductive patterns (Wu, 2012; Hao, 2013; Li, 2017), they can represent the interspecific differences in wild bird intestinal flora.
Therefore, this study used non-invasive sampling methods to collect feces from wild red-crowned and white-naped cranes of the same range. The 16S rRNA high-throughput sequencing technology was applied to analyze the composition and structure of crane intestinal flora. The study hypothesized that: (1) The gut microbial communities of cranes from the same habitats domain have similar composition. (2) The crane intestinal microbial communities perform similar functions. (3) The gut microbial communities of cranes from the same species are more similar than between individuals of different species.
Materials and methods
Ethics statement
Fourteen and eighteen fecal samples of wild red-crowned and white-naped cranes were collected from the natural environment using non-contact and non-invasive sampling methods without damaging cranes. The administration of Zhalong National Nature Reserve approved the sample collection.
Sample collection
Test samples were collected from the Zhalong Nature Reserve (46°52′–47°32′N, 123°47′–124°37′E) between March and April 2021. The average annual precipitation and annual temperature in this area are 416 mm and 3.2°C, respectively. The climate during the sampling period is characterized by low precipitation, rapid warming, sandy winds, and partial dryness. Red-crowned cranes and white-naped cranes mainly eats maize (Zea mays), reeds (Phragmites australis), and Suaeda glauca in this period at Zhalong Nature Reserve (Wu, 2012). A preliminary survey revealed red-crowned, and white-naped cranes forage and rest in the farmland near Kachin Lake located within the reserve. Thus, that farmland was selected as the sampling site for this study (Figure 1).
To ensure proper species identification of each sample, we used telescopes to search for red-crowned crane and white-naped crane flocks, waited until the birds flew away and then collected fresh fecal samples. The interval distance for samples was more than 5 m to avoid individual repetition. The fresh and moist feces were kept in a portable liquid nitrogen tank and transported refrigerated to the lab as quickly as possible. The outside of each sample was cut and discarded to avoid contamination; the rest was stored at –20°C, transported to the laboratory, and stored at –80°C until further analysis for DNA extraction that could ensure the success and accuracy of subsequent sequencing (Supplementary Figure 1).
Bacterial DNA extraction and sequencing
Total bacterial DNA was isolated from fecal samples using the MagPure Soil DNA LQ Kit (Producers: Magen, Cat. No.: D6356-02) following the manufacturer’s instructions. DNA integrity was measured using agarose gel electrophoresis. Next, PCR amplification of the V3-V4 hypervariable regions of the bacterial 16S rRNA gene was performed in a 25 μL reaction using universal primer pairs (343F: 5’-TACGGRAGGCAGCAG-3’; 798R: 5’-AGGGTATCTAATCCT-3’). The PCR products were purified with the Agencourt AMPure XP beads (Beckman Coulter Inc., CA, USA) and quantified using the Qubit dsDNA assay kit (Thermo Fisher Scientific, MA, USA). The library concentrations were adjusted accordingly, and libraries were sequenced on an Illumina NovaSeq6000 on two paired-end read cycles of 250 bases (Illumina Inc., CA, USA; OE Biotech Company; Shanghai, China). The raw data were submitted to the Sequence Read Archive at the NCBI database under BioProject accession number PRJNA835714.
Bioinformatic analysis
Paired-end reads were assembled using the FLASH software (Reyon et al., 2012). Sequences with ambiguous and homologous sequences or with < 200 bp were discarded. The QIIME software (version 1.8.0) (Caporaso et al., 2010) retained reads where ≥ 75% of the bases scored Q20. Then, reads with chimeras were removed using VSEARCH (Rognes et al., 2016). Clean reads were subjected to primer sequences removal and clustering to generate operational taxonomic units (OTUs) using the VSEARCH software with a 97% similarity cutoff (Rognes et al., 2016). A representative read was selected from each OTU using the QIIME package. All representative reads were annotated and blasted against the Silva database (Version 132) using the RDP classifier (Wang et al., 2007).
Statistical analysis
An alpha diversity analysis reflected the degree of bacterial diversity within the intestinal environment of the two crane species. The Chao1 (Chao, 1984), Shannon Wiener (Hill et al., 2003), Simpson (Simpson, 1949), and Good’s Coverage (Esty, 1986) indexes were calculated for the two types of crane intestinal flora based on uniform depth. The Wilcoxon sum rank test determined the difference in the alpha diversity indexes between the intestinal flora of the two crane species. Therefore, the number of actual observed OTUs (Observed Species) was used to plot alpha diversity dilution curves. Adonis analysis (i.e., PERMANOVA analysis) was used to test the difference between the different subgroups (P < 0.05 indicates a significant difference between groups, and the R2 value indicates the degree of the explained difference). A PCoA (Principal coordinates analysis) analysis based on a weighted Unifrac distance algorithm (considering the evolutionary relationships and species abundance) was used to compare the beta diversity between the gut of the red-crowned and white-naped cranes. Furthermore, the LEfSe (Linear discriminant analysis coupled with effect size measurements) analysis revealed the intestinal flora composition of the two crane species. All the statistical analyses were completed using the R package.
The functional prediction analysis of gut microbes was based on the 16S sequence data annotated on the Greengenes database (DeSantis et al., 2006). The PICRUSt (Phylogenetic investigation of communities by reconstruction of unobserved states) software predicted the gene functions of known microbial (Langille et al., 2013). The results were compared with the KEGG (Kyoto Encyclopedia of Genes and Genomes) database to find homologous genes and gene copy numbers. Thus, the corresponding KEGG Ortholog (KO) information was used for predicting functions. The abundance of each functional category was calculated based on the KO information and the corresponding OTU abundance. Finally, the Wilcoxon sum rank test established the functional differences between the gut microorganisms of the two crane species.
Results
Sequencing statistics
This study used the Illumina Novoseq platform to amplify and detect 16S rRNA sequences from the fecal microbiota of 14 red-crowned and 18 white-naped cranes. The 32 stool samples yielded 1,965,683 valid tags, including 43,917–69,817 valid tags per sample measuring 404.12–425.34 bp long. The sequences with ≥ 97% QC were classified as one OTU; thus, 5248 OTUs were obtained (651–2,130 OTUs in each sample). The observed OTUs of randomly selected sequences generated the rarefaction curves for the samples. Moreover, increasing the sequencing depth (the number of sampled sequences) increased the rarefaction curves of most samples, indicating a non-significant increase in new OTUs over the increasing sequencing depth (Supplementary Figure 2). The Good’s Coverage Index showed that the sequencing depth covered > 99% of the species in the sample, indicating that the experimental sampling was adequate.
Comparison of the structure of the intestinal flora between the red-crowned and white-naped cranes
The red-crowned and the white-naped crane share 3088 OTUs, accounting for 58.84% of all the sampled OTUs. The red-crowned crane had 10.83% of the bacterial OTUs (375), and the white-naped crane had 57.80% of the endemic OTUs (1785) (Supplementary Figure 3). Simultaneously, we identified 72 OTUs as the core gut microbiota (the taxa) shared by all the individuals in each group of red-crowned and white-naped cranes (Supplementary Figure 4). The shared taxa were annotated to Firmicutes, Bacteroidetes, Proteobacteria, Actinobacteria, Fusobacteria, Epsilonbacteraeota, and Gemmatimonadetes.
Taxonomically, the red-crowned and white-naped cranes had 26 and 25 annotated intestinal phyla, 54 and 61 classes, 136 and 159 bacterial orders, 248 and 275 families, and 523 and 625 bacterial genera. We plotted the relative abundance bar plot of the intestinal flora of both cranes. The phylum-level flora in the intestine of red-crowned cranes consisted mainly of Firmicutes (87.55%), Proteobacteria (7.65%), and Bacteroidetes (2.14%). In contrast, the phylum-level flora in the intestine of white-naped cranes consisted mainly of Firmicutes (66.91%), Proteobacteria (18.83%), Bacteroidetes (4.80%), Epsilonbacteraeota (4.30%), Actinobacteria (3.30%), and Fusobacteria (1.01%) (Figure 2A). Besides, the genus-level flora in the intestine of red-crowned cranes consisted mainly of Catellicoccus (76.10%), Lactobacillus (7.87%), Neisseria (2.75%), and Streptococcus (1.06%). Catellicoccus (43.39%), Lactobacillus (11.05%), Neisseria (8.09%), Campylobacter (4.25%), Streptococcus (3.25%), Anaerobiospirillum (3.03%), Romboutsia (2.04%), Turicibacter (1.49%), Haemophilus (1.40%), and Lautropia (1.08%) dominated the genus-level flora in the intestine of the white-naped crane (Supplementary Figure 5).
Figure 2. (A) Relative abundance bar plot of the relative abundance of the strains. Relative abundance (%) of the two cranes Top10 phyla. Red-crowned crane, A1–A14; white-naped crane, B1–B18. (B) Annotated branching diagram of differential strains. Red nodes indicate species with relatively high abundance in the red group, green nodes indicate species with relatively high abundance in the green group, yellow nodes indicate species with no significant difference in the comparison between the two groups, the node diameter size is proportional to the relative abundance size, each layer of nodes indicates the phylum/class/order/family/genus from the inside out. (A) Red-crowned crane, (B) white-naped crane.
The LEfSe analysis revealed the composition of the strains that differed between the groups and those that contributed more to the observed differences in each group (Figure 2B). In the gut of red-crowned cranes, Firmicutes (P < 0.05) were differentially dominant species at the phylum level, Catellicoccus (P < 0.01) and Rahnella (P < 0.01) at the genus level. However, the differentially dominant species in the gut of white-naped cranes were Actinobacteria (P < 0.05) and Bacteroidetes (P < 0.05) at the phylum level, Clostridium_sensu_ stricto_1 (P < 0.05), Romboutsia (P < 0.001), Lachnospiraceae_NK4A136_group (P < 0.01), Allorhizobium_Neorhizobium_Pararhizobium_Rhizobium (P < 0.001), Actinoplanes (P < 0.05), Terrisporobacter (P < 0.05), Alloprevotella (P < 0.001), and Turicibacter (P < 0.01) at the genus level.
Analysis of the diversity of the intestinal flora of two species of cranes
The Chao index (P < 0.01), Shannon (P < 0.01), and Simpson indexes (P < 0.01) and observed Species (P < 0.05) showed that the gut flora of white-naped cranes was significantly more diverse than that of the red-crowned crane (Figure 3). Moreover, the Adonis analysis revealed significant differences between the intestinal flora of the two cranes species (P < 0.05, R2 = 0.139). PCoA analysis showed that individuals of the same species clustered together, indicating that both species had their own unique gut microbiota composition. Interestingly, most white-naped cranes were more dispersed, suggesting more specificity in the intestinal flora of white-naped cranes (Figure 4). The clustering of some white-naped cranes showed similarity to the red-crowned crane population, possibly due to the presence of both cranes foraging in the same area.
Figure 3. Alpha diversity index of microbial communities in violin map. (A) Representative chao1 Index, (B) representative observed species, (C) representative Shannon Index, (D) representative Simpson index. *Correlation significant at 0.05 level. **Correlation significant at 0.01 level.
Figure 4. PCoA analysis graph of the colony. The same color in the graph is the same grouping. PC1 and PC2 are the first and second principal coordinates with the greatest explanation of the difference between samples.
Functional prediction analysis of the gut microbial community of the two crane species
We predicted functions of gut microorganisms in red-crowned and white-naped cranes by PICRUSt analysis and annotated 41 functional pathways at KEGG level2. The KEGG pathways for membrane transport (16.19%), carbohydrate metabolism (12.92%), replication and repair (8.62%), amino acid metabolism (8.09%), translation (5.77%), and poorly characterized (5.15%) were highly abundant in intestinal microorganisms from red-crowned cranes. In contrast, the gut microbes from white-naped cranes were highly abundant in KEGG pathways for membrane transport (14.15%), carbohydrate metabolism (11.36%), amino acid metabolism (9.10%), replication and repair (8.56%), translation (5.60%), energy metabolism (5.29%), and poorly characterized (5.11%).
Among the 41 functional pathways, cardiovascular diseases (P < 0.01) and cell motility (P < 0.01), excretory system (P < 0.05), immune system (P < 0.05), endocrine system (P < 0.05), cancers (P < 0.05), circulatory system (P < 0.05), and neurodegenerative diseases (P < 0.05) were significantly different between the two crane species (Figure 5). These eight differentially functional pathways were highly abundant in white-naped than red-crowned cranes.
Figure 5. Heat map of differential functional clustering. Horizontal is the sample number; vertical is the functional annotation information; the clustering tree on the left side of the figure is the functional clustering tree; the upper clustering branches represent samples from different subgroups. The red color indicates a higher relative abundance of species, the blue color indicates a lower relative abundance of species. A, Red-crowned crane; B, white-naped crane.
Discussion
The Zhalong Nature Reserve is among the most important breeding sites for wild and red-crowned cranes in China (Hao, 2013; Li, 2017). For the first time, we used high-throughput sequencing to analyze the intestinal flora of red-crowned and white-naped cranes in the wild. This study annotated 26 phyla and 523 genera in the gut of the red-crowned crane and 25 phyla and 625 genera in the gut of the white-naped cranes, which is significantly greater than previous results on the same species (Trevelline et al., 2020). This difference may be because the cranes in this study were wild and lived in a complex, variable environment. In contrast, the cranes in the Trevelline et al. (2020) study were caged and lived in a more homogeneous environment.
Firmicutes had the highest relative abundance in the gut of both crane species at the phylum level. This observation is consistent with that from the gut of 13 other cranes and many waterfowl species such as egrets (Egretta garzetta), black-crowned night heron (Nycticorax nycticorax), cormorant (Phalacrocorax carbo), and black swans (Cygnus atratus). Numerous studies have proven that Firmicutes are the most abundant in the intestines of almost all vertebrates, not only birds (Ding et al., 2017). This phenomenon may be greatly related to the function of Firmicutes, which degrade sugars and fats (Bäckhed et al., 2004) in intestinal food to produce energy and other nutrients for direct use by the host. Thus, they improve food utilization by the host (Flint et al., 2008). In addition, Firmicutes metabolize carbohydrates to produce short-chain fatty acids (SCFAs), which are also closely linked to host immune regulation (Atarashi et al., 2011; Zhang et al., 2015). Proteobacteria and Bacteroidetes were also common dominant groups in the intestines of both cranes, consistent with the dominant phylum in the intestines of chickens (Gallus gallus domesticus), turkeys (Meleagris gallopavo) (Wei et al., 2013), and 15 cranes species (Trevelline et al., 2020). The gut of swans also has similar dominant phyla composition (Wang et al., 2021). Proteobacteria are a very complex bacterial phylum that utilizes carbon-derived substances and is important in energy accumulation (Lu et al., 2012). Moreover, Proteobacteria have various physiological functions that are important for maintaining the homeostasis of their intestinal microbial communities and the healthy growth and development of their hosts (Tenaillon et al., 2010; Samanta et al., 2012). Bacteroidetes mainly ferment carbohydrates, polysaccharides, proteins, and bile acid, perform steroid metabolism, and promote fat accumulation (Salyers, 1984). This phylum produces butyric acid (Pryde et al., 2002), develops the host immune system, enhances host immunity (Hooper et al., 2001; Hooper, 2004), and are important for balancing the intestinal flora (Sears, 2005). In the spring, both red-crowned and white-naped cranes prefer foods with higher crude protein and crude fat and lower crude fiber (Wu, 2012). Moreover, accumulating Firmicutes, Proteobacteria, and Bacteroidetes benefits both cranes by obtaining sufficient energy and nutrients from these foods. These phyla also facilitate quick recovery of the copious energy consumed during the spring migration and prepare birds for the arrival of the breeding season. Genus Catellicoccus had the highest relative abundance in the gut of both crane species, consistent with the gut of the black-headed gull (Chroicocephalus ridibundus) and the whooper swan (Cygnus cygnus) in the Rongcheng area of Shandong Province (Wang et al., 2021). In contrast, the results differed from those of the cormorant, egret, and black-crowned night heron (Laviad-Shitrit et al., 2019). Catellicoccus is present in most waterfowl and can be used as indicator species for fecal contamination in water or environments contaminated by waterfowl feces (Koskey et al., 2014; Byappanahalli et al., 2015; Safaie et al., 2021). Besides, Catellicoccus also function in nutrient transport and bile acid hydrolysis (Góngora et al., 2021). Lactobacillus was highly abundant in the gut microbial communities of both cranes. Similarly, Lactobacillus was highly abundant in the intestines of other birds, including G. Canadensis (Sandhill Crane) (Trevelline et al., 2020), black-headed gulls (Laviad-Shitrit et al., 2019), black-necked cranes (Grus nigricollis) (Wang et al., 2020), and New Zealand parrots (Strigops habroptilus) (Waite et al., 2012). Moreover, Lactobacillus is a beneficial bacterium in the intestine as it synthesizes vitamins and amino acids required by the animal, promotes mineral absorption and regulates intestinal pH value. Lactobacillus inhibits the growth of pathogenic microorganisms, enhances the barrier function of intestinal epithelial cells, reduces inflammation, and enhances the body’s immune capacity (Artis, 2008; Milani et al., 2017).
In the intestinal flora of the two crane species, a LEfSe analysis revealed Firmicutes as the dominant phylum and Catellicoccus as the dominant genus. The differentially dominant phylum in the intestine of the white-naped crane were Actinobacteria and Bacteroidetes. Normally, the intestinal ratio of Firmicutes to Bacteroidetes (F/B) is related to the strength of fat metabolism in the animal body (Ley et al., 2005). In this study, the intestinal F/B in the red-crowned cranes was higher than that of the white-naped cranes, possibly because the red-crowned cranes have a larger body size compared to the white-naped cranes. Furthermore, the red-crowned cranes create larger nests and lay larger eggs during the breeding season (Hao, 2013; Li, 2017). All these differences normally require more energy consumption. Higher F/B also indicates that the cranes have higher fat metabolism and can obtain higher energy to sustain greater energy expenditure. Additionally, Actinobacteria are a dominant flora in the environment (Wink et al., 2017). The significant abundance of actinobacteria in the gut of the white-naped compared to the red-crowned crane may be due to the close association between the intestinal flora of the white-naped cranes and the environment.
The diversity analysis also showed that the alpha diversity of gut flora was significantly lower in red-crowned than white-naped cranes. Moreover, the Adonis and PCoA analyses showed significantly different structures of the intestinal flora of the two cranes (P < 0.05, R2 = 0.139). The intestinal flora diversity of black-necked cranes is significantly different between different regions (Wang et al., 2020), similar to the lesser white-fronted goose (Anser erythropus) (Liu et al., 2020a). During different migration seasons, the gut flora of hooded cranes (Grus monacha) was significantly different (Zhang et al., 2020). We selected wild red-crowned and white-naped cranes from the same distribution area in the spring. Thus, we excluded the influence of spatial and temporal factors on the intestinal flora diversity of the wild cranes. Studies on the lesser white-fronted goose have shown that food habits may influence intestinal flora diversity (Liu et al., 2020a). Wild red-crowned and white-naped cranes in the Zhalong Reserve feed on similar spring food types, and their trophic ecological niche overlap is approximately 90% (Wu, 2012). Therefore, food factors were not considered the main factors influencing the differences in the intestinal flora between the two crane species. Furthermore, Figure 4 showed that white-naped cranes appear to be separated into two clusters, one similar to red-crowned cranes and another separate cluster along the PC1. A possible explanation for this might be that the two species sharing the same foraging ground. According to our observation, white-naped cranes did not migrate to Zhalong wetland at the same time, some of that migrated at different times would forage in groups before the breeding season. Therefore, some individuals of the white-naped crane appear to be clustered and scattered. This result may be explained by the fact that individuals from different small migratory populations arrive at the Zhalong on different dates, while the intestinal flora of the white-naped cranes that arrived later may still retain the environmental adaptation characteristics of some migratory transit sites (Góngora et al., 2021).
By examining the intestinal flora of three species of geese, Yang et al. (2016) suggested that host species may influence the composition of intestinal flora. The gut microbial profile of 59 tropical bird species showed that host species had a greater effect on the gut flora than the host diet and inhabited region (Hird et al., 2015). Both results are consistent with our view that even with closely related species living in the same environment, the gut microbial community structure develops unique characteristics. The top five functional pathways annotated at level 2 contain functional genes that enrich membrane transport, carbohydrate metabolism, amino acid metabolism, and replication and repair. The results generally agree with the predicted gut microbial function of geese and red-crowned cranes in the middle and lower reaches of the Yangtze River (Yang et al., 2016; Dong et al., 2019). Moreover, the results are very similar to those of captive-bred chickens (Zhang, 2021). The predicted gut microbial functions of yak (Bos mutus) and donkey (Equus asinus) are also consistent with the current study (Liu et al., 2019). However, the predicted results of gut microbial functions of the Siberian tiger (Panthera tigris altaica) and tiger frog (Hoplobatrachus rugulosus) differed greatly from those of the crane (Liu et al., 2020b; Ning et al., 2020). Accordingly, we hypothesized that the closer the genetic distance, the more similar the gut microbial function in animals with similar food habits. Besides, the function of gut microbes is related to the adaptation of the host to the external environment (Newell and Douglas, 2014). For example, an enhanced corn (Zea mays L.) diet may promote the gut microbial metabolism of carbohydrates and lipids (Zheng et al., 2019).
In this study, wild red-crowned and white-naped cranes spent a lot of time foraging in the corn field. The high abundance of genes for gut microbial metabolic functions allows both cranes to efficiently obtain numerous nutrients from food to improve their adaptability to food and the environment. The abundance of the main functional groups of microorganisms in the gut of the two cranes was not significantly different (P > 0.05). However, the abundances of the dominant phylum and genera common to the gut of the two crane species were significantly different (P < 0.05). This difference suggests that the gut of different animals with the same domain distribution may be adapted to their respective habits by enriching different species of microorganisms with similar functions.
Conclusion
This study analyzed the gut microbial communities of the red-crowned and white-naped cranes using 16S rRNA high-throughput sequencing. The results showed significant differences in the composition of gut microbial communities between the two crane species in the same range. However, the major functional groups were not significantly different. This study supports the view that host differences are the main reason for the different gut microbial communities. This data will help to develop more scientific conservation plans for cranes.
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 below: https://www.ncbi.nlm.nih.gov/, PRJNA835714.
Ethics statement
This animal study was reviewed and approved by Animal Welfare and Animal Ethics Committee of Northeast Forestry University.
Author contributions
ZG, YW, and HZ designed the experiments. ZG, HS, ZL, and YT completed the field sampling. HD, ZG, and XJ performed the data analysis and prepared the figures. ZG wrote the manuscript. YT and ZL contributed to the revision of manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by the National Nature Science Foundation of China (No. 31770423), the Fundamental Research Funds for the Central Universities (No. 2572020DR11), the Natural Sciences Foundation of Heilongjiang Province (No. LH2021C004), the Heilongjiang Provincial Scientific Research Institutes Scientific Research Business Fund (No. CZKYF2022-1-C022), and the Outstanding Youth Fund of Heilongjiang Academy of Sciences (No. CXJQ2023ZR01).
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.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1007884/full#supplementary-material
References
Artis, D. (2008). Epithelial-cell recognition of commensal bacteria and maintenance of immune homeostasis in the gut. Nat. Rev. Immunol. 8, 411–420. doi: 10.1038/nri2316
Atarashi, K., Tanoue, T., Shima, T., Imaoka, A., Kuwahara, T., Momose, Y., et al. (2011). Induction of colonic regulatory T cells by indigenous Clostridium species. Science 331, 337–341. doi: 10.1126/science.1198469
Bäckhed, F., Ding, H., Wang, T., Hooper, L. V., Koh, G. Y., Nagy, A., et al. (2004). The gut microbiota as an environmental factor that regulates fat storage. Proc. Natl. Acad. Sci. 101, 15718–15723. doi: 10.1073/pnas.0407076101
Byappanahalli, M. N., Nevers, M. B., Whitman, R. L., Ge, Z., Shively, D., Spoljaric, A., et al. (2015). Wildlife, urban inputs, and landscape configuration are responsible for degraded swimming water quality at an embayed beach. J. Great Lakes Res. 41, 156–163. doi: 10.1016/j.jglr.2014.11.027
Canuti, M., Munro, H. J., Robertson, G. J., Kroyer, A. N., Roul, S., Ojkic, D., et al. (2019). New insight into avian papillomavirus ecology and evolution from characterization of novel wild bird papillomaviruses. Front. Microbiol. 10:701. doi: 10.3389/fmicb.2019.00701
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Chao, A. (1984). Non-parametric estimation of the classes in a population. Scand. J. Statist. 11, 265–270. doi: 10.2307/4615964
Chen, S., Si, H., Li, G., and Liu, H. (2018). Research progress of animal gut microbiota and nutrient metabolism. Feed Industry 39, 33–36. doi: 10.13302/j.cnki.fi.2018.02.006
Cui, L. (2019). Comparative Study on the Gut Microbiota Differences of Chickens with Different Gender. China: Huazhong Agricultural University.
DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05
Ding, Y., Wu, Q., Hu, Y., Wang, X., Nie, Y., Wu, X., et al. (2017). Advances and prospects of gut microbiome in wild mammals. Acta Theriol. Sinica 37, 399–406. doi: 10.16829/j.slxb.201704010
Dong, Y., Xiang, X., Zhao, G., Song, Y., and Zhou, L. (2019). Variations in gut bacterial communities of hooded crane (Grus monacha) over spatial-temporal scales. PeerJ 7:e7045. doi: 10.7717/peerj.7045
Esty, W. W. (1986). The efficiency of Good’s nonparametric coverage estimator. Ann. Statist. 14, 1257–1260. doi: 10.1214/AOS/1176350066
Flint, H. J., Bayer, E. A., Rincon, M. T., Lamed, R., and White, B. A. (2008). Polysaccharide utilization by gut bacteria: potential for new insights from genomic analysis. Nat. Rev. Microbiol. 6, 121–131. doi: 10.1038/nrmicro1817
Gabanyi, I., Lepousez, G., Wheeler, R., Vieites-Prado, A., Nissant, A., Wagner, S., et al. (2022). Bacterial sensing via neuronal Nod2 regulates appetite and body temperature. Science 376:eabj3986. doi: 10.1126/science.abj3986
Gao, L., Liu, L., Du, C., and Hou, Q. (2021). Comparative analysis of fecal bacterial microbiota of six bird species. Front. Vet. Sci. 8:791287. doi: 10.3389/fvets.2021.791287
Góngora, E., Elliott, K. H., and Whyte, L. (2021). Gut microbiome is affected by inter-sexual and inter-seasonal variation in diet for thick-billed murres (Uria lomvia). Sci. Rep. 11:1200. doi: 10.1038/s41598-020-80557-x
Grizotte-Lake, M., Zhong, G., Duncan, K., Kirkwood, J., Iyer, N., Smolenski, I., et al. (2018). Commensals suppress intestinal epithelial cell retinoic acid synthesis to regulate interleukin-22 activity and prevent microbial dysbiosis. Immunity 49, 1103–1115.e6. doi: 10.1016/j.immuni.2018.11.018.
Grond, K., Sandercock, B. K., Jumpponen, A., and Zeglin, L. H. (2018). The avian gut microbiota: community, physiology and function in wild birds. J. Avian Biol. 49:e01788. doi: 10.1111/jav.01788
Grond, K., Santo Domingo, J. W., Lanctot, R. B., Jumpponen, A., Bentzen, R. L., Boldenow, M. L., et al. (2019). Composition and drivers of gut microbial communities in arctic-breeding shorebirds. Front. Microbiol. 10:2258. doi: 10.3389/fmicb.2019.02258
Hao, M. (2013). Comparative Research on Behaviors between the Red-Crowned Crane and White-Naped Crane During the Breeding Season in Zhalong National Nature Reserve. Harbin: Northeast Forestry University.
Hill, T. C., Walsh, K. A., Harris, J. A., and Moffett, B. F. (2003). Using ecological diversity measures with bacterial communities. FEMS Microbiol. Ecol. 43, 1–11. doi: 10.1111/j.1574-6941
Hird, S. M., Sánchez, C., Carstens, B. C., and Brumfield, R. T. (2015). Comparative gut microbiota of 59 neotropical bird species. Front. Microbiol. 6:1403. doi: 10.3389/fmicb.2015.01403
Hooper, L. V. (2004). Bacterial contributions to mammalian gut development. Trends Microbiol. 12, 129–134. doi: 10.1016/j.tim.2004.01.001
Hooper, L. V., Wong, M. H., Thelin, A., Hansson, L., Falk, P. G., and Gordon, J. I. (2001). Molecular analysis of commensal host-microbial relationships in the intestine. Science 291, 881–884. doi: 10.1126/science.291.5505.881
Kohl, K. D. (2012). Diversity and function of the avian gut microbiota. J. Comp. Physiol. B 182, 591–602. doi: 10.1007/s00360-012-0645-z
Kohl, K. D., Brun, A., Caviedes-Vidal, E., and Karasov, W. H. (2019). Age-related changes in the gut microbiota of wild House Sparrow nestlings. IBIS 161, 184–191. doi: 10.1111/ibi.12618
Koskey, A. M., Fisher, J. C., Traudt, M. F., Newton, R. J., and McLellan, S. L. (2014). Analysis of the gull fecal microbial community reveals the dominance of Catellicoccus marimammalium in relation to culturable Enterococci. Appl. Environ. Microbiol. 80, 757–765. doi: 10.1128/AEM.02414-13
Langille, M. G., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676
Laviad-Shitrit, S., Izhaki, I., Lalzar, M., and Halpern, M. (2019). Comparative analysis of intestine microbiota of four wild waterbird species. Front. Microbiol. 10:1911. doi: 10.3389/fmicb.2019.01911
Ley, R. E., Bäckhed, F., Turnbaugh, P., Lozupone, C. A., Knight, R. D., and Gordon, J. I. (2005). Obesity alters gut microbial ecology. Proc. Natl. Acad. Sci. U S A. 102, 11070–11075. doi: 10.1073/pnas.0504978102
Li, Q. (2017). Comparison of Habitat Use between Red-crowned Crane and White-naped Crane in Songnen Plain. Harbin: Northeast Forestry University.
Liu, G., Bou, G., Su, S., Xing, J., Qu, H., Zhang, X., et al. (2019). Microbial diversity within the digestive tract contents of Dezhou donkeys. PLoS One 14:e0226186. doi: 10.1371/journal.pone.0226186
Liu, G., Gong, Z., and Li, Q. (2020a). Variations in gut bacterial communities between lesser white-fronted geese wintering at Caizi and Shengjin lakes in China. Microbiologyopen 9:e1037. doi: 10.1002/mbo3.1037
Liu, W., Chen, Y., Guo, X., Tang, Y., Xiao, H., and Li, S. (2020b). Analysis of microbial composition and function in the intestine of farmed tiger frog(Rana rugulosa). Chinese J. Wildlife 41, 171–181.
Lu, H.-P., Wang, Y.-B., Huang, S.-W., Lin, C.-Y., Wu, M., Hsieh, C.-H., et al. (2012). Metagenomic analysis reveals a functional signature for biomass degradation by cecal microbiota in the leaf-eating flying squirrel (Petaurista alborufus lena). BMC Genom. 13:466. doi: 10.1186/1471-2164-13-466
McFall-Ngai, M., Hadfield, M. G., Bosch, T. C., Carey, H. V., Domazet-Lošo, T., Douglas, A. E., et al. (2013). Animals in a bacterial world, a new imperative for the life sciences. Proc. Natl. Acad. Sci. U S A. 110, 3229–3236. doi: 10.1073/pnas.1218525110
Milani, C., Duranti, S., Bottacini, F., Casey, E., Turroni, F., Mahony, J., et al. (2017). The first microbial colonizers of the human gut: composition, activities, and health implications of the infant gut microbiota. Microbiol. Mol. Biol. Rev. 81:e00036-17. doi: 10.1128/MMBR.00036-17
Muegge, B. D., Kuczynski, J., Knights, D., Clemente, J. C., González, A., Fontana, L., et al. (2011). Diet drives convergence in gut microbiome functions across mammalian phylogeny and within humans. Science 332, 970–974. doi: 10.1126/science.1198719
Newell, P. D., and Douglas, A. E. (2014). Interspecies interactions determine the impact of the gut microbiota on nutrient allocation in Drosophila melanogaster. Appl. Environ. Microbiol. 80, 788–796. doi: 10.1128/aem.02742-13
Ning, Y., Qi, J., Dobbins, M. T., Liang, X., Wang, J., Chen, S., et al. (2020). Comparative analysis of microbial community structure and function in the gut of wild and captive Amur tiger. Front. Microbiol. 11:1665. doi: 10.3389/fmicb.2020.01665
Olvera-Rosales, L.-B., Cruz-Guerrero, A.-E., Ramírez-Moreno, E., Quintero-Lira, A., Contreras-López, E., Jaimez-Ordaz, J., et al. (2021). Impact of the gut microbiota balance on the health–disease relationship: the importance of consuming probiotics and prebiotics. Foods 10:1261. doi: 10.3390/foods10061261
Pryde, S. E., Duncan, S. H., Hold, G. L., Stewart, C. S., and Flint, H. J. (2002). The microbiology of butyrate formation in the human colon. FEMS Microbiol. Lett. 217, 133–139. doi: 10.1111/j.1574-6968
Reyon, D., Tsai, S. Q., Khayter, C., Foden, J. A., Sander, J. D., and Joung, J. K. (2012). FLASH assembly of TALENs for high-throughput genome editing. Nat. Biotechnol. 30, 460–465. doi: 10.1038/nbt.2170
Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahé, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi: 10.7717/peerj.2584
Ruiz-Rodríguez, M., Soler, J. J., Lucas, F. S., Heeb, P., Jose Palacios, M., Martín-Gálvez, D., et al. (2009). Bacterial diversity at the cloaca relates to an immune response in magpie Pica pica and to body condition of great spotted cuckoo Clamator glandarius nestlings. J. Avian Biol. 40, 42–48. doi: 10.1111/J.1600-048X
Safaie, A., Weiskerger, C. J., Nevers, M. B., Byappanahalli, M. N., and Phanikumar, M. S. (2021). Evaluating the impacts of foreshore sand and birds on microbiological contamination at a freshwater beach. Water Res. 190:116671. doi: 10.1016/j.watres.2020.116671
Salyers, A. (1984). Bacteroides of the human lower intestinal tract. Ann. Rev. Microbiol. 38, 293–313. doi: 10.1146/annurev.mi.38.100184.001453
Samanta, A. K., Torok, V. A., Percy, N. J., Abimosleh, S. M., and Howarth, G. S. (2012). Microbial fingerprinting detects unique bacterial communities in the faecal microbiota of rats with experimentally-induced colitis. J. Microbiol. 50, 218–225. doi: 10.1007/s12275-012-1362-8
Sears, C. L. (2005). A dynamic partnership: celebrating our gut flora. Anaerobe 11, 247–251. doi: 10.1016/j.anaerobe.2005.05.001
Sommer, F., and Bäckhed, F. (2013). The gut microbiota—masters of host development and physiology. Nat. Rev. Microbiol. 11, 227–238. doi: 10.1038/nrmicro2974
Sun, C. H., Liu, H. Y., Zhang, Y., and Lu, C. H. (2019). Comparative analysis of the gut microbiota of hornbill and toucan in captivity. Microbiologyopen 8:e00786. doi: 10.1002/mbo3.786
Tenaillon, O., Skurnik, D., Picard, B., and Denamur, E. (2010). The population genetics of commensal Escherichia coli. Nat. Rev. Microbiol. 8, 207–217. doi: 10.1038/nrmicro2298
Tong, X., Rehman, M. U., Huang, S., Jiang, X., Zhang, H., and Li, J. (2018). Comparative analysis of gut microbial community in healthy and tibial dyschondroplasia affected chickens by high throughput sequencing. Microbial Pathogenesis 118, 133–139. doi: 10.1016/j.micpath.2018.03.001
Trevelline, B. K., Fontaine, S. S., Hartup, B. K., and Kohl, K. D. (2019). Conservation biology needs a microbial renaissance: a call for the consideration of host-associated microbiota in wildlife management practices. Proc. R. Soc. B 286:20182448. doi: 10.1098/rspb.2018.2448
Trevelline, B. K., Sosa, J., Hartup, B. K., and Kohl, K. D. (2020). A bird’s-eye view of phylosymbiosis: weak signatures of phylosymbiosis among all 15 species of cranes. Proc. R. Soc. B 287:20192988. doi: 10.1098/rspb.2019.2988
Waite, D. W., Deines, P., and Taylor, M. W. (2012). Gut microbiome of the critically endangered New Zealand parrot, the kakapo (Strigops habroptilus). PLoS One 7:e35803. doi: 10.1371/journal.pone.0035803
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, W., Huang, S., Yang, L., and Zhang, G. (2021). Comparative analysis of the fecal bacterial microbiota of wintering whooper swans (Cygnus Cygnus). Front. Vet. Sci. 8:670645. doi: 10.3389/fvets.2021.670645
Wang, W., Wang, F., Li, L., Wang, A., Sharshov, K., Druzyaka, A., et al. (2020). Characterization of the gut microbiome of black-necked cranes (Grus nigricollis) in six wintering areas in China. Arch. Microbiol. 202, 983–993. doi: 10.1007/s00203-019-01802-0
Wei, S., Morrison, M., and Yu, Z. (2013). Bacterial census of poultry intestinal microbiome. Poultry Sci. 92, 671–683. doi: 10.3382/ps.2012-02822
West, A. G., Waite, D. W., Deines, P., Bourne, D. G., Digby, A., McKenzie, V. J., et al. (2019). The microbiome in threatened species conservation. Biol. Conserv. 229, 85–98. doi: 10.1016/j.biocon.2018.11.016
Wienemann, T., Schmitt-Wagner, D., Meuser, K., Segelbacher, G., Schink, B., Brune, A., et al. (2011). The bacterial microbiota in the ceca of Capercaillie (Tetrao urogallus) differs between wild and captive birds. Systematic Appl. Microbiol. 34, 542–551. doi: 10.1016/j.syapm.2011.06.003
Wink, J., Mohammadipanah, F., and Hamedi, J. (2017). Biology and Biotechnology of Actinobacteria. Berlin: Springer.
Wu, Y. (2012). A Comparative Study on Food Hbbits of Red-crowned Crane and White-naped Crane during breeding period in Zhalong Nature Reserve. Harbin: Northeast Forestry University.
Wu, Y., Yang, Y., Cao, L., Yin, H., Xu, M., Wang, Z., et al. (2018). Habitat environments impacted the gut microbiome of long-distance migratory swan geese but central species conserved. Sci. Rep. 8:13314. doi: 10.1038/s41598-018-31731-9
Xie, Y., Xia, P., Wang, H., Yu, H., Giesy, J. P., Zhang, Y., et al. (2016). Effects of captivity and artificial breeding on microbiota in feces of the red-crowned crane (Grus japonensis). Sci. Rep. 6:33350. doi: 10.1038/srep33350
Yang, Y., Deng, Y., and Cao, L. (2016). Characterising the interspecific variations and convergence of gut microbiota in Anseriformes herbivores at wintering areas. Sci. Rep. 6:32655. doi: 10.1038/srep32655
Yang, Z. (2015). Energy Intake Demand of the Captive Red-crowned Crane (Grus japonensis) in Zhalong National Nature Reserve China. Harbin: Northeast Forestry University.
Zhang, F., Xiang, X., Dong, Y., Yan, S., Song, Y., and Zhou, L. (2020). Significant differences in the gut bacterial communities of Hooded Crane (Grus monacha) in different seasons at a stopover site on the flyway. Animals 10:701. doi: 10.3390/ani10040701
Zhang, J., Guo, Z., Xue, Z., Sun, Z., Zhang, M., Wang, L., et al. (2015). A phylo-functional core of gut microbiota in healthy young Chinese cohorts across lifestyles, geography and ethnicities. ISME J. 9, 1979–1990. doi: 10.1038/ismej.2015.11
Zhang, L. (2021). Charaterization of Gastrointestinal Bacterial Community and the Probiotic Properties of Lactobacillus in YAKS. China: Huazhong Agricultural University, doi: 10.27158/d.cnki.ghznu.2021.000006
Zheng, S., Piao, C., Zhang, Z., Liu, T., Wang, J., Wang, L., et al. (2019). Influence of corn-rich diet on mice gut microbiota structure and function. Acta Microbiol. Sinica 59, 1714–1722. doi: 10.13343/j.cnki.wsxb.20190163
Zou, H., Huang, H., Song, Y., and Wu, Q. (2018). Research progress on cranes in songnen plain, China. Chinese J. Wildlife 39, 433–437.
Keywords: red-crowned crane, white-naped crane, intestinal microorganisms, high-throughput sequencing technology, Zhalong Nature Reserve
Citation: Gao Z, Song H, Dong H, Ji X, Lei Z, Tian Y, Wu Y and Zou H (2022) Comparative analysis of intestinal flora between rare wild red-crowned crane and white-naped crane. Front. Microbiol. 13:1007884. doi: 10.3389/fmicb.2022.1007884
Received: 31 July 2022; Accepted: 14 November 2022;
Published: 01 December 2022.
Edited by:
Nucharin Songsasen, Smithsonian Conservation Biology Institute (SI), United StatesReviewed by:
Shrikant Bhute, University of California, San Diego, United StatesLizhi Zhou, Anhui University, China
Copyright © 2022 Gao, Song, Dong, Ji, Lei, Tian, Wu and Zou. 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: Yining Wu, d3V5aW5pbmduZWZ1QDEyNi5jb20=; Hongfei Zou, aG9uZ2ZlaXpvdUAxNjMuY29t
†These authors have contributed equally to this work