ROLL: a method of preparation of gene-specific oligonucleotide libraries

The selection of nucleic acid sequences capable of specifically and efficiently hybridizing to target sequences is crucial to the success of many applications, including microarrays, PCR and other amplification procedures, antisense inhibition, ribozyme-mediated cleavage, and RNA interference (RNAi). Methods of selection using nucleotide sequence libraries have several advantages over rational approaches using defined sequences. However, the high complexity of completely random (degenerate) libraries and their high toxicity in cell-based assays make their use in many applications impractical. Gene-specific oligonucleotide libraries, which contain all possible sequences of a certain length occurring within a given gene, have much lower complexity and, thus, can significantly simplify and accelerate sequence screening. Here, we describe a new method for the preparation of gene-specific libraries using the ligation of randomized oligonucleotide probes hybridized adjacently on target polynucleotide templates followed by PCR amplification. We call this method random oligonucleotide ligated libraries (ROLL).     

In vitro selection using a dual RNA library that allows primerless selection

High affinity target-binding aptamers are identified from random oligonucleotide libraries by an in vitro selection process called Systematic Evolution of Ligands by EXponential enrichment (SELEX). Since the SELEX process includes a PCR amplification step the randomized region of the oligonucleotide libraries need to be flanked by two fixed primer binding sequences. These primer binding sites are often difficult to truncate because they may be necessary to maintain the structure of the aptamer or may even be part of the target binding motif. We designed a novel type of RNA library that carries fixed sequences which constrain the oligonucleotides into a partly double-stranded structure, thereby minimizing the risk that the primer binding sequences become part of the target-binding motif. Moreover, the specific design of the library including the use of tandem RNA Polymerase promoters allows the selection of oligonucleotides without any primer binding sequences. The library was used to select aptamers to the mirror-image peptide of ghrelin. Ghrelin is a potent stimulator of growth-hormone release and food intake. After selection, the identified aptamer sequences were directly synthesized in their mirror-image configuration. The final 44 nt-Spiegelmer, named NOX-B11-3, blocks ghrelin action in a cell culture assay displaying an IC₅₀ of 4.5 nM at 37°C.     

Generation of amplifiable genome-specific oligonucleotide probes and libraries

Here we describe a process for the generation of oligonucleotide libraries representative of a given nucleic acid. Starting from at random pool of DNA oligonucleotides, the technique selects only those that hybridize to the nucleic acid template. This selection yields a highly specific library that represents an oligonucleotide image of the chosen template. The novel quality of this approach is the generation of amplifiable oligonucleotide probes that are of unique length and are easily subjected to differential selection. Here we apply this technique to produce different genomic oligonucleotide libraries and show that these genomic oligonucleotide libraries do not cross-hybridize. Differential selection of these genomic oligonucleotide libraries produces oligonucleotides that can be used in the identification, characterzation, and isolation of nucleic acids.     

Random oligonucleotide mutagenesis: application to a large protein coding sequence of a major histocompatibility complex class I gene, H-2DP.

We have used random oligonucleotide mutagenesis (or saturation mutagenesis) to create a library of point mutations in the alpha 1 protein domain of a Major Histocompatibility Complex (MHC) molecule. This protein domain is critical for T cell and B cell recognition. We altered the MHC class I H-2DP gene sequence such that synthetic mutant alpha 1 exons (270 bp of coding sequence), which contain mutations identified by sequence analysis, can replace the wild type alpha 1 exon. The synthetic exons were constructed from twelve overlapping oligonucleotides which contained an average of 1.3 random point mutations per intact exon. DNA sequence analysis of mutant alpha 1 exons has shown a point mutant distribution that fits a Poisson distribution, and thus emphasizes the utility of this mutagenesis technique to "scan" a large protein sequence for important mutations. We report our use of saturation mutagenesis to scan an entire exon of the H-2DP gene, a cassette strategy to replace the wild type alpha 1 exon with individual mutant alpha 1 exons, and analysis of mutant molecules expressed on the surface of transfected mouse L cells.     

Determination of optimal sites of antisense oligonucleotide cleavage within TNFalpha mRNA

Antisense oligonucleotides provide a powerful tool in order to determine the consequences of the reduced expression of a selected target gene and may include target validation and therapeutic applications. Methods of predicting optimum antisense sites are not always effective. We have compared the efficacy of antisense oligonucleotides, which were selected in vitro using random combinatorial oligonucleotide libraries of differing length and complexity, upon putative target sites within TNFα mRNA. The relationship of specific target site accessibility and oligonucleotide efficacy with respect to these parameters proved to be complex. Modification of the length of the recognition sequence of the oligonucleotide library illustrated that independent target sites demonstrated a preference for antisense oligonucleotides of a defined and independent optimal length. The efficacy of antisense oligonucleotide sequences selected in vitro paralleled that observed in phorbol 12-myristate 13-acetate (PMA)-activated U937 cells. The application of methylphosphonate:phosphodiester chimaeric oligonucleotides to U937 cells reduced mRNA levels to up to 19.8% that of the untreated cell population. This approach provides a predictive means to profile any mRNA of known sequence with respect to the identification and optimisation of sites accessible to antisense oligonucleotide activity.     

Method to protect a targeted amino acid residue during random mutagenesis

To generate a random mutant library that is free from mutation at a particular amino acid residue, we replace the codon of interest with a detachable, short DNA sequence containing a BsaXI recognition site. After PCR mutagenesis, this sequence is removed and intramolecular ligation of the sequences flanking the insert regenerates the gene. The three-base cohesive ends for ligation correspond to the codon for the targeted residue and any sequences with mutations at this site will fail to ligate. As a result, only the variants that are free from mutation at this site are in the proper reading frame. In a random library of C₃₀ carotenoid synthase CrtM, this method was used to exclude readily accessible mutations at position F26, which confer C₄₀ synthase function. This enabled us to identify two additional mutations, W38C and E180G, which confer the same phenotype but are present in the random library at much lower frequencies.     

Constraining protein sequence space: four amino acid alphabets are sufficient to recapitulate lambda repressor multimerization

Nucleic acid polymers selected from random sequence space constitute an enormous array of catalytic, diagnostic and therapeutic molecules. Despite the fact that proteins are robust polymers with far greater chemical and physical diversity, success in unlocking protein sequence space remains elusive. We have devised a combinatorial strategy for accessing nucleic acid sequence space corresponding to proteins comprising selected amino acid alphabets. Using the SynthOMIC approach (synthesis of ORFs by multimerizing in-frame codons), representative libraries comprising four amino acid alphabets were fused in-frame to the lambda repressor DNA-binding domain to provide an in vivo selection for self-interacting proteins that re-constitute lambda repressor function. The frequency of self-interactors as a function of amino acid composition ranged over five orders of magnitude, from ∼6% of clones in a library comprising the amino acid residues LARE to ∼0.6 in 10⁶ in the MASH library. Sequence motifs were evident by inspection in many cases, and individual clones from each library presented substantial sequence identity with translated proteins by BLAST analysis. We posit that the SynthOMIC approach represents a powerful strategy for creating combinatorial libraries of open reading frames that distils protein sequence space on the basis of three inherent properties: it supports the use of selected amino acid alphabets, eliminates redundant sequences and locally constrains amino acids.     

Calculating complexity of large randomized libraries

Randomized libraries are increasingly popular in protein engineering and other biomedical research fields. Statistics of the libraries are useful to guide and evaluate randomized library construction. Previous works only give the mean of the number of unique sequences in the library, and they can only handle equal molar ratio of the four nucleotides at a small number of mutation sites. We derive formulas to calculate the mean and variance of the number of unique sequences in libraries generated by cassette mutagenesis with mixtures of arbitrary nucleotide ratios. Computer program was developed which utilizes arbitrary numerical precision software package to calculate the statistics of large libraries. The statistics of library with mutations in more than 20 amino acids can be calculated easily. Results show that the nucleotide ratios have significant effects on these statistics. The more skewed the ratio, the larger the library size is needed to obtain the same expected number of unique sequences. The program is freely available at http://graphics.med.yale.edu/cgi-bin/lib_comp.pl.     

Cloning of the thermostable phytase gene (phy) from Bacillus sp. DS11 and its overexpression in Escherichia coli

Phytase hydrolyzes phytate to release inorganic phosphate, which would decrease the addition of phosphorus to feedstuffs for monogastric animals and thus reduce environmental pollution. The gene encoding phytase from Bacillus sp. DS11 was cloned in Escherichia coli and its sequence determined. A 560-bp DNA fragment was used as a probe to screen the genomic library. It was obtained through PCR of Bacillus sp. DS11 chromosomal DNA and two oligonucleotide primers based on N-terminal amino acid sequences of the purified protein and the cyanogen bromide-cleaved 21-kDa fragment. The phy cloned was encoded by a 2.2-kb fragment. This gene comprises 1152 nucleotides and encodes a polypeptide of 383 amino acids with a deduced molecular mass of 41 808 Da. Phytase was produced to 20% content of total soluble proteins in E. coli BL21 (DE3) using the pET22b(+) vector with the inducible T7 promoter. This is the first nucleic sequence report on phytase from a bacterial strain.     

Low free energy cost of very long loop insertions in proteins

Long insertions into a loop of a folded host protein are expected to have destabilizing effects because of the entropic cost associated with loop closure unless the inserted sequence adopts a folded structure with amino- and carboxy-termini in close proximity. A loop entropy reduction screen based on this concept was used in an attempt to retrieve folded sequences from random sequence libraries. A library of long random sequences was inserted into a loop of the SH2 domain, displayed on the surface of M13 phage, and the inserted sequences that did not disrupt SH2 function were retrieved by panning using beads coated with a phosphotyrosine containing SH2 peptide ligand. Two sequences of a library of 2 x 10(8) sequences were isolated after multiple rounds of panning, and were found to have recovery levels similar to the wild-type SH2 domain and to be relatively intolerant to further mutation in PCR mutagenesis experiments. Surprisingly, although these inserted sequences exhibited little nonrandom structure, they do not significantly destabilize the host SH2 domain. Additional insertion variants recovered at lower levels in the panning experiments were also found to have a minimal effect on the stability and peptide-binding function of the SH2 domain. The additional level of selection present in the panning experiments is likely to involve in vivo folding and assembly, as there was a rough correlation between recovery levels in the phage-panning experiments and protein solubility. The finding that loop insertions of 60-80 amino acids have minimal effects on SH2 domain stability suggests that the free energy cost of inserting long loops may be considerably less than polymer theory estimates based on the entropic cost of loop closure, and, hence, that loop insertion may have provided an evolutionary route to multidomain protein structures.     

东方田鼠部分BAC文库的微卫星筛选

目的 从东方田鼠的部分BAC文库中筛选微卫星.方法 应用非放射性的菌落杂交方法和磁珠富集法从东方田鼠的BAC文库中筛选高质量的微卫星标记.结果 以地高辛标记的寡聚核苷酸(CA)20为探针,通过菌落杂交法从136个东方田鼠BAC克隆中筛选出杂交信号最强的20个阳性克隆.再将这20个阳性克隆分别通过链霉亲和素磁珠法构建亚克隆文库,从中选取400个经PCR鉴定为阳性的亚克隆进一步测序分析,共得到220个微卫星序列,阳性率55％.选取重复次数高,侧翼序列完整的微卫星序列设计74对引物,共有35对引物能扩增出清晰的条带,其中16对引物具有多态性.结论 成功且高效地从阳性BAC克隆中筛选出微卫星序列,这些微卫星和阳性BAC克隆可用于后续的定位研究.     

Selection of zinc fingers that bind single-stranded telomeric DNA in the G-quadruplex conformation

There is considerable interest in molecules that bind to telomeric DNA sequences and G-quadruplexes with specificity. Such molecules would be useful to test hypotheses for telomere length regulation, and may have therapeutic potential. The versatility and modular nature of the zinc finger motif makes it an ideal candidate for engineering G-quadruplex-binding proteins. Phage display technology has previously been widely used to screen libraries of zinc fingers for binding to novel duplex DNA sequences. In this study, a three-finger library has been screened for clones that bind to an oligonucleotide containing the human telomeric repeat sequence folded in the G-quadruplex conformation. The selected clones show a strong amino acid consensus, suggesting analogous modes of binding. Binding was found to be both sequence dependent and structure specific. This is the first example of an engineered protein that binds to G-quadruplex DNA, and represents a new type of binding interaction for a zinc finger protein.     

Use of short oligonucleotides to screen cosmid libraries for clones containing G/C-rich sequences

We have developed a method to identify clones containing recognition sequences for enzymes that cut mammalian genomes infrequently by direct screening of genomic libraries. The degenerate oligonucleotide NNGCGGCCGCNN, in which the internal 8 bases correspond to the recognition sequence of Not I, was used to screen a cosmid library, and it led to a greater than 10-fold enrichment in the number of clones containing Not I sites. This technique permits the efficient identification of sufficient clones from a chromosome-specific library to allow the construction of a complete pulsed-field map of that chromosome and to assist in finding genes in genomic DNA.     

A method for construction of long randomized open reading frames and polypeptides.

A method is presented for construction of randomized open reading frame sequences (ORFs) and gene libraries containing them. The building blocks for the ORFs were 75 bp long DNA fragments generated by cloning sequences from a single synthetic oligonucleotide preparation by bridge mutagenesis. The fragments had the property that, regardless of their orientation in the ligated product, the ORF of the construct was maintained. The heterogeneity of the ORFs resulted from the random ligation of 2000 different DNA fragments. The randomized ORFs were cloned downstream from the lac promoter in a multicopy plasmid in Escherichia coli. To test the method, a library of 10(6) clones was constructed.     

Design of nucleic acid sequences for DNA computing based on a thermodynamic approach

We have developed an algorithm for designing multiple sequences of nucleic acids that have a uniform melting temperature between the sequence and its complement and that do not hybridize non-specifically with each other based on the minimum free energy (ΔG_min). Sequences that satisfy these constraints can be utilized in computations, various engineering applications such as microarrays, and nano-fabrications. Our algorithm is a random generate-and-test algorithm: it generates a candidate sequence randomly and tests whether the sequence satisfies the constraints. The novelty of our algorithm is that the filtering method uses a greedy search to calculate ΔG_min. This effectively excludes inappropriate sequences before ΔG_min is calculated, thereby reducing computation time drastically when compared with an algorithm without the filtering. Experimental results in silico showed the superiority of the greedy search over the traditional approach based on the hamming distance. In addition, experimental results in vitro demonstrated that the experimental free energy (ΔG_exp) of 126 sequences correlated well with ΔG_min (|R| = 0.90) than with the hamming distance (|R| = 0.80). These results validate the rationality of a thermodynamic approach. We implemented our algorithm in a graphic user interface-based program written in Java.     

Prediction of common folding structures of homologous RNAs.

We have developed an algorithm and a computer program for simultaneously folding homologous RNA sequences. Given an alignment of M homologous sequences of length N, the program performs phylogenetic comparative analysis and predicts a common secondary structure conserved in the sequences. When the structure is not uniquely determined, it infers multiple structures which appear most plausible. This method is superior to energy minimization methods in the sense that it is not sensitive to point mutation of a sequence. It is also superior to usual phylogenetic comparative methods in that it does not require manual scrutiny for covariation or secondary structures. The most plausible 1-5 structures are produced in O(MN2 + N3) time and O(N2) space, which are the same requirements as those of widely used dynamic programs based on energy minimization for folding a single sequence. This is the first algorithm probably practical both in terms of time and space for finding secondary structures of homologous RNA sequences. The algorithm has been implemented in C on a Sun SparcStation, and has been verified by testing on tRNAs, 5S rRNAs, 16S rRNAs, TAR RNAs of human immunodeficiency virus type 1 (HIV-1), and RRE RNAs of HIV-1. We have also applied the program to cis-acting packaging sequences of HIV-1, for which no generally accepted structures yet exist, and propose potentially stable structures. Simulation of the program with random sequences with the same base composition and the same degree of similarity as the above sequences shows that structures common to homologous sequences are very unlikely to occur by chance in random sequences.     

Production of complex nucleic acid libraries using highly parallel in situ oligonucleotide synthesis

Generation of complex libraries of defined nucleic acid sequences can greatly aid the functional analysis of protein and gene function. Previously, such studies relied either on individually synthesized oligonucleotides or on cellular nucleic acids as the starting material. As each method has disadvantages, we have developed a rapid and cost-effective alternative for construction of small-fragment DNA libraries of defined sequences. This approach uses in situ microarray DNA synthesis for generation of complex oligonucleotide populations. These populations can be recovered and either used directly or immortalized by cloning. From a single microarray, a library containing thousands of unique sequences can be generated. As an example of the potential applications of this technology, we have tested the approach for the production of plasmids encoding short hairpin RNAs (shRNAs) targeting numerous human and mouse genes. We achieved high-fidelity clone retrieval with a uniform representation of intended library sequences.     

Rapid cloning and bioinformatic analysis of spinach Y chromosome-specific EST sequences

The genome of spinach single chromosome complement is about 1000 Mbp, which is the model material to study the molecular mechanisms of plant sex differentiation. The cytological study showed that the biggest spinach chromosome (chromosome 1) was taken as spinach sex chromosome. It had three alleles of sex-related X, X ^m and Y. Many researchers have been trying to clone the sex-determining genes and investigated the molecular mechanism of spinach sex differentiation. However, there are no successful cloned reports about these genes. A new technology combining chromosome microdissection with hybridization-specific amplification (HSA) was adopted. The spinach Y chromosome degenerate oligonucleotide primed-PCR (DOP-PCR) products were hybridized with cDNA of the male spinach flowers in florescence. The female spinach genome was taken as blocker and cDNA library specifically expressed in Y chromosome was constructed. Moreover, expressed sequence tag (EST) sequences in cDNA library were cloned, sequenced and bioinformatics was analysed. There were 63 valid EST sequences obtained in this study. The fragment size was between 53 and 486 bp. BLASTn homologous alignment indicated that 12 EST sequences had homologous sequences of nucleic acids, the rest were new sequences. BLASTx homologous alignment indicated that 16 EST sequences had homologous protein-encoding nucleic acid sequence. The spinach Y chromosome-specific EST sequences laid the foundation for cloning the functional genes, specifically expressed in spinach Y chromosome. Meanwhile, the establishment of the technology system in the research provided a reference for rapid cloning of other biological sex chromosome-specific EST sequences.