Inferring horizontal gene transfer

Characteristics and phylogenetic distribution of megaplasmids and prediction of a putative chromid in Pseudomonas aeruginosa

Research on megaplasmids that contribute to the spread of antimicrobial resistance (AMR) in Pseudomonas aeruginosa strains has grown in recent years due to the now widely used technologies allowing long-read sequencing. Here, we systematically analyzed distinct and consistent genetic characteristics of megaplasmids found in P. aeruginosa. Our data provide information on their phylogenetic distribution and hypotheses tracing the potential evolutionary paths of megaplasmids. Most of the megaplasmids we found belong to the IncP-2-type, with conserved and syntenic genetic backbones carrying modules of genes associated with chemotaxis apparatus, tellurite resistance and plasmid replication, segregation, and transmission. Extensively variable regions harbor abundant AMR genes, especially those encoding β-lactamases such as VIM-2, IMP-45, and KPC variants, which are high-risk elements in nosocomial infection. IncP-2 megaplasmids act as effective vehicles transmitting AMR genes to diverse regions. One evolutionary model of the origin of megaplasmids claims that chromids can develop from megaplasmids. These chromids have been characterized as an intermediate between a megaplasmid and a chromosome, also containing core genes that can be found on the chromosome but not on the megaplasmid. Using in silico prediction, we identified the "PABCH45 unnamed replicon" as a putative chromid in P. aeruginosa, which shows a much higher similarity and closer phylogenetic relationship to chromosomes than to megaplasmids while also encoding plasmid-like partition genes. We propose that such a chromid could facilitate genome expansion, allowing for more rapid adaptations to novel ecological niches or selective conditions, in comparison to megaplasmids.

Novel misos shape distinct microbial ecologies: opportunities for flavourful sustainable food innovation

Fermentation is resurgent around the world as people seek healthier, more sustainable, and tasty food options. This study explores the microbial ecology of miso, a traditional Japanese fermented paste, made with novel regional substrates to develop new plant-based foods. Eight novel miso varieties were developed using different protein-rich substrates: yellow peas, Gotland lentils, and fava beans (each with two treatments: standard and nixtamalisation), as well as rye bread and soybeans. The misos were produced at Noma, a restaurant in Copenhagen, Denmark. Samples were analysed with biological and technical triplicates at the beginning and end of fermentation. We also incorporated in this study six samples of novel misos produced following the same recipe at Inua, a former affiliate restaurant of Noma in Tokyo, Japan. To analyse microbial community structure and diversity, metabarcoding (16S and ITS) and shotgun metagenomic analyses were performed. The misos contain a greater range of microbes than is currently described for miso in the literature. The composition of the novel yellow pea misos was notably similar to the traditional soybean ones, suggesting they are a good alternative, which supports our culinary collaborators' sensory conclusions. For bacteria, we found that overall substrate had the strongest effect, followed by time, treatment (nixtamalisation), and geography. For fungi, there was a slightly stronger effect of geography and a mild effect of substrate, and no significant effects for treatment or time. Based on an analysis of metagenome-assembled genomes (MAGs), strains of Staphylococccus epidermidis differentiated according to substrate. Carotenoid biosynthesis genes in these MAGs appeared in strains from Japan but not from Denmark, suggesting a possible gene-level geographical effect. The benign and possibly functional presence of S. epidermidis in these misos, a species typically associated with the human skin microbiome, suggests possible adaptation to the miso niche, and the flow of microbes between bodies and foods in certain fermentation as more common than is currently recognised. This study improves our understanding of miso ecology, highlights the potential for developing novel misos using diverse local ingredients, and suggests how fermentation innovation can contribute to studies of microbial ecology and evolution.

In Silico Analysis of the Phylogenetic and Physiological Characteristics of Sphingobium indicum B90A: A Hexachlorocyclohexane-Degrading Bacterium

The study focuses on the in silico genomic characterization of Sphingobium indicum B90A, revealing a wealth of genes involved in stress response, carbon monoxide oxidation, β-carotene biosynthesis, heavy metal resistance, and aromatic compound degradation, suggesting its potential as a bioremediation agent. Furthermore, genomic adaptations among nine Sphingomonad strains were explored, highlighting shared core genes via pangenome analysis, including those related to the shikimate pathway and heavy metal resistance. The majority of genes associated with aromatic compound degradation, heavy metal resistance, and stress response were found within genomic islands across all strains. Sphingobium indicum UT26S exhibited the highest number of genomic islands, while Sphingopyxis alaskensis RB2256 had the maximum fraction of its genome covered by genomic islands. The distribution of lin genes varied among the strains, indicating diverse genetic responses to environmental pressures. Additionally, in silico evidence of horizontal gene transfer (HGT) between plasmids pSRL3 and pISP3 of the Sphingobium and Sphingomonas genera, respectively, has been provided. The manuscript offers novel insights into strain B90A, highlighting its role in horizontal gene transfer and refining evolutionary relationships among Sphingomonad strains. The discovery of stress response genes and the czcABCD operon emphasizes the potential of Sphingomonads in consortia development, supported by genomic island analysis.

Phylogenetic analyses reveal insights into interdomain horizontal gene transfer of microbial lipases

Microbial lipases play a pivotal role in a wide range of biotechnological processes and in the human skin microbiome. However, their evolution remains poorly understood. Accessing the evolutionary process of lipases could contribute to future applications in health and biotechnology. We investigated genetic events associated with the evolutionary trajectory of the microbial family LIP lipases. Using phylogenetic analysis, we identified two distinct horizontal gene transfer (HGT) events from Bacteria to Fungi. Further analysis of human cutaneous mycobiome members such as the lipophilic Malassezia yeasts and CUG-Ser-1 clade (including Candida sp. and other microorganisms associated with cutaneous mycobiota) revealed recent evolutionary processes, with multiple gene duplication events. The Lid region of fungal lipases, crucial for substrate interaction, exhibits varying degrees of conservation among different groups. Our findings suggest the adaptability of the fungal LIP family in various genetic and metabolic contexts and its potential role in niche exploration.

A global survey of prokaryotic genomes reveals the eco-evolutionary pressures driving horizontal gene transfer

Horizontal gene transfer, the exchange of genetic material through means other than reproduction, is a fundamental force in prokaryotic genome evolution. Genomic persistence of horizontally transferred genes has been shown to be influenced by both ecological and evolutionary factors. However, there is limited availability of ecological information about species other than the habitats from which they were isolated, which has prevented a deeper exploration of ecological contributions to horizontal gene transfer. Here we focus on transfers detected through comparison of individual gene trees to the species tree, assessing the distribution of gene-exchanging prokaryotes across over a million environmental sequencing samples. By analysing detected horizontal gene transfer events, we show distinct functional profiles for recent versus old events. Although most genes transferred are part of the accessory genome, genes transferred earlier in evolution tend to be more ubiquitous within present-day species. We find that co-occurring, interacting and high-abundance species tend to exchange more genes. Finally, we show that host-associated specialist species are most likely to exchange genes with other host-associated specialist species, whereas species found across different habitats have similar gene exchange rates irrespective of their preferred habitat. Our study covers an unprecedented scale of integrated horizontal gene transfer and environmental information, highlighting broad eco-evolutionary trends.

Modeling the mosaic structure of bacterial genomes to infer their evolutionary history

The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.

Predicting horizontal gene transfers with perfect transfer networks

BACKGROUND

Horizontal gene transfer inference approaches are usually based on gene sequences: parametric methods search for patterns that deviate from a particular genomic signature, while phylogenetic methods use sequences to reconstruct the gene and species trees. However, it is well-known that sequences have difficulty identifying ancient transfers since mutations have enough time to erase all evidence of such events. In this work, we ask whether character-based methods can predict gene transfers. Their advantage over sequences is that homologous genes can have low DNA similarity, but still have retained enough important common motifs that allow them to have common character traits, for instance the same functional or expression profile. A phylogeny that has two separate clades that acquired the same character independently might indicate the presence of a transfer even in the absence of sequence similarity.

OUR CONTRIBUTIONS

We introduce perfect transfer networks, which are phylogenetic networks that can explain the character diversity of a set of taxa under the assumption that characters have unique births, and that once a character is gained it is rarely lost. Examples of such traits include transposable elements, biochemical markers and emergence of organelles, just to name a few. We study the differences between our model and two similar models: perfect phylogenetic networks and ancestral recombination networks. Our goals are to initiate a study on the structural and algorithmic properties of perfect transfer networks. We then show that in polynomial time, one can decide whether a given network is a valid explanation for a set of taxa, and show how, for a given tree, one can add transfer edges to it so that it explains a set of taxa. We finally provide lower and upper bounds on the number of transfers required to explain a set of taxa, in the worst case.

The Theory of Gene Family Histories

Most genes are part of larger families of evolutionary-related genes. The history of gene families typically involves duplications and losses of genes as well as horizontal transfers into other organisms. The reconstruction of detailed gene family histories, i.e., the precise dating of evolutionary events relative to phylogenetic tree of the underlying species has remained a challenging topic despite their importance as a basis for detailed investigations into adaptation and functional evolution of individual members of the gene family. The identification of orthologs, moreover, is a particularly important subproblem of the more general setting considered here. In the last few years, an extensive body of mathematical results has appeared that tightly links orthology, a formal notion of best matches among genes, and horizontal gene transfer. The purpose of this chapter is to broadly outline some of the key mathematical insights and to discuss their implication for practical applications. In particular, we focus on tree-free methods, i.e., methods to infer orthology or horizontal gene transfer as well as gene trees, species trees, and reconciliations between them without using a priori knowledge of the underlying trees or statistical models for the inference of phylogenetic trees. Instead, the initial step aims to extract binary relations among genes.

Dancing the Nanopore limbo - Nanopore metagenomics from small DNA quantities for bacterial genome reconstruction

BACKGROUND

While genome-resolved metagenomics has revolutionized our understanding of microbial and genetic diversity in environmental samples, assemblies of short-reads often result in incomplete and/or highly fragmented metagenome-assembled genomes (MAGs), hampering in-depth genomics. Although Nanopore sequencing has increasingly been used in microbial metagenomics as long reads greatly improve the assembly quality of MAGs, the recommended DNA quantity usually exceeds the recoverable amount of DNA of environmental samples. Here, we evaluated lower-than-recommended DNA quantities for Nanopore library preparation by determining sequencing quality, community composition, assembly quality and recovery of MAGs.

RESULTS

We generated 27 Nanopore metagenomes using the commercially available ZYMO mock community and varied the amount of input DNA from 1000 ng (the recommended minimum) down to 1 ng in eight steps. The quality of the generated reads remained stable across all input levels. The read mapping accuracy, which reflects how well the reads match a known reference genome, was consistently high across all libraries. The relative abundance of the species in the metagenomes was stable down to input levels of 50 ng. High-quality MAGs (> 95% completeness, ≤ 5% contamination) could be recovered from metagenomes down to 35 ng of input material. When combined with publicly available Illumina reads for the mock community, Nanopore reads from input quantities as low as 1 ng improved the quality of hybrid assemblies.

CONCLUSION

Our results show that the recommended DNA amount for Nanopore library preparation can be substantially reduced without any adverse effects to genome recovery and still bolster hybrid assemblies when combined with short-read data. We posit that the results presented herein will enable studies to improve genome recovery from low-biomass environments, enhancing microbiome understanding.

Cyanobacteria as a critical reservoir of the environmental antimicrobial resistome

Antimicrobial resistance (AMR) is predicted to cause a worldwide annual toll of 10 million deaths by 2050. This looming public health threat has been linked to antibiotic overuse and pollution, which places selective pressures on AMR maintenance and transfer in and between microbial populations. We examined the distribution, diversity and potential mobility of AMR genes in cyanobacteria. While cyanobacteria are not pathogenic, we hypothesised that they could be a major environmental reservoir for AMR genes. Genes encoding AMR to seven antimicrobial drug classes were found in 10% of cyanobacterial genomes. AMR genes were found in 13% of freshwater, 19% of terrestrial, 34% of symbiotic, 2% of thermal spring, and 3% of marine genomes. AMR genes were found in five cyanobacterial orders with 23% of Nostocales and 8% of Oscillatoriales strains containing AMR genes. The most frequently observed alleles were ansamycin resistance genes, which were present in 7% of strains. AMR genes responsible for resistance to broad-spectrum β-lactams, chloramphenicols, tetracyclines, macrolides, and aminoglycosides were associated with mobile genetic elements or plasmid replicons or both. These results suggest that cyanobacteria are an extensive reservoir, and potential vector, for AMR genes in diverse terrestrial and aquatic habitats.

A mesophilic phosphotriesterase-like lactonase shows high stability and proficiency as quorum quenching enzyme

The problem of biofilm formation is a serious concern under various pathological conditions such as extensive burns, wounds in diabetic patients, bedsores, cystic fibrosis, nosocomial infections from implantable medical devices such as catheters, valves, etc. Environmental diffusion of biofilm (in pools, wet floors, industrial food plants) that could represent a reservoir of antibiotic resistant bacteria constitues an additional issue. In this work is described a lactonase from Rhodococcus erythropolis, a phosphotriesterase-like lactonase (PLL) enzyme, which has already been studied in the past and can be used for containment of biofilm formation. The protein is 28% and 40% identical with respect to the Pseudomonas diminuta PTE and the thermostable Saccharolobus solfataricus SsoPox respectively. The protein was obtained starting from a synthetic His-tagged gene, expressed in E. coli, purified and further characterized. New properties, not previously known or deducible from its sequence, have been highlighted. These properties are: the enzyme is thermophilic and thermostable even though it originates from a mesophilic bacterium; the enzyme has a long (months) shelf life at 4 °C; the enzyme is not only stable to low concentrations of the oxidant H2O2 but even activated by it at high concentrations; the enzyme proved to be a proficient quorum quenching enzyme, able to hydrolase acyl-homoserine lactones 3oxoC12-HSL and C4-HSL, and can inhibit up to 60% the formation of Pseudomonas aeruginosa (PAO1) biofilm. These different properties make the lactonase useful to fight resistant bacteria that induce inflammatory and infectious processes mediated by the quorum sensing mechanism.

Insights into early evolutionary adaptations of the Akkermansia genus to the vertebrate gut

Akkermansia, a relevant mucin degrader from the vertebrate gut microbiota, is a member of the deeply branched Verrucomicrobiota, as well as the only known member of this phylum to be described as inhabitants of the gut. Only a few Akkermansia species have been officially described so far, although there is genomic evidence addressing the existence of more species-level variants for this genus. This niche specialization makes Akkermansia an interesting model for studying the evolution of microorganisms to their adaptation to the gastrointestinal tract environment, including which kind of functions were gained when the Akkermansia genus originated or how the evolutionary pressure functions over those genes. In order to gain more insight into Akkermansia adaptations to the gastrointestinal tract niche, we performed a phylogenomic analysis of 367 high-quality Akkermansia isolates and metagenome-assembled genomes, in addition to other members of Verrucomicrobiota. This work was focused on three aspects: the definition of Akkermansia genomic species clusters and the calculation and functional characterization of the pangenome for the most represented species; the evolutionary relationship between Akkermansia and their closest relatives from Verrucomicrobiota, defining the gene families which were gained or lost during the emergence of the last Akkermansia common ancestor (LAkkCA) and; the evaluation of the evolutionary pressure metrics for each relevant gene family of main Akkermansia species. This analysis found 25 Akkermansia genomic species clusters distributed in two main clades, divergent from their non-Akkermansia relatives. Pangenome analyses suggest that Akkermansia species have open pangenomes, and the gene gain/loss model indicates that genes associated with mucin degradation (both glycoside hydrolases and peptidases), (micro)aerobic metabolism, surface interaction, and adhesion were part of LAkkCA. Specifically, mucin degradation is a very ancestral innovation involved in the origin of Akkermansia. Horizontal gene transfer detection suggests that Akkermansia could receive genes mostly from unknown sources or from other Gram-negative gut bacteria. Evolutionary metrics suggest that Akkemansia species evolved differently, and even some conserved genes suffered different evolutionary pressures among clades. These results suggest a complex evolutionary landscape of the genus and indicate that mucin degradation could be an essential feature in Akkermansia evolution as a symbiotic species.

Genomic characterization of the antiviral arsenal of Actinobacteria

Phages are ubiquitous in nature, and bacteria with very different genomics, metabolisms, and lifestyles are subjected to their predation. Yet, the defence systems that allow bacteria to resist their phages have rarely been explored experimentally outside a very limited number of model organisms. Actinobacteria (Actinomycetota) are a phylum of GC-rich Gram-positive bacteria, which often produce an important diversity of secondary metabolites. Despite being ubiquitous in a wide range of environments, from soil to fresh and sea water but also the gut microbiome, relatively little is known about the anti-phage arsenal of Actinobacteria. In this work, we used DefenseFinder to systematically detect 131 anti-phage defence systems in 22803 fully sequenced prokaryotic genomes, among which are 2253 Actinobacteria of more than 700 species. We show that, like other bacteria, Actinobacteria encode many diverse anti-phage systems that are often encoded on mobile genetic elements. We further demonstrate that most detected defence systems are absent or rarer in Actinobacteria than in other bacteria, while a few rare systems are enriched (notably gp29-gp30 and Wadjet). We characterize the spatial distribution of anti-phage systems on Streptomyces chromosomes and show that some defence systems (e.g. RM systems) tend to be encoded in the core region, while others (e.g. Lamassu and Wadjet) are enriched towards the extremities. Overall, our results suggest that Actinobacteria might be a source of novel anti-phage systems and provide clues to characterize mechanistic aspects of known anti-phage systems.

Evolutionary barriers to horizontal gene transfer in macrophage-associated Salmonella

Horizontal gene transfer (HGT) is a powerful evolutionary force facilitating bacterial adaptation and emergence of novel phenotypes. Several factors, including environmental ones, are predicted to restrict HGT, but we lack systematic and experimental data supporting these predictions. Here, we address this gap by measuring the relative fitness of 44 genes horizontally transferred from Escherichia coli to Salmonella enterica in infection-relevant environments. We estimated the distribution of fitness effects in each environment and identified that dosage-dependent effects across different environments are a significant barrier to HGT. The majority of genes were found to be deleterious. We also found longer genes had stronger negative fitness consequences than shorter ones, showing that gene length was negatively associated with HGT. Furthermore, fitness effects of transferred genes were found to be environmentally dependent. In summary, a substantial fraction of transferred genes had a significant fitness cost on the recipient, with both gene characteristics and the environment acting as evolutionary barriers to HGT.

Virus-like transposons cross the species barrier and drive the evolution of genetic incompatibilities

Horizontal gene transfer, the movement of genetic material between species, has been reported across all major eukaryotic lineages. However, the underlying mechanisms of transfer and their impact on genome evolution are still poorly understood. While studying the evolutionary origin of a selfish element in the nematode Caenorhabditis briggsae, we discovered that Mavericks, ancient virus-like transposons related to giant viruses and virophages, are one of the long-sought vectors of horizontal gene transfer. We found that Mavericks gained a novel herpesvirus-like fusogen in nematodes, leading to the widespread exchange of cargo genes between extremely divergent species, bypassing sexual and genetic barriers spanning hundreds of millions of years. Our results show how the union between viruses and transposons causes horizontal gene transfer and ultimately genetic incompatibilities in natural populations.

Potential mobile units drive the horizontal transfer of phytoplasma effector phyllogen genes

Phytoplasmas are obligate intracellular plant pathogenic bacteria that can induce phyllody, which is a type of abnormal floral organ development. Phytoplasmas possess phyllogens, which are effector proteins that cause phyllody in plants. Phylogenetic comparisons of phyllogen and 16S rRNA genes have suggested that phyllogen genes undergo horizontal transfer between phytoplasma species and strains. However, the mechanisms and evolutionary implications of this horizontal gene transfer are unclear. Here, we analyzed synteny in phyllogen flanking genomic regions from 17 phytoplasma strains that were related to six 'Candidatus' species, including three strains newly sequenced in this study. Many of the phyllogens were flanked by multicopy genes within potential mobile units (PMUs), which are putative transposable elements found in phytoplasmas. The multicopy genes exhibited two distinct patterns of synteny that correlated with the linked phyllogens. The low level of sequence identities and partial truncations found among these phyllogen flanking genes indicate that the PMU sequences are deteriorating, whereas the highly conserved sequences and functions (e.g., inducing phyllody) of the phyllogens suggest that the latter are important for phytoplasma fitness. Furthermore, although their phyllogens were similar, PMUs in strains related to 'Ca. P. asteris' were often located in different regions of the genome. These findings strongly indicate that PMUs drive the horizontal transfer of phyllogens among phytoplasma species and strains. These insights improve our understanding of how symptom-determinant genes have been shared among phytoplasmas.

Bacterial origin of a key innovation in the evolution of the vertebrate eye

The vertebrate eye was described by Charles Darwin as one of the greatest potential challenges to a theory of natural selection by stepwise evolutionary processes. While numerous evolutionary transitions that led to the vertebrate eye have been explained, some aspects appear to be vertebrate specific with no obvious metazoan precursor. One critical difference between vertebrate and invertebrate vision hinges on interphotoreceptor retinoid-binding protein (IRBP, also known as retinol-binding protein, RBP3), which enables the physical separation and specialization of cells in the vertebrate visual cycle by promoting retinoid shuttling between cell types. While IRBP has been functionally described, its evolutionary origin has remained elusive. Here, we show that IRBP arose via acquisition of novel genetic material from bacteria by interdomain horizontal gene transfer (iHGT). We demonstrate that a gene encoding a bacterial peptidase was acquired prior to the radiation of extant vertebrates >500 Mya and underwent subsequent domain duplication and neofunctionalization to give rise to vertebrate IRBP. Our phylogenomic analyses on >900 high-quality genomes across the tree of life provided the resolution to distinguish contamination in genome assemblies from true instances of horizontal acquisition of IRBP and led us to discover additional independent transfers of the same bacterial peptidase gene family into distinct eukaryotic lineages. Importantly, this work illustrates the evolutionary basis of a key transition that led to the vertebrate visual cycle and highlights the striking impact that acquisition of bacterial genes has had on vertebrate evolution.

The structure and diversity of microbial communities in Paederus fuscipes (Coleoptera: Staphylinidae): from ecological paradigm to pathobiome

BACKGROUND

Paederus fuscipes is medically the most famous rove beetle, which causes dermatitis or conjunctivitis in humans, as well as gastrointestinal toxicosis in livestock, via releasing toxic hemolymph containing pederin. Pedrin biosynthesis genes have been identified in uncultured Pseudomonas-like endosymbionts that are speculated to be acquired through a horizontal transfer. However, the composition of the P. fuscipes microbial community, especially of the gut and genital microbiome, remains unclear. This study was aimed to characterize the structure and diversity of P. fuscipes-associated bacterial communities in terms of gender, organ, and location using the Illumina HiSeq platform in the southern littorals of Caspian Sea.

RESULTS

The OTUs identified from P. fuscipes specimens were collapsed into 40 phyla, 112 classes, 249 orders, 365 families, 576 genera, and 106 species. The most abundant families were Pseudomonadaceae, Spiroplasmataceae, Weeksellaceae, Enterococcaceae, and Rhizobiaceae, respectively. Thirty top genera made up > 94% of the P. fuscipes microbiome, with predominating Pseudomonas, followed by the Spiroplasma, Apibacter, Enterococcus, Dysgonomonas, Sebaldella, Ruminococcus, and Wolbachia. Interesting dissimilarities were also discovered within and between the beetle microbiomes in terms of genders and organs. Analyses showed that Spiroplasma / Apibacter as well as Pseudomonas / Pseudomonas were the most abundant in the genitals / intestines of male and female beetles, respectively. Bacterial richness did not display any significant difference in the three provinces but was higher in male beetles than in females and more in the genitals than intestines.

CONCLUSIONS

The present study identified Pseudomonas-like endobacterium as a common symbiont of P. fuscipes beetles; this bacterium begins its journey from gut and genitalia of females to reach the male rove beetles. Additionally, male and female rove beetles were characterized by distinctive microbiota in different organs, likely reflecting different functions and/or adaptation processes. Evidence of the extension of P. fuscipes microbiome from the environmental paradigm to the pathobiome was also presented herein. A comprehensive survey of P. fuscipes microbiome components may eventually lead to ecological insights into the production and utilization of defensive compound of pederin and also the management of linear dermatitis with the use of available antibiotics against bacterial pathogens released by the beetles. Video Abstract.

Two novel genes identified by large-scale transcriptomic analysis are essential for biofilm and rugose colony development of Vibrio vulnificus

Many pathogenic bacteria form biofilms to survive under environmental stresses and host immune defenses. Differential expression (DE) analysis of the genes in biofilm and planktonic cells under a single condition, however, has limitations to identify the genes essential for biofilm formation. Independent component analysis (ICA), a machine learning algorithm, was adopted to comprehensively identify the biofilm genes of Vibrio vulnificus, a fulminating human pathogen, in this study. ICA analyzed the large-scale transcriptome data of V. vulnificus cells under various biofilm and planktonic conditions and then identified a total of 72 sets of independently co-regulated genes, iModulons. Among the three iModulons specifically activated in biofilm cells, BrpT-iModulon mainly consisted of known genes of the regulon of BrpT, a transcriptional regulator controlling biofilm formation of V. vulnificus. Interestingly, the BrpT-iModulon additionally contained two novel genes, VV1_3061 and VV2_1694, designated as cabH and brpN, respectively. cabH and brpN were shared in other Vibrio species and not yet identified by DE analyses. Genetic and biochemical analyses revealed that cabH and brpN are directly up-regulated by BrpT. The deletion of cabH and brpN impaired the robust biofilm and rugose colony formation. CabH, structurally similar to the previously known calcium-binding matrix protein CabA, was essential for attachment to the surface. BrpN, carrying an acyltransferase-3 domain as observed in BrpL, played an important role in exopolysaccharide production. Altogether, ICA identified two novel genes, cabH and brpN, which are regulated by BrpT and essential for the development of robust biofilms and rugose colonies of V. vulnificus.

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.

Botrytis cinerea identifies host plants via the recognition of antifungal capsidiol to induce expression of a specific detoxification gene

The gray mold pathogen Botrytis cinerea has a broad host range, causing disease in >400 plant species, but it is not known how this pathogen evolved this polyxenous nature. Botrytis cinerea can metabolize a wide range of phytoalexins, including the stilbenoid resveratrol in grape, and the sesquiterpenoids capsidiol in tobacco and rishitin in potato and tomato. In this study, we analyzed the metabolism of sesquiterpenoid phytoalexins by B. cinerea. Capsidiol was dehydrogenated to capsenone, which was then further oxidized, while rishitin was directly oxidized to epoxy- or hydroxyrishitins, indicating that B. cinerea has separate mechanisms to detoxify structurally similar sesquiterpenoid phytoalexins. RNA-seq analysis revealed that a distinct set of genes were induced in B. cinerea when treated with capsidiol or rishitin, suggesting that B. cinerea can distinguish structurally similar phytoalexins to activate appropriate detoxification mechanisms. The gene most highly upregulated by capsidiol treatment encoded a dehydrogenase, designated Bccpdh. Heterologous expression of Bccpdh in a capsidiol-sensitive plant symbiotic fungus, Epichloë festucae, resulted in an acquired tolerance of capsidiol and the ability to metabolize capsidiol to capsenone, while B. cinerea Δbccpdh mutants became relatively sensitive to capsidiol. The Δbccpdh mutant showed reduced virulence on the capsidiol producing Nicotiana and Capsicum species but remained fully pathogenic on potato and tomato. Homologs of Bccpdh are found in taxonomically distant Ascomycota fungi but not in related Leotiomycetes species, suggesting that B. cinerea acquired the ancestral Bccpdh by horizontal gene transfer, thereby extending the pathogenic host range of this polyxenous pathogen to capsidiol-producing plant species.

Unexpected absence of ribosomal protein genes from metagenome-assembled genomes

Metagenome-assembled genomes (MAGs) have revealed the hidden diversity and functions of uncultivated microbes, but their reconstruction from metagenomes remains a computationally difficult task. Repetitive or exogenous sequences, such as ribosomal RNA and horizontally transferred genes, are frequently absent from MAGs because of misassembly and binning errors. Here, we report that ribosomal protein genes are also often absent from MAGs, although they are neither repetitive nor exogenous. Comprehensive analyses of more than 190,000 MAGs revealed that these genes could be missing in more than 20-40% of near-complete (i.e., with completeness of 90% or higher) MAGs. While some uncultivated environmental microbes intrinsically lack some ribosomal protein genes, we found that this unexpected absence is largely due to special evolutionary patterns of codon usage bias in ribosomal protein genes and algorithmic characteristics of metagenomic binning, which is dependent on tetranucleotide frequencies of contigs. This problem reflects the microbial life-history strategy. Fast-growing microbes tend to have this difficulty, likely because of strong evolutionary pressures on ribosomal protein genes toward the efficient assembly of ribosomes. Our observations caution those who study genomics and phylogeny of uncultivated microbes, the diversity and evolution of microbial genes in the central dogma, and bioinformatics in metagenomics.

hgtseq: A Standard Pipeline to Study Horizontal Gene Transfer

Horizontal gene transfer (HGT) is well described in prokaryotes: it plays a crucial role in evolution, and has functional consequences in insects and plants. However, less is known about HGT in humans. Studies have reported bacterial integrations in cancer patients, and microbial sequences have been detected in data from well-known human sequencing projects. Few of the existing tools for investigating HGT are highly automated. Thanks to the adoption of Nextflow for life sciences workflows, and to the standards and best practices curated by communities such as nf-core, fully automated, portable, and scalable pipelines can now be developed. Here we present nf-core/hgtseq to facilitate the analysis of HGT from sequencing data in different organisms. We showcase its performance by analysing six exome datasets from five mammals. Hgtseq can be run seamlessly in any computing environment and accepts data generated by existing exome and whole-genome sequencing projects; this will enable researchers to expand their analyses into this area. Fundamental questions are still open about the mechanisms and the extent or role of horizontal gene transfer: by releasing hgtseq we provide a standardised tool which will enable a systematic investigation of this phenomenon, thus paving the way for a better understanding of HGT.

Target-enriched long-read sequencing (TELSeq) contextualizes antimicrobial resistance genes in metagenomes

BACKGROUND

Metagenomic data can be used to profile high-importance genes within microbiomes. However, current metagenomic workflows produce data that suffer from low sensitivity and an inability to accurately reconstruct partial or full genomes, particularly those in low abundance. These limitations preclude colocalization analysis, i.e., characterizing the genomic context of genes and functions within a metagenomic sample. Genomic context is especially crucial for functions associated with horizontal gene transfer (HGT) via mobile genetic elements (MGEs), for example antimicrobial resistance (AMR). To overcome this current limitation of metagenomics, we present a method for comprehensive and accurate reconstruction of antimicrobial resistance genes (ARGs) and MGEs from metagenomic DNA, termed target-enriched long-read sequencing (TELSeq).

RESULTS

Using technical replicates of diverse sample types, we compared TELSeq performance to that of non-enriched PacBio and short-read Illumina sequencing. TELSeq achieved much higher ARG recovery (>1,000-fold) and sensitivity than the other methods across diverse metagenomes, revealing an extensive resistome profile comprising many low-abundance ARGs, including some with public health importance. Using the long reads generated by TELSeq, we identified numerous MGEs and cargo genes flanking the low-abundance ARGs, indicating that these ARGs could be transferred across bacterial taxa via HGT.

CONCLUSIONS

TELSeq can provide a nuanced view of the genomic context of microbial resistomes and thus has wide-ranging applications in public, animal, and human health, as well as environmental surveillance and monitoring of AMR. Thus, this technique represents a fundamental advancement for microbiome research and application. Video abstract.

Dissemination of metaldehyde catabolic pathways is driven by mobile genetic elements in Proteobacteria

Bioremediation of metaldehyde from drinking water using metaldehyde-degrading strains has recently emerged as a promising alternative. Whole-genome sequencing was used to obtain full genomes for metaldehyde degraders Acinetobacter calcoaceticus E1 and Sphingobium CMET-H. For the former, the genetic context of the metaldehyde-degrading genes had not been explored, while for the latter, none of the degrading genes themselves had been identified. In A. calcoaceticus E1, IS91 and IS6-family insertion sequences (ISs) were found surrounding the metaldehyde-degrading gene cluster located in plasmid pAME76. This cluster was located in closely-related plasmids and associated to identical ISs in most metaldehyde-degrading β- and γ-Proteobacteria, indicating horizontal gene transfer (HGT). For Sphingobium CMET-H, sequence analysis suggested a phytanoyl-CoA family oxygenase as a metaldehyde-degrading gene candidate due to its close homology to a previously identified metaldehyde-degrading gene known as mahX. Heterologous gene expression in Escherichia coli alongside degradation tests verified its functional significance and the degrading gene homolog was henceforth called mahS. It was found that mahS is hosted within the conjugative plasmid pSM1 and its genetic context suggested a crossover between the metaldehyde and acetoin degradation pathways. Here, specific replicons and ISs responsible for maintaining and dispersing metaldehyde-degrading genes in α, β and γ-Proteobacteria through HGT were identified and described. In addition, a homologous gene implicated in the first step of metaldehyde utilisation in an α-Proteobacteria was uncovered. Insights into specific steps of this possible degradation pathway are provided.

Non-synonymous to synonymous substitutions suggest that orthologs tend to keep their functions, while paralogs are a source of functional novelty

Orthologs separate after lineages split from each other and paralogs after gene duplications. Thus, orthologs are expected to remain more functionally coherent across lineages, while paralogs have been proposed as a source of new functions. Because protein functional divergence follows from non-synonymous substitutions, we performed an analysis based on the ratio of non-synonymous to synonymous substitutions (dN/dS), as proxy for functional divergence. We used five working definitions of orthology, including reciprocal best hits (RBH), among other definitions based on network analyses and clustering. The results showed that orthologs, by all definitions tested, had values of dN/dS noticeably lower than those of paralogs, suggesting that orthologs generally tend to be more functionally stable than paralogs. The differences in dN/dS ratios remained suggesting the functional stability of orthologs after eliminating gene comparisons with potential problems, such as genes with high codon usage biases, low coverage of either of the aligned sequences, or sequences with very high similarities. Separation by percent identity of the encoded proteins showed that the differences between the dN/dS ratios of orthologs and paralogs were more evident at high sequence identity, less so as identity dropped. The last results suggest that the differences between dN/dS ratios were partially related to differences in protein identity. However, they also suggested that paralogs undergo functional divergence relatively early after duplication. Our analyses indicate that choosing orthologs as probably functionally coherent remains the right approach in comparative genomics.

Horizontal gene transfer from genetically modified plants - Regulatory considerations

Gene technology regulators receive applications seeking permission for the environmental release of genetically modified (GM) plants, many of which possess beneficial traits such as improved production, enhanced nutrition and resistance to drought, pests and diseases. The regulators must assess the risks to human and animal health and to the environment from releasing these GM plants. One such consideration, of many, is the likelihood and potential consequence of the introduced or modified DNA being transferred to other organisms, including people. While such gene transfer is most likely to occur to sexually compatible relatives (vertical gene transfer), horizontal gene transfer (HGT), which is the acquisition of genetic material that has not been inherited from a parent, is also a possibility considered during these assessments. Advances in HGT detection, aided by next generation sequencing, have demonstrated that HGT occurrence may have been previously underestimated. In this review, we provide updated evidence on the likelihood, factors and the barriers for the introduced or modified DNA in GM plants to be horizontally transferred into a variety of recipients. We present the legislation and frameworks the Australian Gene Technology Regulator adheres to with respect to the consideration of risks posed by HGT. Such a perspective may generally be applicable to regulators in other jurisdictions as well as to commercial and research organisations who develop GM plants.

Integrated Metabolomic-Genomic Workflows Accelerate Microbial Natural Product Discovery

The pairing of analytical chemistry with genomic techniques represents a new wave in natural product chemistry. With an increase in the availability of sequencing and assembly of microbial genomes, interrogation into the biosynthetic capability of producers with valuable secondary metabolites is possible. However, without the development of robust, accessible, and medium to high throughput tools, the bottleneck in pairing metabolic potential and compound isolation will continue. Several innovative approaches have proven useful in the nascent stages of microbial genome-informed drug discovery. Here, we consider a number of these approaches which have led to prioritization of strain targets and have mitigated rediscovery rates. Likewise, we discuss integration of principles of comparative evolutionary studies and retrobiosynthetic predictions to better understand biosynthetic mechanistic details and link genome sequence to structure. Lastly, we discuss advances in engineering, chemistry, and molecular networking and other computational approaches that are accelerating progress in the field of omic-informed natural product drug discovery. Together, these strategies enhance the synergy between cutting edge omics, chemical characterization, and computational technologies that pitch the discovery of natural products with pharmaceutical and other potential applications to the crest of the wave where progress is ripe for rapid advances.

A novel Tetrahymena thermophila sterol C-22 desaturase belongs to the Fatty Acid Hydroxylase/Desaturase superfamily

Sterols in eukaryotic cells play important roles in modulating membrane fluidity and in cell signaling and trafficking. During evolution, a combination of gene losses and acquisitions gave rise to an extraordinary diversity of sterols in different organisms. The sterol C-22 desaturase identified in plants and fungi as a cytochrome P-450 monooxygenase evolved from the first eukaryotic cytochrome P450 and was lost in many lineages. Although the ciliate Tetrahymena thermophila desaturates sterols at the C-22 position, no cytochrome P-450 orthologs are present in the genome. Here, we aim to identify the genes responsible for the desaturation as well as their probable origin. We used gene knockout and yeast heterologous expression approaches to identify two putative genes, retrieved from a previous transcriptomic analysis, as sterol C-22 desaturases. Furthermore, we demonstrate using bioinformatics and evolutionary analyses that both genes encode a novel type of sterol C-22 desaturase that belongs to the large fatty acid hydroxylase/desaturase superfamily and the genes originated by genetic duplication prior to functional diversification. These results stress the widespread existence of nonhomologous isofunctional enzymes among different lineages of the tree of life as well as the suitability for the use of T. thermophila as a valuable model to investigate the evolutionary process of large enzyme families.

Mapping Strengths and Weaknesses of Different Clustering Approaches to Deciphering Bacterial Chimerism

Bacterial genomes are chimeras of DNA of different ancestries. Deconstructing chimeric genomes is central to understanding the evolutionary trajectories of their disparate components and thus the organisms as a whole in the light of their evolutionary contexts. Of specific interest is to delineate and quantify native (vertically inherited) and alien (horizontally acquired) components of bacterial genomes and also specify genomic fractions that represent different donor sources. An agglomerative clustering procedure that prioritizes grouping of proximal similar genomic segments has previously been invoked for this purpose in conjunction with a recursive segmentation procedure. Surprisingly, however, the relative strengths and weaknesses of different clustering approaches to deciphering bacterial chimerism have not yet been investigated, despite the need to robustly interpret tens of thousands of completely sequenced bacterial genomes and nearly complete genome assemblies available in the public databases. To bridge this knowledge gap and develop more robust approaches, we assessed different clustering methods, including segment order based (proximal) clustering, hierarchical clustering, affinity propagation clustering, and a novel network clustering approach on chimeric genomes modeled after bacterial genomes representing a broad spectrum of compositional complexity. Although segment order-based clustering and network clustering compared favorably with the other approaches in discriminating between native and alien DNA at genome optimized settings, network clustering did consistently better than other methods at parametric settings optimized on all test genomes together. Segment order-based clustering and hierarchical clustering outperformed other methods in alien DNA identification while preserving donor identity in the genomes. Our study highlights the strengths and weaknesses of different approaches and suggests how this can be leveraged to achieve a more robust deconstruction of bacterial chimerism.

Current Methods for Recombination Detection in Bacteria

The role of genetic exchanges, i.e., homologous recombination (HR) and horizontal gene transfer (HGT), in bacteria cannot be overestimated for it is a pivotal mechanism leading to their evolution and adaptation, thus, tracking the signs of recombination and HGT events is importance both for fundamental and applied science. To date, dozens of bioinformatics tools for revealing recombination signals are available, however, their pros and cons as well as the spectra of solvable tasks have not yet been systematically reviewed. Moreover, there are two major groups of software. One aims to infer evidence of HR, while the other only deals with horizontal gene transfer (HGT). However, despite seemingly different goals, all the methods use similar algorithmic approaches, and the processes are interconnected in terms of genomic evolution influencing each other. In this review, we propose a classification of novel instruments for both HR and HGT detection based on the genomic consequences of recombination. In this context, we summarize available methodologies paying particular attention to the type of traceable events for which a certain program has been designed.

Phosphonate production by marine microbes: Exploring new sources and potential function

Phosphonates are a class of phosphorus metabolites characterized by a highly stable C-P bond. Phosphonates accumulate to high concentrations in seawater, fuel a large fraction of marine methane production, and serve as a source of phosphorus to microbes inhabiting nutrient-limited regions of the oligotrophic ocean. Here, we show that 15% of all bacterioplankton in the surface ocean have genes phosphonate synthesis and that most belong to the abundant groups Prochlorococcus and SAR11. Genomic and chemical evidence suggests that phosphonates are incorporated into cell-surface phosphonoglycoproteins that may act to mitigate cell mortality by grazing and viral lysis. These results underscore the large global biogeochemical impact of relatively rare but highly expressed traits in numerically abundant groups of marine bacteria.

Phosphonates are organophosphorus metabolites with a characteristic C-P bond. They are ubiquitous in the marine environment, their degradation broadly supports ecosystem productivity, and they are key components of the marine phosphorus (P) cycle. However, the microbial producers that sustain the large oceanic inventory of phosphonates as well as the physiological and ecological roles of phosphonates are enigmatic. Here, we show that phosphonate synthesis genes are rare but widely distributed among diverse bacteria and archaea, including Prochlorococcus and SAR11, the two major groups of bacteria in the ocean. In addition, we show that Prochlorococcus can allocate over 40% of its total cellular P-quota toward phosphonate production. However, we find no evidence that Prochlorococcus uses phosphonates for surplus P storage, and nearly all producer genomes lack the genes necessary to degrade and assimilate phosphonates. Instead, we postulate that phosphonates are associated with cell-surface glycoproteins, suggesting that phosphonates mediate ecological interactions between the cell and its surrounding environment. Our findings indicate that the oligotrophic surface ocean phosphonate pool is sustained by a relatively small fraction of the bacterioplankton cells allocating a significant portion of their P quotas toward secondary metabolism and away from growth and reproduction.

Habitat Adaptation Drives Speciation of a Streptomyces Species with Distinct Habitats and Disparate Geographic Origins

Microbial diversification is driven by geographic and ecological factors, but how the relative importance of these factors varies among species, geographic scales, and habitats remains unclear. Streptomyces, a genus of antibiotic-producing, spore-forming, and widespread bacteria, offers a robust model for identifying the processes underlying population differentiation. We examined the population structure of 37 Streptomyces olivaceus strains isolated from various sources, showing that they diverged into two habitat-associated (free-living and insect-associated) and geographically disparate lineages. More frequent gene flow within than between the lineages confirmed genetic isolation in S. olivaceus. Geographic isolation could not explain the genetic isolation; instead, habitat type was a strong predictor of genetic distance when controlling for geographic distance. The identification of habitat-specific genetic variations, including genes involved in regulation, resource use, and secondary metabolism, suggested a significant role of habitat adaptation in the diversification process. Physiological assays revealed fitness trade-offs under different environmental conditions in the two lineages. Notably, insect-associated isolates could outcompete free-living isolates in a free-iron-deficient environment. Furthermore, substrate (e.g., sialic acid and glycogen) utilization but not thermal traits differentiated the two lineages. Overall, our results argue that adaptive processes drove ecological divergence among closely related streptomycetes, eventually leading to dispersal limitation and gene flow barriers between the lineages. S. olivaceus may best be considered a species complex consisting of two cryptic species. IMPORTANCE Both isolation by distance and isolation by environment occur in bacteria, and different diversification patterns may apply to different species. Streptomyces species, typified by producing useful natural products, are widespread in nature and possess high genetic diversity. However, the ecological processes and evolutionary mechanisms that shape their distribution are not well understood. Here, we show that the population structure of a ubiquitous Streptomyces species complex matches its habitat distribution and can be defined by gene flow discontinuities. Using comparative genomics and physiological assays, we reveal that gains and losses of specific genomic traits play a significant role in the transition between free-living and host-associated lifestyles, driving speciation of the species. These results provide new insights into the evolutionary trajectory of Streptomyces and the notion of species.

RecPD: A Recombination-aware measure of phylogenetic diversity

A critical step in studying biological features (e.g., genetic variants, gene families, metabolic capabilities, or taxa) is assessing their diversity and distribution among a sample of individuals. Accurate assessments of these patterns are essential for linking features to traits or outcomes of interest and understanding their functional impact. Consequently, it is of crucial importance that the measures employed for quantifying feature diversity can perform robustly under any evolutionary scenario. However, the standard measures used for quantifying and comparing the distribution of features, such as prevalence, phylogenetic diversity, and related approaches, either do not take into consideration evolutionary history, or assume strictly vertical patterns of inheritance. Consequently, these approaches cannot accurately assess diversity for features that have undergone recombination or horizontal transfer. To address this issue, we have devised RecPD, a novel recombination-aware phylogenetic-diversity statistic for measuring the distribution and diversity of features under all evolutionary scenarios. RecPD utilizes ancestral-state reconstruction to map the presence / absence of features onto ancestral nodes in a species tree, and then identifies potential recombination events in the evolutionary history of the feature. We also derive several related measures from RecPD that can be used to assess and quantify evolutionary dynamics and correlation of feature evolutionary histories. We used simulation studies to show that RecPD reliably reconstructs feature evolutionary histories under diverse recombination and loss scenarios. We then applied RecPD in two diverse real-world scenarios including a preliminary study type III effector protein families secreted by the plant pathogenic bacterium Pseudomonas syringae and growth phenotypes of the Pseudomonas genus and demonstrate that prevalence is an inadequate measure that obscures the potential impact of recombination. We believe RecPD will have broad utility for revealing and quantifying complex evolutionary processes for features at any biological level.

Phylogenetic diversity is an important concept utilized in evolutionary ecology which has extensive applications in population genetics to help us understand how evolutionary processes have distributed genetic variation among individuals of a species, and how this impacts phenotypic diversification over time. However, existing approaches for studying phylogenetic diversity largely assume that the genetic features follow vertical inheritance, which is frequently violated in the case of microbial genomes due to horizontal transfer. To address this shortcoming, we present RecPD, a recombination-aware phylogenetic diversity measure, which incorporates ancestral state reconstruction to quantify the phylogenetic diversity of genetic features mapped onto a species phylogeny. Through simulation experiments we show that RecPD robustly reconstructs the evolutionary histories of features evolving under various scenarios of recombination and loss. When applied to a real-world example of type III secreted effector protein families from the plant pathogenic bacterium Pseudomonas syringae, RecPD reveals that horizontal transfer has played an important role in shaping the phylogenetic distributions of a substantial proportion of families across the P. syringae species complex. Furthermore, we demonstrate that the traditional measures of feature prevalence are unsuitable as a measure for comparing feature diversity. We also provide a R package implementation of RecPD for public use: https://github.com/cedatorma/recpd.

A Large-Scale Multiple Genome Comparison of Acidophilic Archaea (pH ≤ 5.0) Extends Our Understanding of Oxidative Stress Responses in Polyextreme Environments

Acidophilic archaea thrive in anaerobic and aerobic low pH environments (pH < 5) rich in dissolved heavy metals that exacerbate stress caused by the production of reactive oxygen species (ROS) such as hydrogen peroxide (H2O2), hydroxyl radical (OH) and superoxide (O2−). ROS react with lipids, proteins and nucleic acids causing oxidative stress and damage that can lead to cell death. Herein, genes and mechanisms potentially involved in ROS mitigation are predicted in over 200 genomes of acidophilic archaea with sequenced genomes. These organisms are often be subjected to simultaneous multiple stresses such as high temperature, high salinity, low pH and high heavy metal loads. Some of the topics addressed include: (1) the phylogenomic distribution of these genes and what this can tell us about the evolution of these mechanisms in acidophilic archaea; (2) key differences in genes and mechanisms used by acidophilic versus non-acidophilic archaea and between acidophilic archaea and acidophilic bacteria and (3) how comparative genomic analysis predicts novel genes or pathways involved in oxidative stress responses in archaea and likely horizontal gene transfer (HGT) events.

Bioinformatic and Mechanistic Analysis of the Palmerolide PKS-NRPS Biosynthetic Pathway From the Microbiome of an Antarctic Ascidian

Complex interactions exist between microbiomes and their hosts. Increasingly, defensive metabolites that have been attributed to host biosynthetic capability are now being recognized as products of host-associated microbes. These unique metabolites often have bioactivity targets in human disease and can be purposed as pharmaceuticals. Polyketides are a complex family of natural products that often serve as defensive metabolites for competitive or pro-survival purposes for the producing organism, while demonstrating bioactivity in human diseases as cholesterol lowering agents, anti-infectives, and anti-tumor agents. Marine invertebrates and microbes are a rich source of polyketides. Palmerolide A, a polyketide isolated from the Antarctic ascidian Synoicum adareanum, is a vacuolar-ATPase inhibitor with potent bioactivity against melanoma cell lines. The biosynthetic gene clusters (BGCs) responsible for production of secondary metabolites are encoded in the genomes of the producers as discrete genomic elements. A candidate palmerolide BGC was identified from a S. adareanum microbiome-metagenome based on a high degree of congruence with a chemical structure-based retrobiosynthetic prediction. Protein family homology analysis, conserved domain searches, active site and motif identification were used to identify and propose the function of the ∼75 kbp trans-acyltransferase (AT) polyketide synthase-non-ribosomal synthase (PKS-NRPS) domains responsible for the stepwise synthesis of palmerolide A. Though PKS systems often act in a predictable co-linear sequence, this BGC includes multiple trans-acting enzymatic domains, a non-canonical condensation termination domain, a bacterial luciferase-like monooxygenase (LLM), and is found in multiple copies within the metagenome-assembled genome (MAG). Detailed inspection of the five highly similar pal BGC copies suggests the potential for biosynthesis of other members of the palmerolide chemical family. This is the first delineation of a biosynthetic gene cluster from an Antarctic microbial species, recently proposed as Candidatus Synoicihabitans palmerolidicus. These findings have relevance for fundamental knowledge of PKS combinatorial biosynthesis and could enhance drug development efforts of palmerolide A through heterologous gene expression.

Contribution of a ZIP-family protein to manganese uptake and infective endocarditis virulence in Streptococcus sanguinis

Streptococcus sanguinis is an important cause of infective endocarditis. In strain SK36, the ABC‐family manganese transporter, SsaACB, is essential for virulence. We have now identified a ZIP‐family protein, TmpA, as a secondary manganese transporter. A tmpA mutant had no phenotype, but a ΔssaACB ΔtmpA mutant was more attenuated for serum growth and for virulence in a rabbit model than its ΔssaACB parent. The growth of both mutants was restored by supplemental manganese, but the ΔssaACB ΔtmpA mutant required twenty‐fold more and accumulated less. Although ZIP‐family proteins are known for zinc and iron transport, TmpA‐mediated transport of either metal was minimal. While ssaACB appears ubiquitous in St. sanguinis, tmpA was present in a majority of strains and a mntH gene encoding an NRAMP‐family transporter was identified in relatively few, including VMC66. As in SK36, deletion of ssaACB greatly diminished VMC66 endocarditis virulence and serum growth, and deletion of tmpA from this mutant diminished virulence further. Virulence was not significantly altered by deletion of mntH from either VMC66 or its ΔssaACB mutant. This and the accompanying paper together suggest that SsaACB is of primary importance for endocarditis virulence while secondary transporters TmpA and MntH contribute to growth under differing conditions.

Closed Genome Sequence of Yavru, a Novel Arthrobacter globiformis Phage

We characterized the complete genome sequence of actinobacteriophage Yavru (Siphoviridae), a cluster FE bacteriophage infecting Arthrobacter globiformis NRRL B-2979; it was 89.5% identical to cluster FE phage Whytu, with a capsid width of 50 nm and a tail length of 90 nm. The genome was 15,193 bp in length, with 23 predicted protein-coding genes.

Cultivation of a vampire: 'Candidatus Absconditicoccus praedator'

Halorhodospira halophila, one of the most-xerophilic halophiles, inhabits biophysically stressful and energetically expensive, salt-saturated alkaline brines. Here, we report an additional stress factor that is biotic: a diminutive Candidate-Phyla-Radiation bacterium, that we named 'Ca. Absconditicoccus praedator' M39-6, which predates H. halophila M39-5, an obligately photosynthetic, anaerobic purple-sulfur bacterium. We cultivated this association (isolated from the hypersaline alkaline Lake Hotontyn Nur, Mongolia) and characterized their biology. 'Ca. Absconditicoccus praedator' is the first stably cultivated species from the candidate class-level lineage Gracilibacteria (order-level lineage Absconditabacterales). Its closed-and-curated genome lacks genes for the glycolytic, pentose phosphate- and Entner-Doudoroff pathways which would generate energy/reducing equivalents and produce central carbon currencies. Therefore, 'Ca. Absconditicoccus praedator' is dependent on host-derived building blocks for nucleic acid-, protein-, and peptidoglycan synthesis. It shares traits with (the uncultured) 'Ca. Vampirococcus lugosii', which is also of the Gracilibacteria lineage. These are obligate parasitic lifestyle, feeding on photosynthetic anoxygenic Gammaproteobacteria, and absorption of host cytoplasm. Commonalities in their genomic composition and structure suggest that the entire Absconditabacterales lineage consists of predatory species which act to cull the populations of their respective host bacteria. Cultivation of vampire : host associations can shed light on unresolved aspects of their metabolism and ecosystem dynamics at life-limiting extremes.

GFICLEE: ultrafast tree-based phylogenetic profile method inferring gene function at the genomic-wide level

Background

Phylogenetic profiling is widely used to predict novel members of large protein complexes and biological pathways. Although methods combined with phylogenetic trees have significantly improved prediction accuracy, computational efficiency is still an issue that limits its genome-wise application.

Results

Here we introduce a new tree-based phylogenetic profiling algorithm named GFICLEE, which infers common single and continuous loss (SCL) events in the evolutionary patterns. We validated our algorithm with human pathways from three databases and compared the computational efficiency with current tree-based with 10 different scales genome dataset. Our algorithm has a better predictive performance with high computational efficiency.

Conclusions

The GFICLEE is a new method to infers genome-wide gene function. The accuracy and computational efficiency of GFICLEE make it possible to explore gene functions at the genome-wide level on a personal computer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08070-7.

Evolutionary Divergence of the Wsp Signal Transduction Systems in Beta- and Gammaproteobacteria

Bacteria rapidly adapt to their environment by integrating external stimuli through diverse signal transduction systems. Pseudomonas aeruginosa, for example, senses surface contact through the Wsp signal transduction system to trigger the production of cyclic di-GMP. Diverse mutations in wsp genes that manifest enhanced biofilm formation are frequently reported in clinical isolates of P. aeruginosa and in biofilm studies of Pseudomonas spp. and Burkholderia cenocepacia. In contrast to the convergent phenotypes associated with comparable wsp mutations, we demonstrate that the Wsp system in B. cenocepacia does not impact intracellular cyclic di-GMP levels, unlike that in Pseudomonas spp. Our current mechanistic understanding of the Wsp system is based entirely on the study of four Pseudomonas spp., and its phylogenetic distribution remains unknown. Here, we present a broad phylogenetic analysis to show that the Wsp system originated in the betaproteobacteria and then horizontally transferred to Pseudomonas spp., the sole member of the gammaproteobacteria. Alignment of 794 independent Wsp systems with reported mutations from the literature identified key amino acid residues that fall within and outside annotated functional domains. Specific residues that are highly conserved but uniquely modified in B. cenocepacia likely define mechanistic differences among Wsp systems. We also find the greatest sequence variation in the extracellular sensory domain of WspA, indicating potential adaptations to diverse external stimuli beyond surface contact sensing. This study emphasizes the need to better understand the breadth of functional diversity of the Wsp system as a major regulator of bacterial adaptation beyond B. cenocepacia and select Pseudomonas spp. IMPORTANCE The Wsp signal transduction system serves as an important model system for studying how bacteria adapt to living in densely structured communities known as biofilms. Biofilms frequently cause chronic infections and environmental fouling, and they are very difficult to eradicate. In Pseudomonas aeruginosa, the Wsp system senses contact with a surface, which in turn activates specific genes that promote biofilm formation. We demonstrate that the Wsp system in Burkholderia cenocepacia regulates biofilm formation uniquely from that in Pseudomonas species. Furthermore, a broad phylogenetic analysis reveals the presence of the Wsp system in diverse bacterial species, and sequence analyses of 794 independent systems suggest that the core signaling components function similarly but with key differences that may alter what or how they sense. This study shows that Wsp systems are highly conserved and more broadly distributed than previously thought, and their unique differences likely reflect adaptations to distinct environments.

Comprehensive Genomic Analysis of the Endophytic Bacillus altitudinis Strain GLB197, a Potential Biocontrol Agent of Grape Downy Mildew

Bacillus has been extensively studied for agricultural application as a biocontrol agent. B. altitudinis GLB197, an endophytic bacterium isolated from grape leaves, exhibits distinctive inhibition to grape downy mildew based on unknown mechanisms. To determine the genetic traits involved in the mechanism of biocontrol and host-interaction traits, the genome sequence of GLB197 was obtained and further analyzed. The genome of B. altitudinis GLB197 consisted of one plasmid and a 3,733,835-bp circular chromosome with 41.56% G + C content, containing 3,770 protein-coding genes. Phylogenetic analysis of 17 Bacillus strains using the concatenated 1,226 single-copy core genes divided into different clusters was conducted. In addition, average nucleotide identity (ANI) values indicate that the current taxonomy of some B. pumilus group strains is incorrect. Comparative analysis of B. altitudinis GLB197 proteins with other B. altitudinis strains identified 3,157 core genes. Furthermore, we found that the pan-genome of B. altitudinis is open. The genome of B. altitudinis GLB197 contains one nonribosomal peptide synthetase (NRPS) gene cluster which was annotated as lichenysin. Interestingly, the cluster in B. altitudinis has two more genes than other Bacillus strains (lgrD and lgrB). The two genes were probably obtained via horizontal gene transfer (HGT) during the evolutionary process from Brevibacillus. Taken together, these observations enable the future application of B. altitudinis GLB197 as a biocontrol agent for control of grape downy mildew and promote our understanding of the beneficial interactions between B. altitudinis GLB197 and plants.

Aurantimicrobium photophilum sp. nov., a non-photosynthetic bacterium adjusting its metabolism to the diurnal light cycle and reclassification of Cryobacterium mesophilum as Terrimesophilobacter mesophilus gen. nov., comb. nov

The aerobic primarily chemoorganotrophic actinobacterial strain MWH-Mo1^T was isolated from a freshwater lake and is characterized by small cell lengths of less than 1 µm, small cell volumes of 0.05-0.06 µm³ (ultramicrobacterium), a small genome size of 1.75 Mbp and, at least for an actinobacterium, a low DNA G+C content of 54.6 mol%. Phylogenetic analyses based on concatenated amino acid sequences of 116 housekeeping genes suggested the type strain of Aurantimicrobium minutum affiliated with the family Microbacteriaceae as its closest described relative. Strain MWH-Mo1^T shares with the type strain of that species a 16S rRNA gene sequence similarity of 99.6 % but the genomes of the two strains share an average nucleotide identity of only 79.3 %. Strain MWH-Mo1^T is in many genomic, phenotypic and chemotaxonomic characteristics quite similar to the type strain of A. minutum. Previous intensive investigations revealed two unusual traits of strain MWH-Mo1^T. Although the strain is not known to be phototrophic, the metabolism is adjusted to the diurnal light cycle by up- and down-regulation of genes in light and darkness. This results in faster growth in the presence of light. Additionally, a cell size-independent protection against predation by bacterivorous flagellates, most likely mediated by a proteinaceous cell surface structure, was demonstrated. For the previously intensively investigated aerobic chemoorganotrophic actinobacterial strain MWH-Mo1^T (=CCUG 56426^T=DSM 107758^T), the establishment of the new species Aurantimicrobium photophilum sp. nov. is proposed.

Indirect identification of horizontal gene transfer

Several implicit methods to infer horizontal gene transfer (HGT) focus on pairs of genes that have diverged only after the divergence of the two species in which the genes reside. This situation defines the edge set of a graph, the later-divergence-time (LDT) graph, whose vertices correspond to genes colored by their species. We investigate these graphs in the setting of relaxed scenarios, i.e., evolutionary scenarios that encompass all commonly used variants of duplication-transfer-loss scenarios in the literature. We characterize LDT graphs as a subclass of properly vertex-colored cographs, and provide a polynomial-time recognition algorithm as well as an algorithm to construct a relaxed scenario that explains a given LDT. An edge in an LDT graph implies that the two corresponding genes are separated by at least one HGT event. The converse is not true, however. We show that the complete xenology relation is described by an rs-Fitch graph, i.e., a complete multipartite graph satisfying constraints on the vertex coloring. This class of vertex-colored graphs is also recognizable in polynomial time. We finally address the question “how much information about all HGT events is contained in LDT graphs” with the help of simulations of evolutionary scenarios with a wide range of duplication, loss, and HGT events. In particular, we show that a simple greedy graph editing scheme can be used to efficiently detect HGT events that are implicitly contained in LDT graphs.

Identical sequences found in distant genomes reveal frequent horizontal transfer across the bacterial domain

Horizontal gene transfer (HGT) is an essential force in microbial evolution. Despite detailed studies on a variety of systems, a global picture of HGT in the microbial world is still missing. Here, we exploit that HGT creates long identical DNA sequences in the genomes of distant species, which can be found efficiently using alignment-free methods. Our pairwise analysis of 93,481 bacterial genomes identified 138,273 HGT events. We developed a model to explain their statistical properties as well as estimate the transfer rate between pairs of taxa. This reveals that long-distance HGT is frequent: our results indicate that HGT between species from different phyla has occurred in at least 8% of the species. Finally, our results confirm that the function of sequences strongly impacts their transfer rate, which varies by more than three orders of magnitude between different functional categories. Overall, we provide a comprehensive view of HGT, illuminating a fundamental process driving bacterial evolution.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

SSG-LUGIA: Single Sequence based Genome Level Unsupervised Genomic Island Prediction Algorithm

BACKGROUND

Genomic Islands (GIs) are clusters of genes that are mobilized through horizontal gene transfer. GIs play a pivotal role in bacterial evolution as a mechanism of diversification and adaptation to different niches. Therefore, identification and characterization of GIs in bacterial genomes is important for understanding bacterial evolution. However, quantifying GIs is inherently difficult, and the existing methods suffer from low prediction accuracy and precision-recall trade-off. Moreover, several of them are supervised in nature, and thus, their applications to newly sequenced genomes are riddled with their dependency on the functional annotation of existing genomes.

RESULTS

We present SSG-LUGIA, a completely automated and unsupervised approach for identifying GIs and horizontally transferred genes. SSG-LUGIA is a novel method based on unsupervised anomaly detection technique, accompanied by further refinement using cues from signal processing literature. SSG-LUGIA leverages the atypical compositional biases of the alien genes to localize GIs in prokaryotic genomes. SSG-LUGIA was assessed on a large benchmark dataset `IslandPick' and on a set of 15 well-studied genomes in the literature and followed by a thorough analysis on the well-understood Salmonella typhi CT18 genome. Furthermore, the efficacy of SSG-LUGIA in identifying horizontally transferred genes was evaluated on two additional bacterial genomes, namely, those of Corynebacterium diphtheria NCTC13129 and Pseudomonas aeruginosa LESB58. SSG-LUGIA was examined on draft genomes and was demonstrated to be efficient as an ensemble method.

CONCLUSIONS

Our results indicate that SSG-LUGIA achieved superior performance in comparison to frequently used existing methods. Importantly, it yielded a better trade-off between precision and recall than the existing methods. Its nondependency on the functional annotation of genomes makes it suitable for analyzing newly sequenced, yet uncharacterized genomes. Thus, our study is a significant advance in identification of GIs and horizontally transferred genes. SSG-LUGIA is available as an open source software at https://nibtehaz.github.io/SSG-LUGIA/.

Comparative genome analysis of Corynebacterium species: The underestimated pathogens with high virulence potential

Non-diphtherial Corynebacterium species or diphtheroids were previously considered as the mere contaminants of clinical samples. Of late, they have been reckoned as the formidable infection causing agents of various diseases. While the scientific database is filled with articles that document whole genome analysis of individual isolates, a comprehensive comparative genomic analysis of diphtheroids alongside Corynebacterium diphtheriae is expected to enable us in understanding their genomic as well as evolutionary divergence. Here, we have analysed the whole genome sequences of forty strains that were selected from a range of eleven Corynebacterium species (pathogenic and non-pathogenic). A statistical analysis of the pan and core genomes revealed that even though the core genome is saturated, the pan genome is yet open rendering scope for newer gene families to be accumulated in the course of evolution that might further change the pathogenic behavior of these species. Every strain had bacteriophage components integrated in its genome and some of them were intact and consisted of toxins. The presence of diversified genomic islands was observed across the dataset and most of them consisted of genes for virulence and multidrug resistance. Moreover, the phylogenetic analysis showed that a diphtheroid is the last common ancestor of all the Corynebacterium species. The current study is a compilation of genomic features of pathogenic as well as non-pathogenic Corynebacterium species which provides insights into their virulence potential in the times to come.

Comprehensive and Functional Analysis of Horizontal Gene Transfer Events in Diatoms

Diatoms are a diverse group of mainly photosynthetic algae, responsible for 20% of worldwide oxygen production, which can rapidly respond to favorable conditions and often outcompete other phytoplankton. We investigated the contribution of horizontal gene transfer (HGT) to its ecological success. A large-scale phylogeny-based prokaryotic HGT detection procedure across nine sequenced diatoms showed that 3-5% of their proteome has a horizontal origin and a large influx occurred at the ancestor of diatoms. More than 90% of HGT genes are expressed, and species-specific HGT genes in Phaeodactylum tricornutum undergo strong purifying selection. Genes derived from HGT are implicated in several processes including environmental sensing and expand the metabolic toolbox. Cobalamin (vitamin B12) is an essential cofactor for roughly half of the diatoms and is only produced by bacteria. Five consecutive genes involved in the final synthesis of the cobalamin biosynthetic pathway, which could function as scavenging and repair genes, were detected as HGT. The full suite of these genes was detected in the cold-adapted diatom Fragilariopsis cylindrus. This might give diatoms originating from the Southern Ocean, a region typically depleted in cobalamin, a competitive advantage. Overall, we show that HGT is a prevalent mechanism that is actively used in diatoms to expand its adaptive capabilities.

Evidence for Horizontal and Vertical Transmission of Mtr-Mediated Extracellular Electron Transfer among the Bacteria

Some bacteria and archaea have evolved the means to use extracellular electron donors and acceptors for energy metabolism, a phenomenon broadly known as extracellular electron transfer (EET). One such EET mechanism is the transmembrane electron conduit MtrCAB, which has been shown to transfer electrons derived from metabolic substrates to electron acceptors, like Fe(III) and Mn(IV) oxides, outside the cell. Although most studies of MtrCAB-mediated EET have been conducted in Shewanella oneidensis MR-1, recent investigations in Vibrio and Aeromonas species have revealed that the electron-donating proteins that support MtrCAB in Shewanella are not as representative as previously thought. This begs the question of how widespread the capacity for MtrCAB-mediated EET is, the changes it has accrued in different lineages, and where these lineages persist today. Here, we employed a phylogenetic and comparative genomics approach to identify the MtrCAB system across all domains of life. We found mtrCAB in the genomes of numerous diverse Bacteria from a wide range of environments, and the patterns therein strongly suggest that mtrCAB was distributed through both horizontal and subsequent vertical transmission, and with some cases indicating downstream modular diversification of both its core and accessory components. Our data point to an emerging evolutionary story about metal-oxidizing and -reducing metabolism, demonstrates that this capacity for EET has broad relevance to a diversity of taxa and the biogeochemical cycles they drive, and lays the foundation for further studies to shed light on how this mechanism may have coevolved with Earth's redox landscape. IMPORTANCE While many metabolisms make use of soluble, cell-permeable substrates like oxygen or hydrogen, there are other substrates, like iron or manganese, that cannot be brought into the cell. Some bacteria and archaea have evolved the means to directly "plug in" to such environmental electron reservoirs in a process known as extracellular electron transfer (EET), making them powerful agents of biogeochemical change and promising vehicles for bioremediation and alternative energy. Yet the diversity, distribution, and evolution of EET mechanisms are poorly constrained. Here, we present findings showing that the genes encoding one such EET system (mtrCAB) are present in a broad diversity of bacteria found in a wide range of environments, emphasizing the ubiquity and potential impact of EET in our biosphere. Our results suggest that these genes have been disseminated largely through horizontal transfer, and the changes they have accrued in these lineages potentially reflect adaptations to changing environments.