Gene content dissimilarity for subclassification of highly similar microbial strains

Diverse microbiome functions, limited temporal variation and substantial genomic conservation within sedimentary and granite rock deep underground research laboratories

BACKGROUND

Underground research laboratories (URLs) provide a window on the deep biosphere and enable investigation of potential microbial impacts on nuclear waste, CO2 and H2 stored in the subsurface. We carried out the first multi-year study of groundwater microbiomes sampled from defined intervals between 140 and 400 m below the surface of the Horonobe and Mizunami URLs, Japan.

RESULTS

We reconstructed draft genomes for > 90% of all organisms detected over a four year period. The Horonobe and Mizunami microbiomes are dissimilar, likely because the Mizunami URL is hosted in granitic rock and the Horonobe URL in sedimentary rock. Despite this, hydrogen metabolism, rubisco-based CO2 fixation, reduction of nitrogen compounds and sulfate reduction are well represented functions in microbiomes from both URLs, although methane metabolism is more prevalent at the organic- and CO2-rich Horonobe URL. High fluid flow zones and proximity to subsurface tunnels select for candidate phyla radiation bacteria in the Mizunami URL. We detected near-identical genotypes for approximately one third of all genomically defined organisms at multiple depths within the Horonobe URL. This cannot be explained by inactivity, as in situ growth was detected for some bacteria, albeit at slow rates. Given the current low hydraulic conductivity and groundwater compositional heterogeneity, ongoing inter-site strain dispersal seems unlikely. Alternatively, the Horonobe URL microbiome homogeneity may be explained by higher groundwater mobility during the last glacial period. Genotypically-defined species closely related to those detected in the URLs were identified in three other subsurface environments in the USA. Thus, dispersal rates between widely separated underground sites may be fast enough relative to mutation rates to have precluded substantial divergence in species composition. Species overlaps between subsurface locations on different continents constrain expectations regarding the scale of global subsurface biodiversity.

CONCLUSIONS

Our analyses reveal microbiome stability in the sedimentary rocks and surprising microbial community compositional and genotypic overlap over sites separated by hundreds of meters of rock, potentially explained by dispersal via slow groundwater flow or during a prior hydrological regime. Overall, microbiome and geochemical stability over the study period has important implications for underground storage applications.

Discovery of novel treponemes associated with pododermatitis in elk (Cervus canadensis)

Pododermatitis, also known as treponeme-associated hoof disease (TAHD), presents a significant challenge to elk (Cervus canadensis) populations in the northwestern USA, with Treponema spp. consistently implicated in the lesion development. However, identifying species-specific Treponema strains from these lesions is hindered by its culture recalcitrance and limited genomic information. This study utilized shotgun sequencing, in silico genome reconstruction, and comparative genomics as a culture-independent approach to identify metagenome-assembled Treponema genomes (MATGs) from skin scraping samples collected from captive elk experimentally challenged with TAHD. The genomic analysis revealed 10 new MATGs, with 6 representing novel genomospecies associated with pododermatitis in elk and 4 corresponding to previously identified species-Treponema pedis and Treponema phagedenis. Importantly, genomic signatures of novel genomospecies identified in this study were consistently detected in biopsy samples of free-ranging elk diagnosed with TAHD, indicating a potential etiologic association. Comparative metabolic profiling of the MATGs against other Treponema genomes showed a distinct metabolic profile, suggesting potential host adaptation or geographic uniqueness of these newly identified genomospecies. The discovery of novel Treponema genomospecies enhances our understanding of the pathogenesis of pododermatitis and lays the foundation for the development of improved molecular surveillance tools to monitor and manage the disease in free-ranging elk.IMPORTANCETreponema spp. play an important role in the development of pododermatitis in free-ranging elk; however, the species-specific detection of Treponema from pododermatitis lesions is challenging due to culture recalcitrance and limited genomic information. The study utilized shotgun sequencing and in silico genome reconstruction to identify novel Treponema genomospecies from elk with pododermatitis. The discovery of the novel Treponema species opens new avenues to develop molecular diagnostic and epidemiologic tools for the surveillance of pododermatitis in elk. These findings significantly enhance our understanding of the genomic landscape of the Treponemataceae consortium while offering valuable insights into the etiology and pathogenesis of emerging pododermatitis in elk populations.

Mangrovimonas cancribranchiae sp. nov., a novel bacterial species associated with the gills of the fiddler crab Cranuca inversa (Brachyura, Ocypodidae) from Red Sea mangroves

Two bacteria, UG2_1T and UG2_2, were isolated from the gill tissues of the mangrove fiddler crab Cranuca inversa collected on the east coast of the Red Sea (Thuwal, Saudi Arabia). The cells are Gram-negative, rod-shaped, orange-pigmented, motile by gliding with no flagella, strictly aerobic, and grow at 20-37 °C (optimum, 28-35 °C), at pH 5.0-9.0 (optimum, pH 6.0-7.0), and with 1-11 % (w/v) NaCl (optimum, 2-4 %). They were positive for oxidase and catalase activity. Phylogenetic analysis based on 16S rRNA gene sequences indicated that isolates UG2_1T and UG2_2 belong to the genus Mangrovimonas, showing the highest similarity to Mangrovimonas spongiae HN-E26T (99.4 %). Phylogenomic analysis based on the whole genomes, independently using 49 and 120 concatenated genes, showed that strains UG2_1T and UG2_2 formed a monophyletic lineage in a different cluster from other type strain species within the genus Mangrovimonas. The genome sizes were 3.08 and 3.07 Mbp for UG2_1T and UG2_2, respectively, with a G+C content of 33.8 mol% for both strains. Values of average nucleotide identity and digital DNA-DNA hybridization between the strains and closely related species were 91.0 and 43.5 %, respectively. Chemotaxonomic analysis indicated that both strains had iso-C15 : 0 and iso-C15 : 1 G as dominant fatty acids, and the primary respiratory quinone was identified as MK-6. The major polar lipids comprised phosphatidylethanolamine, one unidentified glycolipid, one unidentified phospholipid, two unidentified aminolipids, and four unidentified lipids. Based on phylogenetic, phylogenomic, genome relatedness, phenotypic, and chemotaxonomical data, the two isolates represent a novel species within the genus Mangrovimonas, with the proposed name Mangrovimonas cancribranchiae sp. nov., and the type strain UG2_1T (=KCTC 102158T=DSM 117025T).

Identification of an Aerococcus urinaeequi isolate by Whole Genome Sequencing and Average Nucleotide Identity analysis

OBJECTIVES

Identification and classification of microorganisms is one of the most important but difficult and challenging issues in microbiology. Whole genome sequencing (WGS), which can give a thorough understanding for the genome of bacteria strain, has been universally used for studying bacterial classification, evolution, and drug-related resistant genes. We in this study aimed to identify a gram-positive, microaerophilic, catalase-negative cocci strain named AV208, which has shown resistance to vancomycin, by whole genome's average nucleotide identity (ANI) and high-throughput sequencing technology.

METHODS

The AV208 strain was identified by following commercially available identification systems, including API 20 Strep system and Vitek 2 Compact gram-positive identification system for biochemical phenotypic test. Matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF-MS) and 16S rRNA gene sequencing were used for confirmation identification. The whole genome of AV208 was sequenced by using high throughput sequencing technology and ANI between AV208, and its phylogenetic neighbors were analyzed by the Orthologous Average Nucleotide Identity Tool (OAT) software. Polymerase chain reaction (PCR) and DNA sequencing were used to investigate the potential molecular mechanism for vancomycin resistance.

RESULTS

The AV208 strain was isolated from an ascites sample from a patient with chronic kidney disease who showed extensive resistance to the drugs detected, such as vancomycin with MIC > 256 μg/ml. With combination of biochemical phenotypic test, MALDI-TOF-MS and 16S rRNA gene sequencing, the AV208 strain was tentatively identified as an Aercoccus viridans. By using complete genome sequence, we found a 96.24% ANI between strain AV208 and Aerococcus urinaeequi CCUG 28094^T, which was higher than that with A. viridans CCUG4311^T (94.9%). The consistency of 16S rRNA sequence of strain AV208 was 100% with A. urinaeequi CCUG 28094^T and 99.9% with A. viridans CCUG4311^T, with only one base difference between them. PCR and sequencing for van genes revealed that AV208 was positive for the vanA gene. A Tn1546 transposon-like structure with vanA gene was found in the genome, which was predicted locating in plasmid, causing vancomycin resistance phenotypes.

CONCLUSIONS

Average nucleotide identity analysis based on whole genome sequence is an accurate and effective method for identification of bacteria, especially for strains that are not discernible by existing methods such as Aerococcus.

Comparative evaluation of Vibrio delineation methodologies in post-genomic era

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

On a Non-Discrete Concept of Prokaryotic Species

The taxonomic concept of species has received continuous attention. A microbial species as a discrete box contains a limited number of highly similar microorganisms assigned to that taxon, following a polyphasic approach. In the 21st Century, with the advancements of sequencing technologies and genomics, the existence of a huge prokaryotic diversity has become well known. At present, the prokaryotic species might no longer have to be understood as discrete values (such as 1 or 2, by homology to Natural numbers); rather, it is expected that some microorganisms could be potentially distributed (according to their genome features and phenotypes) in between others (such as decimal numbers between 1 and 2; real numbers). We propose a continuous species concept for microorganisms, which adapts to the current knowledge on the huge diversity, variability and heterogeneity existing among bacteria and archaea. Likely, this concept could be extended to eukaryotic microorganisms. The continuous species concept considers a species to be delimited by the distance between a range of variable features following a Gaussian-type distribution around a reference organism (i.e., its type strain). Some potential pros and cons of a continuous concept are commented on, offering novel perspectives on our understanding of the highly diversified prokaryotic world, thus promoting discussion and further investigation in the field.

Phylogenomics Reveals Clear Cases of Misclassification and Genus-Wide Phylogenetic Markers for Acinetobacter

The Gram-negative Acinetobacter genus has several species of clear medical relevance. Many fully sequenced genomes belonging to the genus have been published in recent years; however, there has not been a recent attempt to infer the evolutionary history of Acinetobacter with that vast amount of information. Here, through a phylogenomic approach, we established the most up-to-date view of the evolutionary relationships within this genus and highlighted several cases of poor classification, especially for the very closely related species within the Acinetobacter calcoaceticus-Acinetobacter baumannii complex (Acb complex). Furthermore, we determined appropriate phylogenetic markers for this genus and showed that concatenation of the top 13 gives a very decent reflection of the evolutionary relationships for the genus Acinetobacter. The intersection between our top markers and previously defined universal markers is very small. In general, our study shows that, although there seems to be hardly any universal markers, bespoke phylogenomic approaches can be used to infer the phylogeny of different bacterial genera. We expect that ad hoc phylogenomic approaches will be the standard in the years to come and will provide enough information to resolve intricate evolutionary relationships like those observed in the Acb complex.

Genome-Based Comparison of Clostridioides difficile: Average Amino Acid Identity Analysis of Core Genomes

Infections due to Clostridioides difficile (previously known as Clostridium difficile) are a major problem in hospitals, where cases can be caused by community-acquired strains as well as by nosocomial spread. Whole genome sequences from clinical samples contain a lot of information but that needs to be analyzed and compared in such a way that the outcome is useful for clinicians or epidemiologists. Here, we compare 663 public available complete genome sequences of C. difficile using average amino acid identity (AAI) scores. This analysis revealed that most of these genomes (640, 96.5%) clearly belong to the same species, while the remaining 23 genomes produce four distinct clusters within the Clostridioides genus. The main C. difficile cluster can be further divided into sub-clusters, depending on the chosen cutoff. We demonstrate that MLST, either based on partial or full gene-length, results in biased estimates of genetic differences and does not capture the true degree of similarity or differences of complete genomes. Presence of genes coding for C. difficile toxins A and B (ToxA/B), as well as the binary C. difficile toxin (CDT), was deduced from their unique PfamA domain architectures. Out of the 663 C. difficile genomes, 535 (80.7%) contained at least one copy of ToxA or ToxB, while these genes were missing from 128 genomes. Although some clusters were enriched for toxin presence, these genes are variably present in a given genetic background. The CDT genes were found in 191 genomes, which were restricted to a few clusters only, and only one cluster lacked the toxin A/B genes consistently. A total of 310 genomes contained ToxA/B without CDT (47%). Further, published metagenomic data from stools were used to assess the presence of C. difficile sequences in blinded cases of C. difficile infection (CDI) and controls, to test if metagenomic analysis is sensitive enough to detect the pathogen, and to establish strain relationships between cases from the same hospital. We conclude that metagenomics can contribute to the identification of CDI and can assist in characterization of the most probable causative strain in CDI patients.

GET_PHYLOMARKERS, a Software Package to Select Optimal Orthologous Clusters for Phylogenomics and Inferring Pan-Genome Phylogenies, Used for a Critical Geno-Taxonomic Revision of the Genus Stenotrophomonas

The massive accumulation of genome-sequences in public databases promoted the proliferation of genome-level phylogenetic analyses in many areas of biological research. However, due to diverse evolutionary and genetic processes, many loci have undesirable properties for phylogenetic reconstruction. These, if undetected, can result in erroneous or biased estimates, particularly when estimating species trees from concatenated datasets. To deal with these problems, we developed GET_PHYLOMARKERS, a pipeline designed to identify high-quality markers to estimate robust genome phylogenies from the orthologous clusters, or the pan-genome matrix (PGM), computed by GET_HOMOLOGUES. In the first context, a set of sequential filters are applied to exclude recombinant alignments and those producing anomalous or poorly resolved trees. Multiple sequence alignments and maximum likelihood (ML) phylogenies are computed in parallel on multi-core computers. A ML species tree is estimated from the concatenated set of top-ranking alignments at the DNA or protein levels, using either FastTree or IQ-TREE (IQT). The latter is used by default due to its superior performance revealed in an extensive benchmark analysis. In addition, parsimony and ML phylogenies can be estimated from the PGM. We demonstrate the practical utility of the software by analyzing 170 Stenotrophomonas genome sequences available in RefSeq and 10 new complete genomes of Mexican environmental S. maltophilia complex (Smc) isolates reported herein. A combination of core-genome and PGM analyses was used to revise the molecular systematics of the genus. An unsupervised learning approach that uses a goodness of clustering statistic identified 20 groups within the Smc at a core-genome average nucleotide identity (cgANIb) of 95.9% that are perfectly consistent with strongly supported clades on the core- and pan-genome trees. In addition, we identified 16 misclassified RefSeq genome sequences, 14 of them labeled as S. maltophilia, demonstrating the broad utility of the software for phylogenomics and geno-taxonomic studies. The code, a detailed manual and tutorials are freely available for Linux/UNIX servers under the GNU GPLv3 license at https://github.com/vinuesa/get_phylomarkers. A docker image bundling GET_PHYLOMARKERS with GET_HOMOLOGUES is available at https://hub.docker.com/r/csicunam/get_homologues/, which can be easily run on any platform.

Phenetic Comparison of Prokaryotic Genomes Using k-mers

Bacterial genomics studies are getting more extensive and complex, requiring new ways to envision analyses. Using the Ray Surveyor software, we demonstrate that comparison of genomes based on their k-mer content allows reconstruction of phenetic trees without the need of prior data curation, such as core genome alignment of a species. We validated the methodology using simulated genomes and previously published phylogenomic studies of Streptococcus pneumoniae and Pseudomonas aeruginosa. We also investigated the relationship of specific genetic determinants with bacterial population structures. By comparing clusters from the complete genomic content of a genome population with clusters from specific functional categories of genes, we can determine how the population structures are correlated. Indeed, the strain clustering based on a subset of k-mers allows determination of its similarity with the whole genome clusters. We also applied this methodology on 42 species of bacteria to determine the correlational significance of five important bacterial genomic characteristics. For example, intrinsic resistance is more important in P. aeruginosa than in S. pneumoniae, and the former has increased correlation of its population structure with antibiotic resistance genes. The global view of the pangenome of bacteria also demonstrated the taxa-dependent interaction of population structure with antibiotic resistance, bacteriophage, plasmid, and mobile element k-mer data sets.

Antibacterial Activities of Bacteria Isolated from the Marine Sponges Isodictya compressa and Higginsia bidentifera Collected from Algoa Bay, South Africa

Due to the rise in multi-drug resistant pathogens and other diseases, there is renewed interest in marine sponge endosymbionts as a rich source of natural products (NPs). The South African marine environment is rich in marine biota that remains largely unexplored and may represent an important source for the discovery of novel NPs. We first investigated the bacterial diversity associated with five South African marine sponges, whose microbial populations had not previously been investigated, and select the two sponges (Isodictya compressa and Higginsia bidentifera) with highest species richness to culture bacteria. By employing 33 different growth conditions 415 sponge-associated bacterial isolates were cultured and screened for antibacterial activity. Thirty-five isolates showed antibacterial activity, twelve of which exhibited activity against the multi-drug resistant Escherichia coli 1699, implying that some of the bioactive compounds could be novel. Genome sequencing of two of these isolates confirmed that they harbour uncharacterized biosynthetic pathways that may encode novel chemical structures.