Diversity of 16S rRNA genes within individual prokaryotic genomes

Genomic signatures of cold adaptation in the family Colwelliaceae

Psychrophily is a phenotype describing microbial growth at low temperatures; elucidating the biomolecular and genomic adaptations necessary for survival in the cold is important for understanding life in extreme environments on Earth and in outer space. We used comparative genomics and temperature growth experiments of bacteria from the family Colwelliaceae to identify genomic factors correlated with optimal growth temperature (OGT). A phylogenomic analysis of 67 public and 39 newly sequenced strains revealed three main clades of Colwelliaceae. Temperature growth experiments revealed significant differences in mean OGT by clade, wherein strains of Colwelliaceae had similar growth rates at -1 °C but varied in their ability to tolerate 17 °C. Using amino acid compositional indices, a multiple linear regression model was constructed to predict the OGT of these organisms (RMSE 5.2 °C). Investigation of Colwelliaceae functional genes revealed a putative cold-adaptive gene cassette that was present in psychrophilic strains but absent in a closely related strain with a significantly higher OGT. This study also presents genomic evidence suggesting that the clade of Colwelliaceae containing Colwellia hornerae should be investigated as a new genus. These contributions offer key insights into the psychrophily phenotype and its underlying genomic foundation in the family Colwelliaceae.

Pseudomonas fragariae sp. nov., a novel bacterial species causing leaf spots on strawberry (Fragaria×ananassa)

In Florida, angular leaf spot, caused by Xanthomonas fragariae, was the only known bacterial disease in strawberry, which is sporadic and affects the foliage and calyx. However, from the 2019-2020 to 2023-2024 Florida strawberry seasons, unusual bacterial-like symptoms were observed in commercial farms, with reports of up to 30 % disease incidence. Typical lesions were water-soaked and angular in early stages that later became necrotic with a circular-ellipsoidal purple halo, and consistently yielded colonies resembling Pseudomonas on culture media. Strains were pathogenic on strawberry, fluorescent, oxidase- and arginine-dihydrolase-negative, elicited a hypersensitive reaction on tobacco, and lacked pectolytic activity. Although phenotypic assays, such as fatty acid methyl profiles and Biolog protocols, placed the strains into the Pseudomonas group, there was a low similarity at the species level. Further analysis using 16S rRNA genes, housekeeping genes, and whole genome sequencing showed that the strains cluster into the Pseudomonas group but do not share more than 95 % average nucleotide identity compared to representative members. Therefore, the genomic and phenotypic analysis confirm that the strains causing bacterial spot in strawberry represent a new plant pathogenic bacterial species for which we propose the name Pseudomonas fragariae sp. nov. with 20-417T (17T=LMG 32456T=DSM 113340 T) as the type strain, in relation to Fragaria×ananassa, the plant species from which the pathogen was first isolated. Future work is needed to assess the epidemiology, cultivar susceptibility, chemical sensitivity, and disease management of this possible new emerging strawberry pathogen.

The past, present and future of polymicrobial infection research: Modelling, eavesdropping, terraforming and other stories

Over the last two centuries, great advances have been made in microbiology as a discipline. Much of this progress has come about as a consequence of studying the growth and physiology of individual microbial species in well-defined laboratory media; so-called "axenic growth". However, in the real world, microbes rarely live in such "splendid isolation" (to paraphrase Foster) and more often-than-not, share the niche with a plethora of co-habitants. The resulting interactions between species (and even between kingdoms) are only very poorly understood, both on a theoretical and experimental level. Nevertheless, the last few years have seen significant progress, and in this review, we assess the importance of polymicrobial infections, and show how improved experimental traction is advancing our understanding of these. A particular focus is on developments that are allowing us to capture the key features of polymicrobial infection scenarios, especially as those associated with the human airways (both healthy and diseased).

Pangenome Analysis of Helicobacter pylori Isolates from Selected Areas of Africa Indicated Diverse Antibiotic Resistance and Virulence Genes

The challenge facing Helicobacter pylori (H. pylori) infection management in some parts of Africa is the evolution of drug-resistant species, the lack of gold standard in diagnostic methods, and the ineffectiveness of current vaccines against the bacteria. It is being established that even though clinical consequences linked to the bacteria vary geographically, there is rather a generic approach to treatment. This situation has remained problematic in the successful fight against the bacteria in parts of Africa. As a result, this study compared the genomes of selected H. pylori isolates from selected areas of Africa and evaluated their virulence and antibiotic drug resistance, those that are highly pathogenic and are associated with specific clinical outcomes and those that are less virulent and rarely associated with clinical outcomes. 146 genomes of H. pylori isolated from selected locations of Africa were sampled, and bioinformatic tools such as Abricate, CARD RGI, MLST, Prokka, Roary, Phandango, Google Sheets, and iTOLS were used to compare the isolates and their antibiotic resistance or susceptibility. Over 20 k virulence and AMR genes were observed. About 95% of the isolates were genetically diverse, 90% of the isolates harbored shell genes, and 50% harbored cloud and core genes. Some isolates did not retain the cagA and vacA genes. Clarithromycin, metronidazole, amoxicillin, and tinidazole were resistant to most AMR genes (vacA, cagA, oip, and bab). Conclusion. This study found both virulence and AMR genes in all H. pylori strains in all the selected geographies around Africa with differing quantities. MLST, Pangenome, and ORF analyses showed disparities among the isolates. This in general could imply diversities in terms of genetics, evolution, and protein production. Therefore, generic administration of antibiotics such as clarithromycin, amoxicillin, and erythromycin as treatment methods in the African subregion could be contributing to the spread of the bacterium's antibiotic resistance.

Examination of Large Chromosomal Inversions in the Genome of Erwinia amylovora Strains Reveals Worldwide Distribution and North America-Specific Types

Erwinia amylovora is a relatively homogeneous species with low genetic diversity at the nucleotide level. However, phenotypic differences and genomic structural variations among E. amylovora strains have been documented. In this study, we identified 10 large chromosomal inversion (LCI) types in the Spiraeoideae-infecting (SI) E. amylovora strains by combining whole genome sequencing and PCR-based molecular markers. It was found that LCIs were mainly caused by homologous recombination events among seven rRNA operons (rrns) in SI E. amylovora strains. Although ribotyping results identified inter- and intra-variations in the internal transcribed spacer (ITS1 and ITS2) regions among rrns, LCIs tend to occur between rrns transcribed in the opposite directions and with the same tRNA content (tRNA-Glu or tRNA-Ile/Ala) in ITS1. Based on the LCI types, physical/estimated replichore imbalance (PRI/ERI) was examined and calculated. Among the 117 SI strains evaluated, the LCI types of Ea1189, CFBP1430, and Ea273 were the most common, with ERI values at 1.31, 7.87, and 4.47°, respectively. These three LCI types had worldwide distribution, whereas the remaining seven LCI types were restricted to North America (or certain regions of the United States). Our results indicated ongoing chromosomal recombination events in the SI E. amylovora population and showed that LCI events are mostly symmetrical, keeping the ERI less than 15°. These findings provide initial evidence about the prevalence of certain LCI types in E. amylovora strains, how LCI occurs, and its potential evolutionary advantage and history, which might help track the movement of the pathogen.

Genomic insights and comparative analysis of Flavobacterium bizetiae HJ-32-4 isolated from soil

In the late 1970s, Flavobacterium bizetiae was first isolated from diseased fish in Canada. After four decades of preservation, it was reported as a novel species in 2020. Here, we report the first complete genome sequence of HJ-32-4, a novel strain of F. bizetiae. Interestingly, HJ-32-4 was isolated from soil in Gangwon-do, Republic of Korea, unlike the other two previously reported F. bizetiae strains which were isolated from fish. We generated a single circular chromosome of HJ-32-4, comprising 5,745,280 bp with a GC content of 34.2%. The average nucleotide identity (ANI) value of 96.2% indicated that HJ-32-4 belongs to F. bizetiae CIP 105534T. The virulence factor was not detected in the genome. Comparative genomic analysis of F. bizetiae and major flavobacterial pathogens revealed that F. bizetiae had a larger genome size and the ratio of peptidases (PEP) and glycoside hydrolase (GH) genes of F. bizetiae were lower than those of the rest strains, implying that F. bizetiae exhibits similar characteristics with non-pathogenic strains from a genomic point of view. However, further experimental verification is required to ensure these in silico predictions. This study will provide insight into the overall characteristics of HJ-32-4 compared to other strains.

UACG: Up-to-Date Archaeal Core Genes and Software for Phylogenomic Tree Reconstruction

In the post-genomic era, phylogenomics is a powerful and routinely-used tool to discover evolutionary relationships between microorganisms. Inferring phylogenomic trees by concatenating core gene sequences into a supermatrix is the standard method. The previously released up-to-date bacterial core gene (UBCG) tool provides a pipeline to infer phylogenomic trees using single-copy core genes for the Bacteria domain. In this study, we established up-to-date archaeal core gene (UACG), comprising 128 genes suitable for inferring archaeal phylogenomic trees. To test the gene set, we selected the Haloarcula genus and scrutinized its phylogeny. The phylogeny inferred using the UACG tool was consistent with the orthoANIu dendrogram, whereas the 16S rRNA gene phylogeny showed high intragenomic heterogeneity resulting in phylogenetic discrepancies. The software tool using the UACG set is available at https://www.ezbiocloud.net/tools/uacg .

Ensifer canadensis sp. nov. strain T173T isolated from Melilotus albus (sweet clover) in Canada possesses recombinant plasmid pT173b harbouring symbiosis and type IV secretion system genes apparently acquired from Ensifer medicae

A bacterial strain, designated T173^T, was previously isolated from a root-nodule of a Melilotus albus plant growing in Canada and identified as a novel Ensifer lineage that shared a clade with the non-symbiotic species, Ensifer adhaerens. Strain T173^T was also previously found to harbour a symbiosis plasmid and to elicit root-nodules on Medicago and Melilotus species but not fix nitrogen. Here we present data for the genomic and taxonomic description of strain T173^T. Phylogenetic analyses including the analysis of whole genome sequences and multiple locus sequence analysis (MLSA) of 53 concatenated ribosome protein subunit (rps) gene sequences confirmed placement of strain T173^T in a highly supported lineage distinct from named Ensifer species with E. morelensis Lc04^T as the closest relative. The highest digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) values of genome sequences of strain T173^T compared with closest relatives (35.7 and 87.9%, respectively) are well below the respective threshold values of 70% and 95-96% for bacterial species circumscription. The genome of strain T173^T has a size of 8,094,229 bp with a DNA G + C content of 61.0 mol%. Six replicons were detected: a chromosome (4,051,102 bp) and five plasmids harbouring plasmid replication and segregation (repABC) genes. These plasmids were also found to possess five apparent conjugation systems based on analysis of TraA (relaxase), TrbE/VirB4 (part of the Type IV secretion system (T4SS)) and TraG/VirD4 (coupling protein). Ribosomal RNA operons encoding 16S, 23S, and 5S rRNAs that are usually restricted to bacterial chromosomes were detected on plasmids pT173d and pT173e (946,878 and 1,913,930 bp, respectively) as well as on the chromosome of strain T173^T. Moreover, plasmid pT173b (204,278 bp) was found to harbour T4SS and symbiosis genes, including nodulation (nod, noe, nol) and nitrogen fixation (nif, fix) genes that were apparently acquired from E. medicae by horizontal transfer. Data for morphological, physiological and symbiotic characteristics complement the sequence-based characterization of strain T173^T. The data presented support the description of a new species for which the name Ensifer canadensis sp. nov. is proposed with strain T173^T (= LMG 32374^T = HAMBI 3766^T) as the species type strain.

Diagnosis and Treatment of Gordonia Species Infection in a Peach-Faced Lovebird (Agapornis roseicollis)

Respiratory distress is a common presentation for avian species. A 9-week-old peach-faced lovebird (Agapornis roseicollis) was presented with a 2-week history of progressive dyspnea. Computed tomographic (CT) images were suggestive of splenomegaly and bilateral granulomatous pulmonary disease. Polymerase chain reaction testing of samples from the choana, cloaca, and distal tracheal/syringeal area were positive for Mycobacterium species hsp65. A comparison search of the 400 base pair sequence in the NCBI/BLAST/blastn database revealed a best match of 93% similarity to Gordonia species and 91% similarity to Gordonia bronchialis. Gordonia is a genus in the phylum Actinomycetota, the same lineage that includes Mycobacterium species. Gordonia species can be mistaken for Mycobacterium species unless more definitive diagnostic testing is pursued. Infection caused by Gordonia species is rare in humans. Reports commonly cite infection of immunocompromised patients, and to our knowledge, no reports of treatment have been published in the veterinary literature. After the test results were obtained, the patient was treated with azithromycin and pradofloxacin for 3 months. The lovebird was presented for reexamination when the antibiotic treatment was complete. When reexamined, and a second series of CT images evaluated, it was determined that the treatment achieved clinical resolution of signs and lesions.

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.

Evaluation of DNA extraction methods and direct PCR in metabarcoding of mock and marine bacterial communities

Recent advances in new molecular biology methods and next-generation sequencing (NGS) technologies have revolutionized metabarcoding studies investigating complex microbial communities from various environments. The inevitable first step in sample preparation is DNA extraction which introduces its own set of biases and considerations. In this study, we assessed the influence of five DNA extraction methods [B1: phenol/chloroform/isoamyl extraction, B2 and B3: isopropanol and ethanol precipitations, respectively-both modifications of B1, K1: DNeasy PowerWater Kit (QIAGEN), K2: modified DNeasy PowerWater Kit (QIAGEN) and direct PCR approach (P) that completely circumvents this step on community composition and DNA yield of mock and marine sample communities from the Adriatic Sea]. B1-B3 methods generally produced higher DNA yields and more similar microbial communities, but with higher interindividual variability. Each method demonstrated significant differences in a specific community structure, where rare taxa seem to play a crucial role. There was not one superior method closest to the theoretically expected mock community composition, they all demonstrated skewed ratios, but in a similar way which might be attributed to other factors, such as primer bias or 16S rRNA gene count for specific taxa. Direct PCR represents an interesting approach when high throughput in sample processing is required. We emphasize the importance of making a cautious decision about the choice of the extraction method or direct PCR approach, but even more importantly its consistent application throughout the study.

Extraction of High-Quality RNA from S. aureus Internalized by Endothelial Cells

Staphylococcus aureus evades antibiotic therapy and antimicrobial defenses by entering human host cells. Bacterial transcriptomic analysis represents an invaluable tool to unravel the complex interplay between host and pathogen. Therefore, the extraction of high-quality RNA from intracellular S. aureus lays the foundation to acquire meaningful gene expression data. In this study, we present a novel and straightforward strategy to isolate RNA from internalized S. aureus after 90 min, 24 h, and 48 h postinfection. Real-time PCR data were obtained for the target genes agrA and fnba, which play major roles during infection. The commonly used reference genes gyrB, aroE, tmRNA, gmk, and hu were analyzed under different conditions: bacteria from culture (condition I), intracellular bacteria (condition II), and across both conditions I and II. The most stable reference genes were used for the normalization of agrA and fnbA. Delta C_q (quantification cycle) values had a relatively low variability and thus demonstrated the high quality of the extracted RNA from intracellular S. aureus during the early phase of infection. The established protocol allows the extraction and purification of intracellular staphylococcal RNA while minimizing the amount of host RNA in the sample. This approach can leverage reproducible gene expression data to study host-pathogen interactions.

Mock community as an in situ positive control for amplicon sequencing of microbiotas from the same ecosystem

Metataxonomy has become the standard for characterizing the diversity and composition of microbial communities associated with multicellular organisms and their environment. Currently available protocols for metataxonomy assume a uniform DNA extraction, amplification and sequencing efficiency for all sample types and taxa. It has been suggested that the addition of a mock community (MC) to biological samples before the DNA extraction step could aid identification of technical biases during processing and support direct comparisons of microbiota composition, but the impact of MC on diversity estimates of samples is unknown. Here, large and small aliquots of pulverized bovine fecal samples were extracted with no, low or high doses of MC, characterized using standard Illumina technology for metataxonomics, and analysed with custom bioinformatic pipelines. We demonstrated that sample diversity estimates were distorted only if MC dose was high compared to sample mass (i.e. when MC > 10% of sample reads). We also showed that MC was an informative in situ positive control, permitting an estimation of the sample 16S copy number, and detecting sample outliers. We tested this approach on a range of sample types from a terrestrial ecosystem, including rhizosphere soil, whole invertebrates, and wild vertebrate fecal samples, and discuss possible clinical applications.

Optimization of Molecular Methods for Detecting Duckweed-Associated Bacteria

The bacterial colonization dynamics of plants can differ between phylogenetically similar bacterial strains and in the context of complex bacterial communities. Quantitative methods that can resolve closely related bacteria within complex communities can lead to a better understanding of plant-microbe interactions. However, current methods often lack the specificity to differentiate phylogenetically similar bacterial strains. In this study, we describe molecular strategies to study duckweed-associated bacteria. We first systematically optimized a bead-beating protocol to co-isolate nucleic acids simultaneously from duckweed and bacteria. We then developed a generic fingerprinting assay to detect bacteria present in duckweed samples. To detect specific duckweed-bacterium associations, we developed a genomics-based computational pipeline to generate bacterial strain-specific primers. These strain-specific primers differentiated bacterial strains from the same genus and enabled the detection of specific duckweed-bacterium associations present in a community context. Moreover, we used these strain-specific primers to quantify the bacterial colonization of duckweed by normalization to a plant reference gene and revealed differences in colonization levels between strains from the same genus. Lastly, confocal microscopy of inoculated duckweed further supported our PCR results and showed bacterial colonization of the duckweed root-frond interface and root interior. The molecular methods introduced in this work should enable the tracking and quantification of specific plant-microbe associations within plant-microbial communities.

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.

StrainSelect: A novel microbiome reference database that disambiguates all bacterial strains, genome assemblies and extant cultures worldwide

Motivation: Microbial metagenomic profiling software and databases are advancing rapidly for development of novel disease biomarkers and therapeutics yet three problems impede analyses: 1) the conflation of "genome assembly" and "strain" in reference databases; 2) difficulty connecting DNA biomarkers to a procurable strain for laboratory experimentation; and 3) absence of a comprehensive and unified strain-resolved reference database for integrating both shotgun metagenomics and 16S rRNA gene data. Results: We demarcated 681,087 strains, the largest collection of its kind, by filtering public data into a knowledge graph of vertices representing contiguous DNA sequences, genome assemblies, strain monikers and bio-resource center (BRC) catalog numbers then adding inter-vertex edges only for synonyms or direct derivatives. Surprisingly, for 10,043 important strains, we found replicate RefSeq genome assemblies obstructing interpretation of database searches. We organized each strain into eight taxonomic ranks with bootstrap confidence inversely correlated with genome assembly contamination. The StrainSelect database is suited for applications where a taxonomic, functional or procurement reference is needed for shotgun or amplicon metagenomics since 636,568 strains have at least one 16S rRNA gene, 245,005 have at least one annotated genome assembly, and 36,671 are procurable from at least one BRC. The database overcomes all three aforementioned problems since it disambiguates strains from assemblies, locates strains at BRCs, and unifies a taxonomic reference for both 16S rRNA and shotgun metagenomics. Availability: The StrainSelect database is available in igraph and tabular vertex-edge formats compatible with Neo4J. Dereplicated MinHash and fasta databases are distributed for sourmash and usearch pipelines at http://strainselect.secondgenome.com. Contact:todd.desantis@gmail.com. Supplementary information: Supplementary data are available online.

Spatial and Seasonal Patterns of Sediment Bacterial Communities in Large River Cascade Reservoirs: Drivers, Assembly Processes, and Co-occurrence Relationship

Sediment bacteria play an irreplaceable role in promoting the function and biogeochemical cycle of the freshwater ecosystem; however, little is known about their biogeographical patterns and community assembly mechanisms in large river suffering from cascade development. Here, we investigated the spatiotemporal distribution patterns of bacterial communities employing next-generation sequencing analysis and multivariate statistical analyses from the Lancang River cascade reservoirs during summer and winter. We found that sediment bacterial composition has a significant seasonal turnover due to the modification of cascade reservoirs operation mode, and the spatial consistency of biogeographical models (including distance-decay relationship and covariation of community composition with geographical distance) also has subtle changes. The linear regression between the dissimilarity of bacterial communities in sediments, geographical and environmental distance showed that the synergistic effects of geographical and environmental factors explained the influence on bacterial communities. Furthermore, the environmental difference explained little variations (19.40%) in community structure, implying the homogeneity of environmental conditions across the cascade reservoirs of Lancang River. From the quantification of the ecological process, the homogeneous selection was recognized as the dominating factor of bacterial community assembly. The co-occurrence topological network analyses showed that the key genera were more important than the most connected genera. In general, the assembly of bacterial communities in sediment of cascade reservoirs was mediated by both deterministic and stochastic processes and is always dominated by homogeneous selection with the seasonal switching, but the effects of dispersal limitation and ecological drift cannot be ignored.

The Microbiome's Influence on Head and Neck Cancers

PURPOSE OF REVIEW

Head and neck tumors (HNC) rank sixth among cancers worldwide. Due to their late diagnosis and poor prognosis, they are a clinical challenge. However, recent years have seen a dynamic development of science on the microbiome. The aim of the study is to discuss the role of the microbiome in HNC, the impact of the microbiome on oncogenesis, the course of the disease, as well as on treatment, and its toxicity.

RECENT FINDINGS

The microbiome's influence on oncogenesis, the course of the disease, and the effectiveness of oncological treatment have been confirmed in cancers of the colon, pancreas, lungs, and prostate. There is no solid literature on HNC. Many studies indicate disruption of the oral microbiome and periodontal disease as potential cancer risk factors. Disruption of the microbiome increases radiotherapy's toxicity, intensifying radiation reactions. The microbiome plays an important role in cancer. It is a new target in research into new therapies. It may also be a prognostic marker of cancer development. Changes in the composition of the microbiome modulate the effectiveness of oncological treatment. More research is needed on the microbiome and its effects on HNC.

Oral Microbiome in Nonsmoker Patients with Oral Cavity Squamous Cell Carcinoma, Defined by Metagenomic Shotgun Sequencing

Objectives: Smoking is the commonest cause of oral cavity squamous cell carcinoma (OC-SCC), but the etiology of OC-SCC in nonsmokers is unknown. Our primary goal was to use metagenomic shotgun sequencing (MSS) to define the taxonomic composition and functional potential of oral metagenome in nonsmokers with OC-SCC. Methods: We conducted a case-control study with 42 OC-SCC case and 45 control nonsmokers. MSS was performed on DNA extracted from mouthwash samples. Taxonomic analysis and pathway analysis were done using MetaPhlAn2 and HUMAnN2, respectively. Statistical difference was determined using the Mann-Whitney test controlling false discovery rate. Results: There was no significant difference in age, sex, race, or alcohol consumption between OC-SCC and control patients. There was a significant difference in beta diversity between OC-SCC and controls. At the phylum level, Bacteroidetes and Synergistetes were overly represented in OC-SCC while Actinobacteria and Firmicutes were overly represented in controls. At the genus level, Fusobacterium was overly represented in OC-SCC compared with controls, while Corynebacterium, Streptococcus, Actinomyces, Cryptobacterium, and Selenomonas were overly represented in controls. Bacterial pathway analysis identified overrepresentation in OC-SCC of pathways related to metabolism of flavin, biotin, thiamin, heme, sugars, fatty acids, peptidoglycans, and tRNA and overrepresentation of nucleotides and essential amino acids in controls. Conclusions: The oral microbiome in nonsmoker patients with OC-SCC is significantly different from that of nonsmoker control patients in taxonomic compositions and functional potentials. Our study's MSS findings matched with previous 16S-based methods in taxonomic differentiation but varied greatly in functional differentiation of microbiomes in OC-SCC and controls.

Pseudomonas californiensis sp. nov. and Pseudomonas quasicaspiana sp. nov., isolated from ornamental crops in California

Five bacterial strains were isolated from symptomatic leaves of Achillea millefolium, Delphinium sp. and Hydrangea sp. in California. Colonies isolated on King’s medium B (KMB) appeared white, mucoid and round, similar to Pseudomonas species. Phylogenetic analyses based on 16S rRNA, rpoB, rpoD and gyrB genes placed the bacteria into three distinct groups within Pseudomonas that were most closely related to Pseudomonas viridiflava , Pseudomonas cichorii or Pseudomonas caspiana . To further characterize the strains, phenotypic analyses and the following tests were performed: fatty acid methyl ester composition, LOPAT, fluorescence on KMB, Biolog assay, and transmission electron microscopy. Finally, whole genome sequencing of the strains was conducted, and the sequences were compared with reference genomes of Pseudomonas species based on average nucleotide identity (ANI). The first group, which consists of three strains isolated from delphinium, hydrangea and achillea, had 95.6–96.9 % pairwise ANI between each other; the second group consists of two strains isolated from delphinium that had 100 % pairwise ANI. Although comparisons of the two groups with publicly available genomes revealed closest relationships with P. viridiflava (91.6 %), P. caspiana (88.3 %) and P. asturiensis (86.7 %), ANI values were less than 95 % compared to all validly published pseudomonads. Combining genomic and phenotypic data, we conclude that these strains represent two new species and the names proposed are Pseudomonas quasicaspiana sp. nov. (type strain DSMZ 11 30 42T=LMG 32 434T) for the strains isolated from delphinium, achillea and hydrangea and Pseudomonas californiensis sp. nov. (DSMZ 11 30 43T=LMG 32 432T) for the two strains isolated from delphinium. The specific epithets quasicaspiana and californiensis were selected based on the close phylogenetic relationship of strains with P. caspiana and on the geographic location of isolation, respectively.

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.

Multi-omics analysis: Paving the path toward achieving precision medicine in cancer treatment and immuno-oncology

The acceleration of large-scale sequencing and the progress in high-throughput computational analyses, defined as omics, was a hallmark for the comprehension of the biological processes in human health and diseases. In cancerology, the omics approach, initiated by genomics and transcriptomics studies, has revealed an incredible complexity with unsuspected molecular diversity within a same tumor type as well as spatial and temporal heterogeneity of tumors. The integration of multiple biological layers of omics studies brought oncology to a new paradigm, from tumor site classification to pan-cancer molecular classification, offering new therapeutic opportunities for precision medicine. In this review, we will provide a comprehensive overview of the latest innovations for multi-omics integration in oncology and summarize the largest multi-omics dataset available for adult and pediatric cancers. We will present multi-omics techniques for characterizing cancer biology and show how multi-omics data can be combined with clinical data for the identification of prognostic and treatment-specific biomarkers, opening the way to personalized therapy. To conclude, we will detail the newest strategies for dissecting the tumor immune environment and host–tumor interaction. We will explore the advances in immunomics and microbiomics for biomarker identification to guide therapeutic decision in immuno-oncology.

RiboTaxa: combined approaches for rRNA genes taxonomic resolution down to the species level from metagenomics data revealing novelties

Metagenomic classifiers are widely used for the taxonomic profiling of metagenomics data and estimation of taxa relative abundance. Small subunit rRNA genes are a gold standard for phylogenetic resolution of microbiota, although the power of this marker comes down to its use as full-length. We aimed at identifying the tools that can efficiently lead to taxonomic resolution down to the species level. To reach this goal, we benchmarked the performance and accuracy of rRNA-specialized versus general-purpose read mappers, reference-targeted assemblers and taxonomic classifiers. We then compiled the best tools (BBTools, FastQC, SortMeRNA, MetaRib, EMIRGE, VSEARCH, BBMap and QIIME 2’s Sklearn classifier) to build a pipeline called RiboTaxa. Using metagenomics datasets, RiboTaxa gave the best results compared to other tools (i.e. Kraken2, Centrifuge, METAXA2, phyloFlash, SPINGO, BLCA, MEGAN) with precise taxonomic identification and relative abundance description without false positive detection (F-measure of 100% and 83.7% at genus level and species level, respectively). Using real datasets from various environments (i.e. ocean, soil, human gut) and from different approaches (e.g. metagenomics and gene capture by hybridization), RiboTaxa revealed microbial novelties not discerned by current bioinformatics analysis opening new biological perspectives in human and environmental health.

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.

Distribution, Characterization and the Commercialization of Elite Rhizobia Strains in Africa

Grain legumes play a significant role in smallholder farming systems in Africa because of their contribution to nutrition and income security and their role in fixing nitrogen. Biological Nitrogen Fixation (BNF) serves a critical role in improving soil fertility for legumes. Although much research has been conducted on rhizobia in nitrogen fixation and their contribution to soil fertility, much less is known about the distribution and diversity of the bacteria strains in different areas of the world and which of the strains achieve optimal benefits for the host plants under specific soil and environmental conditions. This paper reviews the distribution, characterization, and commercialization of elite rhizobia strains in Africa.

Microbiome differences between wild and aquarium whitespotted eagle rays (Aetobatus narinari)

Background

Animal-associated microbiomes can be influenced by both host and environmental factors. Comparing wild animals to those in zoos or aquariums can help disentangle the effects of host versus environmental factors, while also testing whether managed conditions foster a ‘natural’ host microbiome. Focusing on an endangered elasmobranch species—the whitespotted eagle ray Aetobatus narinari—we compared the skin, gill, and cloaca microbiomes of wild individuals to those at Georgia Aquarium. Whitespotted eagle ray microbiomes from Georgia Aquarium were also compared to those of cownose rays (Rhinoptera bonasus) in the same exhibit, allowing us to explore the effect of host identity on the ray microbiome.

Results

Long-term veterinary monitoring indicated that the rays in managed care did not have a history of disease and maintained health parameters consistent with those of wild individuals, with one exception. Aquarium whitespotted eagle rays were regularly treated to control parasite loads, but the effects on animal health were subclinical. Microbiome α- and β-diversity differed between wild versus aquarium whitespotted eagle rays at all body sites, with α-diversity significantly higher in wild individuals. β-diversity differences in wild versus aquarium whitespotted eagle rays were greater for skin and gill microbiomes compared to those of the cloaca. At each body site, we also detected microbial taxa shared between wild and aquarium eagle rays. Additionally, the cloaca, skin, and gill microbiomes of aquarium eagle rays differed from those of cownose rays in the same exhibit. Potentially pathogenic bacteria were at low abundance in all wild and aquarium rays.

Conclusion

For whitespotted eagle rays, managed care was associated with a microbiome differing significantly from that of wild individuals. These differences were not absolute, as the microbiome of aquarium rays shared members with that of wild counterparts and was distinct from that of a cohabitating ray species. Eagle rays under managed care appear healthy, suggesting that their microbiomes are not associated with compromised host health. However, the ray microbiome is dynamic, differing with both environmental factors and host identity. Monitoring of aquarium ray microbiomes over time may identify taxonomic patterns that co-vary with host health.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-022-00187-8.

New Insights into the Taxonomy of Bacteria in the Genomic Era and a Case Study with Rhizobia

Since early studies, the history of prokaryotes taxonomy has dealt with many changes driven by the development of new and more robust technologies. As a result, the number of new taxa descriptions is exponentially increasing, while an increasing number of others has been subject of reclassification, demanding from the taxonomists more effort to maintain an organized hierarchical system. However, expectations are that the taxonomy of prokaryotes will acquire a more stable status with the genomic era. Other analyses may continue to be necessary to determine microbial features, but the use of genomic data might be sufficient to provide reliable taxa delineation, helping taxonomy to reach the goal of correct classification and identification. Here we describe the evolution of prokaryotes' taxonomy until the genomic era, emphasizing bacteria and taking as an example the history of rhizobia taxonomy. This example was chosen because of the importance of the symbiotic nitrogen fixation of legumes with rhizobia to the nitrogen input to both natural ecosystems and agricultural crops. This case study reports the technological advances and the methodologies used to classify and identify bacterial species and indicates the actual rules required for an accurate description of new taxa.

RiboGrove: a database of full-length prokaryotic 16S rRNA genes derived from completely assembled genomes

16S rRNA gene is frequently used for the identification of prokaryotic organisms and for phylogeny inference. Several specialized public databases exist that contain complete and partial sequences of 16S rRNA genes. In this paper, we present RiboGrove: the first publicly available database that comprises only full-length sequences of 16S rRNA genes originating from completely assembled prokaryotic genomes deposited in RefSeq. Despite being strongly biased towards frequently sequenced genomes, RiboGrove is a useful complement to existing 16S rRNA resources and allows for analyses that would not be possible using amplicon-derived gene sequences. For instance, the absence of partial gene sequences in RiboGrove allowed us to make a summary of prokaryotic organisms, which lack core anti-Shine-Dalgarno sequence in their 16S rRNA genes. In this study, we describe the collected sequence data and present the results of exploratory data analysis of 16S rRNA gene sequences.

Gut microbiome alterations in the crustacean Pacifastacus leniusculus exposed to environmental concentrations of antibiotics and effects on susceptibility to bacteria challenges

Gut-associated microbiota in crustaceans are recognized as a key element for maintaining homeostasis and health in the animal. Since the richness of these microbial communities is strongly influenced by the local environment, especially in aquatic organisms, it is important to address to what extent environmental variations can affect these communities. In the present study, we used high-throughput 16S rRNA sequencing technology to study the composition of gut-associated microbiota of the crayfish Pacifastacus leniusculus after exposure to environmentally-relevant concentrations of an antibiotic, namely sulfamethoxazole. Also, we examined if alterations of microbiota caused by environmentally-relevant concentrations of this antibiotic affected the host susceptibility to bacterial diseases, including Vibrio species. As a result, we found high individual variability of bacterial abundance and composition in the intestinal microbiome of crayfish, in both antibiotic-exposed and antibiotic-free crayfish. However, an increase of chitinolytic bacteria including Vibrio spp. was detected in some animals exposed to the antibiotic. Moreover, when crayfish susceptibility to bacterial infections was tested, the antibiotic-exposed crayfish survived longer than the control crayfish group. This study represents the first approach for investigating the interplay between crayfish and intestinal bacteria during antibiotic-pollution scenarios. Results herein should be considered by scientists before planning experiments under laboratory conditions, especially to study environmental effects on aquatic animals' intestinal health and immune status.

Improved detection of disease-associated gut microbes using 16S sequence-based biomarkers

Background

Sequencing partial 16S rRNA genes is a cost effective method for quantifying the microbial composition of an environment, such as the human gut. However, downstream analysis relies on binning reads into microbial groups by either considering each unique sequence as a different microbe, querying a database to get taxonomic labels from sequences, or clustering similar sequences together. However, these approaches do not fully capture evolutionary relationships between microbes, limiting the ability to identify differentially abundant groups of microbes between a diseased and control cohort. We present sequence-based biomarkers (SBBs), an aggregation method that groups and aggregates microbes using single variants and combinations of variants within their 16S sequences. We compare SBBs against other existing aggregation methods (OTU clustering and Microphenoor DiTaxa features) in several benchmarking tasks: biomarker discovery via permutation test, biomarker discovery via linear discriminant analysis, and phenotype prediction power. We demonstrate the SBBs perform on-par or better than the state-of-the-art methods in biomarker discovery and phenotype prediction.

Results

On two independent datasets, SBBs identify differentially abundant groups of microbes with similar or higher statistical significance than existing methods in both a permutation-test-based analysis and using linear discriminant analysis effect size. . By grouping microbes by SBB, we can identify several differentially abundant microbial groups (FDR <.1) between children with autism and neurotypical controls in a set of 115 discordant siblings. Porphyromonadaceae, Ruminococcaceae, and an unnamed species of Blastocystis were significantly enriched in autism, while Veillonellaceae was significantly depleted. Likewise, aggregating microbes by SBB on a dataset of obese and lean twins, we find several significantly differentially abundant microbial groups (FDR<.1). We observed Megasphaera andSutterellaceae highly enriched in obesity, and Phocaeicola significantly depleted. SBBs also perform on bar with or better than existing aggregation methods as features in a phenotype prediction model, predicting the autism phenotype with an ROC-AUC score of .64 and the obesity phenotype with an ROC-AUC score of .84.

Conclusions

SBBs provide a powerful method for aggregating microbes to perform differential abundance analysis as well as phenotype prediction. Our source code can be freely downloaded from http://github.com/briannachrisman/16s_biomarkers.

A deep-sea bacterium related to coastal marine pathogens

Evolution of virulence traits from adaptation to environmental niches other than the host is probably a common feature of marine microbial pathogens, whose knowledge might be crucial to understand their emergence and pathogenetic potential. Here, we report genome sequence analysis of a novel marine bacterial species, Vibrio bathopelagicus sp. nov., isolated from warm bathypelagic waters (3309 m depth) of the Mediterranean Sea. Interestingly, V. bathopelagicus sp. nov. is closely related to coastal Vibrio strains pathogenic to marine bivalves. V. bathopelagicus sp. nov. genome encodes genes involved in environmental adaptation to the deep-sea but also in virulence, such as the R5.7 element, MARTX toxin cluster, Type VI secretion system and zinc-metalloprotease, previously associated with Vibrio infections in farmed oysters. The results of functional in vitro assays on immunocytes (haemocytes) of the Mediterranean mussel Mytilus galloprovincialis and the Pacific oyster Crassostrea gigas, and of the early larval development assay in Mytilus support strong toxicity of V. bathopelagicus sp. nov. towards bivalves. V. bathopelagicus sp. nov., isolated from a remote Mediterranean bathypelagic site, is an example of a planktonic marine bacterium with genotypic and phenotypic traits associated with animal pathogenicity, which might have played an evolutionary role in the origin of coastal marine pathogens.

Amplicon Sequence Variants Artificially Split Bacterial Genomes into Separate Clusters

Amplicon sequencing variants (ASVs) have been proposed as an alternative to operational taxonomic units (OTUs) for analyzing microbial communities. ASVs have grown in popularity, in part because of a desire to reflect a more refined level of taxonomy since they do not cluster sequences based on a distance-based threshold. However, ASVs and the use of overly narrow thresholds to identify OTUs increase the risk of splitting a single genome into separate clusters. To assess this risk, I analyzed the intragenomic variation of 16S rRNA genes from the bacterial genomes represented in an rrn copy number database, which contained 20,427 genomes from 5,972 species. As the number of copies of the 16S rRNA gene increased in a genome, the number of ASVs also increased. There was an average of 0.58 ASVs per copy of the 16S rRNA gene for full-length 16S rRNA genes. It was necessary to use a distance threshold of 5.25% to cluster full-length ASVs from the same genome into a single OTU with 95% confidence for genomes with 7 copies of the 16S rRNA, such as Escherichia coli. This research highlights the risk of splitting a single bacterial genome into separate clusters when ASVs are used to analyze 16S rRNA gene sequence data. Although there is also a risk of clustering ASVs from different species into the same OTU when using broad distance thresholds, these risks are of less concern than artificially splitting a genome into separate ASVs and OTUs.

IMPORTANCE 16S rRNA gene sequencing has engendered significant interest in studying microbial communities. There has been tension between trying to classify 16S rRNA gene sequences to increasingly lower taxonomic levels and the reality that those levels were defined using more sequence and physiological information than is available from a fragment of the 16S rRNA gene. Furthermore, the naming of bacterial taxa reflects the biases of those who name them. One motivation for the recent push to adopt ASVs in place of OTUs in microbial community analyses is to allow researchers to perform their analyses at the finest possible level that reflects species-level taxonomy. The current research is significant because it quantifies the risk of artificially splitting bacterial genomes into separate clusters. Far from providing a better representation of bacterial taxonomy and biology, the ASV approach can lead to conflicting inferences about the ecology of different ASVs from the same genome.

MYCOBACTERIA IN SKIN LESIONS AND THE HABITAT OF THE ENDANGERED HOUSTON TOAD (ANAXYRUS HOUSTONENSIS)

Head-starting of the federally endangered Houston toad (Anaxyrus houstonensis), that is, the release of egg strands, tadpoles, and metamorphic juveniles produced in captivity into the original breeding ponds, requires assessment of potential threats for the transmission of pathogens from captive to free-ranging toads. We used Illumina-based 16S rRNA V3 amplicon sequencing to investigate the community structure of bacteria from skin lesions of captive Houston toad and habitat (pond) samples. Proteobacteria, alone or together with Actinobacteria and, in some samples, Cyanobacteria represented virtually all reads in tissue lesion samples, whereas pond samples were much more diverse, with Acidobacteria, Actinobacteria, Bacteriodetes, Chloroflexi, Cyanobacteria, Firmicutes, Planctomycetes, Proteobacteria, and Verrucomicrobia present with little variation between samples. If present in lesions, Actinobacteria were largely represented by Mycobacteriaceae, and here mainly by one sequence identical to sequences of members of the Mycobacterium chelonae-abscessus complex. In pond samples, mycobacteria represented only a small portion of the actinobacteria, although at higher diversity with six distinct reads. Sequences for reads obtained from pond samples were identical to those representing the M. chelonae-abscessus complex, a group with Mycobacterium marinum, Mycobacterium kansasii, Mycobacterium avium, a group with Mycobacterium vaccae, Mycobacterium fortuitum, Mycobacterium poriferae, and a group with Mycobacterium elephantis and Mycobacterium celeriflavum, whereas sequences of high similarity were detected for reads related to those of Mycobacterium holsaticum, Mycobacterium pallens, and Mycobacterium obuense, and Mycobacterium goodii. Our results indicated that lesions observed on the Houston toad in captivity are not the result of mycobacteria in every case, and that the presence of mycobacteria in the captive colony does not represent a novel pathogen threat to the wild populations because such bacteria are also seen in the natural pond habitats for the Houston toad.

General Unified Microbiome Profiling Pipeline (GUMPP) for Large Scale, Streamlined and Reproducible Analysis of Bacterial 16S rRNA Data to Predicted Microbial Metagenomes, Enzymatic Reactions and Metabolic Pathways

General Unified Microbiome Profiling Pipeline (GUMPP) was developed for large scale, streamlined and reproducible analysis of bacterial 16S rRNA data and prediction of microbial metagenomes, enzymatic reactions and metabolic pathways from amplicon data. GUMPP workflow introduces reproducible data analyses at each of the three levels of resolution (genus; operational taxonomic units (OTUs); amplicon sequence variants (ASVs)). The ability to support reproducible analyses enables production of datasets that ultimately identify the biochemical pathways characteristic of disease pathology. These datasets coupled to biostatistics and mathematical approaches of machine learning can play a significant role in extraction of truly significant and meaningful information from a wide set of 16S rRNA datasets. The adoption of GUMPP in the gut-microbiota related research enables focusing on the generation of novel biomarkers that can lead to the development of mechanistic hypotheses applicable to the development of novel therapies in personalized medicine.

Bioinformatic pipelines combining denoising and clustering tools allow for more comprehensive prokaryotic and eukaryotic metabarcoding

Environmental DNA metabarcoding is a powerful tool for studying biodiversity. However, bioinformatic approaches need to adjust to the diversity of taxonomic compartments targeted as well as to each barcode gene specificities. We built and tested a pipeline based on read correction with DADA2 allowing analysing metabarcoding data from prokaryotic (16S) and eukaryotic (18S, COI) life compartments. We implemented the option to cluster amplicon sequence variants (ASVs) into operational taxonomic units (OTUs) with swarm, a network-based clustering algorithm, and the option to curate ASVs/OTUs using LULU. Finally, taxonomic assignment was implemented via the Ribosomal Database Project Bayesian classifier (RDP) and BLAST. We validated this pipeline with ribosomal and mitochondrial markers using metazoan mock communities and 42 deep-sea sediment samples. The results show that ASVs and OTUs describe different levels of biotic diversity, the choice of which depends on the research questions. They underline the advantages and complementarity of clustering and LULU-curation for producing metazoan biodiversity inventories at a level approaching the one obtained using morphological criteria. While clustering removes intraspecific variation, LULU effectively removes spurious clusters, originating from errors or intragenomic variability. Swarm clustering affected alpha and beta diversity differently depending on genetic marker. Specifically, d-values > 1 appeared to be less appropriate with 18S for metazoans. Similarly, increasing LULU's minimum ratio level proved essential to avoid losing species in sample-poor data sets. Comparing BLAST and RDP underlined that accurate assignments of deep-sea species can be obtained with RDP, but highlighted the need for a concerted effort to build comprehensive, ecosystem-specific databases.

Shifts in the Skin Microbiota after UVB Treatment in Adult Atopic Dermatitis

BACKGROUND

The pathophysiology in atopic dermatitis (AD) is not fully understood, but immune dysfunction, skin barrier defects, and alterations of the skin microbiota are thought to play important roles. AD skin is frequently colonized with Staphylococcus aureus (S. aureus) and microbial diversity on lesional skin (LS) is reduced compared to on healthy skin. Treatment with narrow-band ultraviolet B (nb-UVB) leads to clinical improvement of the eczema and reduced abundance of S. aureus. However, in-depth knowledge of the temporal dynamics of the skin microbiota in AD in response to nb-UVB treatment is lacking and could provide important clues to decipher whether the microbial changes are primary drivers of the disease, or secondary to the inflammatory process.

OBJECTIVES

To map the temporal shifts in the microbiota of the skin, nose, and throat in adult AD patients after nb-UVB treatment.

METHODS

Skin swabs were taken from lesional AD skin (n = 16) before and after 3 treatments of nb-UVB, and after 6-8 weeks of full-body treatment. We also obtained samples from non-lesional skin (NLS) and from the nose and throat. All samples were characterized by 16S rRNA gene sequencing.

RESULTS

We observed shifts towards higher diversity in the microbiota of lesional AD skin after 6-8 weeks of treatment, while the microbiota of NLS and of the nose/throat remained unchanged. After only 3 treatments with nb-UVB, there were no significant changes in the microbiota.

CONCLUSION

Nb-UVB induces changes in the skin microbiota towards higher diversity, but the microbiota of the nose and throat are not altered.

Metabacillus schmidteae sp. nov., Cultivated from Planarian Schmidtea mediterranea Microbiota

Taxonogenomics combines phenotypic assays and genomic analysis as a means of characterizing novel strains. We used this strategy to study Marseille-P9898T strain, an aerobic, motile, Gram-negative, spore-forming, and rod-shaped bacterium isolated from planarian Schmidtea mediterranea. Marseille-P9898T is catalase-positive and oxidase-negative. The major fatty acids detected are 12-methyl-tetradecanoic acid, 13-methyl-tetradecanoic acid, and hexadecanoic acid. Marseille-P9898T strain shared more than 98% sequence similarity with the Metabacillus niabensis strain 4T19T (98.99%), Metabacillus halosaccharovorans strain E33T (98.75%), Metabacillus malikii strain NCCP-662T (98.19%), and Metabacillus litoralis strain SW-211T (97.15%). Marseille-P9898 strain belongs to Metabacillus genus. Genomic analysis revealed the highest similarities with Ortho-ANI and dDDH, 85.76% with Metabacillus halosaccharovorans, and 34.20% with Bacillus acidicola, respectively. These results show that the Marseille-P9898T strain is a novel bacterial species from Metabacillus genus, for which we propose the name of Metabacillus schmidteae sp. nov. (Type strain Marseille-P9898T = CSUR P9898T = DSM 111480T).

Pedobacter ghigonii sp. nov., Isolated from the Microbiota of the Planarian Schmidtea mediterranea

The planarian S. mediterranea is a platyhelminth with worldwide distribution that can regenerate any part of its body after amputation and has the capacity to eliminate a large spectrum of human bacterial pathogens. Surprisingly, the microbiota of S. mediterranea remains poorly investigated. Using the culturomics strategy to study the bacterial component of planarians, we isolated a new bacterial strain, Marseille-Q2390, which we characterized with the taxono-genomic approach that associates phenotypic assays and genome sequencing and analysis. Strain Marseille-Q2390 exhibited a 16S rRNA sequence similarity of 99.36% with Pedobacter kyungheensis strain THG-T17T, the closest phylogenetic neighbor. It is a white-pigmented, Gram-negative, and rod-shaped bacterium. It grows in aerobic conditions and belongs to the family Sphingobacteriaceae. The genome of strain Marseille-Q2390 is 5,919,359 bp-long, with a G + C content of 40.3%. By comparing its genome with other closely related strains, the highest Orthologous Average Nucleotide Identity (Ortho-ANI) and digital DNA-DNA hybridization (dDDH) values were 85.71% and 30.50%, respectively, which were found with Pedobacter soli strain 15-51T. We conclude that strain Marseille-Q2390T is sufficiently different from other nearby species to be classified within a new species for which we propose the name Pedobacter ghigonii sp. nov.

Genomic analysis of dibenzofuran-degrading Pseudomonas veronii strain Pvy reveals its biodegradative versatility

Certain industrial chemicals accumulate in the environment due to their recalcitrant properties. Bioremediation uses the capability of some environmental bacteria to break down these chemicals and attenuate the pollution. One such bacterial strain, designated Pvy, was isolated from sediment samples from a lagoon in Romania located near an oil refinery due to its capacity to degrade dibenzofuran (DF). The genome sequence of the Pvy strain was obtained using an Oxford Nanopore MiniION platform. According to the consensus 16S rRNA gene sequence that was compiled from six 16S rRNA gene copies contained in the genome and orthologous average nucleotide identity (OrthoANI) calculation, the Pvy strain was identified as Pseudomonas veronii, which confirmed the identification obtained with the aid of MALDI-TOF mass spectrometry and MALDI BioTyper. The genome was analyzed with respect to enzymes responsible for the overall biodegradative versatility of the strain. The Pvy strain was able to derive carbon from naphthalene (NP) and several aromatic compounds of natural origin, including salicylic, protocatechuic, p-hydroxybenzoic, trans-cinnamic, vanillic, and indoleacetic acids or vanillin, and was shown to degrade but not utilize DF. In total seven loci were found in the Pvy genome, which enables the strain to participate in the degradation of these aromatic compounds. Our experimental data also indicate that the transcription of the NP-dioxygenase α-subunit gene (ndoB), carried by the plasmid of the Pvy strain, is inducible by DF. These features make the Pvy strain a potential candidate for various bioremediation applications.

The Subgingival Microbiome in Patients with Down Syndrome and Periodontitis

Objective: To describe the subgingival microbiome of individuals with Down syndrome (DS). Methods: We conducted a cross-sectional observational study that obtained bacterial DNA samples from 50 patients with DS, 25 with periodontitis (PDS) and 25 with a healthy periodontal condition (HDS). The samples were analyzed by sequencing the 16S rRNA gene V3–V4 hypervariable region using the MiSeq System. Taxonomic affiliations were assigned using the naïve Bayesian classifier integrated in QIIME2 plugins. We evaluated the difference in bacteria abundance between the sample groups using Wilcoxon and Kruskal–Wallis tests. We evaluated the alpha diversity of the identified species using the Observed, Chao1metric, ACE and Shannon indices and evaluated beta diversity with principal coordinate analysis (registration code: 2018/510). Results: Twenty-one genera and 39 bacterial species showed a significantly different abundance between the study groups. Among the genera, Porphyromonas, Treponema, Tannerella and Aggregatibacter were more abundant in the PDS group than in the HDS group, as were the less commonly studied Filifactor, Fretibacterium and Desulfobulbus genera. Among the species, Porphyromonas spp. and Tannerella spp. were the most abundant in the PDS group; the most abundant species in the HDS group were Pseudomonas spp., Granulicatella spp. and Gemella spp. Conclusion: Well-recognized periodontal pathogens and newly proposed pathogenic taxa were associated with periodontitis in patients with DS.

Development of a new TaqMan-based PCR assay for the specific detection and quantification of Simkania negevensis

Simkania negevensis is an emerging Chlamydia-like bacterium related to human respiratory diseases. An early and accurate detection of this pathogen could be useful to monitor the potential infectious risks and to set suitable outbreak control measures. In Tunisia, distribution and abundance of S. negevensis remain until now largely unknown. In the present work, a qPCR assay, targeting the 16S rRNA gene, for fast detection and quantification of S. negevensis was developed and validated. A high specificity for S. negevensis detection displaying no cross-reaction with the closely related Chlamydia spp. or the other tested microorganisms was noticed. qPCR assay performance was considered very satisfying with detection limits of 5 DNA copies per reaction. qPCR assay validation was performed by screening 37 clinical specimens and 35 water samples. S. negevensis wasn't detected in respiratory samples, but it was found in four cases of water samples. We suggest that the qPCR assay developed in this study could be considered sufficiently characterized to initiate the quantification of S. negevensis in environmental samples.

Species-specific genomic sequences for classification of bacteria

Modern bacterial classification relies on genomic relatedness. Genetic variation in bacterial populations present a big challenge for taxonomic classification and recently several bacterial species have been reclassified based on the intra-species genome comparison. These were facilitated by next generation sequencing technologies and advances in genome comparison approaches which led to the rearrangement of diverse bacterial species and revolution in the microbial classification system. One of the outcome of these studies is the development of suitable DNA barcodes as reliable and cost-effective method for identifying various bacterial genera. Towards refining this further, we have applied a genome comparison approach in 1104 bacterial genome assemblies (excluding plasmids) to identify unique genomic segments among intra-species genome assemblies. Using extensive bioinformatics analysis, we have identified species-specific genomic regions and designed unique primers for 100 different species (belonging to 62 genera) which includes 62 pathogenic and 13 opportunistic pathogenic bacterial species and built a database (http://slsdb.manipal.edu/bact/). These species-specific genomic regions will have a major impact on in silico and molecular methods aimed at bacterial classification and identification. These may also serve as better DNA barcodes than the markers currently used for delineation of bacteria and may also find application in various translational research programs.

Dispersal mitigates bacterial dominance over microalgal competitor in metacommunities

Ecological theory suggests that a combination of local and regional factors regulate biodiversity and community functioning in metacommunities. The relative importance of different factors structuring communities likely changes over successional time, but to date this concept is scarcely documented. In addition, the few studies describing successional dynamics in metacommunity regulation have only focused on a single group of organisms. Here, we report results of an experimental study testing the effect size of initial local community composition and dispersal between local patches on community dynamics of benthic microalgae and their associated bacteria over community succession. Our results show that over time dispersal outweighed initial effects of community composition on microalgal evenness and biomass, microalgal β-diversity, and the ratio of bacteria to microalgae. At the end of the experiment (ca. 20 microalgae generations), dispersal significantly decreased microalgal evenness and β-diversity by promoting one regionally superior competitor. Dispersal also decreased the ratio of bacteria to microalgae, while it significantly increased microalgal biomass. These results suggest that the dispersal-mediated establishment of a dominant and superior microalgae species prevented bacteria from gaining competitive advantage over the autotrophs in these metacommunities, ultimately maintaining the provision of autotrophic biomass. Our study emphasizes the importance of time for dispersal to be a relevant community-structuring mechanism. Moreover, we highlight the need for considering multiple competitors in complex metacommunity systems to properly pinpoint the consequences of local change in dominance through dispersal for metacommunity function.

Phylogenetic comparison between Type IX Secretion System (T9SS) protein components suggests evidence of horizontal gene transfer

Porphyromonas gingivalis is one of the major bacteria that causes periodontitis. Chronic periodontitis is a severe form of periodontal disease that ultimately leads to tooth loss. Virulence factors that contribute to periodontitis are secreted by Type IX Secretion System (T9SS). There are aspects of T9SS protein components that have yet to be characterised. Thus, the aim of this study is to investigate the phylogenetic relationship between members of 20 T9SS component protein families. The Bayesian Inference (BI) trees for 19 T9SS protein components exhibit monophyletic clades for all major classes under Bacteroidetes with strong support for the monophyletic clades or its subclades that is consistent with phylogeny exhibited by the constructed BI tree of 16S rRNA. The BI tree of PorR is different from the 19 BI trees of T9SS protein components as it does not exhibit monophyletic clades for all major classes under Bacteroidetes. There is strong support for the phylogeny exhibited by the BI tree of PorR which deviates from the phylogeny based on 16S rRNA. Hence, it is possible that the porR gene is subjected to horizontal transfer as it is known that virulence factor genes could be horizontally transferred. Seven genes (porR included) that are involved in the biosynthesis of A-LPS are found to be flanked by insertion sequences (IS5 family transposons). Therefore, the intervening DNA segment that contains the porR gene might be transposed and subjected to conjugative transfer. Thus, the seven genes can be co-transferred via horizontal gene transfer. The BI tree of UgdA does not exhibit monophyletic clades for all major classes under Bacteroidetes which is similar to the BI tree of PorR (both are a part of the seven genes). Both BI trees also exhibit similar topology as the four identified clusters with strong support and have similar relative positions to each other in both BI trees. This reinforces the possibility that porR and the other six genes might be horizontally transferred. Other than the BI tree of PorR, the 19 other BI trees of T9SS protein components also exhibit evidence of horizontal gene transfer. However, their genes might undergo horizontal gene transfer less frequently compared to porR because the intervening DNA segment that contains porR is easily exchanged between bacteria under Bacteroidetes due to the presence of insertion sequences (IS5 family transposons) that flank it. In conclusion, this study can provide a better understanding about the phylogeny of T9SS protein components.

Endosymbiotic adaptations in three new bacterial species associated with Dictyostelium discoideum: Paraburkholderia agricolaris sp. nov., Paraburkholderia hayleyella sp. nov., and Paraburkholderia bonniea sp. nov

Here we give names to three new species of Paraburkholderia that can remain in symbiosis indefinitely in the spores of a soil dwelling eukaryote, Dictyostelium discoideum. The new species P. agricolaris sp. nov., P. hayleyella sp. nov., and P. bonniea sp. nov. are widespread across the eastern USA and were isolated as internal symbionts of wild-collected D. discoideum. We describe these sp. nov. using several approaches. Evidence that they are each a distinct new species comes from their phylogenetic position, average nucleotide identity, genome-genome distance, carbon usage, reduced length, cooler optimal growth temperature, metabolic tests, and their previously described ability to invade D. discoideum amoebae and form a symbiotic relationship. All three of these new species facilitate the prolonged carriage of food bacteria by D. discoideum, though they themselves are not food. Further studies of the interactions of these three new species with D. discoideum should be fruitful for understanding the ecology and evolution of symbioses.

Isolation, characterization and bioaugmentation of an acidotolerant 1,2-dichloroethane respiring Desulfitobacterium species from a low pH aquifer

A Desulfitobacterium sp. strain AusDCA of the Peptococcaceae family capable of respiring 1,2-dichloroethane (1,2-DCA) to ethene anaerobically with ethanol or hydrogen as electron donor at pH 5.0 with optimal range between pH 6.5-7.5 was isolated from an acidic aquifer near Sydney, Australia. Strain AusDCA is distant (94% nucleotide identity) from its nearest phylogenetic neighbor, D. metallireducens, and could represent a new species. Reference gene-based quantification of growth indicated a doubling time of 2 days in cultures buffered at pH 7.2, and a yield of 7.66 (± 4.0) × 106 cells µmol-1 of 1,2-DCA. A putative 1,2-DCA reductive dehalogenase was translated from a dcaAB locus and had high amino acid identity (97.3% for DcaA and 100% for DcaB) to RdhA1B1 of the 1,2-DCA respiring Dehalobacter strain WL. Proteomic analysis confirmed DcaA expression in the pure culture. Dehalogenation of 1,2-DCA (1.6 mM) was observed in batch cultures established from groundwater at pH 5.5 collected 38 days after in situ bioaugmentation but not in cultures established with groundwater collected at the same time from wells not receiving bioaugmentation. Overall, strain AusDCA can tolerate lower pH than previously characterized organohalide respiring bacteria and remained viable in groundwater at pH 5.5.

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.