Horizontal gene transfer in prokaryotes: quantification and classification

Hot Spots of Site-Specific Integration into the Sinorhizobium meliloti Chromosome

The diversity of phage-related sequences (PRSs) and their site-specific integration into the genomes of nonpathogenic, agriculturally valuable, nitrogen-fixing root nodule bacteria, such as Sinorhizobium meliloti, were evaluated in this study. A total of 314 PRSs, ranging in size from 3.24 kb to 88.98 kb, were identified in the genomes of 27 S. meliloti strains. The amount of genetic information foreign to S. meliloti accumulated in all identified PRSs was 6.30 Mb. However, more than 53% of this information was contained in prophages (Phs) and genomic islands (GIs) integrated into genes encoding tRNAs (tRNA genes) located on the chromosomes of the rhizobial strains studied. It was found that phiLM21-like Phs were predominantly abundant in the genomes of S. meliloti strains of distant geographical origin, whereas RR1-A- and 16-3-like Phs were much less common. In addition, GIs predominantly contained fragments of phages infecting bacteria of distant taxa, while rhizobiophage-like sequences were unique. A site-specific integration analysis revealed that not all tRNA genes in S. meliloti are integration sites, but among those in which integration occurred, there were "hot spots" of integration into which either Phs or GIs were predominantly inserted. For the first time, it is shown that at these integration "hot spots", not only is the homology of attP and attB strictly preserved, but integrases in PRSs similar to those of phages infecting the Proteobacteria genera Azospirillum or Pseudomonas are also present. The data presented greatly expand the understanding of the fate of phage-related sequences in host bacterial genomes and also raise new questions about the role of phages in bacterial-phage coevolution.

Discovering genomic islands in unannotated bacterial genomes using sequence embedding

Motivation

Genomic islands (GEIs) are clusters of genes in bacterial genomes that are typically acquired by horizontal gene transfer. GEIs play a crucial role in the evolution of bacteria by rapidly introducing genetic diversity and thus helping them adapt to changing environments. Specifically of interest to human health, many GEIs contain pathogenicity and antimicrobial resistance genes. Detecting GEIs is, therefore, an important problem in biomedical and environmental research. There have been many previous studies for computationally identifying GEIs. Still, most of these studies rely on detecting anomalies in the unannotated nucleotide sequences or on a fixed set of known features on annotated nucleotide sequences.

Results

Here, we present TreasureIsland, which uses a new unsupervised representation of DNA sequences to predict GEIs. We developed a high-precision boundary detection method featuring an incremental fine-tuning of GEI borders, and we evaluated the accuracy of this framework using a new comprehensive reference dataset, Benbow. We show that TreasureIsland's accuracy rivals other GEI predictors, enabling efficient and faster identification of GEIs in unannotated bacterial genomes.

Availability and implementation

TreasureIsland is available under an MIT license at: https://github.com/FriedbergLab/GenomicIslandPrediction.

Epistasis, core-genome disharmony, and adaptation in recombining bacteria

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome's co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background-facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.IMPORTANCEGenetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter, an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.

Horizontal gene transfer in eukaryotes: aligning theory with data

Horizontal gene transfer (HGT), or lateral gene transfer, is the non-sexual movement of genetic information between genomes. It has played a pronounced part in bacterial and archaeal evolution, but its role in eukaryotes is less clear. Behaviours unique to eukaryotic cells - phagocytosis and endosymbiosis - have been proposed to increase the frequency of HGT, but nuclear genomes encode fewer HGTs than bacteria and archaea. Here, I review the existing theory in the context of the growing body of data on HGT in eukaryotes, which suggests that any increased chance of acquiring new genes through phagocytosis and endosymbiosis is offset by a reduced need for these genes in eukaryotes, because selection in most eukaryotes operates on variation not readily generated by HGT.

The contours of evolution: In defence of Darwin's tree of life paradigm

Both the concept of a Darwinian tree of life (TOL) and the possibility of its accurate reconstruction have been much criticized. Criticisms mostly revolve around the extensive occurrence of lateral gene transfer (LGT), instances of uptake of complete organisms to become organelles (with the associated subsequent gene transfer to the nucleus), as well as the implications of more subtle aspects of the biological species concept. Here we argue that none of these criticisms are sufficient to abandon the valuable TOL concept and the biological realities it captures. Especially important is the need to conceptually distinguish between organismal trees and gene trees, which necessitates incorporating insights of widely occurring LGT into modern evolutionary theory. We demonstrate that all criticisms, while based on important new findings, do not invalidate the TOL. After considering the implications of these new insights, we find that the contours of evolution are best represented by a TOL.

Integrating taxonomic signals from MAGs and contigs improves read annotation and taxonomic profiling of metagenomes

Metagenomic analysis typically includes read-based taxonomic profiling, assembly, and binning of metagenome-assembled genomes (MAGs). Here we integrate these steps in Read Annotation Tool (RAT), which uses robust taxonomic signals from MAGs and contigs to enhance read annotation. RAT reconstructs taxonomic profiles with high precision and sensitivity, outperforming other state-of-the-art tools. In high-diversity groundwater samples, RAT annotates a large fraction of the metagenomic reads, calling novel taxa at the appropriate, sometimes high taxonomic ranks. Thus, RAT integrative profiling provides an accurate and comprehensive view of the microbiome from shotgun metagenomics data. The package of Contig Annotation Tool (CAT), Bin Annotation Tool (BAT), and RAT is available at https://github.com/MGXlab/CAT_pack (from CAT pack v6.0). The CAT pack now also supports Genome Taxonomy Database (GTDB) annotations.

Bacterial and viral co-infections in aquaculture under climate warming: co-evolutionary implications, diagnosis, and treatment

Climate change and the associated environmental temperature fluctuations are contributing to increases in the frequency and severity of disease outbreaks in both wild and farmed aquatic species. This has a significant impact on biodiversity and also puts global food production systems, such as aquaculture, at risk. Most infections are the result of complex interactions between multiple pathogens, and understanding these interactions and their co-evolutionary mechanisms is crucial for developing effective diagnosis and control strategies. In this review, we discuss current knowledge on bacteria-bacteria, virus-virus, and bacterial and viral co-infections in aquaculture as well as their co-evolution in the context of global warming. We also propose a framework and different novel methods (e.g. advanced molecular tools such as digital PCR and next-generation sequencing) to (1) precisely identify overlooked co-infections, (2) gain an understanding of the co-infection dynamics and mechanisms by knowing species interactions, and (3) facilitate the development multi-pathogen preventive measures such as polyvalent vaccines. As aquaculture disease outbreaks are forecasted to increase both due to the intensification of practices to meet the protein demand of the increasing global population and as a result of global warming, understanding and treating co-infections in aquatic species has important implications for global food security and the economy.

Structure, function and substrate preferences of archaeal S-adenosyl-L-homocysteine hydrolases

S-Adenosyl-L-homocysteine hydrolase (SAHH) reversibly cleaves S-adenosyl-L-homocysteine, the product of S-adenosyl-L-methionine-dependent methylation reactions. The conversion of S-adenosyl-L-homocysteine into adenosine and L-homocysteine plays an important role in the regulation of the methyl cycle. An alternative metabolic route for S-adenosyl-L-methionine regeneration in the extremophiles Methanocaldococcus jannaschii and Thermotoga maritima has been identified, featuring the deamination of S-adenosyl-L-homocysteine to S-inosyl-L-homocysteine. Herein, we report the structural characterisation of different archaeal SAHHs together with a biochemical analysis of various SAHHs from all three domains of life. Homologues deriving from the Euryarchaeota phylum show a higher conversion rate with S-inosyl-L-homocysteine compared to S-adenosyl-L-homocysteine. Crystal structures of SAHH originating from Pyrococcus furiosus in complex with SLH and inosine as ligands, show architectural flexibility in the active site and offer deeper insights into the binding mode of hypoxanthine-containing substrates. Altogether, the findings of our study support the understanding of an alternative metabolic route for S-adenosyl-L-methionine and offer insights into the evolutionary progression and diversification of SAHHs involved in methyl and purine salvage pathways.

Evolution and variation in amide aminoacyl-tRNA synthesis

The amide proteogenic amino acids, asparagine and glutamine, are two of the twenty amino acids used in translation by all known life. The aminoacyl-tRNA synthetases for asparagine and glutamine, asparaginyl-tRNA synthetase and glutaminyl tRNA synthetase, evolved after the split in the last universal common ancestor of modern organisms. Before that split, life used two-step indirect pathways to synthesize asparagine and glutamine on their cognate tRNAs to form the aminoacyl-tRNA used in translation. These two-step pathways were retained throughout much of the bacterial and archaeal domains of life and eukaryotic organelles. The indirect routes use non-discriminating aminoacyl-tRNA synthetases (non-discriminating aspartyl-tRNA synthetase and non-discriminating glutamyl-tRNA synthetase) to misaminoacylate the tRNA. The misaminoacylated tRNA formed is then transamidated into the amide aminoacyl-tRNA used in protein synthesis by tRNA-dependent amidotransferases (GatCAB and GatDE). The enzymes and tRNAs involved assemble into complexes known as transamidosomes to help maintain translational fidelity. These pathways have evolved to meet the varied cellular needs across a diverse set of organisms, leading to significant variation. In certain bacteria, the indirect pathways may provide a means to adapt to cellular stress by reducing the fidelity of protein synthesis. The retention of these indirect pathways versus acquisition of asparaginyl-tRNA synthetase and glutaminyl tRNA synthetase in lineages likely involves a complex interplay of the competing uses of glutamine and asparagine beyond translation, energetic costs, co-evolution between enzymes and tRNA, and involvement in stress response that await further investigation.

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.

Gene Transfer-Based Phylogenetics: Analytical Expressions and Additivity via Birth-Death Theory

The genomic era has opened up vast opportunities in molecular systematics, one of which is deciphering the evolutionary history in fine detail. Under this mass of data, analyzing the point mutations of standard markers is often too crude and slow for fine-scale phylogenetics. Nevertheless, genome dynamics (GD) events provide alternative, often richer information. The synteny index (SI) between a pair of genomes combines gene order and gene content information, allowing the comparison of genomes of unequal gene content, together with order considerations of their common genes. Recently, genome dynamics has been modeled as a continuous-time Markov process, and gene distance in the genome as a birth-death-immigration process. Nevertheless, due to complexities arising in this setting, no precise and provably consistent estimators could be derived, resulting in heuristic solutions. Here, we extend this modeling approach by using techniques from birth-death theory to derive explicit expressions of the system's probabilistic dynamics in the form of rational functions of the model parameters. This, in turn, allows us to infer analytically accurate distances between organisms based on their SI. Subsequently, we establish additivity of this estimated evolutionary distance (a desirable property yielding phylogenetic consistency). Applying the new measure in simulation studies shows that it provides accurate results in realistic settings and even under model extensions such as gene gain/loss or over a tree structure. In the real-data realm, we applied the new formulation to unique data structure that we constructed-the ordered orthology DB-based on a new version of the EggNOG database, to construct a tree with more than 4.5K taxa. To the best of our knowledge, this is the largest gene-order-based tree constructed and it overcomes shortcomings found in previous approaches. Constructing a GD-based tree allows to confirm and contrast findings based on other phylogenetic approaches, as we show.

Sperm-dependent asexual species and their role in ecology and evolution

Sexual reproduction is the primary mode of reproduction in eukaryotes, but some organisms have evolved deviations from classical sex and switched to asexuality. These asexual lineages have sometimes been viewed as evolutionary dead ends, but recent research has revealed their importance in many areas of general biology. Our review explores the understudied, yet important mechanisms by which sperm-dependent asexuals that produce non-recombined gametes but rely on their fertilization, can have a significant impact on the evolution of coexisting sexual species and ecosystems. These impacts are concentrated around three major fields. Firstly, sperm-dependent asexuals can potentially impact the gene pool of coexisting sexual species by either restricting their population sizes or by providing bridges for interspecific gene flow whose type and consequences substantially differ from gene flow mechanisms expected under sexual reproduction. Secondly, they may impact on sexuals' diversification rates either directly, by serving as stepping-stones in speciation, or indirectly, by promoting the formation of pre- and postzygotic reproduction barriers among nascent species. Thirdly, they can potentially impact on spatial distribution of species, via direct or indirect (apparent) types of competition and Allee effects. For each such mechanism, we provide empirical examples of how natural sperm-dependent asexuals impact the evolution of their sexual counterparts. In particular, we highlight that these broad effects may last beyond the tenure of the individual asexual lineages causing them, which challenges the traditional perception that asexual lineages are short-lived evolutionary dead ends and minor sideshows. Our review also proposes new research directions to incorporate the aforementioned impacts of sperm-dependent asexuals. These research directions will ultimately enhance our understanding of the evolution of genomes and biological interactions in general.

Revisiting the Role of Genetic Variation in Adaptation

AbstractEvolutionary biologists have thought about the role of genetic variation during adaptation for a very long time-before we understood the organization of the genetic code, the provenance of genetic variation, and how such variation influenced the phenotypes on which natural selection acts. Half a century after the discovery of the structure of DNA and the unraveling of the genetic code, we have a rich understanding of these problems and the means to both delve deeper and widen our perspective across organisms and natural populations. The 2022 Vice Presidential Symposium of the American Society of Naturalists highlighted examples of recent insights into the role of genetic variation in adaptive processes, which are compiled in this special section. The work was conducted in different parts of the world, included theoretical and empirical studies with diverse organisms, and addressed distinct aspects of how genetic variation influences adaptation. In our introductory article to the special section, we discuss some important recent insights about the generation and maintenance of genetic variation, its impacts on phenotype and fitness, its fate in natural populations, and its role in driving adaptation. By placing the special section articles in the broader context of recent developments, we hope that this overview will also serve as a useful introduction to the field.

Stepwise pathway for early evolutionary assembly of dissimilatory sulfite and sulfate reduction

Microbial dissimilatory sulfur metabolism utilizing dissimilatory sulfite reductases (Dsr) influenced the biochemical sulfur cycle during Earth's history and the Dsr pathway is thought to be an ancient metabolic process. Here we performed comparative genomics, phylogenetic, and synteny analyses of several Dsr proteins involved in or associated with the Dsr pathway across over 195,000 prokaryotic metagenomes. The results point to an archaeal origin of the minimal DsrABCMK(N) protein set, having as primordial function sulfite reduction. The acquisition of additional Dsr proteins (DsrJOPT) increased the Dsr pathway complexity. Archaeoglobus would originally possess the archaeal-type Dsr pathway and the archaeal DsrAB proteins were replaced with the bacterial reductive-type version, possibly at the same time as the acquisition of the QmoABC and DsrD proteins. Further inventions of two Qmo complex types, which are more spread than previously thought, allowed microorganisms to use sulfate as electron acceptor. The ability to use the Dsr pathway for sulfur oxidation evolved at least twice, with Chlorobi and Proteobacteria being extant descendants of these two independent adaptations.

Unrecognized risk of perfluorooctane sulfonate in promoting conjugative transfers of bacterial antibiotic resistance genes

Antibiotic resistance is a major global health crisis facing humanity, with horizontal gene transfer (HGT) as a principal dissemination mechanism in the natural and clinical environments. Perfluoroalkyl substances (PFASs) are emerging contaminants of global concern due to their high persistence in the environment and adverse effects on humans. However, it is unknown whether PFASs affect the HGT of bacterial antibiotic resistance. Using a genetically engineered Escherichia coli MG1655 as the donor of plasmid-encoded antibiotic resistance genes (ARGs), E. coli J53 and soil bacterial community as two different recipients, this study demonstrated that the conjugation frequency of ARGs between two E. coli strains was (1.45 ± 0.17) × 10-5 and perfluorooctane sulfonate (PFOS) at environmentally relevant concentrations (2-50 μg L-1) increased conjugation transfer between E. coli strains by up to 3.25-fold. Increases in reactive oxygen species production, cell membrane permeability, biofilm formation capacity, and cell contact in two E. coli strains were proposed as major promotion mechanisms from PFOS exposure. Weighted gene co-expression network analysis of transcriptome data identified a series of candidate genes whose expression changes could contribute to the increase in conjugation transfer induced by PFOS. Furthermore, PFOS also generally increased the ARG transfer into the studied soil bacterial community, although the uptake ability of different community members of the plasmid either increased or decreased upon PFOS exposure depending on specific bacterial taxa. Overall, this study reveals an unrecognized risk of PFOS in accelerating the dissemination of antibiotic resistance. IMPORTANCE Perfluoroalkyl substances (PFASs) are emerging contaminants of global concern due to their high persistence in the environment and adverse health effects. Although the influence of environmental pollutants on the spread of antibiotic resistance, one of the biggest threats to global health, has attracted increasing attention in recent years, it is unknown whether environmental residues of PFASs affect the dissemination of bacterial antibiotic resistance. Considering PFASs, often called "forever" compounds, have significantly higher environmental persistence than most emerging organic contaminants, exploring the effect of PFASs on the spread of antibiotic resistance is more environmentally relevant and has essential ecological and health significance. By systematically examining the influence of perfluorooctane sulfonate on the antibiotic resistance gene conjugative transfer, not only at the single-strain level but also at the community level, this study has uncovered an unrecognized risk of PFASs in promoting conjugative transfers of bacterial antibiotic resistance genes, which could be incorporated into the risk assessment framework of PFASs.

Genomic and clinical characterization of Klebsiella pneumoniae carrying the pks island

Background

The pks island and its production of the bacterial secondary metabolite genotoxin, colibactin, have attracted increasing attention. However, genomic articles focusing on pks islands in Klebsiella pneumoniae, as well as comparative genomic studies of mobile genetic elements, such as prophages, plasmids, and insertion sequences, are lacking. In this study, a large-scale analysis was conducted to understand the prevalence and evolution of pks islands, differences in mobile genetic elements between pks-negative and pks-positive K. pneumoniae, and clinical characteristics of infection caused by pks-positive K. pneumoniae.

Methods

The genomes of 2,709 K. pneumoniae were downloaded from public databases, among which, 1,422 were from NCBI and 1,287 were from the China National GeneBank DataBase (CNGBdb). Screening for virulence and resistance genes, phylogenetic tree construction, and pan-genome analysis were performed. Differences in mobile genetic elements between pks-positive and pks-negative strains were compared. The clinical characteristics of 157 pks-positive and 157 pks-negative K. pneumoniae infected patients were investigated.

Results

Of 2,709 K. pneumoniae genomes, 245 pks-positive genomes were screened. The four siderophores, type VI secretion system, and nutritional factor genes were present in at least 77.9% (191/245), 66.9% (164/245), and 63.3% (155/245) of pks-positive strains, respectively. The number and fragment length of prophage were lower in pks-positive strains than in pks-negative strains (p < 0.05). The prevalence of the IS6 family was higher in pks-negative strains than in pks-positive strains, and the prevalence of multiple plasmid replicon types differed between the pks-positive and pks-negative strains (p < 0.05). The detection rate of pks-positive K. pneumoniae in abscess samples was higher than that of pks-negative K. pneumoniae (p < 0.05).

Conclusion

The pks-positive strains had abundant virulence genes. There were differences in the distribution of mobile genetic elements between pks-positive and pks-negative isolates. Further analysis of the evolutionary pattern of pks island and epidemiological surveillance in different populations are needed.

BLASTphylo: An Interactive Web Tool for Taxonomic and Phylogenetic Analysis of Prokaryotic Genes

Identifying a protein's function is crucial to reveal its role in the cellular complex. Computationally, the most common approach is to search for homologous proteins in a large database of proteins of known function using BLAST. One goal of such an analysis is the identification and visualization of the protein in the taxonomy of interest. Another goal is the reconstruction of the phylogenetic history of the protein. However, the BLAST result provides information about the occurrence of the protein in the taxonomy and its putative function mainly in a tabular format. This requires manual interventions and makes the taxonomic identification laborious. Although various tools exist to visualize and annotate large-scale trees, none of them intuitively and interactively visualizes the protein's occurrence in the taxonomy for different taxonomic ranks. To target this gap, we developed BLASTphylo, a web tool that combines BLAST with automatic taxonomic mapping and phylogenetic analysis and provides the results in interactive visualizations. We demonstrate the functionalities of BLASTphylo in two case studies.

Horizontally Transferred DNA in the Genome of the Fungus Pyricularia oryzae is Associated With Repressive Histone Modifications

Horizontal gene transfer (HGT) is a means of exchanging genetic material asexually. The process by which horizontally transferred genes are domesticated by the host genome is of great interest but is not well understood. In this study, we determined the telomere-to-telomere genome sequence of the wheat-infecting Pyricularia oryzae strain Br48. SNP analysis indicated that the Br48 strain is a hybrid of wheat- and Brachiaria-infecting strains by a sexual or parasexual cross. Comparative genomic analysis identified several megabase-scale "insertions" in the Br48 genome, some of which were possibly gained by HGT-related events from related species, such as P. pennisetigena or P. grisea. Notably, the mega-insertions often contained genes whose phylogeny is not congruent with the species phylogeny. Moreover, some of the genes have a close homolog even in distantly related organisms, such as basidiomycetes or prokaryotes, implying the involvement of multiple HGT events. Interestingly, the levels of the silent epigenetic marks H3K9me3 and H3K27me3 in a genomic region tended to be negatively correlated with the phylogenetic concordance of genes in the same region, suggesting that horizontally transferred DNA is preferentially targeted for epigenetic silencing. Indeed, the putative HGT-derived genes were activated when MoKmt6, the gene responsible for H3K27me3 modification, was deleted. Notably, these genes also tended to be up-regulated during infection, suggesting that they are now under host control and have contributed to establishing a fungal niche. In conclusion, this study suggests that epigenetic modifications have played an important role in the domestication of HGT-derived genes in the P. oryzae genome.

Insertion sequence transposition inactivates CRISPR-Cas immunity

CRISPR-Cas immunity systems safeguard prokaryotic genomes by inhibiting the invasion of mobile genetic elements. Here, we screened prokaryotic genomic sequences and identified multiple natural transpositions of insertion sequences (ISs) into cas genes, thus inactivating CRISPR-Cas defenses. We then generated an IS-trapping system, using Escherichia coli strains with various ISs and an inducible cas nuclease, to monitor IS insertions into cas genes following the induction of double-strand DNA breakage as a physiological host stress. We identified multiple events mediated by different ISs, especially IS1 and IS10, displaying substantial relaxed target specificity. IS transposition into cas was maintained in the presence of DNA repair machinery, and transposition into other host defense systems was also detected. Our findings highlight the potential of ISs to counter CRISPR activity, thus increasing bacterial susceptibility to foreign DNA invasion.

A fast and scalable method for inferring phylogenetic networks from trees by aligning lineage taxon strings

The reconstruction of phylogenetic networks is an important but challenging problem in phylogenetics and genome evolution, as the space of phylogenetic networks is vast and cannot be sampled well. One approach to the problem is to solve the minimum phylogenetic network problem, in which phylogenetic trees are first inferred, and then the smallest phylogenetic network that displays all the trees is computed. The approach takes advantage of the fact that the theory of phylogenetic trees is mature, and there are excellent tools available for inferring phylogenetic trees from a large number of biomolecular sequences. A tree-child network is a phylogenetic network satisfying the condition that every nonleaf node has at least one child that is of indegree one. Here, we develop a new method that infers the minimum tree-child network by aligning lineage taxon strings in the phylogenetic trees. This algorithmic innovation enables us to get around the limitations of the existing programs for phylogenetic network inference. Our new program, named ALTS, is fast enough to infer a tree-child network with a large number of reticulations for a set of up to 50 phylogenetic trees with 50 taxa that have only trivial common clusters in about a quarter of an hour on average.

Microbial applications for sustainable space exploration beyond low Earth orbit

With the construction of the International Space Station, humans have been continuously living and working in space for 22 years. Microbial studies in space and other extreme environments on Earth have shown the ability for bacteria and fungi to adapt and change compared to "normal" conditions. Some of these changes, like biofilm formation, can impact astronaut health and spacecraft integrity in a negative way, while others, such as a propensity for plastic degradation, can promote self-sufficiency and sustainability in space. With the next era of space exploration upon us, which will see crewed missions to the Moon and Mars in the next 10 years, incorporating microbiology research into planning, decision-making, and mission design will be paramount to ensuring success of these long-duration missions. These can include astronaut microbiome studies to protect against infections, immune system dysfunction and bone deterioration, or biological in situ resource utilization (bISRU) studies that incorporate microbes to act as radiation shields, create electricity and establish robust plant habitats for fresh food and recycling of waste. In this review, information will be presented on the beneficial use of microbes in bioregenerative life support systems, their applicability to bISRU, and their capability to be genetically engineered for biotechnological space applications. In addition, we discuss the negative effect microbes and microbial communities may have on long-duration space travel and provide mitigation strategies to reduce their impact. Utilizing the benefits of microbes, while understanding their limitations, will help us explore deeper into space and develop sustainable human habitats on the Moon, Mars and beyond.

A host of armor: Prokaryotic immune strategies against mobile genetic elements

Prokaryotic adaptation is strongly influenced by the horizontal acquisition of beneficial traits via mobile genetic elements (MGEs), such as viruses/bacteriophages and plasmids. However, MGEs can also impose a fitness cost due to their often parasitic nature and differing evolutionary trajectories. In response, prokaryotes have evolved diverse immune mechanisms against MGEs. Recently, our understanding of the abundance and diversity of prokaryotic immune systems has greatly expanded. These defense systems can degrade the invading genetic material, inhibit genome replication, or trigger abortive infection, leading to population protection. In this review, we highlight these strategies, focusing on the most recent discoveries. The study of prokaryotic defenses not only sheds light on microbial evolution but also uncovers novel enzymatic activities with promising biotechnological applications.

Distribution and potential ecological risks of microplastics in Zhushan Bay, China

The interaction between microplastics (MPs) and microorganisms may alter the distribution of antibiotic resistance genes (ARGs) in water and increase the ecological risk of drinking water sources. To investigate the characteristics of MPs geographical distribution and its potential ecological risk in typical urban water, this study was conducted in Zhushan Bay, and we carried out a combination of tests to analyze the distribution of MPs and the migration changes of their surface microbial community composition and ARGs in different media by 16S rRNA gene high-throughput sequencing, non-targeted metabolomics and qPCR genomics in the near-shore (I), middle area (Ⅱ) and near-lake (Ⅲ) of Zhushan Bay. The results showed that MPs in fibrous form were dominant in the aquatic environment of Zhushan Bay; Polyurethane (PU) and Silicone were the main MPs types in Zhushan Bay. The abundance of MPs in the water of Zhushan Bay was winter > summer > autumn > spring; and in the sediment was winter > summer > autumn > spring, respectively. The distribution results of MPs in geographical location are as follows: In the water I > Ⅱ > Ⅲ, sediment exhibited Ⅱ > Ⅲ > I. The results indicate that physicochemical factors will affect the geographical distribution of MPs and their surface microbial community composition in the aquatic environment of Zhushan Bay. More cooperative behaviors and increased metabolically important pathways occurred in the microbial network on water-MPs compared to sediment-MPs. However, the microbial community in the sediment-MPs was more stable and had higher abundance of mobile genetic elements (MGEs). A total of 362 differential metabolites were detected, of which 193 were up-regulated and 19 down-regulated differential metabolites. blaTEM, Sul, and inti1 were prevalent in both the water and sediments of Zhushan Bay. Sul1 was most contaminated in ARGs. This study provides the latest field data and insights into MPs pollution in key aquatic environments.

Effects of cadmium and copper mixtures on antibiotic resistance genes in rhizosphere soil

The evolvement and development of antibiotic resistance in microorganisms may be influenced by metals; however, it is still unclear how cadmium (Cd) and copper (Cu) combined affect the distribution and presence of antibiotic-resistance genes (ARGs) in rhizosphere soil. The aims of this research were to (1) compare the distribution patterns of bacterial communities and ARGs in response to the effects of Cd and Cu both separately and combined; (2) explore the possible mechanisms underlying the variation in soil bacterial communities and ARGs in addition to the combined effects of Cd, Cu, and various environmental variables (nutrients, pH, etc.); and (3) provide a reference for assessing the risks of metals (Cd and Cu) and ARGs. The findings showed that the multidrug resistance genes acrA and acrB and the transposon gene intI-1 were present in high relative abundance in bacterial communities. Cadmium and Cu had a substantial interaction effect on the abundance of acrA, whereas Cu had a notable main effect on the abundance of intI-1. According to the network analysis, the strong links between bacterial taxa and specific ARGs revealed that most ARGs were hosted by Proteobacteria, Actinobacteria, and Bacteroidetes. According to structural equation modeling, Cd had a larger effect on ARGs than Cu. Compared to previous analyses of ARGs, bacterial community diversity had little effect on ARGs in this study. Overall, the results may have important consequences for determining the possible hazard of soil metals and extend the understanding of how Cd and Cu co-select ARGs in rhizosphere soils.

Compensatory relationship between low-complexity regions and gene paralogy in the evolution of prokaryotes

The evolution of genomes in all life forms involves two distinct, dynamic types of genomic changes: gene duplication (and loss) that shape families of paralogous genes and extension (and contraction) of low-complexity regions (LCR), which occurs through dynamics of short repeats in protein-coding genes. Although the roles of each of these types of events in genome evolution have been studied, their co-evolutionary dynamics is not thoroughly understood. Here, by analyzing a wide range of genomes from diverse bacteria and archaea, we show that LCR and paralogy represent two distinct routes of evolution that are inversely correlated. The emergence of LCR is a prominent evolutionary mechanism in fast evolving, young protein families, whereas paralogy dominates the comparatively slow evolution of old protein families. The analysis of multiple prokaryotic genomes shows that the formation of LCR is likely a widespread, transient evolutionary mechanism that temporally and locally affects also ancestral functions, but apparently, fades away with time, under mutational and selective pressures, yielding to gene paralogy. We propose that compensatory relationships between short-term and longer-term evolutionary mechanisms are universal in the evolution of life.

Renewing Linnaean taxonomy: a proposal to restructure the highest levels of the Natural System

During the last century enormous progress has been made in the understanding of biological diversity, involving a dramatic shift from macroscopic to microscopic organisms. The question now arises as to whether the Natural System introduced by Carl Linnaeus, which has served as the central system for organizing biological diversity, can accommodate the great expansion of diversity that has been discovered. Important discoveries regarding biological diversity have not been fully integrated into a formal, coherent taxonomic system. In addition, because of taxonomic challenges and conflicts, various proposals have been made to abandon key aspects of the Linnaean system. We review the current status of taxonomy of the living world, focussing on groups at the taxonomic level of phylum and above. We summarize the main arguments against and in favour of abandoning aspects of the Linnaean system. Based on these considerations, we conclude that retaining the Linnaean Natural System provides important advantages. We propose a relatively small number of amendments for extending this system, particularly to include the named rank of world (Latin alternative mundis) formally to include non-cellular entities (viruses), and the named rank of empire (Latin alternative imperium) to accommodate the depth of diversity in (unicellular) eukaryotes that has been uncovered. We argue that in the case of both the eukaryotic domain and the viruses the cladistic approach intrinsically fails. However, the resulting semi-cladistic system provides a productive way forward that can help resolve taxonomic challenges. The amendments proposed allow us to: (i) retain named taxonomic levels and the three-domain system, (ii) improve understanding of the main eukaryotic lineages, and (iii) incorporate viruses into the Natural System. Of note, the proposal described herein is intended to serve as the starting point for a broad scientific discussion regarding the modernization of the Linnaean system.

Screening for Primordial RNA–Peptide Interactions Using High-Density Peptide Arrays

RNA–peptide interactions are an important factor in the origin of the modern mechanism of translation and the genetic code. Despite great progress in the bioinformatics of RNA–peptide interactions due to the rapid growth in the number of known RNA–protein complexes, there is no comprehensive experimental method to take into account the influence of individual amino acids on non-covalent RNA–peptide bonds. First, we designed the combinatorial libraries of primordial peptides according to the combinatorial fusion rules based on Watson–Crick mutations. Next, we used high-density peptide arrays to investigate the interaction of primordial peptides with their cognate homo-oligonucleotides. We calculated the interaction scores of individual peptide fragments and evaluated the influence of the peptide length and its composition on the strength of RNA binding. The analysis shows that the amino acids phenylalanine, tyrosine, and proline contribute significantly to the strong binding between peptides and homo-oligonucleotides, while the sum charge of the peptide does not have a significant effect. We discuss the physicochemical implications of the combinatorial fusion cascade, a hypothesis that follows from the amino acid partition used in the work.

How has horizontal gene transfer shaped the evolution of insect genomes?

As the most diverse group of animals on Earth, insects are key organisms in ecosystems. Horizontal gene transfer (HGT) refers to the transfer of genetic material between species by non-reproductive means. HGT is a major evolutionary force in prokaryotic genome evolution, but its importance in different eukaryotic groups, such as insects, has only recently begun to be understood. Genomic data from hundreds of insect species have enabled the detection of large numbers of HGT events and the elucidation of the functions of some of these foreign genes. Although quantification of the extent of HGT in insects broadens our understanding of its role in insect evolution, the scope of its influence and underlying mechanism(s) of its occurrence remain open questions for the field.

Exploring the expressiveness of abstract metabolic networks

Metabolism is characterised by chemical reactions linked to each other, creating a complex network structure. The whole metabolic network is divided into pathways of chemical reactions, such that every pathway is a metabolic function. A simplified representation of metabolism, which we call an abstract metabolic network, is a graph in which metabolic pathways are nodes and there is an edge between two nodes if their corresponding pathways share one or more compounds. The abstract metabolic network of a given organism results in a small network that requires low computational power to be analysed and makes it a suitable model to perform a large-scale comparison of organisms' metabolism. To explore the potentials and limits of such a basic representation, we considered a comprehensive set of KEGG organisms, represented through their abstract metabolic network. We performed pairwise comparisons using graph kernel methods and analyse the results through exploratory data analysis and machine learning techniques. The results show that abstract metabolic networks discriminate macro evolutionary events, indicating that they are expressive enough to capture key steps in metabolism evolution.

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.

Structural and evolutionary insights into astacin metallopeptidases

The astacins are a family of metallopeptidases (MPs) that has been extensively described from animals. They are multidomain extracellular proteins, which have a conserved core architecture encompassing a signal peptide for secretion, a prodomain or prosegment and a zinc-dependent catalytic domain (CD). This constellation is found in the archetypal name-giving digestive enzyme astacin from the European crayfish Astacus astacus. Astacin catalytic domains span ∼200 residues and consist of two subdomains that flank an extended active-site cleft. They share several structural elements including a long zinc-binding consensus sequence (HEXXHXXGXXH) immediately followed by an EXXRXDRD motif, which features a family-specific glutamate. In addition, a downstream SIMHY-motif encompasses a "Met-turn" methionine and a zinc-binding tyrosine. The overall architecture and some structural features of astacin catalytic domains match those of other more distantly related MPs, which together constitute the metzincin clan of metallopeptidases. We further analysed the structures of PRO-, MAM, TRAF, CUB and EGF-like domains, and described their essential molecular determinants. In addition, we investigated the distribution of astacins across kingdoms and their phylogenetic origin. Through extensive sequence searches we found astacin CDs in > 25,000 sequences down the tree of life from humans beyond Metazoa, including Choanoflagellata, Filasterea and Ichtyosporea. We also found < 400 sequences scattered across non-holozoan eukaryotes including some fungi and one virus, as well as in selected taxa of archaea and bacteria that are pathogens or colonizers of animal hosts, but not in plants. Overall, we propose that astacins originate in the root of Holozoa consistent with Darwinian descent and that the latter genes might be the result of horizontal gene transfer from holozoan donors.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Evolutionary Overview of Aquaporin Superfamily

Aquaporins (AQPs) are present not only in three domains of life, bacteria, eukaryotes, and archaea, but also in viruses. With the accumulating arrays of AQP superfamily, the evolutional relationship has attracted much attention with multiple publications on "the genome-wide identification and phylogenetic analysis" of AQP superfamily. A pair of NPA boxes forming a pore is highly conserved throughout the evolution and renders key residues for the classification of AQP superfamily into four groups: AQP1-like, AQP3-like, AQP8-like, and AQP11-like. The complexity of AQP family has mostly been achieved in nematodes and subsequent evolution has been directed toward increasing the number of AQPs through whole-genome duplications (WGDs) to extend the tissue specific expression and regulation. The discovery of the intracellular AQP (iAQP: AQP8-like and AQP11-like) and substrate transports by the plasma membrane AQP (pAQP: AQP1-like and AQP3-like) have accelerated the AQP research much more toward the transport of substrates with complex profiles. This evolutionary overview based on a simple classification of AQPs into four subfamilies will provide putative structural, functional, and localization information and insights into the role of AQP as well as clues to understand the complex diversity of AQP superfamily.

Mobilome of the Rhus Gall Aphid Schlechtendalia chinensis Provides Insight into TE Insertion-Related Inactivation of Functional Genes

Transposable elements (TEs) comprise a considerable proportion of insect genomic DNA; how they contribute to genome structure and organization is still poorly understood. Here, we present an analysis of the TE repertoire in the chromosome-level genome assembly of Rhus gall aphid Schlechtendalia chinensis. The TE fractions are composed of at least 32 different superfamilies and many TEs from different families were transcriptionally active in the S. chinensis genome. Furthermore, different types of transposase-derived proteins were also found in the S. chinensis genome. We also provide insight into the TEs related insertional inactivation, and exogenization of TEs in functional genes. We considered that the presence of TE fragments in the introns of functional genes could impact the activity of functional genes, and a large number of TE fragments in introns could lead to the indirect inactivation of functional genes. The present study will be beneficial in understanding the role and impact of TEs in genomic evolution of their hosts.

Gene transferability and sociality do not correlate with gene connectivity

The connectivity of a gene, defined as the number of interactions a gene's product has with other genes' products, is a key characteristic of a gene. In prokaryotes, the complexity hypothesis predicts that genes which undergo more frequent horizontal transfer will be less connected than genes which are only very rarely transferred. We tested the role of horizontal gene transfer, and other potentially important factors, by examining the connectivity of chromosomal and plasmid genes, across 134 diverse prokaryotic species. We found that (i) genes on plasmids were less connected than genes on chromosomes; (ii) connectivity of plasmid genes was not correlated with plasmid mobility; and (iii) the sociality of genes (cooperative or private) was not correlated with gene connectivity.

The Role of the Environment in Horizontal Gene Transfer

Gene-by-environment interactions play a crucial role in horizontal gene transfer by affecting how the transferred genes alter host fitness. However, how the environment modulates the fitness effect of transferred genes has not been tested systematically in an experimental study. We adapted a high-throughput technique for obtaining very precise estimates of bacterial fitness, in order to measure the fitness effects of 44 orthologs transferred from Salmonella Typhimurium to Escherichia coli in six physiologically relevant environments. We found that the fitness effects of individual genes were highly dependent on the environment, while the distributions of fitness effects across genes were not, with all tested environments resulting in distributions of same shape and spread. Furthermore, the extent to which the fitness effects of a gene varied between environments depended on the average fitness effect of that gene across all environments, with nearly neutral and nearly lethal genes having more consistent fitness effects across all environments compared to deleterious genes. Put together, our results reveal the unpredictable nature of how environmental conditions impact the fitness effects of each individual gene. At the same time, distributions of fitness effects across environments exhibit consistent features, pointing to the generalizability of factors that shape horizontal gene transfer of orthologous genes.

Zymogenic latency in an ∼250-million-year-old astacin metallopeptidase

The horseshoe crab Limulus polyphemus is one of few extant Limulus species, which date back to ∼250 million years ago under the conservation of a common Bauplan documented by fossil records. It possesses the only proteolytic blood-coagulation and innate immunity system outside vertebrates and is a model organism for the study of the evolution and function of peptidases. The astacins are a family of metallopeptidases that share a central ∼200-residue catalytic domain (CD), which is found in >1000 species across holozoans and, sporadically, bacteria. Here, the zymogen of an astacin from L. polyphemus was crystallized and its structure was solved. A 34-residue, mostly unstructured pro-peptide (PP) traverses, and thus blocks, the active-site cleft of the CD in the opposite direction to a substrate. A central `PP motif' (F35-E-G-D-I39) adopts a loop structure which positions Asp38 to bind the catalytic metal, replacing the solvent molecule required for catalysis in the mature enzyme according to an `aspartate-switch' mechanism. Maturation cleavage of the PP liberates the cleft and causes the rearrangement of an `activation segment'. Moreover, the mature N-terminus is repositioned to penetrate the CD moiety and is anchored to a buried `family-specific' glutamate. Overall, this mechanism of latency is reminiscent of that of the other three astacins with known zymogenic and mature structures, namely crayfish astacin, human meprin β and bacterial myroilysin, but each shows specific structural characteristics. Remarkably, myroilysin lacks the PP motif and employs a cysteine instead of the aspartate to block the catalytic metal.

Modulating Gene Expression within a Microbiome Based on Computational Models

Simple Summary

Development of computational biology methodologies has provided comprehensive understanding of the complexity of microbiomes, and the extensive ways in which they influence their environment. This has awakened a new research goal, aiming to not only understand the mechanisms in which microbiomes function, but to actively modulate and engineer them for various purposes. However, current microbiome engineering techniques are usually manually tailored for a specific system and neglect the different interactions between the new genetic information and the bacterial population, turning a blind eye to processes such as horizontal gene transfer, mutations, and other genetic alterations. In this work, we developed a generic computational method to automatically tune the expression of heterologous genes within a microbiome according to given preferences, to allow the functionality of the engineering process to propagate in longer periods of time. This goal was achieved by treating each part of the gene individually and considering long term fitness effects on the environment, providing computational and experimental evidence for this approach.

Abstract

Recent research in the field of bioinformatics and molecular biology has revealed the immense complexity and uniqueness of microbiomes, while also showcasing the impact of the symbiosis between a microbiome and its host or environment. A core property influencing this process is horizontal gene transfer between members of the bacterial community used to maintain genetic variation. The essential effect of this mechanism is the exposure of genetic information to a wide array of members of the community, creating an additional “layer” of information in the microbiome named the “plasmidome”. From an engineering perspective, introduction of genetic information to an environment must be facilitated into chosen species which will be able to carry out the desired effect instead of competing and inhibiting it. Moreover, this process of information transfer imposes concerns for the biosafety of genetic engineering of microbiomes as exposure of genetic information into unwanted hosts can have unprecedented ecological impacts. Current technologies are usually experimentally developed for a specific host/environment, and only deal with the transformation process itself at best, ignoring the impact of horizontal gene transfer and gene-microbiome interactions that occur over larger periods of time in uncontrolled environments. The goal of this research was to design new microbiome-specific versions of engineered genetic information, providing an additional layer of compatibility to existing engineering techniques. The engineering framework is entirely computational and is agnostic to the selected microbiome or gene by reducing the problem into the following set up: microbiome species can be defined as wanted or unwanted hosts of the modification. Then, every element related to gene expression (e.g., promoters, coding regions, etc.) and regulation is individually examined and engineered by novel algorithms to provide the defined expression preferences. Additionally, the synergistic effect of the combination of engineered gene blocks facilitates robustness to random mutations that might occur over time. This method has been validated using both computational and experimental tools, stemming from the research done in the iGEM 2021 competition, by the TAU group.

Distribution, Polymorphism and Function Characteristics of the GST-Encoding Fhb7 in Triticeae

Encoding a glutathione S-transferase (GST) and conferring resistance to Fusarium head blight (FHB), Fhb7 was successfully isolated from the newly assembled Thinopyrum elongatum genome by researchers, with blasting searches revealing that Thinopyrum gained Fhb7 through horizontal gene transfer from an endophytic Epichloë species. On the contrary, our molecular evidence reveals that the homologs of Fhb7 are distributed commonly in Triticeae. Other than Thinopyrum, the Fhb7 homologs were also detected in four other genera, Elymus, Leymus, Roegneria and Pseudoroegneria, respectively. Sequence comparisons revealed that the protein sequences were at least 94% identical across all of the Fhb7 homologs in Triticeae plants, which in turn suggested that the horizontal gene transfer of the Fhb7 might have occurred before Triticeae differentiation instead of Thinopyrum. The multiple Fhb7 homologs detected in some Triticeae accessions and wheat-Thinopyrum derivatives might be attributed to the alloploid nature and gene duplication during evolution. In addition, we discovered that some wheat-Thinopyrum derivatives carrying the Fhb7 homologs had a completely different reaction to Fusarium head blight, which made us question the ability of the GST-encoding Fhb7 to resist FHB.

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.

Revisiting the Tigger Transposon Evolution Revealing Extensive Involvement in the Shaping of Mammal Genomes

Simple Summary

Despite the discovery of the Tigger family of pogo transposons in the mammalian genome, the evolution profile of this family is still incomplete. Here, we conducted a systematic evolution analysis for Tigger in nature. The data revealed that Tigger was found in a broad variety of animals, and extensive invasion of Tigger was observed in mammal genomes. Common horizontal transfer events of Tigger elements were observed across different lineages of animals, including mammals, that may have led to their widespread distribution, while parasites and invasive species may have promoted Tigger HT events. Our results also indicate that the activity of Tigger transposons tends to be low in vertebrates; only one mammalian genome and fish genome may harbor active Tigger.

Abstract

The data of this study revealed that Tigger was found in a wide variety of animal genomes, including 180 species from 36 orders of invertebrates and 145 species from 29 orders of vertebrates. An extensive invasion of Tigger was observed in mammals, with a high copy number. Almost 61% of those species contain more than 50 copies of Tigger; however, 46% harbor intact Tigger elements, although the number of these intact elements is very low. Common HT events of Tigger elements were discovered across different lineages of animals, including mammals, that may have led to their widespread distribution, whereas Helogale parvula and arthropods may have aided Tigger HT incidences. The activity of Tigger seems to be low in the kingdom of animals, most copies were truncated in the mammal genomes and lost their transposition activity, and Tigger transposons only display signs of recent and current activities in a few species of animals. The findings suggest that the Tigger family is important in structuring mammal genomes.

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.

Enhancing the Phytoremediation of Heavy Metals by Combining Hyperaccumulator and Heavy Metal-Resistant Plant Growth-Promoting Bacteria

Heavy metals (HMs) have become a major environmental pollutant threatening ecosystems and human health. Although hyperaccumulators provide a viable alternative for the bioremediation of HMs, the potential of phytoremediation is often limited by the small biomass and slow growth rate of hyperaccumulators and HM toxicity to plants. Here, plant growth-promoting bacteria (PGPB)-assisted phytoremediation was used to enhance the phytoremediation of HM-contaminated soils. A PGPB with HM-tolerant (HMT-PGPB), Bacillus sp. PGP15 was isolated from the rhizosphere of a cadmium (Cd) hyperaccumulator, Solanum nigrum. Pot experiments demonstrated that inoculation with strain PGP15 could significantly increase the growth of S. nigrum. More importantly, strain PGP15 markedly improved Cd accumulation in S. nigrum while alleviating Cd-induced stress in S. nigrum. Specifically, PGP15 inoculation significantly decreased the contents of H₂O₂, MDA, and O 2 · - in S. nigrum, while the activities (per gram plant fresh weight) of SOD, APX, and CAT were significantly increased in the PGP15-inoculated plants compared with the control sample. These results suggested that the interactions between strain PGP15 and S. nigrum could overcome the limits of phytoremediation alone and highlighted the promising application potential of the PGPB-hyperaccumulator collaborative pattern in the bioremediation of HM-contaminated soils. Furthermore, the PGP15 genome was sequenced and compared with other strains to explore the mechanisms underlying plant growth promotion by HMT-PGPB. The results showed that core genes that define the fundamental metabolic capabilities of strain PGP15 might not be necessary for plant growth promotion. Meanwhile, PGP15-specific genes, including many transposable elements, played a crucial role in the adaptive evolution of HM resistance. Overall, our results improve the understanding of interactions between HMT-PGPB and plants and facilitate the application of HMT-PGPB in the phytoremediation of HM-contaminated soils.

Phthalate hydrolase: distribution, diversity and molecular evolution

The alpha/beta-fold superfamily of hydrolases is rapidly becoming one of the largest groups of structurally related enzymes with diverse catalytic functions. In this superfamily of enzymes, esterase deserves special attention because of their wide distribution in biological systems and importance towards environmental and industrial applications. Among various esterases, phthalate hydrolases are the key alpha/beta enzymes involved in the metabolism of structurally diverse estrogenic phthalic acid esters, ubiquitously distributed synthetic chemicals, used as plasticizer in plastic manufacturing processes. Although they vary both at the sequence and functional levels, these hydrolases use a similar acid-base-nucleophile catalytic mechanism to catalyse reactions on structurally different substrates. The current review attempts to present insights on phthalate hydrolases, describing their sources, structural diversities, phylogenetic affiliations and catalytically different types or classes of enzymes, categorized as diesterase, monoesterase and diesterase-monoesterase, capable of hydrolysing phthalate diester, phthalate monoester and both respectively. Furthermore, available information on in silico analyses and site-directed mutagenesis studies revealing structure-function integrity and altered enzyme kinetics have been highlighted along with the possible scenario of their evolution at the molecular level.

Complex and flexible catabolism in Aromatoleum aromaticum pCyN1

Large quantities of organic matter are continuously deposited, and (a)biotic gradients intersect in the soil-rhizosphere, where biodegradation contributes to the global cycles of elements. The betaproteobacterial genus Aromatoleum comprises cosmopolitan, facultative denitrifying degradation specialists. A. aromaticum pCyN1 stands out for anaerobically decomposing plant-derived monoterpenes in addition to monoaromatic hydrocarbons, polar aromatics and aliphatics. The catabolic network's structure and flexibility in A. aromaticum pCyN1 was studied across 34 growth conditions by superimposing proteome profiles onto the manually annotated 4.37 Mbp genome. Strain pCyN1 employs three fundamentally different enzymes for C-H-bond cleavage at the methyl groups of p-cymene/4-ethyltoluene, toluene and p-cresol, respectively. Regulation of degradation modules displayed substrate specificities ranging from narrow (toluene and cyclohexane carboxylate) via medium-wide (one module shared by p-cymene, 4-ethyltoluene, α-phellandrene, α-terpinene, γ-terpinene and limonene) to broad (central benzoyl-CoA pathway serving 16 aromatic substrates). Remarkably, three variants of ATP-dependent (class I) benzoyl-CoA reductase and four different β-oxidation routes establish a degradation hub that accommodates the substrate diversity. The respiratory system displayed several conspicuous profiles, e.g., the presence of nitrous oxide reductase under oxic and of low-affinity oxidase under anoxic conditions. Overall, nutritional versatility in conjunction with network regulation endow A. aromaticum pCyN1 with broad adaptability. This article is protected by copyright. All rights reserved.

The Notable Achievements and the Prospects of Bacterial Pathogen Genomics

Throughout the entirety of human history, bacterial pathogens have played an important role and even shaped the fate of civilizations. The application of genomics within the last 27 years has radically changed the way we understand the biology and evolution of these pathogens. In this review, we discuss how the short- (Illumina) and long-read (PacBio, Oxford Nanopore) sequencing technologies have shaped the discipline of bacterial pathogen genomics, in terms of fundamental research (i.e., evolution of pathogenicity), forensics, food safety, and routine clinical microbiology. We have mined and discuss some of the most prominent data/bioinformatics resources such as NCBI pathogens, PATRIC, and Pathogenwatch. Based on this mining, we present some of the most popular sequencing technologies, hybrid approaches, assemblers, and annotation pipelines. A small number of bacterial pathogens are of very high importance, and we also present the wealth of the genomic data for these species (i.e., which ones they are, the number of antimicrobial resistance genes per genome, the number of virulence factors). Finally, we discuss how this discipline will probably be transformed in the near future, especially by transitioning into metagenome-assembled genomes (MAGs), thanks to long-read sequencing.

Importance of RpoD- and Non-RpoD-Dependent Expression of Horizontally Acquired Genes in Cupriavidus metallidurans

The genome of the metal-resistant, hydrogen-oxidizing bacterium Cupriavidus metallidurans contains a large number of horizontally acquired plasmids and genomic islands that were integrated into its chromosome or chromid. For the C. metallidurans CH34 wild-type strain growing under nonchallenging conditions, 5,763 transcriptional starting sequences (TSSs) were determined. Using a custom-built motif discovery software based on hidden Markov models, patterns upstream of the TSSs were identified. The pattern TTGACA, −35.6 ± 1.6 bp upstream of the TSSs, in combination with a TATAAT sequence 15.8 ± 1.4 bp upstream occurred frequently, especially upstream of the TSSs for 48 housekeeping genes, and these were assigned to promoters used by RNA polymerase containing the main housekeeping sigma factor RpoD. From patterns upstream of the housekeeping genes, a score for RpoD-dependent promoters in C. metallidurans was derived and applied to all 5,763 TSSs. Among these, 2,572 TSSs could be associated with RpoD with high probability, 373 with low probability, and 2,818 with no probability. In a detailed analysis of horizontally acquired genes involved in metal resistance and not involved in this process, the TSSs responsible for the expression of these genes under nonchallenging conditions were assigned to RpoD- or non-RpoD-dependent promoters. RpoD-dependent promoters occurred frequently in horizontally acquired metal resistance and other determinants, which should allow their initial expression in a new host. However, other sigma factors and sense/antisense effects also contribute—maybe to mold in subsequent adaptation steps the assimilated gene into the regulatory network of the cell.

IMPORTANCE In their natural environment, bacteria are constantly acquiring genes by horizontal gene transfer. To be of any benefit, these genes should be expressed. We show here that the main housekeeping sigma factor RpoD plays an important role in the expression of horizontally acquired genes in the metal-resistant hydrogen-oxidizing bacterium C. metallidurans. By conservation of the RpoD recognition consensus sequence, a newly arriving gene has a high probability to be expressed in the new host cell. In addition to integrons and genes travelling together with that for their sigma factor, conservation of the RpoD consensus sequence may be an important contributor to the overall evolutionary success of horizontal gene transfer in bacteria. Using C. metallidurans as an example, this publication sheds some light on the fate and function of horizontally acquired genes in bacteria.

Comparative Genomics Reveals Genetic Diversity and Metabolic Potentials of the Genus Qipengyuania and Suggests Fifteen Novel Species

Members of the genus Qipengyuania are heterotrophic bacteria frequently isolated from marine environments with great application potential in areas such as carotenoid production. However, the genomic diversity, metabolic function, and adaption of this genus remain largely unclear. Here, 16 isolates related to the genus Qipengyuania were recovered from coastal samples and their genomes were sequenced. The phylogenetic inference of these isolates and reference type strains of this genus indicated that the 16S rRNA gene was insufficient to distinguish them at the species level; instead, the phylogenomic reconstruction could provide the reliable phylogenetic relationships and confirm 15 new well-supported branches, representing 15 putative novel genospecies corroborated by the digital DNA-DNA hybridization and average nucleotide identity analyses. Comparative genomics revealed that the genus Qipengyuania had an open pangenome and possessed multiple conserved genes and pathways related to metabolic functions and environmental adaptation, despite the presence of divergent genomic features and specific metabolic potential. Genetic analysis and pigment detection showed that the members of this genus were identified as carotenoid producers, while some proved to be potentially aerobic anoxygenic photoheterotrophs. Collectively, the first insight into the genetic diversity and metabolic potentials of the genus Qipengyuania will contribute to better understanding of the speciation and adaptive evolution in natural environments.

IMPORTANCE The deciphering of the phylogenetic diversity and metabolic features of the abundant bacterial taxa is critical for exploring their ecological importance and application potential. Qipengyuania is a genus of frequently isolated heterotrophic microorganisms with great industrial application potential. Numerous strains related to the genus Qipengyuania have been isolated from diverse environments, but their genomic diversity and metabolic functions remain unclear. Our study revealed a high degree of genetic diversity, metabolic versatility, and environmental adaptation of the genus Qipengyuania using comparative genomics. Fifteen novel species of this genus have been established using a polyphasic taxonomic approach, expanding the number of described species to almost double. This study provided an overall view of the genus Qipengyuania at the genomic level and will enable us to better uncover its ecological roles and evolutionary history.

Recruitment of Mobile Genetic Elements for Diverse Cellular Functions in Prokaryotes

Prokaryotic genomes are replete with mobile genetic elements (MGE) that span a continuum of replication autonomy. On numerous occasions during microbial evolution, diverse MGE lose their autonomy altogether but, rather than being quickly purged from the host genome, assume a new function that benefits the host, rendering the immobilized MGE subject to purifying selection, and resulting in its vertical inheritance. This mini-review highlights the diversity of the repurposed (exapted) MGE as well as the plethora of cellular functions that they perform. The principal contribution of the exaptation of MGE and their components is to the prokaryotic functional systems involved in biological conflicts, and in particular, defense against viruses and other MGE. This evolutionary entanglement between MGE and defense systems appears to stem both from mechanistic similarities and from similar evolutionary predicaments whereby both MGEs and defense systems tend to incur fitness costs to the hosts and thereby evolve mechanisms for survival including horizontal mobility, causing host addiction, and exaptation for functions beneficial to the host. The examples discussed demonstrate that the identity of an MGE, overall mobility and relationship with the host cell (mutualistic, symbiotic, commensal, or parasitic) are all factors that affect exaptation.