DarkHorse: a method for genome-wide prediction of horizontal gene transfer

Zinc stimulation of phytoplankton in a low carbon dioxide, coastal Antarctic environment

Zinc (Zn) is a key micronutrient used by phytoplankton for carbon (C) acquisition, yet there have been few observations of its influence on natural oceanic phytoplankton populations. In this study, we observed Zn limitation of growth in the natural phytoplankton community of Terra Nova Bay, Antarctica, due to low (~220 μatm) pCO2 conditions, in addition to primary iron (Fe) limitation. Shipboard incubation experiments amended with Zn and Fe resulted in significantly higher chlorophyll a content and dissolved inorganic carbon drawdown compared to Fe addition alone. Zn and Fe response proteins detected in incubation and environmental biomass provided independent verification of algal co-stress for these micronutrients. These observations of Zn limitation under low pCO2 conditions demonstrate Zn can influence coastal primary productivity. Yet, as surface ocean pCO2 rises with continued anthropogenic emissions, the occurrence of Zn/C co-limitation will become rarer, impacting the biogeochemical cycling of Zn and other trace metal micronutrients.

Database Bias in the Detection of Interdomain Horizontal Gene Transfer Events in Pezizomycotina

Horizontal gene transfer (HGT) is a widely acknowledged phenomenon in prokaryotes for generating genetic diversity. However, the impact of this process in eukaryotes, particularly interdomain HGT, is a topic of debate. Although there have been observed biases in interdomain HGT detection, little exploration has been conducted on the effects of imbalanced databases. In our study, we conducted experiments to assess how different databases affect the detection of interdomain HGT using proteomes from the Pezizomycotina fungal subphylum as our focus group. Our objective was to simulate the database imbalance commonly found in public biological databases, where bacterial and eukaryotic sequences are unevenly represented, and demonstrate that an increase in uploaded eukaryotic sequences leads to a decrease in predicted HGTs. For our experiments, four databases with varying proportions of eukaryotic sequences but consistent proportions of bacterial sequences were utilized. We observed a significant reduction in detected interdomain HGT candidates as the proportion of eukaryotes increased within the database. Our data suggest that the imbalance in databases bias the interdomain HGT detection and highlights challenges associated with confirming the presence of interdomain HGT among Pezizomycotina fungi and potentially other groups within Eukarya.

Iron limitation of heterotrophic bacteria in the California Current System tracks relative availability of organic carbon and iron

Iron is an essential nutrient for all microorganisms of the marine environment. Iron limitation of primary production has been well documented across a significant portion of the global surface ocean, but much less is known regarding the potential for iron limitation of the marine heterotrophic microbial community. In this work, we characterize the transcriptomic response of the heterotrophic bacterial community to iron additions in the California Current System, an eastern boundary upwelling system, to detect in situ iron stress of heterotrophic bacteria. Changes in gene expression in response to iron availability by heterotrophic bacteria were detected under conditions of high productivity when carbon limitation was relieved but when iron availability remained low. The ratio of particulate organic carbon to dissolved iron emerged as a biogeochemical proxy for iron limitation of heterotrophic bacteria in this system. Iron stress was characterized by high expression levels of iron transport pathways and decreased expression of iron-containing enzymes involved in carbon metabolism, where a majority of the heterotrophic bacterial iron requirement resides. Expression of iron stress biomarkers, as identified in the iron-addition experiments, was also detected insitu. These results suggest iron availability will impact the processing of organic matter by heterotrophic bacteria with potential consequences for the marine biological carbon pump.

Massive horizontal gene transfer and the evolution of nematomorph-driven behavioral manipulation of mantids

To complete their life cycle, a wide range of parasites must manipulate the behavior of their hosts.1 This manipulation is a well-known example of the "extended phenotype,2" where genes in one organism have phenotypic effects on another organism. Recent studies have explored the parasite genes responsible for such manipulation of host behavior, including the potential molecular mechanisms.3,4 However, little is known about how parasites have acquired the genes involved in manipulating phylogenetically distinct hosts.4 In a fascinating example of the extended phenotype, nematomorph parasites have evolved the ability to induce their terrestrial insect hosts to enter bodies of water, where the parasite then reproduces. Here, we comprehensively analyzed nematomorphs and their mantid hosts, focusing on the transcriptomic changes associated with host manipulations and sequence similarity between host and parasite genes to test molecular mimicry. The nematomorph's transcriptome changed during host manipulation, whereas no distinct changes were found in mantids. We then discovered numerous possible host-derived genes in nematomorphs, and these genes were frequently up-regulated during host manipulation. Our findings suggest a possible general role of horizontal gene transfer (HGT) in the molecular mechanisms of host manipulation, as well as in the genome evolution of manipulative parasites. The evidence of HGT between multicellular eukaryotes remains scarce but is increasing and, therefore, elucidating its mechanisms will advance our understanding of the enduring influence of HGT on the evolution of the web of life.

Short-term acidification promotes diverse iron acquisition and conservation mechanisms in upwelling-associated phytoplankton

Coastal upwelling regions are among the most productive marine ecosystems but may be threatened by amplified ocean acidification. Increased acidification is hypothesized to reduce iron bioavailability for phytoplankton thereby expanding iron limitation and impacting primary production. Here we show from community to molecular levels that phytoplankton in an upwelling region respond to short-term acidification exposure with iron uptake pathways and strategies that reduce cellular iron demand. A combined physiological and multi-omics approach was applied to trace metal clean incubations that introduced 1200 ppm CO2 for up to four days. Although variable, molecular-level responses indicate a prioritization of iron uptake pathways that are less hindered by acidification and reductions in iron utilization. Growth, nutrient uptake, and community compositions remained largely unaffected suggesting that these mechanisms may confer short-term resistance to acidification; however, we speculate that cellular iron demand is only temporarily satisfied, and longer-term acidification exposure without increased iron inputs may result in increased iron stress.

Profiling novel lateral gene transfer events in the human microbiome

Lateral gene transfer (LGT) is an important mechanism for genome diversification in microbial populations, including the human microbiome. While prior work has surveyed LGT events in human-associated microbial isolate genomes, the scope and dynamics of novel LGT events arising in personal microbiomes are not well understood, as there are no widely adopted computational methods to detect, quantify, and characterize LGT from complex microbial communities. We addressed this by developing, benchmarking, and experimentally validating a computational method (WAAFLE) to profile novel LGT events from assembled metagenomes. Applying WAAFLE to >2K human metagenomes from diverse body sites, we identified >100K putative high-confidence but previously uncharacterized LGT events (~2 per assembled microbial genome-equivalent). These events were enriched for mobile elements (as expected), as well as restriction-modification and transport functions typically associated with the destruction of foreign DNA. LGT frequency was quantifiably influenced by biogeography, the phylogenetic similarity of the involved taxa, and the ecological abundance of the donor taxon. These forces manifest as LGT networks in which hub species abundant in a community type donate unequally with their close phylogenetic neighbors. Our findings suggest that LGT may be a more ubiquitous process in the human microbiome than previously described. The open-source WAAFLE implementation, documentation, and data from this work are available at http://huttenhower.sph.harvard.edu/waafle.

Leveraging comparative genomics to uncover alien genes in bacterial genomes

A significant challenge in bacterial genomics is to catalogue genes acquired through the evolutionary process of horizontal gene transfer (HGT). Both comparative genomics and sequence composition-based methods have often been invoked to quantify horizontally acquired genes in bacterial genomes. Comparative genomics methods rely on completely sequenced genomes and therefore the confidence in their predictions increases as the databases become more enriched in completely sequenced genomes. Recent developments including in microbial genome sequencing call for reassessment of alien genes based on information-rich resources currently available. We revisited the comparative genomics approach and developed a new algorithm for alien gene detection. Our algorithm compared favourably with the existing comparative genomics-based methods and is capable of detecting both recent and ancient transfers. It can be used as a standalone tool or in concert with other complementary algorithms for comprehensively cataloguing alien genes in bacterial genomes.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Reduced and Nonreduced Genomes in Paraburkholderia Symbionts of Social Amoebas

The social amoeba Dictyostelium discoideum is a predatory soil protist frequently used for studying host-pathogen interactions. A subset of D. discoideum strains isolated from soil persistently carry symbiotic Paraburkholderia, recently formally described as P. agricolaris, P. bonniea, and P. hayleyella. The three facultative symbiont species of D. discoideum present a unique opportunity to study a naturally occurring symbiosis in a laboratory model protist. There is a large difference in genome size between P. agricolaris (8.7 million base pairs [Mbp]) versus P. hayleyella and P. bonniea (4.1 Mbp). We took a comparative genomics approach and compared the three genomes of D. discoideum symbionts to 12 additional Paraburkholderia genomes to test for genome evolution patterns that frequently accompany host adaptation. Overall, P. agricolaris is difficult to distinguish from other Paraburkholderia based on its genome size and content, but the reduced genomes of P. bonniea and P. hayleyella display characteristics indicative of genome streamlining rather than deterioration during adaptation to their protist hosts. In addition, D. discoideum-symbiont genomes have increased secretion system and motility genes that may mediate interactions with their host. Specifically, adjacent BurBor-like type 3 and T6SS-5-like type 6 secretion system operons shared among all three D. discoideum-symbiont genomes may be important for host interaction. Horizontal transfer of these secretion system operons within the amoeba host environment may have contributed to the unique ability of these symbionts to establish and maintain a symbiotic relationship with D. discoideum. IMPORTANCE Protists are a diverse group of typically single cell eukaryotes. Bacteria and archaea that form long-term symbiotic relationships with protists may evolve in additional ways than those in relationships with multicellular eukaryotes such as plants, animals, or fungi. Social amoebas are a predatory soil protist sometimes found with symbiotic bacteria living inside their cells. They present a unique opportunity to explore a naturally occurring symbiosis in a protist frequently used for studying host-pathogen interactions. We show that one amoeba-symbiont species is similar to other related bacteria in genome size and content, while the two reduced-genome-symbiont species show characteristics of genome streamlining rather than deterioration during adaptation to their host. We also identify sets of genes present in all three amoeba-symbiont genomes that are potentially used for host-symbiont interactions. Because the amoeba symbionts are distantly related, the amoeba host environment may be where these genes were shared among symbionts.

Proteomic traits vary across taxa in a coastal Antarctic phytoplankton bloom

Production and use of proteins is under strong selection in microbes, but it is unclear how proteome-level traits relate to ecological strategies. We identified and quantified proteomic traits of eukaryotic microbes and bacteria through an Antarctic phytoplankton bloom using in situ metaproteomics. Different taxa, rather than different environmental conditions, formed distinct clusters based on their ribosomal and photosynthetic proteomic proportions, and we propose that these characteristics relate to ecological differences. We defined and used a proteomic proxy for regulatory cost, which showed that SAR11 had the lowest regulatory cost of any taxa we observed at our summertime Southern Ocean study site. Haptophytes had lower regulatory cost than diatoms, which may underpin haptophyte-to-diatom bloom progression in the Ross Sea. We were able to make these proteomic trait inferences by assessing various sources of bias in metaproteomics, providing practical recommendations for researchers in the field. We have quantified several proteomic traits (ribosomal and photosynthetic proteomic proportions, regulatory cost) in eukaryotic and bacterial taxa, which can then be incorporated into trait-based models of microbial communities that reflect resource allocation strategies.

Evidence for the Rapid and Divergent Evolution of Mycoplasmas: Structural and Phylogenetic Analysis of Enolases

Mycoplasmas are a group of prokaryotes without cell walls that have evolved through several rounds of degenerative evolution. With a low cell DNA G + C content and definitively long genetic lineages, mycoplasmas are thought to be in a state of rapid evolution. However, little associated evidence has been provided. Enolase is a key enzyme in glycolysis that is widely found in all species from the three domains, and it is evolutionarily conserved. In our previous studies, enolase acted as a virulence factor and participated in cell-surface adhesion in Mycoplasma hyopneumoniae. Furthermore, unique loop regions were first found in the crystal structure of Mhp Eno. Here, enolase structures from Mycoplasma pneumoniae and Mycoplasma bovis were determined. An extra helix 7 is specific and conservatively found in almost all mycoplasma enolases, as confirmed by crystal structures and sequence alignment. Particular motifs for helix 7, which is composed of F-K/G-K-L/F-K-X-A-I, have been proposed and could be regarded as molecular markers. To our surprise, the genetic distances between any two mycoplasma enolases were obviously longer than those between the two corresponding species themselves, indicating divergent evolution of mycoplasma enolases, whereas no horizontal gene transfer was detected in mycoplasma enolase genens. Furthermore, different evolutionary patterns were adopted by different loop regions of mycoplasma enolase. Enolases from different Mycoplasma species also showed different affinities for PLG and fibronectin. Our results indicate the rapid and divergent evolution of mycoplasma enolase and mycoplasmas. This study will also aid understanding the independent evolution of Mycoplasma species after separation from their common ancestor.

Diversity, multifaceted evolution, and facultative saprotrophism in the European Batrachochytrium salamandrivorans epidemic

While emerging fungi threaten global biodiversity, the paucity of fungal genome assemblies impedes thoroughly characterizing epidemics and developing effective mitigation strategies. Here, we generate de novo genomic assemblies for six outbreaks of the emerging pathogen Batrachochytrium salamandrivorans (Bsal). We reveal the European epidemic currently damaging amphibian populations to comprise multiple, highly divergent lineages demonstrating isolate-specific adaptations and metabolic capacities. In particular, we show extensive gene family expansions and acquisitions, through a variety of evolutionary mechanisms, and an isolate-specific saprotrophic lifecycle. This finding both explains the chytrid's ability to divorce transmission from host density, producing Bsal's enigmatic host population declines, and is a key consideration in developing successful mitigation measures.

The transposable element-rich genome of the cereal pest Sitophilus oryzae

BACKGROUND

The rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

RESULTS

We sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino acids required for insect development and cuticle biosynthesis.

CONCLUSIONS

Here we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.

A pan-genome method to determine core regions of the Bacillus subtilis and Escherichia coli genomes

Background: Synthetic engineering of bacteria to produce industrial products is a burgeoning field of research and application. In order to optimize genome design, designers need to understand which genes are essential, which are optimal for growth, and locations in the genome that will be tolerated by the organism when inserting engineered cassettes.

Methods: We present a pan-genome based method for the identification of core regions in a genome that are strongly conserved at the species level.

Results: We show that the core regions determined by our method contain all or almost all essential genes. This demonstrates the accuracy of our method as essential genes should be core genes. We show that we outperform previous methods by this measure. We also explain why there are exceptions to this rule for our method.

Conclusions: We assert that synthetic engineers should avoid deleting or inserting into these core regions unless they understand and are manipulating the function of the genes in that region. Similarly, if the designer wishes to streamline the genome, non-core regions and in particular low penetrance genes would be good targets for deletion. Care should be taken to remove entire cassettes with similar penetrance of the genes within cassettes as they may harbor toxin/antitoxin genes which need to be removed in tandem. The bioinformatic approach introduced here saves considerable time and effort relative to knockout studies on single isolates of a given species and captures a broad understanding of the conservation of genes that are core to a species.

Bacterial O-GlcNAcase genes abundance decreases in ulcerative colitis patients and its administration ameliorates colitis in mice

OBJECTIVE

O-linked N-acetylglucosaminylation (O-GlcNAcylation), controlled by O-GlcNAcase (OGA) and O-GlcNAc transferase (OGT), is an important post-translational modification of eukaryotic proteins and plays an essential role in regulating gut inflammation. Gut microbiota encode various enzymes involved in O-GlcNAcylation. However, the characteristics, abundance and function of these enzymes are unknown.

DESIGN

We first investigated the structure and taxonomic distribution of bacterial OGAs and OGTs. Then, we performed metagenomic analysis to explore the OGA genes abundance in health samples and different diseases. Finally, we employed in vitro and in vivo experiments to determine the effects and mechanisms of bacterial OGAs to hydrolyse O-GlcNAcylated proteins in host cells and suppress inflammatory response in the gut.

RESULTS

We found OGAs, instead of OGTs, are enriched in Bacteroidetes and Firmicutes, the major bacterial divisions in the human gut. Most bacterial OGAs are secreted enzymes with the same conserved catalytic domain as human OGAs. A pooled analysis on 1999 metagenomic samples encompassed six diseases revealed that bacterial OGA genes were conserved in healthy human gut with high abundance, and reduced exclusively in ulcerative colitis. In vitro studies showed that bacterial OGAs could hydrolyse O-GlcNAcylated proteins in host cells, including O-GlcNAcylated NF-κB-p65 subunit, which is important for activating NF-κB signalling. In vivo studies demonstrated that gut bacteria-derived OGAs could protect mice from chemically induced colonic inflammation through hydrolysing O-GlcNAcylated proteins.

CONCLUSION

Our results reveal a previously unrecognised enzymatic activity by which gut microbiota influence intestinal physiology and highlight bacterial OGAs as a promising therapeutic strategy in colonic inflammation.

A sugar utilization phenotype contributes to the formation of genetic exchange communities in lactic acid bacteria

In prokaryotes, a major contributor to genomic evolution is the exchange of genes via horizontal gene transfer (HGT). Areas with a high density of HGT networks are defined as genetic exchange communities (GECs). Although some phenotypes associated with specific ecological niches are linked to GECs, little is known about the phenotypic influences on HGT in bacterial groups within a taxonomic family. Thanks to the published genome sequences and phenotype data of lactic acid bacteria (LAB), it is now possible to obtain more detailed information about the phenotypes that affect GECs. Here, we have investigated the relationship between HGT and internal and external environmental factors for 178 strains from 24 genera in the Lactobacillaceae family. We found a significant correlation between strains with high utilization of sugars and HGT bias. The result suggests that the phenotype of the utilization of a variety of sugars is key to the construction of GECs in this family. This feature is consistent with the fact that the Lactobacillaceae family contributes to the production of a wide variety of fermented foods by sharing niches such as those in vegetables, dairy products and brewing-related environments. This result provides the first evidence that phenotypes associated with ecological niches contribute to form GECs in the LAB family.

Diploid genomic architecture of Nitzschia inconspicua, an elite biomass production diatom

A near-complete diploid nuclear genome and accompanying circular mitochondrial and chloroplast genomes have been assembled from the elite commercial diatom species Nitzschia inconspicua. The 50 Mbp haploid size of the nuclear genome is nearly double that of model diatom Phaeodactylum tricornutum, but 30% smaller than closer relative Fragilariopsis cylindrus. Diploid assembly, which was facilitated by low levels of allelic heterozygosity (2.7%), included 14 candidate chromosome pairs composed of long, syntenic contigs, covering 93% of the total assembly. Telomeric ends were capped with an unusual 12-mer, G-rich, degenerate repeat sequence. Predicted proteins were highly enriched in strain-specific marker domains associated with cell-surface adhesion, biofilm formation, and raphe system gliding motility. Expanded species-specific families of carbonic anhydrases suggest potential enhancement of carbon concentration efficiency, and duplicated glycolysis and fatty acid synthesis pathways across cytosolic and organellar compartments may enhance peak metabolic output, contributing to competitive success over other organisms in mixed cultures. The N. inconspicua genome delivers a robust new reference for future functional and transcriptomic studies to illuminate the physiology of benthic pennate diatoms and harness their unique adaptations to support commercial algae biomass and bioproduct production.

Molecular underpinnings and biogeochemical consequences of enhanced diatom growth in a warming Southern Ocean

Phytoplankton contribute to the Southern Ocean’s (SO) ability to absorb atmospheric CO₂ and shape the stoichiometry of northward macronutrient delivery. Climate change is altering the SO environment, yet we know little about how resident phytoplankton will react to these changes. Here, we studied a natural SO community and compared responses of two prevalent, bloom-forming diatom groups to changes in temperature and iron that are projected to occur by 2100 to 2300. We found that one group, Pseudo-nitzschia, grows better under warmer low-iron conditions by managing cellular iron demand and efficiently increasing photosynthetic capacity. This ability to grow and draw down nutrients in the face of warming, regardless of iron availability, has major implications for ocean ecosystems and global nutrient cycles.

The Southern Ocean (SO) harbors some of the most intense phytoplankton blooms on Earth. Changes in temperature and iron availability are expected to alter the intensity of SO phytoplankton blooms, but little is known about how these changes will influence community composition and downstream biogeochemical processes. We performed light-saturated experimental manipulations on surface ocean microbial communities from McMurdo Sound in the Ross Sea to examine the effects of increased iron availability (+2 nM) and warming (+3 and +6 °C) on nutrient uptake, as well as the growth and transcriptional responses of two dominant diatoms, Fragilariopsis and Pseudo-nitzschia. We found that community nutrient uptake and primary productivity were elevated under both warming conditions without iron addition (relative to ambient −0.5 °C). This effect was greater than additive under concurrent iron addition and warming. Pseudo-nitzschia became more abundant under warming without added iron (especially at 6 °C), while Fragilariopsis only became more abundant under warming in the iron-added treatments. We attribute the apparent advantage Pseudo-nitzschia shows under warming to up-regulation of iron-conserving photosynthetic processes, utilization of iron-economic nitrogen assimilation mechanisms, and increased iron uptake and storage. These data identify important molecular and physiological differences between dominant diatom groups and add to the growing body of evidence for Pseudo-nitzschia’s increasingly important role in warming SO ecosystems. This study also suggests that temperature-driven shifts in SO phytoplankton assemblages may increase utilization of the vast pool of excess nutrients in iron-limited SO surface waters and thereby influence global nutrient distribution and carbon cycling.

Genome of the fatal tapeworm Sparganum proliferum uncovers mechanisms for cryptic life cycle and aberrant larval proliferation

The cryptic parasite Sparganum proliferum proliferates in humans and invades tissues and organs. Only scattered cases have been reported, but S. proliferum infection is always fatal. However, S. proliferum’s phylogeny and life cycle remain enigmatic. To investigate the phylogenetic relationships between S. proliferum and other cestode species, and to examine the mechanisms underlying pathogenicity, we sequenced the entire genomes of S. proliferum and a closely related non–life-threatening tapeworm Spirometra erinaceieuropaei. Additionally, we performed larvae transcriptome analyses of S. proliferum plerocercoid to identify genes involved in asexual reproduction in the host. The genome sequences confirmed that the S. proliferum has experienced a clearly distinct evolutionary history from S. erinaceieuropaei. Moreover, we found that nonordinal extracellular matrix coordination allows asexual reproduction in the host, and loss of sexual maturity in S. proliferum are responsible for its fatal pathogenicity to humans. Our high-quality reference genome sequences should be valuable for future studies of pseudophyllidean tapeworm biology and parasitism.

Kikuchi et al. sequence the genome of the fatal tapeworm Sparganum proliferum and a closely related non–life-threatening tapeworm Spirometra erinaceieuropaei, and describe its genomic features suggesting the natural history and molecular mechanisms underlying pathogenicity. Their findings indicate that nonordinal extracellular matrix coordination is important for its asexual reproduction, and suggest that loss of sexual maturity contributes to the fatal pathogenicity of S. proliferum to humans.

A pan-genome method to determine core regions of the Bacillus subtilis and Escherichia coli genomes

Background: Synthetic engineering of bacteria to produce industrial products is a burgeoning field of research and application. In order to optimize genome design, designers need to understand which genes are essential, which are optimal for growth, and locations in the genome that will be tolerated by the organism when inserting engineered cassettes. Methods: We present a pan-genome based method for the identification of core regions in a genome that are strongly conserved at the species level. Results: We show that the core regions determined by our method contain all or almost all essential genes. This demonstrates the accuracy of our method as essential genes should be core genes. We show that we outperform previous methods by this measure. We also explain why there are exceptions to this rule for our method. Conclusions: We assert that synthetic engineers should avoid deleting or inserting into these core regions unless they understand and are manipulating the function of the genes in that region. Similarly, if the designer wishes to streamline the genome, non-core regions and in particular low penetrance genes would be good targets for deletion. Care should be taken to remove entire cassettes with similar penetrance of the genes within cassettes as they may harbor toxin/antitoxin genes which need to be removed in tandem. The bioinformatic approach introduced here saves considerable time and effort relative to knockout studies on single isolates of a given species and captures a broad understanding of the conservation of genes that are core to a species.

A Multi-Omics Characterization of the Natural Product Potential of Tropical Filamentous Marine Cyanobacteria

Microbial natural products are important for the understanding of microbial interactions, chemical defense and communication, and have also served as an inspirational source for numerous pharmaceutical drugs. Tropical marine cyanobacteria have been highlighted as a great source of new natural products, however, few reports have appeared wherein a multi-omics approach has been used to study their natural products potential (i.e., reports are often focused on an individual natural product and its biosynthesis). This study focuses on describing the natural product genetic potential as well as the expressed natural product molecules in benthic tropical cyanobacteria. We collected from several sites around the world and sequenced the genomes of 24 tropical filamentous marine cyanobacteria. The informatics program antiSMASH was used to annotate the major classes of gene clusters. BiG-SCAPE phylum-wide analysis revealed the most promising strains for natural product discovery among these cyanobacteria. LCMS/MS-based metabolomics highlighted the most abundant molecules and molecular classes among 10 of these marine cyanobacterial samples. We observed that despite many genes encoding for peptidic natural products, peptides were not as abundant as lipids and lipopeptides in the chemical extracts. Our results highlight a number of highly interesting biosynthetic gene clusters for genome mining among these cyanobacterial samples.

Antimicrobial resistance in enteric bacteria: current state and next-generation solutions

Antimicrobial resistance is one of the largest threats to global health and imposes substantial burdens in terms of morbidity, mortality, and economic costs. The gut is a key conduit for the genesis and spread of antimicrobial resistance in enteric bacterial pathogens. Distinct bacterial species that cause enteric disease can exist as invasive enteropathogens that immediately evoke gastrointestinal distress, or pathobionts that can arise from established bacterial commensals to inflict dysbiosis and disease. Furthermore, various environmental reservoirs and stressors facilitate the evolution and transmission of resistance. In this review, we present a comprehensive discussion on circulating resistance profiles and gene mobilization strategies of the most problematic species of enteric bacterial pathogens. Importantly, we present emerging approaches toward surveillance of pathogens and their resistance elements as well as promising treatment strategies that can circumvent common resistance mechanisms.

Genomic Signatures of Honey Bee Association in an Acetic Acid Symbiont

Recent declines in the health of the honey bee have startled researchers and lay people alike as honey bees are agriculture's most important pollinator. Honey bees are important pollinators of many major crops and add billions of dollars annually to the US economy through their services. One factor that may influence colony health is the microbial community. Indeed, the honey bee worker digestive tract harbors a characteristic community of bee-specific microbes, and the composition of this community is known to impact honey bee health. However, the honey bee is a superorganism, a colony of eusocial insects with overlapping generations where nestmates cooperate, building a hive, gathering and storing food, and raising brood. In contrast to what is known regarding the honey bee worker gut microbiome, less is known of the microbes associated with developing brood, with food stores, and with the rest of the built hive environment. More recently, the microbe Bombella apis was identified as associated with nectar, with developing larvae, and with honey bee queens. This bacterium is related to flower-associated microbes such as Saccharibacter floricola and other species in the genus Saccharibacter, and initial phylogenetic analyses placed it as sister to these environmental bacteria. Here, we used comparative genomics of multiple honey bee-associated strains and the nectar-associated Saccharibacter to identify genomic changes that may be associated with the ecological transition to honey bee association. We identified several genomic differences in the honey bee-associated strains, including a complete CRISPR/Cas system. Many of the changes we note here are predicted to confer upon Bombella the ability to survive in royal jelly and defend themselves against mobile elements, including phages. Our results are a first step toward identifying potential function of this microbe in the honey bee superorganism.

Structural basis for carotenoid cleavage by an archaeal carotenoid dioxygenase

Apocarotenoids are important signaling molecules generated from carotenoids through the action of carotenoid cleavage dioxygenases (CCDs). These enzymes have a remarkable ability to cleave carotenoids at specific alkene bonds while leaving chemically similar sites within the polyene intact. Although several bacterial and eukaryotic CCDs have been characterized, the long-standing goal of experimentally visualizing a CCD-carotenoid complex at high resolution to explain this exquisite regioselectivity remains unfulfilled. CCD genes are also present in some archaeal genomes, but the encoded enzymes remain uninvestigated. Here, we address this knowledge gap through analysis of a metazoan-like archaeal CCD from Candidatus Nitrosotalea devanaterra (NdCCD). NdCCD was active toward β-apocarotenoids but did not cleave bicyclic carotenoids. It exhibited an unusual regiospecificity, cleaving apocarotenoids solely at the C14'-C13' alkene bond to produce β-apo-14'-carotenals. The structure of NdCCD revealed a tapered active site cavity markedly different from the broad active site observed for the retinal-forming Synechocystis apocarotenoid oxygenase (SynACO) but similar to the vertebrate retinoid isomerase RPE65. The structure of NdCCD in complex with its apocarotenoid product demonstrated that the site of cleavage is defined by interactions along the substrate binding cleft as well as selective stabilization of reaction intermediates at the scissile alkene. These data on the molecular basis of CCD catalysis shed light on the origins of the varied catalytic activities found in metazoan CCDs, opening the possibility of modifying their activity through rational chemical or genetic approaches.

Microbial genomic island discovery, visualization and analysis

Horizontal gene transfer (also called lateral gene transfer) is a major mechanism for microbial genome evolution, enabling rapid adaptation and survival in specific niches. Genomic islands (GIs), commonly defined as clusters of bacterial or archaeal genes of probable horizontal origin, are of particular medical, environmental and/or industrial interest, as they disproportionately encode virulence factors and some antimicrobial resistance genes and may harbor entire metabolic pathways that confer a specific adaptation (solvent resistance, symbiosis properties, etc). As large-scale analyses of microbial genomes increases, such as for genomic epidemiology investigations of infectious disease outbreaks in public health, there is increased appreciation of the need to accurately predict and track GIs. Over the past decade, numerous computational tools have been developed to tackle the challenges inherent in accurate GI prediction. We review here the main types of GI prediction methods and discuss their advantages and limitations for a routine analysis of microbial genomes in this era of rapid whole-genome sequencing. An assessment is provided of 20 GI prediction software methods that use sequence-composition bias to identify the GIs, using a reference GI data set from 104 genomes obtained using an independent comparative genomics approach. Finally, we present guidelines to assist researchers in effectively identifying these key genomic regions.

Understanding Horizontal Gene Transfer network in human gut microbiota

Background

Horizontal Gene Transfer (HGT) is the process of transferring genetic materials between species. Through sharing genetic materials, microorganisms in the human microbiota form a network. The network can provide insights into understanding the microbiota. Here, we constructed the HGT networks from the gut microbiota sequencing data and performed network analysis to characterize the HGT networks of gut microbiota.

Results

We constructed the HGT network and perform the network analysis to two typical gut microbiota datasets, a 283-sample dataset of Mother-to-Child and a 148-sample dataset of longitudinal inflammatory bowel disease (IBD) metagenome. The results indicated that (1) the HGT networks are scale-free. (2) The networks expand their complexities, sizes, and edge numbers, accompanying the early stage of lives; and microbiota established in children shared high similarity as their mother (p-value = 0.0138), supporting the transmission of microbiota from mother to child. (3) Groups harbor group-specific network edges, and network communities, which can potentially serve as biomarkers. For instances, IBD patient group harbors highly abundant communities of Proteobacteria (p-value = 0.0194) and Actinobacteria (p-value = 0.0316); children host highly abundant communities of Proteobacteria (p-value = 2.8785\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^{-5}$$\end{document}e-5) and Actinobacteria (p-value = 0.0015), and the mothers host highly abundant communities of Firmicutes (p-value = 8.0091\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^{-7}$$\end{document}e-7). IBD patient networks contain more HGT edges in pathogenic genus, including Mycobacterium, Sutterella, and Pseudomonas. Children’s networks contain more edges from Bifidobacterium and Escherichia.

Conclusion

Hence, we proposed the HGT network constructions from the gut microbiota sequencing data. The HGT networks capture the host state and the response of microbiota to the environmental and host changes, and they are essential to understand the human microbiota.

ShadowCaster: Compositional Methods under the Shadow of Phylogenetic Models to Detect Horizontal Gene Transfers in Prokaryotes

Horizontal gene transfer (HGT) plays an important role for evolutionary innovations within prokaryotic communities and is a crucial event for their survival. Several computational approaches have arisen to identify HGT events in recipient genomes. However, this has been proven to be a complex task due to the generation of a great number of false positives and the prediction disagreement among the existing methods. Phylogenetic reconstruction methods turned out to be the most reliable ones, but they are not extensible to all genes/species and are computationally demanding when dealing with large datasets. In contrast, the so-called surrogate methods that use heuristic solutions either based on nucleotide composition patterns or phyletic distribution of BLAST hits can be applied easily to the genomic scale, but they fail in identifying common HGT events. Here, we present ShadowCaster, a hybrid approach that sequentially combines nucleotide composition-based predictions by support vector machines (SVMs) under the shadow of phylogenetic models independent of tree reconstruction, to improve the detection of HGT events in prokaryotes. ShadowCaster successfully predicted close and distant HGT events in both artificial and bacterial genomes. ShadowCaster detected HGT related to heavy metal resistance in the genome of Rhodanobacter denitrificans with higher accuracy than the most popular state-of-the-art computational approaches, encompassing most of the predicted cases made by other methods. ShadowCaster is released at the GitHub platform as an open-source software under the GPLv3 license.

A genomic view of trophic and metabolic diversity in clade-specific Lamellodysidea sponge microbiomes

BACKGROUND

Marine sponges and their microbiomes contribute significantly to carbon and nutrient cycling in global reefs, processing and remineralizing dissolved and particulate organic matter. Lamellodysidea herbacea sponges obtain additional energy from abundant photosynthetic Hormoscilla cyanobacterial symbionts, which also produce polybrominated diphenyl ethers (PBDEs) chemically similar to anthropogenic pollutants of environmental concern. Potential contributions of non-Hormoscilla bacteria to Lamellodysidea microbiome metabolism and the synthesis and degradation of additional secondary metabolites are currently unknown.

RESULTS

This study has determined relative abundance, taxonomic novelty, metabolic capacities, and secondary metabolite potential in 21 previously uncharacterized, uncultured Lamellodysidea-associated microbial populations by reconstructing near-complete metagenome-assembled genomes (MAGs) to complement 16S rRNA gene amplicon studies. Microbial community compositions aligned with sponge host subgroup phylogeny in 16 samples from four host clades collected from multiple sites in Guam over a 3-year period, including representatives of Alphaproteobacteria, Gammaproteobacteria, Oligoflexia, and Bacteroidetes as well as Cyanobacteria (Hormoscilla). Unexpectedly, microbiomes from one host clade also included Cyanobacteria from the prolific secondary metabolite-producer genus Prochloron, a common tunicate symbiont. Two novel Alphaproteobacteria MAGs encoded pathways diagnostic for methylotrophic metabolism as well as type III secretion systems, and have been provisionally assigned to a new order, designated Candidatus Methylospongiales. MAGs from other taxonomic groups encoded light-driven energy production pathways using not only chlorophyll, but also bacteriochlorophyll and proteorhodopsin. Diverse heterotrophic capabilities favoring aerobic versus anaerobic conditions included pathways for degrading chitin, eukaryotic extracellular matrix polymers, phosphonates, dimethylsulfoniopropionate, trimethylamine, and benzoate. Genetic evidence identified an aerobic catabolic pathway for halogenated aromatics that may enable endogenous PBDEs to be used as a carbon and energy source.

CONCLUSIONS

The reconstruction of high-quality MAGs from all microbial taxa comprising greater than 0.1% of the sponge microbiome enabled species-specific assignment of unique metabolic features that could not have been predicted from taxonomic data alone. This information will promote more representative models of marine invertebrate microbiome contributions to host bioenergetics, the identification of potential new sponge parasites and pathogens based on conserved metabolic and physiological markers, and a better understanding of biosynthetic and degradative pathways for secondary metabolites and halogenated compounds in sponge-associated microbiota. Video Abstract.

Genome sequencing of Rigidoporus microporus provides insights on genes important for wood decay, latex tolerance and interspecific fungal interactions

Fungal plant pathogens remain a serious threat to the sustainable agriculture and forestry, despite the extensive efforts undertaken to control their spread. White root rot disease is threatening rubber tree (Hevea brasiliensis) plantations throughout South and Southeast Asia and Western Africa, causing tree mortality and severe yield losses. Here, we report the complete genome sequence of the basidiomycete fungus Rigidoporus microporus, a causative agent of the disease. Our phylogenetic analysis confirmed the position of R. microporus among the members of Hymenochaetales, an understudied group of basidiomycetes. Our analysis further identified pathogen's genes with a predicted role in the decay of plant cell wall polymers, in the utilization of latex components and in interspecific interactions between the pathogen and other fungi. We also detected putative horizontal gene transfer events in the genome of R. microporus. The reported first genome sequence of a tropical rubber tree pathogen R. microporus should contribute to the better understanding of how the fungus is able to facilitate wood decay and nutrient cycling as well as tolerate latex and utilize resinous extractives.

Multi-Omic Profiling of Melophlus Sponges Reveals Diverse Metabolomic and Microbiome Architectures that Are Non-overlapping with Ecological Neighbors

Marine sponge holobionts, defined as filter-feeding sponge hosts together with their associated microbiomes, are prolific sources of natural products. The inventory of natural products that have been isolated from marine sponges is extensive. Here, using untargeted mass spectrometry, we demonstrate that sponges harbor a far greater diversity of low-abundance natural products that have evaded discovery. While these low-abundance natural products may not be feasible to isolate, insights into their chemical structures can be gleaned by careful curation of mass fragmentation spectra. Sponges are also some of the most complex, multi-organismal holobiont communities in the oceans. We overlay sponge metabolomes with their microbiome structures and detailed metagenomic characterization to discover candidate gene clusters that encode production of sponge-derived natural products. The multi-omic profiling strategy for sponges that we describe here enables quantitative comparison of sponge metabolomes and microbiomes to address, among other questions, the ecological relevance of sponge natural products and for the phylochemical assignment of previously undescribed sponge identities.

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

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

Inferring Horizontal Gene Transfer with DarkHorse, Phylomizer, and ETE Toolkits

In this chapter, we describe how to use DarkHorse2.0 to search for xenologs in the genome of the cyanobacterium Synechococcus elongatus PCC 7942. DarkHorse is an implicit phylogenetic method that uses BLAST searches to identify proteins having close homologs of unexpected taxonomic affiliation. Once a set of putative xenologs are identified, Phylomizer is used to reconstruct phylogenetic trees. Phylomizer reproduces all the necessary steps to perform a basic phylogenetic analysis. The combined use of DarkHorse and Phylomizer allows the identification of genes incorporated into a given genome by HGT.

Current and Promising Approaches to Identify Horizontal Gene Transfer Events in Metagenomes

High-throughput shotgun metagenomics sequencing has enabled the profiling of myriad natural communities. These data are commonly used to identify gene families and pathways that were potentially gained or lost in an environment and which may be involved in microbial adaptation. Despite the widespread interest in these events, there are no established best practices for identifying gene gain and loss in metagenomics data. Horizontal gene transfer (HGT) represents several mechanisms of gene gain that are especially of interest in clinical microbiology due to the rapid spread of antibiotic resistance genes in natural communities. Several additional mechanisms of gene gain and loss, including gene duplication, gene loss-of-function events, and de novo gene birth are also important to consider in the context of metagenomes but have been less studied. This review is largely focused on detecting HGT in prokaryotic metagenomes, but methods for detecting these other mechanisms are first discussed. For this article to be self-contained, we provide a general background on HGT and the different possible signatures of this process. Lastly, we discuss how improved assembly of genomes from metagenomes would be the most straight-forward approach for improving the inference of gene gain and loss events. Several recent technological advances could help improve metagenome assemblies: long-read sequencing, determining the physical proximity of contigs, optical mapping of short sequences along chromosomes, and single-cell metagenomics. The benefits and limitations of these advances are discussed and open questions in this area are highlighted.

Diel transcriptional response of a California Current plankton microbiome to light, low iron, and enduring viral infection

Phytoplankton and associated microbial communities provide organic carbon to oceanic food webs and drive ecosystem dynamics. However, capturing those dynamics is challenging. Here, an in situ, semi-Lagrangian, robotic sampler profiled pelagic microbes at 4 h intervals over ~2.6 days in North Pacific high-nutrient, low-chlorophyll waters. We report on the community structure and transcriptional dynamics of microbes in an operationally large size class (>5 μm) predominantly populated by dinoflagellates, ciliates, haptophytes, pelagophytes, diatoms, cyanobacteria (chiefly Synechococcus), prasinophytes (chiefly Ostreococcus), fungi, archaea, and proteobacteria. Apart from fungi and archaea, all groups exhibited 24-h periodicity in some transcripts, but larger portions of the transcriptome oscillated in phototrophs. Periodic photosynthesis-related transcripts exhibited a temporal cascade across the morning hours, conserved across diverse phototrophic lineages. Pronounced silica:nitrate drawdown, a high flavodoxin to ferredoxin transcript ratio, and elevated expression of other Fe-stress markers indicated Fe-limitation. Fe-stress markers peaked during a photoperiodically adaptive time window that could modulate phytoplankton response to seasonal Fe-limitation. Remarkably, we observed viruses that infect the majority of abundant taxa, often with total transcriptional activity synchronized with putative hosts. Taken together, these data reveal a microbial plankton community that is shaped by recycled production and tightly controlled by Fe-limitation and viral activity.

Understanding the Genetic Code

The universal triple-nucleotide genetic code is often viewed as a given, randomly selected through evolution. However, as summarized in this article, many observations and deductions within structural and thermodynamic frameworks help to explain the forces that must have shaped the code during the early evolution of life on Earth.

Comparative Genomics of Cyanobacterial Symbionts Reveals Distinct, Specialized Metabolism in Tropical Dysideidae Sponges

Marine sponges are recognized as valuable sources of bioactive metabolites and renowned as petri dishes of the sea, providing specialized niches for many symbiotic microorganisms. Sponges of the family Dysideidae are well documented to be chemically talented, often containing high levels of polyhalogenated compounds, terpenoids, peptides, and other classes of bioactive small molecules. This group of tropical sponges hosts a high abundance of an uncultured filamentous cyanobacterium, Hormoscilla spongeliae Here, we report the comparative genomic analyses of two phylogenetically distinct Hormoscilla populations, which reveal shared deficiencies in essential pathways, hinting at possible reasons for their uncultivable status, as well as differing biosynthetic machinery for the production of specialized metabolites. One symbiont population contains clustered genes for expanded polybrominated diphenylether (PBDE) biosynthesis, while the other instead harbors a unique gene cluster for the biosynthesis of the dysinosin nonribosomal peptides. The hybrid sequencing and assembly approach utilized here allows, for the first time, a comprehensive look into the genomes of these elusive sponge symbionts.IMPORTANCE Natural products provide the inspiration for most clinical drugs. With the rise in antibiotic resistance, it is imperative to discover new sources of chemical diversity. Bacteria living in symbiosis with marine invertebrates have emerged as an untapped source of natural chemistry. While symbiotic bacteria are often recalcitrant to growth in the lab, advances in metagenomic sequencing and assembly now make it possible to access their genetic blueprint. A cell enrichment procedure, combined with a hybrid sequencing and assembly approach, enabled detailed genomic analysis of uncultivated cyanobacterial symbiont populations in two chemically rich tropical marine sponges. These population genomes reveal a wealth of secondary metabolism potential as well as possible reasons for historical difficulties in their cultivation.

Horizontal gene transfer allowed the emergence of broad host range entomopathogens

The emergence of new pathogenic fungi has profoundly impacted global biota, but the underlying mechanisms behind host shifts remain largely unknown. The endophytic insect pathogen Metarhizium robertsii evolved from fungi that were plant associates, and entomopathogenicity is a more recently acquired adaptation. Here we report that the broad host-range entomopathogen M. robertsii has 18 genes that are derived via horizontal gene transfer (HGT). The necessity of degrading insect cuticle served as a major selective pressure to retain these genes, as 12 are up-regulated during penetration; 6 were confirmed to have a role in penetration, and their collective actions are indispensable for infection. Two lipid-carrier genes are involved in utilizing epicuticular lipids, and a third (MrNPC2a) facilitates hemocoel colonization. Three proteases degraded the procuticular protein matrix, which facilitated up-regulation of other cuticle-degrading enzymes. The three lipid carriers and one of the proteases are present in all analyzed Metarhizium species and are essential for entomopathogenicity. Acquisition of another protease (MAA_01413) in an ancestor of broad host-range lineages contributed to their host-range expansion, as heterologous expression in the locust specialist Metarhizium acridum enabled it to kill caterpillars. Our work reveals that HGT was a key mechanism in the emergence of entomopathogenicity in Metarhizium from a plant-associated ancestor and in subsequent host-range expansion by some Metarhizium lineages.

MetaCHIP: community-level horizontal gene transfer identification through the combination of best-match and phylogenetic approaches

BACKGROUND

Metagenomic datasets provide an opportunity to study horizontal gene transfer (HGT) on the level of a microbial community. However, current HGT detection methods cannot be applied to community-level datasets or require reference genomes. Here, we present MetaCHIP, a pipeline for reference-independent HGT identification at the community level.

RESULTS

Assessment of MetaCHIP's performance on simulated datasets revealed that it can predict HGTs with various degrees of genetic divergence from metagenomic datasets. The results also indicated that the detection of very recent gene transfers (i.e. those with low levels of genetic divergence) from metagenomics datasets is largely affected by the read assembly step. Comparison of MetaCHIP with a previous analysis on soil bacteria showed a high level of consistency for the prediction of recent HGTs and revealed a large number of additional non-recent gene transfers, which can provide new biological and ecological insight. Assessment of MetaCHIP's performance on real metagenomic datasets confirmed the role of HGT in the spread of genes related to antibiotic resistance in the human gut microbiome. Further testing also showed that functions related to energy production and conversion as well as carbohydrate transport and metabolism are frequently transferred among free-living microorganisms.

CONCLUSION

MetaCHIP provides an opportunity to study HGTs among members of a microbial community and therefore has several applications in the field of microbial ecology and evolution. MetaCHIP is implemented in Python and freely available at https://github.com/songweizhi/MetaCHIP .

Pangenomic comparison of globally distributed Poribacteria associated with sponge hosts and marine particles

Candidatus Poribacteria is a little-known bacterial phylum, previously characterized by partial genomes from a single sponge host, but never isolated in culture. We have reconstructed multiple genome sequences from four different sponge genera and compared them to recently reported, uncharacterized Poribacteria genomes from the open ocean, discovering shared and unique functional characteristics. Two distinct, habitat-linked taxonomic lineages were identified, designated Entoporibacteria (sponge-associated) and Pelagiporibacteria (free-living). These lineages differed in flagellar motility and chemotaxis genes unique to Pelagiporibacteria, and highly expanded families of restriction endonucleases, DNA methylases, transposases, CRISPR repeats, and toxin-antitoxin gene pairs in Entoporibacteria. Both lineages shared pathways for facultative anaerobic metabolism, denitrification, fermentation, organosulfur compound utilization, type IV pili, cellulosomes, and bacterial proteosomes. Unexpectedly, many features characteristic of eukaryotic host association were also shared, including genes encoding the synthesis of eukaryotic-like cell adhesion molecules, extracellular matrix digestive enzymes, phosphoinositol-linked membrane glycolipids, and exopolysaccharide capsules. Complete Poribacteria 16S rRNA gene sequences were found to contain multiple mismatches to "universal" 16S rRNA gene primer sets, substantiating concerns about potential amplification failures in previous studies. A newly designed primer set corrects these mismatches, enabling more accurate assessment of Poribacteria abundance in diverse marine habitats where it may have previously been overlooked.

Collection, Culturing, and Genome Analyses of Tropical Marine Filamentous Benthic Cyanobacteria

Decreasing sequencing costs has sparked widespread investigation of the use of microbial genomics to accelerate the discovery and development of natural products for therapeutic uses. Tropical marine filamentous cyanobacteria have historically produced many structurally novel natural products, and therefore present an excellent opportunity for the systematic discovery of new metabolites via the information derived from genomics and molecular genetics. Adequate knowledge transfer and institutional know-how are important to maintain the capability for studying filamentous cyanobacteria due to their unusual microbial morphology and characteristics. Here, we describe workflows, procedures, and commentary on sample collection, cultivation, genomic DNA generation, bioinformatics tools, and biosynthetic pathway analysis concerning filamentous cyanobacteria.

HgtSIM: a simulator for horizontal gene transfer (HGT) in microbial communities

The development and application of metagenomic approaches have provided an opportunity to study and define horizontal gene transfer (HGT) on the level of microbial communities. However, no current metagenomic data simulation tools offers the option to introduce defined HGT within a microbial community. Here, we present HgtSIM, a pipeline to simulate HGT event among microbial community members with user-defined mutation levels. It was developed for testing and benchmarking pipelines for recovering HGTs from complex microbial datasets. HgtSIM is implemented in Python3 and is freely available at: https://github.com/songweizhi/HgtSIM.

Alienness: Rapid Detection of Candidate Horizontal Gene Transfers across the Tree of Life

Horizontal gene transfer (HGT) is the transmission of genes between organisms by other means than parental to offspring inheritance. While it is prevalent in prokaryotes, HGT is less frequent in eukaryotes and particularly in Metazoa. Here, we propose Alienness, a taxonomy-aware web application available at http://alienness.sophia.inra.fr. Alienness parses BLAST results against public libraries to rapidly identify candidate HGT in any genome of interest. Alienness takes as input the result of a BLAST of a whole proteome of interest against any National Center for Biotechnology Information (NCBI) protein library. The user defines recipient (e.g., Metazoa) and donor (e.g., bacteria, fungi) branches of interest in the NCBI taxonomy. Based on the best BLAST E-values of candidate donor and recipient taxa, Alienness calculates an Alien Index (AI) for each query protein. An AI > 0 indicates a better hit to candidate donor than recipient taxa and a possible HGT. Higher AI represent higher gap of E-values between candidate donor and recipient and a more likely HGT. We confirmed the accuracy of Alienness on phylogenetically confirmed HGT of non-metazoan origin in plant-parasitic nematodes. Alienness scans whole proteomes to rapidly identify possible HGT in any species of interest and thus fosters exploration of HGT more easily and largely across the tree of life.

Effect of the environment on horizontal gene transfer between bacteria and archaea

Background

Horizontal gene transfer, the transfer and incorporation of genetic material between different species of organisms, has an important but poorly quantified role in the adaptation of microbes to their environment. Previous work has shown that genome size and the number of horizontally transferred genes are strongly correlated. Here we consider how genome size confuses the quantification of horizontal gene transfer because the number of genes an organism accumulates over time depends on its evolutionary history and ecological context (e.g., the nutrient regime for which it is adapted).

Results

We investigated horizontal gene transfer between archaea and bacteria by first counting reciprocal BLAST hits among 448 bacterial and 57 archaeal genomes to find shared genes. Then we used the DarkHorse algorithm, a probability-based, lineage-weighted method (Podell & Gaasterland, 2007), to identify potential horizontally transferred genes among these shared genes. By removing the effect of genome size in the bacteria, we have identified bacteria with unusually large numbers of shared genes with archaea for their genome size. Interestingly, archaea and bacteria that live in anaerobic and/or high temperature conditions are more likely to share unusually large numbers of genes. However, high salt was not found to significantly affect the numbers of shared genes. Numbers of shared (genome size-corrected, reciprocal BLAST hits) and transferred genes (identified by DarkHorse) were strongly correlated. Thus archaea and bacteria that live in anaerobic and/or high temperature conditions are more likely to share horizontally transferred genes. These horizontally transferred genes are over-represented by genes involved in energy conversion as well as the transport and metabolism of inorganic ions and amino acids.

Conclusions

Anaerobic and thermophilic bacteria share unusually large numbers of genes with archaea. This is mainly due to horizontal gene transfer of genes from the archaea to the bacteria. In general, these transfers are from archaea that live in similar oxygen and temperature conditions as the bacteria that receive the genes. Potential hotspots of horizontal gene transfer between archaea and bacteria include hot springs, marine sediments, and oil wells. Cold spots for horizontal transfer included dilute, aerobic, mesophilic environments such as marine and freshwater surface waters.

The metabolic potential of the single cell genomes obtained from the Challenger Deep, Mariana Trench within the candidate superphylum Parcubacteria (OD1)

Candidate phyla (CP) are broad phylogenetic clusters of organisms that lack cultured representatives. Included in this fraction is the candidate Parcubacteria superphylum. Specific characteristics that have been ascribed to the Parcubacteria include reduced genome size, limited metabolic potential and exclusive reliance on fermentation for energy acquisition. The study of new environmental niches, such as the marine versus terrestrial subsurface, often expands the understanding of the genetic potential of taxonomic groups. For this reason, we analyzed 12 Parcubacteria single amplified genomes (SAGs) from sediment samples collected within the Challenger Deep of the Mariana Trench, obtained during the Deepsea Challenge (DSC) Expedition. Many of these SAGs are closely related to environmental sequences obtained from deep-sea environments based on 16S rRNA gene similarity and BLAST matches to predicted proteins. DSC SAGs encode features not previously identified in Parcubacteria obtained from other habitats. These include adaptation to oxidative stress, polysaccharide modification and genes associated with respiratory nitrate reduction. The DSC SAGs are also distinguished by relative greater abundance of genes for nucleotide and amino acid biosynthesis, repair of alkylated DNA and the synthesis of mechanosensitive ion channels. These results present an expanded view of the Parcubacteria, among members residing in an ultra-deep hadal environment.

The diversity of CO2-concentrating mechanisms in marine diatoms as inferred from their genetic content

Marine diatoms are one of the most ecologically significant primary producers in the ocean. Most diatoms use a CO2-concentrating mechanism (CCM) to overcome the scarcity of CO2 in the ocean and limitations of the carbon-fixing enzyme Rubisco. However, the CCMs in model diatoms differ substantially in their genetic make-up and structural organization. To assess the extent of CCM diversity in marine diatoms more generally, we analyzed genome and transcriptome data from 31 diatom strains to identify putative CCM genes, examine the overall CCM architecture, and study CCM development in the context of the evolutionary history of these diatoms. Key CCM genes [carbonic anhydrases (CAs) and solute carrier 4 (SLC4) bicarbonate transporters] identified in the diatoms were placed into groups of likely orthologs by sequence similarity (OrthoMCL) and phylogenetic methods. These analyses indicated that diatoms seem to share similar HCO3- transporters, but possess a variety of CAs that have either undergone extensive diversification within the diatom lineage or have been acquired through horizontal gene transfer. Hierarchical clustering of the diatom species based on their CCM gene content suggests that CCM development is largely congruent with evolution of diatom species, despite some notable differences in CCM genes even among closely related species.

Metagenomic discovery of polybrominated diphenyl ether biosynthesis by marine sponges

Naturally produced polybrominated diphenyl ethers (PBDEs) pervade the marine environment and structurally resemble toxic man-made brominated flame retardants. PBDEs bioaccumulate in marine animals and are likely transferred to the human food chain. However, the biogenic basis for PBDE production in one of their most prolific sources, marine sponges of the order Dysideidae, remains unidentified. Here, we report the discovery of PBDE biosynthetic gene clusters within sponge-microbiome-associated cyanobacterial endosymbionts through the use of an unbiased metagenome-mining approach. Using expression of PBDE biosynthetic genes in heterologous cyanobacterial hosts, we correlate the structural diversity of naturally produced PBDEs to modifications within PBDE biosynthetic gene clusters in multiple sponge holobionts. Our results establish the genetic and molecular foundation for the production of PBDEs in one of the most abundant natural sources of these molecules, further setting the stage for a metagenomic-based inventory of other PBDE sources in the marine environment.