VizBin - an application for reference-independent visualization and human-augmented binning of metagenomic data

Simultaneous Single-Cell Genome and Transcriptome Sequencing of Termite Hindgut Protists Reveals Metabolic and Evolutionary Traits of Their Endosymbionts

Some of the protist species which colonize the hindguts of wood-feeding Reticulitermes termites are associated with endosymbiotic bacteria belonging to the genus Endomicrobium. In this study, we focused on the endosymbionts of three protist species from Reticulitermes flavipes, as follows: Pyrsonympha vertens, Trichonympha agilis, and Dinenympha species II. Since these protist hosts represented members of different taxa which colonize separate niches within the hindguts of their termite hosts, we investigated if these differences translated to differential gene content and expression in their endosymbionts. Following assembly and comparative genome and transcriptome analyses, we discovered that these endosymbionts differed with respect to some possible niche-specific traits, such as carbon metabolism. Our analyses suggest that species-specific genes related to carbon metabolism were acquired by horizontal gene transfer (HGT) and may have come from taxa which are common in the termite hind gut. In addition, our analyses suggested that these endosymbionts contain and express genes related to natural transformation (competence) and recombination. Taken together, the presence of genes acquired by HGT and a putative competence pathway suggest that these endosymbionts are not cut off from gene flow and that competence may be a mechanism by which members of Endomicrobium can acquire new traits.

IMPORTANCE The composition and structure of wood, which contains cellulose, hemicellulose, and lignin, prevent most organisms from using this common food source. Termites are a rare exception among animals, and they rely on a complex microbiota housed in their hindguts to use wood as a source of food. The lower termite, Reticulitermes flavipes, houses a variety of protists and prokaryotes that are the key players in the disassembly of lignocellulose. Here, we describe the genomes and the gene expression profiles of five Endomicrobium endosymbionts living inside three different protist species from R. flavipes. Data from these genomes suggest that these Endomicrobium species have different mechanisms for using carbon. In addition, they harbor genes that may be used to import DNA from their environment. This process of DNA uptake may contribute to the high levels of horizontal gene transfer noted previously in Endomicrobium species.

De Novo Assembly and Annotation of the Vaginal Metatranscriptome Associated with Bacterial Vaginosis

The vaginal microbiome plays an important role in women’s health and disease. Here we reanalyzed 40 vaginal transcriptomes from a previous study of de novo assembly (metaT-Assembly) followed by functional annotation. We identified 286,293 contigs and further assigned them to 25 phyla, 209 genera, and 339 species. Lactobacillus iners and Lactobacillus crispatus dominated the microbiome of non-bacterial vaginosis (BV) samples, while a complex of microbiota was identified from BV-associated samples. The metaT-Assembly identified a higher number of bacterial species than the 16S rRNA amplicon and metaT-Kraken methods. However, metaT-Assembly and metaT-Kraken exhibited similar major bacterial composition at the species level. Binning of metatranscriptome data resulted in 176 bins from major known bacteria and several unidentified bacteria in the vagina. Functional analyses based on Clusters of Orthologous Genes (COGs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways suggested that a higher number of transcripts were expressed by the microbiome complex in the BV-associated samples than in non-BV-associated samples. The KEGG pathway analysis with an individual bacterial genome identified specific functions of the identified bacterial genome. Taken together, we demonstrated that the metaT-Assembly approach is an efficient tool to understand the dynamic microbial communities and their functional roles associated with the human vagina.

Microbiomes of air dust collected during the ground-based closed bioregenerative life support experiment "Lunar Palace 365"

BACKGROUND

Understanding the dynamics of airborne microbial communities and antibiotic resistance genes (ARGs) in space life support systems is important because potential pathogens and antibiotic resistance pose a health risk to crew that can lead to mission failure. There have been few reports on the distribution patterns of microbiomes and ARGs in space life support systems. In particular, there have been no detailed investigations of microbiomes and/or antibiotic resistance based on molecular methods in long-term confined bioregenerative life support systems (BLSSs). Therefore, in the present study, we collected air dust samples from two crew shifts, different areas, and different time points in the "Lunar Palace 365" experiment. We evaluated microbial diversity, species composition, functional potential, and antibiotic resistance by combining cultivation-independent analyses (amplicon, shot-gun sequencing, and qPCR).

RESULTS

We found that the bacterial community diversity in the Lunar Palace1 (LP1) system was higher than that in a controlled environment but lower than that in an open environment. Personnel exchange led to significant differences in bacterial community diversity, and source tracking analysis revealed that most bacteria in the air derived from the cabin crew and plants, but no differences in microbial function or antibiotic resistance were observed. Thus, human presence had the strongest effect on the succession of microbial diversity in the BLSSs.

CONCLUSIONS

Our results highlight that microbial diversity in BLSSs is heavily influenced by changes in crew and is unique from other open and controlled environments. Our findings can be used to help develop safe, enclosed BLSS that meet the requirements of human survival and habitation in outer space. In addition, our results can further enhance our understanding of the indoor air microbial community and effectively maintain a safe working and living environment, including plant growth.

Maternal and early life exposures and their potential to influence development of the microbiome

At the dawn of the twentieth century, the medical care of mothers and children was largely relegated to family members and informally trained birth attendants. As the industrial era progressed, early and key public health observations among women and children linked the persistence of adverse health outcomes to poverty and poor nutrition. In the time hence, numerous studies connecting genetics ("nature") to public health and epidemiologic data on the role of the environment ("nurture") have yielded insights into the importance of early life exposures in relation to the occurrence of common diseases, such as diabetes, allergic and atopic disease, cardiovascular disease, and obesity. As a result of these parallel efforts in science, medicine, and public health, the developing brain, immune system, and metabolic physiology are now recognized as being particularly vulnerable to poor nutrition and stressful environments from the start of pregnancy to 3 years of age. In particular, compelling evidence arising from a diverse array of studies across mammalian lineages suggest that modifications to our metagenome and/or microbiome occur following certain environmental exposures during pregnancy and lactation, which in turn render risk of childhood and adult diseases. In this review, we will consider the evidence suggesting that development of the offspring microbiome may be vulnerable to maternal exposures, including an analysis of the data regarding the presence or absence of a low-biomass intrauterine microbiome.

Sulfate differentially stimulates but is not respired by diverse anaerobic methanotrophic archaea

Sulfate-coupled anaerobic oxidation of methane (AOM) is a major methane sink in marine sediments. Multiple lineages of anaerobic methanotrophic archaea (ANME) often coexist in sediments and catalyze this process syntrophically with sulfate-reducing bacteria (SRB), but the potential differences in ANME ecophysiology and mechanisms of syntrophy remain unresolved. A humic acid analog, anthraquinone 2,6-disulfonate (AQDS), could decouple archaeal methanotrophy from bacterial sulfate reduction and serve as the terminal electron acceptor for AOM (AQDS-coupled AOM). Here in sediment microcosm experiments, we examined variations in physiological response between two co-occurring ANME-2 families (ANME-2a and ANME-2c) and tested the hypothesis of sulfate respiration by ANME-2. Sulfate concentrations as low as 100 µM increased AQDS-coupled AOM nearly 2-fold matching the rates of sulfate-coupled AOM. However, the SRB partners remained inactive in microcosms with sulfate and AQDS and neither ANME-2 families respired sulfate, as shown by their cellular sulfur contents and anabolic activities measured using nanoscale secondary ion mass spectrometry. ANME-2a anabolic activity was significantly higher than ANME-2c, suggesting that ANME-2a was primarily responsible for the observed sulfate stimulation of AQDS-coupled AOM. Comparative transcriptomics showed significant upregulation of ANME-2a transcripts linked to multiple ABC transporters and downregulation of central carbon metabolism during AQDS-coupled AOM compared to sulfate-coupled AOM. Surprisingly, genes involved in sulfur anabolism were not differentially expressed during AQDS-coupled AOM with and without sulfate amendment. Collectively, this data indicates that ANME-2 archaea are incapable of respiring sulfate, but sulfate availability differentially stimulates the growth and AOM activity of different ANME lineages.

Reconstruction of Archaeal Genomes from Short-Read Metagenomes

As the majority of biological diversity remains unexplored and uncultured, investigating it requires culture-independent approaches. Archaea in particular suffer from a multitude of issues that make their culturing problematic, from them being frequently members of the rare biosphere, to low growth rates, to them thriving under very specific and often extreme environmental and community conditions that are difficult to replicate. OMICs techniques are state of the art approaches that allow direct high-throughput investigations of environmental samples at all levels from nucleic acids to proteins, lipids, and secondary metabolites. Metagenomics, as the foundation for other OMICs techniques, facilitates the identification and functional characterization of the microbial community members and can be combined with other methods to provide insights into the microbial activities, both on the RNA and protein levels. In this chapter, we provide a step-by-step workflow for the recovery of archaeal genomes from metagenomes, starting from raw short-read sequences. This workflow can be applied to recover bacterial genomes as well.

Comparative genomics reveals electron transfer and syntrophic mechanisms differentiating methanotrophic and methanogenic archaea

The anaerobic oxidation of methane coupled to sulfate reduction is a microbially mediated process requiring a syntrophic partnership between anaerobic methanotrophic (ANME) archaea and sulfate-reducing bacteria (SRB). Based on genome taxonomy, ANME lineages are polyphyletic within the phylum Halobacterota, none of which have been isolated in pure culture. Here, we reconstruct 28 ANME genomes from environmental metagenomes and flow sorted syntrophic consortia. Together with a reanalysis of previously published datasets, these genomes enable a comparative analysis of all marine ANME clades. We review the genomic features that separate ANME from their methanogenic relatives and identify what differentiates ANME clades. Large multiheme cytochromes and bioenergetic complexes predicted to be involved in novel electron bifurcation reactions are well distributed and conserved in the ANME archaea, while significant variations in the anabolic C1 pathways exists between clades. Our analysis raises the possibility that methylotrophic methanogenesis may have evolved from a methanotrophic ancestor.

A comparative genomics study of anaerobic methanotrophic (ANME) archaea reveals the genetic "parts list" associated with the repeated evolutionary transition between methanogenic and methanotrophic metabolism in the archaeal domain of life.

Isolation from a fish kill and transcriptomic characterization of Gyrodinium jinhaense off Long Island Sound

First found in Korean coastal water, the dinoflagellate Gyrodinium jinhaense is a recently established species with unclear global distribution and unexplored genomic characteristics. From a laboratory fish mortality event off Long Island Sound, USA, we isolated a dinoflagellate, and by microscopic and molecular (18S rRNA gene; >99% identical) analyses found that it resembles G. jinhaense, hence named G. jinhaense strain AP17. Towards developing a genetic database for this dinoflagellate, a transcriptome of this species was sequenced using RNA-seq, producing 6 Gbp of data that was assembled into over 70,000 unigenes. The assembled transcriptome GC content was approximately 56% and the total Benchmarking Universal Single-Copy Orthologs for Eukaryota and Alveolata databases were approximately 50% and 57%, respectfully. Genes involved in grazing, energy generation, genome architecture, and protein synthesis, processing, and degradation were highly represented in the transcriptome. Moreover, fragments of polyketide synthase and saxitoxin genes were found but saxitoxins were not detected in high performance liquid chromatography measurements. With the first reported transcriptome for the Gyrodinium genus, this study will serve as a baseline for future Gyrodinium genomics and toxicological studies.

Standard Sample Storage Conditions Have an Impact on Inferred Microbiome Composition and Antimicrobial Resistance Patterns

Storage of biological specimens is crucial in the life and medical sciences. Storage conditions for samples can be different for a number of reasons, and it is unclear what effect this can have on the inferred microbiome composition in metagenomics analyses. Here, we assess the effect of common storage temperatures (deep freezer, -80°C; freezer, -20°C; refrigerator, 5°C; room temperature, 22°C) and storage times (immediate sample processing, 0 h; next day, 16 h; over weekend, 64 h; longer term, 4, 8, and 12 months) as well as repeated sample freezing and thawing (2 to 4 freeze-thaw cycles). We examined two different pig feces and sewage samples, unspiked and spiked with a mock community, in triplicate, respectively, amounting to a total of 438 samples (777 Gbp; 5.1 billion reads). Storage conditions had a significant and systematic effect on the taxonomic and functional composition of microbiomes. Distinct microbial taxa and antimicrobial resistance classes were, in some situations, similarly affected across samples, while others were not, suggesting an impact of individual inherent sample characteristics. With an increasing number of freeze-thaw cycles, an increasing abundance of Firmicutes, Actinobacteria, and eukaryotic microorganisms was observed. We provide recommendations for sample storage and strongly suggest including more detailed information in the metadata together with the DNA sequencing data in public repositories to better facilitate meta-analyses and reproducibility of findings. IMPORTANCE Previous research has reported effects of DNA isolation, library preparation, and sequencing technology on metagenomics-based microbiome composition; however, the effect of biospecimen storage conditions has not been thoroughly assessed. We examined the effect of common sample storage conditions on metagenomics-based microbiome composition and found significant and, in part, systematic effects. Repeated freeze-thaw cycles could be used to improve the detection of microorganisms with more rigid cell walls, including parasites. We provide a data set that could also be used for benchmarking algorithms to identify and correct for unwanted batch effects. Overall, the findings suggest that all samples of a microbiome study should be stored in the same way. Furthermore, there is a need to mandate more detailed information about sample storage and processing be published together with DNA sequencing data at the International Nucleotide Sequence Database Collaboration (ENA/EBI, NCBI, DDBJ) or other repositories.

Genome Streamlining, Plasticity, and Metabolic Versatility Distinguish Co-occurring Toxic and Nontoxic Cyanobacterial Strains of Microcoleus

Harmful cyanobacterial bloom occurrences have increased worldwide due to climate change and eutrophication, causing nuisance and animal deaths. Species from the benthic cyanobacterial genus Microcoleus are ubiquitous and form thick mats in freshwater systems, such as rivers, that are sometimes toxic due to the production of potent neurotoxins (anatoxins). Anatoxin-producing (toxic) strains typically coexist with non-anatoxin-producing (nontoxic) strains in mats, although the reason for this is unclear. To determine the genetic mechanisms differentiating toxic and nontoxic Microcoleus, we sequenced and assembled genomes from 11 cultures and compared these to another 31 Microcoleus genomes. Average nucleotide identities (ANI) indicate that toxic and nontoxic strains are distinct species (ANI, <95%), and only 6% of genes are shared across all 42 genomes, suggesting a high level of genetic divergence among Microcoleus strains. Comparative genomics showed substantial genome streamlining in toxic strains and a potential dependency on external sources for thiamine and sucrose. Toxic and nontoxic strains are further differentiated by an additional set of putative nitrate transporter (nitrogen uptake) and cyanophycin (carbon and nitrogen storage) genes, respectively. These genes likely confer distinct competitive advantages based on nutrient availability and suggest nontoxic strains are more robust to nutrient fluctuations. Nontoxic strains also possess twice as many transposable elements, potentially facilitating greater genetic adaptation to environmental changes. Our results offer insights into the divergent evolution of Microcoleus strains and the potential for cooperative and competitive interactions that contribute to the co-occurrence of toxic and nontoxic species within mats. IMPORTANCE Microcoleus autumnalis, and closely related Microcoleus species, compose a geographically widespread group of freshwater benthic cyanobacteria. Canine deaths due to anatoxin-a poisoning, following exposure to toxic proliferations, have been reported globally. While Microcoleus proliferations are on the rise, the mechanisms underpinning competition between, or coexistence of, toxic and nontoxic strains are unknown. This study identifies substantial genetic differences between anatoxin-producing and non-anatoxin-producing strains, pointing to reduced metabolic flexibility in toxic strains, and potential dependence on cohabiting nontoxic strains. Results provide insights into the metabolic and evolutionary differences between toxic and nontoxic Microcoleus, which may assist in predicting and managing aquatic proliferations.

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."

Expanding the repertoire of electron acceptors for the anaerobic oxidation of methane in carbonates in the Atlantic and Pacific Ocean

Authigenic carbonates represent a significant microbial sink for methane, yet little is known about the microbiome responsible for the methane removal. We identify carbonate microbiomes distributed over 21 locations hosted by seven different cold seeps in the Pacific and Atlantic Oceans by carrying out a gene-based survey using 16S rRNA- and mcrA gene sequencing coupled with metagenomic analyses. Based on 16S rRNA gene amplicon analyses, these sites were dominated by bacteria affiliated to the Firmicutes, Alpha- and Gammaproteobacteria. ANME-1 and -2 archaeal clades were abundant in the carbonates yet their typical syntrophic partners, sulfate-reducing bacteria, were not significantly present. Based on mcrA amplicon analyses, the Candidatus Methanoperedens clades were also highly abundant. Our metagenome analysis indicated that methane oxidizers affiliated to the ANME-1 and -2, may be capable of performing complete methane- and potentially short-chain alkane oxidation independently using oxidized sulfur and nitrogen compounds as terminal electron acceptors. Gammaproteobacteria are hypothetically capable of utilizing oxidized nitrogen compounds and may be involved in syntrophy with methane-oxidizing archaea. Carbonate structures represent a window for a more diverse utilization of electron acceptors for anaerobic methane oxidation along the Atlantic and Pacific Margin.

Genomic adaptations enabling Acidithiobacillus distribution across wide-ranging hot spring temperatures and pHs

BACKGROUND

Terrestrial hot spring settings span a broad spectrum of physicochemistries. Physicochemical parameters, such as pH and temperature, are key factors influencing differences in microbial composition across diverse geothermal areas. Nonetheless, analysis of hot spring pools from the Taupo Volcanic Zone (TVZ), New Zealand, revealed that some members of the bacterial genus, Acidithiobacillus, are prevalent across wide ranges of hot spring pHs and temperatures. To determine the genomic attributes of Acidithiobacillus that inhabit such diverse conditions, we assembled the genomes of 19 uncultivated hot spring Acidithiobacillus strains from six geothermal areas and compared these to 37 publicly available Acidithiobacillus genomes from various habitats.

RESULTS

Analysis of 16S rRNA gene amplicons from 138 samples revealed that Acidithiobacillus comprised on average 11.4 ± 16.8% of hot spring prokaryotic communities, with three Acidithiobacillus amplicon sequence variants (ASVs) (TVZ_G1, TVZ_G2, TVZ_G3) accounting for > 90% of Acidithiobacillus in terms of relative abundance, and occurring in 126 out of 138 samples across wide ranges of temperature (17.5-92.9 °C) and pH (1.0-7.5). We recovered 19 environmental genomes belonging to each of these three ASVs, as well as a fourth related group (TVZ_G4). Based on genome average nucleotide identities, the four groups (TVZ_G1-TVZ_G4) constitute distinct species (ANI < 96.5%) of which three are novel Acidithiobacillus species (TVZ_G2-TVZ_G4) and one belongs to Acidithiobacillus caldus (TVZ_G1). All four TVZ Acidithiobacillus groups were found in hot springs with temperatures above the previously known limit for the genus (up to 40 °C higher), likely due to significantly higher proline and GC contents than other Acidithiobacillus species, which are known to increase thermostability. Results also indicate hot spring-associated Acidithiobacillus have undergone genome streamlining, likely due to thermal adaptation. Moreover, our data suggest that Acidithiobacillus prevalence across varied hot spring pHs is supported by distinct strategies, whereby TVZ_G2-TVZ_G4 regulate pH homeostasis mostly through Na⁺/H⁺ antiporters and proton-efflux ATPases, whereas TVZ_G1 mainly relies on amino acid decarboxylases.

CONCLUSIONS

This study provides insights into the distribution of Acidithiobacillus species across diverse hot spring physichochemistries and determines genomic features and adaptations that potentially enable Acidithiobacillus species to colonize a broad range of temperatures and pHs in geothermal environments. Video Abstract.

Acquisition and Adaptation of Ultra-small Parasitic Reduced Genome Bacteria to Mammalian Hosts

The first cultivated representative of the enigmatic phylum Saccharibacteria (formerly TM7) was isolated from humans and revealed an ultra-small cell size (200-300 nm), a reduced genome with limited biosynthetic capabilities, and a unique parasitic lifestyle. TM7x was the only cultivated member of the candidate phyla radiation (CPR), estimated to encompass 26% of the domain Bacteria. Here we report on divergent genomes from major lineages across the Saccharibacteria phylum in humans and mammals, as well as from ancient dental calculus. These lineages are present at high prevalence within hosts. Direct imaging reveals that all groups are ultra-small in size, likely feeding off commensal bacteria. Analyses suggest that multiple acquisition events in the past led to the current wide diversity, with convergent evolution of key functions allowing Saccharibacteria from the environment to adapt to mammals. Ultra-small, parasitic CPR bacteria represent a relatively unexplored paradigm of prokaryotic interactions within mammalian microbiomes.

Discovering multiple types of DNA methylation from bacteria and microbiome using nanopore sequencing

Bacterial DNA methylation occurs at diverse sequence contexts and plays important functional roles in cellular defense and gene regulation. Existing methods for detecting DNA modification from nanopore sequencing data do not effectively support de novo study of unknown bacterial methylomes. In this work, we observed that a nanopore sequencing signal displays complex heterogeneity across methylation events of the same type. To enable nanopore sequencing for broadly applicable methylation discovery, we generated a training dataset from an assortment of bacterial species and developed a method, named nanodisco ( https://github.com/fanglab/nanodisco ), that couples the identification and fine mapping of the three forms of methylation into a multi-label classification framework. We applied it to individual bacteria and the mouse gut microbiome for reliable methylation discovery. In addition, we demonstrated the use of DNA methylation for binning metagenomic contigs, associating mobile genetic elements with their host genomes and identifying misassembled metagenomic contigs.

Persistence of birth mode-dependent effects on gut microbiome composition, immune system stimulation and antimicrobial resistance during the first year of life

Caesarean section delivery (CSD) disrupts mother-to-neonate transmission of specific microbial strains and functional repertoires as well as linked immune system priming. Here we investigate whether differences in microbiome composition and impacts on host physiology persist at 1 year of age. We perform high-resolution, quantitative metagenomic analyses of the gut microbiomes of infants born by vaginal delivery (VD) or by CSD, from immediately after birth through to 1 year of life. Several microbial populations show distinct enrichments in CSD-born infants at 1 year of age including strains of Bacteroides caccae, Bifidobacterium bifidum and Ruminococcus gnavus, whereas others are present at higher levels in the VD group including Faecalibacterium prausnitizii, Bifidobacterium breve and Bifidobacterium kashiwanohense. The stimulation of healthy donor-derived primary human immune cells with LPS isolated from neonatal stool samples results in higher levels of tumour necrosis factor alpha (TNF-α) in the case of CSD extracts over time, compared to extracts from VD infants for which no such changes were observed during the first year of life. Functional analyses of the VD metagenomes at 1 year of age demonstrate a significant increase in the biosynthesis of the natural antibiotics, carbapenem and phenazine. Concurrently, we find antimicrobial resistance (AMR) genes against several classes of antibiotics in both VD and CSD. The abundance of AMR genes against synthetic (including semi-synthetic) agents such as phenicol, pleuromutilin and diaminopyrimidine are increased in CSD children at day 5 after birth. In addition, we find that mobile genetic elements, including phages, encode AMR genes such as glycopeptide, diaminopyrimidine and multidrug resistance genes. Our results demonstrate persistent effects at 1 year of life resulting from birth mode-dependent differences in earliest gut microbiome colonisation.

SLIDR and SLOPPR: flexible identification of spliced leader trans-splicing and prediction of eukaryotic operons from RNA-Seq data

BACKGROUND

Spliced leader (SL) trans-splicing replaces the 5' end of pre-mRNAs with the spliced leader, an exon derived from a specialised non-coding RNA originating from elsewhere in the genome. This process is essential for resolving polycistronic pre-mRNAs produced by eukaryotic operons into monocistronic transcripts. SL trans-splicing and operons may have independently evolved multiple times throughout Eukarya, yet our understanding of these phenomena is limited to only a few well-characterised organisms, most notably C. elegans and trypanosomes. The primary barrier to systematic discovery and characterisation of SL trans-splicing and operons is the lack of computational tools for exploiting the surge of transcriptomic and genomic resources for a wide range of eukaryotes.

RESULTS

Here we present two novel pipelines that automate the discovery of SLs and the prediction of operons in eukaryotic genomes from RNA-Seq data. SLIDR assembles putative SLs from 5' read tails present after read alignment to a reference genome or transcriptome, which are then verified by interrogating corresponding SL RNA genes for sequence motifs expected in bona fide SL RNA molecules. SLOPPR identifies RNA-Seq reads that contain a given 5' SL sequence, quantifies genome-wide SL trans-splicing events and predicts operons via distinct patterns of SL trans-splicing events across adjacent genes. We tested both pipelines with organisms known to carry out SL trans-splicing and organise their genes into operons, and demonstrate that (1) SLIDR correctly detects expected SLs and often discovers novel SL variants; (2) SLOPPR correctly identifies functionally specialised SLs, correctly predicts known operons and detects plausible novel operons.

CONCLUSIONS

SLIDR and SLOPPR are flexible tools that will accelerate research into the evolutionary dynamics of SL trans-splicing and operons throughout Eukarya and improve gene discovery and annotation for a wide range of eukaryotic genomes. Both pipelines are implemented in Bash and R and are built upon readily available software commonly installed on most bioinformatics servers. Biological insight can be gleaned even from sparse, low-coverage datasets, implying that an untapped wealth of information can be retrieved from existing RNA-Seq datasets as well as from novel full-isoform sequencing protocols as they become more widely available.

A protocol combining breath testing and ex vivo fermentations to study the human gut microbiome

This protocol describes the application of breath testing and ex vivo fermentations to study the association between breath methane and the composition and functionality of the gut microbiome. The protocol provides a useful systems biology approach for studying the gut microbiome in humans, which combines standardized methods in human breath testing and fecal sampling. The model described is accessible and easy to repeat, but its relative simplicity means that it can deviate from human physiological conditions.

Assembly of Bacterial Genome Sequences from Metagenomes of Spacecraft Assembly Cleanrooms

Characterizing the microbiome of spacecraft assembly cleanrooms is important for planetary protection. We report two bacterial metagenome-assembled genome sequences (MAGs) reconstructed from metagenomes produced from cleanroom samples from the Kennedy Space Center’s Payload Hazardous Servicing Facility (KSC-PHSF) during the handling of the Phoenix spacecraft. Characterization of these MAGs will enable identification of the strategies underpinning their survival.

Characterizing the microbiome of spacecraft assembly cleanrooms is important for planetary protection. We report two bacterial metagenome-assembled genomes (MAGs) reconstructed from metagenomes produced from cleanroom samples from the Kennedy Space Center’s Payload Hazardous Servicing Facility (KSC-PHSF) during the handling of the Phoenix spacecraft. Characterization of these MAGs will enable identification of the strategies underpinning their survival.

Genome-scale phylogenetic analyses confirm Olpidium as the closest living zoosporic fungus to the non-flagellated, terrestrial fungi

The zoosporic obligate endoparasites, Olpidium, hold a pivotal position to the reconstruction of the flagellum loss in fungi, one of the key morphological transitions associated with the colonization of land by the early fungi. We generated genome and transcriptome data from non-axenic zoospores of Olpidium bornovanus and used a metagenome approach to extract phylogenetically informative fungal markers. Our phylogenetic reconstruction strongly supported Olpidium as the closest zoosporic relative of the non-flagellated terrestrial fungi. Super-alignment analyses resolved Olpidium as sister to the non-flagellated terrestrial fungi, whereas a super-tree approach recovered different placements of Olpidium, but without strong support. Further investigations detected little conflicting signal among the sampled markers but revealed a potential polytomy in early fungal evolution associated with the branching order among Olpidium, Zoopagomycota and Mucoromycota. The branches defining the evolutionary relationships of these lineages were characterized by short branch lengths and low phylogenetic content and received equivocal support for alternative phylogenetic hypotheses from individual markers. These nodes were marked by important morphological innovations, including the transition to hyphal growth and the loss of flagellum, which enabled early fungi to explore new niches and resulted in rapid and temporally concurrent Precambrian diversifications of the ancestors of several phyla of fungi.

Reconstructing Genomes of Carbon Monoxide Oxidisers in Volcanic Deposits Including Members of the Class Ktedonobacteria

Microorganisms can potentially colonise volcanic rocks using the chemical energy in reduced gases such as methane, hydrogen (H₂) and carbon monoxide (CO). In this study, we analysed soil metagenomes from Chilean volcanic soils, representing three different successional stages with ages of 380, 269 and 63 years, respectively. A total of 19 metagenome-assembled genomes (MAGs) were retrieved from all stages with a higher number observed in the youngest soil (1640: 2 MAGs, 1751: 1 MAG, 1957: 16 MAGs). Genomic similarity indices showed that several MAGs had amino-acid identity (AAI) values >50% to the phyla Actinobacteria, Acidobacteria, Gemmatimonadetes, Proteobacteria and Chloroflexi. Three MAGs from the youngest site (1957) belonged to the class Ktedonobacteria (Chloroflexi). Complete cellular functions of all the MAGs were characterised, including carbon fixation, terpenoid backbone biosynthesis, formate oxidation and CO oxidation. All 19 environmental genomes contained at least one gene encoding a putative carbon monoxide dehydrogenase (CODH). Three MAGs had form I coxL operon (encoding the large subunit CO-dehydrogenase). One of these MAGs (MAG-1957-2.1, Ktedonobacterales) was highly abundant in the youngest soil. MAG-1957-2.1 also contained genes encoding a [NiFe]-hydrogenase and hyp genes encoding accessory enzymes and proteins. Little is known about the Ktedonobacterales through cultivated isolates, but some species can utilise H₂ and CO for growth. Our results strongly suggest that the remote volcanic sites in Chile represent a natural habitat for Ktedonobacteria and they may use reduced gases for growth.

Optimised biomolecular extraction for metagenomic analysis of microbial biofilms from high-mountain streams

Glacier-fed streams (GFS) are harsh ecosystems dominated by microbial life organized in benthic biofilms, yet the biodiversity and ecosystem functions provided by these communities remain under-appreciated. To better understand the microbial processes and communities contributing to GFS ecosystems, it is necessary to leverage high throughput sequencing. Low biomass and high inorganic particle load in GFS sediment samples may affect nucleic acid extraction efficiency using extraction methods tailored to other extreme environments such as deep-sea sediments. Here, we benchmarked the utility and efficacy of four extraction protocols, including an up-scaled phenol-chloroform protocol. We found that established protocols for comparable sample types consistently failed to yield sufficient high-quality DNA, delineating the extreme character of GFS. The methods differed in the success of downstream applications such as library preparation and sequencing. An adapted phenol-chloroform-based extraction method resulted in higher yields and better recovered the expected taxonomic profile and abundance of reconstructed genomes when compared to commercially-available methods. Affordable and straight-forward, this method consistently recapitulated the abundance and genomes of a mock community, including eukaryotes. Moreover, by increasing the amount of input sediment, the protocol is readily adjustable to the microbial load of the processed samples without compromising protocol efficiency. Our study provides a first systematic and extensive analysis of the different options for extraction of nucleic acids from glacier-fed streams for high-throughput sequencing applications, which may be applied to other extreme environments.

Termite mounds contain soil-derived methanotroph communities kinetically adapted to elevated methane concentrations

Termite mounds have recently been confirmed to mitigate approximately half of termite methane (CH4) emissions, but the aerobic CH4 oxidising bacteria (methanotrophs) responsible for this consumption have not been resolved. Here, we describe the abundance, composition and CH4 oxidation kinetics of the methanotroph communities in the mounds of three distinct termite species sampled from Northern Australia. Results from three independent methods employed show that methanotrophs are rare members of microbial communities in termite mounds, with a comparable abundance but distinct composition to those of adjoining soil samples. Across all mounds, the most abundant and prevalent methane monooxygenase sequences were affiliated with upland soil cluster α (USCα), with sequences homologous to Methylocystis and tropical upland soil cluster (TUSC) also detected. The reconstruction of a metagenome-assembled genome of a mound USCα representative highlighted the metabolic capabilities of this group of methanotrophs. The apparent Michaelis-Menten kinetics of CH4 oxidation in mounds were estimated from in situ reaction rates. Methane affinities of the communities were in the low micromolar range, which is one to two orders of magnitude higher than those of upland soils, but significantly lower than those measured in soils with a large CH4 source such as landfill cover soils. The rate constant of CH4 oxidation, as well as the porosity of the mound material, were significantly positively correlated with the abundance of methanotroph communities of termite mounds. We conclude that termite-derived CH4 emissions have selected for distinct methanotroph communities that are kinetically adapted to elevated CH4 concentrations. However, factors other than substrate concentration appear to limit methanotroph abundance and hence these bacteria only partially mitigate termite-derived CH4 emissions. Our results also highlight the predominant role of USCα in an environment with elevated CH4 concentrations and suggest a higher functional diversity within this group than previously recognised.

Methanogens Within a High Salinity Oil Reservoir From the Gulf of Mexico

Oil reservoirs contain microbial populations that are both autochthonously and allochthonously introduced by industrial development. These microbial populations are greatly influenced by external factors including, but not limited to, salinity and temperature. In this study, we used metagenomics to examine the microbial populations within five wells of the same hydrocarbon reservoir system in the Gulf of Mexico. These elevated salinity (149–181 ppt salinity, 4–5× salinity of seawater) reservoirs have limited taxonomic and functional microbial diversity dominated by methanogens, Halanaerobium and other Firmicutes lineages, and contained less abundant lineages such as Deltaproteobacteria. Metagenome assembled genomes (MAGs) were generated and analyzed from the various wells. Methanogen MAGs were closely related to Methanohalophilus euhalobius, a known methylotrophic methanogen from a high salinity oil environment. Based on metabolic reconstruction of genomes, the Halanaerobium perform glycine betaine fermentation, potentially produced by the methanogens. Industrial introduction of methanol to prevent methane hydrate formation to this environment is likely to be consumed by these methanogens. As such, this subsurface oil population may represent influences from industrial processes.

Characterization of sponge-associated Verrucomicrobia: microcompartment-based sugar utilization and enhanced toxin-antitoxin modules as features of host-associated Opitutales

Bacteria of the phylum Verrucomicrobia are ubiquitous in marine environments and can be found as free-living organisms or as symbionts of eukaryotic hosts. Little is known about host-associated Verrucomicrobia in the marine environment. Here we reconstructed two genomes of symbiotic Verrucomicrobia from bacterial metagenomes derived from the Atlanto-Mediterranean sponge Petrosia ficiformis and three genomes from strains that we isolated from offshore seawater of the Eastern Mediterranean Sea. Phylogenomic analysis of these five strains indicated that they are all members of Verrucomicrobia subdivision 4, order Opitutales. We compared these novel sponge-associated and seawater-isolated genomes to closely related Verrucomicrobia. Genomic analysis revealed that Planctomycetes-Verrucomicrobia microcompartment gene clusters are enriched in the genomes of symbiotic Opitutales including sponge symbionts but not in free-living ones. We hypothesize that in sponge symbionts these microcompartments are used for degradation of l-fucose and l-rhamnose, which are components of algal and bacterial cell walls and therefore may be found at high concentrations in the sponge tissue. Furthermore, we observed an enrichment of toxin-antitoxin modules in symbiotic Opitutales. We suggest that, in sponges, verrucomicrobial symbionts utilize these modules as a defence mechanism against antimicrobial activity deriving from the abundant microbial community co-inhabiting the host.

Gammaproteobacteria mediating utilization of methyl-, sulfur- and petroleum organic compounds in deep ocean hydrothermal plumes

Deep-sea hydrothermal plumes are considered natural laboratories for understanding ecological and biogeochemical interactions. Previous studies focused on interactions between microorganisms and inorganic, reduced hydrothermal inputs including sulfur, hydrogen, iron, and manganese. However, little is known about transformations of organic compounds, especially methylated, sulfur-containing compounds, and petroleum hydrocarbons. Here, we reconstructed nine gammaproteobacterial metagenome-assembled genomes, affiliated with Methylococcales, Methylophaga, and Cycloclasticus, from three hydrothermal ecosystems. We present evidence that these three groups have high transcriptional activities of genes encoding cycling of C₁-compounds, petroleum hydrocarbons, and organic sulfur in hydrothermal plumes. This includes oxidation of methanethiol, the simplest thermochemically-derived organic sulfur, for energy metabolism in Methylococcales and Cycloclasticus. Together with active transcription of genes for thiosulfate and methane oxidation in Methylococcales, these results suggest an adaptive strategy of versatile and simultaneous use of multiple available electron donors. Meanwhile, the first near-complete MAG of hydrothermal Methylophaga aminisulfidivorans and its transcriptional profile point to active chemotaxis targeting small organic compounds. Petroleum hydrocarbon-degrading Cycloclasticus are abundant and active in plumes of oil spills as well as deep-sea vents, suggesting that they are indigenous and effectively respond to stimulus of hydrocarbons in the deep sea. These findings suggest that these three groups of Gammaproteobacteria transform organic carbon and sulfur compounds via versatile and opportunistic metabolism and modulate biogeochemistry in plumes of hydrothermal systems as well as oil spills, thus contributing broad ecological impact to the deep ocean globally.

Tools for successful proliferation: diverse strategies of nutrient acquisition by a benthic cyanobacterium

Freshwater cyanobacterial blooms have increased worldwide, channeling organic carbon into these systems, and threatening animal health through the production of cyanotoxins. Both toxic and nontoxic Microcoleus proliferations usually occur when there are moderate concentrations of dissolved inorganic nitrogen, but when phosphorus is scarce. In order to understand how Microcoleus establishes thick biofilms (or mats) on riverbeds under phosphorus-limiting conditions, we collected Microcoleus-dominated biofilms over a 19-day proliferation event for proteogenomics. A single pair of nitrogen-dependent Microcoleus species were consistently present in relatively high abundance, although each followed a unique metabolic trajectory. Neither possessed anatoxin gene clusters, and only very low concentrations of anatoxins (~2 µg kg-1) were detected, likely originating from rarer Microcoleus species also present. Proteome allocations were dominated by photosynthesizing cyanobacteria and diatoms, and data indicate biomass was actively recycled by Bacteroidetes and Myxococcales. Microcoleus likely acquired nutrients throughout the proliferation event by uptake of nitrate, urea, and inorganic and organic phosphorus. Both species also harbored genes that could be used for inorganic phosphate solubilization with pyrroloquinoline quinone cofactors produced by cohabiting Proteobacteria. Results indicate that Microcoleus are equipped with diverse mechanisms for nitrogen and phosphorus acquisition, enabling them to proliferate and out-compete others in low-phosphorus waters.

SolidBin: improving metagenome binning with semi-supervised normalized cut

MOTIVATION

Metagenomic contig binning is an important computational problem in metagenomic research, which aims to cluster contigs from the same genome into the same group. Unlike classical clustering problem, contig binning can utilize known relationships among some of the contigs or the taxonomic identity of some contigs. However, the current state-of-the-art contig binning methods do not make full use of the additional biological information except the coverage and sequence composition of the contigs.

RESULTS

We developed a novel contig binning method, Semi-supervised Spectral Normalized Cut for Binning (SolidBin), based on semi-supervised spectral clustering. Using sequence feature similarity and/or additional biological information, such as the reliable taxonomy assignments of some contigs, SolidBin constructs two types of prior information: must-link and cannot-link constraints. Must-link constraints mean that the pair of contigs should be clustered into the same group, while cannot-link constraints mean that the pair of contigs should be clustered in different groups. These constraints are then integrated into a classical spectral clustering approach, normalized cut, for improved contig binning. The performance of SolidBin is compared with five state-of-the-art genome binners, CONCOCT, COCACOLA, MaxBin, MetaBAT and BMC3C on five next-generation sequencing benchmark datasets including simulated multi- and single-sample datasets and real multi-sample datasets. The experimental results show that, SolidBin has achieved the best performance in terms of F-score, Adjusted Rand Index and Normalized Mutual Information, especially while using the real datasets and the single-sample dataset.

AVAILABILITY AND IMPLEMENTATION

https://github.com/sufforest/SolidBin.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Bacterial community in commercial airliner cabins in China

Considering the amount of time that crew members and passengers spend on airliners and the potential health impact of pathogenic bacteria, it is important to understand the population of bacteria inside airliners and the factors affecting the bacterial concentration. This study recorded the species of airborne and cabin surface culturable bacteria inside various airliner. Seven flights ranging from 3 to 5 hours in duration on different types of airliner were chosen. Multiple species of bacteria in the air of the airliners, such as Brachybacterium paraconglomeratum, were identified by means of the 16S ribosomal RNA gene sequencing method, and most of the bacteria were Gram-positive. This study found that the bacterial concentration in the airliners decreases as the relative humidity increases. The decrease in the number of airborne bacteria may be the reason for the reduced occurrence of unwanted symptoms exhibited by passengers.

Applications of weighted association networks applied to compositional data in biology

Next-generation sequencing technologies have generated, and continue to produce, an increasingly large corpus of biological data. The data generated are inherently compositional as they convey only relative information dependent upon the capacity of the instrument, experimental design and technical bias. There is considerable information to be gained through network analysis by studying the interactions between components within a system. Network theory methods using compositional data are powerful approaches for quantifying relationships between biological components and their relevance to phenotype, environmental conditions or other external variables. However, many of the statistical assumptions used for network analysis are not designed for compositional data and can bias downstream results. In this mini-review, we illustrate the utility of network theory in biological systems and investigate modern techniques while introducing researchers to frameworks for implementation. We overview (1) compositional data analysis, (2) data transformations and (3) network theory along with insight on a battery of network types including static-, temporal-, sample-specific- and differential-networks. The intention of this mini-review is not to provide a comprehensive overview of network methods, rather to introduce microbiology researchers to (semi)-unsupervised data-driven approaches for inferring latent structures that may give insight into biological phenomena or abstract mechanics of complex systems.

Novel clostridial lineages recovered from metagenomes of a hot oil reservoir

Oil reservoirs have been shown to house numerous microbial lineages that differ based on the in-situ pH, salinity and temperature of the subsurface environment. Lineages of Firmicutes, including Clostridiales, have been frequently detected in oil reservoirs, but are typically not considered impactful or relevant due to their spore-forming nature. Here we show, using metagenomics, a high temperature oil reservoir of marine salinity contains a microbial population that is predominantly from within the Order Clostridiales. These organisms form an oil-reservoir specific clade based on the phylogenies of both 16S rRNA genes and ribosomal proteins, which we propose to name ^UPetromonas tenebris, meaning they are single-celled organisms from dark rocks. Metagenome-assembled genomes (MAGs) of these Petromonas sp. were obtained and used to determine that these populations, while capable of spore-formation, were also likely replicating in situ in the reservoir. We compared these MAGs to closely related genomes and show that these subsurface Clostridiales differ, from the surface derived genomes, showing signatures of the ability to degrade plant-related compounds, whereas subsurface genomes only show the ability to process simple sugars. The estimation of in-situ replication from genomic data suggest that ^UPetromonas tenebris lineages are functional in-situ and may be specifically adapted to inhabit oil reservoirs.

Genome diversification in globally distributed novel marine Proteobacteria is linked to environmental adaptation

Proteobacteria constitute one of the most diverse and abundant groups of microbes on Earth. In productive marine environments like deep-sea hydrothermal systems, Proteobacteria are implicated in autotrophy coupled to sulfur, methane, and hydrogen oxidation, sulfate reduction, and denitrification. Beyond chemoautotrophy, little is known about the ecological significance of poorly studied Proteobacteria lineages that are globally distributed and active in hydrothermal systems. Here we apply multi-omics to characterize 51 metagenome-assembled genomes from three hydrothermal vent plumes in the Pacific and Atlantic Oceans that are affiliated with nine Proteobacteria lineages. Metabolic analyses revealed these organisms to contain a diverse functional repertoire including chemolithotrophic ability to utilize sulfur and C₁ compounds, and chemoorganotrophic ability to utilize environment-derived fatty acids, aromatics, carbohydrates, and peptides. Comparative genomics with marine and terrestrial microbiomes suggests that lineage-associated functional traits could explain niche specificity. Our results shed light on the ecological functions and metabolic strategies of novel Proteobacteria in hydrothermal systems and beyond, and highlight the relationship between genome diversification and environmental adaptation.

Temperature and Nutrient Levels Correspond with Lineage-Specific Microdiversification in the Ubiquitous and Abundant Freshwater Genus Limnohabitans

Most freshwater bacterial communities are characterized by a few dominant taxa that are often ubiquitous across freshwater biomes worldwide. Our understanding of the genomic diversity within these taxonomic groups is limited to a subset of taxa. Here, we investigated the genomic diversity that enables Limnohabitans, a freshwater genus key in funneling carbon from primary producers to higher trophic levels, to achieve abundance and ubiquity. We reconstructed eight putative Limnohabitans metagenome-assembled genomes (MAGs) from stations located along broad environmental gradients existing in Lake Michigan, part of Earth's largest surface freshwater system. De novo strain inference analysis resolved a total of 23 strains from these MAGs, which strongly partitioned into two habitat-specific clusters with cooccurring strains from different lineages. The largest number of strains belonged to the abundant LimB lineage, for which robust in situ strain delineation had not previously been achieved. Our data show that temperature and nutrient levels may be important environmental parameters associated with microdiversification within the Limnohabitans genus. In addition, strains predominant in low- and high-phosphorus conditions had larger genomic divergence than strains abundant under different temperatures. Comparative genomics and gene expression analysis yielded evidence for the ability of LimB populations to exhibit cellular motility and chemotaxis, a phenotype not yet associated with available Limnohabitans isolates. Our findings broaden historical marker gene-based surveys of Limnohabitans microdiversification and provide in situ evidence of genome diversity and its functional implications across freshwater gradients.IMPORTANCE Limnohabitans is an important bacterial taxonomic group for cycling carbon in freshwater ecosystems worldwide. Here, we examined the genomic diversity of different Limnohabitans lineages. We focused on the LimB lineage of this genus, which is globally distributed and often abundant, and its abundance has shown to be largely invariant to environmental change. Our data show that the LimB lineage is actually comprised of multiple cooccurring populations for which the composition and genomic characteristics are associated with variations in temperature and nutrient levels. The gene expression profiles of this lineage suggest the importance of chemotaxis and motility, traits that had not yet been associated with the Limnohabitans genus, in adapting to environmental conditions.

Gut microbial dysbiosis in individuals with Sjögren's syndrome

BACKGROUND

Autoimmune diseases have been associated with changes in the gut microbiome. In this study, the gut microbiome was evaluated in individuals with dry eye and bacterial compositions were correlated to dry eye (DE) measures. We prospectively included 13 individuals with who met full criteria for Sjögren's (SDE) and 8 individuals with features of Sjögren's but who did not meet full criteria (NDE) for a total of 21 cases as compared to 21 healthy controls. Stool was analyzed by 16S pyrosequencing, and associations between bacterial classes and DE symptoms and signs were examined.

RESULTS

Results showed that Firmicutes was the dominant phylum in the gut, comprising 40-60% of all phyla. On a phyla level, subjects with DE (SDE and NDE) had depletion of Firmicutes (1.1-fold) and an expansion of Proteobacteria (3.0-fold), Actinobacteria (1.7-fold), and Bacteroidetes (1.3-fold) compared to controls. Shannon's diversity index showed no differences between groups with respect to the numbers of different operational taxonomic units (OTUs) encountered (diversity) and the instances these unique OTUs were sampled (evenness). On the other hand, Faith's phylogenetic diversity showed increased diversity in cases vs controls, which reached significance when comparing SDE and controls (13.57 ± 0.89 and 10.96 ± 0.76, p = 0.02). Using Principle Co-ordinate Analysis, qualitative differences in microbial composition were noted with differential clustering of cases and controls. Dimensionality reduction and clustering of complex microbial data further showed differences between the three groups, with regard to microbial composition, association and clustering. Finally, differences in certain classes of bacteria were associated with DE symptoms and signs.

CONCLUSIONS

In conclusion, individuals with DE had gut microbiome alterations as compared to healthy controls. Certain classes of bacteria were associated with DE measures.

Computational profiling of the gut-brain axis: microflora dysbiosis insights to neurological disorders

Almost 2500 years after Hippocrates' observations on health and its direct association to the gastrointestinal tract, a paradigm shift has recently occurred, making the gut and its symbionts (bacteria, fungi, archaea and viruses) a point of convergence for studies. It is nowadays well established that the gut microflora's compositional diversity regulates via its genes (the microbiome) the host's health and provides preliminary insights into disease progression and regulation. The microbiome's involvement is evident in immunological and physiological studies that link changes in its biodiversity to its contributions to the host's phenotype but also in neurological investigations, substantiating the aptly named gut-brain axis. The definitive mechanisms of this last bidirectional interaction will be our main focus because it presents researchers with a new conundrum. In this review, we prospect current literature for computational analysis methodologies that accommodate the need for better understanding of the microbiome-gut-brain interactions and neurological disorder onset and progression, through cross-disciplinary systems biology applications. We will present bioinformatics tools used in exploring these synergies that help build and interpret microbial 16S ribosomal RNA data sets, produced by shotgun and high-throughput sequencing of healthy and neurological disorder samples stored in biological databases. These approaches provide alternative means for researchers to form hypotheses to their inquests faster, cheaper and swith precision. The goal of these studies relies on the integration of combined metagenomics and metabolomics assessments. An accurate characterization of the microbiome and its functionality can support new diagnostic, prognostic and therapeutic strategies for neurological disorders, customized for each individual host.

SCGid: a consensus approach to contig filtering and genome prediction from single-cell sequencing libraries of uncultured eukaryotes

MOTIVATION

Whole-genome sequencing of uncultured eukaryotic genomes is complicated by difficulties in acquiring sufficient amounts of tissue. Single-cell genomics (SCG) by multiple displacement amplification provides a technical workaround, yielding whole-genome libraries which can be assembled de novo. Downsides of multiple displacement amplification include coverage biases and exacerbation of contamination. These factors affect assembly continuity and fidelity, complicating discrimination of genomes from contamination and noise by available tools. Uncultured eukaryotes and their relatives are often underrepresented in large sequence data repositories, further impairing identification and separation.

RESULTS

We compare the ability of filtering approaches to remove contamination and resolve eukaryotic draft genomes from SCG metagenomes, finding significant variation in outcomes. To address these inconsistencies, we introduce a consensus approach that is codified in the SCGid software package. SCGid parallelly filters assemblies using different approaches, yielding three intermediate drafts from which consensus is drawn. Using genuine and mock SCG metagenomes, we show that our approach corrects for variation among draft genomes predicted by individual approaches and outperforms them in recapitulating published drafts in a fast and repeatable way, providing a useful alternative to available methods and manual curation.

AVAILABILITY AND IMPLEMENTATION

The SCGid package is implemented in python and R. Source code is available at http://www.github.com/amsesk/SCGid under the GNU GPL 3.0 license.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Accurate and complete genomes from metagenomes

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

Metabolic potentials of archaeal lineages resolved from metagenomes of deep Costa Rica sediments

Numerous archaeal lineages are known to inhabit marine subsurface sediments, although their distributions, metabolic capacities, and interspecies interactions are still not well understood. Abundant and diverse archaea were recently reported in Costa Rica (CR) margin subseafloor sediments recovered during IODP Expedition 334. Here, we recover metagenome-assembled genomes (MAGs) of archaea from the CR margin and compare them to their relatives from shallower settings. We describe 31 MAGs of six different archaeal lineages (Lokiarchaeota, Thorarchaeota, Heimdallarchaeota, Bathyarcheota, Thermoplasmatales, and Hadesarchaea) and thoroughly analyze representative MAGs from the phyla Lokiarchaeota and Bathyarchaeota. Our analysis suggests the potential capability of Lokiarchaeota members to anaerobically degrade aliphatic and aromatic hydrocarbons. We show it is genetically possible and energetically feasible for Lokiarchaeota to degrade benzoate if they associate with organisms using nitrate, nitrite, and sulfite as electron acceptors, which suggests a possibility of syntrophic relationships between Lokiarchaeota and nitrite and sulfite reducing bacteria. The novel Bathyarchaeota lineage possesses an incomplete methanogenesis pathway lacking the methyl coenzyme M reductase complex and encodes a noncanonical acetogenic pathway potentially coupling methylotrophy to acetogenesis via the methyl branch of Wood-Ljungdahl pathway. These metabolic characteristics suggest the potential of this Bathyarchaeota lineage to be a transition between methanogenic and acetogenic Bathyarchaeota lineages. This work expands our knowledge about the metabolic functional repertoire of marine benthic archaea.

Genome- and Community-Level Interaction Insights into Carbon Utilization and Element Cycling Functions of Hydrothermarchaeota in Hydrothermal Sediment

Hydrothermal vents release reduced compounds and small organic carbon compounds into the surrounding seawater, providing essential substrates for microbial growth and bioenergy transformations. Despite the wide distribution of the marine benthic group E archaea (referred to as Hydrothermarchaeota) in the hydrothermal environment, little is known about their genomic repertoires and biogeochemical significance. Here, we studied four highly complete (>80%) metagenome-assembled genomes (MAGs) from a black smoker chimney and the surrounding sulfur-rich sediments on the South Atlantic Mid-Ocean Ridge and publicly available data sets (the Integrated Microbial Genomes system of the U.S. Department of Energy-Joint Genome Institute and NCBI SRA data sets). Genomic analysis suggested a wide carbon metabolic diversity of Hydrothermarchaeota members, including the utilization of proteins, lactate, and acetate; the anaerobic degradation of aromatics; the oxidation of C₁ compounds (CO, formate, and formaldehyde); the utilization of methyl compounds; CO₂ incorporation by the tetrahydromethanopterin-based Wood-Ljungdahl pathway; and participation in the type III ribulose-1,5-bisphosphate carboxylase/oxygenase-based Calvin-Benson-Bassham cycle. These microbes also potentially oxidize sulfur, arsenic, and hydrogen and engage in anaerobic respiration based on sulfate reduction and denitrification. Among the 140 MAGs reconstructed from the black smoker chimney microbial community (including Hydrothermarchaeota MAGs), community-level metabolic predictions suggested a redundancy of carbon utilization and element cycling functions and interactive syntrophic and sequential utilization of substrates. These processes might make various carbon and energy sources widely accessible to the microorganisms. Further, the analysis suggested that Hydrothermarchaeota members contained important functional components obtained from the community via lateral gene transfer, becoming a distinctive clade. This might serve as a niche-adaptive strategy for metabolizing heavy metals, C₁ compounds, and reduced sulfur compounds. Collectively, the analysis provides comprehensive metabolic insights into the Hydrothermarchaeota IMPORTANCE This study provides comprehensive metabolic insights into the Hydrothermarchaeota from comparative genomics, evolution, and community-level perspectives. Members of the Hydrothermarchaeota synergistically participate in a wide range of carbon-utilizing and element cycling processes with other microorganisms in the community. We expand the current understanding of community interactions within the hydrothermal sediment and chimney, suggesting that microbial interactions based on sequential substrate metabolism are essential to nutrient and element cycling.

Genome-inferred spatio-temporal resolution of an uncultivated Roizmanbacterium reveals its ecological preferences in groundwater

Subsurface ecosystems like groundwater harbour diverse microbial communities, including small-sized, putatively symbiotic organisms of the Candidate Phyla Radiation, yet little is known about their ecological preferences and potential microbial partners. Here, we investigated a member of the superphylum Microgenomates (Cand. Roizmanbacterium ADI133) from oligotrophic groundwater using mini-metagenomics and monitored its spatio-temporal distribution using 16S rRNA gene analyses. A Roizmanbacteria-specific quantitative PCR assay allowed us to track its abundance over the course of 1 year within eight groundwater wells along a 5.4 km hillslope transect, where Roizmanbacteria reached maximum relative abundances of 2.3%. In-depth genomic analyses suggested that Cand. Roizmanbacterium ADI133 is a lactic acid fermenter, potentially able to utilize a range of complex carbon substrates, including cellulose. We hypothesize that it attaches to host cells using a trimeric autotransporter adhesin and inhibits their cell wall biosynthesis using a toxin-antitoxin system. Network analyses based on correlating Cand. Roizmanbacterium ADI133 abundances with amplicon sequencing-derived microbial community profiles suggested one potential host organism, classified as a member of the class Thermodesulfovibrionia (Nitrospirae). By providing lactate as an electron donor Cand. Roizmanbacterium ADI133 potentially mediates the transfer of carbon to other microorganisms and thereby is an important connector in the microbial community.

Light-driven anaerobic microbial oxidation of manganese

Oxygenic photosynthesis supplies organic carbon to the modern biosphere, but it is uncertain when this metabolism originated. It has previously been proposed^1,2 that photosynthetic reaction centres capable of splitting water arose by about 3 billion years ago on the basis of the inferred presence of manganese oxides in Archaean sedimentary rocks. However, this assumes that manganese oxides can be produced only in the presence of molecular oxygen³, reactive oxygen species^4,5 or by high-potential photosynthetic reaction centres^6,7. Here we show that communities of anoxygenic photosynthetic microorganisms biomineralize manganese oxides in the absence of molecular oxygen and high-potential photosynthetic reaction centres. Microbial oxidation of Mn(II) under strictly anaerobic conditions during the Archaean eon would have produced geochemical signals identical to those used to date the evolution of oxygenic photosynthesis before the Great Oxidation Event^1,2. This light-dependent process may also produce manganese oxides in the photic zones of modern anoxic water bodies and sediments.

MetaBMF: a scalable binning algorithm for large-scale reference-free metagenomic studies

MOTIVATION

Metagenomics studies microbial genomes in an ecosystem such as the gastrointestinal tract of a human. Identification of novel microbial species and quantification of their distributional variations among different samples that are sequenced using next-generation-sequencing technology hold the key to the success of most metagenomic studies. To achieve these goals, we propose a simple yet powerful metagenomic binning method, MetaBMF. The method does not require prior knowledge of reference genomes and produces highly accurate results, even at a strain level. Thus, it can be broadly used to identify disease-related microbial organisms that are not well-studied.

RESULTS

Mathematically, we count the number of mapped reads on each assembled genomic fragment cross different samples as our input matrix and propose a scalable stratified angle regression algorithm to factorize this count matrix into a product of a binary matrix and a nonnegative matrix. The binary matrix can be used to separate microbial species and the nonnegative matrix quantifies the species distributions in different samples. In simulation and empirical studies, we demonstrate that MetaBMF has a high binning accuracy. It can not only bin DNA fragments accurately at a species level but also at a strain level. As shown in our example, we can accurately identify the Shiga-toxigenic Escherichia coli O104: H4 strain which led to the 2011 German E.coli outbreak. Our efforts in these areas should lead to (i) fundamental advances in metagenomic binning, (ii) development and refinement of technology for the rapid identification and quantification of microbial distributions and (iii) finding of potential probiotics or reliable pathogenic bacterial strains.

AVAILABILITY AND IMPLEMENTATION

The software is available at https://github.com/didi10384/MetaBMF.

A review of methods and databases for metagenomic classification and assembly

Microbiome research has grown rapidly over the past decade, with a proliferation of new methods that seek to make sense of large, complex data sets. Here, we survey two of the primary types of methods for analyzing microbiome data: read classification and metagenomic assembly, and we review some of the challenges facing these methods. All of the methods rely on public genome databases, and we also discuss the content of these databases and how their quality has a direct impact on our ability to interpret a microbiome sample.

A review of methods and databases for metagenomic classification and assembly

Microbiome research has grown rapidly over the past decade, with a proliferation of new methods that seek to make sense of large, complex data sets. Here, we survey two of the primary types of methods for analyzing microbiome data: read classification and metagenomic assembly, and we review some of the challenges facing these methods. All of the methods rely on public genome databases, and we also discuss the content of these databases and how their quality has a direct impact on our ability to interpret a microbiome sample.

Taxonomic features and comparisons of the gut microbiome from two edible fungus-farming termites (Macrotermes falciger; M. natalensis) harvested in the Vhembe district of Limpopo, South Africa

Background

Termites are an important food resource for many human populations around the world, and are a good supply of nutrients. The fungus-farming ‘higher’ termite members of Macrotermitinae are also consumed by modern great apes and are implicated as critical dietary resources for early hominins. While the chemical nutritional composition of edible termites is well known, their microbiomes are unexplored in the context of human health. Here we sequenced the V4 region of the 16S rRNA gene of gut microbiota extracted from the whole intestinal tract of two Macrotermes sp. soldiers collected from the Limpopo region of South Africa.

Results

Major and minor soldier subcastes of M. falciger exhibit consistent differences in taxonomic representation, and are variable in microbial presence and abundance patterns when compared to another edible but less preferred species, M. natalensis. Subcaste differences include alternate patterns in sulfate-reducing bacteria and methanogenic Euryarchaeota abundance, and differences in abundance between Alistipes and Ruminococcaceae. M. falciger minor soldiers and M. natalensis soldiers have similar microbial profiles, likely from close proximity to the termite worker castes, particularly during foraging and fungus garden cultivation. Compared with previously published termite and cockroach gut microbiome data, the taxonomic representation was generally split between termites that directly digest lignocellulose and humic substrates and those that consume a more distilled form of nutrition as with the omnivorous cockroaches and fungus-farming termites. Lastly, to determine if edible termites may point to a shared reservoir for rare bacterial taxa found in the gut microbiome of humans, we focused on the genus Treponema. The majority of Treponema sequences from edible termite gut microbiota most closely relate to species recovered from other termites or from environmental samples, except for one novel OTU strain, which clustered separately with Treponema found in hunter-gatherer human groups.

Conclusions

Macrotermes consumed by humans display special gut microbial arrangements that are atypical for a lignocellulose digesting invertebrate, but are instead suited to the simplified nutrition in the fungus-farmer diet. Our work brings to light the particular termite microbiome features that should be explored further as avenues in human health, agricultural sustainability, and evolutionary research.

Electronic supplementary material

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

Contrasting Winter Versus Summer Microbial Communities and Metabolic Functions in a Permafrost Thaw Lake

Permafrost thawing results in the formation of thermokarst lakes, which are biogeochemical hotspots in northern landscapes and strong emitters of greenhouse gasses to the atmosphere. Most studies of thermokarst lakes have been in summer, despite the predominance of winter and ice-cover over much of the year, and the microbial ecology of these waters under ice remains poorly understood. Here we first compared the summer versus winter microbiomes of a subarctic thermokarst lake using DNA- and RNA-based 16S rRNA amplicon sequencing and qPCR. We then applied comparative metagenomics and used genomic bin reconstruction to compare the two seasons for changes in potential metabolic functions in the thermokarst lake microbiome. In summer, the microbial community was dominated by Actinobacteria and Betaproteobacteria, with phototrophic and aerobic pathways consistent with the utilization of labile and photodegraded substrates. The microbial community was strikingly different in winter, with dominance of methanogens, Planctomycetes, Chloroflexi and Deltaproteobacteria, along with various taxa of the Patescibacteria/Candidate Phyla Radiation (Parcubacteria, Microgenomates, Omnitrophica, Aminicenantes). The latter group was underestimated or absent in the amplicon survey, but accounted for about a third of the metagenomic reads. The winter lineages were associated with multiple reductive metabolic processes, fermentations and pathways for the mobilization and degradation of complex organic matter, along with a strong potential for syntrophy or cross-feeding. The results imply that the summer community represents a transient stage of the annual cycle, and that carbon dioxide and methane production continue through the prolonged season of ice cover via a taxonomically distinct winter community and diverse mechanisms of permafrost carbon transformation.

Genomic ecology of Marine Group II, the most common marine planktonic Archaea across the surface ocean

Planktonic Archaea have been detected in all the world's oceans and are found from surface waters to the deep sea. The two most common Archaea phyla are Thaumarchaeota and Euryarchaeota. Euryarchaeota are generally more common in surface waters, but very little is known about their ecology and their potential metabolisms. In this study, we explore the genomic ecology of the Marine Group II (MGII), the main marine planktonic Euryarchaeota, and test if it is composed of different ecologically relevant units. We re‐analyzed Tara Oceans metagenomes from the photic layer and the deep ocean by annotating sequences against a custom MGII database and by mapping gene co‐occurrences. Our data provide a global view of the distribution of Euryarchaeota, and more specifically of MGII subgroups, and reveal their association to a number of gene‐coding sequences. In particular, we show that MGII proteorhodopsins were detected in both the surface and the deep chlorophyll maximum layer and that different clusters of these light harvesting proteins were present. Our approach helped describing the set of genes found together with specific MGII subgroups. We could thus define genomic environments that could theoretically describe ecologically meaningful units and the ecological niche that they occupy.