acdc - Automated Contamination Detection and Confidence estimation for single-cell genome data

Gene inversion led to the emergence of brackish archaeal heterotrophs in the aftermath of the Cryogenian Snowball Earth

Land-ocean interactions greatly impact the evolution of coastal life on earth. However, the ancient geological forces and genetic mechanisms that shaped evolutionary adaptations and allowed microorganisms to inhabit coastal brackish waters remain largely unexplored. In this study, we infer the evolutionary trajectory of the ubiquitous heterotrophic archaea Poseidoniales (Marine Group II archaea) presently occurring across global aquatic habitats. Our results show that their brackish subgroups had a single origination, dated to over 600 million years ago, through the inversion of the magnesium transport gene corA that conferred osmotic-stress tolerance. The subsequent loss and gain of corA were followed by genome-wide adjustment, characterized by a general two-step mode of selection in microbial speciation. The coastal family of Poseidoniales showed a rapid increase in the evolutionary rate during and in the aftermath of the Cryogenian Snowball Earth (∼700 million years ago), possibly in response to the enhanced phosphorus supply and the rise of algae. Our study highlights the close interplay between genetic changes and ecosystem evolution that boosted microbial diversification in the Neoproterozoic continental margins, where the Cambrian explosion of animals soon followed.

First Report of IMI-2-Producing Enterobacter bugandensis and CTX-M-55-Producing Escherichia coli isolated from Healthy Volunteers in Tunisia

The aim of this study was to characterize the prevalence of fecal carriage of extended-spectrum beta-lactamases and carbapenemase-producing Gram-negative bacteria among healthy humans in Tunisia. Fifty-one rectal swabs of healthy volunteers were plated on MacConkey agar plates supplemented with cefotaxime or imipenem. The occurrences of resistance genes, integrons, and phylogroup typing were investigated using PCR and sequencing. The genetic relatedness of isolates was determined by pulsed-field-gel-electrophoresis (PFGE) and multilocus-sequence-typing (MLST). Whole-genome-sequencing (WGS) was performed for the carbapenem-resistant isolate. Sixteen ESBL-producing Escherichia coli isolates and one carbapenem-resistant Enterobacter bugandensis were detected out of the fifty-one fecal samples. The ESBL-producing E. coli strains contained genes encoding CTX-M-15 (n = 9), CTX-M-1 (n = 3), CTX-M-27 (n = 3), and CTX-M-55 (n = 1). Three CTX-M-1-producers were of lineages ST131, ST7366, and ST1158; two CTX-M-15-producers belonged to lineage ST925 and ST5100; one CTX-M-27-producer belonged to ST2887, and one CTX-M-15-producer belonged to ST744. Six isolates contained class 1 integrons with the following four gene cassette arrangements: dfrA5 (two isolates), dfrA12-orf-aadA2 (two isolates), dfrA17-aadA5 (one isolate), and aadA1 (one isolate). E. bugandensis belonged to ST1095, produced IMI-2 carbapenemase, and contained qnrE1 and fosA genes. A genome-sequence analysis of the E. bugandensis strain revealed new mutations in the bla_ACT and qnr genes. Our results reveal an alarming rate of ESBL-E. coli in healthy humans in Tunisia and the first description of IMI-2 in E. bugandensis.

On clustering for cell-phenotyping in multiplex immunohistochemistry (mIHC) and multiplexed ion beam imaging (MIBI) data

OBJECTIVE

Multiplex immunohistochemistry (mIHC) and multiplexed ion beam imaging (MIBI) images are usually phenotyped using a manual thresholding process. The thresholding is prone to biases, especially when examining multiple images with high cellularity.

RESULTS

Unsupervised cell-phenotyping methods including PhenoGraph, flowMeans, and SamSPECTRAL, primarily used in flow cytometry data, often perform poorly or need elaborate tuning to perform well in the context of mIHC and MIBI data. We show that, instead, semi-supervised cell clustering using Random Forests, linear and quadratic discriminant analysis are superior. We test the performance of the methods on two mIHC datasets from the University of Colorado School of Medicine and a publicly available MIBI dataset. Each dataset contains a bunch of highly complex images.

Music of metagenomics-a review of its applications, analysis pipeline, and associated tools

This humble effort highlights the intricate details of metagenomics in a simple, poetic, and rhythmic way. The paper enforces the significance of the research area, provides details about major analytical methods, examines the taxonomy and assembly of genomes, emphasizes some tools, and concludes by celebrating the richness of the ecosystem populated by the "metagenome."

Single-Cell Genomics of Novel Actinobacteria With the Wood-Ljungdahl Pathway Discovered in a Serpentinizing System

Serpentinite-hosted systems represent modern-day analogs of early Earth environments. In these systems, water-rock interactions generate highly alkaline and reducing fluids that can contain hydrogen, methane, and low-molecular-weight hydrocarbons-potent reductants capable of fueling microbial metabolism. In this study, we investigated the microbiota of Hakuba Happo hot springs (∼50°C; pH∼10.5-11), located in Nagano (Japan), which are impacted by the serpentinization process. Analysis of the 16S rRNA gene amplicon sequences revealed that the bacterial community comprises Nitrospirae (47%), "Parcubacteria" (19%), Deinococcus-Thermus (16%), and Actinobacteria (9%), among others. Notably, only 57 amplicon sequence variants (ASV) were detected, and fifteen of these accounted for 90% of the amplicons. Among the abundant ASVs, an early-branching, uncultivated actinobacterial clade identified as RBG-16-55-12 in the SILVA database was detected. Ten single-cell genomes (average pairwise nucleotide identity: 0.98-1.00; estimated completeness: 33-93%; estimated genome size: ∼2.3 Mb) that affiliated with this clade were obtained. Taxonomic classification using single copy genes indicates that the genomes belong to the actinobacterial class-level clade UBA1414 in the Genome Taxonomy Database. Based on metabolic pathway predictions, these actinobacteria are anaerobes, capable of glycolysis, dissimilatory nitrate reduction and CO₂ fixation via the Wood-Ljungdahl (WL) pathway. Several other genomes within UBA1414 and two related class-level clades also encode the WL pathway, which has not yet been reported for the Actinobacteria phylum. For the Hakuba actinobacterium, the energy metabolism related to the WL pathway is likely supported by a combination of the Rnf complex, group 3b and 3d [NiFe]-hydrogenases, [FeFe]-hydrogenases, and V-type (H⁺/Na⁺ pump) ATPase. The genomes also harbor a form IV ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) complex, also known as a RubisCO-like protein, and contain signatures of interactions with viruses, including clustered regularly interspaced short palindromic repeat (CRISPR) regions and several phage integrases. This is the first report and detailed genome analysis of a bacterium within the Actinobacteria phylum capable of utilizing the WL pathway. The Hakuba actinobacterium is a member of the clade UBA1414/RBG-16-55-12, formerly within the group "OPB41." We propose to name this bacterium 'Candidatus Hakubanella thermoalkaliphilus.'

Draft Genome Sequence of Haloferax volcanii SS0101, Isolated from Salt Farms in Samut Sakhon, Thailand

Haloferax volcanii SS0101 is a halophilic archaeon isolated from salt farms in Thailand. The genome sequence of H. volcanii SS0101 contains a gene encoding capreomycidine synthase, a key enzyme for capreomycidine biosynthesis. This 3.8-Mb draft genome sequence of H. volcanii SS0101 will provide the tools for investigating genes involved in capeomycidine production in haloarchaea.

Haloferax volcanii SS0101 is a halophilic archaeon isolated from salt farms in Thailand. The genome sequence of H. volcanii SS0101 contains a gene encoding capreomycidine synthase, a key enzyme for capreomycidine biosynthesis. This 3.8-Mb draft genome sequence of H. volcanii SS0101 will provide the tools for investigating genes involved in capeomycidine production in haloarchaea.

Primer-free FISH probes from metagenomics/metatranscriptomics data permit the study of uncharacterised taxa in complex microbial communities

Methods for the study of member species in complex microbial communities remain a high priority, particularly for rare and/or novel member species that might play an important ecological role. Specifically, methods that link genomic information of member species with its spatial structure are lacking. This study adopts an integrative workflow that permits the characterisation of previously unclassified bacterial taxa from microbiomes through: (1) imaging of the spatial structure; (2) taxonomic classification and (3) genome recovery. Our study attempts to bridge the gaps between metagenomics/metatranscriptomics and high-resolution biomass imaging methods by developing new fluorescence in situ hybridisation (FISH) probes—termed as R-Probes—from shotgun reads that harbour hypervariable regions of the 16S rRNA gene. The sample-centric design of R-Probes means that probes can directly hybridise to OTUs as detected in shotgun sequencing surveys. The primer-free probe design captures larger microbial diversity as compared to canonical probes. R-Probes were designed from deep-sequenced RNA-Seq datasets for both FISH imaging and FISH–Fluorescence activated cell sorting (FISH–FACS). FISH–FACS was used for target enrichment of previously unclassified bacterial taxa prior to downstream multiple displacement amplification (MDA), genomic sequencing and genome recovery. After validation of the workflow on an axenic isolate of Thauera species, the techniques were applied to investigate two previously uncharacterised taxa from a tropical full-scale activated sludge community. In some instances, probe design on the hypervariable region allowed differentiation to the species level. Collectively, the workflow can be readily applied to microbiomes for which shotgun nucleic acid survey data is available.

Antarctic Streptomyces fildesensis So13.3 strain as a promising source for antimicrobials discovery

Antarctic have been suggested as an attractive source for antibiotics discovery and members of Streptomyces genus have historically been studied as natural producers of antimicrobial metabolites. Nonetheless, our knowledge on antibiotic-producing Streptomyces from Antarctic is very limited. In this study, the antimicrobial activity of organic extracts from Antarctic Streptomyces strains was evaluated by disk diffusion assays and minimum inhibitory concentration. The strain Streptomyces sp. So13.3 showed the greatest antibiotic activity (MIC = 15.6 μg/mL) against Gram-positive bacteria and growth reduction of Gram‒negative pathogens. The bioactive fraction in the crude extract was revealed by TLC‒bioautography at R_f = 0.78 with molecular weight between 148 and 624 m/z detected by LC-ESI-MS/MS. The strain So13.3 was taxonomically affiliated as Streptomyces fildesensis. Whole genome sequencing and analysis suggested a 9.47 Mb genome size with 42 predicted biosynthetic gene clusters (BGCs) and 56 putative clusters representing a 22% of total genome content. Interestingly, a large number of them (11 of 42 BGCs and 40 of 56 putative BGCs), did not show similarities with other known BGCs. Our results highlight the potential of the Antarctic Streptomyces strains as a promising source of novel antimicrobials, particularly the strain Streptomyces fildesensis So13.3, which first draft genome is reported in this work.

Defending Our Public Biological Databases as a Global Critical Infrastructure

Progress in modern biology is being driven, in part, by the large amounts of freely available data in public resources such as the International Nucleotide Sequence Database Collaboration (INSDC), the world's primary database of biological sequence (and related) information. INSDC and similar databases have dramatically increased the pace of fundamental biological discovery and enabled a host of innovative therapeutic, diagnostic, and forensic applications. However, as high-value, openly shared resources with a high degree of assumed trust, these repositories share compelling similarities to the early days of the Internet. Consequently, as public biological databases continue to increase in size and importance, we expect that they will face the same threats as undefended cyberspace. There is a unique opportunity, before a significant breach and loss of trust occurs, to ensure they evolve with quality and security as a design philosophy rather than costly "retrofitted" mitigations. This Perspective surveys some potential quality assurance and security weaknesses in existing open genomic and proteomic repositories, describes methods to mitigate the likelihood of both intentional and unintentional errors, and offers recommendations for risk mitigation based on lessons learned from cybersecurity.

Niche Differentiation among Three Closely Related Competibacteraceae Clades at a Full-Scale Activated Sludge Wastewater Treatment Plant and Putative Linkages to Process Performance

Multiple clades within a microbial taxon often coexist within natural and engineered environments. Because closely related clades have similar metabolic potential, it is unclear how diversity is sustained and what factors drive niche differentiation. In this study, we retrieved three near-complete Competibacter lineage genomes from activated sludge metagenomes at a full-scale pure oxygen activated sludge wastewater treatment plant. The three genomes represent unique taxa within the Competibacteraceae A comparison of the genomes revealed differences in capacity for exopolysaccharide (EPS) biosynthesis, glucose fermentation to lactate, and motility. Using quantitative PCR (qPCR), we monitored these clades over a 2-year period. The clade possessing genes for motility and lacking genes for EPS biosynthesis (CPB_P15) was dominant during periods of suspended solids in the effluent. Further analysis of operational parameters indicate that the dominance of the CPB_P15 clade is associated with low-return activated sludge recycle rates and low wasting rates, conditions that maintain relatively high levels of biomass within the system.IMPORTANCE Members of the Competibacter lineage are relevant in biotechnology as glycogen-accumulating organisms (GAOs). Here, we document the presence of three Competibacteraceae clades in a full-scale activated sludge wastewater treatment plant and their linkage to specific operational conditions. We find evidence for niche differentiation among the three clades with temporal variability in clade dominance that correlates with operational changes at the treatment plant. Specifically, we observe episodic dominance of a likely motile clade during periods of elevated effluent turbidity, as well as episodic dominance of closely related nonmotile clades that likely enhance floc formation during periods of low effluent turbidity.

Consensus assessment of the contamination level of publicly available cyanobacterial genomes

Publicly available genomes are crucial for phylogenetic and metagenomic studies, in which contaminating sequences can be the cause of major problems. This issue is expected to be especially important for Cyanobacteria because axenic strains are notoriously difficult to obtain and keep in culture. Yet, despite their great scientific interest, no data are currently available concerning the quality of publicly available cyanobacterial genomes. As reliably detecting contaminants is a complex task, we designed a pipeline combining six methods in a consensus strategy to assess the contamination level of 440 genome assemblies of Cyanobacteria. Two methods are based on published reference databases of ribosomal genes (SSU rRNA 16S and ribosomal proteins), one is indirectly based on a reference database of marker genes (CheckM), and three are based on complete genome analysis. Among those genome-wide methods, Kraken and DIAMOND blastx share the same reference database that we derived from Ensembl Bacteria, whereas CONCOCT does not require any reference database, instead relying on differences in DNA tetramer frequencies. Given that all the six methods appear to have their own strengths and limitations, we used the consensus of their rankings to infer that >5% of cyanobacterial genome assemblies are highly contaminated by foreign DNA (i.e., contaminants were detected by 5 or 6 methods). Our results will help researchers to check the quality of publicly available genomic data before use in their own analyses. Moreover, we argue that journals should make mandatory the submission of raw read data along with genome assemblies in order to facilitate the detection of contaminants in sequence databases.

Testing culture purity in prokaryotes: criteria and challenges

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

Characterization of Bathyarchaeota genomes assembled from metagenomes of biofilms residing in mesophilic and thermophilic biogas reactors

BACKGROUND

Previous studies on the Miscellaneous Crenarchaeota Group, recently assigned to the novel archaeal phylum Bathyarchaeota, reported on the dominance of these Archaea within the anaerobic carbohydrate cycle performed by the deep marine biosphere. For the first time, members of this phylum were identified also in mesophilic and thermophilic biogas-forming biofilms and characterized in detail.

RESULTS

Metagenome shotgun libraries of biofilm microbiomes were sequenced using the Illumina MiSeq system. Taxonomic classification revealed that between 0.1 and 2% of all classified sequences were assigned to Bathyarchaeota. Individual metagenome assemblies followed by genome binning resulted in the reconstruction of five metagenome-assembled genomes (MAGs) of Bathyarchaeota. MAGs were estimated to be 65-92% complete, ranging in their genome sizes from 1.1 to 2.0 Mb. Phylogenetic classification based on core gene sets confirmed their placement within the phylum Bathyarchaeota clustering as a separate group diverging from most of the recently known Bathyarchaeota clusters. The genetic repertoire of these MAGs indicated an energy metabolism based on carbohydrate and amino acid fermentation featuring the potential for extracellular hydrolysis of cellulose, cellobiose as well as proteins. In addition, corresponding transporter systems were identified. Furthermore, genes encoding enzymes for the utilization of carbon monoxide and/or carbon dioxide via the Wood-Ljungdahl pathway were detected.

CONCLUSIONS

For the members of Bathyarchaeota detected in the biofilm microbiomes, a hydrolytic lifestyle is proposed. This is the first study indicating that Bathyarchaeota members contribute presumably to hydrolysis and subsequent fermentation of organic substrates within biotechnological biogas production processes.

Homoacetogenesis in Deep-Sea Chloroflexi, as Inferred by Single-Cell Genomics, Provides a Link to Reductive Dehalogenation in Terrestrial Dehalococcoidetes

The deep marine subsurface is one of the largest unexplored biospheres on Earth and is widely inhabited by members of the phylum Chloroflexi. In this report, we investigated genomes of single cells obtained from deep-sea sediments of the Peruvian Margin, which are enriched in such Chloroflexi. 16S rRNA gene sequence analysis placed two of these single-cell-derived genomes (DscP3 and Dsc4) in a clade of subphylum I Chloroflexi which were previously recovered from deep-sea sediment in the Okinawa Trough and a third (DscP2-2) as a member of the previously reported DscP2 population from Peruvian Margin site 1230. The presence of genes encoding enzymes of a complete Wood-Ljungdahl pathway, glycolysis/gluconeogenesis, a Rhodobacter nitrogen fixation (Rnf) complex, glyosyltransferases, and formate dehydrogenases in the single-cell genomes of DscP3 and Dsc4 and the presence of an NADH-dependent reduced ferredoxin:NADP oxidoreductase (Nfn) and Rnf in the genome of DscP2-2 imply a homoacetogenic lifestyle of these abundant marine Chloroflexi. We also report here the first complete pathway for anaerobic benzoate oxidation to acetyl coenzyme A (CoA) in the phylum Chloroflexi (DscP3 and Dsc4), including a class I benzoyl-CoA reductase. Of remarkable evolutionary significance, we discovered a gene encoding a formate dehydrogenase (FdnI) with reciprocal closest identity to the formate dehydrogenase-like protein (complex iron-sulfur molybdoenzyme [CISM], DET0187) of terrestrial Dehalococcoides/Dehalogenimonas spp. This formate dehydrogenase-like protein has been shown to lack formate dehydrogenase activity in Dehalococcoides/Dehalogenimonas spp. and is instead hypothesized to couple HupL hydrogenase to a reductive dehalogenase in the catabolic reductive dehalogenation pathway. This finding of a close functional homologue provides an important missing link for understanding the origin and the metabolic core of terrestrial Dehalococcoides/Dehalogenimonas spp. and of reductive dehalogenation, as well as the biology of abundant deep-sea Chloroflexi.

The deep marine subsurface is one of the largest unexplored biospheres on Earth and is widely inhabited by members of the phylum Chloroflexi. In this report, we investigated genomes of single cells obtained from deep-sea sediments and provide evidence for a homacetogenic lifestyle of these abundant marine Chloroflexi. Moreover, genome signature and key metabolic genes indicate an evolutionary relationship between these deep-sea sediment microbes and terrestrial, reductively dehalogenating Dehalococcoides.

Genomic Variation and Evolution of Vibrio parahaemolyticus ST36 over the Course of a Transcontinental Epidemic Expansion

Vibrio parahaemolyticus is the leading cause of seafood-related infections with illnesses undergoing a geographic expansion. In this process of expansion, the most fundamental change has been the transition from infections caused by local strains to the surge of pandemic clonal types. Pandemic clone sequence type 3 (ST3) was the only example of transcontinental spreading until 2012, when ST36 was detected outside the region where it is endemic in the U.S. Pacific Northwest causing infections along the U.S. northeast coast and Spain. Here, we used genome-wide analyses to reconstruct the evolutionary history of the V. parahaemolyticus ST36 clone over the course of its geographic expansion during the previous 25 years. The origin of this lineage was estimated to be in ~1985. By 1995, a new variant emerged in the region and quickly replaced the old clone, which has not been detected since 2000. The new Pacific Northwest (PNW) lineage was responsible for the first cases associated with this clone outside the Pacific Northwest region. After several introductions into the northeast coast, the new PNW clone differentiated into a highly dynamic group that continues to cause illness on the northeast coast of the United States. Surprisingly, the strains detected in Europe in 2012 diverged from this ancestral group around 2000 and have conserved genetic features present only in the old PNW lineage. Recombination was identified as the major driver of diversification, with some preliminary observations suggesting a trend toward a more specialized lifestyle, which may represent a critical element in the expansion of epidemics under scenarios of coastal warming.IMPORTANCEVibrio parahaemolyticus and Vibrio cholerae represent the only two instances of pandemic expansions of human pathogens originating in the marine environment. However, while the current pandemic of V. cholerae emerged more than 50 years ago, the global expansion of V. parahaemolyticus is a recent phenomenon. These modern expansions provide an exceptional opportunity to study the evolutionary process of these pathogens at first hand and gain an understanding of the mechanisms shaping the epidemic dynamics of these diseases, in particular, the emergence, dispersal, and successful introduction in new regions facilitating global spreading of infections. In this study, we used genomic analysis to examine the evolutionary divergence that has occurred over the course of the most recent transcontinental expansion of a pathogenic Vibrio, the spreading of the V. parahaemolyticus sequence type 36 clone from the region where it is endemic on the Pacific coast of North America to the east coast of the United States and finally to the west coast of Europe.

Analysis of single-cell genome sequences of bacteria and archaea

Single-cell genome sequencing of individual archaeal and bacterial cells is a vital approach to decipher the genetic makeup of uncultured microorganisms. With this review, we describe single-cell genome analysis with a focus on the unique properties of single-cell sequence data and with emphasis on quality assessment and assurance.

Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea

Standards for sequencing the microbial 'uncultivated majority', namely bacterial and archaeal single-cell genome sequences, and genome sequences from metagenomic datasets, are proposed.

We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.