HAPIscreen,a method for high-throughput aptamer identification

Background

Aptamers are oligonucleotides displaying specific binding properties for a predetermined target. They are selected from libraries of randomly synthesized candidates through an in vitro selection process termed SELEX (Systematic Evolution of Ligands by EXponential enrichment) alternating selection and amplification steps. SELEX is followed by cloning and sequencing of the enriched pool of oligonucleotides to enable comparison of the selected sequences. The most represented candidates are then synthesized and their binding properties are individually evaluated thus leading to the identification of aptamers. These post-selection steps are time consuming and introduce a bias to the expense of poorly amplified binders that might be of high affinity and are consequently underrepresented. A method that would circumvent these limitations would be highly valuable.

Results

We describe a novel homogeneous solution-based method for screening large populations of oligonucleotide candidates generated from SELEX. This approach, based on the AlphaScreen^® technology, is carried out on the exclusive basis of the binding properties of the selected candidates without the needs of performing a priori sequencing. It therefore enables the functional identification of high affinity aptamers. We validated the HAPIscreen (High throughput APtamer Identification screen) methodology using aptamers targeted to RNA hairpins, previously identified in our laboratory. We then screened pools of candidates issued from SELEX rounds in a 384 well microplate format and identify new RNA aptamers to pre-microRNAs.

Conclusions

HAPIscreen, an Alphascreen^®-based methodology for the identification of aptamers is faster and less biased than current procedures based on sequence comparison of selected oligonucleotides and sampling binding properties of few individuals. Moreover this methodology allows for screening larger number of candidates. Used here for selecting anti-premiR aptamers, HAPIscreen can be adapted to any type of tagged target and is fully amenable to automation.

     

High-throughput identification of genetic markers using representational oligonucleotide microarray analysis

We describe a novel approach for high-throughput development of genetic markers using representational oligonucleotide microarray analysis. We test the performance of the method in sugar beet (Beta vulgaris L.) as a model for crop plants with little sequence information available. Genomic representations of both parents of a mapping population were hybridized on microarrays containing in total 146,554 custom made oligonucleotides based on sugar beet bacterial artificial chromosome (BAC) end sequences and expressed sequence tags (ESTs). Oligonucleotides showing a signal with one parental line only, were selected as potential marker candidates and placed onto an array, designed for genotyping of 184 F₂ individuals from the mapping population. Utilizing known co-dominant anchor markers we obtained 511 new dominant markers (392 derived from BAC end sequences, and 119 from ESTs) distributed over all nine sugar beet linkage groups and calculated genetic maps. Further improvements for large-scale application of the approach are discussed and its feasibility for the cost-effective and flexible generation of genetic markers is presented.     

An improved hypergeometric probability method for identification of functionally linked proteins using phylogenetic profiles

Predicting functions of proteins and alternatively spliced isoforms encoded in a genome is one of the important applications of bioinformatics in the post-genome era. Due to the practical limitation of experimental characterization of all proteins encoded in a genome using biochemical studies, bioinformatics methods provide powerful tools for function annotation and prediction. These methods also help minimize the growing sequence-to-function gap. Phylogenetic profiling is a bioinformatics approach to identify the influence of a trait across species and can be employed to infer the evolutionary history of proteins encoded in genomes. Here we propose an improved phylogenetic profile-based method which considers the co-evolution of the reference genome to derive the basic similarity measure, the background phylogeny of target genomes for profile generation and assigning weights to target genomes. The ordering of genomes and the runs of consecutive matches between the proteins were used to define phylogenetic relationships in the approach. We used Escherichia coli K12 genome as the reference genome and its 4195 proteins were used in the current analysis. We compared our approach with two existing methods and our initial results show that the predictions have outperformed two of the existing approaches. In addition, we have validated our method using a targeted protein-protein interaction network derived from protein-protein interaction database STRING. Our preliminary results indicates that improvement in function prediction can be attained by using coevolution-based similarity measures and the runs on to the same scale instead of computing them in different scales. Our method can be applied at the whole-genome level for annotating hypothetical proteins from prokaryotic genomes.     

An improved method for the determination of oligonucleotide chain length using phosphodiesterase hydrolysis and ion-pair high-performance liquid chromatography

An improved method for the determination of chain length, leading and/or terminating nucleotide, and base composition of the deoxyribo- and ribooligonucleotides was developed, based upon the phosphodiesterase-catalyzed hydrolysis of the oligonucleotide, followed by ion-pair (hetaeric) high-performance liquid chromatography using tetrabutyl ammonium phosphate as the hetaeron. Chains containing terminal phosphate groups were first dephosphorylated using alkaline phosphatase. Total analysis times involved a 1.5-h phosphodiesterase incubation, followed by a 40-min chromatographic separation. Minimal or no sample preparation is involved. An analysis of the propagation of errors indicated that chains of up to 36 nucleotides in length could be analyzed with less than a 0.5-base-unit error. Ultimate sensitivity will depend upon the particular high-performance liquid chromatography system, and will be limited by the detection limits for a single nucleoside which, for the described system, was on the order of 50 pmol of oligomer.     

An improved single-cell cDNA amplification method for efficient high-density oligonucleotide microarray analysis

A systems-level understanding of a small but essential population of cells in development or adulthood (e.g. somatic stem cells) requires accurate quantitative monitoring of genome-wide gene expression, ideally from single cells. We report here a strategy to globally amplify mRNAs from single cells for highly quantitative high-density oligonucleotide microarray analysis that combines a small number of directional PCR cycles with subsequent linear amplification. Using this strategy, both the representation of gene expression profiles and reproducibility between individual experiments are unambiguously improved from the original method, along with high coverage and accuracy. The immediate application of this method to single cells in the undifferentiated inner cell masses of mouse blastocysts at embryonic day (E) 3.5 revealed the presence of two populations of cells, one with primitive endoderm (PE) expression and the other with pluripotent epiblast-like gene expression. The genes expressed differentially between these two populations were well preserved in morphologically differentiated PE and epiblast in the embryos one day later (E4.5), demonstrating that the method successfully detects subtle but essential differences in gene expression at the single-cell level among seemingly homogeneous cell populations. This study provides a strategy to analyze biophysical events in medicine as well as in neural, stem cell and developmental biology, where small numbers of distinctive or diseased cells play critical roles.     

Sex identification of owls (family Strigidae) using oligonucleotide microarrays

Molecular sexing of the diversified avian family Strigidae is difficult. Sex identification using the intron length difference between W and Z chromosomal CHD1 genes, as visualized by agarose gel electrophoreses, often produces ambiguous results. Here we describe a simple method for sexing a variety of Strigidae species using oligonucleotide microarrays, on which several sex-specific probes operated complementarily or in concert. The sex of 8 owl species was identified clearly on the microarrays through sequence recognition. This sequence-directed method can be easily applied to a wider range of Strigidae species.     

High-resolution identification of chromosomal abnormalities using oligonucleotide arrays containing 116,204 SNPs

Mutation of the human genome ranges from single base-pair changes to whole-chromosome aneuploidy. Karyotyping, fluorescence in situ hybridization, and comparative genome hybridization are currently used to detect chromosome abnormalities of clinical significance. These methods, although powerful, suffer from limitations in speed, ease of use, and resolution, and they do not detect copy-neutral chromosomal aberrations--for example, uniparental disomy (UPD). We have developed a high-throughput approach for assessment of DNA copy-number changes, through use of high-density synthetic oligonucleotide arrays containing 116,204 single-nucleotide polymorphisms, spaced at an average distance of 23.6 kb across the genome. Using this approach, we analyzed samples that failed conventional karyotypic analysis, and we detected amplifications and deletions across a wide range of sizes (1.3-145.9 Mb), identified chromosomes containing anonymous chromatin, and used genotype data to determine the molecular origin of two cases of UPD. Furthermore, our data provided independent confirmation for a case that had been misinterpreted by karyotype analysis. The high resolution of our approach provides more-precise breakpoint mapping, which allows subtle phenotypic heterogeneity to be distinguished at a molecular level. The accurate genotype information provided on these arrays enables the identification of copy-neutral loss-of-heterozygosity events, and the minimal requirement of DNA (250 ng per array) allows rapid analysis of samples without the need for cell culture. This technology overcomes many limitations currently encountered in routine clinical diagnostic laboratories tasked with accurate and rapid diagnosis of chromosomal abnormalities.     

Single nucleotide polymorphism discrimination assisted by improved base stacking hybridization using oligonucleotide microarrays

Efficiencies of mismatch discrimination using size-varied capture probes were examined at various hybridization temperatures. The probes were 17, 15, 13, 11, 9, and 7 nucleotides long and contained single-base mismatches at their 3' ends. The optimal signal intensity and efficiency of base stacking hybridization on mismatch discrimination were observed for capture probes with a melting temperature (Tm) value of 36 degrees C, in the detection of DNA sequence variations at 40 degrees C. We employed asymmetric PCR to prepare single-stranded target DNA labeled with a fluorescent dye, and the PCR product was hybridized on the DNA microarray with no further purification. Our efforts have enhanced the sensitivity and simplified the procedures of base stacking hybridization on mismatch discrimination. As a model experiment, this improved technology was used to identify plasmid templates of human leukocyte antigen (HLA)-A alleles 2601, 2902, and 0206 on oligonucleotide microarrays. It is now possible to apply this simple, rapid, sensitive, and reliable base stacking hybridization technology to detect DNA sequence variations on microarrays in clinical diagnosis and other applications.     

Rapid identification of improved protein ligands using peptoid microarrays

A rapid array-based protocol is presented by which a modest affinity protein-binding small molecule can be appended to a library of peptoids via click chemistry. The array can then be screened for improved ligands that exhibit a higher affinity for the protein target.     

Dynamic diversity and small-molecule evolution: a new paradigm for ligand identification.

A longstanding goal of organic, medicinal and bioorganic chemists has been the discovery of efficient methods for designing or identifying biologically active compounds. Recently, several groups have reported using the directed evolution of combinatorial libraries as a new method of identifying compounds capable of binding tightly to a target molecule. Although significant development remains to be done, the initial results suggest that dynamic diversity and associated selection methods will prove to be valuable additions to the drug-discovery process.     

An improved method for the identification and quantitation of biological lipids by HPLC using laser light-scattering detection

We have developed a simplified and improved high performance liquid chromatography (HPLC) method for the detection and quantitation of tissue lipid using a new laser light-scattering detector (Varex model ELSD II). This detector has a limit of sensitivity of 50 ng for neutral lipid and 200 ng for most phospholipids with excellent reproducibility. By coupling the ELSD II with a ternary gradient normal phase HPLC system, we were able to separate and quantify the major lipid constituents of extracted tissue. This system was used to profile and quantitate the major lipids from rat brain, liver, and heart with greater sensitivity than other available techniques, with the exception of high performance thin-layer chromatography (HPTLC). However, the convenience of HPLC allows for a significant improvement in analysis time with only a threefold reduction in sensitivity when compared to HPTLC.     

Identifying high-affinity aptamer ligands with defined cross-reactivity using high-throughput guided systematic evolution of ligands by exponential enrichment

Oligonucleotide aptamers represent a novel platform for creating ligands with desired specificity, and they offer many potentially significant advantages over monoclonal antibodies in terms of feasibility, cost, and clinical applicability. However, the isolation of high-affinity aptamer ligands from random oligonucleotide pools has been challenging. Although high-throughput sequencing (HTS) promises to significantly facilitate systematic evolution of ligands by exponential enrichment (SELEX) analysis, the enormous datasets generated in the process pose new challenges for identifying those rare, high-affinity aptamers present in a given pool. We show that emulsion PCR preserves library diversity, preventing the loss of rare high-affinity aptamers that are difficult to amplify. We also demonstrate the importance of using reference targets to eliminate binding candidates with reduced specificity. Using a combination of bioinformatics and functional analyses, we show that the rate of amplification is more predictive than prevalence with respect to binding affinity and that the mutational landscape within a cluster of related aptamers can guide the identification of high-affinity aptamer ligands. Finally, we demonstrate the power of this selection process for identifying cross-species aptamers that can bind human receptors and cross-react with their murine orthologs.     

Rapid identification of Aerococcus viridans using the polymerase chain reaction and an oligonucleotide probe

A polymerase chain reaction/oligonucleotide probe method was developed for the specific identification of the Gram-positive bacterium Aerococcus viridans. Primers for the enzymatic amplification reaction were designed from specific sequences within the 16S rRNA. The method was also highly sensitive and 10 cfu of A. viridans could be detected in 5 h although the reliability of detection was poor in mixed cultures with Escherichia coli.     

A flexible method for the fabrication of gold nanostructures using oligonucleotide derivatives

Several linear and branched DNA structures from 80-200 nm with a biotine molecule in the middle have been prepared. These structures have been decorated by addition of positively charged gold nanoparticles carrying 4-(dimethylamino)pyridine ligands. Streptavidin binds to the central biotine molecule introducing a 20 nm gap in the structure in which a biotinylated nanoparticle can be introduced. The simplest structure (80 nm, linear) is formed by 4 oligonucleotides. By changing some of these components changes on length, shape, and recognition system easily can be introduced.     

A reliable method for northern blot analysis using synthetic oligonucleotide probes.

We have developed a method for using short (30-42 base pair) synthetic oligonucleotide DNA probes in Northern blot assays. The method involves labeling the probes to high specific activity, very stringent hybridization and wash conditions, and the presence of several inhibitors of nonspecific binding in the hybridization buffer. We have tested this method with several probes obtained from local and commercial sources. The results with every probe used were high signal-to-noise ratios in an exposure time range of 30 min to 7 days.     

Prediction of phosphotyrosine signaling networks using a scoring matrix-assisted ligand identification approach

Systematic identification of binding partners for modular domains such as Src homology 2 (SH2) is important for understanding the biological function of the corresponding SH2 proteins. We have developed a worldwide web-accessible computer program dubbed SMALI for scoring matrix-assisted ligand identification for SH2 domains and other signaling modules. The current version of SMALI harbors 76 unique scoring matrices for SH2 domains derived from screening oriented peptide array libraries. These scoring matrices are used to search a protein database for short peptides preferred by an SH2 domain. An experimentally determined cut-off value is used to normalize an SMALI score, therefore allowing for direct comparison in peptide-binding potential for different SH2 domains. SMALI employs distinct scoring matrices from Scansite, a popular motif-scanning program. Moreover, SMALI contains built-in filters for phosphoproteins, Gene Ontology (GO) correlation and colocalization of subject and query proteins. Compared to Scansite, SMALI exhibited improved accuracy in identifying binding peptides for SH2 domains. Applying SMALI to a group of SH2 domains identified hundreds of interactions that overlap significantly with known networks mediated by the corresponding SH2 proteins, suggesting SMALI is a useful tool for facile identification of signaling networks mediated by modular domains that recognize short linear peptide motifs.     

Detection of known mutations in hypertrophic cardiomyopathy using oligonucleotide microarrays assisted by improved base stacking hybridization

With the assistance of improved base stacking hybridization, a low-density microarray, containing 12 capture probes, was used to identify 7 known hypertrophic cardiomyopathy-related mutations in the gene of MYH7 (-myosin heavy chain). The hybridization targets, amplified from 11 plasmids containing wild type or mutation sequences of MYH7 and healthy genomic DNA, were prepared by single-step fluorescence labeled asymmetric PCR. Six single base substitutions and a trinucleotide deletion were identified unambiguously, and the average discrimination ratio (Q_mut) for artificial heterozygous samples was as high as 16.2.