Single sample resolution of rare microbial dark matter in a marine invertebrate metagenome

SoxY gene family expansion underpins adaptation to diverse hosts and environments in symbiotic sulfide oxidizers

Sulfur-oxidizing bacteria (SOB) have developed distinct ecological strategies to obtain reduced sulfur compounds for growth. These range from specialists that can only use a limited range of reduced sulfur compounds to generalists that can use many different forms as electron donors. Forming intimate symbioses with animal hosts is another highly successful ecological strategy for SOB, as animals, through their behavior and physiology, can enable access to sulfur compounds. Symbioses have evolved multiple times in a range of animal hosts and from several lineages of SOB. They have successfully colonized a wide range of habitats, from seagrass beds to hydrothermal vents, with varying availability of symbiont energy sources. Our extensive analyses of sulfur transformation pathways in 234 genomes of symbiotic and free-living SOB revealed widespread conservation in metabolic pathways for sulfur oxidation in symbionts from different host species and environments, raising the question of how they have adapted to such a wide range of distinct habitats. We discovered a gene family expansion of soxY in these genomes, with up to five distinct copies per genome. Symbionts harboring only the "canonical" soxY were typically ecological "specialists" that are associated with specific host subfamilies or environments (e.g., hydrothermal vents, mangroves). Conversely, symbionts with multiple divergent soxY genes formed versatile associations across diverse hosts in various marine environments. We hypothesize that expansion and diversification of the soxY gene family could be one genomic mechanism supporting the metabolic flexibility of symbiotic SOB enabling them and their hosts to thrive in a range of different and dynamic environments.IMPORTANCESulfur metabolism is thought to be one of the most ancient mechanisms for energy generation in microorganisms. A diverse range of microorganisms today rely on sulfur oxidation for their metabolism. They can be free-living, or they can live in symbiosis with animal hosts, where they power entire ecosystems in the absence of light, such as in the deep sea. In the millions of years since they evolved, sulfur-oxidizing bacteria have adopted several highly successful strategies; some are ecological "specialists," and some are "generalists," but which genetic features underpin these ecological strategies are not well understood. We discovered a gene family that has become expanded in those species that also seem to be "generalists," revealing that duplication, repurposing, and reshuffling existing genes can be a powerful mechanism driving ecological lifestyle shifts.

Physiological changes during torpor favor association with Endozoicomonas endosymbionts in the urochordate Botrylloides leachii

Environmental perturbations evoke down-regulation of metabolism in some multicellular organisms, leading to dormancy, or torpor. Colonies of the urochordate Botrylloides leachii enter torpor in response to changes in seawater temperature and may survive for months as small vasculature remnants that lack feeding and reproductive organs but possess torpor-specific microbiota. Upon returning to milder conditions, the colonies rapidly restore their original morphology, cytology and functionality while harboring re-occurring microbiota, a phenomenon that has not been described in detail to date. Here we investigated the stability of B. leachii microbiome and its functionality in active and dormant colonies, using microscopy, qPCR, in situ hybridization, genomics and transcriptomics. A novel lineage of Endozoicomonas, proposed here as Candidatus Endozoicomonas endoleachii, was dominant in torpor animals (53-79% read abundance), and potentially occupied specific hemocytes found only in torpid animals. Functional analysis of the metagenome-assembled genome and genome-targeted transcriptomics revealed that Endozoicomonas can use various cellular substrates, like amino acids and sugars, potentially producing biotin and thiamine, but also expressing various features involved in autocatalytic symbiosis. Our study suggests that the microbiome can be linked to the metabolic and physiological states of the host, B. leachii, introducing a model organism for the study of symbioses during drastic physiological changes, such as torpor.

Bacterial Microbiota of Ostreobium, the Coral-Isolated Chlorophyte Ectosymbiont, at Contrasted Salinities

Microscopic filaments of the siphonous green algae Ostreobium (Ulvophyceae, Bryopsidales) colonize and dissolve the calcium carbonate skeletons of coral colonies in reefs of contrasted salinities. Here, we analyzed their bacterial community's composition and plasticity in response to salinity. Multiple cultures of Pocillopora coral-isolated Ostreobium strains from two distinct rbcL lineages representative of IndoPacific environmental phylotypes were pre-acclimatized (>9 months) to three ecologically relevant reef salinities: 32.9, 35.1, and 40.2 psu. Bacterial phylotypes were visualized for the first time at filament scale by CARD-FISH in algal tissue sections, within siphons, at their surface or in their mucilage. Ostreobium-associated microbiota, characterized by bacterial 16S rDNA metabarcoding of cultured thalli and their corresponding supernatants, were structured by host genotype (Ostreobium strain lineage), with dominant Kiloniellaceae or Rhodospirillaceae (Alphaproteobacteria, Rhodospirillales) depending on Ostreobium lineage, and shifted Rhizobiales' abundances in response to the salinity increase. A small core microbiota composed of seven ASVs (~1.5% of thalli ASVs, 19-36% cumulated proportions) was persistent across three salinities in both genotypes, with putative intracellular Amoebophilaceae and Rickettsiales_AB1, as well as Hyphomonadaceae and Rhodospirillaceae also detected within environmental (Ostreobium-colonized) Pocillopora coral skeletons. This novel knowledge on the taxonomic diversity of Ostreobium bacteria paves the way to functional interaction studies within the coral holobiont.

Hyperactive nanobacteria with host-dependent traits pervade Omnitrophota

Candidate bacterial phylum Omnitrophota has not been isolated and is poorly understood. We analysed 72 newly sequenced and 349 existing Omnitrophota genomes representing 6 classes and 276 species, along with Earth Microbiome Project data to evaluate habitat, metabolic traits and lifestyles. We applied fluorescence-activated cell sorting and differential size filtration, and showed that most Omnitrophota are ultra-small (~0.2 μm) cells that are found in water, sediments and soils. Omnitrophota genomes in 6 classes are reduced, but maintain major biosynthetic and energy conservation pathways, including acetogenesis (with or without the Wood-Ljungdahl pathway) and diverse respirations. At least 64% of Omnitrophota genomes encode gene clusters typical of bacterial symbionts, suggesting host-associated lifestyles. We repurposed quantitative stable-isotope probing data from soils dominated by andesite, basalt or granite weathering and identified 3 families with high isotope uptake consistent with obligate bacterial predators. We propose that most Omnitrophota inhabit various ecosystems as predators or parasites.

Seasonal dynamics of a complex cheilostome bryozoan symbiosis: vertical transfer challenged

Symbiotic associations are dynamic systems influenced by both intrinsic and extrinsic factors. Here we describe for the first time the developmental and seasonal changes of the funicular bodies in the bryozoan Dendrobeania fruticosa, which are unique temporary organs of cheilostome bryozoans containing prokaryotic symbionts. Histological and ultrastructural studies showed that these organs undergo strong seasonal modification in the White Sea during the ice-free period. Initially (in June) they play a trophic function and support the development of a large population of bacteria. From June to September, both funicular bodies and bacteria show signs of degradation accompanied by development of presumed virus-like particles (VLPs); these self-organize to hollow spheres inside bacteria and are also detected outside of them. Although the destruction of bacteria coincides with the development of VLPs and spheres, the general picture differs considerably from the known instances of bacteriophagy in bryozoans. We broadly discuss potential routes of bacterial infection in Bryozoa and question the hypothesis of vertical transfer, which, although widely accepted in the literature, is contradicted by molecular, morphological and ecological evidence.

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.

Functional characterization of prokaryotic dark matter: the road so far and what lies ahead

Eight-hundred thousand to one trillion prokaryotic species may inhabit our planet. Yet, fewer than two-hundred thousand prokaryotic species have been described. This uncharted fraction of microbial diversity, and its undisclosed coding potential, is known as the "microbial dark matter" (MDM). Next-generation sequencing has allowed to collect a massive amount of genome sequence data, leading to unprecedented advances in the field of genomics. Still, harnessing new functional information from the genomes of uncultured prokaryotes is often limited by standard classification methods. These methods often rely on sequence similarity searches against reference genomes from cultured species. This hinders the discovery of unique genetic elements that are missing from the cultivated realm. It also contributes to the accumulation of prokaryotic gene products of unknown function among public sequence data repositories, highlighting the need for new approaches for sequencing data analysis and classification. Increasing evidence indicates that these proteins of unknown function might be a treasure trove of biotechnological potential. Here, we outline the challenges, opportunities, and the potential hidden within the functional dark matter (FDM) of prokaryotes. We also discuss the pitfalls surrounding molecular and computational approaches currently used to probe these uncharted waters, and discuss future opportunities for research and applications.

Coral holobiont cues prime Endozoicomonas for a symbiotic lifestyle

Endozoicomonas are prevalent, abundant bacterial associates of marine animals, including corals. Their role in holobiont health and functioning, however, remains poorly understood. To identify potential interactions within the coral holobiont, we characterized the novel isolate Endozoicomonas marisrubri sp. nov. 6c and assessed its transcriptomic and proteomic response to tissue extracts of its native host, the Red Sea coral Acropora humilis. We show that coral tissue extracts stimulated differential expression of genes putatively involved in symbiosis establishment via the modulation of the host immune response by E. marisrubri 6c, such as genes for flagellar assembly, ankyrins, ephrins, and serpins. Proteome analyses revealed that E. marisrubri 6c upregulated vitamin B1 and B6 biosynthesis and glycolytic processes in response to holobiont cues. Our results suggest that the priming of Endozoicomonas for a symbiotic lifestyle involves the modulation of host immunity and the exchange of essential metabolites with other holobiont members. Consequently, Endozoicomonas may play an important role in holobiont nutrient cycling and may therefore contribute to coral health, acclimatization, and adaptation.

Non-essential ribosomal proteins in bacteria and archaea identified using COGs

Ribosomal proteins (RPs) are highly conserved across the bacterial and archaeal domains. Although many RPs are essential for survival, genome analysis demonstrates the absence of some RP genes in many bacterial and archaeal genomes. Furthermore, global transposon mutagenesis and/or targeted deletion showed that elimination of some RP genes had only a moderate effect on the bacterial growth rate. Here, we systematically analyze the evolutionary conservation of RPs in prokaryotes by compiling the list of the ribosomal genes that are missing from one or more genomes in the recently updated version of the Clusters of Orthologous Genes (COG) database. Some of these absences occurred because the respective genes carried frameshifts, presumably, resulting from sequencing errors, while others were overlooked and not translated during genome annotation. Apart from these annotation errors, we identified multiple genuine losses of RP genes in a variety of bacteria and archaea. Some of these losses are clade-specific, whereas others occur in symbionts and parasites with dramatically reduced genomes. The lists of computationally and experimentally defined non-essential ribosomal genes show a substantial overlap, revealing a common trend in prokaryote ribosome evolution that could be linked to the architecture and assembly of the ribosomes. Thus, RPs that are located at the surface of the ribosome and/or are incorporated at a late stage of ribosome assembly are more likely to be non-essential and to be lost during microbial evolution, particularly, in the course of genome compaction.IMPORTANCEIn many prokaryote genomes, one or more ribosomal protein (RP) genes are missing. Analysis of 1,309 prokaryote genomes included in the COG database shows that only about half of the RPs are universally conserved in bacteria and archaea. In contrast, up to 16 other RPs are missing in some genomes, primarily, tiny (<1 Mb) genomes of host-associated bacteria and archaea. Ten universal and nine archaea-specific ribosomal proteins show clear patterns of lineage-specific gene loss. Most of the RPs that are frequently lost from bacterial genomes are located on the ribosome periphery and are non-essential in Escherichia coli and Bacillus subtilis These results reveal general trends and common constraints in the architecture and evolution of ribosomes in prokaryotes.

First evidence of virus-like particles in the bacterial symbionts of Bryozoa

Bacteriophage communities associated with humans and vertebrate animals have been extensively studied, but the data on phages living in invertebrates remain scarce. In fact, they have never been reported for most animal phyla. Our ultrastructural study showed for the first time a variety of virus-like particles (VLPs) and supposed virus-related structures inside symbiotic bacteria in two marine species from the phylum Bryozoa, the cheilostomes Bugula neritina and Paralicornia sinuosa. We also documented the effect of VLPs on bacterial hosts: we explain different bacterial 'ultrastructural types' detected in bryozoan tissues as stages in the gradual destruction of prokaryotic cells caused by viral multiplication during the lytic cycle. We speculate that viruses destroying bacteria regulate symbiont numbers in the bryozoan hosts, a phenomenon known in some insects. We develop two hypotheses explaining exo- and endogenous circulation of the viruses during the life-cycle of B. neritina. Finally, we compare unusual 'sea-urchin'-like structures found in the collapsed bacteria in P. sinuosa with so-called metamorphosis associated contractile structures (MACs) formed in the cells of the marine bacterium Pseudoalteromonas luteoviolacea which are known to trigger larval metamorphosis in a polychaete worm.

Synthesizing Systems Biology Knowledge from Omics Using Genome-Scale Models

Omic technologies have enabled the complete readout of the molecular state of a cell at different biological scales. In principle, the combination of multiple omic data types can provide an integrated view of the entire biological system. This integration requires appropriate models in a systems biology approach. Here, genome-scale models (GEMs) are focused upon as one computational systems biology approach for interpreting and integrating multi-omic data. GEMs convert the reactions (related to metabolism, transcription, and translation) that occur in an organism to a mathematical formulation that can be modeled using optimization principles. A variety of genome-scale modeling methods used to interpret multiple omic data types, including genomics, transcriptomics, proteomics, metabolomics, and meta-omics are reviewed. The ability to interpret omics in the context of biological systems has yielded important findings for human health, environmental biotechnology, bioenergy, and metabolic engineering. The authors find that concurrent with advancements in omic technologies, genome-scale modeling methods are also expanding to enable better interpretation of omic data. Therefore, continued synthesis of valuable knowledge, through the integration of omic data with GEMs, are expected.

Accurate and complete genomes from metagenomes

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

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.

Computational Modeling of the Human Microbiome

The impact of microorganisms on human health has long been acknowledged and studied, but recent advances in research methodologies have enabled a new systems-level perspective on the collections of microorganisms associated with humans, the human microbiome. Large-scale collaborative efforts such as the NIH Human Microbiome Project have sought to kick-start research on the human microbiome by providing foundational information on microbial composition based upon specific sites across the human body. Here, we focus on the four main anatomical sites of the human microbiome: gut, oral, skin, and vaginal, and provide information on site-specific background, experimental data, and computational modeling. Each of the site-specific microbiomes has unique organisms and phenomena associated with them; there are also high-level commonalities. By providing an overview of different human microbiome sites, we hope to provide a perspective where detailed, site-specific research is needed to understand causal phenomena that impact human health, but there is equally a need for more generalized methodology improvements that would benefit all human microbiome research.

Autometa: automated extraction of microbial genomes from individual shotgun metagenomes

Shotgun metagenomics is a powerful, high-resolution technique enabling the study of microbial communities in situ. However, species-level resolution is only achieved after a process of 'binning' where contigs predicted to originate from the same genome are clustered. Such culture-independent sequencing frequently unearths novel microbes, and so various methods have been devised for reference-free binning. As novel microbiomes of increasing complexity are explored, sometimes associated with non-model hosts, robust automated binning methods are required. Existing methods struggle with eukaryotic contamination and cannot handle highly complex single metagenomes. We therefore developed an automated binning pipeline, termed 'Autometa', to address these issues. This command-line application integrates sequence homology, nucleotide composition, coverage and the presence of single-copy marker genes to separate microbial genomes from non-model host genomes and other eukaryotic contaminants, before deconvoluting individual genomes from single metagenomes. The method is able to effectively separate over 1000 genomes from a metagenome, allowing the study of previously intractably complex environments at the level of single species. Autometa is freely available at https://bitbucket.org/jason_c_kwan/autometa and as a docker image at https://hub.docker.com/r/jasonkwan/autometa under the GNU Affero General Public License 3 (AGPL 3).

Diversity and Dynamics of "Candidatus Endobugula" and Other Symbiotic Bacteria in Chinese Populations of the Bryozoan, Bugula neritina

Bugula neritina is a common invasive cosmopolitan bryozoan that harbors (like many sessile marine invertebrates) a symbiotic bacterial (SB) community. Among the SB of B. neritina, "Candidatus Endobugula sertula" continues to receive the greatest attention, because it is the source of bryostatins. The bryostatins are potent bioactive polyketides, which have been investigated for their therapeutic potential to treat various cancers, Alzheimer's disease, and AIDS. In this study, we compare the metagenomics sequences for the 16S ribosomal RNA gene of the SB communities from different geographic and life cycle samples of Chinese B. neritina. Using a variety of approaches for estimating alpha/beta diversity and taxonomic abundance, we find that the SB communities vary geographically with invertebrate and fish mariculture and with latitude and environmental temperature. During the B. neritina life cycle, we find that the diversity and taxonomic abundances of the SB communities change with the onset of host metamorphosis, filter feeding, colony formation, reproduction, and increased bryostatin production. "Ca. Endobugula sertula" is confirmed as the symbiont of the Chinese "Ca. Endobugula"/B. neritina symbiosis. Our study extends our knowledge about B. neritina symbiosis from the New to the Old World and offers new insights into the environmental and life cycle factors that can influence its SB communities, "Ca. Endobugula," and bryostatins more globally.

Bacterial spores, from ecology to biotechnology

The production of a highly specialized cell structure called a spore is a remarkable example of a survival strategy displayed by bacteria in response to challenging environmental conditions. The detailed analysis and description of the process of sporulation in selected model organisms have generated a solid background to understand the cellular processes leading to the formation of this specialized cell. However, much less is known regarding the ecology of spore-formers. This research gap needs to be filled as the feature of resistance has important implications not only on the survival of spore-formers and their ecology, but also on the use of spores for environmental prospection and biotechnological applications.

Ca. Endozoicomonas cretensis: A Novel Fish Pathogen Characterized by Genome Plasticity

Endozoicomonas bacteria are generally beneficial symbionts of diverse marine invertebrates including reef-building corals, sponges, sea squirts, sea slugs, molluscs, and Bryozoans. In contrast, the recently reported Ca. Endozoicomonas cretensis was identified as a vertebrate pathogen, causing epitheliocystis in fish larvae resulting in massive mortality. Here, we described the Ca. E. cretensis draft genome, currently undergoing genome decay as evidenced by massive insertion sequence (IS element) expansion and pseudogene formation. Many of the insertion sequences are also predicted to carry outward-directed promoters, implying that they may be able to modulate the expression of neighbouring coding sequences (CDSs). Comparative genomic analysis has revealed many Ca. E. cretensis-specific CDSs, phage integration and novel gene families. Potential virulence related CDSs and machineries were identified in the genome, including secretion systems and related effector proteins, and systems related to biofilm formation and directed cell movement. Mucin degradation would be of importance to a fish pathogen, and many candidate CDSs associated with this pathway have been identified. The genome may reflect a bacterium in the process of changing niche from symbiont to pathogen, through expansion of virulence genes and some loss of metabolic capacity.

Microbial-derived products as potential new antimicrobials

Due to the continuing global concerns involving antibiotic resistance, there is a need for scientific forums to assess advancements in the development of antimicrobials and their alternatives that might reduce development and spread of antibiotic resistance among bacterial pathogens. The objectives of the 2nd International Symposium on Alternatives to Antibiotics were to highlight promising research results and novel technologies that can provide alternatives to antibiotics for use in animal health and production, assess challenges associated with their authorization and commercialization for use, and provide actionable strategies to support their development. The session on microbial-derived products was directed at presenting novel technologies that included exploiting CRISPR-Cas nucleases to produce sequence-specific antimicrobials, probiotics development via fecal microbiome transplants among monogastric production animals such as chickens and mining microbial sources such as bacteria or yeast to identify new antimicrobial compounds. Other research has included continuing development of antimicrobial peptides such as newly discovered bacteriocins as alternatives to antibiotics, use of bacteriophages accompanied by development of unique lytic proteins with specific cell-wall binding domains and novel approaches such as microbial-ecology guided discovery of anti-biofilm compounds discovered in marine environments. The symposium was held at the Headquarters of the World Organisation for Animal Health (OIE) in Paris, France during 12-15 December 2016.

Chemoenzymatic Dissection of Polyketide β-Branching in the Bryostatin Pathway

β-Branching is an expansion upon canonical polyketide synthase extension that allows for the installation of diverse chemical moieties in several natural products. Several of these moieties are unique among natural products, including the two vinyl methylesters found in the core structure of bryostatins. This family of molecules is derived from an obligate bacterial symbiont of a sessile marine bryozoan, Bugula neritina. Within this family, bryostatin 1 has been investigated as an anticancer, neuroprotective, and immunomodulatory compound. We have turned to the biosynthetic gene cluster within the bacterial symbiont to investigate the biosynthesis of bryostatins. Recent sequencing efforts resulted in the annotation of two missing genes: bryT and bryU. Using novel chemoenzymatic techniques, we have validated these as the missing enoyl-CoA hydratase and donor acyl carrier protein, essential components of the β-branching cassette of the bryostatin pathway. Together, this cassette installs the vinyl methylester moieties essential to the activity of bryostatins.

The Who, Why, and How of Small-Molecule Production in Invertebrate Microbiomes: Basic Insights Fueling Drug Discovery

Bacteria have supplied us with many bioactive molecules for use in medicine and agriculture. However, rates of discovery have decreased as the biosynthetic capacity of the culturable biosphere has been continuously mined for many decades.

Bacteria have supplied us with many bioactive molecules for use in medicine and agriculture. However, rates of discovery have decreased as the biosynthetic capacity of the culturable biosphere has been continuously mined for many decades. The as-yet-uncultured biosphere is likely to hold far greater biosynthetic potential, especially where ecological niches favor the selection of therapeutically useful bioactivities. I outline here how metagenomics and other systems biology approaches can be used to gain insight into small-molecule biosynthesis and the selective forces which shape it. I also argue that we need a greater understanding of the function of small molecules in complex microbiomes and rational synthetic biology methods to functionally reconstruct large biosynthetic pathways in heterologous hosts.

Increased Biosynthetic Gene Dosage in a Genome-Reduced Defensive Bacterial Symbiont

Secondary metabolites, which are small-molecule organic compounds produced by living organisms, provide or inspire drugs for many different diseases. These natural products have evolved over millions of years to provide a survival benefit to the producing organism and often display potent biological activity with important therapeutic applications. For instance, defensive compounds in the environment may be cytotoxic to eukaryotic cells, a property exploitable for cancer treatment. Here, we describe the genome of an uncultured symbiotic bacterium that makes such a cytotoxic metabolite. This symbiont is losing genes that do not endow a selective advantage in a hospitable host environment. Secondary metabolism genes, however, are repeated multiple times in the genome, directly demonstrating their selective advantage. This finding shows the strength of selective forces in symbiotic relationships and suggests that uncultured bacteria in such relationships should be targeted for drug discovery efforts.

A symbiotic lifestyle frequently results in genome reduction in bacteria; the isolation of small populations promotes genetic drift and the fixation of deletions and deleterious mutations over time. Transitions in lifestyle, including host restriction or adaptation to an intracellular habitat, are thought to precipitate a wave of sequence degradation events and consequent proliferation of pseudogenes. We describe here a verrucomicrobial symbiont of the tunicate Lissoclinum sp. that appears to be undergoing such a transition, with low coding density and many identifiable pseudogenes. However, despite the overall drive toward genome reduction, this symbiont maintains seven copies of a large polyketide synthase (PKS) pathway for the mandelalides (mnd), cytotoxic compounds that likely constitute a chemical defense for the host. There is evidence of ongoing degradation in a small number of these repeats—including variable borders, internal deletions, and single nucleotide polymorphisms (SNPs). However, the gene dosage of most of the pathway is increased at least 5-fold. Correspondingly, this single pathway accounts for 19% of the genome by length and 25.8% of the coding capacity. This increased gene dosage in the face of generalized sequence degradation and genome reduction suggests that mnd genes are under strong purifying selection and are important to the symbiotic relationship.

IMPORTANCE Secondary metabolites, which are small-molecule organic compounds produced by living organisms, provide or inspire drugs for many different diseases. These natural products have evolved over millions of years to provide a survival benefit to the producing organism and often display potent biological activity with important therapeutic applications. For instance, defensive compounds in the environment may be cytotoxic to eukaryotic cells, a property exploitable for cancer treatment. Here, we describe the genome of an uncultured symbiotic bacterium that makes such a cytotoxic metabolite. This symbiont is losing genes that do not endow a selective advantage in a hospitable host environment. Secondary metabolism genes, however, are repeated multiple times in the genome, directly demonstrating their selective advantage. This finding shows the strength of selective forces in symbiotic relationships and suggests that uncultured bacteria in such relationships should be targeted for drug discovery efforts.

Author Video: An author video summary of this article is available.

Interpreting Microbial Biosynthesis in the Genomic Age: Biological and Practical Considerations

Genome mining has become an increasingly powerful, scalable, and economically accessible tool for the study of natural product biosynthesis and drug discovery. However, there remain important biological and practical problems that can complicate or obscure biosynthetic analysis in genomic and metagenomic sequencing projects. Here, we focus on limitations of available technology as well as computational and experimental strategies to overcome them. We review the unique challenges and approaches in the study of symbiotic and uncultured systems, as well as those associated with biosynthetic gene cluster (BGC) assembly and product prediction. Finally, to explore sequencing parameters that affect the recovery and contiguity of large and repetitive BGCs assembled de novo, we simulate Illumina and PacBio sequencing of the Salinispora tropica genome focusing on assembly of the salinilactam (slm) BGC.

metaSPAdes: a new versatile metagenomic assembler

While metagenomics has emerged as a technology of choice for analyzing bacterial populations, the assembly of metagenomic data remains challenging, thus stifling biological discoveries. Moreover, recent studies revealed that complex bacterial populations may be composed from dozens of related strains, thus further amplifying the challenge of metagenomic assembly. metaSPAdes addresses various challenges of metagenomic assembly by capitalizing on computational ideas that proved to be useful in assemblies of single cells and highly polymorphic diploid genomes. We benchmark metaSPAdes against other state-of-the-art metagenome assemblers and demonstrate that it results in high-quality assemblies across diverse data sets.

Lack of Overt Genome Reduction in the Bryostatin-Producing Bryozoan Symbiont "Candidatus Endobugula sertula"

The uncultured bacterial symbiont "Candidatus Endobugula sertula" is known to produce cytotoxic compounds called bryostatins, which protect the larvae of its host, Bugula neritina The symbiont has never been successfully cultured, and it was thought that its genome might be significantly reduced. Here, we took a shotgun metagenomics and metatranscriptomics approach to assemble and characterize the genome of "Ca Endobugula sertula." We found that it had specific metabolic deficiencies in the biosynthesis of certain amino acids but few other signs of genome degradation, such as small size, abundant pseudogenes, and low coding density. We also identified homologs to genes associated with insect pathogenesis in other gammaproteobacteria, and these genes may be involved in host-symbiont interactions and vertical transmission. Metatranscriptomics revealed that these genes were highly expressed in a reproductive host, along with bry genes for the biosynthesis of bryostatins. We identified two new putative bry genes fragmented from the main bry operon, accounting for previously missing enzymatic functions in the pathway. We also determined that a gene previously assigned to the pathway, bryS, is not expressed in reproductive tissue, suggesting that it is not involved in the production of bryostatins. Our findings suggest that "Ca Endobugula sertula" may be able to live outside the host if its metabolic deficiencies are alleviated by medium components, which is consistent with recent findings that it may be possible for "Ca Endobugula sertula" to be transmitted horizontally.

IMPORTANCE

The bryostatins are potent protein kinase C activators that have been evaluated in clinical trials for a number of indications, including cancer and Alzheimer's disease. There is, therefore, considerable interest in securing a renewable supply of these compounds, which is currently only possible through aquaculture of Bugula neritina and total chemical synthesis. However, these approaches are labor-intensive and low-yielding and thus preclude the use of bryostatins as a viable therapeutic agent. Our genome assembly and transcriptome analysis for "Ca Endobugula sertula" shed light on the metabolism of this symbiont, potentially aiding isolation and culturing efforts. Our identification of additional bry genes may also facilitate efforts to express the complete pathway heterologously.