A free-standing homing endonuclease targets an intron insertion site in the psbA gene of cyanophages

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.

Abundant and cosmopolitan lineage of cyanopodoviruses lacking a DNA polymerase gene

Cyanopodoviruses affect the mortality and population dynamics of the unicellular picocyanobacteria Prochlorococcus and Synechococcus, the dominant primary producers in the oceans. Known cyanopodoviruses all contain the DNA polymerase gene (DNA pol) that is important for phage DNA replication and widely used in field quantification and diversity studies. However, we isolated 18 cyanopodoviruses without identifiable DNA pol. They form a new MPP-C clade that was separated from the existing MPP-A, MPP-B, and P-RSP2 clades. The MPP-C phages have the smallest genomes (37.3-37.9 kb) among sequenced cyanophages, and show longer latent periods than the MPP-B phages. Metagenomic reads of both clades are highly abundant in surface waters, but the MPP-C phages show higher relative abundance in surface waters than in deeper waters, while MPP-B phages have higher relative abundance in deeper waters. Our study reveals that cyanophages with distinct genomic contents and infection kinetics can exhibit different depth profiles in the oceans.

The Coevolution of Fungal Mitochondrial Introns and Their Homing Endonucleases (GIY-YIG and LAGLIDADG)

Fungal mitochondrial (mt) genomes exhibit great diversity in size which is partially attributed to their variable intergenic regions and most importantly to the inclusion of introns within their genes. These introns belong to group I or II, and both of them are self-splicing. The majority of them carry genes encoding homing endonucleases, either LAGLIDADG or GIY-YIG. In this study, it was found that these intronic homing endonucleases genes (HEGs) may originate from mt free-standing open reading frames which can be found nowadays in species belonging to Early Diverging Fungi as “living fossils.” A total of 487 introns carrying HEGs which were located in the publicly available mt genomes of representative species belonging to orders from all fungal phyla was analyzed. Their distribution in the mt genes, their insertion target sequence, and the phylogenetic analyses of the HEGs showed that these introns along with their HEGs form a composite structure in which both selfish elements coevolved. The invasion of the ancestral free-standing HEGs in the introns occurred through a perpetual mechanism, called in this study as “aenaon” hypothesis. It is based on recombination, transpositions, and horizontal gene transfer events throughout evolution. HEGs phylogenetically clustered primarily according to their intron hosts and secondarily to the mt genes carrying the introns and their HEGs. The evolutionary models created revealed an “intron-early” evolution which was enriched by “intron-late” events through many different independent recombinational events which resulted from both vertical and horizontal gene transfers.

Recent and Ongoing Horizontal Transfer of Mitochondrial Introns Between Two Fungal Tree Pathogens

Two recently introduced fungal plant pathogens (Ceratocystis lukuohia and Ceratocystis huliohia) are responsible for Rapid ‘ōhi‘a Death (ROD) in Hawai‘i. Despite being sexually incompatible, the two pathogens often co-occur in diseased ‘ōhi‘a sapwood, where genetic interaction is possible. We sequenced and annotated 33 mitochondrial genomes of the two pathogens and related species, and investigated 35 total Ceratocystis mitogenomes. Ten mtDNA regions [one group I intron, seven group II introns, and two autonomous homing endonuclease (HE) genes] were heterogeneously present in C. lukuohia mitogenomes, which were otherwise identical. Molecular surveys with specific primers showed that the 10 regions had uneven geographic distribution amongst populations of C. lukuohia. Conversely, identical orthologs of each region were present in every studied isolate of C. huliohia regardless of geographical origin. Close relatives of C. lukuohia lacked or, rarely, had few and dissimilar orthologs of the 10 regions, whereas most relatives of C. huliohia had identical or nearly identical orthologs. Each region included or worked in tandem with HE genes or reverse transcriptase/maturases that could facilitate interspecific horizontal transfers from intron-minus to intron-plus alleles. These results suggest that the 10 regions originated in C. huliohia and are actively moving to populations of C. lukuohia, perhaps through transient cytoplasmic contact of hyphal tips (anastomosis) in the wound surface of ‘ōhi‘a trees. Such contact would allow for the transfer of mitochondria followed by mitochondrial fusion or cytoplasmic exchange of intron intermediaries, which suggests that further genomic interaction may also exist between the two pathogens.

The mitochondrial genome of Endoconidiophora resinifera is intron rich

Endoconidiophora resinifera (=Ceratocystis resinifera) is a blue-stain fungus that occurs on conifers. The data showed that the Endoconidiophora resinifera mitochondrial genome is one of the largest mitochondrial genomes (>220 kb) so far reported among members of the Ascomycota. An exceptional large number of introns (81) were noted and differences among the four strains were restricted to minor variations in intron numbers and a few indels and single nucleotide polymorphisms. The major differences among the four strains examined are due to size polymorphisms generated by the absence or presence of mitochondrial introns. Also, these mitochondrial genomes encode the largest cytochrome oxidase subunit 1 gene (47.5 kb) reported so far among the fungi. The large size for this gene again can be attributed to the large number of intron insertions. This study reports the first mitochondrial genome for the genus Endoconidiophora, previously members of this genus were assigned to Ceratocystis. The latter genus has recently undergone extensive taxonomic revisions and the mitochondrial genome might provide loci that could be applied as molecular markers assisting in the identification of taxa within this group of economically important fungi. The large mitochondrial genome also may provide some insight on mechanisms that can lead to mitochondrial genome expansion.

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.

Energy limitation of cyanophage development: implications for marine carbon cycling

Marine cyanobacteria are responsible for ~25% of the fixed carbon that enters the ocean biosphere. It is thought that abundant co-occurring viruses play an important role in regulating population dynamics of cyanobacteria and thus the cycling of carbon in the oceans. Despite this, little is known about how viral infections ‘play-out’ in the environment, particularly whether infections are resource or energy limited. Photoautotrophic organisms represent an ideal model to test this since available energy is modulated by the incoming light intensity through photophosphorylation. Therefore, we exploited phototrophy of the environmentally relevant marine cyanobacterium Synechococcus and monitored growth of a cyanobacterial virus (cyanophage). We found that light intensity has a marked effect on cyanophage infection dynamics, but that this is not manifest by a change in DNA synthesis. Instead, cyanophage development appears energy limited for the synthesis of proteins required during late infection. We posit that acquisition of auxiliary metabolic genes (AMGs) involved in light-dependent photosynthetic reactions acts to overcome this limitation. We show that cyanophages actively modulate expression of these AMGs in response to light intensity and provide evidence that such regulation may be facilitated by a novel mechanism involving light-dependent splicing of a group I intron in a photosynthetic AMG. Altogether, our data offers a mechanistic link between diurnal changes in irradiance and observed community level responses in metabolism, i.e., through an irradiance-dependent, viral-induced release of dissolved organic matter (DOM).

Comparative Functional Genomic Analysis of Two Vibrio Phages Reveals Complex Metabolic Interactions with the Host Cell

Sequencing and annotation was performed for two large double stranded DNA bacteriophages, φGrn1 and φSt2 of the Myoviridae family, considered to be of great interest for phage therapy against Vibrios in aquaculture live feeds. In addition, phage–host metabolic interactions and exploitation was studied by transcript profiling of selected viral and host genes. Comparative genomic analysis with other large Vibrio phages was also performed to establish the presence and location of homing endonucleases highlighting distinct features for both phages. Phylogenetic analysis revealed that they belong to the “schizoT4like” clade. Although many reports of newly sequenced viruses have provided a large set of information, basic research related to the shift of the bacterial metabolism during infection remains stagnant. The function of many viral protein products in the process of infection is still unknown. Genome annotation identified the presence of several viral open reading frames (ORFs) participating in metabolism, including a Sir2/cobB (sirtuin) protein and a number of genes involved in auxiliary NAD⁺ and nucleotide biosynthesis, necessary for phage DNA replication. Key genes were subsequently selected for detail study of their expression levels during infection. This work suggests a complex metabolic interaction and exploitation of the host metabolic pathways and biochemical processes, including a possible post-translational protein modification, by the virus during infection.

Programmed translational bypassing elements in mitochondria: structure, mobility, and evolutionary origin

Programmed translational bypassing enables ribosomes to 'ignore' a precise mRNA interval of several dozen nucleotides. Well-characterized bypassed sequences include hop and byp elements, present in bacteriophage T4 and mitochondria of the yeast Magnusiomyces capitatus, respectively. The bypassing mechanism of byps is probably similar to that of hop, yet the former appears more effective and less constrained as to sequence context. Furthermore, both elements are mobile but hop moves as part of a cassette including a homing endonuclease, whereas byps seem to spread like miniature DNA transposable elements known as GC clusters. Here, we argue that hop and byps arose independently by convergent evolution, and that byps evolved in magnusiomycete mitochondria due to (as yet unknown) alterations of the mitochondrial translation machinery.

Survey of chimeric IStron elements in bacterial genomes: multiple molecular symbioses between group I intron ribozymes and DNA transposons

IStrons are chimeric genetic elements composed of a group I intron associated with an insertion sequence (IS). The group I intron is a catalytic RNA providing the IStron with self-splicing ability, which renders IStron insertions harmless to the host genome. The IS element is a DNA transposon conferring mobility, and thus allowing the IStron to spread in genomes. IStrons are therefore a striking example of a molecular symbiosis between unrelated genetic elements endowed with different functions. In this study, we have conducted the first comprehensive survey of IStrons in sequenced genomes that provides insights into the distribution, diversity, origin and evolution of IStrons. We show that IStrons have a restricted phylogenetic distribution limited to two bacterial phyla, the Firmicutes and the Fusobacteria. Nevertheless, diverse IStrons representing two major groups targeting different insertion site motifs were identified. This taken with the finding that while the intron components of all IStrons belong to the same structural class, they are fused to different IS families, indicates that multiple intron–IS symbioses have occurred during evolution. In addition, introns and IS elements related to those that were at the origin of IStrons were also identified.

Mapping homing endonuclease cleavage sites using in vitro generated protein

Mapping the precise position of endonucleolytic cleavage sites is a fundamental experimental technique used to describe the function of a homing endonuclease. However, these proteins are often recalcitrant to cloning and over-expression in biological systems because of toxicity induced by spurious DNA cleavage events. In this chapter we outline the steps to successfully express a homing endonuclease in vitro and use this product in nucleotide-resolution cleavage assays.

Horizontal transfer and gene conversion as an important driving force in shaping the landscape of mitochondrial introns

Group I introns are highly dynamic and mobile, featuring extensive presence-absence variation and widespread horizontal transfer. Group I introns can invade intron-lacking alleles via intron homing powered by their own encoded homing endonuclease gene (HEG) after horizontal transfer or via reverse splicing through an RNA intermediate. After successful invasion, the intron and HEG are subject to degeneration and sequential loss. It remains unclear whether these mechanisms can fully address the high dynamics and mobility of group I introns. Here, we found that HEGs undergo a fast gain-and-loss turnover comparable with introns in the yeast mitochondrial 21S-rRNA gene, which is unexpected, as the intron and HEG are generally believed to move together as a unit. We further observed extensively mosaic sequences in both the introns and HEGs, and evidence of gene conversion between HEG-containing and HEG-lacking introns. Our findings suggest horizontal transfer and gene conversion can accelerate HEG/intron degeneration and loss, or rescue and propagate HEG/introns, and ultimately result in high HEG/intron turnover rate. Given that up to 25% of the yeast mitochondrial genome is composed of introns and most mitochondrial introns are group I introns, horizontal transfer and gene conversion could have served as an important mechanism in introducing mitochondrial intron diversity, promoting intron mobility and consequently shaping mitochondrial genome architecture.

Bacterial group I introns: mobile RNA catalysts

Group I introns are intervening sequences that have invaded tRNA, rRNA and protein coding genes in bacteria and their phages. The ability of group I introns to self-splice from their host transcripts, by acting as ribozymes, potentially renders their insertion into genes phenotypically neutral. Some group I introns are mobile genetic elements due to encoded homing endonuclease genes that function in DNA-based mobility pathways to promote spread to intronless alleles. Group I introns have a limited distribution among bacteria and the current assumption is that they are benign selfish elements, although some introns and homing endonucleases are a source of genetic novelty as they have been co-opted by host genomes to provide regulatory functions. Questions regarding the origin and maintenance of group I introns among the bacteria and phages are also addressed.

Riboregulators and the role of Hfq in photosynthetic bacteria

Anoxygenic and oxygenic bacteria directly convert solar energy into biomass using photosynthesis. The formation and composition of photosynthetic complexes has to be tightly controlled in response to environmental conditions, as exposure to sunlight can be harmful due to the generation of reactive oxygen species and the damaging effects of UV irradiation. Therefore, photosynthetic bacteria are exposed to a particular set of regulatory challenges in addition to those that also affect other bacteria, requiring sophisticated regulatory systems. Indeed, hundreds of potential regulatory RNAs have been identified in photosynthetic model bacteria as well as antisense RNAs (asRNAs) of up to several kb in length that protect certain mRNAs from degradation. The trans-acting small non-coding RNAs (sRNAs), PcrZ and PsrR1, control pigment and photosystem biogenesis in Rhodobacter sphaeroides and cyanobacteria, respectively. The asRNAs IsrR and As1_flv4 act as negative regulators and the asRNAs PsbA2R and PsbA3R as positive effectors of photosynthesis gene expression in Synechocystis 6803.

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.

The global virome: not as big as we thought?

Viruses likely infect all organisms, serving to unknown extent as genetic vectors in complex networks of organisms. Environmental virologists have revealed that these abundant nanoscale entities are global players with critical roles in every ecosystem investigated. Curiously, novel genes dominate viral genomes and metagenomes, which has led to the suggestion that viruses represent the largest reservoir of unexplored genetic material on Earth with literature estimates, extrapolating from 14 mycobacteriophage genomes, suggesting that two billion phage-encoded ORFs remain to be discovered. Here we examine (meta)genomic data available in the decade since this provocative assertion, and use 'protein clusters' to evaluate whether sampling technologies have advanced to the point that we may be able to sample 'all' of viral diversity in nature.

I-PfoP3I: a novel nicking HNH homing endonuclease encoded in the group I intron of the DNA polymerase gene in Phormidium foveolarum phage Pf-WMP3

Homing endonucleases encoded in a group I self-splicing intron in a protein-coding gene in cyanophage genomes have not been reported, apart from some free-standing homing edonucleases. In this study, a nicking DNA endonuclease, I-PfoP3I, encoded in a group IA2 intron in the DNA polymerase gene of a T7-like cyanophage Pf-WMP3, which infects the freshwater cyanobacterium Phormidium foveolarum is described. The Pf-WMP3 intron splices efficiently in vivo and self-splices in vitro simultaneously during transcription. I-PfoP3I belongs to the HNH family with an unconventional C-terminal HNH motif. I-PfoP3I nicks the intron-minus Pf-WMP3 DNA polymerase gene more efficiently than the Pf-WMP4 DNA polymerase gene that lacks any intervening sequence in vitro, indicating the variable capacity of I-PfoP3I. I-PfoP3I cleaves 4 nt upstream of the intron insertion site on the coding strand of EXON 1 on both intron-minus Pf-WMP3 and Pf-WMP4 DNA polymerase genes. Using an in vitro cleavage assay and scanning deletion mutants of the intronless target site, the minimal recognition site was determined to be a 14 bp region downstream of the cut site. I-PfoP3I requires Mg²⁺, Ca²⁺ or Mn²⁺ for nicking activity. Phylogenetic analysis suggests that the intron and homing endonuclease gene elements might be inserted in Pf-WMP3 genome individually after differentiation from Pf-WMP4. To our knowledge, this is the first report of the presence of a group I self-splicing intron encoding a functional homing endonuclease in a protein-coding gene in a cyanophage genome.

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.

Home and away- the evolutionary dynamics of homing endonucleases

Background

Homing endonucleases (HEases) are a large and diverse group of site-specific DNAases. They reside within self-splicing introns and inteins, and promote their horizontal dissemination. In recent years, HEases have been the focus of extensive research due to their promising potential use in gene targeting procedures for the treatment of genetic diseases and for the genetic engineering of crop, animal models and cell lines.

Results

Using mathematical analysis and computational modeling, we present here a novel account for the evolution and population dynamics of HEase genes (HEGs). We describe HEGs as paradoxical selfish elements whose long-term persistence in a single population relies on low transmission rates and a positive correlation between transmission efficiency and toxicity.

Conclusion

Plausible conditions allow HEGs to sustain at high frequency through long evolutionary periods, with the endonuclease frequency being either at equilibrium or periodically oscillating. The predictions of our model may prove important not only for evolutionary theory but also for gene therapy and bio-engineering applications of HEases.

A homing endonuclease and the 50-nt ribosomal bypass sequence of phage T4 constitute a mobile DNA cassette

Since its initial description more than two decades ago, the ribosome bypass (or "hop") sequence of phage T4 stands out as a uniquely extreme example of programmed translational frameshifting. The gene for a DNA topoisomerase subunit of T4 has been split by a 1-kb insertion into two genes that retain topoisomerase function. A second 50-nt insertion, beginning with an in-phase stop codon, is inserted near the start of the newly created downstream gene 60. Instead of terminating at this stop codon, approximately half of the ribosomes skip 50 nucleotides and continue translation in a new reading frame. However, no functions, regulatory or otherwise, have been imputed for the truncated peptide that results from termination at codon 46 or for the bypass sequence itself. Moreover, how this unusual mRNA organization arose and why it is maintained have never been explained. We show here that a homing endonuclease (MobA) is encoded in the insertion that created gene 60, and the mobA gene together with the bypass sequence constitute a mobile DNA cassette. The bypass sequence provides protection against self-cleavage by the nuclease, whereas the nuclease promotes horizontal spread of the entire cassette to related bacteriophages. Group I introns frequently provide protection against self-cleavage by associated homing endonucleases. We present a scenario by which the bypass sequence, which is otherwise a unique genetic element, might have been derived from a degenerate group I intron.

cis-antisense RNA, another level of gene regulation in bacteria

A substantial amount of antisense transcription is a hallmark of gene expression in eukaryotes. However, antisense transcription was first demonstrated in bacteria almost 50 years ago. The transcriptomes of bacteria as different as Helicobacter pylori, Bacillus subtilis, Escherichia coli, Synechocystis sp. strain PCC6803, Mycoplasma pneumoniae, Sinorhizobium meliloti, Geobacter sulfurreducens, Vibrio cholerae, Chlamydia trachomatis, Pseudomonas syringae, and Staphylococcus aureus have now been reported to contain antisense RNA (asRNA) transcripts for a high percentage of genes. Bacterial asRNAs share functional similarities with trans-acting regulatory RNAs, but in addition, they use their own distinct mechanisms. Among their confirmed functional roles are transcription termination, codegradation, control of translation, transcriptional interference, and enhanced stability of their respective target transcripts. Here, we review recent publications indicating that asRNAs occur as frequently in simple unicellular bacteria as they do in higher organisms, and we provide a comprehensive overview of the experimentally confirmed characteristics of asRNA actions and intimately linked quantitative aspects. Emerging functional data suggest that asRNAs in bacteria mediate a plethora of effects and are involved in far more processes than were previously anticipated. Thus, the functional impact of asRNAs should be considered when developing new strategies against pathogenic bacteria and when optimizing bacterial strains for biotechnology.

Recurrent insertion of 5'-terminal nucleotides and loss of the branchpoint motif in lineages of group II introns inserted in mitochondrial preribosomal RNAs

A survey of sequence databases revealed 10 instances of subgroup IIB1 mitochondrial ribosomal introns with 1 to 33 additional nucleotides inserted between the 5' exon and the consensus sequence at the intron 5' end. These 10 introns depart further from the IIB1 consensus in their predicted domain VI structure: In contrast to its basal helix and distal GNRA terminal loop, the middle part of domain VI is highly variable and lacks the bulging A that serves as the branchpoint in lariat formation. In vitro experiments using two closely related IIB1 members inserted at the same ribosomal RNA site in the basidiomycete fungi Grifola frondosa and Pycnoporellus fulgens revealed that both ribozymes are capable of efficient self-splicing. However, whereas the Grifola intron was excised predominantly as a lariat, the Pycnoporellus intron, which possesses six additional nucleotides at the 5' end, yielded only linear products, consistent with its predicted domain VI structure. Strikingly, all of the introns with 5' terminal insertions lack the EBS2 exon-binding site. Moreover, several of them are part of the small subset of group II introns that encode potentially functional homing endonucleases of the LAGLIDADG family rather than reverse transcriptases. Such coincidences suggest causal relationships between the shift to DNA-based mobility, the loss of one of the two ribozyme sites for binding the 5' exon, and the exclusive use of hydrolysis to initiate splicing.

Genomic analysis of oceanic cyanobacterial myoviruses compared with T4-like myoviruses from diverse hosts and environments

T4-like myoviruses are ubiquitous, and their genes are among the most abundant documented in ocean systems. Here we compare 26 T4-like genomes, including 10 from non-cyanobacterial myoviruses, and 16 from marine cyanobacterial myoviruses (cyanophages) isolated on diverse Prochlorococcus or Synechococcus hosts. A core genome of 38 virion construction and DNA replication genes was observed in all 26 genomes, with 32 and 25 additional genes shared among the non-cyanophage and cyanophage subsets, respectively. These hierarchical cores are highly syntenic across the genomes, and sampled to saturation. The 25 cyanophage core genes include six previously described genes with putative functions (psbA, mazG, phoH, hsp20, hli03, cobS), a hypothetical protein with a potential phytanoyl-CoA dioxygenase domain, two virion structural genes, and 16 hypothetical genes. Beyond previously described cyanophage-encoded photosynthesis and phosphate stress genes, we observed core genes that may play a role in nitrogen metabolism during infection through modulation of 2-oxoglutarate. Patterns among non-core genes that may drive niche diversification revealed that phosphorus-related gene content reflects source waters rather than host strain used for isolation, and that carbon metabolism genes appear associated with putative mobile elements. As well, phages isolated on Synechococcus had higher genome-wide %G+C and often contained different gene subsets (e.g. petE, zwf, gnd, prnA, cpeT) than those isolated on Prochlorococcus. However, no clear diagnostic genes emerged to distinguish these phage groups, suggesting blurred boundaries possibly due to cross-infection. Finally, genome-wide comparisons of both diverse and closely related, co-isolated genomes provide a locus-to-locus variability metric that will prove valuable for interpreting metagenomic data sets.

Evolutionary dynamics of the mS952 intron: a novel mitochondrial group II intron encoding a LAGLIDADG homing endonuclease gene

Examination of the mitochondrial small subunit ribosomal RNA (rns) gene of five species of the fungal genus Leptographium revealed that the gene has been invaded at least once at position 952 by a group II intron encoding a LAGLIDADG homing endonuclease gene. Phylogenetic analyses of the intron and homing endonuclease sequences indicated that each element in Leptographium species forms a single clade and is closely related to the group II intron/homing endonuclease gene composite element previously reported at position 952 of the mitochondrial rns gene of Cordyceps species and of Cryphonectria parasitica. The results of an intron survey of the mt rns gene of Leptographium species superimposed onto the phylogenetic analysis of the host organisms suggest that the composite element was transmitted vertically in Leptographium lundbergii. However, its stochastic distribution among strains of L. wingfieldii, L. terebrantis, and L. truncatum suggests that it has been horizontally transmitted by lateral gene transfer among these species, although the random presence of the intron may reflect multiple random loss events. A model is proposed describing the initial invasion of the group II intron in the rns gene of L. lundbergii by a LAGLIDADG homing endonuclease gene and subsequent evolution of this gene to recognize a novel DNA target site, which may now promote the mobility of the intron and homing endonuclease gene as a composite element.

Evolution of introns in the archaeal world

The self-splicing group I introns are removed by an autocatalytic mechanism that involves a series of transesterification reactions. They require RNA binding proteins to act as chaperones to correctly fold the RNA into an active intermediate structure in vivo. Pre-tRNA introns in Bacteria and in higher eukaryote plastids are typical examples of self-splicing group I introns. By contrast, two striking features characterize RNA splicing in the archaeal world. First, self-splicing group I introns cannot be found, to this date, in that kingdom. Second, the RNA splicing scenario in Archaea is uniform: All introns, whether in pre-tRNA or elsewhere, are removed by tRNA splicing endonucleases. We suggest that in Archaea, the protein recruited for splicing is the preexisting tRNA splicing endonuclease and that this enzyme, together with the ligase, takes over the task of intron removal in a more efficient fashion than the ribozyme. The extinction of group I introns in Archaea would then be a consequence of recruitment of the tRNA splicing endonuclease. We deal here with comparative genome analysis, focusing specifically on the integration of introns into genes coding for 23S rRNA molecules, and how this newly acquired intron has to be removed to regenerate a functional RNA molecule. We show that all known oligomeric structures of the endonuclease can recognize and cleave a ribosomal intron, even when the endonuclease derives from a strain lacking rRNA introns. The persistence of group I introns in mitochondria and chloroplasts would be explained by the inaccessibility of these introns to the endonuclease.

Regulatory RNAs in cyanobacteria: developmental decisions, stress responses and a plethora of chromosomally encoded cis-antisense RNAs

Cyanobacteria are the only prokaryotes which directly convert solar energy into biomass using oxygenic photosynthesis. Therefore, these bacteria are of interest for the production of biofuels, biotechnology and are of tremendous relevance for primary carbon fixation in many ecosystems. Mechanisms controlling gene expression cannot be understood entirely without information on the numbers and functions of regulatory RNAs. In cyanobacteria, non-coding RNAs have been characterized from simple unicellular species such as Prochlorococcus up to complex species such as Anabaena. Several of these RNAs function in the control of stress responses, photosynthesis, outer cell membrane protein biosynthesis and the differentiation of cells.

Homing endonucleases residing within inteins: evolutionary puzzles awaiting genetic solutions

Inteins are selfish genetic elements that disrupt the sequence of protein-coding genes and are excised post-translationally. Most inteins also contain a HEN (homing endonuclease) domain, which is important for their horizontal transmission. The present review focuses on the evolution of inteins and their nested HENs, and highlights several unsolved questions that could benefit from molecular genetic approaches. Such approaches can be well carried out in halophilic archaea, which are naturally intein-rich and have highly developed genetic tools for their study. In particular, the fitness effects of habouring an intein/HEN can be tested in direct competition assays, providing additional insights that will improve current evolutionary models.

Horizontal gene transfers with or without cell fusions in all categories of the living matter

This article reviews the history of widespread exchanges of genetic segments initiated over 3 billion years ago, to be part of their life style, by sphero-protoplastic cells, the ancestors of archaea, prokaryota, and eukaryota. These primordial cells shared a hostile anaerobic and overheated environment and competed for survival. "Coexist with, or subdue and conquer, expropriate its most useful possessions, or symbiose with it, your competitor" remain cellular life's basic rules. This author emphasizes the role of viruses, both in mediating cell fusions, such as the formation of the first eukaryotic cell(s) from a united crenarchaeon and prokaryota, and the transfer of host cell genes integrated into viral (phages) genomes. After rising above the Darwinian threshold, rigid rules of speciation and vertical inheritance in the three domains of life were established, but horizontal gene transfers with or without cell fusions were never abolished. The author proves with extensive, yet highly selective documentation, that not only unicellular microorganisms, but the most complex multicellular entities of the highest ranks resort to, and practice, cell fusions, and donate and accept horizontally (laterally) transferred genes. Cell fusions and horizontally exchanged genetic materials remain the fundamental attributes and inherent characteristics of the living matter, whether occurring accidentally or sought after intentionally. These events occur to cells stagnating for some 3 milliard years at a lower yet amazingly sophisticated level of evolution, and to cells achieving the highest degree of differentiation, and thus functioning in dependence on the support of a most advanced multicellular host, like those of the human brain. No living cell is completely exempt from gene drains or gene insertions.

Genomes of the T4-related bacteriophages as windows on microbial genome evolution

The T4-related bacteriophages are a group of bacterial viruses that share morphological similarities and genetic homologies with the well-studied Escherichia coli phage T4, but that diverge from T4 and each other by a number of genetically determined characteristics including the bacterial hosts they infect, the sizes of their linear double-stranded (ds) DNA genomes and the predicted compositions of their proteomes. The genomes of about 40 of these phages have been sequenced and annotated over the last several years and are compared here in the context of the factors that have determined their diversity and the diversity of other microbial genomes in evolution. The genomes of the T4 relatives analyzed so far range in size between ~160,000 and ~250,000 base pairs (bp) and are mosaics of one another, consisting of clusters of homology between them that are interspersed with segments that vary considerably in genetic composition between the different phage lineages. Based on the known biological and biochemical properties of phage T4 and the proteins encoded by the T4 genome, the T4 relatives reviewed here are predicted to share a genetic core, or "Core Genome" that determines the structural design of their dsDNA chromosomes, their distinctive morphology and the process of their assembly into infectious agents (phage morphogenesis). The Core Genome appears to be the most ancient genetic component of this phage group and constitutes a mere 12-15% of the total protein encoding potential of the typical T4-related phage genome. The high degree of genetic heterogeneity that exists outside of this shared core suggests that horizontal DNA transfer involving many genetic sources has played a major role in diversification of the T4-related phages and their spread to a wide spectrum of bacterial species domains in evolution. We discuss some of the factors and pathways that might have shaped the evolution of these phages and point out several parallels between their diversity and the diversity generally observed within all groups of interrelated dsDNA microbial genomes in nature.

A group II intron encodes a functional LAGLIDADG homing endonuclease and self-splices under moderate temperature and ionic conditions

A group II intron encoding a protein belonging to the LAGLIDADG family of homing endonucleases was identified in the mitochondrial rns gene of the filamentous fungus Leptographium truncatum, and the catalytic activities of both the intron and its encoded protein were characterized. A model of the RNA secondary structure indicates that the intron is a member of the IIB1 subclass and the open reading frame is inserted in ribozyme domain III. In vitro assays carried out with two versions of the intron, one in which the open reading frame was removed and the other in which it was present, demonstrate that both versions of the intron readily self-splice at 37 degrees C and at a concentration of MgCl(2) as low as 6 mM. The open reading frame encodes a functional LAGLIDADG homing endonuclease that cleaves 2 (top strand) and 6 (bottom strand) nucleotides (nt) upstream of the intron insertion site, generating 4 nt 3' OH overhangs. In vitro splicing assays carried out in the absence and presence of the intron-encoded protein indicate that the protein does not enhance intron splicing, and RNA-binding assays show that the protein does not appear to bind to the intron RNA precursor transcript. These findings raise intriguing questions concerning the functional and evolutionary relationships of the two components of this unique composite element.

An antisense RNA in a lytic cyanophage links psbA to a gene encoding a homing endonuclease

Cyanophage genomes frequently possess the psbA gene, encoding the D1 polypeptide of photosystem II. This protein is believed to maintain host photosynthetic capacity during infection and enhance phage fitness under high-light conditions. Although the first documented cyanophage-encoded psbA gene contained a group I intron, this feature has not been widely reported since, despite a plethora of new sequences becoming available. In this study, we show that in cyanophage S-PM2, this intron is spliced during the entire infection cycle. Furthermore, we report the widespread occurrence of psbA introns in marine metagenomic libraries, and with psbA often adjacent to a homing endonuclease (HE). Bioinformatic analysis of the intergenic region between psbA and the adjacent HE gene F-CphI in S-PM2 showed the presence of an antisense RNA (asRNA) connecting these two separate genetic elements. The asRNA is co-regulated with psbA and F-CphI, suggesting its involvement with their expression. Analysis of scaffolds from global ocean survey datasets shows this asRNA to be commonly associated with the 3' end of cyanophage psbA genes, implying that this potential mechanism of regulating marine 'viral' photosynthesis is evolutionarily conserved. Although antisense transcription is commonly found in eukaryotic and increasingly also in prokaryotic organisms, there has been no indication for asRNAs in lytic phages so far. We propose that this asRNA also provides a means of preventing the formation of mobile group I introns within cyanophage psbA genes.

Phage T4 mobE promotes trans homing of the defunct homing endonuclease I-TevIII

Homing endonucleases are site-specific DNA endonucleases that typically function as mobile genetic elements by introducing a double-strand break (DSB) in genomes that lack the endonuclease, resulting in a unidirectional gene conversion event that mobilizes the homing endonuclease gene and flanking DNA. Here, we characterize phage T4-encoded mobE, a predicted free-standing HNH family homing endonuclease. We show that mobE is promoterless and dependent on upstream transcription for expression, and that an internal intrinsic terminator regulates mobE transcript levels. Crucially, in vivo mapping experiments revealed a MobE-dependent, strand-specific nick in the non-coding strand of the nrdB gene of phage T2. An internal deletion of the predicted HNH catalytic motif of MobE abolishes nicking, and reduces high-frequency inheritance of mobE. Sequence polymorphisms of progeny phage that inherit mobE are consistent with DSB repair pathways. Significantly, we found that mobility of the neighboring I-TevIII, a defunct homing endonuclease encoded within a group I intron interrupting the nrdB gene of phage T4, was dependent on an intact mobE gene. Thus, our data indicate that the stagnant nrdB intron and I-TevIII are mobilized in trans as a consequence of a MobE-dependent gene conversion event, facilitating persistence of genetic elements that have no inherent means of promoting their own mobility.

Conservation of intron and intein insertion sites: implications for life histories of parasitic genetic elements

BACKGROUND

Inteins and introns are genetic elements that are removed from proteins and RNA after translation or transcription, respectively. Previous studies have suggested that these genetic elements are found in conserved parts of the host protein. To our knowledge this type of analysis has not been done for group II introns residing within a gene. Here we provide quantitative statistical support from an analyses of proteins that host inteins, group I introns, group II introns and spliceosomal introns across all three domains of life.

RESULTS

To determine whether or not inteins, group I, group II, and spliceosomal introns are found preferentially in conserved regions of their respective host protein, conservation profiles were generated and intein and intron positions were mapped to the profiles. Fisher's combined probability test was used to determine the significance of the distribution of insertion sites across the conservation profile for each protein. For a subset of studied proteins, the conservation profile and insertion positions were mapped to protein structures to determine if the insertion sites correlate to regions of functional activity. All inteins and most group I introns were found to be preferentially located within conserved regions; in contrast, a bacterial intein-like protein, group II and spliceosomal introns did not show a preference for conserved sites.

CONCLUSIONS

These findings demonstrate that inteins and group I introns are found preferentially in conserved regions of their respective host proteins. Homing endonucleases are often located within inteins and group I introns and these may facilitate mobility to conserved regions. Insertion at these conserved positions decreases the chance of elimination, and slows deletion of the elements, since removal of the elements has to be precise as not to disrupt the function of the protein. Furthermore, functional constrains on the targeted site make it more difficult for hosts to evolve immunity to the homing endonuclease. Therefore, these elements will better survive and propagate as molecular parasites in conserved sites. In contrast, spliceosomal introns and group II introns do not show significant preference for conserved sites and appear to have adopted a different strategy to evade loss.

Genetic insertions and diversification of the PolB-type DNA polymerase (gp43) of T4-related phages

In Escherichia coli phage T4 and many of its phylogenetic relatives, gene 43 consists of a single cistron that encodes a PolB family (PolB-type) DNA polymerase. We describe the divergence of this phage gene and its protein product (gp43) (gene product 43) among 26 phylogenetic relatives of T4 and discuss our observations in the context of diversity among the widely distributed PolB enzymes in nature. In two T4 relatives that grow in Aeromonas salmonicida phages 44RR and 25, gene 43 is fragmented by different combinations of three distinct types of DNA insertion elements: (a) a short intercistronic untranslated sequence (IC-UTS) that splits the polymerase gene into two cistrons, 43A and 43B, corresponding to N-terminal (gp43A) and C-terminal (gp43B) protein products; (b) a freestanding homing endonuclease gene (HEG) inserted between the IC-UTS and the 43B cistron; and (c) a group I intron in the 43B cistron. Phage 25 has all three elements, whereas phage 44RR has only the IC-UTS. We present evidence that (a) the split gene of phage 44RR encodes a split DNA polymerase consisting of a complex between gp43A and gp43B subunits; (b) the putative HEG encodes a double-stranded DNA endonuclease that specifically cleaves intron-free homologues of the intron-bearing 43B site; and (c) the group I intron is a self-splicing RNA. Our results suggest that some freestanding HEGs can mediate the homing of introns that do not encode their own homing enzymes. The results also suggest that different insertion elements can converge on a polB gene and evolve into a single integrated system for lateral transfer of polB genetic material. We discuss the possible pathways for the importation of such insertion elements into the genomes of T4-related phages.

A likely pathway for formation of mobile group I introns

Mobile group I introns are RNA splicing elements that have been invaded by endonuclease genes. These endonucleases facilitate intron mobility by a unidirectional, duplicative gene-conversion process known as homing [1]. Survival of the invading endonuclease depends upon its ability to promote intron mobility. Therefore, the endonuclease must either quickly change its cleavage specificity to match the site of intron insertion, or it must already be preadapted to cleave this sequence. Here we show that the group I intron in the DNA polymerase gene of T7-like bacteriophage PhiI is mobile, dependent upon its intronic HNH homing endonuclease gene, I-TslI. We also show that gene 5.3 of phage T3, located adjacent to its intronless DNA polymerase gene, is a homologous homing endonuclease gene whose protein product initiates efficient spread of gene 5.3 into empty sites in related phages. Both of these endonucleases cleave intronless DNA polymerase genes at identical positions. This shared feature between an intronic and free-standing endonuclease is unprecedented. Based on this evidence, we propose that introns and their homing endonucleases evolve separately to target the same highly conserved sequences, uniting afterwards to create a composite mobile element.

Selfish DNA: homing endonucleases find a home

Self-splicing group I introns come in two flavours - those with a homing endonuclease to promote mobility of the intron, and those without an endonuclease. How homing endonucleases and self-splicing introns associate to form a composite selfish genetic element is a question of long-standing interest. Recent work has revealed that a shared characteristic of both introns and endonucleases, the targeting of conserved sequences, may provide the impetus for the evolution of composite mobile genetic elements.