Use of extremely short F?rster resonance energy transfer probes in real-time polymerase chain reaction

Described in the article is a new approach for the sequence-specific detection of nucleic acids in real-time polymerase chain reaction (PCR) using fluorescently labeled oligonucleotide probes. The method is based on the production of PCR amplicons, which fold into dumbbell-like secondary structures carrying a specially designed ‘probe-luring’ sequence at their 5′ ends. Hybridization of this sequence to a complementary ‘anchoring’ tail introduced at the 3′ end of a fluorescent probe enables the probe to bind to its target during PCR, and the subsequent probe cleavage results in the florescence signal. As it has been shown in the study, this amplicon-endorsed and guided formation of the probe-target duplex allows the use of extremely short oligonucleotide probes, up to tetranucleotides in length. In particular, the short length of the fluorescent probes makes possible the development of a ‘universal’ probe inventory that is relatively small in size but represents all possible sequence variations. The unparalleled cost-effectiveness of the inventory approach is discussed. Despite the short length of the probes, this new method, named Angler real-time PCR, remains highly sequence specific, and the results of the study indicate that it can be effectively used for quantitative PCR and the detection of polymorphic variations.     

Optimization of the BLASTN substitution matrix for prediction of non-specific DNA microarray hybridization

DNA microarray measurements are susceptible to error caused by non-specific hybridization between a probe and a target (cross-hybridization), or between two targets (bulk-hybridization). Search algorithms such as BLASTN can quickly identify potentially hybridizing sequences. We set out to improve BLASTN accuracy by modifying the substitution matrix and gap penalties. We generated gene expression microarray data for samples in which 1 or 10% of the target mass was an exogenous spike of known sequence. We found that the 10% spike induced 2-fold intensity changes in 3% of the probes, two-third of which were decreases in intensity likely caused by bulk-hybridization. These changes were correlated with similarity between the spike and probe sequences. Interestingly, even very weak similarities tended to induce a change in probe intensity with the 10% spike. Using this data, we optimized the BLASTN substitution matrix to more accurately identify probes susceptible to non-specific hybridization with the spike. Relative to the default substitution matrix, the optimized matrix features a decreased score for A–T base pairs relative to G–C base pairs, resulting in a 5–15% increase in area under the ROC curve for identifying affected probes. This optimized matrix may be useful in the design of microarray probes, and in other BLASTN-based searches for hybridization partners.     

Targeted resequencing of HIV variants by microarray thermodynamics

Within a single infected individual, a virus population can have a high genomic variability. In the case of HIV, several mutations can be present even in a small genomic window of 20–30 nucleotides. For diagnostics purposes, it is often needed to resequence genomic subsets where crucial mutations are known to occur. In this article, we address this issue using DNA microarrays and inputs from hybridization thermodynamics. Hybridization signals from multiple probes are analysed, including strong signals from perfectly matching (PM) probes and a large amount of weaker cross-hybridization signals from mismatching (MM) probes. The latter are crucial in the data analysis. Seven coded clinical samples (HIV-1) are analyzed, and the microarray results are in full concordance with Sanger sequencing data. Moreover, the thermodynamic analysis of microarray signals resolves inherent ambiguities in Sanger data of mixed samples and provides additional clinically relevant information. These results show the reliability and added value of DNA microarrays for point-of-care diagnostic purposes.     

Novel DNA probes with low background and high hybridization-triggered fluorescence

Novel fluorogenic DNA probes are described. The probes (called Pleiades) have a minor groove binder (MGB) and a fluorophore at the 5′-end and a non-fluorescent quencher at the 3′-end of the DNA sequence. This configuration provides surprisingly low background and high hybridization-triggered fluorescence. Here, we comparatively study the performance of such probes, MGB-Eclipse probes, and molecular beacons. Unlike the other two probe formats, the Pleiades probes have low, temperature-independent background fluorescence and excellent signal-to-background ratios. The probes possess good mismatch discrimination ability and high rates of hybridization. Based on the analysis of fluorescence and absorption spectra we propose a mechanism of action for the Pleiades probes. First, hydrophobic interactions between the quencher and the MGB bring the ends of the probe and, therefore, the fluorophore and the quencher in close proximity. Second, the MGB interacts with the fluorophore and independent of the quencher is able to provide a modest (2–4-fold) quenching effect. Joint action of the MGB and the quencher is the basis for the unique quenching mechanism. The fluorescence is efficiently restored upon binding of the probe to target sequence due to a disruption in the MGB–quencher interaction and concealment of the MGB moiety inside the minor groove.     

Characterization of mismatch and high-signal intensity probes associated with Affymetrix genechips

MOTIVATION: For Affymetrix microarray platforms, gene expression is determined by computing the difference in signal intensities between perfect match (PM) and mismatch (MM) probesets. Although the use of PM is not controversial, MM probesets have been associated with variance and ultimately inaccurate gene expression calls. A principal focus of this study was to investigate the nature of the MM signal intensities and demonstrate its contribution to the experimental results. RESULTS: While most MM intensities were likely associated with random noise, a subset of approximately 20% (99,485) of the MM probes displayed relatively high signal intensities to the corresponding PM probes (MM > PM) in a non-random fashion; 13,440 of these probes demonstrated exceptionally high 'outlier' intensities. About 15,938 PM probes also demonstrated exceptionally high outlier intensities consistently across all hybridizations. About 92% of the MM > PM probes had either a dThymidine (dT) or a dCytidine (dC) at the 13th position of the probe sequence. MM and PM probes displaying extremely high outlier intensities contained high dC rich nucleotides, and low dA contents at other nucleotides positions along the 25mer probe sequence. Differentially expressed genes generated using Genechip Operating System (GCOS) or modified PM-only methods were also examined. Of those candidate genes identified in the PM-only method, 157 of them were designated by GCOS as absent across all datasets and many others contained probes with MM > PM signal intensities. Our data suggests that MM intensity from PM signal can be a major source of error analysis, leading to fewer potentially biologically important candidate genes. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

AffyMAPSDetector: a software tool to characterize Affymetrix GeneChip™ expression arrays with respect to SNPs

Background  

Affymetrix gene expression arrays incorporate paired perfect match (PM) and mismatch (MM) probes to distinguish true signals from those arising from cross-hybridization events. A MM signal often shows greater intensity than a PM signal; we propose that one underlying cause is the presence of allelic variants arising from single nucleotide polymorphisms (SNPs). To annotate and characterize SNP contributions to anomalous probe binding behavior we have developed a software tool called AffyMAPSDetector.     

Bacillus subtilis RNase J1 endonuclease and 5′ exonuclease activities in the turnover of ΔermC mRNA

RNase J1, a ribonuclease with 5′ exonuclease and endonuclease activities, is an important factor in Bacillus subtilis mRNA decay. A model for RNase J1 endonuclease activity in mRNA turnover has RNase J1 binding to the 5′ end and tracking to a target site downstream, where it makes a decay-initiating cleavage. The upstream fragment from this cleavage is degraded by 3′ exonucleases; the downstream fragment is degraded by RNase J1 5′ exonuclease activity. Previously, ΔermC mRNA was used to show 5′-end dependence of mRNA turnover. Here we used ΔermC mRNA to probe RNase J1-dependent degradation, and the results were consistent with aspects of the model. ΔermC mRNA showed increased stability in a mutant strain that contained a reduced level of RNase J1. In agreement with the tracking concept, insertion of a strong stem–loop structure at +65 resulted in increased stability. Weakening this stem–loop structure resulted in reversion to wild-type stability. RNA fragments containing the 3′ end were detected in a strain with reduced RNase J1 expression, but were undetectable in the wild type. The 5′ ends of these fragments mapped to the upstream side of predicted stem–loop structures, consistent with an impediment to RNase J1 5′ exonuclease processivity. A ΔermC mRNA deletion analysis suggested that decay-initiating endonuclease cleavage could occur at several sites near the 3′ end. However, even in the absence of these sites, stability was further increased in a strain with reduced RNase J1, suggesting alternate pathways for decay that could include exonucleolytic decay from the 5′ end.     

Redefinition of Affymetrix probe sets by sequence overlap with cDNA microarray probes reduces cross-platform inconsistencies in cancer-associated gene expression measurements

Background  

Comparison of data produced on different microarray platforms often shows surprising discordance. It is not clear whether this discrepancy is caused by noisy data or by improper probe matching between platforms. We investigated whether the significant level of inconsistency between results produced by alternative gene expression microarray platforms could be reduced by stringent sequence matching of microarray probes. We mapped the short oligo probes of the Affymetrix platform onto cDNA clones of the Stanford microarray platform. Affymetrix probes were reassigned to redefined probe sets if they mapped to the same cDNA clone sequence, regardless of the original manufacturer-defined grouping. The NCI-60 gene expression profiles produced by Affymetrix HuFL platform were recalculated using these redefined probe sets and compared to previously published cDNA measurements of the same panel of RNA samples.     

A new label technology for the detection of specific polymerase chain reaction products in a closed tube

A novel signal generation principle suitable for real time and end-point detection of specific PCR products in a closed tube is described. Linear DNA probes were labeled at their 5′-ends with a stable, fluorescent terbium chelate. The fluorescence intensity of this chelate is lower when it is coupled to single-stranded DNA than when the chelate is free in solution. The synthesized probes were used in the real time monitoring of PCR using a prototype instrument that consisted of a fluorometer coupled to a thermal cycler. When the probe anneals to a complementary target amplicon, the 5′→3′ exonucleolytic activity of DNA polymerase detaches the label from the probe. This results in an enhanced terbium fluorescence signal. Since terbium has a long excited state lifetime, its fluorescence can be measured in a time-resolved manner, which results in a low background fluorescence and a 1000-fold signal amplification. The detection method is quantitative over an extremely wide linear range (at least 10–10⁷ initial template molecules). The label strategy can easily be combined with existing label technologies, such as TaqMan 5′-exonuclease assays, in order to carry out multiplex assays that do not suffer from overlapping emission peaks of the fluorophores.     

Roles for TbDSS-1 in RNA surveillance and decay of maturation by-products from the 12S rRNA locus

The Trypanosoma brucei exoribonuclease, TbDSS-1, has been implicated in multiple aspects of mitochondrial RNA metabolism. Here, we investigate the role of TbDSS-1 in RNA processing and surveillance by analyzing 12S rRNA processing intermediates in TbDSS-1 RNAi cells. RNA fragments corresponding to leader sequence upstream of 12S rRNA accumulate upon TbDSS-1 depletion. The 5′ extremity of 12S rRNA is generated by endonucleolytic cleavage, and TbDSS-1 degrades resulting upstream maturation by-products. RNAs with 5′ ends at position −141 and 3′ ends adjacent to the mature 5′ end of 12S rRNA are common and invariably possess oligo(U) tails. 12S rRNAs with mature 3′ ends and unprocessed 5′ ends also accumulate in TbDSS-1 depleted cells, suggesting that these RNAs represent dead-end products normally destined for decay by TbDSS-1 in an RNA surveillance pathway. Together, these data indicate dual roles for TbDSS-1 in degradation of 12S rRNA maturation by-products and as part of a mitochondrial RNA surveillance pathway that eliminates stalled 12S processing intermediates. We further provide evidence that TbDSS-1 degrades RNAs originating upstream of the first gene on the minor strand of the mitochondrial maxicircle suggesting that TbDSS-1 also removes non-functional RNAs generated from other regions of the mitochondrial genome.     

New approach to real-time nucleic acids detection: folding polymerase chain reaction amplicons into a secondary structure to improve cleavage of F?rster resonance energy transfer probes in 5′-nuclease assays

The article describes a new technology for real-time polymerase chain reaction (PCR) detection of nucleic acids. Similar to Taqman, this new method, named Snake, utilizes the 5′-nuclease activity of Thermus aquaticus (Taq) DNA polymerase that cleaves dual-labeled Förster resonance energy transfer (FRET) probes and generates a fluorescent signal during PCR. However, the mechanism of the probe cleavage in Snake is different. In this assay, PCR amplicons fold into stem–loop secondary structures. Hybridization of FRET probes to one of these structures leads to the formation of optimal substrates for the 5′-nuclease activity of Taq. The stem–loop structures in the Snake amplicons are introduced by the unique design of one of the PCR primers, which carries a special 5′-flap sequence. It was found that at a certain length of these 5′-flap sequences the folded Snake amplicons have very little, if any, effect on PCR yield but benefit many aspects of the detection process, particularly the signal productivity. Unlike Taqman, the Snake system favors the use of short FRET probes with improved fluorescence background. The head-to-head comparison study of Snake and Taqman revealed that these two technologies have more differences than similarities with respect to their responses to changes in PCR protocol, e.g. the variations in primer concentration, annealing time, PCR asymmetry. The optimal PCR protocol for Snake has been identified. The technology’s real-time performance was compared to a number of conventional assays including Taqman, 3′-MGB-Taqman, Molecular Beacon and Scorpion primers. The test trial showed that Snake supersedes the conventional assays in the signal productivity and detection of sequence variations as small as single nucleotide polymorphisms. Due to the assay’s cost-effectiveness and simplicity of design, the technology is anticipated to quickly replace all known conventional methods currently used for real-time nucleic acid detection.     

Hybridization of single-stranded DNA targets to immobilized complementary DNA probes: comparison of hairpin versus linear capture probes

A microtiter-based assay system is described in which DNA hairpin probes with dangling ends and single-stranded, linear DNA probes were immobilized and compared based on their ability to capture single-strand target DNA. Hairpin probes consisted of a 16 bp duplex stem, linked by a T₂-biotin·dT-T₂ loop. The third base was a biotinylated uracil (U_B) necessary for coupling to avidin coated microtiter wells. The capture region of the hairpin was a 3′ dangling end composed of either 16 or 32 bases. Fundamental parameters of the system, such as probe density and avidin adsorption capacity of the plates were characterized. The target DNA consisted of 65 bases whose 3′ end was complementary to the dangling end of the hairpin or to the linear probe sequence. The assay system was employed to measure the time dependence and thermodynamic stability of target hybridization with hairpin and linear probes. Target molecules were labeled with either a 5′-FITC, or radiolabeled with [γ-³³P]ATP and captured by either linear or hairpin probes affixed to the solid support. Over the range of target concentrations from 10 to 640 pmol hybridization rates increased with increasing target concentration, but varied for the different probes examined. Hairpin probes displayed higher rates of hybridization and larger equilibrium amounts of captured targets than linear probes. At 25 and 45°C, rates of hybridization were better than twice as great for the hairpin compared with the linear capture probes. Hairpin–target complexes were also more thermodynamically stable. Binding free energies were evaluated from the observed equilibrium constants for complex formation. Results showed the order of stability of the probes to be: hairpins with 32 base dangling ends > hairpin probes with l6 base dangling ends > 16 base linear probes > 32 base linear probes. The physical characteristics of hairpins could offer substantial advantages as nucleic acid capture moieties in solid support based hybridization systems.     

Isoenergetic penta- and hexanucleotide microarray probing and chemical mapping provide a secondary structure model for an RNA element orchestrating R2 retrotransposon protein function

LNA (locked nucleic acids, i.e. oligonucleotides with a methyl bridge between the 2′ oxygen and 4′ carbon of ribose) and 2,6-diaminopurine were incorporated into 2′-O-methyl RNA pentamer and hexamer probes to make a microarray that binds unpaired RNA approximately isoenergetically. That is, binding is roughly independent of target sequence if target is unfolded. The isoenergetic binding and short probe length simplify interpretation of binding to a structured RNA to provide insight into target RNA secondary structure. Microarray binding and chemical mapping were used to probe the secondary structure of a 323 nt segment of the 5′ coding region of the R2 retrotransposon from Bombyx mori (R2Bm 5′ RNA). This R2Bm 5′ RNA orchestrates functioning of the R2 protein responsible for cleaving the second strand of DNA during insertion of the R2 sequence into the genome. The experimental results were used as constraints in a free energy minimization algorithm to provide an initial model for the secondary structure of the R2Bm 5′ RNA.     

Light-directed 5′→3′ synthesis of complex oligonucleotide microarrays

Light-directed synthesis of high-density microarrays is currently performed in the 3′→5′ direction due to constraints in existing synthesis chemistry. This results in the probes being unavailable for many common types of enzymatic modification. Arrays that are synthesized in the 5′→3′ direction could be utilized to perform parallel genotyping and resequencing directly on the array surface, dramatically increasing the throughput and reducing the cost relative to existing techniques. In this report we demonstrate the use of photoprotected phosphoramidite monomers for light-directed array synthesis in the 5′→3′ direction, using maskless array synthesis technology. These arrays have a dynamic range of >2.5 orders of magnitude, sensitivity below 1 pM and a coefficient of variance of <10% across the array surface. Arrays containing >150 000 probe sequences were hybridized to labeled mouse cRNA producing highly concordant data (average R² = 0.998). We have also shown that the 3′ ends of array probes are available for sequence-specific primer extension and ligation reactions.     

Pat1 contributes to the RNA binding activity of the Lsm1-7–Pat1 complex

A major mRNA decay pathway in eukaryotes is initiated by deadenylation followed by decapping of the oligoadenylated mRNAs and subsequent 5′-to-3′ exonucleolytic degradation of the capless mRNA. In this pathway, decapping is a rate-limiting step that requires the hetero-octameric Lsm1-7–Pat1 complex to occur at normal rates in vivo. This complex is made up of the seven Sm-like proteins, Lsm1 through Lsm7, and the Pat1 protein. It binds RNA and has a unique binding preference for oligoadenylated RNAs over polyadenylated RNAs. Such binding ability is crucial for its mRNA decay function in vivo. In order to determine the contribution of Pat1 to the function of the Lsm1-7–Pat1 complex, we compared the RNA binding properties of the Lsm1-7 complex purified from pat1Δ cells and purified Pat1 fragments with that of the wild-type Lsm1-7–Pat1 complex. Our studies revealed that both the Lsm1-7 complex and purified Pat1 fragments have very low RNA binding activity and are impaired in the ability to recognize the oligo(A) tail on the RNA. However, reconstitution of the Lsm1-7–Pat1 complex from these components restored these abilities. We also observed that Pat1 directly contacts RNA in the context of the Lsm1-7–Pat1 complex. These studies suggest that the unique RNA binding properties and the mRNA decay function of the Lsm1-7–Pat1 complex involve cooperation of residues from both Pat1 and the Lsm1-7 ring. Finally our studies also revealed that the middle domain of Pat1 is essential for the interaction of Pat1 with the Lsm1-7 complex in vivo.     

A robust method for estimating gene expression states using Affymetrix microarray probe level data

Background  

Microarray technology is a high-throughput method for measuring the expression levels of thousand of genes simultaneously. The observed intensities combine a non-specific binding, which is a major disadvantage with microarray data. The Affymetrix GeneChip assigned a mismatch (MM) probe with the intention of measuring non-specific binding, but various opinions exist regarding usefulness of MM measures. It should be noted that not all observed intensities are associated with expressed genes and many of those are associated with unexpressed genes, of which measured values express mere noise due to non-specific binding, cross-hybridization, or stray signals. The implicit assumption that all genes are expressed leads to poor performance of microarray data analyses. We assume two functional states of a gene - expressed or unexpressed - and propose a robust method to estimate gene expression states using an order relationship between PM and MM measures.     

Tests of rRNA hybridization to microarrays suggest that hybridization characteristics of oligonucleotide probes for species discrimination cannot be predicted

Hybridization of rRNAs to microarrays is a promising approach for prokaryotic and eukaryotic species identification. Typically, the amount of bound target is measured by fluorescent intensity and it is assumed that the signal intensity is directly related to the target concentration. Using thirteen different eukaryotic LSU rRNA target sequences and 7693 short perfect match oligonucleotide probes, we have assessed current approaches for predicting signal intensities by comparing Gibbs free energy (ΔG°) calculations to experimental results. Our evaluation revealed a poor statistical relationship between predicted and actual intensities. Although signal intensities for a given target varied up to 70-fold, none of the predictors were able to fully explain this variation. Also, no combination of different free energy terms, as assessed by principal component and neural network analyses, provided a reliable predictor of hybridization efficiency. We also examined the effects of single-base pair mismatch (MM) (all possible types and positions) on signal intensities of duplexes. We found that the MM effects differ from those that were predicted from solution-based hybridizations. These results recommend against the application of probe design software tools that use thermodynamic parameters to assess probe quality for species identification. Our results imply that the thermodynamic properties of oligonucleotide hybridization are by far not yet understood.