Abstract: Manipulating microorganisms to increase soil organic carbon (SOC) in croplands remains a challenge. Soil microbes are important drivers of SOC sequestration, especially via their necromass accumulation. However, microbial parameters are rarely used to predict cropland SOC stocks, possibly due to uncertainties regarding the relationships between microbial carbon pools, community properties and SOC. Herein we evaluated the microbial community properties (diversity and network complexity), microbial carbon pools (biomass and necromass carbon) and SOC in 468 cropland soils across northeast China. We found that not only microbial necromass carbon but also microbial community properties (diversity and network complexity) and biomass carbon were correlated with SOC. Microbial biomass carbon and diversity played more important role in predicting SOC for maize, while microbial network complexity was more important for rice. Models to predict SOC performed better when the microbial community and microbial carbon pools were included simultaneously. Taken together our results suggest that microbial carbon pools and community properties influence SOC accumulation in croplands, and management practices that improve these microbial parameters may increase cropland SOC levels.

Abstract: Abstract Ecosystem functions and services are under threat from anthropogenic global change at a planetary scale. Microorganisms are the dominant drivers of nearly all ecosystem functions and therefore ecosystem-scale responses are dependent on responses of resident microbial communities. However, the specific characteristics of microbial communities that contribute to ecosystem stability under anthropogenic stress are unknown. We evaluated bacterial drivers of ecosystem stability by generating wide experimental gradients of bacterial diversity in soils, applying stress to the soils, and measuring responses of several microbial-mediated ecosystem processes, including C and N cycling rates and soil enzyme activities. Some processes (e.g., C mineralization) exhibited positive correlations with bacterial diversity and losses of diversity resulted in reduced stability of nearly all processes. However, comprehensive evaluation of all potential bacterial drivers of the processes revealed that bacterial α diversity per se was never among the most important predictors of ecosystem functions. Instead, key predictors included total microbial biomass, 16S gene abundance, bacterial ASV membership, and abundances of specific prokaryotic taxa and functional groups (e.g., nitrifying taxa). These results suggest that bacterial α diversity may be a useful indicator of soil ecosystem function and stability, but that other characteristics of bacterial communities are stronger statistical predictors of ecosystem function and better reflect the biological mechanisms by which microbial communities influence ecosystems. Overall, our results provide insight into the role of microorganisms in supporting ecosystem function and stability by identifying specific characteristics of bacterial communities that are critical for understanding and predicting ecosystem responses to global change.

Abstract: Developing low-cadmium (Cd) rice cultivars has emerged as a promising avenue for food safety in Cd-contaminated farmlands. The root-associated microbiomes of rice have been shown to enhance rice growth and alleviate Cd stress. However, the microbial taxon-specific Cd resistance mechanisms underlying different Cd accumulation characteristics between different rice cultivars remain largely unknown. This study compared low-Cd cultivar XS14 and hybrid rice cultivar YY17 for Cd accumulation with five soil amendments. The results showed that XS14 was characterized by more variable community structures and stable co-occurrence networks in the soil-root continuum compared to YY17. The stronger stochastic processes in assembly of the XS14 (~25%) rhizosphere community than that of YY17 (~12%) suggested XS14 may have higher resistance to changes in soil properties. Microbial co-occurrence networks and machine learning models jointly identified keystone indicator microbiota, such as Desulfobacteria in XS14 and Nitrospiraceae in YY17. Meanwhile, genes involved in sulfur cycling and nitrogen cycling were observed among the root-associated microbiome of these two cultivars, respectively. Microbiomes in the rhizosphere and root of XS14 showed a higher diversity in functioning, with the significant enrichment of functional genes related to amino acid and carbohydrate transport and metabolism, and sulfur cycling. Our findings revealed differences and similarities in the microbial communities associated with two rice cultivars, as well as bacterial biomarkers predictive of Cd-accumulation capacity. Thus, we provide new insights into taxon-specific recruitment strategies of two rice cultivars under Cd stress and highlight the utility of biomarkers in offering clues for enhancing crop resilience to Cd stresses in the future.

Abstract: Oxygen is a key regulatory factor of anaerobic ammonium oxidation (anammox). Although the inhibitory effect of oxygen is evident, a wide range of oxygen sensitivities of anammox bacteria have been reported so far, which makes it difficult to model the marine nitrogen loss and design anammox-based technologies. Here, oxygen tolerance and detoxification mechanisms of four genera of anammox bacteria; one marine species ("Ca. Scalindua sp.") and four freshwater anammox species ("Ca. Brocadia sinica", "Ca. Brocadia sapporoensis", "Ca. Jettenia caeni", and "Ca. Kuenenia stuttgartiensis") were determined and then related to the activities of anti-oxidative enzymes. Highly enriched planktonic anammox cells were exposed to various levels of oxygen, and oxygen inhibition kinetics (50% inhibitory concentration (IC50) and upper O2 limits (DOmax) of anammox activity) were quantitatively determined. A marine anammox species, "Ca. Scalindua sp.", exhibited much higher oxygen tolerance capability (IC50 = 18.0 µM and DOmax = 51.6 µM) than freshwater species (IC50 = 2.7-4.2 µM and DOmax = 10.9-26.6 µM). The upper DO limit of "Ca. Scalindua sp." was much higher than the values reported so far (~20 µM). Furthermore, the oxygen inhibition was reversible even after exposed to ambient air for 12-24 h. The comparative genome analysis confirmed that all anammox species commonly possess the genes considered to function for reduction of O2, superoxide anion (O2•-), and H2O2. However, the superoxide reductase (Sor)-peroxidase dependent detoxification system alone may not be sufficient for cell survival under microaerobic conditions. Despite the fact that anaerobes normally possess no or little superoxide dismutase (Sod) or catalase (Cat), only Scalindua exhibited high Sod activity of 22.6 ± 1.9 U/mg-protein with moderate Cat activity of 1.6 ± 0.7 U/mg-protein, which was consistent with the genome sequence analysis. This Sod-Cat dependent detoxification system could be responsible for the higher O2 tolerance of Scalindua than other freshwater anammox species lacking the Sod activity.

Abstract: Abstract Inferring the functional capabilities of bacteria from metagenome-assembled genomes (MAGs) is becoming a central process in microbiology. Here we show that the completeness of genomes has a significant impact on the recovered functional signal, spanning all domains of metabolic functions. We identify factors that affect this relationship between genome completeness and function fullness, and provide baseline knowledge to guide efforts to correct for this overlooked bias in metagenomic functional inference.

Abstract: Crop breeding has traditionally ignored the plant-associated microbial communities. Considering the interactions between plant genotype and associated microbiota is of value since different genotypes of the same crop often harbor distinct microbial communities which can influence the plant phenotype. However, recent studies have reported contrasting results, which led us to hypothesize that the effect of genotype is constrained by growth stages, sampling year and plant compartment. To test this hypothesis, we sampled bulk soil, rhizosphere soil and roots of 10 field-grown wheat genotypes, twice per year, for 4 years. DNA was extracted and regions of the bacterial 16 S rRNA and CPN60 genes and the fungal ITS region were amplified and sequenced. The effect of genotype was highly contingent on the time of sampling and on the plant compartment sampled. Only for a few sampling dates, were the microbial communities significantly different across genotypes. The effect of genotype was most often significant for root microbial communities. The three marker genes used provided a highly coherent picture of the effect of genotype. Taken together, our results confirm that microbial communities in the plant environment strongly vary across compartments, growth stages, and years, and that this can mask the effect of genotype.

Abstract: Stony coral tissue loss disease (SCTLD) has been causing significant whole colony mortality on reefs in Florida and the Caribbean. The cause of SCTLD remains unknown, with the limited concurrence of SCTLD-associated bacteria among studies. We conducted a meta-analysis of 16S ribosomal RNA gene datasets generated by 16 field and laboratory SCTLD studies to find consistent bacteria associated with SCTLD across disease zones (vulnerable, endemic, and epidemic), coral species, coral compartments (mucus, tissue, and skeleton), and colony health states (apparently healthy colony tissue (AH), and unaffected (DU) and lesion (DL) tissue from diseased colonies). We also evaluated bacteria in seawater and sediment, which may be sources of SCTLD transmission. Although AH colonies in endemic and epidemic zones harbor bacteria associated with SCTLD lesions, and aquaria and field samples had distinct microbial compositions, there were still clear differences in the microbial composition among AH, DU, and DL in the combined dataset. Alpha-diversity between AH and DL was not different; however, DU showed increased alpha-diversity compared to AH, indicating that, prior to lesion formation, corals may undergo a disturbance to the microbiome. This disturbance may be driven by Flavobacteriales, which were especially enriched in DU. In DL, Rhodobacterales and Peptostreptococcales-Tissierellales were prominent in structuring microbial interactions. We also predict an enrichment of an alpha-toxin in DL samples which is typically found in Clostridia. We provide a consensus of SCTLD-associated bacteria prior to and during lesion formation and identify how these taxa vary across studies, coral species, coral compartments, seawater, and sediment.

Abstract: Abstract Soil microbiomes are sensitive to current and previous soil conditions, and bacterial ‘bioindicators’ of biological, physical, and chemical soil properties have considerable potential for soil health assessment. However, the lack of ecological or physiological information for most soil microorganisms limits our ability to interpret the associations of bioindicators and, thus, their utility for guiding management. We identified bioindicators of tillage intensity and twelve soil properties used to rate soil health using a 16S rRNA gene-based survey of farmland across North America. We then inferred the genomic traits of bioindicators and evaluated their environment-wide associations (EWAS) with respect to agricultural management practice, disturbance, and plant associations with 89 studies from agroecosystems. Most bioindicators were either positively correlated with biological properties (e.g., organic matter) or negatively correlated with physical and chemical properties. Higher soil health ratings corresponded with smaller genome size and higher coding density, while lower ratings corresponded with larger genomes and higher rrn copy number. Community-weighted genome size explained most variation in health ratings. EWAS linked prominent bioindicators with the impacts of environmental disturbances. Our findings provide ecological insights into bioindicators of soil properties relevant to soil health management, illustrating the tight coupling of microbiome and soil function.

Abstract: Emiliania huxleyi is a unicellular micro-alga that forms massive oceanic blooms and plays key roles in global biogeochemical cycles. Mounting studies demonstrate various stimulatory and inhibitory influences that bacteria have on the E. huxleyi physiology. To investigate these algal-bacterial interactions, laboratory co-cultures have been established by us and by others. Owing to these co-cultures, various mechanisms of algal-bacterial interactions have been revealed, many involving bacterial pathogenicity towards algae. However, co-cultures represent a significantly simplified system, lacking the complexity of bacterial communities. In order to investigate bacterial pathogenicity within an ecologically relevant context, it becomes imperative to enhance the microbial complexity of co-culture setups. Phaeobacter inhibens bacteria are known pathogens that cause the death of E. huxleyi algae in laboratory co-culture systems. The bacteria depend on algal exudates for growth, but when algae senesce, bacteria switch to a pathogenic state and induce algal death. Here we investigate whether P. inhibens bacteria can induce algal death in the presence of a complex bacterial community. We show that an E. huxleyi-associated bacterial community protects the alga from the pathogen, although the pathogen occurs within the community. To study how the bacterial community regulates pathogenicity, we reduced the complex bacterial community to a five-member synthetic community (syncom). The syncom is comprised of a single algal host and five isolated bacterial species, which represent major bacterial groups that are naturally associated with E. huxleyi. We discovered that a single bacterial species in the reduced community, Sulfitobacter pontiacus, protects the alga from the pathogen. We further found that algal protection from P. inhibens pathogenicity is a shared trait among several Sulfitobacter species. Algal protection by bacteria might be a common phenomenon with ecological significance, which is overlooked in reduced co-culture systems.

Abstract: The environmental impacts of genetically modified (GM) plants remain a controversial global issue. To address these issues, comprehensive environmental risk assessments of GM plants is critical for the sustainable development and application of transgenic technology. In this paper, significant differences were not observed between microbial metagenomic and metabolomic profiles in surface waters of the Bt rice (T1C-1, the transgenic line) and non-Bt cultivars (Minghui 63 (the isogenic line) and Zhonghua 11 (the conventional japonica cultivar)). In contrast, differences in these profiles were apparent in the rhizospheres. T1C-1 planting increased soil microbiome diversity and network stability, but did not significantly alter the abundances of potential probiotic or phytopathogenic microorganisms compared with Minghui 63 and Zhonghua 11, which revealed no adverse effects of T1C-1 on soil microbial communities. T1C-1 planting could significantly alter soil C and N, probably via the regulation of the abundances of enzymes related to soil C and N cycling. In addition, integrated multi-omic analysis of root exudate metabolomes and soil microbiomes showed that the abundances of various metabolites released as root exudates were significantly correlated with subsets of microbial populations including the Acidobacteria, Actinobacteria, Chloroflexi, and Gemmatimonadetes that were differentially abundant in T1C-1 and Mnghui 63 soils. Finally, the potential for T1C-1-associated root metabolites to exert growth effects on T1C-1-associated species was experimentally validated by analysis of bacterial cultures, revealing that Bt rice planting could selectively modulate specific root microbiota. Overall, this study indicate that Bt rice can directly modulate rhizosphere microbiome assemblages by altering the metabolic compositions of root exudates that then alters soil metabolite profiles and physiochemical properties. This study unveils the mechanistic associations of Bt plant-microorganism-environment, which provides comprehensive insights into the potential ecological impacts of GM plants.

Abstract: For a deeper and comprehensive understanding of the composition and function of rhizosphere microbiomes, we need to focus at the scale of individual roots in standardized growth containers. Root exudation patterns are known to vary along distinct parts of the root even in juvenile plants giving rise to spatially distinct microbial niches. To address this, we analyzed the microbial community from two spatially distinct zones of the developing primary root (tip and base) in young Brachypodium distachyon grown in natural soil using standardized fabricated ecosystems known as EcoFABs as well as in more conventional pot and tubes. 16S rRNA based community analysis showed a strong rhizosphere effect resulting in significant enrichment of several OTUs belonging to Actinobacteria, Bacteroidetes, Firmicutes and Proteobacteria. However, microbial community composition did not differ between root tips and root base or across different growth containers. Functional analysis of bulk metagenomics revealed significant differences between root tips and bulk soil. The genes associated with different metabolic pathways and root colonization were enriched in root tips. On the other hand, genes associated with nutrient-limitation and environmental stress were prominent in the bulk soil compared to root tips, implying the absence of easily available, labile carbon and nutrients in bulk soil relative to roots. Such insights into the relationships between developing root and microbial communities are critical for judicious understanding of plant-microbe interactions in early developmental stages of plants.

Abstract: Marine sediments comprise one of the largest environments on the planet, and their microbial inhabitants are significant players in global carbon and nutrient cycles. Recent studies using metagenomic techniques have shown the complexity of these communities and identified novel microorganisms from the ocean floor. Here, we obtained 77 metagenome-assembled genomes (MAGs) from the bacterial phylum Armatimonadota in the Guaymas Basin, Gulf of California, and the Bohai Sea, China. These MAGs comprise two previously undescribed classes within Armatimonadota, which we propose naming Hebobacteria and Zipacnadia. They are globally distributed in hypoxic and anoxic environments and are dominant members of deep-sea sediments (up to 1.95% of metagenomic raw reads). The classes described here also have unique metabolic capabilities, possessing pathways to reduce carbon dioxide to acetate via the Wood-Ljungdahl pathway (WLP) and generating energy through the oxidative branch of glycolysis using carbon dioxide as an electron sink, maintaining the redox balance using the WLP. Hebobacteria may also be autotrophic, not previously identified in Armatimonadota. Furthermore, these Armatimonadota may play a role in sulfur and nitrogen cycling, using the intermediate compounds hydroxylamine and sulfite. Description of these MAGs enhances our understanding of diversity and metabolic potential within anoxic habitats worldwide.

Abstract: Forests are composed of various plant species, and rhizosphere soil microbes are driven by root exudates. However, the interplay between root exudates, microbial communities in the rhizosphere soil of canopy trees, understory shrubs, grasses, and their responses to nitrogen (N) deposition remains unclear. Pinus tabulaeformis, Rosa xanthina, and Carex lancifolia were used to investigate root exudates, rhizosphere soil microbial communities, and their responses to N application in forest ecosystem. Root exudate abundances of P. tabulaeformis were significantly higher than that of R. xanthina and C. lancifolia, with carbohydrates dominating P. tabulaeformis and R. xanthina root exudates, fatty acids prevailing in C. lancifolia root exudates. Following N application, root exudate abundances of P. tabulaeformis and R. xanthina initially increased before decreasing, whereas those of C. lancifolia decreased. Microbial number of rhizosphere soil of C. lancifolia was higher than that of P. tabulaeformis and R. xanthina, but there was insignificant variation of rhizosphere soil microbial diversity among plant species. N application exerted promotional and inhibitory impacts on bacterial and fungal numbers, respectively, while bacterial and fungal diversities were increased by N application. Overall, N application had negative effects on root exudates of P. tabulaeformis, inhibiting rhizosphere soil microbial populations. N application suppressed rhizosphere soil microbial populations by increasing root exudates of R. xanthina. Conversely, N application elevated nutrient content in the rhizosphere soil of C. lancifolia, reducing root exudates and minimally promoting microbial populations. This study highlights the importance of understory vegetation in shaping soil microbial communities within forests under N deposition.

Abstract: Acinetobacter is present in the livestock environment, but little is known about their antibiotic resistance and pathogenic species in the farm groundwater. Here we investigated antibiotic resistance of Acinetobacter in the swine farm groundwater (JZPG) and residential groundwater (JZG) of a swine farming village, in comparison to a nearby (3.5 km) non-farming village (WTG) using metagenomic and culture-based approaches. Results showed that the abundance of antibiotic resistome in some JZG and all JZPG (~3.4 copies/16S rRNA gene) was higher than that in WTG (~0.7 copies/16S rRNA gene), indicating the influence of farming activities on both groundwater types. Acinetobacter accounted for ~95.7% of the bacteria in JZG and JZPG, but only ~8.0% in WTG. They were potential hosts of ~95.6% of the resistome in farm affected groundwater, which includes 99 ARG subtypes against 23 antibiotic classes. These ARGs were associated with diverse intrinsic and acquired resistance mechanisms, and the predominant ARGs were tetracyclines and fluoroquinolones resistance genes. Metagenomic binning analysis elucidated that non-baumannii Acinetobacter including A. oleivorans, A. beijerinckii, A. seifertii, A. bereziniae and A. modestus might pose environmental risks because of multidrug resistance, pathogenicity and massive existence in the groundwater. Antibiotic susceptibility tests showed that the isolated strains were resistant to multiple antibiotics including sulfamethoxazole (resistance ratio: 96.2%), levofloxacin (42.5%), gatifloxacin (39.0%), ciprofloxacin (32.6%), tetracycline (32.0%), doxycycline (29.0%) and ampicillin (12.0%) as well as last-resort polymyxin B (31.7%), colistin (24.1%) and tigecycline (4.1%). The findings highlight potential prevalence of groundwater-borne antibiotic-resistant pathogenic Acinetobacter in the livestock environment.

Abstract: Microbiome analysis through 16S rRNA gene sequencing is a crucial tool for understanding the microbial ecology of any habitat or ecosystem. However, workflows require large equipment, stable internet, and extensive computing power such that most of the work is performed far away from sample collection in both space and time. Performing amplicon sequencing and analysis at sample collection would have positive implications in many instances including remote fieldwork and point-of-care medical diagnoses. Here we present SituSeq, an offline and portable workflow for the sequencing and analysis of 16S rRNA gene amplicons using Nanopore sequencing and a standard laptop computer. SituSeq was validated by comparing Nanopore 16S rRNA gene amplicons, Illumina 16S rRNA gene amplicons, and Illumina metagenomes, sequenced using the same environmental DNA. Comparisons revealed consistent community composition, ecological trends, and sequence identity across platforms. Correlation between the abundance of taxa in each taxonomic level in Illumina and Nanopore data sets was high (Pearson's r > 0.9), and over 70% of Illumina 16S rRNA gene sequences matched a Nanopore sequence with greater than 97% sequence identity. On board a research vessel on the open ocean, SituSeq was used to analyze amplicon sequences from deep sea sediments less than 2 h after sequencing, and 8 h after sample collection. The rapidly available results informed decisions about subsequent sampling in near real-time while the offshore expedition was still underway. SituSeq is a portable and user-friendly workflow that helps to bring the power of microbial genomics and diagnostics to many more researchers and situations.

Abstract: Abstract 16S rRNA gene copy number (16S GCN) varies among bacterial species and this variation introduces potential biases to microbial diversity analyses using 16S rRNA read counts. To correct the biases, methods have been developed to predict 16S GCN. A recent study suggests that the prediction uncertainty can be so great that copy number correction is not justified in practice. Here we develop RasperGade16S, a novel method and software to better model and capture the inherent uncertainty in 16S GCN prediction. RasperGade16S implements a maximum likelihood framework of pulsed evolution model and explicitly accounts for intraspecific GCN variation and heterogeneous GCN evolution rates among species. Using cross-validation, we show that our method provides robust confidence estimates for the GCN predictions and outperforms other methods in both precision and recall. We have predicted GCN for 592605 OTUs in the SILVA database and tested 113842 bacterial communities that represent an exhaustive and diverse list of engineered and natural environments. We found that the prediction uncertainty is small enough for 99% of the communities that 16S GCN correction should improve their compositional and functional profiles estimated using 16S rRNA reads. On the other hand, we found that GCN variation has limited impacts on beta-diversity analyses such as PCoA, NMDS, PERMANOVA and random-forest test.

Abstract: Aquifers are populated by highly diverse microbial communities, including unusually small bacteria and archaea. The recently described Patescibacteria (or Candidate Phyla Radiation) and DPANN radiation are characterized by ultra-small cell and genomes sizes, resulting in limited metabolic capacities and probable dependency on other organisms to survive. We applied a multi-omics approach to characterize the ultra-small microbial communities over a wide range of aquifer groundwater chemistries. Results expand the known global range of these unusual organisms, demonstrate the wide geographical range of over 11,000 subsurface-adapted Patescibacteria, Dependentiae and DPANN archaea, and indicate that prokaryotes with ultra-small genomes and minimalistic metabolism are a characteristic feature of the terrestrial subsurface. Community composition and metabolic activities were largely shaped by water oxygen content, while highly site-specific relative abundance profiles were driven by a combination of groundwater physicochemistries (pH, nitrate-N, dissolved organic carbon). We provide insights into the activity of ultra-small prokaryotes with evidence that they are major contributors to groundwater community transcriptional activity. Ultra-small prokaryotes exhibited genetic flexibility with respect to groundwater oxygen content, and transcriptionally distinct responses, including proportionally greater transcription invested into amino acid and lipid metabolism and signal transduction in oxic groundwater, along with differences in taxa transcriptionally active. Those associated with sediments differed from planktonic counterparts in species composition and transcriptional activity, and exhibited metabolic adaptations reflecting a surface-associated lifestyle. Finally, results showed that groups of phylogenetically diverse ultra-small organisms co-occurred strongly across sites, indicating shared preferences for groundwater conditions.

Abstract: Abstract Viruses of the phylum Nucleocytoviricota are ubiquitous in ocean waters and play important roles in shaping the dynamics of marine ecosystems. In this study, we leveraged the bioGEOTRACES metagenomic dataset collected across the Atlantic and Pacific Oceans to investigate the biogeography of these viruses in marine environments. We identified 330 viral genomes, including 212 in the order Imitervirales and 54 in the order Algavirales . We found that most viruses appeared to be prevalent in shallow waters (<150 m), and that viruses of the Mesomimiviridae ( Imitervirales) and Prasinoviridae ( Algavirales) are by far the most abundant and diverse groups in our survey. Five mesomimiviruses and one prasinovirus are particularly widespread in oligotrophic waters; annotation of these genomes revealed common stress response systems, photosynthesis-associated genes, and oxidative stress modulation genes that may be key to their broad distribution in the pelagic ocean. We identified a latitudinal pattern in viral diversity in one cruise that traversed the North and South Atlantic Ocean, with viral diversity peaking at high latitudes of the northern hemisphere. Community analyses revealed three distinct Nucleocytoviricota communities across latitudes, categorized by latitudinal distance towards the equator. Our results contribute to the understanding of the biogeography of these viruses in marine systems.

Abstract: Abstract The application of meta-barcoding, qPCR, and metagenomics to aquatic eukaryotic microbial communities requires knowledge of genomic copy number variability (CNV). CNV may be particularly relevant to functional genes, impacting dosage and expression, yet little is known of the scale and role of CNV in microbial eukaryotes. Here, we quantify CNV of rRNA and a gene involved in Paralytic Shellfish Toxin (PST) synthesis ( sxtA4 ), in 51 strains of 4 Alexandrium (Dinophyceae) species. Genomes varied up to threefold within species and ~7-fold amongst species, with the largest ( A. pacificum , 130 ± 1.3 pg cell −1 /~127 Gbp) in the largest size category of any eukaryote. Genomic copy numbers (GCN) of rRNA varied by 6 orders of magnitude amongst Alexandrium (10 2 – 10 8 copies cell −1 ) and were significantly related to genome size. Within the population CNV of rRNA was 2 orders of magnitude (10 5 – 10 7 cell −1 ) in 15 isolates from one population, demonstrating that quantitative data based on rRNA genes needs considerable caution in interpretation, even if validated against locally isolated strains. Despite up to 30 years in laboratory culture, rRNA CNV and genome size variability were not correlated with time in culture. Cell volume was only weakly associated with rRNA GCN (20–22% variance explained across dinoflagellates, 4% in Gonyaulacales). GCN of sxtA4 varied from 0–10 2 copies cell −1 , was significantly related to PSTs (ng cell −1 ), displaying a gene dosage effect modulating PST production. Our data indicate that in dinoflagellates, a major marine eukaryotic group, low-copy functional genes are more reliable and informative targets for quantification of ecological processes than unstable rRNA genes.

Abstract: The oscillating redox conditions that characterize coastal sandy sediments foster microbial communities capable of respiring oxygen and nitrate simultaneously, thereby increasing the potential for organic matter remineralization, nitrogen (N)-loss and emissions of the greenhouse gas nitrous oxide. It is unknown to what extent these conditions also lead to overlaps between dissimilatory nitrate and sulfate respiration. Here, we show that sulfate and nitrate respiration co-occur in the surface sediments of an intertidal sand flat. Furthermore, we found strong correlations between dissimilatory nitrite reduction to ammonium (DNRA) and sulfate reduction rates. Until now, the nitrogen and sulfur cycles were assumed to be mainly linked in marine sediments by the activity of nitrate-reducing sulfide oxidisers. However, transcriptomic analyses revealed that the functional marker gene for DNRA (nrfA) was more associated with microorganisms known to reduce sulfate rather than oxidise sulfide. Our results suggest that when nitrate is supplied to the sediment community upon tidal inundation, part of the sulfate reducing community may switch respiratory strategy to DNRA. Therefore increases in sulfate reduction rate in-situ may result in enhanced DNRA and reduced denitrification rates. Intriguingly, the shift from denitrification to DNRA did not influence the amount of N2O produced by the denitrifying community. Our results imply that microorganisms classically considered as sulfate reducers control the potential for DNRA within coastal sediments when redox conditions oscillate and therefore retain ammonium that would otherwise be removed by denitrification, exacerbating eutrophication.

Abstract: Carbohydrates are chemically and structurally diverse, represent a substantial fraction of marine organic matter and are key substrates for heterotrophic microbes. Studies on carbohydrate utilisation by marine microbes have been centred on phytoplankton blooms in temperate regions, while far less is known from high-latitude waters and during later seasonal stages. Here, we combine glycan microarrays and analytical chromatography with metagenomics and metatranscriptomics to show the spatial heterogeneity in glycan distribution and potential carbohydrate utilisation by microbes in Atlantic waters of the Arctic. The composition and abundance of monomers and glycan structures in POM varied with location and depth. Complex fucose-containing sulfated polysaccharides, known to accumulate in the ocean, were consistently detected, while the more labile β-1,3-glucan exhibited a patchy distribution. Through 'omics analysis, we identify variations in the abundance and transcription of carbohydrate degradation-related genes across samples at the community and population level. The populations contributing the most to transcription were taxonomically related to those known as primary responders and key carbohydrate degraders in temperate ecosystems, such as NS4 Marine Group and Formosa. The unique transcription profiles for these populations suggest distinct substrate utilisation potentials, with predicted glycan targets corresponding to those structurally identified in POM from the same sampling sites. By combining cutting-edge technologies and protocols, we provide insights into the carbohydrate component of the carbon cycle in the Arctic during late summer and present a high-quality dataset that will be of great value for future comparative analyses.

Abstract: Heterotrophic bacteria in the ocean invest carbon, nitrogen, and energy in extracellular enzymes to hydrolyze large substrates to smaller sizes suitable for uptake. Since hydrolysis products produced outside of a cell may be lost to diffusion, the return on this investment is uncertain. Selfish bacteria change the odds in their favor by binding, partially hydrolyzing, and transporting polysaccharides into the periplasmic space without loss of hydrolysis products. We expected selfish bacteria to be most common in the upper ocean, where phytoplankton produce abundant fresh organic matter, including complex polysaccharides. We, therefore, sampled water in the western North Atlantic Ocean at four depths from three stations differing in physiochemical conditions; these stations and depths also differed considerably in microbial community composition. To our surprise, we found that selfish bacteria are common throughout the water column of the ocean, including at depths greater than 5500 m. Selfish uptake as a strategy thus appears to be geographically-and phylogenetically-widespread. Since processing and uptake of polysaccharides require enzymes that are highly sensitive to substrate structure, the activities of these bacteria might not be reflected by measurements relying on uptake only of low molecular weight substrates. Moreover, even at the bottom of the ocean, the supply of structurally-intact polysaccharides, and therefore the return on enzymatic investment, must be sufficient to maintain these organisms.

Abstract: Abstract Marine protists are major components of the oceanic microbiome that remain largely unrepresented in culture collections and genomic reference databases. The exploration of this uncharted protist diversity in oceanic communities relies essentially on studying genetic markers from the environment as taxonomic barcodes. Here we report that across 6 large scale spatio-temporal planktonic surveys, half of the genetic barcodes remain taxonomically unassigned at the genus level, preventing a fine ecological understanding for numerous protist lineages. Among them, parasitic Syndiniales (Dinoflagellata) appear as the least described protist group. We have developed a computational workflow, integrating diverse 18S rDNA gene metabarcoding datasets, in order to infer large-scale ecological patterns at 100% similarity of the genetic marker, overcoming the limitation of taxonomic assignment. From a spatial perspective, we identified 2171 unassigned clusters, i.e., Syndiniales sequences with 100% similarity, exclusively shared between the Tropical/Subtropical Ocean and the Mediterranean Sea among all Syndiniales orders and 25 ubiquitous clusters shared within all the studied marine regions. From a temporal perspective, over 3 time-series, we highlighted 39 unassigned clusters that follow rhythmic patterns of recurrence and are the best indicators of parasite community’s variation. These clusters withhold potential as ecosystem change indicators, mirroring their associated host community responses. Our results underline the importance of Syndiniales in structuring planktonic communities through space and time, raising questions regarding host-parasite association specificity and the trophic mode of persistent Syndiniales, while providing an innovative framework for prioritizing unassigned protist taxa for further description.

Abstract: Cyanobacterial mats inhabit a variety of aquatic habitats, including the most extreme environments on Earth. They can thrive in a wide range of phosphorus (P) levels and are thus important players for ecosystem primary production and P cycling at the sediment-water interface. Polyphosphate (polyP), the major microbial P storage molecule, is assigned a critical role in compensating for phosphate fluctuations in planktonic cyanobacteria, but little is known about potentially analogous mechanisms of mat-forming cyanobacteria. To investigate acclimation strategies of cyanobacterial mats to fluctuating phosphate concentrations, laboratory batch experiments were conducted, in which the cosmopolitan mat-forming, marine cyanobacterium Sodalinema stali was exposed to low dissolved P concentrations, followed by a P pulse. Our results show that the cyanobacteria dynamically adjusted cellular P content to ambient phosphate concentrations and that they had accumulated polyP during periods of high phosphate availability, which was subsequently recycled to sustain growth during phosphate scarcity. However, following the depletion of dispensable cellular P sources, including polyP, we observed a reallocation of P contained in DNA into polyP, accompanied by increasing alkaline phosphatase activity. This suggests a change of the metabolic focus from growth towards maintenance and the attempt to acquire organic P, which would be naturally contained in the sediment. P overplus uptake following a simulated P pulse further suggests that Sodalinema-dominated mats exhibit elaborated mechanisms to cope with severe P fluctuations to overcome unfavourable environmental conditions, and potentially modulate critical P fluxes in the aquatic cycle.

Abstract: The oceanic igneous crust is a vast reservoir for microbial life, dominated by diverse and active bacteria, archaea, and fungi. Archaeal and bacterial viruses were previously detected in oceanic crustal fluids at the Juan de Fuca Ridge (JdFR). Here we report the discovery of two eukaryotic Nucleocytoviricota genomes from the same crustal fluids by sorting and sequencing single virions. Both genomes have a tRNATyr gene with an intron (20 bps) at the canonical position between nucleotide 37 and 38, a common feature in eukaryotic and archaeal tRNA genes with short introns (<100 bps), and fungal genes acquired through horizontal gene transfer (HGT) events. The dominance of Ascomycota fungi as the main eukaryotes in crustal fluids and the evidence for HGT point to these fungi as the putative hosts, making these the first putative fungi-Nucleocytoviricota specific association. Our study suggests active host-viral dynamics for the only eukaryotic group found in the subsurface oceanic crust and raises important questions about the impact of viral infection on the productivity and biogeochemical cycling in this ecosystem.

Abstract: Abstract The inter-individual variations of gut microbiome contribute to the different responses toward drug therapy among populations, developing a reliable ex vivo culture method for mixed bacteria is the urgent need for predicting personal reaction to drug therapy. Unfortunately, very few attentions have been paid to the bias that could be introduced during the culture process for mixed bacteria. Here we systemically evaluated the factors that may affect the outcomes of cultured bacteria from human feces. We demonstrated that inter-individual difference of host gut microbiome was the main factor affecting the outcomes of cultured bacteria, followed by the culture medium and time point. We further optimized a new medium termed GB based on our established multi-dimensional evaluation method, which could mimic the status of in situ host gut microbiome to the highest extent. Finally, we assessed the inter-individual metabolism by host gut microbiome from 10 donors on three frequently used clinical drugs (aspirin, levodopa and doxifluridine) based on the optimized GB medium. Our results revealed obvious variation in drug metabolism by microbiome from different donors, especially levodopa and doxifluridine. This work suggested the optimized culture medium had the potential for exploring the inter-individual impacts of host gut microbiome on drug metabolism.

Abstract: Due to the wide use of antibiotics, intensive aquaculture farms have been recognized as a significant reservoir of antibiotic resistomes. Although the prevalence of colistin resistance genes and multidrug-resistant bacteria (MDRB) has been documented, empirical evidence for the transmission of colistin and multidrug resistance between bacterial communities in aquaculture farms through horizontal gene transfer (HGT) is lacking. Here, we report the prevalence and transmission risk of colistin and multidrug resistance in 27 aquaculture water samples from 9 aquaculture zones from over 5000 km of subtropical coastlines in southern China. The colistin resistance gene mcr-1, mobile genetic element (MGE) intl1 and 13 typical antibiotic resistance genes (ARGs) were prevalent in all the aquaculture water samples. Most types of antibiotic (especially colistin) resistance are transmissible in bacterial communities based on evidence from laboratory conjugation and transformation experiments. Diverse MDRB were detected in most of the aquaculture water samples, and a strain with high-level colistin resistance, named Ralstonia pickettii MCR, was isolated. The risk of horizontal transfer of the colistin resistance of R. pickettii MCR through conjugation and transformation was low, but the colistin resistance could be steadily transmitted to offspring through vertical transfer. The findings have important implications for the future regulation of antibiotic use in aquaculture farms globally to address the growing threat posed by antibiotic resistance to human health.