- Open Access
Composition and diversity of gut microbiota in Pomacea canaliculata in sexes and between developmental stages
BMC Microbiology volume 21, Article number: 200 (2021)
The apple snail, Pomacea canaliculata, is one of the world’s 100 worst invasive alien species and vector of some pathogens relevant to human health.
On account of the importance of gut microbiota to the host animals, we compared the communities of the intestinal microbiota from P. canaliculata collected at different developmental stages (juvenile and adult) and different sexes by using high-throughput sequencing.
The core bacteria phyla of P. canaliculata gut microbiota included Tenericutes (at an average relative abundance of 45.7 %), Firmicutes (27.85 %), Proteobacteria (11.86 %), Actinobacteria (4.45 %), and Cyanobacteria (3.61 %). The female group possessed the highest richness values, whereas the male group possessed the lowest bacterial richness and diversity compared with the female and juvenile group. Both the developmental stages and sexes had important effects on the composition of the intestinal microbiota of P. canaliculata. By LEfSe analysis, microbes from the phyla Proteobacteria and Actinobacteria were enriched in the female group, phylum Bacteroidetes was enriched in the male group, family Mycoplasmataceae and genus Leuconostoc were enriched in the juvenile group. PICRUSt analysis predicted twenty-four metabolic functions in all samples, including general function prediction, amino acid transport and metabolism, transcription, replication, recombination and repair, carbohydrate transport and metabolism, etc.
This study provided a general understanding of the diversity characteristics of intestinal microbial communities of P. canaliculata, and indicated that developmental stage and gender could both influence the intestinal microbes of P. canaliculata. Further study may focus on the interaction between the gut microbiota and their host.
Pomacea is freshwater gastropods native to South America . Several Pomacea species such as Pomacea canaliculata, P. maculata have been introduced and become severe invasive pests in many parts of the world including Asian countries, North America, islands of the Pacific, and Europe [2, 3]. P. canaliculata is listed as “100 world’ s worst invasive alien species” by the International Union for Conservative of Nature and the Invasive Species Specialist Group . Due to its high fecundity (a female snail in its lifetime could averagely spawn 13,764 eggs and reproduce 6070 young snails ), fast growth and voracious appetite for vegetation, P. canaliculata is a notorious pest, causing serious damages to aquatic crops such as rice and lotus, as well as to wetland floral diversity and ecosystem functioning . P. canaliculata is also a vector of a parasitic nematode Angiostrongylus cantonensis, which causes eosinophilic meningitis in humans . P. canaliculata was introduced to China for commercial purposes in the early 1980 s, and widely distributed in most areas of southern part of China at present .
Recent studies have highlighted the ability of the gut bacteria of animals in multiple physiological processes of their hosts, primarily including digestion, nutrition, development, reproduction, immunity, and environmental resistance [8, 9]. Bacteria in the intestine of fish can produce various enzymes, such as protease, amylase, and lipase, for digestion . Furthermore, the intestinal microorganisms of snail giant African snail (Achatina fulica) play an important role in cellulose decomposition . It is increasingly recognized on the importance of intestinal microbial structure and function, especially the potential contribution to nutrition utilization of host .
Gut microbiomes can be influenced by host development and growth stages in some animals, such as zebrafish (Danio rerio), bovine and Atlantic salmon (Salmo salar) [13,14,15]. Moreover, gender has a significant correlation with the intestinal microbiota. The alpha diversity of the gut microbiome in three-spined stickleback (Gasterosteus aculeatus) and European perch (Perca fluviatilis) is different between sex . The significant separation is found in the microbial communities between male and female fathead minnow (Pimephales promelas) based on beta diversity metrics .
Only limited studies have focused on the gut microbiota of P. canaliculata. Cheng et al.  investigated intestinal bacterial communities in P. canaliculata by denaturing gradient gel electrophoresis (DGGE). The results showed that the microflora of the male and female P. canaliculata gastric and intestinal contents was the same, containing 22 species of bacteria in all of the samples. Li et al.  compared the gut microbes of P. canaliculata at different gut sections via high-throughput sequencing of the 16 S rRNA gene. A total of 29 phyla and 111 genera of bacteria were identified in all of the samples. High-throughput sequencing opens up the possibility of conducting large-scale studies analyzing thousands of samples simultaneously to survey microbial communities at an unprecedented spatial and temporal resolution .
The differences related to developmental stages and gender in the gut microbiota of P. canaliculata have not been fully investigated. Diverse factors, such as diet, age, antibiotics, stress, psychological factors, mode of delivery, environmental factor, and exercise, can influence the status of the gut microbiota . The influence of sex on the gut microbiota is not so clear when compared with other factors such as diet and medication . Nevertheless, the effect of different genders on the gut microbiota and their interactions with other factors should be routinely analyzed. To better understand the bacterial community in the gut of this invasive snail and provide insight into their adaptive strategies in the environments, we investigated the differences related to developmental stages and gender in the gut microbiota of P. canaliculata by high-throughput sequencing in this study. Our study provided a framework for characterizing age-bacterial and sex-bacterial communities in P. canaliculata and supplemented information on intestinal content microbiota in P. canaliculata.
Materials and methods
P. canaliculata individuals were collected from an artificial pond located in Lishui City (E 119º54’, N 28º28’), Zhejiang province, China, in August 2019, including 11 adults (six males and five females, shell length ranged from 48.74mm to 57.76mm) and 15 juvenile snails (shell length ranged from 17.51mm to 19.71mm). All the testing snails sampled from artificial pond were preliminarily discerned by shell morphological analysis and using primers LCO1490/HCO2198 to amplify cytochrome c oxidase subunit I (COI) gene to identify P. canaliculata which could be used for the experiments [23, 24]. PCR amplification products were sequenced by Shanghai Sangon Biological Engineering Technology & Services Co., Ltd. P. canaliculata samples were divided into adult (≥ 30 mm shell length) and juvenile groups (10–25 mm shell length) . In addition, the adult P. canaliculata were divided into male and female groups. The shell was removed from each snail after being wiped by 75 % ethanol three times and rinsing it twice in distilled water. Dissection was performed on ice using sterilized tools. Because the coiled gut of juvenile individuals was too small, samples of three snails were collected for sequencing . Coiled gut content from three juvenile individuals were dissected and pooled to compose one juvenile biological sample and one adult individual were dissected and pooled to compose one adult biological sample . Coiled gut contents were extracted carefully to avoid rupturing the gut wall. There were five and eleven samples representing the juvenile and the adult group for intestine microbial analysis respectively. The intestinal content of each sample was placed in sterile tubes and stored at -80 °C until DNA extraction.
DNA extraction and high-throughput 16S rRNA gene sequencing
Extraction of microbial DNA from approximately 0.02 g of coiled gut contents extracted from each P. canaliculata sample was performed using the FastDNA® Kit (MP Biomedicals, CA, USA). The V3-V4 hyper-variable regions of the 16S rRNA genes were amplified by using the specific primer 338F (5’-ACTCCTACGGGAGGCAGCAG-3’) and 806R (5’-GGACTACHVGGGTWTCTAAT-3’). The PCR procedure was as follows: denaturation at 95 °C for 3 min, 27 cycles of 95 °C denaturation for 30 s, 55 °C annealing for 30 s, and 72 °C extension for 45 s, and a final extension at 72 °C for 10 min. PCR products were detected by 2 % agarose gel and further purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) and quantified using QuantiFluor™-ST (Promega, USA). Purified amplicons were pooled and sequenced according to the standard protocols using the Illumina MiSeq platform (Majorbio Bio-Pharm Technology Co. Ltd., Shanghai, China).
Statistical and bioinformatics analyses
The raw sequences were processed by Quantitative Insights into Microbial Ecology (QIIME) software (Version1.9.1, http://qiime.org/install/index.html). Raw fastq files were quality-filtered by Trimmomatic and merged by FLASH according to the following criteria: (i) The reads were truncated at any site receiving an average quality score < 20 over a 50 bp sliding window. (ii) Sequences overlapped longer than 10 bp were merged according to their overlap with mismatch no more than 2 bp. (iii) Sequences of each sample were separated according to barcodes (exactly matching) and primers (allowing 2 nucleotide mismatching), and reads containing ambiguous bases were removed. All of the remaining high-quality sequences were clustered into operational taxonomic units (OTUs) at a 97 % identity threshold using USEARCH (Version 7.0, http://drive5.com/uparse/) . Rarefaction curves were generated based on normalized OTU numbers using Mothur software (Version 1.30.2, https://www.mothur.org/wiki/Download_mothur) . Using RDP Native Bayesian Classifier (Version 2.11), the taxonomy of each 16S rRNA gene sequence was analyzed against the Silva 16S rRNA database (Version 132, http://www.arb-silva.de) . Subsequently, BLAST search was further carried out on NCBI nucleotide collection (nr/nt) using the BLASTn algorithm to obtain information on taxonomic identity of the noranked OTUs on genus level, and the closest matches to bacterial strains were obtained.
Community diversity was estimated using alpha diversity indices including the Chao and Shannon indices using Mothur (Version 1.30.2). To display any discrepancy among the three groups (female; male; juvenile), nonmetric multidimensional scaling (NMDS) analysis was used based on the unweighted and weighted Unifrac distances and analysis of similarity (ANOSIM) based on 999 permutations. Based on Bray-Curtis distance, analysis of molecular variance (AMOVA), was performed to analyze the differences between the female and male group, and between the juvenile and adult group by Mothur software (Version 1.30.2). The differences in intestinal microbial composition among the three groups (female; male; juvenile) were analyzed by linear discriminant analysis coupled with effect size (LEfSe) analysis (LDA score > 4) and Kruskal-Wallis H test . Based on the evolutionary genealogy of genes from Non-supervised Orthologous Groups (EggNOG) database and the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, PICRUSt software was used to predict microbial functions and functional pathways respectively . LEfSe analysis (LDA score > 2) was used to identify significantly enriched functions and KEGG pathways among the three groups (female; male; juvenile). To identify statistically significant differences of alpha diversity between the groups, Wilcoxon rank-sum test was used by SPSS 19.0 software (IBM, Armonk, USA). P-value < 0.05 was considered statistically significant.
Sequencing depth and alpha diversity indices
DNA extracted from sixteen P. canaliculata samples were amplified successfully, and 559,712 valid sequences were obtained. P. canaliculata yielded 2,098 valid OTUs at a 97 % identity. By comparing the OTUs with SILVA 132 database for species annotation, 2098 OTUs were annotated to 30 bacterial phyla, 69 classes, 197 orders, 352 families, 715 genera and 1173 species. Noranked OTUs on genus level were showed in Supplementary file 1 (Datafile S1). The rarefaction curve of the Shannon index on OTU level reached asymptote (Supplementary file 3: Fig. S1), which indicated that the sequencing depth was sufficient to represent the majority of species richness in each sample.
The alpha diversity measures (Shannon index and Chao diversity index) were calculated in each group to examine whether the female, male and juvenile P. canaliculata had differences in alpha diversity. The female group had the highest Shannon and Chao richness values, whereas the male group had the lowest bacterial richness and diversity. With reference to Shannon indices, the female group (3.35 ± SD 0.43) was significantly higher than the juvenile group (1.99 ± 0.71) (P < 0.05) and the male group (1.89 ± 0.40) (P < 0.01), while the Shannon indices of the juvenile group was not significantly different from the male group (P > 0.05) (Fig. 1a). With reference to Chao indices, both female group (1106.90 ± 133.68) and juvenile group (1068.7 ± 192.51) were significantly higher than the male group (723.44 ± 224.25) (P < 0.05). While the Chao indices between female group and juvenile group had no significant difference (P > 0.05) (Fig. 1b).
Taxonomic composition and beta diversity analysis
In terms of prominent bacteria (relative abundance > 1 % based on classifications in any group), P. canaliculata harbored bacteria from seven phyla, including Tenericutes, with an average relative abundance of 45.7 %, Firmicutes (27.85 %), Proteobacteria (11.86 %), Actinobacteria (4.45 %), Cyanobacteria (3.61 %), Chloroflexi (3.4 %), Bacteroidetes (2.25 %) (Fig. 2a) and eleven genera, including Leuconostoc (45.68 %), Lactococcus (20.64 %), Bacillus (4.19 %), Enterobacter (3.25 %), Cloacibacterium (2.18 %), Mycobacterium (1.77 %), Caldilinea (1.52 %), Litorilinea (1.51 %), unclassified genus from order Rhizobiales (1.37 %), Nordella (1.22 %), Methylocystis (1.1 %) (Fig. 2b).
To evaluate the overall difference in the beta diversity, NMDS analysis based on the unweighted unifrac distances showed that the three groups (female; male; juvenile) could be separated from each other. The ANOSIM analysis revealed significant differences in the structure of gut microbiota among different groups (R = 0.6282, P = 0.001; Fig. 3a). However, when we used weighted UniFrac distance based ANOSIM analysis to account for the abundance information, the microbiota of the male and juvenile group clustered closely and could be separated from the female group (R = 0.4131, P = 0.002; Fig. 3b).
AMOVA analysis was performed to further compare the intestinal microbial communities between different sex and between different developmental stages. The result showed that the intestinal microbial communities were significantly different between female and male group based on Bray-Curtis distance analyses (P < 0.01). The intestinal microbial communities were also significantly different between juvenile and adult group (including female and male group) based on Bray-Curtis distance analyses (P < 0.05).
Microbial community similarities and differences
A Venn diagram showed that a total of 696 OTUs were identified as core bacterial OTUs in all groups. The number of unique OTUs in each group was 672 (Female), 304 (Male), and 158 (Juvenile), respectively (Fig. 4a).
Female, male and juvenile P. canaliculata exhibited different bacterial taxa in their coiled guts (Fig. 3). According to the Kruskal-Wallis H test, we found some significant differences in the abundances of the bacterial communities among the samples from female, male and juvenile at phylum and genus levels (Supplementary file 2: Datafile S2). There were significant differences in the abundance of Tenericutes [H(3,N = 16) = 9.706, P = 0.008], Proteobacteria [H(3,N = 16) = 9.11, P = 0.011], Actinobacteria [H(3,N = 16) = 8.169, P = 0.017] and Bacteroidetes [H(3,N = 16) = 6.11, P = 0.047] in female, male and juvenile group at phylum level (P < 0.05). The 17 microbial genera with statistically significant differences in their abundances (P < 0.05) mainly included Leuconostoc [H(3,N = 16) = 9.706, P = 0.008], Caldilinea [H(3,N = 16) = 9.934, P = 0.007], Lysinibacillus [H(3,N = 16) = 9.993, P = 0.007], Lactobacillus [H(3,N = 16) = 11.812, P = 0.003]. We also performed linear discriminant analysis effect size (LEfSe) tests on the samples to detect any relative abundance differences (average relative abundance > 1 %) in the bacterial taxon (including phylum, class, order, and family). The LEfSe results were similar to the Kruskal-Wallis H test. LEfSe identified 18 discriminative features (LDA score > 4) whose relative abundance varied significantly between groups (Fig. 4b). Microbiotas of the female group were enriched in phyla Proteobacteria, Actinobacteria and Chloroflexi, classes Alphaproteobacteria, Actinobacteria and Anaerolineae, orders Rhizobiales, SBR1031, family Leptolyngbyaceae and genus Caldilinea. Microbiotas of the male group were enriched in phylum Bacteroidetes and class Bacteroidia. Microbiotas of the juvenile group were enriched in phylum Tenericutes, class Mollicutes, order Mycoplasmatales, family Mycoplasmataceae and genus Leuconostoc.
Prediction of bacterial functions in P. canaliculata
The result of PICRUSt analysis based on the 16S rRNA composition data of each sample from EggNOG database showed a total of 24 metabolic functions predicted in all samples (Fig. 5). Function unknown accounted for the highest proportion (10.2 %), followed by general function prediction (8.69 %), amino acid transport and metabolism (8.6 %), transcription (6.82 %), replication, recombination and repair (6.68 %), carbohydrate transport and metabolism (6.65 %), energy production and conversion (6.53 %), cell wall/membrane/envelope biogenesis (6.19 %), inorganic ion transport and metabolism (6.15 %), translation, ribosomal structure and biogenesis (5.59 %), signal transduction mechanisms (5.04 %) and lipid transport and metabolism (4.34 %), etc. To better understand the functional differences, LEfSe test (LDA > 2) was performed in the P. canaliculata gut microbiotas from the three groups (female; male; juvenile). The result showed that general function prediction was enriched in the male group, while intracellular trafficking secretion and vesicular transport were enriched in the juvenile group (Supplementary file 4: Fig. S2).
In addition, predicted functional pathways were annotated using the KEGG database and a total of 5416 KEGG Orthology (KOs) were mapped to 299 level 3 KEGG pathways and were then classified into 41 level 2 KEGG pathways. Of these 41 secondary KEGG pathways, 12 pathways were involved in metabolism, 8 pathways involved in organismal systems, 6 pathways involved in human diseases, 4 pathways involved in cellular processes, 4 pathways involved in genetic information processing, 4 pathways involved in unclassified functions, and 3 pathways involved in environmental information processing. LEfSe identified 3 pathways (LDA score > 2) whose relative abundance varied significantly between groups. Environmental information processing (level 1) and secretion system (level 3) were enriched in the juvenile group, carbohydrate metabolism (level 2) was enriched in the male group (Supplementary file 5: Fig. S3).
The gut microbiota of animals plays an important role in food ingestion, digestion and nutrient absorption of the host . During past years, only limited studies have focused on the gut microbiota of P. canaliculata [18, 19].
In this study, there were no dramatic differences in the bacterial composition of the female, male and juvenile group at the phylum level. The top 3 phyla in the three groups were Tenericutes, Firmicutes and Proteobacteria. The fourth abundant phylum in the intestinal contents of P. canaliculata was Actinobacteria by using high-throughput sequencing. Similar findings were also obtained in other snails. Proteobacteria (the relative abundance of juvenile was 36.0 %, the relative abundance of adult was 31.6 %), Firmicutes (Juvenile: 14.4 %, Adult: 6.7 %), Actinobacteria (Juvenile: 8.2 %, Adult: 12.6 %), and Tenericutes (Juvenile: 7.3 %, Adult: 6.2 %) were the predominant phyla of gut microbiota in the snail Radix auricularia . Proteobacteria, Firmicutes, and Actinobacteria were also identified as the dominant bacterial taxa at the phylum level in the gut microbiota of adult Oncomelania hupensis .
The spatial structure of the microbiota in the intestine of female P. canaliculata showed that the abundance of Tenericutes was higher in the intestine than in the buccal mass and stomach . Tenericutes was also the most abundant microbiota in the gut content of P. canaliculata in the present study. However, at the genus level, there were some differences, which might be attributed to different sample collection methods.
Firmicutes was one of the most abundant bacterial phyla in the gut content of P. canaliculata, which may play an important role in protein degradation . P. canaliculata has a putative symbiont (plausibly a cyanobacterium) in the digestive gland, which travels from the stomach to the coiled gut [37,38,39]. The cyanobacterium produces a protease of 30 kDa that participates in the digestion of proteins . In addition, Firmicutes has been reported to be able to promote preservation of gut homeostasis and host immunity development . Pomacea has a well-developed innate cellular system and it has organs (lung and kidney) that can act as barriers against different antigens [40, 41]. In this study, some gut bacteria of P. canaliculata may play important roles in immunity. For example, Lactococcus (phylum Firmicutes; family Streptococcaceae), one of the most abundant bacterial genera in the intestinal content of P. canaliculata, was identified as dominant genus in lactic acid microflora of the gut of snails Cornu aspersum and Oncomelania hupensis [35, 42]. In addition, L. lactis subsp. lactis from the Arapaima gigas fish displayed in vitro antagonism against all 10 pathogens tested, including Citrobacter freundii, Pseudomonas sp., Citrobacter freundii, Enterobacter sp., Pseudomonas stutzeri, Pseudomonas sp., Micrococcus luteus, Aeromonas hydrophila ATCC 7966, Staphylococcus agalactiae and Aeromonas hydrophila DRM CPQBA 228-08 . The genus Bacillus (family Bacillaceae, phylum Firmicutes) was identified as the core element of the intestinal microbiota of snail Oncomelania hupensis , which was also found in this study. Bacillus cereus is common in soil and food, and some strains are pathogenic to humans, while others act as probiotics for animals . Bacillus cereus var. toyoi has been reported to be able to promote the growth and immune status of livestock and fish .
Proteobacteria, which was identified as the dominant phylum in the gut bacterial community of P. canaliculata, was also reported in the gut bacteria of other snails, such as two different terrestrial snails, Achatina fulica, Helix pomatia [46, 47], and three freshwater snails, Biomphalaria pfeifferi, Bulinus africanus and Helisoma duryi . In previous studies, 70 % of the isolated cellulolytic bacteria from the gut of Holotrichia parallela larvae were Proteobacteria . The genus Enterobacter (family Enterobacteriaceae, phylum Proteobacteria) was found in the gut of P. canaliculata. Enterobacter has been identified as cellulolytic species and reported as predominant in A. fulica and H. pomatia [46, 47], associated with carboxymethyl cellulase (CMCase) activity [46, 49, 50]. Enterobacter was isolated from the coiled gut of female P. canaliculata. Uricase specific activity could be determined in Enterobacter. The uric acid-degrading bacterium may participate in recycling of combined nitrogen in P. canaliculata . Under a long period of hypometabolism, P. canaliculata using uric acid as a POS (Preparation for Oxidative Stress) like strategy to tolerate prolonged drought and endure low temperatures .
P. canaliculata is a voracious herbivorous snail with great environmental and ecological importance . Its capacity to process a broad variety of vegetable organic matter is due to the presence of cellulolytic enzymes, both from the gut microbiome and stomach of the snail [53, 54]. The bacterial communities inside the gut of the snail may have crucial importance in the digestion of cellulose and other plant wall components . Higher proportions of Proteobacteria and Firmicutes were often associated with diets containing plant ingredients . Based on a culture-dependent cellulose-degrading bacteria screening method, Actinobacteria representatives could be easily recovered from the intestinal tract of A. fulica and could be cultivated to produce a wide range of glycoside hydrolases . In addition, cellulases A0A2T7PPN6 (GHF9), A0A2T7NYY1(GHF10), A0A2T7NZR0 (GHF10), and Pc89752 (GHF10) were uniquely found in the digestive contents indicating a possible origin in unicellular glands of the gut or in commensal symbiotic organisms . A bacterial cellulase secreted by Bacillus sp. from Ampullaria crossean has been characterized at the molecular and biochemical level [57, 58]. Detection of these and other cellulase-producing bacteria in the present study is consistent with the voracious herbivorous diet of P. canaliculata. These bacteria could help P. canaliculata to digest and absorb nutrients from the plant materials in their diets via the microbe production of various digestive enzymes to decompose various cellulose and hemicellulose.
There were distinct differences in the richness and diversity of the microbial community of P. canaliculata from the female, male and juvenile groups. Male group showed the lowest richness and diversity values and the female group showed the highest richness and diversity values.
The gut microbial analysis revealed significant differences related to developmental stages and gender indicated by NMDS analysis based on unweighted UniFrac distance and the AMOVA analysis (P < 0.05), indicating that sex and developmental stage have effect on gut bacterial community. In previous studies, higher ingestion rates of juvenile snails in relation to those adults have been reported in P. canaliculata . P. canaliculata of small size have higher foraging and competitive abilities than big size . There were some differences related to diet between juvenile and adult P. canaliculata. Adult P. canaliculata consume Eichhornia crassipes at a high rate while juvenile P. canaliculata snails do not consume this macrophyte at all [61, 62]. In addition, snails of different sizes could also use different strategies when fed on the same macrophyte species . Bacterial microbiota from juvenile and adult snails were significantly different in precious reports in which the changes of the microbiota were observed at different developmental stages of the aquatic animals, such as southern catfish (Silurus meridionalis) and white cachama (Piaractus brachypomus) [63, 64], proposing that as the aquatic animal transits through youthhood to adulthood, the microbiota drifts, even when the environment remained constant.
In addition, higher growth rates, ingestion rates and growth efficiencies of female snails in relation to those of male counterparts have been reported in P. canaliculata, probably resulting from a higher waste of undigested food related to smaller mid-gut gland in male, which in turn results in lower food assimilation . Females of P. canaliculata reach larger sizes than males, both in natural habitats and reared in the laboratory, a pattern that results from higher growth rates rather than from different survivorship rates [66, 67]. An indoor incubation test found that at both hungry and full states, the females had a significantly higher foraging rate than the males . The higher abundance of Proteobacteria and Actinobacteria in the female group was found by LEfSe analysis. Previous studies have shown that Proteobacteria (67.13 %) and Actinobacteria (23.15 %) are representative of cellulolytic bacterial community in the hindgut of Holotrichia parallela larvae . The phyla Proteobacteria and Actinobacteria play important roles in digesting plant and providing energy and nutrients for the host [11, 49], which may be related to higher growth rates of the female.
However, NMDS analysis based on weighted UniFrac distance showed that juvenile and male groups cluster together. This is at least in part due to abundance information, which can obscure significant patterns of variation in the taxa that are present [69, 70], indicating that taking the abundance of bacterial taxon into account reveals similarities between juvenile and male populations.
Leuconostoc (family Streptocnccaceae) was the most abundant genus in the gut of P. canaliculata and it was enriched in the gut microbiotas of the juvenile group by LEfSe analysis. Leuconostoc is present in many plant materials such as vegetables, silage and fermented food products and is highly beneficial to the host by fermenting various dairy products , which is related to the diet of P. canaliculata . Leuconostoc bacteria are gram-positive, which may play important roles in the production of polysaccharides, mannitol, vitamins-K, bacteriocins and the hydrolysis of α-galactosides, which may be related to higher ingestion rates for juvenile snails than adult ones [71,72,73,74]. Recent studies have been reported that Leuconostoc pseudomesenteroides could significantly restore intestinal disorder caused by a high-fat diet .
The family Mycoplasmataceae, was also found to be enriched in the gut microbiotas of the juvenile group, related to the changes of the host environment. For example, the gut microbiome of marine S. salar with a higher abundance of Mycoplasmataceae was less rich and diverse than that of freshwater juveniles . Moreover, KEGG pathways of environmental information processing (level 1) and secretion system (level 3), function prediction of intracellular trafficking secretion and vesicular transport were enriched in the juvenile, which might be related to the differences between juvenile and adult and could be very important in allowing juvenile to adapt to a complex digestive environment.
In the present study, Bacteroidetes was found to be enriched in the gut microbiotas of the male group. Bacteroidetes is known for fermentative metabolism and degradation of oligosaccharides derived from plant material . Thus, the presence of Bacteroidetes might make an effect on the male group and allowed them to maximize the energy. Moreover, KEGG pathways of carbohydrate metabolism (level 2) were enriched in the male group. However, further study is needed to investigate the underlying mechanism on the increased Bacteroidetes in the male group. The results in our study may serve as a preliminary indication for the function of bacteria communities. Further analyses of metagenomic and metatranscriptomic approaches are required to illustrate the interactions of hosts and microbiota, as well as microbial structure and function in P. canaliculata.
Our study provided a general understanding of the diversity characteristics of the intestinal microbial communities of P. canaliculata. It indicated that developmental stage and gender could both influence the intestinal microbes of P. canaliculata. We also found some gut microbiome associated with diets containing plant ingredients. This study can provide insight into their adaptive strategies in the environments. Further study will focus on the interaction between the gut microbiota and the host.
Availability of data and materials
The raw data are available from the SRA database (PRJNA673528).
Operational taxonomic unit
Nonmetric multidimensional scaling analysis
Analysis of similarity
Analysis of molecular variance
Linear discriminant analysis effect size
Non-supervised Orthologous Groups
Encyclopedia of Genes and Genomes
Hayes KA, Burks RL, Castro-Vazquez A, Darby PC, Heras H, Martín PR et al. Insights from an integrated view of the biology of apple snails (Caenogastropoda: Ampullariidae). Malacologia. 2015;58(1–2):245–302.
Hayes KA, Joshi RC, Tiengo SC, Cowie RH. Out of South America: multiple origins of non-native apple snails in Asia. Diver Distrib. 2008;14(4):701–12.
Rawlings TA, Hayes KA, Cowie RH, Collins TM. The identity, distribution, and impacts of non-native apple snails in the continental United States. BMC Evol Biol. 2007;7(1):97.
Lowe S, Browne M, Boudjelas S, De Poorter M. 100 of the world’s worst invasive alien species a selection from the Global Invasive Species Database. The Invasive Species Specialist Group of the Species Survival Commission of the World Conservation Union. 2000; 12 pp.
Liu J, He YJ, Tan JC, Xu CX, Zhong L, Wang ZG, Liao QG. Characteristics of Pomacea canaliculata reproduction under natural conditions. Chin J Appl Ecol. 2012;23(2):559–65.
Carlsson NOL, Bronmark C, Hansson LA. Invading herbivory: the golden apple snail alters ecosystem functioning in Asian wetlands. Ecology. 2004;85(6):1575–80.
Lv S, Zhang Y, Liu HX, Hu L, Liu Q, Wei FR, Guo YH, Steinmann P, Hu W, Zhou XN, Utzinger J. Phylogenetic evidence for multiple and secondary introductions of invasive snails: Pomacea species in the People’s Republic of China. Diver Distrib. 2013;19(2):147–56.
Nayak SK. Role of gastrointestinal microbiota in fish. Aquacult Res. 2010;41(11):1553–73.
Sommer F, Bäckhed F. The gut microbiota-masters of host development and physiology. Nat Rev Microbiol. 2013;11(4):227–38.
Bairagi A, Ghosh KS, Sen SK, Ray AK. Enzyme producing bacterial flora isolated from fish digestive tracts. Aquacult Int. 2002;10:109–21.
Pinheiro GL, Correa RF, Cunha RS, Cardoso AM, Chaia C, Clementino MM, Garcia ES, De Souza W, Frasés S. Isolation of aerobic cultivable cellulolytic bacteria from different regions of the gastrointestinal tract of giant land snail Achatina fulica. Front Microbiol. 2015;6:860.
Chaiyapechara S, Rungrassamee W, Suriyachay I, Kuncharin Y, Klanchui A, Karoonuthaisiri N, Jiravanichpaisal P. Bacterial community associated with the intestinal tract of P. monodon in commercial farms. Microb Ecol. 2012;63(4):938–53.
Llewellyn MS, McGinnity P, Dionne M, Letourneau J, Thonier F, Carvalho GR, Creer S, Derome N. The biogeography of the atlantic salmon (Salmo salar) gut microbiome. ISME J. 2016;10(5):1280–4.
Jami E, Israel A, Kotser A, Mizrahi I. Exploring the bovine rumen bacterial community from birth to adulthood. ISME J. 2013;7(6):1069–79.
Stephens WZ, Burns AR, Stagaman K, Wong S, Rawls JF, Guillemin K, Bohannan BJM. The composition of the zebrafish intestinal microbial community varies across development. ISME J. 2015;10(3):644–54.
Bolnick DI, Snowberg LK, Hirsch PE, Lauber CL, Org E, Parks B, Lusis AJ, Knight R, Caporaso JG, Svanbäck R. Individual diet has sex-dependent effects on vertebrate gut microbiota. Nat Commun. 2014;5(1):4500.
Debofsky A, Xie Y, Grimard C, Alcaraz AJ, Brinkmann M, Hecker M, Giesya JP. Differential responses of gut microbiota of male and female fathead minnow (Pimephales promelas) to a short-term environmentally-relevant, aqueous exposure to benzo[a]pyrene. Chemosphere. 2020;252:126461.
Cheng TY, Qu LX, Wang XJ. The flora of the Pomacea canaliculata gastric and intestinal contents with PCR-DGGE. Chin Agricult Sci Bull. 2015;11:98–102.
Li LH, Lv S, Lu Y, Bi DQ, Guo YH, Wu JT, Yue ZY, Mao GY, Guo ZX, Zhang Y, Tang YF. Spatial structure of the microbiome in the gut of Pomacea canaliculata. BMC Microbiol. 2019;19(1):273.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, Fierer N, Knight R. Global patterns of 16s rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci USA. 2011;108(Supplement 1):4516–22.
Osadchiy V, Martin CR, Mayer EA. The gut-brain axis and the microbiome: mechanisms and clinical implications. Clin Gastroenterol Hepatol. 2019;17(2):322–32.
Kim YS, Unno T, Kim BY, Park MS. Sex differences in gut microbiota. World J Mens Health. 2020;38(1):48–60.
Hayes KA, Cowie RH, Thiengo SC, Strong EE. Comparing apples with apples: clarifying the identities of two highly invasive Neotropical Ampwlarida& (Gaenogastonoda). Zool J Linnean Soc. 2012;166(4):723–53.
Yang Q, Liu S, He C, Cowie RH, Yua XP, Hayes KA. Invisible apple snail invasions: importance of continued vigilance and rigorous taxonomic assessments. Pest Manage Sci. 2019;75(5):1277–86.
Bernatis JL, Mcgaw IJ, Cross CL. Abiotic tolerances in different life stages of apple snails Pomacea canaliculata and Pomacea maculata and the implications for distribution. J Shellfish Res. 2016;35(4):1013–25.
Hu Z, Tong Q, Chang J, Yu J, Li S, Niu H, Ma D. Gut bacterial communities in the freshwater snail Planorbella trivolvis and their modification by a non-herbivorous diet. PeerJ. 2021;9:e10716.
Godoy MS, Castro-Vazquez A, Vega IA. Correction: Endosymbiotic and host proteases in the digestive tract of the invasive snail Pomacea canaliculata: Diversity, origin and characterization. PLoS One. 2013;8(7):10.
Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10(10):996–8.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(D1):D590–6.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6):R60.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40:D109–14.
Pawar KD, Banskar S, Rane SD, Charan SS, Kulkarni GJ, Sawant SS, Ghate HV, Patole MS, Shouche YS. Bacterial diversity in different regions of gastrointestinal tract of giant African snail (Achatina fulica). Microbiologyopen. 2012;1(4):415–26.
Hu Z, Chen X, Chang J, Yu J, Tong Q, Li S, Niu H. Compositional and predicted functional analysis of the gut microbiota of Radix auricularia (Linnaeus) via high-throughput Illumina sequencing. PeerJ. 2018;6(4):e5537.
Hao Y, Guan W, Wu H, Li L, Abe EM, Xue J, Qin Z, Wang Q, Lv S, Xu J, Wang W, Li S. Intestinal microbiome profiles in Oncomelania hupensis in mainland China. Acta Trop. 2020;201:105202.
Colston TJ, Jackson CR. Microbiome evolution along divergent branches of the vertebrate tree of life: What is known and unknown. Mol Ecol. 2016;25(16):3776–800.
Vega IA, Damborenea MC, Gamarra-Luques C, Koch E, Cueto JA, Castro-Vazquez A. Facultative and obligate symbiotic associations of Pomacea canaliculata (Caenogastropoda, Ampullariidae). Biocell. 2006;30(2):367–75.
Dellagnola FA, Rodriguez C, Castro-Vazquez A, Vega IA. A multiple comparative study of putative endosymbionts in three coexisting apple snail species. PeerJ. 2019;7:e8125.
Vega IA, Gamarra-Luques C, Koeh E, Bussmann LE, Castro-Vazquez A. A study of corpuscular DNA and midgut gland occupancy by putative symbiotic elements in Pomacea canaliculata (Caenogastropoda, Ampullariidae). Symbiosis. 2005;39(1):37–45.
Cueto JA, Rodriguez C, Vega IA, CastroVazquez A. Immune defenses of the invasive apple snail Pomacea canaliculata (Caenogastropoda, Ampullariidae): phagocytic hemocytes in the circulation and the kidney. PLoS One. 2015;10(4):e0123964.
Rodriguez C, Prieto GI, Vega IA, Castro-Vazquez A. Assessment of the kidney and lung as immune barriers and hematopoietic sites in the invasive apple snail Pomacea canaliculata. PeerJ. 2018;6:e5789.
Koleva Z, Dedov I, Kizheva J, Lipovanska R, Moncheva P, Hristova P. Lactic acid microflora of the gut of snail Cornu aspersum. Biotechnol Biotechnol Equip. 2014;28(4):627–34.
Gabriella DVP, Cunha DGD, Pedreira Mourino JL, Rodiles A, Jaramillo-Torres A, Merrifield DL. Characterization of microbiota in Arapaima gigas intestine and isolation of potential probiotic bacteria. J Appl Microbiol. 2017;123(5):1298–311.
Gao T, Ding Y, Wu Q, Wang J, Zhang J, Yu S, Yu P, Liu C, Kong L, Feng Z, Chen M, Wu S, Zeng H, Wu H. Prevalence, virulence genes, antimicrobial susceptibility, and genetic diversity of Bacillus cereus isolated from pasteurized milk in China. Front Microbiol. 2018;9:533.
Gisbert E, Castillo M, Skalli A, Andree KB, Badiola I. Bacillus cereus var. Toyoi promotes growth, affects the histological organization and microbiota of the intestinal mucosa in rainbow trout fingerlings. J Animal Sci. 2013;91(6):2766–74.
Pawar KD, Dar MA, Rajput BP, Kulkarni GJ. Enrichment and identification of cellulolytic bacteria from the gastrointestinal tract of giant African snail, Achatina fulica. Appl Biochemist Biotechnol. 2015;175(4):1971–80.
Nicolai A, Rouland-Lefevre C, Ansart A, Filser J, Lenz R, Pando A, Charrier M. Inter-population differences and seasonal dynamic of the bacterial gut community in the endangered land snail Helix pomatia (gastropoda: helicidae). Malacologia. 2015;59(1):177–90.
Van Horn DJ, Garcia JR, Loker ES, Mitchell KR, Mkoji GM, Adema CM, Takacs-Vesbach CD. Complex intestinal bacterial communities in three species of planorbid snails. J Molluscan Stud. 2012;78(1):74–80.
Huang S, Sheng P, Zhang H. Isolation and identification of cellulolytic bacteria from the gut of Holotrichia parallela larvae (Coleoptera: Scarabaeidae). Int J Mol Sci. 2012;13(3):2563–77.
Saha S, Roy RN, Sen SK, Ray AK. Characterization of cellulase-producing bacteria from the digestive tract of tilapia, Oreochromis mossambica (Peters) and grass carp, Ctenopharyngodon idella (Valenciennes). Aquacult Res. 2006;37(4):380–8.
Koch E, Lozada M, Dionisi H, Castro-Vazquez A. Uric acid-degrading bacteria in the gut of the invading apple snail Pomacea canaliculata and their possible symbiotic significance. Symbiosis. 2014;63:149–55.
Giraud-Billoud M, Rivera-Ingraham GA, Moreira DC, Burmester T, Castro-Vazquez A, Carvajalino-Fernández JM, Dafre A, Niu C, Tremblay N, Paital B, Rosa R, Storey JM, Vega IA, Zhang W, Yepiz-Plascencia G, Zenteno-Savin T, Storey KB, Hermes-Lima M. Twenty years of the‘Preparation for Oxidative Stress’ (POS) theory: Ecophysiological advantages and molecular strategies. Comp Biochem Physiol A Mol Integr Physiol. 2019;234:36–49.
Imjongjirak C, Amparyup P, Sittipraneed S. Cloning, genomic organization and expression of two glycosyl hydrolase family 10 (GHF10) genes from golden apple snail (Pomacea canaliculata). DNA Seq. 2008;19(3):224–36.
Liu C, Zhang Y, Ren Y, Wang H, Li S, Jiang F, Yin L, Qiao X, Zhang G, Qian W, Liu B, Fan W. The genome of the golden apple snail Pomacea canaliculata provides insight into stress tolerance and invasive adaptation. Gigascience. 2018;7(9):giy101.
Escobar-Correas S, Mendoza-Porras O, Dellagnola FA, Colgrave ML, Vega IA. Integrative proteomic analysis of digestive tract glycosidases from the invasive golden apple snail, Pomacea canaliculata. J Proteome Res. 2019;18(9):3342–52.
Rimoldi S, Terova G, Ascione C, Giannico R, Brambilla F. Next generation sequencing for gut microbiome characterization in rainbow trout (Oncorhynchus mykiss) fed animal by-product meals as an alternative to fishmeal protein sources. PLoS One. 2018;13(3):e0193652.
Zhang S, Yin QY, Li YH, Ding M, Xu GJ, Zhao FK. Molecular and biochemical characterization of Ba-EGA, a cellulase secreted by Bacillus sp. AC-1 from Ampullaria crosseans. Appl Microbiol Biotechnol. 2007;75(6):1327–34.
Li YH, Ding M, Wang J, Xu GJ, Zhao F. A novel thermoacidophilic endoglucanase, Ba-EGA, from a new cellulosedegrading bacterium, Bacillus sp. AC-1. Appl Microbiol Biotechnol. 2006;70(4):430–6.
Boland BB, Meerhoff M, Fosalba C, Mazzeo N, Barnes MA, Burks RL. Juvenile snails, adult appetites: contrasting resource consumption between two species of apple snails (Pomacea). J Molluscan Stud. 2008;74(1):47–54.
Carlsson NOL, Brnmark C. Size dependent effects of an invasive herbivorous snail (Pomacea canaliculata) on macrophytes and periphyton in Asian wetlands. Freshwater Biol. 2006;51(4):695–704.
Carlsson NOL, Lacoursière JO. Herbivory on aquatic vascular plants by the introduced golden apple snail (Pomacea canaliculata) in Lao PDR. Biol Invasions. 2005;7:233–41.
Lach L, Britton DK, Rundell RJ, Cowie RH. Food preference and reproductive plasticity in an invasive freshwater snails. Biol Invasions. 2000;2:279–88.
Zhang Z, Li D, Refaey MM, Xu W, Tang R, Li L. Host age affects the development of Southern catfish gut bacterial community divergent from that in the food and rearing water. Front Microbiol. 2018;9:495.
Castañeda-Monsalve VA, Junca H, García-Bonilla E, Montoya-Campuzano OI, Moreno-Herrera CX. Characterization of the gastrointestinal bacterial microbiome of farmed juvenile and adult white Cachama (Piaractus brachypomus). Aquaculture. 2019;521(15):734325.
Tamburi NE, Martín PR. Feeding rates and food conversion efficiencies in the apple snail Pomacea canaliculata (Caenogastropoda: Ampullariidae). Malacologia. 2009;51(2):221–32.
Martín PR, Estebenet AL. Interpopulation variation in life-history traits of Pomacea canaliculata (Gastropoda: Ampullariidae) in Southwestern Buenos Aires Province, Argentina. Malacologia. 2002;44(1):153–63.
Estebenet AL, Martin PR, Burela S. Conchological variation in Pomacea canaliculata and other South American Ampullariidae (Caenogastropoda, Architaenioglossa). Biocell. 2006;30:329–35.
Xu WB, Zhong QH, Li LF, Zhang JE, Luo MZ, Zhuo BL, Qin Z. Differences of female and male Pomacea canaliculata in foraging and heterosexual choice. Chin J Ecol. 2011;30(11):2528–33.
Lozupone CA, Hamady M, Kelley ST, Knight R. Quantitative and qualitative diversity measures lead to different insights into factors that structure microbial communities. Appl Environ Microbiol. 2007;73(5):1576–85.
Chang Q, Luan Y, Sun F. Variance adjusted weighted UniFrac: a powerful beta diversity measure for comparing communities based on phylogeny. BMC Bioinformatics. 2011;12(1):118.
Hemme D, Foucaud-Scheunemann C. Leuconostoc, characteristics, use in dairy technology and prospects in functional foods. Int Dairy J. 2004;14(6):467–94.
Stiles ME. Bacteriocins produced by Leuconostoc species. J Dairy Sci. 1994;77(9):2718–24.
Sybesma W, Starrenburg M, Tijsseling L, Hoefnagel MHN, Hugenholtz J. Effect of cultivation conditions on folate production by lactic bacteria. Appl Environ Microbiol. 2003;69(8):452–8.
Weymarn N, Hujanen M, Leisola M. Production of D-mannitol by heterofermentatative lactic acid bacteria. Process Biochem. 2002;37(11):1207–13.
Sun M, Wang Q, Zhang MM, Zhang GH, Wu T, Liu R, Sui W, Zhang JJ, Yin JJ, Zhang M. Leuconostoc pseudomesenteroides improves microbiota dysbiosis and liver metabolism imbalance and ameliorates the correlation between dihydroceramide and strains of Firmicutes and Proteobacteria in high fat diet obese mice. Food Function. 2020;11:6855–65.
Ghanbari M, Kneifel W, Domig KJ. A new view of the fish gut microbiome: advances from next-generation sequencing. Aquaculture. 2015;448:464–75.
We thank Guohua Ding for help with snail collection.
This research was sponsored by the Natural Science Foundation of China (No.31770402 and 31670422), the Natural Science Foundation of Jiangsu Province (No. BK20171407), “333 high level talents project” and Qing Lan Project of Jiangsu province to Lian Chen, Shanghai Agricultural Science Committee Youth Talents Development Plan No. 2018 (1–32), Biodiversity Investigation, Observation and Assessment Program of Ministry of Ecology and Environment of China.
Ethics approval and consent to participate
This research was approved by the ethics committee of Nanjing Normal University.
Consent for publication
All authors have no conflict of interest related to this study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Noranked OTUs on genus level using the BLASTn on NCBI nucleotide collection.
Abundance differences in the gut microbiota communities by Kruskal-Wallis H test, among the samples from female, male and juvenile at phylum and genus levels.
The rarefaction curve of Shannon index on OTU level.
LEfSe analysis of gut microbiota predictive COG functions in each group (LDA score>2).
LEfSe analysis of gut microbiota predictive KEGG functional pathways in each group (LDA score>2)
About this article
Cite this article
Chen, L., Li, S., Xiao, Q. et al. Composition and diversity of gut microbiota in Pomacea canaliculata in sexes and between developmental stages. BMC Microbiol 21, 200 (2021). https://doi.org/10.1186/s12866-021-02259-2
- Pomacea canaliculata
- Gut microbiota
- 16S rRNA gene
- High-throughput sequencing