Ocean plankton. Determinants of community structure in the global plankton interactome

Biogeographic patterns of meio- and micro-eukaryotic communities in dam-induced river-reservoir systems

Although the Three Gorges Dam (TGD) is the world's largest hydroelectric dam, little is known about the spatial-temporal patterns and community assembly mechanisms of meio- and micro-eukaryotes and its two subtaxa (zooplankton and zoobenthos). This knowledge gap is particularly evident across various habitats and during different water-level periods, primarily arising from the annual regular dam regulation. To address this inquiry, we employed mitochondrial cytochrome c oxidase I (COI) gene-based environmental DNA (eDNA) metabarcoding technology to systematically analyze the biogeographic pattern of the three communities within the Three Gorges Reservoir (TGR). Our findings reveal distinct spatiotemporal characteristics and complementary patterns in the distribution of meio- and micro-eukaryotes. The three communities showed similar biogeographic patterns and assembly processes. Notably, the diversity of these three taxa gradually decreased along the river. Their communities were less shaped by stochastic processes, which gradually decreased along the longitudinal riverine-transition-lacustrine gradient. Hence, deterministic factors, such as seasonality, environmental, and spatial variables, along with species interactions, likely play a pivotal role in shaping these communities. Environmental factors primarily drive seasonal variations in these communities, while hydrological conditions, represented as spatial distance, predominantly influence spatial variations. These three communities followed the distance-decay pattern. In winter, compared to summer, both the decay and species interrelationships are more pronounced. Taken together, this study offers fresh insights into the composition and diversity patterns of meio- and micro-eukaryotes at the spatial-temporal level. It also uncovers the mechanisms behind community assembly in various environmental niches within the dam-induced river-reservoir systems. KEY POINTS: • Distribution and diversity of meio- and micro-eukaryotes exhibit distinct spatiotemporal patterns in the TGR. • Contribution of stochastic processes in community assembly gradually decreases along the river. • Deterministic factors and species interactions shape meio- and micro-eukaryotic community.

Differences in archaeal diversity and potential ecological functions between saline and hypersaline lakes on Qinghai-Tibet Plateau were driven by multiple environmental and non-environmental factors beyond the salinity

BACKGROUND

Saline lakes are home to various archaea that play special and crucial roles in the global biogeochemical cycle. The Qinghai-Tibet Plateau hosts a large number of lakes with diverse salinity ranging from 0.1 to over 400 g/L, harboring complex and diverse archaea. To the best of our knowledge, the formation mechanisms and potential ecological roles of archaea in Qinghai-Tibetan Plateau saline lakes remain largely unknown.

RESULTS

Using High-throughput Illumina sequencing, we uncovered the vastly distinct archaea communities between two typical saline lakes with significant salinity differences on the Qinghai Tibet Plateau (Qinghai saline lake and Chaka hypersaline lake) and suggested archaea played different important roles in methanogenesis-related and nitrate reduction-related functions of these two lakes, respectively. Rather than the individual effect of salinity, the composite effect of salinity with diverse environmental parameters (e.g., temperature, chlorophyll a, total nitrogen, and total phosphorus) dominated the explanation of the variations in archaeal community structure in different habitats. Based on the network analysis, we further found the correlations between dominant archaeal OTUs were tight but significantly different between the two habitats, implying that archaeal interactions may also largely determine the shape of archaeal communities.

CONCLUSION

The present study improved our understanding of the structure and function of archaea in different saline lakes on the Qinghai-Tibet Plateau and provided a new perspective on the mechanisms underlying shaping their communities.

Enhancing soil health and nutrient cycling through soil amendments: Improving the synergy of bacteria and fungi

The synergy between bacteria and fungi is a key determinant of soil health and have a positive effect on plant development under drought conditions, with the potentially enhancing the sustainability of amending soil with natural materials. However, identifying how soil amendments influence plant growth is often difficult due to the complexity of microorganisms and their links with different soil amendment types and environmental factors. To address this, we conducted a field experiment to examine the impact of soil amendments (biochar, Bacillus mucilaginosus, Bacillus subtilis and super absorbent polymer) on plant growth. We also assessed variations in microbial community, links between fungi and bacteria, and soil available nutrients, while exploring how the synergistic effects between fungus and bacteria influenced the response of soil amendments to plant growth. This study revealed that soil amendments reduced soil bacterial diversity but increased the proportion of the family Enterobacteriaceae, Nitrosomonadaceae, and also increased soil fungal diversity and the proportion of the sum of the family Lasiosphaeriaceae, Chaetomiaceae, Pleosporaceae. Changes in soil microbial communities lead to increase the complexity of microbial co-occurrence networks. Furthermore, this heightened network complexity enhanced the synergy of soil bacteria and fungi, supporting bacterial functions related to soil nutrient cycling, such as metabolic functions and genetic, environmental, and cellular processes. Hence, the BC and BS had 3.0-fold and 0.5-fold greater root length densities than CK and apple tree shoot growth were increased by 62.14 %,50.53 % relative to CK, respectively. In sum, our results suggest that the synergistic effect of bacteria and fungi impacted apple tree growth indirectly by modulating soil nutrient cycling. These findings offer a new strategy for enhancing the quality of arable land in arid and semi-arid regions.

Using Network Analysis and Predictive Functional Analysis to Explore the Fluorotelomer Biotransformation Potential of Soil Microbial Communities

Microbial transformation of per- and polyfluoroalkyl substances (PFAS), including fluorotelomer-derived PFAS, by native microbial communities in the environment has been widely documented. However, few studies have identified the key microorganisms and their roles during the PFAS biotransformation processes. This study was undertaken to gain more insight into the structure and function of soil microbial communities that are relevant to PFAS biotransformation. We collected 16S rRNA gene sequencing data from 8:2 fluorotelomer alcohol and 6:2 fluorotelomer sulfonate biotransformation studies conducted in soil microcosms under various redox conditions. Through co-occurrence network analysis, several genera, including Variovorax, Rhodococcus, and Cupriavidus, were found to likely play important roles in the biotransformation of fluorotelomers. Additionally, a metagenomic prediction approach (PICRUSt2) identified functional genes, including 6-oxocyclohex-1-ene-carbonyl-CoA hydrolase, cyclohexa-1,5-dienecarbonyl-CoA hydratase, and a fluoride-proton antiporter gene, that may be involved in defluorination. This study pioneers the application of these bioinformatics tools in the analysis of PFAS biotransformation-related sequencing data. Our findings serve as a foundational reference for investigating enzymatic mechanisms of microbial defluorination that may facilitate the development of efficient microbial consortia and/or pure microbial strains for PFAS biotransformation.

Biogeographical and biodiversity patterns of planktonic microeukaryotes along the tropical western to eastern Pacific Ocean transect revealed by metabarcoding

Microeukaryotic plankton (0.2-200 µm), which are morphologically and genetically highly diverse, play a crucial role in ocean productivity and carbon consumption. The Pacific Ocean (PO), one of the world's largest oligotrophic regions, remains largely unexplored in terms of the biogeography and biodiversity of microeukaryotes based on large-scale sampling. We investigated the horizontal distribution of microeukaryotes along a 16,000 km transect from the west to the east of the PO. The alpha diversity indices showed a distinct decreasing trend from west to east, which was highly correlated with water temperature. The microeukaryotic community, which was clustered into the western, central, and eastern PO groups, displayed a significant distance-decay relationship. Syndiniales, a lineage of parasitic dinoflagellates, was ubiquitously distributed along the transect and dominated the community in terms of both sequence and zero-radius operational taxonomic unit (ZOTU) proportions. The prevailing dominance of Syndiniales-affiliated ZOTUs and their close associations with dinoflagellates, diatoms, and radiolarians, as revealed by SparCC correlation analysis, suggested that parasitism may be an important trophic strategy in the surface waters of the PO. Geographical distance and temperature were the most important environmental factors that significantly correlated with community structure. Overall, our study sheds more light on the distribution pattern of both alpha and beta diversities of microeukaryotic communities and highlighted the importance of parasitisms by Syndiniales across the tropical PO.IMPORTANCEUnderstanding the biogeographical and biodiversity patterns of microeukaryotic communities is essential to comprehending their roles in biogeochemical cycling. In this study, planktonic microeukaryotes were collected along a west-to-east Pacific Ocean transect (ca. 16,000 km). Our study revealed that the alpha diversity indices were highly correlated with water temperature, and the microeukaryotic communities displayed a distinct geographical distance-driven pattern. The predominance of the parasitic dinoflagellate lineage Syndiniales and their close relationship with other microeukaryotic groups suggest that parasitism may be a crucial survival strategy for microeukaryotes in the surface waters of the Pacific Ocean. Our findings expand our understanding of the biodiversity and biogeographical pattern of microeukaryotes and highlight the significance of parasitic Syndiniales in the surface ocean.

New plastids, old proteins: repeated endosymbiotic acquisitions in kareniacean dinoflagellates

Dinoflagellates are a diverse group of ecologically significant micro-eukaryotes that can serve as a model system for plastid symbiogenesis due to their susceptibility to plastid loss and replacement via serial endosymbiosis. Kareniaceae harbor fucoxanthin-pigmented plastids instead of the ancestral peridinin-pigmented ones and support them with a diverse range of nucleus-encoded plastid-targeted proteins originating from the haptophyte endosymbiont, dinoflagellate host, and/or lateral gene transfers (LGT). Here, we present predicted plastid proteomes from seven distantly related kareniaceans in three genera (Karenia, Karlodinium, and Takayama) and analyze their evolutionary patterns using automated tree building and sorting. We project a relatively limited ( ~ 10%) haptophyte signal pointing towards a shared origin in the family Chrysochromulinaceae. Our data establish significant variations in the functional distributions of these signals, emphasizing the importance of micro-evolutionary processes in shaping the chimeric proteomes. Analysis of plastid genome sequences recontextualizes these results by a striking finding the extant kareniacean plastids are in fact not all of the same origin, as two of the studied species (Karlodinium armiger, Takayama helix) possess plastids from different haptophyte orders than the rest.

Differential effects of warming on the complexity and stability of the microbial network in Phragmites australis and Spartina alterniflora wetlands in Yancheng, Jiangsu Province, China

The impact of climate warming on soil microbial communities can significantly influence the global carbon cycle. Coastal wetlands, in particular, are susceptible to changes in soil microbial community structure due to climate warming and the presence of invasive plant species. However, there is limited knowledge about how native and invasive plant wetland soil microbes differ in their response to warming. In this study, we investigated the temporal dynamics of soil microbes (prokaryotes and fungi) under experimental warming in two coastal wetlands dominated by native Phragmites australis (P. australis) and invasive Spartina alterniflora (S. alterniflora). Our research indicated that short-term warming had minimal effects on microbial abundance, diversity, and composition. However, it did accelerate the succession of soil microbial communities, with potentially greater impacts on fungi than prokaryotes. Furthermore, in the S. alterniflora wetland, experimental warming notably increased the complexity and connectivity of the microbial networks. While in the P. australis wetland, it decreased these factors. Analysis of robustness showed that experimental warming stabilized the co-occurrence network of the microbial community in the P. australis wetland, but destabilized it in the S. alterniflora wetland. Additionally, the functional prediction analysis using the Faprotax and FunGuild databases revealed that the S. alterniflora wetland had a higher proportion of saprotrophic fungi and prokaryotic OTUs involved in carbon degradation (p < 0.05). With warming treatments, there was an increasing trend in the proportion of prokaryotic OTUs involved in carbon degradation, particularly in the S. alterniflora wetland. Therefore, it is crucial to protect native P. australis wetlands from S. alterniflora invasion to mitigate carbon emissions and preserve the health of coastal wetland ecosystems under future climate warming in China.

Herptile gut microbiomes: a natural system to study multi-kingdom interactions between filamentous fungi and bacteria

Reptiles and amphibians (herptiles) are some of the most endangered and threatened species on the planet and numerous conservation strategies are being implemented with the goal of ensuring species recovery. Little is known, however, about the gut microbiome of wild herptiles and how it relates to the health of these populations. Here, we report results from the gut microbiome characterization of both a broad survey of herptiles, and the correlation between the fungus Basidiobolus, and the bacterial community supported by a deeper, more intensive sampling of Plethodon glutinosus, known as slimy salamanders. We demonstrate that bacterial communities sampled from frogs, lizards, and salamanders are structured by the host taxonomy and that Basidiobolus is a common and natural component of these wild gut microbiomes. Intensive sampling of multiple hosts across the ecoregions of Tennessee revealed that geography and host:geography interactions are strong predictors of distinct Basidiobolus operational taxonomic units present within a given host. Co-occurrence analyses of Basidiobolus and bacterial community diversity support a correlation and interaction between Basidiobolus and bacteria, suggesting that Basidiobolus may play a role in structuring the bacterial community. We further the hypothesis that this interaction is advanced by unique specialized metabolism originating from horizontal gene transfer from bacteria to Basidiobolus and demonstrate that Basidiobolus is capable of producing a diversity of specialized metabolites including small cyclic peptides.IMPORTANCEThis work significantly advances our understanding of biodiversity and microbial interactions in herptile microbiomes, the role that fungi play as a structural and functional members of herptile gut microbiomes, and the chemical functions that structure microbiome phenotypes. We also provide an important observational system of how the gut microbiome represents a unique environment that selects for novel metabolic functions through horizontal gene transfer between fungi and bacteria. Such studies are needed to better understand the complexity of gut microbiomes in nature and will inform conservation strategies for threatened species of herpetofauna.

Deviation from neutral species abundance distributions unveils geographical differences in the structure of diatom communities

In recent years, the application of metagenomics techniques has advanced our understanding of plankton communities and their global distribution. Despite this progress, the relationship between the abundance distribution of diatom species and varying marine environmental conditions remains poorly understood. This study, leveraging data from the Tara Oceans expedition, tests the hypothesis that diatoms in sampled stations display a consistent species abundance distribution structure, as though they were sampled from a single ocean-wide metacommunity. Using a neutral sampling theory, we thus develop a framework to estimate the structure and diversity of diatom communities at each sampling station given the shape of the species abundance distribution of the metacommunity and the information of a reference station. Our analysis reveals a substantial temperature gradient in the discrepancies between predicted and observed biodiversity across the sampled stations. These findings challenge the hypothesis of a single neutral metacommunity, indicating that environmental differences substantially influence both the composition and structure of diatom communities.

Algal blooms in the ocean: hot spots for chemically mediated microbial interactions

The cycling of major nutrients in the ocean is affected by large-scale phytoplankton blooms, which are hot spots of microbial life. Diverse microbial interactions regulate bloom dynamics. At the single-cell level, interactions between microorganisms are mediated by small molecules in the chemical crosstalk that determines the type of interaction, ranging from mutualism to pathogenicity. Algae interact with viruses, bacteria, parasites, grazers and other algae to modulate algal cell fate, and these interactions are dependent on the environmental context. Recent advances in mass spectrometry and single-cell technologies have led to the discovery of a growing number of infochemicals - metabolites that convey information - revealing the ability of algal cells to govern biotic interactions in the ocean. The diversity of infochemicals seems to account for the specificity in cellular response during microbial communication. Given the immense impact of algal blooms on biogeochemical cycles and climate regulation, a major challenge is to elucidate how microscale interactions control the fate of carbon and the recycling of major elements in the ocean. In this Review, we discuss microbial interactions and the role of infochemicals in algal blooms. We further explore factors that can impact microbial interactions and the available tools to decipher them in the natural environment.

Spatiotemporal Changes of Antibiotic Resistance, Potential Pathogens, and Health Risk in Kindergarten Dust

The microorganisms present in kindergartens are extremely important for children's health during their three-year preschool education. To assess the risk of outdoor dust in kindergartens, the antibiotic resistome and potential pathogens were investigated in dust samples collected from 59 kindergartens in Xiamen, southeast China in both the winter and summer. Both high-throughput quantitative PCR and metagenome analysis revealed a higher richness and abundance of antibiotic resistance genes (ARGs) in winter (P < 0.05). Besides, the bloom of ARGs and potential pathogens was evident in the urban kindergartens. The co-occurrence patterns among ARGs, mobile genetic elements (MGEs), and potential pathogens suggested some bacterial pathogens were potential hosts of ARGs and MGEs. We found a large number of high-risk ARGs in the dust; the richness and abundance of high-risk ARGs were higher in winter and urban kindergartens compared to in summer and peri-urban kindergartens, respectively. The results of the co-occurrence patterns and high-risk ARGs jointly reveal that urbanization will significantly increase the threat of urban dust to human beings and their risks will be higher in winter. This study unveils the close association between ARGs/mobile ARGs and potential pathogens and emphasizes that we should pay more attention to the health risks induced by their combination.

Positive associations fuel soil biodiversity and ecological networks worldwide

Microbial interactions are key to maintaining soil biodiversity. However, whether negative or positive associations govern the soil microbial system at a global scale remains virtually unknown, limiting our understanding of how microbes interact to support soil biodiversity and functions. Here, we explored ecological networks among multitrophic soil organisms involving bacteria, protists, fungi, and invertebrates in a global soil survey across 20 regions of the planet and found that positive associations among both pairs and triads of soil taxa governed global soil microbial networks. We further revealed that soil networks with greater levels of positive associations supported larger soil biodiversity and resulted in lower network fragility to withstand potential perturbations of species losses. Our study provides unique evidence of the widespread positive associations between soil organisms and their crucial role in maintaining the multitrophic structure of soil biodiversity worldwide.

Depth significantly affects plastisphere microbial evenness, assembly and co-occurrence pattern but not richness and composition

Microplastics have become one of the hot concerns of global marine pollution. In recent years, diversity and abiotic influence factors of plastisphere microbial communities were well documented, but our knowledge of their assembly mechanisms and co-occurrence patterns remains unclear, especially the effects of depth on them. Here, we collected microorganisms on microplastics to investigate how ocean depth affects on microbial diversity, community composition, assembly processes and co-occurrence patterns. Our results indicated that there were similar microbial richness and community compositions but microbial evenness and unique microbes were obviously different in different ocean layers. Our findings also demonstrated that deterministic processes played dominant roles in the assembly of the mesopelagic plastisphere microbial communities, while the bathypelagic microbial community assembly was mainly shaped by stochastic processes. In addition, the co-occurrence networks suggested that the relationships between microorganisms in the mesopelagic layer were more complex and stable than those in the bathypelagic layer. Simultaneously, we also found that Proteobacteria and Actinobacteriota were the most abundant keystones which played important roles in microbial co-occurrence networks at both layers. This study enhanced our understanding of microbial diversity, assembly mechanism, and co-occurrence pattern on plastisphere surfaces, and provided useful insights into microorganisms capable of degrading plastics and microbial remediation.

Dynamics of bacterioplankton communities in the estuary areas of the Taihu Lake: Distinct ecological mechanisms of abundant and rare communities

As the crucial confluences of rivers and lakes, the estuary areas with varied hydrodynamic exchanges intensively affect the bacterioplankton communities, whereas the ecological characteristics of the bacterioplankton in the areas have not been well understood. Here, the distribution patterns and assembly mechanisms of bacterioplankton communities in the estuary areas of the Taihu Lake were investigated using high-throughput sequencing and multivariate statistical analyses. Our results showed obvious seasonal variations in bacterioplankton diversity and community composition, which had significant correlations with water temperature. Neutral and null models together revealed that stochastic processes (especially dispersal limitation) were the major processes in shaping the communities across different seasons. By contrast, heterogeneous selection in deterministic processes exhibited increased impacts on community assembly during summer and autumn, which was significantly related to the comprehensive water quality index (WQI) rather than any single factor. In this study, rare communities displayed more pronounced seasonal dynamics compared to abundant communities, likely due to their sensitivity towards environmental factors. Accordingly, the heterogeneous selection of deterministic processes largely shaped the rare communities. These results enriched our understanding of the assembly mechanisms of bacterioplankton communities in estuary areas and emphasized the specific co-occurrence patterns of abundant and rare communities.

An Assessment of Human Opportunistic Pathogenic Bacteria on Daily Necessities in Nanjing City during Plum Rain Season

The plum rain season is a special climatic phenomenon in east Asia, which is characterized by persistent rainfall, a high temperature, and humidity, providing suitable environmental conditions for certain pathogenic bacteria, thus increasing the incidence of respiratory, gastrointestinal, and urinary diseases. However, studies on human opportunistic pathogenic bacteria communities during the plum rain season are still limited. In this study, the characteristics of human opportunistic pathogenic bacterial communities on daily necessities during the non-plum and plum rain seasons were investigated using high-throughput sequencing technology. The results revealed that the relative abundance of human opportunistic pathogenic bacteria was higher in the plum rain season (cotton cloth: 2.469%, electric bicycles: 0.724%, rice: 3.737%, and washbasins: 5.005%) than in the non-plum rain season (cotton cloth: 1.425%, electric bicycles: 0.601%, rice: 2.426%, and washbasins: 4.801%). Both temperature and relative humidity affected human opportunistic pathogenic bacterial communities. Stochastic processes dominated the assembly process of human opportunistic pathogenic bacterial communities, and undominated processes prevailed. The stability of the co-occurrence network was higher in the non-plum rain season than that in the plum rain season. In addition, the proportion of deterministic processes showed the same trend as the complexity of the co-occurrence network.

Depth shapes microbiome assembly and network stability in the Mariana Trench

IMPORTANCE

Exploring microbial interactions and their stability/resilience from the surface to the hadal ocean is critical for further understanding of the microbiome structure and ecosystem function in the Mariana Trench. Vertical gradients did not destabilize microbial communities after long-term evolution and adaption. The uniform niche breadth, diversity, community complexity, and stability of microbiomes in both upper bathypelagic and hadal waters suggest the consistent roles of microbiomes in elemental cycling and adaptive strategies to overcome extreme environmental conditions. Compared with microeukaryotes, bacteria and archaea play a pivotal role in shaping the stability of the hadal microbiome. The consistent co-occurrence stability of microbiomes across vertical gradients was observed in the Mariana Trench. These results illuminate a key principle of microbiomes inhabiting the deepest trench: although distinct microbial communities occupy specific habitats, the interactions within microbial communities remain consistently stable from the upper bathypelagic to the hadal waters.

Structure and assembly process of fungal communities in the Yangtze River Estuary

Marine fungi are essential for the ecological function of estuarine ecosystems. However, limited studies have reported on the structure and assembly pattern of the fungal communities in estuaries. The purpose of this study is to reveal the structure and the ecological process of the fungal community in the Yangtze River Estuary (YRE) by using the amplicon sequencing method. Phyla of Ascomycota, Basidiomycota, and Chytridiomycota were dominant in the seawater and sediment samples from YRE. The null model analysis, community-neutral community model (NCM), and phylogenetic normalized stochasticity ratio (pNST) showed that the stochastic process dominated the assembly of fungal communities in YRE. Drift and homogeneous dispersal were the predominant stochastic processes for the fungal community assembly in seawater and sediment samples, respectively. The co-occurrence network analysis showed that fungal communities were more complex and closely connected in the sediment than in the seawater samples. Phyla Ascomycota, Basidiomycota, and Mucoromycota were the potential keystone taxa in the network. These findings demonstrated the importance of stochastic processes for the fungal community assembly, thereby widening our knowledge of the community structure and dynamics of fungi for future study and utilization in the YRE ecosystem.

Disentangling microbial networks across pelagic zones in the tropical and subtropical global ocean

Microbial interactions are vital in maintaining ocean ecosystem function, yet their dynamic nature and complexity remain largely unexplored. Here, we use association networks to investigate possible ecological interactions in the marine microbiome among archaea, bacteria, and picoeukaryotes throughout different depths and geographical regions of the tropical and subtropical global ocean. Our findings reveal that potential microbial interactions change with depth and geographical scale, exhibiting highly heterogeneous distributions. A few potential interactions were global, meaning they occurred across regions at the same depth, while 11-36% were regional within specific depths. The bathypelagic zone had the lowest proportion of global associations, and regional associations increased with depth. Moreover, we observed that most surface water associations do not persist in deeper ocean layers despite microbial vertical dispersal. Our work contributes to a deeper understanding of the tropical and subtropical global ocean interactome, which is essential for addressing the challenges posed by global change.

Deciphering microbiomes dozens of meters under our feet and their edaphoclimatic and spatial drivers

Microbes inhabiting deep soil layers are known to be different from their counterpart in topsoil yet remain under investigation in terms of their structure, function, and how their diversity is shaped. The microbiome of deep soils (>1 m) is expected to be relatively stable and highly independent from climatic conditions. Much less is known, however, on how these microbial communities vary along climate gradients. Here, we used amplicon sequencing to investigate bacteria, archaea, and fungi along fifteen 18-m depth profiles at 20-50-cm intervals across contrasting aridity conditions in semi-arid forest ecosystems of China's Loess Plateau. Our results showed that bacterial and fungal α diversity and bacterial and archaeal community similarity declined dramatically in topsoil and remained relatively stable in deep soil. Nevertheless, deep soil microbiome still showed the functional potential of N cycling, plant-derived organic matter degradation, resource exchange, and water coordination. The deep soil microbiome had closer taxa-taxa and bacteria-fungi associations and more influence of dispersal limitation than topsoil microbiome. Geographic distance was more influential in deep soil bacteria and archaea than in topsoil. We further showed that aridity was negatively correlated with deep-soil archaeal and fungal richness, archaeal community similarity, relative abundance of plant saprotroph, and bacteria-fungi associations, but increased the relative abundance of aerobic ammonia oxidation, manganese oxidation, and arbuscular mycorrhizal in the deep soils. Root depth, complexity, soil volumetric moisture, and clay play bridging roles in the indirect effects of aridity on microbes in deep soils. Our work indicates that, even microbial communities and nutrient cycling in deep soil are susceptible to changes in water availability, with consequences for understanding the sustainability of dryland ecosystems and the whole-soil in response to aridification. Moreover, we propose that neglecting soil depth may underestimate the role of soil moisture in dryland ecosystems under future climate scenarios.

Role of Syndiniales parasites in depth-specific networks and carbon flux in the oligotrophic ocean

Microbial associations that result in phytoplankton mortality are important for carbon transport in the ocean. This includes parasitism, which in microbial food webs is dominated by the marine alveolate group, Syndiniales. Parasites are expected to contribute to carbon recycling via host lysis; however, knowledge on host dynamics and correlation to carbon export remain unclear and limit the inclusion of parasitism in biogeochemical models. We analyzed a 4-year 18S rRNA gene metabarcoding dataset (2016-19), performing network analysis for 12 discrete depths (1-1000 m) to determine Syndiniales-host associations in the seasonally oligotrophic Sargasso Sea. Analogous water column and sediment trap data were included to define environmental drivers of Syndiniales and their correlation with particulate carbon flux (150 m). Syndiniales accounted for 48-74% of network edges, most often associated with Dinophyceae and Arthropoda (mainly copepods) at the surface and Rhizaria (Polycystinea, Acantharea, and RAD-B) in the aphotic zone. Syndiniales were the only eukaryote group to be significantly (and negatively) correlated with particulate carbon flux, indicating their contribution to flux attenuation via remineralization. Examination of Syndiniales amplicons revealed a range of depth patterns, including specific ecological niches and vertical connection among a subset (19%) of the community, the latter implying sinking of parasites (infected hosts or spores) on particles. Our findings elevate the critical role of Syndiniales in marine microbial systems and reveal their potential use as biomarkers for carbon export.

Sea-ice melt determines seasonal phytoplankton dynamics and delimits the habitat of temperate Atlantic taxa as the Arctic Ocean atlantifies

The Arctic Ocean is one of the regions where anthropogenic environmental change is progressing most rapidly and drastically. The impact of rising temperatures and decreasing sea ice on Arctic marine microbial communities is yet not well understood. Microbes form the basis of food webs in the Arctic Ocean, providing energy for larger organisms. Previous studies have shown that Atlantic taxa associated with low light are robust to more polar conditions. We compared to which extent sea ice melt influences light-associated phytoplankton dynamics and biodiversity over two years at two mooring locations in the Fram Strait. One mooring is deployed in pure Atlantic water, and the second in the intermittently ice-covered Marginal Ice Zone. Time-series analysis of amplicon sequence variants abundance over a 2-year period, allowed us to identify communities of co-occurring taxa that exhibit similar patterns throughout the annual cycle. We then examined how alterations in environmental conditions affect the prevalence of species. During high abundance periods of diatoms, polar phytoplankton populations dominated, while temperate taxa were weakly represented. Furthermore, we found that polar pelagic and ice-associated taxa, such as Fragilariopsis cylindrus and Melosira arctica, were more common in Atlantic conditions, while temperate taxa, such as Odontella aurita and Proboscia alata, were less abundant under polar conditions. This suggests that sea ice melt may act as a barrier to the northward expansion of temperate phytoplankton, preventing their dominance in regions still strongly influenced by polar conditions. Our findings highlight the complex interactions between sea ice melt, phytoplankton dynamics, and biodiversity in the Arctic.

Environmental heterogeneity associated with boat activity shapes bacteria and microeukaryotic communities with discrepant response patterns

With the development of global tourism, tourist boats, a significant form of anthropogenic disturbance, are having an increasingly serious impact on the structure and function of aquatic ecosystems. In this study, the effects of different intensities of tourist boat activities on the microbial communities of West lake, were investigated by high-throughput sequencing. The results showed significant differences in the composition of bacterioplankton and microeukaryotic communities between the high-intensity boat activity (HIBA) area and low-intensity boat activity (LIBA) area. Variation partitioning analysis showed that environmental factors contributed the most to microbial community variation, and the effect of boat activities on microbial communities mainly occurred through coupling with environmental factors. The contribution of boat activity to microbial community changes occupies the second place, the first being environmental factors. Co-occurrence network analyses showed that microbial communities in the HIBA area had more nodes and edges, higher connectivity and lower modularity than in the LIBA area, suggesting a more complex and stable network. Networks of associations between potential keystone taxa and environmental factors reveal the way in which boat activity affects microbial communities. The bacterial community responded strongly to environmental factors associated with boat activities, whereas the microeukaryotic community was more likely to be regulated by interspecific interactions. This also suggests that when faced with disturbances from the boat activity, microeukaryotes might exert a stronger direct resistance effect compared to bacterioplankton. These findings imply that bacterioplankton and microeukaryotes demonstrate distinct response patterns in the presence of disturbance caused by boat activity. Our research expand our understanding of the effects of boat activities on aquatic ecosystems and provide further insights into the assessment of anthropogenic disturbances in aquatic ecosystems.

Globally-distributed microbial eukaryotes exhibit endemism at deep-sea hydrothermal vents

Single-celled microbial eukaryotes inhabit deep-sea hydrothermal vent environments and play critical ecological roles in the vent-associated microbial food web. 18S rRNA amplicon sequencing of diffuse venting fluids from four geographically- and geochemically-distinct hydrothermal vent fields was applied to investigate community diversity patterns among protistan assemblages. The four vent fields include Axial Seamount at the Juan de Fuca Ridge, Sea Cliff and Apollo at the Gorda Ridge, all in the NE Pacific Ocean, and Piccard and Von Damm at the Mid-Cayman Rise in the Caribbean Sea. We describe species diversity patterns with respect to hydrothermal vent field and sample type, identify putative vent endemic microbial eukaryotes, and test how vent fluid geochemistry may influence microbial community diversity. At a semi-global scale, microbial eukaryotic communities at deep-sea vents were composed of similar proportions of dinoflagellates, ciliates, Rhizaria, and stramenopiles. Individual vent fields supported distinct and highly diverse assemblages of protists that included potentially endemic or novel vent-associated strains. These findings represent a census of deep-sea hydrothermal vent protistan communities. Protistan diversity, which is shaped by the hydrothermal vent environment at a local scale, ultimately influences the vent-associated microbial food web and the broader deep-sea carbon cycle.

Land use conversion increases network complexity and stability of soil microbial communities in a temperate grassland

Soils harbor highly diverse microbial communities that are critical to soil health, but agriculture has caused extensive land use conversion resulting in negative effects on critical ecosystem processes. However, the responses and adaptations of microbial communities to land use conversion have not yet been understood. Here, we examined the effects of land conversion for long-term crop use on the network complexity and stability of soil microbial communities over 19 months. Despite reduced microbial biodiversity in comparison with native tallgrass prairie, conventionally tilled (CT) cropland significantly increased network complexity such as connectivity, connectance, average clustering coefficient, relative modularity, and the number of species acting at network hubs and connectors as well as resulted in greater temporal variation of complexity indices. Molecular ecological networks under CT cropland became significantly more robust and less vulnerable, overall increasing network stability. The relationship between network complexity and stability was also substantially strengthened due to land use conversion. Lastly, CT cropland decreased the number of relationships between network structure and environmental properties instead being strongly correlated to management disturbances. These results indicate that agricultural disturbance generally increases the complexity and stability of species "interactions", possibly as a trade-off for biodiversity loss to support ecosystem function when faced with frequent agricultural disturbance.

Influence of phycospheric bacterioplankton disruption or removal on algae growth and survival

Phycospheric bacteria play a crucial role in the survival of microalgae. However, the potential of using the growth regulation and community structure modulation of phycospheric bacteria to prevent the occurrence of blooms is yet to be verified. The phycospheric bacterioplankton of Cyclotella sp. can be categorized into HNA (high nucleic acid) bacteria and LNA (low nucleic acid) bacteria. 16S rRNA sequencing showed that the HNA bacteria exhibited higher α-diversity compared to the LNA bacteria, and the microbial community composition also exhibited variations. Metagenomic sequencing further indicated the distinct ecological functions between HNA and LNA bacteria. Furthermore, the study showcased the restorative capacity of the phycospheric bacterioplankton. Biomass analysis revealed that the recovery of phycospheric bacterioplankton positively influenced the microalgae growth, thus affirming the significance of phycospheric bacterioplankton to microalgae. The community structure of phycospheric bacterioplankton demonstrated a notable decrease in the abundance of restored LNA core bacteria. Additionally, the restored phycospheric bacterioplankton exhibited a more complex co-occurrence network structure, resulting in decreased resistance and sensitivity of microalgae to adverse environments. The presence of phycospheric bacterioplankton provides a protective shield for microalgae, and thus destabilizing or removing phycospheric bacterioplankton may effectively inhibit growth of microalgae.

Seasonal dynamics of bacterial community and co-occurrence with eukaryotic phytoplankton in the Pearl River Estuary

In this study, we investigated the taxonomic composition of the bacteria and phytoplankton communities in the Pearl River Estuary (PRE) through Illumina sequencing of the V3-V4 region of the 16 S rRNA gene. Furthermore, their relationships as well as recorded environmental variables were explored by co-occurrence networks. Bacterial community composition was different in two size fractions, as well as along the salinity gradient across two seasons. Free-living (FL) communities were dominated by pico-sized Cyanobacteria (Synechococcus CC9902) while Exiguobacterium, Halomonas and Pseudomonas were predominantly associated with particle-associated (PA) lifestyle, and Cyanobium PCC-6307 exhibited seasonal shifts in lifestyles in different seasons. In wet season, bacterial community composition was characterized by abundance of Cyanobacteria, Actinobacteria, and Bacteroidetes, which were tightly linked with high riverine inflow. While in dry season, Proteobacteria increased in prevalence, especially for Psychrobacter, NOR5/OM60 clade and Pseudomonas, which were thrived in lower water temperature and higher salinity. Moreover, we discovered that differences between PA and FL composition were more significant in the wet season than in the dry season, which may be due to better nutritional conditions of particles (indicated by POC%) in the wet season and then attract more diverse PA populations. Based on the analysis of plastidial 16 S rRNA genes, abundant small-sized mixotrophic phytoplankton (Dinophyceae, Euglenida and Haptophyta) were identified in the PRE. The complexity of co-occurrence network increased from FL to PA fractions in both seasons, which suggested that suspended particles can provide ecological niches for particle-associated colonizers contributing to the maintenance of a more stable community structure. In addition, the majority of phytoplankton species exhibited positive co-occurrences with both other phytoplankton species and bacterial counterparts, indicating the mutual cooperation between phytoplankton assemblages and specific bacterial populations e likely benefited from phytoplankton-derived organic compounds. This study enhances our understanding of the seasonal and spatial dynamics of bacterial communities and their potential relationship with phytoplankton assembly in estuarine waters.

Patterns of bacterial distance decay and community assembly in different land-use types as influenced by tillage management and soil layers

Land use and cover change are major factors driving global change and greatly impact terrestrial organisms, especially soil microbial diversity. Little is known, however, about bacterial diversity, distribution patterns and assembly processes across different land use types. In this study, therefore, we conducted a large-scale field survey of 48 sampling sites, encompassing different land use types in Xuancheng city, China, with different degrees of soil disturbance and different soil horizons. The distance-decay relationships (DDRs), assembly processes and the spatial patterns of soil bacterial communities were investigated based on high-throughput sequencing data. We found that the DDRs might be weakened by anthropogenic disturbances, which were not observed in tilled soils, while a decreasing trend was observed along the soil horizon in untilled soils. The relative importance of environmental factors and geographic distance varied with soil tillage. Specifically, bacterial communities in tilled soils were driven by non-spatially autocorrelated environmental factors, while untilled soils were more susceptible to geographic distance. In addition, the heterogeneity of soil properties, as well as the differences in soil bacterial niche width and niche overlap, determined the assembly processes of the bacterial community, resulting in opposite trends along the soil layers in tilled and untilled soils. These findings expand the current understanding of the biogeography of soil bacterial communities across different land use types.

Integrated Biogeography and Assembly Mechanisms of Microeukaryotic Communities in Coastal Waters Near Shellfish Cultivation

The Lianjiang coast in the East China Sea is a typical subtropical marine ecosystem, and shellfish cultivation occupies almost all of the available tidal flats. Many studies have investigated the effects of shellfish cultivation on benthic organisms and sediments, while the impact of shellfish cultivation on plankton ecosystems is still poorly understood. This study investigated the biogeographical patterns of microeukaryotic communities from Lianjiang coastal waters in four seasons using 18S ribosomal RNA gene amplicon sequencing. Microeukaryotes were mainly comprised of Dinoflagellata, Diatomea, Arthropoda, Ciliophora, Chlorophyta, Protalveolata, Cryptophyceae, and Ochrophyta, and presented significant differences in three habitats (the aquaculture area, confluent area, and offshore area) and four seasons. Similarity percentage analysis revealed that Paracalanus parvus, Heterocapsa rotundata, Bestiolina similis, and five additional key taxa contributed to spatio-temporal differences. Seasonal environmental and spatial factors explained 27.47% of microeukaryotic community variation on average, with 11.11% of the variation shared. Environmental variables, particularly depth, pH, and nitrite concentration, were strongly associated with the microeukaryotic community compositions. The neutral community model further demonstrated that stochastic processes were sufficient in shaping substantial variation in microeukaryotic communities across four seasons, which may reveal the remaining unexplained microeukaryotic community variation. We further divided four seasons into the aquaculture stages and non-aquaculture stages, and speculated that aquaculture activities may increase the dispersal limitation of microeukaryotes in coastal waters, especially for the big bodied-microbes like Arthropoda. The results provide a better understanding of the biogeographical patterns, processes, and mechanisms of microeukaryotic communities near shellfish cultivation.

Dispersal Limitation Controlling the Assembly of the Fungal Community in Karst Caves

As a unique ecosystem, Karst caves harbor an impressive diversity of specific fungi. However, the factors and mechanisms that shape fungal biodiversity in caves remain elusive. In this study, we explored the assembly patterns of fungal communities based on our previous research in eight representative Karst caves in Southwest China. Our results indicated that dispersal limitation plays a crucial role in shaping the overall fungal community as well as specific communities in rock, sediment, and water samples. However, "Undominated" processes contributed more than dispersal limitation in air samples. Interestingly, the dominant assembly processes varied between caves. Consistently, environmental selection had a minor impact on the assembly of fungal communities. Among the examined spatial and environmental variables, latitude, longitude, altitude, and temperature were found to significantly influence fungal communities irrespective of substrate type. These findings provide valuable insights into the ecological factors governing fungal community assembly in Karst caves.

Distinct patterns of distribution, community assembly and cross-domain co-occurrence of planktonic archaea in four major estuaries of China

BACKGROUND

Archaea are key mediators of estuarine biogeochemical cycles, but comprehensive studies comparing archaeal communities among multiple estuaries with unified experimental protocols during the same sampling periods are scarce. Here, we investigated the distribution, community assembly, and cross-domain microbial co-occurrence of archaea in surface waters across four major estuaries (Yellow River, Yangtze River, Qiantang River, and Pearl River) of China cross climatic zones (~ 1,800 km) during the winter and summer cruises.

RESULTS

The relative abundance of archaea in the prokaryotic community and archaeal community composition varied with estuaries, seasons, and stations (reflecting local environmental changes such as salinity). Archaeal communities in four estuaries were overall predominated by ammonia-oxidizing archaea (AOA) (aka. Marine Group (MG) I; primarily Nitrosopumilus), while the genus Poseidonia of Poseidoniales (aka. MGII) was occasionally predominant in Pearl River estuary. The cross-estuary dispersal of archaea was largely limited and the assembly mechanism of archaea varied with estuaries in the winter cruise, while selection governed archaeal assembly in all estuaries in the summer cruise. Although the majority of archaea taxa in microbial networks were peripherals and/or connectors, extensive and distinct cross-domain associations of archaea with bacteria were found across the estuaries, with AOA as the most crucial archaeal group. Furthermore, the expanded associations of MGII taxa with heterotrophic bacteria were observed, speculatively indicating the endogenous demand for co-processing high amount and diversity of organic matters in the estuarine ecosystem highly impacted by terrestrial/anthropogenic input, which is worthy of further study.

CONCLUSIONS

Our results highlight the lack of common patterns in the dynamics of estuarine archaeal communities along the geographic gradient, expanding the understanding of roles of archaea in microbial networks of this highly dynamic ecosystem.

Unraveling antibiotic resistomes associated with bacterial and viral communities in intertidal mudflat aquaculture area

The extensive use of antibiotics in intertidal mudflat aquaculture area has substantially increased the dissemination risk of antibiotic resistance genes (ARGs). As hosts of ARGs, bacteria and virus exert vital effects on ARG dissemination. However, the insights for the interrelationships among ARGs, bacteria, and virus have not been thoroughly explored in intertidal mudflat. Therefore, this study attempts to unravel the occurrence, dissemination, evolution, and driving mechanisms of ARGs associated with bacterial and viral communities using metagenomic sequencing in a typical intertidal mudflat. Abundant and diverse ARGs (22 types and 437 subtypes) were identified and those of ARGs were higher in spring than in autumn. It is worthy noted that virus occupied a more essential position than bacteria for ARGs dissemination through network analysis. Meanwhile, nitrogen exerted indirect effect on ARG profiles by shaping viral and bacterial diversity. According to the results of neutral and null models, deterministic processes dominated the ARG community assembly by controlling sediment nitrogen and antibiotics. Homogeneous and variable selection dominated phylogenetic turnover of ARG community, contributing 46.15% and 45.90% of the total processes, respectively. This study can hence theoretically support for the ARG pollution control and management in intertidal mudflat aquaculture area.

A Metagenomic and Amplicon Sequencing Combined Approach Reveals the Best Primers to Study Marine Aerobic Anoxygenic Phototrophs

Studies based on protein-coding genes are essential to describe the diversity within bacterial functional groups. In the case of aerobic anoxygenic phototrophic (AAP) bacteria, the pufM gene has been established as the genetic marker for this particular functional group, although available primers are known to have amplification biases. We review here the existing primers for pufM gene amplification, design new ones, and evaluate their phylogenetic coverage. We then use samples from contrasting marine environments to evaluate their performance. By comparing the taxonomic composition of communities retrieved with metagenomics and with different amplicon approaches, we show that the commonly used PCR primers are biased towards the Gammaproteobacteria phylum and some Alphaproteobacteria clades. The metagenomic approach, as well as the use of other combinations of the existing and newly designed primers, show that these groups are in fact less abundant than previously observed, and that a great proportion of pufM sequences are affiliated to uncultured representatives, particularly in the open ocean. Altogether, the framework developed here becomes a better alternative for future studies based on the pufM gene and, additionally, serves as a reference for primer evaluation of other functional genes.

Epilithic Bacterial Assemblages on Subtidal Rocky Reefs: Variation Among Alternative Habitats at Ambient and Enhanced Nutrient Levels

Temperate rocky reefs often support mosaics of alternative habitats such as macroalgal forests, algal turfs and sea urchin barrens. Although the composition of epilithic microbial biofilms (EMBs) is recognized as a major determinant of macroalgal recruitment, their role in regulating the stability of alternative habitats on temperate rocky reefs remains unexplored. On shallow rocky reefs of the Island of Capraia (NW Mediterranean), we compared EMB structure among canopy stands formed by the fucoid Ericaria brachycarpa, algal turfs, and urchin barrens under ambient versus experimentally enhanced nutrient levels. The three habitats shared a core microbial community consisting of 21.6 and 25.3% of total ASVs under ambient and enhanced nutrient conditions, respectively. Although Gammaproteobacteria, Alphaproteobacteria and Flavobacteriia were the most abundant classes across habitats, multivariate analyses at the ASV level showed marked differences in EMB composition among habitats. Enhancing nutrient level had no significant effect on EMBs, although it increased their similarity between macroalgal canopy and turf habitats. At both ambient and enriched nutrient levels, ASVs mostly belonging to Proteobacteria and Bacteroidetes were more abundant in EMBs from macroalgal canopies than barrens. In contrast, ASVs belonging to the phylum of Proteobacteria and, in particular, to the families of Rhodobacteraceae and Flavobacteriaceae at ambient nutrient levels and of Rhodobacteraceae and Bacteriovoracaceae at enhanced nutrient levels were more abundant in turf than canopy habitats. Our results show that primary surfaces from alternative habitats that form mosaics on shallow rocky reefs in oligotrophic areas host distinct microbial communities that are, to some extent, resistant to moderate nutrient enhancement. Understanding the role of EMBs in generating reinforcing feedback under different nutrient loading regimes appears crucial to advance our understanding of the mechanisms underpinning the stability of habitats alternative to macroalgal forests as well as their role in regulating reverse shifts.

Water table level controls methanogenic and methanotrophic communities and methane emissions in a Sphagnum-dominated peatland

Peatlands are important sources of the greenhouse gas methane emissions equipoised by methanogens and methanotrophs. However, knowledge about how microbial functional groups associated with methane production and oxidation respond to water table fluctuations has been limited to date. Here, methane-related microbial communities and the potentials of methane production and oxidation were determined along sectioned peat layers in a subalpine peatland across four Sphagnum-dominated sites with different water table levels. Methane fluxes were also monitored at these sites. The results showed that mcrA gene copies for methanogens were the highest in the 10- to 15-cm peat layer, which was also characterized by the maximum potential methane production (24.53 ± 1.83 nmol/g/h). Copy numbers of the pmoA gene for type Ia and Ib methanotrophs were enriched in the 0-5 cm peat layer with the highest potential methane oxidation (43.09 ± 3.44 nmol/g/h). For the type II methanotrophs, the pmoA gene copies were higher in the 10- to 15-cm peat layer. Hydrogenotrophic methanogens and type II methanotrophs dominated the methane functional groups. Deterministic process contributed more to methanogenic and methanotrophic community assemblages in comparison with stochastic process. The level of water table significantly shaped methanogenic and methanotrophic community structures and regulated methane fluxes. Compared with vascular plants, Sphagnum mosses significantly reduced the methane emissions in peatlands. Collectively, these findings enhance a comprehensive understanding of the effect of the water table level on methane functional groups, with consequential implications for reducing methane emissions within peatland ecosystems.IMPORTANCEThe water table level is recognized as a critical factor in regulating methane emissions, which are largely dependent on the balance of methanogens and methanotrophs. Previous studies on peat methane emissions have been mostly focused on spatial-temporal variations and the relationship with meteorological conditions. However, the role of the water table level in methane emissions remains unknown. In this work, four representative microhabitats along a water table gradient in a Sphagnum-dominated peatland were sampled to gain an insight into methane functional communities and methane emissions as affected by the water table level. The changes in methane-related microbial community structure and assembly were used to characterize the response to the water table level. This study improves the understanding of the changes in methane-related microbial communities and methane emissions with water table levels in peatlands.

Unraveling the roles of coastal bacterial consortia in degradation of various lignocellulosic substrates

Lignocellulose, as the most abundant natural organic carbon on earth, plays a key role in regulating the global carbon cycle, but there have been only few studies in marine ecosystems. Little information is available about the extant lignin-degrading bacteria in coastal wetlands, limiting our understanding of their ecological roles and traits in lignocellulose degradation. We utilized in situ lignocellulose enrichment experiments coupled with 16S rRNA amplicon and shotgun metagenomics sequencing to identify and characterize bacterial consortia attributed to different lignin/lignocellulosic substrates in the southern-east intertidal zone of East China Sea. We found the consortia enriched on woody lignocellulose showed higher diversity than those on herbaceous substrate. This also revealed substrate-dependent taxonomic groups. A time-dissimilarity pattern with increased alpha diversity over time was observed. Additionally, this study identified a comprehensive set of genes associated with lignin degradation potential, containing 23 gene families involved in lignin depolymerization, and 371 gene families involved in aerobic/anaerobic lignin-derived aromatic compound pathways, challenging the traditional view of lignin recalcitrance within marine ecosystems. In contrast to similar cellulase genes among the lignocellulose substrates, significantly different ligninolytic gene groups were observed between consortia under woody and herbaceous substrates. Importantly, we not only observed synergistic degradation of lignin and hemi-/cellulose, but also pinpointed the potential biological actors at the levels of taxa and functional genes, which indicated that the alternation of aerobic and anaerobic catabolism could facilitate lignocellulose degradation. Our study advances the understanding of coastal bacterial community assembly and metabolic potential for lignocellulose substrates. IMPORTANCE It is essential for the global carbon cycle that microorganisms drive lignocellulose transformation, due to its high abundance. Previous studies were primarily constrained to terrestrial ecosystems, with limited information about the role of microbes in marine ecosystems. Through in situ lignocellulose enrichment experiment coupled with high-throughput sequencing, this study demonstrated different impacts that substrates and exposure times had on long-term bacterial community assembly and pinpointed comprehensive, yet versatile, potential decomposers at the levels of taxa and functional genes in response to different lignocellulose substrates. Moreover, the links between ligninolytic functional traits and taxonomic groups of substrate-specific populations were revealed. It showed that the synergistic effect of lignin and hemi-/cellulose degradation could enhance lignocellulose degradation under alternation of aerobic and anaerobic conditions. This study provides valuable taxonomic and genomic insights into coastal bacterial consortia for lignocellulose degradation.

Inferring microbial co-occurrence networks from amplicon data: a systematic evaluation

Microbes commonly organize into communities consisting of hundreds of species involved in complex interactions with each other. 16S ribosomal RNA (16S rRNA) amplicon profiling provides snapshots that reveal the phylogenies and abundance profiles of these microbial communities. These snapshots, when collected from multiple samples, can reveal the co-occurrence of microbes, providing a glimpse into the network of associations in these communities. However, the inference of networks from 16S data involves numerous steps, each requiring specific tools and parameter choices. Moreover, the extent to which these steps affect the final network is still unclear. In this study, we perform a meticulous analysis of each step of a pipeline that can convert 16S sequencing data into a network of microbial associations. Through this process, we map how different choices of algorithms and parameters affect the co-occurrence network and identify the steps that contribute substantially to the variance. We further determine the tools and parameters that generate robust co-occurrence networks and develop consensus network algorithms based on benchmarks with mock and synthetic data sets. The Microbial Co-occurrence Network Explorer, or MiCoNE (available at https://github.com/segrelab/MiCoNE) follows these default tools and parameters and can help explore the outcome of these combinations of choices on the inferred networks. We envisage that this pipeline could be used for integrating multiple data sets and generating comparative analyses and consensus networks that can guide our understanding of microbial community assembly in different biomes. IMPORTANCE Mapping the interrelationships between different species in a microbial community is important for understanding and controlling their structure and function. The surge in the high-throughput sequencing of microbial communities has led to the creation of thousands of data sets containing information about microbial abundances. These abundances can be transformed into co-occurrence networks, providing a glimpse into the associations within microbiomes. However, processing these data sets to obtain co-occurrence information relies on several complex steps, each of which involves numerous choices of tools and corresponding parameters. These multiple options pose questions about the robustness and uniqueness of the inferred networks. In this study, we address this workflow and provide a systematic analysis of how these choices of tools affect the final network and guidelines on appropriate tool selection for a particular data set. We also develop a consensus network algorithm that helps generate more robust co-occurrence networks based on benchmark synthetic data sets.

MArVD2: a machine learning enhanced tool to discriminate between archaeal and bacterial viruses in viral datasets

Our knowledge of viral sequence space has exploded with advancing sequencing technologies and large-scale sampling and analytical efforts. Though archaea are important and abundant prokaryotes in many systems, our knowledge of archaeal viruses outside of extreme environments is limited. This largely stems from the lack of a robust, high-throughput, and systematic way to distinguish between bacterial and archaeal viruses in datasets of curated viruses. Here we upgrade our prior text-based tool (MArVD) via training and testing a random forest machine learning algorithm against a newly curated dataset of archaeal viruses. After optimization, MArVD2 presented a significant improvement over its predecessor in terms of scalability, usability, and flexibility, and will allow user-defined custom training datasets as archaeal virus discovery progresses. Benchmarking showed that a model trained with viral sequences from the hypersaline, marine, and hot spring environments correctly classified 85% of the archaeal viruses with a false detection rate below 2% using a random forest prediction threshold of 80% in a separate benchmarking dataset from the same habitats.

Ocean-wide comparisons of mesopelagic planktonic community structures

For decades, marine plankton have been investigated for their capacity to modulate biogeochemical cycles and provide fishery resources. Between the sunlit (epipelagic) layer and the deep dark waters, lies a vast and heterogeneous part of the ocean: the mesopelagic zone. How plankton composition is shaped by environment has been well-explored in the epipelagic but much less in the mesopelagic ocean. Here, we conducted comparative analyses of trans-kingdom community assemblages thriving in the mesopelagic oxygen minimum zone (OMZ), mesopelagic oxic, and their epipelagic counterparts. We identified nine distinct types of intermediate water masses that correlate with variation in mesopelagic community composition. Furthermore, oxygen, NO3- and particle flux together appeared as the main drivers governing these communities. Novel taxonomic signatures emerged from OMZ while a global co-occurrence network analysis showed that about 70% of the abundance of mesopelagic plankton groups is organized into three community modules. One module gathers prokaryotes, pico-eukaryotes and Nucleo-Cytoplasmic Large DNA Viruses (NCLDV) from oxic regions, and the two other modules are enriched in OMZ prokaryotes and OMZ pico-eukaryotes, respectively. We hypothesize that OMZ conditions led to a diversification of ecological niches, and thus communities, due to selective pressure from limited resources. Our study further clarifies the interplay between environmental factors in the mesopelagic oxic and OMZ, and the compositional features of communities.

Untangling microbiota diversity and assembly patterns in the world's longest underground culvert water diversion canal

The long-distance underground box culvert water transport system (LUBWT) is a crucial link between the source of drinking water and the consumers. It must ensure the stability of water quality during transportation. However, uncontrollable microbial growth can develop in the water delivery system during the long delivery process, posing a risk to health and safety. Therefore, we applied 16 s and 18 s gene sequence analysis in order to study microbial communities in box culvert waters sampled in 2021, as well as a molecular ecological network-based approach to decipher microbial interactions and stability. Our findings revealed that, in contrast to natural freshwater ecosystems, micro-eukaryotes in LUBWT have complex interactions such as predation, parasitism, and symbiosis due to their semi-enclosed box culvert environment. Total nitrogen may be the primary factor affecting bacterial community interactions in addition to temperature. Moreover, employing stability indicators such as robustness and vulnerability, we also found that microbial stability varied significantly from season to season, with summer having the higher stability of microbial communities. Not only that but also the stability of the micronuclei also varied greatly during water transport, which might also be related to the complex interactions among the micro-eukaryotes. To summarize, our study reveals the microbial interactions and stability in LUBWT, providing essential ecological knowledge to ensure the safety of LUBWT's water quality.

Modules in robust but low-efficiency phyllosphere fungal networks drive saponin accumulation in leaves of different Panax species

BACKGROUND

The phyllosphere mycobiome plays a crucial role in plant fitness and ecosystem functions. The complex microbial ecological networks (MEN) formed by these fungi remain poorly understood, particularly with regard to their organization strategy and their contributions to plant secondary metabolites such as saponin.

RESULTS

In this study, we constructed six MENs from leaf epiphytic and endophytic mycobiomes of three Panax species distributed in the northeast and southwest ends of mainland China. Hub nodes were absent in these MENs, which were significantly more complex, robust, and less efficient compared to random graphs (P < 0.05), indicating a hub-independent high-robustness strategy to maintain structural homeostasis. The important roles of specific MEN modules in shaping leaf saponin profiles of each Panax species were proved by multiple machine learning algorithms. Positive regulation modules (PRMs) of total saponin content were further identified, which exhibited more deterministic ecological assembly and comprised of highly connected nodes as well as higher proportion of plant-associated fungal guilds compared to other network members, indicating their tight links with host plant. The significant and direct effects (P < 0.05) of PRMs on total saponin accumulation were validated by well-fitted structural equation models (χ² < 0.3, P > 0.5). Taxonomic analysis revealed that Pleosporales and Chaetothyriales were significantly overrepresented by positive regulation taxa (PRT) of total saponin content (FDR < 0.05). Across PRT identified in three Panax species, Epicoccum and Coniothyrium were conservatively present, while species-specific taxa such as Agaricales were also found, indicating the conservatism and specificity of plant-fungi interactions associated with leaf saponin accumulation in Panax genus.

CONCLUSIONS

These findings provide a foundation for understanding mechanisms maintaining the steady state of phyllosphere mycobiome in healthy plant, and offer clues for engineering phyllosphere mycobiome to improve the accumulation of bioactive secondary metabolites on the basis of network modules.

Temperature dependence of parasitoid infection and abundance of a diatom revealed by automated imaging and classification

Diatoms are a group of phytoplankton that contribute disproportionately to global primary production. Traditional paradigms that suggest diatoms are consumed primarily by larger zooplankton are challenged by sporadic parasitic "epidemics" within diatom populations. However, our understanding of diatom parasitism is limited by difficulties in quantifying these interactions. Here, we observe the dynamics of Cryothecomonas aestivalis (a protist) infection of an important diatom on the Northeast U.S. Shelf (NES), Guinardia delicatula, with a combination of automated imaging-in-flow cytometry and a convolutional neural network image classifier. Application of the classifier to >1 billion images from a nearshore time series and >20 survey cruises across the broader NES reveals the spatiotemporal gradients and temperature dependence of G. delicatula abundance and infection dynamics. Suppression of parasitoid infection at temperatures <4 °C drives annual cycles in both G. delicatula infection and abundance, with an annual maximum in infection observed in the fall-winter preceding an annual maximum in host abundance in the winter-spring. This annual cycle likely varies spatially across the NES in response to variable annual cycles in water temperature. We show that infection remains suppressed for ~2 mo following cold periods, possibly due to temperature-induced local extinctions of the C. aestivalis strain(s) that infect G. delicatula. These findings have implications for predicting impacts of a warming NES surface ocean on G. delicatula abundance and infection dynamics and demonstrate the potential of automated plankton imaging and classification to quantify phytoplankton parasitism in nature across unprecedented spatiotemporal scales.

Deterministic processes dominate archaeal community assembly from the Pearl River to the northern South China Sea

Archaea play a significant role in the biogeochemical cycling of nutrients in estuaries. However, comprehensive researches about their assembly processes remain notably insufficient. In this study, we systematically examined archaeal community dynamics distinguished between low-salinity and high-salinity groups in water and surface sediments over a 600-kilometer range from the upper Pearl River (PR) to the northern South China Sea (NSCS). Neutral community model analysis together with null model analysis showed that their C-score values were greater than 2, suggesting that deterministic processes could dominate the assembly of those planktonic or benthic archaeal communities at both the low-salinity and high-salinity sites. And deterministic processes contributed more in the low-salinity than high-salinity environments from the PR to the NSCS. Furthermore, through the co-occurrence network analysis, we found that the archaeal communities in the low-salinity groups possessed closer interactions and higher proportions of negative interactions than those in the high-salinity groups, which might be due to the larger environmental heterogeneities reflected by the nutrient concentrations of those low-salinity samples. Collectively, our work systematically investigated the composition and co-occurrence networks of archaeal communities in water as well as sediments from the PR to the NSCS, yielding new insights into the estuary's archaeal community assembly mechanisms.

Deciphering Interactions Within a 4-Strain Riverine Bacterial Community

The dynamics of a community of four planktonic bacterial strains isolated from river water was followed in R2 broth for 72 h in batch experiments. These strains were identified as Janthinobacterium sp., Brevundimonas sp., Flavobacterium sp. and Variovorax sp. 16S rRNA gene sequencing and flow cytometry analyses were combined to monitor the change in abundance of each individual strain in bi-cultures and quadri-culture. Two interaction networks were constructed that summarize the impact of the strains on each other's growth rate in exponential phase and carrying capacity in stationary phase. The networks agree on the absence of positive interactions but also show differences, implying that ecological interactions can be specific to particular growth phases. Janthinobacterium sp. was the fastest growing strain and dominated the co-cultures. However, its growth rate was negatively affected by the presence of other strains 10 to 100 times less abundant than Janthinobacterium sp. In general, we saw a positive correlation between growth rate and carrying capacity in this system. In addition, growth rate in monoculture was predictive of carrying capacity in co-culture. Taken together, our results highlight the necessity to take growth phases into account when measuring interactions within a microbial community. In addition, evidence that a minor strain can greatly influence the dynamics of a dominant one underlines the necessity to choose population models that do not assume a linear dependency of interaction strength to abundance of other species for accurate parameterization from such empirical data.

Diversity and putative interactions of parasitic alveolates belonging to Syndiniales at a coastal Pacific site

The dinoflagellate lineage Syndiniales currently consists entirely of parasitic species that fall into five well-supported clades. Environmental sequencing studies worldwide have found an abundance of Syndiniales in a variety of marine ecosystems, but very little is known about the majority of Syndiniales species including two entire clades which have only been observed in sequence data. Syndiniales are known to have a wide range of hosts, but only a few dozen interactions have been confirmed through observation of actual infections. Here, we describe the diversity of Syndiniales found at the Scripps Institution of Oceanography pier over the course of a year based on 18S sequencing. We find Syndiniales to be the most species (amplicon sequence variant)-rich taxonomic group and for its members to be present and abundant throughout the year. We used several analytical techniques to identify potential parasite-host interactions which we were then able to visualize over time. Using mock communities and size fractionation of seawater, we suggest that the majority of Syndiniales sequences that are found in environmental studies belong to the free-living dinospore stage rather than representing active infections.

Abundance does not predict extinction risk in the fossil record of marine plankton

A major premise of ecological neutral theory is that population size is inversely related to extinction risk. This idea is central to modern biodiversity conservation efforts, which often rely on abundance metrics to partially determine species extinction risk. However, limited empirical studies have tested whether extinction is indeed more probable for species with low abundances. Here we use the fossil record of Neogene radiolaria to test the relationship between relative abundance and longevity (time from first to last occurrence). Our dataset includes abundance histories for 189 polycystine radiolarian species from the Southern Ocean, and 101 species from the tropical Pacific. Using linear regression analyses, we show that neither maximum nor average relative abundance are significant predictors of longevity in either oceanographic region. This suggests that neutral theory fails to explain the plankton ecological-evolutionary dynamics we observe. Extrinsic factors are likely more important than neutral dynamics in controlling radiolarian extinction.

The underappreciated role of biotic factors in controlling the bloom ecology of potentially harmful microalgae in the Hudson-Raritan Bay

Despite widespread distribution of harmful algal blooms (HABs) and new and improved methods for detecting and quantifying them, no unifying ecological explanation has been found. Improved understanding depends upon local, ecological studies that include analysis of phytoplankton species data in relation to both abiotic and biotic parameters. Ecological network analysis was used to detect co-occurrence patterns among abiotic and biotic parameters in a long-term monitoring dataset (i.e., 2010-2021) from the eutrophic Hudson-Raritan Estuary (HRE) between the states of New York and New Jersey. The regular co-occurrence of potentially harmful bloom-forming species with companion species observed through microscopy was supported by the results of ecological network analysis, which showed that there were far more associations between HAB species and biotic parameters (∼95%) than abiotic parameters (∼5%). Temperature was the environmental variable that was most associated with HAB species throughout the estuary. The numerous network associations of HAB species with one another and with diatoms, dinoflagellates, and zooplankton highlight the complexity of planktonic food webs and interactions. Results also suggest that some taxa may play a central role in structuring the HRE plankton communities. These findings demonstrate that biotic associations play an underappreciated role in plankton structure and the value of examining the ecology of HAB species within the breadth of their biological communities. While network analysis does not fully explain and confirm complex associations among species, it does provide fresh insights and testable hypotheses to strengthen understanding and improve prediction.

Stronger Geographic Limitations Shape a Rapid Turnover and Potentially Highly Connected Network of Core Bacteria on Microplastics

Core microbiota is shared microbial taxa within the same habitat, which is important for understanding the stable and consistent components of the complex microbial assembly. However, information on the microplastic core bacteria from the river ecosystems is poor. Here, we investigated the composition and function of microplastic core bacteria from the Three Gorges Reservoir area along the approximate 662 km of the Yangtze River via full-length 16S rRNA gene sequencing, compared with those in water, sediment, and soil. The results showed that the spatial turnover of bacterial communities in four habitats supported deterministic processes dominated by niche differentiation, which shaped their core bacteria. The composition and function of microplastic core bacteria were significantly different from those in the other three habitats. Rhodobacteraceae was the main component of microplastic core bacteria, while the main component of core bacteria in water, sediment, and soil were Burkholderiaceae (21.90%), Burkholderiaceae (5.01%), Nitrosomonadaceae (4.61%), respectively. Furthermore, microplastic core bacteria had stronger geographic limitations along the Yangtze River in the Three Gorges Reservoir area. Stronger geographic limitations shaped the rapid community turnover and a potentially more connected network for the microplastic core bacteria than water, sediment, and soil. More importantly, microplastic core bacteria had strong potential functions of drug resistance and could cause risks to ecosystems and human health. Microplastic core bacteria were mainly influenced by sediment core bacteria, although the bacteria colonizing on microplastics could be from all the contact environments and original sources. These findings provide important insights into the composition, function, and association of microplastic core bacteria with their surrounding environment.

Resource-dependent biodiversity and potential multi-trophic interactions determine belowground functional trait stability

BACKGROUND

For achieving long-term sustainability of intensive agricultural practices, it is pivotal to understand belowground functional stability as belowground organisms play essential roles in soil biogeochemical cycling. It is commonly believed that resource availability is critical for controlling the soil biodiversity and belowground organism interactions that ultimately lead to the stabilization or collapse of terrestrial ecosystem functions, but evidence to support this belief is still limited. Here, we leveraged field experiments from the Chinese National Ecosystem Research Network (CERN) and two microcosm experiments mimicking high and low resource conditions to explore how resource availability mediates soil biodiversity and potential multi-trophic interactions to control functional trait stability.

RESULTS

We found that agricultural practice-induced higher resource availability increased potential cross-trophic interactions over 316% in fields, which in turn had a greater effect on functional trait stability, while low resource availability made the stability more dependent on the potential within trophic interactions and soil biodiversity. This large-scale pattern was confirmed by fine-scale microcosm systems, showing that microcosms with sufficient nutrient supply increase the proportion of potential cross-trophic interactions, which were positively associated with functional stability. Resource-driven belowground biodiversity and multi-trophic interactions ultimately feedback to the stability of plant biomass.

CONCLUSIONS

Our results indicated the importance of potential multi-trophic interactions in supporting belowground functional trait stability, especially when nutrients are sufficient, and also suggested the ecological benefits of fertilization programs in modern agricultural intensification. Video Abstract.

Significant Effects of Associated Microorganisms on the Community of Photosynthetic Picoeukaryotes

Photosynthetic picoeukaryotes (PPEs) form associations with other microorganisms. However, whether and how the associated microbes affect PPE communities remain unknown. We used flow cytometric cell sorting combined with parallel high-throughput sequencing of the 18S and 16S rRNA genes to simultaneously investigate PPEs and their associated microbial communities in the Yangtze-connected Lake Dongting. The lake harbors a great diversity of PPEs. PPE communities exhibited significant temporal rather than spatial variations. Two distinct PPE taxa affiliated with Discostella nipponica and Poterioochromonas malhamensis were dominant during winter/spring and summer, respectively. Parallel high-throughput sequencing revealed a great diversity of associated bacteria and non-pigmented eukaryotes (NPEs) in PPEs sorts. Proteobacteria, Actinobacteria, Bacteroidetes, and Cyanobacteria among the associated bacteria and fungi among the associated NPEs were dominant. PPEs were more apparently associated with bacteria than with NPEs. The co-occurrence network of PPEs and associated microbes formed five major modules, which exhibited distinct temporal patterns, being specific to a certain period. Variations in PPEs communities were significantly correlated with both environmental factors and associated microbial communities. In variation partitioning analysis, the associated bacteria explained the greatest variations in PPE communities, and associated bacteria and NPEs co-explained a large portion of environmental effects on PPE communities. Our results highlight the significance of associated microbes in shaping PPE communities.