Inteins, often referred to as protein introns, are highly mobile genetic elements that invade conserved genes throughout the tree of life. Inteins have been found to invade a wide variety of key genes within actinophages. While in the process of conducting a survey of these inteins in actinophages, we discovered that one protein family of methylases contained a putative intein, and two other unique insertion elements. These methylases are known to occur commonly in phages as orphan methylases (possibly as a form of resistance to restriction–modification systems). We found that the methylase family is not conserved within phage clusters and has a disparate distribution across divergent phage groups. We determined that two of the three insertion elements have a patchy distribution within the methylase protein family. Additionally, we found that the third insertion element is likely a second homing endonuclease, and that all three elements (the intein, the homing endonuclease, and what we refer to as the ShiLan domain) have different insertion sites that are conserved in the methylase gene family. Furthermore, we find strong evidence that both the intein and ShiLan domain are partaking in long-distance horizontal gene transfer events between divergent methylases in disparate phage hosts within the already dispersed methylase distribution. The reticulate evolutionary history of methylases and their insertion elements reveals high rates of gene transfer and within-gene recombination in actinophages.

Horizontal gene transfer (HGT) is fundamental to microbial evolution and adaptation. When a gene is horizontally transferred, it may either add itself as a new gene to the recipient genome (possibly displacing nonhomologous genes) or replace an existing homologous gene. Currently, studies do not usually distinguish between “additive” and “replacing” HGTs, and their relative frequencies, integration mechanisms, and specific roles in microbial evolution are poorly understood. In this work, we develop a novel computational framework for large-scale classification of HGTs as either additive or replacing. Our framework leverages recently developed phylogenetic approaches for HGT detection and classifies HGTs inferred between terminal edges based on gene orderings along genomes and phylogenetic relationships between the microbial species under consideration. The resulting method, called DART, is highly customizable and scalable and can classify a large fraction of inferred HGTs with high confidence and statistical support. Our application of DART to a large dataset of thousands of gene families from 103 Aeromonas genomes provides insights into the relative frequencies, functional biases, and integration mechanisms of additive and replacing HGTs. Among other results, we find that (i) the relative frequency of additive HGT increases with increasing phylogenetic distance, (ii) replacing HGT dominates at shorter phylogenetic distances, (iii) additive and replacing HGTs have strikingly different functional profiles, (iv) homologous recombination in flanking regions of a novel gene may be a frequent integration mechanism for additive HGT, and (v) phages and mobile genetic elements likely play an important role in facilitating additive HGT.

Background/Objectives: The ShiLan domain was previously identified as an insertion sequence in a phage DNA methylase gene that exhibited similar evolutionary patterns to that of an active intein or self-splicing intron but could not be identified as either. It produces no internal stop codons when read in frame with its host methylase gene, leading to the thought that it may not be an intron and rather be an abnormal type of intein. However, the sequence has no detectable self-splicing domains, which are essential for intein persistence, as preventing an intein from successfully splicing is often detrimental to proper host protein function. Methods: The analysis of alternate open reading frames for the full nucleotide sequence of this insertion element revealed the insertion to be an out-of-frame histidine-asparagine-histidine (HNH) endonuclease. A GTG start codon is located 18 bp into the insertion, and a TAA stop codon within the last four bases of the insertion (TAAC). When this frame is read, an HNH endonuclease is revealed. In-depth computational analysis could not retrieve support for this element being any known type of self-splicing element, neither intein nor intron. When read in-frame with the methylase gene, this insertion is predicted to take on a looping structure that may be able to avoid interference with the DNA methylase activity. We performed searches for sequences similar in nature to the inserted out-of-frame HNH and found several in other phages and prokaryotes. We present our survey of these out-of-frame endonuclease insertion elements as well as some speculation on how these endonucleases are getting translated to facilitate their homing activity. Conclusions: These findings expand our understanding of the possible arrangements for and prevalence of unorthodox mobile genetic elements and overlapping open reading frames in phages.

Intein presence within Actinobacteriophages (within PhagesDB) was last surveyed in 2016, and despite a 5-fold increase in the size of the database, has not been updated since. To address this, we present a modern survey of the current iteration of the PhagesDB database. We developed a new algorithm — Iterative Cluster Expansion BLAST (ICE-BLAST) — to expand our search to more divergent sequences. Nearly 800 inteins were retrieved through this process; the majority of which were previously unreported. We describe the nature of these inteins, their classes, integration target sites, distribution within phage clusters, and explore the geographical location of nearly identical intein sequences found in divergent exteins. Our findings suggest that these inteins recently invaded local phage populations. We also find two instances of a Cas4 exonuclease intein evolving from a terminase large subunit intein, and propose a model by which one of these inteins was able to utilize sequence similarity conferred by a shared nucleotide binding site to jump between genes. Additionally, we find inteins with never-before-reported homing endonucleases, and inteins with homing endonucleases encoded in a reading frame separate from that which encodes the extein and the intein’s self-splicing domain. We provide predicted structures for these elements and hypothesize on their evolution and relation to free-standing homing endonucleases within phage genomes. Finally, we provide evidence that these “non-canonical” inteins are still transferring between host genomes, in a fashion similar to other inteins with canonical homing endonucleases within the dataset.

Many double-stranded DNA viruses, including tailed bacteriophages (phages) and herpesviruses, use the HK97-fold in their major capsid protein to make the capsomers of the icosahedral viral capsid. Following the genome packaging at near-crystalline densities, the capsid is subjected to a major expansion and stabilization step that allows it to withstand environmental stresses and internal high pressure. Several different mechanisms for stabilizing the capsid have been structurally characterized, but how these mechanisms have evolved is still not understood. Using cryo-EM structure determination, structural comparisons, phylogenetic analyses, and Alphafold predictions, we have constructed a detailed structural dendrogram describing the evolution of capsid structural stability within the actinobacteriophages. The cryo-EM reconstructions of ten capsids solved to resolutions between 2.2 and 4 Ångstroms revealed that eight of them exhibit major capsid proteins that are linked by a covalent cross-linking (isopeptide bond) between subunits that was first described in the HK97 phage. Those covalent interactions ultimately lead to the formation of mutually interlinked capsomers that has been compared to the structure of chain mail. However, three of the closely related phages do not exhibit such an isopeptide bond as demonstrated by both our cryo-EM maps and the lack of the required residue. This work raises questions about the importance of previously described capsid stabilization mechanisms.

Inteins, often referred to as protein introns, are highly mobile genetic elements that invade conserved genes throughout the tree of life. Inteins have been found to invade a wide variety of key genes within actinophages. While in the process of conducting a survey of these inteins in actinophages we discovered that one protein family of methylases contained a putative intein, and two other unique insertion elements. These methylases are known to occur commonly in phages as orphan methylases (possibly as a form of resistance to restriction-modification systems). We found that the methylase family is not conserved within phage clusters and has a disparate distribution across divergent phage groups. We determined that two of the three insertion elements have a patchy distribution within the methylase protein family. We also found that the third insertion element is likely a second homing endonuclease, and that all three elements (the intein, the homing endonuclease, and what we refer to as the ShiLan domain) all have different insertion sites that are conserved in the methylase gene family. Furthermore, we find strong evidence that both the intein and ShiLan domain are partaking in long distance horizontal gene transfer events between divergent methylases in disparate phage hosts within the already dispersed methylase distribution. The reticulate evolutionary history of methylases and their insertion elements reveals high rates of gene transfer and within-gene recombination in actinophages.