Creating New Genes by Plasmid Recombination in Escherichia coli and Bacillus subtilis

Gene shuffling is a way of creating proteins with interesting new characteristics, starting from diverged sequences. We tested an alternative to gene shuffling based on plasmid recombination and found that Bacillus subtilis efficiently recombines sequences with 4% divergence, and Escherichia coli mutS is more appropriate for sequences with 22% divergence.     

Nucleotide sequences of bovine alpha S1- and kappa-casein cDNAs

The nucleotide sequences corresponding to bovine alpha S1- and kappa-casein mRNAs are presented. An unusual alpha S1-casein cDNA has been characterised whose 5' end commences upstream from its putative TATA box. The alpha S1-casein mRNA is compared to rat alpha-casein mRNA and two components of divergence are identified. Firstly, the two sequences have diverged at a high point mutation rate and the rate of amino acid replacement by this mechanism is at least as great as the rate of divergence of any other part of the mRNAs. Secondly, the protein coding sequence has been subjected to several insertion/deletion events, one of which may be an example of exon shuffling . The kappa-casein mRNA sequence verifies the proposition that it has arisen from a different ancestral gene to the other caseins.     

An effective family shuffling method using single-stranded DNA

Family shuffling, which is one of the most powerful techniques for in vitro protein evolution, always involves the problem of reassembling the gene fragments into parental gene sequences, because such a process prevents the formation of chimeric sequences. In order to improve the efficiency of hybrid formation in family shuffling, single-stranded DNAs (ssDNAs) were used as templates. The ssDNAs of two catechol 2,3-dioxygenase genes, nahH and xylE, were prepared, the xylE strand being complementary to the nahH strand. When these ssDNAs were digested by DNase I and reassembled, chimeric genes were obtained at a rate of 14%, which was much higher than the rate of less than 1% obtained by shuffling with double-stranded DNAs. Chimeric catechol 2,3-dioxygenases that were more thermally stable than the parental enzymes, XylE and NahH, were obtained by this ssDNA-based DNA shuffling.     

Integrated Mimicry of B Cell Antibody Mutagenesis Using Yeast Homologous Recombination

Antibody affinity maturation proceeds in vivo via a combination of point mutations, insertions, deletions, and combinatorial shuffling of light chains or portions of the heavy chain, thereby reducing the probability of trapping in local affinity optima in sequence space. In vivo homologous recombination in yeast can be exploited to mimic the broad spectrum of mutational types deployed by B cells, incorporating both receptor revision and receptor editing together with polymerase-directed point mutagenesis. This method was used to effect a 10,000-fold affinity improvement in an anti-peptide single-chain antibody in three rounds of mutagenesis and screening, and a 1,000-fold affinity improvement in an anti-protein single-chain antibody in a single round. When recombinational mutagenesis (CDR or chain shuffling) was directly compared to error-prone PCR, the recombinational approach yielded greater affinity improvement with substantially reduced divergence from germline sequences, demonstrating an advantage of simultaneously testing a broad range of mutational strategies.     

Versatile capacity of shuffled cytochrome P450s for dye production

DNA family shuffling is a relatively new method of directed evolution used to create novel enzymes in order to improve their existing properties or to develop new features. This method of evolution in vitro has one basic requirement: a high similarity of initial parental sequences. Cytochrome P450 enzymes are relatively well conserved in their amino acid sequences. Members of the same family can have more than 40% of sequence identity at the protein level and are therefore good candidates for DNA family shuffling. These xenobiotic-metabolising enzymes have an ability to metabolise a wide range of chemicals and produce a variety of products including blue pigments such as indigo. By applying the specifically designed DNA family shuffling approach, catalytic properties of cytochrome P450 enzymes were further extended in the chimeric progeny to include a new range of blue colour formations. This mini-review evokes the possibility of exploiting directed evolution of cytochrome P450s and the novel enzymes created by DNA family shuffling for the production of new dyes.     

Directed evolution of proteins by exon shuffling

Evolution of eukaryotes is mediated by sexual recombination of parental genomes. Crossovers occur in random, but homologous, positions at a frequency that depends on DNA length. As exons occupy only 1% of the human genome and introns about 24%, by far most of the crossovers occur between exons, rather than inside. The natural process of creating new combinations of exons by intronic recombination is called exon shuffling. Our group is developing in vitro formats for exon shuffling and applying these to the directed evolution of proteins. Based on the splice frame junctions, nine classes of exons and three classes of introns can be distinguished. Splice frame diagrams of natural genes show how the splice frame rules govern exon shuffling. Here, we review various approaches to constructing libraries of exon-shuffled genes. For example, exon shuffling of human pharmaceutical proteins can generate libraries in which all of the sequences are fully human, without the point mutations that raise concerns about immunogenicity.     

Correcting errors in synthetic DNA through consensus shuffling

Although efficient methods exist to assemble synthetic oligonucleotides into genes and genomes, these suffer from the presence of 1–3 random errors/kb of DNA. Here, we introduce a new method termed consensus shuffling and demonstrate its use to significantly reduce random errors in synthetic DNA. In this method, errors are revealed as mismatches by re-hybridization of the population. The DNA is fragmented, and mismatched fragments are removed upon binding to an immobilized mismatch binding protein (MutS). PCR assembly of the remaining fragments yields a new population of full-length sequences enriched for the consensus sequence of the input population. We show that two iterations of consensus shuffling improved a population of synthetic green fluorescent protein (GFPuv) clones from ~60 to >90% fluorescent, and decreased errors 3.5- to 4.3-fold to final values of ~1 error per 3500 bp. In addition, two iterations of consensus shuffling corrected a population of GFPuv clones where all members were non-functional, to a population where 82% of clones were fluorescent. Consensus shuffling should facilitate the rapid and accurate synthesis of long DNA sequences.     

Frequent appearance of novel protein-coding sequences by frameshift translation

Genomic duplication, followed by divergence, contributes to organismal evolution. Several mechanisms, such as exon shuffling and alternative splicing, are responsible for novel gene functions, but they generate homologous domains and do not usually lead to drastic innovation. Major novelties can potentially be introduced by frameshift mutations and this idea can explain the creation of novel proteins. Here, we employ a strategy using simulated protein sequences and identify 470 human and 108 mouse frameshift events that originate new gene segments. No obvious interspecies overlap was observed, suggesting high rates of acquisition of evolutionary events. This inference is supported by a deficiency of TpA dinucleotides in the protein-coding sequences, which decreases the occurrence of translational termination, even on the complementary strand. Increased usage of the TGA codon as the termination signal in newer genes also supports our inference. This suggests that tolerated frameshift changes are a prevalent mechanism for the rapid emergence of new genes and that protein-coding sequences can be derived from existing or ancestral exons rather than from events that result in noncoding sequences becoming exons.     

Chimeric gene library construction by a simple and highly versatile method using recombination-dependent exponential amplification

A simple and efficient method for the construction of chimeric gene libraries termed RDA-PCR (recombination-dependent exponential amplification polymerase chain reaction) was developed by modifying polymerase chain reaction. A chimeric gene library is generated from homologous parental genes with additional primer-annealing sequences at their "heads" and "tails". Two primers ("skew primers") are designed to exclusively anneal to either the heads of maternal genes or the tails of paternal genes. During the RDA-PCR, short annealing/extension periods facilitate homologous recombination. The chimeric sequences can be exponentially amplified to form the chimeric gene library, whereas parental sequences without crossovers are not amplified. As a model, we constructed a chimeric gene library of yellow and green fluorescent protein (yfp and gfp, respectively). The crossover point profile of RDA-PCR clones was compared with those obtained by (modified) family shuffling. PCR restriction fragment polymorphism (PCR-RFLP) analysis of the RDA-PCR clones showed a high content of chimeric genes in the library, whereas family shuffling required the modification using skew primers for selective enrichment of chimeric sequences. PCR-RFLP analysis also indicated that the crossover points of RDA-PCR chimeras were distributed over the entire protein-coding region. Moreover, as few as 2 bp of the continual identity of nucleotides were found at the crossover points at high frequency (30% of the tested clones), suggesting that RDA-PCR resulted in a higher diversity in crossover points than family shuffling.     

Digestion-after-shuffling: a method for combinatorial library construction

We have developed a convenient method for family shuffling of amino acid sequences, termed digestion-after-shuffling. After DNA shuffling of homologous genes, plasmid mixture is extracted from a library and used for several double digestions with restriction enzymes. For each double digestion, two restriction enzymes are selected, corresponding to the single restriction sites of different parental genes. After digestions, fragments with expected sizes are obtained by gel purification and religated to construct recombinant plasmids. Thus, the obtained genes should be chimeras and have at least two restriction sites originating from different parental sequences.     

Prospecting metagenomic enzyme subfamily genes for DNA family shuffling by a novel PCR-based approach

DNA family shuffling is a powerful method for enzyme engineering, which utilizes recombination of naturally occurring functional diversity to accelerate laboratory-directed evolution. However, the use of this technique has been hindered by the scarcity of family genes with the required level of sequence identity in the genome database. We describe here a strategy for collecting metagenomic homologous genes for DNA shuffling from environmental samples by truncated metagenomic gene-specific PCR (TMGS-PCR). Using identified metagenomic gene-specific primers, twenty-three 921-bp truncated lipase gene fragments, which shared 64-99% identity with each other and formed a distinct subfamily of lipases, were retrieved from 60 metagenomic samples. These lipase genes were shuffled, and selected active clones were characterized. The chimeric clones show extensive functional and genetic diversity, as demonstrated by functional characterization and sequence analysis. Our results indicate that homologous sequences of genes captured by TMGS-PCR can be used as suitable genetic material for DNA family shuffling with broad applications in enzyme engineering.     

DNA shuffling method for generating estrogen receptor alpha and beta chimeras in yeast

To facilitate our study of the molecular basis for the estrogen receptor (ER) subtype selectivity of novel ligands, we used DNA shuffling to construct chimeric ERs having ligand binding domains derived from both ER alpha and ER beta. The efficiency of chimera generation was low with traditional DNA shuffling protocols. Furthermore, ER ligand binding domain sequences lack convenient restriction sites for introducing chimeric ligand binding domain sequences into expression vectors. To overcome these problems, we developed a modified strategy whereby chimeric sequences were exclusively amplified from among the reassembled products from DNA shuffling using a special pair of PCR primers whose 3' ends specifically match the alpha and beta sequences, respectively, and whose 5' ends match sequences outside the ER beta ligand binding domain. When chimeric ligand binding domain DNA sequences, amplified with these primers, were co-transformed into a yeast strain with a linearized expression vector for ER beta, an active expression vector was produced by homologous recombination. Twenty-two different crossover sites were found; most occurred when there was a stretch of eight or more identical base pairs in both sequences, and many were concentrated in the regions important for studying ligand binding and transactivation. This method should prove to be useful for generating chimeric gene products from parent templates that share relatively low sequence identity.     

过氧化氢酶基因定向进化文库的构建及其评价

采用易错聚合酶链反应和DNA改组技术构建野生型梅花鹿过氧化氢酶(CAT)基因的突变体文库,并随机对两种方法所得产物各5个样品做序列测定。序列分析结果表明突变率分别为0.329%和27.58%,易错聚合酶链反应体系的错配率可以比普通PCR体系提高约10倍,DNA改组的突变率则更高,但是难以避免由于突变率太高造成的目的基因无法正确翻译这一情况。另外,应用邻接法(neighbor-joining, NJ)对随机选择的过氧化氢酶基因突变体序列和野生型序列做核酸和蛋白质序列的NJ进化树,进化关系与突变率分析基本一致。     

Unraveling the sequence dynamics of the formation of genus-specific satellite DNAs in the family solanaceae

Tandemly repeated DNAs, referred to as satellite DNAs, often occur in a genome in a genus-specific manner. However, the mechanisms for generation and evolution for these sequences are largely unknown because of the uncertain origins of the satellite DNAs. We found highly divergent genus-specific satellite DNAs that showed sequence similarity with genus-specific intergenic spacers (IGSs) in the family Solanaceae, which includes the genera Nicotiana, Solanum and Capsicum. The conserved position of the IGS between 25S and 18S rDNA facilitates comparison of IGS sequences across genera, even in the presence of very low sequence similarity. Sequence comparison of IGS may elucidate the procedure of the genesis of complex monomer units of the satellite DNAs. Within the IGS of Capsicum species, base substitutions and copy number variation of subrepeat monomers were causes of monomer divergence in IGS sequences. At the level of inter-generic IGS sequences of the family Solanaceae, however, genus-specific motif selection, motif shuffling between subrepeats and differential amplification among motifs were involved in formation of genus-specific IGS. Therefore, the genus-specific satellite DNAs in Solanaceae plants can be generated from differentially organized repeat monomers of the IGS rather than by accumulation of mutations from pre-existent satellite DNAs.     

Functional promoter analysis using an approach based on an in vitro evolution strategy

In vitro evolution imitates the natural evolution of genes and has been very successfully applied to the modification of coding sequences, but it has not yet been applied to promoter sequences. We propose an alternative method for functional promoter analysis by applying an in vitro evolution scheme consisting of rounds of error-prone PCR, followed by DNA shuffling and selection of mutant promoter activities. We modified the activity in embryogenic sugarcane cells of the promoter region of the "Goldfinger" isolate of banana streak virus and obtained mutant promoter sequences that showed an average mutation rate of 2.5% after applying one round of error-prone PCR and DNA shuffling. Selection and sequencing of promoter sequences with decreased or unaltered activity allowed us to rapidly map the position of one cis-acting element that influenced promoter activity in embryogenic sugarcane cells and to discover neutral mutations that did not affect promoter function. The "selective-shotgun" approach of this promoter analysis method immediately after the promoter boundaries have been defined by 5' deletion analysis dramatically reduces the labor associated with traditional "linker-scanning" deletion analysis to reveal the position of functional promoter domains. Furthermore, this method allows the entire promoter to be investigated at once, rather than selected domains or nucleotides, increasing the prospect of identifying interacting promoter regions.     

indel-Seq-Gen: a new protein family simulator incorporating domains, motifs, and indels

Reconstructing the evolutionary history of protein sequences will provide a better understanding of divergence mechanisms of protein superfamilies and their functions. Long-term protein evolution often includes dynamic changes such as insertion, deletion, and domain shuffling. Such dynamic changes make reconstructing protein sequence evolution difficult and affect the accuracy of molecular evolutionary methods, such as multiple alignments and phylogenetic methods. Unfortunately, currently available simulation methods are not sufficiently flexible and do not allow biologically realistic dynamic protein sequence evolution. We introduce a new method, indel-Seq-Gen (iSG), that can simulate realistic evolutionary processes of protein sequences with insertions and deletions (indels). Unlike other simulation methods, iSG allows the user to simulate multiple subsequences according to different evolutionary parameters, which is necessary for generating realistic protein families with multiple domains. iSG tracks all evolutionary events including indels and outputs the "true" multiple alignment of the simulated sequences. iSG can also generate a larger sequence space by allowing the use of multiple related root sequences. With all these functions, iSG can be used to test the accuracy of, for example, multiple alignment methods, phylogenetic methods, evolutionary hypotheses, ancestral protein reconstruction methods, and protein family classification methods. We empirically evaluated the performance of iSG against currently available methods by simulating the evolution of the G protein-coupled receptor and lipocalin protein families. We examined their true multiple alignments, reconstruction of the transmembrane regions and beta-strands, and the results of similarity search against a protein database using the simulated sequences. We also presented an example of using iSG for examining how phylogenetic reconstruction is affected by high indel rates.     

A complete mitochondrial genome sequence of Ogura-type male-sterile cytoplasm and its comparative analysis with that of normal cytoplasm in radish (Raphanus sativus L.)

ABSTRACT: BACKGROUND: Plant mitochondrial genome has unique features such as large size, frequent recombination and incorporation of foreign DNA. Cytoplasmic male sterility (CMS) is caused by rearrangement of the mitochondrial genome, and a novel chimeric open reading frame (ORF) created by shuffling of endogenous sequences is often responsible for CMS. The Ogura-type male-sterile cytoplasm is one of the most extensively studied cytoplasms in Brassicaceae. Although the gene orf138 has been isolated as a determinant of Ogura-type CMS, no homologous sequence to orf138 has been found in public databases. Therefore, how orf138 sequence was created is a mystery. In this study, we determined the complete nucleotide sequence of two radish mitochondrial genomes, namely, Ogura- and normal-type genomes, and analyzed them to reveal the origin of the gene orf138. RESULTS: Ogura- and normal-type mitochondrial genomes were assembled to 258,426-bp and 244,036-bp circular sequences, respectively. Normal-type mitochondrial genome contained 33 protein-coding and three rRNA genes, which are well conserved with the reported mitochondrial genome of rapeseed. Ogura-type genomes contained same genes and additional atp9. As for tRNA, normal-type contained 17 tRNAs, while Ogura type contained 17 tRNAs and one additional trnfM. The gene orf138 was specific to Ogura-type mitochondrial genome, and no sequence homologous to it was found in normal-type genome. Comparative analysis of the two genomes revealed that radish mitochondrial genome consists of 11 syntenic regions (length >3kb, similarity >99.9%). It was shown that short repeats and overlapped repeats present in the edge of syntenic regions were involved in recombination events during evolution to interconvert two types of mitochondrial genome. Ogura-type mitochondrial genome has four unique regions (2,803 bp, 1,601 bp, 451 bp and 15,255 bp in size) that are non-syntenic to normal-type genome, and the gene orf138 was found to be located at the edge of the largest unique region. Blast analysis performed to assign the unique regions showed that about 80% of the region was covered by short homologous sequences to the mitochondrial sequences of normal-type radish or other reported Brassicaceae species, although no homology was found for the remaining 20% of sequences. CONCLUSIONS: Ogura-type mitochondrial genome was highly rearranged compared with the normal-type genome by recombination through one large repeat and multiple short repeats. The rearrangement has produced four unique regions in Ogura-type mitochondrial genome, and most of the unique regions are composed of known Brassicaceae mitochondrial sequences. This suggests that the regions unique to the Ogura-type genome were generated by integration and shuffling of pre-existing mitochondrial sequences during the evolution of Brassicaceae, and novel genes such as orf138 could have been created by the shuffling process of mitochondrial genome.     

Divergence of unique DNA sequences of bivalve mollusks from Mytilinae subfamily (Bivalvia mytilidae)

The DNA-DNA hybridization method was used to determine the divergence degree of unique sequences of five bivalvia species belonging to the Mytilinae subfamily. The matrix of delta Tm values for heteroduplexes of unique sequences was found which made it possible to define three phylogenetic branches within the subfamily. Under non-stringent hybridization conditions (55 degrees C, 0.5M PB) the divergence between species of any two branches was about 14% of nucleotide substitutions. The hybridization of [3H] unique sequences of Mytilinae with DNA fragments of Modiolus modiolus, a representative of the closest relative Modiolinae subfamily, showed that the divergence rate of unique sequences in two phylogenetic lines of the Mytilus genusis higher than in the line of the Crenomytilus genus I0.35-0.23% and 0.1% of nucleotide substitutions per one million years, respectively). According to the matrix of delta Tm values for five species of Mytilinae and Modiolus modiolus, a phylogenetic tree was built reflecting the differences between the divergence rates in different branches.