Insights into the identification of bitter peptides from Jinhua ham and its taste mechanism by molecular docking and transcriptomics analysis.

In order to identify the peptides responsible for bitter defects and to understand the mechanism of bitterness in dry-cured ham, the peptides were identified by LC-MS/MS, and the interaction between bitter peptides and receptor proteins were evaluated by molecular docking and molecular dynamics simulation; the signal transduction mechanism of bitter peptides was investigated using the model of HEK-293T cells by calcium imaging and transcriptomics analysis. The results of LC-MS/MS showed that 11 peptides were identified from the high bitterness fraction of defective ham; peptides PKAPPAK, VTDTTR and YIIEK derived from titin showed the highest bitterness values compared with other peptides. The results of molecular docking showed that lower CDOCKER energy was observed in the interaction between these peptides and hT2R16 in comparison with these receptors of hT2R1, hT2R4, hT2R5, hT2R8 and hT2R14, and the interaction of hT2R16 and peptides was stabilized by hydrophobic interaction and hydrogen bond. The average RMSF values of VTDTTR were higher than that of YIIEK and PKAPPAK, while EC50 values of VTDTTR were lower compared with PKAPPAK and YIIEK. Transcriptomics analysis showed that 529 differentially expressed genes were identified in HEK-293T cells during the stimulating by VTDTTR and were mainly enriched into neuroactive ligand-receptor interaction, MAPK pathway, cAMP pathway and calcium signaling pathway, which were mainly responsible for the bitter signal transduction of VTDTTR. These results could provide evidence for understanding the bitter defects of dry-cured ham and the taste mechanism of bitter peptide.

Genomic and Proteomic Analyses of Extracellular Products Reveal Major Virulence Factors Likely Accounting for Differences in Pathogenicity to Bivalves between Vibrio mediterranei Strains.

Vibrio mediterranei, a bacterial pathogen of bivalves, has exhibited strain-dependent virulence. The mechanisms behind the variations in bivalve pathogenicity between V. mediterranei strains have remained unclear. However, a preliminary analysis of the extracellular product (ECP) proteomes has revealed differences in protein compositions between low- and high-virulence strains; in addition to 1265 shared proteins, 127 proteins have been identified to be specific to one low-virulence strain and 95 proteins to be specific to two high-virulence strains. We further studied the ECP proteins of the three V. mediterranei strains from functional perspectives using integrated genomics and proteomics approaches. The results showed that lipid metabolism, transporter activity and membrane transporter pathways were more enriched in the ECPs of the two high-virulence strains than in those of the low-virulence strain. Additionally, 73 of the 95 high-virulence strain-specific proteins were found to have coding genes in the genome but were not expressed in the low-virulence strain. Moreover, comparisons with known virulence factors in the Virulence Factor Database (VFDB) and the Pathogen-Host Interactions Database (PHI-base) allowed us to predict more than 10 virulence factors in the categories of antimicrobial activity/competitive advantage, the effector delivery system and immune modulation, and the high-virulence strain-specific ECP proteins consisted of a greater percentage of known virulence factors than the low-virulence strain. Particularly, two virulence factors, MtrC and KatG, were identified in the ECPs of the two high-virulence strains but not in those of the low-virulence strain. Most coding genes of the ECP proteins including known virulence factors were identified on chromosome 1 of V. mediterranei. Our findings indicate that variations in virulence factor composition in the bacterial ECPs may partially account for the differences in the bivalve pathogenicity between V. mediterranei strains.

Acute Salinity Stress Disrupts Gut Microbiota Homeostasis and Reduces Network Connectivity and Cooperation in Razor Clam Sinonovacula constricta.

Accumulating evidence demonstrates that it is of great importance to maintain a stable and functional gut microbial community for host's growth and health. However, gut microenvironment is constantly affected by diverse environmental factors. Salinity can cause stress, including hypersaline or hyposaline stress to aquatic species, thereby affecting their growth conditions. Razor clam (Sinonovacula constricta), an economically important bivalve species, inhabits in intertidal and estuarine zones and constantly experiences salinity stress. Yet little is known about how and to what extent clam gut microbiota is affected by salinity stress, while this knowledge is fundamental for clam aquaculture health management. To address this concern, this study compared the temporal differences of gut bacterial signatures and community assembly of S. constricta under normal salinity (NS), low salinity (LS), and high salinity (HS) conditions. Acute salinity stress affected the compositions, structures, and functional potentials of clam gut microbial community, of which salinity stress, hours post stress, and their interaction respectively constrained 7.6%, 16.4%, and 7.9% of community variation. Phylogenetic bin-based null model result revealed that the gut bacterial assembly of three salinity groups seemed to be largely driven by stochastic processes. Network analysis indicated that gut bacterial interspecies interaction exhibited less connected and lower cooperative activity under the conditions of LS and HS compared with NS. Notably, some pathogenic bacteria, including Vibrio and Pseudoalteromonas, were identified as keystone taxa of gut microbial networks in LS and HS groups. Above findings suggest that the clams under LS and HS conditions might be at a higher risk of developing disease. Our findings enhance the mechanism understanding of gut microbial assembly in S. constricta under abiotic factor challenge, which has important implications for clam health control from a microbial ecological perspective.

Expression and Functional Analysis of AMT1 Gene Responding to High Ammonia Stress in Razor Clam (Sinonovacula constricta).

Ammonium transporter 1 (AMT1), a member of ammonia (NH3/NH4+) transport proteins, has been found to have ammonia transport activity in plants and microorganisms. However, the functional characteristics and molecular mechanisms of AMT1 in mollusks remain unclear. The razor clam (Sinonovacula constricta) is a suitable model species to explore the molecular mechanism of ammonia excretion because of the high concentration of ambient ammonia it is exposed to in the clam-fish-shrimp polyculture system. Here, the expression of AMT1 in S. constricta (Sc-AMT1) in response to high ammonia (12.85 mmol/L NH4Cl) stress was identified by real-time quantitative PCR (qRT-PCR), Western blotting, RNA interference, and immunofluorescence analysis. Additionally, the association between the SNP_g.15211125A > T linked with Sc-AMT1 and ammonia tolerance was validated by kompetitive allele-specific PCR (KASP). A significant upregulated expression of Sc-AMT1 was observed during ammonia exposure, and Sc-AMT1 was found to be localized in the flat cells of gill. Moreover, the interference with Sc-AMT1 significantly upregulated the hemolymph ammonia levels, accompanied by the increased mRNA expression of Rhesus glycoprotein (Rh). Taken together, our findings imply that AMT1 may be a primary contributor to ammonia excretion in S. constricta, which is the basis of their ability to inhabit benthic water with high ammonia levels.

Characterization of Vibrio mediterranei Isolates as Causative Agents of Vibriosis in Marine Bivalves.

Marine bivalves include species important globally for aquaculture and estuary ecology. However, epizootics of variable etiologies often pose a threat to the marine fishery industry and ecosystem by causing significant mortalities in related species. One of such diseases is larval vibriosis caused by bacteria of the genus Vibrio, which frequently occurs and causes mass mortalities in bivalve larvae and juveniles in hatcheries. During a mass mortality of razor clam, Sinonovacula constricta, juveniles in a shellfish hatchery in 2019, Vibrio mediterranei was identified as a dominant bacterial species in diseased animals and their rearing water. In this study, we selected and characterized 11 V. mediterranei isolates and studied their pathogenicity to the larvae and juveniles of S. constricta and Crossostrea sikamea. We found that V. mediterranei isolates showed various degrees of pathogenicity to the experimental animals by immersion. Injection of the extracellular products (ECPs) of the strains into clam juveniles resulted in similar pathogenicity with strain immersion. Furthermore, the measurements of enzyme activity suggested the existence of virulence factors in the ECPs of disease-causing V. mediterranei strains. Additionally, proteomic analysis revealed that more than 700 differentially expressed proteins were detected in the ECPs among V. mediterranei strains with different levels of virulence, and the higher expressed proteins in the ECPs of highly virulent strains were involved mainly in the virulence-related pathways. This research represented the first characterization of the V. mediterranei strains as causative agents for larval bivalve vibriosis. The mechanisms underlying the pathogenicity and related strain variability are under further study. IMPORTANCE In the marine environment, Vibrio members have a significant impact on aquatic organisms. Larval vibriosis, caused by bacteria of the genus Vibrio, often poses a threat to the marine fishery industry and ecosystem by causing the mortality of bivalves. However, the emerging pathogens of larval vibriosis in bivalves have not been explored fully. Vibrio mediterranei, the dominant bacterium isolated from moribund clam juveniles in a mortality event, may be responsible for the massive mortality of bivalve juveniles and vibriosis occurrence. Thus, it is necessary to study the pathogenic mechanisms of V. mediterranei to bivalve larvae. We found that V. mediterranei was the pathogen of larval bivalve vibriosis, and its extracellular products contributed a critical role for virulence in juveniles. This research is the first report of V. mediterranei as a causative agent for vibriosis in bivalve juveniles. Our results provide valuable information for understanding the pathogenic mechanism of V. mediterranei to bivalve larvae.

Changes in Bacterial Communities of Kumamoto Oyster Larvae During Their Early Development and Following Vibrio Infection Resulting in a Mass Mortality Event.

Vibrio and Ostreid herpesvirus 1 are responsible for mass mortalities of oyster larvae in hatcheries. Relevant works have focused on their relationships with the disease when larval mortality occurs. On the contrary, little is known about how the resident microbiota in oyster larvae responds to Vibrio-infected disease causing mortality as the disease progressed, whereas this knowledge is fundamental to unveil the etiology of the disease. Here, we analyzed the temporal succession of the microbiome of Kumamoto oyster (Crassostrea sikamea) larvae during their early development, accompanied by a Vibrio-caused mortality event that occurred at the post D-stage of larval development in a shellfish hatchery in Ningbo, China, on June 2020. The main causative agent of larval mortality was attributable to Vibrio infection, which was confirmed by linearly increased Vibrio abundance over disease progression. Larval bacterial communities dramatically changed over host development and disease progression, as highlighted by reduced α-diversity and less diverse core taxa when the disease occurred. Null model and phylogenetic-based mean nearest taxon distance analyses showed that the relative importance of deterministic processes governing larval bacterial assembly initially increased over host development, whereas this dominance was depleted over disease progression. Furthermore, we screened the disease-discriminatory taxa with a significant change in their relative abundances, which could be indicative of disease progression. In addition, network analysis revealed that disease occurrence remodeled the co-occurrence patterns and niche characteristics of larval microbiota. Our findings demonstrate that the dysbiosis of resident bacterial communities and the shift of microecological mechanisms in the larval microbiome may contribute to mortality during oyster early development.

PDIA6 involves the thermal stress response of razor clam, Sinonovacula constricta.

Protein disulfide isomerases A6 (PDIA6), an oxidoreductase and isomerase, catalyzes the oxidation reduction and isomerization of disulfide bonds, and serves as molecular chaperone to prevent the buildup of misfolded proteins under various environmental insults. However, the role of PDIA6 in mollusks remains largely obscure, although its multifunctional protein has been reported in other species under adverse conditions. To fill this gap, we identified PDIA6 from the razor clam Sinonovacula constricta (ScPDIA6) and investigated its expression patterns in response to thermal stress. Tissue distribution showed that the mRNA transcript of ScPDIA6 was ubiquitously expressed in nine tested tissues. Temporal expression profiles by qPCR revealed that ScPDIA6 in the gill and mantle was significantly increased by hyper-thermic treatment. Further, Western blot and immunofluorescence indicated that ScPDIA6 was significantly upregulated by thermal treatment at the protein level. Additionally, the survival test demonstrated that the viability of E. coli cells expressing recombinant ScPDIA6 protein increased at 42 °C compared with empty vector. Overall, these findings suggested that ScPDIA6 may play a pivotal role in counteracting thermal stress. This study will provide valuable reference data resource for understanding the potential role of PDIA6 in mollusks.

Bacterial Community Dynamics in Kumamoto Oyster Crassostrea sikamea Hatchery During Larval Development.

Increasing evidence indicates that microbes colonized in early life stages have a long-term effect on animal wellbeing in later life stages. Related research is still limited in aquatic animals, particularly in bivalve mollusks. In this study, we analyzed the dynamics of the bacterial composition of the pelagic larval stages (fertilized egg, trochophore, D-stage, veliger, and pediveliger) and the sessile postlarval stage (spat) of Kumamoto oyster (Crassostrea sikamea) and their relationships with the rearing water bacterioplankton in a hatchery by using Illumina sequencing of bacterial 16S rRNA gene. Both bacterioplankton and larval bacterial communities changed greatly over larval development, and the two communities remarkably differed (r = 0.956, P < 0.001), as highlighted by the differences in the dominant taxa and bacterial diversity. Ecological processes of larval bacterial communities were measured by abundance-unweighted and abundance-weighted standardized effect sizes of the mean nearest taxon distance (ses.MNTD). The unweighted ses.MNTD analysis revealed that the deterministic process constrained the larval bacterial assembly, whereas the weighted ses.MNTD analysis showed that larval bacterial composition was initially governed by stochasticity and then gradually by determinism in the later stages. SourceTracker analysis revealed that the larval bacteria were primarily derived from an internal source, mainly from larvae at the present stage. Additionally, the abundances of larval bacterial-mediated functional pathways that were involved in the amino acid, energy, lipid and carbohydrate metabolisms significantly altered with the larval development. These findings suggest that bacteria assemble into distinct communities in larvae and rearing water in the hatchery system, and the dynamics of bacterial community composition in larvae is likely associated with larval developmental stages.

Characteristics of glutathione peroxidase gene and its responses to ammonia-N stress in razor clam Sinonovacula constricta.

Glutathione peroxidase (GPX) is a crucial enzyme in the antioxidant defense system. However, the previous studies on the structure and functions of mollusk GPX genes are still very limited. Here, we investigated the GPX gene from Sinonovacula constricta (Sc-GPX), and its expression profiles, protein localization, gene function and association with ammonia tolerance. The full length of sequence of Sc-GPX was 1781 bp, containing an open reading frame (ORF) of 588 bp encoding 195 amino acids. Quantitative expression of seven adult tissues showed that Sc-GPX was most abundant in hepatopancreas, followed by gills. Furthermore, the enzyme activity of Sc-GPX in hepatopancreas increased significantly under different ammonia concentrations (100, 140, and 180 mg/L) (P < 0.01). Further, we explored the mRNA expression level, histological structure and histo-cellular localization in gills and hepatopancreas of Sc-GPX under 140 mg/L ammonia stress. The mRNA expression level in gills and hepatopancreas of Sc-GPX increased significantly (P < 0.05) and immunohistochemistry results suggested that the columnar cells of gills filaments and the endothelial cells of hepatopancreas were the major sites for the action of Sc-GPX protein. In addition, we performed western blotting (WB), RNA interference (RNAi) and single nucleotide polymorphisms (SNPs) in the hepatopancreas of Sc-GPX under ammonia stress (140 mg/L). WB results indicated that the protein expression of Sc-GPX increased significantly (P < 0.01) after ammonia challenge. In addition, expression of Sc-GPX mRNA were significantly downregulated at 24 and 48 h after RNAi (P < 0.01). The association analysis between ammonia-tolerance group and control group identified six SNPs in coding sequence (CDS) of Sc-GPX from 449 individuals. Among them, c.162A > C was missense mutation, which lead to the amino acid change from Lys to Asn. These findings revealed that Sc-GPX may play a critical role in clam ammonia detoxification.

Expression and functional characterization of peptidoglycan recognition protein-S6 involved in antibacterial responses in the razor clam Sinonovacula constricta.

It has been recognized that peptidoglycan recognition proteins (PGRPs), structurally conserved molecules, play crucial roles in the innate immunity of invertebrate. However, few studies have been taken to explore their potential functions. In this study, a novel PGRP from the razor clam Sinonovacula constrict designated as ScPGRP-S6 was identified and characterized. The open reading frame (ORF) of ScPGRP-S6 was 666 bp in length, encoding a protein of 221 amino acid with a signal peptide (1-30) and a typical PGRP domain (39-187). The sequence alignment combined with phylogenetic analysis collectively confirmed that ScPGRP-S6 was a novel member belonging to PGRP-S family. The mRNA transcript of ScPGRP-S6 in the hepatopancreases was significantly up-regulated after peptidoglycan (PGN) stimulation, while it was moderately up-regulated after lipopolysaccharide (LPS) stimulation. The result of immunofluorescence detection demonstrated that the positive signal enhanced obviously after Vibrio parahaemolyticus challenge. Notably, the recombinant protein of ScPGRP-S6 (designed as rScPGRP-S6) exhibited high agglutination activity towards V. parahaemolyticus but weak to Staphylococcus aureus. Furthermore, rScPGRP-S6 showed strong amidase and antibacterial activity in the presence of Zn2+. Collectively, our results manifested that ScPGRP-S6 could act as a scavenger in the innate immune response of S. constricta.

Effect of Rhizophora apiculata plantation for improving water quality, growth, and health of mud crab.

A deteriorated water quality is closely associated with disease outbreaks in aquaculture, where microorganisms play indispensable roles in improving water quality and aquatic animals' health. Mangrove is known to be a natural water quality filter and microbiological buffer of pathogen and prebiotics. However, it is unclear how and to what extent Rhizophora apiculata plantation is of benefits to the gut microbiota and growth over mud crab (Scylla paramamosain) aging. To address these concerns, we explored the bacterial communities in mud crab gut and rearing water at 45, 114, and 132 days after incubation, roughly corresponding to juvenile, pre-adult, and adult stages of mud crab. Results showed that 1-year R. apiculata plantation slightly increased the body weight of mud crab and improved water quality to a certain extent. Both bacterioplankton and gut bacterial communities were highly temporal dynamic, while the two communities were significantly distinct (ANOSIM r = 0.90, P = 0.0001). Relative abundances of dominant taxa in water and gut significantly varied between the plantation and the control conditions over mud crab aging. R. apiculata plantation promoted the stability of gut microbiota, as evidenced by more diverse core species. Furthermore, R. apiculata plantation led to the dominance of Verrucomicrobiae species in water and probiotic Bacteroidetes and Lactobacillales taxa in gut. A structural equation model revealed that water variables directly constrained gut microbiota, which in turn affected the body weight of mud crab (r = 0.52, P < 0.001). In addition, functional pathways facilitating immunity and lipid metabolism significantly increased in mud crab gut under the plantation, while those involved in infectious diseases exhibited the opposing trend. These findings greatly expand our understanding of the R. apiculata plantation effects on water quality, gut microbiota, and growth feature of mud crab. Overall, R. apiculata plantation is beneficial for mud crab growth and health. KEY POINTS: • A short-term R. apiculata plantation could potentially improve water quality. • Bacterioplankton is more sensitive than mud crab gut microbiota in response to R. apiculata plantation. • R. apiculata plantation enhances mud crab resistance against pathogen invasion. • R. apiculata plantation alters mud crab gut microbiota, which in turn promotes their body weight.

Concept of microbial gatekeepers: Positive guys?

Microorganisms contribute diverse and fundamental roles in biogeochemical processes. In a given microbial community, individuals interact with one another to form complex regulatory networks in which gatekeepers contribute disproportional roles in sustaining stability, dynamics, and function. Owing to the ecological and functional importance of microbial gatekeeper, this review provides an overview on its history, identification, roles, application in biological sciences, and clinical diagnostics. The roles of microbial gatekeepers can be beneficial or detrimental, which depends on our purpose. As the field is rather new, some limitations are raised, and further efforts devoted to solving these concerns are proposed. Collectively, gatekeepers provide promising targets for sustaining and re-establishing a desired microbial community. However, substantial obstacles, such as factors governing gatekeepers, must be overcome to manipulate gatekeepers as positive guys.

Response of host-bacterial colonization in shrimp to developmental stage, environment and disease.

The host-associated microbiota is increasingly recognized to facilitate host fitness, but the understanding of the underlying ecological processes that govern the host-bacterial colonization over development and, particularly, under disease remains scarce. Here, we tracked the gut microbiota of shrimp over developmental stages and in response to disease. The stage-specific gut microbiotas contributed parallel changes to the predicted functions, while shrimp disease decoupled this intimate association. After ruling out the age-discriminatory taxa, we identified key features indicative of shrimp health status. Structural equation modelling revealed that variations in rearing water led to significant changes in bacterioplankton communities, which subsequently affected the shrimp gut microbiota. However, shrimp gut microbiotas are not directly mirrored by the changes in rearing bacterioplankton communities. A neutral model analysis showed that the stochastic processes that govern gut microbiota tended to become more important as healthy shrimp aged, with 37.5% stochasticity in larvae linearly increasing to 60.4% in adults. However, this defined trend was skewed when disease occurred. This departure was attributed to the uncontrolled growth of two candidate pathogens (over-represented taxa). The co-occurrence patterns provided novel clues on how the gut commensals interact with candidate pathogens in sustaining shrimp health. Collectively, these findings offer updated insight into the ecological processes that govern the host-bacterial colonization in shrimp and provide a pathological understanding of polymicrobial infections.

Quantitative PCR Analysis of Gut Disease-Discriminatory Phyla for Determining Shrimp Disease Incidence.

There is evidence that gut microbial signatures are indicative of host health status. However, few efforts have been devoted to establishing an applicable technique for determining disease incidence by using gut microbial signatures. Herein, we established a quantitative PCR (qPCR)-based approach to detect the relative abundances of gut disease-discriminatory phyla, which in turn afforded independent variables for quantitatively determining the incidence of shrimp disease. Given the temporal dynamics of gut bacterial communities as healthy shrimp aged, we identified disease-discriminatory phyla after ruling out age-discriminatory phyla. The top 10 disease-discriminatory phyla contributed to an overall 93.2% accuracy in diagnosis (n = 103 samples from shrimp that were determined with high confidence to be healthy or that exhibited apparent disease symptoms and subsequent death), with 70% diagnosis accuracy at the disease onset stage, when symptoms or signs of disease were not apparent. 16S rRNA gene-targeted group-specific primers of five disease-discriminatory phyla were then designed according to their compositions within shrimp gut microbiota, and other primers were borrowed from previous studies. The relative abundances of the 10 disease-discriminatory phyla assayed by qPCR exhibited a high consistency (r = 0.946, P < 0.001) with those detected by Illumina sequencing. Notably, using the profiles of disease-discriminatory phyla assayed by qPCR and the corresponding weight coefficients as independent variables, we were able to accurately estimate the incidences of future disease outcome. This work establishes an applicable technique to quantitatively determine the incidence and onset of shrimp disease, which is a valuable attempt to translate scientific research into a practical application.IMPORTANCE Current studies have identified gut microbial signatures of host health using high-throughput sequencing (HTS) techniques. However, HTS is still expensive and time-consuming and requires a high technical ability, thereby impeding its application in routine monitoring in aquaculture. Hence, it is necessary to seek an alternative strategy to overcome these shortcomings. Herein, we establish a qPCR-based approach to detect the relative abundances of gut disease-discriminatory phyla, which in turn afford independent variables to quantitatively determine the incidence and onset of shrimp disease. Notably, there is a high consistency between the accuracies of disease diagnosis achieved by qPCR and HTS. This applicable technique makes important progress toward defining a diseased state in shrimp and toward solving an important animal health management-driven economic problem.

[Effects of Organic Pollutants on the Bacterioplankton Community in Hangzhou Bay].

Hangzhou Bay suffers from intensive anthropogenic disturbances and a huge amount of terrestrial inputs, and thus has become one of the most seriously contaminated coastal zones in China. There is evidence that microbes play a dominant role in pollutant biodegradation and serve as biomarkers for pollution levels. However, it remains unclear how the bacterioplankton communities respond to organic contaminants. To fill this knowledge gap, we collected surface water samples (0.5 m below the surface layer) from 13 sites across Hangzhou Bay and 8 control sites across its adjacent offshore areas. Using Illumina sequencing based on analysis of the bacterial 16S rRNA gene, we explored the effects of increasing organic pollution levels on the bacterioplankton community compositions (BCCs). The results revealed that the organic pollution level (A) in Hangzhou Bay (13.2±1.6) was significantly (P<0.001) higher than in the control zone (5.4±3.0). The distribution and diversity of bacterioplankton communities were significantly distinct between the two zones. The dominant bacterioplankton lineages in Hangzhou Bay were γ-Proteobacteria (24.4%±5.5%), α-Proteobacteria (16.5%±7.7%), and Planctomycetes (13.9%±8.6%), whereas those in the adjacent zones were Cyanobacteria (20.1%±7.5%), Bacteroidetes (18.4%±1.5%), Actinobacteria (17.5%±4.2%), γ-Proteobacteria (16.6%±1.2%), and α-Proteobacteria (14.3%±1.7%). Multivariate regression tree (MRT) analysis showed that the bacterioplankton community diversity was primarily affected by suspended particulates (SP), nitrite, oil, and organic pollutants, which respectively explained 22.0%, 6.5%, 6.0%, and 5.5% of the variance in diversity. Redundancy analysis (RDA) illustrated that the bacterioplankton community distribution was controlled by organic pollutants, COD, Chla, TN, nitrate, and salinity, which cumulatively governed 71.0% of the variation in BCCs. Organic pollutants alone controlled 6.5% variance, which was higher than any other single factor. Additionally, 35 sensitive species were identified via the indicator value method and their relative abundances were significantly associated (P<0.05 in each case) with the organic pollution level, thereby indicating their potential for evaluating coastal pollution. Collectively, our work demonstrates that BCCs are sensitive to coastal pollution and provides biomarkers for elevated pollution levels.

Starvation stress affects the interplay among shrimp gut microbiota, digestion and immune activities.

Aquatic animals are frequently suffered from starvation due to restricted food availability or deprivation. It is currently known that gut microbiota assists host in nutrient acquisition. Thus, exploring the gut microbiota responses would improve our understanding on physiological adaptation to starvation. To achieve this, we investigated how the gut microbiota and shrimp digestion and immune activities were affected under starvation stress. The results showed that the measured digestion activities in starved shrimp were significantly lower than in normal cohorts; while the measured immune activities exhibited an opposite trend. A structural equation modeling (SEM) revealed that changes in the gut bacterial community were directly related to digestive and immune enzyme activities, which in turn markedly affected shrimp growth traits. Notably, several gut bacterial indicators that characterized the shrimp nutrient status were identified, with more abundant opportunistic pathogens in starved shrimp, although there were no statistical differences in the overall diversity and the structures of gut bacterial communities between starved and normal shrimp. Starved shrimp exhibited less connected and cooperative interspecies interaction as compared with normal cohorts. Additionally, the functional pathways involved in carbohydrate and protein digestion, glycan biosynthesis, lipid and enzyme metabolism remarkably decreased in starved shrimp. These attenuations could increase the susceptibility of starved shrimp to pathogens infection. In summary, this study provides novel insights into the interplay among shrimp digestion, immune activities and gut microbiota in response to starvation stress.

Quantitative prediction of shrimp disease incidence via the profiles of gut eukaryotic microbiota.

One common notion is emerging that gut eukaryotes are commensal or beneficial, rather than detrimental. To date, however, surprisingly few studies have been taken to discern the factors that govern the assembly of gut eukaryotes, despite growing interest in the dysbiosis of gut microbiota-disease relationship. Herein, we firstly explored how the gut eukaryotic microbiotas were assembled over shrimp postlarval to adult stages and a disease progression. The gut eukaryotic communities changed markedly as healthy shrimp aged, and converged toward an adult-microbiota configuration. However, the adult-like stability was distorted by disease exacerbation. A null model untangled that the deterministic processes that governed the gut eukaryotic assembly tended to be more important over healthy shrimp development, whereas this trend was inverted as the disease progressed. After ruling out the baseline of gut eukaryotes over shrimp ages, we identified disease-discriminatory taxa (species level afforded the highest accuracy of prediction) that characteristic of shrimp health status. The profiles of these taxa contributed an overall 92.4% accuracy in predicting shrimp health status. Notably, this model can accurately diagnose the onset of shrimp disease. Interspecies interaction analysis depicted how the disease-discriminatory taxa interacted with one another in sustaining shrimp health. Taken together, our findings offer novel insights into the underlying ecological processes that govern the assembly of gut eukaryotes over shrimp postlarval to adult stages and a disease progression. Intriguingly, the established model can quantitatively and accurately predict the incidences of shrimp disease.

Integrating molecular and ecological approaches to identify potential polymicrobial pathogens over a shrimp disease progression.

It is now recognized that some gut diseases attribute to polymicrobial pathogens infections. Thus, traditional isolation of single pathogen from disease subjects could bias the identification of causal agents. To fill this gap, using Illumina sequencing of the bacterial 16S rRNA gene, we explored the dynamics of gut bacterial communities over a shrimp disease progression. The results showed significant differences in the gut bacterial communities between healthy and diseased shrimp. Potential pathogens were inferred by a local pathogens database, of which two OTUs (affiliated with Vibrio tubiashii and Vibrio harveyi) exhibited significantly higher abundances in diseased shrimp as compared to healthy subjects. The two OTUs cumulatively contributed 64.5% dissimilarity in the gut microbiotas between shrimp health status. Notably, the random Forest model depicted that profiles of the two OTUs contributed 78.5% predicted accuracy of shrimp health status. Removal of the two OTUs from co-occurrence networks led to network fragmentation, suggesting their gatekeeper features. For these evidences, the two OTUs were inferred as candidate pathogens. Three virulence genes (bca, tlpA, and fdeC) that were coded by the two candidate pathogens were inferred by a virulence factor database, which were enriched significantly (P < 0.05 in the three cases, as validated by qPCR) in diseased shrimp as compared to healthy ones. The two candidate pathogens were repressed by Flavobacteriaceae, Garvieae, and Photobacrerium species in healthy shrimp, while these interactions shifted into synergy in disease cohorts. Collectively, our findings offer a frame to identify potential polymicrobial pathogen infections from an ecological perspective.

The gut eukaryotic microbiota influences the growth performance among cohabitating shrimp.

Increasing evidence has revealed a close interplay between the gut bacterial communities and host growth performance. However, until recently, studies generally ignored the contribution of eukaryotes, endobiotic organisms. To fill this gap, we used Illumina sequencing technology on eukaryotic 18S rRNA gene to compare the structures of gut eukaryotic communities among cohabitating retarded, overgrown, and normal shrimp obtained from identically managed ponds. Results showed that a significant difference between gut eukaryotic communities differed significantly between water and intestine and among three shrimp categories. Structural equation modeling revealed that changes in the gut eukaryotic community were positively related to digestive enzyme activities, which in turn influenced shrimp growth performance (λ = 0.97, P < 0.001). Overgrown shrimp exhibited a more complex and cooperative gut eukaryotic interspecies interaction than retarded and normal shrimp, which may facilitate their nutrient acquisition efficiency. Notably, the distribution of dominant eukaryotic genera and shifts in keystone species were closely concordant with shrimp growth performance. In summary, this study provides an integrated overview on direct roles of gut eukaryotic communities in shrimp growth performance instead of well-studied bacterial assembly.

Bacterioplankton assembly and interspecies interaction indicating increasing coastal eutrophication.

Anthropogenic perturbations impose negative effects on coastal ecosystems, such as increasing levels of eutrophication. Given the biogeochemical significance of microorganisms, understanding the processes and mechanisms underlying their spatial distribution under changing environmental conditions is critical. To address this question, we examined how coastal bacterioplankton communities respond to increasing eutrophication levels created by anthropogenic perturbations. The results showed that the magnitude of changes in the bacterioplankton community compositions (BCCs) and the importance of deterministic processes that constrained bacterial assembly were closely associated with eutrophication levels. Moreover, increasing eutrophication significantly (P < 0.001) attenuated the distance decay rate, with a random spatial distribution of BCCs in the undisturbed location. In contrast, the complexity of interspecies interaction was enhanced under moderate eutrophication levels but declined under heavy eutrophication. Changes in the relative abundances of 27 bacterial families were significantly correlated with eutrophication levels. Notably, the pattern of enrichment or decrease for a given bacterial family was consistent with its known ecological functions. Our findings demonstrate that the magnitude of changes in BCCs and underlying determinism are dependent on eutrophication levels. However, the buffer capacity of bacterioplankton community is limited, with disrupted interspecies interaction occurring under heavy eutrophication. As such, bacterial assemblages are sensitive to changes in environmental conditions and could thus potentially serve as bio-indicators for increasing eutrophication.