Activity, specificity and structure of I-Bth0305I: a representative of a new homing endonuclease family

Nuclease genes occupy boundaries of genetic exchange between bacteriophages

Homing endonuclease genes (HEGs) are ubiquitous selfish elements that generate targeted double-stranded DNA breaks, facilitating the recombination of the HEG DNA sequence into the break site and contributing to the evolutionary dynamics of HEG-encoding genomes. Bacteriophages (phages) are well-documented to carry HEGs, with the paramount characterization of HEGs being focused on those encoded by coliphage T4. Recently, it has been observed that the highly sampled vibriophage, ICP1, is similarly enriched with HEGs distinct from T4's. Here, we examined the HEGs encoded by ICP1 and diverse phages, proposing HEG-driven mechanisms that contribute to phage evolution. Relative to ICP1 and T4, we found a variable distribution of HEGs across phages, with HEGs frequently encoded proximal to or within essential genes. We identified large regions (> 10kb) of high nucleotide identity flanked by HEGs, deemed HEG islands, which we hypothesize to be mobilized by the activity of flanking HEGs. Finally, we found examples of domain swapping between phage-encoded HEGs and genes encoded by other phages and phage satellites. We anticipate that HEGs have a larger impact on the evolutionary trajectory of phages than previously appreciated and that future work investigating the role of HEGs in phage evolution will continue to highlight these observations.

Nuclease genes occupy boundaries of genetic exchange between bacteriophages

Homing endonuclease genes (HEGs) are ubiquitous selfish elements that generate targeted double-stranded DNA breaks, facilitating the recombination of the HEG DNA sequence into the break site and contributing to the evolutionary dynamics of HEG-encoding genomes. Bacteriophages (phages) are well-documented to carry HEGs, with the paramount characterization of HEGs being focused on those encoded by coliphage T4. Recently, it has been observed that the highly sampled vibriophage, ICP1, is similarly enriched with HEGs distinct from T4’s. Here, we examined the HEGs encoded by ICP1 and diverse phages, proposing HEG-driven mechanisms that contribute to phage evolution. Relative to ICP1 and T4, we found a variable distribution of HEGs across phages, with HEGs frequently encoded proximal to or within essential genes. We identified large regions (> 10kb) of high nucleotide identity flanked by HEGs, deemed HEG islands, which we hypothesize to be mobilized by the activity of flanking HEGs. Finally, we found examples of domain swapping between phage-encoded HEGs and genes encoded by other phages and phage satellites. We anticipate that HEGs have a larger impact on the evolutionary trajectory of phages than previously appreciated and that future work investigating the role of HEGs in phage evolution will continue to highlight these observations.

Comparative analysis reveals distinctive genomic features of Taiwan hot-spring cyanobacterium Thermosynechococcus sp. TA-1

Thermosynechococcus is a genus of thermophilic unicellular cyanobacteria that dominates microbial mats in Asian non-acidic hot springs. These cyanobacteria are the major primary producers in their ecological niches and are promising sources of thermostable enzymes for biotechnology applications. To improve our understanding of these organisms, we conducted whole-genome sequencing of a novel strain for comparative analysis with other representatives in the same genus. This newly characterized strain, Thermosynechococcus sp. TA-1, was isolated from the Taian hot springs in Taiwan. Analyses based on average nucleotide identity (ANI) and genome-scale phylogeny suggested that TA-1 and another Taiwanese strain CL-1 belong to a novel species-level taxon. Two metagenome-assembled genomes (MAGs) originated from India represent the sister group, and Thermosynechococcus elongatus PKUAC-SCTE542 from China is the next closest lineage. All cultivated strains and MAGs from Japan form a separate monophyletic clade and could be classified into two species-level taxa. Intriguingly, although TA-1 and CL-1 share 97.0% ANI, the genome alignment identified at least 16 synteny breakpoints that are mostly associated with transposase genes, which illustrates the dynamic nature of their chromosomal evolution. Gene content comparisons identified multiple features distinct at species- or strain-level among these Thermosynechococcus representatives. Examples include genes involved in bicarbonate transportation, nitric oxide protection, urea utilization, kanamycin resistance, restriction-modification system, and chemotaxis. Moreover, we observed the insertion of type II inteins in multiple genes of the two Taiwanese strains and inferred putative horizontal transfer of an asparagine synthase gene (asnB) associated with exopolysaccharides gene cluster. Taken together, while previous work suggested that strains in this genus share a highly conserved genomic core and no clear genetic differentiation could be linked to environmental factors, we found that the overall pattern of gene content divergence is largely congruent with core genome phylogeny. However, it is difficult to distinguish between the roles of phylogenetic relatedness and geographic proximity in shaping the genetic differentiation. In conclusion, knowledge of the genomic differentiation among these strains provides valuable resources for future functional characterization.

Comparative Genomics of Three Novel Jumbo Bacteriophages Infecting Staphylococcus aureus

The majority of previously described Staphylococcus aureus bacteriophages belong to three major groups: P68-like podophages, Twort-like or K-like myophages, and a more diverse group of temperate siphophages. Here we present three novel S. aureus "jumbo" phages: MarsHill, Madawaska, and Machias. These phages were isolated from swine production environments in the United States and represent a novel clade of S. aureus myophage. The average genome size for these phages is ∼269 kb with each genome encoding ∼263 predicted protein-coding genes. Phage genome organization and content is similar to known jumbo phages of Bacillus, including AR9 and vB_BpuM-BpSp. All three phages possess genes encoding complete virion and non-virion RNA polymerases, multiple homing endonucleases, and a retron-like reverse transcriptase. Like AR9, all of these phages are presumed to have uracil-substituted DNA which interferes with DNA sequencing. These phages are also able to transduce host plasmids, which is significant as these phages were found circulating in swine production environments and can also infect human S. aureus isolates. Importance of work: This study describes the comparative genomics of three novel S. aureus jumbo phages: MarsHill, Madawaska, and Machias. These three S. aureus myophages represent an emerging class of S. aureus phage. These genomes contain abundant introns which show a pattern consistent with repeated acquisition rather than vertical inheritance, suggesting intron acquisition and loss is an active process in the evolution of these phages. These phages have presumably hypermodified DNA which inhibits sequencing by several different common platforms. Therefore, these phages also represent potential genomic diversity that has been missed due to the limitations of standard sequencing techniques. In particular, such hypermodified genomes may be missed by metagenomic studies due to their resistance to standard sequencing techniques. Phage MarsHill was found to be able to transduce host DNA at levels comparable to that found for other transducing S. aureus phages, making them a potential vector for horizontal gene transfer in the environment.

Metal effect on intein splicing: A review

Inteins are intervening polypeptides that interrupt the functional domains of several important proteins across the three domains of life. Inteins excise themselves from the precursor protein, ligating concomitant extein residues in a process called protein splicing. Post-translational auto-removal of inteins remain critical for the generation of active proteins. The perspective of inteins in science is a robust field of research, however fundamental studies centralized upon splicing regulatory mechanism are imperative for addressing more intricate issues. Controlled engineering of intein splicing has many applications; intein inhibition can facilitate novel drug design, while activation of intein splicing is exploited in protein purification. This paper provides a comprehensive review of the past and recent advances in the splicing regulation via metal-intein interaction. We compare the behavior of different metal ions on diverse intein systems. Though metals such as Zn, Cu, Pt, Cd, Co, Ni exhibit intein inhibitory effect heterogeneously on different inteins, divalent metal ions such as Ca and Mg fail to do so. The observed diversity in the metal-intein interaction arises mostly due to intein polymorphism and variations in atomic structure of metals. A mechanistic understanding of intein regulation by metals in native as well as synthetically engineered intein systems may yield potent intein inhibitors via direct or indirect approach.

Engineering and functionalization of large circular tandem repeat protein nanoparticles

Protein engineering has enabled the design of molecular scaffolds that display a wide variety of sizes, shapes, symmetries and subunit compositions. Symmetric protein-based nanoparticles that display multiple protein domains can exhibit enhanced functional properties due to increased avidity and improved solution behavior and stability. Here we describe the creation and characterization of a computationally designed circular tandem repeat protein (cTRP) composed of 24 identical repeated motifs, which can display a variety of functional protein domains (cargo) at defined positions around its periphery. We demonstrate that cTRP nanoparticles can self-assemble from smaller individual subunits, can be produced from prokaryotic and human expression platforms, can employ a variety of cargo attachment strategies and can be used for applications (such as T-cell culture and expansion) requiring high-avidity molecular interactions on the cell surface.

F-CphI represents a new homing endonuclease family using the Endo VII catalytic motif

Background

There are six known families of homing endonucleases, LAGLIDADG, GIY-YIG, HNH, His-Cys box, PD-(D/E)-XK, and EDxHD, which are characterized by their conserved residues. Previously, we discovered a novel homing endonuclease F-CphI encoded by ORF177 of cyanophage S-PM2. F-CphI does not resemble any characterized homing endonucleases. Instead, the C-terminus of F-CphI aligns well with the N-terminal catalytic domain of a Holliday junction DNA resolvase, phage T4 endonuclease VII (Endo VII).

Results

A PSI-BLAST search resulted in a total of 313 Endo VII motif–containing sequences in sequenced genomes. Multiple sequence alignment showed that the catalytically important residues of T4 Endo VII were all well conserved in these proteins. Our site-directed mutagenesis studies further confirmed that the catalytically important residues of T4 Endo VII were also essential for F-CphI activity, and thus F-CphI might use a similar protein fold as Endo VII for DNA cleavage. A phylogenetic tree of the Endo VII motif–containing sequences showed that putative resolvases grouped into one clade while putative homing endonucleases and restriction endonucleases grouped into another clade.

Conclusions

Based on the unique conserved residues, we proposed that F-CphI represents a new homing endonuclease family, which was named the DHHRN family. Our phylogenetic analysis could be used to predict the functions of many previously unknown proteins.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0132-5) contains supplementary material, which is available to authorized users.

Resilience of biochemical activity in protein domains in the face of structural divergence

Recent studies point to the prevalence of the evolutionary phenomenon of drastic structural transformation of protein domains while continuing to preserve their basic biochemical function. These transformations span a wide spectrum, including simple domains incorporated into larger structural scaffolds, changes in the structural core, major active site shifts, topological rewiring and extensive structural transmogrifications. Proteins from biological conflict systems, such as toxin-antitoxin, restriction-modification, CRISPR/Cas, polymorphic toxin and secondary metabolism systems commonly display such transformations. These include endoDNases, metal-independent RNases, deaminases, ADP ribosyltransferases, immunity proteins, kinases and E1-like enzymes. In eukaryotes such transformations are seen in domains involved in chromatin-related peptide recognition and protein/DNA-modification. Intense selective pressures from 'arms-race'-like situations in conflict and macromolecular modification systems could favor drastic structural divergence while preserving function.

Structure of 5-hydroxymethylcytosine-specific restriction enzyme, AbaSI, in complex with DNA

AbaSI, a member of the PvuRts1I-family of modification-dependent restriction endonucleases, cleaves deoxyribonucleic acid (DNA) containing 5-hydroxymethylctosine (5hmC) and glucosylated 5hmC (g5hmC), but not DNA containing unmodified cytosine. AbaSI has been used as a tool for mapping the genomic locations of 5hmC, an important epigenetic modification in the DNA of higher organisms. Here we report the crystal structures of AbaSI in the presence and absence of DNA. These structures provide considerable, although incomplete, insight into how this enzyme acts. AbaSI appears to be mainly a homodimer in solution, but interacts with DNA in our structures as a homotetramer. Each AbaSI subunit comprises an N-terminal, Vsr-like, cleavage domain containing a single catalytic site, and a C-terminal, SRA-like, 5hmC-binding domain. Two N-terminal helices mediate most of the homodimer interface. Dimerization brings together the two catalytic sites required for double-strand cleavage, and separates the 5hmC binding-domains by ∼70 Å, consistent with the known activity of AbaSI which cleaves DNA optimally between symmetrically modified cytosines ∼22 bp apart. The eukaryotic SET and RING-associated (SRA) domains bind to DNA containing 5-methylcytosine (5mC) in the hemi-methylated CpG sequence. They make contacts in both the major and minor DNA grooves, and flip the modified cytosine out of the helix into a conserved binding pocket. In contrast, the SRA-like domain of AbaSI, which has no sequence specificity, contacts only the minor DNA groove, and in our current structures the 5hmC remains intra-helical. A conserved, binding pocket is nevertheless present in this domain, suitable for accommodating 5hmC and g5hmC. We consider it likely, therefore, that base-flipping is part of the recognition and cleavage mechanism of AbaSI, but that our structures represent an earlier, pre-flipped stage, prior to actual recognition.

Homing endonucleases from mobile group I introns: discovery to genome engineering

Homing endonucleases are highly specific DNA cleaving enzymes that are encoded within genomes of all forms of microbial life including phage and eukaryotic organelles. These proteins drive the mobility and persistence of their own reading frames. The genes that encode homing endonucleases are often embedded within self-splicing elements such as group I introns, group II introns and inteins. This combination of molecular functions is mutually advantageous: the endonuclease activity allows surrounding introns and inteins to act as invasive DNA elements, while the splicing activity allows the endonuclease gene to invade a coding sequence without disrupting its product. Crystallographic analyses of representatives from all known homing endonuclease families have illustrated both their mechanisms of action and their evolutionary relationships to a wide range of host proteins. Several homing endonucleases have been completely redesigned and used for a variety of genome engineering applications. Recent efforts to augment homing endonucleases with auxiliary DNA recognition elements and/or nucleic acid processing factors has further accelerated their use for applications that demand exceptionally high specificity and activity.

Homing endonucleases: from genetic anomalies to programmable genomic clippers

Homing endonucleases are strong drivers of genetic exchange and horizontal transfer of both their own genes and their local genetic environment. The mechanisms that govern the function and evolution of these genetic oddities have been well documented over the past few decades at the genetic, biochemical, and structural levels. This wealth of information has led to the manipulation and reprogramming of the endonucleases and to their exploitation in genome editing for use as therapeutic agents, for insect vector control and in agriculture. In this chapter we summarize the molecular properties of homing endonucleases and discuss their strengths and weaknesses in genome editing as compared to other site-specific nucleases such as zinc finger endonucleases, TALEN, and CRISPR-derived endonucleases.

Group I introns in Staphylococcus bacteriophages

The abundance of group I introns, intragenic RNA sequences capable of self-splicing, in Gram-positive bacteriophage genomes, is illustrated by various new group I introns recently described in Staphylococcus phage genomes. These introns were found to interrupt DNA metabolism genes as well as late genes. These group I introns often code for homing endonucleases, which promote lateral transfer of group I introns, thereby enabling spread through a population. Homing endonucleases encoded by group I introns in Staphylococcus phage genomes were predicted to belong to the GIY-YIG, LAGLIDADG, HNH or EDxHD family of endonucleases. The group I intron distribution in Staphylococcus phage genomes exemplifies the homology between these introns as well as the encoded endonucleases. Despite several suggested functions, the role of group I introns in bacteriophages remains unclear or might be nonexistent. However, transcriptome analysis might provide additional information to elucidate the possible purpose of group I introns in phage genomes.

Romulus and Remus, two phage isolates representing a distinct clade within the Twortlikevirus genus, display suitable properties for phage therapy applications

The renewed interest in controlling Staphylococcus aureus infections using their natural enemies, bacteriophages, has led to the isolation of a limited number of virulent phages so far. These phages are all members of the Twortlikevirus, displaying little variance. We present two novel closely related (95.9% DNA homology) lytic myoviruses, Romulus and Remus, with double-stranded DNA (dsDNA) genomes of 131,333 bp and 134,643 bp, respectively. Despite their relatedness to Staphylococcus phages K, G1, ISP, and Twort and Listeria phages A511 and P100, Romulus and Remus can be proposed as isolates of a new species within the Twortlikevirus genus. A distinguishing feature for these phage genomes is the unique distribution of group I introns compared to that in other staphylococcal myoviruses. In addition, a hedgehog/intein domain was found within their DNA polymerase genes, and an insertion sequence-encoded transposase exhibits splicing behavior and produces a functional portal protein. From a phage therapy application perspective, Romulus and Remus infected approximately 70% of the tested S. aureus isolates and displayed promising lytic activity against these isolates. Furthermore, both phages showed a rapid initial adsorption and demonstrated biofilm-degrading capacity in a proof-of-concept experiment.

Natural zinc ribbon HNH endonucleases and engineered zinc finger nicking endonuclease

Many bacteriophage and prophage genomes encode an HNH endonuclease (HNHE) next to their cohesive end site and terminase genes. The HNH catalytic domain contains the conserved catalytic residues His-Asn-His and a zinc-binding site [CxxC]₂. An additional zinc ribbon (ZR) domain with one to two zinc-binding sites ([CxxxxC], [CxxxxH], [CxxxC], [HxxxH], [CxxC] or [CxxH]) is frequently found at the N-terminus or C-terminus of the HNHE or a ZR domain protein (ZRP) located adjacent to the HNHE. We expressed and purified 10 such HNHEs and characterized their cleavage sites. These HNHEs are site-specific and strand-specific nicking endonucleases (NEase or nickase) with 3- to 7-bp specificities. A minimal HNH nicking domain of 76 amino acid residues was identified from Bacillus phage γ HNHE and subsequently fused to a zinc finger protein to generate a chimeric NEase with a new specificity (12–13 bp). The identification of a large pool of previously unknown natural NEases and engineered NEases provides more ‘tools’ for DNA manipulation and molecular diagnostics. The small modular HNH nicking domain can be used to generate rare NEases applicable to targeted genome editing. In addition, the engineered ZF nickase is useful for evaluation of off-target sites in vitro before performing cell-based gene modification.

A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity

Clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated (Cas) systems provide bacteria and archaea with adaptive immunity against viruses and plasmids by using CRISPR RNAs (crRNAs) to guide the silencing of invading nucleic acids. We show here that in a subset of these systems, the mature crRNA that is base-paired to trans-activating crRNA (tracrRNA) forms a two-RNA structure that directs the CRISPR-associated protein Cas9 to introduce double-stranded (ds) breaks in target DNA. At sites complementary to the crRNA-guide sequence, the Cas9 HNH nuclease domain cleaves the complementary strand, whereas the Cas9 RuvC-like domain cleaves the noncomplementary strand. The dual-tracrRNA:crRNA, when engineered as a single RNA chimera, also directs sequence-specific Cas9 dsDNA cleavage. Our study reveals a family of endonucleases that use dual-RNAs for site-specific DNA cleavage and highlights the potential to exploit the system for RNA-programmable genome editing.

A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity

Clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated (Cas) systems provide bacteria and archaea with adaptive immunity against viruses and plasmids by using CRISPR RNAs (crRNAs) to guide the silencing of invading nucleic acids. We show here that in a subset of these systems, the mature crRNA that is base-paired to trans-activating crRNA (tracrRNA) forms a two-RNA structure that directs the CRISPR-associated protein Cas9 to introduce double-stranded (ds) breaks in target DNA. At sites complementary to the crRNA-guide sequence, the Cas9 HNH nuclease domain cleaves the complementary strand, whereas the Cas9 RuvC-like domain cleaves the noncomplementary strand. The dual-tracrRNA:crRNA, when engineered as a single RNA chimera, also directs sequence-specific Cas9 dsDNA cleavage. Our study reveals a family of endonucleases that use dual-RNAs for site-specific DNA cleavage and highlights the potential to exploit the system for RNA-programmable genome editing.

Homing endonucleases: DNA scissors on a mission

Buried within the genomes of many microorganisms are genetic elements that encode rare-cutting homing endonucleases that assist in the mobility of the elements that encode them, such as the self-splicing group I and II introns and in some cases inteins. There are several different families of homing endonucleases and their ability to initiate and target specific sequences for lateral transfers makes them attractive reagents for gene targeting. Homing endonucleases have been applied in promoting DNA modification or genome editing such as gene repair or "gene knockouts". This review examines the categories of homing endonucleases that have been described so far and their possible applications to biotechnology. Strategies to engineer homing endonucleases to alter target site specificities will also be addressed. Alternatives to homing endonucleases such as zinc finger nucleases, transcription activator-like effector nucleases, triplex forming oligonucleotide nucleases, and targetrons are also briefly discussed.

Sequence, structure and functional diversity of PD-(D/E)XK phosphodiesterase superfamily

Proteins belonging to PD-(D/E)XK phosphodiesterases constitute a functionally diverse superfamily with representatives involved in replication, restriction, DNA repair and tRNA-intron splicing. Their malfunction in humans triggers severe diseases, such as Fanconi anemia and Xeroderma pigmentosum. To date there have been several attempts to identify and classify new PD-(D/E)KK phosphodiesterases using remote homology detection methods. Such efforts are complicated, because the superfamily exhibits extreme sequence and structural divergence. Using advanced homology detection methods supported with superfamily-wide domain architecture and horizontal gene transfer analyses, we provide a comprehensive reclassification of proteins containing a PD-(D/E)XK domain. The PD-(D/E)XK phosphodiesterases span over 21,900 proteins, which can be classified into 121 groups of various families. Eleven of them, including DUF4420, DUF3883, DUF4263, COG5482, COG1395, Tsp45I, HaeII, Eco47II, ScaI, HpaII and Replic_Relax, are newly assigned to the PD-(D/E)XK superfamily. Some groups of PD-(D/E)XK proteins are present in all domains of life, whereas others occur within small numbers of organisms. We observed multiple horizontal gene transfers even between human pathogenic bacteria or from Prokaryota to Eukaryota. Uncommon domain arrangements greatly elaborate the PD-(D/E)XK world. These include domain architectures suggesting regulatory roles in Eukaryotes, like stress sensing and cell-cycle regulation. Our results may inspire further experimental studies aimed at identification of exact biological functions, specific substrates and molecular mechanisms of reactions performed by these highly diverse proteins.

Structural, functional and evolutionary relationships between homing endonucleases and proteins from their host organisms

Homing endonucleases (HEs) are highly specific DNA-cleaving enzymes that are encoded by invasive DNA elements (usually mobile introns or inteins) within the genomes of phage, bacteria, archea, protista and eukaryotic organelles. Six unique structural HE families, that collectively span four distinct nuclease catalytic motifs, have been characterized to date. Members of each family display structural homology and functional relationships to a wide variety of proteins from various organisms. The biological functions of those proteins are highly disparate and include non-specific DNA-degradation enzymes, restriction endonucleases, DNA-repair enzymes, resolvases, intron splicing factors and transcription factors. These relationships suggest that modern day HEs share common ancestors with proteins involved in genome fidelity, maintenance and gene expression. This review summarizes the results of structural studies of HEs and corresponding proteins from host organisms that have illustrated the manner in which these factors are related.