Quantification of PCR bias caused by a single nucleotide polymorphism in SMN gene dosage analysis

Multi-template polymerase chain reaction

PCR is a formidable and potent technology that serves as an indispensable tool in a wide range of biological disciplines. However, due to the ease of use and often lack of rigorous standards many PCR applications can lead to highly variable, inaccurate, and ultimately meaningless results. Thus, rigorous method validation must precede its broad adoption to any new application. Multi-template samples possess particular features, which make their PCR analysis prone to artifacts and biases: multiple homologous templates present in copy numbers that vary within several orders of magnitude. Such conditions are a breeding ground for chimeras and heteroduplexes. Differences in template amplification efficiencies and template competition for reaction compounds undermine correct preservation of the original template ratio. In addition, the presence of inhibitors aggravates all of the above-mentioned problems. Inhibitors might also have ambivalent effects on the different templates within the same sample. Yet, no standard approaches exist for monitoring inhibitory effects in multitemplate PCR, which is crucial for establishing compatibility between samples.

Length and GC-biases during sequencing library amplification: a comparison of various polymerase-buffer systems with ancient and modern DNA sequencing libraries

High-throughput sequencing technologies frequently necessitate the use of PCR for sequencing library amplification. PCR is a sometimes enigmatic process and is known to introduce biases. Here we perform a simple amplification-sequencing assay using 10 commercially available polymerase-buffer systems to amplify libraries prepared from both modern and ancient DNA. We compare the performance of the polymerases with respect to a previously uncharacterized template length bias, as well as GC-content bias, and find that simply avoiding certain polymerase can dramatically decrease the occurrence of both. For amplification of ancient DNA, we found that some commonly used polymerases strongly bias against amplification of endogenous DNA in favor of GC-rich microbial contamination, in our case reducing the fraction of endogenous sequences to almost half.

Nested PCR bias: a case study of Pseudomonas spp. in soil microcosms

Nested PCR methods combined with denaturing gradient gel electrophoresis (DGGE) are widely used for the detection of low copy number species or for the analysis of group-specific community profiles. With an appropriate number of PCR cycles during the first round of amplification, initial differences in the copy numbers of different DNA fragments that are targeted can be maintained during the second round without significant bias. However, if an excessive number of cycles in used in the first round, relative differences in the copy numbers of the targeted sequences can be obscured. Here we demonstrate the effect of "nested PCR bias" in a case study with PCR-DGGE of 16S rRNA gene sequences targeting Pseudomonas spp. following exposure of soil to naphthalene vapors. Our results demonstrate artifacts caused by nested PCR bias can be substantially minimized by calibrating the number of first round PCR cycles, thereby preserving the ability to obtain semiquantitative data for evaluating changes in gene copy numbers over time.

Amplification-free Illumina sequencing-library preparation facilitates improved mapping and assembly of (G+C)-biased genomes

Amplification artifacts introduced during library preparation for the Illumina Genome Analyzer increase the likelihood that an appreciable proportion of these sequences will be duplicates, and cause an uneven distribution of read coverage across the targeted sequencing regions. As a consequence, these unfavorable features result in difficulties in genome assembly and variation analysis from the short reads, particularly when the sequences are from genomes with base compositions at the extremes of high or low GC content. Here we present an amplification-free method of library preparation, in which the cluster amplification step, rather than the polymerase chain reaction, enriches for fully ligated template strands, reducing the incidence of duplicate sequences, improving read mapping and SNP calling and aiding de novo assembly. We illustrate this by generating and analysing DNA sequences from extremely GC-poor (Plasmodium falciparum), GC-neutral (Escherichia coli) and high GC (Bordetella pertussis) genomes.

Occurrence of genetic bottlenecks during citrus tristeza virus acquisition by Toxoptera citricida under field conditions

In this study, we address the involvement of T. citricida in strain segregation and genetic bottleneck events by comparing the nucleotide diversity of CTV coat protein (CP) gene variants present in field-grown trees with that of variants retrieved from single apterous aphids. Plant material and aphids were collected in orange orchards in the northern part of Portugal. Shoots from two trees that were found to be positive using ELISA and twenty-four apterous aphids from these same trees were selected for individual molecular assays. CTV was detected in 60% of the aphids by amplification of a 417-bp fragment of the CP gene. Analysis of molecular variance (AMOVA) of this fragment revealed that most of the variation of the virus was found among individual aphids (FSC: 0.766) within each location. Nucleotide diversity comparison between the pool of sequences obtained from a given shoot and sequences obtained from individual aphids present on that shoot showed a reduction of more than one order of magnitude in most cases. Computer simulations of random virus acquisition by single aphids showed that in 54% of the cases only a single CP gene phylogenetic group was acquired. However, a small number of aphids (e.g. 6) was enough to acquire the full complement of phylogenetic groups present.

Robust quantification of the SMN gene copy number by real-time TaqMan PCR

Spinal muscular atrophy (SMA) is an autosomal recessive disease caused by mutation or deletion of the survival motor neuron gene 1 (SMN1). The highly homologous gene, SMN2, is present in all patients, but it cannot compensate for loss of SMN1. SMN2 differs from SMN1 by a few nucleotide changes, but a C --> T transition in exon 7 leads to exon skipping. As a result, most transcripts from the SMN2 gene lack exon 7. Although SMN1 is the disease-determining gene, the number of SMN2 copies appears to modulate SMA clinical phenotypes. Thus, determining the SMN copy number is important for clinical diagnosis and prognosis. We have developed a quantitative real-time TaqMan polymerase chain reaction assay for both the SMN1 and SMN2 genes, in which reliable copy number determination was possible on deoxyribonucleic acid samples obtained by two different isolation methods and from two different sources (human blood and skin fibroblasts). For SMN1, allele specificity was attained solely by addition of an allele-specific forward primer and, for SMN2, by addition of a specific forward primer and a nonextending oligonucleotide (SMN1 blocker) that reduced nonspecific amplification from SMN1 to a negligible level. We validated the reliability of this real-time polymerase chain reaction approach and found that the coefficient of variation for all the gene copy number measurements was below 10%. Quantitative analysis of the SMN copy number in SMA fibroblasts by this approach showed deletion of SMN1 and an inverse correlation between the SMN2 copy number and severity of the disease.

Harnessing asymmetrical substrate recognition by thermostable EndoV to achieve balanced linear amplification in multiplexed SNP typing

Multiplexed amplification of specific DNA sequences, by PCR or by strand-displacement amplification, is an intrinsically biased process. The relative abundance of amplified DNA can be altered significantly from the original representation and, in extreme cases, allele dropout can occur. In this paper, we present a method of linear amplification of DNA that relies on the cooperative, sequence-dependent functioning of the DNA mismatch-repair enzyme endonuclease V (EndoV) from Thermotoga maritima (Tma) and Bacillus stearothermophilus (Bst) DNA polymerase. Tma EndoV can nick one strand of unmodified duplex DNA, allowing extension by Bst polymerase. By controlling the bases surrounding a mismatch and the mismatch itself, the efficiency of nicking by EndoV and extension by Bst polymerase can be controlled. The method currently allows 100-fold multiplexed amplification of target molecules to be performed isothermally, with an average change of <1.3-fold in their original representation. Because only a single primer is necessary, primer artefacts and nonspecific amplification products are minimized.

Precision and performance characteristics of bisulfite conversion and real-time PCR (MethyLight) for quantitative DNA methylation analysis

Assays to measure DNA methylation, which are important in epigenetic research and clinical diagnostics, typically rely on conversion of unmethylated cytosine to uracil by sodium bisulfite. However, no study has comprehensively evaluated the precision and performance characteristics of sodium bisulfite conversion and subsequent quantitative methylation assay. We developed quantitative real-time polymerase chain reaction (MethyLight) to measure percentage of methylated reference (PMR, ie, degree of methylation) for the MGMT, MLH1, and CDKN2A (p16) promoters. To measure the precision of bisulfite conversion, we bisulfite-treated seven different aliquots of DNA from each of four paraffin-embedded colon cancer samples. To assess run-to-run variation, we repeated MethyLight five times. Bisulfite-to-bisulfite coefficient of variation (CV) of PMR ranged from 0.10 to 0.38 (mean, 0.21), and run-to-run CV of PMR ranged from 0.046 to 0.60 (mean, 0.31). Interclass correlation coefficients were 0.74 to 0.84 for the three loci, indicating good reproducibility. DNA mixing study with methylated and unmethylated DNA showed good linearity of the assay. Of 272 colorectal cancers evaluated, most showed PMR either <1 or >10, and promoter methylation (PMR >4) was tightly associated with loss of respective protein expression (P < 10(-16)). In conclusion, sodium bisulfite conversion and quantitative MethyLight assays have good precision and linearity and can be effectively used for high-throughput DNA methylation analysis on paraffin-embedded tissue.

Sensitive sequencing method for KRAS mutation detection by Pyrosequencing

Both benign and malignant tumors represent heterogenous tissue containing tumor cells and non-neoplastic mesenchymal and inflammatory cells. To detect a minority of mutant KRAS alleles among abundant wild-type alleles, we developed a sensitive DNA sequencing assay using Pyrosequencing, ie, nucleotide extension sequencing with an allele quantification capability. We designed our Pyrosequencing assay for use with whole-genome-amplified DNA from paraffin-embedded tissue. Assessing various mixtures of DNA from mutant KRAS cell lines and DNA from a wild-type KRAS cell line, we found that mutation detection rates for Pyrosequencing were superior to dideoxy sequencing. In addition, Pyrosequencing proved superior to dideoxy sequencing in the detection of KRAS mutations from DNA mixtures of paraffin-embedded colon cancer and normal tissue as well as from paraffin-embedded pancreatic cancers. Quantification of mutant alleles by Pyrosequencing was precise and useful for assay validation, monitoring, and quality assurance. Our Pyrosequencing method is simple, robust, and sensitive, with a detection limit of approximately 5% mutant alleles. It is particularly useful for tumors containing abundant non-neoplastic cells. In addition, the applicability of this assay for DNA amplified by whole-genome amplification technique provides an expanded source of DNA for large-scale studies.

Evaluation of temperature gradient capillary electrophoresis for detection of the Factor V Leiden mutation: coincident identification of a novel polymorphism in Factor V

AIM

The Factor V Leiden mutation (G1691A) is a clinically important polymorphism that results in an increased risk of thrombosis. The goal of this study was to compare a temperature gradient capillary electrophoresis (TGCE) platform for the detection of Factor V gene mutations to a conventional restriction fragment length polymorphism (RFLP) assay.

METHODS

Three hundred and four samples were analyzed by both TGCE and a common clinical Mnl I/RFLP assay. Concordance of results between the two assays was observed for 302/304 (99.3%) of the samples.

RESULTS

All of the Leiden mutants (23/23, 100%) were identified by TGCE. Of the two discrepant results, one was caused by low peak heights in the TGCE output data and was easily rectified by the addition of a minimum peak height threshold. The second discrepancy resulted from the presence of a G-->A transition 95 bp downstream of the Leiden mutation site. This polymorphism represents a previously unreported alteration of the Factor V gene.

CONCLUSIONS

The TGCE assay is less labor-intensive and has a higher throughput capacity than the Mnl I/RFLP assay. TGCE is a less specific assay than the Mnl I/RFLP assay that allows for the detection of novel polymorphisms, but also creates the need for all positive TGCE results to be confirmed by an alternate method such as sequencing. Our results demonstrate that TGCE is a highly sensitive method for mutation detection and has utility for mutation discovery analysis.

DFold: PCR design that minimizes secondary structure and optimizes downstream genotyping applications

Secondary structures in polymerase chain reaction (PCR) target sequences have a negative impact on amplification success rates and on downstream uses of PCR products. For example, signal strength and allele discrimination in single nucleotide polymorphism (SNP) genotyping methods can be compromised by allele-biased amplification and/or by PCR product folding that limits access of interrogating probes. To increase the fidelity and robustness of PCR, and to aid follow-on applications, we have developed DFold (http://dfold.cgb.ki.se)-a generalized software solution that creates PCR oligonucleotide primer designs devoid of stable secondary structures. We demonstrate the effectiveness of the tool by applying it to a range of dynamic allele-specific hybridization (DASH) assay designs, many of which we evaluate in the laboratory. We further consider how the system throughput may be made sufficiently high for use upon millions of target sequences in order to support whole-genome analyses.

New insights on the evolution of the SMN1 and SMN2 region: simulation and meta-analysis for allele and haplotype frequency calculations

Most spinal muscular atrophy patients lack both copies of SMN1. Loss of SMN1 ('0-copy alleles') can occur by gene deletion or SMN1-to-SMN2 gene conversion. Despite worldwide efforts to map the segmental duplications within the SMN region, most assemblies do not correctly delineate these genes. A near pericentromeric location provides impetus for the strong evidence that SMN1 and SMN2 arose from a primate-specific paralogous gene duplication. Here we meta-analyzed our recent laboratory results together with available published data, in order to calculate new mutation rates and allele/haplotype frequencies in this recalcitrant and highly unstable region of the human genome. Based on our tested assumption of compliance with Hardy-Weinberg equilibrium, we conclude that the SMN1 allele frequencies are: '0-copy disease alleles,' 0.013; '1-copy normal alleles,' 0.95; '2-copy normal alleles (ie, two copies of SMN1 on one chromosome),' 0.038; and '1(D) disease alleles (SMN1 with a small intragenic mutation),' 0.00024. The SMN1 haplotype ['(SMN1 copy number)-(SMN2 copy number)'] frequencies are: '0-0,' 0.00048; '0-1,' 0.0086; '0-2,' 0.0042; '1-0,' 0.27; '1-1,' 0.66; '1-2,' 0.015; '2-0,' 0.027; and '2-1,' 0.012. Paternal and maternal de novo mutation rates are 2.1 x 10(-4) and 4.2 x 10(-5), respectively. Our data provide the basis for the most accurate genetic risk calculations, as well as new insights on the evolution of the SMN region, with evidence that nucleotide position 840 (where a transition 840C>T functionally distinguishes SMN2 from SMN1) constitutes a mutation hotspot. Our data also suggest selection of the 1-1 haplotype and the presence of rare chromosomes with three copies of SMN1.

Rapid detection of single nucleotide polymorphisms associated with spinal muscular atrophy by use of a reusable fibre-optic biosensor

Rapid (<2 min) and quantitative genotyping for single nucleotide polymorphisms (SNPs) associated with spinal muscular atrophy was done using a reusable (approximately 80 cycles of application) fibre-optic biosensor over a clinically relevant range (0-4 gene copies). Sensors were functionalized with oligonucleotide probes immobilized at high density (approximately 7 pmol/cm2) to impart enhanced selectivity for SNP discrimination and used in a total internal reflection fluorescence detection motif to detect 202 bp PCR amplicons from patient samples. Real-time detection may be done over a range of ionic strength conditions (0.1-1.0 M) without stringency rinsing to remove non-selectively bound materials and without loss of selectivity, permitting a means for facile sample preparation. By using the time-derivative of fluorescence intensity as the analytical parameter, linearity of response may be maintained while allowing for significant reductions in analysis time (10-100-fold), permitting for the completion of measurements in under 1 min.

Inverse correlation between SMN1 and SMN2 copy numbers: evidence for gene conversion from SMN2 to SMN1

Most carriers of autosomal recessive spinal muscular atrophy (SMA) have only one copy of SMN1 because of SMN1 gene deletions or gene conversions from SMN1 to SMN2, which has only one base difference in coding sequence from SMN1. Using SMN gene dosage analysis, we determined the copy numbers of SMN1 and SMN2 in the general population as well as in SMA patients and carriers. Increased SMN1 copy number is associated with decreased SMN2 copy number in the general population; that is, SMN2 copy number was decreased to one or zero copies in 11 of 13 individuals with three or four copies of SMN1, whereas only 71 of 164 individuals with two copies of SMN1 had one or zero copies of SMN2 (P<0.01). SMN2 copy number was increased to three or four in a subset of SMN1 deletion/conversion carriers, and in most SMA patients with a milder phenotype. In conclusion, our data provide evidence that gene conversion from SMN2 to SMN1 occurs, and that SMN1 converted from SMN2 is present in the general population.