Abstract
Microbiomes can modulate plant and animal responses to environmental fluctuations, revealing their central role in global biogeochemical cycles and the physiology of individual organisms, and in how these fluctuations are modulated. In marine ecosystems, Eastern Boundary Upwelling Systems (EBUS) allow us to examine host-microbiome dynamic interactions that respond to environmental variation, offering critical insights into the potential impacts of global environmental change. This study investigates the microbiota response of two marine invertebrates that co-occur at sites with strong synoptic and seasonal fluctuations in temperature and dissolved oxygen concentration between winter and summer. Specimens of the sessile marine tunicate Pyura chilensis (“Piure”) and the porcelain crab Allopetrolisthes punctatus (“Changai”) were collected at about 10 m depth under contrasting oxygen regimes in summer and winter at Chome, central Chile. Microbiota from gills and digestive glands were assessed via 16S rRNA Nanopore sequencing, which showed clear differences in bacterial community structure and composition between species and across summer and winter surveys. Functional predictions suggest that shifts in the microbiota are primarily associated with variations in exposure to hypoxic waters induced by intensified upwelling. This suggests that the microbiomes may positively influence host performance under environmental stress, enhancing resistance to low dissolved oxygen conditions.
Introduction
Thanks to advances in molecular and data analysis techniques, the last decade has witnessed a revolution in our understanding of the critical role microbial communities play in the biogeochemical cycles, on which all life depends (; ), and in the physiology and homeostasis of individual organisms. The discovery that microbial communities inhabiting other organisms (“internal microbiota”) can modulate plant and animal responses to environmental fluctuations, once thought to be largely independent of microbes, is reshaping our conception of what constitutes an individual, how species’ niches should be defined, and the units upon which natural selection operates (; ). We are only beginning to understand how natural environmental variability influences interactions among microbial communities and their hosts. Understanding how microbiotas respond to the environment is essential for anticipating the potential effects of anthropogenic and climatic disturbances. In marine ecosystems, key environmental variables such as seawater temperature, pH, and dissolved oxygen are being altered by climate change and other global drivers. In this context, Eastern Boundary Upwelling Systems (EBUS), where the wind-driven onshore transport of subsurface waters produces large natural fluctuations in these biologically critical variables, offer an exceptional opportunity to assess the impacts of environmental variability on microbiota communities. Indeed, besides temperature and nutrient levels, large fluctuations in pH and dissolved oxygen concentration characterize upwelling shores (; ; ) and provide clues to global phenomena such as ocean acidification and deoxygenation (; ; ).
The coastal ocean is a heterogeneous and temporally variable environment, with complex circulation patterns and oceanographic conditions, strongly influenced by coastal bathymetry and the interaction with wind forcing and riverine inputs (; ). Physical variables such as temperature, salinity, dissolved oxygen (DO), and pH fluctuate widely over temporal scales that ranging from seasons to days, hours, and even minutes (; ; ; ; ). Over the past decade, extreme environmental variability in the coastal ocean has been increasingly reported, including in upwelling nearshore ecosystems (; ). While temperature has remained relatively stable along these shores (; ; ), several environmental changes associated with climate change have been reported, including shifts in primary productivity and in the intensity and frequency of hypoxic events (; ; ). The ecosystem response to hypoxia depends on the timing, predictability, and intensity of oxygen depletion. Sustained hypoxia in nearshore areas is generally accompanied by a loss of biodiversity and ecosystem functions, and by increased mortality of benthic organisms, followed by some degree of renewed colonization upon return to normoxia (; ).
The response of marine metazoans to oxygen depletion has been studied in terms of behavioral, physiological, and metabolic changes (). Various taxa exhibit different levels of tolerance to low-oxygen conditions (; ). Among large metazoans, cnidarians and annelids are generally considered most tolerant to hypoxia (; ), whereas crustaceans and fish are known to be more sensitive to low-oxygen conditions (; ). Although an organism’s survival under environmental stress, such as a decline in oxygen availability, has been typically viewed as the result of its own capacity to respond or acclimate, it is now known that its microbiota plays a role in responding to environmental changes (; ). Changes in water quality can reshape microbial communities and may contribute to organism health if there is a significant alteration in their microbiota structure and abundance (; ). Given the influence of the microbiota on the physiology and metabolism of marine organisms, microbes are expected to contribute to their adaptation to a changing ocean environment (; ). For instance, exposure to hypoxia has been shown to alter the microbiota of mussels (; ) and fish (; ; ). Reporting changes in microbiota richness, and an enrichment of microbiota with functions related to methane metabolism, carbon fixation, and key metabolic pathways has been observed in fish exposed to hypoxia (; ; ).
In the mussel, Mytilus coruscus, changes in their gut microbiota have been triggered by changes in DO levels (). A recent study of the blue mussel Mytilus chilensis gut exposed to hypoxic conditions revealed significant alterations in microbial composition, with a shift toward facultative anaerobes (). Functional analysis indicated a decline in critical microbial functions associated with nutrient metabolism and immune support, potentially jeopardizing the health and survival of the host (). Highlighting the importance of understanding the role of host-microbiota interactions in increasingly altered marine environments, especially in regions impacted by coastal hypoxia. This study focuses on the microbiota of benthic invertebrates from the upwelling coast of central Chile, which experiences large seasonal and shorter-term changes in the oxygen conditions of near-bottom waters (). We selected two species of filter feeders with contrasting behavior: the sessile tunicate Pyura chilensis, widely distributed along the Chilean coast (), a conspicuous species considered a bioengineer, a dominant competitor of hard substrates, and a species with economic relevance for small-scale aquaculture activity in Chile (; ). The other specie evaluated in this study is the porcelanid crab Allopetrolisthes punctatus, which is highly mobile, locally abundant, but patchily distributed between northern and south-central Chile (18–40°S) (). This crab shows gregarious behavior, forming dense aggregations that form multi-layered clusters (). We characterized the gill and gut microbiotas of both species and compared them with the ambient microbiota’s composition and changes between the upwelling season in summer–characterized by low DO levels–and the non-upwelling period in winter, which has higher DO levels. We hypothesized that, like their hosts, the microbiota of mobile and sessile marine invertebrates varied in response to environmental stress (i.e., DO), thereby influencing host adaptation capacity.
Materials and methods
Study area
The study was conducted at a nearshore site on the central Chile coast - Chome, 36° 46′ 24′′ S, 73° 12′ 38′′ W (Figure 1A). This small fishing cove is located just north of the Biobío River mouth; its near-bottom conditions are strongly influenced by the seasonality of coastal upwelling and by the interaction between tidal currents and the Biobío submarine canyon (; ; ). Temporal variability in near-bottom hydrographic conditions, temperature and dissolved oxygen, at Chome was characterized from continuous records gathered from July 2024 to May 2025 using a suite of environmental sensors moored at ca. 2 m above the bottom at a nominal depth of 30 m. The sensor array included loggers for dissolved oxygen and temperature (MiniDOT, PME, USA), pressure (HOBO U20, Onset Computer, USA), and conductivity (HOBO U24, Onset Computer, USA). All loggers were programmed to record data at 10-min intervals and were serviced every 3–4 months. The data gathered from the loggers were quality controlled, curated, and converted into hourly time series prior to the analysis. To visualize the seasonal variability in exposure to hypoxic conditions, the fraction of recorded values below 2 mg/L was computed for each month. The seasonal variability in water temperature, and its co-variation with DO variability, were assessed by computing monthly distributions from all recorded values. Additionally, the RMS (root mean squared) of first-differenced hourly time series was computed for each month as an index of shorter-term variability in near-bottom temperatures. Higher short-term variability in temperature should indicate a stronger influence of advective processes in changing both temperature and dissolved oxygen conditions in bottom waters.
FIGURE 1
Biological sampling
Thirty adults - of combined male and female specimens - of both the tunicate Pyura chilensis (Piure) and the porcelain crab Allopetrolisthes punctatus (Changai) were collected from two different sampling locations through scuba diving at about 10 m depth in July (winter) of 2024 February 2025 (summer). Tissue samples were extracted from the digestive gland and gills of each specimen, then stored at −80 °C until DNA extraction. Additionally, to evaluate the ambient microbiota, two samples of water from the same location as the organisms was collected and filtered through a 0.22 μm membrane filter (Isopore™). DNA extraction was performed using the PureLink™ Microbiome DNA Purification Kit (Invitrogen), following the manufacturer’s instructions. DNA quality and concentrations were assessed using a Qubit 4 fluorometer (ThermoScientific, USA). DNA of all individuals was diluted to the same concentration, 100 ng/ul, to perform a two-pool for tissue and organisms of 10 individuals to continue with the 16S rRNA amplification. The sampling design for this study adhered to the three Rs (3Rs) guidelines for animal testing.
Nanopore library synthesis and sequencing
The complete 16S rRNA gene was amplified from two replicates per sample using the universal primers 27F and 1492R. LongAmp Taq DNA polymerase (New England Biolabs, USA) was utilized for amplification at a reaction volume of 15 μL under the following conditions: 95 °C for 1 min, followed by 30 cycles at 95 °C for 20 s, 56 °C for 30 s, and 65 °C for 1 min, concluding with a final extension at 65 °C for 5 min. The efficacy of the PCR was assessed through 1% agarose gel electrophoresis. PCR products were purified with a 1:2 sample-to-bead ratio using Agencourt AMPure XP beads (Beckman Coulter, USA) for Nanopore library synthesis. The samples were placed in a magnetic stand for washing with freshly prepared 70% ethanol. After rinsing in ultrapure sterile water, the product was separated from the beads by placing it in a magnetic stand. The purified amplicon was quantified again using the Qubit 4 fluorometer (ThermoScientific, USA) and utilized as a template for library synthesis with the 16S rDNA Barcoding Kit (SQK-RAB204, Oxford Nanopore Technologies) following the manufacturer’s instructions. Briefly, the purified 16S rDNA amplicon from a pool of ten individuals was mixed with barcodes. Then, a PCR was performed as stated in the previous step, with a final reaction volume of 50 μL, using the LongAmp Taq Polymerase (New England Biolabs). The PCR product was incubated for 5 min at room temperature using a HulaMixer (ThermoScientific, USA) and subsequently cleaned with Agencourt AMPure XP beads. The purified product was eluted in 10 μL of the elution buffer (10 mM Tris-HCl, pH 8.0, with 50 mM NaCl) following the procedure for washing the magnetic beads. The final concentration of the library was measured using a Qubit 4 fluorometer (Thermo Fisher Scientific, USA), and libraries were tested using High Sensitivity ScreenTape D5000 (Agilent, USA) according to the manufacturer’s instructions on a TapeStation Bioanalyzer 2100 (Agilent, USA). Following the methodology of Oxford Nanopore Technologies, the libraries were pooled in multiplex mode and placed into the flowcell MK1 Spot-ON FLO-MIN107-R9. The ZymoBIOMICS Microbial Community DNA Standard (Cat. No. D6305) was included in the runs. Two biological samples of each tissue, organism, and sample point were sequenced. Sequencing efficiency was monitored using the software MinKNOW 2.0 from Oxford Nanopore Technologies.
Bioinformatic analysis
Fast5 files were base-called using Guppy (version 3.2.2, Oxford Nanopore Technologies, UK), and a filtering step was applied to retain only sequences with a Q-score of ≥7 (quality filter). The sequencing performance of summer and winter samples was analyzed separately using pycoQC (). The reads demultiplex, primers, and adapter trimming were performed with Porechop (). Taxonomic classification of reads was performed with EMU v.3.4.5 (). The minimum abundance per sample was filtered at 0.1% using the R package dplyr (). Supplementary Table 1 reports all data used for the analysis (count, taxonomy, abundance, and metadata tables). Diversity patterns and ordination analyses were determined using the R packages phyloseq, vegan, and ape (). Taxonomic distribution and distance were analyzed with phyloseq and fantaxtic, while the statistical comparison was performed with the R package ggpubr. The taxonomic abundance comparison among bacterial communities present in each evaluated tissue and organism was performed by heat-trees using metacoder R package (). The bacterial community composition was displayed by the Krona visualization tool running in Galaxy ().
Functional MetaCyc gene pathways prediction (, ) of bacterial communities among organisms, between tissues, and between seasons was performed using PICRUSt2 (). Read abundances were normalized by 16S rRNA gene copy number, and gene abundances were estimated by multiplying the normalized read counts by the predicted gene copy numbers. Differentially abundant metabolic pathways were analyzed with the STAMP software, using Welch’s t-test with a p-value cutoff of 0.05 ().
Community interactions among bacterial species were inferred using the network approach (R package SpiecEasi; Sparse Inverse Covariance Estimation for Ecological Association Inference) (), and the igraph software (). The correlation matrix obtained was filtered using an absolute correlation score greater than or equal to 0.5, where each node represents one species and each edge represents the correlations between the species abundances. Finally, relationships among bacterial communities were identified in gill samples, and their putative functional pathways were determined using the R package mixOmics ().
Results
Oxygen and temperature variability in bottom waters
The time series of temperature and dissolved oxygen concentrations revealed a clear and strong seasonal signal in both variables, with substantial differences in the oxygenation regime of bottom waters between winter and summer months, accompanied by a seasonal change in not only the mean values of temperature but also in its short-term (hourly) variability (Figures 1B,C). From December to March (austral summer), DO levels remained below 2 mg/L for over 60% of the time (Figure 1B). Throughout this period, the frequency distributions of DO were skewed toward low and very low values. Contrasting distributions of DO were observed for the May-August period (autumn-winter), and especially in the June-July window (austral winter), with less than 4% of the recorded DO values below 2 mg/L. Wide and symmetrical distributions of DO values were observed in the transition months of September and April-May (Figure 1B). Exposure to severe hypoxia was highest in December-February, when >75% of the DO values recorded was below 2 mg/L (Figure 1B, red line and symbols). Strong seasonal pattern was also observed for temperature (Figure 1C). However, the minimum and maximum values were not observed in the same months as dissolved oxygen (Figure 1C). The lowest temperatures were recorded in November (spring) with 10.8 °C, and the warmest month was May with 12.5 °C (late autumn). The narrowest range of temperature was observed in June-July (2 °C), when the short-term (hourly) variability of temperature was also minimal (Figure 1C, blue line and symbols). The largest amplitudes in recorded temperature values were observed in December and March (>5 °C), whereas the short-term variability of temperature was highest in February-March (ca. 4 times as in July-August).
Invertebrates’ microbiota community structure under oxic and hypoxic conditions
A total of 2,659,031 reads were obtained after trimming and filtering, comprising 1,254,397 for the summertime, hypoxia-dominated conditions and 1,404,634 reads for the wintertime oxygenated conditions (Supplementary Figure 1 and Supplementary Tables 2, 3). Using the data of both tissues per species alpha diversity analyses revealed significant differences among the animals’ microbiota and the environment, where is observed more diversity with significant differences according to the Shannon and Fisher index (Figure 2A). Relate with the oxygenation condition, significant differences in Shannon index were observed between the crab and the tunicate under hypoxic conditions (Figure 2A). Simpson index also demonstrates significant differences between the crab and the tunicate under oxic and hypoxic condition. Also, the microbiota diversity in the tunicate is highly differences when is compare oxic and hypoxic condition according Simpson index (Figure 2A). A Nonmetric Multidimensional Scaling (NMSD) analysis indicated differences in the community structure of bacteria in the microbiota and those in the ambient microbial communities (ANOSIM R = 0.9, p < 0.002, Figure 2B), consistent with differences between species (Figure 2B). Similarly, there were differences in the bacterial community structure across seasons (Principal Coordinates Analysis; PCoA, Figure 2C) and organism (Figure 2D). However, there were no apparent shifts in the microbiota of tunicates and crabs when comparing between oxygenation regimes (Supplementary Figure 2).
FIGURE 2
Compositional changes in the bacterial community during oxic and hypoxic conditions
Differences in microbiota communities’ abundance were observed in the digestive gland and gill tissues between the tunicate and the crab (Figure 3). In the tunicate P. chilensis, the bacterial community at family level exhibited three clades: two for hypoxia, one each for gills and digestive gland bacteria, and another clade associated with oxic conditions that also separates both tissues (Figure 3A). Under both oxygenation regimes, the most abundant bacterial family in the tunicate’s microbiota is Nostocaceae from the Cyanobacteria order (Figure 3A). As for the crab A. punctatus, its bacterial community at family level separated into two main clades: hypoxic and oxic, also exhibiting differentiation between digestive gland and gill tissues (Figure 3A). The family Ahrensiaceae becomes prominent under hypoxic conditions while the family Gemmateceae dominates under oxic regime (Figure 3A). Notably, bacteria from the families Methylothermaceae and Thiotrichaceae were identified in both species under both oxygenation conditions.
FIGURE 3
A total of 33 bacterial species constitute the digestive gland core, while 10 bacterial species were identified for the gill core (Figure 3B). The tunicate’s digestive gland revealed 90 and 54 exclusive species during the hypoxic and oxic periods, respectively. Samples from the crab’s digestive gland 37 and 64 exclusive species in hypoxic and oxic conditions, respectively (Figure 3B). In the gills, the tunicate exhibited the highest number of exclusive bacterial species during the hypoxia-dominated period (264 species, Supplementary Table 4). In contrast, during the oxygenated season, there were only 79 exclusive bacterial species in the tunicate’s microbiota (Figure 3B). In the crab’s gills, 87 and 43 exclusive bacteria were identified under oxic and hypoxic conditions, respectively (Figure 3B). A core of bacteria was identified among species and oxygen regimes for each tissue with 33 bacterial species in the digestive gland core, and 10 in the gills core, with abundance differences among the oxygen condition and organism (Figure 3C).
Shifts in bacterial community interactions during oxic and hypoxic conditions
Network analyses were conducted to unveil bacterial co-occurrence. The network displayed greater complexity when animals were exposed to oxic compared to hypoxic conditions (Figure 4). Under hypoxia-dominated conditions, there were five hubs (i.e., key taxa by phylum) involving interaction among Actinomycetota, Bacteroidota, Cyanobacteriota, Bacillota, Pseudomonadota (Figure 4A). The taxa with the highest interaction (i.e., closeness) were Bacteroidota and Cyanobacteriota (Figure 4B). Under oxygenated conditions, six hubs were observed: the phyla Actinomycetota, Cyanobacteriota, Bacillota, Planctomycetota, Pseudomonadota, and Verrucomicrobiota (Figure 4C). Interestingly, Cyanobacteriota exhibited the highest node influence in the connected network (i.e., eigenvector) for both oxygenation regimes (Figures 4B,D). At the species level, Salmonella enterica, Escherichia coli, and Pseudahrensia todarodis exhibited the highest closeness levels under hypoxia. Meanwhile, Pseudomonata species such as Pluralibacter pyrinus and Helicobacter pylori displayed high values in oxic conditions.
FIGURE 4
Using the date of both organisms, oxygen condition-specific microbiota was identified in gills and digestive glands (Supplementary Figure 3 and Supplementary Table 5). Under hypoxic conditions, 56 species were annotated, representing bacterial classes such as Alphaproteobacteria, Bacteroidia, and Flavobacteriia (Supplementary Figure 3 and Supplementary Table 6). Under oxygenated conditions, 58 bacterial species were present in both tissues; the most abundant bacterial class was Cyanobacteria, primarily represented by the genus Synechococcus, which dominated under this regime (Supplementary Figure 3 and Supplementary Table 5).
Functional prediction of the host microbiota
Using the data of both tissues per species, a functional prediction of MetaCyc gene pathways was conducted from the identified bacterial species (Nanopore data using Emu; Supplementary Table 6). In terms of putative microbiota function in hypoxic conditions, ten gene pathways were significantly different between the tunicate and the crab under hypoxic conditions, including amino acid degradation, cell structure biosynthesis, and carbohydrate biosynthesis/degradation (Figure 5A). Under oxic conditions 20 pathways were identified with significant differences between the tunicate and the crab (Figure 5B). Among the highly represented gene pathways are Nucleoside and Nucleotide Biosynthesis, Amino Acid Biosynthesis, and Cofactor-Prosthetics Group-electron Carrier Cofactor-and Vitamin Biosynthesis (Figure 5B).
FIGURE 5
Notably, the fermentation of short-chain fatty acids appeared in a higher proportion in the tunicate’s tissues during hypoxia, compared to the crab under the same conditions (Figure 6A). Furthermore, the tunicate exhibited a greater potential for the fermentation of acetate, butanoate, propanoate, and pyruvate to ethanol during hypoxia, whereas the crab displayed a higher potential for these metabolic processes under oxic conditions (Figure 6A). A high proportion of bacteria associated fermentation of short-chain fatty acids were observed in both organism under hypoxia and oxic conditions, for instance, Planctomycetota, Verrucomicrobiota, Actinomycetota, and Pseudomonadota which are involved in butanoate fermentation in the gill tissues of both species (Figure 6B).
FIGURE 6
Finally, a network analysis was conducted to explore the relationships between the most abundant gill microbiota families and the MetaCyc gene pathways. For the crab’s microbiota, Vibriocaceae found under hypoxic conditions were positively correlated only with Glycolysis and negatively correlated with most of the MetaCyc gene pathways, including the TCA cycle, fermentation, carbohydrate biosynthesis, and inorganic nutrient metabolism (Supplementary Figure 4). In contrast, the Enterobacteriaceae, Roseobacteraceae, Helicobacteraceae, and Bartonellaceae families were positively correlated with most of the MetaCyc gene pathways (Supplementary Figure 4). For the microbiota of tunicate gills, the bacterial families observed under hypoxic and oxic conditions were positively correlated with MetaCyc gene pathways (Supplementary Figure 4). For instance, the bacteria of the Roseobacteraceae family under hypoxic conditions are positively correlated with amino acid degradation. In the case of microbiota under oxic conditions, the Clostridiaceae, Carnobacteriaceae, Corynebacteriaceae, and Xanthomonadaceae families exhibited a positive correlation with Carbohydrate Degradation, Secondary Metabolite Biosynthesis, Cofactor Prosthetic Groups, Electron Carriers, and Vitamin Biosynthesis.
Discussion
Microbial communities play a central role in marine biogeochemical cycles. Their composition responds to environmental variability over various spatial and temporal scales. Although efforts have been made to understand spatial/biogeographical patterns (), successional processes (), and responses to hypoxia in environmental microbial communities (), little is known about the responses of microbiota in marine organisms to changes in critical environmental variables. The ubiquity of microorganisms in the environments in which animals are immersed suggests a functional connection between the ambient microbial community and the well-being of various hosts. For instance, juvenile invertebrates or fish larvae deprived of their gut microbiota experience adverse effects on their growth and survival (; ). The role of microbial communities in marine ecosystem functioning is unequivocal, with microorganisms being key drivers of the ocean’s biogeochemistry and having a critical role in maintaining animal health ().
Dissolved oxygen is a primary driver of coastal animal benthic community structure, behavior, physiology, and ultimately species survival (). In recent years, ocean deoxygenation has increased, undermining the adaptation and survival of marine organisms worldwide (; ; ). The Chilean coast experiences seasonal upwelling, leading to recurring environmental hypoxia, especially during the summer (; ). Our study area was chosen as a representative site for the seasonal and shorter-term variation in dissolved oxygen levels that are forced mainly by wind-driven upwelling, but also by the interaction of coastal flows - including tidal currents - and the local topography (). Thus, conditions at our study site ranged between moderate to severe hypoxia during the austral spring-summer season and well-oxygenated conditions in the austral winter. Here, variations in the microbiota of model benthic organisms, such as tunicates and decapod crustaceans, were evaluated to elucidate their role in marine organism adaptation to environmental hypoxia. Upwelling-driven hypoxia is coupled with shifts in primary productivity and turbidity. For filter-feeding tunicates like P. chilensis, the distribution of primary producers likely drives Cyanobacteria dominance during high light transmittance. Conversely, the influx of Ahrensiaceae in crabs during hypoxia could be explained to sediment resuspension. Thus, while feeding strategies and particle distribution play a role in the microbiota, dissolved oxygen remains the primary physiological driver.
Microbiota studies in ascidians () and crabs () have reported species-specific microbiota composition. Also, microbiota variation has been reported in different marine organisms in response to low-oxygen levels (; ; ). Here, the microbiota diversity in the Tunicate is more affected by oxygen conditions than in the Crab. Interestingly, in the Tunicate, microbiota diversity increases under hypoxic conditions, which suggests a role for this microbiota in the Tunicate’s adaptation to low oxygen availability. In the Chinese mitten crab (Eriocheir sinensis), changes in intestinal microbiota have been reported in response to hypoxia (), with a high dominance of the Bacillota and Bacteroidota phyla. At the class level, the authors reported a decrease in Gammaproteobacteria in the hypoxia-exposed group compared with the control (). In our study, the class Gammaproteobacteria was more abundant in samples under hypoxic conditions than under oxic conditions. In another study, performed in the oyster Crassostrea hongkongensis, a significant increase in Pseudomonadota abundance was reported under low-oxygen conditions (). Interestingly, the microbiota in the tunicate’s gill displayed a high abundance of Bacteroidota bacteria under hypoxic conditions, compared to what has been reported in the Chinese mitten crab. Moreover, Bacilli bacteria have been associated with hypoxia in Chinese mitten crab (), and oysters (), suggesting that those bacterial groups are common among marine invertebrates exposed to hypoxia. At the tissue level, it is important to note that, while the gill tissue is permanently exposed to ambient waters and, therefore, the associated changes in oxygenation levels, the digestive gland is contained and, therefore, under permanent hypoxic to anoxic conditions. Consistently, in our study, Bacillota were observed in the tunicate’s digestive gland in both seasons (hypoxia-dominated and oxic conditions), similar to what was reported in the mussels Mytilus couscous (M. couscous) () and M. chilensis ().
Microbiota interactions play a crucial role in maintaining ecological balance, and a robust microbial network (i.e., high connectivity and sub-nodes) is associated with a balanced microbiota with a positive impact on the host (). For example, a study on the fish Bostrichthys sinensis showed less robust microbiota networks under hypoxic conditions than those exposed to oxic conditions (). Overall, our model species exhibited more nodes and sub-nodes in oxic than in hypoxic conditions, suggesting that hypoxic stress affects both the host and the microbiota.
Hypoxic stress triggers a cascade of metabolic effects leading to a high energy demand, which can induce changes in the organism’s response to stressors, including pathogens or even affect digestion regulation (; ). Due to the microbiota’s control over several host processes at the metabolic level (; ; ), understanding its function is essential for grasping the regulation of host metabolism. Even though assessing gene expression is beyond the scope of this study, the functional prediction analysis of the microbiota revealed the potential for various metabolic functions. In this context, environmental hypoxia appeared to promote processes such as carbohydrate biosynthesis, degradation, and amino acid degradation. This suggests that reduced oxygen levels have an impact beyond the known physiological stress of marine animals subjected to hypoxia, but also at the microbiota level in terms of composition and functioning. This aligns with observations conducted on the benthic fish Bostrichthys sinensis, in which similar functional network differences were reported comparing control (oxic) to hypoxic conditions ().
In humans, the short-chain fatty acids (SCFAs) produced by gut microbiota play multiple roles (). Despite the taxonomic dominance of Cyanobacteria in the tunicate under oxic conditions, our functional predictions - normalized for gene copy number - distinctly identified an enrichment of SCFA fermentation pathways during hypoxia. This suggests that the predicted metabolic shift toward anaerobic functions (i.e., butanoate and acetate production), is a response to hypoxia independently of the dominant aerobic taxa. Moreover, SCFAs have been associated with stabilizing the hypoxia-inducible factor (HIF) gene in epithelia with healthy microbiota (). This occurs under limited oxygen concentrations, triggering processes such as energy metabolism, inflammation, and apoptosis (; ). Although this mechanism has not been reported for invertebrates, high abundances of bacteria associated with butyrate fermentation were found under hypoxia within the tunicate’s digestive gland, suggesting a similar regulatory role of HIF in invertebrates. This type of response would probably prevent cell apoptosis and aid the host organism’s adaptation to prolonged hypoxia. This aligns with the role of butyrate as a histone deacetylase inhibitor and regulator of intestinal genes (; ). Indeed, butyrate inhibits the PHDs (pyruvate dehydrogenase) gene in human microbiota and stabilizes HIF (). Altogether, this suggests a crucial symbiotic relationship where the microbiota’s metabolic byproducts, such as butyrate, are key to the host’s physiological adaptation to low-oxygen environments.
Although a stable core of species persists across seasons reaching 33 in the digestive gland, 10 in the gills, the microbiota is primarily characterized by large-scale shifts in overall community structure and the recruitment of exclusive species during hypoxia. Thus, the non-core community reacts to a changing environment. Upwelling-driven environmental variability in coastal waters, particularly the seasonal hypoxia that characterizes the shores of central Chile and other Eastern Boundary regions, profoundly impacts animal physiology at both the host and microbiota levels. These impacts lead to a drastic reduction in microbial diversity and significant shifts in community composition, favoring anaerobic and facultative anaerobic taxa in the exposed gills and, interestingly, also in the permanently hypoxic micro-environment of the digestive glands. Crucially, oxic conditions host a different microbial community and function, indicating persistent alterations rather than transient perturbations. These microbial shifts are likely to alter metabolic functions, especially within the host gut, with considerable implications for the host’s health. While microbiota compositional shifts are clearly identified, potential functions are predicted with their limitations (). Thus, future studies are required, for instance, to use a common-garden design to evaluate, under controlled conditions, changes in the microbiome, and to include metatranscriptome and metabolome analyses to determine whether the functional consequences inferred in this study are consistent. Finally, understanding the effects of cumulative stress, particularly the response to recurrent environmental factors such as hypoxia, is vital for accurately predicting animal adaptation and the long-term resilience and health of aquatic ecosystems. Here, model invertebrate species from subtidal waters along the Chilean coast, and their microbiotas, can shed light on how coastal organisms are coping with the local manifestations of climate change and anthropogenic pollution.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Ethics statement
The animal study was approved by Ethic Committee of the University of Concepcion, Chile. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
VV-M: Formal analysis, Methodology, Conceptualization, Writing – original draft, Funding acquisition, Writing – review & editing, Investigation. FT: Conceptualization, Funding acquisition, Writing – review & editing, Formal analysis, Writing – original draft. CG-E: Writing – original draft, Conceptualization, Writing – review & editing. MM: Formal analysis, Conceptualization, Funding acquisition, Visualization, Writing – review & editing. SN: Conceptualization, Writing – review & editing. MG: Conceptualization, Writing – review & editing. MM-R: Formal analysis, Writing – review & editing, Data curation. GJ: Conceptualization, Writing – original draft, Methodology, Writing – review & editing, Funding acquisition, Formal analysis.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This research work was supported by COPAS Coastal (ANID FB210021), at Universidad de Concepción.
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2026.1764313/full#supplementary-material
Supplementary Figure 1Rarefaction curve of raw data was obtained for each sample.
Supplementary Figure 2Microbiota diversity for the Crab and Tunicate gills and digestive gland was determined by Alpha diversity, which was assessed using Shannon, Simpson, and Fisher diversity index.
Supplementary Figure 3Microbiota presented in gills and digestive gland tissues during hypoxic and oxic conditions. (A) The upset plot showed the number of bacterial species present and shared in both tissues. (B) Bacteria were identified at the genus level in gills and the digestive gland under each oxygen condition.
Supplementary Figure 4Network analysis of bacterial families and the Top MetaCyc gene pathways designated for the gill microbiota of both analyzed organisms. Red: bacteria families in a hypoxic condition; Green: bacteria families in an oxic condition. The red lines indicate negative correlations, and the green lines indicate positive correlations.
References
1
AdamsK. A.BarthJ. A.ChanF. (2013). Temporal variability of near-bottom dissolved oxygen during upwelling off central Oregon.J. Geophys. Res. Oceans1184839–4854. 10.1002/jgrc.20361
2
AguirreC.García-LoyolaS.TestaG.SilvaD.FaríasL. (2018). Insight into anthropogenic forcing on coastal upwelling off south-central Chile.Elem. Sci. Anth.6:59. 10.1525/elementa.314
3
ApprillA. (2017). Marine animal microbiomes: Toward understanding host–microbiome interactions in a changing ocean.Front. Mar. Sci.4:222. 10.3389/fmars.2017.00222
4
BacciG.MeriggiN.ChengC. L. Y.NgK. H.IannucciA.MengoniA.et al. (2023). Species-specific gill’s microbiome of eight crab species with different breathing adaptations.Sci. Rep.13:21033. 10.1038/s41598-023-48308-w
5
BäckhedF. (2011). Programming of host metabolism by the gut microbiota.Ann. Nutr. Metab.5844–52. 10.1159/000328042
6
BanerjeeS.SchlaeppiK.van der HeijdenM. G. A. (2018). Keystone taxa as drivers of microbiome structure and functioning.Nat. Rev. Microbiol.16567–576. 10.1038/s41579-018-0024-1
7
BarthJ. A.PierceS. D.CarterB. R.ChanF.ErofeevA. Y.FisherJ. L.et al. (2024). Widespread and increasing near-bottom hypoxia in the coastal ocean off the United States Pacific Northwest.Sci. Rep.14:3798. 10.1038/s41598-024-54476-0
8
BreitburgD.LevinL. A.OschliesA.GrégoireM.ChavezF. P.ConleyD. J.et al. (2018). Declining oxygen in the global ocean and coastal waters.Science359:eaam7240. 10.1126/science.aam7240
9
BremnerJ.RogersS.FridC. (2006). Matching biological traits to environmental conditions in marine benthic ecosystems.J. Mar. Syst.60302–316. 10.1016/j.jmarsys.2006.02.004
10
CárcamoP. J.Hernández-MirandaE.VeasR.QuiñonesR. A. (2017). Macrofaunal community structure in Bahía Concepción (Chile) before and after the 8.8 Mw Maule mega-earthquake and tsunami.Mar. Environ. Res.130233–247. 10.1016/j.marenvres.2017.07.022
11
CaspiR.BillingtonR.FulcherC. A.KeselerI. M.KothariA.KrummenackerM.et al. (2018). The MetaCyc database of metabolic pathways and enzymes.Nucleic Acids Res.46D633–D639. 10.1093/nar/gkx935
12
CaspiR.BillingtonR.KeselerI. M.KothariA.KrummenackerM.MidfordP. E.et al. (2020). The MetaCyc database of metabolic pathways and enzymes - a 2019 update.Nucleic Acids Res.48D445–D453. 10.1093/nar/gkz862
13
ChanF.BarthJ. A.BlanchetteC. A.ByrneR. H.ChavezF.CheritonO.et al. (2017). Persistent spatial structuring of coastal ocean acidification in the California Current System.Sci. Rep.7:2526. 10.1038/s41598-017-02777-y
14
ChenJ.HuangJ.PengJ.YangC.LiaoY.LiJ.et al. (2022). Effects of hypoxic stress on the digestion, energy metabolism, oxidative stress regulation, and immune function of the pearl oyster (Pinctada fucata martensii).Aquaculture Rep.25:101246. 10.1016/j.aqrep.2022.101246
15
CodelloA.HoseG. C.CharitonA. (2022). Microbial co-occurrence networks as a biomonitoring tool for aquatic environments: A review.Mar. Freshwater Res.74409–422. 10.1071/MF22045
16
ConnollyN.CrosslandM.PearsonR. (2004). Effect of low dissolved oxygen on survival, emergence, and drift of tropical stream macroinvertebrates.J. North Am. Benthol. Soc.23251–270. 10.1899/0887-35932004023<0251:EOLDOO<2.0.CO;2
17
CsardiG.NepuszT. (2006). The igraph software package for complex network research.InterJ. Complex Syst.16951–9. 10.5281/zenodo.12578402
18
CurryK. D.WangQ.NuteM. G.TyshaievaA.ReevesE.SorianoS.et al. (2022). Emu: Species-level microbial community profiling of full-length 16S rRNA Oxford Nanopore sequencing data.Nat. Methods19845–853. 10.1038/s41592-022-01520-4
19
da Silva OliveiraF. A.ColaresG. B.HissaD. C.AngelimA. L.MeloV. M. M.LotufoT. M. C. (2013). Microbial epibionts of the colonial ascidians Didemnum galacteum and Cystodytes sp.Symbiosis5957–63. 10.1007/s13199-012-0210-2
20
DaiW.YeJ.LiuS.ChangG.XuH.LinZ.et al. (2022). Bacterial community dynamics in kumamoto oyster crassostrea sikamea hatchery during larval development.Front. Microbiol.13:933941. 10.3389/fmicb.2022.933941
21
DavieJ. R. (2003). Inhibition of histone deacetylase activity by butyrate.J. Nutr.133(7 Suppl), 2485S–2493S. 10.1093/jn/133.7.2485S
22
De La MazaL.FaríasL. (2023). The intensification of coastal hypoxia off central Chile: Long term and high frequency variability.Front. Earth Sci.10:929271. 10.3389/feart.2022.929271
23
DengD.HuS.LinZ.GengJ.QianZ.ZhangK.et al. (2024). High temperature aggravated hypoxia-induced intestine toxicity on juvenile Chinese mitten crab (Eriocheir sinensis).Comp. Biochem. Physiol. Part D: Genom. Proteom.52:101288. 10.1016/j.cbd.2024.101288
24
DiazR. J.RosenbergR. (2008). Spreading dead zones and consequences for marine ecosystems.Science321926–929. 10.1126/science.1156401
25
DixonP. (2003). VEGAN, a package of R functions for community ecology.J. Vegetation Sci.14927–930. 10.1658/1100-92332003014[0927:VAPORF]2.0.CO;2
26
DoneyS. C.BuschD. S.CooleyS. R.KroekerK. J. (2020). The impacts of ocean acidification on marine ecosystems and reliant human communities.Ann. Rev. Environ. Resources4583–112. 10.1146/annurev-environ-012320-083019
27
DouglasG. M.MaffeiV. J.ZaneveldJ. R.YurgelS. N.BrownJ. R.TaylorC. M.et al. (2020). PICRUSt2 for prediction of metagenome functions.Nat. Biotechnol.38685–688. 10.1038/s41587-020-0548-6
28
FanS.LiH.ZhaoR. (2020). Effects of normoxic and hypoxic conditions on the immune response and gut microbiota of Bostrichthys sinensis.Aquaculture525:735336. 10.1016/j.aquaculture.2020.735336
29
FosterZ. S.SharptonT. J.GrünwaldN. J. (2017). Metacoder: An R package for visualization and manipulation of community taxonomic diversity data.PLoS Comput. Biol.13:e1005404. 10.1371/journal.pcbi.1005404
30
FusiM.StephensonF.NavarreteS. A.TapiaF. J.LargierJ. L.MarascoR.et al. (2025). The ecology of the oxyscape in coastal ecosystems.Trends Ecol. Evol.40791–804. 10.1016/j.tree.2025.06.008
31
García-ReyesM.ThompsonS. A.Rogers-BennettL.SydemanW. J. (2022). Winter oceanographic conditions predict summer bull kelp canopy cover in northern California.PLoS One17:e0267737. 10.1371/journal.pone.0267737
32
GranthamB. A.ChanF.NielsenK. J.FoxD. S.BarthJ. A.HuyerA.et al. (2004). Upwelling-driven nearshore hypoxia signals ecosystem and oceanographic changes in the northeast Pacific.Nature429749–754. 10.1038/nature02605
33
GreijerA. E.Van der WallE. (2004). The role of hypoxia inducible factor 1 (HIF-1) in hypoxia induced apoptosis.J. Clin. Pathol.571009–1014. 10.1136/jcp.2003.015032
34
HallinS.BodelierP. L. E. (2020). Grand challenges in terrestrial microbiology: Moving on from a decade of progress in microbial biogeochemistry.Front. Microbiol.11:981. 10.3389/fmicb.2020.00981
35
HayeP. A.TuronX.SegoviaN. I. (2021). Time or space? Relative importance of geographic distribution and interannual variation in three lineages of the ascidian Pyura chilensis in the Southeast Pacific Coast.Front. Mar. Sci.8:657411. 10.3389/fmars.2021.657411
36
HernándezA.TapiaF. J. (2021). Connecting spatial structure in subtidal benthic communities with temporal variability in bottom temperature and dissolved oxygen along an upwelling coast.Estuarine Coastal Shelf Sci.250:107166. 10.1016/j.ecss.2021.107166
37
Hernández-MirandaE.VeasR.LabraF. A.SalamancaM.QuiñonesR. A. (2012). Response of the epibenthic macrofaunal community to a strong upwelling-driven hypoxic event in a shallow bay of the southern Humboldt Current System.Mar. Environ. Res.7916–28. 10.1016/j.marenvres.2012.04.004
38
JessenG. L.LichtschlagA.RametteA.PantojaS.RosselP. E.SchubertC. J.et al. (2017). Hypoxia causes preservation of labile organic matter and changes seafloor microbial community composition (Black Sea).Sci. Adv.3:e1601897. 10.1126/sciadv.1601897
39
JohnsonN. C.MarínC. (2025). Functional team selection as a framework for local adaptation in plants and their belowground microbiomes.ISME J.19:wraf137. 10.1093/ismejo/wraf137
40
KhanB.ClintonS. M.HampT. J.OliverJ. D.RingwoodA. H. (2018). Potential impacts of hypoxia and a warming ocean on oyster microbiomes.Mar. Environ. Res.13927–34. 10.1016/j.marenvres.2018.04.018
41
KirincichA. R.BarthJ. A.GranthamB. A.MengeB. A.LubchencoJ. (2005). Wind-driven inner-shelf circulation off central Oregon during summer.J. Geophys. Res. Oceans110C10S03.1–C10S03.17. 10.1029/2004JC002611
42
KnightR.VrbanacA.TaylorB. C.AksenovA.CallewaertC.DebeliusJ.et al. (2018). Best practices for analysing microbiomes.Nat. Rev. Microbiol.16410–422. 10.1038/s41579-018-0029-9
43
KurtzZ. D.MüllerC. L.MiraldiE. R.LittmanD. R.BlaserM. J.BonneauR. A. (2015). Sparse and compositionally robust inference of microbial ecological networks.PLoS Comput. Biol.11:e1004226. 10.1371/journal.pcbi.1004226
44
LagosM. E.BarnecheD. R.WhiteC. R.MarshallD. J. (2017). Do low oxygen environments facilitate marine invasions? Relative tolerance of native and invasive species to low oxygen conditions.Glob. Chang Biol.232321–2330. 10.1111/gcb.13668
45
LeeP.ChandelN. S.SimonM. C. (2020). Cellular adaptation to hypoxia through hypoxia inducible factors and beyond.Nat. Rev. Mol. Cell Biol.21268–283. 10.1038/s41580-020-0227-y
46
LeeY.ByeonE.KimD. H.MaszczykP.WangM.WuR. S. S.et al. (2023). Hypoxia in aquatic invertebrates: Occurrence and phenotypic and molecular responses.Aquat. Toxicol.263:106685. 10.1016/j.aquatox.2023.106685
47
LegerA.LeonardiT. (2019). pycoQC, interactive quality control for Oxford Nanopore Sequencing.J. Open Source Software4:1236. 10.21105/joss.01236
48
LevinL. A. (2003). Oxygen Minimum Zone Benthos: Adaptation and Community Response to Hypoxia, Oceanography and Marine Biology, An Annual Review.Boca Raton, FL: CRC Press.
49
LevinL.EkauW.GoodayA.JorissenF.MiddelburgJ.NaqviS.et al. (2009). Effects of natural and human-induced hypoxia on coastal benthos.Biogeosciences62063–2098. 10.5194/bgd-6-3563-2009c
50
LindsayE. C.MetcalfeN. B.LlewellynM. S. (2020). The potential role of the gut microbiota in shaping host energetics and metabolic rate.J. Anim. Ecol.892415–2426. 10.1111/1365-2656.13327
51
LiuY.WangJ.DingJ.ZhangY.HouC.ShenW.et al. (2024). Effects of hypoxia stress on oxidative stress, apoptosis and microorganisms in the intestine of large yellow croaker (Larimichthys crocea).Aquaculture581:740444. 10.1016/j.aquaculture.2023.740444
52
ManríquezP. H.CastillaJ. C. (2005). Self-fertilization as an alternative mode of reproduction in the solitary tunicate Pyura chilensis.Mar. Ecol. Prog. Ser.305113–125. 10.3354/meps305113
53
Montúfar-RomeroM.Valenzuela-MirandaD.Valenzuela-MuñozV.Morales-RiveraM. F.Gallardo-EscárateC. (2025). Microbiota dysbiosis in mytilus chilensis is induced by hypoxia, leading to molecular and functional consequences.Microorganisms13:825. 10.3390/microorganisms13040825
54
MuñozR.TapiaF. J.SobarzoM. (2025). Temporal variability of bottom hypoxia in open and semi-enclosed coastal areas in an upwelling region.Prog. Oceanography233:103435. 10.1016/j.pocean.2025.103435
55
OndovB. D.BergmanN. H.PhillippyA. M. (2011). Interactive metagenomic visualization in a Web browser.BMC Bioinform.12:385. 10.1186/1471-2105-12-385
56
ParksD. H.TysonG. W.HugenholtzP.BeikoR. G. (2014). STAMP: Statistical analysis of taxonomic and functional profiles.Bioinformatics303123–3124. 10.1093/bioinformatics/btu494
57
Pérez-ValdésM. P.Figueroa-AguileraD.Rojas-PérezC. (2017). Ciclo reproductivo de la ascidia Pyura chilensis (Urochordata: Ascidiacea) procedente de líneas de cultivo de mitílidos.Rev. Biol. Mar. Oceanografía52333–342. 10.4067/S0718-19572017000200012
58
RawlsJ. F.SamuelB. S.GordonJ. I. (2004). Gnotobiotic zebrafish reveal evolutionarily conserved responses to the gut microbiota.Proc. Natl. Acad. Sci. U. S. A.1014596–4601. 10.1073/pnas.0400706101
59
RohartF.GautierB.SinghA.Lê CaoK. A. (2017). mixOmics: An R package for ’omics feature selection and multiple data integration.PLoS Comput. Biol.13:e1005752. 10.1371/journal.pcbi.1005752
60
Santa CruzF.RetamalM. A. (2018). Biology, distribution, and abundance of Allopetrolisthes punctatus (Guérin, 1835)(Decapoda, Anomura, Porcellanidae), in the rocky subtidal of Chome, central Chile.Latin Am. J. Aquat. Res.46383–391. 10.3856/vol46-issue2-fulltext-14
61
SehnalL.Brammer-RobbinsE.WormingtonA. M.BlahaL.BisesiJ.LarkinI.et al. (2021). Microbiome composition and function in aquatic vertebrates: Small organisms making big impacts on aquatic animal health.Front. Microbiol.12:567408. 10.3389/fmicb.2021.567408
62
Sison-MangusM. P.MushegianA. A.EbertD. (2015). Water fleas require microbiota for survival, growth and reproduction.ISME J.959–67. 10.1038/ismej.2014.116
63
SobarzoM.SaldíasG. S.TapiaF. J.BravoL.MoffatC.LargierJ. L. (2016). On subsurface cooling associated with the Biobio River Canyon (Chile).J. Geophys. Res. Oceans1214568–4584. 10.1002/2016JC011796
64
SunF.WangY.WangC.ZhangL.TuK.ZhengZ. (2019). Insights into the intestinal microbiota of several aquatic organisms and association with the surrounding environment.Aquaculture507196–202. 10.1016/j.aquaculture.2019.04.026
65
SunS.JonesR. B.FodorA. A. (2020). Inference-based accuracy of metagenome prediction tools varies across sample types and functional categories.Microbiome8:46. 10.1186/s40168-020-00815-y
66
SydemanW. J.ThompsonS. A.García-ReyesM.KahruM.PetersonW. T.LargierJ. L. (2014). Multivariate ocean-climate indicators (MOCI) for the central California Current: Environmental change, 1990–2010.Prog. Oceanography120352–369. 10.1016/j.pocean.2013.10.017
67
TapiaF. J.LargierJ. L.CastilloM.WietersE. A.NavarreteS. A. (2014). Latitudinal discontinuity in thermal conditions along the nearshore of central-northern Chile.PLoS One9:e110841. 10.1371/journal.pone.0110841
68
TeelingH.FuchsB. M.BecherD.KlockowC.GardebrechtA.BennkeC. M.et al. (2012). Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom.Science336608–611. 10.1126/science.1218344
69
ToppingD. L.CliftonP. M. (2001). Short-chain fatty acids and human colonic function: Roles of resistant starch and nonstarch polysaccharides.Physiol. Rev.811031–1064. 10.1152/physrev.2001.81.3.1031
70
TremaroliV.BäckhedF. (2012). Functional interactions between the gut microbiota and host metabolism.Nature489242–249. 10.1038/nature11552
71
Ullah KhanF.ShangY.ChangX.KongH.ZuberiA.FangJ. K. H.et al. (2021). Effects of ocean acidification, hypoxia, and warming on the gut microbiota of the thick shell mussel mytilus coruscus through 16S rRNA gene sequencing.Front. Mar. Sci.8:736338. 10.3389/fmars.2021.736338
72
Vaquer-SunyerR.DuarteC. M. (2008). Thresholds of hypoxia for marine biodiversity.Proc. Natl. Acad. Sci. U. S. A.10515452–15457. 10.1073/pnas.0803833105
73
VivianiC. A.HillerA.WerdingB. (2010). Swarming in open space in the rocky intertidal: A new population-settlement strategy in the Eastern Pacific porcellanid crab, Allopetrolisthes punctatus (Decapoda, Anomura, Porcellanidae).Crustaceana83435–442. 10.1163/001121610X489395
74
WalterR. K.WoodsonC. B.LearyP. R.MonismithS. G. (2014). Connecting wind-driven upwelling and offshore stratification to nearshore internal bores and oxygen variability.J. Geophys. Res. Oceans1193517–3534. 10.1002/2014JC009998
75
WangR. X.HenenM. A.LeeJ. S.VögeliB.ColganS. P. (2021). Microbiota-derived butyrate is an endogenous HIF prolyl hydroxylase inhibitor.Gut Microbes13:1938380. 10.1080/19490976.2021.1938380
76
WangR. X.LeeJ. S.CampbellE. L.ColganS. P. (2020). Microbiota-derived butyrate dynamically regulates intestinal homeostasis through regulation of actin-associated protein synaptopodin.Proc. Natl. Acad. Sci. U. S. A.11711648–11657. 10.1073/pnas.1917597117
77
WangW.-Z.HuangJ.-S.ZhangJ.-D.WangZ.-L.LiH.-J.AmenyogbeE.et al. (2021). Effects of hypoxia stress on the intestinal microflora of juvenile of cobia (Rachycentron canadum).Aquaculture536:736419. 10.1016/j.aquaculture.2021.736419
78
WickR. R.JuddL. M.GorrieC. L.HoltK. E. (2017). Completing bacterial genome assemblies with multiplex MinION sequencing.Microb. Genom.3:e000132. 10.1099/mgen.0.000132
79
WickhamH.FrançoisR.HenryL.MüllerK.VaughanD. (2023). dplyr: A Grammar of Data Manipulation. R package version 1.1.4.San Francisco, CA: GitHub.
80
XieZ.LiY.XiongK.TuZ.WaihoK.YangC.et al. (2023). Combined effect of salinity and hypoxia on digestive enzymes and intestinal microbiota in the oyster Crassostrea hongkongensis.Environ. Pollut.331(Pt 2):121921. 10.1016/j.envpol.2023.121921
81
ZhengX.FuD.ChengJ.TangR.ChuM.ChuP.et al. (2021). Effects of hypoxic stress and recovery on oxidative stress, apoptosis, and intestinal microorganisms in Pelteobagrus vachelli.Aquaculture543:736945. 10.1016/j.aquaculture.2021.736945
82
ZingerL.Amaral-ZettlerL. A.FuhrmanJ. A.Horner-DevineM. C.HuseS. M.WelchD. B.et al. (2011). Global patterns of bacterial beta-diversity in seafloor and seawater ecosystems.PLoS One6:e24570. 10.1371/journal.pone.0024570
Summary
Keywords
16S rRNA, adaptation, crustacea, dissolved oxygen, microbiota, Pacific Ocean, tunicate
Citation
Valenzuela-Muñoz V, Tapia FJ, Gallardo-Escárate C, Mestre M, Navarrete SA, Gutiérrez MH, Morales-Rivera MF and Jessen GL (2026) Microbiota responses to environmental stress in mobile and sessile marine invertebrates: evidence for the effect of dissolved oxygen variations. Front. Microbiol. 17:1764313. doi: 10.3389/fmicb.2026.1764313
Received
09 December 2025
Revised
10 April 2026
Accepted
24 April 2026
Published
13 May 2026
Volume
17 - 2026
Edited by
Marcia Saraiva, University of Porto, Portugal
Reviewed by
Kae Kyoung Kwon, Korea Institute of Ocean Science and Technology (KIOST), Republic of Korea
Raphael Timotheus Koll, University of Hamburg, Germany
Updates
Copyright
© 2026 Valenzuela-Muñoz, Tapia, Gallardo-Escárate, Mestre, Navarrete, Gutiérrez, Morales-Rivera and Jessen.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Valentina Valenzuela-Muñoz, valevalenzuela@udec.clGerdhard L. Jessen, gerdhard.jessen@uach.cl
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.