Sequence analysis of tyrosine recombinases allows annotation of mobile genetic elements in prokaryotic genomes

Genomic analysis of hyperparasitic viruses associated with entomopoxviruses

Polinton-like viruses (PLVs) are a diverse group of small integrative dsDNA viruses that infect diverse eukaryotic hosts. Many PLVs are hypothesized to parasitize viruses in the phylum Nucleocytoviricota for their own propagation and spread. Here, we analyze the genomes of novel PLVs associated with the occlusion bodies of entomopoxvirus (EPV) infections of two separate lepidopteran hosts. The presence of these elements within EPV occlusion bodies suggests that they are the first known hyperparasites of poxviruses. We find that these PLVs belong to two distinct lineages that are highly diverged from known PLVs. These PLVs possess mosaic genomes, and some essential genes share homology with mobile genes within EPVs. Based on this homology and observed PLV mosaicism, we propose a mechanism to explain the turnover of PLV replication and integration genes.

Systematic identification of cargo-mobilizing genetic elements reveals new dimensions of eukaryotic diversity

Cargo-mobilizing mobile elements (CMEs) are genetic entities that faithfully transpose diverse protein coding sequences. Although common in bacteria, we know little about eukaryotic CMEs because no appropriate tools exist for their annotation. For example, Starships are giant fungal CMEs whose functions are largely unknown because they require time-intensive manual curation. To address this knowledge gap, we developed starfish, a computational workflow for high-throughput eukaryotic CME annotation. We applied starfish to 2 899 genomes of 1 649 fungal species and found that starfish recovers known Starships with 95% combined precision and recall while expanding the number of annotated elements ten-fold. Extant Starship diversity is partitioned into 11 families that differ in their enrichment patterns across fungal classes. Starship cargo changes rapidly such that elements from the same family differ substantially in their functional repertoires, which are predicted to contribute to diverse biological processes such as metabolism. Many elements have convergently evolved to insert into 5S rDNA and AT-rich sequence while others integrate into random locations, revealing both specialist and generalist strategies for persistence. Our work establishes a framework for advancing mobile element biology and provides the means to investigate an emerging dimension of eukaryotic genetic diversity, that of genomes within genomes.

Redefining the bacteriophage mv4 site-specific recombination system and the sequence specificity of its attB and core-attP sites

Through their involvement in the integration and excision of a large number of mobile genetic elements, such as phages and integrative and conjugative elements (ICEs), site-specific recombination systems based on heterobivalent tyrosine recombinases play a major role in genome dynamics and evolution. However, despite hundreds of these systems having been identified in genome databases, very few have been described in detail, with none from phages that infect Bacillota (formerly Firmicutes). In this study, we reanalyzed the recombination module of Lactobacillus delbrueckii subsp. bulgaricus phage mv4, previously considered atypical compared with classical systems. Our results reveal that mv4 integrase is a 369 aa protein with all the structural hallmarks of recombinases from the Tn916 family and that it cooperatively interacts with its recombination sites. Using randomized DNA libraries, NGS sequencing, and other molecular approaches, we show that the 21-bp core-attP and attB sites have structural similarities to classical systems only if considering the nucleotide degeneracy, with two 7-bp inverted regions corresponding to mv4Int core-binding sites surrounding a 7-bp strand-exchange region. We also examined the different compositional constraints in the core-binding regions, which define the sequence space of permissible recombination sites.

The evolution of antibiotic resistance islands occurs within the framework of plasmid lineages

Bacterial pathogens carrying multidrug resistance (MDR) plasmids are a major threat to human health. The acquisition of antibiotic resistance genes (ARGs) in plasmids is often facilitated by mobile genetic elements that copy or translocate ARGs between DNA molecules. The agglomeration of mobile elements in plasmids generates resistance islands comprising multiple ARGs. However, whether the emergence of resistance islands is restricted to specific MDR plasmid lineages remains understudied. Here we show that the agglomeration of ARGs in resistance islands is biased towards specific large plasmid lineages. Analyzing 6784 plasmids in 2441 Escherichia, Salmonella, and Klebsiella isolates, we quantify that 84% of the ARGs in MDR plasmids are found in resistance islands. We furthermore observe rapid evolution of ARG combinations in resistance islands. Most regions identified as resistance islands are shared among closely related plasmids but rarely among distantly related plasmids. Our results suggest the presence of barriers for the dissemination of ARGs between plasmid lineages, which are related to plasmid genetic properties, host range and the plasmid evolutionary history. The agglomeration of ARGs in plasmids is attributed to the workings of mobile genetic elements that operate within the framework of existing plasmid lineages.

An integrated in-silico approach for drug target identification in human pathogen Shigella dysenteriae

Shigella dysenteriae, is a Gram-negative bacterium that emerged as the second most significant cause of bacillary dysentery. Antibiotic treatment is vital in lowering Shigella infection rates, yet the growing global resistance to broad-spectrum antibiotics poses a significant challenge. The persistent multidrug resistance of S. dysenteriae complicates its management and control. Hence, there is an urgent requirement to discover novel therapeutic targets and potent medications to prevent and treat this disease. Therefore, the integration of bioinformatics methods such as subtractive and comparative analysis provides a pathway to compute the pan-genome of S. dysenteriae. In our study, we analysed a dataset comprising 27 whole genomes. The S. dysenteriae strain SD197 was used as the reference for determining the core genome. Initially, our focus was directed towards the identification of the proteome of the core genome. Moreover, several filters were applied to the core genome, including assessments for non-host homology, protein essentiality, and virulence, in order to prioritize potential drug targets. Among these targets were Integration host factor subunit alpha and Tyrosine recombinase XerC. Furthermore, four drug-like compounds showing potential inhibitory effects against both target proteins were identified. Subsequently, molecular docking analysis was conducted involving these targets and the compounds. This initial study provides the list of novel targets against S. dysenteriae. Conclusively, future in vitro investigations could validate our in-silico findings and uncover potential therapeutic drugs for combating bacillary dysentery infection.

Bacteriophage lambda site-specific recombination

The site-specific recombination pathway of bacteriophage λ encompasses isoenergetic but highly directional and tightly regulated integrative and excisive reactions that integrate and excise the vial chromosome into and out of the bacterial chromosome. The reactions require 240 bp of phage DNA and 21 bp of bacterial DNA comprising 16 protein binding sites that are differentially used in each pathway by the phage-encoded Int and Xis proteins and the host-encoded integration host factor and factor for inversion stimulation proteins. Structures of higher-order protein-DNA complexes of the four-way Holliday junction recombination intermediates provided clarifying insights into the mechanisms, directionality, and regulation of these two pathways, which are tightly linked to the physiology of the bacterial host cell. Here we review our current understanding of the mechanisms responsible for regulating and executing λ site-specific recombination, with an emphasis on key studies completed over the last decade.

Geographical distribution of mobile genetic elements in microbial communities along the Yucatan coast

Horizontal gene transfer (HGT) is a well-documented strategy used by bacteria to enhance their adaptability to challenging environmental conditions. Through HGT, a group of conserved genetic elements known as mobile genetic elements (MGEs) is disseminated within bacterial communities. MGEs offer numerous advantages to the host, increasing its fitness by acquiring new functions that help bacteria contend with adverse conditions, including exposure to heavy metal and antibiotics. This study explores MGEs within microbial communities along the Yucatan coast using a metatranscriptomics approach. Prior to this research, nothing was known about the coastal Yucatan's microbial environmental mobilome and HGT processes between these bacterial communities. This study reveals a positive correlation between MGEs and antibiotic resistance genes (ARGs) along the Yucatan coast, with higher MGEs abundance in more contaminated sites. The Proteobacteria and Firmicutes groups exhibited the highest number of MGEs. It's important to highlight that the most abundant classes of MGEs might not be the ones most strongly linked to ARGs, as observed for the recombination/repair class. This work presents the first geographical distribution of the environmental mobilome in Yucatan Peninsula mangroves.

Horizontal transfer potential of antibiotic resistance genes in wastewater treatment plants unraveled by microfluidic-based mini-metagenomics

Wastewater treatment plants (WWTPs) are known to harbor antibiotic resistance genes (ARGs), which can potentially spread to the environment and human populations. However, the extent and mechanisms of ARG transfer in WWTPs are not well understood due to the high microbial diversity and limitations of molecular techniques. In this study, we used a microfluidic-based mini-metagenomics approach to investigate the transfer potential and mechanisms of ARGs in activated sludge from WWTPs. Our results show that while diverse ARGs are present in activated sludge, only a few highly similar ARGs are observed across different taxa, indicating limited transfer potential. We identified two ARGs, ermF and tla-1, which occur in a variety of bacterial taxa and may have high transfer potential facilitated by mobile genetic elements. Interestingly, genes that are highly similar to the sequences of these two ARGs, as identified in this study, display varying patterns of abundance across geographic regions. Genes similar to ermF found are widely found in Asia and the Americas, while genes resembling tla-1 are primarily detected in Asia. Genes similar to both genes are barely detected in European WWTPs. These findings shed light on the limited horizontal transfer potential of ARGs in WWTPs and highlight the importance of monitoring specific ARGs in different regions to mitigate the spread of antibiotic resistance.

Large language models improve annotation of prokaryotic viral proteins

Viral genomes are poorly annotated in metagenomic samples, representing an obstacle to understanding viral diversity and function. Current annotation approaches rely on alignment-based sequence homology methods, which are limited by the paucity of characterized viral proteins and divergence among viral sequences. Here we show that protein language models can capture prokaryotic viral protein function, enabling new portions of viral sequence space to be assigned biologically meaningful labels. When applied to global ocean virome data, our classifier expanded the annotated fraction of viral protein families by 29%. Among previously unannotated sequences, we highlight the identification of an integrase defining a mobile element in marine picocyanobacteria and a capsid protein that anchors globally widespread viral elements. Furthermore, improved high-level functional annotation provides a means to characterize similarities in genomic organization among diverse viral sequences. Protein language models thus enhance remote homology detection of viral proteins, serving as a useful complement to existing approaches.

Integrase-associated niche differentiation of endogenous large DNA viruses in crustaceans

IMPORTANCE

Crustacean genomes harbor sequences originating from a family of large DNA viruses called nimaviruses, but it is unclear why they are present. We show that endogenous nimaviruses selectively insert into repetitive sequences within the host genome, and this insertion specificity was correlated with different types of integrases, which are DNA recombination enzymes encoded by the nimaviruses themselves. This suggests that endogenous nimaviruses have colonized various genomic niches through the acquisition of integrases with different insertion specificities. Our results point to a novel survival strategy of endogenous large DNA viruses colonizing the host genomes. These findings may clarify the evolution and spread of nimaviruses in crustaceans and lead to measures to control and prevent the spread of pathogenic nimaviruses in aquaculture settings.

Complete genomes of Clostridium botulinum type B(F) isolates associated with a 1995 foodborne botulism outbreak from commercial pâté reveals a recombination event disrupting the ntnh gene

Foodborne botulism is a neuroparalytic disease caused by ingestion of foods contaminated with botulinum neurotoxin (BoNT), produced by Clostridium botulinum. In 1995 a husband and wife from Québec, Canada, were hospitalized for several months with prolonged muscle paralysis after ingesting a commercial pâté de campagne. Examination of faecal samples from both patients and the pâté produced viable Group I (proteolytic) C. botulinum type B from each of the three samples. Whole genome sequencing revealed that all three isolates contain identical bont/B5 and bont/F2 genes encoded on a plasmid. Both faecal isolate genomes were identical in chromosome and plasmid length, as well as gene content. The genome of the pâté isolate was nearly identical to that of the faecal isolates with the notable difference of a missing 13-gene insertion on the bont/B5 cluster disrupting the ntnh gene. Examination of the insertion revealed several mobile genetic elements that participate in recombination.

Annotation and Comparative Genomics of Prokaryotic Transposable Elements

The data generated in nearly 30 years of bacterial genome sequencing has revealed the abundance of transposable elements (TE) and their importance in genome and transcript remodeling through the mediation of DNA insertions and deletions, structural rearrangements, and regulation of gene expression. Furthermore, what we have learned from studying transposition mechanisms and their regulation in bacterial TE is fundamental to our current understanding of TE in other organisms because much of what has been observed in bacteria is conserved in all domains of life. However, unlike eukaryotic TE, prokaryotic TE sequester and transmit important classes of genes that impact host fitness, such as resistance to antibiotics and heavy metals and virulence factors affecting animals and plants, among other acquired traits. This provides dynamism and plasticity to bacteria, which would otherwise be propagated clonally. The insertion sequences (IS), the simplest form of prokaryotic TE, are autonomous and compact mobile genetic elements. These can be organized into compound transposons, in which two similar IS can flank any DNA segment and render it transposable. Other more complex structures, called unit transposons, can be grouped into four major families (Tn3, Tn7, Tn402, Tn554) with specific genetic characteristics. This chapter will revisit the prominent structural features of these elements, focusing on a genomic annotation framework and comparative analysis. Relevant aspects of TE will also be presented, stressing their key position in genome impact and evolution, especially in the emergence of antimicrobial resistance and other adaptive traits.

Meeting report: transposable elements at the crossroads of evolution, health and disease 2023

The conference "Transposable Elements at the Crossroads of Evolution, Health and Disease" was hosted by Keystone Symposia in Whistler, British Columbia, Canada, on September 3-6, 2023, and was organized by Kathleen Burns, Harmit Malik and Irina Arkhipova. The central theme of the meeting was the incredible diversity of ways in which transposable elements (TEs) interact with the host, from disrupting the existing genes and pathways to creating novel gene products and expression patterns, enhancing the repertoire of host functions, and ultimately driving host evolution. The meeting was organized into six plenary sessions and two afternoon workshops with a total of 50 invited and contributed talks, two poster sessions, and a career roundtable. The topics ranged from TE roles in normal and pathological processes to restricting and harnessing TE activity based on mechanistic insights gained from genetic, structural, and biochemical studies.

A systematic approach to classify and characterize genomic islands driven by conjugative mobility using protein signatures

Genomic islands (GIs) play a crucial role in the spread of antibiotic resistance, virulence factors and antiviral defense systems in a broad range of bacterial species. However, the characterization and classification of GIs are challenging due to their relatively small size and considerable genetic diversity. Predicting their intercellular mobility is of utmost importance in the context of the emerging crisis of multidrug resistance. Here, we propose a large-scale classification method to categorize GIs according to their mobility profile and, subsequently, analyze their gene cargo. We based our classification decision scheme on a collection of mobility protein motif definitions available in publicly accessible databases. Our results show that the size distribution of GI classes correlates with their respective structure and complexity. Self-transmissible GIs are usually the largest, except in Bacillota and Actinomycetota, accumulate antibiotic and phage resistance genes, and favour the use of a tyrosine recombinase to insert into a host's replicon. Non-mobilizable GIs tend to use a DDE transposase instead. Finally, although tRNA genes are more frequently targeted as insertion sites by GIs encoding a tyrosine recombinase, most GIs insert in a protein-encoding gene. This study is a stepping stone toward a better characterization of mobile GIs in bacterial genomes and their mechanism of mobility.

For Someone, You Are the Whole World: Host-Specificity of Salmonella enterica

Salmonella enterica is a bacterial pathogen known to cause gastrointestinal infections in diverse hosts, including humans and animals. Despite extensive knowledge of virulence mechanisms, understanding the factors driving host specificity remains limited. In this study, we performed a comprehensive pangenome-wide analysis of S. enterica to identify potential loci determining preference towards certain hosts. We used a dataset of high-quality genome assemblies grouped into 300 reference clusters with a special focus on four host groups: humans, pigs, cattle, and birds. The reconstructed pangenome was shown to be open and enriched with the accessory component implying high genetic diversity. Notably, phylogenetic inferences did not correspond to the distribution of affected hosts, as large compact phylogenetic groups were absent. By performing a pangenome-wide association study, we identified potential host specificity determinants. These included multiple genes encoding proteins involved in distinct infection stages, e.g., secretion systems, surface structures, transporters, transcription regulators, etc. We also identified antibiotic resistance loci in host-adapted strains. Functional annotation corroborated the results obtained with significant enrichments related to stress response, antibiotic resistance, ion transport, and surface or extracellular localization. We suggested categorizing the revealed specificity factors into three main groups: pathogenesis, resistance to antibiotics, and propagation of mobile genetic elements (MGEs).

Modularity and diversity of target selectors in Tn7 transposons

To spread, transposons must integrate into target sites without disruption of essential genes while avoiding host defense systems. Tn7-like transposons employ multiple mechanisms for target-site selection, including protein-guided targeting and, in CRISPR-associated transposons (CASTs), RNA-guided targeting. Combining phylogenomic and structural analyses, we conducted a broad survey of target selectors, revealing diverse mechanisms used by Tn7 to recognize target sites, including previously uncharacterized target-selector proteins found in newly discovered transposable elements (TEs). We experimentally characterized a CAST I-D system and a Tn6022-like transposon that uses TnsF, which contains an inactivated tyrosine recombinase domain, to target the comM gene. Additionally, we identified a non-Tn7 transposon, Tsy, encoding a homolog of TnsF with an active tyrosine recombinase domain, which we show also inserts into comM. Our findings show that Tn7 transposons employ modular architecture and co-opt target selectors from various sources to optimize target selection and drive transposon spread.

Large language models improve annotation of viral proteins

Viral sequences are poorly annotated in environmental samples, a major roadblock to understanding how viruses influence microbial community structure. Current annotation approaches rely on alignment-based sequence ho-mology methods, which are limited by available viral sequences and sequence divergence in viral proteins. Here, we show that protein language model representations capture viral protein function beyond the limits of remote sequence homology by targeting two axes of viral sequence annotation: systematic labeling of protein families and function identification for biologic discovery. Protein language model representations capture protein functional properties specific to viruses and expand the annotated fraction of ocean virome viral protein sequences by 37%. Among unannotated viral protein families, we identify a novel DNA editing protein family that defines a new mobile element in marine picocyanobacteria. Protein language models thus significantly enhance remote homology detection of viral proteins and can be utilized to enable new biological discovery across diverse functional categories.

Starships are active eukaryotic transposable elements mobilized by a new family of tyrosine recombinases

Transposable elements in eukaryotic organisms have historically been considered "selfish," at best conferring indirect benefits to their host organisms. The Starships are a recently discovered feature in fungal genomes that are, in some cases, predicted to confer beneficial traits to their hosts and also have hallmarks of being transposable elements. Here, we provide experimental evidence that Starships are indeed autonomous transposons, using the model Paecilomyces variotii, and identify the HhpA "Captain" tyrosine recombinase as essential for their mobilization into genomic sites with a specific target site consensus sequence. Furthermore, we identify multiple recent horizontal gene transfers of Starships, implying that they jump between species. Fungal genomes have mechanisms to defend against mobile elements, which are frequently detrimental to the host. We discover that Starships are also vulnerable to repeat-induced point mutation defense, thereby having implications on the evolutionary stability of such elements.

A Maverick-like cluster in the genome of a pathogenic, moderately virulent strain of Gallibacterium anatis, ESV200, a transient biofilm producer

Introduction

Gallibacterium anatis causes gallibacteriosis in birds. These bacteria produce biofilms and secrete several fimbrial appendages as tools to cause disease in animals. G. anatis strains contain up to three types of fimbriae. Complete genome sequencing is the strategy currently used to determine variations in the gene content of G. anatis, although today only the completely circularized genome of G. anatis UMN179 is available.

Methods

The appearance of growth of various strains of G. anatis in liquid culture medium was studied. Biofilm production and how the amount of biofilm was affected by DNase, Proteinase K, and Pronase E enzymes were analyzed. Fimbrial gene expression was performed by protein analysis and qRT-PCR. In an avian model, the pathogenesis generated by the strains G. anatis ESV200 and 12656-12 was investigated. Using bioinformatic tools, the complete genome of G. anatis ESV200 was comparatively studied to search for virulence factors that would help explain the pathogenic behavior of this strain.

Results and Discussion

G. anatis ESV200 strain differs from the 12656-12 strain because it produces a biofilm at 20%. G. anatis ESV200 strain express fimbrial genes and produces biofilm but with a different structure than that observed for strain 12656-12. ESV200 and 12656-12 strains are pathogenic for chickens, although the latter is the most virulent. Here, we show that the complete genome of the ESV200 strain is similar to that of the UNM179 strain. However, these strains have evolved with many structural rearrangements; the most striking chromosomal arrangement is a Maverick-like element present in the ESV200 strain.

Novel integrative elements and genomic plasticity in ocean ecosystems

Horizontal gene transfer accelerates microbial evolution. The marine picocyanobacterium Prochlorococcus exhibits high genomic plasticity, yet the underlying mechanisms are elusive. Here, we report a novel family of DNA transposons-"tycheposons"-some of which are viral satellites while others carry cargo, such as nutrient-acquisition genes, which shape the genetic variability in this globally abundant genus. Tycheposons share distinctive mobile-lifecycle-linked hallmark genes, including a deep-branching site-specific tyrosine recombinase. Their excision and integration at tRNA genes appear to drive the remodeling of genomic islands-key reservoirs for flexible genes in bacteria. In a selection experiment, tycheposons harboring a nitrate assimilation cassette were dynamically gained and lost, thereby promoting chromosomal rearrangements and host adaptation. Vesicles and phage particles harvested from seawater are enriched in tycheposons, providing a means for their dispersal in the wild. Similar elements are found in microbes co-occurring with Prochlorococcus, suggesting a common mechanism for microbial diversification in the vast oligotrophic oceans.

Deciphering a novel chloramphenicols resistance mechanism: Oxidative inactivation of the propanediol pharmacophore

Expanding knowledge about new types of antibiotic resistance genes is of great significance in dealing with the global antibiotic resistance crisis. Herein, a novel oxidoreductase capO was discovered to be responsible for oxidative inactivation of chloramphenicol and thiamphenicol. The antibiotic resistance mechanism was comprehensively deciphered using multi-omics and multiscale computational approaches. A 66,383 bp DNA fragment carrying capO was shared among four chloramphenicol-resistant strains, and the co-occurrence of capO with a mobile genetic element cluster revealed its potential mobility among different taxa. Metagenomic analysis of 772 datasets indicated that chloramphenicol was the crucial driving factor for the development and accumulation of capO in activated sludge bioreactors treating antibiotic production wastewater. Therefore, we should pay sufficient attention to its possible prevalence and transfer to pathogens, especially in some hotspot environments contaminated with high concentrations of chloramphenicols. This finding significantly expands our knowledge boundary about chloramphenicols resistance mechanisms.

Generating a Small Shuttle Vector for Effective Genetic Engineering of Methanosarcina mazei Allowed First Insights in Plasmid Replication Mechanism in the Methanoarchaeon

Due to their role in methane production, methanoarchaea are of high ecological relevance and genetic systems have been ever more established in the last two decades. The system for protein expression in Methanosarcina using a comprehensive shuttle vector is established; however, details about its replication mechanism in methanoarchaea remain unknown. Here, we report on a significant optimisation of the rather large shuttle vector pWM321 (8.9 kbp) generated by Metcalf through a decrease in its size by about 35% by means of the deletion of several non-coding regions and the ssrA gene. The resulting plasmid (pRS1595) still stably replicates in M. mazei and—most likely due to its reduced size—shows a significantly higher transformation efficiency compared to pWM321. In addition, we investigate the essential gene repA, coding for a rep type protein. RepA was heterologously expressed in Escherichia coli, purified and characterised, demonstrating the significant binding and nicking activity of supercoiled plasmid DNA. Based on our findings we propose that the optimised shuttle vector replicates via a rolling circle mechanism with RepA as the initial replication protein in Methanosarcina. On the basis of bioinformatic comparisons, we propose the presence and location of a double-strand and a single-strand origin, which need to be further verified.

mobileOG-db: a Manually Curated Database of Protein Families Mediating the Life Cycle of Bacterial Mobile Genetic Elements

Bacterial mobile genetic elements (MGEs) encode functional modules that perform both core and accessory functions for the element, the latter of which are often only transiently associated with the element. The presence of these accessory genes, which are often close homologs to primarily immobile genes, incur high rates of false positives and, therefore, limits the usability of these databases for MGE annotation. To overcome this limitation, we analyzed 10,776,849 protein sequences derived from eight MGE databases to compile a comprehensive set of 6,140 manually curated protein families that are linked to the “life cycle” (integration/excision, replication/recombination/repair, transfer, stability/transfer/defense, and phage-specific processes) of plasmids, phages, integrative, transposable, and conjugative elements. We overlay experimental information where available to create a tiered annotation scheme of high-quality annotations and annotations inferred exclusively through bioinformatic evidence. We additionally provide an MGE-class label for each entry (e.g., plasmid or integrative element), and assign to each entry a major and minor category. The resulting database, mobileOG-db (for mobile orthologous groups), comprises over 700,000 deduplicated sequences encompassing five major mobileOG categories and more than 50 minor categories, providing a structured language and interpretable basis for an array of MGE-centered analyses. mobileOG-db can be accessed at mobileogdb.flsi.cloud.vt.edu/, where users can select, refine, and analyze custom subsets of the dynamic mobilome.

IMPORTANCE The analysis of bacterial mobile genetic elements (MGEs) in genomic data is a critical step toward profiling the root causes of antibiotic resistance, phenotypic or metabolic diversity, and the evolution of bacterial genera. Existing methods for MGE annotation pose high barriers of biological and computational expertise to properly harness. To bridge this gap, we systematically analyzed 10,776,849 proteins derived from eight databases of MGEs to identify 6,140 MGE protein families that can serve as candidate hallmarks, i.e., proteins that can be used as “signatures” of MGEs to aid annotation. The resulting resource, mobileOG-db, provides a multilevel classification scheme that encompasses plasmid, phage, integrative, and transposable element protein families categorized into five major mobileOG categories and more than 50 minor categories. mobileOG-db thus provides a rich resource for simple and intuitive element annotation that can be integrated seamlessly into existing MGE detection pipelines and colocalization analyses.

Identification and Characterisation of pST1023 A Mosaic, Multidrug-Resistant and Mobilisable IncR Plasmid

We report the identification and characterisation of a mosaic, multidrug-resistant and mobilisable IncR plasmid (pST1023) detected in Salmonella ST1023, a monophasic variant 4,[5],12:i: strain of widespread pandemic lineage, reported as a Southern European clone. pST1023 contains exogenous DNA regions, principally gained from pSLT-derivatives and IncI1 plasmids. Acquisition from IncI1 included oriT and nikAB and these conferred the ability to be mobilisable in the presence of a helper plasmid, as we demonstrated with the conjugative plasmids pST1007-1D (IncFII) or pVC1035 (IncC). A sul3-associated class 1 integron, conferring resistance to aminoglycosides, chloramphenicol and trimethoprim-sulphonamides, was also embedded in the acquired IncI1 DNA segment. pST1023 also harboured an additional site-specific recombination system (rfsF/rsdB) and IS elements of the IS1, IS5 (IS903 group) and IS6 families. Four of the six IS26 elements present constituted two pseudo-compound-transposons, named PCT-sil and PCT-Tn10 (identified here for the first time). The study further highlighted the mosaic genetic architecture and the clinical importance of IncR plasmids. Moreover, it provides the first experimental data on the ability of IncR plasmids to be mobilised and their potential role in the horizontal spread of antimicrobial-resistant genes.

The First Cbk-Like Phage Infecting Erythrobacter, Representing a Novel Siphoviral Genus

Erythrobacter is an important and widespread bacterial genus in the ocean. However, our knowledge about their phages is still rare. Here, a novel lytic phage vB_EliS-L02, infecting Erythrobacter litoralis DSM 8509, was isolated and purified from Sanggou Bay seawater, China. Morphological observation revealed that the phage belonged to Cbk-like siphovirus, with a long prolate head and a long tail. The host range test showed that phage vB_EliS-L02 could only infect a few strains of Erythrobacter, demonstrating its potential narrow-host range. The genome size of vB_EliS-L02 was 150,063 bp with a G+C content of 59.43%, encoding 231 putative open reading frames (ORFs), but only 47 were predicted to be functional domains. Fourteen auxiliary metabolic genes were identified, including phoH that may confer vB_EliS-L02 the advantage of regulating phosphate uptake and metabolism under a phosphate-limiting condition. Genomic and phylogenetic analyses indicated that vB_EliS-L02 was most closely related to the genus Lacusarxvirus with low similarity (shared genes < 30%, and average nucleotide sequence identity < 70%), distantly from other reported phages, and could be grouped into a novel viral genus cluster, in this study as Eliscbkvirus. Meanwhile, the genus Eliscbkvirus and Lacusarxvirus stand out from other siphoviral genera and could represent a novel subfamily within Siphoviridae, named Dolichocephalovirinae-II. Being a representative of an understudied viral group with manifold adaptations to the host, phage vB_EliS-L02 could improve our understanding of the virus–host interactions and provide reference information for viral metagenomic analysis in the ocean.

Giant Starship Elements Mobilize Accessory Genes in Fungal Genomes

Accessory genes are variably present among members of a species and are a reservoir of adaptive functions. In bacteria, differences in gene distributions among individuals largely result from mobile elements that acquire and disperse accessory genes as cargo. In contrast, the impact of cargo-carrying elements on eukaryotic evolution remains largely unknown. Here, we show that variation in genome content within multiple fungal species is facilitated by Starships, a newly discovered group of massive mobile elements that are 110 kb long on average, share conserved components, and carry diverse arrays of accessory genes. We identified hundreds of Starship-like regions across every major class of filamentous Ascomycetes, including 28 distinct Starships that range from 27 to 393 kb and last shared a common ancestor ca. 400 Ma. Using new long-read assemblies of the plant pathogen Macrophomina phaseolina, we characterize four additional Starships whose activities contribute to standing variation in genome structure and content. One of these elements, Voyager, inserts into 5S rDNA and contains a candidate virulence factor whose increasing copy number has contrasting associations with pathogenic and saprophytic growth, suggesting Voyager's activity underlies an ecological trade-off. We propose that Starships are eukaryotic analogs of bacterial integrative and conjugative elements based on parallels between their conserved components and may therefore represent the first dedicated agents of active gene transfer in eukaryotes. Our results suggest that Starships have shaped the content and structure of fungal genomes for millions of years and reveal a new concerted route for evolution throughout an entire eukaryotic phylum.

Landscape of mobile genetic elements and their antibiotic resistance cargo in prokaryotic genomes

Prokaryotic Mobile Genetic Elements (MGEs) such as transposons, integrons, phages and plasmids, play important roles in prokaryotic evolution and in the dispersal of cargo functions like antibiotic resistance. However, each of these MGE types is usually annotated and analysed individually, hampering a global understanding of phylogenetic and environmental patterns of MGE dispersal. We thus developed a computational framework that captures diverse MGE types, their cargos and MGE-mediated horizontal transfer events, using recombinases as ubiquitous MGE marker genes and pangenome information for MGE boundary estimation. Applied to ∼84k genomes with habitat annotation, we mapped 2.8 million MGE-specific recombinases to six operational MGE types, which together contain on average 13% of all the genes in a genome. Transposable elements (TEs) dominated across all taxa (∼1.7 million occurrences), outnumbering phages and phage-like elements (<0.4 million). We recorded numerous MGE-mediated horizontal transfer events across diverse phyla and habitats involving all MGE types, disentangled and quantified the extent of hitchhiking of TEs (17%) and integrons (63%) with other MGE categories, and established TEs as dominant carriers of antibiotic resistance genes. We integrated all these findings into a resource (proMGE.embl.de), which should facilitate future studies on the large mobile part of genomes and its horizontal dispersal.

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.

IntegronFinder 2.0: Identification and Analysis of Integrons across Bacteria, with a Focus on Antibiotic Resistance in Klebsiella

Integrons are flexible gene-exchanging platforms that contain multiple cassettes encoding accessory genes whose order is shuffled by a specific integrase. Integrons embedded within mobile genetic elements often contain multiple antibiotic resistance genes that they spread among nosocomial pathogens and contribute to the current antibiotic resistance crisis. However, most integrons are presumably sedentary and encode a much broader diversity of functions. IntegronFinder is a widely used software to identify novel integrons in bacterial genomes, but has aged and lacks some useful functionalities to handle very large datasets of draft genomes or metagenomes. Here, we present IntegronFinder version 2. We have updated the code, improved its efficiency and usability, adapted the output to incomplete genome data, and added a few novel functions. We describe these changes and illustrate the relevance of the program by analyzing the distribution of integrons across more than 20,000 fully sequenced genomes. We also take full advantage of its novel capabilities to analyze close to 4000 Klebsiella pneumoniae genomes for the presence of integrons and antibiotic resistance genes within them. Our data show that K. pneumoniae has a large diversity of integrons and the largest mobile integron in our database of plasmids. The pangenome of these integrons contains a total of 165 different gene families with most of the largest families being related with resistance to numerous types of antibiotics. IntegronFinder is a free and open-source software available on multiple public platforms.

Towards the biogeography of prokaryotic genes

Microbial genes encode the majority of the functional repertoire of life on earth. However, despite increasing efforts in metagenomic sequencing of various habitats^1-3, little is known about the distribution of genes across the global biosphere, with implications for human and planetary health. Here we constructed a non-redundant gene catalogue of 303 million species-level genes (clustered at 95% nucleotide identity) from 13,174 publicly available metagenomes across 14 major habitats and use it to show that most genes are specific to a single habitat. The small fraction of genes found in multiple habitats is enriched in antibiotic-resistance genes and markers for mobile genetic elements. By further clustering these species-level genes into 32 million protein families, we observed that a small fraction of these families contain the majority of the genes (0.6% of families account for 50% of the genes). The majority of species-level genes and protein families are rare. Furthermore, species-level genes, and in particular the rare ones, show low rates of positive (adaptive) selection, supporting a model in which most genetic variability observed within each protein family is neutral or nearly neutral.

Conformational dynamics promotes disordered regions from function-dispensable to essential in evolved site-specific DNA recombinases

•

New functional regions emerging in evolution of DNA site-specific recombinase tails.

•

Transient structural nucleation promotes function-dispensable regions to essential.

•

Molecular dynamics reveals conformational diversity and its functional implications.

•

Evolved disordered molecular mechanisms of N-term tails for protein stability.

•

Structural disorder-based link between protein evolution, stability and function.

Protein intrinsically disordered regions (IDRs) play pivotal roles in molecular recognition and regulatory processes through structural disorder-to-order transitions. To understand and exploit the distinctive functional implications of IDRs and to unravel the underlying molecular mechanisms, structural disorder-to-function relationships need to be deciphered. The DNA site-specific recombinase system Cre/loxP represents an attractive model to investigate functional molecular mechanisms of IDRs. Cre contains a functionally dispensable disordered N-terminal tail, which becomes indispensable in the evolved Tre/loxLTR recombinase system. The difficulty to experimentally obtain structural information about this tail has so far precluded any mechanistic study on its involvement in DNA recombination. Here, we use in vitro and in silico evolution data, conformational dynamics, AI-based folding simulations, thermodynamic stability calculations, mutagenesis and DNA recombination assays to investigate how evolution and the dynamic behavior of this IDR may determine distinct functional properties. Our studies suggest that partial conformational order in the N-terminal tail of Tre recombinase and its packing to a conserved hydrophobic surface on the protein provide thermodynamic stability. Based on our results, we propose a link between protein stability and function, offering new plausible atom-detailed mechanistic insights into disorder-function relationships. Our work highlights the potential of N-terminal tails to be exploited for regulation of the activity of Cre-like tyrosine-type SSRs, which merits future investigations and could be of relevance in future rational engineering for their use in biotechnology and genomic medicine.