Anvi'o: an advanced analysis and visualization platform for 'omics data

Pan-genome Analysis of Ancient and Modern Salmonella enterica Demonstrates Genomic Stability of the Invasive Para C Lineage for Millennia

Salmonella enterica serovar Paratyphi C causes enteric (paratyphoid) fever in humans. Its presentation can range from asymptomatic infections of the blood stream to gastrointestinal or urinary tract infection or even a fatal septicemia [1]. Paratyphi C is very rare in Europe and North America except for occasional travelers from South and East Asia or Africa, where the disease is more common [2, 3]. However, early 20^th-century observations in Eastern Europe [3, 4] suggest that Paratyphi C enteric fever may once have had a wide-ranging impact on human societies. Here, we describe a draft Paratyphi C genome (Ragna) recovered from the 800-year-old skeleton (SK152) of a young woman in Trondheim, Norway. Paratyphi C sequences were recovered from her teeth and bones, suggesting that she died of enteric fever and demonstrating that these bacteria have long caused invasive salmonellosis in Europeans. Comparative analyses against modern Salmonella genome sequences revealed that Paratyphi C is a clade within the Para C lineage, which also includes serovars Choleraesuis, Typhisuis, and Lomita. Although Paratyphi C only infects humans, Choleraesuis causes septicemia in pigs and boar [5] (and occasionally humans), and Typhisuis causes epidemic swine salmonellosis (chronic paratyphoid) in domestic pigs [2, 3]. These different host specificities likely evolved in Europe over the last ∼4,000 years since the time of their most recent common ancestor (tMRCA) and are possibly associated with the differential acquisitions of two genomic islands, SPI-6 and SPI-7. The tMRCAs of these bacterial clades coincide with the timing of pig domestication in Europe [6].

•

Salmonella enterica aDNA sequences were found within 800-year-old teeth and bone

•

The invasive Para C lineage was defined from 50,000 modern S. enterica genomes

•

The Para C lineage includes Ragna, the aDNA genome, and human and swine pathogens

•

Only few genomic changes occurred in the Para C lineage over its 3,000-year history

Metagenomics reveals functional synergy and novel polysaccharide utilization loci in the Castor canadensis fecal microbiome

The North American beaver (Castor canadensis) has long been considered an engineering marvel, transforming landscapes and shaping biological diversity through its dam building behavior. While the beaver possesses conspicuous morphological features uniquely adapted for the use of woody plants as construction materials and dietary staples, relatively little is known about the specialized microorganisms inhabiting the beaver gastrointestinal tract and their functional roles in determining host nutrition. Here we use a combination of shotgun metagenomics, functional screening and carbohydrate biochemistry to chart the community structure and metabolic power of the beaver fecal microbiome. We relate this information to the metabolic capacity of other wood feeding and hindgut fermenting organisms and profile the functional repertoire of glycoside hydrolase (GH) families distributed among and between population genome bins. Metagenomic screening revealed novel mechanisms of xylan oligomer degradation involving GH43 enzymes from uncharacterized subfamilies and divergent polysaccharide utilization loci, indicating the potential for synergistic biomass deconstruction. Together, these results open a functional metagenomic window on less conspicuous adaptations enabling the beaver microbiome to efficiently convert woody plants into host nutrition and point toward rational design of enhanced enzyme mixtures for biorefining process streams.

Transient Osmotic Perturbation Causes Long-Term Alteration to the Gut Microbiota

Osmotic diarrhea is a prevalent condition in humans caused by food intolerance, malabsorption, and widespread laxative use. Here, we assess the resilience of the gut ecosystem to osmotic perturbation at multiple length and timescales using mice as model hosts. Osmotic stress caused reproducible extinction of highly abundant taxa and expansion of less prevalent members in human and mouse microbiotas. Quantitative imaging revealed decimation of the mucus barrier during osmotic perturbation, followed by recovery. The immune system exhibited temporary changes in cytokine levels and a lasting IgG response against commensal bacteria. Increased osmolality prevented growth of commensal strains in vitro, revealing one mechanism contributing to extinction. Environmental availability of microbiota members mitigated extinction events, demonstrating how species reintroduction can affect community resilience. Our findings (1) demonstrate that even mild osmotic diarrhea can cause lasting changes to the microbiota and host and (2) lay the foundation for interventions that increase system-wide resilience.

Nitrogen-fixing populations of Planctomycetes and Proteobacteria are abundant in surface ocean metagenomes

Nitrogen fixation in the surface ocean impacts global marine nitrogen bioavailability and thus microbial primary productivity. Until now, cyanobacterial populations have been viewed as the main suppliers of bioavailable nitrogen in this habitat. Although PCR amplicon surveys targeting the nitrogenase reductase gene have revealed the existence of diverse non-cyanobacterial diazotrophic populations, subsequent quantitative PCR surveys suggest that they generally occur in low abundance. Here, we use state-of-the-art metagenomic assembly and binning strategies to recover nearly one thousand non-redundant microbial population genomes from the TARA Oceans metagenomes. Among these, we provide the first genomic evidence for non-cyanobacterial diazotrophs inhabiting surface waters of the open ocean, which correspond to lineages within the Proteobacteria and, most strikingly, the Planctomycetes. Members of the latter phylum are prevalent in aquatic systems, but have never been linked to nitrogen fixation previously. Moreover, using genome-wide quantitative read recruitment, we demonstrate that the discovered diazotrophs were not only widespread but also remarkably abundant (up to 0.3% of metagenomic reads for a single population) in both the Pacific Ocean and the Atlantic Ocean northwest. Our results extend decades of PCR-based gene surveys, and substantiate the importance of heterotrophic bacteria in the fixation of nitrogen in the surface ocean.

Gut microbes contribute to variation in solid organ transplant outcomes in mice

BACKGROUND

Solid organ transplant recipients show heterogeneity in the occurrence and timing of acute rejection episodes. Understanding the factors responsible for such variability in patient outcomes may lead to improved diagnostic and therapeutic approaches. Rejection kinetics of transplanted organs mainly depends on the extent of genetic disparities between donor and recipient, but a role for environmental factors is emerging. We have recently shown that major alterations of the microbiota following broad-spectrum antibiotics, or use of germ-free animals, promoted longer skin graft survival in mice. Here, we tested whether spontaneous differences in microbial colonization between genetically similar individuals can contribute to variability in graft rejection kinetics.

RESULTS

We compared rejection kinetics of minor mismatched skin grafts in C57BL/6 mice from Jackson Laboratory (Jax) and Taconic Farms (Tac), genetically similar animals colonized by different commensal microbes. Female Tac mice rejected skin grafts from vendor-matched males more quickly than Jax mice. We observed prolonged graft survival in Tac mice when they were exposed to Jax mice microbiome through co-housing or fecal microbiota transplantation (FMT) by gastric gavage. In contrast, exposure to Tac mice did not change graft rejection kinetics in Jax mice, suggesting a dominant suppressive effect of Jax microbiota. High-throughput sequencing of 16S rRNA gene amplicons from Jax and Tac mice fecal samples confirmed a convergence of microbiota composition after cohousing or fecal transfer. Our analysis of amplicon data associated members of a single bacterial genus, Alistipes, with prolonged graft survival. Consistent with this finding, members of the genus Alistipes were absent in a separate Tac cohort, in which fecal transfer from Jax mice failed to prolong graft survival.

CONCLUSIONS

These results demonstrate that differences in resident microbiome in healthy individuals may translate into distinct kinetics of graft rejection, and contribute to interpersonal variability in graft outcomes. The association between Alistipes and prolonged skin graft survival in mice suggests that members of this genus might affect host physiology, including at sites distal to the gastrointestinal tract. Overall, these findings allude to a potential therapeutic role for specific gut microbes to promote graft survival through the administration of probiotics, or FMT.

Pan-genomic approach shows insight of genetic divergence and pathogenic-adaptation of Pasteurella multocida

Pasteurella multocida is a gram-negative, non-motile bacterial pathogen, which is associated with chronic and acute infections as snuffles, pneumonia, atrophic rhinitis, fowl cholera and hemorrhagic septicemia. These diseases affect a wide range of domestic animals, leading to significant morbidity and mortality and causing significant economic losses worldwide. Due to the interest in deciphering the genetic diversity and process adaptive between P. multocida strains, this work aimed was to perform a pan-genome analysis to evidence horizontal gene transfer and positive selection among 23 P. multocida strains isolated from distinct diseases and hosts. The results revealed an open pan-genome containing 3585 genes and an accessory genome presenting 1200 genes. The phylogenomic analysis based on the presence/absence of genes and islands exhibit high levels of plasticity, which reflects a high intraspecific diversity and a possible adaptive mechanism responsible for the specific disease manifestation between the established groups (pneumonia, fowl cholera, hemorrhagic septicemia and snuffles). Additionally, we identified differences in accessory genes among groups, which are involved in sugar metabolism and transport systems, virulence-related genes and a high concentration of hypothetical proteins. However, there was no specific indispensable functional mechanism to decisively correlate the presence of genes and their adaptation to a specific host/disease. Also, positive selection was found only for two genes from sub-group hemorrhagic septicemia, serotype B. This comprehensive comparative genome analysis will provide new insights of horizontal gene transfers that play an essential role in the diversification and adaptation mechanism into P. multocida species to a specific disease.

SAMSA2: a standalone metatranscriptome analysis pipeline

Background

Complex microbial communities are an area of growing interest in biology. Metatranscriptomics allows researchers to quantify microbial gene expression in an environmental sample via high-throughput sequencing. Metatranscriptomic experiments are computationally intensive because the experiments generate a large volume of sequence data and each sequence must be compared with reference sequences from thousands of organisms.

Results

SAMSA2 is an upgrade to the original Simple Annotation of Metatranscriptomes by Sequence Analysis (SAMSA) pipeline that has been redesigned for standalone use on a supercomputing cluster. SAMSA2 is faster due to the use of the DIAMOND aligner, and more flexible and reproducible because it uses local databases. SAMSA2 is available with detailed documentation, and example input and output files along with examples of master scripts for full pipeline execution.

Conclusions

SAMSA2 is a rapid and efficient metatranscriptome pipeline for analyzing large RNA-seq datasets in a supercomputing cluster environment. SAMSA2 provides simplified output that can be examined directly or used for further analyses, and its reference databases may be upgraded, altered or customized to fit the needs of any experiment.

Profiling microbial strains in urban environments using metagenomic sequencing data

Background

The microbial communities populating human and natural environments have been extensively characterized with shotgun metagenomics, which provides an in-depth representation of the microbial diversity within a sample. Microbes thriving in urban environments may be crucially important for human health, but have received less attention than those of other environments. Ongoing efforts started to target urban microbiomes at a large scale, but the most recent computational methods to profile these metagenomes have never been applied in this context. It is thus currently unclear whether such methods, that have proven successful at distinguishing even closely related strains in human microbiomes, are also effective in urban settings for tasks such as cultivation-free pathogen detection and microbial surveillance. Here, we aimed at a) testing the currently available metagenomic profiling tools on urban metagenomics; b) characterizing the organisms in urban environment at the resolution of single strain and c) discussing the biological insights that can be inferred from such methods.

Results

We applied three complementary methods on the 1614 metagenomes of the CAMDA 2017 challenge. With MetaMLST we identified 121 known sequence-types from 15 species of clinical relevance. For instance, we identified several Acinetobacter strains that were close to the nosocomial opportunistic pathogen A. nosocomialis. With StrainPhlAn, a generalized version of the MetaMLST approach, we inferred the phylogenetic structure of Pseudomonas stutzeri strains and suggested that the strain-level heterogeneity in environmental samples is higher than in the human microbiome. Finally, we also probed the functional potential of the different strains with PanPhlAn. We further showed that SNV-based and pangenome-based profiling provide complementary information that can be combined to investigate the evolutionary trajectories of microbes and to identify specific genetic determinants of virulence and antibiotic resistances within closely related strains.

Conclusion

We show that strain-level methods developed primarily for the analysis of human microbiomes can be effective for city-associated microbiomes. In fact, (opportunistic) pathogens can be tracked and monitored across many hundreds of urban metagenomes. However, while more effort is needed to profile strains of currently uncharacterized species, this work poses the basis for high-resolution analyses of microbiomes sampled in city and mass transportation environments.

Reviewers

This article was reviewed by Alexandra Bettina Graf, Daniel Huson and Trevor Cickovski.

Electronic supplementary material

The online version of this article (10.1186/s13062-018-0211-z) contains supplementary material, which is available to authorized users.

Microbiome Dynamics in a Large Artificial Seawater Aquarium

Artificial habitats for animals have high commercial and societal value. Microbial communities (microbiomes) in such habitats may play ecological roles similar to those in nature. However, this hypothesis remains largely untested. Georgia Aquarium's Ocean Voyager (OV) exhibit is a closed-system aquatic habitat that mimics the oligotrophic open ocean and houses thousands of large marine animals, including fish, sea turtles, and whale sharks. We present a 14-month time series characterizing the OV water column microbiome. The composition and stability of the microbiome differed from those of natural marine environments with similar chemical features. The composition shifted dramatically over the span of 2 weeks and was characterized by bloom events featuring members of two heterotrophic bacterial lineages with cosmopolitan distributions in the oceans. The relative abundances of these lineages were inversely correlated, suggesting an overlap in ecological niches. Transcript mapping to metagenome-assembled genomes (MAGs) of these taxa identified unique characteristics, including the presence and activity of genes for the synthesis and degradation of cyanophycin, an amino acid polymer linked to environmental stress and found frequently in cyanobacteria but rarely in heterotrophic bacteria. The dominant MAGs also contained and transcribed plasmid-associated sequences, suggesting a role for conjugation in adaptation to the OV environment. These findings indicate a highly dynamic microbiome despite the stability of the physical and chemical parameters of the water column. Characterizing how such fluctuations affect microbial function may inform our understanding of animal health in closed aquaculture systems.

IMPORTANCE Public aquariums play important societal roles, for example, by promoting science education and helping conserve biodiversity. The health of aquarium animals depends on interactions with the surrounding microbiome. However, the extent to which aquariums recreate a stable and natural microbial ecosystem is uncertain. This study describes the taxonomic composition of the water column microbiome over 14 months in a large indoor aquatic habitat, the Ocean Voyager exhibit at the Georgia Aquarium. Despite stable water column conditions, the exhibit experienced blooms in which the abundance of a single bacterial strain increased to over 65% of the community. Genome analysis indicated that the OV's dominant strains share unique adaptations, notably genes for storage polymers associated with environmental stress. These results, interpreted alongside data from natural ocean systems and another artificial seawater aquarium, suggest a highly dynamic aquarium microbiome and raise questions of how microbiome stability may affect the ecological health of the habitat.

Peripartum Antibiotics Promote Gut Dysbiosis, Loss of Immune Tolerance, and Inflammatory Bowel Disease in Genetically Prone Offspring

Factors affecting the developing neonatal gut microbiome and immune networks may increase the risk of developing complex immune disorders such as inflammatory bowel diseases (IBD). In particular, peripartum antibiotics have been suggested as risk factors for human IBD, although direct evidence is lacking. Therefore, we examined the temporal impact of the commonly used antibiotic cefoperazone on both maternal and offspring microbiota when administered to dams during the peripartum period in the IL-10-deficient murine colitis model. By rigorously controlling for cage, gender, generational, and murine pathobiont confounders, we observed that offspring from cefoperazone-exposed dams develop a persistent gut dysbiosis into adulthood associated with skewing of the host immune system and increased susceptibility to spontaneous and chemically dextran sodium sulfate (DSS)-induced colitis. Thus, early life exposure to antibiotic-induced maternal dysbiosis during a critical developmental window for gut microbial assemblage and immune programming elicits a lasting impact of increased IBD risk on genetically susceptible offspring.

Microbial community and metagenome dynamics during biodegradation of dispersed oil reveals potential key-players in cold Norwegian seawater

Oil biodegradation as a weathering process has been extensively investigated over the years, especially after the Deepwater Horizon blowout. In this study, we performed microcosm experiments at 5 °C with chemically dispersed oil in non-amended seawater. We link biodegradation processes with microbial community and metagenome dynamics and explain the succession based on substrate specialization. Reconstructed genomes and 16S rRNA gene analysis revealed that Bermanella and Zhongshania were the main contributors to initial n-alkane breakdown, while subsequent abundances of Colwellia and microorganisms closely related to Porticoccaceae were involved in secondary n‑alkane breakdown and beta‑oxidation. Cycloclasticus, Porticoccaceae and Spongiiabcteraceae were associated with degradation of mono- and poly-cyclic aromatics. Successional pattern of genes coding for hydrocarbon degrading enzymes at metagenome level, and reconstructed genomic content, revealed a high differentiation of bacteria involved in hydrocarbon biodegradation. A cooperation among oil degrading microorganisms is thus needed for the complete substrate transformation.

On the Road to Strain-Resolved Comparative Metagenomics

Metagenomics has transformed microbiology, but its potential has not been fully expressed yet. From computational methods for digging deeper into metagenomes to study designs for addressing specific hypotheses, the Segata Lab is pursuing an integrative metagenomic approach to describe and model human-associated microbial communities as collections of strains.

Metagenomics has transformed microbiology, but its potential has not been fully expressed yet. From computational methods for digging deeper into metagenomes to study designs for addressing specific hypotheses, the Segata Lab is pursuing an integrative metagenomic approach to describe and model human-associated microbial communities as collections of strains. Linking strain variants to host phenotypes and performing cultivation-free population genomics require large cohorts and meta-analysis strategies to synthesize available cohorts but can revolutionize our understanding of the personalized host-microbiome interface which is at the base of human health.

Comparative genomics sheds light on niche differentiation and the evolutionary history of comammox Nitrospira

The description of comammox Nitrospira spp., performing complete ammonia-to-nitrate oxidation, and their co-occurrence with canonical β-proteobacterial ammonia oxidizing bacteria (β-AOB) in the environment, calls into question the metabolic potential of comammox Nitrospira and the evolutionary history of their ammonia oxidation pathway. We report four new comammox Nitrospira genomes, constituting two novel species, and the first comparative genomic analysis on comammox Nitrospira. Unlike canonical Nitrospira, comammox Nitrospira genomes lack genes for assimilatory nitrite reduction, suggesting that they have lost the potential to use external nitrite nitrogen sources. By contrast, compared to canonical Nitrospira, comammox Nitrospira harbor a higher diversity of urea transporters and copper homeostasis genes and lack cyanate hydratase genes. Additionally, the two comammox clades differ in their ammonium uptake systems. Contrary to β-AOB, comammox Nitrospira genomes have single copies of the two central ammonia oxidation pathway operons. Similar to ammonia oxidizing archaea and some oligotrophic AOB strains, they lack genes involved in nitric oxide reduction. Furthermore, comammox Nitrospira genomes encode genes that might allow efficient growth at low oxygen concentrations. Regarding the evolutionary history of comammox Nitrospira, our analyses indicate that several genes belonging to the ammonia oxidation pathway could have been laterally transferred from β-AOB to comammox Nitrospira. We postulate that the absence of comammox genes in other sublineage II Nitrospira genomes is the result of subsequent loss.

Herbivorous turtle ants obtain essential nutrients from a conserved nitrogen-recycling gut microbiome

Nitrogen acquisition is a major challenge for herbivorous animals, and the repeated origins of herbivory across the ants have raised expectations that nutritional symbionts have shaped their diversification. Direct evidence for N provisioning by internally housed symbionts is rare in animals; among the ants, it has been documented for just one lineage. In this study we dissect functional contributions by bacteria from a conserved, multi-partite gut symbiosis in herbivorous Cephalotes ants through in vivo experiments, metagenomics, and in vitro assays. Gut bacteria recycle urea, and likely uric acid, using recycled N to synthesize essential amino acids that are acquired by hosts in substantial quantities. Specialized core symbionts of 17 studied Cephalotes species encode the pathways directing these activities, and several recycle N in vitro. These findings point to a highly efficient N economy, and a nutritional mutualism preserved for millions of years through the derived behaviors and gut anatomy of Cephalotes ants.

A software tool 'CroCo' detects pervasive cross-species contamination in next generation sequencing data

Background

Multiple RNA samples are frequently processed together and often mixed before multiplex sequencing in the same sequencing run. While different samples can be separated post sequencing using sample barcodes, the possibility of cross contamination between biological samples from different species that have been processed or sequenced in parallel has the potential to be extremely deleterious for downstream analyses.

Results

We present CroCo, a software package for identifying and removing such cross contaminants from assembled transcriptomes. Using multiple, recently published sequence datasets, we show that cross contamination is consistently present at varying levels in real data. Using real and simulated data, we demonstrate that CroCo detects contaminants efficiently and correctly. Using a real example from a molecular phylogenetic dataset, we show that contaminants, if not eliminated, can have a decisive, deleterious impact on downstream comparative analyses.

Conclusions

Cross contamination is pervasive in new and published datasets and, if undetected, can have serious deleterious effects on downstream analyses. CroCo is a database-independent, multi-platform tool, designed for ease of use, that efficiently and accurately detects and removes cross contamination in assembled transcriptomes to avoid these problems. We suggest that the use of CroCo should become a standard cleaning step when processing multiple samples for transcriptome sequencing.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0486-7) contains supplementary material, which is available to authorized users.

Genome-resolved metagenomics identifies genetic mobility, metabolic interactions, and unexpected diversity in perchlorate-reducing communities

Dissimilatory perchlorate reduction is an anaerobic respiratory pathway that in communities might be influenced by metabolic interactions. Because the genes for perchlorate reduction are horizontally transferred, previous studies have been unable to identify uncultivated perchlorate-reducing populations. Here we recovered metagenome-assembled genomes from perchlorate-reducing sediment enrichments and employed a manual scaffolding approach to reconstruct gene clusters for perchlorate reduction found within mobile genetic elements. De novo assembly and binning of four enriched communities yielded 48 total draft genomes. In addition to canonical perchlorate reduction gene clusters and taxa, a new type of gene cluster with an alternative perchlorate reductase was identified. Phylogenetic analysis indicated past exchange between these gene clusters, and the presence of plasmids with either gene cluster shows that the potential for gene transfer via plasmid persisted throughout enrichment. However, a majority of genomes in each community lacked perchlorate reduction genes. Putative chlorate-reducing or sulfur-reducing populations were dominant in most communities, supporting the hypothesis that metabolic interactions might result from perchlorate reduction intermediates and byproducts. Other populations included a novel phylum-level lineage (Ca. Muirbacteria) and epibiotic prokaryotes with no known role in perchlorate reduction. These results reveal unexpected genetic diversity, suggest that perchlorate-reducing communities involve substantial metabolic interactions, and encourage expanded strategies to further understand the evolution and ecology of this metabolism.

Genome-resolved metagenomics of sugarcane vinasse bacteria

BACKGROUND

The production of 1 L of ethanol from sugarcane generates up to 12 L of vinasse, which is a liquid waste containing an as-yet uncharacterized microbial assemblage. Most vinasse is destined for use as a fertilizer on the sugarcane fields because of the high organic and K content; however, increased N₂O emissions have been observed when vinasse is co-applied with inorganic N fertilizers. Here we aimed to characterize the microbial assemblage of vinasse to determine the gene potential of vinasse microbes for contributing to negative environmental effects during fertirrigation and/or to the obstruction of bioethanol fermentation.

RESULTS

We measured chemical characteristics and extracted total DNA from six vinasse batches taken over 1.5 years from a bioethanol and sugar mill in Sao Paulo State. The vinasse microbial assemblage was characterized by low alpha diversity with 5-15 species across the six vinasses. The core genus was Lactobacillus. The top six represented bacterial genera across the samples were Lactobacillus, Megasphaera and Mitsuokella (Phylum Firmicutes, 35-97% of sample reads); Arcobacter and Alcaligenes (Phylum Proteobacteria, 0-40%); Dysgonomonas (Phylum Bacteroidetes, 0-53%); and Bifidobacterium (Phylum Actinobacteria, 0-18%). Potential genes for denitrification but not nitrification were identified in the vinasse metagenomes, with putative nirK and nosZ genes the most represented. Binning resulted in 38 large bins with between 36.0 and 99.3% completeness, and five small mobile element bins. Of the large bins, 53% could be classified at the phylum level as Firmicutes, 15% as Proteobacteria, 13% as unknown phyla, 13% as Bacteroidetes and 6% as Actinobacteria. The large bins spanned a range of potential denitrifiers; moreover, the genetic repertoires of all the large bins included the presence of genes involved in acetate, CO₂, ethanol, H₂O₂, and lactose metabolism; for many of the large bins, genes related to the metabolism of mannitol, xylose, butyric acid, cellulose, sucrose, "3-hydroxy" fatty acids and antibiotic resistance were present based on the annotations. In total, 21 vinasse bacterial draft genomes were submitted to the genome repository.

CONCLUSIONS

Identification of the gene repertoires of vinasse bacteria and assemblages supported the idea that organic carbon and nitrogen present in vinasse together with microbiological variation of vinasse might lead to varying patterns of N₂O emissions during fertirrigation. Furthermore, we uncovered draft genomes of novel strains of known bioethanol contaminants, as well as draft genomes unknown at the phylum level. This study will aid efforts to improve bioethanol production efficiency and sugarcane agriculture sustainability.

Comparative genomic inference suggests mixotrophic lifestyle for Thorarchaeota

Thorarchaeota are a new archaeal phylum within the Asgard superphylum, whose ancestors have been proposed to play possible ecological roles in cellular evolution. However, little is known about the lifestyles of these uncultured archaea. To provide a better resolution of the ecological roles and metabolic capacity of Thorarchaeota, we obtained Thorarchaeota genomes reconstructed from metagenomes of different depth layers in mangrove and mudflat sediments. These genomes from deep anoxic layers suggest the presence of Thorarchaeota with the potential to degrade organic matter, fix inorganic carbon, reduce sulfur/sulfate and produce acetate. In particular, Thorarchaeota may be involved in ethanol production, nitrogen fixation, nitrite reduction, and arsenic detoxification. Interestingly, these Thorarchaeotal genomes are inferred to contain the tetrahydromethanopterin and tetrahydrofolate Wood-Ljungdahl (WL) pathways for CO₂ reduction, and the latter WL pathway appears to have originated from bacteria. These archaea are predicted to be able to use various inorganic and organic carbon sources, possessing genes inferred to encode ribulose bisphosphate carboxylase-like proteins (normally without RuBisCO activity) and a near-complete Calvin-Benson-Bassham cycle. The existence of eukaryotic selenocysteine insertion sequences and many genes for proteins previously considered eukaryote-specific in Thorarchaeota genomes provide new insights into their evolutionary roles in the origin of eukaryotic cellular complexity. Resolving the metabolic capacities of these enigmatic archaea and their origins will enhance our understanding of the origins of eukaryotes and their roles in ecosystems.

Systems biology of the human microbiome

Recent research has shown that the microbiome-a collection of microorganisms, including bacteria, fungi, and viruses, living on and in a host-are of extraordinary importance in human health, even from conception and development in the uterus. Therefore, to further our ability to diagnose disease, to predict treatment outcomes, and to identify novel therapeutics, it is essential to include microbiome and microbial metabolic biomarkers in Systems Biology investigations. In clinical studies or, more precisely, Systems Medicine approaches, we can use the diversity and individual characteristics of the personal microbiome to enhance our resolution for patient stratification. In this review, we explore several Systems Medicine approaches, including Microbiome Wide Association Studies to understand the role of the human microbiome in health and disease, with a focus on 'preventive medicine' or P4 (i.e., personalized, predictive, preventive, participatory) medicine.

Diversity of Rare and Abundant Prokaryotic Phylotypes in the Prony Hydrothermal Field and Comparison with Other Serpentinite-Hosted Ecosystems

The Bay of Prony, South of New Caledonia, represents a unique serpentinite-hosted hydrothermal field due to its coastal situation. It harbors both submarine and intertidal active sites, discharging hydrogen- and methane-rich alkaline fluids of low salinity and mild temperature through porous carbonate edifices. In this study, we have extensively investigated the bacterial and archaeal communities inhabiting the hydrothermal chimneys from one intertidal and three submarine sites by 16S rRNA gene amplicon sequencing. We show that the bacterial community of the intertidal site is clearly distinct from that of the submarine sites with species distribution patterns driven by only a few abundant populations, affiliated to the Chloroflexi and Proteobacteria phyla. In contrast, the distribution of archaeal taxa seems less site-dependent, as exemplified by the co-occurrence, in both submarine and intertidal sites, of two dominant phylotypes of Methanosarcinales previously thought to be restricted to serpentinizing systems, either marine (Lost City Hydrothermal Field) or terrestrial (The Cedars ultrabasic springs). Over 70% of the phylotypes were rare and included, among others, all those affiliated to candidate divisions. We finally compared the distribution of bacterial and archaeal phylotypes of Prony Hydrothermal Field with those of five previously studied serpentinizing systems of geographically distant sites. Although sensu stricto no core microbial community was identified, a few uncultivated lineages, notably within the archaeal order Methanosarcinales and the bacterial class Dehalococcoidia (the candidate division MSBL5) were exclusively found in a few serpentinizing systems while other operational taxonomic units belonging to the orders Clostridiales, Thermoanaerobacterales, or the genus Hydrogenophaga, were abundantly distributed in several sites. These lineages may represent taxonomic signatures of serpentinizing ecosystems. These findings extend our current knowledge of the microbial diversity inhabiting serpentinizing systems and their biogeography.

Draft Genome Sequences of Novel Pseudomonas, Flavobacterium, and Sediminibacterium [corrected] Strains from a Freshwater Ecosystem

Freshwater ecosystems represent 0.01% of the water on Earth, but they support 6% of global biodiversity that is still mostly uncharacterized. Here, we describe the genome sequences of three strains belonging to novel species in the Pseudomonas, Flavobacterium, and Sediminibacterium genera recovered from a water sample of Lake Garda, Italy.

Metagenomic investigation of vestimentiferan tubeworm endosymbionts from Mid-Cayman Rise reveals new insights into metabolism and diversity

BACKGROUND

The microbial endosymbionts of two species of vestimentiferan tubeworms (Escarpia sp. and Lamellibrachia sp.2) collected from an area of low-temperature hydrothermal diffuse vent flow at the Mid-Cayman Rise (MCR) in the Caribbean Sea were characterized using microscopy, phylogenetic analyses, and a metagenomic approach.

RESULTS

Bacteria, with a typical Gram negative cell envelope contained within membrane-bound vacuoles, were observed within the trophosome of both tubeworm species. Phylogenetic analysis of the 16S rRNA gene and ITS region suggested MCR individuals harbored highly similar endosymbionts that were > 98% identical, with the exception of two symbionts that showed a 60 bp insertion within the ITS region. All sequences from MCR endosymbionts formed a separate well-supported clade that diverged from those of symbionts of seep and vent vestimentiferans from the Pacific, Gulf of Mexico, and Mediterranean Sea. The metagenomes of the symbionts of two specimens of each tubeworm species were sequenced, and two distinct Gammaproteobacteria metagenome-assembled genomes (MAGs) of more than 4 Mbp assembled. An Average Nucleotide Identity (ANI) of 86.5% between these MAGs, together with distinct 16S rRNA gene and ITS sequences, indicate the presence of multiple endosymbiont phylotypes at the MCR, with one MAG shared between one Escarpia and two Lamellibrachia individuals, indicating these endosymbionts are not specific to either host species. Genes for sulfur and hydrogen oxidation, nitrate reduction (assimilatory and dissimilatory), glycolysis and the Krebs cycle, peptide, sugar, and lipid transporters, and both rTCA and CBB carbon fixation cycles were detected in the MAGs, highlighting key and shared functions with symbiont metagenomes of the vestimentiferans Riftia, Tevnia, and Ridgeia from the Pacific. The potential for a second hydrogen oxidation pathway (via a bidirectional hydrogenase), formate dehydrogenase, a catalase, and several additional peptide transporters were found exclusively in the MCR endosymbiont MAGs.

CONCLUSIONS

The present study adds new evidence that tubeworm endosymbionts can potentially switch from autotrophic to heterotrophic metabolism, or may be mixotrophic, presumably while free-living, and also suggests their versatile metabolic potential may enable both the host and symbionts to exploit a wide range of environmental conditions. Together, the marked gene content and sequence dissimilarity at the rRNA operon and whole genome level between vent and seep symbionts suggest these newly described endosymbionts from the MCR belong to a novel tubeworm endosymbiont genera, introduced as Candidatus Vondammii.

Linking pangenomes and metagenomes: the Prochlorococcus metapangenome

Pangenomes offer detailed characterizations of core and accessory genes found in a set of closely related microbial genomes, generally by clustering genes based on sequence homology. In comparison, metagenomes facilitate highly resolved investigations of the relative distribution of microbial genomes and individual genes across environments through read recruitment analyses. Combining these complementary approaches can yield unique insights into the functional basis of microbial niche partitioning and fitness, however, advanced software solutions are lacking. Here we present an integrated analysis and visualization strategy that provides an interactive and reproducible framework to generate pangenomes and to study them in conjunction with metagenomes. To investigate its utility, we applied this strategy to a Prochlorococcus pangenome in the context of a large-scale marine metagenomic survey. The resulting Prochlorococcus metapangenome revealed remarkable differential abundance patterns between very closely related isolates that belonged to the same phylogenetic cluster and that differed by only a small number of gene clusters in the pangenome. While the relationships between these genomes based on gene clusters correlated with their environmental distribution patterns, phylogenetic analyses using marker genes or concatenated single-copy core genes did not recapitulate these patterns. The metapangenome also revealed a small set of core genes that mostly occurred in hypervariable genomic islands of the Prochlorococcus populations, which systematically lacked read recruitment from surface ocean metagenomes. Notably, these core gene clusters were all linked to sugar metabolism, suggesting potential benefits to Prochlorococcus from a high sequence diversity of sugar metabolism genes. The rapidly growing number of microbial genomes and increasing availability of environmental metagenomes provide new opportunities to investigate the functioning and the ecology of microbial populations, and metapangenomes can provide unique insights for any taxon and biome for which genomic and sufficiently deep metagenomic data are available.

The reconstruction of 2,631 draft metagenome-assembled genomes from the global oceans

Microorganisms play a crucial role in mediating global biogeochemical cycles in the marine environment. By reconstructing the genomes of environmental organisms through metagenomics, researchers are able to study the metabolic potential of Bacteria and Archaea that are resistant to isolation in the laboratory. Utilizing the large metagenomic dataset generated from 234 samples collected during the Tara Oceans circumnavigation expedition, we were able to assemble 102 billion paired-end reads into 562 million contigs, which in turn were co-assembled and consolidated in to 7.2 million contigs ≥2 kb in length. Approximately 1 million of these contigs were binned to reconstruct draft genomes. In total, 2,631 draft genomes with an estimated completion of ≥50% were generated (1,491 draft genomes >70% complete; 603 genomes >90% complete). A majority of the draft genomes were manually assigned phylogeny based on sets of concatenated phylogenetic marker genes and/or 16S rRNA gene sequences. The draft genomes are now publically available for the research community at-large.

Draft Genome Sequence of the Predatory Marine Bacterium Halobacteriovorax sp. Strain JY17

A draft genome sequence of Halobacteriovorax sp. strain JY17 was assembled from a metagenomic data set. The 3.47-Mbp genome of this unusual predatory bacterium contains 3,263 protein-coding sequences, 33 tRNAs, and 2 copies each of the 16S, 23S, and 5S rRNA genes. This is only the third sequenced representative of this genus.

Metagenomics-guided analysis of microbial chemolithoautotrophic phosphite oxidation yields evidence of a seventh natural CO2 fixation pathway

Dissimilatory phosphite oxidation (DPO), a microbial metabolism by which phosphite (HPO32-) is oxidized to phosphate (PO43-), is the most energetically favorable chemotrophic electron-donating process known. Only one DPO organism has been described to date, and little is known about the environmental relevance of this metabolism. In this study, we used 16S rRNA gene community analysis and genome-resolved metagenomics to characterize anaerobic wastewater treatment sludge enrichments performing DPO coupled to CO2 reduction. We identified an uncultivated DPO bacterium, Candidatus Phosphitivorax (Ca. P.) anaerolimi strain Phox-21, that belongs to candidate order GW-28 within the Deltaproteobacteria, which has no known cultured isolates. Genes for phosphite oxidation and for CO2 reduction to formate were found in the genome of Ca. P. anaerolimi, but it appears to lack any of the known natural carbon fixation pathways. These observations led us to propose a metabolic model for autotrophic growth by Ca. P. anaerolimi whereby DPO drives CO2 reduction to formate, which is then assimilated into biomass via the reductive glycine pathway.

Reef-Building Corals as a Tool for Climate Change Research in the Genomics Era

Coral reef ecosystems are among the most biodiverse habitats in the marine realm. They not only contribute with a plethora of ecosystem services, but they also are beneficial to humankind via nurturing marine fisheries and sustaining recreational activities. We will discuss the biology of coral reefs and their ecophysiology including the complex bacterial microbiota associated with them.

Decontaminating eukaryotic genome assemblies with machine learning

BACKGROUND

High-throughput sequencing has made it theoretically possible to obtain high-quality de novo assembled genome sequences but in practice DNA extracts are often contaminated with sequences from other organisms. Currently, there are few existing methods for rigorously decontaminating eukaryotic assemblies. Those that do exist filter sequences based on nucleotide similarity to contaminants and risk eliminating sequences from the target organism.

RESULTS

We introduce a novel application of an established machine learning method, a decision tree, that can rigorously classify sequences. The major strength of the decision tree is that it can take any measured feature as input and does not require a priori identification of significant descriptors. We use the decision tree to classify de novo assembled sequences and compare the method to published protocols.

CONCLUSIONS

A decision tree performs better than existing methods when classifying sequences in eukaryotic de novo assemblies. It is efficient, readily implemented, and accurately identifies target and contaminant sequences. Importantly, a decision tree can be used to classify sequences according to measured descriptors and has potentially many uses in distilling biological datasets.

Bacteriophage Distributions and Temporal Variability in the Ocean's Interior

Bacteriophages are numerically the most abundant DNA-containing entities in the oligotrophic ocean, yet how specific phage populations vary over time and space remains to be fully explored. Here, we conducted a metagenomic time-series survey of double-stranded DNA phages throughout the water column in the North Pacific Subtropical Gyre, encompassing 1.5 years from depths of 25 to 1,000 m. Viral gene sequences were identified in assembled metagenomic samples, yielding an estimated 172,385 different viral gene families. Viral marker gene distributions suggested that lysogeny was more prevalent at mesopelagic depths than in surface waters, consistent with prior prophage induction studies using mitomycin C. A total of 129 ALOHA viral genomes and genome fragments from 20 to 108 kbp were selected for further study, which represented the most abundant phages in the water column. Phage genotypes displayed discrete population structures. Most phages persisted throughout the time-series and displayed a strong depth structure that mirrored the stratified depth distributions of co-occurring bacterial taxa in the water column. Mesopelagic phages were distinct from surface water phages with respect to diversity, gene content, putative life histories, and temporal persistence, reflecting depth-dependent differences in host genomic architectures and phage reproductive strategies. The spatiotemporal distributions of the most abundant open-ocean bacteriophages that we report here provide new insight into viral temporal persistence, life history, and virus-host-environment interactions throughout the open-ocean water column.

The North Pacific Subtropical Gyre represents one of the largest biomes on the planet, where microbial communities are central mediators of ecosystem dynamics and global biogeochemical cycles. Critical members of these communities are the viruses of marine bacteria, which can alter microbial metabolism and significantly influence their survival and productivity. To better understand these viral assemblages, we conducted genomic analyses of planktonic viruses over a seasonal cycle to ocean depths of 1,000 m. We identified 172,385 different viral gene families and 129 unique virus genotypes in this open-ocean setting. The spatiotemporal distributions of the most abundant open-ocean viruses that we report here provide new insights into viral temporal variability, life history, and virus-host-environment interactions throughout the water column.

Draft Genome Sequence of the Planktic Cyanobacterium Tychonema bourrellyi, Isolated from Alpine Lentic Freshwater

We describe here the draft genome sequence of the cyanobacterium Tychonema bourrellyi, assembled from a metagenome of a nonaxenic culture. The strain (FEM_GT703) was isolated from a freshwater sample taken from Lake Garda, Italy. The draft genome sequence represents the first assembled T. bourrellyi strain.

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.

Comparative Pan-Genome Analysis of Piscirickettsia salmonis Reveals Genomic Divergences within Genogroups

Piscirickettsia salmonis is the etiological agent of salmonid rickettsial septicemia, a disease that seriously affects the salmonid industry. Despite efforts to genomically characterize P. salmonis, functional information on the life cycle, pathogenesis mechanisms, diagnosis, treatment, and control of this fish pathogen remain lacking. To address this knowledge gap, the present study conducted an in silico pan-genome analysis of 19 P. salmonis strains from distinct geographic locations and genogroups. Results revealed an expected open pan-genome of 3,463 genes and a core-genome of 1,732 genes. Two marked genogroups were identified, as confirmed by phylogenetic and phylogenomic relationships to the LF-89 and EM-90 reference strains, as well as by assessments of genomic structures. Different structural configurations were found for the six identified copies of the ribosomal operon in the P. salmonis genome, indicating translocation throughout the genetic material. Chromosomal divergences in genomic localization and quantity of genetic cassettes were also found for the Dot/Icm type IVB secretion system. To determine divergences between core-genomes, additional pan-genome descriptions were compiled for the so-termed LF and EM genogroups. Open pan-genomes composed of 2,924 and 2,778 genes and core-genomes composed of 2,170 and 2,228 genes were respectively found for the LF and EM genogroups. The core-genomes were functionally annotated using the Gene Ontology, KEGG, and Virulence Factor databases, revealing the presence of several shared groups of genes related to basic function of intracellular survival and bacterial pathogenesis. Additionally, the specific pan-genomes for the LF and EM genogroups were defined, resulting in the identification of 148 and 273 exclusive proteins, respectively. Notably, specific virulence factors linked to adherence, colonization, invasion factors, and endotoxins were established. The obtained data suggest that these genes could be directly associated with inter-genogroup differences in pathogenesis and host-pathogen interactions, information that could be useful in designing novel strategies for diagnosing and controlling P. salmonis infection.

Genomic variation in microbial populations inhabiting the marine subseafloor at deep-sea hydrothermal vents

Little is known about evolutionary drivers of microbial populations in the warm subseafloor of deep-sea hydrothermal vents. Here we reconstruct 73 metagenome-assembled genomes (MAGs) from two geochemically distinct vent fields in the Mid-Cayman Rise to investigate patterns of genomic variation within subseafloor populations. Low-abundance populations with high intra-population diversity coexist alongside high-abundance populations with low genomic diversity, with taxonomic differences in patterns of genomic variation between the mafic Piccard and ultramafic Von Damm vent fields. Populations from Piccard are significantly enriched in nonsynonymous mutations, suggesting stronger purifying selection in Von Damm relative to Piccard. Comparison of nine Sulfurovum MAGs reveals two high-coverage, low-diversity MAGs from Piccard enriched in unique genes related to the cellular membrane, suggesting these populations were subject to distinct evolutionary pressures that may correlate with genes related to nutrient uptake, biofilm formation, or viral invasion. These results are consistent with distinct evolutionary histories between geochemically different vent fields, with implications for understanding evolutionary processes in subseafloor microbial populations.

Emerging Diagnostic and Therapeutic Potentials of Human Hair Proteomics

The use of noninvasive human substrates to interrogate pathophysiological conditions has become essential in the post- Human Genome Project era. Due to its high turnover rate, and its long term capability to incorporate exogenous and endogenous substances from the circulation, hair testing is emerging as a key player in monitoring long term drug compliance, chronic alcohol abuse, forensic toxicology, and biomarker discovery, among other things. Novel high-throughput 'omics based approaches like proteomics have been underutilized globally in comprehending human hair morphology and its evolving use as a diagnostic testing substrate in the era of precision medicine. There is paucity of scientific evidence that evaluates the difference in drug incorporation into hair based on lipid content, and very few studies have addressed hair growth rates, hair forms, and the biological consequences of hair grooming or bleaching. It is apparent that protein-based identification using the human hair proteome would play a major role in understanding these parameters akin to DNA single nucleotide polymorphism profiling, up to single amino acid polymorphism resolution. Hence, this work seeks to identify and discuss the progress made thus far in the field of molecular hair testing using proteomic approaches, and identify ways in which proteomics would improve the field of hair research, considering that the human hair is mostly composed of proteins. Gaps in hair proteomics research are identified and the potential of hair proteomics in establishing a historic medical repository of normal and disease-specific proteome is also discussed.

Acetoclastic Methanosaeta are dominant methanogens in organic-rich Antarctic marine sediments

Despite accounting for the majority of sedimentary methane, the physiology and relative abundance of subsurface methanogens remain poorly understood. We combined intact polar lipid and metagenome techniques to better constrain the presence and functions of methanogens within the highly reducing, organic-rich sediments of Antarctica's Adélie Basin. The assembly of metagenomic sequence data identified phylogenic and functional marker genes of methanogens and generated the first Methanosaeta sp. genome from a deep subsurface sedimentary environment. Based on structural and isotopic measurements, glycerol dialkyl glycerol tetraethers with diglycosyl phosphatidylglycerol head groups were classified as biomarkers for active methanogens. The stable carbon isotope (δ¹³C) values of these biomarkers and the Methanosaeta partial genome suggest that these organisms are acetoclastic methanogens and represent a relatively small (0.2%) but active population. Metagenomic and lipid analyses suggest that Thaumarchaeota and heterotrophic bacteria co-exist with Methanosaeta and together contribute to increasing concentrations and δ¹³C values of dissolved inorganic carbon with depth. This study presents the first functional insights of deep subsurface Methanosaeta organisms and highlights their role in methane production and overall carbon cycling within sedimentary environments.

Pathogenicity of Human ST23 Streptococcus agalactiae to Fish and Genomic Comparison of Pathogenic and Non-pathogenic Isolates

Streptococcus agalactiae, or Group B Streptococcus (GBS), is a major pathogen causing neonatal sepsis and meningitis, bovine mastitis, and fish meningoencephalitis. CC23, including its namesake ST23, is not only the predominant GBS strain derived from human and cattle, but also can infect a variety of homeothermic and poikilothermic species. However, it has never been characterized in fish. This study aimed to determine the pathogenicity of ST23 GBS to fish and explore the mechanisms causing the difference in the pathogenicity of ST23 GBS based on the genome analysis. Infection of tilapia with 10 human-derived ST23 GBS isolates caused tissue damage and the distribution of pathogens within tissues. The mortality rate of infection was ranged from 76 to 100%, and it was shown that the mortality rate caused by only three human isolates had statistically significant difference compared with fish-derived ST7 strain (P < 0.05), whereas the mortality caused by other seven human isolates did not show significant difference compared with fish-derived ST7 strain. The genome comparison and prophage analysis showed that the major genome difference between virulent and non-virulent ST23 GBS was attributed to the different prophage sequences. The prophage in the P1 region contained about 43% GC and encoded 28–39 proteins, which can mediate the acquisition of YafQ/DinJ structure for GBS by phage recombination. YafQ/DinJ belongs to one of the bacterial toxin–antitoxin (TA) systems and allows cells to cope with stress. The ST23 GBS strains carrying this prophage were not pathogenic to tilapia, but the strains without the prophage or carrying the pophage that had gene mutation or deletion, especially the deletion of YafQ/DinJ structure, were highly pathogenic to tilapia. In conclusion, human ST23 GBS is highly pathogenic to fish, which may be related to the phage recombination.

Shotgun metagenomics, from sampling to analysis

Diverse microbial communities of bacteria, archaea, viruses and single-celled eukaryotes have crucial roles in the environment and in human health. However, microbes are frequently difficult to culture in the laboratory, which can confound cataloging of members and understanding of how communities function. High-throughput sequencing technologies and a suite of computational pipelines have been combined into shotgun metagenomics methods that have transformed microbiology. Still, computational approaches to overcome the challenges that affect both assembly-based and mapping-based metagenomic profiling, particularly of high-complexity samples or environments containing organisms with limited similarity to sequenced genomes, are needed. Understanding the functions and characterizing specific strains of these communities offers biotechnological promise in therapeutic discovery and innovative ways to synthesize products using microbial factories and can pinpoint the contributions of microorganisms to planetary, animal and human health.

DESMAN: a new tool for de novo extraction of strains from metagenomes

We introduce DESMAN for De novo Extraction of Strains from Metagenomes. Large multi-sample metagenomes are being generated but strain variation results in fragmentary co-assemblies. Current algorithms can bin contigs into metagenome-assembled genomes but are unable to resolve strain-level variation. DESMAN identifies variants in core genes and uses co-occurrence across samples to link variants into haplotypes and abundance profiles. These are then searched for against non-core genes to determine the accessory genome of each strain. We validated DESMAN on a complex 50-species 210-genome 96-sample synthetic mock data set and then applied it to the Tara Oceans microbiome.

Microbial diversity in the hypersaline Lake Meyghan, Iran

Lake Meyghan is one of the largest and commercially most important salt lakes in Iran. Despite its inland location and high altitude, Lake Meyghan has a thalassohaline salt composition suggesting a marine origin. Inputs of fresh water by rivers and rainfall formed various basins characterized by different salinities. We analyzed the microbial community composition of three basins by isolation and culturing of microorganisms and by analysis of the metagenome. The basins that were investigated comprised a green ~50 g kg⁻¹ salinity brine, a red ~180 g kg⁻¹ salinity brine and a white ~300 g kg⁻¹ salinity brine. Using different growth media, 57 strains of Bacteria and 48 strains of Archaea were isolated. Two bacterial isolates represent potential novel species with less than 96% 16S rRNA gene sequence identity to known species. Abundant isolates were also well represented in the metagenome. Bacteria dominated the low salinity brine, with Alteromonadales (Gammaproteobacteria) as a particularly important taxon, whereas the high salinity brines were dominated by haloarchaea. Although the brines of Lake Meyghan differ in geochemical composition, their ecosystem function appears largely conserved amongst each other while being driven by different microbial communities.

Capturing the genetic makeup of the active microbiome in situ

More than any other technology, nucleic acid sequencing has enabled microbial ecology studies to be complemented with the data volumes necessary to capture the extent of microbial diversity and dynamics in a wide range of environments. In order to truly understand and predict environmental processes, however, the distinction between active, inactive and dead microbial cells is critical. Also, experimental designs need to be sensitive toward varying population complexity and activity, and temporal as well as spatial scales of process rates. There are a number of approaches, including single-cell techniques, which were designed to study in situ microbial activity and that have been successively coupled to nucleic acid sequencing. The exciting new discoveries regarding in situ microbial activity provide evidence that future microbial ecology studies will indispensably rely on techniques that specifically capture members of the microbiome active in the environment. Herein, we review those currently used activity-based approaches that can be directly linked to shotgun nucleic acid sequencing, evaluate their relevance to ecology studies, and discuss future directions.

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.

Integrated meta-omic analyses of the gastrointestinal tract microbiome in patients undergoing allogeneic hematopoietic stem cell transplantation

In patients undergoing allogeneic hematopoietic stem cell transplantation (allo-HSCT), treatment-induced changes to the gastrointestinal tract (GIT) microbiome have been linked to adverse outcomes, most notably graft-versus-host disease (GvHD). However, it is presently unknown whether this relationship is causal or consequential. Here, we performed an integrated meta-omic analysis to probe deeper into the GIT microbiome changes during allo-HSCT and its accompanying treatments. We used 16S and 18S rRNA gene amplicon sequencing to resolve archaea, bacteria, and eukaryotes within the GIT microbiomes of 16 patients undergoing allo-HSCT for the treatment of hematologic malignancies. These results revealed a major shift in the GIT microbiome after allo-HSCT including a marked reduction in bacterial diversity, accompanied by only limited changes in eukaryotes and archaea. An integrated analysis of metagenomic and metatranscriptomic data was performed on samples collected from a patient before and after allo-HSCT for acute myeloid leukemia. This patient developed severe GvHD, leading to death 9 months after allo-HSCT. In addition to drastically decreased bacterial diversity, the post-treatment microbiome showed a higher overall number and higher expression levels of antibiotic resistance genes (ARGs). One specific Escherichia coli strain causing a paravertebral abscess was linked to GIT dysbiosis, suggesting loss of intestinal barrier integrity. The apparent selection for bacteria expressing ARGs suggests that prophylactic antibiotic administration may adversely affect the overall treatment outcome. We therefore assert that such analyses including information about the selection of pathogenic bacteria expressing ARGs may assist clinicians in "personalizing" regimens for individual patients to improve overall outcomes.

RNA-Based Amplicon Sequencing Reveals Microbiota Development during Ripening of Artisanal versus Industrial Lard d'Arnad

Valle d'Aosta Lard d'Arnad is a protected designation of origin (PDO) product produced from fat of the shoulder and back of heavy pigs. Its manufacturing process can be very diverse, especially regarding the maturation temperature and the NaCl concentration used for the brine; thereby, the main goal of this study was to investigate the impact of those parameters on the microbiota developed during curing and ripening. Three farms producing Lard d'Arnad were selected. Two plants, reflecting the industrial process characterized either by low maturation temperature (plant A [10% NaCl, 2°C]) or by using a low NaCl concentration (plant B [2.5% NaCl, 4°C]), were selected, while the third was characterized by an artisanal process (plant C [30% NaCl, 8°C]). Lard samples were obtained at time 0 and after 7, 15, 30, 60, and 90 days of maturation. From each plant, 3 independent lots were analyzed. The diversity of live microbiota was evaluated by using classical plate counts and amplicon target sequencing of small subunit (SSU) rRNA. The main taxa identified by sequencing were Acinetobacter johnsonii, Psychrobacter, Staphylococcus equorum, Staphylococcus sciuri, Pseudomonas fragi, Brochothrix, Halomonas, and Vibrio, and differences in their relative abundances distinguished samples from the individual plants. The composition of the microbiota was more similar among plants A and B, and it was characterized by the higher presence of taxa recognized as undesired bacteria in food-processing environments. Oligotype analysis of Halomonas and Acinetobacter revealed the presence of several characteristic oligotypes associated with A and B samples.IMPORTANCE Changes in the food production process can drastically affect the microbial community structure, with a possible impact on the final characteristics of the products. The industrial processes of Lard d'Arnad production are characterized by a reduction in the salt concentration in the brines to address a consumer demand for less salty products; this can negatively affect the dynamics and development of the live microbiota and, as a consequence, can negatively impact the quality of the final product due to the higher abundance of spoilage bacteria. This study is an overview of the live microbiota that develop during lard manufacturing, and it highlights the importance of the use of traditional process to produce PDO from a spoilage perspective.

BlobTools: Interrogation of genome assemblies

The goal of many genome sequencing projects is to provide a complete representation of a target genome (or genomes) as underpinning data for further analyses. However, it can be problematic to identify which sequences in an assembly truly derive from the target genome(s) and which are derived from associated microbiome or contaminant organisms.  We present BlobTools, a modular command-line solution for visualisation, quality control and taxonomic partitioning of genome datasets. Using guanine+cytosine content of sequences, read coverage in sequencing libraries and taxonomy of sequence similarity matches, BlobTools can assist in primary partitioning of data, leading to improved assemblies, and screening of final assemblies for potential contaminants.  Through simulated paired-end read dataset,s containing a mixture of metazoan and bacterial taxa, we illustrate the main BlobTools workflow and suggest useful parameters for taxonomic partitioning of low-complexity metagenome assemblies.