Carbohydrate catabolic capability of a Flavobacteriia bacterium isolated from hadal water

Eliminating the ecological hazards of Heterosigma akashiwo bloom by a microbial algicide: removal of nitrite contamination, redirection of carbon flow and restoration of metabolic generalists

Harmful algal blooms (HABs) attracted much attention due to their extensive ecological hazards and the increasing influences on global biogeochemical cycles with the intensification of human impact and global warming. Lysing algal cells with species-specific microbial algicide seemed to be promising to eliminate HABs, but the potential ecotoxicity was rarely studied. In this study, microcosms simulating Heterosigma akashiwo blooms were established to reveal the influences of a microbial algicide from Streptomyces sp. U3 on the biological, physicochemical parameters and bacterial community. The results showed that H. akashiwo bloom accumulated nitrite to a lethal dose, produced bio-labile DOM with widespread influences and enriched pathogenic Coxiella to a high abundance. Lysing H. akashiwo cells by microbial algicide induced a bacterial bloom, eliminated nitrite contamination, enhanced the recalcitrance of DOM, and restored bacterial population from a Gammaproteobacteria-dominant community during bloom back to an Alphaproteobacteria-dominant community similar to the non-bloom seawater. Succession of bacterial genera further suggested that the variation from algal exudates to lysates promoted the restoration of metabolic generalists, which redirected the carbon flow to a less ecologically impactive path. This study revealed the benefits of using microbial algicide to remediate the ecological hazards of HABs, which provided references for future application.

Spatial variation and metabolic diversity of microbial communities in the surface sediments of the Mariana Trench

Mariana Trench represents the deepest and one of least explored biosphere on Earth, and its carbon sources include euphotic sinking, lateral transportation and diffusion from underlying crust, etc. By far the spatial variation of microbial community with associated organic carbon degradation potential in the surface sediments of the Mariana Trench were still largely unknown. Based on the high-throughput 16S rRNA amplicon sequencing, significantly different microbial community structure was overserved between the shallow (<10,000 m) and deep stations (>10,000 m), which could be explained by spatial variation of Chloroflexi, Proteobacteria and Crenarchaeota, with sampling depth and total organic carbon (TOC) content as the environmental driving forces. During the 109-day incubation with Biolog EcoPlate™ microplate, polymers and carbohydrates were preferentially used, followed by amino acids and carboxylic acids, and microbial metabolic diversity was significantly different between the shallow and deep stations. The metabolic diversity of microorganisms at most shallow stations was significantly lower than that at deep stations. This could potentially be attributed the metabolic capabilities of different microbial groups with varied ecological niches, and reflected the initial preference of carbon source by the nature microbes as well. Our study obtained a rough assessment of physiological and taxonomic characteristics of the trench sediment microbial community with polyphasic approaches. Distinct microbial structure and potential carbon metabolic functions in different sampling depths might led to the differentiation of ecological niches, which enable various microorganisms to make full use of the limited resources in the deep sea, and provided a research basis for further exploration of the carbon cycle in different deep-sea regions.

Spatial Variability of Bacterial Community Compositions in the Mariana Trench

Hadal microorganisms play an important role in the biogeochemical processes in marine ecosystems and act as a valuable resource for industrial applications. This paper presents the bacterial community analysis of samples taken from the Challenger Deep within the Mariana Trench, which is the deepest site in the ocean. High-throughput 16S rRNA gene amplicon sequencing was used to reveal that the vertically sampled bacterial populations at eight stations varied at the surface to 10 km depth. The surface water samples harbored a distinct bacterial assemblage, while the mesopelagic and bathyal samples manifested different bacterial community composition, which was not consistent with previous studies. Gammaproteobacteria was the most abundant bacteria in the bathyal and hadal water. The hadal bacterial community consisted mostly of Alteromonadales and Oceanospirillales. The former was widely spread in the water column, which might suggest habitat partitioning at the genus and OTU levels, while the latter might represent hadal-enriched hydrocarbon degraders. The present work complements the current knowledge and understanding of the bathyal and hadal bacterial communities of the Mariana Trench.

Coastal Microbial Communities Disrupted During the 2018 Hurricane Season in Outer Banks, North Carolina

Hurricane frequencies and intensities are expected to increase under warming climate scenarios, increasing potential to disrupt microbial communities from steady-state conditions and alter ecosystem function. This study shows the impact of hurricane season on microbial community dynamics within the barrier island system of Outer Banks, North Carolina. We found that the passage of two sequential energetic hurricanes in 2018 (Florence and Michael) were correlated with shifts in total and active (DNA and RNA) portions of bacterial communities but not in archaeal communities, and within surface waters but not within the sediment. These microbial community shifts were distinct from non-hurricane season conditions, suggesting significant implications for nutrient cycling in nearshore and offshore environments. Hurricane-influenced marine sites in the coastal North Atlantic region had lower microbial community evenness and Shannon diversity, in addition to increased relative abundance of copiotrophic microbes compared to non-hurricane conditions. The abundance of functional genes associated with carbon and nitrogen cycling pathways were also correlated with the storm season, potentially shifting microbial communities at offshore sites from autotroph-dominated to heterotroph-dominated and leading to impacts on local carbon budgets. Understanding the geographic- and system-dependent responses of coastal microbial communities to extreme storm disturbances is critical for predicting impacts to nutrient cycling and ecosystem stability in current and future climate scenarios.

Diversity of epiphytic bacterial communities on male and female Porphyra haitanensis

Purpose

To study the structure of the epiphytic bacterial community of the male and female Porphyra haitanensis, in order to explore the similarities and differences of epiphytic bacterial community structure between dioecious macroalgae.

Methods

Collection of male and female Porphyra haitanensis from the intertidal zone of Niushan Island, Fujian, China. Epiphytic bacteria were collected and studied, and the community composition and diversity of epiphytic bacteria were explored using high-throughput sequencing technology.

Results

There was no significant difference between male and female Porphyra haitanensis on α-diversity and β-diversity. Proteobacteria and Bacteroidetes were the core microbiota in male and female Porphyra haitanensis. Bacteria from the Maribacter (male 14.87%, female 1.66%) and the Tenacibaculum (male 1.44%, female 25.78%) were the most indicative epiphytic bacterial taxa on male and female Porphyra haitanensis.

Conclusions

Sex differences have some influence on the construction of epiphytic bacterial communities in Porphyra haitanensis, but they are not the decisive factors affecting the construction of epiphytic bacterial communities in Porphyra haitanensis.

Phage-encoded ribosomal protein S21 expression is linked to late-stage phage replication

The ribosomal protein S21 (bS21) gene has been detected in diverse viruses with a large range of genome sizes, yet its in situ expression and potential significance have not been investigated. Here, we report five closely related clades of bacteriophages (phages) represented by 47 genomes (8 curated to completion and up to 331 kbp in length) that encode a bS21 gene. The bS21 gene is on the reverse strand within a conserved region that encodes the large terminase, major capsid protein, prohead protease, portal vertex proteins, and some hypothetical proteins. Based on CRISPR spacer targeting, the predominance of bacterial taxonomic affiliations of phage genes with those from Bacteroidetes, and the high sequence similarity of the phage bS21 genes and those from Bacteroidetes classes of Flavobacteriia, Cytophagia and Saprospiria, these phages are predicted to infect diverse Bacteroidetes species that inhabit a range of depths in freshwater lakes. Thus, bS21 phages have the potential to impact microbial community composition and carbon turnover in lake ecosystems. The transcriptionally active bS21-encoding phages were likely in the late stage of replication when collected, as core structural genes and bS21 were highly expressed. Thus, our analyses suggest that the phage bS21, which is involved in translation initiation, substitutes into the Bacteroidetes ribosomes and selects preferentially for phage transcripts during the late-stage replication when large-scale phage protein production is required for assembly of phage particles.

Scientific and technological progress in the microbial exploration of the hadal zone

The hadal zone is the deepest point in the ocean with a depth that exceeds 6000 m. Exploration of the biological communities in hadal zone began in the 1950s (the first wave of hadal exploration) and substantial advances have been made since the turn of the twenty-first century (the second wave of hadal exploration), resulting in a focus on the hadal sphere as a research hotspot because of its unique physical and chemical conditions. A variety of prokaryotes are found in the hadal zone. The mechanisms used by these prokaryotes to manage the high hydrostatic pressures and acquire energy from the environment are of substantial interest. Moreover, the symbioses between microbes and hadal animals have barely been studied. In addition, equipment has been developed that can now mimic hadal environments in the laboratory and allow cultivation of microbes under simulated in situ pressure. This review provides a brief summary of recent progress in the mechanisms by which microbes adapt to high hydrostatic pressures, manage limited energy resources and coexist with animals in the hadal zone, as well as technical developments in the exploration of hadal microbial life.

Description of Aureibaculum luteum sp. nov. and Aureibaculum flavum sp. nov. isolated from Antarctic intertidal sediments

Two Gram-stain-negative, aerobic, non-motile, and rod-shaped bacterial strains, designated SM1352^T and A20^T, were isolated from intertidal sediments collected from King George Island, Antarctic. They shared 99.8% 16S rRNA gene sequence similarity with each other and had the highest sequence similarity of 98.1% to type strain of Aureibaculum marinum but < 93.4% sequence similarity to those of other known bacterial species. The genomes of strains SM1352^T and A20^T consisted of 5,108,092 bp and 4,772,071 bp, respectively, with the G + C contents both being 32.0%. They respectively encoded 4360 (including 37 tRNAs and 6 rRNAs) and 4032 (including 36 tRNAs and 5 rRNAs) genes. In the phylogenetic trees based on 16S rRNA gene and single-copy orthologous clusters (OCs), both strains clustered with Aureibaculum marinum and together formed a separate branch within the family Flavobacteriaceae. The ANI and DDH values between the two strains and Aureibaculum marinum BH-SD17^T were all below the thresholds for species delineation. The major cellular fatty acids (> 10%) of the two strains included iso-C_15:0, iso-C_15:1 G, iso-C_17:0 3-OH. Their polar lipids predominantly included phosphatidylethanolamine, one unidentified aminophospholipid, one unidentified aminolipid, and two unidentified lipids. Genomic comparison revealed that both strains possessed much more glycoside hydrolases and sulfatase-rich polysaccharide utilization loci (PULs) than Aureibaculum marinum BH-SD17^T. Based on the above polyphasic evidences, strains SM1352^T and A20^T represent two novel species within the genus Aureibaculum, for which the names Aureibaculum luteum sp. nov. and Aureibaculum flavum sp. nov. are proposed. The type strains are SM1352^T (= CCTCC AB 2014243 ^T = JCM 30335 ^T) and A20^T (= CCTCC AB 2020370 ^T = KCTC 82503 ^T), respectively.

Ammonia stress disrupts intestinal microbial community and amino acid metabolism of juvenile yellow catfish (Pelteobagrus fulvidraco)

Ammonia has adverse effects on aquatic animals, which is also widely distributed in natural aquatic environments and intensive aquaculture systems. The intestine is a primary defensive line for aquatic animals, the accumulation of ammonia in the aquatic environment can cause irreversible damage to intestinal function. In this study, we investigated the effects of acute ammonia stress on the reaction characteristics of digestive function, amino acid metabolism, and the variation in the intestinal microbiota of juvenile yellow catfish (Pelteobagrus fulvidraco). Thus, the yellow catfish was placed in water with the addition of ammonia at 0 (control), 14.6, and 146 mg/L total ammonia nitrogen for 96-h. The present study observed that ammonia accumulated in the intestine and muscle (ammonia contents in the intestine and muscle increased) and induced the activities of protein digestive enzymes dysfunction (pepsin increased while trypsin decreased). Ammonia stress changed various amino acids composition (proline, arginine, lysine, histidine, phenylalanine, tyrosine, leucine, isoleucine, valine, alanine, glutamic acid, tyrosine, and aspartic acid contents were increased in muscle) and increased the activities of alanine aminotransferase and aspartate aminotransferase in muscle. Furthermore, through 16 S rRNA gene analysis, ammonia stress-induced reduction in diversity, richness, and evenness and structure of microbiota alteration in the intestine. At the phylum level, the abundance of Fusobacteria increased while Firmicutes and Actinobacteria decreased significantly. At the genus level, the abundance of beneficial microbiota Cetobacterium significantly increased after ammonia stress. In conclusion, activation of amino acid synthesis in muscle may be involved in ammonia detoxification after severe ammonia stress. The accumulation of ammonia can disrupt the intestinal digestive function and intestinal microbiota community. The Cetobacterium may be a new potential positive factor in the resistance of ammonia toxicity.

Effects of salinity shock on simultaneous nitrification and denitrification by a membrane bioreactor: Performance, sludge activity, and functional microflora

Physical and chemical treatments of Tungsten smelting wastewater, with high salt content and low C/N ratio, are often tedious. As a solution, this study suggested a simultaneous nitrification and denitrification membrane bioreactor (SND-MBR) for salinity gradient domestication. During the salinity acclimation period, we observed 20% and 11% removal of NH₄⁺-N and Chemical Oxygen Demand (COD), respectively. However, the SND efficiency reached 95.55% after stable operation at 3.0% salinity. Through stoichiometric and kinetic analyses, we confirmed that increased salinity significantly inhibited electron transport system activity, nitrification, and denitrification, evidenced by the extremely low ammonia monooxygenase and nitrite reductase activities. Further high-throughput sequencing showed that Nitrosomonas dominated the functional microbial flora succession and denitrification in high salinity environments. In comparison with a control, the Kyoto Encyclopedia of Genes and Genomes analysis showed that wastewater salinity weakened the functional gene level of MBR microbial flora, and the enzyme key to the assimilation nitrate reduction changed from nitrate reductase to assimilation nitrate reductase.

Comparative genomic analysis of obligately piezophilic Moritella yayanosii DB21MT-5 reveals bacterial adaptation to the Challenger Deep, Mariana Trench

Hadal trenches are the deepest but underexplored ecosystems on the Earth. Inhabiting the trench bottom is a group of micro-organisms termed obligate piezophiles that grow exclusively under high hydrostatic pressures (HHP). To reveal the genetic and physiological characteristics of their peculiar lifestyles and microbial adaptation to extreme high pressures, we sequenced the complete genome of the obligately piezophilic bacterium Moritella yayanosii DB21MT-5 isolated from the deepest oceanic sediment at the Challenger Deep, Mariana Trench. Through comparative analysis against pressure sensitive and deep-sea piezophilic Moritella strains, we identified over a hundred genes that present exclusively in hadal strain DB21MT-5. The hadal strain encodes fewer signal transduction proteins and secreted polysaccharases, but has more abundant metal ion transporters and the potential to utilize plant-derived saccharides. Instead of producing osmolyte betaine from choline as other Moritella strains, strain DB21MT-5 ferments on choline within a dedicated bacterial microcompartment organelle. Furthermore, the defence systems possessed by DB21MT-5 are distinct from other Moritella strains but resemble those in obligate piezophiles obtained from the same geographical setting. Collectively, the intensive comparative genomic analysis of an obligately piezophilic strain Moritella yayanosii DB21MT-5 demonstrates a depth-dependent distribution of energy metabolic pathways, compartmentalization of important metabolism and use of distinct defence systems, which likely contribute to microbial adaptation to the bottom of hadal trench.

Enhanced treatment of shale gas fracturing waste fluid through plant-microbial synergism

Cost-efficient and environmentally friendly treatment of hydraulic fracturing effluents is of great significance for the sustainable development of shale gas exploration. We investigated the synergistic effects of plant-microbial treatment of shale gas fracturing waste fluid. The results showed that illumination wavelength and temperature are direct drivers for microbial treatment effects of COD_Cr and BOD₅, while exhibit little effects on nitrogen compounds, TDS, EC, and SS removals as well as microbial species and composition. Plant-microbial synergism could significantly enhance the removal of pollutants compared with removal efficiency without plant enhancement. Additionally, the relative abundance and structure of microorganisms in the hydraulic fracturing effluents greatly varied with the illumination wavelength and temperature under plant-microbial synergism. 201.24 g water dropwort and 435 mg/L activated sludge with illumination of 450-495 nm (blue) at 25 °C was proved as the best treatment condition for shale gas fracturing waste fluid samples, which showed the highest removal efficiency of pollutants and the lowest algal toxicity in treated hydraulic fracturing effluents. The microbial community composition (36.73% Flavobacteriia, 25.01% Gammaproteobacteria, 18.55% Bacteroidia, 9.3% Alphaproteobacteria, 4.1% Cytophagia, and 2.83% Clostridia) was also significantly different from other treatments. The results provide a potential technical solution for improved treatment of shale gas hydraulic fracturing effluents.

Performance and enhancement mechanism of corncob guiding chromium (VI) bioreduction

Chromium-contaminated groundwater has drawn extensive attention due to its high toxicity and wide application. Although bioremediation is considered to be an effective approach for Cr(VI) removal, a better method is still urgently needed. In this study, corncob-guided Cr(VI) reduction achieved the highest removal efficiency due to the highest amount of total carbon and available carbon emissions. After verifying the sustainability and operational feasibility of this approach, the broad-spectrum applicability of corncob to guide Cr(VI) bioreduction was further explored under various operating conditions. In addition, it suggested that the carrier effect, nutrient element release and electron shuttle effect were the main mechanisms enhancing the reduction process, with approximate contribution rates of 12.5%, 7.5% and 75%, respectively. Microbiological analysis demonstrated that the addition of solid-phase carbon sources increased the abundance of microbes related to carbon metabolism and promoted the expression of glycolytic metabolic pathway. Furthermore, the addition of corncob led to an elevation of expression level of the electron transport pathway, which is consistent with the function of the electron shuttle. This study provides theoretical and technical support for the practical application of corncob-mediated Cr(VI) bioreduction.

Improvement in nitrogen removal and changes in community structure in a sequencing batch reactor bioaugmented with P. stutzeri strain XL-2

The aim of this work was to study the bioaugmentation of P. stutzeri strain XL-2 in activated sludge to improve nitrogn removal from wastewater with the guide of growth kinetics. When 4250 mg/L COD and 80 mg/L NH₄⁺-N were applied, the TN removal efficiency in a bioaugmented sequencing batch reactor (SBR_XL) achieved 95%, while that in the control reactor (SBR_C) without strain XL-2 was only 84% (P < 0.05). The microbial community analysis demonstrated that strain XL-2 was successfully bioaugmented in SBR_XL, and increasing influent COD concentration promoted its abundance. Influent COD concentration played a dominant role in affecting community structure, while the bioaugmentation of strain XL-2 had much less impact on the community structure. Combined with principal coordinates analysis, redundancy analysis and FAPROTAX, the improvement of TN removal was mainly achieved by the bioaugmentation of strain XL-2, which played a major role in promoting aerobic denitrification.

Effects of Microcystis aeruginosa and microcystin-LR on intestinal histology, immune response, and microbial community in Litopenaeus vannamei

Microcystis aeruginosa (MA) is a primary hazardous cyanobacteria species in aquatic ecosystems that can produce microcystin-LR (MC-LR), which harms aquatic animals. The intestine is an important target tissue for MA and MC-LR. In this study, we investigated the effects of MA and MC-LR exposure on the intestinal microbiota variation and immune responses of Litopenaeus vannamei. Shrimp were experimentally exposed to MA and MC-LR for 72 h. The results showed that both MA and MC-LR exposure caused marked histological variation and apoptosis characteristics and increased oxidative stress in the intestine. Furthermore, the relative expression levels of antimicrobial peptide genes (ALF, Crus, Pen-3) decreased, while those of pro-inflammatory cytokines (MyD88, Rel, TNF-a), a pattern-recognition receptor (TLR4) and a mediator of apoptosis (Casp-3) increased. MA and MC-LR exposure also caused intestinal microbiota variation, including decreasing microbial diversity and disturbing microbial composition. Specifically, the relative abundance of Proteobacteria decreased in the two stress groups; that of Bacteroidetes decreased in the MA group but increased in the MC-LR group, while Tenericutes varied inversely with Bacteroidetes. Our results indicate that MA and MC-LR exposure causes intestinal histopathological and microbiota variations and induces oxidative stress and immune responses in L. vannamei. In conclusion, this study reveals the negative effects of MA and MC-LR on the intestinal health of shrimp, which should be considered in aquaculture.

Insights into the Vertical Stratification of Microbial Ecological Roles across the Deepest Seawater Column on Earth

The Earth's oceans are a huge body of water with physicochemical properties and microbial community profiles that change with depth, which in turn influences their biogeochemical cycling potential. The differences between microbial communities and their functional potential in surface to hadopelagic water samples are only beginning to be explored. Here, we used metagenomics to investigate the microbial communities and their potential to drive biogeochemical cycling in seven different water layers down the vertical profile of the Challenger Deep (0-10,500 m) in the Mariana Trench, the deepest natural point in the Earth's oceans. We recovered 726 metagenome-assembled genomes (MAGs) affiliated to 27 phyla. Overall, biodiversity increased in line with increased depth. In addition, the genome size of MAGs at ≥4000 m layers was slightly larger compared to those at 0-2000 m. As expected, surface waters were the main source of primary production, predominantly from Cyanobacteria. Intriguingly, microbes conducting an unusual form of nitrogen metabolism were identified in the deepest waters (>10,000 m), as demonstrated by an enrichment of genes encoding proteins involved in dissimilatory nitrate to ammonia conversion (DNRA), nitrogen fixation and urea transport. These likely facilitate the survival of ammonia-oxidizing archaea α lineage, which are typically present in environments with a high ammonia concentration. In addition, the microbial potential for oxidative phosphorylation and the glyoxylate shunt was enhanced in >10,000 m waters. This study provides novel insights into how microbial communities and their genetic potential for biogeochemical cycling differs through the Challenger deep water column, and into the unique adaptive lifestyle of microbes in the Earth's deepest seawater.

Comparative genomic analysis of Flavobacteriaceae: insights into carbohydrate metabolism, gliding motility and secondary metabolite biosynthesis

BACKGROUND

Members of the bacterial family Flavobacteriaceae are widely distributed in the marine environment and often found associated with algae, fish, detritus or marine invertebrates. Yet, little is known about the characteristics that drive their ubiquity in diverse ecological niches. Here, we provide an overview of functional traits common to taxonomically diverse members of the family Flavobacteriaceae from different environmental sources, with a focus on the Marine clade. We include seven newly sequenced marine sponge-derived strains that were also tested for gliding motility and antimicrobial activity.

RESULTS

Comparative genomics revealed that genome similarities appeared to be correlated to 16S rRNA gene- and genome-based phylogeny, while differences were mostly associated with nutrient acquisition, such as carbohydrate metabolism and gliding motility. The high frequency and diversity of genes encoding polymer-degrading enzymes, often arranged in polysaccharide utilization loci (PULs), support the capacity of marine Flavobacteriaceae to utilize diverse carbon sources. Homologs of gliding proteins were widespread among all studied Flavobacteriaceae in contrast to members of other phyla, highlighting the particular presence of this feature within the Bacteroidetes. Notably, not all bacteria predicted to glide formed spreading colonies. Genome mining uncovered a diverse secondary metabolite biosynthesis arsenal of Flavobacteriaceae with high prevalence of gene clusters encoding pathways for the production of antimicrobial, antioxidant and cytotoxic compounds. Antimicrobial activity tests showed, however, that the phenotype differed from the genome-derived predictions for the seven tested strains.

CONCLUSIONS

Our study elucidates the functional repertoire of marine Flavobacteriaceae and highlights the need to combine genomic and experimental data while using the appropriate stimuli to unlock their uncharted metabolic potential.

Evidence in the Japan Sea of microdolomite mineralization within gas hydrate microbiomes

Over the past 15 years, massive gas hydrate deposits have been studied extensively in Joetsu Basin, Japan Sea, where they are associated primarily with active gas chimney structures. Our research documents the discovery of spheroidal microdolomite aggregates found in association with other impurities inside of these massive gas hydrates. The microdolomites are often conjoined and show dark internal cores occasionally hosting saline fluid inclusions. Bacteroidetes sp. are concentrated on the inner rims of microdolomite grains, where they degrade complex petroleum-macromolecules present as an impurity within yellow methane hydrate. These oils show increasing biodegradation with depth which is consistent with the microbial activity of Bacteroidetes. Further investigation of these microdolomites and their contents can potentially yield insight into the dynamics and microbial ecology of other hydrate localities. If microdolomites are indeed found to be ubiquitous in both present and fossil hydrate settings, the materials preserved within may provide valuable insights into an unusual microhabitat which could have once fostered ancient life.

Transcriptomic and microbiota response on Litopenaeus vannamei intestine subjected to acute sulfide exposure

Harmful effects of water pollutants are myriad. Sulfide from water bodies affects the aquatic animals. Intestine barrier function serves as the front-line of animals defense. Our previous study confirmed the toxic effect of sulfide on intestine immune response of Litopenaeus vannamei, but the underlying mechanisms remained elusive. Therefore, in this study, we investigated the transcriptomic and microbiota responses of the L. vannamei intestine subjected to acute sulfide exposure. Sulfide decreased bacterial richness and altered the intestine microbiota composition. Specifically, sulfide increased the abundances of Bacteroidetes and Actinobacteria, but decreased the abundance of Proteobacteria. At the genus level, sulfide increased typical cellulolytic characteristics bacteria, such as Formosa, Sphingomonas, and Demequina. RNA-seq analysis identified differential expression of 1799 genes (701 up-regulated and 1098 down-regulated) were grouped into 267 pathways. The most enriched pathway 'amoebiasis' was related to the intestine mucus homeostasis. A number of immune-related genes associated with antimicrobial, antioxidant, pathogen attachment and recognition, and apoptosis processes in contrasting accessions; they were correlated with the abundance of intestine bacterial at the phylum level. This study provides an insight into the mechanisms associated with molecular and microbiota response and processes involved in adaptation strategies towards sulfide stress.