Rapid screening of complex DNA samples by single-molecule amplification and sequencing

Microbial cloning makes Sanger sequencing of complex DNA samples possible but is labor intensive. We present a simple, rapid and robust method that enables laboratories without special equipment to perform single-molecule amplicon sequencing, although in a low-throughput manner, from sub-picogram quantities of DNA. The method can also be used for quick quality control of next-generation sequencing libraries, as was demonstrated for a metagenomic sample.     

Global single-cell cDNA amplification to provide a template for representative high-density oligonucleotide microarray analysis

We describe here a protocol for the representative amplification of global mRNAs from typical single mammalian cells to provide a template for high-density oligonucleotide microarray analysis. A single cell is lysed in a tube without purification and first-strand cDNAs are synthesized using a poly(dT)-tailed primer. Unreacted primer is specifically eliminated by exonuclease treatment and second strands are generated with a second poly(dT)-tailed primer after poly(dA) tailing of the first-strand cDNAs. The cDNAs are split into four tubes, which are independently directionally amplified by PCR, and then recombined. The amplified products (approximately 100 ng) show superior representation and reproducibility of original gene expression, especially for genes expressed in more than 20 copies per cell, compared with those obtained by a conventional PCR protocol, and can effectively be used for quantitative PCR and EST analyses. The cDNAs are then subjected to another PCR amplification with primers bearing the T7 promoter sequence. The resultant cDNA products are gel purified, amplified by one final cycle and used for isothermal linear amplification by T7 RNA polymerase to synthesize cRNAs for microarray hybridization. This protocol yields cDNA templates sufficient for more than 80 microarray hybridizations from a single cell, and can be completed in 5-6 days.     

Kinesin's biased stepping mechanism: amplification of neck linker zippering

A physically motivated model of kinesin's motor function is developed within the framework of rectified Brownian motion. The model explains how the amplification of neck linker zippering arises naturally through well-known formulae for overdamped dynamics, thereby providing a means to understand how weakly-favorable zippering leads to strongly favorable plus-directed binding of a free kinesin head to microtubule. Additional aspects of kinesin's motion, such as head coordination and rate-limiting steps, are directly related to the force-dependent inhibition of ATP binding to a microtubule bound head. The model of rectified Brownian motion is presented as an alternative to power stroke models and provides an alternative interpretation for the significance of ATP hydrolysis in the kinesin stepping cycle.     

Thermodynamic and kinetic analysis of sensitivity amplification in biological signal transduction

Based on a thermodynamic analysis of the kinetic model for the protein phosphorylation-dephosphorylation cycle, we study the ATP (or GTP) energy utilization of this ubiquitous biological signal transduction process. It is shown that the free energy from hydrolysis inside cells, DeltaG (phosphorylation potential), controls the amplification and sensitivity of the switch-like cellular module; the response coefficient of the sensitivity amplification approaches the optimal 1 and the Hill coefficient increases with increasing DeltaG. We discover that zero-order ultrasensitivity is mathematically equivalent to allosteric cooperativity. Furthermore, we show that the high amplification in ultrasensitivity is mechanistically related to the proofreading kinetics for protein biosynthesis. Both utilize multiple kinetic cycles in time to gain temporal cooperativity, in contrast to allosteric cooperativity that utilizes multiple subunits in a protein.     

Whole genome amplification: abundant supplies of DNA from precious samples or clinical specimens

Whole genome amplification methods are used to generate the large amounts of DNA that are required for genetic testing. Preparation of genomic DNA from clinical samples is a bottleneck in high-throughput genotyping and is frequently limited by the amount of specimen available. Precious DNA collections used for association and linkage analysis can be a nonrenewable resource available to only a limited number of laboratories. To address these needs, amplified DNA is now being used in a growing number of research and diagnostic applications. A new method, called multiple displacement amplification, dramatically improves the high-fidelity reproduction of genomic DNA, with 10–100 kb amplified DNA products providing uniform coverage of genes.     

Segmentation of complex cell clusters in microscopic images: application to bone marrow samples.

BACKGROUND: Morphologic examination of bone marrow and peripheral blood samples continues to be the cornerstone in diagnostic hematology. In recent years, interest in automatic leukocyte classification using image analysis has increased rapidly. Such systems collect a series of images in which each cell must be segmented accurately to be classified correctly. Although segmentation algorithms have been developed for sparse cells in peripheral blood, the problem of segmenting the complex cell clusters characterizing bone marrow images is harder and has not been addressed previously. METHODS: We present a novel algorithm for segmenting clusters of any number of densely packed cells. The algorithm first oversegments the image into cell subparts. These parts are then assembled into complete cells by solving a combinatorial optimization problem in an efficient way. RESULTS: Our experimental results show that the algorithm succeeds in correctly segmenting densely clustered leukocytes in bone marrow images. CONCLUSIONS: The presented algorithm enables image analysis-based analysis of bone marrow samples for the first time and may also be adopted for other digital cytometric applications where separation of complex cell clusters is required.     

Plant pyruvate dehydrogenase complex: analysis of the kinetic properties and metabolite regulation

The pyruvate dehydrogenase complex was isolated from the mitochondria of broccoli florets and shown to be similar in its reaction mechanism to the complexes from other sources. Three families of parallel lines were obtained for the initial velocity patterns, indicating a multisite ping-pong mechanism. The apparent K_m values obtained were 321 ± 18, 148 ± 13, and 7.2 ± 0.51 μm for pyruvate, NAD⁺, and CoA, respectively. Product inhibition studies using acetyl-CoA and NADH yielded results which were in agreement with those predicted by the multisite ping-pong mechanism. Acetyl-CoA and NADH were found to be competitive inhibitors versus CoA and NAD⁺, respectively. All other substrate-product combinations showed uncompetitive inhibition patterns, except for acetyl-CoA versus NAD⁺. Among various metabolites tested, only hydroxypyruvate (K_i = 0.11 mM) and glyoxylate (K_i = 3.27 mM) were found to be capable of inhibiting the broccoli enzyme to a significant degree. Initial velocity patterns using Mg²⁺? or Ca²⁺-thiamine pyrophosphate and pyruvate as the variable substrate were found to be consistent with an equilibrium ordered mechanism where Mg? or Ca-thiamine pyrophosphate bind first, with dissociation constants of 33.8 and 3 μm, respectively. The Mg- or Ca-thiamine pyrophosphate complexes also dissociated rapidly from the enzyme complex.     

Non-cognate enzyme-DNA complex: structural and kinetic analysis of EcoRV endonuclease bound to the EcoRI recognition site GAATTC

The crystal structure of EcoRV endonuclease bound to non-cognate DNA at 2.0 angstroms resolution shows that very small structural adaptations are sufficient to ensure the extreme sequence specificity characteristic of restriction enzymes. EcoRV bends its specific GATATC site sharply by 50 degrees into the major groove at the center TA step, generating unusual base-base interactions along each individual DNA strand. In the symmetric non-cognate complex bound to GAATTC, the center step bend is relaxed to avoid steric hindrance caused by the different placement of the exocyclic thymine methyl groups. The decreased base-pair unstacking in turn leads to small conformational rearrangements in the sugar-phosphate backbone, sufficient to destabilize binding of crucial divalent metal ions in the active site. A second crystal structure of EcoRV bound to the base-analog GAAUTC site shows that the 50 degrees center-step bend of the DNA is restored. However, while divalent metals bind at high occupancy in this structure, one metal ion shifts away from binding at the scissile DNA phosphate to a position near the 3'-adjacent phosphate group. This may explain why the 10(4)-fold attenuated cleavage efficiency toward GAATTC is reconstituted by less than tenfold toward GAAUTC. Examination of DNA binding and bending by equilibrium and stopped-flow florescence quenching and fluorescence resonance energy transfer (FRET) methods demonstrates that the capacity of EcoRV to bend the GAATTC non-cognate site is severely limited, but that full bending of GAAUTC is achieved at only a threefold reduced rate compared with the cognate complex. Together, the structural and biochemical data demonstrate the existence of distinct mechanisms for ensuring specificity at the bending and catalytic steps, respectively. The limited conformational rearrangements observed in the EcoRV non-cognate complex provide a sharp contrast to the extensive structural changes found in a non-cognate BamHI-DNA crystal structure, thus demonstrating a diversity of mechanisms by which restriction enzymes are able to achieve specificity.     

Linkage analysis of complex diseases using microsatellites and single-nucleotide polymorphisms: application to alcoholism

The efficacy of linkage studies using microsatellites and single-nucleotide polymorphisms (SNPs) was evaluated. Analyzed data were supplied by the Collaborative Study on the Genetics of Alcoholism (COGA). Alcoholism was analyzed together with a simulated trait caused by a gene of known position, through a nonparametric linkage test (NPL). For the alcoholism trait, four densities of SNPs (1 SNP per 0.2 cM, 0.5 cM, 1 cM and 2 cM) showed higher peaks of NPL z scores and smaller significant p-values than the usual 10-cM density of microsatellites. However, the two highest densities of SNPs had unstable z score signals, and therefore were difficult to interpret. Analyzing a simulated trait with the same markers in the same pedigrees, we confirmed the higher power of all four densities of SNPs compared to the 10-cM microsatellites panel, although the existence of other confounding peaks was confirmed for maps that are denser than 1 SNP/cM. We further showed that estimating the gene position using SNPs is far less biased than using the usual panel of microsatellites (biases of 0-2 cM for SNPs vs. 8.9 cM for microsatellites). We conclude that using dense maps of SNPs in linkage analysis is more powerful and less biased than using the 10-cM maps of microsatellites. However, linkage signals can be unstable and difficult to interpret when several SNPs are genotyped per centimorgan. The power and accuracy of 1 SNP/cM or 1 SNP/2 cM may be sufficient in a genome-wide linkage scan while denser maps may be most useful in fine-gene mapping studies exploiting linkage disequilibrium.     

Pyrophosphorolysis-activated polymerization (PAP): application to allele-specific amplification

To measure mutation load or to detect minimal residual disease, a robust method for identifying one mutant allele in the range of 10(6)-10(9) wild-type alleles would be advantageous. Herein, we present evidence that pyrophosphorolysis-activated polymerization (PAP) has the potential to provide a highly specific and robust method of allele-specific amplification if DNA polymerases with higher pyrophosphorolysis activity can be found or engineered. In PAP, pyrophosphorolysis and polymerization by DNA polymerase are coupled serially by utilizing a pyrophosphorolysis-activatable oligonucleotide (P*). P*, which is an allele-specific oligonucleotide with a dideoxynucleotide at the 3' terminus, can be activated by pyrophosphorolysis to remove the 3' terminal dideoxynucleotide in the presence of pyrophosphate (PPi) and the complementary strand of the allelic template; then the activated P* can be extended by DNA polymerization. Specificity results from both pyrophosphorolysis and polymerization because significant nonspecific amplification requires the combination of mismatch pyrophosphorolysis and misincorporation by the DNA polymerase, which is an extremely rare event. Proof of principle has been achieved with a polymorphic site within the human D1 dopamine receptor gene. The effects of the dideoxyoligonucleotide sequences, DNA polymerases, PPi concentrations, allele-specific templates, pH and dNTP concentrations were examined.     

gamma-Glutamylcyclotransferase: inhibition by D-beta-aminoglutaryl-L-alanine and analysis of the solvent kinetic isotope effect

Spin-echo NMR spectroscopy was used to record the cleavage of a gamma-glutamyl--amino-acid by (5-L-glutamyl)-L-amino-acid 5-glutamyltransferase (cyclizing) (gamma-glutamylcyclotransferase) in human erythrocyte hemolysates. The Michaelis-Menten steady-state kinetic parameters were obtained by fitting the integrated Michaelis-Menten equation to the reaction time curves. The product, L-5-oxoproline, was shown to be an inhibitor of the reaction. The active site of the enzyme was probed by studies of the inhibition by D- and L-beta-aminoglutaryl-L-alanine which are the beta-amino-acid isomers of D- and L-gamma-glutamyl-L-alanine (the latter being a natural substrate of the enzyme); the D-isomer was the more potent inhibitor (Ki = 0.30 +/- 0.02 mmol/l water). When the alanyl alpha-carboxyl of the inhibitor was reduced to a hydroxyl (i.e. to give D-beta-aminoglutaryl-L-alaninol) the potency of inhibition was reduced. The previously reported kinetic isotope effect of solvent 2H2O on the enzyme-catalyzed reaction has been further studied using a proton inventory. We propose that the solvent kinetic isotope effect is due to an intramolecular proton transfer between the glutamyl amino group and the peptide bond nitrogen.     

An effective method to extract DNA from environmental samples for polymerase chain reaction amplification and DNA fingerprint analysis

A rapid direct-extraction method was used to obtain DNA from environmental soil samples. Heat, enzymes, and guanidine isothiocyanate were utilized to lyse cells. The DNA was purified by agarose gel electrophoresis, amplified with 16S rRNA-based primers by use of the polymerase chain reaction, and then digested with the restriction endonucleasePalI. The extraction method was used to obtain DNA from a variety of plants, bacteria, and fungi includingGossypium hirsucum (cotton),Pseudomonas, Bacillus, Streptomyces, andColletotrichum. Up to 100 g DNA/g (wet weight) of soil and 400 g DNA/g of plant material were recovered. Restriction endonuclease analysis patterns of amplified rDNA from pure microbial cultures and plant species contained three to five different DNA fragments. Amplified rDNA of mixed population DNA extracts from soil samples, digested with the restriction endonucleasePalI, contained 12–20 DNA fragments, appearing as sample fingerprints. Results from eight environmental soil samples that were analyzed suggest that the amplified rDNA fingerprints can be used to help characterize the genetic and biological diversity of the microbial populations in these samples.     

Use of whole genome amplification to rescue DNA from plasma samples

While DNA of good quality and sufficient amount can be obtained easily from whole blood, buccal swabs, surgical specimens, or cell lines, these DNA-rich sources are not always available. This is particularly the case in studies for which biological specimens were collected when genotyping assays were not widely available. In those studies, serum or plasma is often the only source of DNA. Newly developed whole genome amplification (WGA) methods, based on phi29 polymerase, may play a significant role in recovering DNA in such instances. We tested a total of 528 plasma samples kept in storage at -40 degrees C for approximately 10 years for 8 single nucleotide polymorphisms (SNPs) using the 5' exonuclease (TaqMan) assay. These specimens yielded undetectable levels of DNA following extraction with an affinity column but produced an average 52.7 microg (standard deviation of 31.2 microg) of DNA when column-extracted DNA was used as a template for WGA. This increased the genotyping success rate from 54% to 93%. There were only 3 disagreements out of 364 paired genotyping results for pre- and post-WGA DNAs, indicating an error rate of 0.82%. These results are encouraging for expanding the use of poor DNA resources in genotyping studies.     

Recovering motifs from biased genomes: application of signal correction

A significant problem in biological motif analysis arises when the background symbol distribution is biased (e.g. high/low GC content in the case of DNA sequences). This can lead to overestimation of the amount of information encoded in a motif. A motif can be depicted as a signal using information theory (IT). We apply two concepts from IT, distortion and patterned interference (a type of noise), to model genomic and codon bias respectively. This modeling approach allows us to correct a raw signal to recover signals that are weakened by compositional bias. The corrected signal is more likely to be discriminated from a biased background by a macromolecule. We apply this correction technique to recover ribosome-binding site (RBS) signals from available sequenced and annotated prokaryotic genomes having diverse compositional biases. We observed that linear correction was sufficient for recovering signals even at the extremes of these biases. Further comparative genomics studies were made possible upon correction of these signals. We find that the average Euclidian distance between RBS signal frequency matrices of different genomes can be significantly reduced by using the correction technique. Within this reduced average distance, we can find examples of class-specific RBS signals. Our results have implications for motif-based prediction, particularly with regards to the estimation of reliable inter-genomic model parameters.     

Efficient procedure for purification of obligate intracellular Wolbachia pipientis and representative amplification of its genome by multiple-displacement amplification

Bacteria belonging to the genus Wolbachia are obligatory microendocytobionts that infect a variety of arthropods and a majority of filarial nematode species, where they induce reproductive alterations or establish a mutualistic symbiosis. Although two whole genome sequences of Wolbachia pipientis, for strain wMel from Drosophila melanogaster and strain wBm from Brugia malayi, have been fully completed and six other genome sequencing projects are ongoing (http://www.genomesonline.org/index.cgi?want=Prokaryotic+Ongoin), genetic analyses of these bacteria are still scarce, mainly due to the inability to cultivate them outside of eukaryotic cells. Usually, a large amount of host tissue (a thousand individuals, or about 10 g) is required in order to purify Wolbachia and extract its DNA, which is often recovered in small amounts and contaminated by host cell DNA, thus hindering genomic studies. In this report, we describe an efficient and reliable procedure to representatively amplify the Wolbachia genome by multiple-displacement amplification from limited infected host tissue (0.2 g or 2 x 10(7) cells). We obtained sufficient amounts (8 to 10 microg) of DNA of suitable quality for genomic studies, and we demonstrated that the amplified DNA contained all of the Wolbachia loci targeted. In addition, our data indicated that the genome of strain wRi, an obligatory endosymbiont of Drosophila simulans, shares a similar overall architecture with its relative strain wMel.