New types of conserved sequence domains in DNA-binding regions of homing endonucleases

Prevalence of Group II Introns in Phage Genomes

Although bacteriophage genomes are under strong selective pressure for high coding density, they are still frequently invaded by mobile genetic elements (MGEs). Group II introns are MGEs that reduce host burden by autocatalytically splicing out of RNA before translation. While widely known in bacterial, archaeal, and eukaryotic organellar genomes, group II introns have been considered absent in phage. Identifying group II introns in genome sequences has previously been challenging because of their lack of primary sequence similarity. Advances in RNA structure-based homology searches using covariance models has provided the ability to identify the conserved secondary structures of group II introns. Here, we discover that group II introns are widely prevalent in phages from diverse phylogenetic backgrounds, from endosymbiont phage to jumbophage.

Functional Validation of Different Alternative Splicing Variants of the Chrysanthemum lavandulifolium ClNUM1 Gene in Tobacco

The Asteraceae are widely distributed throughout the world, with diverse functions and large genomes. Many of these genes remain undiscovered and unstudied. In this study, we discovered a new gene ClNUM1 in Chrysanthemum lavandulifolium and studied its function. In this study, bioinformatics, RT-qPCR, paraffin sectioning, and tobacco transgenics were utilized to bioinformatically analyze and functionally study the three variable splice variants of the unknown gene ClNUM1 cloned from C. lavandulifolium. The results showed that ClNUM1.1 and ClNUM1.2 had selective 3' splicing and selective 5' splicing, and ClNUM1.3 had selective 5' splicing. When the corresponding transgenic tobacco plants were subjected to abiotic stress treatment, in the tobacco seedlings, the ClNUM1.1 gene and the ClNUM1.2 gene enhanced salt and low-temperature tolerance and the ClNUM1.3 gene enhanced low-temperature tolerance; in mature tobacco plants, the ClNUM1.1 gene was able to enhance salt and low-temperature tolerance, and the ClNUM1.2 and ClNUM1.3 genes were able to enhance low-temperature tolerance. In summary, there are differences in the functions of the different splice variants and the different seedling stages of transgenic tobacco, but all of them enhanced the resistance of tobacco to a certain extent. The analysis and functional characterization of the ClNUM1 gene provided new potential genes and research directions for abiotic resistance breeding in Chrysanthemum.

DNA-binding and protein structure of nuclear factors likely acting in genetic information processing in the Paulinella chromatophore

The chromatophores in Paulinella are evolutionary-early-stage photosynthetic organelles. Biological processes in chromatophores depend on a combination of chromatophore and nucleus-encoded proteins. Interestingly, besides proteins carrying chromatophore-targeting signals, a large arsenal of short chromatophore-targeted proteins (sCTPs; <90 amino acids) without recognizable targeting signals were found in chromatophores. This situation resembles endosymbionts in plants and insects that are manipulated by host-derived antimicrobial peptides. Previously, we identified an expanded family of sCTPs of unknown function, named here "DNA-binding (DB)-sCTPs". DB-sCTPs contain a ~45 amino acid motif that is conserved in some bacterial proteins with predicted functions in DNA processing. Here, we explored antimicrobial activity, DNA-binding capacity, and structures of three purified recombinant DB-sCTPs. All three proteins exhibited antimicrobial activity against bacteria involving membrane permeabilization, and bound to bacterial lipids in vitro. A combination of in vitro assays demonstrated binding of recombinant DB-sCTPs to chromatophore-derived genomic DNA sequences with an affinity in the low nM range. Additionally, we report the 1.2 Å crystal structure of one DB-sCTP. In silico docking studies suggest that helix α2 inserts into the DNA major grove and the exposed residues, that are highly variable between different DB-sCTPs, confer interaction with the DNA bases. Identification of photosystem II subunit CP43 as a potential interaction partner of one DB-sCTP, suggests DB-sCTPs to be involved in more complex regulatory mechanisms. We hypothesize that membrane binding of DB-sCTPs is related to their import into chromatophores. Once inside, they interact with the chromatophore genome potentially providing nuclear control over genetic information processing.

The complete genome and comparative analysis of the phage phiC120 infecting multidrug-resistant Escherichia coli and Salmonella strains

Phages infecting Salmonella and Escherichia coli are promising agents for therapeutics and biological control of these foodborne pathogens, in particular those strains with resistance to several antibiotics. In an effort to assess the potential of the phage phiC120, a virulent phage isolated from horse feces in Mexico, we characterized its morphology, host range and complete genome. Herein, we showed that phiC120 possesses strong lytic activity against several multidrug-resistant E. coli O157: H7 and Salmonella strains, and its morphology indicated that is a member of Myoviridae family. The phiC120 genome is double-stranded DNA and consists of 186,570 bp in length with a 37.6% G + C content. A total of 281 putative open reading frames (ORFs) and two tRNAs were found, where 150 ORFs encoded hypothetical proteins with unknown function. Comparative analysis showed that phiC120 shared high similarity at nucleotide and protein levels with coliphages RB69 and phiE142. Detailed phiC120 analysis revealed that ORF 94 encodes a putative depolymerase, meanwhile genes encoding factors associated with lysogeny, toxins, and antibiotic resistance were absent; however, ORF 95 encodes a putative protein with potential allergenic and pro-inflammatory properties, making needed further studies to guarantee the safety of phiC120 for human use. The characterization of phiC120 expands our knowledge about the biology of coliphages and provides novel insights supporting its potential for the development of phage-based applications to control unwanted bacteria.

Two New Lytic Bacteriophages of the Myoviridae Family Against Carbapenem-Resistant Acinetobacter baumannii

Two lytic bacteriophages, WCHABP1 and WCHABP12, were recovered from hospital sewage and were able to infect 9 and 12 out of 18 carbapenem-resistant Acinetobacter baumannii clinical strains, which belonged to different clones. Electron microscopy scan showed that both bacteriophages had the similar morphology as those of the Myoviridae family. Whole genomic sequencing revealed 45.4- or 45.8-kb genome with a 37.6% GC content for WCHABP1 and WCHABP12, both of which showed significant DNA sequence similarity with bacteriophages of the Ap22virus genus within the Myoviridae family. Taxonomic analysis was therefore performed following the proposal approved by the International Committee on Taxonomy of Viruses, which confirmed that WCHABP1 and WCHABP12 represented two new species of the Ap22virus genus. No tRNAs but 88 and 89 open reading frames (ORFs) were predicted for the two bacteriophages, among which 22 and 21 had known function and encoded proteins for morphogenesis, packaging, lysis, and nucleiotide metabolism. The C-terminal amino acids of the large unit of fiber tail proteins varied between the bacteriophages, which may explain their different host ranges. For most lytic bacteriophages, a set of holin and endolysin are required for lysis. However, no known holin-encoding genes were identified in WCHABP1 and WCHABP12, suggesting that they may use alternative, yet-to-be-identified, novel holins for host cell membrane lysis. To test the efficacy of the bacteriophages in protecting against A. baumannii infection, a Galleria mellonella larva model was used. Only <20% G. mellonella larvae survived at 96 h after being infected by carbapenem-resistant A. baumannii strains, from which the two bacteriophages were recovered. With the administration of WCHABP1 and WCHABP12, the survival of larvae increased to 75%, while the treatment of polymyxin B only slightly increased the survival rate to 25%. The isolation of two new lytic bacteriophages in this study could expand our sight on Acinetobacter bacteriophages and may offer new potential therapeutic alternatives against A. baumannii.

Systematic classification of the His-Me finger superfamily

The His-Me finger endonucleases, also known as HNH or ββα-metal endonucleases, form a large and diverse protein superfamily. The His-Me finger domain can be found in proteins that play an essential role in cells, including genome maintenance, intron homing, host defense and target offense. Its overall structural compactness and non-specificity make it a perfectly-tailored pathogenic module that participates on both sides of inter- and intra-organismal competition. An extremely low sequence similarity across the superfamily makes it difficult to identify and classify new His-Me fingers. Using state-of-the-art distant homology detection methods, we provide an updated and systematic classification of His-Me finger proteins. In this work, we identified over 100 000 proteins and clustered them into 38 groups, of which three groups are new and cannot be found in any existing public domain database of protein families. Based on an analysis of sequences, structures, domain architectures, and genomic contexts, we provide a careful functional annotation of the poorly characterized members of this superfamily. Our results may inspire further experimental investigations that should address the predicted activity and clarify the potential substrates, to provide more detailed insights into the fundamental biological roles of these proteins.

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.

EheA from Exiguobacterium sp. yc3 is a novel thermostable DNase belonging to HNH endonuclease superfamily

The HNH endonuclease superfamily usually contains a conserved HNH motif in the sequence, and the second subfamily of it uses N to replace the second H in the HNH motif. A bacterium with extracellular thermostable DNase was isolated and identified as Exiguobacterium sp. yc3. A 20 kDa putative DNase was later purified and the encoding gene of it was amplified and sequenced, the deduced amino acid sequence analysis showed that the protein belongs to the HNH endonuclease superfamily, and therefore it was named as EheA ( E: xiguobacterium H: NH E: ndonuclease). Characterization of the recombinant EheA confirmed that EheA is a DNase. By site-directed mutation method, H116, N141 and N156 were proved to be essential for the DNase activity. EheA is the first experimentally determined bacterial source endonuclease belonging to the second subfamily of HNH superfamily. Further bioinformatic analysis showed that EheA homologue genes are conserved in the Exiguobacterium species, which suggests their possible important functions for Exiguobacterium species. And as a thermostable DNase, EheA also has a promising future in many application fields.

Function of the N-terminal segment of the RecA-dependent nuclease Ref

The bacteriophage P1 Ref (recombination enhancement function) protein is a RecA-dependent, HNH endonuclease. It can be directed to create targeted double-strand breaks within a displacement loop formed by RecA. The 76 amino acid N-terminal region of Ref is positively charged (25/76 amino acid residues) and inherently unstructured in solution. Our investigation of N-terminal truncation variants shows this region is required for DNA binding, contains a Cys involved in incidental dimerization and is necessary for efficient Ref-mediated DNA cleavage. Specifically, Ref N-terminal truncation variants lacking between 21 and 47 amino acids are more effective RecA-mediated targeting nucleases. We propose a more refined set of options for the Ref-mediated cleavage mechanism, featuring the N-terminal region as an anchor for at least one of the DNA strand cleavage events.

LTR retroelements are intrinsic components of transcriptional networks in frogs

Background

LTR retroelements (LTR REs) constitute a major group of transposable elements widely distributed in eukaryotic genomes. Through their own mechanism of retrotranscription LTR REs enrich the genomic landscape by providing genetic variability, thus contributing to genome structure and organization. Nonetheless, transcriptomic activity of LTR REs still remains an obscure domain within cell, developmental, and organism biology.

Results

Here we present a first comparative analysis of LTR REs for anuran amphibians based on a full depth coverage transcriptome of the European pool frog, Pelophylax lessonae, the genome of the African clawed frog, Silurana tropicalis (release v7.1), and additional transcriptomes of S. tropicalis and Cyclorana alboguttata. We identified over 1000 copies of LTR REs from all four families (Bel/Pao, Ty1/Copia, Ty3/Gypsy, Retroviridae) in the genome of S. tropicalis and discovered transcripts of several of these elements in all RNA-seq datasets analyzed. Elements of the Ty3/Gypsy family were most active, especially Amn-san elements, which accounted for approximately 0.27% of the genome in Silurana. Some elements exhibited tissue specific expression patterns, for example Hydra1.1 and MuERV-like elements in Pelophylax. In S. tropicalis considerable transcription of LTR REs was observed during embryogenesis as soon as the embryonic genome became activated, i.e. at midblastula transition. In the course of embryonic development the spectrum of transcribed LTR REs changed; during gastrulation and neurulation MuERV-like and SnRV like retroviruses were abundantly transcribed while during organogenesis transcripts of the XEN1 retroviruses became much more active.

Conclusions

The differential expression of LTR REs during embryogenesis in concert with their tissue-specificity and the protein domains they encode are evidence for the functional roles these elements play as integrative parts of complex regulatory networks. Our results support the meanwhile widely accepted concept that retroelements are not simple “junk DNA” or “harmful genomic parasites” but essential components of the transcriptomic machinery in vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-626) contains supplementary material, which is available to authorized users.

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

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

The monomeric GIY-YIG homing endonuclease I-BmoI uses a molecular anchor and a flexible tether to sequentially nick DNA

The GIY-YIG nuclease domain is found within protein scaffolds that participate in diverse cellular pathways and contains a single active site that hydrolyzes DNA by a one-metal ion mechanism. GIY-YIG homing endonucleases (GIY-HEs) are two-domain proteins with N-terminal GIY-YIG nuclease domains connected to C-terminal DNA-binding and they are thought to function as monomers. Using I-BmoI as a model GIY-HE, we test mechanisms by which the single active site is used to generate a double-strand break. We show that I-BmoI is partially disordered in the absence of substrate, and that the GIY-YIG domain alone has weak affinity for DNA. Significantly, we show that I-BmoI functions as a monomer at all steps of the reaction pathway and does not transiently dimerize or use sequential transesterification reactions to cleave substrate. Our results are consistent with the I-BmoI DNA-binding domain acting as a molecular anchor to tether the GIY-YIG domain to substrate, permitting rotation of the GIY-YIG domain to sequentially nick each DNA strand. These data highlight the mechanistic differences between monomeric GIY-HEs and dimeric or tetrameric GIY-YIG restriction enzymes, and they have implications for the use of the GIY-YIG domain in genome-editing applications.

Correlation between structure, protein composition, morphogenesis and cytopathology of Glossina pallidipes salivary gland hypertrophy virus

The Glossina pallidipes salivary gland hypertrophy virus (GpSGHV) is a dsDNA virus with rod-shaped, enveloped virions. Its 190 kb genome contains 160 putative protein-coding ORFs. Here, the structural components, protein composition and associated aspects of GpSGHV morphogenesis and cytopathology were investigated. Four morphologically distinct structures: the nucleocapsid, tegument, envelope and helical surface projections, were observed in purified GpSGHV virions by electron microscopy. Nucleocapsids were present in virogenic stroma within the nuclei of infected salivary gland cells, whereas enveloped virions were located in the cytoplasm. The cytoplasm of infected cells appeared disordered and the plasma membranes disintegrated. Treatment of virions with 1 % NP-40 efficiently partitioned the virions into envelope and nucleocapsid fractions. The fractions were separated by SDS-PAGE followed by in-gel trypsin digestion and analysis of the tryptic peptides by liquid chromatography coupled to electrospray and tandem mass spectrometry. Using the MaxQuant program with Andromeda as a database search engine, a total of 45 viral proteins were identified. Of these, ten and 15 were associated with the envelope and the nucleocapsid fractions, respectively, whilst 20 were detected in both fractions, most likely representing tegument proteins. In addition, 51 host-derived proteins were identified in the proteome of the virus particle, 13 of which were verified to be incorporated into the mature virion using a proteinase K protection assay. This study provides important information about GpSGHV biology and suggests options for the development of future anti-GpSGHV strategies by interfering with virus-host interactions.

Genome sequence of the phage clP1, which infects the beer spoilage bacterium Pediococcus damnosus

Pediococcus damnosus (P. damnosus) bacteriophage (phage) clP1 is a novel virulent phage isolated from a municipal sewage sample collected in Southern Ireland. This phage infects the beer spoilage strain P. damnosus P82 which was isolated from German breweries. Sequencing of the phage has revealed a linear double stranded DNA genome of 38,013 base pairs (bp) with an overall GC content of 47.6%. Fifty seven open reading frames (ORFs) were identified of which 30 showed homology to previously sequenced proteins, and as a consequence 20 of these were assigned predicted functions. The majority of genes displayed homology with genes from the Lactobacillus plantarum phage phiJL-1. All genes were located on the same coding strand and in the same orientation. Morphological characterisation placed phage clP1 as a member of the Siphoviridae family with an isometric head (59 nm diameter) and non-contractile tail (length 175 nm; diameter 10nm. Interestingly, the phage clP1 genome was found to share very limited identity with other phage genome sequences in the database, and was hence considered unique. This was highlighted by the genome organisation which differed slightly to the consensus pattern of genomic organisation usually found in Siphoviridae phages. With the genetic machinery present for a lytic lifecycle and the absence of potential endotoxin factors, this phage may have applications in the biocontrol of beer spoilage bacteria. To our knowledge, this study represents the first reported P. damnosus phage genome sequence.

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

Novel family of putative homing endonuclease genes was recently discovered during analyses of metagenomic and genomic sequence data. One such protein is encoded within a group I intron that resides in the recA gene of the Bacillus thuringiensis 0305ϕ8–36 bacteriophage. Named I-Bth0305I, the endonuclease cleaves a DNA target in the uninterrupted recA gene at a position immediately adjacent to the intron insertion site. The enzyme displays a multidomain, homodimeric architecture and footprints a DNA region of ∼60 bp. Its highest specificity corresponds to a 14-bp pseudopalindromic sequence that is directly centered across the DNA cleavage site. Unlike many homing endonucleases, the specificity profile of the enzyme is evenly distributed across much of its target site, such that few single base pair substitutions cause a significant decrease in cleavage activity. A crystal structure of its C-terminal domain confirms a nuclease fold that is homologous to very short patch repair (Vsr) endonucleases. The domain architecture and DNA recognition profile displayed by I-Bth0305I, which is the prototype of a homing lineage that we term the ‘EDxHD’ family, are distinct from previously characterized homing endonucleases.

Genomic analysis of bacteriophage ϕAB1, a ϕKMV-like virus infecting multidrug-resistant Acinetobacter baumannii

We present the complete genomic sequence of a lytic bacteriophage ϕAB1 which can infect many clinical isolates of multidrug-resistant Acinetobacter baumannii. The recently isolated bacteriophage displays morphology resembling Podoviridae family. The ϕAB1 genome is a linear double-stranded DNA of 41,526 bp containing 46 possible open reading frames (ORFs). The majority of the predicted structural proteins were identified as part of the phage particle by mass spectrometry analysis. According to the virion morphology, overall genomic structure, and the phylogenetic tree of RNA polymerase, we propose that ϕAB1 is a new member of the ϕKMV-like phages. Additionally, we identified four ORFs encoding putative HNH endonucleases, one of which is presumed to integrate and create a genes-in-pieces DNA polymerase. Also, a potential lysis cassette was identified in the late genome. The lytic power of this bacteriophage combined with its specificity for A. baumannii makes ϕAB1 an attractive agent for therapeutic or disinfection applications.

A unified genetic, computational and experimental framework identifies functionally relevant residues of the homing endonuclease I-BmoI

Insight into protein structure and function is best obtained through a synthesis of experimental, structural and bioinformatic data. Here, we outline a framework that we call MUSE (mutual information, unigenic evolution and structure-guided elucidation), which facilitated the identification of previously unknown residues that are relevant for function of the GIY-YIG homing endonuclease I-BmoI. Our approach synthesizes three types of data: mutual information analyses that identify co-evolving residues within the GIY-YIG catalytic domain; a unigenic evolution strategy that identifies hyper- and hypo-mutable residues of I-BmoI; and interpretation of the unigenic and co-evolution data using a homology model. In particular, we identify novel positions within the GIY-YIG domain as functionally important. Proof-of-principle experiments implicate the non-conserved I71 as functionally relevant, with an I71N mutant accumulating a nicked cleavage intermediate. Moreover, many additional positions within the catalytic, linker and C-terminal domains of I-BmoI were implicated as important for function. Our results represent a platform on which to pursue future studies of I-BmoI and other GIY-YIG-containing proteins, and demonstrate that MUSE can successfully identify novel functionally critical residues that would be ignored in a traditional structure-function analysis within an extensively studied small domain of approximately 90 amino acids.

Homing endonucleases: from basics to therapeutic applications

Homing endonucleases (HE) are double-stranded DNAses that target large recognition sites (12-40 bp). HE-encoding sequences are usually embedded in either introns or inteins. Their recognition sites are extremely rare, with none or only a few of these sites present in a mammalian-sized genome. However, these enzymes, unlike standard restriction endonucleases, tolerate some sequence degeneracy within their recognition sequence. Several members of this enzyme family have been used as templates to engineer tools to cleave DNA sequences that differ from their original wild-type targets. These custom HEs can be used to stimulate double-strand break homologous recombination in cells, to induce the repair of defective genes with very low toxicity levels. The use of tailored HEs opens up new possibilities for gene therapy in patients with monogenic diseases that can be treated ex vivo. This review provides an overview of recent advances in this field.

Structure of bacteriophage T4 endonuclease II mutant E118A, a tetrameric GIY-YIG enzyme

Coliphage T4 endonuclease II (EndoII), encoded by gene denA, is a small (16 kDa, 136 aa) enzyme belonging to the GIY-YIG family of endonucleases, which lacks a C-terminal domain corresponding to that providing most of the binding energy in the structurally characterized GIY-YIG endonucleases, I-TevI and UvrC. In vivo, it is involved in degradation of host DNA, permitting scavenging of host-derived nucleotides for phage DNA synthesis. EndoII primarily catalyzes single-stranded nicking of DNA; 5- to 10-fold less frequently double-stranded breaks are produced. The Glu118Ala mutant of EndoII was crystallized in space group P2(1) with four monomers in the asymmetric unit. The fold of the EndoII monomer is similar to that of the catalytic domains of UvrC and I-TevI. In contrast to these enzymes, EndoII forms a striking X-shaped tetrameric structure composed as a dimer of dimers, with a protruding hairpin domain not present in UvrC or I-TevI providing most of the dimerization and tetramerization interfaces. A bound phosphate ion in one of the four active sites of EndoII likely mimics the scissile phosphate in a true substrate complex. In silico docking experiments showed that a protruding loop containing a nuclease-associated modular domain 3 element is likely to be involved in substrate binding, as well as residues forming a separate nucleic acid binding surface adjacent to the active site. The positioning of these sites within the EndoII primary dimer suggests that the substrate would bind to a primary EndoII dimer diagonally over the active sites, requiring significant distortion of the enzyme or the substrate DNA, or both, for simultaneous nicking of both DNA strands. The scarcity of potential nucleic acid binding residues between the active sites indicates that EndoII may bind its substrate inefficiently across the two sites in the dimer, offering a plausible explanation for the catalytic preponderance of single-strand nicks. Mutations analyzed in earlier functional studies are discussed in their structural context.

Biochemical and mutagenic analysis of I-CreII reveals distinct but important roles for both the H-N-H and GIY-YIG motifs

Homing endonucleases typically contain one of four conserved catalytic motifs, and other elements that confer tight DNA binding. I-CreII, which catalyzes homing of the Cr.psbA4 intron, is unusual in containing two potential catalytic motifs, H-N-H and GIY-YIG. Previously, we showed that cleavage by I-CreII leaves ends (2-nt 3' overhangs) that are characteristic of GIY-YIG endonucleases, yet it has a relaxed metal requirement like H-N-H enzymes. Here we show that I-CreII can bind DNA without an added metal ion, and that it binds as a monomer, akin to GIY-YIG enzymes. Moreover, cleavage of supercoiled DNA, and estimates of strand-specific cleavage rates, suggest that I-CreII uses a sequential cleavage mechanism. Alanine substitution of a number of residues in the GIY-YIG motif, however, did not block cleavage activity, although DNA binding was substantially reduced in several variants. Substitution of conserved histidines in the H-N-H motif resulted in variants that did not promote DNA cleavage, but retained high-affinity DNA binding-thus identifying it as the catalytic motif. Unlike the non-specific H-N-H colicins, however; substitution of the conserved asparagine substantially reduced DNA binding (though not the ability to promote cleavage). These results indicate that, in I-CreII, two catalytic motifs have evolved to play important roles in specific DNA binding. The data also indicate that only the H-N-H motif has retained catalytic ability.

Amino acid residues in the GIY-YIG endonuclease II of phage T4 affecting sequence recognition and binding as well as catalysis

Phage T4 endonuclease II (EndoII), a GIY-YIG endonuclease lacking a carboxy-terminal DNA-binding domain, was subjected to site-directed mutagenesis to investigate roles of individual amino acids in substrate recognition, binding, and catalysis. The structure of EndoII was modeled on that of UvrC. We found catalytic roles for residues in the putative catalytic surface (G49, R57, E118, and N130) similar to those described for I-TevI and UvrC; in addition, these residues were found to be important for substrate recognition and binding. The conserved glycine (G49) and arginine (R57) were essential for normal sequence recognition. Our results are in agreement with a role for these residues in forming the DNA-binding surface and exposing the substrate scissile bond at the active site. The conserved asparagine (N130) and an adjacent proline (P127) likely contribute to positioning the catalytic domain correctly. Enzymes in the EndoII subfamily of GIY-YIG endonucleases share a strongly conserved middle region (MR, residues 72 to 93, likely helical and possibly substituting for heterologous helices in I-TevI and UvrC) and a less strongly conserved N-terminal region (residues 12 to 24). Most of the conserved residues in these two regions appeared to contribute to binding strength without affecting the mode of substrate binding at the catalytic surface. EndoII K76, part of a conserved NUMOD3 DNA-binding motif of homing endonucleases found to overlap the MR, affected both sequence recognition and catalysis, suggesting a more direct involvement in positioning the substrate. Our data thus suggest roles for the MR and residues conserved in GIY-YIG enzymes in recognizing and binding the substrate.

Stabilization of pSW100 from Pantoea stewartii by the F conjugation system

Plasmid pSW100 is 1 of the 13 plasmids from Pantoea stewartii subsp. stewartii SW2 which has a replicon that resembles that of ColE1. This work uses a pSW100 derivative, pSW140K, to study how the pSW100 replicon is stably maintained in its hosts. Our results indicate that although pSW140K is stable in Escherichia coli HB101, the plasmid is rapidly lost in another E. coli strain, DH5alpha, indicating that the genetic background of an E. coli strain affects the stability of pSW140K. Mutagenesis of E. coli HB101 with EZ::TN <DHFR-1> revealed that mutations in traC, traF, traG, traN, and traV, which encode the components of the sex pilus assembly, reduce plasmid stability. Furthermore, this work identified that a 38-bp region located immediately upstream of the RNAII promoter is critical to the maintenance of plasmid stability in E. coli HB101. TraC binds to the region, and in addition, deleting the region destabilizes the plasmid. Furthermore, inserting this 38-bp fragment into a plasmid that contains the minimal replicon from pSW200 stabilizes the plasmid in E. coli HB101. Fluorescence in situ hybridization and immunofluorescence staining also revealed that derivatives of pSW100, pSW128A, and TraC are colocalized in cells, suggesting that pSW100 may use the sex pilus assembly as a partition apparatus to ensure the even distribution of the plasmid during cell division, which may thus maintain the plasmid's stability.

Phage T4 SegB protein is a homing endonuclease required for the preferred inheritance of T4 tRNA gene region occurring in co-infection with a related phage

Homing endonucleases initiate nonreciprocal transfer of DNA segments containing their own genes and the flanking sequences by cleaving the recipient DNA. Bacteriophage T4 segB gene, which is located in a cluster of tRNA genes, encodes a protein of unknown function, homologous to homing endonucleases of the GIY-YIG family. We demonstrate that SegB protein is a site-specific endonuclease, which produces mostly 3′ 2-nt protruding ends at its DNA cleavage site. Analysis of SegB cleavage sites suggests that SegB recognizes a 27-bp sequence. It contains 11-bp conserved sequence, which corresponds to a conserved motif of tRNA TψC stem-loop, whereas the remainder of the recognition site is rather degenerate. T4-related phages T2L, RB1 and RB3 contain tRNA gene regions that are homologous to that of phage T4 but lack segB gene and several tRNA genes. In co-infections of phages T4 and T2L, segB gene is inherited with nearly 100% of efficiency. The preferred inheritance depends absolutely on the segB gene integrity and is accompanied by the loss of the T2L tRNA gene region markers. We suggest that SegB is a homing endonuclease that functions to ensure spreading of its own gene and the surrounding tRNA genes among T4-related phages.

Occurrence of protein structure elements in conserved sequence regions

BACKGROUND

Conserved protein sequence regions are extremely useful for identifying and studying functionally and structurally important regions. By means of an integrated analysis of large-scale protein structure and sequence data, structural features of conserved protein sequence regions were identified.

RESULTS

Helices and turns were found to be underrepresented in conserved regions, while strands were found to be overrepresented. Similar numbers of loops were found in conserved and random regions.

CONCLUSION

These results can be understood in light of the structural constraints on different secondary structure elements, and their role in protein structural stabilization and topology. Strands can tolerate fewer sequence changes and nonetheless keep their specific shape and function. They thus tend to be more conserved than helices, which can keep their shape and function with more changes. Loop behavior can be explained by the presence of both constrained and freely changing loops in proteins. Our detailed statistical analysis of diverse proteins links protein evolution to the biophysics of protein thermodynamic stability and folding. The basic structural features of conserved sequence regions are also important determinants of protein structure motifs and their function.

Phylogenomic analysis of the GIY-YIG nuclease superfamily

Background

The GIY-YIG domain was initially identified in homing endonucleases and later in other selfish mobile genetic elements (including restriction enzymes and non-LTR retrotransposons) and in enzymes involved in DNA repair and recombination. However, to date no systematic search for novel members of the GIY-YIG superfamily or comparative analysis of these enzymes has been reported.

Results

We carried out database searches to identify all members of known GIY-YIG nuclease families. Multiple sequence alignments together with predicted secondary structures of identified families were represented as Hidden Markov Models (HMM) and compared by the HHsearch method to the uncharacterized protein families gathered in the COG, KOG, and PFAM databases. This analysis allowed for extending the GIY-YIG superfamily to include members of COG3680 and a number of proteins not classified in COGs and to predict that these proteins may function as nucleases, potentially involved in DNA recombination and/or repair. Finally, all old and new members of the GIY-YIG superfamily were compared and analyzed to infer the phylogenetic tree.

Conclusion

An evolutionary classification of the GIY-YIG superfamily is presented for the very first time, along with the structural annotation of all (sub)families. It provides a comprehensive picture of sequence-structure-function relationships in this superfamily of nucleases, which will help to design experiments to study the mechanism of action of known members (especially the uncharacterized ones) and will facilitate the prediction of function for the newly discovered ones.

Distance determination by GIY-YIG intron endonucleases: discrimination between repression and cleavage functions

GIY-YIG homing endonucleases are modular proteins, with conserved N-terminal catalytic domains connected by linkers to C-terminal DNA-binding domains. I-TevI, the T4 phage GIY-YIG intron endonuclease, functions both in promoting td intron homing, and in acting as a transcriptional autorepressor. Repression is achieved by binding to an operator, which is cleaved at 100-fold reduced efficiency relative to the intronless homing site. The linker includes a zinc finger, which functions in distance determination, to constrain the catalytic domain to cleave the homing site at a fixed position. Here we show that I-BmoI, a related GIY-YIG endonuclease lacking a zinc finger, also possesses some cleavage distance discrimination. Furthermore, hybrid endonucleases constructed by swapping the domains of I-BmoI and I-TevI are active, precise and demonstrate that features other than the zinc finger facilitate distance determination. Most importantly, I-TevI zinc finger mutants cleave the operator more efficiently than the homing site, the converse of wild-type protein. These results are consistent with the zinc finger acting as a measuring device, directing efficient cleavage of the homing site to promote intron mobility, while reducing cleavage at the operator to ensure transcriptional autorepression and phage viability.

Genome Comparison of Pseudomonas aeruginosa Large Phages

Pseudomonas aeruginosa phage EL is a dsDNA phage related to the giant phiKZ-like Myoviridae. The EL genome sequence comprises 211,215 bp and has 201 predicted open reading frames (ORFs). The EL genome does not share DNA sequence homology with other viruses and micro-organisms sequenced to date. However, one-third of the predicted EL gene products (gps) shares similarity (Blast alignments of 17-55% amino acid identity) with phiKZ proteins. Comparative EL and phiKZ genomics reveals that these giant phages are an example of substantially diverged genetic mosaics. Based on the position of similar EL and phiKZ predicted gene products, five genome regions can be delineated in EL, four of which are relatively conserved between EL and phiKZ. Region IV, a 17.7 kb genome region with 28 predicted ORFs, is unique to EL. Fourteen EL ORFs have been assigned a putative function based on protein similarity. Assigned proteins are involved in DNA replication and nucleotide metabolism (NAD+-dependent DNA ligase, ribonuclease HI, helicase, thymidylate kinase), host lysis and particle structure. EL-gp146 is the first chaperonin GroEL sequence identified in a viral genome. Besides a putative transposase, EL harbours predicted mobile endonucleases related to H-N-H and LAGLIDADG homing endonucleases associated with group I intron and intein intervening sequences.

One-Block CYRCA: an automated procedure for identifying multiple-block alignments from single block queries

One-Block CYRCA is an automated procedure for identifying multiple-block alignments from single block queries (). It is based on the LAMA and CYRCA block-to-block alignment methods. The procedure identifies whether the query blocks can form new multiple-block alignments (block sets) with blocks from a database or join pre-existing database block sets. Using pre-computed LAMA block alignments and CYRCA sets from the Blocks database reduces the computation time. LAMA and CYRCA are highly sensitive and selective methods that can augment many other sequence analysis approaches.

Discovery of the principal specific transcription factors of Apicomplexa and their implication for the evolution of the AP2-integrase DNA binding domains

The comparative genomics of apicomplexans, such as the malarial parasite Plasmodium, the cattle parasite Theileria and the emerging human parasite Cryptosporidium, have suggested an unexpected paucity of specific transcription factors (TFs) with DNA binding domains that are closely related to those found in the major families of TFs from other eukaryotes. This apparent lack of specific TFs is paradoxical, given that the apicomplexans show a complex developmental cycle in one or more hosts and a reproducible pattern of differential gene expression in course of this cycle. Using sensitive sequence profile searches, we show that the apicomplexans possess a lineage-specific expansion of a novel family of proteins with a version of the AP2 (Apetala2)-integrase DNA binding domain, which is present in numerous plant TFs. About 20–27 members of this apicomplexan AP2 (ApiAP2) family are encoded in different apicomplexan genomes, with each protein containing one to four copies of the AP2 DNA binding domain. Using gene expression data from Plasmodium falciparum, we show that guilds of ApiAP2 genes are expressed in different stages of intraerythrocytic development. By analogy to the plant AP2 proteins and based on the expression patterns, we predict that the ApiAP2 proteins are likely to function as previously unknown specific TFs in the apicomplexans and regulate the progression of their developmental cycle. In addition to the ApiAP2 family, we also identified two other novel families of AP2 DNA binding domains in bacteria and transposons. Using structure similarity searches, we also identified divergent versions of the AP2-integrase DNA binding domain fold in the DNA binding region of the PI-SceI homing endonuclease and the C-terminal domain of the pleckstrin homology (PH) domain-like modules of eukaryotes. Integrating these findings, we present a reconstruction of the evolutionary scenario of the AP2-integrase DNA binding domain fold, which suggests that it underwent multiple independent combinations with different types of mobile endonucleases or recombinases. It appears that the eukaryotic versions have emerged from versions of the domain associated with mobile elements, followed by independent lineage-specific expansions, which accompanied their recruitment to transcription regulation functions.

Homing into the origin of the AP2 DNA binding domain

The AP2 DNA binding domain was thought to be plant specific because of its presence in plant, but not animal, transcriptional regulators, particularly members of the AP2/ERF family. Two recent studies have identified the AP2 domain in bacteria, bacteriophage and a ciliate as part of proteins that also encode site-specific endonucleases. The association of AP2 with an enzyme known to catalyze its own movement within populations and between species explains the unusual distribution of AP2 and, as such, adds to a growing list of phenomena where mobile DNA has promoted evolutionary novelty.

DNA binding and cleavage by the HNH homing endonuclease I-HmuI

The structure of I-HmuI, which represents the last family of homing endonucleases without a defining crystallographic structure, has been determined in complex with its DNA target. A series of diverse protein structural domains and motifs, contacting sequential stretches of nucleotide bases, are distributed along the DNA target. I-HmuI contains an N-terminal domain with a DNA-binding surface found in the I-PpoI homing endonuclease and an associated HNH/N active site found in the bacterial colicins, and a C-terminal DNA-binding domain previously observed in the I-TevI homing endonuclease. The combination and exchange of these features between protein families indicates that the genetic mobility associated with homing endonucleases extends to the level of independent structural domains. I-HmuI provides an unambiguous structural connection between the His-Cys box endonucleases and the bacterial colicins, supporting the hypothesis that these enzymes diverged from a common ancestral nuclease.

From endonucleases to transcription factors: evolution of the AP2 DNA binding domain in plants

All members of the AP2/ERF family of plant transcription regulators contain at least one copy of a DNA binding domain called the AP2 domain. The AP2 domain has been considered plant specific. Here, we show that homologs are present in the cyanobacterium Trichodesmium erythraeum, the ciliate Tetrahymena thermophila, and the viruses Enterobacteria phage Rb49 and Bacteriophage Felix 01. We demonstrate that the T. erythraeum AP2 domain selectively binds stretches of poly(dG)/poly(dC) showing functional conservation with plant AP2/ERF proteins. The newly discovered nonplant proteins bearing an AP2 domain are predicted to be HNH endonucleases. Sequence conservation extends outside the AP2 domain to include part of the endonuclease domain for the T. erythraeum protein and some plant AP2/ERF proteins. Our phylogenetic analysis of the broader family of AP2 domains supports the possibility of lateral gene transfer. We hypothesize that a horizontal transfer of an HNH-AP2 endonuclease from bacteria or viruses into plants may have led to the origin of the AP2/ERF family of transcription factors via transposition and homing processes.