Autometa: automated extraction of microbial genomes from individual shotgun metagenomes

From pollen to putrid: Comparative metagenomics reveals how microbiomes support dietary specialization in vulture bees

For most animals, the microbiome is key for nutrition and pathogen defence, and is often shaped by diet. Corbiculate bees, including honey bees, bumble bees, and stingless bees, share a core microbiome that has been shaped, at least in part, by the challenges associated with pollen digestion. However, three species of stingless bees deviate from the general rule of bees obtaining their protein exclusively from pollen (obligate pollinivores) and instead consume carrion as their sole protein source (obligate necrophages) or consume both pollen and carrion (facultative necrophages). These three life histories can provide missing insights into microbiome evolution associated with extreme dietary transitions. Here, we investigate, via shotgun metagenomics, the functionality of the microbiome across three bee diet types: obligate pollinivory, obligate necrophagy, and facultative necrophagy. We find distinct differences in microbiome composition and gene functional profiles between the diet types. Obligate necrophages and pollinivores have more specialized microbes, whereas facultative necrophages have a diversity of environmental microbes associated with several dietary niches. Our study suggests that necrophagous bee microbiomes may have evolved to overcome cellular stress and microbial competition associated with carrion. We hypothesize that the microbiome evolved social phenotypes, such as biofilms, that protect the bees from opportunistic pathogens present on carcasses, allowing them to overcome novel nutritional challenges. Whether specific microbes enabled diet shifts or diet shifts occurred first and microbial evolution followed requires further research to disentangle. Nonetheless, we find that necrophagous microbiomes, vertebrate and invertebrate alike, have functional commonalities regardless of their taxonomy.

Prevalence and diversity of TAL effector-like proteins in fungal endosymbiotic Mycetohabitans spp

Endofungal Mycetohabitans (formerly Burkholderia) spp. rely on a type III secretion system to deliver mostly unidentified effector proteins when colonizing their host fungus, Rhizopus microsporus. The one known secreted effector family from Mycetohabitans consists of homologues of transcription activator-like (TAL) effectors, which are used by plant pathogenic Xanthomonas and Ralstonia spp. to activate host genes that promote disease. These 'Burkholderia TAL-like (Btl)' proteins bind corresponding specific DNA sequences in a predictable manner, but their genomic target(s) and impact on transcription in the fungus are unknown. Recent phenotyping of Btl mutants of two Mycetohabitans strains revealed that the single Btl in one Mycetohabitans endofungorum strain enhances fungal membrane stress tolerance, while others in a Mycetohabitans rhizoxinica strain promote bacterial colonization of the fungus. The phenotypic diversity underscores the need to assess the sequence diversity and, given that sequence diversity translates to DNA targeting specificity, the functional diversity of Btl proteins. Using a dual approach to maximize capture of Btl protein sequences for our analysis, we sequenced and assembled nine Mycetohabitans spp. genomes using long-read PacBio technology and also mined available short-read Illumina fungal-bacterial metagenomes. We show that btl genes are present across diverse Mycetohabitans strains from Mucoromycota fungal hosts yet vary in sequences and predicted DNA binding specificity. Phylogenetic analysis revealed distinct clades of Btl proteins and suggested that Mycetohabitans might contain more species than previously recognized. Within our data set, Btl proteins were more conserved across M. rhizoxinica strains than across M. endofungorum, but there was also evidence of greater overall strain diversity within the latter clade. Overall, the results suggest that Btl proteins contribute to bacterial-fungal symbioses in myriad ways.

Short macrocyclic peptides in sponge genomes

Sponges (Porifera) contain many peptide-specialized metabolites with potent biological activities and significant roles in shaping marine ecology. It is well established that symbiotic bacteria produce bioactive "sponge" peptides, both on the ribosome (RiPPs) and nonribosomally. Here, we demonstrate that sponges themselves also produce many bioactive macrocyclic peptides, such as phakellistatins and related proline-rich macrocyclic peptides (PRMPs). Using the Stylissa carteri sponge transcriptome, methods were developed to find sequences encoding 46 distinct RiPP-type core peptides, of which ten encoded previously identified PRMP sequences. With this basis set, the genome and transcriptome of the sponge Axinella corrugata was interrogated to find 35 PRMP precursor peptides encoding 31 unique core peptide sequences. At least 11 of these produced cyclic peptides that were present in the sponge and could be characterized by mass spectrometry, including stylissamides A-D and seven previously undescribed compounds. Precursor peptides were encoded in the A. corrugata genome, confirming their animal origin. The peptides contained signal peptide sequences and highly repetitive recognition sequence-core peptide elements with up to 25 PRMP copies in a single precursor. In comparison to sponges without PRMPs, PRMP sponges are incredibly enriched in potentially secreted polypeptides, with >23,000 individual signal peptide encoding genes found in a single transcriptome. The similarities between PRMP biosynthetic genes and neuropeptides in terms of their biosynthetic logic suggest a fundamental biology linked to circular peptides, possibly indicating a widespread and underappreciated diversity of signaling peptide post-translational modifications across the animal kingdom.

Deficient butyrate metabolism in the intestinal microbiome is a potential risk factor for recurrent kidney stone disease

Intestinal microbiome dysbiosis is a known risk factor for recurrent kidney stone disease (KSD) with prior data suggesting a role for dysfunctional metabolic pathways other than those directly utilizing oxalate. To identify alternative mechanisms, the current study analyzed differences in the metabolic potential of intestinal microbiomes of patients (n = 17) and live-in controls (n = 17) and determined their relevance to increased risk for KSD using shotgun metagenomic sequencing. We found no differences in the abundance of genes associated with known oxalate degradation pathways, supporting the notion that dysfunction in other metabolic pathways plays a role in KSD. Further analysis showed decreased abundance of key enzymes involved in butyrate biosynthesis in patient intestinal microbiomes. Furthermore, de novo construction of microbial genomes showed that the majority of genes significantly enriched in non-stone formers are affiliated with Faecalibacterium prausnitzii, a major butyrate producer. Specifically pertaining to butyrate metabolism, the majority of abundant genes mapped back to F. prausnitzii, Alistipes spp., and Akkermansia muciniphila. No differences were observed in ascorbate or glyoxylate metabolic pathways. Collectively, these data suggest that impaired bacterial-associated butyrate metabolism may be an oxalate-independent mechanism that contributes to an increased risk for recurrent KSD. This indicates that the role of the intestinal microbiome in recurrent KSD is multi-factorial, which is representative of the highly intertwined metabolic nature of this complex environment. Future bacteria-based treatments must not be restricted to targeting only oxalate metabolism.

Methylotrophic methanogenesis in the Archaeoglobi revealed by cultivation of Ca. Methanoglobus hypatiae from a Yellowstone hot spring

Over the past decade, environmental metagenomics and polymerase chain reaction-based marker gene surveys have revealed that several lineages beyond just a few well-established groups within the Euryarchaeota superphylum harbor the genetic potential for methanogenesis. One of these groups are the Archaeoglobi, a class of thermophilic Euryarchaeota that have long been considered to live non-methanogenic lifestyles. Here, we enriched Candidatus Methanoglobus hypatiae, a methanogen affiliated with the family Archaeoglobaceae, from a hot spring in Yellowstone National Park. The enrichment is sediment-free, grows at 64-70°C and a pH of 7.8, and produces methane from mono-, di-, and tri-methylamine. Ca. M. hypatiae is represented by a 1.62 Mb metagenome-assembled genome with an estimated completeness of 100% and accounts for up to 67% of cells in the culture according to fluorescence in situ hybridization. Via genome-resolved metatranscriptomics and stable isotope tracing, we demonstrate that Ca. M. hypatiae expresses methylotrophic methanogenesis and energy-conserving pathways for reducing monomethylamine to methane. The detection of Archaeoglobi populations related to Ca. M. hypatiae in 36 geochemically diverse geothermal sites within Yellowstone National Park, as revealed through the examination of previously published gene amplicon datasets, implies a previously underestimated contribution to anaerobic carbon cycling in extreme ecosystems.

The polyketide to fatty acid transition in the evolution of animal lipid metabolism

Animals synthesize simple lipids using a distinct fatty acid synthase (FAS) related to the type I polyketide synthase (PKS) enzymes that produce complex specialized metabolites. The evolutionary origin of the animal FAS and its relationship to the diversity of PKSs remain unclear despite the critical role of lipid synthesis in cellular metabolism. Recently, an animal FAS-like PKS (AFPK) was identified in sacoglossan molluscs. Here, we explore the phylogenetic distribution of AFPKs and other PKS and FAS enzymes across the tree of life. We found AFPKs widely distributed in arthropods and molluscs (>6300 newly described AFPK sequences). The AFPKs form a clade with the animal FAS, providing an evolutionary link bridging the type I PKSs and the animal FAS. We found molluscan AFPK diversification correlated with shell loss, suggesting AFPKs provide a chemical defense. Arthropods have few or no PKSs, but our results indicate AFPKs contributed to their ecological and evolutionary success by facilitating branched hydrocarbon and pheromone biosynthesis. Although animal metabolism is well studied, surprising new metabolic enzyme classes such as AFPKs await discovery.

Resistance mechanisms for Gram-negative bacteria-specific lipopeptides, turnercyclamycins, differ from that of colistin

IMPORTANCE

Bacterial resistance to antibiotics is a crisis. Acinetobacter baumannii is among the CDC urgent threat pathogens in part for this reason. Lipopeptides known as turnercyclamycins are produced by symbiotic bacteria that normally live in marine mollusks, where they may be involved in shaping their symbiotic niche. Turnercyclamycins killed Gram-negative pathogens including drug-resistant Acinetobacter, but how do the mechanisms of resistance compare to other lipopeptide drugs? Here, we define resistance from a truncation of MlaA, a protein involved in regulating bacterial membrane phospholipids. Intriguingly, this resistance mechanism only affected one turnercyclamycin variant, which differed only in two atoms in the lipid tail of the compounds. We could not obtain significant resistance to the second turnercyclamycin variant, which was also effective in an infection model. This study reveals an unexpected subtlety in resistance to lipopeptide antibiotics, which may be useful in the design and development of antibiotics to combat drug resistance.

Artificial intelligence for natural product drug discovery

Developments in computational omics technologies have provided new means to access the hidden diversity of natural products, unearthing new potential for drug discovery. In parallel, artificial intelligence approaches such as machine learning have led to exciting developments in the computational drug design field, facilitating biological activity prediction and de novo drug design for molecular targets of interest. Here, we describe current and future synergies between these developments to effectively identify drug candidates from the plethora of molecules produced by nature. We also discuss how to address key challenges in realizing the potential of these synergies, such as the need for high-quality datasets to train deep learning algorithms and appropriate strategies for algorithm validation.

Diversity and function of methyl-coenzyme M reductase-encoding archaea in Yellowstone hot springs revealed by metagenomics and mesocosm experiments

Metagenomic studies on geothermal environments have been central in recent discoveries on the diversity of archaeal methane and alkane metabolism. Here, we investigated methanogenic populations inhabiting terrestrial geothermal features in Yellowstone National Park (YNP) by combining amplicon sequencing with metagenomics and mesocosm experiments. Detection of methyl-coenzyme M reductase subunit A (mcrA) gene amplicons demonstrated a wide diversity of Mcr-encoding archaea inhabit geothermal features with differing physicochemical regimes across YNP. From three selected hot springs we recovered twelve Mcr-encoding metagenome assembled genomes (MAGs) affiliated with lineages of cultured methanogens as well as Candidatus (Ca.) Methanomethylicia, Ca. Hadesarchaeia, and Archaeoglobi. These MAGs encoded the potential for hydrogenotrophic, aceticlastic, hydrogen-dependent methylotrophic methanogenesis, or anaerobic short-chain alkane oxidation. While Mcr-encoding archaea represent minor fractions of the microbial community of hot springs, mesocosm experiments with methanogenic precursors resulted in the stimulation of methanogenic activity and the enrichment of lineages affiliated with Methanosaeta and Methanothermobacter as well as with uncultured Mcr-encoding archaea including Ca. Korarchaeia, Ca. Nezhaarchaeia, and Archaeoglobi. We revealed that diverse Mcr-encoding archaea with the metabolic potential to produce methane from different precursors persist in the geothermal environments of YNP and can be enriched under methanogenic conditions. This study highlights the importance of combining environmental metagenomics with laboratory-based experiments to expand our understanding of uncultured Mcr-encoding archaea and their potential impact on microbial carbon transformations in geothermal environments and beyond.

Changes in community structures and functions of the gut microbiomes of deep-sea cold seep mussels during in situ transplantation experiment

BACKGROUND

Many deep-sea invertebrates largely depend on chemoautotrophic symbionts for energy and nutrition, and some of them have reduced functional digestive tracts. By contrast, deep-sea mussels have a complete digestive system although symbionts in their gills play vital roles in nutrient supply. This digestive system remains functional and can utilise available resources, but the roles and associations among gut microbiomes in these mussels remain unknown. Specifically, how the gut microbiome reacts to environmental change is unclear.

RESULTS

The meta-pathway analysis showed the nutritional and metabolic roles of the deep-sea mussel gut microbiome. Comparative analyses of the gut microbiomes of original and transplanted mussels subjected to environmental change revealed shifts in bacterial communities. Gammaproteobacteria were enriched, whereas Bacteroidetes were slightly depleted. The functional response for the shifted communities was attributed to the acquisition of carbon sources and adjusting the utilisation of ammonia and sulphide. Self-protection was observed after transplantation.

CONCLUSION

This study provides the first metagenomic insights into the community structure and function of the gut microbiome in deep-sea chemosymbiotic mussels and their critical mechanisms for adapting to changing environments and meeting of essential nutrient demand.

Genetic mapping of microbial and host traits reveals production of immunomodulatory lipids by Akkermansia muciniphila in the murine gut

The molecular bases of how host genetic variation impacts the gut microbiome remain largely unknown. Here we used a genetically diverse mouse population and applied systems genetics strategies to identify interactions between host and microbe phenotypes including microbial functions, using faecal metagenomics, small intestinal transcripts and caecal lipids that influence microbe-host dynamics. Quantitative trait locus (QTL) mapping identified murine genomic regions associated with variations in bacterial taxa; bacterial functions including motility, sporulation and lipopolysaccharide production and levels of bacterial- and host-derived lipids. We found overlapping QTL for the abundance of Akkermansia muciniphila and caecal levels of ornithine lipids. Follow-up in vitro and in vivo studies revealed that A. muciniphila is a major source of these lipids in the gut, provided evidence that ornithine lipids have immunomodulatory effects and identified intestinal transcripts co-regulated with these traits including Atf3, which encodes for a transcription factor that plays vital roles in modulating metabolism and immunity. Collectively, these results suggest that ornithine lipids are potentially important for A. muciniphila-host interactions and support the role of host genetics as a determinant of responses to gut microbes.

Cyclitol metabolism is a central feature of Burkholderia leaf symbionts

The symbioses between plants of the Rubiaceae and Primulaceae families with Burkholderia bacteria represent unique and intimate plant-bacterial relationships. Many of these interactions have been identified through PCR-dependent typing methods, but there is little information available about their functional and ecological roles. We assembled 17 new endophyte genomes representing endophytes from 13 plant species, including those of two previously unknown associations. Genomes of leaf endophytes belonging to Burkholderia s.l. show extensive signs of genome reduction, albeit to varying degrees. Except for one endophyte, none of the bacterial symbionts could be isolated on standard microbiological media. Despite their taxonomic diversity, all endophyte genomes contained gene clusters linked to the production of specialized metabolites, including genes linked to cyclitol sugar analog metabolism and in one instance non-ribosomal peptide synthesis. These genes and gene clusters are unique within Burkholderia s.l. and are likely horizontally acquired. We propose that the acquisition of secondary metabolite gene clusters through horizontal gene transfer is a prerequisite for the evolution of a stable association between these endophytes and their hosts.

Taxonomic Assignment-Based Genome Reconstruction from Apical Periodontal Metagenomes to Identify Antibiotic Resistance and Virulence Factors

Primary apical periodontitis occurs due to various insults to the dental pulp including microbial infections, physical and iatrogenic trauma, whereas inadequate elimination of intraradicular infection during root canal treatment may lead to secondary apical periodontitis. We explored the complex intra-radicular microbial communities and their functional potential through genome reconstruction. We applied shotgun metagenomic sequencing, binning and functional profiling to identify the significant contributors to infection at the acute and chronic apical periodontal lesions. Our analysis revealed the five classified clusters representing Enterobacter, Enterococcus, Lacticaseibacillus, Pseudomonas, Streptococcus and one unclassified cluster of contigs at the genus level. Of them, the major contributors were Pseudomonas, with 90.61% abundance in acute conditions, whereas Enterobacter followed by Enterococcus with 69.88% and 15.42% abundance, respectively, in chronic conditions. Enterobacter actively participated in antibiotic target alteration following multidrug efflux-mediated resistance mechanisms, predominant in the chronic stage. The prediction of pathways involved in the destruction of the supportive tissues of the tooth in Enterobacter and Pseudomonas support their crucial role in the manifestation of respective disease conditions. This study provides information about the differential composition of the microbiome in chronic and acute apical periodontitis. It takes a step to interpret the role of a single pathogen, solely or predominantly, in establishing endodontic infection types through genome reconstruction following high throughput metagenomic DNA analysis. The resistome prediction sheds a new light on the therapeutic treatment guidelines for endodontists. However, it needs further conclusive research to support this outcome using a larger number of samples with similar etiological conditions, but different demographic origin.

Mycoparasites, Gut Dwellers, and Saprotrophs: Phylogenomic Reconstructions and Comparative Analyses of Kickxellomycotina Fungi

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

GutAlive® enables DNA-based microbiome analysis without disrupting the original composition and diversity

Introduction

A precise fecal microbiome analysis requires normalized methods for microbiome sampling, transport and manipulation in order to obtain a representative snapshot of the microbial community. GutAlive^® is the unique stool collection kit that generates an anaerobic atmosphere enabling oxygen sensitive bacteria to survive, maintaining the original microbiome composition and diversity.

Methods

Five stool samples from different donors were collected using two different sampling devices, GutAlive^® and Zymo DNA/RNA Shield^®, and processed at four different time points. Shotgun metagenomics was used to evaluate the influence of the device and the processing timing on the microbial populations to unravel the potential fluctuations on the composition and diversity of the fecal microbiome and the metabolic pathways profiling. Additionally, RT-qPCR was used to quantify bacterial cell viability for downstream applications of microbiota samples beyond metagenomics.

Results

Our results show that GutAlive^® enables bacterial cell viability overtime preserving DNA integrity, obtaining high-quantity and high-quality DNA to perform microbiome analysis using shotgun metagenomics. Based on the taxonomic profiling, metabolic pathways analysis, phylogeny and metagenome-assembled genomes, GutAlive^® displayed greater performance without significant variability over time, showcasing the stabilization of the microbiome preserving the original composition and diversity. Indeed, this DNA stabilization is enabled with the preservation of bacterial viability on an anaerobic environment inside of the sampling device, without the addition of any reagents that interact directly with sample.

Conclusion

All the above makes GutAlive^® an user-friendly kit for self-collection of biological samples, suitable for microbiome analysis, diagnostics, fecal microbiota transplant and bacterial isolation, maintaining the stability and bacterial viability over time, preserving the original composition and diversity of the microbiome.

Host innate immune responses and microbiome profile of neonatal calves challenged with Cryptosporidium parvum and the effect of bovine colostrum supplementation

Introduction

Calves are highly susceptible to gastrointestinal infection with Cryptosporidium parvum (C. parvum), which can result in watery diarrhea and eventually death or impaired development. With little to no effective therapeutics, understanding the host's microbiota and pathogen interaction at the mucosal immune system has been critical to identify and test novel control strategies.

Methods

Herein, we used an experimental model of C. parvum challenge in neonatal calves to describe the clinical signs and histological and proteomic profiling of the mucosal innate immunity and microbiota shifts by metagenomics in the ileum and colon during cryptosporidiosis. Also, we investigated the impact of supplemental colostrum feeding on C. parvum infection.

Results

We showed that C. parvum challenged calves experienced clinical signs including pyrexia and diarrhea 5 days post challenge. These calves showed ulcerative neutrophil ileitis with a proteomic signature driven by inflammatory effectors, including reactive oxygen species and myeloperoxidases. Colitis was also noticed with an aggravated mucin barrier depletion and incompletely filled goblet cells. The C. parvum challenged calves also displayed a pronounced dysbiosis with a high prevalence of Clostridium species (spp.) and number of exotoxins, adherence factors, and secretion systems related to Clostridium spp. and other enteropathogens, including Campylobacter spp., Escherichia sp., Shigella spp., and Listeria spp. Daily supplementation with a high-quality bovine colostrum product mitigated some of the clinical signs and modulated the gut immune response and concomitant microbiota to a pattern more similar to that of healthy unchallenged calves.

Discussion

C. parvum infection in neonatal calves provoked severe diarrheic neutrophilic enterocolitis, perhaps augmented due to the lack of fully developed innate gut defenses. Colostrum supplementation showed limited effect mitigating diarrhea but demonstrated some clinical alleviation and specific modulatory influence on host gut immune responses and concomitant microbiota.

Exploring microbial functional biodiversity at the protein family level-From metagenomic sequence reads to annotated protein clusters

Metagenomics has enabled accessing the genetic repertoire of natural microbial communities. Metagenome shotgun sequencing has become the method of choice for studying and classifying microorganisms from various environments. To this end, several methods have been developed to process and analyze the sequence data from raw reads to end-products such as predicted protein sequences or families. In this article, we provide a thorough review to simplify such processes and discuss the alternative methodologies that can be followed in order to explore biodiversity at the protein family level. We provide details for analysis tools and we comment on their scalability as well as their advantages and disadvantages. Finally, we report the available data repositories and recommend various approaches for protein family annotation related to phylogenetic distribution, structure prediction and metadata enrichment.

Whole-Genome Sequencing of Bacterial Endophytes From Fresh and Preserved Plant Specimens

Many plants harbor symbiotic bacteria in their leaves, sometimes within structures visible with the naked eye. These bacteria play critical roles for host development and defense, but are often not amenable to culture. Gaining insight into the functions of these obligate endophytic bacteria hinges on culture-independent omics approaches, which have seen tremendous development in recent years. We describe in this chapter a set of protocols for the extraction and bioinformatic analysis of bacterial genomic DNA from leaf samples of various origins, including fresh, silica-preserved, or herbarium specimens.

Binning Metagenomic Contigs Using Unsupervised Clustering and Reference Databases

Metagenomics can directly extract the genetic material of all microorganisms from the environment, and obtain metagenomic samples with a large number of unknown DNA sequences. Binning of metagenomic contigs is a hot topic in metagenomics research. There are two key challenges for the current unsupervised metagenomic clustering algorithms. First, unsupervised metagenomic clustering methods rarely use reference databases, causing a certain waste of resources. Second, unsupervised metagenomic clustering methods are restricted by the characteristics of the sequences and the clustering algorithms, and the binning effect is limited. Therefore, a new binning method for metagenomic contigs using unsupervised clustering methods and reference databases is proposed to address these challenges, to make full use of the advantages of unsupervised clustering methods and reference databases constructed by scientists to improve the overall binning effect. This method uses the integrated SVM classification model to further bin the unsupervised clustering parts that do not perform well. Our proposed method was tested on simulated datasets and a real dataset and compared with other state-of-the-art metagenomic clustering methods including CONCOCT, Metabin2.0, Autometa, and MetaBAT. The results show that our method can achieve higher precision rate and improve the binning effect.

Constructing metagenome-assembled genomes for almost all components in a real bacterial consortium for binning benchmarking

Background

So far, a lot of binning approaches have been intensively developed for untangling metagenome-assembled genomes (MAGs) and evaluated by two main strategies. The strategy by comparison to known genomes prevails over the other strategy by using single-copy genes. However, there is still no dataset with all known genomes for a real (not simulated) bacterial consortium yet.

Results

Here, we continue investigating the real bacterial consortium F1RT enriched and sequenced by us previously, considering the high possibility to unearth all MAGs, due to its low complexity. The improved F1RT metagenome reassembled by metaSPAdes here utilizes about 98.62% of reads, and a series of analyses for the remaining reads suggests that the possibility of containing other low-abundance organisms in F1RT is greatly low, demonstrating that almost all MAGs are successfully assembled. Then, 4 isolates are obtained and individually sequenced. Based on the 4 isolate genomes and the entire metagenome, an elaborate pipeline is then in-house developed to construct all F1RT MAGs. A series of assessments extensively prove the high reliability of the herein reconstruction. Next, our findings further show that this dataset harbors several properties challenging for binning and thus is suitable to compare advanced binning tools available now or benchmark novel binners. Using this dataset, 8 advanced binning algorithms are assessed, giving useful insights for developing novel approaches. In addition, compared with our previous study, two novel MAGs termed FC8 and FC9 are discovered here, and 7 MAGs are solidly unearthed for species without any available genomes.

Conclusion

To our knowledge, it is the first time to construct a dataset with almost all known MAGs for a not simulated consortium. We hope that this dataset will be used as a routine toolkit to complement mock datasets for evaluating binning methods to further facilitate binning and metagenomic studies in the future.

Uncovering Lasonolide A Biosynthesis Using Genome-Resolved Metagenomics

Invertebrates, particularly sponges, have been a dominant source of new marine natural products. For example, lasonolide A (LSA) is a potential anticancer molecule isolated from the marine sponge Forcepia sp., with nanomolar growth inhibitory activity and a unique cytotoxicity profile against the National Cancer Institute 60-cell-line screen. Here, we identified the putative biosynthetic pathway for LSA. Genomic binning of the Forcepia sponge metagenome revealed a Gram-negative bacterium belonging to the phylum Verrucomicrobia as the candidate producer of LSA. Phylogenetic analysis showed that this bacterium, here named "Candidatus Thermopylae lasonolidus," only has 88.78% 16S rRNA identity with the closest relative, Pedosphaera parvula Ellin514, indicating that it represents a new genus. The lasonolide A (las) biosynthetic gene cluster (BGC) was identified as a trans-acyltransferase (AT) polyketide synthase (PKS) pathway. Compared with its host genome, the las BGC exhibits a significantly different GC content and pentanucleotide frequency, suggesting a potential horizontal acquisition of the gene cluster. Furthermore, three copies of the putative las pathway were identified in the candidate producer genome. Differences between the three las repeats were observed, including the presence of three insertions, two single-nucleotide polymorphisms, and the absence of a stand-alone acyl carrier protein in one of the repeats. Even though the verrucomicrobial producer shows signs of genome reduction, its genome size is still fairly large (about 5 Mbp), and, compared to its closest free-living relative, it contains most of the primary metabolic pathways, suggesting that it is in the early stages of reduction. IMPORTANCE While sponges are valuable sources of bioactive natural products, a majority of these compounds are produced in small quantities by uncultured symbionts, hampering the study and clinical development of these unique compounds. Lasonolide A (LSA), isolated from marine sponge Forcepia sp., is a cytotoxic molecule active at nanomolar concentrations, which causes premature chromosome condensation, blebbing, cell contraction, and loss of cell adhesion, indicating a novel mechanism of action and making it a potential anticancer drug lead. However, its limited supply hampers progression to clinical trials. We investigated the microbiome of Forcepia sp. using culture-independent DNA sequencing, identified genes likely responsible for LSA synthesis in an uncultured bacterium, and assembled the symbiont's genome. These insights provide future opportunities for heterologous expression and cultivation efforts that may minimize LSA's supply problem.

A global phylogenomic and metabolic reconstruction of the large intestine bacterial community of domesticated cattle

BACKGROUND

The large intestine is a colonization site of beneficial microbes complementing the nutrition of cattle but also of zoonotic and animal pathogens. Here, we present the first global gene catalog of cattle fecal microbiomes, a proxy of the large intestine microbiomes, from 436 metagenomes from six countries.

RESULTS

Phylogenomics suggested that the reconstructed genomes and their close relatives form distinct branches and produced clustering patterns that were reminiscent of the metagenomics sample origin. Bacterial taxa had distinct metabolic profiles, and complete metabolic pathways were mainly linked to carbohydrates and amino acids metabolism. Dietary changes affected the community composition, diversity, and potential virulence. However, predicted enzymes, which were part of complete metabolic pathways, remained present, albeit encoded by different microbes.

CONCLUSIONS

Our findings provide a global insight into the phylogenetic relationships and the metabolic potential of a rich yet understudied bacterial community and suggest that it provides valuable services to the host. However, we tentatively infer that members of that community are not irreplaceable, because similar to previous findings, symbionts of complex bacterial communities of mammals are expendable if there are substitutes that can perform the same task. Video Abstract.

Culexarchaeia, a novel archaeal class of anaerobic generalists inhabiting geothermal environments

Geothermal environments, including terrestrial hot springs and deep-sea hydrothermal sediments, often contain many poorly understood lineages of archaea. Here, we recovered ten metagenome-assembled genomes (MAGs) from geothermal sediments and propose that they constitute a new archaeal class within the TACK superphylum, "Candidatus Culexarchaeia", named after the Culex Basin in Yellowstone National Park. Culexarchaeia harbor distinct sets of proteins involved in key cellular processes that are either phylogenetically divergent or are absent from other closely related TACK lineages, with a particular divergence in cell division and cytoskeletal proteins. Metabolic reconstruction revealed that Culexarchaeia have the capacity to metabolize a wide variety of organic and inorganic substrates. Notably, Culexarchaeia encode a unique modular, membrane associated, and energy conserving [NiFe]-hydrogenase complex that potentially interacts with heterodisulfide reductase (Hdr) subunits. Comparison of this [NiFe]-hydrogenase complex with similar complexes from other archaea suggests that interactions between membrane associated [NiFe]-hydrogenases and Hdr may be more widespread than previously appreciated in both methanogenic and non-methanogenic lifestyles. The analysis of Culexarchaeia further expands our understanding of the phylogenetic and functional diversity of lineages within the TACK superphylum and the ecology, physiology, and evolution of these organisms in extreme environments.

The Role of Feeding Characteristics in Shaping Gut Microbiota Composition and Function of Ensifera (Orthoptera)

Feeding habits were the primary factor affecting the gut bacterial communities in Ensifera. However, the interaction mechanism between the gut microbiota and feeding characteristics is not precisely understood. Here, the gut microbiota of Ensifera with diverse feeding habits was analyzed by shotgun metagenomic sequencing to further clarify the composition and function of the gut microbiota and its relationship with feeding characteristics. Our results indicate that under the influence of feeding habits, the gut microbial communities of Ensifera showed specific characteristics. Firstly, the gut microbial communities of the Ensifera with different feeding habits differed significantly, among which the gut microbial diversity of the herbivorous Mecopoda niponensis was the highest. Secondly, the functional genes related to feeding habits were in high abundance. Thirdly, the specific function of the gut microbial species in the omnivorous Gryllotalpa orientalis showed that the more diverse the feeding behavior of Ensifera, the worse the functional specificity related to the feeding characteristics of its gut microbiota. However, feeding habits were not the only factors affecting the gut microbiota of Ensifera. Some microorganisms' genes, whose functions were unrelated to feeding characteristics but were relevant to energy acquisition and nutrient absorption, were detected in high abundance. Our results were the first to report on the composition and function of the gut microbiota of Ensifera based on shotgun metagenomic sequencing and to explore the potential mechanism of the gut microbiota's association with diverse feeding habits.

Ancient plant-like terpene biosynthesis in corals

Octocorals are major contributors of terpenoid chemical diversity in the ocean. Natural products from other sessile marine animals are primarily biosynthesized by symbiotic microbes rather than by the host. Here, we challenge this long-standing paradigm by describing a monophyletic lineage of animal-encoded terpene cyclases (TCs) ubiquitous in octocorals. We characterized 15 TC enzymes from nine genera, several of which produce precursors of iconic coral-specific terpenoids, such as pseudopterosin, lophotoxin and eleutherobin. X-ray crystallography revealed that coral TCs share conserved active site residues and structural features with bacterial TCs. The identification of coral TCs enabled the targeted identification of the enzyme that constructs the coral-exclusive capnellane scaffold. Several TC genes are colocalized with genes that encode enzymes known to modify terpenes. This work presents an example of biosynthetic capacity in the kingdom Animalia that rivals the chemical complexity generated by plants, unlocking the biotechnological potential of octocorals for biomedical applications.

Comparative functional analysis of the urinary tract microbiome for individuals with or without calcium oxalate calculi

Individuals with urinary stone disease (USD) exhibit dysbiosis in the urinary tract and the loss of Lactobacillus that promote urinary tract health. However, the microbial metabolic functions that differentiate individuals with USD from healthy individuals are unknown. The objective of the current study was to determine the microbial functions across prokaryotic, viral, fungal, and protozoan domains that are associated with calcium oxalate (CaOx) stone formers through comparative shotgun metagenomics of midstream, voided urine samples for a small number of patients (n = 5 CaOx stone formers, n = 5 healthy controls). Results revealed that CaOx stone formers had reduced levels of genes associated with oxalate metabolism, as well as transmembrane transport, proteolysis, and oxidation-reduction processes. From 17 draft genomes extracted from the data and > 42,000 full length reference genomes, genes enriched in the Control group mapped overwhelming to Lactobacillus crispatus and those associated with CaOx mapped to Pseudomonas aeruginosa and Burkholderia sp. The microbial functions that differentiated the clinical cohorts are associated with known mechanisms of stone formation. While the prokaryotes most differentiated the CaOx and Control groups, a diverse, trans-domain microbiome was apparent. While our sample numbers were small, results corroborate previous studies and suggest specific microbial metabolic pathways in the urinary tract that modulate stone formation. Future studies that target these metabolic pathways as well as the influence of viruses, fungi, and protozoa on urinary tract physiology is warranted.

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 Provides Insight into the Function of Broad-Host Range Sponge Symbionts

The fossil record indicates that the earliest evidence of extant marine sponges (phylum Porifera) existed during the Cambrian explosion and that their symbiosis with microbes may have begun in their extinct ancestors during the Precambrian period. Many symbionts have adapted to their sponge host, where they perform specific, specialized functions. There are also widely distributed bacterial taxa such as Poribacteria, SAUL, and Tethybacterales that are found in a broad range of invertebrate hosts. Here, we added 11 new genomes to the Tethybacterales order, identified a novel family, and show that functional potential differs between the three Tethybacterales families. We compare the Tethybacterales with the well-characterized Entoporibacteria and show that these symbionts appear to preferentially associate with low-microbial abundance (LMA) and high-microbial abundance (HMA) sponges, respectively. Within these sponges, we show that these symbionts likely perform distinct functions and may have undergone multiple association events, rather than a single association event followed by coevolution. IMPORTANCE Marine sponges often form symbiotic relationships with bacteria that fulfil a specific need within the sponge holobiont, and these symbionts are often conserved within a narrow range of related taxa. To date, there exist only three known bacterial taxa (Entoporibacteria, SAUL, and Tethybacterales) that are globally distributed and found in a broad range of sponge hosts, and little is known about the latter two. We show that the functional potential of broad-host range symbionts is conserved at a family level and that these symbionts have been acquired several times over evolutionary history. Finally, it appears that the Entoporibacteria are associated primarily with high-microbial abundance sponges, while the Tethybacterales associate with low-microbial abundance sponges.

Hydrodynamic disturbance controls microbial community assembly and biogeochemical processes in coastal sediments

The microbial community composition and biogeochemical dynamics of coastal permeable (sand) sediments differs from cohesive (mud) sediments. Tide- and wave-driven hydrodynamic disturbance causes spatiotemporal variations in oxygen levels, which select for microbial generalists and disrupt redox cascades. In this work, we profiled microbial communities and biogeochemical dynamics in sediment profiles from three sites varying in their exposure to hydrodynamic disturbance. Strong variations in sediment geochemistry, biogeochemical activities, and microbial abundance, composition, and capabilities were observed between the sites. Most of these variations, except for microbial abundance and diversity, significantly correlated with the relative disturbance level of each sample. In line with previous findings, metabolically flexible habitat generalists (e.g., Flavobacteriaceae, Woeseaiceae, Rhodobacteraceae) dominated in all samples. However, we present evidence that aerobic specialists such as ammonia-oxidizing archaea (Nitrosopumilaceae) were more abundant and active in more disturbed samples, whereas bacteria capable of sulfate reduction (e.g., uncultured Desulfobacterales), dissimilatory nitrate reduction to ammonium (DNRA; e.g., Ignavibacteriaceae), and sulfide-dependent chemolithoautotrophy (e.g., Sulfurovaceae) were enriched and active in less disturbed samples. These findings are supported by insights from nine deeply sequenced metagenomes and 169 derived metagenome-assembled genomes. Altogether, these findings suggest that hydrodynamic disturbance is a critical factor controlling microbial community assembly and biogeochemical processes in coastal sediments. Moreover, they strengthen our understanding of the relationships between microbial composition and biogeochemical processes in these unique environments.

Arms race in a cell: genomic, transcriptomic, and proteomic insights into intracellular phage-bacteria interplay in deep-sea snail holobionts

BACKGROUND

Deep-sea animals in hydrothermal vents often form endosymbioses with chemosynthetic bacteria. Endosymbionts serve essential biochemical and ecological functions, but the prokaryotic viruses (phages) that determine their fate are unknown.

RESULTS

We conducted metagenomic analysis of a deep-sea vent snail. We assembled four genome bins for Caudovirales phages that had developed dual endosymbiosis with sulphur-oxidising bacteria (SOB) and methane-oxidising bacteria (MOB). Clustered regularly interspaced short palindromic repeat (CRISPR) spacer mapping, genome comparison, and transcriptomic profiling revealed that phages Bin1, Bin2, and Bin4 infected SOB and MOB. The observation of prophages in the snail endosymbionts and expression of the phage integrase gene suggested the presence of lysogenic infection, and the expression of phage structural protein and lysozyme genes indicated active lytic infection. Furthermore, SOB and MOB appear to employ adaptive CRISPR-Cas systems to target phage DNA. Additional expressed defence systems, such as innate restriction-modification systems and dormancy-inducing toxin-antitoxin systems, may co-function and form multiple lines for anti-viral defence. To counter host defence, phages Bin1, Bin2, and Bin3 appear to have evolved anti-restriction mechanisms and expressed methyltransferase genes that potentially counterbalance host restriction activity. In addition, the high-level expression of the auxiliary metabolic genes narGH, which encode nitrate reductase subunits, may promote ATP production, thereby benefiting phage DNA packaging for replication.

CONCLUSIONS

This study provides new insights into phage-bacteria interplay in intracellular environments of a deep-sea vent snail. Video Abstract.

A high-quality carabid genome assembly provides insights into beetle genome evolution and cold adaptation

The hyperdiverse order Coleoptera comprises a staggering ~25% of known species on Earth. Despite recent breakthroughs in next generation sequencing, there remains a limited representation of beetle diversity in assembled genomes. Most notably, the ground beetle family Carabidae, comprising more than 40,000 described species, has not been studied in a comparative genomics framework using whole genome data. Here we generate a high-quality genome assembly for Nebria riversi, to examine sources of novelty in the genome evolution of beetles, as well as genetic changes associated with specialization to high-elevation alpine habitats. In particular, this genome resource provides a foundation for expanding comparative molecular research into mechanisms of insect cold adaptation. Comparison to other beetles shows a strong signature of genome compaction, with N. riversi possessing a relatively small genome (~147 Mb) compared to other beetles, with associated reductions in repeat element content and intron length. Small genome size is not, however, associated with fewer protein-coding genes, and an analysis of gene family diversity shows significant expansions of genes associated with cellular membranes and membrane transport, as well as protein phosphorylation and muscle filament structure. Finally, our genomic analyses show that these high-elevation beetles have endosymbiotic Spiroplasma, with several metabolic pathways (e.g., propanoate biosynthesis) that might complement N. riversi, although its role as a beneficial symbiont or as a reproductive parasite remains equivocal.

Bacteriome depiction and the trophic status of the largest Northern highland lake from Andes system: Lago de Tota, Boyacá, Colombia

Lago de Tota is the largest highland lake in Colombia and one of the most remarkable of Northern Andean Mountain range. This lake is under an anthropogenic-based eutrophication process as a consequence of non-sustainable agriculture practices developing nearby. Notable relationship between the trophic status and Bacteriome loop dynamics has been increasingly disclosed in lakes worldwide. We performed a 16S sequencing analysis to depict the bacterial community present and we inferred its potential gene function in Lago de Tota. Parameters for determining current trophic condition such as total nitrogen (TN), dissolved carbon (DOC), particulate organic matter (POM), and chlorophyll-a (chl-a) were measured. A total of 440 Operational Taxonomic Units (OTUs) arranged into 50 classes were identified based on V3-V4 regions of the 16S rRNA gene, harboring high-frequent likely found environmental classes such as Actinobacteria, Gammaproteobacteria, Bacteroidia, Acidimicrobia, and Verrucomicrobiae. A total of 26 bacterial classes configure most abundant predicted functional processes involved in organic matter decomposition (i.e., carbohydrate metabolism, amino acid metabolism, xenobiotic biodegradation, and energy metabolism). In general, Actinobacteria, Alphaproteobacteria, and Gammaproteobacteria show the highest potential gene functional contributors, although other low-frequent classes OTUs are also relevant in processes of carbohydrate metabolism, xenobiotic biodegradation, and energy metabolism. The Trophic State Index indicates an oligo-mesotrophic status, and additional variables measured (i.e., POM, DOC) suggest the increasing carbon accumulation. Results provide preliminary evidence for several bacteria groups related to eutrophication of Lago de Tota. Under this picture, we suggest that further studies for Bacteriome loop spatial-temporal description are essential to inform local water quality monitoring strategies.

First finding of free-living representatives of Prokinetoplastina and their nuclear and mitochondrial genomes

Kinetoplastids are heterotrophic flagellated protists, including important parasites of humans and animals (trypanosomatids), and ecologically important free-living bacterial consumers (bodonids). Phylogenies have shown that the earliest-branching kinetoplastids are all parasites or obligate endosymbionts, whose highly-derived state makes reconstructing the ancestral state of the group challenging. We have isolated new strains of unusual free-living flagellates that molecular phylogeny shows to be most closely related to endosymbiotic and parasitic Perkinsela and Ichthyobodo species that, together with unidentified environmental sequences, form the clade at the base of kinetoplastids. These strains are therefore the first described free-living prokinetoplastids, and potentially very informative in understanding the evolution and ancestral states of morphological and molecular characteristics described in other kinetoplastids. Overall, we find that these organisms morphologically and ultrastructurally resemble some free-living bodonids and diplonemids, and possess nuclear genomes with few introns, polycistronic mRNA expression, high coding density, and derived traits shared with other kinetoplastids. Their genetic repertoires are more diverse than the best-studied free-living kinetoplastids, which is likely a reflection of their higher metabolic potential. Mitochondrial RNAs of these new species undergo the most extensive U insertion/deletion editing reported so far, and limited deaminative C-to-U and A-to-I editing, but we find no evidence for mitochondrial trans-splicing.

Genomic profiling of bacterial and fungal communities and their predictive functionality during pulque fermentation by whole-genome shotgun sequencing

Pulque is a culturally important 4,000-year-old traditional Mexican fermented drink. Pulque is produced by adding fresh aguamiel (agave sap) to mature pulque, resulting in a mixture of microbial communities and chemical compositions. We performed shotgun metagenomic sequencing of five stages of pulque fermentation to characterize organismal and functional diversity. We identified 6 genera (Acinetobacter, Lactobacillus, Lactococcus, Leuconostoc, Saccharomyces and Zymomonas) and 10 species (Acinetobacter boissieri, Acinetobacter nectaris, Lactobacillus sanfranciscensis, Lactococcus lactis, Lactococcus piscium, Lactococcus plantarum, Leuconostoc citreum, Leuconostoc gelidum, Zymomonas mobilis and Saccharomyces cerevisiae) that were present ≥ 1% in at least one stage of pulque fermentation. The abundance of genera and species changed during fermentation and was associated with a decrease in sucrose and increases in ethanol and lactic acid, suggesting that resource competition shapes organismal diversity. We also predicted functional profiles, based on organismal gene content, for each fermentation stage and identified an abundance of genes associated with the biosynthesis of folate, an essential B-vitamin. Additionally, we investigated the evolutionary relationships of S. cerevisiae and Z. mobilis, two of the major microbial species found in pulque. For S. cerevisiae, we used a metagenomics assembly approach to identify S. cerevisiae scaffolds from pulque, and performed phylogenetic analysis of these sequences along with a collection of 158 S. cerevisiae strains. This analysis suggests that S. cerevisiae from pulque is most closely related to Asian strains isolated from sake and bioethanol. Lastly, we isolated and sequenced the whole-genomes of three strains of Z. mobilis from pulque and compared their relationship to seven previously sequenced isolates. Our results suggest pulque strains may represent a distinct lineage of Z. mobilis.

Metagenomic approaches in microbial ecology: an update on whole-genome and marker gene sequencing analyses

Metagenomics and marker gene approaches, coupled with high-throughput sequencing technologies, have revolutionized the field of microbial ecology. Metagenomics is a culture-independent method that allows the identification and characterization of organisms from all kinds of samples. Whole-genome shotgun sequencing analyses the total DNA of a chosen sample to determine the presence of micro-organisms from all domains of life and their genomic content. Importantly, the whole-genome shotgun sequencing approach reveals the genomic diversity present, but can also give insights into the functional potential of the micro-organisms identified. The marker gene approach is based on the sequencing of a specific gene region. It allows one to describe the microbial composition based on the taxonomic groups present in the sample. It is frequently used to analyse the biodiversity of microbial ecosystems. Despite its importance, the analysis of metagenomic sequencing and marker gene data is quite a challenge. Here we review the primary workflows and software used for both approaches and discuss the current challenges in the field.

Evaluating metagenomics tools for genome binning with real metagenomic datasets and CAMI datasets

Background

Shotgun metagenomics based on untargeted sequencing can explore the taxonomic profile and the function of unknown microorganisms in samples, and complement the shortage of amplicon sequencing. Binning assembled sequences into individual groups, which represent microbial genomes, is the key step and a major challenge in metagenomic research. Both supervised and unsupervised machine learning methods have been employed in binning. Genome binning belonging to unsupervised method clusters contigs into individual genome bins by machine learning methods without the assistance of any reference databases. So far a lot of genome binning tools have emerged. Evaluating these genome tools is of great significance to microbiological research. In this study, we evaluate 15 genome binning tools containing 12 original binning tools and 3 refining binning tools by comparing the performance of these tools on chicken gut metagenomic datasets and the first CAMI challenge datasets.

Results

For chicken gut metagenomic datasets, original genome binner MetaBat, Groopm2 and Autometa performed better than other original binner, and MetaWrap combined the binning results of them generated the most high-quality genome bins. For CAMI datasets, Groopm2 achieved the highest purity (> 0.9) with good completeness (> 0.8), and reconstructed the most high-quality genome bins among original genome binners. Compared with Groopm2, MetaBat2 had similar performance with higher completeness and lower purity. Genome refining binners DASTool predicated the most high-quality genome bins among all genomes binners. Most genome binner performed well for unique strains. Nonetheless, reconstructing common strains still is a substantial challenge for all genome binner.

Conclusions

In conclusion, we tested a set of currently available, state-of-the-art metagenomics hybrid binning tools and provided a guide for selecting tools for metagenomic binning by comparing range of purity, completeness, adjusted rand index, and the number of high-quality reconstructed bins. Furthermore, available information for future binning strategy were concluded.

MinYS: mine your symbiont by targeted genome assembly in symbiotic communities

Most metazoans are associated with symbionts. Characterizing the effect of a particular symbiont often requires getting access to its genome, which is usually done by sequencing the whole community. We present MinYS, a targeted assembly approach to assemble a particular genome of interest from such metagenomic data. First, taking advantage of a reference genome, a subset of the reads is assembled into a set of backbone contigs. Then, this draft assembly is completed using the whole metagenomic readset in a de novo manner. The resulting assembly is output as a genome graph, enabling different strains with potential structural variants coexisting in the sample to be distinguished. MinYS was applied to 50 pea aphid resequencing samples, with variable diversity in symbiont communities, in order to recover the genome sequence of its obligatory bacterial symbiont, Buchnera aphidicola. It was able to return high-quality assemblies (one contig assembly in 90% of the samples), even when using increasingly distant reference genomes, and to retrieve large structural variations in the samples. Because of its targeted essence, it outperformed standard metagenomic assemblers in terms of both time and assembly quality.

MetaBCC-LR: metagenomics binning by coverage and composition for long reads

MOTIVATION

Metagenomics studies have provided key insights into the composition and structure of microbial communities found in different environments. Among the techniques used to analyse metagenomic data, binning is considered a crucial step to characterize the different species of micro-organisms present. The use of short-read data in most binning tools poses several limitations, such as insufficient species-specific signal, and the emergence of long-read sequencing technologies offers us opportunities to surmount them. However, most current metagenomic binning tools have been developed for short reads. The few tools that can process long reads either do not scale with increasing input size or require a database with reference genomes that are often unknown. In this article, we present MetaBCC-LR, a scalable reference-free binning method which clusters long reads directly based on their k-mer coverage histograms and oligonucleotide composition.

RESULTS

We evaluate MetaBCC-LR on multiple simulated and real metagenomic long-read datasets with varying coverages and error rates. Our experiments demonstrate that MetaBCC-LR substantially outperforms state-of-the-art reference-free binning tools, achieving ∼13% improvement in F1-score and ∼30% improvement in ARI compared to the best previous tools. Moreover, we show that using MetaBCC-LR before long-read assembly helps to enhance the assembly quality while significantly reducing the assembly cost in terms of time and memory usage. The efficiency and accuracy of MetaBCC-LR pave the way for more effective long-read-based metagenomics analyses to support a wide range of applications.

AVAILABILITY AND IMPLEMENTATION

The source code is freely available at: https://github.com/anuradhawick/MetaBCC-LR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Secondary Metabolism in the Gill Microbiota of Shipworms (Teredinidae) as Revealed by Comparison of Metagenomes and Nearly Complete Symbiont Genomes

We define a system in which the major symbionts that are important to host biology and to the production of secondary metabolites can be cultivated. We show that symbiotic bacteria that are critical to host nutrition and lifestyle also have an immense capacity to produce a multitude of diverse and likely novel bioactive secondary metabolites that could lead to the discovery of drugs and that these pathways are found within shipworm gills. We propose that, by shaping associated microbial communities within the host, the compounds support the ability of shipworms to degrade wood in marine environments. Because these symbionts can be cultivated and genetically manipulated, they provide a powerful model for understanding how secondary metabolism impacts microbial symbiosis.

Shipworms play critical roles in recycling wood in the sea. Symbiotic bacteria supply enzymes that the organisms need for nutrition and wood degradation. Some of these bacteria have been grown in pure culture and have the capacity to make many secondary metabolites. However, little is known about whether such secondary metabolite pathways are represented in the symbiont communities within their hosts. In addition, little has been reported about the patterns of host-symbiont co-occurrence. Here, we collected shipworms from the United States, the Philippines, and Brazil and cultivated symbiotic bacteria from their gills. We analyzed sequences from 22 shipworm gill metagenomes from seven shipworm species and from 23 cultivated symbiont isolates. Using (meta)genome sequencing, we demonstrate that the cultivated isolates represent all the major bacterial symbiont species and strains in shipworm gills. We show that the bacterial symbionts are distributed among shipworm hosts in consistent, predictable patterns. The symbiotic bacteria harbor many gene cluster families (GCFs) for biosynthesis of bioactive secondary metabolites, only <5% of which match previously described biosynthetic pathways. Because we were able to cultivate the symbionts and to sequence their genomes, we can definitively enumerate the biosynthetic pathways in these symbiont communities, showing that ∼150 of ∼200 total biosynthetic gene clusters (BGCs) present in the animal gill metagenomes are represented in our culture collection. Shipworm symbionts occur in suites that differ predictably across a wide taxonomic and geographic range of host species and collectively constitute an immense resource for the discovery of new biosynthetic pathways corresponding to bioactive secondary metabolites.

IMPORTANCE We define a system in which the major symbionts that are important to host biology and to the production of secondary metabolites can be cultivated. We show that symbiotic bacteria that are critical to host nutrition and lifestyle also have an immense capacity to produce a multitude of diverse and likely novel bioactive secondary metabolites that could lead to the discovery of drugs and that these pathways are found within shipworm gills. We propose that, by shaping associated microbial communities within the host, the compounds support the ability of shipworms to degrade wood in marine environments. Because these symbionts can be cultivated and genetically manipulated, they provide a powerful model for understanding how secondary metabolism impacts microbial symbiosis.

Genetic and Biochemical Reconstitution of Bromoform Biosynthesis in Asparagopsis Lends Insights into Seaweed Reactive Oxygen Species Enzymology

Marine macroalgae, seaweeds, are exceptionally prolific producers of halogenated natural products. Biosynthesis of halogenated molecules in seaweeds is inextricably linked to reactive oxygen species (ROS) signaling as hydrogen peroxide serves as a substrate for haloperoxidase enzymes that participate in the construction these halogenated molecules. Here, using red macroalga Asparagopsis taxiformis, a prolific producer of the ozone depleting molecule bromoform, we provide the discovery and biochemical characterization of a ROS-producing NAD(P)H oxidase from seaweeds. This discovery was enabled by our sequencing of Asparagopsis genomes, in which we find the gene encoding the ROS-producing enzyme to be clustered with genes encoding bromoform-producing haloperoxidases. Biochemical reconstitution of haloperoxidase activities establishes that fatty acid biosynthesis can provide viable hydrocarbon substrates for bromoform production. The ROS production haloperoxidase enzymology that we describe here advances seaweed biology and biochemistry by providing the molecular basis for decades worth of physiological observations in ROS and halogenated natural product biosyntheses.

Metagenomic Exploration of the Marine Sponge Mycale hentscheli Uncovers Multiple Polyketide-Producing Bacterial Symbionts

Marine sponges have been a prolific source of unique bioactive compounds that are presumed to act as a deterrent to predation. Many of these compounds have potential therapeutic applications; however, the lack of efficient and sustainable synthetic routes frequently limits clinical development. Here, we describe a metagenomic investigation of Mycale hentscheli, a chemically gifted marine sponge that possesses multiple distinct chemotypes. We applied shotgun metagenomic sequencing, hybrid assembly of short- and long-read data, and metagenomic binning to obtain a comprehensive picture of the microbiome of five specimens, spanning three chemotypes. Our data revealed multiple producing species, each having relatively modest secondary metabolomes, that contribute collectively to the chemical arsenal of the holobiont. We assembled complete genomes for multiple new genera, including two species that produce the cytotoxic polyketides pateamine and mycalamide, as well as a third high-abundance symbiont harboring a proteusin-type biosynthetic pathway that appears to encode a new polytheonamide-like compound. We also identified an additional 188 biosynthetic gene clusters, including a pathway for biosynthesis of peloruside. These results suggest that multiple species cooperatively contribute to defensive symbiosis in M. hentscheli and reveal that the taxonomic diversity of secondary-metabolite-producing sponge symbionts is larger and richer than previously recognized.IMPORTANCE Mycale hentscheli is a marine sponge that is rich in bioactive small molecules. Here, we use direct metagenomic sequencing to elucidate highly complete and contiguous genomes for the major symbiotic bacteria of this sponge. We identify complete biosynthetic pathways for the three potent cytotoxic polyketides which have previously been isolated from M. hentscheli Remarkably, and in contrast to previous studies of marine sponges, we attribute each of these metabolites to a different producing microbe. We also find that the microbiome of M. hentscheli is stably maintained among individuals, even over long periods of time. Collectively, our data suggest a cooperative mode of defensive symbiosis in which multiple symbiotic bacterial species cooperatively contribute to the defensive chemical arsenal of the holobiont.

Horizontal Gene Transfer to a Defensive Symbiont with a Reduced Genome in a Multipartite Beetle Microbiome

Symbiotic mutualisms of bacteria and animals are ubiquitous in nature, running a continuum from facultative to obligate from the perspectives of both partners. The loss of functions required for living independently but not within a host gives rise to reduced genomes in many symbionts. Although the phenomenon of genome reduction can be explained by existing evolutionary models, the initiation of the process is not well understood. Here, we describe the microbiome associated with the eggs of the beetle Lagria villosa, consisting of multiple bacterial symbionts related to Burkholderia gladioli, including a reduced-genome symbiont thought to be the exclusive producer of the defensive compound lagriamide. We show that the putative lagriamide-producing symbiont is the only member of the microbiome undergoing genome reduction and that it has already lost the majority of its primary metabolism and DNA repair pathways. The key step preceding genome reduction in the symbiont was likely the horizontal acquisition of the putative lagriamide lga biosynthetic gene cluster. Unexpectedly, we uncovered evidence of additional horizontal transfers to the symbiont's genome while genome reduction was occurring and despite a current lack of genes needed for homologous recombination. These gene gains may have given the genome-reduced symbiont a selective advantage in the microbiome, especially given the maintenance of the large lga gene cluster despite ongoing genome reduction.IMPORTANCE Associations between microorganisms and an animal, plant, or fungal host can result in increased dependence over time. This process is due partly to the bacterium not needing to produce nutrients that the host provides, leading to loss of genes that it would need to live independently and to a consequent reduction in genome size. It is often thought that genome reduction is aided by genetic isolation-bacteria that live in monocultures in special host organs, or inside host cells, have less access to other bacterial species from which they can obtain genes. Here, we describe exposure of a genome-reduced beetle symbiont to a community of related bacteria with nonreduced genomes. We show that the symbiont has acquired genes from other bacteria despite going through genome reduction, suggesting that isolation has not yet played a major role in this case of genome reduction, with horizontal gene gains still offering a potential route for adaptation.

Trophic Status Is Associated With Community Structure and Metabolic Potential of Planktonic Microbiota in Plateau Lakes

Microbes in various aquatic ecosystems play a key role in global energy fluxes and biogeochemical processes. However, the detailed patterns on the functional structure and the metabolic potential of microbial communities in freshwater lakes with different trophic status remain to be understood. We employed a metagenomics workflow to analyze the correlations between trophic status and planktonic microbiota in freshwater lakes on Yun-Gui Plateau, China. Our results revealed that microbial communities in the eutrophic and mesotrophic-oligotrophic lake ecosystems harbor distinct community structure and metabolic potential. Cyanobacteria were dominant in the eutrophic ecosystems, mainly driving the processes of aerobic respiration, fermentation, nitrogen assimilation, nitrogen mineralization, assimilatory sulfate reduction and sulfur mineralization in this ecosystem group. Actinobacteria, Proteobacteria (Alpha-, Beta-, and Gammaproteobacteria), Verrucomicrobia and Planctomycetes, occurred more often in the mesotrophic-oligotrophic ecosystems than those in the eutrophic ecosystems, and these taxa potentially mediate the above metabolic processes. In these two groups of ecosystems, a difference in the abundance of functional genes involved in carbohydrate metabolism, energy metabolism, glycan biosynthesis and metabolism, and metabolism of cofactors and vitamins significantly contribute to the distinct functional structure of microbiota from surface water. Furthermore, the microbe-mediated metabolic potentials for carbon, nitrogen and sulfur transformation showed differences in the two ecosystem groups. Compared with the mesotrophic-oligotrophic ecosystems, planktonic microbial communities in the eutrophic ecosystems showed higher potential for aerobic carbon fixation, fermentation, methanogenesis, anammox, denitrification, and sulfur mineralization, but they showed lower potential for aerobic respiration, CO oxidation, nitrogen fixation, and assimilatory sulfate reduction. This study offers insights into the relationships of trophic status to planktonic microbial community structure and its metabolic potential, and identifies the main taxa responsible for the biogeochemical cycles of carbon, nitrogen and sulfur in freshwater lake environments.