Neutral community model explains the bacterial community assembly in freshwater lakes

Effects of reduced flow gradient on benthic biofilm communities' ecological network and community assembly

The intensification of human activities has led to flow reduction and cut-off in most global rivers, seriously affecting riverine organisms and the biogeochemical processes. As key indicators of river ecosystems' structure and function, benthic biofilms play a critical role in driving primary production and material cycling in rivers. This research aimed to investigate the characteristics of microbial communities' complexity and stability during river flow reduction. Benthic biofilms were grown in artificial channels and subjected to eight gradients of flow reduction (represented by flow velocity from 0.4 to 110 cm/s). Biofilms' biodiversity, ecological networks and community assembly of bacteria, fungi and algae were investigated by high-throughput sequencing. Results showed significant differences in community composition and structure under different flow conditions. The eight flow gradients' microbial communities were divided into three groups: low, medium and high flows. The flow reduction led to significant decreases in bacterial and fungal communities' Chao1 index. Low flow conditions enriched the bacterial phyla Oxyphotobacteria, Alphaproteobacteria and Mollicutes, but significantly decreased the fungal phylum Chytridiomycota. Lowering flow reduced the fungal network's number of nodes and increased the algal network's number of edges. Cross-domain interactions network analysis showed a gradual increase in node and edge numbers with decreasing flow, while decreasing average path length. The neutral model predicted stochastic processes primarily drove biofilm community assembly, and that model's explanations decreased as the flow gradient decreased. The null model analysis revealed diffusion limitation as the most common stochastic ecological process for bacterial and algal communities, with reduced flow reducing heterogeneous selection and increasing diffusion-limited processes. This study provides an in-depth analysis of flow reduction's effects on biofilm communities' ecological networks and community assembly.

Analysis of bacterial community structure, functional variation, and assembly mechanisms in multi-media habitats of lakes during the frozen period

Ice, water, and sediment represent three interconnected habitats in lake ecosystems, and bacteria are crucial for maintaining ecosystem equilibrium and elemental cycling across these habitats. However, the differential characteristics and driving mechanisms of bacterial community structures in the ice, water, and sediments of seasonally frozen lakes remain unclear. In this study, high-throughput sequencing technology was used to analyze and compare the structure, function, network characteristics, and assembly mechanisms of bacterial communities in the ice, water, and sediment of Wuliangsuhai, a typical cold region in Inner Mongolia. The results showed that the bacterial communities in the ice and water phases had similar diversity and composition, with Proteobacteria, Bacteroidota, Actinobacteria, Campilobacterota, and Cyanobacteria as dominant phyla. The bacterial communities in sediments displayed significant differences from ice and water, with Chloroflexi, Proteobacteria, Firmicutes, Desulfobacterota, and Acidobacteriota being the dominant phyla. Notably, the bacterial communities in water exhibited higher spatial variability in their distribution than those in ice and sediment. This study also revealed that during the frozen period, the bacterial community species in the ice, water, and sediment media were dominated by cooperative relationships. Community assembly was primarily influenced by stochastic processes, with dispersal limitation and drift identified as the two most significant factors within this process. However, heterogeneous selection also played a significant role in the community composition. Furthermore, functions related to nitrogen, phosphorus, sulfur, carbon, and hydrogen cycling vary among bacterial communities in ice, water, and sediment. These findings elucidate the intrinsic mechanisms driving variability in bacterial community structure and changes in water quality across different media phases (ice, water, and sediment) in cold-zone lakes during the freezing period, offering new insights for water environmental protection and ecological restoration efforts in such environments.

Effect of Altitude Gradients on the Spatial Distribution Mechanism of Soil Bacteria in Temperate Deciduous Broad-Leaved Forests

Soil bacteria are an important part of the forest ecosystem, and they play a crucial role in driving energy flow and material circulation. Currently, many uncertainties remain about how the composition and distribution patterns of bacterial communities change along altitude gradients, especially in forest ecosystems with strong altitude gradients in climate, vegetation, and soil properties. Based on dynamic site monitoring of the Baiyun Mountain Forest National Park (33°38'-33°42' N, 111°47'-111°51' E), this study used Illumina technology to sequence 120 soil samples at the site and explored the spatial distribution mechanisms and ecological processes of soil bacteria under different altitude gradients. Our results showed that the composition of soil bacterial communities varied significantly between different altitude gradients, affecting soil bacterial community building by influencing the balance between deterministic and stochastic processes; in addition, bacterial communities exhibited broader ecological niche widths and a greater degree of stochasticity under low-altitude conditions, implying that, at lower altitudes, community assembly is predominantly influenced by stochastic processes. Light was the dominant environmental factor that influenced variation in the entire bacterial community as well as other taxa across different altitude gradients. Moreover, changes in the altitude gradient could cause significant differences in the diversity and community composition of bacterial taxa. Our study revealed significant differences in bacterial community composition in the soil under different altitude gradients. The bacterial communities at low elevation gradients were mainly controlled by stochasticity processes, and bacterial community assembly was strongly influenced by deterministic processes at middle altitudes. Furthermore, light was an important environmental factor that affects differences. This study revealed that the change of altitude gradient had an important effect on the development of the soil bacterial community and provided a theoretical basis for the sustainable development and management of soil bacteria.

DNA metabarcoding reveals ecological patterns and driving mechanisms of archaeal, bacterial, and eukaryotic communities in sediments of the Sansha Yongle Blue Hole

The Sansha Yongle Blue Hole (SYBH) is the world's deepest marine blue hole with unique physicochemical characteristics. However, our knowledge of the biodiversity and community structure in SYBH sediments remains limited, as past studies have mostly focused on microbial communities in the water column. Here, we collected sediment samples from the aerobic zone (3.1 to 38.6 m) and the deep anaerobic zone (150 m, 300 m) of the SYBH and extracted DNA to characterize the archaeal, bacterial, and eukaryotic communities inhabiting these sediments. Our results showed that the archaeal and bacterial communities were dominated by Thaumarchaeota and Proteobacteria, respectively. The dominant taxa of eukaryotes in different sites varied greatly, mainly including Phaeophyceae, Annelida, Diatomea and Arthropoda. All three examined domains showed clear vertical distributions and significant differences in community composition between the aerobic and anaerobic zones. Sulfide played a prominent role in structuring the three domains, followed by salinity, nitrous oxide, pH, temperature and dissolved oxygen, all of which were positively correlated with the turnover component, the main contributor to beta diversity. Neutral community model revealed that stochastic processes contributed to more than half of the community variations across the three domains. Co-occurrence network showed an equal number of positive and negative interactions in the archaeal network, while positive interactions accounted for ~ 80% in the bacterial and eukaryotic networks. Our findings reveal the ecological features of prokaryotes and eukaryotes in SYBH sediments and shed new light on community dynamics and survival strategies in the special environment of marine blue holes.

Phyllosphere bacterial community dynamics in response to bacterial wildfire disease: succession and interaction patterns

Plants interact with complex microbial communities in which microorganisms play different roles in plant development and health. While certain microorganisms may cause disease, others promote nutrient uptake and resistance to stresses through a variety of mechanisms. Developing plant protection measures requires a deeper comprehension of the factors that influence multitrophic interactions and the organization of phyllospheric communities. High-throughput sequencing was used in this work to investigate the effects of climate variables and bacterial wildfire disease on the bacterial community's composition and assembly in the phyllosphere of tobacco (Nicotiana tabacum L.). The samples from June (M1), July (M2), August (M3), and September (M4) formed statistically separate clusters. The assembly of the whole bacterial population was mostly influenced by stochastic processes. PICRUSt2 predictions revealed genes enriched in the M3, a period when the plant wildfire disease index reached climax, were associated with the development of the wildfire disease (secretion of virulence factor), the enhanced metabolic capacity and environmental adaption. The M3 and M4 microbial communities have more intricate molecular ecological networks (MENs), bursting with interconnections within a densely networked bacterial population. The relative abundances of plant-beneficial and antagonistic microbes Clostridiales, Bacillales, Lactobacillales, and Sphingobacteriales, showed significant decrease in severally diseased sample (M3) compared to the pre-diseased samples (M1/M2). Following the results of MENs, we further test if the correlating bacterial pairs within the MEN have the possibility to share functional genes and we have unraveled 139 entries of such horizontal gene transfer (HGT) events, highlighting the significance of HGT in shaping the adaptive traits of plant-associated bacteria across the MENs, particularly in relation to host colonization and pathogenicity.

Temperature and nutrients alter the relative importance of stochastic and deterministic processes in the coastal macroinvertebrates biodiversity assembly on long-time scales

Macroinvertebrates play a vital role in coastal ecosystems and are an important indicator of ecosystem quality. Both anthropogenic activity and environmental changes may lead to significant changes in the marine macroinvertebrate community. However, the assembly process of benthic biodiversity and its mechanism driven by environmental factors at large scales remains unclear. Here, using the benthic field survey data of 15 years at large spatial and temporal scales from the Yellow Sea Large Marine Ecosystem, we investigated the relative importance of environmental selection, dispersal processes, random-deterministic processes of macroinvertebrates community diversity assembly, and the responses of this relative importance driven by temperature and nutrients. Results showed that the macroinvertebrates community diversity is mainly affected by dispersal. Nitrogen and phosphorus are the most important negative factors among environmental variables, while geographical distance is the main limiting factor of β diversity. Within the range of 0.35-0.70 mg/L of nutrients, increasing nutrient concentration can significantly facilitate the contribution of the decay effect to β diversity. Within the temperature range studied (15.0-18.0°C), both warming and cooling can lead to a greater tendency for species diversity assembly processes to be dominated by deterministic processes. The analysis contributes to a better understanding of the assembly process of the diversity of coastal marine macroinvertebrates communities and how they adapt to global biogeochemical processes.

Effects of lactic acid bacteria on rearing water bacterial community in Eriocheir sinensis culture

The Chinese mitten crab (CMC, Eriocheir sinensi) culture in ponds is a unique aquaculture system. Probiotics are commonly used in the maintenance of the health of pond-cultured CMCs. However, the effects of probiotics on the bacterial community of CMC-culturing water remain unclear. This study utilized 16S rRNA gene amplicon sequencing to assess changes in the bacterial community composition, diversity, assembly, and co-occurrence patterns in CMC-culturing water following probiotic application. The results indicate that the α-diversity of the bacterial community in CMC-culturing water varied with time following probiotic application. The addition of probiotics to the water resulted in an increase in the occurrence of new operational taxonomic units (OTUs). The bacterial community assembly in the CMC-culturing water was shaped by a balance between deterministic and stochastic processes, while commercial probiotics enhanced the proportion of heterogeneous selection. In addition, including OTU2953 (Burkholderiaceae) and OTU3005 (Lactobacillaceae), from the commercial probiotics served as keystone species in the bacterial network of CMC-culturing water. Overall, probiotic application had a significant impact on the bacterial ecology of CMC-culturing water.

Effects of Cascade Reservoirs on Spatiotemporal Dynamics of the Sedimentary Bacterial Community: Co-occurrence Patterns, Assembly Mechanisms, and Potential Functions

Dam construction as an important anthropogenic activity significantly influences ecological processes in altered freshwater bodies. However, the effects of multiple cascade dams on microbial communities have been largely overlooked. In this study, the spatiotemporal distribution, co-occurrence relationships, assembly mechanisms, and functional profiles of sedimentary bacterial communities were systematically investigated in 12 cascade reservoirs across two typical karst basins in southwest China over four seasons. A significant spatiotemporal heterogeneity was observed in bacterial abundance and diversity. Co-occurrence patterns in the Wujiang Basin exhibited greater edge counts, graph density, average degree, robustness, and reduced modularity, suggesting more intimate and stronger ecological interactions among species than in the Pearl River Basin. Furthermore, Armatimonadota and Desulfobacterota, identified as keystone species, occupied a more prominent niche than the dominant species. A notable distance-decay relationship between geographical distance and community dissimilarities was identified in the Pearl River Basin. Importantly, in the Wujiang Basin, water temperature emerged as the primary seasonal variable steering the deterministic process of bacterial communities, whereas 58.5% of the explained community variance in the neutral community model (NCM) indicated that stochastic processes governed community assembly in the Pearl River Basin. Additionally, principal component analysis (PCA) revealed more pronounced seasonal dynamics in nitrogen functional compositions than spatial variation in the Wujiang Basin. Redundancy analysis (RDA) results indicated that in the Wujiang Basin, environmental factors and in Pearl River Basin, geographical distance, reservoir age, and hydraulic retention time (HRT), respectively, influenced the abundance of nitrogen-related genes. Notably, these findings offer novel insights: building multiple cascade reservoirs could lead to a cascading decrease in biodiversity and resilience in the river-reservoir ecosystem.

Microplastic polymer properties as deterministic factors driving terrestrial plastisphere microbiome assembly and succession in the field

Environmental microplastic (MP) is ubiquitous in aquatic and terrestrial ecosystems providing artificial habitats for microbes. Mechanisms of MP colonization, MP polymer impacts, and effects on soil microbiomes are largely unknown in terrestrial systems. Therefore, we experimentally tested the hypothesis that MP polymer type is an important deterministic factor affecting MP community assembly by incubating common MP polymer types in situ in landfill soil for 14 months. 16S rRNA gene amplicon sequencing indicated that MP polymers have specific impacts on plastisphere microbiomes, which are subsets of the soil microbiome. Chloroflexota, Gammaproteobacteria, certain Nitrososphaerota, and Nanoarchaeota explained differences among MP polymers and time points. Plastisphere microbial community composition derived from different MP diverged over time and was enriched in potential pathogens. PICRUSt predictions of pathway abundances and quantitative PCR of functional marker genes indicated that MP polymers exerted an ambivalent effect on genetic potentials of biogeochemical cycles. Overall, the data indicate that (i) polymer type as deterministic factor rather than stochastic factors drives plastisphere community assembly, (ii) MP impacts greenhouse gas metabolism, xenobiotic degradation and pathogen distribution, and (iii) MP serves as an ideal model system for studying fundamental questions in microbial ecology such as community assembly mechanisms in terrestrial environments.

Planktonic microbial community and biological metabolism in a subtropical drinking water river-reservoir system

To understand the dynamics of planktonic microbial community and its metabolism processes in subtropical drinking water river-reservoir system with lower man-made pollution loading, this study selected Dongzhen river-reservoir system in Mulan Creek as object to investigate spatial-temporal characteristics of community profile and functional genes involved in biological metabolism, and to analyze the influence of environmental factors. The results indicated that Proteobacteria and Actinobacteria were the most diverse phyla with proportion ranges of 9%-80% in target system, and carbohydrate metabolism (5.76-7.12 × 10-2), amino acid metabolism (5.78-7.21 × 10-2) and energy metabolism (4.07-5.17 × 10-2) were found to be the dominant pathways of biological metabolism. Although there were variations in biological properties both spatially and temporally, seasonal variation had a greater influence on microbial community and biological metabolism, than locational differences. Regarding the role of environmental factors, this study revealed that microbial diversity could be affected by multiple abiotic factors, with total organic carbon, total phosphorus and temperature being more influential (absolute value of standardized regression weights >2.13). Stochastic processes dominated the microbial community assembly (R2 of neutral community model = 0.645), while niche-based processes differences represented by nutrients, temperature and pH level played secondary roles (R > 0.388, P < 0.01). Notably, the synergistic influences among the environmental factors accounted for the higher percentages of community variation (maximum proportion up to 17.6%). Additionally, pH level, temperature, and concentrations of dissolved oxygen, carbon and nitrogen were found to be the significant factors affecting carbon metabolism pathways (P < 0.05), yet only total organic carbon significantly affected on nitrogen transformation (P < 0.05). In summary, the microbial profile in reservoir is not completely dominated by that in feeding river, and planktonic microbial community and its metabolism in subtropical drinking water river-reservoir system are shaped by multiple abiotic and biotic factors with underlying interactions.

Composition and assembly mechanisms of prokaryotic communities in wetlands, and their relationships with different vegetation and reclamation methods

Coastal wetlands are undergoing substantial transformations globally as a result of increased human activities. However, compared to other ecosystems, diversity and functional characteristics of microbial communities in reclaimed coastal wetlands are not well studied compared to other ecosystems. This is important because it is known that microorganisms can play a crucial role in biogeochemical cycling within coastal wetland ecosystems. Hence, this study utilized the high-throughput sequencing technique to investigate the structure and assembly processes of microbial communities in reclaimed coastal wetlands. The results revealed a substantial change in soil properties following coastal wetland reclamation. Remarkably, the reclaimed soil exhibited significantly lower pH, soil organic carbon (SOC), and total salinity (TS) values (p < 0.05). The dominant phyla included Proteobacteria, Chloroflexi, Bacteroidetes, Acidobacteria, and Planctomycetes among study sites. However, the relative abundance of Proteobacteria increased from un-reclaimed coastal wetlands to reclaimed ones. The Proteobacteria, Chloroflexi, and Acidobacteria showed higher relative abundance in vegetated soil compared to bare soil, while Bacteroidetes and Planctomycetes exhibited the opposite trend. Notably, vegetation types exerted the strongest influence on microbial diversity, surpassing the effects of soil types and depth (F = 34.49, p < 0.001; F = 25.49, p < 0.001; F = 3.173, p < 0.078, respectively). Stochastic assembly processes dominated in un-reclaimed soil, whereas deterministic processes governed the assembly in artificial sea embankment wetlands (SEW). The presence of Spartina alterniflora in all soil types (except SEW soils) indicated stochastic assembly, while Phragmites australis in reclaimed soils pointed toward deterministic microbial assembly. Furthermore, environmental factors such as pH, soil water content (SWC), SOC, total carbon (TC), total nitrogen (TN), total phosphorus (TP), NH4+-N, vegetation types, soil depth, and geographic distance exhibited significant effects on microbial beta diversity indices. Co-occurrence network analysis revealed a stronger association between taxa in SEW compared to land reclaimed from wetlands (LRW) and natural coastal wetlands (NCW). The bottom soil layer exhibited more complex network interactions than the topsoil layer. Besides soil parameters, reclamation and varieties of vegetation were also substantial factors influencing the composition, diversity, and assembly processes of microbial communities in coastal wetlands.

Contrasting Mechanisms Determine the Microeukaryotic and Syndiniales Community Assembly in a Eutrophic bay

Syndiniales is a diverse parasitic group, increasingly gaining attention owing to its high taxonomic diversity in marine ecosystems and inhibitory effects on the dinoflagellate blooms. However, their seasonal dynamics, host interactions, and mechanisms of community assembly are largely unknown, particularly in eutrophic waters. Here, using 18S rRNA gene amplicon sequencing, we intended to elucidate the interactions between Syndiniales and microeukaryotes, as well as community assembly processes in a eutrophic bay. The results showed that Syndiniales group II was dominating throughout the year, with substantially higher abundance in the winter and spring, whereas Syndiniales group I was more abundant in the summer and autumn. Temperature and Dinoflagellata were the most important abiotic and biotic factors driving variations of the Syndiniales community, respectively. The assembly processes of microeukaryotes and Syndiniales were completely different, with the former being controlled by a balance between homogeneous selection and drift and the latter being solely governed by drift. Network analysis revealed that Syndiniales group II had the largest number of interactions with microeukaryotes, and they primarily associated with Dinoflagellata in the winter, while interactions with Chlorophyta and Bacillariophyta increased dramatically in summer and autumn. These findings provide significant insights in understanding the interactions and assembly processes of Syndiniales throughout the year, which is critical in revealing the roles of single-celled parasites in driving protist dynamics in eutrophic waters.

Genome-resolved metagenomics provides insights into the ecological roles of the keystone taxa in heavy-metal-contaminated soils

Microorganisms that exhibit resistance to environmental stressors, particularly heavy metals, have the potential to be used in bioremediation strategies. This study aimed to explore and identify microorganisms that are resistant to heavy metals in soil environments as potential candidates for bioremediation. Metagenomic analysis was conducted using microbiome metagenomes obtained from the rhizosphere of soil contaminated with heavy metals and mineral-affected soil. The analysis resulted in the recovery of a total of 175 metagenome-assembled genomes (MAGs), 73 of which were potentially representing novel taxonomic levels beyond the genus level. The constructed ecological network revealed the presence of keystone taxa, including Rhizobiaceae, Xanthobacteraceae, Burkholderiaceae, and Actinomycetia. Among the recovered MAGs, 50 were associated with these keystone taxa. Notably, these MAGs displayed an abundance of genes conferring resistance to heavy metals and other abiotic stresses, particularly those affiliated with the keystone taxa. These genes were found to combat excessive accumulation of zinc/manganese, arsenate/arsenite, chromate, nickel/cobalt, copper, and tellurite. Furthermore, the keystone taxa were found to utilize both organic and inorganic energy sources, such as sulfur, arsenic, and carbon dioxide. Additionally, these keystone taxa exhibited the ability to promote vegetation development in re-vegetated mining areas through phosphorus solubilization and metabolite secretion. In summary, our study highlights the metabolic adaptability and ecological significance of microbial keystone taxa in mineral-affected soils. The MAGs associated with keystone taxa exhibited a markedly higher number of genes related to abiotic stress resistance and plant growth promotion compared to non-keystone taxa MAGs.

Expanding known viral diversity in the healthy infant gut

The gut microbiome is shaped through infancy and impacts the maturation of the immune system, thus protecting against chronic disease later in life. Phages, or viruses that infect bacteria, modulate bacterial growth by lysis and lysogeny, with the latter being especially prominent in the infant gut. Viral metagenomes (viromes) are difficult to analyse because they span uncharted viral diversity, lacking marker genes and standardized detection methods. Here we systematically resolved the viral diversity in faecal viromes from 647 1-year-olds belonging to Copenhagen Prospective Studies on Asthma in Childhood 2010, an unselected Danish cohort of healthy mother-child pairs. By assembly and curation we uncovered 10,000 viral species from 248 virus family-level clades (VFCs). Most (232 VFCs) were previously unknown, belonging to the Caudoviricetes viral class. Hosts were determined for 79% of phage using clustered regularly interspaced short palindromic repeat spacers within bacterial metagenomes from the same children. Typical Bacteroides-infecting crAssphages were outnumbered by undescribed phage families infecting Clostridiales and Bifidobacterium. Phage lifestyles were conserved at the viral family level, with 33 virulent and 118 temperate phage families. Virulent phages were more abundant, while temperate ones were more prevalent and diverse. Together, the viral families found in this study expand existing phage taxonomy and provide a resource aiding future infant gut virome research.

Multi-group biodiversity distributions and drivers of metacommunity organization along a glacial-fluvial-limnic pathway on the Tibetan plateau

Extensive global glacial retreats are threatening cryosphere ecosystem functioning and the associated biota in glacier-fed water systems. Understanding multi-group biodiversity distributions and compositional variation across diverse but hydrologically linked habitats under varying glacial influences will help explain the mechanisms underlying glacial community organization and ecosystem processes. However, such data are generally lacking due to the difficulty of obtaining biodiversity information across wide taxonomic ranges. Here, we used a multi-marker environmental DNA metabarcoding approach to simultaneously investigate the spatial patterns of community compositions and assembly mechanisms of four taxonomic groups (cyanobacteria, diatoms, invertebrates, and vertebrates) along the flowpaths of a tributary of Lake Nam Co on the Tibetan Plateau-from its glacier headwaters, through its downstream river and wetlands, to its estuary. We detected 869 operational taxonomic units: 119 cyanobacterial, 395 diatom, 269 invertebrate, and 86 vertebrate. Taxonomic richnesses consistently increased from upstream to downstream, and although all groups showed community similarity distance decay patterns, the trend for vertebrates was the weakest. Cyanobacteria, diatom, and invertebrate community compositions were significantly correlated with several environmental factors, while the vertebrate community was only correlated with waterway width. Variation partitioning analysis indicated that varying extents of environmental conditions and spatial factors affected community organizations for different groups. Furthermore, stochastic processes contributed prominently to the microorganisms' community assembly (Sloan's neutral model R² = 0.77 for cyanobacteria and 0.73 for diatoms) but were less important for macroorganisms (R² = 0.21 for invertebrates and 0.15 for vertebrates). That trend was further substantiated by modified stochasticity ratio analyses. This study provides the first holistic picture of the diverse biotic communities residing in a series of hydrologically connected glacier-influenced habitats. Our results both uncovered the distinct mechanisms that underlie the metacommunity organizations of different glacial organisms and helped comprehensively predict the ecological impacts of the world's melting glaciers.

Responses of the bacterial community of tobacco phyllosphere to summer climate and wildfire disease

Both biotic and abiotic factors continually affect the phyllospheric ecology of plants. A better understanding of the drivers of phyllospheric community structure and multitrophic interactions is vital for developing plant protection strategies. In this study, 16S rRNA high-throughput sequencing was applied to study how summer climatic factors and bacterial wildfire disease have affected the composition and assembly of the bacterial community of tobacco (Nicotiana tabacum L.) phyllosphere. Our results indicated that three time series groups (T1, T2 and T3) formed significantly distinct clusters. The neutral community model (NCM) and beta nearest taxon index (betaNTI) demonstrated that the overall bacterial community assembly was predominantly driven by stochastic processes. Variance partitioning analysis (VPA) further showed that the complete set of the morbidity and climatic variables together could explain 35.7% of the variation of bacterial communities. The node numbers of the molecular ecological networks (MENs) showed an overall uptrend from T1 to T3. Besides, Pseudomonas is the keystone taxa in the MENs from T1 to T3. PICRUSt2 predictions revealed significantly more abundant genes of osmoprotectant biosynthesis/transport in T2, and more genes for pathogenicity and metabolizing organic substrate in T3. Together, this study provides insights into spatiotemporal patterns, processes and response mechanisms underlying the phyllospheric bacterial community.

Salinity and seasonality shaping free-living and particle-associated bacterioplankton community assembly in lakeshores of the northeastern Qinghai-Tibet Plateau

Microorganisms in lakeshore zones are essential for pollution interception and biodiversity maintenance. However, the biogeographic patterns of bacterioplankton communities in lakeshore zones and the mechanisms that driving them are poorly understood. We analyzed the 16 S rRNA gene sequences of particle-associated (PA) and free-living (FL) bacterioplankton communities in the lakeshore zones of 14 alpine lakes in two seasons on Qinghai-Tibet Plateau to investigate the bacterial diversity, composition and assembly processes. Our results revealed that PA and FL bacterioplankton communities were driven by both seasonality and salinity in the lakeshores on Qinghai-Tibet Plateau. Compared to FL bacterioplankton, PA bacterioplankton communities were more susceptible to seasonality than spatial salinity. FL bacterioplankton communities were more salinity constrained than the PA counterpart. Besides, the Stegen null model analyses have validated a quantitative bias on stochastic processes at different spatial scales. At a regional scale, stochasticity was the predominant assembly process in both PA and FL bacterioplankton. While at a subregional scale, dispersal limitation was the main contributor of stochastic processes for PA bacterioplankton in summer and heterogeneous selection was the dominant deterministic processes in winter, whereas the community assembly of FL bacterioplankton was more stochastic processes (i.e., dispersal limitation) dominated in the freshwater type but deterministic process (i.e., heterogeneous selection) increased with increasing salinity. Our study provides new insights into both significant spatiotemporal patterns and distinct assembly processes of PA and FL bacterioplankton in alpine lakeshores on the northeastern Qinghai-Tibet Plateau.

Geographic pattern of phytoplankton community and their drivers in lakes of middle and lower reaches of Yangtze River floodplain, China

Disentangling the relative contributions of deterministic and stochastic processes was critical to compressive understanding of underlying mechanism governing geographic pattern and assembly of phytoplankton community, while it was seldom performed in connected lakes under human pressure. Here, we investigated phytoplankton community pattern in relation to environmental and spatial factors over 81 lakes located in the middle and lower reaches of Yangtze River (MLYR) floodplain, where many lakes suffered from eutrophication and cyanobacterial blooms. A majority of MLYR lakes had higher phytoplankton abundance surpassing 107 cells/L and were dominated by common bloom-forming cyanobacterial genera, including Pseudanabaena, Microcystis, Merismopedia, Dolichospermum, Limnothrix, and Raphidiopsis. Phytoplankton community exhibited a striking geographical pattern both for taxonomic and functional compositions, while functional groups were less sensitive, and dissimilarity in communities displayed no significant increases with increasing geographical distance. Further, species richness explained much higher percentage of community variations than species turnover, indicating a reduced effect of environmental filtering of phytoplankton species with tolerance to similar environments in connected MLYR lakes. Both deterministic and stochastic processes governed assembly and biogeographic of phytoplankton community. Variation partition analysis showed that spatial factors exhibited greater influence on phytoplankton community compared to environmental variables. The stronger influence of spatial factors was further demonstrated by Mantel test and neutral community model. These findings indicate that deterministic and stochastic processes exhibited similar biogeographic patterns for phytoplankton community in MLYR lakes, but stochastic process was overwhelmingly dominated. Moreover, a large proportion of unexplained variation implies that complex interactions exist to shape assembly mechanism of phytoplankton community in MLYR lakes.

Community assembly of microbial habitat generalists and specialists in urban aquatic ecosystems explained more by habitat type than pollution gradient

Urban freshwater ecosystems have important ecosystem functions, provide habitats for diverse microbial communities and are susceptible to multiple interconnected factors such as environmental pollution. Despite the ecological significance of bacteria and microeukaryotes, little is known about how their community assembly responds to various environmental factors across water and sediment habitats and ecological processes shaping them. Here, environmental DNA-based approaches were used to investigate the community assembly processes of bacteria and microeukaryotes (including habitat generalists and specialists) in urban water and sediment across an urban-pollution gradient in Wuhan, central China. The diversity, community composition and potential function of bacteria and microeukaryotes showed significantly stronger variation between water and sediment than across an urban pollution gradient. Although, bacterial and microeukaryotic community assemblies were dominated by strong selection processes in both water and sediment habitats, but a contrasting community assembly mechanism was identified between habitat generalists and specialists. Bacterial and microeukaryotic communities showed a greater response to physicochemical variability in water, while a strong distance-decay relationship was found in sediment. Further, cross-kingdom microbial network analysis revealed strong modular associations of bacteria and microeukaryotes, meanwhile, microeukaryotic habitat specialists might be keystone, but generalists have higher proportion of connections in the networks. This study provides significant insights into the response of bacteria and microeukaryotes to different urban pollutions between water and sediment, and the ecological processes structuring microbial community dynamics across habitat types under anthropogenic disturbances.

Robust host source tracking building on the divergent and non-stochastic assembly of gut microbiomes in wild and farmed large yellow croaker

BACKGROUND

Given the lack of genetic background, the source tracking unknown individuals of fish species with both farmed and wild populations often cannot be robustly achieved. The gut microbiome, which is shaped by both deterministic and stochastic processes, can serve as a molecular marker of fish host source tracking, particularly as an alternative to the yet-to-be-established host genetic marker. A candidate for testing the feasibility is the large yellow croaker, Larimichthys crocea, which is carnivorous and ranks the top mariculture fish in China. Wild resource of this fish was depleted decades ago and might have potential problematic estimation because of escaping of farmed individuals.

RESULTS

The rectums of wild (n = 212) and farmed (n = 79) croakers from multiple batches were collected for the profiling of their gut bacterial communities. The farmed individuals had a higher alpha diversity and lower bacterial load than the wild individuals. The gut microbiota of the two sources exhibited divergence and high inter-batch variation, as featured by the dominance of Psychrobacter spp. in the wild group. Predicted functional capacity of the gut microbiome and representative isolates showed differences in terms of host source. This difference can be linked to the potential diet divergence between farmed and wild fishes. The non-stochastic distribution pattern of the core gut microbiota of the wild and farmed individuals supports the feasibility of microbiota-based host source tracking via the machine learning algorithm. A random forest classifier based on the divergence and non-stochastic assembly of the gut microbiome was robust in terms of host source tracking the individuals from all batches of croaker, including a newly introduced batch.

CONCLUSIONS

Our study revealed the divergence of gut microbiota and related functional profiles between wild and farmed croakers. For the first time, with representative datasets and non-stochastic patterns, we have verified that gut microbiota can be robustly applied to the tracking of host source even in carnivorous fish. Video abstract.

Distribution and Assembly Processes of Soil Fungal Communities along an Altitudinal Gradient in Tibetan Plateau

In soil ecosystems, fungi exhibit diverse biodiversity and play an essential role in soil biogeochemical cycling. Fungal diversity and assembly processes across soil strata along altitudinal gradients are still unclear. In this study, we investigated the structure and abundance of soil fungal communities among soil strata and elevational gradients on the Tibetan Plateau using Illumina MiSeq sequencing of internal transcribed spacer1 (ITS1). The contribution of neutral and niche ecological processes were quantified using a neutral community model and a null model-based methodology. Our results showed that fungal gene abundance increased along altitudinal gradients, while decreasing across soil strata. Along with altitudinal gradients, fungal α-diversity (richness) decreased from surface to deeper soil layers, while β-diversity showed weak correlations with elevations. The neutral community model showed an excellent fit for neutral processes and the lowest migration rate (R² = 0.75). The null model showed that stochastic processes dominate in all samples (95.55%), dispersal limitations were dominated at the surface layer and decreased significantly with soil strata, while undominated processes (ecological drift) show a contrary trend. The log-normal model and the null model (βNTI) correlation analysis also neglect the role of niche-based processes. We conclude that stochastic dispersal limitations, together with ecological drifts, drive fungal communities.

Diversity Distribution, Driving Factors and Assembly Mechanisms of Free-Living and Particle-Associated Bacterial Communities at a Subtropical Marginal Sea

Free-living (FL) and particle-associated (PA) bacterioplankton communities play critical roles in biogeochemical cycles in the ocean. However, their community composition, assembly process and functions in the continental shelf and slope regions are poorly understood. Based on 16S rRNA gene amplicon sequencing, we investigated bacterial communities’ driving factors, assembly processes and functional potentials at a subtropical marginal sea. The bacterioplankton community showed specific distribution patterns with respect to lifestyle (free living vs. particle associated), habitat (slope vs. shelf) and depth (surface vs. DCM and Bottom). Salinity and water temperature were the key factors modulating turnover in the FL community, whereas nitrite, silicate and phosphate were the key factors for the PA community. Model analyses revealed that stochastic processes outweighed deterministic processes and had stronger influences on PA than FL. Homogeneous selection (Hos) was more responsible for the assembly and turnover of FL, while drift and dispersal limitation contributed more to the assembly of PA. Importantly, the primary contributor to Hos in PA was Gammaproteobacteria:Others, whereas that in FL was Cyanobacteria:Bin6. Finally, the PICRUSt2 analysis indicated that the potential metabolisms of carbohydrates, cofactors, amino acids, terpenoids, polyketides, lipids and antibiotic resistance were markedly enriched in PA than FL.

Integrating Stochastic and Deterministic Process in the Biogeography of N2-Fixing Cyanobacterium Candidatus Atelocyanobacterium Thalassa

UCYN-A is one of the most widespread and important marine diazotrophs. Its unusual distribution in both cold/warm and coastal/oceanic waters challenges current understanding about what drives the biogeography of diazotrophs. This study assessed the community assembly processes of the nitrogen-fixing cyanobacterium UCYN-A, developing a framework of assembly processes underpinning the microbial biogeography and diversity. High-throughput sequencing and a qPCR approach targeting the nifH gene were used to investigate three tropical seas: the Bay of Bengal, the Western Pacific Ocean, and the South China Sea. Based on the neutral community model and two types of null models calculating the β-nearest taxon index and the normalized stochasticity ratio, we found that stochastic assembly processes could explain 66-92% of the community assembly; thus, they exert overwhelming influence on UCYN-A biogeography and diversity. Among the deterministic processes, temperature and coastal/oceanic position appeared to be the principal environmental factors driving UCYN-A diversity. In addition, a close linkage between assembly processes and UCYN-A abundance/diversity/drivers can provide clues for the unusual global distribution of UCYN-A.

Assembly Processes and Co-occurrence Patterns of Abundant and Rare Bacterial Community in the Eastern Indian Ocean

Microbial communities are composed of many rare species and a few abundant species. Considering the disproportionate importance of rare species for ecosystem functioning, it is important to understand the mechanisms structuring the rare and abundant components of a diverse community in response to environmental changes. Here, we used a 16S ribosomal RNA gene sequencing approach to investigate the bacterial community diversity in the Eastern Indian Ocean (EIO) during the monsoon and intermonsoon. We employed a phylogenetic null model and network analysis to evaluate the assembly processes and co-occurrence pattern of the microbial community. We found that higher bacterial diversity was detected in the intermonsoon with high temperature and low Chlorophyll a concentrations and N/P ratios. The balance between ecological deterministic processes and stochastic processes varied with seasons in the EIO. Meanwhile, conditionally rare taxa (CRT) were more likely modulated by variable selection processes than always rare taxa (ART) and abundant taxa (AT) (CRT > ART > AT). By linking assembly process and species co-occurrence, we demonstrated that the microbial co-occurrence associations tended to be higher when deterministic processes (mainly variable selection) were weaker. This negative trend was observed in rare species rather than abundant species. The linkage could enhance our understanding of the underlying mechanisms underpinning the generation and maintenance of microbial community diversity.

Changes in protist communities in drainages across the Pearl River Delta under anthropogenic influence

Drainages in the Pearl River Delta urban agglomeration (PRDUA) host vital aquatic ecosystems and face enormous pressures from human activities in one of the largest urban agglomerations in the world. Despite being crucial components of aquatic ecosystems, the interactions and assembly processes of the protistan community are rarely explored in areas with serious anthropogenic disturbance. To elucidate the mechanisms of these processes, we used environmental DNA sequencing of 18S rDNA to investigate the influence of environmental factors and species interactions on the protistan community and its assembly in drainages of the PRDUA during summer. The protistan community showed a high level of diversity and a marked spatial pattern in this region. Community assembly was driven primarily by stochastic processes based on the Sloan neutral community model, explaining 74.28%, 75.82%, 73.67%, 74.40% and 51.24% of community variations in the BJ (Beijiang), XJ (Xijiang), PRD (Pearl River Delta), PRE (Pearl River Estuary) areas and in total, respectively. Meanwhile, environmental variables including temperature, pH, dissolved oxygen, transparency, nutrients and land use were strongly correlated with the composition and assembly of the protistan community, explaining 40.40% of variation in the protistan community. Furthermore, the bacterial community was simultaneously analysed by the 16S rDNA sequencing. Co-occurrence network analysis revealed that species interactions within bacteria (81.41% positive) or protists (82.80% positive), and those between bacteria and protists (50% positive and 50% negative) impacted the protistan community assembly. In summary, stochastic processes dominated, whereas species interactions and environmental factors also played important roles in shaping the protistan communities in drainages across the PRDUA. This study provides insights into the ecological patterns, assembly processes and species interactions underlying protistan dynamics in urban aquatic ecosystems experiencing serious anthropogenic disturbance.

Transient-rare Bacterial Taxa Are Assembled Neutrally across Temporal Scales

Despite the importance of microbial communities in ecosystem functions, the mechanisms underlying the assembly of rare taxa over time are poorly understood. It remains largely unknown whether rare taxa exhibit similar assembly processes to common taxa in local communities. We herein retrieved the 16S rRNA sequences of bacteria collected bimonthly for 2 years from the Pohang wastewater treatment plant. The transient-rare taxa showed different abundance distributions from the common taxa. Transient-rare taxon assemblages also exhibited higher temporal variations than common taxon assemblages, suggesting the distinct ecological patterns of the two assemblages. A multivariate analysis revealed that environmental parameters accounted for 25.3 and 61.6% of temporal variations in the transient-rare and common taxon assemblages, respectively. The fitting of all observed taxa to a neutral community model revealed that 96.4% of the transient-rare taxa (relative abundance, 71.4%) and 73.3% of the common taxa (relative abundance, 45.6%) followed the model, suggesting that stochastic mechanisms were more important than deterministic ones in the assembly of the transient-rare taxa. Collectively, the present results indicate that the transient-rare bacterial taxa at the Pohang wastewater treatment plant differed from the common taxa in ecological patterns, suggesting that dispersal is a key process in their assembly.

Bacterial Community Spacing Is Mainly Shaped by Unique Species in the Subalpine Natural Lakes of China

Bacterial communities have been described as early indicators of both regional and global climatic change and play a critical role in the global biogeochemical cycle. Exploring the mechanisms that determine the diversity patterns of bacterial communities and how they share different habitats along environmental gradients are, therefore, a central theme in microbial ecology research. We characterized the diversity patterns of bacterial communities in Pipahai Lake (PPH), Mayinghai Lake (MYH), and Gonghai Lake (GH), three subalpine natural lakes in Ningwu County, Shanxi, China, and analyzed the distribution of their shared and unique taxa (indicator species). Results showed that the species composition and structure of bacterial communities were significantly different among the three lakes. Both the structure of the entire bacterial community and the unique taxa were significantly influenced by the carbon content (TOC and IC) and space distance; however, the structure of the shared taxa was affected by conductivity (EC), pH, and salinity. The structure of the entire bacterial community and unique taxa were mainly affected by the same factors, suggesting that unique taxa may be important in maintaining the spatial distribution diversity of bacterial communities in subalpine natural freshwater lakes. Our results provide new insights into the diversity maintenance patterns of the bacterial communities in subalpine lakes, and suggest dispersal limitation on bacterial communities between adjacent lakes, even in a small local area. We revealed the importance of unique taxa in maintaining bacterial community structure, and our results are important in understanding how bacterial communities in subalpine lakes respond to environmental change in local habitats.

Responses of abundant and rare bacterioplankton to temporal change in a subtropical urban reservoir

Investigation of bacterial community dynamics across different time scales is important for understanding how environmental conditions drive community change over time. Bacterioplankton from the surface waters of a subtropical urban reservoir in southeast China were analyzed through high-frequency sampling over 13 months to compare patterns and ecological processes between short (0‒8 weeks), medium (9‒24 weeks) and long (25‒53 weeks) time intervals. We classified the bacterial community into different subcommunities: abundant taxa (AT); conditionally rare taxa (CRT); rare taxa (RT). CRT contributed > 65% of the alpha-diversity, and temporal change of beta-diversities was more pronounced for AT and CRT than RT. The bacterial community exhibited a directional change in the short- and medium-time intervals and a convergent dynamic during the long-time interval due to a seasonal cycle. Cyanobacteria exhibited a strong succession pattern than other phyla. CRT accounted for > 76% of the network nodes in three stations. The bacteria-environment relationship and deterministic processes were stronger for large sample size at station G (n = 116) than small sample size at stations C (n = 12) and L (n = 22). These findings suggest that a high-frequency sampling approach can provide a better understanding on the time scales at which bacterioplankton can change fast between being abundant or rare, thus providing the facts about environmental factors driving microbial community dynamics. Patterns and processes in alpha- and beta-diversities and community assembly of bacterioplankton differ among different time intervals (short-, medium- and long-time intervals) and different subcommunities (abundant, conditionally rare and rare taxa) in a subtropical urban reservoir, demonstrating the importance of temporal scale and high-frequency sampling in microbial community ecology.

Spatial and seasonal variations of sediment bacterial communities in a river-bay system in South China

River-bay systems are transitional areas that hold important roles in biogeochemical processes between continents and oceans. However, composition and structure of microbial communities shaped by such environments have not been clear yet. In this study, we used high-throughput sequencing of 16S rRNA genes to analyze the diversity and composition of sediment bacterial communities from the Shenzhen river-bay system during dry and wet seasons. The results showed that sediment bacterial community structure was varied according to habitats (river vs. estuary) and seasons (wet season vs. dry season). The alpha diversity of sediment bacterial community was significantly higher in the dry season than in the wet season, while no significant difference of alpha diversity was found between river and estuary. Neutral community model revealed a significant influence of stochastic processes on sediment bacterial community assembly, especially in the wet season. However, the beta nearest-taxon index indicated that deterministic processes were more responsible for the assembly of sediment bacterial community. Additionally, redundancy analysis suggested strong links between sediment bacterial communities and environmental factors in Shenzhen river-bay system, with the environmental factors explaining 63.5% of the bacterial community variation. Specifically, NH₄⁺, pH, and salinity were the three most important contributing factors that shaped the sediment bacterial communities. Overall, this study provides a valuable reference to get insights into the spatiotemporal pattern of sediment bacterial communities in a typical river-bay system. KEY POINTS: • Stochastic processes contribute sediment bacterial community assembly. • Deterministic processes dominate sediment bacterial community assembly. • Environmental factors shape sediment bacterial communities.

The balance between deterministic and stochastic processes in structuring lake bacterioplankton community over time

One major goal in microbial ecology is to establish the importance of deterministic and stochastic processes for community assembly. This is relevant to explain and predict how diversity changes at different temporal scales. However, understanding of the relative quantitative contribution of these processes and particularly of how they may change over time is limited. Here, we assessed the importance of deterministic and stochastic processes based on the analysis of the bacterial microbiome in one alpine oligotrophic and in one subalpine mesotrophic lake, which were sampled over two consecutive years at different time scales. We found that in both lakes, homogeneous selection (i.e., a deterministic process) was the main assembly process at the annual scale and explained 66.7% of the bacterial community turnover, despite differences in diversity and temporal variability patterns between ecosystems. However, in the alpine lake, homogenizing dispersal (i.e., a stochastic process) was the most important assembly process at the short‐term (daily and weekly) sampling scale and explained 55% of the community turnover. Alpha diversity differed between lakes, and seasonal stability of the bacterial community was more evident in the oligotrophic lake than in the mesotrophic one. Our results demonstrate how important forces that govern temporal changes in bacterial communities act at different time scales. Overall, our study validates on a quantitative basis, the importance and dominance of deterministic processes in structuring bacterial communities in freshwater environments over long time scales.

Species Sorting and Neutral Theory Analyses Reveal Archaeal and Bacterial Communities Are Assembled Differently in Hot Springs

Although the recognition of archaea as one of the three kingdoms in the tree of life has been nearly a half-century long, the comparative investigations on their ecological adaptations with bacteria have been limited. The mechanisms of their community assembly and diversity maintenance in hot springs have not been addressed. The mechanistic study is critical not only for understanding the hot-spring microbiome structure and dynamics, but also for shedding light on their evolutionary adaptations. We applied the neutral theory model and species sorting paradigm of metacommunity theory to investigate how hot-spring microbial communities were assembled, how their diversities were maintained, and how the temperature and pH influence these mechanisms. Through rigorous statistical tests based on the neutral theory and species sorting paradigm, we found (i) According to the neutral theory, archaeal and bacterial communities are assembled differently, with stochastic neutral force playing a more significant role in archaeal communities than in bacterial communities (neutrality-rate = 52.9 vs. 15.8%, p-value < 0.05). (ii) Temperature and pH account for rather limited (<10%) variations in hot-spring microbiomes based on the species sorting paradigm. The pH has more significant influences than temperature on archaeal communities, and both pH and temperature have similarly low influences on bacterial community structure. (iii) We postulate that the differences between archaea and bacteria are likely due to the longer evolutionary history and better adaptation of archaea to host spring environments.

Dispersal limitation drives biogeographical patterns of anammox bacterial communities across the Yangtze River

Over the past few decades, anaerobic ammonium oxidation (anammox) has been extensively documented at different scales in natural ecosystems. Previous studies have stated that the community composition of anammox bacteria is shaped mainly by environmental factors, whereas spatial factors have been largely overlooked. This study investigated biogeographical patterns of anammox bacterial communities using 42 sediment samples along a 4300-km stretch of the Yangtze River, the longest river in Asia. A significant distance-decay relationship was observed for anammox bacterial community similarity, which was significantly influenced by mean dendritic distance rather than environmental factors. This implied that dispersal limitation plays an important role in shaping biogeographical pattern of anammox community. Furthermore, our results revealed that neutral processes played vital role in shaping community assembly of anammox bacteria, and their communities were seriously dispersal limited. These findings contrast with previous studies on community similarities between broad taxonomic groups, which are mainly determined by niche-based selection owing to greater niche distances within broad taxonomic groups than in anammox bacteria. Importantly, the slope of the distance-decay curve was much steeper than previously reported for whole bacteria, which indicating the species turnover rate of anammox bacterial community was significantly higher than that of the whole bacterial community. Anammox bacteria harbor strong adhesion ability and low dispersal potential, and ultimately exhibited a high species turnover rate. Together, in the context of biogeography, our results highlight the importance of dispersal limitation in shaping the biogeographical pattern of anammox bacterial community.

Identification and characterization of core sludge and biofilm microbiota in anaerobic membrane bioreactors

An analysis of sludge (i.e., 63 samples) and biofilm (i.e., 79 samples) sampled from 13 anaerobic membrane bioreactors (AnMBR) was conducted. Predominant microbial community identification and multivariate analysis indicate that these reactors showed different microbial community structure, but these differences had no impact on the overall AnMBR performance. Instead, core microbial genera which occurred in ≥90% of sludge (20 genera) and biofilm (12 genera) samples could potentially account for the AnMBR performance. A further calculation on net growth rate (NGR) of core genera in sludge suggested distribution into two main groups (i.e., I: low relative abundance and NGR, II: high relative abundance or high NGR). Consistent positive correlations between bacterial genera were observed among those that exhibited either high relative abundance or high NGR. The anaerobic microbial consortium in both sludge and biofilm were largely affected by stochastic dispersal and migration processes (i.e., neutral assembly). However, Acinetobacter spp. and Methanobacterium spp. occurred consistently in higher frequency in the biofilm but in lower occurrence frequency in the AnMBR permeate. Findings from this study suggest first, specific core microorganisms exist in the sludge regardless of the operating conditions of the AnMBRs, and second, prevention of biofoulant layer on anaerobic membranes can be devised by minimizing attachment of microbes on surfaces in a non-selective manner.

Stochastic processes shape microeukaryotic community assembly in a subtropical river across wet and dry seasons

BACKGROUND

The deep mechanisms (deterministic and/or stochastic processes) underlying community assembly are a central challenge in microbial ecology. However, the relative importance of these processes in shaping riverine microeukaryotic biogeography is still poorly understood. Here, we compared the spatiotemporal and biogeographical patterns of microeukaryotic community using high-throughput sequencing of 18S rRNA gene and multivariate statistical analyses from a subtropical river during wet and dry seasons.

RESULTS

Our results provide the first description of biogeographical patterns of microeukaryotic communities in the Tingjiang River, the largest river in the west of Fujian province, southeastern China. The results showed that microeukaryotes from both wet and dry seasons exhibited contrasting community compositions, which might be owing to planktonic microeukaryotes having seasonal succession patterns. Further, all components of the microeukaryotic communities (including total, dominant, always rare, and conditionally rare taxa) exhibited a significant distance-decay pattern in both seasons, and these communities had a stronger distance-decay relationship during the dry season, especially for the conditionally rare taxa. Although several variables had a significant influence on the microeukaryotic communities, the environmental and spatial factors showed minor roles in shaping the communities. Importantly, these microeukaryotic communities were strongly driven by stochastic processes, with 89.9%, 88.5%, and 89.6% of the community variation explained by neutral community model during wet, dry, and both seasons, respectively. The neutral community model also explained a large fraction of the community variation across different taxonomic groups and levels. Additionally, the microeukaryotic taxa, which were above and below the neutral prediction, were ecologically and taxonomically distinct groups, which might be interactively structured by deterministic and stochastic processes.

CONCLUSIONS

This study demonstrated that stochastic processes are sufficient in shaping substantial variation in river microeukaryotic metacommunity across different hydrographic regimes, thereby providing a better understanding of spatiotemporal patterns, processes, and mechanisms of microeukaryotic community in waters.

A comprehensive census of lake microbial diversity on a global scale

Despite recent interest in microbial diversity and community structure of lakes across various spatial scales, a global biogeographic distribution pattern and its controlling factors have not been fully disclosed. Here, we compiled and analyzed 88,334,735 environmental 16S rRNA sequences from 431 lakes across a wide range of geographical distance and environmental conditions (in particular, salinity, 0-373.3 g L^-1). Our results showed that lake sediments inhabit significantly (ANOVA: P<0.001) more diverse microbial communities than lake waters. Non-metric dimensional scaling (NMDS) ordinations indicated that microbial community compositions differed distinctly among sample types (freshwater vs. saline, water vs. sediment) and geographic locations. Mantel and partial Mantel tests showed that microbial community composition in lake water was significantly (P=0.001) correlated with geographic distance, salinity, and pH. Statistical analyses based on neutral community and null models indicated that stochastic processes may play predominant roles in shaping the microbial biogeographic distribution patterns in the studied global lake waters. The dispersal-related stochasticity (e.g., homogenizing dispersal) exhibited a stronger influence on the distribution of microbial community in freshwater lakes than in saline lakes. Overall, this work expands our understanding of the impact of geographic distance, environmental conditions, and stochastic processes on microbial distribution in global lakes.

Biogeographic patterns of abundant and rare bacterioplankton in three subtropical bays resulting from selective and neutral processes

Unraveling the relative importance of ecological processes regulating microbial community structure is a central goal in microbial ecology. Here, we used high-throughput sequencing to examine the relative contribution of selective and neutral processes in the assembly of abundant and rare subcommunities from three subtropical bays of China. We found that abundant and rare bacterial taxa were distinctly different in diversity, despite the similar biogeographic patterns and strong distance-decay relationships, but the dispersal of rare bacterial taxa was more limited than that of abundant taxa. Furthermore, the environmental (selective processes) and spatial (neutral processes) factors seemed to govern the assembly and biogeography of abundant and rare bacterial subcommunities, although both factors explained only a small fraction of variation within the rare subcommunity. More importantly, variation partitioning (based on adjusted R2 in redundancy analysis) showed that spatial factors exhibited a slightly greater influence on both abundant and rare subcommunities compared to environmental selection; however, the abundant subcommunity had a much stronger response to spatial factors (17.3% of pure variance was explained) than that shown by the rare bacteria (3.5%). These results demonstrate that environmental selection and neutral processes explained the similar biogeographic patterns of abundant and rare subcommunities, but a large proportion of unexplained variation in the rare taxa (91.1%) implies that more complex assembly mechanisms may exist to shape the rare bacterial assemblages in the three subtropical bays.

The antibiotic resistome of free-living and particle-attached bacteria under a reservoir cyanobacterial bloom

In freshwater systems, both antibiotic resistance genes (ARGs) and cyanobacterial blooms attract global public health concern. Cyanobacterial blooms can greatly impact bacterial taxonomic communities, but very little is known about the influence of the blooms on antibiotic resistance functional community. In this study, the ARGs in both free-living (FL) and particle-attached (PA) bacteria under bloom and non-bloom conditions were simultaneously investigated in a subtropical reservoir using high-throughput approaches. In total, 145 ARGs and 9 mobile genetic elements (MGEs) were detected. The most diverse and dominant of which (68.93%) were multidrug resistance genes and efflux pump mechanism. The richness of ARGs in both FL and PA bacteria was significantly lower during the bloom period compared with non-bloom period. The abundance of ARGs in FL bacteria was significantly lower under bloom condition than in the non-bloom period, but the abundance of ARGs in PA bacteria stayed constant. More importantly, the resistant functional community in PA bacteria was more strongly influenced by the cyanobacterial bloom than in the FL bacteria, although >96% ARGs were shared in both FL and PA bacteria or both bloom and non-bloom periods. We also compared the community compositions between taxonomy and function, and found antibiotic resistant communities were highly variable and exhibited lower similarity between bloom and non-bloom periods than seen in the taxonomic composition, with an exception of FL bacteria. Altogether, cyanobacterial blooms appear to have stronger inhibitory effect on ARG abundance in FL bacteria, and stronger influence on antibiotic resistant community composition in PA bacteria. Our results further suggested that both neutral and selective processes interactively affected the ARG composition dynamics of the FL and PA bacteria. However, the antibiotic resistant community of FL bacteria exhibited a higher level of temporal stochasticity following the bloom event than PA bacteria. Therefore, we emphasized the bacterial lifestyles as an important mechanism, giving rise to different responses of antibiotic resistant community to the cyanobacterial bloom.

Importance of Local and Regional Scales in Shaping Mycobacterial Abundance in Freshwater Lakes

Biogeographical studies considering the entire bacterial community may underestimate mechanisms of bacterial assemblages at lower taxonomic levels. In this context, the study aimed to identify factors affecting the spatial and temporal dynamic of the Mycobacterium, a genus widespread in aquatic ecosystems. Nontuberculous mycobacteria (NTM) density variations were quantified in the water column of freshwater lakes at the regional scale (annual monitoring of 49 lakes in the Paris area) and at the local scale (2-year monthly monitoring in Créteil Lake) by real-time quantitative PCR targeting the atpE gene. At the regional scale, mycobacteria densities in water samples ranged from 6.7 × 10³ to 1.9 × 10⁸ genome units per liter. Density variations were primarily explained by water pH, labile iron, and dispersal processes through the connection of the lakes to a river. In Créteil Lake, no spatial variation of mycobacterial densities was noticed over the 2-year monthly survey, except after large rainfall events. Indeed, storm sewer effluents locally and temporarily increased NTM densities in the water column. The temporal dynamic of the NTM densities in Créteil Lake was associated with suspended solid concentrations. No clear seasonal variation was noticed despite a shift in NTM densities observed over the 2012-2013 winter. Temporal NTM densities fluctuations were well predicted by the neutral community model, suggesting a random balance between loss and gain of mycobacterial taxa within Créteil Lake. This study highlights the importance of considering multiple spatial scales for understanding the spatio-temporal dynamic of bacterial populations in natural environments.

Ecological patterns and adaptability of bacterial communities in alkaline copper mine drainage

Environmental gradient have strong effects on community assembly processes. In order to reveal the effects of alkaline mine drainage (AlkMD) on bacterial and denitrifying bacterial community compositions and diversity in tailings reservoir, here we conducted an experiment to examine all and core bacterial taxa and denitrifying functional genes's (nirS, nirK, nosZ_Ι) abundance along a chemical gradient in tailings water in Shibahe copper tailings in Zhongtiaoshan, China. Differences in bacterial and denitrifying bacterial community compositions in different habitats and their relationships with environmental parameters were analyzed. The results showed that the richness and diversity of bacterial community in downstream seeping water (SDSW) were the largest, while that in upstream tailings water (STW1) were the lowest. The diversity and abundance of bacterial communities tended to increase from STW1 to SDSW. The variation of bacterial community diversity was significantly related to electroconductibility (EC), nitrate (NO₃^-), nitrite (NO₂^-), total carbon (TC), inorganic carbon (IC) and sulfate (SO₄^2-), but was not correlated with geographic distance in local scale. Core taxa from class to genus were all significantly related to NO₃^- and NO₂^-. Core taxa Rhodobacteraceae, Rhodobacter, Acinetobacter and Hydrogenophaga were typical denitrifying bacteria. The variation trends of these groups were consistent with the copy number of nirS, nirK and nosZ_Ι, demonstrating their importance in the process of nitrogen reduction. The copy number of nirK, nosZ_Ι and nirS/16S rDNA, nirK/16Sr DNA correlated strongly with NO₃^-, NO₂^- and IC, but nirS and nosZ_I/16SrDNA had no significant correlation with NO₃^- and NO₂^-. The copy numbers of denitrifying functional genes (nirS, nirK and nosZ_Ι) were negatively correlated with heavy metal plumbum (Pb) and zinc (Zn). It showed that heavy metal contamination was an important factor affecting the structure of denitrifying bacterial community in AlkMD. In this study we have identified the distribution pattern of bacterial community along physiochemical gradients in alkaline tailings reservoir and displayed the driving force of shaping the structure of bacterial community. The influence of NO₃^-, NO₂^-, IC and heavy metal Pb and Zn on bacterial community might via their influence on the functional groups involving nitrogen, carbon and metal metabolisms.

Patterns and Processes in Marine Microeukaryotic Community Biogeography from Xiamen Coastal Waters and Intertidal Sediments, Southeast China

Microeukaryotes play key roles in the structure and functioning of marine ecosystems. Little is known about the relative importance of the processes that drive planktonic and benthic microeukaryotic biogeography in subtropical offshore areas. This study compares the microeukaryotic community compositions (MCCs) from offshore waters (n = 12) and intertidal sediments (n = 12) around Xiamen Island, southern China, using high-throughput sequencing of 18S rDNA. This work further quantifies the relative contributions of spatial and environmental variables on the distribution of marine MCCs (including total, dominant, rare and conditionally rare taxa). Our results showed that planktonic and benthic MCCs were significantly different, and the benthic richness (6627 OTUs) was much higher than that for plankton (4044 OTUs) with the same sequencing effort. Further, we found that benthic MCCs exhibited a significant distance-decay relationship, whereas the planktonic communities did not. After removing two unique sites (N2 and N3), however, 72% variation in planktonic community was explained well by stochastic processes. More importantly, both the environmental and spatial factors played significant roles in influencing the biogeography of total and dominant planktonic and benthic microeukaryotic communities, although their relative effects on these community variations were different. However, a high proportion of unexplained variation in the rare taxa (78.1–97.4%) and conditionally rare taxa (49.0–81.0%) indicated that more complex mechanisms may influence the assembly of the rare subcommunity. These results demonstrate that patterns and processes in marine microeukaryotic community assembly differ among the different habitats (coastal water vs. intertidal sediment) and different communities (total, dominant, rare and conditionally rare microeukaryotes), and provide novel insight on the microeukaryotic biogeography and ecological mechanisms in coastal waters and intertidal sediments at local scale.

Shotgun mitogenomics across body size classes in a local assemblage of tropical Diptera: Phylogeny, species diversity and mitochondrial abundance spectrum

Mitochondrial genomes can be assembled readily from shotgun-sequenced DNA mixtures of mass-trapped arthropods ("mitochondrial metagenomics"), speeding up the taxonomic characterization. Bulk sequencing was conducted on some 800 individuals of Diptera obtained by canopy fogging of a single tree in Borneo dominated by small (<1.5 mm) individuals. Specimens were split into five body size classes for DNA extraction, to equalize read numbers across specimens and to study how body size, a key ecological trait, interacts with species and phylogenetic diversity. Genome assembly produced 304 orthologous mitochondrial contigs presumed to each represent a different species. The small-bodied fraction was the by far most species-rich (187 contigs). Identification of contigs was through phylogenetic analysis together with 56 reference mitogenomes, which placed most of the Bornean community into seven clades of small-bodied species, indicating phylogenetic conservation of body size. Mapping of shotgun reads against the mitogenomes showed wide ranges of read abundances within each size class. Ranked read abundance plots were largely log-linear, indicating a uniformly filled abundance spectrum, especially for small-bodied species. Small-bodied species differed greatly from other size classes in neutral metacommunity parameters, exhibiting greater levels of immigration, besides greater total community size. We suggest that the established uses of mitochondrial metagenomics for analysis of species and phylogenetic diversity can be extended to parameterize recent theories of community ecology and biodiversity, and by focusing on the number mitochondria, rather than individuals, a new theoretical framework for analysis of mitochondrial abundance spectra can be developed that incorporates metabolic activity approximated by the count of mitochondria.

Inferring neutral biodiversity parameters using environmental DNA data sets

The DNA present in the environment is a unique and increasingly exploited source of information for conducting fast and standardized biodiversity assessments for any type of organisms. The datasets resulting from these surveys are however rarely compared to the quantitative predictions of biodiversity models. In this study, we simulate neutral taxa-abundance datasets, and artificially noise them by simulating noise terms typical of DNA-based biodiversity surveys. The resulting noised taxa abundances are used to assess whether the two parameters of Hubbell's neutral theory of biodiversity can still be estimated. We find that parameters can be inferred provided that PCR noise on taxa abundances does not exceed a certain threshold. However, inference is seriously biased by the presence of artifactual taxa. The uneven contribution of organisms to environmental DNA owing to size differences and barcode copy number variability does not impede neutral parameter inference, provided that the number of sequence reads used for inference is smaller than the number of effectively sampled individuals. Hence, estimating neutral parameters from DNA-based taxa abundance patterns is possible but requires some caution. In studies that include empirical noise assessments, our comprehensive simulation benchmark provides objective criteria to evaluate the robustness of neutral parameter inference.