Modular networks and cumulative impact of lateral transfer in prokaryote genome evolution

Chemical fertilizers promote dissemination of ARGs in maize rhizosphere: An overlooked risk revealed after 37-year traditional agriculture practice

Bacterial communities in soil and rhizosphere maintain a large collection of antibiotic resistance genes (ARGs). However, few of these ARGs and antibiotic resistant bacteria (ARB) are well-characterized under traditional farming practices. Here we compared the ARG profiles of maize rhizosphere and their bulk soils using metagenomic analysis to identify the ARG dissemination and explored the potential impact of chemical fertilization on ARB. Results showed a relatively lower abundance but higher diversity of ARGs under fertilization than straw-return. Moreover, the abundance and diversity of MGEs were significantly promoted by chemical fertilizer inputs in the rhizosphere compared to bulk soil. Machine learning and bipartite networks identified three bacterial genera (Pseudomonas, Bacillus and Streptomyces) as biomarkers for ARG accumulation. Thus we cultured 509 isolates belonging to these three genera from the rhizosphere and tested their antimicrobial susceptibility, and found that multi-resistance was frequently observed among Pseudomonas isolates. Assembly-based tracking explained that ARGs and four class I integrons (LR134330, LS998783, CP065848, LT883143) were co-occurred among contigs from Pseudomonas sp. Chemical fertilizers may shape the resistomes of maize rhizosphere, highlighting that rhizosphere carried multidrug-resistant Pseudomonas isolates, which may pose a risk to animal and human health. This study adds knowledge of long-term chemical fertilization on ARG dissemination in farmland systems and provides information for decision-making in agricultural production and monitoring.

The radical impact of oxygen on prokaryotic evolution-enzyme inhibition first, uninhibited essential biosyntheses second, aerobic respiration third

Molecular oxygen is a stable diradical. All O2-dependent enzymes employ a radical mechanism. Generated by cyanobacteria, O2 started accumulating on Earth 2.4 billion years ago. Its evolutionary impact is traditionally sought in respiration and energy yield. We mapped 365 O2-dependent enzymatic reactions of prokaryotes to phylogenies for the corresponding 792 protein families. The main physiological adaptations imparted by O2-dependent enzymes were not energy conservation, but novel organic substrate oxidations and O2-dependent, hence O2-tolerant, alternative pathways for O2-inhibited reactions. Oxygen-dependent enzymes evolved in ancestrally anaerobic pathways for essential cofactor biosynthesis including NAD+, pyridoxal, thiamine, ubiquinone, cobalamin, heme, and chlorophyll. These innovations allowed prokaryotes to synthesize essential cofactors in O2-containing environments, a prerequisite for the later emergence of aerobic respiratory chains.

Tandem gene duplication selected by activation of horizontally transferred gene in bacteria

Horizontal gene transfer occurs frequently in bacteria, but the mechanism driving activation and optimization of the expression of horizontally transferred genes (HTGs) in new recipient strains is not clear. Our previous study found that spontaneous tandem DNA duplication resulted in rapid activation of HTGs. Here, we took advantage of this finding to develop a novel technique for tandem gene duplication, named tandem gene duplication selected by activation of horizontally transferred gene in bacteria (TDAH), in which tandem duplication was selected by the activation of horizontally transferred selectable marker gene. TDAH construction does not contain any reported functional elements based on homologous or site-specific recombination and DNA amplification. TDAH only contains an essential selectable marker for copy number selection and 9-bp-microhomology border sequences for precise illegitimate recombination. One transformation and 3 days were enough to produce a high-copy strain, so its procedure is simple and fast. Without subsequent knockout of the endogenous recombination system, TDAH could also generate the relatively stable high-copy tandem duplication for plasmid-carried and genome-integrated DNA. TDAH also showed an excellent capacity for increase gene expression and worked well in different industrial bacteria. We also applied TDAH to select the optimal high copy number of ribA for vitamin B2 production in E. coli; the yield was improved by 3.5 times and remained stable even after 12 subcultures. TDAH is a useful tool for recombinant protein production and expression optimization of biosynthetic pathways. KEY POINTS: • We develop a novel and efficient technique (TDAH) for tandem gene duplication in bacterium. TDAH is based on the mechanism of HTG rapid activation. TDAH does not contain any reported functional elements based on homologous recombination and DNA amplification. TDAH only contains an essential selectable marker for copy number selection, so its construction and procedure are very simple and fast. • TDAH is the first reported selected and stable tandem-gene-duplication technique in which the selected high-copy plasmid-carried and genome-integrated DNA could remain stable without the subsequent knockout of recombination system. • TDAH showed an excellent capacity for regulating gene expression and worked well in different industrial bacteria, indicating it is a useful tool for recombinant protein production and expression optimization of biosynthetic pathways. • TDAH was applied to select the optimal high copy number of ribA for vitamin B2 production in E. coli; the yield was improved by 3.5-fold and remained stable even after 12 subcultures.

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.

Differential responses of the gut microbiome and resistome to antibiotic exposures in infants and adults

Despite their crucial importance for human health, there is still relatively limited knowledge on how the gut resistome changes or responds to antibiotic treatment across ages, especially in the latter case. Here, we use fecal metagenomic data from 662 Danish infants and 217 young adults to fill this gap. The gut resistomes are characterized by a bimodal distribution driven by E. coli composition. The typical profile of the gut resistome differs significantly between adults and infants, with the latter distinguished by higher gene and plasmid abundances. However, the predominant antibiotic resistance genes (ARGs) are the same. Antibiotic treatment reduces bacterial diversity and increased ARG and plasmid abundances in both cohorts, especially core ARGs. The effects of antibiotic treatments on the gut microbiome last longer in adults than in infants, and different antibiotics are associated with distinct impacts. Overall, this study broadens our current understanding of gut resistome dynamics and the impact of antibiotic treatment across age groups.

Elucidating applied voltage on the fate of antibiotic resistance genes in microbial electrolysis cell: Focusing on its transmission between anolyte and biofilm microbes

Microbial electrolysis cell (MEC) system to treat wastewater containing antibiotics has been researched actively in past years. However, the fate of antibiotic resistant genes (ARGs) in MEC is not fully revealed. The effect of applied voltage on the migration of ARGs between anolyte and biofilm microbes via examining the microbial physiology and abundances of macrolide resistance genes (MRGs) and mobile genetic elements (MGEs) was elucidated in this research. Results showed that the abundance of MRGs and MGEs was decreased in the anolyte, but their abundances were increased on the electrode biofilm, indicating their transmission from anolyte to biofilm microbes. Increased applied voltage enhanced adenosine triphosphate (ATP), reactive oxygen species (ROS), and cell membrane permeability of electrode microorganisms. The structure of the electrode microbial community was shifted through applied voltage, and the abundance of electroactive microorganisms (Geobacter, Azospirillum and Dechlorobacter) was significantly improved. Network analysis revealed that Geobacter and Geothrix were potential hosts for MRGs. Therefore, the horizontal and vertical gene transfer of ARGs could be increased by the applied voltage, leading to the enriched ARGs at the electrode biofilm. This study provides evidence and insights into the transmission of ARGs between anolyte and biofilm microbes in MEC system. SYNOPSIS: This study revealed the effect of applied voltage on ARGs in MEC and the potential migration mechanism of ARGs.

Extracellular Polysaccharide Receptor and Receptor-Binding Proteins of the Rhodobacter capsulatus Bacteriophage-like Gene Transfer Agent RcGTA

A variety of prokaryotes produce a bacteriophage-like gene transfer agent (GTA), and the alphaproteobacterial Rhodobacter capsulatus RcGTA is a model GTA. Some environmental isolates of R. capsulatus lack the ability to acquire genes transferred by the RcGTA (recipient capability). In this work, we investigated the reason why R. capsulatus strain 37b4 lacks recipient capability. The RcGTA head spike fiber and tail fiber proteins have been proposed to bind extracellular oligosaccharide receptors, and strain 37b4 lacks a capsular polysaccharide (CPS). The reason why strain 37b4 lacks a CPS was unknown, as was whether the provision of a CPS to 37b4 would result in recipient capability. To address these questions, we sequenced and annotated the strain 37b4 genome and used BLAST interrogations of this genome sequence to search for homologs of genes known to be needed for R. capsulatus recipient capability. We also created a cosmid-borne genome library from a wild-type strain, mobilized the library into 37b4, and used the cosmid-complemented strain 37b4 to identify genes needed for a gain of function, allowing for the acquisition of RcGTA-borne genes. The relative presence of CPS around a wild-type strain, 37b4, and cosmid-complemented 37b4 cells was visualized using light microscopy of stained cells. Fluorescently tagged head spike fiber and tail fiber proteins of the RcGTA particle were created and used to measure the relative binding to wild-type and 37b4 cells. We found that strain 37b4 lacks recipient capability because of an inability to bind RcGTA; the reason it is incapable of binding is that it lacks CPS, and the absence of CPS is due to the absence of genes previously shown to be needed for CPS production in another strain. In addition to the head spike fiber, we found that the tail fiber protein also binds to the CPS.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

Adding context to the pneumococcal core genes using bioinformatic analysis of the intergenic pangenome of Streptococcus pneumoniae

Introduction: Whole genome sequencing offers great opportunities for linking genotypes to phenotypes aiding in our understanding of human disease and bacterial pathogenicity. However, these analyses often overlook non-coding intergenic regions (IGRs). By disregarding the IGRs, crucial information is lost, as genes have little biological function without expression. Methods/Results: In this study, we present the first complete pangenome of the important human pathogen Streptococcus pneumoniae (pneumococcus), spanning both the genes and IGRs. We show that the pneumococcus species retains a small core genome of IGRs that are present across all isolates. Gene expression is highly dependent on these core IGRs, and often several copies of these core IGRs are found across each genome. Core genes and core IGRs show a clear linkage as 81% of core genes are associated with core IGRs. Additionally, we identify a single IGR within the core genome that is always occupied by one of two highly distinct sequences, scattered across the phylogenetic tree. Discussion: Their distribution indicates that this IGR is transferred between isolates through horizontal regulatory transfer independent of the flanking genes and that each type likely serves different regulatory roles depending on their genetic context.

The ESKAPE mobilome contributes to the spread of antimicrobial resistance and CRISPR-mediated conflict between mobile genetic elements

Mobile genetic elements (MGEs) mediate the shuffling of genes among organisms. They contribute to the spread of virulence and antibiotic resistance (AMR) genes in human pathogens, such as the particularly problematic group of ESKAPE pathogens. Here, we performed the first systematic analysis of MGEs, including plasmids, prophages, and integrative and conjugative/mobilizable elements (ICEs/IMEs), across all ESKAPE pathogens. We found that different MGE types are asymmetrically distributed across these pathogens, and that most horizontal gene transfer (HGT) events are restricted by phylum or genus. We show that the MGEs proteome is involved in diverse functional processes and distinguish widespread proteins within the ESKAPE context. Moreover, anti-CRISPRs and AMR genes are overrepresented in the ESKAPE mobilome. Our results also underscore species-specific trends shaping the number of MGEs, AMR, and virulence genes across pairs of conspecific ESKAPE genomes with and without CRISPR-Cas systems. Finally, we observed that CRISPR spacers found on prophages, ICEs/IMEs, and plasmids have different targeting biases: while plasmid and prophage CRISPRs almost exclusively target other plasmids and prophages, respectively, ICEs/IMEs CRISPRs preferentially target prophages. Overall, our study highlights the general importance of the ESKAPE mobilome in contributing to the spread of AMR and mediating conflict among MGEs.

Improved Fixed-Parameter Algorithm for the Tree Containment Problem on Unrooted Phylogenetic Network

Phylogenetic trees are unable to represent the evolutionary process for a collection of species if reticulation events happened, and a generalized model named phylogenetic network was introduced consequently. However, the representation of the evolutionary process for one gene is actually a phylogenetic tree that is "contained" in the phylogenetic network for the considered species containing the gene. Thus a fundamental computational problem named Tree Containment problem arises, which asks whether a phylogenetic tree is contained in a phylogenetic network. The previous research on the problem mainly focused on its rooted version of which the considered tree and network are rooted, and several algorithms were proposed when the considered network is binary or structure-restricted. There is almost no algorithm for its unrooted version except the recent fixed-parameter algorithm with runtime O(4^kn²), where k and n are the reticulation number and size of the considered unrooted binary phylogenetic network N, respectively. As the runtime is a little expensive when considering big values of k, we aim to improve it and successfully propose a fixed-parameter algorithm with runtime O(2.594^kn²) in the paper. Additionally, we experimentally show its effectiveness on biological data and simulated data.

A Novel Module Promotes Horizontal Gene Transfer in Azorhizobium caulinodans ORS571

Azorhizobium caulinodans ORS571 contains an 87.6 kb integrative and conjugative element (ICE^Ac) that conjugatively transfers symbiosis genes to other rhizobia. Many hypothetical redundant gene fragments (rgfs) are abundant in ICE^Ac, but their potential function in horizontal gene transfer (HGT) is unknown. Molecular biological methods were employed to delete hypothetical rgfs, expecting to acquire a minimal ICE^Ac and consider non-functional rgfs as editable regions for inserting genes related to new symbiotic functions. We determined the significance of rgf4 in HGT and identified the physiological function of genes designated rihF1a (AZC_3879), rihF1b (AZC_RS26200), and rihR (AZC_3881). In-frame deletion and complementation assays revealed that rihF1a and rihF1b work as a unit (rihF1) that positively affects HGT frequency. The EMSA assay and lacZ-based reporter system showed that the XRE-family protein RihR is not a regulator of rihF1 but promotes the expression of the integrase (intC) that has been reported to be upregulated by the LysR-family protein, AhaR, through sensing host’s flavonoid. Overall, a conservative module containing rihF1 and rihR was characterized, eliminating the size of ICE^Ac by 18.5%. We propose the feasibility of constructing a minimal ICE^Ac element to facilitate the exchange of new genetic components essential for symbiosis or other metabolic functions between soil bacteria.

Rapid acquisition of microorganisms and microbial genes can help explain punctuated evolution

The punctuated mode of evolution posits that evolution occurs in rare bursts of rapid evolutionary change followed by long periods of genetic stability (stasis). The accepted cause for the rapid changes in punctuated evolution is special ecological circumstances – selection forces brought about by changes in the environment. This article presents a complementary explanation for punctuated evolution by the rapid formation of genetic variants in animals and plants by the acquisition of microorganisms from the environment into microbiomes and microbial genes into host genomes by horizontal gene transfer. Several examples of major evolutionary events driven by microorganisms are discussed, including the formation of the first eukaryotic cell, the ability of some animals to digest cellulose and other plant cell-wall complex polysaccharides, dynamics of root system architecture, and the formation of placental mammals. These changes by cooperation were quantum leaps in the evolutionary development of complex bilolgical systems and can contribute to an understanding of the mechanisms underlying punctuated evolution.

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.

Realistic Gene Transfer to Gene Duplication Ratios Identify Different Roots in the Bacterial Phylogeny Using a Tree Reconciliation Method

The rooting of phylogenetic trees permits important inferences about ancestral states and the polarity of evolutionary events. Recently, methods that reconcile discordance between gene-trees and species-trees—tree reconciliation methods—are becoming increasingly popular for rooting species trees. Rooting via reconciliation requires values for a particular parameter, the gene transfer to gene duplication ratio (T:D), which in current practice is estimated on the fly from discordances observed in the trees. To date, the accuracy of T:D estimates obtained by reconciliation analyses has not been compared to T:D estimates obtained by independent means, hence the effect of T:D upon inferences of species tree roots is altogether unexplored. Here we investigated the issue in detail by performing tree reconciliations of more than 10,000 gene trees under a variety of T:D ratios for two phylogenetic cases: a bacterial (prokaryotic) tree with 265 species and a fungal-metazoan (eukaryotic) tree with 31 species. We show that the T:D ratios automatically estimated by a current tree reconciliation method, ALE, generate virtually identical T:D ratios across bacterial genes and fungal-metazoan genes. The T:D ratios estimated by ALE differ 10- to 100-fold from robust, ALE-independent estimates from real data. More important is our finding that the root inferences using ALE in both datasets are strongly dependent upon T:D. Using more realistic T:D ratios, the number of roots inferred by ALE consistently increases and, in some cases, clearly incorrect roots are inferred. Furthermore, our analyses reveal that gene duplications have a far greater impact on ALE’s preferences for phylogenetic root placement than gene transfers or gene losses do. Overall, we show that obtaining reliable species tree roots with ALE is only possible when gene duplications are abundant in the data and the number of falsely inferred gene duplications is low. Finding a sufficient sample of true gene duplications for rooting species trees critically depends on the T:D ratios used in the analyses. T:D ratios, while being important parameters of genome evolution in their own right, affect the root inferences with tree reconciliations to an unanticipated degree.

A network approach to decipher the dynamics of Lysobacteraceae plasmid gene sharing

Plasmids provide an efficient vehicle for gene sharing among bacterial populations, playing a key role in bacterial evolution. Network approaches are particularly suitable to represent multipartite relationships and are useful tools to characterize plasmid-mediated gene sharing events. The Lysobacteraceae bacterial family gathers plant commensal, plant pathogenic and opportunistic human pathogens for which plasmid mediated adaptation was reported. We searched for homologues of plasmid gene sequences from this family in all the diversity of available bacterial genome sequences and built a network of plasmid gene sharing from the results. While plasmid genes are openly shared between the bacteria of the Lysobacteraceae family, taxonomy strongly defined the boundaries of these exchanges, that only barely reached other families. Most inferred plasmid gene sharing events involved a few genes only, and evidence of full plasmid transfers were restricted to taxonomically close taxon. We detected multiple plasmid-chromosome gene transfers, among which the otherwise known sharing of a heavy metal resistance transposon. In the network, bacterial lifestyles shaped sub-structures of isolates colonizing specific ecological niches and harboring specific types of resistance genes. Genes associated to pathogenicity or antibiotic and metal resistance were among those that most importantly structured the network, highlighting the imprints of human-mediated selective pressure on pathogenic populations. A massive sequencing effort on environmental Lysobacteraceae is therefore required to refine our understanding on how this reservoir fuels the emergence and the spread of genes amongst this family and its potential impact on plant, animal and human health.

Major episodes of horizontal gene transfer drove the evolution of land plants

How horizontal gene transfer (HGT) has contributed to the evolution of animals and plants remains a major puzzle. Despite recent progress, defining the overall scale and pattern of HGT events in land plants has been largely elusive. In this study, we performed systematic analyses for acquired genes in different plant groups and throughout land plant evolution. We found that relatively recent HGT events occurred in charophytes and all major land plant groups, but their frequency declined rapidly in seed plants. Two major episodes of HGT events occurred in land plant evolution, corresponding to the early evolution of streptophytes and the origin of land plants, respectively. Importantly, a vast majority of the genes acquired in the two episodes have been retained in descendant groups, affecting numerous activities and processes of land plants. We analyzed some of the acquired genes involved in stress responses, ion and metabolite transport, growth and development, and specialized metabolism, and further assessed the cumulative effects of HGT in land plants.

The gossip paradox: Why do bacteria share genes?

Bacteria, in contrast to eukaryotic cells, contain two types of genes: chromosomal genes that are fixed to the cell, and plasmids, smaller loops of DNA capable of being passed from one cell to another. The sharing of plasmid genes between individual bacteria and between bacterial lineages has contributed vastly to bacterial evolution, allowing specialized traits to 'jump ship' between one lineage or species and the next. The benefits of this generosity from the point of view of both recipient cell and plasmid are generally understood: plasmids receive new hosts and ride out selective sweeps across the population, recipient cells gain new traits (such as antibiotic resistance). Explaining this behavior from the point of view of donor cells is substantially more difficult. Donor cells pay a fitness cost in order to share plasmids, and run the risk of sharing advantageous genes with their competition and rendering their own lineage redundant, while seemingly receiving no benefit in return. Using both compartment based models and agent based simulations we demonstrate that 'secretive' genes which restrict horizontal gene transfer are favored over a wide range of models and parameter values, even when sharing carries no direct cost. 'Generous' chromosomal genes which are more permissive of plasmid transfer are found to have neutral fitness at best, and are generally disfavored by selection. Our findings lead to a peculiar paradox: given the obvious benefits of keeping secrets, why do bacteria share information so freely?

Evolving perspective on the origin and diversification of cellular life and the virosphere

The tree of life (TOL) is a powerful framework to depict the evolutionary history of cellular organisms through time, from our microbial origins to the diversification of multicellular eukaryotes that shape the visible biosphere today. During the past decades, our perception of the TOL has fundamentally changed, in part, due to profound methodological advances, which allowed a more objective approach to studying organismal and viral diversity and led to the discovery of major new branches in the TOL as well as viral lineages. Phylogenetic and comparative genomics analyses of these data have, among others, revolutionized our understanding of the deep roots and diversity of microbial life, the origin of the eukaryotic cell, eukaryotic diversity, as well as the origin, and diversification of viruses. In this review, we provide an overview of some of the recent discoveries on the evolutionary history of cellular organisms and their viruses and discuss a variety of complementary techniques that we consider crucial for making further progress in our understanding of the TOL and its interconnection with the virosphere.

Population Genomics of Microbial Biostalactites: Non-recombinogenic Genome Islands and Microdiversification by Transposons

Intrapopulation genetic variability in prokaryotes is receiving increasing attention thanks to improving sequencing methods; however, the ability to distinguish intrapopulation variability from species clusters or initial stages of gene flow barrier development remains insufficient. To overcome this limitation, we took advantage of the lifestyle of Ferrovum myxofaciens, a species that may represent 99% of prokaryotic microbiome of biostalactites growing at acid mine drainage springs. We gained four complete and one draft metagenome-assembled F. myxofaciens genomes using Oxford Nanopore and Illumina sequencing and mapped the reads from each sample on the reference genomes to assess the intrapopulation variability. We observed two phenomena associated with intrapopulation variability: hypervariable regions affected by mobilome expansion called “scrapyards,” and variability in gene disruptions caused by transposons within each population. Both phenomena were previously described in prokaryotes. However, we present here for the first time scrapyard regression and the development of a new one. Nearly complete loss of intrapopulation short sequence variability in the old scrapyard and high variability in the new one suggest that localized gene flow suppression is necessary for scrapyard formation. Concerning the variable gene disruptions, up to 9 out of 41 occurrences per sample were located in highly conserved diguanylate cyclases/phosphodiesterases. We propose that microdiversification of life strategies may be an adaptive outcome of random diguanylate cyclase elimination. The mine biostalactites thus proved as a unique model system for describing genomic intrapopulation processes, as they offer easily sampleable units enriched in a single microbial species.

Contamination detection in genomic data: more is not enough

The decreasing cost of sequencing and concomitant augmentation of publicly available genomes have created an acute need for automated software to assess genomic contamination. During the last 6 years, 18 programs have been published, each with its own strengths and weaknesses. Deciding which tools to use becomes more and more difficult without an understanding of the underlying algorithms. We review these programs, benchmarking six of them, and present their main operating principles. This article is intended to guide researchers in the selection of appropriate tools for specific applications. Finally, we present future challenges in the developing field of contamination detection.

Modular evolution of secretion systems and virulence plasmids in a bacterial species complex

BACKGROUND

Many named species as defined in current bacterial taxonomy correspond to species complexes. Uncertainties regarding the organization of their genetic diversity challenge research efforts. We utilized the Agrobacterium tumefaciens species complex (a.k.a. Agrobacterium biovar 1), a taxon known for its phytopathogenicity and applications in transformation, as a study system and devised strategies for investigating genome diversity and evolution of species complexes.

RESULTS

We utilized 35 genome assemblies, including 14 newly generated ones, to achieve a phylogenetically balanced sampling of A. tumefaciens. Our genomic analysis suggested that the 10 genomospecies described previously are distinct biological species and supported a quantitative guideline for species delineation. Furthermore, our inference of gene content and core-genome phylogeny allowed for investigations of genes critical in fitness and ecology. For the type VI secretion system (T6SS) involved in interbacterial competition and thought to be conserved, we detected multiple losses and one horizontal gene transfer. For the tumor-inducing plasmids (pTi) and pTi-encoded type IV secretion system (T4SS) that are essential for agrobacterial phytopathogenicity, we uncovered novel diversity and hypothesized their involvement in shaping this species complex. Intriguingly, for both T6SS and T4SS, genes encoding structural components are highly conserved, whereas extensive diversity exists for genes encoding effectors and other proteins.

CONCLUSIONS

We demonstrate that the combination of a phylogeny-guided sampling scheme and an emphasis on high-quality assemblies provides a cost-effective approach for robust analysis in evolutionary genomics. We show that the T6SS VgrG proteins involved in specific effector binding and delivery can be classified into distinct types based on domain organization. The co-occurrence patterns of VgrG-associated domains and the neighboring genes that encode different chaperones/effectors can be used to infer possible interacting partners. Similarly, the associations between plant host preference and the pTi type among these strains can be used to infer phenotype-genotype correspondence. Our strategies for multi-level investigations at scales that range from whole genomes to intragenic domains and phylogenetic depths from between- to within-species are applicable to other bacteria. Furthermore, modularity observed in the molecular evolution of genes and domains is useful for inferring functional constraints and informing experimental works.

Assembling a Reference Phylogenomic Tree of Bacteria and Archaea by Summarizing Many Gene Phylogenies

Phylogenomics is the inference of phylogenetic trees based on multiple marker genes sampled in the genomes of interest. An important challenge in phylogenomics is the potential incongruence among the evolutionary histories of individual genes, which can be widespread in microorganisms due to the prevalence of horizontal gene transfer. This protocol introduces the procedures for building a phylogenetic tree of a large number of microbial genomes using a broad sampling of marker genes that are representative of whole-genome evolution. The protocol highlights the use of a gene tree summary method, which can effectively reconstruct the species tree while accounting for the topological conflicts among individual gene trees. The pipeline described in this protocol is scalable to tens of thousands of genomes while retaining high accuracy. We discussed multiple software tools, libraries, and scripts to enable convenient adoption of the protocol. The protocol is suitable for microbiology and microbiome studies based on public genomes and metagenomic data.

Metaplasmidome-encoded functions of Siberian low-centered polygonal tundra soils

Plasmids have the potential to transfer genetic traits within bacterial communities and thereby serve as a crucial tool for the rapid adaptation of bacteria in response to changing environmental conditions. Our knowledge of the environmental pool of plasmids (the metaplasmidome) and encoded functions is still limited due to a lack of sufficient extraction methods and tools for identifying and assembling plasmids from metagenomic datasets. Here, we present the first insights into the functional potential of the metaplasmidome of permafrost-affected active-layer soil-an environment with a relatively low biomass and seasonal freeze-thaw cycles that is strongly affected by global warming. The obtained results were compared with plasmid-derived sequences extracted from polar metagenomes. Metaplasmidomes from the Siberian active layer were enriched via cultivation, which resulted in a longer contig length as compared with plasmids that had been directly retrieved from the metagenomes of polar environments. The predicted hosts of plasmids belonged to Moraxellaceae, Pseudomonadaceae, Enterobacteriaceae, Pectobacteriaceae, Burkholderiaceae, and Firmicutes. Analysis of their genetic content revealed the presence of stress-response genes, including antibiotic and metal resistance determinants, as well as genes encoding protectants against the cold.

Altered gut ecosystems plus the microbiota's potential for rapid evolution: A recipe for inevitable change with unknown consequences

In a single human gut, which is estimated to produce 1000-times more bacteria in a single day than the entire human population on Earth as of 2020, the potential for evolution is vast. In addition to the sheer volume of reproductive events, prokaryotes can transfer most genes horizontally, greatly accelerating their potential to evolve. In the face of this evolutionary potential, Westernization has led to profound changes in the ecosystem of the gut, including increased chronic inflammation in many individuals and dramatically reduced fiber consumption and decreased seasonal variation in the diet of most individuals. Experimental work using a variety of model systems has shown that bacteria will evolve within days to weeks when faced with substantial environmental changes. However, studies evaluating the effects of inflammation of the gut on the microbiota are still in their infancy and generally confounded by the effects of the microbiota on the immune system. At the same time, experimental data indicate that complete loss of fiber from the diet constitutes an extinction-level event for the gut microbiota. However, these studies evaluating diet may not apply to Westernized humans who typically have reduced but not absent levels of fiber in their diet. Thus, while it is expected that the microbiota will evolve rapidly in the face of Westernization, experimental studies that address the magnitude of that evolution are generally lacking, and it remains unknown to what extent this evolutionary process affects disease and the ability to treat the disease state.

Gene Duplications Are At Least 50 Times Less Frequent than Gene Transfers in Prokaryotic Genomes

The contribution of gene duplications to the evolution of eukaryotic genomes is well studied. By contrast, studies of gene duplications in prokaryotes are scarce and generally limited to a handful of genes or careful analysis of a few prokaryotic lineages. Systematic broad-scale studies of prokaryotic genomes that sample available data are lacking, leaving gaps in our understanding of the contribution of gene duplications as a source of genetic novelty in the prokaryotic world. Here, we report conservative and robust estimates for the frequency of recent gene duplications within prokaryotic genomes relative to recent lateral gene transfer (LGT), as mechanisms to generate multiple copies of related sequences in the same genome. We obtain our estimates by focusing on evolutionarily recent events among 5,655 prokaryotic genomes, thereby avoiding vagaries of deep phylogenetic inference and confounding effects of ancient events and differential loss. We find that recent, genome-specific gene duplications are at least 50 times less frequent and probably 100 times less frequent than recent, genome-specific, gene acquisitions via LGT. The frequency of gene duplications varies across lineages and functional categories. The findings improve our understanding of genome evolution in prokaryotes and have far-reaching implications for evolutionary models that entail LGT to gene duplications ratio as a parameter.

Photosynthetic Systems Suggest an Evolutionary Pathway to Diderms

Bacteria are divided primarily into monoderms (with one cell membrane, and usually Gram-positive, due to a thick peptidoglycan layer) and diderms (with two cell membranes, and mostly Gram-negative, due to a thin peptidoglycan layer sandwiched between the two membranes). Photosynthetic species are spread among the taxonomic groups, some having type I reaction centers (RCI in monoderm phylum Firmicutes; and diderm phyla Acidobacteria and Chlorobi), others with type II reaction centers (RCII in monoderm phylum Chloroflexi; and diderm taxa Gemmatimonadetes, and alpha-, beta-, and gamma-Proteobacteria), and some containing both (RCI and RCII, only in diderm phylum Cyanobacteria). In most bacterial phylograms, photosystem types and diderm taxa are polyphyletic. A more parsimonious arrangement, which is supported by photosystem evolution, as well as additional sets of molecular characters, suggests that endosymbiotic events resulted in the formation of the diderms. In the model presented, monoderms readily form a monophyletic group, while diderms are produced by at least two endosymbiotic events, followed by additional evolutionary changes.

Promoter Architecture Differences among Alphaproteobacteria and Other Bacterial Taxa

Much of our knowledge of bacterial transcription initiation has been derived from studying the promoters of Escherichia coli and Bacillus subtilis. Given the expansive diversity across the bacterial phylogeny, it is unclear how much of this knowledge can be applied to other organisms. Here, we report on bioinformatic analyses of promoter sequences of the primary σ factor (σ70) by leveraging publicly available transcription start site (TSS) sequencing data sets for nine bacterial species spanning five phyla. This analysis identifies previously unreported differences in the -35 and -10 elements of σ70-dependent promoters in several groups of bacteria. We found that Actinobacteria and Betaproteobacteria σ70-dependent promoters lack the TTG triad in their -35 element, which is predicted to be conserved across the bacterial phyla. In addition, the majority of the Alphaproteobacteria σ70-dependent promoters analyzed lacked the thymine at position -7 that is highly conserved in other phyla. Bioinformatic examination of the Alphaproteobacteria σ70-dependent promoters identifies a significant overrepresentation of essential genes and ones encoding proteins with common cellular functions downstream of promoters containing an A, C, or G at position -7. We propose that transcription of many σ70-dependent promoters in Alphaproteobacteria depends on the transcription factor CarD, which is an essential protein in several members of this phylum. Our analysis expands the knowledge of promoter architecture across the bacterial phylogeny and provides new information that can be used to engineer bacteria for use in medical, environmental, agricultural, and biotechnological processes. IMPORTANCE Transcription of DNA to RNA by RNA polymerase is essential for cells to grow, develop, and respond to stress. Understanding the process and control of transcription is important for health, disease, the environment, and biotechnology. Decades of research on a few bacteria have identified promoter DNA sequences that are recognized by the σ subunit of RNA polymerase. We used bioinformatic analyses to reveal previously unreported differences in promoter DNA sequences across the bacterial phylogeny. We found that many Actinobacteria and Betaproteobacteria promoters lack a sequence in their -35 DNA recognition element that was previously assumed to be conserved and that Alphaproteobacteria lack a thymine residue at position -7, also previously assumed to be conserved. Our work reports important new information about bacterial transcription, illustrates the benefits of studying bacteria across the phylogenetic tree, and proposes new lines of future investigation.

Prokaryotic taxonomy and nomenclature in the age of big sequence data

The classification of life forms into a hierarchical system (taxonomy) and the application of names to this hierarchy (nomenclature) is at a turning point in microbiology. The unprecedented availability of genome sequences means that a taxonomy can be built upon a comprehensive evolutionary framework, a longstanding goal of taxonomists. However, there is resistance to adopting a single framework to preserve taxonomic freedom, and ever increasing numbers of genomes derived from uncultured prokaryotes threaten to overwhelm current nomenclatural practices, which are based on characterised isolates. The challenge ahead then is to reach a consensus on the taxonomic framework and to adapt and scale the existing nomenclatural code, or create a new code, to systematically incorporate uncultured taxa into the chosen framework.

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.

Microbial driven genetic variation in holobionts

Genetic variation in holobionts, (host and microbiome), occurring by changes in both host and microbiome genomes, can be observed from two perspectives: observable variations and the processes that bring about the variation. The observable includes the enormous genetic diversity of prokaryotes, which gave rise to eukaryotic organisms. Holobionts then evolved a rich microbiome with a stable core containing essential genes, less so common taxa, and a more diverse non-core enabling considerable genetic variation. The result being that, the human gut microbiome, for example, contains 1,000 times more unique genes than are present in the human genome. Microbial driven genetic variation processes in holobionts include: (1) Acquisition of novel microbes from the environment, which bring in multiple genes in one step, (2) amplification/reduction of certain microbes in the microbiome, that contribute to holobiont` s adaptation to changing conditions, (3) horizontal gene transfer between microbes and between microbes and host, (4) mutation, which plays an important role in optimizing interactions between different microbiota and between microbiota and host. We suggest that invertebrates and plants, where microbes can live intracellularly, have a greater chance of genetic exchange between microbiota and host, thus a greater chance of vertical transmission and a greater effect of microbiome on evolution of host than vertebrates. However, even in vertebrates the microbiome can aid in environmental fluctuations by amplification/reduction and by acquisition of novel microorganisms.

Horizontal Gene Transfer Involving Chloroplasts

Horizontal gene transfer (HGT)- is defined as the acquisition of genetic material from another organism. However, recent findings indicate a possible role of HGT in the acquisition of traits with adaptive significance, suggesting that HGT is an important driving force in the evolution of eukaryotes as well as prokaryotes. It has been noted that, in eukaryotes, HGT is more prevalent than originally thought. Mitochondria and chloroplasts lost a large number of genes after their respective endosymbiotic events occurred. Even after this major content loss, organelle genomes still continue to lose their own genes. Many of these are subsequently acquired by intracellular gene transfer from the original plastid. The aim of our review was to elucidate the role of chloroplasts in the transfer of genes. This review also explores gene transfer involving mitochondrial and nuclear genomes, though recent studies indicate that chloroplast genomes are far more active in HGT as compared to these other two DNA-containing cellular compartments.

Extensive genome-wide duplications in the eastern oyster (Crassostrea virginica)

Genomic structural variation is an important source of genetic and phenotypic diversity, playing a critical role in evolution. The recent availability of a high-quality reference genome for the eastern oyster, Crassostrea virginica, and whole-genome sequence data of samples from across the species range in the USA, provides an opportunity to explore structural variation across the genome of this species. Our analysis shows significantly greater individual-level duplications of regions across the genome than that of most model vertebrate species. Duplications are widespread across all ten chromosomes with variation in frequency per chromosome. The eastern oyster shows a large interindividual variation in duplications as well as particular chromosomal regions with a higher density of duplications. A high percentage of duplications seen in C. virginica lie completely within genes and exons, suggesting the potential for impacts on gene function. These results support the hypothesis that structural changes may play a significant role in standing genetic variation in C. virginica, and potentially have a role in their adaptive and evolutionary success. Altogether, these results suggest that copy number variation plays an important role in the genomic variation of C. virginica. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Ten Years of Collaborative Progress in the Quest for Orthologs

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

The metabolic network of the last bacterial common ancestor

Bacteria are the most abundant cells on Earth. They are generally regarded as ancient, but due to striking diversity in their metabolic capacities and widespread lateral gene transfer, the physiology of the first bacteria is unknown. From 1089 reference genomes of bacterial anaerobes, we identified 146 protein families that trace to the last bacterial common ancestor, LBCA, and form the conserved predicted core of its metabolic network, which requires only nine genes to encompass all universal metabolites. Our results indicate that LBCA performed gluconeogenesis towards cell wall synthesis, and had numerous RNA modifications and multifunctional enzymes that permitted life with low gene content. In accordance with recent findings for LUCA and LACA, analyses of thousands of individual gene trees indicate that LBCA was rod-shaped and the first lineage to diverge from the ancestral bacterial stem was most similar to modern Clostridia, followed by other autotrophs that harbor the acetyl-CoA pathway.

Contribution of horizontal gene transfer to the functionality of microbial biofilm on a macroalgae

Horizontal gene transfer (HGT) is thought to be an important driving force for microbial evolution and niche adaptation and has been show in vitro to occur frequently in biofilm communities. However, the extent to which HGT takes place and what functions are being transferred in more complex and natural biofilm systems remains largely unknown. To address this issue, we investigated here HGT and enrichment of gene functions in the biofilm community of the common kelp (macroalgae) Ecklonia radiata in comparison to microbial communities in the surrounding seawater. We found that HGTs in the macroalgal biofilms were dominated by transfers between bacterial members of the same class or order and frequently involved genes for nutrient transport, sugar and phlorotannin degradation as well as stress responses, all functions that would be considered beneficial for bacteria living in this particular niche. HGT did not appear to be driven by mobile gene elements, indicating rather an involvement of unspecific DNA uptake (e.g. natural transformation). There was also a low overlap between the gene functions subject to HGT and those enriched in the biofilm community in comparison to planktonic community members. This indicates that much of the functionality required for bacteria to live in an E. radiata biofilm might be derived from vertical or environmental transmissions of symbionts. This study enhances our understanding of the relative role of evolutionary and ecological processes in driving community assembly and genomic diversity of biofilm communities.

Deep learning for HGT insertion sites recognition

Background

Horizontal Gene Transfer (HGT) refers to the sharing of genetic materials between distant species that are not in a parent-offspring relationship. The HGT insertion sites are important to understand the HGT mechanisms. Recent studies in main agents of HGT, such as transposon and plasmid, demonstrate that insertion sites usually hold specific sequence features. This motivates us to find a method to infer HGT insertion sites according to sequence features.

Results

In this paper, we propose a deep residual network, DeepHGT, to recognize HGT insertion sites. To train DeepHGT, we extracted about 1.55 million sequence segments as training instances from 262 metagenomic samples, where the ratio between positive instances and negative instances is about 1:1. These segments are randomly partitioned into three subsets: 80% of them as the training set, 10% as the validation set, and the remaining 10% as the test set. The training loss of DeepHGT is 0.4163 and the validation loss is 0.423. On the test set, DeepHGT has achieved the area under curve (AUC) value of 0.8782. Furthermore, in order to further evaluate the generalization of DeepHGT, we constructed an independent test set containing 689,312 sequence segments from another 147 gut metagenomic samples. DeepHGT has achieved the AUC value of 0.8428, which approaches the previous test AUC value. As a comparison, the gradient boosting classifier model implemented in PyFeat achieve an AUC value of 0.694 and 0.686 on the above two test sets, respectively. Furthermore, DeepHGT could learn discriminant sequence features; for example, DeepHGT has learned a sequence pattern of palindromic subsequences as a significantly (P-value=0.0182) local feature. Hence, DeepHGT is a reliable model to recognize the HGT insertion site.

Conclusion

DeepHGT is the first deep learning model that can accurately recognize HGT insertion sites on genomes according to the sequence pattern.

A spectrum of verticality across genes

Lateral gene transfer (LGT) has impacted prokaryotic genome evolution, yet the extent to which LGT compromises vertical evolution across individual genes and individual phyla is unknown, as are the factors that govern LGT frequency across genes. Estimating LGT frequency from tree comparisons is problematic when thousands of genomes are compared, because LGT becomes difficult to distinguish from phylogenetic artefacts. Here we report quantitative estimates for verticality across all genes and genomes, leveraging a well-known property of phylogenetic inference: phylogeny works best at the tips of trees. From terminal (tip) phylum level relationships, we calculate the verticality for 19,050,992 genes from 101,422 clusters in 5,655 prokaryotic genomes and rank them by their verticality. Among functional classes, translation, followed by nucleotide and cofactor biosynthesis, and DNA replication and repair are the most vertical. The most vertically evolving lineages are those rich in ecological specialists such as Acidithiobacilli, Chlamydiae, Chlorobi and Methanococcales. Lineages most affected by LGT are the α-, β-, γ-, and δ- classes of Proteobacteria and the Firmicutes. The 2,587 eukaryotic clusters in our sample having prokaryotic homologues fail to reject eukaryotic monophyly using the likelihood ratio test. The low verticality of α-proteobacterial and cyanobacterial genomes requires only three partners-an archaeal host, a mitochondrial symbiont, and a plastid ancestor-each with mosaic chromosomes, to directly account for the prokaryotic origin of eukaryotic genes. In terms of phylogeny, the 100 most vertically evolving prokaryotic genes are neither representative nor predictive for the remaining 97% of an average genome. In search of factors that govern LGT frequency, we find a simple but natural principle: Verticality correlates strongly with gene distribution density, LGT being least likely for intruding genes that must replace a preexisting homologue in recipient chromosomes. LGT is most likely for novel genetic material, intruding genes that encounter no competing copy.

phyloFlash: Rapid Small-Subunit rRNA Profiling and Targeted Assembly from Metagenomes

To track organisms across all domains of life, the SSU rRNA gene is the gold standard. Many environmental microbes are known only from high-throughput sequence data, but the SSU rRNA gene, the key to visualization by molecular probes and link to existing literature, is often missing from metagenome-assembled genomes (MAGs). The easy-to-use phyloFlash software suite tackles this gap with rapid, SSU rRNA-centered taxonomic classification, targeted assembly, and graph-based linking to MAGs. Starting from a cleaned reference database, phyloFlash profiles the taxonomic diversity and assembles the sorted SSU rRNA reads. The phyloFlash design is domain agnostic and covers eukaryotes, archaea, and bacteria alike. phyloFlash also provides utilities to visualize multisample comparisons and to integrate the recovered SSU rRNAs in a metagenomics workflow by linking them to MAGs using assembly graph parsing.

The small-subunit rRNA (SSU rRNA) gene is the key marker in molecular ecology for all domains of life, but it is largely absent from metagenome-assembled genomes that often are the only resource available for environmental microbes. Here, we present phyloFlash, a pipeline to overcome this gap with rapid, SSU rRNA-centered taxonomic classification, targeted assembly, and graph-based binning of full metagenomic assemblies. We show that a cleanup of artifacts is pivotal even with a curated reference database. With such a filtered database, the general-purpose mapper BBmap extracts SSU rRNA reads five times faster than the rRNA-specialized tool SortMeRNA with similar sensitivity and higher selectivity on simulated metagenomes. Reference-based targeted assemblers yielded either highly fragmented assemblies or high levels of chimerism, so we employ the general-purpose genomic assembler SPAdes. Our optimized implementation is independent of reference database composition and has satisfactory levels of chimera formation. phyloFlash quickly processes Illumina (meta)genomic data, is straightforward to use, even as part of high-throughput quality control, and has user-friendly output reports. The software is available at https://github.com/HRGV/phyloFlash (GPL3 license) and is documented with an online manual.

IMPORTANCE To track organisms across all domains of life, the SSU rRNA gene is the gold standard. Many environmental microbes are known only from high-throughput sequence data, but the SSU rRNA gene, the key to visualization by molecular probes and link to existing literature, is often missing from metagenome-assembled genomes (MAGs). The easy-to-use phyloFlash software suite tackles this gap with rapid, SSU rRNA-centered taxonomic classification, targeted assembly, and graph-based linking to MAGs. Starting from a cleaned reference database, phyloFlash profiles the taxonomic diversity and assembles the sorted SSU rRNA reads. The phyloFlash design is domain agnostic and covers eukaryotes, archaea, and bacteria alike. phyloFlash also provides utilities to visualize multisample comparisons and to integrate the recovered SSU rRNAs in a metagenomics workflow by linking them to MAGs using assembly graph parsing.

Whitefly genomes contain ribotoxin coding genes acquired from plants

Ribosome inactivating proteins (RIPs) are RNA N-glycosidases that depurinate a specific adenine residue in the conserved sarcin/ricin loop of 28S rRNA. These enzymes are widely distributed among plants and bacteria. Previously, we have described for the first time RIP genes in mosquitoes belonging to the Culicidae family. We showed that these genes are derived from a single event of horizontal gene transfer (HGT) from a prokaryotic donor. Mosquito RIP genes are evolving under purifying selection, strongly suggesting that these toxins have acquired a functional role. In this work, we show the existence of two RIP encoding genes in the genome of the whitefly Bemisia tabaci, a hemiptera species belonging to the Aleyrodidae family distantly related to mosquitoes. Contamination artifacts were ruled out analyzing three independent B. tabaci genome databases. In contrast to mosquito RIPs, whitefly genes harbor introns and according to transcriptomic evidence are transcribed and spliced. Phylogeny and the taxonomic distribution strongly support that whitefly RIP genes are derived from an independent HGT event from a plant source. These results, along with our previous description of RIPs in Diptera, suggest that the acquired genes are functional in these insects and confer some fitness advantage.

Batch-Learning Self-Organizing Map Identifies Horizontal Gene Transfer Candidates and Their Origins in Entire Genomes

Horizontal gene transfer (HGT) has been widely suggested to play a critical role in the environmental adaptation of microbes; however, the number and origin of the genes in microbial genomes obtained through HGT remain unknown as the frequency of detected HGT events is generally underestimated, particularly in the absence of information on donor sequences. As an alternative to phylogeny-based methods that rely on sequence alignments, we have developed an alignment-free clustering method on the basis of an unsupervised neural network “Batch-Learning Self-Organizing Map (BLSOM)” in which sequence fragments are clustered based solely on oligonucleotide similarity without taxonomical information, to detect HGT candidates and their origin in entire genomes. By mapping the microbial genomic sequences on large-scale BLSOMs constructed with nearly all prokaryotic genomes, HGT candidates can be identified, and their origin assigned comprehensively, even for microbial genomes that exhibit high novelty. By focusing on two types of Alphaproteobacteria, specifically psychrotolerant Sphingomonas strains from an Antarctic lake, we detected HGT candidates using BLSOM and found higher proportions of HGT candidates from organisms belonging to Betaproteobacteria in the genomes of these two Antarctic strains compared with those of continental strains. Further, an origin difference was noted in the HGT candidates found in the two Antarctic strains. Although their origins were highly diversified, gene functions related to the cell wall or membrane biogenesis were shared among the HGT candidates. Moreover, analyses of amino acid frequency suggested that housekeeping genes and some HGT candidates of the Antarctic strains exhibited different characteristics to other continental strains. Lys, Ser, Thr, and Val were the amino acids found to be increased in the Antarctic strains, whereas Ala, Arg, Glu, and Leu were decreased. Our findings strongly suggest a low-temperature adaptation process for microbes that may have arisen convergently as an independent evolutionary strategy in each Antarctic strain. Hence, BLSOM analysis could serve as a powerful tool in not only detecting HGT candidates and their origins in entire genomes, but also in providing novel perspectives into the environmental adaptations of microbes.

Exploration of space to achieve scientific breakthroughs

Living organisms adapt to changing environments using their amazing flexibility to remodel themselves by a process called evolution. Environmental stress causes selective pressure and is associated with genetic and phenotypic shifts for better modifications, maintenance, and functioning of organismal systems. The natural evolution process can be used in complement to rational strain engineering for the development of desired traits or phenotypes as well as for the production of novel biomaterials through the imposition of one or more selective pressures. Space provides a unique environment of stressors (e.g., weightlessness and high radiation) that organisms have never experienced on Earth. Cells in the outer space reorganize and develop or activate a range of molecular responses that lead to changes in cellular properties. Exposure of cells to the outer space will lead to the development of novel variants more efficiently than on Earth. For instance, natural crop varieties can be generated with higher nutrition value, yield, and improved features, such as resistance against high and low temperatures, salt stress, and microbial and pest attacks. The review summarizes the literature on the parameters of outer space that affect the growth and behavior of cells and organisms as well as complex colloidal systems. We illustrate an understanding of gravity-related basic biological mechanisms and enlighten the possibility to explore the outer space environment for application-oriented aspects. This will stimulate biological research in the pursuit of innovative approaches for the future of agriculture and health on Earth.

Characterization of horizontally acquired ribotoxin encoding genes and their transcripts in Aedes aegypti

Ribosome Inactivating Proteins (RIPs) are RNA N-glycosidases that depurinate a specific adenine residue in the conserved sarcin/ricin loop of the 28S rRNA. The occurrence of RIP genes has been described in a wide range of plant taxa, as well as in several species of bacteria and fungi. A remarkable case is the presence of these genes in metazoans belonging to the Culicinae subfamily. We reported that these genes are derived from a single horizontal gene transfer event, most likely from a bacterial donor species. Moreover, we have shown evidence that mosquito RIP genes are evolving under purifying selection, suggesting that these toxins have acquired a functional role in these organisms. In the present work, we characterized the intra-specific sequence variability of Aedes aegypti RIP genes (RIPAe1, RIPAe2, and RIPAe3) and tested their expression at the mRNA level. Our results show that RIPAe2 and RIPAe3 are transcribed and polyadenylated, and their expression levels are modulated across the developmental stages. Varibility among genes was observed, including the existence of null alleles for RIPAe1 and RIPAe2, with variants showing partial deletions. These results further support the existence of a physiological function for these foreign genes in mosquitoes. The possible nature of this functionality is discussed.

Large-scale network analysis captures biological features of bacterial plasmids

Many bacteria can exchange genetic material through horizontal gene transfer (HGT) mediated by plasmids and plasmid-borne transposable elements. Here, we study the population structure and dynamics of over 10,000 bacterial plasmids, by quantifying their genetic similarities and reconstructing a network based on their shared k-mer content. We use a community detection algorithm to assign plasmids into cliques, which correlate with plasmid gene content, bacterial host range, GC content, and existing classifications based on replicon and mobility (MOB) types. Further analysis of plasmid population structure allows us to uncover candidates for yet undescribed replicon genes, and to identify transposable elements as the main drivers of HGT at broad phylogenetic scales. Our work illustrates the potential of network-based analyses of the bacterial 'mobilome' and opens up the prospect of a natural, exhaustive classification framework for bacterial plasmids.

Delineation of a Subgroup of the Genus Paraburkholderia, Including P. terrae DSM 17804T, P. hospita DSM 17164T, and Four Soil-Isolated Fungiphiles, Reveals Remarkable Genomic and Ecological Features-Proposal for the Definition of a P. hospita Species Cluster

The fungal-interactive (fungiphilic) strains BS001, BS007, BS110, and BS437 have previously been preliminarily assigned to the species Paraburkholderia terrae. However, in the (novel) genus Paraburkholderia, an as-yet unresolved subgroup exists, that clusters around Paraburkholderia hospita (containing the species P. terrae, P. hospita, and Paraburkholderia caribensis). To shed light on the precise relationships across the respective type strains and the novel fungiphiles, we here compare their genomic and ecophysiological features. To reach this goal, the genomes of the three type strains, with sizes ranging from 9.0 to 11.5 Mb, were de novo sequenced and the high-quality genomes analyzed. Using whole-genome, ribosomal RNA and marker-gene-concatenate analyses, close relationships between P. hospita DSM 17164T and P. terrae DSM 17804T, versus more remote relationships to P. caribensis DSM 13236T, were found. All four fungiphilic strains clustered closely to the two-species cluster. Analyses of average nucleotide identities (ANIm) and tetranucleotide frequencies (TETRA) confirmed the close relationships between P. hospita DSM 17164T and P. terrae DSM 17804T (ANIm = 95.42; TETRA = 0.99784), as compared with the similarities of each one of these strains to P. caribensis DSM 13236T. A species cluster was thus proposed. Furthermore, high similarities of the fungiphilic strains BS001, BS007, BS110, and BS437 with this cluster were found, indicating that these strains also make part of it, being closely linked to P. hospita DSM 17164T (ANIm = 99%; TETRA = 0.99). We propose to coin this cluster the P. hospita species cluster (containing P. hospita DSM 17164T, P. terrae DSM 17804T, and strains BS001, BS007, BS110, and BS437), being clearly divergent from the closely related species P. caribensis (type strain DSM 13236T). Moreover, given their close relatedness to P. hospita DSM 17164T within the cluster, we propose to rename the four fungiphilic strains as members of P. hospita. Analysis of migratory behavior along with fungal growth through soil revealed both P. terrae DSM 17804T and P. hospita DSM 17164T (next to the four fungiphilic strains) to be migration-proficient, whereas P. caribensis DSM 13236T was a relatively poor migrator. Examination of predicted functions across the genomes of the seven investigated strains, next to several selected additional ones, revealed the common presence of features in the P. hospita cluster strains that are potentially important in interactions with soil fungi. Thus, genes encoding specific metabolic functions, biofilm formation (pelABCDEFG, pgaABCD, alginate-related genes), motility/chemotaxis, type-4 pili, and diverse secretion systems were found.

Horizontal Gene Transfer and Endophytes: An Implication for the Acquisition of Novel Traits

Horizontal gene transfer (HGT), an important evolutionary mechanism observed in prokaryotes, is the transmission of genetic material across phylogenetically distant species. In recent years, the availability of complete genomes has facilitated the comprehensive analysis of HGT and highlighted its emerging role in the adaptation and evolution of eukaryotes. Endophytes represent an ecologically favored association, which highlights its beneficial attributes to the environment, in agriculture and in healthcare. The HGT phenomenon in endophytes, which features an important biological mechanism for their evolutionary adaptation within the host plant and simultaneously confers "novel traits" to the associated microbes, is not yet completely understood. With a focus on the emerging implications of HGT events in the evolution of biological species, the present review discusses the occurrence of HGT in endophytes and its socio-economic importance in the current perspective. To our knowledge, this review is the first report that provides a comprehensive insight into the impact of HGT in the adaptation and evolution of endophytes.

The evolution of MarR family transcription factors as counter-silencers in regulatory networks

Gene duplication facilitates the evolution of biological complexity, as one copy of a gene retains its original function while a duplicate copy can acquire mutations that would otherwise diminish fitness. Duplication has played a particularly important role in the evolution of regulatory networks by permitting novel regulatory interactions and responses to stimuli. The diverse MarR family of transcription factors (MFTFs) illustrate this concept, ranging from highly specific repressors of single operons to pleiotropic global regulators controlling hundreds of genes. MFTFs are often genetically and functionally linked to antimicrobial efflux systems. However, the SlyA MFTF lineage in the Enterobacteriaceae plays little or no role in regulating efflux but rather functions as transcriptional counter-silencers, which alleviate xenogeneic silencing of horizontally acquired genes and facilitate bacterial evolution by horizontal gene transfer. This review will explore recent advances in our understanding of MFTF traits that have contributed to their functional evolution.