Molecular evolution of bacteriophages: evidence of selection against the recognition sites of host restriction enzymes

Restriction enzymes produced by bacteria serve as a defense against invading bacteriophages, and so phages without other protection would be expected to undergo selection to eliminate recognition sites for these enzymes from their genomes. The observed frequencies of all restriction sites in the genomes of all completely sequenced DNA phages (T7, lambda, phi X174, G4, M13, f1, fd, and IKe) have been compared to expected frequencies derived from trinucleotide frequencies. Attention was focused on 6-base palindromes since they comprise the typical recognition sites for type II restriction enzymes. All of these coliphages, with the exception of lambda and G4, exhibit significant avoidance of the particular sequences that are enterobacterial restriction sites. As expected, the sequenced fraction of the genome of phi 29, a Bacillus subtilis phage, lacks Bacillus restriction sites. By contrast, the RNA phage MS2, several viruses that infect eukaryotes (EBV, adenovirus, papilloma, and SV40), and three mitochondrial genomes (human, mouse, and cow) were found not to lack restriction sites. Because the particular palindromes avoided correspond closely with the recognition sites for host enzymes and because other viruses and small genomes do not show this avoidance, it is concluded that the effect indeed results from natural selection.      

Using shotgun sequence data to find active restriction enzyme genes

Whole genome shotgun sequence analysis has become the standard method for beginning to determine a genome sequence. The preparation of the shotgun sequence clones is, in fact, a biological experiment. It determines which segments of the genome can be cloned into Escherichia coli and which cannot. By analyzing the complete set of sequences from such an experiment, it is possible to identify genes lethal to E. coli. Among this set are genes encoding restriction enzymes which, when active in E. coli, lead to cell death by cleaving the E. coli genome at the restriction enzyme recognition sites. By analyzing shotgun sequence data sets we show that this is a reliable method to detect active restriction enzyme genes in newly sequenced genomes, thereby facilitating functional annotation. Active restriction enzyme genes have been identified, and their activity demonstrated biochemically, in the sequenced genomes of Methanocaldococcus jannaschii, Bacillus cereus ATCC 10987 and Methylococcus capsulatus.     

Common patterns in type II restriction enzyme binding sites

Restriction enzymes are among the best studied examples of DNA binding proteins. In order to find general patterns in DNA recognition sites, which may reflect important properties of protein–DNA interaction, we analyse the binding sites of all known type II restriction endonucleases. We find a significantly enhanced GC content and discuss three explanations for this phenomenon. Moreover, we study patterns of nucleotide order in recognition sites. Our analysis reveals a striking accumulation of adjacent purines (R) or pyrimidines (Y). We discuss three possible reasons: RR/YY dinucleotides are characterized by (i) stronger H-bond donor and acceptor clusters, (ii) specific geometrical properties and (iii) a low stacking energy. These features make RR/YY steps particularly accessible for specific protein–DNA interactions. Finally, we show that the recognition sites of type II restriction enzymes are underrepresented in host genomes and in phage genomes.     

Development of primer sets designed for use with the PCR to amplify conserved genes from filamentous ascomycetes.

We constructed nine sets of oligonucleotide primers on the basis of the results of DNA hybridization of cloned genes from Neurospora crassa and Aspergillus nidulans to the genomes of select filamentous ascomycetes and deuteromycetes (with filamentous ascomycete affiliations). Nine sets of primers were designed to amplify segments of DNA that span one or more introns in conserved genes. PCR DNA amplification with the nine primer sets with genomic DNA from ascomycetes, deuteromycetes, basidiomycetes, and plants revealed that five of the primer sets amplified a product only from DNA of the filamentous ascomycetes and deuteromycetes. The five primer sets were constructed from the N. crassa genes for histone 3, histone 4, beta-tubulin, and the plasma membrane ATPase. With these five primer sets, polymorphisms were observed in both the size of and restriction enzyme sites in the amplified products from the filamentous ascomycetes. The primer sets described here may provide useful tools for phylogenetic studies and genome analyses in filamentous ascomycetes and deuteromycetes (with ascomycete affiliations), as well as for the rapid differentiation of fungal species by PCR.     

Rapid and reliable identification of rice genomes by RFLP analysis of PCR-amplified Adh genes.

The rice genus (Oryza L.) consists of 24 species with 10 recognized genome types. With the realization of many useful genes in species of wild rice, continuous efforts have been made to understand their genomic composition and relationships. However, the identification of rice genomes has often been difficult owing to complex morphological variation and formation of allotetraploids. Here we propose a rapid and reliable method for identifying rice genomes based on the restriction sites of PCR-amplified Adh genes. The experimental procedure was as follows: (i) amplify a portion of Adh1 and Adh2 genes with the locus-specific PCR primers; (ii) digest PCR products with restriction enzymes that distinguish different genomes; and (iii) run the digested products on 1.4% agarose gel, and photograph. Using various combinations of restriction digestion of the two Adh genes, all of the rice genomes can be identified.     

Comparison of Methods of Detection of Exceptional Sequences in Prokaryotic Genomes

Many proteins need recognition of specific DNA sequences for functioning. The number of recognition sites and their distribution along the DNA might be of biological importance. For example, the number of restriction sites is often reduced in prokaryotic and phage genomes to decrease the probability of DNA cleavage by restriction endonucleases. We call a sequence an exceptional one if its frequency in a genome significantly differs from one predicted by some mathematical model. An exceptional sequence could be either under- or over-represented, depending on its frequency in comparison with the predicted one. Exceptional sequences could be considered biologically meaningful, for example, as targets of DNA-binding proteins or as parts of abundant repetitive elements. Several methods to predict frequency of a short sequence in a genome, based on actual frequencies of certain its subsequences, are used. The most popular are methods based on Markov chain models. But any rigorous comparison of the methods has not previously been performed. We compared three methods for the prediction of short sequence frequencies: the maximum-order Markov chain model-based method, the method that uses geometric mean of extended Markovian estimates, and the method that utilizes frequencies of all subsequences including discontiguous ones. We applied them to restriction sites in complete genomes of 2500 prokaryotic species and demonstrated that the results depend greatly on the method used: lists of 5% of the most under-represented sites differed by up to 50%. The method designed by Burge and coauthors in 1992, which utilizes all subsequences of the sequence, showed a higher precision than the other two methods both on prokaryotic genomes and randomly generated sequences after computational imitation of selective pressure. We propose this method as the first choice for detection of exceptional sequences in prokaryotic genomes.     

Fighting Fire with Fire: Endogenous Retrovirus Envelopes as Restriction Factors

A considerable portion of vertebrate genomes are made up of endogenous retroviruses (ERVs). While aberrant or uncontrolled ERV expression has been perceived as a potential cause of disease, there is mounting evidence that some ERVs have become integral components of normal host development and physiology. Here, we revisit the longstanding concept that some of the gene products encoded by ERVs and other endogenous viral elements may offer to the host protection against viral infection. Notably, proteins produced from envelope (env) genes have been shown to act as restriction factors against related exogenous retroviruses in chickens, sheep, mice, and cats. Based on the proposed mode of restriction and the domain architecture of known antiretroviral env, we argue that many more env gene-derived restriction factors await discovery in vertebrate genomes, including the human genome.     

Genome majority vote improves gene predictions

Recent studies have noted extensive inconsistencies in gene start sites among orthologous genes in related microbial genomes. Here we provide the first documented evidence that imposing gene start consistency improves the accuracy of gene start-site prediction. We applied an algorithm using a genome majority vote (GMV) scheme to increase the consistency of gene starts among orthologs. We used a set of validated Escherichia coli genes as a standard to quantify accuracy. Results showed that the GMV algorithm can correct hundreds of gene prediction errors in sets of five or ten genomes while introducing few errors. Using a conservative calculation, we project that GMV would resolve many inconsistencies and errors in publicly available microbial gene maps. Our simple and logical solution provides a notable advance toward accurate gene maps.     

Massively parallel characterization of restriction endonucleases

Restriction endonucleases are highly specific in recognizing the particular DNA sequence they act on. However, their activity is affected by sequence context, enzyme concentration and buffer composition. Changes in these factors may lead to either ineffective cleavage at the cognate restriction site or relaxed specificity allowing cleavage of degenerate ‘star’ sites. Additionally, uncharacterized restriction endonucleases and engineered variants present novel activities. Traditionally, restriction endonuclease activity is assayed on simple substrates such as plasmids and synthesized oligonucleotides. We present and use high-throughput Illumina sequencing-based strategies to assay the sequence specificity and flanking sequence preference of restriction endonucleases. The techniques use fragmented DNA from sequenced genomes to quantify restriction endonuclease cleavage on a complex genomic DNA substrate in a single reaction. By mapping millions of restriction site–flanking reads back to the Escherichia coli and Drosophila melanogaster genomes we were able to quantitatively characterize the cognate and star site activity of EcoRI and MfeI and demonstrate genome-wide decreases in star activity with engineered high-fidelity variants EcoRI-HF and MfeI-HF, as well as quantify the influence on MfeI cleavage conferred by flanking nucleotides. The methods presented are readily applicable to all type II restriction endonucleases that cleave both strands of double-stranded DNA.     

Identification of protein binding sites in genomic DNA by two-dimensional gel electrophoresis.

We describe a simple two-dimensional electrophoresis procedure to identify the recognition sites of DNA-binding proteins within large DNA molecules. Using this approach, we have mapped E. coli IHF (Integration Host Factor) binding sites within phage Lambda (48 kb) and phage Mu (39 kb) DNA. We are also able to visualize IHF binding sites in E. coli chromosomal DNA (4,700 kb). We present an extension of this technique using direct amplification by PCR of the isolated restriction fragments, which should permit the cloning of a collection of recognition sequences for DNA binding proteins in complex genomes.     

Dissection of the Salmonella typhimurium genome by use of a Tn5 derivative carrying rare restriction sites.

A polylinker with rare restriction sites was introduced into a mini-Tn5 derivative. These sites include M.XbaI-DpnI (TCTAGATCTAGA), which is rare in most bacterial genomes, SwaI (ATTTAAAT) and PacI (TTAATTAA), which are rare in G+C-rich genomes, NotI (GCGGCCGC) and SfiI (GGCCN5GGCC), which are rare in A+T-rich genomes, and BlnI (CCTAGG), SpeI (ACTAGT), and XbaI (TCTAGA), which are rare in the genomes of many gram-negative bacteria. This Tn5(pfm) (pulsed-field mapping) transposon carries resistance to chloramphenicol and kanamycin to allow selection in a wide variety of background genomes. This Tn5(pfm) was integrated randomly into the Salmonella typhimurium and Serratia marcescens genomes. Integration of the new rare SwaI, PacI, BlnI, SpeI, and XbaI sites was assayed by restriction digestion and pulsed-field gel electrophoresis. Tn5(pfm) constructs could be valuable tools for pulsed-field mapping of gram-negative bacterial genomes by assisting in the production of physical maps and restriction fragment catalogs. For the first applications of a Tn5(pfm), we bisected five of the six largest BlnI fragments in the S. typhimurium genome, bisected the linearized 90-kb pSLT plasmid, and used Tn5(pfm) and Tn10 to trisect the largest BlnI fragment.     

Emergent neutrality in adaptive asexual evolution

In nonrecombining genomes, genetic linkage can be an important evolutionary force. Linkage generates interference interactions, by which simultaneously occurring mutations affect each other's chance of fixation. Here, we develop a comprehensive model of adaptive evolution in linked genomes, which integrates interference interactions between multiple beneficial and deleterious mutations into a unified framework. By an approximate analytical solution, we predict the fixation rates of these mutations, as well as the probabilities of beneficial and deleterious alleles at fixed genomic sites. We find that interference interactions generate a regime of emergent neutrality: all genomic sites with selection coefficients smaller in magnitude than a characteristic threshold have nearly random fixed alleles, and both beneficial and deleterious mutations at these sites have nearly neutral fixation rates. We show that this dynamic limits not only the speed of adaptation, but also a population's degree of adaptation in its current environment. We apply the model to different scenarios: stationary adaptation in a time-dependent environment and approach to equilibrium in a fixed environment. In both cases, the analytical predictions are in good agreement with numerical simulations. Our results suggest that interference can severely compromise biological functions in an adapting population, which sets viability limits on adaptive evolution under linkage.     

Conservation and variation of nucleotide sequences within related bacterial genomes: Escherichia coli strains.

Changes in the patterns produced by annealing restriction endonuclease digests of bacterial genomes with probe deoxyribonucleic acids (DNAs) containing small portions of a bacterial genome provide sensitive indicator of the degree of nucleotide sequence relatedness that exists in localized regions of the genomes of closely related bacteria. We have used five probe DNAs to explore the relatedness of parts of the genomes of six laboratory Escherichi coli strains. A range in in the amount of variability in the positions of restriction enzyme cleavage sites in the selected portions of the genomes was found. Portions of the genome that are believed to be inacative were more variable than portions that contained functional genes: the sites in and near regions of homology to phage lambda DNA in the genome showed the greatest variability. These regions probably represent remnants of cryptic prophages. Variability was assessed pairwise among four of the E. coli strains and ranged from 5 to > 25% base pair substitutions in the lambda-related regions. In contrast, the endonuclease cleavage sites in the trp, tna, lac, thy regions, and one other as-yet-unidentified segment of the genome were more highly conserved. It seems likely that these sites lie in genetic locations that are subject to functional constraints.     

Amplification of plant genomic DNA by Phi29 DNA polymerase for use in physical mapping of the hypermethylated genomic region

Plant genomes contain a heavily methylated region in which cytosines are methylated in both the symmetrical and asymmetrical sequences. The physical mapping of such a hypermethylated region is difficult because many restriction enzymes are sensitive to methylated cytosine residues in their recognition sites. The Phi29 DNA polymerase provides an efficient and representative amplification of the genomic DNA that is methylation-free. Using this amplified genomic DNA, we were able to show that a heavily methylated genomic DNA region becomes amenable to physical mapping with any restriction enzymes. This protocol will be especially useful for analysis of the heavily methylated region of plant genomes.     

Achilles' heel cleavage: creation of rare restriction sites in lambda phage genomes and evaluation of additional operators, repressors and restriction/modification systems

A novel technique for the creation of rare restriction sites was described by Koob et al. [Science 241 (1988) 1084-1086]. This technique, Achilles' heel cleavage (AC), relies on the use of a bound repressor molecule to protect only one of many identical restriction sites from a modification methyltransferase that inactivates all other restriction sites. The technique was applied to a small plasmid and shown to work efficiently with two repressor/operator systems: lac repressor/lacO operator and lambda repressor/lambda oL1 operator. Here, we have extended these results to a lac operator carried by a much larger vector, namely a 44-kb phage lambda construct. In addition, we have evaluated the effect of altering the stability of the lac repressor/lac operator complex by varying both the operator and the repressor. We have also evaluated several more restriction/modification systems (MboI, Dam, MspI and AluI) in addition to HhaI and HaeII used earlier. Finally, we extended the AC technique to a third system, that of the phage 434 repressor and a synthetic 434 operator. From our results we conclude that the AC method should be applicable to the mapping of large genomes and to measuring the strength of operator-repressor interactions. AC could also be applied to identifying and evaluating many different DNA-binding proteins and their sites of action.     

Efficient enumeration of phylogenetically informative substrings.

We study the problem of enumerating substrings that are common amongst genomes that share evolutionary descent. For example, one might want to enumerate all identical (therefore conserved) substrings that are shared between all mammals and not found in non-mammals. Such collection of substrings may be used to identify conserved subsequences or to construct sets of identifying substrings for branches of a phylogenetic tree. For two disjoint sets of genomes on a phylogenetic tree, a substring is called a tag if it is found in all of the genomes of one set and none of the genomes of the other set. We present a near-linear time algorithm that finds all tags in a given phylogeny; and a sublinear space algorithm (at the expense of running time) that is more suited for very large data sets. Under a stochastic model of evolution, we show that a simple process of tag-generation essentially captures all possible ways of generating tags. We use this insight to develop a faster tag discovery algorithm with a small chance of error. However, since tags are not guaranteed to exist in a given data set, we generalize the notion of a tag from a single substring to a set of substrings. We present a linear programming-based approach for finding approximate generalized tag sets. Finally, we use our tag enumeration algorithm to analyze a phylogeny containing 57 whole microbial genomes. We find tags for all nodes in the phylogeny except the root for which we find generalized tag sets.     

Identification of the avian myeloblastosis virus genome. II. Restriction endonuclease analysis of DNA from lambda proviral recombinants and leukemic myeoblast clones.

Two lambda proviral DNA recombinants were characterized with a number of restriction endonucleases. One recombinant contained a complete presumptive avian myeloblastosis virus (AMV) provirus flanked by cellular sequences on either side, and the second recombinant contained 85% of a myeloblastosis-associated virus type 1 (MAV-1)-like provirus with cellular sequences adjacent to the 5' end of the provirus. Comparing the restriction maps for the proviral DNAs contained in each lambda hybrid showed that the putative AMV and MAV-1-like genomes shared identical enzyme sites for 3.6 megadaltons beginning at the 5' termini of the proviruses with respect to viral RNA. Two enzyme sites near the 3'-end of the MAV-1-like provirus were not present in the putative AMV genome. We also examined a number of leukemic myeloblast clones for proviral content and cell-provirus integration sites. The presumptive AMV provirus was present in all the leukemic myeloblast clones regardless of the endogenous proviral content of the target cells or the AMV pseudotype used for conversion. Multiple cellular sites were suitable for integration of the putative AMV genome and the helper genomes. The proviral genomes were all integrated colinearly with respect to linear viral DNA.     

Web-based primer design for single nucleotide polymorphism analysis

The detection of single nucleotide polymorphisms by PCR is necessary for many types of genetic analysis, from mapping genomes to tracking specific mutations. This technique is most commonly used when polymorphisms alter restriction endonuclease recognition sites. Here we describe a web-based program, dCAPS Finder 2.0, that facilitates the design of mismatched PCR primers to create or remove a restriction endonuclease recognition site relative to the polymorphism being analyzed.