Facile SNP detection using bifunctional, cross-linking oligonucleotide probes

A facile, sensitive method for detecting specific sequences of oligonucleotides was developed. Detection of DNA sequences with single nucleotide discrimination is achieved by combining the selectivity of hybridization with an efficient cross-linking reaction. Readily synthesized bifunctional oligonucleotide probes containing a modified pyrimidine that is capable of forming interstrand cross-links under mild oxidative conditions internally, and biotin at their 5′-termini were used to discriminate between 16-nt long sites in plasmid DNA that differ by a single nucleotide. The target sequence was detected via fluorescence spectroscopy by utilizing conjugates of avidin and horseradish peroxidase in a microtiter plate assay. The method is able to detect as little as 250 fmol of target without using PCR and exhibits single nucleotide discrimination that approaches 200:1. In principle, this method is capable of probing any target sequence containing a 2′-deoxyadenosine.     

Thermodynamic stability of base pairs between 2-hydroxyadenine and incoming nucleotides as a determinant of nucleotide incorporation specificity during replication

We investigated the thermodynamic stability of double-stranded DNAs with an oxidative DNA lesion, 2-hydroxyadenine (2-OH-Ade), in two different sequence contexts (5′-GA*C-3′ and 5′-TA*A-3′, A* represents 2-OH-Ade). When an A*–N pair (N, any nucleotide base) was located in the center of a duplex, the thermodynamic stabilities of the duplexes were similar for all the natural bases except A (N = T, C and G). On the other hand, for the duplexes with the A*–N pair at the end, which mimic the nucleotide incorporation step, the stabilities of the duplexes were dependent on their sequence. The order of stability is T > G > C >> A in the 5′-GA*C-3′ sequences and T > A > C > G in the 5′-TA*A-3′ sequences. Because T/G/C and T/A are nucleotides incorporated opposite to 2-OH-Ade in the 5′-GA*C-3′ and 5′-TA*A-3′ sequences, respectively, these results agree with the tendency of mutagenic misincorporation of the nucleotides opposite to 2-OH-Ade in vitro. Thus, the thermodynamic stability of the A*–N base pair may be an important factor for the mutation spectra of 2-OH-Ade.     

Real-time closed tube single nucleotide polymorphism (SNP) quantification in pooled samples by quencher extension (QEXT)

Quencher extension (QEXT) is a novel single step closed tube real-time method to quantify SNPs using reporters and quenchers in combination with primer extension. A probe with a 5′-reporter dye is single base extended with a dideoxy nucleotide containing a quencher dye if the target SNP allele is present. The extension is recorded from the quenching (reduced fluorescence) of the reporter dye. This avoids the influence of the unincorporated dye-labeled nucleotides, resulting in high accuracy and a high signal-to-noise ratio. The relative amount of a specific SNP allele is determined from the nucleotide incorporation rate in a thermo-cycling reaction. We tested the QEXT assay using five SNPs in the Listeria monocytogenes inlA gene as a model system. The presence of the target SNP alleles was determined with high statistical confidence (P < 0.0005). The quantitative detection limits were between 0 and 5% for the targeted SNP alleles on a background of other SNP alleles (P < 0.05). The QEXT method is directly adaptable to current real-time PCR equipment and is thus suited for high throughput and a wide application range.     

The Accumulation of Calcium Ions by Sarcotubular Vesicles

The accumulation of Ca⁺⁺ by microsomal (sarcotubular) preparations of rabbit skeletal muscle in the presence of oxalate, and the concurrent splitting of nucleoside triphosphate, displayed moderate nucleotide specificity in the sequence ATP > GTP, CTP, ITP > UTP > (ADP) > ATetraP for the former, ATP > (ADP) > ITP > GTP > CTP > UTP > ATetraP for the latter process. The "calcium pump" was weakly inhibited by caffeine, and was inhibited together with the ATPase by pyridoxalphosphate. Carnosine had no effect as such nor in the presence of pyridoxalphosphate except at high concentration; thiourea and p-chloromercuribenzoate were inhibiting while iodoacetate was inactive. Ca⁺⁺ accumulation and ATPase were inhibited by atabrine (not tested on ATPase), dinitrophenol, and amytal. High concentrations of oligomycin and rutamycin inhibited Ca⁺⁺ uptake while slightly stimulating ATPase. Antimycin A stimulated the Ca⁺⁺ uptake. These results are discussed in the light of their possible relation to partial reactions in oxidative phosphorylation. The Ca⁺⁺ uptake and relaxing factor activities did not behave identically throughout. This is in part ascribed to changes in reactivity of actomyosin in the relaxation test, in part to the participation of relaxing substances other than the calcium pump.     

APOBEC3G cytosine deamination hotspots are defined by both sequence context and single-stranded DNA secondary structure

Apolipoprotein B mRNA-editing, enzyme-catalytic, polypeptide-like 3G (i.e., APOBEC3G or A3G) is an evolutionarily conserved cytosine deaminase that potently restricts human immunodeficiency virus type 1 (HIV-1), retrotransposons and other viruses. A3G has a nucleotide target site specificity for cytosine dinucleotides, though only certain cytosine dinucleotides are ‘hotspots’ for cytosine deamination, and others experience little or no editing by A3G. The factors that define these critical A3G hotspots are not fully understood. To investigate how A3G hotspots are defined, we used an in vitro fluorescence resonance energy transfer-based oligonucleotide assay to probe the site specificity of A3G. Our findings strongly suggest that the target single-stranded DNA (ssDNA) secondary structure as well as the bases directly 3′ and 5′ of the cytosine dinucleotide are critically important A3G recognition. For instance, A3G cannot readily deaminate a cytosine dinucleotide in ssDNA stem structures or in nucleotide base loops composed of three bases. Single-stranded nucleotide loops up to seven bases in length were poor targets for A3G activity unless cytosine residues flanked the cytosine dinucleotide. Furthermore, we observed that A3G favors adenines, cytosines and thymines flanking the cytosine dinucleotide target in unstructured regions of ssDNA. Low cytosine deaminase activity was detected when guanines flanked the cytosine dinucleotide. Taken together, our findings provide the first demonstration that A3G cytosine deamination hotspots are defined by both the sequence context of the cytosine dinucleotide target as well as the ssDNA secondary structure. This knowledge can be used to better trace the origins of mutations to A3G activity, and illuminate its impact on processes such as HIV-1 genetic variation.     

Hybridization of 2'-O-methyl and 2'-deoxy molecular beacons to RNA and DNA targets

Molecular beacons are stem–loop hairpin oligonucleotide probes labeled with a fluorescent dye at one end and a fluorescence quencher at the other end; they can differentiate between bound and unbound probes in homogeneous hybridization assays with a high signal-to-background ratio and enhanced specificity compared with linear oligonucleotide probes. However, in performing cellular imaging and quantification of gene expression, degradation of unmodified molecular beacons by endogenous nucleases can significantly limit the detection sensitivity, and results in fluorescence signals unrelated to probe/target hybridization. To substantially reduce nuclease degradation of molecular beacons, it is possible to protect the probe by substituting 2′-O-methyl RNA for DNA. Here we report the analysis of the thermodynamic and kinetic properties of 2′-O-methyl and 2′-deoxy molecular beacons in the presence of RNA and DNA targets. We found that in terms of molecular beacon/target duplex stability, 2′-O-methyl/RNA > 2′-deoxy/RNA > 2′-deoxy/DNA > 2′-O-methyl/DNA. The improved stability of the 2′-O-methyl/RNA duplex was accompanied by a slightly reduced specificity compared with the duplex of 2′-deoxy molecular beacons and RNA targets. However, the 2′-O-methyl molecular beacons hybridized to RNA more quickly than 2′-deoxy molecular beacons. For the pairs tested, the 2′-deoxy-beacon/DNA-target duplex showed the fastest hybridization kinetics. These findings have significant implications for the design and application of molecular beacons.     

Fractal landscapes and molecular evolution: modeling the myosin heavy chain gene family.

Mapping nucleotide sequences onto a "DNA walk" produces a novel representation of DNA that can then be studied quantitatively using techniques derived from fractal landscape analysis. We used this method to analyze 11 complete genomic and cDNA myosin heavy chain (MHC) sequences belonging to 8 different species. Our analysis suggests an increase in fractal complexity for MHC genes with evolution with vertebrate > invertebrate > yeast. The increase in complexity is measured by the presence of long-range power-law correlations, which are quantified by the scaling exponent alpha. We develop a simple iterative model, based on known properties of polymeric sequences, that generates long-range nucleotide correlations from an initially noncorrelated coding region. This new model-as well as the DNA walk analysis-both support the intron-late theory of gene evolution.     

Electron transfer in DNA duplexes containing 2-methyl-1,4-naphthoquinone

2-Methyl-1,4-naphthoquinone (menadione, MQ) was linked to synthetic oligonucleotides and exposed to near-UV light to generate base radical cations in DNA. This model system of electron transfer induced alkali-labile breaks at GG doublets, similar to anthraquinone and metallointercalators systems. In sharp contrast to other systems, the photolysis of MQ–DNA duplexes gave interstrand cross-links and alkali-labile breaks at bases on the complementary strand opposite the MQ moiety. For sequences with an internal MQ, the formation of cross-links with A and C opposite the MQ moiety was 2- to 3-fold greater than that with G and T. The yield of cross-links was more than 10-fold greater than that of breaks opposite MQ, which in turn was more than 2-fold greater than breaks at GG doublets. The yield of damage at GG doublets greatly increased for a sequence with a terminal MQ. The distribution of base damage was measured by enzymatic digestion and HPLC analysis (dAdo > dThd > dGuo > dCyd). The formation of novel products in MQ–DNA duplexes was attributed to the ability of excited MQ to generate the radical cations of all four DNA bases; thus, this photochemical reaction provides an ideal model system to study the effects of ionizing radiation and one-electron oxidants.     

Sequence-specific recognition of colicin E5, a tRNA-targeting ribonuclease

Colicin E5 is a novel Escherichia coli ribonuclease that specifically cleaves the anticodons of tRNA^Tyr, tRNA^His, tRNA^Asn and tRNA^Asp. Since this activity is confined to its 115 amino acid long C-terminal domain (CRD), the recognition mechanism of E5-CRD is of great interest. The four tRNA substrates share the unique sequence UQU within their anticodon loops, and are cleaved between Q (modified base of G) and 3′ U. Synthetic minihelix RNAs corresponding to the substrate tRNAs were completely susceptible to E5-CRD and were cleaved in the same manner as the authentic tRNAs. The specificity determinant for E5-CRD was YGUN at −1 to +3 of the ‘anticodon’. The YGU is absolutely required and the extent of susceptibility of minihelices depends on N (third letter of the anticodon) in the order A > C > G > U accounting for the order of susceptibility tRNA^Tyr > tRNA^Asp > tRNA^His, tRNA^Asn. Contrastingly, we showed that GpUp is the minimal substrate strictly retaining specificity to E5-CRD. The effect of contiguous nucleotides is inconsistent between the loop and linear RNAs, suggesting that nucleotide extension on each side of GpUp introduces a structural constraint, which is reduced by a specific loop structure formation that includes a 5′ pyrimidine and 3′ A.     

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.     

Use of Subtractive Hybridization To Design Habitat-Based Oligonucleotide Probes for Investigation of Natural Bacterial Communities

We describe a rapid oligonucleotide probe design strategy based on subtractive hybridization which yields probes for 16S rRNA or rRNA genes of individual members of microbial communities that are specific within the context of those communities. This strategy circumvents the need to sequence many similar or identical clones of dominant members of a community. Radioactively labeled subfragments of a cloned 16S rRNA gene sequence for which a probe is required (target) were hybridized with biotinylated total 16S ribosomal DNA (rDNA) amplified from the microbial community, and the hybrids formed were subsequently discarded. The remaining enriched fragments were used to screen a library consisting of cloned subfragments of the target sequence by colony hybridization in order to identify the variable regions of the 16S rRNA gene with the required specificity. The sequencing of random clones in one 16S rDNA library demonstrated that only those clones with 100% sequence identity with the probe fragment were detected by it. Moreover, sequencing of other, randomly selected, probe-positive clones revealed 100% sequence identity with the probe. Probes developed in this way tended to correspond to more variable regions of the 16S rRNA if the target sequences were similar to the sequences of other clones in the library and to less variable regions if the target sequences were phylogenetically isolated within the clone library. Although the absolute specificity of the latter probes, as assessed by comparison with available database sequences, was lower than the absolute specificity of the probes from the more variable regions, they were specific within the context of the environmental samples from which they were derived.     

Targeted gene modification in mismatch-repair-deficient embryonic stem cells by single-stranded DNA oligonucleotides

Gene targeting through homologous recombination in murine embryonic stem (ES) cells is already strongly suppressed by DNA mismatch-repair (MMR)-dependent anti-recombination when targeting construct and target locus differ at <1% of the nucleotide positions. We demonstrate that MMR activity also raises a strong impediment to gene modification mediated by small synthetic DNA oligonucleotide sequences. In the absence of the DNA MMR gene MSH2, synthetic single-stranded deoxyribo-oligonucleotides can be used to site-specifically modify the ES cell genome. We show that PCR-based procedures can be used to identify and clone modified cells. By this method we have substituted a single codon in the retinoblastoma gene.     

Adenosine triphosphatase from soybean callus and root cells

The ATPase activity of a membrane fraction from soybean (Glycine max L.) root and callus cells, presumed to be enriched in plasma membrane, has been characterized with respect to ion stimulation, pH requirement, and nucleotide specificity. The enzyme from both sources was activated by divalent cations (Mg²⁺ > Mn²⁺ > Zn²⁺ > Ca²⁺ > Sr²⁺) and further stimulated by monovalent salts. Preparations from root cells were stimulated by monovalent ions according to the sequence: K⁺ > Rb⁺ > Choline⁺ > Na⁺ > Li⁺ > NH₄⁺ > Cs⁺ > tris⁺. Membrane preparations from callus cells showed similar stimulatory patterns except for a slight preference for Na⁺ over K⁺. No synergism between K⁺ and Na⁺ was found with preparations from either cell source.     

Influence of the three nucleotides upstream of the initiation codon on expression of the Escherichia coli lacZ gene in Saccharomyces cerevisiae.

By introducing synthetic oligonucleotides into a lacZ-yeast expression vector a set of 47 plasmids (out of 64 possible) was generated, differing only in the three bases immediately upstream of the AUG initiation codon of the Escherichia coli lacZ gene. Expression of the beta-galactosidase fusion protein encoded by the different plasmids was determined in Saccharomyces cerevisiae by immunogel electrophoresis. Among the clones tested we found a factor 3 difference in expression. A slight nucleotide preference was found in positions -3(A > G > C = U) and -2 (G > C = U > A). The choice of the nucleotide at position -1 immediately 5' of the AUG did not effect translation efficiency. Increasing homology to the yeast consensus sequence (AAAAAAAUGUCU) was not concomitant with an increased translation efficiency. Our results indicate that the choice of nucleotides immediately preceding the initiation codon in yeast does not dramatically influence translation efficiency, as in prokaryotes or higher eukaryotes.     

Novel PCR Primers for the Archaeal Phylum Thaumarchaeota Designed Based on the Comparative Analysis of 16S rRNA Gene Sequences

Based on comparative phylogenetic analysis of 16S rRNA gene sequences deposited in an RDP database, we constructed a local database of thaumarchaeotal 16S rRNA gene sequences and developed a novel PCR primer specific for the archaeal phylum Thaumarchaeota. Among 9,727 quality-filtered (chimeral-checked, size >1.2 kb) archaeal sequences downloaded from the RDP database, 1,549 thaumarchaeotal sequences were identified and included in our local database. In our study, Thaumarchaeota included archaeal groups MG-I, SAGMCG-I, SCG, FSCG, RC, and HWCG-III, forming a monophyletic group in the phylogenetic tree. Cluster analysis revealed 114 phylotypes for Thaumarchaeota. The majority of the phylotypes (66.7%) belonged to the MG-I and SCG, which together contained most (93.9%) of the thaumarchaeotal sequences in our local database. A phylum-directed primer was designed from a consensus sequence of the phylotype sequences, and the primer’s specificity was evaluated for coverage and tolerance both in silico and empirically. The phylum-directed primer, designated THAUM-494, showed >90% coverage for Thaumarchaeota and <1% tolerance to non-target taxa, indicating high specificity. To validate this result experimentally, PCRs were performed with THAUM-494 in combination with a universal archaeal primer (ARC917R or 1017FAR) and DNAs from five environmental samples to construct clone libraries. THAUM-494 showed a satisfactory specificity in empirical studies, as expected from the in silico results. Phylogenetic analysis of 859 cloned sequences obtained from 10 clone libraries revealed that >95% of the amplified sequences belonged to Thaumarchaeota. The most frequently sampled thaumarchaeotal subgroups in our samples were SCG, MG-I, and SAGMCG-I. To our knowledge, THAUM-494 is the first phylum-level primer for Thaumarchaeota. Furthermore, the high coverage and low tolerance of THAUM-494 will make it a potentially valuable tool in understanding the phylogenetic diversity and ecological niche of Thaumarchaeota.     

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.     

Characterization of fully 2'-modified oligoribonucleotide hetero- and homoduplex hybridization and nuclease sensitivity.

The nuclease stability and melting temperatures (Tm) were compared for fully modified oligoribonucleotide sequences containing 2'-fluoro, 2'-O-methyl, 2'-O-propyl and 2'-O-pentyl nucleotides. Duplexes formed between 2' modified oligoribonucleotides and RNA have typical A-form geometry as observed by circular dichroism spectroscopy. Modifications, with the exception of 2'-O-pentyl, were observed to increase the Tm of duplexes formed with complementary RNA. Modified homoduplexes showed significantly higher Tms, with the following Tm order: 2'-fluoro:2'fluoro > 2'-O-propyl:2'-O-propyl > 2'-O-methyl:2'-O- methyl > RNA:RNA > DNA:DNA. The nuclease stability of 2'-modified oligoribonucleotides was examined using snake venom phosphodiesterase (SVPD) and nuclease S1. The stability imparted by 2' modifications was observed to correlate with the size of the modification. An additional level of nuclease stability was present in oligoribonucleotides having the potential for forming secondary structure, but only for 2' modified oligoribonucleotides and not for 2'-deoxy oligoribonucleotides.     

Histological aspects of in situ hybridization

Summary This study examined the detection of cellular poly(A) sequences in mouse liver sections by in situ hybridization using a ³H-labelled poly(dT) probe. Parameters examined included possible losses of target poly(A) sequences from sectioned cells, access of probe to target sequences, section thickness, hybridization conditions, autoradioigraphic efficiency, specific activity of probes and specificity of reaction. An improved protocol was devised that resulted in good preservation of histological detail in sectioned tissue blocks, and a calculated hybridization efficiency of 50%–100%. With the use of probes of defined sequence, the protocol should allow detection of unique mRNA sequences within single cells with an estimated sensitivity of 6–12 unique mRNA molecules per sectioned cell.     

A novel endonuclease IV post-PCR genotyping system

Here we describe a novel endonuclease IV (Endo IV) based assay utilizing a substrate that mimics the abasic lesions that normally occur in double-stranded DNA. The three component substrate is characterized by single-stranded DNA target, an oligonucleotide probe, separated from a helper oligonucleotide by a one base gap. The oligonucleotide probe contains a non-fluorescent quencher at the 5′ end and fluorophore attached to the 3′ end through a special rigid linker. Fluorescence of the oligonucleotide probe is efficiently quenched by the interaction of terminal dye and quencher when not hybridized. Upon hybridization of the oligonucleotide probe and helper probe to their complementary target, the phosphodiester linkage between the rigid linker and the 3′ end of the probe is efficiently cleaved, generating a fluorescent signal. In this study, the use of the Endo IV assay as a post-PCR amplification detection system is demonstrated. High sensitivity and specificity are illustrated using single nucleotide polymorphism detection.     

A thermodynamic scale for leucine zipper stability and dimerization specificity: e and g interhelical interactions.

The leucine zipper is a dimeric coiled-coil protein structure composed of two amphipathic alpha-helices with the hydrophobic surfaces interacting to create the dimer interface. This structure has been found to mediate the dimerization of two abundant classes of DNA binding proteins: the bZIP and bHLH-Zip proteins. Several workers have reported that amino acids in the e and g positions of the coiled coil can modulate dimerization stability and specificity. Using the bZIP protein VBP as a host molecule, we report a thermodynamic scale (delta delta G) for 27 interhelical interactions in 35 proteins between amino acids in the g and the following e positions (g<==>e') of a leucine zipper coiled coil. We have examined the four commonly occurring amino acids in the e and g positions of bZIP proteins, lysine (K), arginine (R), glutamine (Q), glutamic acid (E), as well as the only other remaining charged amino acid aspartic acid (D), and finally alanine (A) as a reference amino acid. These results indicate that E<==>R is the most stable interhelical pair, being 0.35 kcal/mol more stable than E<==>K. A thermodynamic cycle analysis shows that the E<==>R pair is 1.33 kcal/mol more stable than A<==>A with -1.14 kcal/mol of coupling energy (delta delta Gint) coming from the interaction of E with R. The E<==>K coupling energy is only -0.14 kcal/mol. E interacts with more specificity than Q. The R<==>R pair is less stable than the K<==>K by 0.24 kcal/mol. R interacts with more specificity than K. Q forms more stable pairs with the basic amino acids K and R rather than with E. Changing amino acids in the e position to A creates bZIP proteins that form tetramers.