Invasion of a multitude of genetic niches by mobile endonuclease genes

A Photoresponsive Homing Endonuclease for Programmed DNA Cleavage

Homing endonucleases are used in a wide range of biotechnological applications including gene editing, in gene drive systems, and for the modification of DNA structures, arrays, and prodrugs. However, controlling nuclease activity and sequence specificity remain key challenges when developing new tools. Here a photoresponsive homing endonuclease was engineered for optical control of DNA cleavage by partitioning DNA binding and nuclease domains of the monomeric homing endonuclease I-TevI into independent polypeptide chains. Use of the Aureochrome1a light-oxygen-voltage domain delivered control of dimerization with light. Illumination reduced the concentration needed to achieve 50% cleavage of the homing target site by 6-fold when compared to the dark state, resulting in an up to 9-fold difference in final yields between cleavage products. I-TevI nucleases with and without a native I-TevI zinc finger motif displayed different nuclease activity and sequence preference impacting the promiscuity of the nuclease domain. By harnessing an alternative DNA binding domain, target preference was reprogrammed only when the nuclease lacked the I-TevI zinc finger motif. This work establishes a first-generation photoresponsive platform for spatiotemporal activation of DNA cleavage.

Synthetic gene drives as an anthropogenic evolutionary force

Genetic drive represents a fundamental evolutionary force that can exact profound change to the genetic composition of populations by biasing allele transmission. Herein I propose that the use of synthetic homing gene drives, the human-mediated analog of endogenous genetic drives, warrants the designation of 'genetic welding' as an anthropogenic evolutionary force. Conceptually, this distinction parallels that of artificial and natural selection. Genetic welding is capable of imposing complex and rapid heritable phenotypic change on entire populations, whether motivated by biodiversity conservation or public health. Unanticipated possible long-term evolutionary outcomes, however, demand further investigation and bioethical consideration. The emerging importance of genetic welding also compels our explicit recognition of genetic drive as an addition to the other four fundamental forces of evolution.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

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.

Comparative mitogenomics of Agaricomycetes: Diversity, abundance, impact and coding potential of putative open-reading frames

The mitochondrion is an organelle found in eukaryote organisms, and it is vital for different cellular pathways. The mitochondrion has its own DNA molecule and, because its genetic content is relatively conserved, despite the variation of size and structure, mitogenome sequences have been widely used as a promising molecular biomarker for taxonomy and evolution in fungi. In this study, the mitogenomes of two fungal species of Agaricomycetes class, Phellinotus piptadeniae and Trametes villosa, were assembled and annotated for the first time. We used these newly sequenced mitogenomes for comparative analyses with other 55 mitogenomes of Agaricomycetes available in public databases. Mitochondrial DNA (mtDNA) size and content are highly variable and non-coding and intronic regions, homing endonucleases (HEGs), and unidentified ORFs (uORFs) significantly contribute to the total size of the mitogenome. Furthermore, accessory genes (most of them as HEGs) are shared between distantly related species, most likely as a consequence of horizontal gene transfer events. Conversely, uORFs are only shared between taxonomically related species, most probably as a result of vertical evolutionary inheritance. Additionally, codon usage varies among mitogenomes and the GC content of mitochondrial features may be used to distinguish coding from non-coding sequences. Our results also indicated that transposition events of mitochondrial genes to the nuclear genome are not common. Despite the variation of size and content of the mitogenomes, mitochondrial genes seemed to be reliable molecular markers in our time-divergence analysis, even though the nucleotide substitution rates of mitochondrial and nuclear genomes of fungi are quite different. We also showed that many events of mitochondrial gene shuffling probably happened amongst the Agaricomycetes during evolution, which created differences in the gene order among species, even those of the same genus. Altogether, our study revealed new information regarding evolutionary dynamics in Agaricomycetes.

Engineering the Composition and Fate of Wild Populations with Gene Drive

Insects play important roles as predators, prey, pollinators, recyclers, hosts, parasitoids, and sources of economically important products. They can also destroy crops; wound animals; and serve as vectors for plant, animal, and human diseases. Gene drive-a process by which genes, gene complexes, or chromosomes encoding specific traits are made to spread through wild populations, even if these traits result in a fitness cost to carriers-provides new opportunities for altering populations to benefit humanity and the environment in ways that are species specific and sustainable. Gene drive can be used to alter the genetic composition of an existing population, referred to as population modification or replacement, or to bring about population suppression or elimination. We describe technologies under consideration, progress that has been made, and remaining technological hurdles, particularly with respect to evolutionary stability and our ability to control the spread and ultimate fate of genes introduced into populations.

The yeast mating-type switching endonuclease HO is a domesticated member of an unorthodox homing genetic element family

The mating-type switching endonuclease HO plays a central role in the natural life cycle of Saccharomyces cerevisiae, but its evolutionary origin is unknown. HO is a recent addition to yeast genomes, present in only a few genera close to Saccharomyces. Here we show that HO is structurally and phylogenetically related to a family of unorthodox homing genetic elements found in Torulaspora and Lachancea yeasts. These WHO elements home into the aldolase gene FBA1, replacing its 3' end each time they integrate. They resemble inteins but they operate by a different mechanism that does not require protein splicing. We show that a WHO protein cleaves Torulaspora delbrueckii FBA1 efficiently and in an allele-specific manner, leading to DNA repair by gene conversion or NHEJ. The DNA rearrangement steps during WHO element homing are very similar to those during mating-type switching, and indicate that HO is a domesticated WHO-like element.

In the same way as a sperm from a male and an egg from a female join together to form an embryo in most animals, yeast cells have two sexes that coordinate how they reproduce. These are called “mating types” and, rather than male or female, an individual yeast cell can either be mating type “a” or “alpha”. Every yeast cell contains the genes for both mating types, and each cell’s mating type is determined by which of those genes it has active.

Only one mating type gene can be ‘on’ at a time, but some yeast species can swap mating type on demand by switching the corresponding genes ‘on’ or ‘off’. This switch is unusual. Rather than simply activate one of the genes it already has, the yeast cell keeps an inactive version of each mating type gene tucked away, makes a copy of the gene it wants to be active and pastes that copy into a different location in its genome. To do all of this yeast need another gene called HO. This gene codes for an enzyme that cuts the DNA at the location of the active mating type gene. This makes an opening that allows the cell to replace the ‘a’ gene with the ‘alpha’ gene, or vice versa. This system allows yeast cells to continue mating even if all the cells in a colony start off as the same mating type. But, cutting into the DNA is risky, and can damage the health of the cell. So, why did yeast cells evolve a system that could cause them harm?

To find out where the HO gene came from, Coughlan et al. searched through all the available genomes from yeast species for other genes with similar sequences and identified a cluster which they nicknamed “weird HO” genes, or WHO genes for short. Testing these genes revealed that they also code for enzymes that make cuts in the yeast genome, but the way the cell repairs the cuts is different. The WHO genes are jumping genes. When the enzyme encoded by a WHO gene makes a cut in the genome, the yeast cell copies the gene into the gap, allowing the gene to ‘jump’ from one part of the genome to another. It is possible that this was the starting point for the evolution of the HO gene. Changes to a WHO gene could have allowed it to cut into the mating type region of the yeast genome, giving the yeast an opportunity to ‘domesticate’ it. Over time, the yeast cell stopped the WHO gene from jumping into the gap and started using the cut to change its mating type.

Understanding how cells adapt genes for different purposes is a key question in evolutionary biology. There are many other examples of domesticated jumping genes in other organisms, including in the human immune system. Understanding the evolution of HO not only sheds light on how yeast mating type switching evolved, but on how other species might harness and adapt their genes.

An Evolutionary Perspective on Yeast Mating-Type Switching

Cell differentiation in yeast species is controlled by a reversible, programmed DNA-rearrangement process called mating-type switching. Switching is achieved by two functionally similar but structurally distinct processes in the budding yeast Saccharomyces cerevisiae and the fission yeast Schizosaccharomyces pombe. In both species, haploid cells possess one active and two silent copies of the mating-type locus (a three-cassette structure), the active locus is cleaved, and synthesis-dependent strand annealing is used to replace it with a copy of a silent locus encoding the opposite mating-type information. Each species has its own set of components responsible for regulating these processes. In this review, we summarize knowledge about the function and evolution of mating-type switching components in these species, including mechanisms of heterochromatin formation, MAT locus cleavage, donor bias, lineage tracking, and environmental regulation of switching. We compare switching in these well-studied species to others such as Kluyveromyces lactis and the methylotrophic yeasts Ogataea polymorpha and Komagataella phaffii. We focus on some key questions: Which cells switch mating type? What molecular apparatus is required for switching? Where did it come from? And what is the evolutionary purpose of switching?

Three new active members of the I-OnuI family of homing endonucleases

In vitro characterization of 3 LAGLIDADG-type homing endonucleases (HEs) (I-CcaI, I-CcaII, and I-AstI) that belong to the I-OnuI family showed that they are functional HEs that cleave their respective cognate target sites. These endonucleases are encoded within group ID introns and appear to be orthologues that have inserted into 3 different mitochondrial genes: rns, rnl, and cox3. The endonuclease activity of I-CcaI was tested using various substrates, and its minimum DNA recognition sequence was estimated to be 26 nt. This set of HEs may provide some insight into how these types of mobile elements can migrate into new locations. This study provides additional endonucleases that can be added to the catalog of currently available HEs that may have various biotechnology applications.

Insertion of Group II Intron-Based Ribozyme Switches into Homing Endonuclease Genes

Fungal mitochondrial genomes act as "reservoirs" for homing endonucleases. These enzymes with their DNA site-specific cleavage activities are attractive tools for genome editing, targeted mutagenesis and gene therapy applications. Herein, we present strategies where homing endonuclease open reading frames (HEases ORFs) are interrupted with group II intron sequences. The ultimate goal is to achieve in vivo expression of HEases that can be regulated by manipulating the splicing efficiency of the HEase ORF-embedded group II introns. That addition of exogenous magnesium chloride (MgCl₂) appears to stimulate splicing of nonnative group II introns in Escherichia coli and the addition of cobalt chloride (CoCl₂) to the growth medium antagonizes the expression of HEase activity (i.e., splicing). Group II introns are potentially autocatalytic self-splicing elements and thus can be used as molecular switches that allow for temporal regulated HEase expression. This should be useful in precision genome engineering, mutagenesis, and minimizing off-target activities.

The diversity of mtDNA rns introns among strains of Ophiostoma piliferum, Ophiostoma pluriannulatum and related species

Background

Based on previous studies, it was suspected that the mitochondrial rns gene within the Ophiostomatales is rich in introns. This study focused on a collection of strains representing Ophiostoma piliferum, Ophiostoma pluriannulatum and related species that cause blue-stain; these fungi colonize the sapwood of trees and impart a dark stain. This reduces the value of the lumber. The goal was to examine the mtDNA rns intron landscape for these important blue stain fungi in order to facilitate future annotation of mitochondrial genomes (mtDNA) and to potentially identify mtDNA introns that can encode homing endonucleases which may have applications in biotechnology.

Results

Comparative sequence analysis identified five intron insertion sites among the ophiostomatoid fungi examined. Positions mS379 and mS952 harbor group II introns, the mS379 intron encodes a reverse transcriptase, and the mS952 intron encodes a potential homing endonuclease. Positions mS569, mS1224, and mS1247 have group I introns inserted and these encode intact or eroded homing endonuclease open reading frames (ORF). Phylogenetic analysis of the intron ORFs showed that they can be found in the same insertion site in closely and distantly related species.

Conclusions

Based on the molecular markers examined (rDNA internal transcribed spacers and rns introns), strains representing O. pilifera, O. pluriannulatum and Ophiostoma novae-zelandiae could not be resolved. Phylogenetic studies suggest that introns are gained and lost and that horizontal transfer could explain the presence of related intron in distantly related fungi. With regard to the mS379 group II intron, this study shows that mitochondrial group II introns and their reverse transcriptases may also follow the life cycle previously proposed for group I introns and their homing endonucleases. This consists of intron invasion, decay of intron ORF, loss of intron, and possible reinvasion.

Electronic supplementary material

The online version of this article (doi:10.1186/s40064-016-3076-6) contains supplementary material, which is available to authorized users.

Using Group II Introns for Attenuating the In Vitro and In Vivo Expression of a Homing Endonuclease

In Chaetomium thermophilum (DSM 1495) within the mitochondrial DNA (mtDNA) small ribosomal subunit (rns) gene a group IIA1 intron interrupts an open reading frame (ORF) encoded within a group I intron (mS1247). This arrangement offers the opportunity to examine if the nested group II intron could be utilized as a regulatory element for the expression of the homing endonuclease (HEase). Constructs were generated where the codon-optimized ORF was interrupted with either the native group IIA1 intron or a group IIB type intron. This study showed that the expression of the HEase (in vivo) in Escherichia coli can be regulated by manipulating the splicing efficiency of the HEase ORF-embedded group II introns. Exogenous magnesium chloride (MgCl₂) stimulated the expression of a functional HEase but the addition of cobalt chloride (CoCl₂) to growth media antagonized the expression of HEase activity. Ultimately the ability to attenuate HEase activity might be useful in precision genome engineering, minimizing off target activities, or where pathways have to be altered during a specific growth phase.

Genomic investigation reveals evolution and lifestyle adaptation of endophytic Staphylococcus epidermidis

Staphylococcus epidermidis is a major human associated bacterium and also an emerging nosocomial pathogen. There are reports of its association to rodents, sheep and plants. However, comparative and evolutionary studies of ecologically diverse strains of S. epidermidis are lacking. Here, we report the whole genome sequences of four S. epidermidis strains isolated from surface sterilized rice seeds along with genome sequence of type strain. Phylogenomic analysis of rice endophytic S. epidermidis (RESE) with “type strain” unequivocally established their species identity. Whole genome based tree of 93 strains of S. epidermidis revealed RESE as distinct sub-lineage which is more related to rodent sub-lineage than to majority of human lineage strains. Furthermore, comparative genomics revealed 20% variable gene-pool in S. epidermidis, suggesting that genomes of ecologically diverse strains are under flux. Interestingly, we were also able to map several genomic regions that are under flux and gave rise to RESE strains. The largest of these genomic regions encodes a cluster of genes unique to RESE that are known to be required for survival and stress tolerance, apart from those required for adaptation to plant habitat. The genomes and genes of RESE represent distinct ecological resource/sequences and provided first evolutionary insights into adaptation of S. epidermidis to plants.

Dynamics differentiate between active and inactive inteins

The balance between stability and dynamics for active enzymes can be somewhat quantified by studies of intein splicing and cleaving reactions. Inteins catalyze the ligation of flanking host exteins while excising themselves. The potential for applications led to engineering of a mini-intein splicing domain, where the homing endonuclease domain of the Mycobacterium tuberculosis RecA (Mtu recA) intein was removed. The remaining domains were linked by several short peptides, but splicing activity in all was substantially lower than the full-length intein. Native splicing activity was restored in some cases by a V67L mutation. Using computations and experiments, we examine the impact of this mutation on the stability and conformational dynamics of the mini-intein splicing domain. Molecular dynamics simulations were used to delineate the factors that determine the active state, including the V67L mini-intein mutant, and peptide linker. We found that (1) the V67L mutation lowers the global fluctuations in all modeled mini-inteins, stabilizing the mini-intein constructs; (2) the connecting linker length affects intein dynamics; and (3) the flexibilities of the linker and intein core are higher in the active structure. We have observed that the interaction of the linker region and a turn region around residues 35-41 provides the pathway for the allostery interaction. Our experiments reveal that intein catalysis is characterized by non-linear Arrhenius plot, confirming the significant contribution of protein conformational dynamics to intein function. We conclude that while the V67L mutation stabilizes the global structure, cooperative dynamics of all intein regions appear more important for intein function than high stability. Our studies suggest that effectively quenching the conformational dynamics of an intein through engineered allosteric interactions could deactivate intein splicing or cleaving.

I-OmiI and I-OmiII: two intron-encoded homing endonucleases within the Ophiostoma minus rns gene

The mitochondrial small subunit ribosomal RNA (rns) gene of the ascomycetous fungus Ophiostoma minus [strain WIN(M)371] was found to contain a group IC2 and a group IIB1 intron at positions mS569 and mS952 respectively. Both introns have open reading frames (ORFs) embedded that encode double motif LAGLIDADG homing endonucleases (I-OmiI and I-OmiII respectively). Codon-optimized versions of I-OmiI and I-OmiII were synthesized for overexpression in Escherichia coli. The in vitro characterization of I-OmiII showed that it is a functional homing endonuclease that cleaves the rns target site two nucleotides upstream (sense strand) of the intron insertion site generating 4 nucleotide 3' overhangs. The endonuclease activity of I-OmiII was tested using linear and circular substrates and cleavage activity was evaluated at various temperatures. The I-OmiI protein was expressed in E. coli, but purification was difficult, thus the endonuclease activity of this protein was tested via in vivo assays. Overall this study showed that there are many native forms of functional homing endonucleases yet to be discovered among fungal mtDNA genomes.

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

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.

A bacterial one-hybrid system to isolate homing endonuclease variants with altered DNA target specificities

Chromosomal cleavage near the site of mutations that cause disease can facilitate the targeted repair of the locus. Gene therapy protocols therefore require the engineering of DNA endonucleases that target specific genomic loci. Here, we describe a bacterial one-hybrid selection system that has been used to isolate derivatives of the I-SceI homing endonuclease from combinatorial libraries that display altered DNA recognition specificities. The construction of plasmid expression libraries, the development of reporter strains, and the utilization of these components in the bacterial one-hybrid system are detailed.

The susceptibility of Pseudomonas aeruginosa strains from cystic fibrosis patients to bacteriophages

Phage therapy may become a complement to antibiotics in the treatment of chronic Pseudomonas aeruginosa infection. To design efficient therapeutic cocktails, the genetic diversity of the species and the spectrum of susceptibility to bacteriophages must be investigated. Bacterial strains showing high levels of phage resistance need to be identified in order to decipher the underlying mechanisms. Here we have selected genetically diverse P. aeruginosa strains from cystic fibrosis patients and tested their susceptibility to a large collection of phages. Based on plaque morphology and restriction profiles, six different phages were purified from "pyophage", a commercial cocktail directed against five different bacterial species, including P. aeruginosa. Characterization of these phages by electron microscopy and sequencing of genome fragments showed that they belong to 4 different genera. Among 47 P. aeruginosa strains, 13 were not lysed by any of the isolated phages individually or by pyophage. We isolated two new phages that could lyse some of these strains, and their genomes were sequenced. The presence/absence of a CRISPR-Cas system (Clustered Regularly Interspaced Short Palindromic Repeats and Crisper associated genes) was investigated to evaluate the role of the system in phage resistance. Altogether, the results show that some P. aeruginosa strains cannot support the growth of any of the tested phages belonging to 5 different genera, and suggest that the CRISPR-Cas system is not a major defence mechanism against these lytic phages.

Flow cytometric assays for interrogating LAGLIDADG homing endonuclease DNA-binding and cleavage properties

A fast, easy, and scalable method to assess the properties of site-specific nucleases is crucial to -understanding their in cellulo behavior in genome engineering or population-level gene drive applications. Here we describe an analytical platform that enables high-throughput, semiquantitative interrogation of the DNA-binding and catalytic properties of LAGLIDADG homing endonucleases (LHEs). Using this platform, natural or engineered LHEs are expressed on the surface of Saccharomyces cerevisiae yeast where they can be rapidly evaluated against synthetic DNA target sequences using flow cytometry.

Comparative analysis of mitochondrial genomes of Rhizophagus irregularis - syn. Glomus irregulare - reveals a polymorphism induced by variability generating elements

Arbuscular mycorrhizal (AM) fungi are involved in one of the most widespread plant-fungus interactions. A number of studies on the population dynamics of AM fungi have used mitochondrial (mt) DNA sequences, and yet mt AM fungus genomes are poorly known. To date, four mt genomes of three species of AM fungi are available, among which are two from Rhizophagus irregularis. In order to study intra- and interstrain mt genome variability of R. irregularis, we sequenced and de novo assembled four additional mt genomes of this species. We used 454 pyrosequencing and Illumina technologies to directly sequence mt genomes from total genomic DNA. The mt genomes are unique within each strain. Interstrain divergences in genome size, as a result of highly polymorphic intergenic and intronic sequences, were observed. The polymorphism is brought about by three types of variability generating element (VGE): homing endonucleases, DNA polymerase domain-containing open reading frames and small inverted repeats. Based on VGE positioning, mt sequences and nuclear markers, two subclades of R. irregularis were characterized. The discovery of VGEs highlights the great intraspecific plasticity of the R. irregularis mt genome. VGEs allow the design of powerful mt markers for the typing and monitoring of R. irregularis strains in genetic and population studies.

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.

General principles of single-construct chromosomal gene drive

Gene drive systems are genetic elements capable of spreading into a population even if they confer a fitness cost to their host. We consider a class of drive systems consisting of a chromosomally located, linked cluster of genes, the presence of which renders specific classes of offspring arising from specific parental crosses unviable. Under permissive conditions, a number of these elements are capable of distorting the offspring ratio in their favor. We use a population genetic framework to derive conditions under which these elements spread to fixation in a population or induce a population crash. Many of these systems can be engineered using combinations of toxin and antidote genes, analogous to Medea, which consists of a maternal toxin and zygotic antidote. The majority of toxin–antidote drive systems require a critical frequency to be exceeded before they spread into a population. Of particular interest, a Z-linked Medea construct with a recessive antidote is expected to induce an all-male population crash for release frequencies above 50%. We suggest molecular tools that may be used to build these systems, and discuss their relevance to the control of a variety of insect pest species, including mosquito vectors of diseases such as malaria and dengue fever.

Restriction endonucleases: natural and directed evolution

Type II restriction endonucleases (REs) are highly sequence-specific compared with other classes of nucleases. PD-(D/E)XK nucleases, initially represented by only type II REs, now comprise a large and extremely diverse superfamily of proteins and, although sharing a structurally conserved core, typically display little or no detectable sequence similarity except for the active site motifs. Sequence similarity can only be observed in methylases and few isoschizomers. As a consequence, REs are classified according to combinations of functional properties rather than on the basis of genetic relatedness. New alignment matrices and classification systems based on structural core connectivity and cleavage mechanisms have been developed to characterize new REs and related proteins. REs recognizing more than 300 distinct specificities have been identified in RE database (REBASE: http://rebase.neb.com/cgi-bin/statlist ) but still the need for newer specificities is increasing due to the advancement in molecular biology and applications. The enzymes have undergone constant evolution through structural changes in protein scaffolds which include random mutations, homologous recombinations, insertions, and deletions of coding DNA sequences but rational mutagenesis or directed evolution delivers protein variants with new functions in accordance with defined biochemical or environmental pressures. Redesigning through random mutation, addition or deletion of amino acids, methylation-based selection, synthetic molecules, combining recognition and cleavage domains from different enzymes, or combination with domains of additional functions change the cleavage specificity or substrate preference and stability. There is a growing number of patents awarded for the creation of engineered REs with new and enhanced properties.

Comprehensive homing endonuclease target site specificity profiling reveals evolutionary constraints and enables genome engineering applications

Homing endonucleases (HEs) promote the evolutionary persistence of selfish DNA elements by catalyzing element lateral transfer into new host organisms. The high site specificity of this lateral transfer reaction, termed homing, reflects both the length (14-40 bp) and the limited tolerance of target or homing sites for base pair changes. In order to better understand molecular determinants of homing, we systematically determined the binding and cleavage properties of all single base pair variant target sites of the canonical LAGLIDADG homing endonucleases I-CreI and I-MsoI. These Chlorophyta algal HEs have very similar three-dimensional folds and recognize nearly identical 22 bp target sites, but use substantially different sets of DNA-protein contacts to mediate site-specific recognition and cleavage. The site specificity differences between I-CreI and I-MsoI suggest different evolutionary strategies for HE persistence. These differences also provide practical guidance in target site finding, and in the generation of HE variants with high site specificity and cleavage activity, to enable genome engineering applications.

Divalent metal ion differentially regulates the sequential nicking reactions of the GIY-YIG homing endonuclease I-BmoI

Homing endonucleases are site-specific DNA endonucleases that function as mobile genetic elements by introducing double-strand breaks or nicks at defined locations. Of the major families of homing endonucleases, the modular GIY-YIG endonucleases are least understood in terms of mechanism. The GIY-YIG homing endonuclease I-BmoI generates a double-strand break by sequential nicking reactions during which the single active site of the GIY-YIG nuclease domain must undergo a substantial reorganization. Here, we show that divalent metal ion plays a significant role in regulating the two independent nicking reactions by I-BmoI. Rate constant determination for each nicking reaction revealed that limiting divalent metal ion has a greater impact on the second strand than the first strand nicking reaction. We also show that substrate mutations within the I-BmoI cleavage site can modulate the first strand nicking reaction over a 314-fold range. Additionally, in-gel DNA footprinting with mutant substrates and modeling of an I-BmoI-substrate complex suggest that amino acid contacts to a critical GC-2 base pair are required to induce a bottom-strand distortion that likely directs conformational changes for reaction progress. Collectively, our data implies mechanistic roles for divalent metal ion and substrate bases, suggesting that divalent metal ion facilitates the re-positioning of the GIY-YIG nuclease domain between sequential nicking reactions.

Transposable elements in yeasts

With the development of new sequencing technologies in the past decade, yeast genomes have been extensively sequenced and their structures investigated. Transposable elements (TEs) are ubiquitous in eukaryotes and constitute a limited part of yeast genomes. However, due to their ability to move in genomes and generate dispersed repeated sequences, they contribute to modeling yeast genomes and thereby induce plasticity. This review assesses the TE contents of yeast genomes investigated so far. Their diversity and abundance at the inter- and intraspecific levels are presented, and their effects on gene expression and genome stability is considered. Recent results concerning TE-host interactions are also analyzed.

Comprehensive computational design of mCreI homing endonuclease cleavage specificity for genome engineering

Homing endonucleases (HEs) cleave long (∼ 20 bp) DNA target sites with high site specificity to catalyze the lateral transfer of parasitic DNA elements. In order to determine whether comprehensive computational design could be used as a general strategy to engineer new HE target site specificities, we used RosettaDesign (RD) to generate 3200 different variants of the mCreI LAGLIDADG HE towards 16 different base pair positions in the 22 bp mCreI target site. Experimental verification of a range of these designs demonstrated that over 2/3 (24 of 35 designs, 69%) had the intended new site specificity, and that 14 of the 15 attempted specificity shifts (93%) were achieved. These results demonstrate the feasibility of using structure-based computational design to engineer HE variants with novel target site specificities to facilitate genome engineering.

Social networking between mobile introns and their host genes

Homing endonucleases have long been known as the orchestrators of intron mobility. However, the extent of their influence on the intron and its genetic and cellular environment is still being elucidated. The accompanying paper emphasizes the importance of temporal control of endonuclease expression on splicing, expression of the host gene and cellular metabolism, while it raises questions to guide future inquiry.

Evolution of I-SceI homing endonucleases with increased DNA recognition site specificity

Elucidating how homing endonucleases undergo changes in recognition site specificity will facilitate efforts to engineer proteins for gene therapy applications. I-SceI is a monomeric homing endonuclease that recognizes and cleaves within an 18-bp target. It tolerates limited degeneracy in its target sequence, including substitution of a C:G(+4) base pair for the wild-type A:T(+4) base pair. Libraries encoding randomized amino acids at I-SceI residue positions that contact or are proximal to A:T(+4) were used in conjunction with a bacterial one-hybrid system to select I-SceI derivatives that bind to recognition sites containing either the A:T(+4) or the C:G(+4) base pairs. As expected, isolates encoding wild-type residues at the randomized positions were selected using either target sequence. All I-SceI proteins isolated using the C:G(+4) recognition site included small side-chain substitutions at G100 and either contained (K86R/G100T, K86R/G100S and K86R/G100C) or lacked (G100A, G100T) a K86R substitution. Interestingly, the binding affinities of the selected variants for the wild-type A:T(+4) target are 4- to 11-fold lower than that of wild-type I-SceI, whereas those for the C:G(+4) target are similar. The increased specificity of the mutant proteins is also evident in binding experiments in vivo. These differences in binding affinities account for the observed ∼36-fold difference in target preference between the K86R/G100T and wild-type proteins in DNA cleavage assays. An X-ray crystal structure of the K86R/G100T mutant protein bound to a DNA duplex containing the C:G(+4) substitution suggests how sequence specificity of a homing enzyme can increase. This biochemical and structural analysis defines one pathway by which site specificity is augmented for a homing endonuclease.

Mitochondrial genome from the facultative anaerobe and petite-positive yeast Dekkera bruxellensis contains the NADH dehydrogenase subunit genes

The progenitor of the Dekkera/Brettanomyces clade separated from the Saccharomyces/Kluyveromyces clade over 200 million years ago. However, within both clades, several lineages developed similar physiological traits. Both Saccharomyces cerevisiae and Dekkera bruxellensis are facultative anaerobes; in the presence of excess oxygen and sugars, they accumulate ethanol (Crabtree effect) and they both spontaneously generate respiratory-deficient mutants (petites). In order to understand the role of respiratory metabolism, the mitochondrial DNA (mtDNA) molecules of two Dekkera/Brettanomyces species were analysed. Dekkera bruxellensis mtDNA shares several properties with S. cerevisiae, such as the large genome size (76 453 bp), and the organization of the intergenic sequences consisting of spacious AT-rich regions containing a number of hairpin GC-rich cluster-like elements. In addition to a basic set of the mitochondrial genes coding for the components of cytochrome oxidase, cytochrome b, subunits of ATPase, two rRNA subunits and 25 tRNAs, D. bruxellensis also carries genes for the NADH dehydrogenase complex. Apparently, in yeast, the loss of this complex is not a precondition to develop a petite-positive, Crabtree-positive and anaerobic nature. On the other hand, mtDNA from a petite-negative Brettanomyces custersianus is much smaller (30 058 bp); it contains a similar gene set and has only short intergenic sequences.

Switching the mechanism of mating type switching: a domesticated transposase supplants a domesticated homing endonuclease

Programmed DNA rearrangements are critical for the development of many organisms and, intriguingly, can be catalyzed by domesticated mobile genetic elements. In this issue of Genes & Development, Barsoum and colleagues (pp. 33-44) demonstrate that, in the budding yeast Kluyveromyces lactis, a DNA rearrangement associated with mating type switching requires a domesticated transposase and occurs through a mechanism distinct from that in the related yeast, Saccharomyces cerevisiae. Thus, mechanisms for mating type switching have evolved multiple times, indicating the relative ease with which mobile genetic elements can be captured.

The Cryphonectria parasitica mitochondrial rns gene: plasmid-like elements, introns and homing endonucleases

The mt-rns gene of Cryphonectria parasitica is 9872bp long and includes two group I and two group II introns. An analysis of intronic protein-encoding sequences revealed that LAGLIDADG ORFs, which usually are associated with group I introns, were transferred at least twice into group II introns. A plasmid-like mitochondrial element (plME) that appears in high amounts in previously mutagen-induced mit1 and mit2 hypovirulent mutants of the Ep155 standard virulent strain of C. parasitica was found to be derived from a short region of the mt-rns gene, including the exon 1 and most of the first intron. The plME is a 4.2-kb circular, multimeric DNA and an autonomously-replicating mtDNA fragment. Although sexual transmission experiments indicate that the plME does not directly cause hypovirulence, its emergence is one manifestation of the many complex molecular and genetic events that appear to underlie this phenotype.

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.

An RNA hairpin sequesters the ribosome binding site of the homing endonuclease mobE gene

Previous transcript mapping of the bacteriophage Aeh1 nrd operon revealed a predicted RNA hairpin upstream of the homing endonuclease mobE gene. We enzymatically mapped the hairpin, showing that the mobE ribosome binding site is sequestered. Cloning of the hairpin upstream of lacZ resulted in reduced beta-galactosidase activity, consistent with translational regulation.

Engineering variants of the I-SceI homing endonuclease with strand-specific and site-specific DNA-nicking activity

The number of strand-specific nicking endonucleases that are currently available for laboratory procedures and applications in vivo is limited, and none is sufficiently specific to nick single target sites within complex genomes. The extreme target specificity of homing endonucleases makes them attractive candidates for engineering high-specificity nicking endonucleases. I-SceI is a monomeric homing enzyme that recognizes an 18 bp asymmetric target sequence, and cleaves both DNA strands to leave 3'-overhangs of 4 bp. In single turnover experiments using plasmid substrates, I-SceI generates transient open circle intermediates during the conversion of supercoiled to linear DNA, indicating that the enzyme cleaves the two DNA strands sequentially. A novel hairpin substrate was used to demonstrate that although wild-type I-SceI cleaves either the top or bottom DNA strand first to generate two nicked DNA intermediates, the enzyme has a preference for cleaving the bottom strand. The kinetics data are consistent with a parallel sequential reaction mechanism. Substitution of two pseudo-symmetric residues, Lys122 and Lys223, markedly reduces top and bottom-strand cleavage, respectively, to generate enzymes with significant strand- and sequence-specific nicking activity. The two active sites are partially interdependent, since alterations to one site affect the second. The kinetics analysis is consistent with X-ray crystal structures of I-SceI/DNA complexes that reveal a role for the lysines in establishing important solvent networks that include nucleophilic water molecules thought to attack the scissile phosphodiester bonds.

Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution

Group I and group II introns are different catalytic self-splicing and mobile RNA elements that contribute to genome dynamics. In this study, we have analyzed their distribution and evolution in 29 sequenced genomes from the Bacillus cereus group of bacteria. Introns were of different structural classes and evolutionary origins, and a large number of nearly identical elements are shared between multiple strains of different sources, suggesting recent lateral transfers and/or that introns are under a strong selection pressure. Altogether, 73 group I introns were identified, inserted in essential genes from the chromosome or newly described prophages, including the first elements found within phages in bacterial plasmids. Notably, bacteriophages are an important source for spreading group I introns between strains. Furthermore, 77 group II introns were found within a diverse set of chromosomal and plasmidic genes. Unusual findings include elements located within conserved DNA metabolism and repair genes and one intron inserted within a novel retroelement. Group II introns are mainly disseminated via plasmids and can subsequently invade the host genome, in particular by coupling mobility with host cell replication. This study reveals a very high diversity and variability of mobile introns in B. cereus group strains.

The sporadic occurrence of a group I intron-like element in the mtDNA rnl gene of Ophiostoma novo-ulmi subsp. americana

The presence of group I intron-like elements within the U7 region of the mtDNA large ribosomal subunit RNA gene (rnl) was investigated in strains of Ophiostoma novo-ulmi subsp. americana from Canada, Europe and Eurasia, and in selected strains of O. ips, O. minus, O. piceae, O. ulmi, and O. himal-ulmi. This insertion is of interest as it has been linked previously to the generation of plasmid-like mtDNA elements in diseased strains of O. novo-ulmi. Among 197 O. novo-ulmi subsp. americana strains tested, 61 contained a 1.6kb insertion within the rnl-U7 region and DNA sequence analysis suggests the presence of a group I intron (IA1 type) that encodes a potential double motif LAGLIDADG homing endonuclease-like gene (HEG). Phylogenetic analysis of rnl-U7 intron encoded HEG-like elements supports the view that double motif HEGs originated from a duplication event of a single-motif HEG followed by a fusion event that combined the two copies into one open reading frame (ORF). The data also show that rnl-U7 intron encoded ORFs belong to a clade that includes ORFs inserted into different types of group I introns, e.g. IB, ID, IC3, IA1, present within a variety of different mtDNA genes, such as the small ribosomal subunit RNA gene (rns), apo-cytochrome b gene (cob), NADH dehydrogenase subunit 5 (nad5), cytochrome oxidase subunit 1 gene (coxI), and ATPase subunit 9 gene (atp9). We also compared the occurrence of the rnl-U7 intron in our collection of 227 strains with the presence of the rnl-U11 group I intron and concluded that the U7 intron appears to be an optional element and the U11 intron is probably essential among the strains tested.

Altered target site specificity variants of the I-PpoI His-Cys box homing endonuclease

We used a yeast one-hybrid assay to isolate and characterize variants of the eukaryotic homing endonuclease I-PpoI that were able to bind a mutant, cleavage-resistant I-PpoI target or ‘homing’ site DNA in vivo. Native I-PpoI recognizes and cleaves a semi-palindromic 15-bp target site with high specificity in vivo and in vitro. This target site is present in the 28S or equivalent large subunit rDNA genes of all eukaryotes. I-PpoI variants able to bind mutant target site DNA had from 1 to 8 amino acid substitutions in the DNA–protein interface. Biochemical characterization of these proteins revealed a wide range of site–binding affinities and site discrimination. One-third of variants were able to cleave target site DNA, but there was no systematic relationship between site-binding affinity and site cleavage. Computational modeling of several variants provided mechanistic insight into how amino acid substitutions that contact, or are adjacent to, specific target site DNA base pairs determine I-PpoI site-binding affinity and site discrimination, and may affect cleavage efficiency.

Multiple controls regulate the expression of mobE, an HNH homing endonuclease gene embedded within a ribonucleotide reductase gene of phage Aeh1

Mobile genetic elements have the potential to influence the expression of genes surrounding their insertion site upon invasion of a genome. Here, we examine the transcriptional organization of a ribonucleotide reductase operon (nrd) that has been invaded by an HNH family homing endonuclease, mobE. In Aeromonas hydrophila phage Aeh1, mobE has inserted into the large-subunit gene (nrdA) of aerobic ribonucleotide reductase (RNR), splitting it into two smaller genes, nrdA-a and nrdA-b. This gene organization differs from that in phages T4, T6, RB2, RB3, RB15, and LZ7, where mobE is inserted in the nrdA-nrdB intergenic region. We present evidence that the expression of Aeh1 mobE is regulated by transcriptional, posttranscriptional, and translational controls. An Aeh1-specific late promoter drives expression of mobE, but strikingly the mobE transcript is processed internally at an RNase E-like site. We also identified a putative stem-loop structure upstream of mobE that sequesters the mobE ribosome binding site, presumably acting to down regulate MobE translation. Moreover, our transcriptional analyses indicate that the surrounding nrd genes of phage Aeh1 are expressed by a different strategy than are the corresponding phage T4 genes and that transcriptional readthrough is the only mechanism by which the promoterless Aeh1 nrdB gene is expressed. We suggest that the occurrence of multiple layers of control to limit the expression of mobE to late in the Aeh1 infection cycle is an adaptation of Aeh1 to reduce any effects on expression of the surrounding nrd genes early in phage infection when RNR function is critical.

Identification of a new subfamily of HNH nucleases and experimental characterization of a representative member, HphI restriction endonuclease

The restriction endonuclease (REase) R. HphI is a Type IIS enzyme that recognizes the asymmetric target DNA sequence 5'-GGTGA-3' and in the presence of Mg(2+) hydrolyzes phosphodiester bonds in both strands of the DNA at a distance of 8 nucleotides towards the 3' side of the target, producing a 1 nucleotide 3'-staggered cut in an unspecified sequence at this position. REases are typically ORFans that exhibit little similarity to each other and to any proteins in the database. However, bioinformatics analyses revealed that R.HphI is a member of a relatively big sequence family with a conserved C-terminal domain and a variable N-terminal domain. We predict that the C-terminal domains of proteins from this family correspond to the nuclease domain of the HNH superfamily rather than to the most common PD-(D/E)XK superfamily of nucleases. We constructed a three-dimensional model of the R.HphI catalytic domain and validated our predictions by site-directed mutagenesis and studies of DNA-binding and catalytic activities of the mutant proteins. We also analyzed the genomic neighborhood of R.HphI homologs and found that putative nucleases accompanied by a DNA methyltransferase (i.e. predicted REases) do not form a single group on a phylogenetic tree, but are dispersed among free-standing putative nucleases. This suggests that nucleases from the HNH superfamily were independently recruited to become REases in the context of RM systems multiple times in the evolution and that members of the HNH superfamily may be much more frequent among the so far unassigned REase sequences than previously thought.