Freshwater gastropods are widely distributed and play an important role in aquatic ecosystems. Symbiotic microorganisms represented by gut microbes can affect the physiological and biochemical activities of their hosts. However, few studies have investigated the response of the gut microbial community of snails to environmental stress. In this study, the dynamics of the gut microbiota of the gastropod Bellamya aeruginosa were tracked to explore their responses in terms of their composition and function to cyanobacterial bloom. Differences in gut microbial community structures during periods of non-cyanobacterial bloom and cyanobacterial bloom were determined. Results showed that the alpha diversity of the gut microbiota exposed to cyanobacterial bloom was lower than that of the gut microbiota exposed to non-cyanobacterial bloom. The main genera differentiating the two periods were Faecalibacterium, Subdoligranulum, Ralstonia, and Pelomonas. Microcystins (MCs) and water temperature (WT) were the primary factors influencing the gut microbial community of B. aeruginosa; between them, the influence of MCs was greater than that of WT. Fourteen pathways (level 2) were notably different between the two periods. The pathways of carbohydrate metabolism, immune system, environmental adaptation, and xenobiotics biodegradation and metabolism in these differential pathways exhibited a strong linear regression relationship with MCs and WT. Changes in the functions of the gut microbiota may help B. aeruginosa meet its immunity and energy needs during cyanobacterial bloom stress. These results provide key information for understanding the response pattern of freshwater snail intestinal flora to cyanobacterial blooms and reveal the underlying environmental adaptation mechanism of gastropods from the perspective of intestinal flora.

Frequent outbreaks of harmful cyanobacterial blooms and the cyanotoxins they produce not only seriously jeopardize the health of freshwater ecosystems but also directly affect the survival of aquatic organisms. In this study, the dynamic characteristics and response patterns of transcriptomes and gut microbiomes in gastropod Bellamya aeruginosa were investigated to explore the underlying response mechanisms to toxic cyanobacterial exposure. The results showed that toxic cyanobacteria exposure induced overall hepatopancreatic transcriptome changes. A total of 2128 differentially expressed genes were identified at different exposure stages, which were mainly related to antioxidation, immunity, and metabolism of energy substances. In the early phase (the first 7 days of exposure), the immune system may notably be the primary means of resistance to toxin stress, and it performs apoptosis to kill damaged cells. In the later phase (the last 7 days of exposure), oxidative stress and the degradation activities of exogenous substances play a dominant role, and nutrient substance metabolism provides energy to the body throughout the process. Microbiomic analysis showed that toxic cyanobacteria increased the diversity of gut microbiota, enhanced interactions between gut microbiota, and altered microbiota function. In addition, the changes in gut microbiota were correlated with the expression levels of antioxidant-, immune-, metabolic-related differentially expressed genes. These results provide a comprehensive understanding of gastropods and intestinal microbiota response to toxic cyanobacterial stress.

Gut microbes play a critical role in helping hosts adapt to external environmental changes and are becoming an important phenotype for evaluating the response of aquatic animals to environmental stresses. However, few studies have reported the role that gut microbes play after the exposure of gastropods to bloom-forming cyanobacteria and toxins. In this study, we investigated the response pattern and potential role of intestinal flora in freshwater gastropod Bellamya aeruginosa when exposed to toxic and non-toxic strains of Microcystis aeruginosa, respectively. Results showed that the composition of the intestinal flora of the toxin-producing cyanobacteria group (T group) changed significantly over time. The concentration of microcystins (MCs) in hepatopancreas tissue decreased from 2.41 ± 0.12 on day 7 to 1.43 ± 0.10 μg·g−1 dry weight on day 14 in the T group. The abundance of cellulase-producing bacteria (Acinetobacter) was significantly higher in the non-toxic cyanobacteria group (NT group) than that in the T group on day 14, whereas the relative abundance of MC-degrading bacteria (Pseudomonas and Ralstonia) was significantly higher in the T group than that in the NT group on day 14. In addition, the co-occurrence networks in the T group were more complex than that in the NT group at day 7 and day 14. Some genera identified as key nodes, such as Acinetobacter, Pseudomonas, and Ralstonia, showed different patterns of variation in the co-occurrence network. Network nodes clustered to Acinetobacter increased in the NT group from day 7 to day 14, whereas the interactions between Pseudomonas and Ralstonia and other bacteria almost changed from positive correlations in the D7T group to negative correlations in the D14T group. These results suggested that these bacteria not only have the ability to improve host resistance to toxic cyanobacterial stress by themselves, but they can also further assist host adaptation to environmental stress by regulating the interaction patterns within the community. This study provides useful information for understanding the role of freshwater gastropod gut flora in response to toxic cyanobacteria and reveals the underlying tolerance mechanisms of B. aeruginosa to toxic cyanobacteria.