Primer selection impacts specific population abundances but not community dynamics in a monthly time-series 16S rRNA gene amplicon analysis of coastal marine bacterioplankton

Multiple cropping effectively increases soil bacterial diversity, community abundance and soil fertility of paddy fields

BACKGROUND

Crop diversification is considered as an imperative approach for synchronizing the plant nutrient demands and soil nutrient availability. Taking two or more crops from the same field in one year is considered as multiple cropping. It improves the diversity and abundance of soil microbes, thereby improving the growth and yield of crops. Therefore, the present study was conducted to explore the effects of different multiple winter cropping on soil microbial communities in paddy fields. In this study, eight rice cropping patterns from two multiple cropping systems with three different winter crops, including Chinese milk vetch (CMV), rape, and wheat were selected. The effects of different multiple winter cropping on soil microbial abundance, community structure, and diversity in paddy fields were studied by 16 S rRNA high-throughput sequencing and real-time fluorescence quantitative polymerase chain reaction (PCR).

RESULTS

The results showed that different multiple winter cropping increased the operational taxonomic units (OTUs), species richness, and community richness index of the bacterial community in 0 ~ 20 cm soil layer. Moreover, soil physical and chemical properties of different multiple cropping patterns also affected the diversity and abundance of microbial bacterial communities. The multiple cropping increased soil potassium and nitrogen content, which significantly affected the diversity and abundance of bacterial communities, and it also increased the overall paddy yield. Moreover, different winter cropping changed the population distribution of microorganisms, and Proteobacteria, Acidobacteria, Nitrospira, and Chloroflexi were identified as the most dominant groups. Multiple winter cropping, especially rape-early rice-late rice (TR) andChinese milk vetch- early rice-late rice (TC) enhanced the abundance of Proteobacteria, Acidobacteria, and Actinobacteria and decreased the relative abundance of Verrucomicrobia and Euryarchaeota.

CONCLUSION

In conclusion, winter cropping of Chinese milk vetch and rape were beneficial to improve the soil fertility, bacteria diversity, abundance and rice yield.

Changes in microbial community structure related to biodegradation of eelgrass (Zostera marina)

Seagrass meadows produce organic carbon and deposit it on the seabed through the decaying process. Microbial activity is closely related to the process of eelgrass death and collapse. We investigated the microbial community structure of eelgrass during the eelgrass decomposition process by using a microcosm containing raw seawater and excised eelgrass leaves collected from a Zostera marina bed in Futtsu, Chiba Prefecture, Japan. The fast-growing microbes (i.e., Alphaproteobacteria, Gammaproteobacteria, and Flavobacteriia) rapidly adhered to the eelgrass leaf surface and proliferated in the first two weeks but gradually decreased the relative abundance as the months moved on. On the other hand, the slow-growing microbes (i.e., Cytophagia, Anaerolineae, Thaumarchaeota, and Actinobacteria) became predominant over the eelgrass surface late in the culture experiment (120, 180 days). The fast-growing groups of Gammaproteobacteria and Flavobacteriia appear to be closely related to the initial decomposition of eelgrass, especially the rapid decomposition of leaf-derived biopolymers. Changes in nitrogen content due to the bacterial rapid consumption of readily degradable organic carbon induced changes in the community structure at the early stage of eelgrass decomposition. In addition, shifts in the C/N ratio were driven by microbial community changes during later decomposition phases.

Gradients of bacteria in the oceanic water column reveal finely-resolved vertical distributions

Bacterial communities directly influence ecological processes in the ocean, and depth has a major influence due to the changeover in primary energy sources between the sunlit photic zone and dark ocean. Here, we examine the abundance and diversity of bacteria in Monterey Bay depth profiles collected from the surface to just above the sediments (e.g., 2000 m). Bacterial abundance in these Pacific Ocean samples decreased by >1 order of magnitude, from 1.22 ±0.69 ×106 cells ml-1 in the variable photic zone to 1.44 ± 0.25 ×105 and 6.71 ± 1.23 ×104 cells ml-1 in the mesopelagic and bathypelagic, respectively. V1-V2 16S rRNA gene profiling showed diversity increased sharply between the photic and mesopelagic zones. Weighted Gene Correlation Network Analysis clustered co-occurring bacterial amplicon sequence variants (ASVs) into seven subnetwork modules, of which five strongly correlated with depth-related factors. Within surface-associated modules there was a clear distinction between a 'copiotrophic' module, correlating with chlorophyll and dominated by e.g., Flavobacteriales and Rhodobacteraceae, and an 'oligotrophic' module dominated by diverse Oceanospirillales (such as uncultured JL-ETNP-Y6, SAR86) and Pelagibacterales. Phylogenetic reconstructions of Pelagibacterales and SAR324 using full-length 16S rRNA gene data revealed several additional subclades, expanding known microdiversity within these abundant lineages, including new Pelagibacterales subclades Ia.B, Id, and IIc, which comprised 4-10% of amplicons depending on the subclade and depth zone. SAR324 and Oceanospirillales dominated in the mesopelagic, with SAR324 clade II exhibiting its highest relative abundances (17±4%) in the lower mesopelagic (300-750 m). The two newly-identified SAR324 clades showed highest relative abundances in the photic zone (clade III), while clade IV was extremely low in relative abundance, but present across dark ocean depths. Hierarchical clustering placed microbial communities from 900 m samples with those from the bathypelagic, where Marinimicrobia was distinctively relatively abundant. The patterns resolved herein, through high resolution and statistical replication, establish baselines for marine bacterial abundance and taxonomic distributions across the Monterey Bay water column, against which future change can be assessed.

Direct observations of microbial community succession on sinking marine particles

Microbial community dynamics on sinking particles control the amount of carbon that reaches the deep ocean and the length of time that carbon is stored, with potentially profound impacts on Earth's climate. A mechanistic understanding of the controls on sinking particle distributions has been hindered by limited depth- and time-resolved sampling and methods that cannot distinguish individual particles. Here, we analyze microbial communities on nearly 400 individual sinking particles in conjunction with more conventional composite particle samples to determine how particle colonization and community assembly might control carbon sequestration in the deep ocean. We observed community succession with corresponding changes in microbial metabolic potential on the larger sinking particles transporting a significant fraction of carbon to the deep sea. Microbial community richness decreased as particles aged and sank; however, richness increased with particle size and the attenuation of carbon export. This suggests that the theory of island biogeography applies to sinking marine particles. Changes in POC flux attenuation with time and microbial community composition with depth were reproduced in a mechanistic ecosystem model that reflected a range of POC labilities and microbial growth rates. Our results highlight microbial community dynamics and processes on individual sinking particles, the isolation of which is necessary to improve mechanistic models of ocean carbon uptake.

Evaluation of multiple displacement amplification for metagenomic analysis of low biomass samples

Combining multiple displacement amplification (MDA) with metagenomics enables the analysis of samples with extremely low DNA concentrations, making them suitable for high-throughput sequencing. Although amplification bias and nonspecific amplification have been reported from MDA-amplified samples, the impact of MDA on metagenomic datasets is not well understood. We compared three MDA methods (i.e. bulk MDA, emulsion MDA, and primase MDA) for metagenomic analysis of two DNA template concentrations (approx. 1 and 100 pg) derived from a microbial community standard "mock community" and two low biomass environmental samples (i.e. borehole fluid and groundwater). We assessed the impact of MDA on metagenome-based community composition, assembly quality, functional profiles, and binning. We found amplification bias against high GC content genomes but relatively low nonspecific amplification such as chimeras, artifacts, or contamination for all MDA methods. We observed MDA-associated representational bias for microbial community profiles, especially for low-input DNA and with the primase MDA method. Nevertheless, similar taxa were represented in MDA-amplified libraries to those of unamplified samples. The MDA libraries were highly fragmented, but similar functional profiles to the unamplified libraries were obtained for bulk MDA and emulsion MDA at higher DNA input and across these MDA libraries for the groundwater sample. Medium to low-quality bins were possible for the high input bulk MDA metagenomes for the most simple microbial communities, borehole fluid, and mock community. Although MDA-based amplification should be avoided, it can still reveal meaningful taxonomic and functional information from samples with extremely low DNA concentration where direct metagenomics is otherwise impossible.

Suboxic DOM is bioavailable to surface prokaryotes in a simulated overturn of an oxygen minimum zone, Devil's Hole, Bermuda

Oxygen minimum zones (OMZs) are expanding due to increased sea surface temperatures, subsequent increased oxygen demand through respiration, reduced oxygen solubility, and thermal stratification driven in part by anthropogenic climate change. Devil's Hole, Bermuda is a model ecosystem to study OMZ microbial biogeochemistry because the formation and subsequent overturn of the suboxic zone occur annually. During thermally driven stratification, suboxic conditions develop, with organic matter and nutrients accumulating at depth. In this study, the bioavailability of the accumulated dissolved organic carbon (DOC) and the microbial community response to reoxygenation of suboxic waters was assessed using a simulated overturn experiment. The surface inoculated prokaryotic community responded to the deep (formerly suboxic) 0.2 μm filtrate with cell densities increasing 2.5-fold over 6 days while removing 5 μmol L-1 of DOC. After 12 days, the surface community began to shift, and DOC quality became less diagenetically altered along with an increase in SAR202, a Chloroflexi that can degrade recalcitrant dissolved organic matter (DOM). Labile DOC production after 12 days coincided with an increase of Nitrosopumilales, a chemoautotrophic ammonia oxidizing archaea (AOA) that converts ammonia to nitrite based on the ammonia monooxygenase (amoA) gene copy number and nutrient data. In comparison, the inoculation of the deep anaerobic prokaryotic community into surface 0.2 μm filtrate demonstrated a die-off of 25.5% of the initial inoculum community followed by a 1.5-fold increase in cell densities over 6 days. Within 2 days, the prokaryotic community shifted from a Chlorobiales dominated assemblage to a surface-like heterotrophic community devoid of Chlorobiales. The DOM quality changed to less diagenetically altered material and coincided with an increase in the ribulose-1,5-bisphosphate carboxylase/oxygenase form I (cbbL) gene number followed by an influx of labile DOM. Upon reoxygenation, the deep DOM that accumulated under suboxic conditions is bioavailable to surface prokaryotes that utilize the accumulated DOC initially before switching to a community that can both produce labile DOM via chemoautotrophy and degrade the more recalcitrant DOM.

Microbial communities from weathered outcrops of a sulfide-rich ultramafic intrusion, and implications for mine waste management

The Duluth Complex (DC) contains sulfide-rich magmatic intrusions that represent one of the largest known economic deposits of copper, nickel, and platinum group elements. Previous work showed that microbial communities associated with experimentally-weathered DC waste rock and tailings were dominated by uncultivated taxa and organisms not typically associated with mine waste. However, those experiments were designed for kinetic testing and do not necessarily represent the conditions expected for long-term environmental weathering. We used 16S rRNA gene methods to characterize the microbial communities present on the surfaces of naturally-weathered and historically disturbed outcrops of DC material. Rock surfaces were dominated by diverse uncultured Ktedonobacteria, Acetobacteria, and Actinobacteria, with abundant algae and other phototrophs. These communities were distinct from microbial assemblages from experimentally-weathered DC rocks, suggesting different energy and nutrient resources in environmental samples. Sulfide mineral incubations performed with and without algae showed that photosynthetic microorganisms could have an inhibitory effect on autotrophic populations, resulting in slightly lower sulfate release and differences in dominant microorganisms. The microbial assemblages from these weathered outcrops show how communities develop during weathering of sulfide-rich DC rocks and represent baseline data that could evaluate the effectiveness of future reclamation of waste produced by large-scale mining operations.

Primer selection impacts the evaluation of microecological patterns in environmental microbiomes

This study revealed that primer selection substantially influences the taxonomic and predicted functional composition and the characterization of microecological patterns, which was not alleviated by close-reference clustering. Biases were relatively consistent across different habitats in community profiling but not in microecological patterns. These primer biases could be attributed to multiple aspects, including taxa specificity, regional hypervariability, and amplification efficiency.

Initial Post-Release Performance of Cultured Cyprinus chilia Juveniles in a Shallow Lake in Southwestern China

The post-release performance of cultured fish is crucial for understanding the viability of cultured fish and assessing the effects of stock enhancement programs. This study aimed to investigate the initial post-release performance of cultured Cyprinus chilia juveniles by examining their movement, spatial distribution, gut fullness, and gut microbiota in nature. In July 2022, a total of 20,000 C. chilia juveniles, tagged with visible implant fluorescence (VIE), were released into Qilu Lake, a shallow lake in southwestern China. Subsequently, continuous recapture was conducted at fixed recapture sites using trap nets during the first 7 days, one month and three months after release. Out of the released fish, 512 were recaptured, resulting in a recapture rate of 2.56%. The recaptured fish had a 100% tag retention rate. The majority (98.05%) of the recaptured fish were found in the recapture sites located on the eastern or western lakeshore, while only 10 fish were recaptured from the recapture sites in the northern lake area. The water depth range where the recaptured fish were found ranged from 190 to 350 cm, with most fish preferring depths less than 300 cm. The majority of the released fish migrated towards the eastern and western lakeshore, with long-distance movement (greater than 100 m) primarily occurring within the first four days after release. The level of gut fullness in the released fish initially decreased and then increased over time following release. Regarding gut microbiota, the dominant phyla observed in most samples were Firmicutes, Proteobacteria, Cyanobacteria, and Fusobacteria. Furthermore, significant variations in the dominant genera were observed across different samples. Principal coordinates analysis (PCoA) revealed clear separation between the microbial communities of pre-release and post-release C. chilia juveniles. This study demonstrated that VIE tagging was a suitable method for short-term marking of C. chilia juveniles. Lakeshores with water depths less than 300 cm were identified as preferred habitats for C. chilia juveniles. The primary adaptation period for cultured C. chilia juveniles released into the natural environment was found to be approximately 4-5 days. These findings contribute to our understanding of the post-release performance of cultured fish and may provide guidance for the management and evaluation of relevant stock enhancement programs.

Identification of rare microbial colonizers of plastic materials incubated in a coral reef environment

Plastic waste accumulation in marine environments has complex, unintended impacts on ecology that cross levels of community organization. To measure succession in polyolefin-colonizing marine bacterial communities, an in situ time-series experiment was conducted in the oligotrophic coastal waters of the Bermuda Platform. Our goals were to identify polyolefin colonizing taxa and isolate bacterial cultures for future studies of the biochemistry of microbe-plastic interactions. HDPE, LDPE, PP, and glass coupons were incubated in surface seawater for 11 weeks and sampled at two-week intervals. 16S rDNA sequencing and ATR-FTIR/HIM were used to assess biofilm community structure and chemical changes in polymer surfaces. The dominant colonizing taxa were previously reported cosmopolitan colonizers of surfaces in marine environments, which were highly similar among the different plastic types. However, significant differences in rare community composition were observed between plastic types, potentially indicating specific interactions based on surface chemistry. Unexpectedly, a major transition in community composition occurred in all material treatments between days 42 and 56 (p < 0.01). Before the transition, Alteromonadaceae, Marinomonadaceae, Saccharospirillaceae, Vibrionaceae, Thalassospiraceae, and Flavobacteriaceae were the dominant colonizers. Following the transition, the relative abundance of these taxa declined, while Hyphomonadaceae, Rhodobacteraceae and Saprospiraceae increased. Over the course of the incubation, 8,641 colonizing taxa were observed, of which 25 were significantly enriched on specific polyolefins. Seven enriched taxa from families known to include hydrocarbon degraders (Hyphomonadaceae, Parvularculaceae and Rhodobacteraceae) and one n-alkane degrader (Ketobacter sp.). The ASVs that exhibited associations with specific polyolefins are targets of ongoing investigations aimed at retrieving plastic-degrading microbes in culture.

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

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

Evaluation of different 16S rRNA gene hypervariable regions and reference databases for profiling engineered microbiota structure and functional guilds in a swine wastewater treatment plant

High-throughput 16S rRNA gene amplicon sequencing technology is widely applied for environmental microbiota structure analysis to derive knowledge that informs microbiome-based surveillance and oriented bioengineering. However, it remains elusive how the selection of 16S rRNA gene hypervariable regions and reference databases affects microbiota diversity and structure profiling. This study systematically evaluated the fitness of different frequently used reference databases (i.e. SILVA 138 SSU, GTDB bact120_r207, Greengenes 13_5 and MiDAS 4.8) and primers of 16S rRNA gene in microbiota profiling of anaerobic digestion and activated sludge collected from a full-scale swine wastewater treatment plant (WWTP). The comparative results showed that MiDAS 4.8 achieved the highest levels of taxonomic diversity and species-level assignment rate. For whichever sample groups, microbiota richness captured by different primers decreased in the following order: V4 > V4-V5 > V3-V4 > V6-V8/V1-V3. Using primer-bias-free metagenomic data results as the judging standard, V4 region also best characterized microbiota structure and well represented typical functional guilds (e.g. methanogens, ammonium oxidizers and denitrifiers), while V6-V8 regions largely overestimated the archaeal methanogens (mainly Methanosarcina) by over 30 times. Therefore, MiDAS 4.8 database and V4 region are recommended for best simultaneous analysis of bacterial and archaeal community diversity and structure of the examined swine WWTP.

Reduced Arbuscular Mycorrhizal Fungi (AMF) Diversity in Light and Moderate Fire Sites in Taiga Forests, Northeast China

Forest fires are an important disturbance factor in forest ecosystems, and obviously change the soil environment. Arbuscular mycorrhizal fungi, as a medium and bridge between vegetation and soil, play a crucial role in mediating plant nutrient uptake and regulating the productivity, stability, and succession of vegetation-soil systems. To investigate the effects of forest fires on the community structure and diversity of arbuscular mycorrhizal fungi in cold-temperate Larix gmelinii forests, we collected soils from light, moderate, and heavy fire disturbance forests and a natural forest as a control forest in Greater Khingan Larix gmelinii forests, in the northeast of China. The community structure and diversity of arbuscular mycorrhizal fungi was sequenced using Illumina MiSeq technology and we analyzed the correlation with the soil physicochemical characteristics. The results showed that the contents of microbial biomass content (MBC), moisture content (MC), total nitrogen (TN), and available phosphors (AP) increased significantly (p < 0.05) with increasing fire intensity (from Light to heavy fire), but available potassium (AK) decreased significantly (p < 0.05). These changes were not significant. A total of 14,554 valid sequences from all sequences were classified into 66 ASVs that belonged into one phylum, one order, four families, and four genera. The genera included Glomus, Ambispora, Paraglomus, and Acaulospora, and Glomus was the dominant genus (the genera with the five most relative abundances) in the control and heavy-fire forests. Non-metric multidimensional scaling (NMDS) analysis showed that forest fires significantly affected the community structure of arbuscular mycorrhizal fungi (p < 0.01). Redundancy analysis (RDA) showed that MBC, SOC, and AP contents significantly affected the composition structure and diversity of arbuscular mycorrhizal fungi communities. This study indicated that forest fires affected the composition and diversity of soil arbuscular mycorrhizal fungi communities through changing the soil physicochemical parameters (MBC, SOC, and AP) in cold-temperate Larix gmelinii forests. The study of soil physicochemical properties and arbuscular mycorrhizal fungi diversity in cold-temperate Larix gmelinii forests in the Greater Khingan Mountains after forest fires provides a reference basis for the revegetation and reconstruction of fire sites.

Heavy Metal Contamination Alters the Co-Decomposition of Leaves of the Invasive Tree Rhus typhina L. and the Native Tree Koelreuteria paniculata Laxm

Invasive and native plants can coexist in the same habitat; however, the decomposition process may be altered by the mixing of invasive and native leaves. Heavy metal contamination may further alter the co-decomposition of both leaf types. This study evaluated the effects of two concentrations (35 mg·L^-1 and 70 mg·L^-1) and three types (Pb, Cu, and combined Pb + Cu) of heavy metal contamination on the co-decomposition of leaves of the invasive tree Rhus typhina L. and the native tree Koelreuteria paniculata Laxm, as well as the mixed effect intensity of the co-decomposition of the mixed leaves. A polyethylene litterbag experiment was performed over six months. The decomposition coefficient of the two trees, mixed effect intensity of the co-decomposition, soil pH and enzymatic activities, soil bacterial alpha diversity, and soil bacterial community structure were determined. A high concentration of Pb and combined Pb + Cu significantly reduced the decomposition rate of R. typhina leaves. A high concentration of Pb or Cu significantly reduced the decomposition rate of the mixed leaves. In general, R. typhina leaves decomposed faster than K. paniculata leaves did. There were synergistic effects observed for the co-decomposition of the mixed leaves treated with combined Pb + Cu, regardless of concentration, but there were antagonistic effects observed for the co-decomposition of the mixed leaves treated with either Pb or Cu, regardless of concentration. A high concentration of Pb or Cu may increase antagonistic effects regarding the co-decomposition of mixed-leaf groups. Thus, heavy metal contamination can significantly affect the intensity of the mixed effect on the co-decomposition of heterogeneous groups of leaves.

Selecting 16S rRNA Primers for Microbiome Analysis in a Host-Microbe System: The Case of the Jellyfish Rhopilema nomadica

Amplicon sequencing of the 16S rRNA gene is extensively used to characterize bacterial communities, including those living in association with eukaryotic hosts. Deciding which region of the 16S rRNA gene to analyze and selecting the appropriate PCR primers remains a major decision when initiating any new microbiome study. Based on a detailed literature survey of studies focusing on cnidarian microbiomes, we compared three commonly used primers targeting different hypervariable regions of the 16S rRNA gene, V1V2, V3V4, and V4V5, using the jellyfish Rhopilema nomadica as a model. Although all primers exhibit a similar pattern in bacterial community composition, the performance of the V3V4 primer set was superior to V1V2 and V4V5. The V1V2 primers misclassified bacteria from the Bacilli class and exhibited low classification resolution for Rickettsiales, which represent the second most abundant 16S rRNA gene sequence in all the primers. The V4V5 primer set detected almost the same community composition as the V3V4, but the ability of these primers to also amplify the eukaryotic 18S rRNA gene may hinder bacterial community observations. However, after overcoming the challenges possessed by each one of those primers, we found that all three of them show very similar bacterial community dynamics and compositions. Nevertheless, based on our results, we propose that the V3V4 primer set is potentially the most suitable for studying jellyfish-associated bacterial communities. Our results suggest that, at least for jellyfish samples, it may be feasible to directly compare microbial community estimates from different studies, each using different primers but otherwise similar experimental protocols. More generally, we recommend specifically testing different primers for each new organism or system as a prelude to large-scale 16S rRNA gene amplicon analyses, especially of previously unstudied host-microbe associations.

The Preventive Effects of Probiotic Prevotella histicola on the Bone Loss of Mice with Ovariectomy-Mediated Osteoporosis

It has been demonstrated that the disturbance of gut microbiota (GM) is closely related to the reduction of bone mass and incidence of osteoporosis (OP). The aim of this study is to investigate whether the supplementation of Prevotella histicola (Ph) can prevent the bone loss in mice with ovariectomy (OVX)-mediated OP, and further explore relevant mechanisms. Regular (once a day for 8 consecutive weeks) and quantitative (200 µL/d) perfusion of Ph (the bacteria that orally gavaged) was conducted starting from 1 week after the construction of mice models. Bone mass and bone microstructure were detected by Micro-computed tomography (Micro-CT). Expressions of intestinal permeability, pro-inflammatory cytokines, and osteogenic and osteoclastic activities of mice were analyzed by histological staining and immunohistochemistry (IHC). 16S rRNA high throughput sequencing technique was applied to analyze the alterations of composition, abundance, and diversity of collected feces. Regular and quantitative perfusion of Ph mitigated the bone loss in mice with OVX-mediated OP. Compared with OVX + PBS group, perfusion of Ph repressed osteoclastogenesis and promoted osteogenesis, reduced release of pro-inflammatory cytokine cytokines (interleukin-1β (IL-1β) and tumor necrosis factor-α (TNF-α)), and reversed expressions of tight junction proteins (zonula occludens protein 1 (ZO-1) and Occludin). Besides, the perfusion of Ph improved the composition, abundance, and diversity of GM. Collectively, this study revealed that regular and quantitative perfusion of Ph can improve the bone loss in mice with OVX-mediated OP by repairing intestinal mucosal barrier damage, optimizing intestinal permeability, inhibiting release of pro-osteoclastogenic cytokines, and improving disturbance of GM.

Seasonal variations of soil bacterial and fungal communities in a subtropical Eucalyptus plantation and their responses to throughfall reduction

Climatic change causes obvious seasonal meteorological drought in southern China, yet there is a lack of comprehensive in situ studies on the effects of drought in Eucalyptus plantations. Here, a 50% throughfall reduction (TR) experiment was conducted to investigate the seasonal variations of soil bacterial and fungal communities and functions in a subtropical Eucalyptus plantation and their responses to TR treatment. Soil samples were collected from control (CK) and TR plots in the dry and rainy seasons and were subjected to high-throughput sequencing analysis. Results showed that TR treatment significantly reduced soil water content (SWC) in the rainy season. In CK and TR treatments, fungal alpha-diversity decreased in the rainy season while bacterial alpha-diversity did not change significantly between dry and rainy seasons. Moreover, bacterial networks were more affected by seasonal variations compared with fungal networks. Redundancy analysis showed that alkali hydrolyzed nitrogen and SWC contributed the most to the bacterial and fungal communities, respectively. Functional prediction indicated that the expression of soil bacterial metabolic functions and symbiotic fungi decreased in the rainy season. In conclusion, seasonal variations have a stronger effect on soil microbial community composition, diversity, and function compared with TR treatment. These findings could be used to develop management practices for subtropical Eucalyptus plantations and help maintain soil microbial diversity to sustain long-term ecosystem function and services in response to future changes in precipitation patterns.

Temporal dynamics of geothermal microbial communities in Aotearoa-New Zealand

Microbial biogeography studies, in particular for geothermal-associated habitats, have focused on spatial patterns and/or individual sites, which have limited ability to describe the dynamics of ecosystem behaviour. Here, we report the first comprehensive temporal study of bacterial and archaeal communities from an extensive range of geothermal features in Aotearoa-New Zealand. One hundred and fifteen water column samples from 31 geothermal ecosystems were taken over a 34-month period to ascertain microbial community stability (control sites), community response to both natural and anthropogenic disturbances in the local environment (disturbed sites) and temporal variation in spring diversity across different pH values (pH 3, 5, 7, 9) all at a similar temperature of 60–70°C (pH sites). Identical methodologies were employed to measure microbial diversity via 16S rRNA gene amplicon sequencing, along with 44 physicochemical parameters from each feature, to ensure confidence in comparing samples across timeframes. Our results indicated temperature and associated groundwater physicochemistry were the most likely parameters to vary stochastically in these geothermal features, with community abundances rather than composition more readily affected by a changing environment. However, variation in pH (pH ±1) had a more significant effect on community structure than temperature (±20°C), with alpha diversity failing to adequately measure temporal microbial disparity in geothermal features outside of circumneutral conditions. While a substantial physicochemical disturbance was required to shift community structures at the phylum level, geothermal ecosystems were resilient at this broad taxonomic rank and returned to a pre-disturbed state if environmental conditions re-established. These findings highlight the diverse controls between different microbial communities within the same habitat-type, expanding our understanding of temporal dynamics in extreme ecosystems.

Microbial Interactions with Dissolved Organic Matter Are Central to Coral Reef Ecosystem Function and Resilience

To thrive in nutrient-poor waters, coral reefs must retain and recycle materials efficiently. This review centers microbial processes in facilitating the persistence and stability of coral reefs, specifically the role of these processes in transforming and recycling the dissolved organic matter (DOM) that acts as an invisible currency in reef production, nutrient exchange, and organismal interactions. The defining characteristics of coral reefs, including high productivity, balanced metabolism, high biodiversity, nutrient retention, and structural complexity, are inextricably linked to microbial processing of DOM. The composition of microbes and DOM in reefs is summarized, and the spatial and temporal dynamics of biogeochemical processes carried out by microorganisms in diverse reef habitats are explored in a variety of key reef processes, including decomposition, accretion, trophictransfer, and macronutrient recycling. Finally, we examine how widespread habitat degradation of reefs is altering these important microbe-DOM interactions, creating feedbacks that reduce reef resilience to global change.

Microbial communities on eelgrass (Zostera marina) thriving in Tokyo Bay and the possible source of leaf-attached microbes

Zostera marina (eelgrass) is classified as one of the marine angiosperms and is widely distributed throughout much of the Northern Hemisphere. The present study investigated the microbial community structure and diversity of Z. marina growing in Futtsu bathing water, Chiba prefecture, Japan. The purpose of this study was to provide new insight into the colonization of eelgrass leaves by microbial communities based on leaf age and to compare these communities to the root-rhizome of Z. marina, and the surrounding microenvironments (suspended particles, seawater, and sediment). The microbial composition of each sample was analyzed using 16S ribosomal gene amplicon sequencing. Each sample type was found to have a unique microbial community structure. Leaf-attached microbes changed in their composition depending on the relative age of the eelgrass leaf. Special attention was given to a potential microbial source of leaf-attached microbes. Microbial communities of marine particles looked more like those of eelgrass leaves than those of water samples. This finding suggests that leaf-attached microbes were derived from suspended particles, which could allow them to go back and forth between eelgrass leaves and the water column.

Particle-attached Microbes in Eelgrass Vegetation Areas Differ in Community Structure Depending on the Distance from the Eelgrass Bed

Zostera marina (eelgrass) is a submerged flowering plant often found in the coastal areas of Japan. Large amounts of suspended particles form in highly productive environments, such as eelgrass beds, and the behavior of these particles is expected to affect the surrounding microbial community. We investigated the microbial community structure of suspended particles in three eelgrass fields (Ikuno-Shima Is., Mutsu Bay, and Nanao Bay) and inferred the formation and dynamics of suspended particles from a microbial community structure ana-lysis. Seawater samples were collected directly above each eelgrass bed (eelgrass-covering) and from locations dozens of meters away from the eelgrass bed (bare-ground). In consideration of the two different lifestyles of marine microbes, microbial communities were obtained from particle-attached (PA) and free-living (FL) states. Differences in microbial diversity and community structures were observed between PA and FL in all eelgrass beds. The FL microbial community was similar between the two sampling points (eelgrass-covering and bare-ground), whereas a significant difference was noted in the microbial community structure of suspended particles between the two sampling points. This difference appeared to be due to the supply of organic matter from the eelgrass sea ground and leaf-attached detritus produced by microbial activity. In addition, the classes Flavobacteriia, Alphaproteobacteria, and Gammaproteobacteria were abundant in the PA and FL fractions. Furthermore, many sequences of the key groups (e.g., Planctomycetes and Verrucomicrobia) were exclusively detected in the PA fraction, in which they may circulate nutrients. The present results provide insights into the microbial communities of suspended particles and provide the first step towards understanding their biogeochemical impact on the eelgrass bed.

Fine-scale evaluation of two standard 16S rRNA gene amplicon primer pairs for analysis of total prokaryotes and archaeal nitrifiers in differently managed soils

The advance of high-throughput molecular biology tools allows in-depth profiling of microbial communities in soils, which possess a high diversity of prokaryotic microorganisms. Amplicon-based sequencing of 16S rRNA genes is the most common approach to studying the richness and composition of soil prokaryotes. To reliably detect different taxonomic lineages of microorganisms in a single soil sample, an adequate pipeline including DNA isolation, primer selection, PCR amplification, library preparation, DNA sequencing, and bioinformatic post-processing is required. Besides DNA sequencing quality and depth, the selection of PCR primers and PCR amplification reactions arguably have the largest influence on the results. This study tested the performance and potential bias of two primer pairs, i.e., 515F (Parada)-806R (Apprill) and 515F (Parada)-926R (Quince) in the standard pipelines of 16S rRNA gene Illumina amplicon sequencing protocol developed by the Earth Microbiome Project (EMP), against shotgun metagenome-based 16S rRNA gene reads. The evaluation was conducted using five differently managed soils. We observed a higher richness of soil total prokaryotes by using reverse primer 806R compared to 926R, contradicting to in silico evaluation results. Both primer pairs revealed various degrees of taxon-specific bias compared to metagenome-derived 16S rRNA gene reads. Nonetheless, we found consistent patterns of microbial community variation associated with different land uses, irrespective of primers used. Total microbial communities, as well as ammonia oxidizing archaea (AOA), the predominant ammonia oxidizers in these soils, shifted along with increased soil pH due to agricultural management. In the unmanaged low pH plot abundance of AOA was dominated by the acid-tolerant NS-Gamma clade, whereas limed agricultural plots were dominated by neutral-alkaliphilic NS-Delta/NS-Alpha clades. This study stresses how primer selection influences community composition and highlights the importance of primer selection for comparative and integrative studies, and that conclusions must be drawn with caution if data from different sequencing pipelines are to be compared.

Adeno-associated virus vector intraperitoneal injection induces colonic mucosa and submucosa transduction and alters the diversity and composition of the faecal microbiota in rats

Background

Viral vector technology, especially recombinant adeno-associated virus vector (rAAV) technology, has shown great promise in preclinical research for clinical applications. Several studies have confirmed that rAAV can successfully transduce the enteric nervous system (ENS), and rAAV gene therapy has been approved by the Food and Drug Administration (FDA) for the treatment of the early childhood blindness disease Leber congenital amaurosis and spinal muscular atrophy (SMA). However, until now, it has not been possible to determine the effect of AAV9 on intestinal microbiota.

Methods

We examined the efficiency of AAV9-mediated ascending colon, transverse colon and descending colon transduction through intraperitoneal (IP) injection, performed 16S rRNA gene amplicon sequencing and analysed specific faecal microbial signatures following AAV9 IP injection via bioinformatics methods in Sprague-Dawley (SD) rats.

Results

Our results showed (1) efficient transduction of the mucosa and submucosa of the ascending, transverse, and descending colon following AAV9 IP injection; (2) a decreased alpha diversity and an altered overall microbial composition following AAV9 IP injection; (3) significant enrichments in a total of 5 phyla, 10 classes, 13 orders, 15 families, 29 genera, and 230 OTUs following AAV9 IP injection; and (4) AAV9 can significantly upregulate the relative abundance of anaerobic microbiota which is one of the seven high-level phenotypes that BugBase could predict.

Conclusion

In summary, these data show that IP injection of AAV9 can successfully induce the transduction of the colonic mucosa and submucosa and alter the diversity and composition of the faecal microbiota in rats.

Diversity and putative metabolic function of prokaryotic communities in tank bromeliads along an elevation gradient in tropical Mexico

Tank bromeliads are unique canopy microhabitats that offer freshwater and organic nutrient-rich substrates in the Neotropics. In them it is possible to thoroughly characterize environmental factors and species composition of terrestrial and aquatic biota. Therefore, these plants have been used as natural models to study how communities are distributed and assembled. Here we used amplicon sequencing of the 16S rRNA gene and their functional annotations to study the diversity and metabolic potential of prokaryotic communities in tank bromeliads in five different forests along an elevation gradient in tropical Mexico. Furthermore, we analyzed the effects of vegetation type and environmental factors inside the tanks on prokaryotic composition. We found a high prokaryotic diversity in tank bromeliads along the elevation gradient. Prokaryotes commonly observed in acidic environments rich in organic carbon, and the potential pathogen Pasteurella multocida, were present in all samples, but few amplicon sequence variants were shared between forests. The prokaryotic composition was affected by forest type, and comparisons against null models suggest that it was shaped by non-neutral processes. Furthermore, prokaryotic community changes significantly covaried with tank water temperature, pH, and inorganic carbon. We found a high diversity of putative metabolic groups dominated by chemoheterotrophs and fermenters, but taxonomic groups involved in nitrogen and sulfur cycling were also present in all samples. These results suggest that tank bromeliads promote taxonomic and metabolic diversity of the prokaryotic community at a local and regional scale and play an important role in the biogeochemistry of forest canopies in the Neotropics.

Effects of green tide on microbial communities in waters of the Jiangsu coastal area, China

Recently, green tide outbreaks have resulted in severe coastal ecology and economic effects in China. Jiangsu coastal areas are usually the site of early green tide outbreaks. To clarify the effects of green tide outbreaks in Jiangsu coastal areas, this study analyzed microbial communities during green tide-free and green tide outbreak periods (May and July, respectively) through 16S rDNA sequencing. Sequences were clustered into 4117 operational taxonomic units (OTUs), 1044 and 3834 of which were obtained from the May and July groups, respectively. Redundancy analysis indicated that green tide occurrence was closely associated with the temperature, pH, and concentrations of various nutrients. Diversity analysis revealed that the July group had a richer microbial community than the May group, in agreement with the results of propagule culture. Moreover, comparative analysis revealed that samples in the May and July groups clustered together. According to Megan analysis, the May group had much more Psychrobacter, Sulfitobacter, and Marinomonas than the July group, whereas the other genera were predominantly found in July, such as Ascidiacerhabitans, Synechococcus Hydrotalea, and Burkholderia-Paraburkholderia. These findings suggest that green tide outbreaks affect marine microbial communities, and detecting the changes in the identified genera during green tide outbreaks may contribute to green tide forecasting. PRACTITIONER POINTS: Jiangsu coastal areas are usually the site of early green tide outbreaks. Green tide occurrence was related to the concentrations of various nutrients. Microbial species and community structure significantly changed after green tide outbreak.

Interaction and spatio-taxonomic patterns of the soil microbiome around oil production wells impacted by petroleum hydrocarbons

Numerous onshore oil production wells currently exist, and the petroleum hydrocarbon contamination of the surrounding soil caused by oil production wells is not well understood. Moreover, the impact of the distribution of the total petroleum hydrocarbons (TPH) in the soil on the microbiota requires further investigation. Accordingly, in this study, the distribution of petroleum hydrocarbons in the soils around oil production wells was investigated, and their alteration of the microbiota was revealed. The results revealed that in the horizontal direction, the heavily TPH-contaminated soils were mainly distributed within a circle with a radius of 200 cm centered on the oil production well; and in the vertical direction, the heavily TPH-contaminated soils were distributed within the 0-50 cm soil layer. A significant positive correlation was found between the microbial abundance and the TPH concentration in the soil with relatively low total carbon contents. Heavy TPH contamination (TPH concentration of >3000 mg/kg) significantly reduced the microbial diversity and altered the microbiota compared with the light TPH contamination (TPH concentration of around 1000 mg/kg). In the heavily TPH-contaminated soils, the relative abundances of the Proteobacteria and Bacteroides increased significantly; the network complexity among the soil microorganisms decreased; and the co-occurrence patterns were altered. In summary, the results of this study have reference value in the remediation of soils around oil production wells and provide guidance for the construction of microbial remediation systems for petroleum contamination.

Dynamic Changes in the Gut Microbial Community and Function during Broiler Growth

During the entire growth process, gut microbiota continues to change and has a certain impact on the performance of broilers. Here, we used 16S rRNA gene sequencing to explore the dynamic changes in the fecal bacterial communities and functions in 120 broilers from 4 to 16 weeks of age. We found that the main phyla (Firmicutes, Fusobacteria, Proteobacteria, and Bacteroides) accounted for more than 93.5% of the total bacteria in the feces. The alpha diversity of the fecal microbiota showed a downward trend with time, and the beta diversity showed significant differences at various time points. Then, the study on the differences of microbiota between high-weight (HW) and low-weight (LW) broilers showed that there were differences in the diversity and composition of microbiota between high- and low-weight broilers. Furthermore, we identified 22 genera that may be related to the weight change of broilers. The analysis of flora function reveals their changes in metabolism, genetic information processing, and environmental information processing. Finally, combined with microbial function and cecal transcriptome results, we speculated that microorganisms may affect the immune level and energy metabolism level of broilers through their own carbohydrate metabolism and lipid metabolism and then affect body weight (BW). Our results will help to expand our understanding of intestinal microbiota and provide guidance for the production of high-quality broilers.

IMPORTANCE The intestinal microbiota has a certain impact on the performance of broilers. However, the change of intestinal microbiota after 4 weeks of age is not clear, and the mechanism of the effect of microorganisms on the weight change of broilers needs more exploration. After 4 weeks of age, the alpha diversity of microorganisms in broiler feces decreased, and the dominant bacteria were Firmicutes, Fusobacteria, Proteobacteria, and Bacteroides. There were differences in microbiota diversity and composition between high- and low-weight broilers. Intestinal microorganisms may affect the immune level and energy metabolism level of broilers through their own carbohydrate metabolism and lipid metabolism and then affect the body weight. The results are helpful to increase the understanding of intestinal microbiota and provide reference for the production of high-quality broilers.

A robust approach to estimate relative phytoplankton cell abundances from metagenomes

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.

Nutrient supply controls the linkage between species abundance and ecological interactions in marine bacterial communities

Nutrient scarcity is pervasive for natural microbial communities, affecting species reproduction and co-existence. However, it remains unclear whether there are general rules of how microbial species abundances are shaped by biotic and abiotic factors. Here we show that the ribosomal RNA gene operon (rrn) copy number, a genomic trait related to bacterial growth rate and nutrient demand, decreases from the abundant to the rare biosphere in the nutrient-rich coastal sediment but exhibits the opposite pattern in the nutrient-scarce pelagic zone of the global ocean. Both patterns are underlain by positive correlations between community-level rrn copy number and nutrients. Furthermore, inter-species co-exclusion inferred by negative network associations is observed more in coastal sediment than in ocean water samples. Nutrient manipulation experiments yield effects of nutrient availability on rrn copy numbers and network associations that are consistent with our field observations. Based on these results, we propose a “hunger games” hypothesis to define microbial species abundance rules using the rrn copy number, ecological interaction, and nutrient availability.

Environmental and biotic factors control ecological communities. Here, the authors study community ribosomal rRNA gene copy number in coastal sediment and ocean bacterial communities, and in microcosm nutrient addition experiments, to propose a conceptual framework of how nutrient supply and ecological interactions shape the community.

Soil Bacterial and Fungal Community Responses to Throughfall Reduction in a Eucalyptus Plantation in Southern China

In subtropical plantations in southern China, how soil microbial communities respond to climate change-induced drought is poorly understood. A field experiment was conducted in a subtropical Eucalyptus plantation to determine the impacts of 50% of throughfall reduction (TR) on soil microbial community composition, function, and soil physicochemical properties. Results showed that TR reduced soil water content (SWC) and soil available phosphorus (AP) content. TR significantly altered 196 bacterial operational taxonomic units (OTUs), most of them belonging to Acidobacteria, Actinobacteria, and Proteobacteria, while there were fewer changes in fungal OTUs. At the phylum level, TR increased the relative abundance of Acidobacteria at 0–20 cm soil depth by 37.18%, but failed to influence the relative abundance of the fungal phylum. Notably, TR did not alter the alpha diversity of the bacterial and fungal communities. The redundancy analysis showed that the bacterial communities were significantly correlated with SWC, and fungal communities were significantly correlated with AP content. According to predictions of bacterial and fungal community functions using PICRUSt2 and FUNGuild platforms, TR had different effects on both bacterial and fungal communities. Overall, SWC and AP decreased during TR, resulting in greater changes in soil bacterial community structure, but did not dramatically change soil fungal community structure.

Long-term and combined effects of N-[2-(nitrooxy)ethyl]-3-pyridinecarboxamide and fumaric acid on methane production, rumen fermentation, and lactation performance in dairy goats

BACKGROUND

In recent years, nitrooxy compounds have been identified as promising inhibitors of methanogenesis in ruminants. However, when animals receive a nitrooxy compound, a high portion of the spared hydrogen is eructated as gas, which partly offsets the energy savings of CH₄ mitigation. The objective of the present study was to evaluate the long-term and combined effects of supplementation with N-[2-(nitrooxy)ethyl]-3-pyridinecarboxamide (NPD), a methanogenesis inhibitor, and fumaric acid (FUM), a hydrogen sink, on enteric CH₄ production, rumen fermentation, bacterial populations, apparent nutrient digestibility, and lactation performance of dairy goats.

RESULTS

Twenty-four primiparous dairy goats were used in a randomized complete block design with a 2 × 2 factorial arrangement of treatments: supplementation without or with FUM (32 g/d) or NPD (0.5 g/d). All samples were collected every 3 weeks during a 12-week feeding experiment. Both FUM and NPD supplementation persistently inhibited CH₄ yield (L/kg DMI, by 18.8% and 18.1%, respectively) without negative influence on DMI or apparent nutrient digestibility. When supplemented in combination, no additive CH₄ suppression was observed. FUM showed greater responses in increasing the molar proportion of propionate when supplemented with NPD than supplemented alone (by 10.2% vs. 4.4%). The rumen microbiota structure in the animals receiving FUM was different from that of the other animals, particularly changed the structure of phylum Firmicutes. Daily milk production and serum total antioxidant capacity were improved by NPD, but the contents of milk fat and protein were decreased, probably due to the bioactivity of absorbed NPD on body metabolism.

CONCLUSIONS

Supplementing NPD and FUM in combination is a promising way to persistently inhibit CH₄ emissions with a higher rumen propionate proportion. However, the side effects of this nitrooxy compound on animals and its residues in animal products need further evaluation before it can be used as an animal feed additive.

Microbial diversity in tropical marine sediments assessed using culture-dependent and culture-independent techniques

The microbial communities associated with marine sediments are critical for ecosystem function yet remain poorly characterized. While culture-independent (CI) techniques capture the broadest perspective on community composition, culture-dependent (CD) methods can select for low abundance taxa that are missed using CI approaches. This study aimed to assess microbial diversity in tropical marine sediments at five shallow-water sites in Belize using both CD and CI techniques. The CD methods captured approximately 3% of the >800 genera detected across all sites using the CI approach. Additionally, 39 genera were only detected in culture, revealing rare taxa that were missed with the CI approach. Significantly different communities were detected across sites, with rare taxa playing an important role in distinguishing among communities. This study provides important baseline data describing shallow-water sediment microbial communities, evidence that standard cultivation techniques may be more effective than previously recognized, and the first steps towards identifying new taxa that are amenable to agar plate cultivation.

Evaluating and Improving Small Subunit rRNA PCR Primer Coverage for Bacteria, Archaea, and Eukaryotes Using Metagenomes from Global Ocean Surveys

Small subunit rRNA (SSU rRNA) amplicon sequencing can quantitatively and comprehensively profile natural microbiomes, representing a critically important tool for studying diverse global ecosystems. However, results will only be accurate if PCR primers perfectly match the rRNA of all organisms present. To evaluate how well marine microorganisms across all 3 domains are detected by this method, we compared commonly used primers with >300 million rRNA gene sequences retrieved from globally distributed marine metagenomes. The best-performing primers compared to 16S rRNA of bacteria and archaea were 515Y/926R and 515Y/806RB, which perfectly matched over 96% of all sequences. Considering cyanobacterial and chloroplast 16S rRNA, 515Y/926R had the highest coverage (99%), making this set ideal for quantifying marine primary producers. For eukaryotic 18S rRNA sequences, 515Y/926R also performed best (88%), followed by V4R/V4RB (18S rRNA specific; 82%)-demonstrating that the 515Y/926R combination performs best overall for all 3 domains. Using Atlantic and Pacific Ocean samples, we demonstrate high correspondence between 515Y/926R amplicon abundances (generated for this study) and metagenomic 16S rRNA (median R² = 0.98, n = 272), indicating amplicons can produce equally accurate community composition data compared with shotgun metagenomics. Our analysis also revealed that expected performance of all primer sets could be improved with minor modifications, pointing toward a nearly completely universal primer set that could accurately quantify biogeochemically important taxa in ecosystems ranging from the deep sea to the surface. In addition, our reproducible bioinformatic workflow can guide microbiome researchers studying different ecosystems or human health to similarly improve existing primers and generate more accurate quantitative amplicon data. IMPORTANCE PCR amplification and sequencing of marker genes is a low-cost technique for monitoring prokaryotic and eukaryotic microbial communities across space and time but will work optimally only if environmental organisms match PCR primer sequences exactly. In this study, we evaluated how well primers match globally distributed short-read oceanic metagenomes. Our results demonstrate that primer sets vary widely in performance, and that at least for marine systems, rRNA amplicon data from some primers lack significant biases compared to metagenomes. We also show that it is theoretically possible to create a nearly universal primer set for diverse saline environments by defining a specific mixture of a few dozen oligonucleotides, and present a software pipeline that can guide rational design of primers for any environment with available meta'omic data.

The Impact of Intensive Fish Farming on Pond Sediment Microbiome and Antibiotic Resistance Gene Composition

Aquaculture is a fast-growing animal food sector, and freshwater fish farming is particularly common in Central and Eastern Europe. As the biodiversity of fishery ponds is changed toward fulfilling the industrial needs, precautions should be taken to keep the system sustainable and protect the adjacent environment from possible damage. Due to risk of infectious diseases, antibiotics are used in aquaculture production systems. The constant exposure to antimicrobials can contribute to the rise of antibiotic resistance in aquaculture products and the adjacent ecosystems, with possibility of dissemination to the wider environment as well as between animals and humans. Even though previous studies have found antibiotic resistance genes in the sediments and water of farming ponds, the tendency and direction of spreading is not clear yet. The objective of this project was to evaluate the influence of intensive fish farming on the condition of water bodies used for the aquaculture and the environment, concentrating on the impact of the aquaculture on the surrounding water ecosystems as well as the possibility of transferring the pollutants and antibiotic resistance genes to both environment and the human hosts. Combined measurement of antibiotic and heavy metal contamination, toxicity assessment, microorganism diversity, and the detection of common antibiotic resistance genes was performed in the sediments of one fishery farm ponds as well as sampling points upstream and downstream. All the tested sediment samples did not show significantly elevated heavy metal concentrations and no substantial veterinary antibiotic pollution. From the antibiotic resistance genes tested, the presence of aminoglycoside and β-lactam resistance determinants as well as the presence of integrons could be of concern for the possibility of transfer to humans. However, despite the lack of heavy metal and antibiotic pollution, the sediments showed toxicity, the cause of which should be explored more.

Organic Matter Composition at Ocean Station Papa Affects Its Bioavailability, Bacterioplankton Growth Efficiency and the Responding Taxa

The bioavailability of organic matter (OM) to marine heterotrophic bacterioplankton is determined by both the chemical composition of OM and the microbial community composition. In the current study, changes in OM bioavailability were identified at Ocean Station Papa as part of the 2018 Export Processes in the Ocean from Remote Sensing (EXPORTS) field study. Removal rates of carbon (C) in controlled experiments were significantly correlated with the initial composition of total hydrolyzable amino acids, and C removal rates were high when the amino acid degradation index suggested a more labile composition. Carbon remineralization rates averaged 0.19 ± 0.08 μmol C L-1 d-1 over 6-10 days while bacterial growth efficiencies averaged 31 ± 7%. Amino acid composition and tandem mass spectrometry analysis of compound classes also revealed transformations to a more degraded OM composition during experiments. There was a log2-fold increase in the relative abundances of 16S rDNA-resolved bacterioplankton taxa in most experiments by members of the Methylophilaceae family (OM43 genus) and KI89A order. Additionally, when OM was more bioavailable, relative abundances increased by at least threefold for the classes Bacteroidetes (Flavobacteriaceae NS2b genus), Alphaproteobacteria (Rhodobacteraceae Sulfitobacter genus), and Gammaproteobacteria (Alteromonadales and Ectothiorhodospiraceae orders). Our data suggest that a diverse group of bacterioplankton was responsible for removing organic carbon and altering the OM composition to a more degraded state. Elevated community diversity, as inferred from the Shannon-Wiener H index, may have contributed to relatively high growth efficiencies by the bacterioplankton. The data presented here shed light on the interconnections between OM bioavailability and key bacterioplankton taxa for the degradation of marine OM.

Comprehensive single-PCR 16S and 18S rRNA community analysis validated with mock communities, and estimation of sequencing bias against 18S

Universal primers for SSU rRNA genes allow profiling of natural communities by simultaneously amplifying templates from Bacteria, Archaea, and Eukaryota in a single PCR reaction. Despite the potential to show relative abundance for all rRNA genes, universal primers are rarely used, due to various concerns including amplicon length variation and its effect on bioinformatic pipelines. We thus developed 16S and 18S rRNA mock communities and a bioinformatic pipeline to validate this approach. Using these mocks, we show that universal primers (515Y/926R) outperformed eukaryote-specific V4 primers in observed versus expected abundance correlations (slope = 0.88 vs. 0.67-0.79), and mock community members with single mismatches to the primer were strongly underestimated (threefold to eightfold). Using field samples, both primers yielded similar 18S beta-diversity patterns (Mantel test, p < 0.001) but differences in relative proportions of many rarer taxa. To test for length biases, we mixed mock communities (16S + 18S) before PCR and found a twofold underestimation of 18S sequences due to sequencing bias. Correcting for the twofold underestimation, we estimate that, in Southern California field samples (1.2-80 μm), there were averages of 35% 18S, 28% chloroplast 16S, and 37% prokaryote 16S rRNA genes. These data demonstrate the potential for universal primers to generate comprehensive microbiome profiles.

Description and comparison of the skin and ear canal microbiota of non-allergic and allergic German shepherd dogs using next generation sequencing

Atopic dermatitis is one of the most common skin diseases in dogs. Pathogenesis is complex and incompletely understood. Skin colonizing bacteria likely play an important role in the severity of this disease. Studying the canine skin microbiota using traditional microbiological methods has many limitations which can be overcome by molecular procedures. The aim of this study was to describe the bacterial microbiota of the skin and ear canals of healthy non-allergic and allergic German shepherd dogs (GSDs) without acute flare or concurrent skin infection and to compare both. Bacterial 16S rRNA gene amplicon sequence data revealed no differences of bacterial community patterns between the different body sites (axilla, front dorsal interdigital skin, groin, and ear canals) in non-allergic dogs. The microbiota at the different body sites of non-allergic GSDs showed no significant differences. Only for the samples obtained from the axilla the bacterial microbiota of allergic dogs was characterized by a lower species richness compared to that of non-allergic dogs and the bacterial community composition of the skin and ear canals of allergic dogs showed body site specific differences compared to non-allergic dogs. Actinobacteria was the most abundant phylum identified from the non-allergic dogs and Proteobacteria from allergic dogs. Macrococcus spp. were more abundant on non-allergic skin while Sphingomonas spp. were more abundant on the allergic skin. Forward step redundancy analysis of metadata indicated that the household the dogs came from had the strongest impact on the composition of the skin microbiome followed by sex, host health status and body site.

In-depth Spatiotemporal Characterization of Planktonic Archaeal and Bacterial Communities in North and South San Francisco Bay

Despite being the largest estuary on the west coast of North America, no in-depth survey of microbial communities in San Francisco Bay (SFB) waters currently exists. In this study, we analyze bacterioplankton and archaeoplankton communities at several taxonomic levels and spatial extents (i.e., North versus South Bay) to reveal patterns in alpha and beta diversity. We assess communities using high-throughput sequencing of the 16S rRNA gene in 177 water column samples collected along a 150-km transect over a 2-year monthly time-series. In North Bay, the microbial community is strongly structured by spatial salinity changes while in South Bay seasonal variations dominate community dynamics. Along the steep salinity gradient in North Bay, we find that operational taxonomic units (OTUs; 97% identity) have higher site specificity than at coarser taxonomic levels and turnover ("species" replacement) is high, revealing a distinct brackish community (in oligo-, meso-, and polyhaline samples) from fresh and marine end-members. At coarser taxonomic levels (e.g., phylum, class), taxa are broadly distributed across salinity zones (i.e., present/abundant in a large number of samples) and brackish communities appear to be a mix of fresh and marine communities. We also observe variations in brackish communities between samples with similar salinities, likely related to differences in water residence times between North and South Bay. Throughout SFB, suspended particulate matter is positively correlated with richness and influences changes in beta diversity. Within several abundant groups, including the SAR11 clade (comprising up to 30% of reads in a sample), OTUs appear to be specialized to a specific salinity range. Some other organisms also showed pronounced seasonal abundance, including Synechococcus, Ca. Actinomarina, and Nitrosopumilus-like OTUs. Overall, this study represents the first in-depth spatiotemporal survey of SFB microbial communities and provides insight into how planktonic microorganisms have specialized to different niches along the salinity gradient.

Terrestrial Inputs Shape Coastal Bacterial and Archaeal Communities in a High Arctic Fjord (Isfjorden, Svalbard)

The Arctic is experiencing dramatic changes including increases in precipitation, glacial melt, and permafrost thaw, resulting in increasing freshwater runoff to coastal waters. During the melt season, terrestrial runoff delivers carbon- and nutrient-rich freshwater to Arctic coastal waters, with unknown consequences for the microbial communities that play a key role in determining the cycling and fate of terrestrial matter at the land-ocean interface. To determine the impacts of runoff on coastal microbial (bacteria and archaea) communities, we investigated changes in pelagic microbial community structure between the early (June) and late (August) melt season in 2018 in the Isfjorden system (Svalbard). Amplicon sequences of the 16S rRNA gene were generated from water column, river and sediment samples collected in Isfjorden along fjord transects from shallow river estuaries and glacier fronts to the outer fjord. Community shifts were investigated in relation to environmental gradients, and compared to river and marine sediment microbial communities. We identified strong temporal and spatial reorganizations in the structure and composition of microbial communities during the summer months in relation to environmental conditions. Microbial diversity patterns highlighted a reorganization from rich communities in June toward more even and less rich communities in August. In June, waters enriched in dissolved organic carbon (DOC) provided a niche for copiotrophic taxa including Sulfitobacter and Octadecabacter. In August, lower DOC concentrations and Atlantic water inflow coincided with a shift toward more cosmopolitan taxa usually associated with summer stratified periods (e.g., SAR11 Clade Ia), and prevalent oligotrophic marine clades (OM60, SAR92). Higher riverine inputs of dissolved inorganic nutrients and suspended particulate matter also contributed to spatial reorganizations of communities in August. Sentinel taxa of this late summer fjord environment included taxa from the class Verrucomicrobiae (Roseibacillus, Luteolibacter), potentially indicative of a higher fraction of particle-attached bacteria. This study highlights the ecological relevance of terrestrial runoff for Arctic coastal microbial communities and how its impacts on biogeochemical conditions may make these communities susceptible to climate change.

Comparison of Two 16S rRNA Primers (V3-V4 and V4-V5) for Studies of Arctic Microbial Communities

Microbial communities of the Arctic Ocean are poorly characterized in comparison to other aquatic environments as to their horizontal, vertical, and temporal turnover. Yet, recent studies showed that the Arctic marine ecosystem harbors unique microbial community members that are adapted to harsh environmental conditions, such as near-freezing temperatures and extreme seasonality. The gene for the small ribosomal subunit (16S rRNA) is commonly used to study the taxonomic composition of microbial communities in their natural environment. Several primer sets for this marker gene have been extensively tested across various sample sets, but these typically originated from low-latitude environments. An explicit evaluation of primer-set performances in representing the microbial communities of the Arctic Ocean is currently lacking. To select a suitable primer set for studying microbiomes of various Arctic marine habitats (sea ice, surface water, marine snow, deep ocean basin, and deep-sea sediment), we have conducted a performance comparison between two widely used primer sets, targeting different hypervariable regions of the 16S rRNA gene (V3-V4 and V4-V5). We observed that both primer sets were highly similar in representing the total microbial community composition down to genus rank, which was also confirmed independently by subgroup-specific catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) counts. Each primer set revealed higher internal diversity within certain bacterial taxonomic groups (e.g., the class Bacteroidia by V3-V4, and the phylum Planctomycetes by V4-V5). However, the V4-V5 primer set provides concurrent coverage of the archaeal domain, a relevant component comprising 10-20% of the community in Arctic deep waters and the sediment. Although both primer sets perform similarly, we suggest the use of the V4-V5 primer set for the integration of both bacterial and archaeal community dynamics in the Arctic marine environment.

Different responses of soil bacterial and fungal communities to nitrogen deposition in a subtropical forest

China has experienced a widespread increase in N deposition due to intensive anthropogenic activities, particularly in the subtropical regions. However, the effects of long-term N deposition on soil bacterial and fungal abundance, diversity, and community composition remain largely unclear. We assessed the effects of N deposition on soil microbial communities in summer and winter, using quantitative polymerase chain reaction and Illumina Miseq sequencing of bacterial 16S rRNA and fungal ITS genes from subtropical natural forest soils. The abundance of both soil bacteria and fungi exhibited a decreasing pattern with increasing N deposition rates. Nitrogen deposition increased bacterial diversity in both summer and winter, whereas fungal diversity was significantly decreased in summer, but greatly increased under the highest level of N deposition (150 kg N ha^-1 yr^-1) in winter. Nitrogen deposition significantly increased the relative abundance of bacterial phyla Actinobacteria, Chloroflexi, and WPS-2, but decreased that of Acidobacteria and Verrucomicrobia. In addition, N deposition significantly decreased the relative abundance of Ascomycetes, but did not exert any significant effect on Basidiomycetes. The bacterial and fungal community compositions were greatly influenced by N deposition, with soil N availability and soil pH identified as the two most influential soil properties. This study demonstrates that the fungal community was more sensitive than the bacterial community to N deposition, and further emphasizes the importance of simultaneously evaluating soil bacterial and fungal communities in response to global environmental changes.

Microbiome Analysis Using 16S Amplicon Sequencing: From Samples to ASVs

In this chapter, we will present an outline of a typical experimental and bioinformatic workflow for identification of bacterial amplicon sequence variants (ASVs) present in a set of samples. This chapter is written from a bioinformatic point of view; therefore, the specific experimental protocols are not detailed, but rather the impact of various experimental decisions on the downstream analysis is described. Emphasis is made on the transition from reads to ASVs, describing the Deblur algorithm.

The Response of Estuarine Ammonia-Oxidizing Communities to Constant and Fluctuating Salinity Regimes

Aerobic nitrification is a fundamental nitrogen biogeochemical process that links the oxidation of ammonia to the removal of fixed nitrogen in eutrophicated water bodies. However, in estuarine environments there is an enormous variability of water physicochemical parameters that can affect the ammonia oxidation biological process. For instance, it is known that salinity can affect nitrification performance, yet there is still a lack of information on the ammonia-oxidizing communities behavior facing daily salinity fluctuations. In this work, laboratory experiments using upstream and downstream estuarine sediments were performed to address this missing gap by comparing the effect of daily salinity fluctuations with constant salinity on the activity and diversity of ammonia-oxidizing microorganisms (AOM). Activity and composition of AOM were assessed, respectively by using nitrogen stable isotope technique and 16S rRNA gene metabarcoding analysis. Nitrification activity was negatively affected by daily salinity fluctuations in upstream sediments while no effect was observed in downstream sediments. Constant salinity regime showed clearly higher rates of nitrification in upstream sediments while a similar nitrification performance between the two salinity regimes was registered in the downstream sediments. Results also indicated that daily salinity fluctuation regime had a negative effect on both ammonia-oxidizing bacteria (AOB) and ammonia-oxidizing archaea (AOA) community’s diversity. Phylogenetically, the estuarine downstream AOM were dominated by AOA (0.92–2.09%) followed by NOB (0.99–2%), and then AOB (0.2–0.32%); whereas NOB dominated estuarine upstream sediment samples (1.4–9.5%), followed by AOA (0.27–0.51%) and AOB (0.01–0.23%). Analysis of variance identified the spatial difference between samples (downstream and upstream) as the main drivers of AOA and AOB diversity. Our study indicates that benthic AOM inhabiting different estuarine sites presented distinct plasticity toward the salinity regimes tested. These findings help to improve our understanding in the dynamics of the nitrogen cycle of estuarine systems by showing the resilience and consequently the impact of different salinity regimes on the diversity and activity of ammonia oxidizer communities.

Bacterial Taxa Migrating from the Mediterranean Sea into the Red Sea Revealed a Higher Prevalence of Anti-Lessepsian Migrations

In 1869, the Suez Canal was opened, which brought the waters of the Mediterranean and the Red Sea into direct contact. Notably, the Suez Canal was constructed for navigation purposes without focusing on the ecological impacts. The Suez Canal paved the way for species migration from the Red Sea to the Mediterranean Sea through Lessepsian migration, named after Ferdinand de Lesseps, while the migration from the Mediterranean Sea to the Red Sea is called the anti-Lessepsian migration. It has been argued in the past that the migrating species had negative consequences for the host environment as well as of humans. Few studies to date have attempted to map the microorganism migration problem because the traditional ways of measuring the community's richness and dissimilarities failed to provide enough detection of the migrating taxa. We collected 22 seawater samples from different locations in Egypt, in relationship to the migration across and to/from the Suez Canal. The V3-V4 regions of 16s genes were amplified and sequenced by the next generation Illumina MiSeq sequencer. Bioinformatics analysis revealed 15 taxa that migrated from the Mediterranean Sea to the Red Sea (i.e., anti-Lessepsian migration) such as the genera Fluvicola, HTCC2207, and Persicirhabdus. The family OCS155 is the only one that migrated from the Red Sea to the Mediterranean Sea (Lessepsian migration). Seven anti-Lessepsian migrants colonized the Suez Canal more than the Mediterranean Sea such as the genera Marinobacter and Halomonas. These findings collectively suggest that the anti-Lessepsian migration is more predominant than the Lessepsian migration in the bacterial community. This study paves the way for future research questions as well. For example, why is the anti-Lessepsian migration more common than the Lessepsian route in bacteria? Why do certain taxa stop migration at the Suez Canal, and why do certain taxa present in higher frequencies in the Suez Canal? Which taxa continue migration to the Indian Ocean and the Atlantic Ocean, and what is the impact of the anti-Lessepsian migration on the bacterial community? Understanding microbial diversity in a context of microorganism migration across seas and oceans remains a prime topic in biodiversity research and systems science.

Stable Isotope Probing Identifies Bacterioplankton Lineages Capable of Utilizing Dissolved Organic Matter Across a Range of Bioavailability

Bacterioplankton consume about half of the dissolved organic matter (DOM) produced by phytoplankton. DOM released from phytoplankton consists of a myriad of compounds that span a range of biological reactivity from labile to recalcitrant. Linking specific bacterioplankton lineages to the incorporation of DOM compounds into biomass is important to understand microbial niche partitioning. We conducted a series of DNA-stable isotope probing (SIP) experiments using ¹³C-labeled substrates of varying lability including amino acids, cyanobacteria lysate, and DOM from diatom and cyanobacteria isolates concentrated on solid phase extraction PPL columns (SPE-DOM). Amendments of substrates into Sargasso Sea bacterioplankton communities were conducted to explore microbial response and DNA-SIP was used to determine which lineages of Bacteria and Archaea were responsible for uptake and incorporation. Greater increases in bacterioplankton abundance and DOC removal were observed in incubations amended with cyanobacteria-derived lysate and amino acids compared to the SPE-DOM, suggesting that the latter retained proportionally more recalcitrant DOM compounds. DOM across a range of bioavailability was utilized by diverse prokaryotic taxa with copiotrophs becoming the most abundant ¹³C-incorporating taxa in the amino acid treatment and oligotrophs becoming the most abundant ¹³C-incorporating taxa in SPE-DOM treatments. The lineages that responded to SPE-DOM amendments were also prevalent in the mesopelagic of the Sargasso Sea, suggesting that PPL extraction of phytoplankton-derived DOM isolates compounds of ecological relevance to oligotrophic heterotrophic bacterioplankton. Our study indicates that DOM quality is an important factor controlling the diversity of the microbial community response, providing insights into the roles of different bacterioplankton in resource exploitation and efficiency of marine carbon cycling.

Naught all zeros in sequence count data are the same

Genomic studies feature multivariate count data from high-throughput DNA sequencing experiments, which often contain many zero values. These zeros can cause artifacts for statistical analyses and multiple modeling approaches have been developed in response. Here, we apply different zero-handling models to gene-expression and microbiome datasets and show models can disagree substantially in terms of identifying the most differentially expressed sequences. Next, to rationally examine how different zero handling models behave, we developed a conceptual framework outlining four types of processes that may give rise to zero values in sequence count data. Last, we performed simulations to test how zero handling models behave in the presence of these different zero generating processes. Our simulations showed that simple count models are sufficient across multiple processes, even when the true underlying process is unknown. On the other hand, a common zero handling technique known as "zero-inflation" was only suitable under a zero generating process associated with an unlikely set of biological and experimental conditions. In concert, our work here suggests several specific guidelines for developing and choosing state-of-the-art models for analyzing sparse sequence count data.

Niche differentiation of denitrifying anaerobic methane oxidizing bacteria and archaea leads to effective methane filtration in a Tibetan alpine wetland

Denitrifying anaerobic methane oxidation (DAMO) is a vital methane sink in wetlands. However, the interactions and niche partitioning of DAMO bacteria and archaea in freshwater wetland soils, in addition to the interactions among microorganisms that couple methane and nitrogen cycling is still unclear, despite that these factors may govern the fate of methane and nitrogen in wetlands. Here, we evaluated the vertical distribution of DAMO bacteria and archaea in soil layers along with the potential interactions among populations in the methane-coupled nitrogen cycling microbial community of Tibetan freshwater wetlands. A combination of molecular biology, stable isotope tracer technology, and microbial bioinformatics was used to evaluate these interrelated dynamics. The abundances and potential methane oxidation rates indicated that DAMO bacteria and archaea differentially occupy surface and subsurface soil layers, respectively. The inferred interactions between DAMO bacteria and nitrogen cycling microorganisms within their communities are complex, DAMO bacteria apparently achieve an advantage in the highly competitive environment of surface soils layers and occupy a specific niche in those environments. Conversely, the apparent relationships between DAMO archaea and nitrogen cycling microorganisms are relatively simple, wherein high levels of cooperation are inferred between DAMO archaea and nitrate-producing organisms in subsurface soils layers. These results suggest that the vertical distribution patterns of DAMO bacteria and archaea enable them to play significant roles in the methane oxidation activity of different soil layers and collectively form an effective methane filtration consortium.

Chinese Liquor Fermentation: Identification of Key Flavor-Producing Lactobacillus spp. by Quantitative Profiling with Indigenous Internal Standards

Identifying the functional microbes in spontaneous food fermentation is important for improving food quality. To identify the key flavor producers in Chinese liquor fermentation, we propose a novel quantitative microbiome profiling method that uses indigenous internal standards to normalize high-throughput amplicon sequencing results. We screened Lactobacillus acetotolerans and Lactobacillus jinshani as indigenous internal standards based on their high distribution frequencies and relative abundances. After determining the absolute abundance of indigenous internal standards using quantitative PCR with species-specific primers, the liquor-fermented bacterial community and its dynamics were better characterized by internal standards normalization. Based on quantitative microbiome profiling, we identified that Lactobacillus was a key flavor producer correlated with eight flavor compounds. Metatranscriptomic analysis indicated that Lactobacillus was active in transcribing genes involving the biosynthesis of flavor compounds and their precursors. This work has developed a novel and extensible absolute quantification method for microbiota that will alleviate concerns in the statistical analyses based on relative microbiome profiling, and shed insights into the function of Lactobacillus in food fermentation. It can potentially be applied to other microbial ecology studies.IMPORTANCE In this study, we developed a novel strategy using indigenous internal standards to normalize the high-throughput amplicon sequencing results. We chose two Lactobacillus species as indigenous internal standards and characterized the absolute abundance of the bacterial community. Further, we identified Lactobacillus as the key flavor producer using quantitative microbiome profiling combined with multivariate statistics and metatranscriptomic analysis. This work developed a novel strategy for absolute quantitative abundance analysis of microbiota and expanded our understanding of the role of Lactobacillus in food fermentation.

Interactions between elevated CO2 levels and floating aquatic plants on the alteration of bacterial function in carbon assimilation and decomposition in eutrophic waters

Elevated atmospheric CO₂ concentration (eCO₂) may have different effects on the bacterial community with regard to C assimilation and decomposition in eutrophic waters compared to that in fresh waters with intermediate levels of nutrients and oceans. Aquatic plant growth under eCO₂ could further modify microbial activities associated with the C cycle in eutrophic waters. Therefore, there is an urgent need to further study how eCO₂ and its interactions with the growth of aquatic plants affect the composition and function of the bacterial community involved in mediating the C cycle in eutrophic waters. Accordingly, we designed a microcosm experiment to investigate the effects of ambient and high CO₂ concentrations on bacterial community composition and function in eutrophic waters with and without the growth of Eichhornia crassipes (Mart.) Solms. The results from 16S rRNA gene sequencing, function prediction, and q-PCR showed that eCO₂ significantly increased the abundance of bacterial and functional genes involved in CO₂ assimilation (photosynthetic bacteria; cbbL IA & IC, cbbL ID, cbbM, pufM) and C decomposition (Acidimicrobiia, Thermoleophilia, Gaiellales; ChiA), illustrating the functional enrichment with photoautotrophy, hydrocarbon degradation, cellulolysis, and aromatic hydrocarbon degradation. However, eCO₂ decreased the abundance of some chemoautotrophic bacteria, including nitrifying bacteria (Nitrospirae, Nitrosomonadaceae). In contrast, the cultivation of E. crassipes decreased the abundance of photosynthetic bacteria but increased the abundance of bacteria involved in complex C decomposition associated with root exudates and degradation, e.g. Fibrobacteres, Sphingobacteriales, Sphingomonadales, and Rhizobiales. eCO₂ and growth of E. crassipes had opposite effects on algal density in eutrophic waters, creating interactive effects that further decreased the diversity of the bacterial community and abundance of some CO₂-assimilating bacteria with nitrifying characteristics (Nitrosomonadaceae) and some C-degrading bacteria (Fibrobacteres) with denitrifying properties (Flavobacteriaceae, Sphingomonadaceae, and Gemmobacter). Therefore, the interactions between aquatic plants and the bacterial community in eutrophic waters under eCO₂ would be beneficial to the environment and help alleviate the greenhouse effect.