A comparison of rpoB and 16S rRNA as markers in pyrosequencing studies of bacterial diversity

Broad-range amplification and sequencing of the rpoB gene: a novel assay for bacterial identification in clinical microbiology

The rpoB gene has been proposed as a promising phylogenetic marker for bacterial identification, providing theoretically improved species-level resolution compared to the 16S rRNA gene for a range of clinically important taxa. However, its utility in diagnostic microbiology has been limited by the lack of broad-range primers allowing for its amplification from most species with a single PCR assay. Here, we present an assay for broad-range partial amplification and Sanger sequencing of the rpoB gene. To reduce cross-reactivity and allow for rpoB amplification directly from patient samples, primers were based on the dual priming oligonucleotide principle. The resulting amplicon is ~550 base pairs in length and appropriate for species-level identification. Systematic in silico evaluation of a wide selection of taxa demonstrated improved resolution within multiple important genera, including Enterococcus, Fusobacterium, Mycobacterium, Streptococcus, and Staphylococcus species and several genera within the Enterobacteriaceae family. Broad-range rpoB amplification and Sanger sequencing of 115 bacterial isolates provided unambiguous species-level identification for 97 (84%) isolates, as compared to 57 (50%) using a clinical 16S rRNA gene assay. Several unresolved taxonomic matters disguised by the low resolution of the 16S rRNA gene were revealed using the rpoB gene. Using a collection of 33 clinical specimens harboring bacteria and assumed to contain high concentrations of human DNA, the rpoB assay identified the pathogen in 29 specimens (88%). Broad-range rpoB amplification and sequencing provides a promising tool for bacterial identification, improving discrimination between closely related species and making it amenable for use in culture-based and culture-independent diagnostic approaches.

Environmental DNA: The next chapter

Molecular tools are an indispensable part of ecology and biodiversity sciences and implemented across all biomes. About a decade ago, the use and implementation of environmental DNA (eDNA) to detect biodiversity signals extracted from environmental samples opened new avenues of research. Initial eDNA research focused on understanding population dynamics of target species. Its scope thereafter broadened, uncovering previously unrecorded biodiversity via metabarcoding in both well-studied and understudied ecosystems across all taxonomic groups. The application of eDNA rapidly became an established part of biodiversity research, and a research field by its own. Here, we revisit key expectations made in a land-mark special issue on eDNA in Molecular Ecology in 2012 to frame the development in six key areas: (1) sample collection, (2) primer development, (3) biomonitoring, (4) quantification, (5) behaviour of DNA in the environment and (6) reference database development. We pinpoint the success of eDNA, yet also discuss shortfalls and expectations not met, highlighting areas of research priority and identify the unexpected developments. In parallel, our retrospective couples a screening of the peer-reviewed literature with a survey of eDNA users including academics, end-users and commercial providers, in which we address the priority areas to focus research efforts to advance the field of eDNA. With the rapid and ever-increasing pace of new technical advances, the future of eDNA looks bright, yet successful applications and best practices must become more interdisciplinary to reach its full potential. Our retrospect gives the tools and expectations towards concretely moving the field forward.

A Comprehensive View of Food Microbiota: Introducing FoodMicrobionet v5

Amplicon-targeted metagenomics is now the standard approach for the study of the composition and dynamics of food microbial communities. Hundreds of papers on this subject have been published in scientific journals and the information is dispersed in a variety of sources, while raw sequences and their metadata are available in public repositories for some, but not all, of the published studies. A limited number of web resources and databases allow scientists to access this wealth of information but their level of annotation on studies and samples varies. Here, we report on the release of FoodMicrobionet v5, a comprehensive database of metataxonomic studies on bacterial and fungal communities of foods. The current version of the database includes 251 published studies (11 focusing on fungal microbiota, 230 on bacterial microbiota, and 10 providing data for both bacterial and fungal microbiota) and 14,035 samples with data on bacteria and 1114 samples with data on fungi. The new structure of the database is compatible with interactive apps and scripts developed for previous versions and allows scientists, R&D personnel in industries and regulators to access a wealth of information on food microbial communities.

Macroevolutionary Dynamics in Micro-organisms: Generalists Give Rise to Specialists Across Biomes in the Ubiquitous Bacterial Phylum Myxococcota

Prokaryotes dominate the Tree of Life, but our understanding of the macroevolutionary processes generating this diversity is still limited. Habitat transitions are thought to be a key driver of prokaryote diversity. However, relatively little is known about how prokaryotes successfully transition and persist across environments, and how these processes might vary between biomes and lineages. Here, we investigate biome transitions and specialization in natural populations of a focal bacterial phylum, the Myxococcota, sampled across a range of replicated soils and freshwater and marine sediments in Cornwall (UK). By targeted deep sequencing of the protein-coding gene rpoB, we found >2,000 unique Myxococcota lineages, with the majority (77%) classified as biome specialists and with only <5% of lineages distributed across the salt barrier. Discrete character evolution models revealed that specialists in one biome rarely transitioned into specialists in another biome. Instead, evolved generalism mediated transitions between biome specialists. State-dependent diversification models found variation in speciation rates across the tree, but this variation was independent of biome association or specialization. Our findings were robust to phylogenetic uncertainty, different levels of species delineation, and different assumed amounts of unsampled diversity resulting in an incomplete phylogeny. Overall, our results are consistent with a "jack-of-all-trades" tradeoff where generalists suffer a cost in any individual environment, resulting in rapid evolution of niche specialists and shed light on how bacteria could transition between biomes.

Comparison of MALDI-TOF mass spectrometry and 16S rDNA sequencing for identification of environmental bacteria: a case study of cave mussel-associated culturable microorganisms

Matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) is routinely used as a rapid and cost-effective method for pathogen identification in clinical settings. In comparison, its performance in other microbiological fields, such as environmental microbiology, is still being tested, although isolates of environmental microbes are essential for in-depth in vivo studies of their biology, including biotechnological applications. We investigated the applicability of MALDI-TOF MS for the identification of bacterial isolates from a highly oligotrophic environment - Dinaric Karst caves, which likely harbor specific microorganisms. We cultured bacteria from the shell surface of the endemic mussel Congeria jalzici, one of the three known cave mussels in the world that lives in the Dinaric karst underground. The bacterial isolates were obtained by swabbing the shell surface of mussels living in microhabitats with different amounts of water: 10 air-exposed mussels, 10 submerged mussels, and 10 mussels in the hygropetric zone. A collection of 87 pure culture isolates was obtained, mostly belonging to the phylum Bacillota (72%), followed by Pseudomonadota (16%), Actinomycetota (11%), and Bacteroidota (1%). We compared the results of MALDI-TOF MS identification (Bruker databases DB-5989 and version 11, v11) with the results of 16S rDNA-based phylogenetic analysis, a standard procedure for bacterial identification. Identification to the genus level based on 16S rDNA was possible for all isolates and clearly outperformed the results from MALDI-TOF MS, although the updated MALDI-TOF MS database v11 gave better results than the DB-5989 version (85% versus 62%). However, identification to the species-level by 16S rDNA sequencing was achieved for only 17% of isolates, compared with 14% and 40% for the MALDI-TOF MS databases DB-5989 and v11 database, respectively. In conclusion, our results suggest that continued enrichment of MALDI-TOF MS libraries will result with this method soon becoming a rapid, accurate, and efficient tool for assessing the diversity of culturable bacteria from different environmental niches.

Enhanced specificity of Bacillus metataxonomics using a tuf-targeted amplicon sequencing approach

Bacillus species are ubiquitous in nature and have tremendous application potential in agriculture, medicine, and industry. However, the individual species of this genus vary widely in both ecological niches and functional phenotypes, which, hence, requires accurate classification of these bacteria when selecting them for specific purposes. Although analysis of the 16S rRNA gene has been widely used to disseminate the taxonomy of most bacterial species, this gene fails proper classification of Bacillus species. To circumvent this restriction, we designed novel primers and optimized them to allow exact species resolution of Bacillus species in both synthetic and natural communities using high-throughput amplicon sequencing. The primers designed for the tuf gene were not only specific for the Bacillus genus but also sufficiently discriminated species both in silico and in vitro in a mixture of 11 distinct Bacillus species. Investigating the primers using a natural soil sample, 13 dominant species were detected including Bacillus badius, Bacillus velezensis, and Bacillus mycoides as primary members, neither of which could be distinguished with 16S rRNA sequencing. In conclusion, a set of high-throughput primers were developed which allows unprecedented species-level identification of Bacillus species and aids the description of the ecological distribution of Bacilli in various natural environment.

Antibiotic resistance genes fate during food waste management - Comparison between thermal treatment, hyperthermophilic composting, and anaerobic membrane bioreactor

The fate of eight different antibiotic resistance genes (ARGs) in food waste (sul1, sul2, tetO, tetW, ermF, ermB, ampC, oxa-1), intI1, and rpoB were monitored during thermal treatment (pyrolysis and incineration), hyperthermophilic composting, and anaerobic membrane bioreactor (AnMBR) treatment. ARGs in food waste ranged from 2.9 × 106 to 3.5 × 109 copies/kg with ampC being the least abundant and sul1 being the most abundant. Thermal treatment achieved removal below detection limits of all ARGs. Only two ARGs (sul1 and ampC) persisted in hyperthermophilic composting. While all genes except for ermB decreased in the AnMBR effluent relative to the food waste feed, sul1 remained at relatively high abundance. Biosolids on the contrary, accumulated tetO, ampC and sul2 in all tested operating conditions. Thermal treatment, despite limited resource recovery, provides the most effective mitigation of ARG risk in food waste.

Microbial Diversity Biased Estimation Caused by Intragenomic Heterogeneity and Interspecific Conservation of 16S rRNA Genes

The 16S rRNA gene has been extensively used as a molecular marker to explore evolutionary relationships and profile microbial composition throughout various environments. Despite its convenience and prevalence, limitations are inevitable. Variable copy numbers, intragenomic heterogeneity, and low taxonomic resolution have caused biases in estimating microbial diversity. Here, analysis of 24,248 complete prokaryotic genomes indicated that the 16S rRNA gene copy number ranged from 1 to 37 in bacteria and 1 to 5 in archaea, and intragenomic heterogeneity was observed in 60% of prokaryotic genomes, most of which were below 1%. The overestimation of microbial diversity caused by intragenomic variation and the underestimation introduced by interspecific conservation were calculated when using full-length or partial 16S rRNA genes. Results showed that, at the 100% threshold, microbial diversity could be overestimated by as much as 156.5% when using the full-length gene. The V4 to V5 region-based analyses introduced the lowest overestimation rate (4.4%) but exhibited slightly lower species resolution than other variable regions under the 97% threshold. For different variable regions, appropriate thresholds rather than the canonical value 97% were proposed for minimizing the risk of splitting a single genome into multiple clusters and lumping together different species into the same cluster. This study has not only updated the 16S rRNA gene copy number and intragenomic variation information for the currently available prokaryotic genomes, but also elucidated the biases in estimating prokaryotic diversity with quantitative data, providing references for choosing amplified regions and clustering thresholds in microbial community surveys. IMPORTANCE Microbial diversity is typically analyzed using marker gene-based methods, of which 16S rRNA gene sequencing is the most widely used approach. However, obtaining an accurate estimation of microbial diversity remains a challenge, due to the intragenomic variation and low taxonomic resolution of 16S rRNA genes. Comprehensive examination of the bias in estimating such prokaryotic diversity using 16S rRNA genes within ever-increasing prokaryotic genomes highlights the importance of the choice of sequencing regions and clustering thresholds based on the specific research objectives.

Modularity and diversity of target selectors in Tn7 transposons

To spread, transposons must integrate into target sites without disruption of essential genes while avoiding host defense systems. Tn7-like transposons employ multiple mechanisms for target-site selection, including protein-guided targeting and, in CRISPR-associated transposons (CASTs), RNA-guided targeting. Combining phylogenomic and structural analyses, we conducted a broad survey of target selectors, revealing diverse mechanisms used by Tn7 to recognize target sites, including previously uncharacterized target-selector proteins found in newly discovered transposable elements (TEs). We experimentally characterized a CAST I-D system and a Tn6022-like transposon that uses TnsF, which contains an inactivated tyrosine recombinase domain, to target the comM gene. Additionally, we identified a non-Tn7 transposon, Tsy, encoding a homolog of TnsF with an active tyrosine recombinase domain, which we show also inserts into comM. Our findings show that Tn7 transposons employ modular architecture and co-opt target selectors from various sources to optimize target selection and drive transposon spread.

Effect of vaginal flora on clinical outcome of frozen embryo transfer

Objective

Vaginal microbiota imbalance is a high risk factor for premature birth, and is closely related to female pelvic inflammation and sexually transmitted diseases. The effect of vaginal microbiota on the outcome of assisted reproductive technology is still unclear. In this study, the vaginal microbial composition and the pregnancy outcome of frozen embryo transfer (FET) was investigated.

Methods

From October 2020 to December 2021, 275 FET cycles were selected from the center of reproductive medicine in Affiliated Hospital of Nantong University. Vaginal secretions were collected on the day of endometrium transformation, and smears were Gram stained. According to the Nugent score they were divided into three groups, including normal group, mild dysbiosis group and sever dysbiosis group. The clinical outcomes of each group were compared.

Results

In 275 FET cycles, the embryo implantation rate, clinical pregnancy rate and ongoing pregnancy rate in the normal group (66.9%,84.3% and 83.1%) were significantly higher than those in the mild dysbiosis group (45.5%, 57.3% and 49.3%) and in sever dysbiosis group (29.6%, 34.2% and 27%). The difference was statistically significant (P<0.01). When compared the preclinical pregnancy loss rate and the miscarriage rate, the normal group (1.3% and 1.3%) was significantly lower than those in the mild dysbiosis group (20.4% and 14.0%) and the sever dysbiosis group (25.5% and 21.1%). The difference was statistically significant (P<0.01), but there was no significant difference between the mild dysbiosis group and sever dysbiosis group (P>0.05).

Conclusion

Nugent score is directly related to the clinical outcome of FET. The Lactobacilli-dominant vaginal flora was a favorable factor for the good clinical outcome of FET, while asymptomatic bacterial vaginosis had a negative correlation with the outcome of FET.

MALDI-TOF MS database expansion for identification of Bacillus and related genera isolated from a pharmaceutical facility

Bacillus and related genera are among the main bacterial groups isolated from pharmaceutical production areas. The identification of Bacillus species and related genera by classical methods is particularly difficult, due to similarities between closely related species. The Matrix-Assisted Laser Desorption Ionization-Time of Flight Mass Spectrometry (MALDI-TOF MS) is one of the most promising techniques for chemotaxonomic characterization of microorganisms, being an alternative to genotypic methods. This study aimed to identify Bacillus strains and related genera isolated from immunobiological production areas by phylogenetic analysis of housekeeping genes and expand the database associated with MALDI-TOF MS to improve their identification. In a previous study, 97 aerobic endospore-forming bacteria isolated from a pharmaceutical facility were analyzed by MALDI-TOF MS and 16S rRNA gene full-length sequencing. All strains were identified as Bacillus and related genera by the latest methodology. Among the 97 strains, 22 were unidentified and 2 strains were misidentified by MALDI-TOF MS. In the present study, these 24 strains were subjected to 16S rRNA gene phylogenetic analysis. Strains not identified at species level by this methodology were submitted to rpoB gene phylogenetic analysis. After identifying the strains, 19 of the 24 strains were incubated for 24, 48, and 72 h on Tryptic Soy Agar and Sheep Blood Agar and subjected to analysis by MALDI-TOF MS. A SuperSpectrum for each strain was created and entered into the equipment database. Finally, the 24 strains were again submitted to proteomic analysis by MALDI-TOF MS, and, at this time, all were correctly identified. The genotypic identification of in-house isolated strains and the introduction of these spectra in MALDI-TOF MS, in order to obtain a customized database, proved to be an extremely effective tool in the identification of Bacillus and related genera from pharmaceutical industry origin.

Housekeeping gene gyrA, a potential molecular marker for Bacillus ecology study

Bacillus is a genus of microorganisms (bacteria) and contains many important commercial species used in industry, agriculture and healthcare. Many different Bacilli are relatively well understood at the single-cell level; however, molecular tools that determine the diversity and ecology of Bacillus community are limited, which limits our understanding of how the Bacillus community works. In the present study, we investigated the potential of the housekeeping gene gyrA as a molecular marker for determining the diversity of Bacillus species. The amplification efficiency for Bacillus species diversity could be greatly improved by primer design. Therefore, we designed a novel primer pair gyrA3 that can detect at least 92 Bacillus species and related species. For B. amyloliquefaciens, B. pumilus, and B. megaterium, we observed that the high variability of the gyrA gene allows for more detailed clustering at the subspecies level that cannot be achieved by the 16S rRNA gene. Since gyrA provides better phylogenetic resolution than 16S rRNA and informs on the diversity of the Bacillus community, we propose that the gyrA gene may have broad application prospects in the study of Bacillus ecology.

Molecular identification of nontuberculous mycobacteria using the rpoB, argH and cya genes analysis

Nontuberculous mycobacterial (NTM) infections are growing worldwide especially in immunocompromised individuals. Since treatment of NTM infections is species-specific, the precise identification of NTM to species level is critical for an optimal treatment. This study was aimed to identify different NTM species by sequencing the rpoB gene and evaluating the effectiveness of argH and cya gene markers. In total 64 clinical isolates suspected to NTM were collected. The identification of the isolates was done by standard conventional methods and PCR-based rpoB gene and sequence analysis. PCR sequencing of argH and cya genes was performed to evaluate the efficacy of these genes in identifying and differentiating different species and subspecies of NTM. Among 64 isolates tested, 51 (79.68%) were detected by conventional tests as NTM. The results of rpoB sequence analysis revealed that the 56 clinical isolates were identified in 10 species of NTM and 8 remaining isolates which showed ambiguous results by rpoB sequencing, application of argH and cya sequencing could detect these isolates. Furthermore, by using cya gene sequencing, M. abscessus subspecies were properly differentiated. Although the rpoB sequencing as a standard method, is beneficial for detecting various species of NTM, however, based on our findings, argH and cya gene markers have a superb ability to discriminate closely related species. Further investigations are required to verify our outcomes.

Relative expression and validation of Aeromonas salmonicida subsp. salmonicida reference genes during ex vivo and in vivo fish infection

The genus Aeromonas is found worldwide in freshwater and marine environments and has been implicated in the etiology of human and animal diseases. In fish, among Aeromonas species, A. salmonicida causes massive mortality and great economic losses in marine and continental aquaculture species. Currently, several aspects of the clinical signs and pathogenesis of this Gram-negative bacterium have been described; however, determination of an appropriate reference gene is essential to normalize cellular mRNA data remain unknown. Here we evaluate the stability of seven candidate reference genes to be used for data normalization during ex vivo and in vivo experiments conducted in Atlantic cod, Atlantic salmon, and lumpfish. To assess this, raw Ct values obtained were evaluated by using geNorm, NormFinder, BestKeeper, Delta Ct comparison, and the comprehensive ranking, through the bioinformatic open-access portal RefFinder. We determined that fabD and era were most suitable reference genes in Atlantic cod primary macrophages, hfq and era in Atlantic salmon primary macrophages, rpoB and fabD in lumpfish head kidney samples, and hfq and era in lumpfish spleen. Our study demonstrates that use of multiple reference genes and its validation before measurements helps to minimize variability arising in qPCR studies that evaluate A. salmonicida gene expression in fish tissues. Overall, this study provided with an expanded list of reliable reference genes for A. salmonicida gene expression using qPCR during fish infection studies.

Analysis of Streptomyces Volatilomes Using Global Molecular Networking Reveals the Presence of Metabolites with Diverse Biological Activities

Streptomyces species produce a wide variety of specialized metabolites, some of which are used for communication or competition for resources in their natural environments. In addition, many natural products used in medicine and industry are derived from Streptomyces, and there has been interest in their capacity to produce volatile organic compounds (VOCs) for different industrial and agricultural applications. Recently, a machine-learning workflow called MSHub/GNPS was developed, which enables auto-deconvolution of gas chromatography-mass spectrometry (GC-MS) data, molecular networking, and library search capabilities, but it has not been applied to Streptomyces volatilomes. In this study, 131 Streptomyces isolates from the island of Newfoundland were phylogenetically typed, and 37 were selected based on their phylogeny and growth characteristics for VOC analysis using both a user-guided (conventional) and an MSHub/GNPS-based approach. More VOCs were annotated by MSHub/GNPS than by the conventional method. The number of unknown VOCs detected by the two methods was higher than those annotated, suggesting that many novel compounds remain to be identified. The molecular network generated by GNPS can be used to guide the annotation of such unknown VOCs in future studies. However, the number of overlapping VOCs annotated by the two methods is relatively small, suggesting that a combination of analysis methods might be required for robust volatilome analysis. More than half of the VOCs annotated with high confidence by the two approaches are plant-associated, many with reported bioactivities such as insect behavior modulation. Details regarding the properties and reported functions of such VOCs are described. IMPORTANCE This study represents the first detailed analysis of Streptomyces volatilomes using MSHub/GNPS, which in combination with a routinely used conventional method led to many annotations. More VOCs could be annotated using MSHub/GNPS as compared to the conventional method, many of which have known antimicrobial, anticancer, and insect behavior-modulating activities. The identification of numerous plant-associated VOCs by both approaches in the current study suggests that their production could be a more widespread phenomenon by members of the genus, highlighting opportunities for their large-scale production using Streptomyces. Plant-associated VOCs with antimicrobial activities, such as 1-octen-3-ol, octanol, and phenylethyl alcohol, have potential applications as fumigants. Furthermore, many of the annotated VOCs are reported to influence insect behavior, alluding to a possible explanation for their production based on the functions of other recently described Streptomyces VOCs in dispersal and nutrient acquisition.

Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies

BACKGROUND

The 16S rRNA gene is used extensively in bacterial phylogenetics, in species delineation, and now widely in microbiome studies. However, the gene suffers from intragenomic heterogeneity, and reports of recombination and an unreliable phylogenetic signal are accumulating. Here, we compare core gene phylogenies to phylogenies constructed using core gene concatenations to estimate the strength of signal for the 16S rRNA gene, its hypervariable regions, and all core genes at the intra- and inter-genus levels. Specifically, we perform four intra-genus analyses (Clostridium, n = 65; Legionella, n = 47; Staphylococcus, n = 36; and Campylobacter, n = 17) and one inter-genus analysis [41 core genera of the human gut microbiome (31 families, 17 orders, and 12 classes), n = 82].

RESULTS

At both taxonomic levels, the 16S rRNA gene was recombinant and subject to horizontal gene transfer. At the intra-genus level, the gene showed one of the lowest levels of concordance with the core genome phylogeny (50.7% average). Concordance for hypervariable regions was lower still, with entropy masking providing little to no benefit. A major factor influencing concordance was SNP count, which showed a positive logarithmic association. Using this relationship, we determined that 690 ± 110 SNPs were required for 80% concordance (average 16S rRNA gene SNP count was 254). We also found a wide range in 16S-23S-5S rRNA operon copy number among genomes (1-27). At the inter-genus level, concordance for the whole 16S rRNA gene was markedly higher (73.8% - 10th out of 49 loci); however, the most concordant hypervariable regions (V4, V3-V4, and V1-V2) ranked in the third quartile (62.5 to 60.0%).

CONCLUSIONS

Ramifications of a poor phylogenetic performance for the 16S rRNA gene are far reaching. For example, in addition to incorrect species/strain delineation and phylogenetic inference, it has the potential to confound community diversity metrics if phylogenetic information is incorporated - for example, with popular approaches such as Faith's phylogenetic diversity and UniFrac. Our results highlight the problematic nature of these approaches and their use (along with entropy masking) is discouraged. Lastly, the wide range in 16S rRNA gene copy number among genomes also has a strong potential to confound diversity metrics. Video Abstract.

Comparison of Actinobacteria communities from human‐impacted and pristine karst caves

Actinobacteria are important cave inhabitants, but knowledge of how anthropization and anthropization‐related visual marks affect this community on cave walls is lacking. We compared Actinobacteria communities among four French limestone caves (Mouflon, Reille, Rouffignac, and Lascaux) ranging from pristine to anthropized, and within Lascaux Cave between marked (wall visual marks) and unmarked areas in different rooms (Sas‐1, Passage, Apse, and Diaclase). In addition to the 16S rRNA gene marker, 441 bp fragments of the hsp65 gene were used and an hsp65‐related taxonomic database was constructed for the identification of Actinobacteria to the species level by Illumina‐MiSeq analysis. The hsp65 marker revealed higher resolution for species and higher richness (99% operational taxonomic units cutoff) versus the 16S rRNA gene; however, more taxa were identified at higher taxonomic ranks. Actinobacteria communities varied between Mouflon and Reille caves (both pristine), and Rouffignac and Lascaux (both anthropized). Rouffignac displayed high diversity of Nocardia, suggesting human inputs, and Lascaux exhibited high Mycobacterium relative abundance, whereas Gaiellales were typical in pristine caves and the Diaclase (least affected area of Lascaux Cave). Within Lascaux, Pseudonocardiaceae dominated on unmarked walls and Streptomycetaceae (especially Streptomyces mirabilis) on marked walls, indicating a possible role in mark formation. A new taxonomic database  was developed. Although not all Actinobacteria species were represented, the use of the hsp65 marker enabled species‐level variations of the Actinobacteria community to be documented based on the extent of anthropogenic pressure. This approach proved effective when comparing different limestone caves or specific conditions within one cave.

Fine-Scale Adaptations to Environmental Variation and Growth Strategies Drive Phyllosphere Methylobacterium Diversity

Methylobacterium is a prevalent bacterial genus of the phyllosphere. Despite its ubiquity, little is known about the extent to which its diversity reflects neutral processes like migration and drift, versus environmental filtering of life history strategies and adaptations. In two temperate forests, we investigated how phylogenetic diversity within Methylobacterium is structured by biogeography, seasonality, and growth strategies. Using deep, culture-independent barcoded marker gene sequencing coupled with culture-based approaches, we uncovered a considerable diversity of Methylobacterium in the phyllosphere. We cultured different subsets of Methylobacterium lineages depending upon the temperature of isolation and growth (20°C or 30°C), suggesting long-term adaptation to temperature. To a lesser extent than temperature adaptation, Methylobacterium diversity was also structured across large (>100 km; between forests) and small (<1.2 km; within forests) geographical scales, among host tree species, and was dynamic over seasons. By measuring the growth of 79 isolates during different temperature treatments, we observed contrasting growth performances, with strong lineage- and season-dependent variations in growth strategies. Finally, we documented a progressive replacement of lineages with a high-yield growth strategy typical of cooperative, structured communities in favor of those characterized by rapid growth, resulting in convergence and homogenization of community structure at the end of the growing season. Together, our results show how Methylobacterium is phylogenetically structured into lineages with distinct growth strategies, which helps explain their differential abundance across regions, host tree species, and time. This work paves the way for further investigation of adaptive strategies and traits within a ubiquitous phyllosphere genus. IMPORTANCE Methylobacterium is a bacterial group tied to plants. Despite the ubiquity of methylobacteria and the importance to their hosts, little is known about the processes driving Methylobacterium community dynamics. By combining traditional culture-dependent and -independent (metabarcoding) approaches, we monitored Methylobacterium diversity in two temperate forests over a growing season. On the surface of tree leaves, we discovered remarkably diverse and dynamic Methylobacterium communities over short temporal (from June to October) and spatial (within 1.2 km) scales. Because we cultured different subsets of Methylobacterium diversity depending on the temperature of incubation, we suspected that these dynamics partly reflected climatic adaptation. By culturing strains under laboratory conditions mimicking seasonal variations, we found that diversity and environmental variations were indeed good predictors of Methylobacterium growth performances. Our findings suggest that Methylobacterium community dynamics at the surface of tree leaves results from the succession of strains with contrasting growth strategies in response to environmental variations.

Applicability of rpoB Gene for PCR-RFLP based Discrimination of Bifidobacterial Species Isolated from Human and Animal Sources

Bifidobacteria are widely used as probiotics for their application in the development of functional food and prophylactic therapy. This has necessitated the development of a molecular approach for the genera to be widely identified up to species and subspecies level. In the current study, PCR-RFLP of the partial RNA polymerase β-subunit (rpoB) gene fragment was evaluated for differential identification of Bifidobacterium species. The rpoB gene partial sequences of 575 bp were amplified from 93 previously identified isolates collected from various sources of human and animal origin along with 12 standard reference strains. The PCR amplified products were digested with three restriction endonucleases HhaI, HinfI and BanI separately. Dendrograms constructed from the patterns of HhaI, were found to be more discriminatory and successfully differentiated all the twelve species and also at sub-species level in between B. longum subsp. longum and B. longum subsp. infantis. However, B. adolescentis and B. pseudocatenulatum group clusters were not separated and represented by one group. The groups were further discriminated by HinfI restriction digestion. A separate combination thereof may be used for inferring the classification of bifidobacterial species targeted on rpoB PCR-RFLP analysis. To our knowledge, this work is the first report based on use of rpoB PCR-RFLP for discrimination of the isolates of genus Bifidobacterium and also provides insights into specific advantages of this method over hsp60 PCR-RFLP in differentiating B. longum subsp. longum and B. longum subsp. infantis.

ddPCR allows 16S rRNA gene amplicon sequencing of very small DNA amounts from low-biomass samples

BACKGROUND

One limiting factor of short amplicon 16S rRNA gene sequencing approaches is the use of low DNA amounts in the amplicon generation step. Especially for low-biomass samples, insufficient or even commonly undetectable DNA amounts can limit or prohibit further analysis in standard protocols.

RESULTS

Using a newly established protocol, very low DNA input amounts were found sufficient for reliable detection of bacteria using 16S rRNA gene sequencing compared to standard protocols. The improved protocol includes an optimized amplification strategy by using a digital droplet PCR. We demonstrate how PCR products are generated even when using very low concentrated DNA, unable to be detected by using a Qubit. Importantly, the use of different 16S rRNA gene primers had a greater effect on the resulting taxonomical profiles compared to using high or very low initial DNA amounts.

CONCLUSION

Our improved protocol takes advantage of ddPCR and allows faithful amplification of very low amounts of template. With this, samples of low bacterial biomass become comparable to those with high amounts of bacteria, since the first and most biasing steps are the same. Besides, it is imperative to state DNA concentrations and volumes used and to include negative controls indicating possible shifts in taxonomical profiles. Despite this, results produced by using different primer pairs cannot be easily compared.

Using amplicon sequencing of rpoB for identification of inoculant rhizobia from peanut nodules

To improve the nitrogen fixation, legume crops are often inoculated with selected effective rhizobia. However, there is large variation in how well the inoculant strains compete with the indigenous microflora in soil. To assess the success of the inoculant, it is necessary to distinguish it from other, closely related strains. Methods used until now have generally been based either on fingerprinting methods or on the use of reporter genes. Nevertheless, these methods have their shortcomings, either because they do not provide sufficiently specific information on the identity of the inoculant strain, or because they use genetically modified organisms that need prior authorization to be applied in the field or other uncontained environments. Another possibility is to target a gene that is naturally present in the bacterial genomes. Here we have developed a method that is based on amplicon sequencing of the bacterial housekeeping gene rpoB, encoding the beta-subunit of the RNA polymerase, which has been proposed as an alternative to the 16S rRNA gene to study the diversity of rhizobial populations in soils. We evaluated the method under laboratory and field conditions. Peanut seeds were inoculated with various Bradyrhizobium strains. After nodule development, DNA was extracted from selected nodules and the nodulating rhizobia were analysed by amplicon sequencing of the rpoB gene. The analyses of the sequence data showed that the method reliably identified bradyrhizobial strains in nodules, at least at the species level, and could be used to assess the competitiveness of the inoculant compared to other bradyrhizobia.

Phylogenomic Analysis Substantiates the gyrB Gene as a Powerful Molecular Marker to Efficiently Differentiate the Most Closely Related Genera Myxococcus, Corallococcus, and Pyxidicoccus

Rapid and accurate strain identification of the most closely related genera Myxococcus, Corallococcus, and Pyxidicoccus can enhance the efficiency of the mining of novel secondary metabolites through dereplication. However, the commonly used 16S rRNA gene sequencing cannot accurately differentiate members of the three genera above, and the whole-genome sequencing is unable to rapidly and inexpensively provide species assignation toward a large number of isolates. To overcome the limitations, the gyrB gene was investigated as a candidate genetic marker for exploring the phylogenetic relationships of bacteria within the three genera and for developing the gyrB-based typing method. Here, the bacterial phylogeny and species affiliations of the three genera were determined based on the phylogenomic reconstruction and the analysis of digital DNA–DNA hybridization values among 90 genomes, further confirming nine novel taxa and assigning over one-third of genomes to defined species. The phylogenetic relationships of these strains based on the gyrB gene sequences were congruent with those based on their genome sequences, allowing the use of the gyrB gene as a molecular marker. The gyrB gene-specific primers for the PCR-amplification and sequencing of bacteria within the three genera were designed and validated for 31 isolates from our group collection. The gyrB-based taxonomic tool proved to be able to differentiate closely related isolates at the species level. Based on the newly proposed 98.6% identity threshold for the 966-bp gyrB gene and the phylogenetic inference, these isolates were assigned into two known species and eight additional putative new species. In summary, this report demonstrated that the gyrB gene is a powerful phylogenetic marker for taxonomy and phylogeny of bacteria within the closely related genera Myxococcus, Corallococcus, and Pyxidicoccus, particularly in the case of hundreds or thousands of isolates in environmental studies.

Nutrient supplementation experiments with saltern microbial communities implicate utilization of DNA as a source of phosphorus

All environments including hypersaline ones harbor measurable concentrations of dissolved extracellular DNA (eDNA) that can be utilized by microbes as a nutrient. However, it remains poorly understood which eDNA components are used, and who in a community utilizes it. For this study, we incubated a saltern microbial community with combinations of carbon, nitrogen, phosphorus, and DNA, and tracked the community response in each microcosm treatment via 16S rRNA and rpoB gene sequencing. We show that microbial communities used DNA only as a phosphorus source, and provision of other sources of carbon and nitrogen was needed to exhibit a substantial growth. The taxonomic composition of eDNA in the water column changed with the availability of inorganic phosphorus or supplied DNA, hinting at preferential uptake of eDNA from specific organismal sources. Especially favored for growth was eDNA from the most abundant taxa, suggesting some haloarchaea prefer eDNA from closely related taxa. The preferential eDNA consumption and differential growth under various nutrient availability regimes were associated with substantial shifts in the taxonomic composition and diversity of microcosm communities. Therefore, we conjecture that in salterns the microbial community assembly is driven by the available resources, including eDNA.

Historical, current, and emerging tools for identification and serotyping of Shigella

The Shigella genus includes serious foodborne disease etiologic agents, with 4 species and 54 serotypes. Identification at species and serotype levels is a crucial task in microbiological laboratories. Nevertheless, the genetic similarity between Shigella spp. and Escherichia coli challenges the correct identification and serotyping of Shigella spp., with subsequent negative repercussions on surveillance, epidemiological investigations, and selection of appropriate treatments. For this purpose, multiple techniques have been developed historically ranging from phenotype-based methods and single or multilocus molecular techniques to whole-genome sequencing (WGS). To facilitate the selection of the most relevant method, we herein provide a global overview of historical and emerging identification and serotyping techniques with a particular focus on the WGS-based approaches. This review highlights the excellent discriminatory power of WGS to more accurately elucidate the epidemiology of Shigella spp., disclose novel promising genomic targets for surveillance methods, and validate previous well-established methods.

Mycorrhizal symbiosis modulates the rhizosphere microbiota to promote rhizobia-legume symbiosis

Plants establish symbioses with mutualistic fungi, such as arbuscular mycorrhizal (AM) fungi, and bacteria, such as rhizobia, to exchange key nutrients and thrive. Plants and symbionts have coevolved and represent vital components of terrestrial ecosystems. Plants employ an ancestral AM signaling pathway to establish intracellular symbioses, including the legume-rhizobia symbiosis, in their roots. Nevertheless, the relationship between the AM and rhizobial symbioses in native soil is poorly understood. Here, we examined how these distinct symbioses affect root-associated bacterial communities in Medicago truncatula by performing quantitative microbiota profiling (QMP) of 16S rRNA genes. We found that M. truncatula mutants that cannot establish AM or rhizobia symbiosis have an altered microbial load (quantitative abundance) in the rhizosphere and roots, and in particular that AM symbiosis is required to assemble a normal quantitative root-associated microbiota in native soil. Moreover, quantitative microbial co-abundance network analyses revealed that AM symbiosis affects Rhizobiales hubs among plant microbiota and benefits the plant holobiont. Through QMP of rhizobial rpoB and AM fungal SSU rRNA genes, we revealed a new layer of interaction whereby AM symbiosis promotes rhizobia accumulation in the rhizosphere of M. truncatula. We further showed that AM symbiosis-conditioned microbial communities within the M. truncatula rhizosphere could promote nodulation in different legume plants in native soil. Given that the AM and rhizobial symbioses are critical for crop growth, our findings might inform strategies to improve agricultural management. Moreover, our work sheds light on the co-evolution of these intracellular symbioses during plant adaptation to native soil conditions.

Differential prevalence and host-association of antimicrobial resistance traits in disinfected and non-disinfected drinking water systems

Antimicrobial resistance (AMR) in drinking water has received less attention than its counterparts in the urban water cycle. While culture-based techniques or gene-centric PCR have been used to probe the impact of treatment approaches (e.g., disinfection) on AMR in drinking water, to our knowledge there is no systematic comparison of AMR trait distribution and prevalence between disinfected and disinfectant residual-free drinking water systems. We used metagenomics to assess the associations between disinfectant residuals and AMR prevalence and its host association in full-scale drinking water distribution systems (DWDSs) with and without disinfectant residuals. While the differences in AMR profiles between DWDSs were associated with the presence or absence of disinfectant, they were also associated with overall water chemistry and more importantly with microbial community structure. AMR genes and mechanisms differentially abundant in disinfected systems were primarily associated with nontuberculous mycobacteria (NTM). Finally, evaluation of metagenome assembled genomes (MAGs) also suggests that NTM possessing AMR genes conferring intrinsic resistance to key antibiotics were prevalent in disinfected systems, whereas such NTM genomes were not detected in disinfectant residual free DWDSs. Altogether, our findings provide insights into the drinking water resistome and its association with potential opportunistic pathogens, particularly in systems with disinfectant residual.

Membrane Fouling Inversely Impacts Intracellular and Extracellular Antibiotic Resistance Gene Abundances in the Effluent of an Anaerobic Membrane Bioreactor

Anaerobic membrane bioreactors (AnMBRs) can significantly reduce the release of antibiotic resistance elements to the environment. The purpose of this study was to elucidate the role of membrane fouling layers (biofilms) in mitigating the release of intracellular and extracellular antibiotic resistance genes (iARGs and eARGs) from an AnMBR. The AnMBR was equipped with three membrane modules, each exhibiting a different level of fouling. Results showed that the absolute abundance of ARGs decreased gradually in the suspended biomass during operation of the AnMBR. Normalized abundances of targeted ARGs and intI1 were found to be significantly higher in the fouling layers compared to the suspended biomass, implying adsorption or an increased potential for horizontal gene transfer of ARGs in the biofilm. Effluent ARG data revealed that the highly fouled (HF) membrane significantly reduced the absolute abundance of eARGs. However, the HF membrane effluent concomitantly had the highest absolute abundance of iARGs. Nevertheless, total ARG abundance (sum of iARG and eARG) in the effluent of the AnMBR was not impacted by the extent of fouling. These results suggest a need for a combination of different treatment technologies to effectively prevent antibiotic resistance proliferation associated with these two ARG fractions.

A Novel Primer Mixture for GH48 Genes: Quantification and Identification of Truly Cellulolytic Bacteria in Biogas Fermenters

Genomic studies revealed the glycoside hydrolases of family 48 (GH48) as a powerful marker for the identification of truly cellulolytic bacteria. Here we report an improved method for detecting cellulolytic bacteria in lab-scale biogas fermenters by using GH48 genes as a molecular marker in DNA and RNA samples. We developed a mixture of primers for the specific amplification of a GH48 gene region in a broad range of bacteria. Additionally, we built a manually curated reference database containing GH48 gene sequences directly linked to the corresponding taxonomic information. Phylogenetic correlation analysis of GH48 to 16S rRNA gene sequences revealed that GH48 gene sequences with 94% identity belong with high confidence to the same genus. Applying this analysis, GH48 amplicon reads revealed that at mesophilic fermenter conditions, 50-99% of the OTUs appear to belong to novel taxa. In contrast, at thermophilic conditions, GH48 gene sequences from the genus Hungateiclostridium dominated with 60-91% relative abundance. The novel primer combinations enabled detection and relative quantification of a wide spectrum of GH48 genes in cellulolytic microbial communities. Deep phylogenetic correlation analysis and a simplified taxonomic identification with the novel database facilitate identification of cellulolytic organisms, including the detection of novel taxa in biogas fermenters.

Challenges in identification of enteroinvasive Escherichia coli and Shigella spp. in Lebanon

This study aimed to evaluate the routine identification tools available in Lebanon for differentiation of Escherichia coli and Shigella spp. The identification of 43 isolates defined as Shigella spp. by Api 20E was accessed using MALDI-TOF, serological testing, duplex PCR targeting ipaH (present in Shigella spp. and enteroinvasive E. coli "EIEC") and lacY (found in E. coli including EIEC but not Shigella spp.) as well as gyrB gene sequencing. Antibiotic susceptibility was investigated as well as Shiga-toxin production. All isolates were identified as E. coli by MALDI-TOF while the PCR showed a disparate group of 26 EIEC, 11 Shigella spp., 5 E. coli and 1 inactive E. coli. However, the sequencing of gyrB gene, which was recently described as a suitable marker for distinguishing E. coli and Shigella spp., identified all isolates as E. coli. Antibiotic resistance was noticeable against ß-lactams, rifampicin, trimethoprim-sulfamethoxazole, gentamicin, and ciprofloxacin. The most common variants of beta-lactamase genes were blaTEM-1, blaCTX-M-15, and blaCTX-M-3. A great discordance between the used methods in identification was revealed herein. An accurate identification technique able to distinguish E. coli from Shigella spp. in routine laboratories is a pressing need in order to select the appropriate treatment and assess the epidemiology of these bacteria.

The Human Oral Microbiome in Health and Disease: From Sequences to Ecosystems

Abstract: The human oral cavity is home to an abundant and diverse microbial community (i.e., the oral microbiome), whose composition and roles in health and disease have been the focus of intense research in recent years. Thanks to developments in sequencing-based approaches, such as 16S ribosomal RNA metabarcoding, whole metagenome shotgun sequencing, or meta-transcriptomics, we now can efficiently explore the diversity and roles of oral microbes, even if unculturable. Recent sequencing-based studies have charted oral ecosystems and how they change due to lifestyle or disease conditions. As studies progress, there is increasing evidence of an important role of the oral microbiome in diverse health conditions, which are not limited to diseases of the oral cavity. This, in turn, opens new avenues for microbiome-based diagnostics and therapeutics that benefit from the easy accessibility of the oral cavity for microbiome monitoring and manipulation. Yet, many challenges remain ahead. In this review, we survey the main sequencing-based methodologies that are currently used to explore the oral microbiome and highlight major findings enabled by these approaches. Finally, we discuss future prospects in the field.

Genome sequencing as a new window into the microbial community of membrane bioreactors - A critical review

Recent developed sequencing techniques have resulted in a new and unprecedented way to study biological wastewater treatment, in which most organisms are uncultivable. This review provides (i) an insight on state-of-the-art sequencing techniques and their limitations; (ii) a critical assessment of the microbial community in biological reactor and biofouling layer in a membrane bioreactor (MBR). The data from high-throughput sequencing has been used to infer microbial growth conditions and metabolisms of microorganisms present in MBRs at the time of sampling. These data shed new insight to two fundamental questions about a microbial community in the MBR process namely the microbial composition (who are they?) and the functions of each specific microbial assemblage (what are their function?). The results to date also highlight the complexity of the microbial community growing on MBRs. Environmental conditions are dynamic and diverse, and can influence the diversity and structural dynamics of any given microbial community for wastewater treatment. The benefits of understanding the structure of microbial communities on three major aspects of the MBR process (i.e. nutrient removal, biofouling control, and micropollutant removal) were symmetrically delineated. This review also indicates that the deployment of microbial community analysis for a practical engineering context, in terms of process design and system optimization, can be further realized.

Comparative Genomic Analysis of Soil Dwelling Bacteria Utilizing a Combinational Codon Usage and Molecular Phylogenetic Approach Accentuating on Key Housekeeping Genes

Soil is a diversified and complex ecological niche, home to a myriad of microorganisms particularly bacteria. The physico-chemical complexities of soil results in a plethora of physiological variations to exist within the different types of soil dwelling bacteria, giving rise to a wide variation in genome structure and complexity. This serves as an attractive proposition to analyze and compare the genome of a large number soil bacteria to comprehend their genome complexity and evolution. In this study a combination of codon usage and molecular phylogenetics of the whole genome and key housekeeping genes like infB (translation initiation factor 2), trpB (tryptophan synthase, beta subunit), atpD (ATP synthase, beta subunit), and rpoB (RNA polymerase, beta subunit) of 92 soil bacterial species spread across the entire eubacterial domain and residing in different soil types was performed. The results indicated the direct relationship of genome size with codon bias and coding frequency in the studied bacteria. The codon usage profile demonstrated by the gene trpB was found to be relatively different from the rest of the housekeeping genes with a large number of bacteria having a greater percentage of genes with Nc values less than the Nc of trpB. The results from the overall codon usage bias profile also depicted that the codon usage bias in the key housekeeping genes of soil bacteria was majorly due to selectional pressure and not mutation. The analysis of hydrophobicity of the gene product encoded by the rpoB coding sequences demonstrated tight clustering across all the soil bacteria suggesting conservation of protein structure for maintenance of form and function. The phylogenetic affinities inferred using 16S rRNA gene and the housekeeping genes demonstrated conflicting signals with trpB gene being the noisiest one. The housekeeping gene atpD was found to depict the least amount of evolutionary change in the soil bacteria considered in this study except in two Clostridium species. The phylogenetic and codon usage analysis of the soil bacteria consistently demonstrated the relatedness of Azotobacter chroococcum with different species of the genus Pseudomonas.

rpoB, a promising marker for analyzing the diversity of bacterial communities by amplicon sequencing

Background

Microbiome composition is frequently studied by the amplification and high-throughput sequencing of specific molecular markers (metabarcoding). Various hypervariable regions of the 16S rRNA gene are classically used to estimate bacterial diversity, but other universal bacterial markers with a finer taxonomic resolution could be employed. We compared specificity and sensitivity between a portion of the rpoB gene and the V3 V4 hypervariable region of the 16S rRNA gene.

Results

We first designed universal primers for rpoB suitable for use with Illumina sequencing-based technology and constructed a reference rpoB database of 45,000 sequences. The rpoB and V3 V4 markers were amplified and sequenced from (i) a mock community of 19 bacterial strains from both Gram-negative and Gram-positive lineages; (ii) bacterial assemblages associated with entomopathogenic nematodes. In metabarcoding analyses of mock communities with two analytical pipelines (FROGS and DADA2), the estimated diversity captured with the rpoB marker resembled the expected composition of these mock communities more closely than that captured with V3 V4. The rpoB marker had a higher level of taxonomic affiliation, a higher sensitivity (detection of all the species present in the mock communities), and a higher specificity (low rates of spurious OTU detection) than V3 V4. We compared the performance of the rpoB and V3 V4 markers in an animal ecosystem model, the infective juveniles of the entomopathogenic nematode Steinernema glaseri carrying the symbiotic bacteria Xenorhabdus poinarii. Both markers showed the bacterial community associated with this nematode to be of low diversity (< 50 OTUs), but only rpoB reliably detected the symbiotic bacterium X. poinarii.

Conclusions

Our results confirm that different microbiota composition data may be obtained with different markers. We found that rpoB was a highly appropriate marker for assessing the taxonomic structure of mock communities and the nematode microbiota. Further studies on other ecosystems should be considered to evaluate the universal usefulness of the rpoB marker. Our data highlight two crucial elements that should be taken into account to ensure more reliable and accurate descriptions of microbial diversity in high-throughput amplicon sequencing analyses: i) the need to include mock communities as controls; ii) the advantages of using a multigenic approach including at least one housekeeping gene (rpoB is a good candidate) and one variable region of the 16S rRNA gene. This study will be useful to the growing scientific community describing bacterial communities by metabarcoding in diverse ecosystems.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1546-z) contains supplementary material, which is available to authorized users.

Validation and Retrospective Clinical Evaluation of a Quantitative 16S rRNA Gene Metagenomic Sequencing Assay for Bacterial Pathogen Detection in Body Fluids

Next-generation sequencing-based 16S rRNA gene metagenomic sequencing (16S MG) technology has tremendous potential for improving diagnosis of bacterial infections given its quantitative capability and culture-independent approach. We validated and used a quantitative 16S MG assay to identify and quantify bacterial species in clinical samples from a wide spectrum of infections, including meningitis, septic arthritis, brain abscess, intra-abdominal abscess, soft tissue abscess, and pneumonia. Twenty clinical samples were tested, and 16S MG identified a total of 34 species, compared with 22 species and three descriptive findings identified by culture. 16S MG results matched culture results in 75% (15/20) of the samples but detected at least one more species in five samples, including one culture-negative cerebrospinal fluid sample that was found to contain Streptococcus intermedius. Shotgun metagenomic sequencing verified the presence of all additional species. The 16S MG assay is highly sensitive, with a limit of detection of 10 to 100 colony-forming units/mL. Other performance characteristics, including linearity, precision, and specificity, all met the requirements for a clinical test. This assay showed the advantages of accurate identification and quantification of bacteria in culture-negative and polymicrobial infections for which conventional microbiology methods are limited. It also showed promises to serve unmet clinical needs for solving difficult infectious diseases cases.

Metagenomic insights into diazotrophic communities across Arctic glacier forefields

Microbial nitrogen fixation is crucial for building labile nitrogen stocks and facilitating higher plant colonisation in oligotrophic glacier forefield soils. Here, the diazotrophic bacterial community structure across four Arctic glacier forefields was investigated using metagenomic analysis. In total, 70 soil metagenomes were used for taxonomic interpretation based on 185 nitrogenase (nif) sequences, extracted from assembled contigs. The low number of recovered genes highlights the need for deeper sequencing in some diverse samples, to uncover the complete microbial populations. A key group of forefield diazotrophs, found throughout the forefields, was identified using a nifH phylogeny, associated with nifH Cluster I and III. Sequences related most closely to groups including Alphaproteobacteria, Betaproteobacteria, Cyanobacteria and Firmicutes. Using multiple nif genes in a Last Common Ancestor analysis revealed a diverse range of diazotrophs across the forefields. Key organisms identified across the forefields included Nostoc, Geobacter, Polaromonas and Frankia. Nitrogen fixers that are symbiotic with plants were also identified, through the presence of root associated diazotrophs, which fix nitrogen in return for reduced carbon. Additional nitrogen fixers identified in forefield soils were metabolically diverse, including fermentative and sulphur cycling bacteria, halophiles and anaerobes.

Diversion and phylogenetic relatedness of filterable bacteria from Norwegian tap and bottled waters

Numerous articles have documented the existence of filterable bacteria. Where filtration is the chosen method of sterilization for medicinal or media components, these bacteria will by definition render products non-sterile. They may further represent a health hazard to the end user. A wide-range of bacterial genera were found in bottled and tap water filtrates from 0.2 μm filters, including genera housing opportunistic pathogens (e.g. Methylobacterium) and endospore formers (Paenibacillus). Two municipal tap water isolates were only distantly related to named species. One of these grew on agar, and could potentially provide hitherto unharvested useful biological products. The other grew only in water, and failed to produce colonies on media targeting either heterotrophs or autotrophs. The present study is one of very few looking at filterable bacteria in bottled waters intended for human consumption and the first identifying the filterable portion. It extends the range of known habitats of filterable bacteria and provides data on two new or novel species.

Deciphering intra-species bacterial diversity of meat and seafood spoilage microbiota using gyrB amplicon sequencing: A comparative analysis with 16S rDNA V3-V4 amplicon sequencing

Meat and seafood spoilage ecosystems harbor extensive bacterial genomic diversity that is mainly found within a small number of species but within a large number of strains with different spoilage metabolic potential. To decipher the intraspecies diversity of such microbiota, traditional metagenetic analysis using the 16S rRNA gene is inadequate. We therefore assessed the potential benefit of an alternative genetic marker, gyrB, which encodes the subunit B of DNA gyrase, a type II DNA topoisomerase. A comparison between 16S rDNA-based (V3-V4) amplicon sequencing and gyrB-based amplicon sequencing was carried out in five types of meat and seafood products, with five mock communities serving as quality controls. Our results revealed that bacterial richness in these mock communities and food samples was estimated with higher accuracy using gyrB than using16S rDNA. However, for Firmicutes species, 35% of putative gyrB reads were actually identified as sequences of a gyrB paralog, parE, which encodes subunit B of topoisomerase IV; we therefore constructed a reference database of published sequences of both gyrB and pare for use in all subsequent analyses. Despite this co-amplification, the deviation between relative sequencing quantification and absolute qPCR quantification was comparable to that observed for 16S rDNA for all the tested species. This confirms that gyrB can be used successfully alongside 16S rDNA to determine the species composition (richness and evenness) of food microbiota. The major benefit of gyrB sequencing is its potential for improving taxonomic assignment and for further investigating OTU richness at the subspecies level, thus allowing more accurate discrimination of samples. Indeed, 80% of the reads of the 16S rDNA dataset were represented by thirteen 16S rDNA-based OTUs that could not be assigned at the species-level. Instead, these same clades corresponded to 44 gyrB-based OTUs, which differentiated various lineages down to the subspecies level. The increased ability of gyrB-based analyses to track and trace phylogenetically different groups of strains will generate improved resolution and more reliable results for studies of the strains implicated in food processes.

Testing culture purity in prokaryotes: criteria and challenges

Reliance on pure cultures was introduced at the beginning of microbiology as a discipline and has remained significant although their adaptive properties are essentially dissimilar from those of mixed cultures and environmental populations. They are needed for (i) taxonomic identification; (ii) diagnostics of pathogens; (iii) virulence and pathogenicity studies; (iv) elucidation of metabolic properties; (v) testing sensitivity to antibiotics; (vi) full-length genome assembly; (vii) strain deposition in microbial collections; and (viii) description of new species with name validation. Depending on the specific task there are alternative claims for culture purity, i.e., when conventional criteria are satisfied or when looking deeper is necessary. Conventional proof (microscopic and plating controls) has a low resolution and depends on the observer's personal judgement. Phenotypic criteria alone cannot prove culture purity and should be complemented with genomic criteria. We consider the possible use of DNA high-throughput culture sequencing data to define criteria for only one genospecies, axenic state detection panel and only one genome. The second and third of these are preferable, although their resolving capacity (depth) is limited. Because minor contaminants may go undetected, even with deep sequencing, the reliably pure culture would be a clonal culture launched from a single cell or trichome (multicellular bacterium). Although this type of culture is associated with technical difficulties and cannot be employed on a large scale (the corresponding inoculums may have low chances of growth when transferred to solid media), it is hoped that the high-throughput culturing methods introduced by 'culturomics' will overcome this obstacle.

Correcting for 16S rRNA gene copy numbers in microbiome surveys remains an unsolved problem

The 16S ribosomal RNA gene is the most widely used marker gene in microbial ecology. Counts of 16S sequence variants, often in PCR amplicons, are used to estimate proportions of bacterial and archaeal taxa in microbial communities. Because different organisms contain different 16S gene copy numbers (GCNs), sequence variant counts are biased towards clades with greater GCNs. Several tools have recently been developed for predicting GCNs using phylogenetic methods and based on sequenced genomes, in order to correct for these biases. However, the accuracy of those predictions has not been independently assessed. Here, we systematically evaluate the predictability of 16S GCNs across bacterial and archaeal clades, based on ∼ 6,800 public sequenced genomes and using several phylogenetic methods. Further, we assess the accuracy of GCNs predicted by three recently published tools (PICRUSt, CopyRighter, and PAPRICA) over a wide range of taxa and for 635 microbial communities from varied environments. We find that regardless of the phylogenetic method tested, 16S GCNs could only be accurately predicted for a limited fraction of taxa, namely taxa with closely to moderately related representatives (≲15% divergence in the 16S rRNA gene). Consistent with this observation, we find that all considered tools exhibit low predictive accuracy when evaluated against completely sequenced genomes, in some cases explaining less than 10% of the variance. Substantial disagreement was also observed between tools (R2<0.5) for the majority of tested microbial communities. The nearest sequenced taxon index (NSTI) of microbial communities, i.e., the average distance to a sequenced genome, was a strong predictor for the agreement between GCN prediction tools on non-animal-associated samples, but only a moderate predictor for animal-associated samples. We recommend against correcting for 16S GCNs in microbiome surveys by default, unless OTUs are sufficiently closely related to sequenced genomes or unless a need for true OTU proportions warrants the additional noise introduced, so that community profiles remain interpretable and comparable between studies.

Multiplex PCR identification and culture-independent quantification of Bacillus licheniformis by qPCR using specific DNA markers

Probiotics benefits in fish farming have been usually inferred appraising the effects observed on the host and not through the direct assessment of probiotic dynamics in the host gut microbiota. To overcome this gap, quantitative PCR (qPCR) can be a powerful approach to study the bacterial dynamics in fish gut microbiota. The presented work proposes four B. licheniformis-specific DNA markers and details a qPCR method to track putative probiotics B. licheniformis on fish gut. The four B. licheniformis-specific DNA markers - BL5B (hypothetical protein BL00303), BL8A (serA2), BL13C (rfaB) and BL18A (ligD) - were selected and validated by PCR and multiplex-PCR with 20 B. licheniformis isolates and a broad range of non-target bacteria. To assess the dynamics of B. licheniformis in the digesta of farmed fish, a qPCR was validated using markers BL8A and BL18A and calibration curves obtained for both markers with digesta samples spiked with B. licheniformis cells showed a high correlation (R² > 0.99) over 6 log units (CFU/reaction), and a limit of detection (LOD) as low as 247 CFUs/reaction. Furthermore, the consistent qPCR repeatability and reproducibility underline the specificity and reliability of the qPCR proposed. Ultimately, the possibility to monitor the dynamics of B. licheniformis probiotics in the gut microbiota of farmed fish might be instrumental to optimize best practices in aquaculture.

Accurate differentiation of Escherichia coli and Shigella serogroups: challenges and strategies

Shigella spp. and Escherichia coli are closely related; both belong to the family Enterobacteriaceae. Phenotypically, Shigella spp. and E. coli share many common characteristics, yet they have separate entities in epidemiology and clinical disease, which poses a diagnostic challenge. We collated information for the best possible approach to differentiate clinically relevant E. coli from Shigella spp. We found that a molecular approach is required for confirmation. High discriminatory potential is seen with whole genome sequencing analysed for k-mers and single nucleotide polymorphism. Among these, identification using single nucleotide polymorphism is easy to perform and analyse, and it thus appears more promising. Among the nonmolecular methods, matrix-assisted desorption ionization–time of flight mass spectrometry may be applicable when data analysis is assisted with advanced analytic tools.

Strategies and Methodologies for Developing Microbial Detoxification Systems to Mitigate Mycotoxins

Mycotoxins, the secondary metabolites of mycotoxigenic fungi, have been found in almost all agricultural commodities worldwide, causing enormous economic losses in livestock production and severe human health problems. Compared to traditional physical adsorption and chemical reactions, interest in biological detoxification methods that are environmentally sound, safe and highly efficient has seen a significant increase in recent years. However, researchers in this field have been facing tremendous unexpected challenges and are eager to find solutions. This review summarizes and assesses the research strategies and methodologies in each phase of the development of microbiological solutions for mycotoxin mitigation. These include screening of functional microbial consortia from natural samples, isolation and identification of single colonies with biotransformation activity, investigation of the physiological characteristics of isolated strains, identification and assessment of the toxicities of biotransformation products, purification of functional enzymes and the application of mycotoxin decontamination to feed/food production. A full understanding and appropriate application of this tool box should be helpful towards the development of novel microbiological solutions on mycotoxin detoxification.

It Is Imperative to Establish a Pellucid Definition of Chimeric RNA and to Clear Up a Lot of Confusion in the Relevant Research

There have been tens of thousands of RNAs deposited in different databases that contain sequences of two genes and are coined chimeric RNAs, or chimeras. However, "chimeric RNA" has never been lucidly defined, partly because "gene" itself is still ill-defined and because the means of production for many RNAs is unclear. Since the number of putative chimeras is soaring, it is imperative to establish a pellucid definition for it, in order to differentiate chimeras from regular RNAs. Otherwise, not only will chimeric RNA studies be misled but also characterization of fusion genes and unannotated genes will be hindered. We propose that only those RNAs that are formed by joining two RNA transcripts together without a fusion gene as a genomic basis should be regarded as authentic chimeras, whereas those RNAs transcribed as, and cis-spliced from, single transcripts should not be deemed as chimeras. Many RNAs containing sequences of two neighboring genes may be transcribed via a readthrough mechanism, and thus are actually RNAs of unannotated genes or RNA variants of known genes, but not chimeras. In today's chimeric RNA research, there are still several key flaws, technical constraints and understudied tasks, which are also described in this perspective essay.

Short-Term Exposure of Paddy Soil Microbial Communities to Salt Stress Triggers Different Transcriptional Responses of Key Taxonomic Groups

Soil salinization due to seawater intrusion along coastal areas is an increasing threat to rice cultivation worldwide. While the detrimental impact on rice growth and yield has been thoroughly studied, little is known about how severe salinity affects structure and function of paddy soil microbial communities. Here, we examined their short-term responses to half- and full-strength seawater salinity in controlled laboratory experiments. Slurry microcosms were incubated under anoxic conditions, with rice straw added as carbon source. Stress exposure time was for 2 days after a pre-incubation period of 7 days. Relative to the control, moderate (300 mM NaCl) and high (600 mM NaCl) salt stress suppressed both net consumption of acetate and methane production by 50% and 70%, respectively. Correspondingly, community-wide mRNA expression decreased by 50–65%, with significant changes in relative transcript abundance of family-level groups. mRNA turnover was clearly more responsive to salt stress than rRNA dynamics. Among bacteria, Clostridiaceae were most abundant and the only group whose transcriptional activity was strongly stimulated at 600 mM NaCl. In particular, clostridial mRNA involved in transcription/translation, fermentation, uptake and biosynthesis of compatible solutes, and flagellar motility was significantly enriched in response salt stress. None of the other bacterial groups were able to compete at 600 mM NaCl. Their responses to 300 mM NaCl were more diverse. Lachnospiraceae increased, Ruminococcaceae maintained, and Peptococcaceae, Veillonellaceae, and Syntrophomonadaceae decreased in relative mRNA abundance. Among methanogens, Methanosarcinaceae were most dominant. Relative to other family-level groups, salt stress induced a significant enrichment of transcripts related to the CO dehydrogenase/acetyl-coenzyme A synthase complex, methanogenesis, heat shock, ammonium uptake, and thermosomes, but the absolute abundance of methanosarcinal mRNA decreased. Most strikingly, the transcriptional activity of the Methanocellaceae was completely suppressed already at 300 mM NaCl. Apparently, the key taxonomic groups involved in the methanogenic breakdown of plant polymers significantly differ in their ability to cope with severe salt stress. Presumably, this different ability is directly linked to differences in their genetic potential and metabolic flexibility to reassign available energy resources for cellular adaptation to salt stress.

Bacterial Human Virulence Genes across Diverse Habitats As Assessed by In silico Analysis of Environmental Metagenomes

The occurrence and distribution of clinically relevant bacterial virulence genes across natural (non-human) environments is not well understood. We aimed to investigate the occurrence of homologs to bacterial human virulence genes in a variety of ecological niches to better understand the role of natural environments in the evolution of bacterial virulence. Twenty four bacterial virulence genes were analyzed in 46 diverse environmental metagenomic datasets, representing various soils, seawater, freshwater, marine sediments, hot springs, the deep-sea, hypersaline mats, microbialites, gutless worms and glacial ice. Homologs to 16 bacterial human virulence genes, involved in urinary tract infections, gastrointestinal diseases, skin diseases, and wound and systemic infections, showed global ubiquity. A principal component analysis did not demonstrate clear trends across the metagenomes with respect to occurrence and frequency of observed gene homologs. Full-length (>95%) homologs of several virulence genes were identified, and translated sequences of the environmental and clinical genes were up to 50-100% identical. Furthermore, phylogenetic analyses indicated deep branching positions of some of the environmental gene homologs, suggesting that they represent ancient lineages in the phylogeny of the clinical genes. Fifteen virulence gene homologs were detected in metatranscriptomes, providing evidence of environmental expression. The ubiquitous presence and transcription of the virulence gene homologs in non-human environments point to an important ecological role of the genes for the activity and survival of environmental bacteria. Furthermore, the high degree of sequence conservation between several of the environmental and clinical genes suggests common ancestral origins.

Marker genes that are less conserved in their sequences are useful for predicting genome-wide similarity levels between closely related prokaryotic strains

BACKGROUND

The 16s rRNA gene is so far the most widely used marker for taxonomical classification and separation of prokaryotes. Since it is universally conserved among prokaryotes, it is possible to use this gene to classify a broad range of prokaryotic organisms. At the same time, it has often been noted that the 16s rRNA gene is too conserved to separate between prokaryotes at finer taxonomic levels.

RESULTS

In this paper, we examine how well levels of similarity of 16s rRNA and 73 additional universal or nearly universal marker genes correlate with genome-wide levels of gene sequence similarity. We demonstrate that the percent identity of 16s rRNA predicts genome-wide levels of similarity very well for distantly related prokaryotes, but not for closely related ones. In closely related prokaryotes, we find that there are many other marker genes for which levels of similarity are much more predictive of genome-wide levels of gene sequence similarity. Finally, we show that the identities of the markers that are most useful for predicting genome-wide levels of similarity within closely related prokaryotic lineages vary greatly between lineages. However, the most useful markers are always those that are least conserved in their sequences within each lineage.

CONCLUSIONS

Our results show that by choosing markers that are less conserved in their sequences within a lineage of interest, it is possible to better predict genome-wide gene sequence similarity between closely related prokaryotes than is possible using the 16s rRNA gene. We point readers towards a database we have created (POGO-DB) that can be used to easily establish which markers show lowest levels of sequence conservation within different prokaryotic lineages.

Yucatán in black and red: Linking edaphic analysis and pyrosequencing-based assessment of bacterial and fungal community structures in the two main kinds of soil of Yucatán State

Yucatán State is dominated by two kinds of soil, named "Black Leptosol" and "Red Leptosol", which are interwoven across the State. In this work, we analyzed the relation between the edaphic characteristics and the bacterial and fungal community structures in these two kinds of Leptosol. The results revealed that Black Leptosol (BlaS) had a higher content of calcium carbonates, organic matter, nitrogen, and phosphorus than Red Leptosol (RedS). The most outstanding difference in the bacterial community structure between BlaS and RedS was that while in BlaS Actinobacteria was the most abundant phylum (43.7%), followed by Acidobacteria (26.9%) and Proteobacteria (23.6%), in RedS the bacterial community was strongly dominated by Acidobacteria (83%). Two fungal phyla were identified in both kinds of soil; Ascomycota, with 77% in BlaS and 56% in RedS, and Basidiomycota, with 22% in RedS and only 0.67% in BlaS. The most relevant difference between the two fungal communities was that excepting for Fusarium sp., all the species they had were different. Thus, in contrast with bacterial communities, where most of the major OTUs were present in both kinds of soil, fungal communities appeared to be unique to each kind of Leptosol.