rrnDB: documenting the number of rRNA and tRNA genes in bacteria and archaea

Residue quality drives SOC sequestration by altering microbial taxonomic composition and ecophysiological function in desert ecosystem

Plant residues are important sources of soil organic carbon in terrestrial ecosystems. The degradation of plant residue by microbes can influence the soil carbon cycle and sequestration. However, little is known about the microbial composition and function, as well as the accumulation of soil organic carbon (SOC) in response to the inputs of different quality plant residues in the desert environment. The present study evaluated the effects of plant residue addition from Pinus sylvestris var. mongolica (Pi), Artemisia desertorum (Ar) and Amorpha fruticosa (Am) on desert soil microbial community composition and function in a field experiment in the Mu Us Desert. The results showed that the addition of the three plant residues with different C/N ratios induced significant variation in soil microbial communities. The Am treatment (low C/N ratio) improved microbial diversity compared with the Ar and Pi treatments (medium and high C/N ratios). The variations in the taxonomic and functional compositions of the dominant phyla Actinobacteria and Proteobacteria were higher than those of the other phyla among the different treatments. Moreover, the network links between Proteobacteria and other phyla and the CAZyme genes abundances from Proteobacteria increased with increasing residue C/N, whereas those decreased for Actinobacteria. The SOC content of the Am, Ar and Pi treatments increased by 45.73%, 66.54% and 107.99%, respectively, as compared to the original soil. The net SOC accumulation was positively correlated with Proteobacteria abundance and negatively correlated with Actinobacteria abundance. These findings showed that changing the initial quality of plant residue from low C/N to high C/N can result in shifts in taxonomic and functional composition from Actinobacteria to Proteobacteria, which favors SOC accumulation. This study elucidates the ecophysiological roles of Actinobacteria and Proteobacteria in the desert carbon cycle, expands our understanding of the potential microbial-mediated mechanisms by which plant residue inputs affect SOC sequestration in desert soils, and provides valuable guidance for species selection in desert vegetation reconstruction.

A longitudinal field study of commercial honey bees shows that non-native probiotics do not rescue antibiotic treatment, and are generally not beneficial

Probiotics are widely used in agriculture including commercial beekeeping, but there is little evidence supporting their effectiveness. Antibiotic treatments can greatly distort the gut microbiome, reducing its protective abilities and facilitating the growth of antibiotic resistant pathogens. Commercial beekeepers regularly apply antibiotics to combat bacterial infections, often followed by an application of non-native probiotics advertised to ease the impact of antibiotic-induced gut dysbiosis. We tested whether probiotics affect the gut microbiome or disease prevalence, or rescue the negative effects of antibiotic induced gut dysbiosis. We found no difference in the gut microbiome or disease markers by probiotic application or antibiotic recovery associated with probiotic treatment. A colony-level application of the antibiotics oxytetracycline and tylosin produced an immediate decrease in gut microbiome size, and over the longer-term, very different and persistent dysbiotic effects on the composition and membership of the hindgut microbiome. Our results demonstrate the lack of probiotic effect or antibiotic rescue, detail the duration and character of dysbiotic states resulting from different antibiotics, and highlight the importance of the gut microbiome for honeybee health.

Bacteria can maintain rRNA operons solely on plasmids for hundreds of millions of years

It is generally assumed that all bacteria must have at least one rRNA operon (rrn operon) on the chromosome, but some strains of the genera Aureimonas and Oecophyllibacter carry their sole rrn operon on a plasmid. However, other related strains and species have chromosomal rrn loci, suggesting that the exclusive presence of rrn operons on a plasmid is rare and unlikely to be stably maintained over long evolutionary periods. Here, we report the results of a systematic search for additional bacteria without chromosomal rrn operons. We find that at least four bacterial clades in the phyla Bacteroidota, Spirochaetota, and Pseudomonadota (Proteobacteria) lost chromosomal rrn operons independently. Remarkably, Persicobacteraceae have apparently maintained this peculiar genome organization for hundreds of millions of years. In our study, all the rrn-carrying plasmids in bacteria lacking chromosomal rrn loci possess replication initiator genes of the Rep_3 family. Furthermore, the lack of chromosomal rrn operons is associated with differences in copy numbers of rrn operons, plasmids, and chromosomal tRNA genes. Thus, our findings indicate that the absence of rrn loci in bacterial chromosomes can be stably maintained over long evolutionary periods.

Characterization of SARS-CoV-2 Distribution and Microbial Succession in a Clinical Microbiology Testing Facility during the SARS-CoV-2 Pandemic

The exchange of microbes between humans and the built environment is a dynamic process that has significant impact on health. Most studies exploring the microbiome of the built environment have been predicated on improving our understanding of pathogen emergence, persistence, and transmission. Previous studies have demonstrated that SARS-CoV-2 presence significantly correlates with the proportional abundance of specific bacteria on surfaces in the built environment. However, in these studies, SARS-CoV-2 originated from infected patients. Here, we perform a similar assessment for a clinical microbiology lab while staff were handling SARS-CoV-2 infected samples. The goal of this study was to understand the distribution and dynamics of microbial population on various surfaces within different sections of a clinical microbiology lab during a short period of 2020 Coronavirus disease (COVID-19) pandemic. We sampled floors, benches, and sinks in 3 sections (bacteriology, molecular microbiology, and COVID) of an active clinical microbiology lab over a 3-month period. Although floor samples harbored SARS-CoV-2, it was rarely identified on other surfaces, and bacterial diversity was significantly greater on floors than sinks and benches. The floors were primarily colonized by bacteria common to natural environments (e.g., soils), and benchtops harbored a greater proportion of human-associated microbes, including Staphylococcus and Streptococcus. Finally, we show that the microbial composition of these surfaces did not change over time and remained stable. Despite finding viruses on the floors, no lab-acquired infections were reported during the study period, which suggests that lab safety protocols and sanitation practices were sufficient to prevent pathogen exposures. IMPORTANCE For decades, diagnostic clinical laboratories have been an integral part of the health care systems that perform diagnostic tests on patient's specimens in bulk on a regular basis. Understanding their microbiota should assist in designing and implementing disinfection, and cleaning regime in more effective way. To our knowledge, there is a lack of information on the composition and dynamics of microbiota in the clinical laboratory environments, and, through this study, we have tried to fill that gap. This study has wider implications as understanding the makeup of microbes on various surfaces within clinical laboratories could help identify any pathogenic bacterial taxa that could have colonized these surfaces, and might act as a potential source of laboratory-acquired infections. Mapping the microbial community within these built environments may also be critical in assessing the reliability of laboratory safety and sanitation practices to lower any potential risk of exposures to health care workers.

Genomic Features Predict Bacterial Life History Strategies in Soil, as Identified by Metagenomic Stable Isotope Probing

Bacteria catalyze the formation and destruction of soil organic matter, but the bacterial dynamics in soil that govern carbon (C) cycling are not well understood. Life history strategies explain the complex dynamics of bacterial populations and activities based on trade-offs in energy allocation to growth, resource acquisition, and survival. Such trade-offs influence the fate of soil C, but their genomic basis remains poorly characterized. We used multisubstrate metagenomic DNA stable isotope probing to link genomic features of bacteria to their C acquisition and growth dynamics. We identify several genomic features associated with patterns of bacterial C acquisition and growth, notably genomic investment in resource acquisition and regulatory flexibility. Moreover, we identify genomic trade-offs defined by numbers of transcription factors, membrane transporters, and secreted products, which match predictions from life history theory. We further show that genomic investment in resource acquisition and regulatory flexibility can predict bacterial ecological strategies in soil. IMPORTANCE Soil microbes are major players in the global carbon cycle, yet we still have little understanding of how the carbon cycle operates in soil communities. A major limitation is that carbon metabolism lacks discrete functional genes that define carbon transformations. Instead, carbon transformations are governed by anabolic processes associated with growth, resource acquisition, and survival. We use metagenomic stable isotope probing to link genome information to microbial growth and carbon assimilation dynamics as they occur in soil. From these data, we identify genomic traits that can predict bacterial ecological strategies which define bacterial interactions with soil carbon.

Impact of 16S rRNA Gene Redundancy and Primer Pair Selection on the Quantification and Classification of Oral Microbiota in Next-Generation Sequencing

This study aimed to evaluate the number of 16S rRNA genes in the complete genomes of the bacterial and archaeal species inhabiting the human mouth and to assess how the use of different primer pairs would affect the detection and classification of redundant amplicons and matching amplicons (MAs) from different taxa. A total of 518 oral-bacterial and 191 oral-archaeal complete genomes were downloaded from the NCBI database, and their complete 16S rRNA genes were extracted. The numbers of genes and variants per genome were calculated. Next, 39 primer pairs were used to search for matches in the genomes and obtain amplicons. For each primer, we calculated the number of gene amplicons, variants, genomes, and species detected and the percentage of coverage at the species level with no MAs (SC-NMA). The results showed that 94.09% of oral bacteria and 52.59% of oral archaea had more than one intragenomic 16S rRNA gene. From 1.29% to 46.70% of bacterial species and from 4.65% to 38.89% of archaea detected by the primers had MAs. The best primers were the following (SC-NMA; region; position for Escherichia coli [GenBank version no. J01859.1]): KP_F048-OP_R030 for bacteria (93.55%; V3 to V7; 342 to 1079), KP_F018-KP_R063 for archaea (89.63%; V3 to V9; undefined to 1506), and OP_F114-OP_R121 for both domains (92.52%; V3 to V9; 340 to 1405). In addition to 16S rRNA gene redundancy, the presence of MAs must be controlled to ensure an accurate interpretation of microbial diversity data. The SC-NMA is a more useful parameter than the conventional coverage percentage for selecting the best primer pairs. The pairs used the most in the oral microbiome literature were not among the best performers. IMPORTANCE Hundreds of publications have studied the oral microbiome through 16S rRNA gene sequencing. However, none have assessed the number of 16S rRNA genes in the genomes of oral microbes, or how the use of primer pairs targeting different regions affects the detection of MAs from different taxa. Here, we found that almost all oral bacteria and more than half of oral archaea have more than one intragenomic 16S rRNA gene. The performance of the primer pairs in not detecting MAs increases as the length of the amplicon augments. As none of those most employed in the oral literature were among the best performers, we selected a series of primers to detect bacteria and/or archaea based on their percentage of species detected without MAs. The intragenomic 16S rRNA gene redundancy and the presence of MAs between distinct taxa need to be considered to ensure an accurate interpretation of microbial diversity data.

Genomic Landscape Highlights Molecular Mechanisms Involved in Silicate Solubilization, Stress Tolerance, and Potential Growth-Promoting Activity of Bacterium Enterobacter sp. LR6

Silicon (Si) is gaining widespread attention due to its prophylactic activity to protect plants under stress conditions. Despite Si's abundance in the earth's crust, most soils do not have enough soluble Si for plants to absorb. In the present study, a silicate-solubilizing bacterium, Enterobacter sp. LR6, was isolated from the rhizospheric soil of rice and subsequently characterized through whole-genome sequencing. The size of the LR6 genome is 5.2 Mb with a GC content of 54.9% and 5182 protein-coding genes. In taxogenomic terms, it is similar to E. hormaechei subsp. xiangfangensis based on average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH). LR6 genomic data provided insight into potential genes involved in stress response, secondary metabolite production, and growth promotion. The LR6 genome contains two aquaporins, of which the aquaglyceroporin (GlpF) is responsible for the uptake of metalloids including arsenic (As) and antimony (Sb). The yeast survivability assay confirmed the metalloid transport activity of GlpF. As a biofertilizer, LR6 isolate has a great deal of tolerance to high temperatures (45 °C), salinity (7%), and acidic environments (pH 9). Most importantly, the present study provides an understanding of plant-growth-promoting activity of the silicate-solubilizing bacterium, its adaptation to various stresses, and its uptake of different metalloids including As, Ge, and Si.

First Record of the Rare Species Aeromonas lusitana from Rainbow Trout (Oncorhynchus mykiss, Walbaum): Comparative Analysis with the Existing Strains

The species Aeromonas lusitana was first described in 2016 with five strains recovered from untreated water and vegetables from Portugal. Since then, no further records exist of this species. During a surveillance study on the presence of Aeromonas in fish farms in Mexico, a new strain (ESV-351) of the mentioned species isolated from a rainbow trout was recovered. It was identified because it clustered phylogenetically with the type strain of A. lusitana based on the analysis of the rpoD gene sequences. In the present study, phenotypic characteristics, antimicrobial resistance profiles, and the presence of putative virulence genes of this novel strain (ESV-351) were determined in parallel to the five isolates from the original species description. Phenotypic differential characteristics exhibited by A. lusitana ESV-351 depicted an evident similarity to the characteristics exhibited by the other evaluated strains. However, the novel strain was positive for the production of indole using conventional methods, while the rest of the strains, including the type strain, were negative for its production. Furthermore, intermediate resistance to ampicillin, amoxicillin-clavulanic acid and cephalothin was detected in both the novel and the type strain. Five different virulence-related genes were detected in the novel strain and in the previously described strains, with the type strain exhibiting the highest number of virulence-related genes. In addition to this, the genome of the novel strain (ESV-351) was sequenced and compared with the genomes from the type strain (A. lusitana CECT 7828T) and other Aeromonas spp. The genomic analysis defined Aeromonas tecta as the closest species to A. lusitana with a highly similar number of predicted proteins. The genomic size, the number of protein-encoding genes and the number of different tRNAs, among other characteristics, make it possible to propose that the ESV-351 strain could potentially have the capacity to adapt to different environments. Genome comparison of the ESV-351 strain with the type strain revealed that both possess a similar sequence of the citrate synthase gene. In addition to this finding, the chromosomal region containing the citrate synthase locus of the novel strain exhibits some similarity to the chromosomal region in the genome of the A. hydrophila type strain and other known human pathogens, such as Vibrio cholerae. This could suggest a possible virulence role for the citrate synthase gene in A. lusitana (ESV-351).

Effect of Intragenomic Sequence Heterogeneity among Multiple 16S rRNA Genes on Species Identification of Elizabethkingia

Accurate identification of Elizabethkingia species mostly requires the use of molecular techniques, and 16S rRNA gene sequencing is generally considered the method of choice. In this study, we evaluated the effect of intraspecific diversity among the multiple copies of the 16S rRNA gene on the accuracy of species identification in the genus Elizabethkingia. Sequences of 16S rRNA genes obtained from the 32 complete whole-genome sequences of Elizabethkingia deposited in GenBank and from 218 clinical isolates collected from 5 hospitals in Taiwan were analyzed. Four or five copies of 16S rRNA were identified in the Elizabethkingia species with complete genome sequences. The dissimilarity among the copies of the16S rRNA gene was <1% in all Elizabethkingia strains. E. meningoseptica demonstrated a significantly higher rate of nucleotide variations in the 16S rRNA than did E. anophelis (P = 0.011). Nucleotide alterations occurred more frequently in regions V2 and V6 than in other hypervariable regions (P < 0.001). E. meningoseptica, E. anophelis, and E. argenteiflava strains were clustered distinctly in the phylogenetic tree inferred from 16S rRNA genes, and the intragenomic variation of gene sequences had no profound effect on the classification of taxa. However, E. miricola, E. bruuniana, E. ursingii, and E. occulta were grouped closely in the phylogenetic analysis, and the variation among the multiple copies of the 16S rRNA in one E. ursingii strain affected species classification. Other marker genes may be required to supplement the species classification of closely related taxa in the genus Elizabethkingia. IMPORTANCE Incorrect identification of bacterial species would influence the epidemiology and clinical analysis of patients infected with Elizabethkingia. The results of the present study suggest that 16S rRNA gene sequencing should not be considered the gold standard for the accurate identification of Elizabethkingia species.

LotuS2: an ultrafast and highly accurate tool for amplicon sequencing analysis

BACKGROUND

Amplicon sequencing is an established and cost-efficient method for profiling microbiomes. However, many available tools to process this data require both bioinformatics skills and high computational power to process big datasets. Furthermore, there are only few tools that allow for long read amplicon data analysis. To bridge this gap, we developed the LotuS2 (less OTU scripts 2) pipeline, enabling user-friendly, resource friendly, and versatile analysis of raw amplicon sequences.

RESULTS

In LotuS2, six different sequence clustering algorithms as well as extensive pre- and post-processing options allow for flexible data analysis by both experts, where parameters can be fully adjusted, and novices, where defaults are provided for different scenarios. We benchmarked three independent gut and soil datasets, where LotuS2 was on average 29 times faster compared to other pipelines, yet could better reproduce the alpha- and beta-diversity of technical replicate samples. Further benchmarking a mock community with known taxon composition showed that, compared to the other pipelines, LotuS2 recovered a higher fraction of correctly identified taxa and a higher fraction of reads assigned to true taxa (48% and 57% at species; 83% and 98% at genus level, respectively). At ASV/OTU level, precision and F-score were highest for LotuS2, as was the fraction of correctly reported 16S sequences.

CONCLUSION

LotuS2 is a lightweight and user-friendly pipeline that is fast, precise, and streamlined, using extensive pre- and post-ASV/OTU clustering steps to further increase data quality. High data usage rates and reliability enable high-throughput microbiome analysis in minutes.

AVAILABILITY

LotuS2 is available from GitHub, conda, or via a Galaxy web interface, documented at http://lotus2.earlham.ac.uk/ . Video Abstract.

Towards quantitative metabarcoding of eukaryotic plankton: an approach to improve 18S rRNA gene copy number bias

Plankton metabarcoding is increasingly implemented in marine ecosystem assessments and is more cost-efficient and less time-consuming than monitoring based on microscopy (morphological). 18S rRNA gene is the most widely used marker for groups’ and species’ detection and classification within marine eukaryotic microorganisms. These datasets have commonly relied on the acquisition of organismal abundances directly from the number of DNA sequences (i.e. reads). Besides the inherent technical biases in metabarcoding, the largely varying 18S rRNA gene copy numbers (GCN) among marine protists (ranging from tens to thousands) is one of the most important biological biases for species quantification. In this work, we present a gene copy number correction factor (CF) for four marine planktonic groups: Bacillariophyta, Dinoflagellata, Ciliophora miscellaneous and flagellated cells. On the basis of the theoretical assumption that ‘1 read’ is equivalent to ‘1 GCN’, we used the GCN median values per plankton group to calculate the corrected cell number and biomass relative abundances. The species-specific absolute GCN per cell were obtained from various studies published in the literature. We contributed to the development of a species-specific 18S rRNA GCN database proposed by previous authors. To assess the efficiency of the correction factor we compared the metabarcoding, morphological and corrected relative abundances (in cell number and biomass) of 15 surface water samples collected in the Belgian Coastal Zone. Results showed that the application of the correction factor over metabarcoding results enables us to significantly improve the estimates of cell abundances for Dinoflagellata, Ciliophora and flagellated cells, but not for Bacillariophyta. This is likely to due to large biovolume plasticity in diatoms not corresponding to genome size and gene copy numbers. C-biomass relative abundance estimations directly from amplicon reads were only improved for Dinoflagellata and Ciliophora. The method is still facing biases related to the low number of species GCN assessed. Nevertheless, the increase of species in the GCN database may lead to the refinement of the proposed correction factor.

Long-term nitrogen deposition enhances microbial capacities in soil carbon stabilization but reduces network complexity

BACKGROUND

Anthropogenic activities have increased the inputs of atmospheric reactive nitrogen (N) into terrestrial ecosystems, affecting soil carbon stability and microbial communities. Previous studies have primarily examined the effects of nitrogen deposition on microbial taxonomy, enzymatic activities, and functional processes. Here, we examined various functional traits of soil microbial communities and how these traits are interrelated in a Mediterranean-type grassland administrated with 14 years of 7 g m^-2 year^-1 of N amendment, based on estimated atmospheric N deposition in areas within California, USA, by the end of the twenty-first century.

RESULTS

Soil microbial communities were significantly altered by N deposition. Consistent with higher aboveground plant biomass and litter, fast-growing bacteria, assessed by abundance-weighted average rRNA operon copy number, were favored in N deposited soils. The relative abundances of genes associated with labile carbon (C) degradation (e.g., amyA and cda) were also increased. In contrast, the relative abundances of functional genes associated with the degradation of more recalcitrant C (e.g., mannanase and chitinase) were either unchanged or decreased. Compared with the ambient control, N deposition significantly reduced network complexity, such as average degree and connectedness. The network for N deposited samples contained only genes associated with C degradation, suggesting that C degradation genes became more intensely connected under N deposition.

CONCLUSIONS

We propose a conceptual model to summarize the mechanisms of how changes in above- and belowground ecosystems by long-term N deposition collectively lead to more soil C accumulation. Video Abstract.

Metagenomic Analysis Using Phylogenetic Placement-A Review of the First Decade

Phylogenetic placement refers to a family of tools and methods to analyze, visualize, and interpret the tsunami of metagenomic sequencing data generated by high-throughput sequencing. Compared to alternative (e. g., similarity-based) methods, it puts metabarcoding sequences into a phylogenetic context using a set of known reference sequences and taking evolutionary history into account. Thereby, one can increase the accuracy of metagenomic surveys and eliminate the requirement for having exact or close matches with existing sequence databases. Phylogenetic placement constitutes a valuable analysis tool per se, but also entails a plethora of downstream tools to interpret its results. A common use case is to analyze species communities obtained from metagenomic sequencing, for example via taxonomic assignment, diversity quantification, sample comparison, and identification of correlations with environmental variables. In this review, we provide an overview over the methods developed during the first 10 years. In particular, the goals of this review are 1) to motivate the usage of phylogenetic placement and illustrate some of its use cases, 2) to outline the full workflow, from raw sequences to publishable figures, including best practices, 3) to introduce the most common tools and methods and their capabilities, 4) to point out common placement pitfalls and misconceptions, 5) to showcase typical placement-based analyses, and how they can help to analyze, visualize, and interpret phylogenetic placement data.

Comparison of an automated DNA extraction and 16S rDNA real time PCR/sequencing diagnostic method using optimized reagents with culture during a 15-month study using specimens from sterile body sites

Background

16S rDNA-PCR for the identification of a bacterial species is an established method. However, the DNA extraction reagents as well as the PCR reagents may contain residual bacterial DNA, which consequently generates false-positive PCR results. Additionally, previously used methods are frequently time-consuming. Here, we describe the results obtained with a new technology that uses DNA-free reagents for automated DNA extraction and subsequent real time PCR using sterile clinical specimens.

Results

In total, we compared 803 clinical specimens using real time PCR and culturing. The clinical specimens were mainly of orthopedic origin received at our diagnostic laboratory. In 595 (74.1%) samples, the results were concordant negative, and in 102 (12.7%) the results were concordant positive. A total of 170 (21.2%) clinical specimens were PCR-positive, of which 62 (36.5% from PCR positive, 7.7% in total) gave an additional benefit to the patient since only the PCR result was positive. Many of these 62 positive specimens were strongly positive based on crossingpoint values (54% < Cp 30), and these 62 positive clinical specimens were diagnosed as medically relevant as well. Thirty-eight (4.2%) clinical specimens were culture-positive (25 of them were only enrichment culture positive) but PCR-negative, mainly for S. epidermidis, S. aureus and C. acnes. The turnaround times for negative specimens were 4 hours (automated DNA extraction and real time PCR) and 1 working day for positive specimens (including Sanger sequencing). Melting-curve analysis of SYBR Green-PCR enables the differentiation of specific and unspecific PCR products. Using Ripseq, even mixed infections of 2 bacterial species could be resolved.

Conclusions

For endocarditis cases, the added benefit of PCR is obvious. The crucial innovations of the technology enable timely reporting of explicit reliable results for adequate treatment of patients. Clinical specimens with truly PCR-positive but culture-negative results represent an additional benefit for patients. Very few results at the detection limit still have to be critically examined.

Antibiotic Exposure, Common Morbidities and Main Intestinal Microbial Groups in Very Preterm Neonates: A Pilot Study

Prematurity exposes newborns to increased risks of infections and it is associated with critical morbidities. Preterm infants often require antibiotic therapies that can affect the correct establishment of gut microbiota. The aim of this study was to investigate targeted intestinal bacteria in preterm neonates with common morbidities and receiving antibiotic treatments of variable duration. Stool samples were collected after birth, at 15, 30 and 90 days of life. qPCR quantification of selected microbial groups (Bifidobacterium spp., Bacteroides fragilis group, Enterobacteriaceae, Clostridium cluster I and total bacteria) was performed and correlation between their levels, the duration of antibiotic treatment and different clinical conditions was studied. An increasing trend over time was observed for all microbial groups, especially for Bifdobacterium spp. Prolonged exposure to antibiotics in the first weeks of life affected Clostridium and B. fragilis levels, but these changes no longer persisted at 90 days of life. Variations of bacterial counts were associated with the length of hospital stay, feeding and mechanical ventilation. Late-onset sepsis and patent ductus arteriosus reduced the counts of Bifidobacterium, whereas B. fragilis was influenced by compromised respiratory conditions. This study can be a start point for the identification of microbial biomarkers associated with some common morbidities and tailored strategies for a healthy microbial development.

Metagenomic community composition and resistome analysis in a full-scale cold climate wastewater treatment plant

BACKGROUND

Wastewater treatment plants are an essential part of maintaining the health and safety of the general public. However, they are also an anthropogenic source of antibiotic resistance genes. In this study, we characterized the resistome, the distribution of classes 1-3 integron-integrase genes (intI1, intI2, and intI3) as mobile genetic element biomarkers, and the bacterial and phage community compositions in the North End Sewage Treatment Plant in Winnipeg, Manitoba. Samples were collected from raw sewage, returned activated sludge, final effluent, and dewatered sludge. A total of 28 bacterial and viral metagenomes were sequenced over two seasons, fall and winter. Integron-integrase genes, the 16S rRNA gene, and the coliform beta-glucuronidase gene were also quantified during this time period.

RESULTS

Bacterial classes observed above 1% relative abundance in all treatments were Actinobacteria (39.24% ± 0.25%), Beta-proteobacteria (23.99% ± 0.16%), Gamma-proteobacteria (11.06% ± 0.09%), and Alpha-proteobacteria (9.18 ± 0.04%). Families within the Caudovirales order: Siphoviridae (48.69% ± 0.10%), Podoviridae (23.99% ± 0.07%), and Myoviridae (19.94% ± 0.09%) were the dominant phage observed throughout the NESTP. The most abundant bacterial genera (in terms of average percent relative abundance) in influent, returned activated sludge, final effluent, and sludge, respectively, includes Mycobacterium (37.4%, 18.3%, 46.1%, and 7.7%), Acidovorax (8.9%, 10.8%, 5.4%, and 1.3%), and Polaromonas (2.5%, 3.3%, 1.4%, and 0.4%). The most abundant class of antibiotic resistance in bacterial samples was tetracycline resistance (17.86% ± 0.03%) followed by peptide antibiotics (14.24% ± 0.03%), and macrolides (10.63% ± 0.02%). Similarly, the phage samples contained a higher prevalence of macrolide (30.12% ± 0.30%), peptide antibiotic (10.78% ± 0.13%), and tetracycline (8.69% ± 0.11%) resistance. In addition, intI1 was the most abundant integron-integrase gene throughout treatment (1.14 × 104 gene copies/mL) followed by intI3 (4.97 × 103 gene copies/mL) while intI2 abundance remained low (6.4 × 101 gene copies/mL).

CONCLUSIONS

Wastewater treatment successfully reduced the abundance of bacteria, DNA phage and antibiotic resistance genes although many antibiotic resistance genes remained in effluent and biosolids. The presence of integron-integrase genes throughout treatment and in effluent suggests that antibiotic resistance genes could be actively disseminating resistance between both environmental and pathogenic bacteria.

[Environmental Monitoring in a Pharmaceutical Manufacturing Facility Using a Culture Independent Approach]

Strict microbial control is required in manufacturing facilities to ensure the quality of pharmaceuticals and foods. Environmental microbial monitoring plays a fundamental role in reducing the risk of microbial contamination. Appropriate microbial control requires an understanding of abundance and community structures of microbes in the target environment. However, most of these microbes are not culturable using conventional methods. In this study, we determined the number of microbial particles and assessed the environmental microbiome in a pharmaceutical manufacturing facility, using high-throughput sequencing of rRNA gene fragments. Our results provide fundamental data for the evaluation and control of microbes in the pharmaceutical and food industries.

Going Beyond Bacteria: Uncovering the Role of Archaeome and Mycobiome in Inflammatory Bowel Disease

Inflammatory Bowel Disease (IBD) is a multifaceted class of relapsing-remitting chronic inflammatory conditions where microbiota dysbiosis plays a key role during its onset and progression. The human microbiota is a rich community of bacteria, viruses, fungi, protists, and archaea, and is an integral part of the body influencing its overall homeostasis. Emerging evidence highlights dysbiosis of the archaeome and mycobiome to influence the overall intestinal microbiota composition in health and disease, including IBD, although they remain some of the least understood components of the gut microbiota. Nonetheless, their ability to directly impact the other commensals, or the host, reasonably makes them important contributors to either the maintenance of the mucosal tissue physiology or to chronic intestinal inflammation development. Therefore, the full understanding of the archaeome and mycobiome dysbiosis during IBD pathogenesis may pave the way to the discovery of novel mechanisms, finally providing innovative therapeutic targets that can soon implement the currently available treatments for IBD patients.

Evolutionary characteristics and phylogeny of cotton chloroplast tRNAs

The novel structural variations were identified in cotton chloroplast tRNAs and gene loss events were more obvious than duplications in chloroplast tRNAs. Transfer RNAs (tRNA) have long been believed an evolutionary-conserved molecular family, which play the key roles in the process of protein biosynthesis in plant life activities. In this study, we detected the evolutionary characteristics and phylogeny of chloroplast tRNAs in cotton plants, an economic and fibered important taxon in the world. We firstly annotated the chloroplast tRNAs of 27 Gossypium species to analyze their genetic composition, structural characteristics and evolution. Compared with the traditional view of evolutionary conservation of tRNA, some novel tRNA structural variations were identified in cotton plants. I.g., tRNA^Val-UAC and tRNA^Ile-GAU only contained one intron in the anti-condon loop region of tRNA secondary structure, respectively. In the variable region, some tRNAs contained a circle structure with a few nucleotides. Interestingly, the calculation result of free energy indicated that the variation of novel tRNAs contributed to the stability of tRNA structure. Phylogenetic analysis suggested that chloroplast tRNAs have evolved from multiple common ancestors, and the tRNA^Met seemed to be an ancestral tRNA, which can be duplicated and diversified to produce other tRNAs. The chloroplast tRNAs contained a group I intron in cotton plants, and the evolutionary analysis of introns indicated that group I intron of chloroplast tRNA originated from cyanobacteria. Analysis of gene duplication and loss events showed that gene loss events were more obvious than duplications in Gossypium chloroplast tRNAs. Additionally, we found that the rate of transition was higher than the ones of transversion in cotton chloroplast tRNAs. This study provided new insights into the structural characteristics and evolution of chloroplast tRNAs in cotton plants.

K-strategy species plays a pivotal role in the natural attenuation of petroleum hydrocarbon pollution in aquifers

The natural attenuation of petroleum hydrocarbons is inseparable from the action of microorganisms, while the degradation methods and ecological strategies of microorganisms in petroleum-contaminated aquifers are still under debate. In the present study, 16 S rRNA sequencing and quantitative real-time polymerase chain reaction were used to assess the potential microbial degradation of petroleum hydrocarbons, and the ecological strategy of microorganisms under petroleum stress was analyzed through a co-occurrence network. The results showed that the microbial community in sediments exhibit higher efficiency and stability and stronger ecological function than that in groundwater. Keystone species coordinated with the community to execute ecosystem processes and tended to choose a K-strategy to survive, with the aquifer sediment being the main site of petroleum hydrocarbon degradation. Under natural conditions, the presence of petroleum hydrocarbons at concentrations higher than 126 μg kg^-1 and 5557 μg kg^-1 was not conducive to the microbial degradation of polycyclic aromatic hydrocarbons and alkanes, respectively. These results can be used as a reference for an enhanced bioremediation of contaminated groundwater. Overall, these findings provide support to managers for developing environmental management strategies.

Microbial generalist or specialist: Intraspecific variation and dormancy potential matter

Microbial generalists and specialists coexist in the soil environment while having distinctive impacts on microbial community dynamics. In microbial ecology, the underlying mechanisms as to why a species is a generalist or a specialist remain ambiguous. Herein, we collected soils across a national scale and identified bacterial generalists and specialists according to niche breadth at the species level (OTU level), and the single-nucleotide differences in each species were measured to investigate intraspecific variation (at zero-radius OTU level). Compared with that of the specialists, the intraspecific variation of the generalists was much higher, which ensured their wider niche breadth and lower variability. The higher asynchrony and different niche preferences of conspecific individuals and the higher dormancy potential within the generalists further contributed to their stability in varying environments. Besides, generalists were less controlled by environmental filtering, which was indicated by the stronger signature of stochastic processes in their assembly, and had higher diversification and transition rates that allowed them to adapt to environmental changes to a greater extent than specialists. Overall, this study provides a new comprehensive understanding of the rules of assembly and the evolutionary roles of bacterial generalists and specialists. It also highlights the importance of intraspecific variation and the dormancy potential in the stability of species.

Restoration, conservation and phytoplankton hysteresis

Phytoplankton growth depends not only upon external factors that are not strongly altered by the presence of phytoplankton, such as temperature, but also upon factors that are strongly influenced by activity of phytoplankton, including photosynthetically active radiation, and the availability of the macronutrients carbon, nitrogen, phosphorus and, for some, silicate. Since phytoplankton therefore modify, and to an extent create, their own habitats, established phytoplankton communities can show resistance and resilience to change, including managed changes in nutrient regimes. Phytoplankton blooms and community structures can be predicted from the overall biogeochemical setting and inputs, but restorations may be influenced by the physiological responses of established phytoplankton taxa to nutrient inputs, temperature, second-order changes in illumination and nutrient recycling. In this review we discuss the contributions of phytoplankton ecophysiology to biogeochemical hysteresis and possible effects on community composition in the face of management, conservation or remediation plans.

Honeybee Exposure to Veterinary Drugs: How Is the Gut Microbiota Affected?

Several studies have outlined that a balanced gut microbiota offers metabolic and protective functions supporting honeybee health and performance. The present work contributes to increasing knowledge on the impact on the honeybee gut microbiota of the three most common veterinary drugs (oxytetracycline, sulfonamides, and tylosin). The study was designed with a semi-field approach in micro-hives containing about 500 honeybees. Micro-hives were located in an incubator during the day and moved outdoors in the late afternoon, considering the restrictions on the use of antibiotics in the open field but allowing a certain freedom to honeybees; 6 replicates were considered for each treatment. The absolute abundance of the major gut microbial taxa in newly eclosed individuals was studied with qPCR and next-generation sequencing. Antimicrobial resistance genes for the target antibiotics were also monitored using a qPCR approach. The results showed that the total amount of gut bacteria was not altered by antibiotic treatment, but qualitative variations were observed. Tylosin treatment determined a significant decrease of α- and β-diversity indices and a strong depletion of the rectum population (lactobacilli and bifidobacteria) while favoring the ileum microorganisms (Gilliamella, Snodgrassella, and Frischella spp.). Major changes were also observed in honeybees treated with sulfonamides, with a decrease in Bartonella and Frischella core taxa and an increase of Bombilactobacillus spp. and Snodgrassella spp. The present study also shows an important effect of tetracycline that is focused on specific taxa with minor impact on alfa and beta diversity. Monitoring of antibiotic resistance genes confirmed that honeybees represent a great reservoir of tetracycline resistance genes. Tetracycline and sulfonamides resistance genes tended to increase in the gut microbiota population upon antibiotic administration. IMPORTANCE This study investigates the impact of the three most widely used antibiotics in the beekeeping sector (oxytetracycline, tylosin, and sulfonamides) on the honeybee gut microbiota and on the spread of antibiotic resistance genes. The research represents an advance to the present literature, considering that the tylosin and sulfonamides effects on the gut microbiota have never been studied. Another original aspect lies in the experimental approach used, as the study looks at the impact of veterinary drugs and feed supplements 24 days after the beginning of the administration, in order to explore perturbations in newly eclosed honeybees, instead of the same treated honeybee generation. Moreover, the study was not performed with cage tests but in micro-hives, thus achieving conditions closer to real hives. The study reaches the conclusion that the most common veterinary drugs determine changes in some core microbiota members and that incidence of resistance genes for tetracycline and sulfonamides increases following antibiotic treatment.

Characterizing Microbiomes via Sequencing of Marker Loci: Techniques To Improve Throughput, Account for Cross-Contamination, and Reduce Cost

New approaches to characterizing microbiomes via high-throughput sequencing provide impressive gains in efficiency and cost reduction compared to approaches that were standard just a few years ago. However, the speed of method development has been such that staying abreast of the latest technological advances is challenging. Moreover, shifting laboratory protocols to include new methods can be expensive and time consuming. To facilitate adoption of new techniques, we provide a guide and review of recent advances that are relevant for single-locus sequence-based study of microbiomes—from extraction to library preparation—including a primer regarding the use of liquid-handling automation in small-scale academic settings. Additionally, we describe several amendments to published techniques to improve throughput, track contamination, and reduce cost. Notably, we suggest adding synthetic DNA molecules to each sample during nucleic acid extraction, thus providing a method of documenting incidences of cross-contamination. We also describe a dual-indexing scheme for Illumina sequencers that allows multiplexing of many thousands of samples with minimal PhiX input. Collectively, the techniques that we describe demonstrate that laboratory technology need not impose strict limitations on the scale of molecular microbial ecology studies.

IMPORTANCE New methods to characterize microbiomes reduce technology-imposed limitations to study design, but many new approaches have not been widely adopted. Here, we present techniques to increase throughput and reduce contamination alongside a thorough review of current best practices.

Microbiome Aggregated Traits and Assembly Are More Sensitive to Soil Management than Diversity

How soil is managed, particularly for agriculture, exerts stresses upon soil microbiomes, resulting in altered community structures and functional states. Understanding how soil microbiomes respond to combined stresses is important for predicting system performance under different land use scenarios, aids in identification of the most environmentally benign managements, and provides insight into how system function can be recovered in degraded soils. We use a long-established field experiment to study the effects of combined chronic (press) disturbance of the magnitude of organic carbon inputs with acute (pulse) effects of physical disturbance by tillage and chemical disturbance due to inorganic fertilization and pesticide application. We show that because of the variety of ways it can be assessed, biodiversity-here based on microbial small subunit rRNA gene phylotypes-does not provide a consistent view of community change. In contrast, aggregated traits associated with soil microbiomes indicate general loss of function, measured as a reduction of average genome lengths, associated with chronic reduction of organic inputs in arable or bare fallow soils and altered growth strategies associated with rRNA operon copy number in prokaryotes, as well as a switch to pathogenicity in fungal communities. In addition, pulse disturbance by soil tillage is associated with an increased influence of stochastic processes upon prokaryote community assembly, but fungicide used in arable soils results in niche assembly of fungal communities compared to untilled grassland. Overall, bacteria, archaea, and fungi do not share a common response to land management change, and estimates of biodiversity do not capture important facets of community adaptation to stresses adequately. IMPORTANCE Changes in soil microbiome diversity and function brought about by land management are predicted to influence a range of environmental services provided by soil, including provision of food and clean water. However, opportunities to compare the long-term effects of combinations of stresses imposed by different management approaches are limited. We exploit a globally unique 50-year field experiment, demonstrating that soil management practices alter microbiome diversity, community traits, and assembly. Grassland soil microbiomes are dominated by fewer-but phylogenetically more diverse-prokaryote phylotypes which sustain larger genomes than microbiomes in arable or bare fallow soil maintained free of plants. Dominant fungi in grassland soils are less phylogenetically diverse than those in arable or fallow soils. Soil tillage increases stochastic processes in microbiome assembly: this, combined with reduced plant biomass, presents opportunities for organisms with a capacity for pathogenesis to become established in stressed soils.

Temperature sensitivity of SOM decomposition is linked with a K-selected microbial community

Temperature sensitivity (Q₁₀ ) of soil organic matter (SOM) decomposition is a crucial parameter to predict the fate of soil carbon (C) under global warming. Nonetheless, the response pattern of Q₁₀ to continuous warming and the underlying mechanisms are still under debate, especially considering the complex interactions between Q₁₀ , SOM quality, and soil microorganisms. We examined the Q₁₀ of SOM decomposition across a mean annual temperature (MAT) gradient from -1.9 to 5.1°C in temperate mixed forest ecosystems in parallel with SOM quality and bioavailability, microbial taxonomic composition, and functional genes responsible for organic carbon decomposition. Within this temperature gradient of 7.0°C, the Q₁₀ values increased with MAT, but decreased with SOM bioavailability. The Q₁₀ values increased with the prevalence of K-strategy of soil microbial community, which was characterized by: (i) high ratios of oligotrophic to copiotrophic taxa, (ii) ectomycorrhizal to saprotrophic fungi, (iii) functional genes responsible for degradation of recalcitrant to that of labile C, and (iv) low average 16S rRNA operon copy number. Because the recalcitrant organic matter was mainly utilized by the K-strategists, these findings independently support the carbon quality-temperature theory from the perspective of microbial taxonomic composition and functions. A year-long incubation experiment was performed to determine the response of labile and recalcitrant C pools to warming based on the two-pool model. The decomposition of recalcitrant SOM was more sensitive to increased temperature in southern warm regions, which might attribute to the dominance of K-selected microbial communities. It implies that climate warming would mobilize the larger recalcitrant pools in warm regions, exacerbating the positive feedback between increased MAT and CO₂ efflux. This is the first attempt to link temperature sensitivity of SOM decomposition with microbial eco-strategies by incorporating the genetic information and disentangling the complex relationship between Q₁₀ and soil microorganisms.

Honeybees Exposure to Natural Feed Additives: How Is the Gut Microbiota Affected?

The role of a balanced gut microbiota to maintain health and prevent diseases is largely established in humans and livestock. Conversely, in honeybees, studies on gut microbiota perturbations by external factors have started only recently. Natural methods alternative to chemical products to preserve honeybee health have been proposed, but their effect on the gut microbiota has not been examined in detail. This study aims to investigate the effect of the administration of a bacterial mixture of bifidobacteria and Lactobacillaceae and a commercial product HiveAlive^TM on honeybee gut microbiota. The study was developed in 18 hives of about 2500 bees, with six replicates for each experimental condition for a total of three experimental groups. The absolute abundance of main microbial taxa was studied using qPCR and NGS. The results showed that the majority of the administered strains were detected in the gut. On the whole, great perturbations upon the administration of the bacterial mixture and the plant-based commercial product were not observed in the gut microbiota. Significant variations with respect to the untreated control were only observed for Snodgrassella sp. for the bacterial mixture, Bartonella sp. in HiveAlive^TM and Bombilactobacillus sp. for both. Therefore, the studied approaches are respectful of the honeybee microbiota composition, conceivably without compromising the bee nutritional, social and ecological functions.

Comparative evaluation of Vibrio delineation methodologies in post-genomic era

Vibrios are widespread in both marine and coastal water environments and are recognized as one of the most important prokaryotic pathogens because they may potentially threaten the health of both aquacultures and human beings. However, owing to highly similar physiological and biochemical properties, accurate classification and identification of Vibrio strains remains challenging. This hampers further research on the physiology, pathogeny, genomics, epidemics, and ecology of vibrios. Here, we comparatively evaluated multiple approaches including 16S rRNA gene identity, average nucleotide identity (ANI), gene content similarity and mutilocus sequence analysis (MLSA) to investigate their ability in delineating Vibrio strains. In addition, we also evaluated the possibility of applying bacterial prophages in classifying and identifying Vibrio strains. Our results showed that MLSA outperformed other methods in discriminating Vibrio species, suggesting that the other four approaches should be used with cautions in Vibrio delineation. Interestingly, we also found that prophages identified in Vibrio strains were highly specific at strain- and species-level, suggesting that prophages held the potential to be used for microbial species, sub-species, and strain-level identifications. This study is expected to provide valuable insights into the taxonomic identification and classification of complex microbial groups in the post-genomic era.

Neonicotinoids in the agroecosystem: In-field long-term assessment on honeybee colony strength and microbiome

Bees can be severely affected by various plant protection products (PPP). Among these, neonicotinoid insecticides are of concern as they have been shown to be responsible for extensive honeybee colonies death when released into the environment. Also, sublethal neonicotinoid doses contaminating single honeybees and their colonies (e.g. through contaminated pollen) are responsible for honeybees physiological alterations with probable implication also on microbiome functionality. Honeybees show symbiotic interactions with specific gut bacteria that can enhance the adult host performances. Among the known mechanisms, the modulation of the immune system, the degradation of recalcitrant secondary plant metabolites, pollen digestion, and hormonal signaling, are the most important functional benefits for the host honeybee. To date, few research efforts have aimed at revealing the impact of PPP on the gut microbial community of managed and wild honeybees. The majority of the existing literature relays on cage or semifield tests of short duration for research investigating neonicotinoids-gut microbiome interactions. This research wanted to unravel the impact of two neonicotinoids (i.e. imidacloprid and thiacloprid) in natural field conditions up to 5 weeks of exposure. A long-term impact of neonicotinoids on gut microbial community of honeybees was observed. The alterations affected several microbial genera and species such as Frischella spp., lactobacilli and bifidobacteria, whose shifting is implicated in intestinal dysbiosis. Long-term impact leading to dysbiosis was detected in case of exposure to imidacloprid, whereas thiacloprid exposure stimulated temporary dysbiosis. Moreover, the microbial diversity was significantly reduced in neonicotinoid-treated groups. Overall, the reported results support a compromised functionality of the gut microbial community, that might reflect a lower efficiency in the ecosystemic functionality of honeybees.

Development of an absolute quantification method for ribosomal RNA gene copy numbers per eukaryotic single cell by digital PCR

Recent increase of Harmful Algal Blooms (HAB) causes world-wide ecological, economical, and health issues, and more attention is paid to frequent coastal monitoring for the early detection of HAB species to prevent or reduce such impacts. Use of molecular tools in addition to traditional microscopy-based observation has become one of the promising methodologies for coastal monitoring. However, as ribosomal RNA (rRNA) genes are commonly targeted in molecular studies, variability in the rRNA gene copy number within and between species must be considered to provide quantitative information in quantitative PCR (qPCR), digital PCR (dPCR), and metabarcoding analyses. Currently, this information is only available for a limited number of species. The present study utilized a dPCR technology to quantify copy numbers of rRNA genes per single cell in 16 phytoplankton species, the majority of which are toxin-producers, using a newly developed universal primer set accompanied by a labeled probe with a fluorophore and a double-quencher. In silico PCR using the newly developed primers allowed the detection of taxa from 8 supergroups, demonstrating universality and broad coverage of the primer set. Chelex buffer was found to be suitable for DNA extraction to obtain DNA fragments with suitable size to avoid underestimation of the copy numbers. The study successfully demonstrated the first comparison of absolute quantification of 18S rRNA copy numbers per cell from 16 phytoplankton species by the dPCR technology.

Lipoptena fortisetosa as a vector of Bartonella bacteria in Japanese sika deer (Cervus nippon)

Background

Two species of deer ked (Lipoptena cervi and L. mazamae) have been identified as vectors of Bartonella bacteria in cervids in Europe and the USA. In an earlier study we showed that Japanese sika deer (Cervus nippon) harbor three Bartonella species, namely B. capreoli (lineage A) and two novel Bartonella species (lineages B and C); however, there is currently no information on the vector of Bartonella bacteria in sika deer. The aim of this study was to clarify potential vectors of Bartonella in Japanese sika deer.

Methods

Thirty-eight wingless deer keds (L. fortisetosa) and 36 ticks (Haemaphysalis and Ixodes species) were collected from sika deer. The prevalence of Bartonella in the arthropods was evaluated by real-time PCR targeting the 16S−23S internal transcribed spacer (ITS) and by culture of the organisms. The total number of Bartonella bacteria were quantified using real-time PCR. The distribution of Bartonella bacteria in deer ked organs was examined by immunofluorescence analysis. The relationship of Bartonella strains isolated from sika deer and arthropods were examined by a phylogenetic analysis based on concatenated sequences of the gltA, rpoB, ftsZ, and ribC genes, followed by a BLAST search for gltA and rpoB.

Results

Bartonella prevalence in deer keds was 87.9% by real-time PCR and 51.5% in culture and that in the ticks was 8.3% by real-time PCR and 2.8% in culture. The mean number of Bartonella bacteria per ked was calculated to be 9.2 × 10⁵ cells. Bartonella aggregates were localized in the midgut of the keds. The phylogenetic analysis and BLAST search showed that both the host deer and the keds harbored two Bartonella species (lineages B and C), while B. capreoli (lineage A) was not detected in the keds. Two novel Bartonella species (lineages D and E) were isolated from one ked.

Conclusions

Lipoptena fortisetosa likely serves as a vector of at least two Bartonella species (lineages B and C), whereas ticks do not seem to play a significant role in the transmission of Bartonella between sika deer based on the lower detection rates of Bartonella in ticks compared to keds. Bartonella species in lineages D and E appear to be L. fortisetosa-specific strains.

The host-associated archaeome

Host-associated microbial communities have an important role in shaping the health and fitness of plants and animals. Most studies have focused on the bacterial, fungal or viral communities, but often the archaeal component has been neglected. The archaeal community, the so-called archaeome, is now increasingly recognized as an important component of host-associated microbiomes. It is composed of various lineages, including mainly Methanobacteriales and Methanomassiliicoccales (Euryarchaeota), as well as representatives of the Thaumarchaeota. Host-archaeome interactions have mostly been delineated from methanogenic archaea in the gastrointestinal tract, where they contribute to substantial methane production and are potentially also involved in disease-relevant processes. In this Review, we discuss the diversity and potential roles of the archaea associated with protists, plants and animals. We also present the current understanding of the archaeome in humans, the specific adaptations involved in interaction with the resident microbial community as well as with the host, and the roles of the archaeome in both health and disease.

Key internal factors leading to the variability in CO2 fixation efficiency of different sulfur-oxidizing bacteria during autotrophic cultivation

Variability in the apparent CO₂ fixation yield of four aerobic sulfur-oxidizing bacteria (Halothiobacillus neapolitanus DSM 15147, Thiobacillus thioparus DSM 505, Thiomonas intermedia DSM 18155, and Starkeya novella DSM 506) in autotrophic culturing was studied, and mutual effects of key intrinsic factors on CO₂ fixation were explored. DSM 15147 and DSM 505 exhibited much higher CO₂ fixation yields than DSM 18155 and DSM 506. The differences in CO₂ fixation yield were determined not only by cbb gene transcription, but also by cell synthesis rate, which was determined by rRNA gene copy number; the rRNA gene copy number had a more significant effect than cbb gene transcription on the apparent CO₂ fixation yield. Moreover, accumulation of EDOC was observed in all four strains during chemoautotrophic cultivation, and the proportion of EDOC accounting for total fixed organic carbon (TOC; EDOC/TOC ratio) was much higher in DSM 18155 and DSM 506 than in DSM 15147 and DSM 505. The accumulation of EDOC led to a significant decrease in the cbb gene transcription efficiency during cultivation, and a further feedback inhibitory effect on CO₂ fixation.

Regeneration of unconventional natural gas by methanogens co-existing with sulfate-reducing prokaryotes in deep shale wells in China

Biogenic methane in shallow shale reservoirs has been proven to contribute to economic recovery of unconventional natural gas. However, whether the microbes inhabiting the deeper shale reservoirs at an average depth of 4.1 km and even co-occurring with sulfate-reducing prokaryote (SRP) have the potential to produce biomethane is still unclear. Stable isotopic technique with culture-dependent and independent approaches were employed to investigate the microbial and functional diversity related to methanogenic pathways and explore the relationship between SRP and methanogens in the shales in the Sichuan Basin, China. Although stable isotopic ratios of the gas implied a thermogenic origin for methane, the decreased trend of stable carbon and hydrogen isotope value provided clues for increasing microbial activities along with sustained gas production in these wells. These deep shale-gas wells harbored high abundance of methanogens (17.2%) with ability of utilizing various substrates for methanogenesis, which co-existed with SRP (6.7%). All genes required for performing methylotrophic, hydrogenotrophic and acetoclastic methanogenesis were present. Methane production experiments of produced water, with and without additional available substrates for methanogens, further confirmed biomethane production via all three methanogenic pathways. Statistical analysis and incubation tests revealed the partnership between SRP and methanogens under in situ sulfate concentration (~ 9 mg/L). These results suggest that biomethane could be produced with more flexible stimulation strategies for unconventional natural gas recovery even at the higher depths and at the presence of SRP.

Complete genome sequence of the biocontrol agent Serratia marcescens strain N4-5 uncovers an assembly artefact

Serratia marcescens are gram-negative bacteria found in several environmental niches, including the plant rhizosphere and patients in hospitals. Here, we present the genome of Serratia marcescens strain N4-5 (=NRRL B-65519), which has a size of 5,074,473 bp (664-fold coverage) and contains 4840 protein coding genes, 21 RNA genes, and an average G + C content of 59.7%. N4-5 harbours a plasmid of 11,089 bp and 43.5% G + C content that encodes six unique CDS repeated 2.5× times totalling 13 CDS. Our genome assembly and manual curation uncovered the insertion of two extra copies of the 5S rRNA gene in the assembled sequence, which was confirmed by PCR and Sanger sequencing to be a misassembly. This artefact was subsequently removed from the final assembly. The occurrence of extra copies of the 5S rRNA gene was also observed in most complete genomes of Serratia spp. deposited in public databases in our comparative analysis. These elements, which also occur naturally, can easily be confused with true genetic variation. Efforts to discover and correct assembly artefacts should be made in order to generate genome sequences that represent the biological truth underlying the studied organism. We present the genome of N4-5 and discuss genes potentially involved in biological control activity against plant pathogens and also the possible mechanisms responsible for the artefact we observed in our initial assembly. This report raises awareness about the extra copies of the 5S rRNA gene in sequenced bacterial genomes as they may represent misassemblies and therefore should be verified experimentally.

Implementation and integration of microbial source tracking in a river watershed monitoring plan

Fecal pollution of water bodies poses a serious threat for public health and ecosystems. Microbial source tracking (MST) is used to track the source of this pollution facilitating better management of pollution at the source. In this study we tested 12 MST markers to track human, ruminant, sheep, horse, pig and gull pollution to assess their usefulness as an effective management tool of water quality. First, the potential of the selected markers to track the source was evaluated using fresh fecal samples. Subsequently, we evaluated their performance in a catchment with different impacts, considering land use and environmental conditions. All MST markers showed high sensitivity and specificity, although none achieved 100% for both. Although some of the MST markers were detected in hosts other than the intended ones, their abundance in the target group was always several orders of magnitude higher than in the non-target hosts, demonstrating their suitability to distinguish between sources of pollution. The MST analysis matched the land use in the watershed allowing an accurate assessment of the main sources of pollution, in this case mainly human and ruminant pollution. Correlating environmental parameters including temperature and rainfall with MST markers provided insight into the dynamics of the pollution in the catchment. The levels of the human marker showed a significant negative correlation with rainfall in human polluted areas suggesting a dilution of the pollution, whereas at agricultural areas the ruminant marker increased with rainfall. There were no seasonal differences in the levels of human marker, indicating human pollution as a constant pressure throughout the year, whereas the levels of the ruminant marker was influenced by the seasons, being more abundant in summer and autumn. MST analysis integrated with land use and environmental data can improve the management of fecal polluted areas and set up best practice.

The quest for absolute abundance: The use of internal standards for DNA-based community ecology

To characterize microbiomes and other ecological assemblages, ecologists routinely sequence and compare loci that differ among focal taxa. Counts of these sequences convey information regarding the occurrence and relative abundances of taxa, but provide no direct measure of their absolute abundances, due to the technical limitations of the sequencing process. The relative abundances in compositional data are inherently constrained and difficult to interpret. The incorporation of internal standards (ISDs; colloquially referred to as 'spike-ins') into DNA pools can ameliorate the problems posed by relative abundance data and allow absolute abundances to be approximated. Unfortunately, many laboratory and sampling biases cause ISDs to underperform or fail. Here, we discuss how careful deployment of ISDs can avoid these complications and be an integral component of well-designed studies seeking to characterize ecological assemblages via sequencing of DNA.

Plant growth promoting rhizobacteria isolated from halophytes and drought-tolerant plants: genomic characterisation and exploration of phyto-beneficial traits

Plant growth promoting rhizobacteria (PGPR) are able to provide cross-protection against multiple stress factors and facilitate growth of their plant symbionts in many ways. The aim of this study was to isolate and characterize rhizobacterial strains under natural conditions, associated with naturally occurring representatives of wild plant species and a local tomato cultivar, growing in differently stressed Mediterranean ecosystems. A total of 85 morphologically different rhizospheric strains were isolated; twenty-five exhibited multiple in vitro PGP-associated traits, including phosphate solubilization, indole-3-acetic acid production, and 1-aminocyclopropane-1-carboxylate deaminase activity. Whole genome analysis was applied to eight selected strains for their PGP potential and assigned seven strains to Gammaproteobacteria, and one to Bacteroidetes. The genomes harboured numerous genes involved in plant growth promotion and stress regulation. They also support the notion that the presence of gene clusters with potential PGP functions is affirmative but not necessary for a strain to promote plant growth under abiotic stress conditions. The selected strains were further tested for their ability to stimulate growth under stress. This initial screening led to the identification of some strains as potential PGPR for increasing crop production in a sustainable manner.

A prospective longitudinal study on the microbiota composition in amyotrophic lateral sclerosis

BACKGROUND

A connection between amyotrophic lateral sclerosis (ALS) and altered gut microbiota composition has previously been reported in animal models. This work is the first prospective longitudinal study addressing the microbiota composition in ALS patients and the impact of a probiotic supplementation on the gut microbiota and disease progression.

METHODS

Fifty patients and 50 matched controls were enrolled. The microbial profile of stool samples from patients and controls was analyzed via PCR-Denaturing Gradient Gel Electrophoresis, and the main microbial groups quantified via qPCR. The whole microbiota was then analyzed via next generation sequencing after amplification of the V3-V4 region of 16S rDNA. Patients were then randomized to receive probiotic treatment or placebo and followed up for 6 months with ALSFRS-R, BMI, and FVC%.

RESULTS

The results demonstrate that the gut microbiota of ALS patients is characterized by some differences with respect to controls, regardless of the disability degree. Moreover, the gut microbiota composition changes during the course of the disease as demonstrated by the significant decrease in the number of observed operational taxonomic unit during the follow-up. Interestingly, an unbalance between potentially protective microbial groups, such as Bacteroidetes, and other with potential neurotoxic or pro-inflammatory activity, such as Cyanobacteria, has been shown. The 6-month probiotic treatment influenced the gut microbial composition; however, it did not bring the biodiversity of intestinal microbiota of patients closer to that of control subjects and no influence on the progression of the disease measured by ALSFRS-R was demonstrated.

CONCLUSIONS

Our study poses the bases for larger clinical studies to characterize the microbiota changes as a novel ALS biomarker and to test new microbial strategy to ameliorate the health status of the gut.

TRIAL REGISTRATION

CE 107/14, approved by the Ethics Committee of the "Maggiore della Carità" University Hospital, Italy.

Ribosome engineering reveals the importance of 5S rRNA autonomy for ribosome assembly

5S rRNA is an indispensable component of cytoplasmic ribosomes in all species. The functions of 5S rRNA and the reasons for its evolutionary preservation as an independent molecule remain unclear. Here we used ribosome engineering to investigate whether 5S rRNA autonomy is critical for ribosome function and cell survival. By linking circularly permutated 5S rRNA with 23S rRNA we generated a bacterial strain devoid of free 5S rRNA. Viability of the engineered cells demonstrates that autonomous 5S rRNA is dispensable for cell growth under standard conditions and is unlikely to have essential functions outside the ribosome. The fully assembled ribosomes carrying 23S-5S rRNA are highly active in translation. However, the engineered cells accumulate aberrant 50S subunits unable to form stable 70S ribosomes. Cryo-EM analysis revealed a malformed peptidyl transferase center in the misassembled 50S subunits. Our results argue that the autonomy of 5S rRNA is preserved due to its role in ribosome biogenesis.

Temporal variations in bacterial community diversity and composition throughout intensive care unit renovations

BACKGROUND

Inanimate surfaces within a hospital serve as a reservoir of microbial life that may colonize patients and ultimately result in healthcare associated infections (HAIs). Critically ill patients in intensive care units (ICUs) are particularly vulnerable to HAIs. Little is known about how the microbiome of the ICU is established or what factors influence its evolution over time. A unique opportunity to bridge the knowledge gap into how the ICU microbiome evolves emerged in our health system, where we were able to characterize microbial communities in an established hospital ICU prior to closing for renovations, during renovations, and then after re-opening.

RESULTS

We collected swab specimens from ICU bedrails, computer keyboards, and sinks longitudinally at each renovation stage, and analyzed the bacterial compositions on these surfaces by 16S rRNA gene sequencing. Specimens collected before ICU closure had the greatest alpha diversity, while specimens collected after the ICU had been closed for over 300 days had the least. We sampled the ICU during the 45 days after re-opening; however, within that time frame, the alpha diversity never reached pre-closure levels. There were clear and significant differences in microbiota compositions at each renovation stage, which was driven by environmental bacteria after closure and human-associated bacteria after re-opening and before closure.

CONCLUSIONS

Overall, we identified significant differences in microbiota diversity and community composition at each renovation stage. These data help to decipher the evolution of the microbiome in the most critical part of the hospital and demonstrate the significant impacts that microbiota from patients and staff have on the evolution of ICU surfaces. Video Abstract.

Microbial Fe(II) oxidation by Sideroxydans lithotrophicus ES-1 in the presence of Schlöppnerbrunnen fen-derived humic acids

Controlled laboratory experiments were combined with field measurements to better understand the interactions between dissolved organic matter (DOM) and reduced iron in organic-rich peatlands. Addition of peat-derived humic acid extract (HA) to Sideroxydans lithotrophicus ES-1 liquid cultures led to higher cell numbers and up to 1.4 times higher Fe(II) oxidation rates compared to chemical controls. This effect was positively correlated with increasing HA concentrations. Similar Fe(III) (oxyhydr)oxide mineralogies were formed both abiotically and biotically irrespective of HA amendment, but minerals formed in the presence of ES-1 and HA were smaller. ES-1 growth with HA promoted aggregation of Fe(III) products in agarose-stabilized gradient tubes as shown by voltammetric profiling. In situ voltammetry in an acidic, iron-rich peatland revealed a gap between oxygen penetration and iron reduction that may reflect active Fe(II)-oxidizing microorganisms. The highest abundance of Fe(II) oxidizers Sideroxydans (4.9 × 107 gene copies gww-1) and Gallionella (1.5 × 107 gene copies gww-1) in the upper peat layer coincided with small-sized minerals resembling nanoparticulate ferrihydrite or goethite. Our results suggest that microbially mediated Fe(II) oxidation dominates in the presence of DOM leading to the formation of nano-sized biogenic Fe(III) (oxyhydr)oxides that might be readily bioavailable and likely important to iron and carbon cycling.

A bird's-eye view of phylosymbiosis: weak signatures of phylosymbiosis among all 15 species of cranes

In numerous animal clades, the evolutionary history of host species drives patterns of gut microbial community structure, resulting in more divergent microbiota with increasing phylogenetic distance between hosts. This phenomenon, termed phylosymbiosis, has been observed in diverse evolutionary lineages, but has been difficult to detect in birds. Previous tests of phylosymbiosis among birds have been conducted using wild individuals, and thus interspecific differences in diet and environment may have masked a phylogenetic signal. Therefore, we tested for phylosymbiosis among all 15 species of cranes (family Gruidae) housed in the same captive environment and maintained on identical diets. 16S rRNA sequencing revealed that crane species harbour distinct gut microbiota. Overall, we detected marginally significant patterns of phylosymbiosis, the strength of which was increased when including the estimates of absolute microbial abundance (rather than relative abundance) derived from microbial densities determined by flow cytometry. Using this approach, we detected the statistically significant signatures of phylosymbiosis only after removing male cranes from our analysis, suggesting that using mixed-sex animal cohorts may prevent the detection of phylosymbiosis. Though weak compared with mammals (and especially insects), these results provide evidence of phylosymbiosis in birds. We discuss the potential differences between birds and mammals, such as transmission routes and host filtering, that may underlie the differences in the strength of phylosymbiosis.

Genome Insights of the Plant-Growth Promoting Bacterium Cronobacter muytjensii JZ38 With Volatile-Mediated Antagonistic Activity Against Phytophthora infestans

Salinity stress is a major challenge to agricultural productivity and global food security in light of a dramatic increase of human population and climate change. Plant growth promoting bacteria can be used as an additional solution to traditional crop breeding and genetic engineering. In the present work, the induction of plant salt tolerance by the desert plant endophyte Cronobacter sp. JZ38 was examined on the model plant Arabidopsis thaliana using different inoculation methods. JZ38 promoted plant growth under salinity stress via contact and emission of volatile compounds. Based on the 16S rRNA and whole genome phylogenetic analysis, fatty acid analysis and phenotypic identification, JZ38 was identified as Cronobacter muytjensii and clearly separated and differentiated from the pathogenic C. sakazakii. Full genome sequencing showed that JZ38 is composed of one chromosome and two plasmids. Bioinformatic analysis and bioassays revealed that JZ38 can grow under a range of abiotic stresses. JZ38 interaction with plants is correlated with an extensive set of genes involved in chemotaxis and motility. The presence of genes for plant nutrient acquisition and phytohormone production could explain the ability of JZ38 to colonize plants and sustain plant growth under stress conditions. Gas chromatography-mass spectrometry analysis of volatiles produced by JZ38 revealed the emission of indole and different sulfur volatile compounds that may play a role in contactless plant growth promotion and antagonistic activity against pathogenic microbes. Indeed, JZ38 was able to inhibit the growth of two strains of the phytopathogenic oomycete Phytophthora infestans via volatile emission. Genetic, transcriptomic and metabolomics analyses, combined with more in vitro assays will provide a better understanding the highlighted genes' involvement in JZ38's functional potential and its interaction with plants. Nevertheless, these results provide insight into the bioactivity of C. muytjensii JZ38 as a multi-stress tolerance promoting bacterium with a potential use in agriculture.

Bacterial and Eukaryotic Small-Subunit Amplicon Data Do Not Provide a Quantitative Picture of Microbial Communities, but They Are Reliable in the Context of Ecological Interpretations

High-throughput sequencing (HTS) of amplified fragments of rRNA genes provides unprecedented insight into the diversity of prokaryotic and eukaryotic microorganisms. Unfortunately, HTS data are prone to quantitative biases, which may lead to an erroneous picture of microbial community composition and thwart efforts to advance its understanding. These concerns motivated us to investigate how accurately HTS data characterize the variability of microbial communities, the relative abundances of specific phylotypes, and their relationships with environmental factors in comparison to an established microscopy-based method. We compared results obtained by HTS and catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) from three independent aquatic time series for both prokaryotic and eukaryotic microorganisms (almost 900 data points, the largest obtained with both methods so far). HTS and CARD-FISH data disagree with regard to relative abundances of bacterial and eukaryotic phylotypes but identify similar environmental drivers shaping bacterial and eukaryotic communities.

High-throughput sequencing (HTS) of gene amplicons is a preferred method of assessing microbial community composition, because it rapidly provides information from a large number of samples at high taxonomic resolution and low costs. However, mock community studies show that HTS data poorly reflect the actual relative abundances of individual phylotypes, casting doubt on the reliability of subsequent statistical analysis and data interpretation. We investigated how accurately HTS data reflect the variability of bacterial and eukaryotic community composition and their relationship with environmental factors in natural samples. For this, we compared results of HTS from three independent aquatic time series (n = 883) with those from an established, quantitative microscopic method (catalyzed reporter deposition-fluorescence in situ hybridization [CARD-FISH]). Relative abundances obtained by CARD-FISH and HTS disagreed for most bacterial and eukaryotic phylotypes. Nevertheless, the two methods identified the same environmental drivers to shape bacterial and eukaryotic communities. Our results show that amplicon data do provide reliable information for their ecological interpretations. Yet, when studying specific phylogenetic groups, it is advisable to combine HTS with quantification using microscopy and/or the addition of internal standards.

IMPORTANCE High-throughput sequencing (HTS) of amplified fragments of rRNA genes provides unprecedented insight into the diversity of prokaryotic and eukaryotic microorganisms. Unfortunately, HTS data are prone to quantitative biases, which may lead to an erroneous picture of microbial community composition and thwart efforts to advance its understanding. These concerns motivated us to investigate how accurately HTS data characterize the variability of microbial communities, the relative abundances of specific phylotypes, and their relationships with environmental factors in comparison to an established microscopy-based method. We compared results obtained by HTS and catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) from three independent aquatic time series for both prokaryotic and eukaryotic microorganisms (almost 900 data points, the largest obtained with both methods so far). HTS and CARD-FISH data disagree with regard to relative abundances of bacterial and eukaryotic phylotypes but identify similar environmental drivers shaping bacterial and eukaryotic communities.

Promiscuous Ribozymes and Their Proposed Role in Prebiotic Evolution

The ability of enzymes, including ribozymes, to catalyze side reactions is believed to be essential to the evolution of novel biochemical activities. It has been speculated that the earliest ribozymes, whose emergence marked the origin of life, were low in activity but high in promiscuity, and that these early ribozymes gave rise to specialized descendants with higher activity and specificity. Here, we review the concepts related to promiscuity and examine several cases of highly promiscuous ribozymes. We consider the evidence bearing on the question of whether de novo ribozymes would be quantitatively more promiscuous than later evolved ribozymes or protein enzymes. We suggest that while de novo ribozymes appear to be promiscuous in general, they are not obviously more promiscuous than more highly evolved or active sequences. Promiscuity is a trait whose value would depend on selective pressures, even during prebiotic evolution.

Functional Genomics Insights Into the Pathogenicity, Habitat Fitness, and Mechanisms Modifying Plant Development of Rhodococcus sp. PBTS1 and PBTS2

Pistachio Bushy Top Syndrome (PBTS) is a recently emerged disease that has strongly impacted the pistachio industry in California, Arizona, and New Mexico. The disease is caused by two bacteria, designated PBTS1 that is related to Rhodococcus corynebacterioides and PBTS2 that belongs to the species R. fascians. Here, we assessed the pathogenic character of the causative agents and examined their chromosomal sequences to predict the presence of particular functions that might contribute to the observed co-occurrence and their effect on plant hosts. In diverse assays, we confirmed the pathogenicity of the strains on "UCB-1" pistachio rootstock and showed that they can also impact the development of tobacco species, but concurrently inconsistencies in the ability to induce symptoms were revealed. We additionally evidence that fas genes are present only in a subpopulation of pure PBTS1 and PBTS2 cultures after growth on synthetic media, that these genes are easily lost upon cultivation in rich media, and that they are enriched for in an in planta environment. Analysis of the chromosomal sequences indicated that PBTS1 and PBTS2 might have complementary activities that would support niche partitioning. Growth experiments showed that the nutrient utilization pattern of both PBTS bacteria was not identical, thus avoiding co-inhabitant competition. PBTS2 appeared to have the potential to positively affect the habitat fitness of PBTS1 by improving its resistance against increased concentrations of copper and penicillins. Finally, mining the chromosomes of PBTS1 and PBTS2 suggested that the bacteria could produce cytokinins, auxins, and plant growth-stimulating volatiles and that PBTS2 might interfere with ethylene levels, in support of their impact on plant development. Subsequent experimentation supported these in silico predictions. Altogether, our data provide an explanation for the observed pathogenic behavior and unveil part of the strategies used by PBTS1 and PBTS2 to interact with plants.

Time Series Genome-Centric Analysis Unveils Bacterial Response to Operational Disturbance in Activated Sludge

Understanding ecosystem response to disturbances and identifying the most critical traits for the maintenance of ecosystem functioning are important goals for microbial community ecology. In this study, we used 16S rRNA amplicon sequencing and metagenomics to investigate the assembly of bacterial populations in a full-scale municipal activated sludge wastewater treatment plant over a period of 3 years, including a 9-month period of disturbance characterized by short-term plant shutdowns. Following the reconstruction of 173 metagenome-assembled genomes, we assessed the functional potential, the number of rRNA gene operons, and the in situ growth rate of microorganisms present throughout the time series. Operational disturbances caused a significant decrease in bacteria with a single copy of the rRNA (rrn) operon. Despite moderate differences in resource availability, replication rates were distributed uniformly throughout time, with no differences between disturbed and stable periods. We suggest that the length of the growth lag phase, rather than the growth rate, is the primary driver of selection under disturbed conditions. Thus, the system could maintain its function in the face of disturbance by recruiting bacteria with the capacity to rapidly resume growth under unsteady operating conditions.IMPORTANCE Disturbance is a key determinant of community assembly and dynamics in natural and engineered ecosystems. Microbiome response to disturbance is thought to be influenced by bacterial growth traits and life history strategies. In this time series observational study, the response to disturbance of microbial communities in a full-scale activated sludge wastewater treatment plant was assessed by computing specific cellular traits of genomes retrieved from metagenomes. It was found that the genomes observed in disturbed periods have more copies of the rRNA operon than genomes observed in stable periods, whereas the in situ mean relative growth rates of bacteria present during stable and disturbed periods were indistinguishable. From these intriguing observations, we infer that the length of the lag phase might be a growth trait that affects the microbial response to disturbance. Further exploration of this hypothesis could contribute to better understanding of the adaptive response of microbiomes to unsteady environmental conditions.

Seasonal development of a coastal microbial mat

Growth and activity of coastal microbial mats is strongly seasonal. The development of these mats starts in early spring and fully maturate during late summer, where after growth ceases and subsequently the mat deteriorates by erosion and decomposition in winter. Here, the composition of the microbial community of three different mats developing along the tidal gradient of the North Sea beach of the Dutch barrier island Schiermonnikoog was analysed. The 16S ribosomal RNA molecules and the associated gene were sequenced in order to obtain the active (RNA) and resident (DNA) community members, respectively. Proteobacteria, Cyanobacteria, and Bacteroidetes dominated the mats during the whole year but considerable differences among these groups were found along the tidal gradient and seasonally when observed at a finer taxonomic resolution. Richness and diversity increased during the year starting from a pioneering community that is gradually succeeded by a more diverse climax community. The initial pioneers consisted of the cold-adapted photoautotrophic cyanobacterium Nodularia sp. and potential cold adapted members of the alphaproteobacterial Loktanella genus. These pioneers were succeeded by, amongst others, cyanobacteria belonging to the genera Leptolyngbya, Lyngbya, and Phormidium. At the upper littoral (Dune site), which was characterized by an extensive salt marsh vegetation, the mats contained a distinct bacterial community that potentially contribute to or benefit from plant decay. This study reports in detail on the seasonal changes and succession of these coastal microbial mat communities and discusses the potential forces that drive these changes.